Fast Radio Bursts (FRBs) are millisecond-long radio pulses of extragalactic origin with peak luminosities far exceeding any Milky Way sources. The prevalent invocation for the FRB origin involves magnetars: young, magnetically powered neutron stars with the strongest magnetic fields in the Universe. A magnetar-defining signature is the emission of bright, hard X-ray bursts of sub-second duration. These occur in isolation or during a burst storm, when several hundred are observed within minutes to hours. On April 27th 2020, the Galactic magnetar SGR J1935+2154 entered an active period, emitting hundreds of X-ray bursts in a few hours. Remarkably, only one of these temporally coincided with an FRB-like radio burst. Here we report on the spectral and temporal analyses of 24 X-ray bursts emitted 13 hours prior to the FRB and seen simultaneously with NASA NICER and Fermi/GBM missions in their combined energy range. We demonstrate that the FRB-associated X-ray burst is very similar temporally, albeit strikingly different spectrally, from the 24 NICER/GBM bursts. If the FRB-associated burst were drawn from this magnetar burst population, its occurrence rate would be at most around 1 in 7000. This rarity combined with the unusual X-ray burst spectrum is perhaps indicative of an uncommon locale for the origin of the FRB-associated burst. We suggest that this unique event originated in quasi-polar open or closed magnetic field lines extending to high altitudes where radio emission can be generated, possibly from a collimated plasma flow.

We report on the timing and spectral properties of the soft X-ray emission from the magnetar 1E 2259+586 from January 2013, $\sim 8$ months after the detection of an anti-glitch, until September 2019, using the Neil Gehrels Swift and NICER observatories. During this time span, we detect two timing discontinuities. The first, occurring around 5 years after the April 2012 anti-glitch, is a relatively large spin-up glitch with a fractional amplitude $\Delta\nu/\nu=1.24(2)\times10^{-6}$. We find no evidence for flux enhancement or change in the spectral or pulse profile shape around the time of this glitch. This is consistent with the picture that a significant number of magnetar spin-up glitches are radiatively-quiet. Approximately 1.5 years later in April 2019, 1E 2259+586 exhibited an anti-glitch with spin-down of a fractional amplitude $\Delta\nu/\nu=-5.8(1)\times10^{-7}$; similar to the fractional change detected in 2012. We do not, however, detect any change to the pulse-profile shape or increase in the rms pulsed flux of the source, nor do we see any possible bursts from its direction around the time of the anti-glitch; all of which occurred during the 2012 event. Hence, similar to spin-up glitches, anti-glitches can occur silently. This may suggest that these phenomena originate in the neutron star interior, and that their locale and triggering mechanism do not necessarily have to be connected to the magnetosphere. Lastly, our observations suggest that the occurrence rate of spin-up and spin-down glitches is about the same in 1E 2259+586, with the former having a larger net fractional change.

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.

Magnetars are young, rotating neutron stars that possess larger magnetic fields ($B$ $\approx$ $10^{13}$-$10^{15}$ G) and longer rotational periods ($P$ $\approx$ 1-12 s) than ordinary pulsars. In contrast to rotation-powered pulsars, magnetar emission is thought to be fueled by the evolution and decay of their powerful magnetic fields. They display highly variable radio and X-ray emission, but the processes responsible for this behavior remain a mystery. We report the discovery of bright, persistent individual X-ray pulses from XTE J1810-197, a transient radio magnetar, using the Neutron star Interior Composition Explorer (NICER) following its recent radio reactivation. Similar behavior has only been previously observed from a magnetar during short time periods following a giant flare. However, the X-ray pulses presented here were detected outside of a flaring state. They are less energetic and display temporal structure that differs from the impulsive X-ray events previously observed from the magnetar class, such as giant flares and short X-ray bursts. Our high frequency radio observations of the magnetar, carried out simultaneously with the X-ray observations, demonstrate that the relative alignment between the X-ray and radio pulses varies on rotational timescales. No correlation was found between the amplitudes or temporal structure of the X-ray and radio pulses. The magnetar's 8.3 GHz radio pulses displayed frequency structure, which was not observed in the pulses detected simultaneously at 31.9 GHz. Many of the radio pulses were also not detected simultaneously at both frequencies, which indicates that the underlying emission mechanism producing these pulses is not broadband. We find that the radio pulses from XTE J1810-197 share similar characteristics to radio bursts detected from fast radio burst (FRB) sources, some of which are now thought to be produced by active magnetars.

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.

In November 2019, MAXI detected an X-ray outburst from the known Be X-ray binary system RX J0209.6-7427 located in the outer wing of the Small Magellanic Cloud. We followed the outburst of the system with NICER which led to the discovery of X-ray pulsations with a period of 9.3 s. We analyzed simultaneous X-ray data obtained with NuSTAR and NICER allowing us to characterize the spectrum and provide an accurate estimate of its bolometric luminosity. During the outburst the maximum broadband X-ray luminosity of the system reached $1-2\times10^{39}$ erg/s, thus exceeding by about one order of magnitude the Eddington limit for a typical 1.4 $M_{\odot}$ mass neutron star (NS). Monitoring observations with Fermi/GBM and NICER allowed us to study the spin evolution of the NS and compare it with standard accretion torque models. We found that the NS magnetic field should be of the order of $3\times10^{12}$ G. We conclude that RX J0209.6-7427 exhibited one of the brightest outbursts observed from a Be X-ray binary pulsar in the Magellanic Clouds, reaching similar luminosity level to the 2016 outburst of SMC X-3. Despite the super-Eddington luminosity of RX J0209.6-7427, the NS appears to have only a moderate magnetic field strength.

We search for extrasolar planets around millisecond pulsars using pulsar timing data and seek to determine the minimum detectable planetary masses as a function of orbital period. Using the 11-year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), we look for variations from our models of pulse arrival times due to the presence of exoplanets. No planets are detected around the millisecond pulsars in the NANOGrav 11-year data set, but taking into consideration the noise levels of each pulsar and the sampling rate of our observations, we develop limits that show we are sensitive to planetary masses as low as that of the moon. We analyzed potential planet periods, P, in the range 7 days < P < 2000 days, with somewhat smaller ranges for some binary pulsars. The planetary mass limit for our median-sensitivity pulsar within this period range is 1 M_moon (P / 100 days)^(-2/3).

We report on the phase-coherent timing analysis of the accreting millisecond X-ray pulsar IGR J17591-2342, using Neutron Star Interior Composition Explorer (NICER) data taken during the outburst of the source between 2018 August 15 and 2018 October 17. We obtain an updated orbital solution of the binary system. We investigate the evolution of the neutron star spin frequency during the outburst, reporting a refined estimate of the spin frequency and the first estimate of the spin frequency derivative ($\dot{\nu} \sim -7\times 10^{-14}$ Hz s$^{-1}$), confirmed independently from the modelling of the fundamental frequency and its first harmonic. We further investigate the evolution of the X-ray pulse phases adopting a physical model that accounts for the accretion material torque as well as the magnetic threading of the accretion disc in regions where the Keplerian velocity is slower than the magnetosphere velocity. From this analysis we estimate the neutron star magnetic field $B_{eq} = 2.8(3)\times10^{8}$ G. Finally, we investigate the pulse profile dependence on energy finding that the observed behaviour of the pulse fractional amplitude and lags as a function of energy are compatible with a thermal Comptonisation of the soft photons emitted from the neutron star caps.

Astrophysics spans an enormous range of questions on scales from individual planets to the entire cosmos. To address the richness of 21st century astrophysics requires a corresponding richness of telescopes spanning all bands and all messengers. Much scientific benefit comes from having the multi-wavelength capability available at the same time. Most of these bands,or measurement sensitivities, require space-based missions. Historically, NASA has addressed this need for breadth with a small number of flagship-class missions and a larger number of Explorer missions. While the Explorer program continues to flourish, there is a large gap between Explorers and strategic missions. A fortunate combination of new astrophysics technologies with new, high capacity, low dollar-per-kg to orbit launchers, and new satellite buses allow for cheaper missions with capabilities approaching strategic mission levels. NASA has recognized these developments by calling for Probe-class mission ideas for mission studies, spanning most of the electromagnetic spectrum from GeV gamma-rays to the far infrared, and the new messengers of neutrinos and ultra-high energy cosmic rays. The key insight from the Probes exercise is that order-of-magnitude advances in science performance metrics are possible across the board for initial total cost estimates in the range 500M-1B dollars.

We report the discovery of PSR~J2251$-$3711, a radio pulsar with a spin period of 12.1 seconds, the second longest currently known. Its timing parameters imply a characteristic age of 15 Myr, a surface magnetic field of $1.3 \times 10^{13}$~G and a spin-down luminosity of $2.9 \times 10^{29}~\mathrm{erg~s}^{-1}$. Its dispersion measure of 12.12(1)~$\mathrm{pc}~\mathrm{cm}^{-3}$ leads to distance estimates of 0.5 and 1.3 kpc according to the NE2001 and YMW16 Galactic free electron density models, respectively. Some of its single pulses show an uninterrupted 180 degree sweep of the phase-resolved polarization position angle, with an S-shape reminiscent of the rotating vector model prediction. However, the fact that this sweep occurs at different phases from one pulse to another is remarkable and without straightforward explanation. Although PSR~J2251$-$3711 lies in the region of the $P-\dot{P}$ parameter space occupied by the X-ray Isolated Neutron Stars (XINS), there is no evidence for an X-ray counterpart in our Swift XRT observation; this places a 99\%-confidence upper bound on its unabsorbed bolometric thermal luminosity of $1.1 \times 10^{31}~(d / 1~\mathrm{kpc})^2~\mathrm{erg/s}$ for an assumed temperature of 85 eV, where $d$ is the distance to the pulsar. Further observations are needed to determine whether it is a rotation-powered pulsar with a true age of at least several Myr, or a much younger object such as an XINS or a recently cooled magnetar. Extreme specimens like PSR J2251$-$3711 help bridge populations in the so-called neutron star zoo in an attempt to understand their origins and evolution.

We present two years of Neutron star Interior Composition Explorer (NICER) X-ray observations of three energetic rotation-powered millisecond pulsars (MSPs): PSRs B1937+21, B1821-24, and J0218+4232. We fit Gaussians and Lorentzians to the pulse profiles for different energy sub-bands of the soft X-ray regime to measure the energy dependence of pulse separation and width. We find that the separation between pulse components of PSR J0218+4232 decreases with increasing energy at $\gt 3\sigma$ confidence. The 95% upper limit on pulse separation evolution for PSRs B1937+21 and B1821-24 is less than 2 milliperiods per keV. Our phase-resolved spectral results provide updated constraints on the non-thermal X-ray emission of these three pulsars. The photon indices of the modeled X-ray emission spectra for each pulse component of PSR B1937+21 are inconsistent with each other at the 90% confidence level, suggesting different emission origins for each pulse. We find that the PSR B1821-24 and PSR J0218+4232 emission spectra are invariant with phase at the 90% confidence level. We describe the implications of our profile and spectral results in the context of equatorial current sheet emission models for these three MSPs with non-thermal, magnetospheric X-ray emission.

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.

We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030$+$0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer (NICER) X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply-connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass $M$ and equatorial radius $R_\mathrm{eq}$ are, respectively, $1.34_{-0.16}^{+0.15}$ M$_{\odot}$ and $12.71_{-1.19}^{+1.14}$ km, whilst the compactness $GM/R_\mathrm{eq}c^2 = 0.156_{-0.010}^{+0.008}$ is more tightly constrained; the credible interval bounds reported here are approximately the $16\%$ and $84\%$ quantiles in marginal posterior mass.

NICER observed several rotation-powered millisecond pulsars to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections ($\ge4.7\sigma$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sample: PSR J0614-3329, PSR J0636+5129, PSR J0751+1807, PSR J1012+5307, and PSR J2241-5236, while PSR J1552+5437 and PSR J1744-1134 remain undetected. Of those, only PSR J0751+1807 and PSR J1012+5307 had pulsations previously detected at the 1.7$\sigma$ and almost 3$\sigma$ confidence levels, respectively, in XMM-Newton data. All detected sources exhibit broad sine-like pulses, which are indicative of surface thermal radiation. As such, these MSPs are promising targets for future X-ray observations aimed at constraining the neutron star mass-radius relation and the dense matter equation of state using detailed pulse profile modeling. Furthermore, we find that three of the detected millisecond pulsars exhibit a significant phase offset between their X-ray and radio pulses.

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.

We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the instrument, observations, and data processing/reduction procedures, as well as the series of investigations conducted to ensure that the properties of the data sets are suitable for parameter estimation analyses to produce reliable constraints on the neutron star mass-radius relation and the dense matter equation of state. We find that the long-term timing and flux behavior and the Fourier-domain properties of the event data do not exhibit any anomalies that could adversely affect the intended measurements. From phase-selected spectroscopy, we find that emission from the individual pulse peaks is well described by a single-temperature hydrogen atmosphere spectrum, with the exception of PSR J0437-4715, for which multiple temperatures are required.

Recent modeling of NICER observations of thermal X-ray pulsations from the surface of the isolated millisecond pulsar PSR J0030+0451 suggests that the hot emitting regions on the pulsar's surface are far from antipodal, which is at odds with the classical assumption that the magnetic field in the pulsar magnetosphere is predominantly that of a centered dipole. Here, we review these results and examine previous attempts to constrain the magnetospheric configuration of PSR J0030+0451. To the best of our knowledge, there is in fact no direct observational evidence that PSR J0030+0451's magnetic field is a centered dipole. Developing models of physically motivated, non-canonical magnetic field configurations and the currents that they can support poses a challenging task. However, such models may have profound implications for many aspects of pulsar research, including pulsar braking, estimates of birth velocities, and interpretations of multi-wavelength magnetospheric emission.

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.

In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux, and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth's orbital period, enabling pulsar timing based astrometry of this accreting millisecond pulsar.

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.

The Low-Frequency Array radio telescope discovered the $707$ Hz binary millisecond pulsar (MSP) J0952$-$0607 in a targeted radio pulsation search of an unidentified $\textit{Fermi}$ gamma-ray source. This source shows a weak energy flux of $F_\gamma = 2.6 \times 10^{-12}\,\text{erg}\,\text{cm}^{-2}\,\text{s}^{-1}$ in the energy range between $100\,\text{MeV}$ and $100\,\text{GeV}$. Here we report the detection of pulsed gamma-ray emission from PSR$\,$J0952$-$0607 in a very sensitive gamma-ray pulsation search. The pulsar's rotational, binary, and astrometric properties are measured over seven years of $\textit{Fermi}$-Large Area Telescope data. For this we take into account the uncertainty on the shape of the gamma-ray pulse profile. We present an updated radio-timing solution now spanning more than two years and show results from optical modeling of the black-widow-type companion based on new multi-band photometric data taken with HiPERCAM on the Gran Telescopio Canarias on La Palma and ULTRACAM on the New Technology Telescope at ESO La Silla. PSR$\,$J0952$-$0607 is now the fastest-spinning pulsar for which the intrinsic spin-down rate has been reliably constrained ($\dot{P}_\text{int} \lesssim 4.6 \times 10^{-21}\,\text{s}\,\text{s}^{-1}$). The inferred surface magnetic field strength of $B_\text{surf} \lesssim 8.2 \times 10^{7}\,\text{G}$ is among the ten lowest of all known pulsars. This discovery is another example of an extremely fast spinning black-widow pulsar hiding within an unidentified $\textit{Fermi} gamma-ray source. In the future such systems might help to pin down the maximum spin frequency and the minimum surface magnetic field strength of MSPs.

NGC300 ULX1 is an ultra-luminous X-ray pulsar, showing an unprecedented spin evolution, from about 126 s to less than 20 s in only 4 years, consistent with steady mass accretion rate. Following its discovery we have been monitoring the system with Swift/XRT and NICER to further study its properties. We found that even though the observed flux of the system dropped by a factor of $\gtrsim$20, the spin-up rate remained almost constant. A possible explanation is that the decrease in the observed flux is a result of increased absorption of obscuring material due to outflows or a precessing accretion disk.

We have been conducting the GMRT High Resolution Southern Sky (GHRSS) survey for the last four years and have discovered 18 pulsars to date. The GHRSS survey is an off-Galactic-plane survey at 322 MHz in a region of the sky (declination range -40 degrees to -54 degrees) complementary to other ongoing low-frequency surveys. In this paper we report the discovery of three pulsars, PSRs J1239-48, J1516-43 and J1726-52. We also present timing solutions for three pulsars previously discovered with the GHRSS survey: PSR J2144-5237, a millisecond pulsar with a period P=5 ms in a 10 day orbit around a < 0.18 Msun companion; PSR J1516-43, a mildly recycled P=36 ms pulsar in a 228 day orbit with a companion of mass ~0.4 Msun; and the P=320 ms PSR J0514-4408 which we show is a source of pulsed $\gamma$-ray emission. We also report radio polarimetric observations of three of the GHRSS discoveries, PSRs J0418-4154, J0514-4408 and J2144-5237.

Thermally emitting X-ray isolated neutron stars represent excellent targets for testing cooling surface emission and atmosphere models, which are used to infer physical parameters of the neutron star. Among the seven known members of this class, RX J1605.3+3249 is the only one that still lacks confirmation of its spin period. Here we analyze NICER and XMM-Newton observations of RX J1605.3+3249, in order to address its timing and spectral behavior. Contrary to a previous tentative detection, but in agreement with the recent work by Pires et al. (2019), we find no significant pulsation with pulsed fraction higher than 1.3% (3{\sigma}) for periods above 150 ms. We also find a limit of 2.6% for periods above 2 ms, despite searches in different energy bands. The X-ray spectrum can be fit by either a double-blackbody model or by a single-temperature magnetized atmosphere model, both modified by a Gaussian absorption line at ~0.44 keV. The origin of the absorption feature as a proton cyclotron line or as an atomic transition in the neutron star atmosphere is discussed. The predictions of the best-fit X-ray models extended to IR, optical and UV bands are compared with archival data. Our results are interpreted in the framework of a fallback disk scenario.

We report the discovery of soft X-ray pulsations from the nearby millisecond pulsar PSR J1231$-$1411 using NICER. The pulsed emission is characterized by a broad and asymmetric main pulse and a much fainter secondary interpulse, with a total pulsed count rate of 0.055 c s$^{-1}$ in the 0.35-1.5 keV band. We analyzed Fermi LAT data to update the pulse timing model covering 10 years of data and used that model to coherently combine NICER data over a year of observations. Spectral modeling suggests that the flux is dominated by thermal emission from a hot spot (or spots) on the neutron star surface. The phase relationship between the X-ray pulse and the radio and $\gamma$ rays provides insight into the geometry of the system.

We present evidence for three spin-down glitches (or `anti-glitches') in the ultraluminous accreting X-ray pulsar NGC 300 ULX-1 in timing observations made with the Neutron Star Interior Composition Explorer (NICER). Our timing analysis reveals three sudden spin-down events of magnitudes $\Delta\nu = -23, -30,$ and $-43 \,\mu$Hz (fractional amplitudes $\Delta\nu/\nu = -4.4, -5.5,$ and $-7.7 \times 10^{-4}$). We determined fully phase-coherent timing solutions through the first two glitches, giving us high confidence in their detection, while the third candidate glitch is somewhat less secure. These are larger in magnitude (and opposite in sign) than any known radio pulsar glitch. This may be caused by the prolonged rapid spin-up of the pulsar causing a sudden transfer of angular momentum between the superfluid and non-superfluid components of the star. We find no evidence for profile or spectral changes at the epochs of the glitches, supporting the conclusion that these are due to the same process as in normal pulsar glitches, but in reverse.

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.

We report on the Neutron Star Interior Composition Explorer (NICER) monitoring campaign of the 468 Hz accreting millisecond X-ray pulsar IGR J17379-3747. From a detailed spectral and timing analysis of the coherent pulsations we find that they show a strong energy dependence, with soft thermal emission lagging about 640 microseconds behind the hard, Comptonized emission. Additionally, we observe uncommonly large pulse fractions, with measured amplitudes in excess of 20% sinusoidal fractional amplitude across the NICER passband and fluctuations of up to ~70%. Based on a phase-resolved spectral analysis, we suggest that these extreme properties might be explained if the source has an unusually favorable viewing geometry with a large magnetic misalignment angle. Due to these large pulse fractions, we were able to detect pulsations down to quiescent luminosities (~5 x 10^33 erg s^-1). We discuss these low-luminosity pulsations in the context of transitional millisecond pulsars.

The first multimessenger observation of a neutron star merger was independently detected in gamma-rays by Fermi-GBM and INTEGRAL SPI-ACS and gravitational waves by Advanced LIGO and Advanced Virgo. Gravitational waves are emitted from systems with accelerating quadrupole moments, and detectable sources are expected to be compact objects. Nearly all distant astrophysical gamma-ray sources are compact objects. Therefore, serendipitous observations of these two messengers will continue to uncover the sources of gravitational waves and gamma-rays, and enable multimessenger science across the Astro2020 thematic areas. This requires upgrades to the ground-based gravitational wave network and ~keV-MeV gamma-ray coverage for observations of neutron star mergers, and broadband coverage in both gravitational waves and gamma-rays to monitor other expected joint sources.

The unknown state of matter at ultra-high density, large proton/neutron number asymmetry, and low temperature is a major long-standing problem in modern physics. Neutron stars provide the only known setting in the Universe where matter in this regime can stably exist. Valuable information about the interior structure of neutron stars can be extracted via sensitive observations of their exteriors. There are several complementary techniques that require different combinations of high time resolution, superb spectral resolution, and high spatial resolution. In the upcoming decade and beyond, measurements of the masses and radii of an ensemble of neutron stars using these techniques, based on data from multiple proposed next-generation X-ray telescopes, can produce definitive empirical constraints on the allowed dense matter equation of state.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

X-ray photometry and optical spectra are presented covering the periastron passage of the highly-eccentric, ~50 year binary system PSR J2032+4127 in November 2017. This system consists of a 143 ms pulsar in orbit around a massive OB star, MT 91-213. The data show dramatic changes during the encounter as the pulsar wind collided with the stellar wind. The X-ray flux rose on the approach to periastron, then underwent a major dip in the few days around periastron, and then gradually declined over the next few weeks. The optical spectroscopy revealed a steady decline in the H? line strength on the approach to periastron (from an Equivalent Width of -15A to -7A) implying a truncation of the OB star's circumstellar disk by the approaching neutron star. Smooth Particle Hydrodynamic (SPH) modelling is used here to model the system within the context of the observed behaviour and predict the geometrical configuration of the circumstellar disk with respect to the pulsar's orbit.

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.

Access to 50 years of data has led to the discovery of pulsar emission and rotation variability on timescales of months and years. Most of this long-term variability has been seen in long-period pulsars, with relatively little focus on recycled millisecond pulsars. We have analyzed a 38-pulsar sub-set of the 45 millisecond pulsars in the NANOGrav 11-year data set, in order to review their pulse profile stability. The most variability, on any timescale, is seen in PSRs J1713+0747, B1937+21 and J2145-0750. The strongest evidence for long-timescale pulse profile changes is seen in PSRs B1937+21 and J1643-1224. We have focused our analyses on these four pulsars in an attempt to elucidate the causes of their profile variability. Effects of scintillation seem to be responsible for the profile modifications of PSR J2145-0750. We see evidence that imperfect polarization calibration contributes to the profile variability of PSRs J1713+0747 and B1937+21, along with radio frequency interference around 2 GHz, but find that propagation effects also have an influence. The changes seen in PSR J1643-1224 have been reported previously, yet elude explanation beyond their astrophysical nature. Regardless of cause, unmodeled pulse profile changes are detrimental to the accuracy of pulsar timing and must be incorporated into the timing models where possible.

We report on the discovery by the Nuclear Spectroscopic Telescope Array (NuSTAR) and the Neutron Star Interior Composition Explorer (NICER) of the accreting millisecond X-ray pulsar IGR J17591-2342, detecting coherent X-ray pulsations around 527.4 Hz (1.9 ms) with a clear Doppler modulation. This implies an orbital period of ~8.8 hours and a projected semi-major axis of ~1.23 lt-s. From the binary mass function, we estimate a minimum companion mass of 0.42 solar masses, obtained assuming a neutron star mass of 1.4 solar masses and an inclination angle lower than 60 degrees, as suggested by the absence of eclipses or dips in the light-curve of the source. The broad-band energy spectrum is dominated by Comptonisation of soft thermal seed photons with a temperature of ~0.7 keV by electrons heated to 21 keV. We also detect black-body-like thermal direct emission compatible with an emission region of a few kilometers and temperature compatible with the seed source of Comptonisation. A weak Gaussian line centered on the iron K-alpha; complex can be interpreted as a signature of disc reflection. A similar spectrum characterises the NICER spectra, measured during the outburst fading.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has observed dozens of millisecond pulsars for over a decade. We have accrued a large collection of dispersion measure (DM) measurements sensitive to the total electron content between Earth and the pulsars at each observation. All lines of sight cross through the solar wind which produces correlated DM fluctuations in all pulsars. We develop and apply techniques for extracting the imprint of the solar wind from the full collection of DM measurements in the recently released NANOGrav 11-yr data set. We filter out long time scale DM fluctuations attributable to structure in the interstellar medium and carry out a simultaneous analysis of all pulsars in our sample that can differentiate the correlated signature of the wind from signals unique to individual lines of sight. When treating the solar wind as spherically symmetric and constant in time, we find the electron number density at 1 A.U. to be $7.9\pm0.2$ cm$^{-3}$. Our data shows little evidence of long-term variation in the density of the wind. We argue that our techniques paired with a high cadence, low radio frequency observing campaign of near-ecliptic pulsars would be capable of mapping out large-scale latitudinal structure in the wind.

We present results of recent Neutron Star Interior Composition Explorer observations of the accreting millisecond X-ray pulsar IGR J17062-6143 that show that it resides in a circular, ultracompact binary with a 38 minute orbital period. NICER observed the source for approximately 26 ksec over a 5.3 day span in 2017 August, and again for 14 and 11 ksec in 2017 October and November, respectively. A power spectral analysis of the August exposure confirms the previous detection of pulsations at 163.656 Hz in Rossi X-ray Timing Explorer data, and reveals phase modulation due to orbital motion of the neutron star. A coherent search for the orbital solution using the Z^2 method finds a best-fitting circular orbit with a period of 2278.21 s (37.97 min), a projected semi-major axis of 0.00390 lt-sec, and a barycentric pulsar frequency of 163.6561105 Hz. This is currently the shortest known orbital period for an AMXP. The mass function is 9.12 e-8} solar masses, presently the smallest known for a stellar binary. The minimum donor mass ranges from about 0.005 - 0.007 solar masses, for a neutron star mass from 1.2 - 2 solar masses. Assuming mass transfer is driven by gravitational radiation, we find donor mass and binary inclination bounds of 0.0175 - 0.0155 solar masses and 19 deg < i < 27.5 deg, where the lower and upper bounds correspond to 1.4 and 2 solar mass neutron stars, respectively. Folding the data accounting for the orbital modulation reveals a sinusoidal profile with fractional amplitude 2.04 +- 0.11 % (0.3 - 3.2 keV).

We report on the coherent timing analysis of the 182 Hz accreting millisecond X-ray pulsar Swift J1756.9$-$2508 during its 2018 outburst as observed with the Neutron Star Interior Composition Explorer (NICER). Combining our NICER observations with Rossi X-ray Timing Explorer observations of the 2007 and 2009 outbursts, we also studied the long-term spin and orbital evolution of this source. We find that the binary system is well described by a constant orbital period model, with an upper limit on the orbital period derivative of $|\dot P_b| < 7.4\times 10^{-13}$ s s$^{-1}$. Additionally, we improve upon the source coordinates through astrometric analysis of the pulse arrival times, finding R.A. = $17^{\rm h}56^{\rm m}57.18^{\rm s}\pm0.08^{\rm s}$ and Decl. = $-25^{\circ}06'27.8''\pm3.5''$, while simultaneously measuring the long-term spin frequency derivative as $\dot\nu = -7.3\times 10^{-16}$ Hz s$^{-1}$. We briefly discuss the implications of these measurements in the context of the wider population of accreting millisecond pulsars.

Three periastron passages of the PSR B1259$-$63/LS 2883 binary system, consisting of a 48 ms rotation-powered pulsar and a $\sim30$ M$_{\odot}$ Be star, have been observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope, in 2010, 2014, and 2017. During the most-recent periastron passage, sustained low-level gamma-ray emission was observed over a $\sim3$-week long interval immediately after periastron, which was followed by an interval of no emission. Sporadic flares were detected starting 40 days post-periastron and lasted approximately 50 days, during which the emission displayed significant spectral curvature, variability on timescales as short as 1.5 minutes, and peak flux levels well in excess of the pulsar spin-down power. By contrast, during the 2010 and 2014 periastron passages, significant gamma-ray emission was not observed with the LAT until 30 and 32 days post-periastron, respectively. The previous flares did not exhibit spectral curvature, showed no short term variability, and did not exceed the pulsar spin-down power. The high flux and short timescales observed in 2017 suggest significant beaming of the emission is required and constrain the size of the emission region. The flares occur long enough after periastron that the neutron star should already have passed through the extended disk-like outflow, thus constraining options for target material and seed photon sources for inverse Compton models.

We describe the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept that will provide an unprecedented view of the X-ray sky, performing timing and spectroscopy over both a broad energy band (0.2-30 keV) and a wide range of timescales from microseconds to years. STROBE-X comprises two narrow-field instruments and a wide field monitor. The soft or low-energy band (0.2-12 keV) is covered by an array of lightweight optics (3-m focal length) that concentrate incident photons onto small solid-state detectors with CCD-level (85-175 eV) energy resolution, 100 ns time resolution, and low background rates. This technology has been fully developed for NICER and will be scaled up to take advantage of the longer focal length of STROBE-X. The higher-energy band (2-30 keV) is covered by large-area, collimated silicon drift detectors that were developed for the European LOFT mission concept. Each instrument will provide an order of magnitude improvement in effective area over its predecessor (NICER in the soft band and RXTE in the hard band). Finally, STROBE-X offers a sensitive wide-field monitor (WFM), both to act as a trigger for pointed observations of X-ray transients and also to provide high duty-cycle, high time-resolution, and high spectral-resolution monitoring of the variable X-ray sky. The WFM will boast approximately 20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger investigations with a large instantaneous field of view. This mission concept will be presented to the 2020 Decadal Survey for consideration.

We present the polarization pulse profiles for 28 pulsars observed with the Arecibo Observatory by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing project at 2.1 GHz, 1.4 GHz, and 430 MHz. These profiles represent some of the most sensitive polarimetric millisecond pulsar profiles to date, revealing the existence of microcomponents (that is, pulse components with peak intensities much lower than the total pulse peak intensity). Although microcomponents have been detected in some pulsars previously, we present microcomponents for PSRs B1937+21, J1713+0747, and J2234+0944 for the first time. These microcomponents can have an impact on pulsar timing, geometry, and flux density determination. We present rotation measures for all 28 pulsars, determined independently at different observation frequencies and epochs, and find the Galactic magnetic fields derived from these rotation measures to be consistent with current models. These polarization profiles were made using measurement equation template matching, which allows us to generate the polarimetric response of the Arecibo Observatory on an epoch-by-epoch basis. We use this method to describe its time variability, and find that the polarimetric responses of the Arecibo Observatory's 1.4 and 2.1 GHz receivers vary significantly with time.

Swift J0243.6+6124 is a newly discovered Galactic Be/X-ray binary, revealed in late September 2017 in a giant outburst with a peak luminosity of 2E+39 (d/7 kpc)^2 erg/s (0.1-10 keV), with no formerly reported activity. At this luminosity, Swift J0243.6+6124 is the first known galactic ultraluminous X-ray pulsar. We describe Neutron star Interior Composition Explorer (NICER)} and Fermi Gamma-ray Burst Monitor (GBM) timing and spectral analyses for this source. A new orbital ephemeris is obtained for the binary system using spin-frequencies measured with GBM and 15-50 keV fluxes measured with the Neil Gehrels Swift Observatory Burst Alert Telescope to model the system's intrinsic spin-up. Power spectra measured with NICER show considerable evolution with luminosity, including a quasi-periodic oscillation (QPO) near 50 mHz that is omnipresent at low luminosity and has an evolving central frequency. Pulse profiles measured over the combined 0.2-100 keV range show complex evolution that is both luminosity and energy dependent. Near the critical luminosity of L~1E+38 erg/s, the pulse profiles transition from single-peaked to double peaked, the pulsed fraction reaches a minimum in all energy bands, and the hardness ratios in both NICER and GBM show a turn-over to softening as the intensity increases. This behavior repeats as the outburst rises and fades, indicating two distinct accretion regimes. These two regimes are suggestive of the accretion structure on the neutron star surface transitioning from a Coulomb collisional stopping mechanism at lower luminosities to a radiation-dominated stopping mechanism at higher luminosities. This is the highest observed (to date) value of the critical luminosity, suggesting a magnetic field of B ~1E+13 G.

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.

Black widows and redbacks are binary systems consisting of a millisecond pulsar in a close binary with a companion having matter driven off of its surface by the pulsar wind. X-rays due to an intra-binary shock have been observed from many of these systems, as well as orbital variations in the optical emission from the companion due to heating and tidal distortion. We have been systematically studying these systems in radio, optical and X-rays. Here we will present an overview of X-ray and optical studies of these systems, including new XMM-Newton and NuStar data obtained from several of them, along with new optical photometry.

We describe an image-based method that uses two radio criteria, compactness and spectral index, to identify promising pulsar candidates among Fermi Large Area Telescope (LAT) unassociated sources. These criteria are applied to those radio sources from the Giant Metrewave Radio Telescope all-sky survey at 150 MHz (TGSS ADR1) found within the error ellipses of unassociated sources from the 3FGL catalog and a preliminary source list based on 7 years of LAT data. After follow-up interferometric observations to identify extended or variable sources, a list of 16 compact, steep-spectrum candidates is generated. An ongoing search for pulsations in these candidates, in gamma rays and radio, has found six millisecond pulsars and one normal pulsar. A comparison of this method with existing selection criteria based on gamma-ray spectral and variability properties suggests that the pulsar discovery space using Fermi may be larger than previously thought. Radio imaging is a hitherto underutilized source selection method that can be used, as with other multi-wavelength techniques, in the search for Fermi pulsars.

We report on the analysis of 13 gamma-ray pulsars discovered in the Einstein@Home blind search survey using Fermi Large Area Telescope (LAT) Pass 8 data. The 13 new gamma-ray pulsars were discovered by searching 118 unassociated LAT sources from the third LAT source catalog (3FGL), selected using the Gaussian Mixture Model (GMM) machine learning algorithm on the basis of their gamma-ray emission properties being suggestive of pulsar magnetospheric emission. The new gamma-ray pulsars have pulse profiles and spectral properties similar to those of previously-detected young gamma-ray pulsars. Follow-up radio observations have revealed faint radio pulsations from two of the newly-discovered pulsars, and enabled us to derive upper limits on the radio emission from the others, demonstrating that they are likely radio-quiet gamma-ray pulsars. We also present results from modeling the gamma-ray pulse profiles and radio profiles, if available, using different geometric emission models of pulsars. The high discovery rate of this survey, despite the increasing difficulty of blind pulsar searches in gamma rays, suggests that new systematic surveys such as presented in this article should be continued when new LAT source catalogs become available.

This chapter describes the development of X-ray Navigation at the Naval Research Laboratory (NRL) within its astrophysics research programs. The prospects for applications emerged from early discoveries of X-ray source classes and their properties. Starting around 1988 some NRL X-ray astronomy programs included navigation as one of the motivations. The USA experiment (1999) was the first flight payload with an explicit X-ray navigation theme. Subsequently, NRL has continued to work in this area through participation in DARPA and NASA programs. Throughout, the general concept of X-ray navigation (XRNAV) has been broad enough to encompass many different uses of X-ray source observations for attitude determination, position determination, and timekeeping. Pulsar-based X-ray navigation (XNAV) is a special case.

The Neutron Star Interior Composition Explorer (NICER) presents an exciting new capability for exploring the modulation properties of X-ray emitting neutron stars, including large area, low background, extremely precise absolute event time stamps, superb low-energy response and flexible scheduling. The Pulsation Searches and Multiwavelength Coordination working group has designed a 2.5 Ms observing program to search for emission and characterize the modulation properties of about 30 known or suspected neutron star sources across a number of source categories. A key early goal will be to search for pulsations from millisecond pulsars that might exhibit thermal pulsations from the surface suitable for pulse profile modeling to constrain the neutron star equation of state. In addition, we will search for pulsations from transitional millisecond pulsars, isolated neutron stars, low-mass X-ray binaries (LMXBs), accretion-powered millisecond pulsars, central compact objects and other sources. We present our science plan and initial results from the first months of the NICER mission, including the discovery of pulsations from the millisecond pulsar J1231-1411.

The Station Explorer for X-ray Timing and Navigation Technology (SEXTANT) is a technology demonstration enhancement to the Neutron Star Interior Composition Explorer (NICER) mission, which is scheduled to launch in 2017 and will be hosted as an externally attached payload on the International Space Station (ISS). During NICER's 18-month baseline science mission to understand ultra-dense matter through observations of neutron stars in the soft X-ray band, SEXTANT will, for the first-time, demonstrate real-time, on-board X-ray pulsar navigation. Using NICER/SEXTANT as an example, we describe the factors that determine the measurement errors on pulse times of arrival, including source and background count rates, and pulse profile shapes. We then describe properties of the SEXTANT navigation pulsar catalog and prospects for growing it once NICER launches. Finally, we describe the factors affecting the prediction of pulse arrival times in advance, including variable interstellar propagation effect and red timing noise. Together, all of these factors determine how well a particular realization of an X-ray pulsar-based navigation system will perform.

The Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X) probes strong gravity for stellar mass to supermassive black holes and ultradense matter with unprecedented effective area, high time-resolution, and good spectral resolution, while providing a powerful time-domain X-ray observatory.

We report the discovery of PSR J0952$-$0607, a 707-Hz binary millisecond pulsar which is now the fastest-spinning neutron star known in the Galactic field (i.e., outside of a globular cluster). PSR J0952$-$0607 was found using LOFAR at a central observing frequency of 135 MHz, well below the 300 MHz to 3 GHz frequencies typically used in pulsar searches. The discovery is part of an ongoing LOFAR survey targeting unassociated Fermi Large Area Telescope $\gamma$-ray sources. PSR J0952$-$0607 is in a 6.42-hr orbit around a very low-mass companion ($M_\mathrm{c}\gtrsim0.02$ M$_\odot$) and we identify a strongly variable optical source, modulated at the orbital period of the pulsar, as the binary companion. The light curve of the companion varies by 1.6 mag from $r^\prime=22.2$ at maximum to $r^\prime>23.8$, indicating that it is irradiated by the pulsar wind. Swift observations place a 3-$\sigma$ upper limit on the $0.3-10$ keV X-ray luminosity of $L_X < 1.1 \times 10^{31}$ erg s$^{-1}$ (using the 0.97 kpc distance inferred from the dispersion measure). Though no eclipses of the radio pulsar are observed, the properties of the system classify it as a black widow binary. The radio pulsed spectrum of PSR J0952$-$0607, as determined through flux density measurements at 150 and 350 MHz, is extremely steep with $\alpha\sim-3$ (where $S \propto \nu^{\alpha}$). We discuss the growing evidence that the fastest-spinning radio pulsars have exceptionally steep radio spectra, as well as the prospects for finding more sources like PSR J0952$-$0607.

In 2012, five high school students involved in the Pulsar Search Collaboratory discovered the millisecond pulsar PSR J1400$-$1431 and initial timing parameters were published in Rosen et al. (2013) a year later. Since then, we have obtained a phase-connected timing solution spanning five years, resolving a significant position discrepancy and measuring $\dot{P}$, proper motion, parallax, and a monotonic slope in dispersion measure over time. Due to PSR J1400$-$1431's proximity and significant proper motion, we use the Shklovskii effect and other priors to determine a 95% confidence interval for PSR J1400$-$1431's distance, $d=270^{+130}_{-80}$ pc. With an improved timing position, we present the first detection of the pulsar's low-mass white dwarf (WD) companion using the Goodman Spectrograph on the 4.1-m SOAR telescope. Deeper imaging suggests that it is a cool DA-type WD with $T_{\rm eff}=3000\pm100$ K and $R/R_\odot=(2.19\pm0.03)\times10^{-2}\,(d/270\,{\rm pc})$. We show a convincing association between PSR J1400$-$1431 and a $\gamma$-ray point source, 3FGL J1400.5$-$1437, but only weak (3.3-$\sigma$) evidence of pulsations after folding $\gamma$-ray photons using our radio timing model. We detect an X-ray counterpart with XMM-Newton but the measured X-ray luminosity ($1\times10^{29}$ ergs s$^{-1}$) makes PSR J1400$-$1431 the least X-ray luminous rotation-powered millisecond pulsar (MSP) detected to date. Together, our findings present a consistent picture of a nearby ($d\approx230$ pc) MSP in a 9.5-day orbit around a cool, $\sim$0.3 M$_\odot$ WD companion, with orbital inclination, $i\gtrsim60^\circ$.

As of today, more than 2500 pulsars have been found, nearly all in the Milky Way, with the exception of ~28 pulsars in the Small and Large Magellanic Clouds. However, there have been few published attempts to search for pulsars deeper in our Galactic neighborhood. Two of the more promising Local Group galaxies are IC 10 and NGC 6822 (also known as Barnard's Galaxy) due to their relatively high star formation rate and their proximity to our galaxy. IC 10 in particular, holds promise as it is the closest starburst galaxy to us and harbors an unusually high number of Wolf-Rayet stars, implying the presence of many neutron stars. We observed IC 10 and NGC 6822 at 820 MHz with the Green Bank Telescope for ~15 and 5 hours respectively, and put a strong upper limit of 0.1 mJy on pulsars in either of the two galaxies. We also performed single pulse searches of both galaxies with no firm detections.

We present an analysis of the gamma-ray measurements by the Large Area Telescope onboard the \textit{Fermi Gamma-ray Space Telescope} in the region of the supernova remnant~(SNR) Monoceros Loop~(G205.5$+$0.5). The brightest gamma-ray peak is spatially correlated with the Rosette Nebula, which is a molecular cloud complex adjacent to the southeast edge of the SNR. After subtraction of this emission by spatial modeling, the gamma-ray emission from the SNR emerges, which is extended and fit by a Gaussian spatial template. The gamma-ray spectra are significantly better reproduced by a curved shape than a simple power law. The luminosities between 0.2--300~GeV are $\sim$~$4 \times 10^{34}$~erg~s$^{-1}$ for the SNR and $\sim$~$3 \times 10^{34}$~erg~s$^{-1}$ for the Rosette Nebula, respectively. We argue that the gamma rays likely originate from the interactions of particles accelerated in the SNR. The decay of neutral pions produced in nucleon-nucleon interactions of accelerated hadrons with interstellar gas provides a reasonable explanation for the gamma-ray emission of both the Rosette Nebula and the Monoceros SNR.

LOFT-P is a concept for a NASA Astrophysics Probe-Class (<$1B) X-ray timing mission, based on the LOFT concept originally proposed to ESAs M3 and M4 calls. LOFT-P requires very large collecting area (>6 m^2, >10x RXTE), high time resolution, good spectral resolution, broad-band spectral coverage (2-30 keV), highly flexible scheduling, and an ability to detect and respond promptly to time-critical targets of opportunity. It addresses science questions such as: What is the equation of state of ultra dense matter? What are the effects of strong gravity on matter spiraling into black holes? It would be optimized for sub-millisecond timing to study phenomena at the natural timescales of neutron star surfaces and black hole event horizons and to measure mass and spin of black holes. These measurements are synergistic to imaging and high-resolution spectroscopy instruments, addressing much smaller distance scales than are possible without very long baseline X-ray interferometry, and using complementary techniques to address the geometry and dynamics of emission regions. A sky monitor (2-50 keV) acts as a trigger for pointed observations, providing high duty cycle, high time resolution monitoring of the X-ray sky with ~20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger studies. A probe-class mission concept would employ lightweight collimator technology and large-area solid-state detectors, technologies which have been recently greatly advanced during the ESA M3 study. Given the large community interested in LOFT (>800 supporters, the scientific productivity of this mission is expected to be very high, similar to or greater than RXTE (~2000 refereed publications). We describe the results of a study, recently completed by the MSFC Advanced Concepts Office, that demonstrates that LOFT-P is feasible within a NASA probe-class mission budget.

PSR J1024$-$0719 is a millisecond pulsar that was long thought to be isolated. However, puzzling results concerning its velocity, distance, and low rotational period derivative have led to reexamination of its properties. We present updated radio timing observations along with new and archival optical data that show PSR J1024$-$0719 is most likely in a long period (2$-$20 kyr) binary system with a low-mass ($\approx 0.4\,M_\odot$) low-metallicity ($Z \approx -0.9\,$ dex) main sequence star. Such a system can explain most of the anomalous properties of this pulsar. We suggest that this system formed through a dynamical exchange in a globular cluster that ejected it into a halo orbit, consistent with the low observed metallicity for the stellar companion. Further astrometric and radio timing observations such as measurement of the third period derivative could strongly constrain the range of orbital parameters.

In a search of the unidentified Fermi gamma-ray source 3FGL J1417.5-4402 with the Parkes radio telescope, we discovered PSR J1417-4402, a 2.66 ms pulsar having the same 5.4 day orbital period as the optical and X-ray binary identified by Strader et al. The existence of radio pulsations implies that the neutron star is currently not accreting. Substantial outflows from the companion render the radio pulsar undetectable for more than half of the orbit, and may contribute to the observed Halpha emission. Our initial pulsar observations, together with the optically inferred orbit and inclination, imply a mass ratio of 0.171+/-0.002, a companion mass of M_2=0.33+/-0.03 Msun, and a neutron star mass in the range 1.77<M_1<2.13 Msun. However, there remains a discrepancy between the distance of 4.4 kpc inferred from the optical properties of the companion and the smaller radio dispersion measure distance of 1.6 kpc. The smaller distance would reduce the inferred Roche-lobe filling factor, increase the inferred inclination angle, and decrease the masses. As a wide binary, PSR J1417-4402 differs from the radio-eclipsing black widow and redback pulsars being discovered in large numbers by Fermi. It is probably a system that began mass transfer onto the neutron star after the companion star left the main sequence. The companion should end its evolution as a He white dwarf in a 6-20 day orbit, i.e., as a typical binary millisecond pulsar companion.

Millisecond pulsars (MSPs) represent nearly half of the more than 160 currently known $\gamma$-ray pulsars detected by the Large Area Telescope on the \textit{Fermi} satellite, and a third of all known MSPs are seen in $\gamma$ rays. The least energetic $\gamma$-ray MSPs enable us to probe the so-called deathline for high-energy emission, i.e., the spin-down luminosity limit under which pulsars (PSRs) cease to produce detectable high-energy radiation. Characterizing the MSP luminosity distribution helps to determine their contribution to the Galactic diffuse $\gamma$-ray emission. We made use of the high-quality pulsar timing data recorded at the Nan\c{c}ay Radio Telescope over several years to characterize the properties of a selection of MSPs. For one of the pulsars, the dataset was complemented with Westerbork Synthesis Radio Telescope observations. The rotation ephemerides derived from this analysis were also used to search the LAT data for new $\gamma$-ray MSPs. For the MSPs considered in this study, we obtained new transverse proper motion measurements or updated the existing ones, and placed new distance constraints for some of them, with four new timing parallax measurements. We discovered significant GeV $\gamma$-ray signals from four MSPs, i.e., PSRs J0740+6620, J0931$-$1902, J1455$-$3330, and J1730$-$2304. The latter is now the least energetic $\gamma$-ray pulsar found to date. Despite the improved $\dot E$ and $L_\gamma$ estimates, the relationship between these two quantities remains unclear, especially at low $\dot E$ values.

We have discovered six radio millisecond pulsars (MSPs) in a search with the Arecibo telescope of 34 unidentified gamma-ray sources from the Fermi Large Area Telescope (LAT) 4-year point source catalog. Among the 34 sources, we also detected two MSPs previously discovered elsewhere. Each source was observed at a center frequency of 327 MHz, typically at three epochs with individual integration times of 15 minutes. The new MSP spin periods range from 1.99 to 4.66 ms. Five of the six pulsars are in interacting compact binaries (period < 8.1 hr), while the sixth is a more typical neutron star-white dwarf binary with an 83-day orbital period. This is a higher proportion of interacting binaries than for equivalent Fermi-LAT searches elsewhere. The reason is that Arecibo's large gain afforded us the opportunity to limit integration times to 15 minutes, which significantly increased our sensitivity to these highly accelerated systems. Seventeen of the remaining 26 gamma-ray sources are still categorized as strong MSP candidates, and will be re-searched.

We report on radio timing and multiwavelength observations of the 4.66 ms redback pulsar J1048+2339, which was discovered in an Arecibo search targeting the Fermi-LAT source 3FGLJ1048.6+2338. Two years of timing allowed us to derive precise astrometric and orbital parameters for the pulsar. PSR J1048+2339 is in a 6-hour binary, and exhibits radio eclipses over half the orbital period and rapid orbital period variations. The companion has a minimum mass of 0.3 solar masses, and we have identified a $V \sim 20$ variable optical counterpart in data from several surveys. The phasing of its $\sim 1$~mag modulation at the orbital period suggests highly efficient and asymmetric heating by the pulsar wind, which may be due to an intrabinary shock that is distorted near the companion, or to the companion's magnetic field channeling the pulsar wind to specific locations on its surface. We also present gamma-ray spectral analysis of the source and preliminary results from searches for gamma-ray pulsations using the radio ephemeris.

We are conducting a survey for pulsars and transients using the Giant Metrewave Radio Telescope (GMRT). The GMRT High Resolution Southern Sky (GHRSS) survey is an off-Galactic-plane (|b|>5) survey in the declination range -40 deg to -54 deg at 322 MHz. With the high time (up to 30.72 micro-sec) and frequency (up to 0.016275 MHz) resolution observing modes, the 5-sigma detection limit is 0.5 mJy for a 2 ms pulsar with 10% duty cycle at 322 MHz. Total GHRSS sky coverage of 2866 square-deg, will result from 1953 pointing, each covering 1.8 square-deg. The 10-sigma detection limit for a 5 milli-sec transient burst is 1.6 Jy for the GHRSS survey. In addition, the GHRSS survey can reveal transient events like rotating radio transients or fast radio bursts. With 35% of the survey completed (i.e. 1000 square-deg), we report the discovery of 10 pulsars, one of which is a millisecond pulsar (MSP), one of the highest pulsar per square degree discovery rates for any off-Galactic plane survey. We re-detected 23 known in-beam pulsars. Utilising the imaging capability of the GMRT we also localised 4 of the GHRSS pulsars (including the MSP) in the gated image plane within +/- 10 arc-second. We demonstrated rapid convergence in pulsar timing with a more precise position than is possible with single dish discoveries. We also exhibited that we can localise the brightest transient sources with simultaneously obtained lower time resolution imaging data, demonstrating a technique that may have application in the SKA.

We have constructed timing solutions for 81 gamma-ray pulsars covering more than five years of Fermi data. The sample includes 37 radio-quiet or radio-faint pulsars which cannot be timed with other telescopes. These timing solutions and the corresponding pulse times of arrival are prerequisites for further study, e.g. phase-resolved spectroscopy or searches for mode switches. Many gamma-ray pulsars are strongly affected by timing noise, and we present a new method for characterizing the noise process and mitigating its effects on other facets of the timing model. We present an analysis of timing noise over the population using a new metric for characterizing its strength and spectral shape, namely its time-domain correlation. The dependence of the strength on spin frequency and spin-down rate is in good agreement with previous studies. We find that noise process power spectra $S(f)$ for unrecycled pulsars are steep, with strong correlations over our entire data set and spectral indices $S(f)\propto f^{-\alpha}$ of 5 to 9. One possible explanation for these results is the occurrence of unmodelled, episodic 'microglitches'. Finally, we show that our treatment of timing noise results in robust parameter estimation, and in particular we measure a precise timing position for each pulsar. We extensively validate our results with multi-wavelength astrometry, and using our updated position, we firmly identify the X-ray counterpart of PSR J1418-6058.

In a search with the Parkes radio telescope of 56 unidentified Fermi-LAT gamma-ray sources, we have detected 11 millisecond pulsars (MSPs), 10 of them discoveries, of which five were reported in Kerr et al. (2012). We did not detect radio pulsations from another six pulsars now known in these sources. We describe the completed survey, which included multiple observations of many targets done to minimize the impact of interstellar scintillation, acceleration effects in binary systems, and eclipses. We consider that 23 of the 39 remaining sources may still be viable pulsar candidates. We present timing solutions and polarimetry for five of the MSPs, and gamma-ray pulsations for PSR J1903-7051 (pulsations for five others were reported in the second Fermi-LAT catalog of gamma-ray pulsars). Two of the new MSPs are isolated and five are in >1 d circular orbits with 0.2-0.3 Msun presumed white dwarf companions. PSR J0955-6150, in a 24 d orbit with a ~0.25 Msun companion but eccentricity of 0.11, belongs to a recently identified class of eccentric MSPs. PSR J1036-8317 is in an 8 hr binary with a >0.14 Msun companion that is probably a white dwarf. PSR J1946-5403 is in a 3 hr orbit with a >0.02 Msun companion with no evidence of radio eclipses.

We present the Robert C. Byrd Green Bank Telescope discovery of the highly eccentric binary millisecond pulsar PSR J1835$-$3259A in the Fermi Large Area Telescope-detected globular cluster NGC 6652. Timing over one orbit yields the pulse period 3.89 ms, orbital period 9.25 d, eccentricity $\sim 0.95$, and an unusually high companion mass of $0.74\,M_{\odot}$ assuming a $1.4\,M_{\odot}$ pulsar. We caution that the lack of data near periastron prevents a precise measurement of the eccentricity, and that further timing is necessary to constrain this and the other orbital parameters. From tidal considerations, we find that the companion must be a compact object. This system likely formed through an exchange encounter in the dense cluster environment. Our initial timing results predict the measurements of at least two post-Keplerian parameters with long-term phase-connected timing: the rate of periastron advance $\dot{\omega} \sim 0.1^{\circ}\,$yr$^{-1}$, requiring 1 yr of phase connection; and the Einstein delay $\gamma_{\mathrm{GR}} \sim 10\,$ms, requiring 2-3 yr of timing. For an orbital inclination $i > 50^{\circ}$, a measurement of $\sin{i}$ is also likely. PSR J1835$-$3259A thus provides an opportunity to measure the neutron star mass with high precision; to probe the cluster environment; and, depending on the nature of the companion, to investigate the limits of general relativity.

The 1.69 ms spin period of PSR J1227-4853 was recently discovered in radio observations of the low-mass X-ray binary XSS J12270-4859 following the announcement of a possible transition to a rotation-powered millisecond pulsar state, inferred from decreases in optical, X-ray, and gamma-ray flux from the source. We report the detection of significant (5$\sigma$) gamma-ray pulsations after the transition, at the known spin period, using ~1 year of data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. The gamma-ray light curve of PSR J1227-4853 can be fit by one broad peak, which occurs at nearly the same phase as the main peak in the 1.4 GHz radio profile. The partial alignment of light-curve peaks in different wavebands suggests that at least some of the radio emission may originate at high altitude in the pulsar magnetosphere, in extended regions co-located with the gamma-ray emission site. We folded the LAT data at the orbital period, both pre- and post-transition, but find no evidence for significant modulation of the gamma-ray flux. Analysis of the gamma-ray flux over the mission suggests an approximate transition time of 2012 November 30. Continued study of the pulsed emission and monitoring of PSR J1227-4853, and other known redback systems, for subsequent flux changes will increase our knowledge of the pulsar emission mechanism and transitioning systems.

We present initial pulsar results from the first station of the Long Wavelength Array (LWA1) obtained during the commissioning period of LWA1 and early science results. We present detections of periodic emission from 44 previously known pulsars, including 3 millisecond pulsars (MSPs). The effects of the interstellar medium on pulsar emission are significantly enhanced at the low frequencies of the LWA1 band (10--88 MHz), making LWA1 a very sensitive instrument for characterizing changes in dispersion measures (DM) and other effects from the interstellar medium. Pulsars also often have significant evolution in their pulse profile at low frequency and a break in their spectral index. We report DM measurements for 44 pulsars, mean flux density measurements for 36 pulsars, and multi-frequency component spacing and widths for 15 pulsars with more than one profile component. For 27 pulsars, we report spectral index measurements within our frequency range. We also introduce the LWA1 Pulsar Data Archive, which stores reduced data products from LWA1 pulsar observations. Reduced data products for the observations presented here can be found on the archive. Reduced data products from future LWA1 pulsar observations will also be made available through the archive.

We consider the X-ray properties of the redback class of eclipsing millisecond pulsars. These are transitional systems between accreting low-mass X-ray binaries and binary millisecond pulsars orbiting white dwarfs, and hence their companions are non-degenerate and nearly Roche-lobe filling. The X-ray luminosity seems to scale with the fraction of the pulsar sky subtended by the companion, suggesting the shock region is not much larger than the companion, which is supported by modeling of the orbital light curves. The typical X-ray photon spectral index is $\sim 1$ and the typical 0.3-8 keV X-ray efficiency, assuming a shock size on the order of the companion's Roche lobe cross-section, is on the order of 10%. We present an overview of previous investigations, and present new observations of two redbacks, a Chandra observation of PSR J1628$-$3205 and a XMM-Newton observation of PSR J2129$-$0429. The latter shows a clearly double peaked orbital light curve with variation of the non-thermal flux by a factor of $\sim 11$, with peaks around orbital phases 0.6 and 0.9. We suggest the magnetic field of the companion plays a significant role in the X-ray emission from intrabinary shocks in redbacks.

PSR J2032+4127 is a gamma-ray and radio-emitting pulsar which has been regarded as a young luminous isolated neutron star. However, its recent spin-down rate has extraordinarily increased by a factor of two. We present evidence that this is due to its motion as a member of a highly-eccentric binary system with a 15-solar-mass Be star, MT91~213. Timing observations show that, not only are the positions of the two stars coincident within 0.4 arcsec, but timing models of binary motion of the pulsar fit the data much better than a model of a young isolated pulsar. MT91~213, and hence the pulsar, lie in the Cyg~OB2 stellar association, which is at a distance of only 1.4-1.7 kpc. The pulsar is currently on the near side of, and accelerating towards, the Be star, with an orbital period of 20-30 years. The next periastron is well-constrained to occur in early 2018, providing an opportunity to observe enhanced high-energy emission as seen in other Be-star binary systems.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of magnetospheres of isolated neutron stars. For a summary, we refer to the paper.

Since their discovery in the late 1960's the population of known neutron stars (NSs) has grown to ~2500. The last five decades of observations have yielded many surprises and demonstrated that the observational properties of NSs are remarkably diverse. The surveys that will be performed with SKA (the Square Kilometre Array) will produce a further tenfold increase in the number of Galactic NSs known. Moreover, the SKA's broad spectral coverage, sub-arraying and multi-beaming capabilities will allow us to characterise these sources with unprecedented efficiency, in turn enabling a giant leap in the understanding of their properties. Here we review the NS population and outline our strategies for studying each of the growing number of diverse classes that are populating the "NS zoo". Some of the main scientific questions that will be addressed by the much larger statistical samples and vastly improved timing efficiency provided by SKA include: (i) the spin period and spin-down rate distributions (and thus magnetic fields) at birth, and the associated information about the SNe wherein they are formed; (ii) the radio pulsar-magnetar connection; (iii) the link between normal radio pulsars, intermittent pulsars and rotating radio transients; (iv) the slowest possible spin period for a radio pulsar (revealing the conditions at the pulsar death-line); (v) proper motions of pulsars (revealing SN kick physics); (vi) the mass distribution of NSs (vii) the fastest possible spin period for a recycled pulsar (constraining magnetosphere-accretion disc interactions, gravitational wave radiation and the equation-of-state); (viii) the origin of high eccentricity millisecond pulsars (MSPs); (ix) the formation channels for recently identified triple systems; and finally (x) how isolated MSPs are formed. We expect that the SKA will break new ground unveiling exotic systems that will challenge... [abridged]

XSS J12270-4859 is an X-ray binary associated with the Fermi LAT gamma-ray source 1FGL J1227.9-4852. In 2012 December, this source underwent a transition where the X-ray and optical luminosity dropped and the spectral signatures of an accretion disc disappeared. We report the discovery of a 1.69 millisecond pulsar (MSP), PSR J1227-4853, at a dispersion measure of 43.4 pc cm$^{-3}$ associated with this source, using the GMRT at 607 MHz. This demonstrates that, post-transition, the system hosts an active radio MSP. This is the third system after PSR J1023+0038 and PSR J1824-2452I showing evidence of state switching between radio MSP and low-mass X-ray binary (LMXB) states. We report timing observations of PSR J1227-4853 with the GMRT and Parkes, which give a precise determination of the rotational and orbital parameters of the system. The companion mass measurement of 0.17 to 0.46 M$_{sun}$ suggests that this is a redback system. PSR J1227-4853 is eclipsed for about 40% of its orbit at 607 MHz; with additional short-duration eclipses at all orbital phases. We also find that the pulsar is very energetic, with a spin-down luminosity of ~ 10$^{35}$ erg s$^{-1}$. We report simultaneous imaging and timing observations with the GMRT, which suggests that eclipses are caused by absorption, rather than dispersion smearing or scattering.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produce upper limits on the emission amplitude. We highlight the results from seven young pulsars with large spin-down luminosities. We reach within a factor of five of the canonical spin-down limit for all seven of these, whilst for the Crab and Vela pulsars we further surpass their spin-down limits. We present new or updated limits for 172 other pulsars (including both young and millisecond pulsars). Now that the detectors are undergoing major upgrades, and, for completeness, we bring together all of the most up-to-date results from all pulsars searched for during the operations of the first-generation LIGO, Virgo and GEO600 detectors. This gives a total of 195 pulsars including the most recent results described in this paper.

Millisecond pulsars (MSPs) are a growing class of gamma-ray emitters. Pulsed gamma-ray signals have been detected from more than 40 MSPs with the Fermi Large Area Telescope (LAT). The wider radio beams and more compact magnetospheres of MSPs enable studies of emission geometries over a broader range of phase space than non-recycled radio-loud gamma-ray pulsars. We have modeled the gamma-ray light curves of 40 LAT-detected MSPs using geometric emission models assuming a vacuum retarded-dipole magnetic field. We modeled the radio profiles using a single-altitude hollow-cone beam, with a core component when indicated by polarimetry; however, for MSPs with gamma-ray and radio light curve peaks occurring at nearly the same rotational phase we assume that the radio emission is co-located with the gamma rays and caustic in nature. The best-fit parameters and confidence intervals are determined using a maximum likelihood technique. We divide the light curves into three model classes, with gamma-ray peaks trailing (Class I), aligned (Class II) or leading (Class III) the radio peaks. Outer gap and slot gap (two-pole caustic) models best fit roughly equal numbers of Class I and II, while Class III are exclusively fit with pair-starved polar cap models. Distinguishing between the model classes based on typical derived parameters is difficult. We explore the evolution of magnetic inclination angle with period and spin-down power, finding possible correlations. While the presence of significant off-peak emission can often be used as a discriminator between outer gap and slot gap models, a hybrid model may be needed.

Eclipsing millisecond pulsars in close ($P_b < 1$~day) binary systems provide a different view of pulsar winds and shocks than do isolated pulsars. Since 2009, the numbers of these systems known in the Galactic field has increased enormously. We have been systematically studying many of these newly discovered systems at multiple wavelengths. Typically, the companion is nearly Roche-lobe filling and heated by the pulsar which drives mass loss from the companion. The pulsar wind shocks with this material just above the surface of the companion. We discuss various observational properties of this shock, including radio eclipses, orbitally modulated X-ray emission, and the potential for $\gamma$-ray emission. Redbacks, whose companions are likely non-degenerate and significantly more massive, generally have more luminous shocks than black widows which have very low mass companions. This is expected since the more massive redback companions intercept a greater fraction of the pulsar wind. We also compare these systems to accreting millisecond pulsars, which may be progenitors of black widows and in some cases can pass back and forth between redback and accretion phases.

We report the discovery of four gamma-ray pulsars, detected in computing-intensive blind searches of data from the Fermi Large Area Telescope (LAT). The pulsars were found using a novel search approach, combining volunteer distributed computing via Einstein@Home and methods originally developed in gravitational-wave astronomy. The pulsars PSRs J0554+3107, J1422-6138, J1522-5735, and J1932+1916 are young and energetic, with characteristic ages between 35 and 56 kyr and spin-down powers in the range $6\times10^{34}$ - $10^{36}$ erg s$^{-1}$. They are located in the Galactic plane and have rotation rates of less than 10 Hz, among which the 2.1 Hz spin frequency of PSR J0554+3107 is the slowest of any known gamma-ray pulsar. For two of the new pulsars, we find supernova remnants coincident on the sky and discuss the plausibility of such associations. Deep radio follow-up observations found no pulsations, suggesting that all four pulsars are radio-quiet as viewed from Earth. These discoveries, the first gamma-ray pulsars found by volunteer computing, motivate continued blind pulsar searches of the many other unidentified LAT gamma-ray sources.

We describe the first X-ray observations of five short orbital period ($P_B < 1$ day), $\gamma$-ray emitting, binary millisecond pulsars. Four of these, PSRs J0023+0923, J1124$-$3653, J1810+1744, and J2256$-$1024 are `black-widow' pulsars, with degenerate companions of mass $\ll0.1 M_{\odot}$, three of which exhibit radio eclipses. The fifth source, PSR J2215+5135, is an eclipsing `redback' with a near Roche-lobe filling $\sim$0.2 solar mass non-degenerate companion. Data were taken using the \textit{Chandra X-Ray Observatory} and covered a full binary orbit for each pulsar. Two pulsars, PSRs J2215+5135 and J2256$-$1024, show significant orbital variability while PSR J1124$-$3653 shows marginal orbital variability. The lightcurves for these three pulsars have X-ray flux minima coinciding with the phases of the radio eclipses. This phenomenon is consistent with an intrabinary shock emission interpretation for the X-rays. The other two pulsars, PSRs J0023+0923 and J1810+1744, are fainter and do not demonstrate variability at a level we can detect in these data. All five spectra are fit with three separate models: a power-law model, a blackbody model, and a combined model with both power-law and blackbody components. The preferred spectral fits yield power-law indices that range from 1.3 to 3.2 and blackbody temperatures in the hundreds of eV. The spectrum for PSR J2215+5135 shows a significant hard X-ray component, with a large number of counts above 2 keV, which is additional evidence for the presence of intrabinary shock emission and is similar to what has been detected in the low-mass X-ray binary to millisecond pulsar transition object PSR J1023+0038.

We report a 5.4\sigma\ detection of pulsed gamma rays from PSR B1821-24 in the globular cluster M28 using ~44 months of Fermi Large Area Telescope (LAT) data that have been reprocessed with improved instrument calibration constants. We constructed a phase-coherent ephemeris, with post-fit residual RMS of 3 \mu s, using radio data spanning ~23.2 years, enabling measurements of the multi-wavelength light curve properties of PSR B1821-24 at the milliperiod level. We fold RXTE observations of PSR B1821-24 from 1996 to 2007 and discuss implications on the emission zones. The gamma-ray light curve consists of two peaks, separated by 0.41$\pm$0.02 in phase, with the first gamma-ray peak lagging the first radio peak by 0.05$\pm$0.02 in phase, consistent with the phase of giant radio pulses. We observe significant emission in the off-peak interval of PSR B1821-24 with a best-fit LAT position inconsistent with the core of M28. We do not detect significant gamma-ray pulsations at the spin or orbital periods from any other known pulsar in M28, and we place limits on the number of energetic pulsars in the cluster. The derived gamma-ray efficiency, ~2%, is typical of other gamma-ray pulsars with comparable spin-down power, suggesting that the measured spin-down rate ($2.2\times10^{36}$ erg s$^{-1}$) is not appreciably distorted by acceleration in the cluster potential. This confirms PSR B1821-24 as the second very energetic millisecond pulsar in a globular cluster and raises the question of whether these represent a separate class of objects that only form in regions of very high stellar density

We present flux density measurements and pulse profiles for the millisecond pulsar PSR J2145-0750 spanning 37 to 81 MHz using data obtained from the first station of the Long Wavelength Array. These measurements represent the lowest frequency detection of pulsed emission from a millisecond pulsar to date. We find that the pulse profile is similar to that observed at 102 MHz. We also find that the flux density spectrum between ~40 MHz to 5 GHz is suggestive of a break and may be better fit by a model that includes spectral curvature with a rollover around 730 MHz rather than a single power law.

The very high energy (VHE; >100 GeV) source HESS J0632+057 has been recently confirmed as a \gamma-ray binary, a subclass of the high mass X-ray binary (HMXB) population, through the detection of an orbital period of 321 days. We performed a deep search for the emission of HESS J0632+057 in the GeV energy range using data from the Fermi Large Area Telescope (LAT). The analysis was challenging due to the source being located in close proximity to the bright \gamma-ray pulsar PSR J0633+0632 and lying in a crowded region of the Galactic plane where there is prominent diffuse emission. We formulated a Bayesian block algorithm adapted to work with weighted photon counts, in order to define the off-pulse phases of PSR J0633+0632. A detailed spectral-spatial model of a 5 deg circular region centred on the known location of HESS J0632+057 was generated to accurately model the LAT data. No significant emission from the location of HESS J0632+057 was detected in the 0.1-100 GeV energy range integrating over ~3.5 years of data; with a 95% flux upper limit of F_{0.1-100 GeV} < 3 x 10-8 ph cm-2 s-1. A search for emission over different phases of the orbit also yielded no significant detection. A search for source emission on shorter timescales (days--months) did not yield any significant detections. We also report the results of a search for radio pulsations using the 100-m Green Bank Telescope (GBT). No periodic signals or individual dispersed bursts of a likely astronomical origin were detected. We estimated the flux density limit of < 90/40 \mu Jy at 2/9 GHz. The LAT flux upper limits combined with the detection of HESS J0632+057 in the 136-400 TeV energy band by the MAGIC collaboration imply that the VHE spectrum must turn over at energies <136 GeV placing constraints on any theoretical models invoked to explain the \gamma-ray emission.

Using the Giant Metrewave Radio Telescope (GMRT) we performed deep observations to search for radio pulsations in the directions of unidentified Fermi Large Area Telescope (LAT) gamma-ray sources. We report the discovery of an eclipsing black-widow millisecond pulsar, PSR J1544+4937, identified with the un-cataloged gamma-ray source Fermi J1544.2+4941. This 2.16 ms pulsar is in a 2.9 hours compact circular orbit with a very low-mass companion (Mc > 0.017 Msun). At 322 MHz this pulsar is found to be eclipsing for 13% of its orbit, whereas at 607 MHz the pulsar is detected throughout the low-frequency eclipse phase. Variations in the eclipse ingress phase are observed, indicating a clumpy and variable eclipsing medium. Moreover, additional short-duration absorption events are observed around the eclipse boundaries. Using the radio timing ephemeris we were able to detect gamma-ray pulsations from this pulsar, confirming it as the source powering the gamma-ray emission.

We are conducting a search for new gamma-ray binaries by making high signal-to-noise light curves of all cataloged Fermi LAT sources and searching for periodic variability using appropriately weighted power spectra. The light curves are created using a variant of aperture photometry where photons are weighted by the probability that they came from the source of interest. From this analysis we find that the light curves of a number of sources near the ecliptic plane are contaminated by gamma-ray emission from the Moon. This shows itself as modulation on the Moon's sidereal period in the power spectra. We demonstrate that this contamination can be removed by excluding times when the Moon was too close to a source. We advocate that this data screening should generally be used when analyzing LAT data from a source located close to the path of the Moon.

In the last few years, over 43 millisecond radio pulsars have been discovered by targeted searches of unidentified gamma-ray sources found by the Fermi Gamma-Ray Space Telescope. A large fraction of these millisecond pulsars are in compact binaries with low-mass companions. These systems often show eclipses of the pulsar signal and are commonly known as black widows and redbacks because the pulsar is gradually destroying its companion. In this paper, we report on the optical discovery of four strongly irradiated millisecond pulsar companions. All four sources show modulations of their color and luminosity at the known orbital periods from radio timing. Light curve modelling of our exploratory data shows that the equilibrium temperature reached on the companion's dayside with respect to their nightside is consistent with about 10-30% of the available spin-down energy from the pulsar being reprocessed to increase the companion's dayside temperature. This value compares well with the range observed in other irradiated pulsar binaries and offers insights about the energetics of the pulsar wind and the production of gamma-ray emission. In addition, this provides a simple way of estimating the brightness of irradiated pulsar companions given the pulsar spin-down luminosity. Our analysis also suggests that two of the four new irradiated pulsar companions are only partially filling their Roche lobe. Some of these sources are relatively bright and represent good targets for spectroscopic follow-up. These measurements could enable, among other things, mass determination of the neutron stars in these systems.

We report on the discovery of gamma-ray pulsations from five millisecond pulsars (MSPs) using the Fermi Large Area Telescope (LAT) and timing ephemerides provided by various radio observatories. We also present confirmation of the gamma-ray pulsations from a sixth source, PSR J2051-0827. Five of these six MSPs are in binary systems: PSRs J1713+0747, J1741+1351, J1600-3053 and the two black widow binary pulsars PSRs J0610-2100 and 2051-0827. The only isolated MSP is the nearby PSR J1024-0719, which is also known to emit X-rays. We present X-ray observations in the direction of PSRs J1600-3053 and J2051-0827. While the latter is firmly detected, we an only give upper limits for the X-ray flux of the former. There are no dedicated X-ray observations available for the other 3 objects. The MSPs mentioned above, together with most of the MSPs detected by Fermi, are used to put together a sample of 30 gamma-ray MSPs which is used to study the morphology and phase connection of radio and gamma-ray pulse profiles. We show that MSPs with pulsed gamma-ray emission which is phase aligned with the radio emission present the steepest radio spectra and the largest magnetic fields at the light cylinder among all MSPs. As well, we also observe a trend towards very low, or undetectable, radio linear polarisation levels. These properties could be attributed to caustic radio emission produced at a range of different altitudes in the magnetosphere. We note that most of these characteristics are also observed in the Crab pulsar, the only other radio pulsar known to exhibit phase-aligned radio and gamma-ray emission.

LWA1 is a new radio telescope operating in the frequency range 10-88 MHz, located in central New Mexico. The telescope consists of 258 pairs of dipole-type antennas whose outputs are individually digitized and formed into beams. Simultaneously, signals from all dipoles can be recorded using one of the instrument's "all dipoles" modes, facilitating all-sky imaging. Notable features of the instrument include high intrinsic sensitivity (about 6 kJy zenith system equivalent flux density), large instantaneous bandwidth (up to 78 MHz), and 4 independently-steerable beams utilizing digital "true time delay" beamforming. This paper summarizes the design of LWA1 and its performance as determined in commissioning experiments. We describe the method currently in use for array calibration, and report on measurements of sensitivity and beamwidth.

Using the 100-m Effelsberg radio telescope operating at 1.36 GHz, we have performed a targeted radio pulsar survey of 289 unassociated gamma-ray sources discovered by the Large Area Telescope (LAT) aboard the Fermi satellite and published in the 1FGL catalogue (Abdo et al., 2010). This survey resulted in the discovery of millisecond pulsar J1745+1017, which resides in a short-period binary system with a low-mass companion, Mmin ~ 0.0137 Msun, indicative of `Black Widow' type systems. A two-year timing campaign has produced a refined radio ephemeris, accurate enough to allow for phase-folding of the LAT photons, resulting in the detection of a dual-peaked gamma-ray light-curve, proving that PSR J1745+1017 is the source responsible for the gamma-ray emission seen in 1FGL J1745.5 + 1018 (2FGL J1745.6+1015; Nolan et al., 2012). We find the gamma-ray spectrum of PSR J1745+1017 to be well modelled by an exponentially-cut-off power law with cut-off energy 3.2 GeV and photon index 1.6. The observed sources are known to contain a further 10 newly discovered pulsars which were undetected in this survey. Our radio observations of these sources are discussed and in all cases limiting flux densities are calculated. The reasons behind the seemingly low yield of discoveries are also discussed.

We report the detection of radio emission from PSR J1311-3430, the first millisecond pulsar discovered in a blind search of Fermi Large Area Telescope (LAT) gamma-ray data. We detected radio pulsations at 2 GHz, visible for <10% of ~4.5-hrs of observations using the Green Bank Telescope (GBT). Observations at 5 GHz with the GBT and at several lower frequencies with Parkes, Nancay, and the Giant Metrewave Radio Telescope resulted in non-detections. We also report the faint detection of a steep spectrum continuum radio source (0.1 mJy at 5 GHz) in interferometric imaging observations with the Jansky Very Large Array. These detections demonstrate that PSR J1311-3430, is not radio quiet and provides additional evidence that the radio beaming fraction of millisecond pulsars is very large. The radio detection yields a distance estimate of 1.4 kpc for the system, yielding a gamma-ray efficiency of 30%, typical of LAT-detected MSPs. We see apparent excess delay in the radio pulsar as the pulsar appears from eclipse and we speculate on possible mechanisms for the non-detections of the pulse at other orbital phases and observing frequencies.

Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

We describe the first X-ray observations of binary millisecond pulsars PSRs J0023+0923, J1810+1744, J2215+5135, and J2256-1024. All four are Fermi gamma-ray sources and three are 'black-widow' pulsars, with companions of mass < 0.1 solar masses. Data were taken using the Chandra X-Ray Observatory and covered a full binary orbit for each pulsar. Two pulsars, PSRs J2215+5135 and J2256-1024, show significant orbital variability and X-ray flux minima at the times of eclipses observed at radio wavelengths. This phenomenon is consistent with intrabinary shock emission characteristic of black-widow pulsars. The other two pulsars, PSRs J0023+0923 and J1810+1744, do not demonstrate significant variability, but are fainter than the other two sources. Spectral fits yield power-law indices that range from 1.4 to 2.3 and blackbody temperatures in the hundreds of eV. The spectrum for PSR J2215+5135 shows a significant hard X-ray component (41% of counts are above 2 keV), which is additional evidence for the presence of intrabinary shock emission.

We describe an "active" antenna system for HF/VHF (long wavelength) radio astronomy that has been successfully deployed 256-fold as the first station (LWA1) of the planned Long Wavelength Array. The antenna system, consisting of crossed dipoles, an active balun/preamp, a support structure, and a ground screen has been shown to successfully operate over at least the band from 20 MHz (15 m wavelength) to 80 MHz (3.75 m wavelength) with a noise figure that is at least 6 dB better than the Galactic background emission noise temperature over that band. Thus, the goal to design and construct a compact, inexpensive, rugged, and easily assembled antenna system that can be deployed many-fold to form numerous large individual "stations" for the purpose of building a large, long wavelength synthesis array telescope for radio astronomical and ionospheric observations was met.

The Long Wavelength Array Software Library (LSL) is a Python module that provides a collection of utilities to analyze and export data collected at the first station of the Long Wavelength Array, LWA1. Due to the nature of the data format and large-N ($\gtrsim$100 inputs) challenges faced by the LWA, currently available software packages are not suited to process the data. Using tools provided by LSL, observers can read in the raw LWA1 data, synthesize a filter bank, and apply incoherent de-dispersion to the data. The extensible nature of LSL also makes it an ideal tool for building data analysis pipelines and applying the methods to other low frequency arrays.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

We report the discovery of PSR J1838-0537, a gamma-ray pulsar found through a blind search of data from the Fermi Large Area Telescope (LAT). The pulsar has a spin frequency of 6.9 Hz and a frequency derivative of -2.2e-11 Hz/s, implying a young characteristic age of 4970 years and a large spin-down power of 5.9e36 erg/s. Follow-up observations with radio telescopes detected no pulsations, thus PSR J1838-0537 appears radio-quiet as viewed from Earth. In September 2009 the pulsar suffered the largest glitch so far seen in any gamma-ray-only pulsar, causing a relative increase in spin frequency of about 5.5e-6. After the glitch, during a putative recovery period, the timing analysis is complicated by the sparsity of the LAT photon data, the weakness of the pulsations, and the reduction in average exposure from a coincidental, contemporaneous change in the LAT's sky-survey observing pattern. The pulsar's sky position is coincident with the spatially extended TeV source HESS J1841-055 detected by the High Energy Stereoscopic System (H.E.S.S.). The inferred energetics suggest that HESS J1841-055 contains a pulsar wind nebula powered by the pulsar.

The first station of the Long Wavelength Array (LWA1) was completed in April 2011 and is currently performing observations resulting from its first call for proposals in addition to a continuing program of commissioning and characterization observations. The instrument consists of 258 dual-polarization dipoles, which are digitized and combined into beams. Four independently-steerable dual-polarization beams are available, each with two "tunings" of 16 MHz bandwidth that can be independently tuned to any frequency between 10 MHz and 88 MHz. The system equivalent flux density for zenith pointing is ~3 kJy and is approximately independent of frequency; this corresponds to a sensitivity of ~5 Jy/beam (5sigma, 1 s); making it one of the most sensitive meter-wavelength radio telescopes. LWA1 also has two "transient buffer" modes which allow coherent recording from all dipoles simultaneously, providing instantaneous all-sky field of view. LWA1 provides versatile and unique new capabilities for Galactic science, pulsar science, solar and planetary science, space weather, cosmology, and searches for astrophysical transients. Results from LWA1 will detect or tightly constrain the presence of hot Jupiters within 50 parsecs of Earth. LWA1 will provide excellent resolution in frequency and in time to examine phenomena such as solar bursts, and pulsars over a 4:1 frequency range that includes the poorly understood turnover and steep-spectrum regimes. Observations to date have proven LWA1's potential for pulsar observing, and just a few seconds with the completed 256-dipole LWA1 provide the most sensitive images of the sky at 23 MHz obtained yet. We are operating LWA1 as an open skies radio observatory, offering ~2000 beam-hours per year to the general community.

We have developed a compact cryogenic system with a pulse tube refrigerator and a coaxial heat exchanger. This liquefaction-purification system not only saves the cooling power used to reach high gaseous recirculation rate, but also reduces the impurity level with high speed. The heat exchanger operates with an efficiency of 99%, which indicates the possibility for fast xenon gas recirculation in a highpressurized large-scale xenon storage with much less thermal losses.

We present a summary of the Fermi Pulsar Search Consortium (PSC), an international collaboration of radio astronomers and members of the Large Area Telescope (LAT) collaboration, whose goal is to organize radio follow-up observations of Fermi pulsars and pulsar candidates among the LAT gamma-ray source population. The PSC includes pulsar observers with expertise using the world's largest radio telescopes that together cover the full sky. We have performed very deep observations of all 35 pulsars discovered in blind frequency searches of the LAT data, resulting in the discovery of radio pulsations from four of them. We have also searched over 300 LAT gamma-ray sources that do not have strong associations with known gamma-ray emitting source classes and have pulsar-like spectra and variability characteristics. These searches have led to the discovery of a total of 43 new radio millisecond pulsars (MSPs) and four normal pulsars. These discoveries greatly increase the known population of MSPs in the Galactic disk, more than double the known population of so-called `black widow' pulsars, and contain many promising candidates for inclusion in pulsar timing arrays.

We report the discovery of the millisecond pulsar PSR J2043+1711 in a search of a Fermi Large Area Telescope (LAT) source with no known associations, with the Nancay Radio Telescope. The new pulsar, confirmed with the Green Bank Telescope, has a spin period of 2.38 ms, is relatively nearby (d <~ 2 kpc), and is in a 1.48 day orbit around a low mass companion, probably a He-type white dwarf. Pulsed gamma-ray emission was detected in the data recorded by the Fermi LAT. The gamma-ray light curve and spectral properties are typical of other gamma-ray millisecond pulsars seen with Fermi. X-ray observations of the pulsar with Suzaku and the Swift/XRT yielded no detection. At 1.4 GHz we observe strong flux density variations because of interstellar diffractive scintillation, however a sharp peak can be observed at this frequency during bright scintillation states. At 327 MHz the pulsar is detected with a much higher signal-to-noise ratio and its flux density is far more steady. However, at that frequency the Arecibo instrumentation cannot yet fully resolve the pulse profile. Despite that, our pulse time-of-arrival measurements have a post-fit residual rms of 2 \mus. This and the expected stability of this system has made PSR J2043+1711 one of the first new Fermi-selected millisecond pulsars to be added to pulsar gravitational wave timing arrays. It has also allowed a significant measurement of relativistic delays in the times of arrival of the pulses due to the curvature of space-time near the companion, but not yet with enough precision to derive useful masses for the pulsar and the companion. A mass for the pulsar between 1.7 and 2.0 solar masses can be derived if a standard millisecond pulsar formation model is assumed. In this article we also present a comprehensive summary of pulsar searches in Fermi LAT sources with the Nancay Radio Telescope to date.

We have discovered five millisecond pulsars (MSPs) in a survey of 14 unidentified Fermi-LAT sources in the southern sky using the Parkes radio telescope. PSRs J0101-6422, J1514-4946, and J1902-5105 reside in binaries, while PSRs J1658-5324 and J1747-4036 are isolated. Using an ephemeris derived from timing observations of PSR J0101-6422 (P =2.57 ms, DM=12 pc cm-3), we have detected {\gamma}-ray pulsations and measured its proper motion. Its {\gamma}-ray spectrum (a power law of {\Gamma} = 0.9 with a cutoff at 1.6 GeV) and efficiency are typical of other MSPs, but its radio and {\gamma}-ray light curves challenge simple geometric models of emission. The high success rate of this survey-enabled by selecting {\gamma}-ray sources based on their detailed spectral characteristics-and other similarly successful searches indicate that a substantial fraction of the local population of MSPs may soon be known.

In a radio search with the Green Bank Telescope of three unidentified low Galactic latitude Fermi-LAT sources, we have discovered the middle-aged pulsar J2030+3641, associated with 1FGL J2030.0+3641 (2FGL J2030.0+3640). Following the detection of gamma-ray pulsations using a radio ephemeris, we have obtained a phase-coherent timing solution based on gamma-ray and radio pulse arrival times that spans the entire Fermi mission. With a rotation period of 0.2 s, spin-down luminosity of 3e34 erg/s, and characteristic age of 0.5 Myr, PSR J2030+3641 is a middle-aged neutron star with spin parameters similar to those of the exceedingly gamma-ray-bright and radio-undetected Geminga. Its gamma-ray flux is 1% that of Geminga, primarily because of its much larger distance, as suggested by the large integrated column density of free electrons, DM=246 pc/cc. We fit the gamma-ray light curve, along with limited radio polarimetric constraints, to four geometrical models of magnetospheric emission, and while none of the fits have high significance some are encouraging and suggest that further refinements of these models may be worthwhile. We argue that not many more non-millisecond radio pulsars may be detected along the Galactic plane that are responsible for LAT sources, but that modified methods to search for gamma-ray pulsations should be productive -- PSR J2030+3641 would have been found blindly in gamma rays if only >0.8 GeV photons had been considered, owing to its relatively flat spectrum and location in a region of high soft background.

We report the detection of pulsed gamma-ray emission from the fast millisecond pulsars (MSPs) B1937+21 (also known as J1939+2134) and B1957+20 (J1959+2048) using 18 months of survey data recorded by the \emph{Fermi} Large Area Telescope (LAT) and timing solutions based on radio observations conducted at the Westerbork and Nan\c{c}ay radio telescopes. In addition, we analyzed archival \emph{RXTE} and \emph{XMM-Newton} X-ray data for the two MSPs, confirming the X-ray emission properties of PSR B1937+21 and finding evidence ($\sim 4\sigma$) for pulsed emission from PSR B1957+20 for the first time. In both cases the gamma-ray emission profile is characterized by two peaks separated by half a rotation and are in close alignment with components observed in radio and X-rays. These two pulsars join PSRs J0034-0534 and J2214+3000 to form an emerging class of gamma-ray MSPs with phase-aligned peaks in different energy bands. The modeling of the radio and gamma-ray emission profiles suggests co-located emission regions in the outer magnetosphere.

We report the discovery of nine previously unknown gamma-ray pulsars in a blind search of data from the Fermi Large Area Telescope (LAT). The pulsars were found with a novel hierarchical search method originally developed for detecting continuous gravitational waves from rapidly rotating neutron stars. Designed to find isolated pulsars spinning at up to kHz frequencies, the new method is computationally efficient, and incorporates several advances, including a metric-based gridding of the search parameter space (frequency, frequency derivative and sky location) and the use of photon probability weights. The nine pulsars have spin frequencies between 3 and 12 Hz, and characteristic ages ranging from 17 kyr to 3 Myr. Two of them, PSRs J1803-2149 and J2111+4606, are young and energetic Galactic-plane pulsars (spin-down power above 6e35 erg/s and ages below 100 kyr). The seven remaining pulsars, PSRs J0106+4855, J0622+3749, J1620-4927, J1746-3239, J2028+3332, J2030+4415, J2139+4716, are older and less energetic; two of them are located at higher Galactic latitudes (|b| > 10 deg). PSR J0106+4855 has the largest characteristic age (3 Myr) and the smallest surface magnetic field (2e11 G) of all LAT blind-search pulsars. PSR J2139+4716 has the lowest spin-down power (3e33 erg/s) among all non-recycled gamma-ray pulsars ever found. Despite extensive multi-frequency observations, only PSR J0106+4855 has detectable pulsations in the radio band. The other eight pulsars belong to the increasing population of radio-quiet gamma-ray pulsars.

We have searched for millisecond pulsars (MSPs) in two globular clusters detected by the Fermi Large Area Telescope. These clusters contained no known MSPs prior to their detections in gamma rays. The discovery of gamma ray emission from many MSPs and the prevalence of MSPs in globular clusters points to a population of MSPs as the likely source of the detected GeV emission, directing our search for new cluster MSPs. We observed NGC 6652 and NGC 6388 at 2 GHz with the Green Bank Ultimate Pulsar Processing Instrument pulsar backend at the Green Bank Telescope using the coherent dedispersion mode. We have discovered one MSP in the gamma ray error circle of NGC 6652. This pulsar is interesting because, while positionally coincident with the GC, it has a much lower dispersion measure than expected from the NE2001 galactic free electron density model. It is unclear whether the MSP is a foreground pulsar or a cluster member, and whether the pulsar, cluster, or both, is responsible for the gamma ray emission. Timing the MSP will give the pulsar position and a solid identification of the pulsar as a cluster member if it is within a few core radii of the cluster center, as well as the opportunity to search for gamma ray pulsations and determine the origin of the GeV emission.

AXTAR is a NASA MIDEX mission concept for X-ray timing of compact objects that combines very large collecting area, broadband spectral coverage, high time resolution, highly flexible scheduling, and an ability to respond promptly to time-critical targets of opportunity. It is optimized for submillisecond timing of bright Galactic X-ray sources in order to study phenomena at the natural time scales of neutron star surfaces and black hole event horizons, thus probing the physics of ultradense matter, strongly curved spacetimes, and intense magnetic fields. AXTAR's main instrument is a collimated, thick Si pixel detector with 2-50 keV coverage and over 3 square meters effective area. For timing observations of accreting neutron stars and black holes, AXTAR provides at least a factor of five improvement in sensitivity over the RXTE PCA. AXTAR also carries a sensitive sky monitor that acts as a trigger for pointed observations of X-ray transients in addition to providing high duty cycle monitoring of the X-ray sky. We review the science goals and design choices that face a next generation timing mission. We then describe the technical concept for AXTAR and summarize a preliminary mission design study at the NASA/MSFC Advanced Concepts Office.

We report the probable identification of the X-ray counterpart to the gamma-ray pulsar PSR J2021+4026 using imaging with the Chandra X-ray Observatory ACIS and timing analysis with the Fermi satellite. Given the statistical and systematic errors, the positions determined by both satellites are coincident. The X-ray source position is R.A. 20h21m30.733s, Decl. +40 deg 26 min 46.04sec (J2000) with an estimated uncertainty of 1.3 arsec combined statistical and systematic error. Moreover, both the X-ray to gamma-ray and the X-ray to optical flux ratios are sensible assuming a neutron star origin for the X-ray flux. The X-ray source has no cataloged infrared-to-visible counterpart and, through new observations, we set upper limits to its optical emission of i' >23.0 mag and r' > 25.2mag. The source exhibits an X-ray spectrum with most likely both a powerlaw and a thermal component. We also report on the X-ray and visible light properties of the 43 other sources detected in our Chandra observation.

Since the launch of the Fermi satellite, the number of known gamma-ray pulsars has increased tenfold. Most gamma-ray detected pulsars are young and energetic, and many are associated with TeV sources. PSR J1357-6429 is a high spin-down power pulsar (Edot = 3.1 * 10^36 erg/s), discovered during the Parkes multibeam survey of the Galactic plane, with significant timing noise typical of very young pulsars. In the very-high-energy domain, H.E.S.S. has reported the detection of the extended source HESS J1356-645 (intrinsic Gaussian width of 12') whose centroid lies 7' from PSR J1357-6429. Using a rotational ephemeris obtained with 74 observations made with the Parkes telescope at 1.4 GHz, we phase-fold more than two years of gamma-ray data acquired by Fermi-LAT as well as those collected with XMM-Newton, and perform gamma-ray spectral modeling. Significant gamma and X-ray pulsations are detected from PSR J1357-6429. The light curve in both bands shows one broad peak. Gamma-ray spectral analysis of the pulsed emission suggests that it is well described by a simple power-law of index 1.5 +/- 0.3stat +/- 0.3syst with an exponential cut-off at 0.8 +/- 0.3stat +/- 0.3syst GeV and an integral photon flux above 100 MeV of (6.5 +/- 1.6stat +/- 2.3syst) * 10^-8 cm^-2 s^-1. The X-ray spectra obtained from the new data provide results consistent with those reported by Zavlin (2007). Upper limits on the gamma-ray emission from its potential pulsar wind nebula (PWN) are also reported. Assuming a distance of 2.4 kpc, the Fermi LAT energy flux yields a gamma-ray luminosity for PSR J1357-6429 of L_gamma = (2.13 +/- 0.25stat +/- 0.83syst) * 10^34 erg/s, consistent with an L_gamma \propto sqrt(Edot) relationship. The Fermi non-detection of the pulsar wind nebula associated with HESS J1356-645 provides new constraints on the electron population responsible for the extended TeV emission.

One of the main results of the Fermi Gamma-Ray Space Telescope is the discovery of {\gamma}-ray selected pulsars. The high magnetic field pulsar, PSR J0007+7303 in CTA1, was the first ever to be discovered through its {\gamma}-ray pulsations. Based on analysis of 2 years of LAT survey data, we report on the discovery of {\gamma}-ray emission in the off-pulse phase interval at the ~ 6{\sigma} level. The flux from this emission in the energy range E \geq 100 MeV is F_100 = (1.73\pm0.40)\times10^(-8) photons/cm^2/s and is best fitted by a power law with a photon index of {\Gamma} = 2.54\pm0.14. The pulsed {\gamma}-ray flux in the same energy range is F_100 = (3.95\pm0.07)\times10^(-7) photons/cm^2/s and is best fitted by an exponentially-cutoff power-law spectrum with a photon index of {\Gamma} = 1.41 \pm 0.23 and a cutoff energy E_c = 4.04 \pm 0.20 GeV. We find no flux variability neither at the 2009 May glitch nor in the long term behavior. We model the {\gamma}-ray light curve with two high-altitude emission models, the outer gap and slot gap, and find that the model that best fits the data depends strongly on the assumed origin of the off-pulse emission. Both models favor a large angle between the magnetic axis and observer line of sight, consistent with the nondetection of radio emission being a geometrical effect. Finally we discuss how the LAT results bear on the understanding of the cooling of this neutron star.

LS I +61 303 and LS 5039 are exceptionally rare examples of HMXBs with MeV-TeV emission, making them two of only five known or proposed "gamma-ray binaries". There has been disagreement within the literature over whether these systems are microquasars, with stellar winds accreting onto a compact object to produce high energy emission and relativistic jets, or whether their emission properties might be better explained by a relativistic pulsar wind colliding with the stellar wind. Here we present an attempt to detect radio pulsars in both systems with the Green Bank Telescope. The upper limits of flux density are between 4.1-14.5 uJy, and we discuss the null results of the search. Our spherically symmetric model of the wind of LS 5039 demonstrates that any pulsar emission will be strongly absorbed by the dense wind unless there is an evacuated region formed by a relativistic colliding wind shock. LS I +61 303 contains a rapidly rotating Be star whose wind is concentrated near the stellar equator. As long as the pulsar is not eclipsed by the circumstellar disk or viewed through the densest wind regions, detecting pulsed emission may be possible during part of the orbit.

We report on the full analysis of a Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) monitoring of the {\gamma}-ray binary system LS I +61 303. The data set covers 42 contiguous cycles of the system orbital motion. Analyzing this X-ray monitoring dataset, the largest to date for this source, we report on the variability of the orbital profile and the spectral distribution, and provide strong evidence for an anti-correlation between flux and spectral index (the higher the flux, the harder the spectral index). Furthermore, we present the analysis of two newly discovered ks-timescale flares, which present significant variability also on shorter timescales, and tend to occur at orbital phases between 0.6-0.9. However, a detailed timing analysis of the flares does not show any coherent or quasi-coherent (QPO) structure in their power spectra. We also investigated the possible appearance of the radio super-orbital modulation at X-rays energies, but we could not unambiguously detect such modulation in the system flux history, nor in the evolution of its orbital modulation fraction.

We report the discovery of two millisecond pulsars in a search for radio pulsations at the positions of \emph{Fermi Large Area Telescope} sources with no previously known counterparts, using the Nan\c{c}ay radio telescope. The two millisecond pulsars, PSRs J2017+0603 and J2302+4442, have rotational periods of 2.896 and 5.192 ms and are both in binary systems with low-eccentricity orbits and orbital periods of 2.2 and 125.9 days respectively, suggesting long recycling processes. Gamma-ray pulsations were subsequently detected for both objects, indicating that they power the associated \emph{Fermi} sources in which they were found. The gamma-ray light curves and spectral properties are similar to those of previously-detected gamma-ray millisecond pulsars. Detailed modeling of the observed radio and gamma-ray light curves shows that the gamma-ray emission seems to originate at high altitudes in their magnetospheres. Additionally, X-ray observations revealed the presence of an X-ray source at the position of PSR J2302+4442, consistent with thermal emission from a neutron star. These discoveries along with the numerous detections of radio-loud millisecond pulsars in gamma rays suggest that many \emph{Fermi} sources with no known counterpart could be unknown millisecond pulsars.

We present precise phase-connected pulse timing solutions for 16 gamma-ray-selected pulsars recently discovered using the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope plus one very faint radio pulsar (PSR J1124-5916) that is more effectively timed with the LAT. We describe the analysis techniques including a maximum likelihood method for determining pulse times of arrival from unbinned photon data. A major result of this work is improved position determinations, which are crucial for multi-wavelength follow up. For most of the pulsars, we overlay the timing localizations on X-ray images from Swift and describe the status of X-ray counterpart associations. We report glitches measured in PSRs J0007+7303, J1124-5916, and J1813-1246. We analyze a new 20 ks Chandra ACIS observation of PSR J0633+0632 that reveals an arcminute-scale X-ray nebula extending to the south of the pulsar. We were also able to precisely localize the X-ray point source counterpart to the pulsar and find a spectrum that can be described by an absorbed blackbody or neutron star atmosphere with a hard powerlaw component. Another Chandra ACIS image of PSR J1732-3131 reveals a faint X-ray point source at a location consistent with the timing position of the pulsar. Finally, we present a compilation of new and archival searches for radio pulsations from each of the gamma-ray-selected pulsars as well as a new Parkes radio observation of PSR J1124-5916 to establish the gamma-ray to radio phase offset.

Using the Parkes radio telescope we have carried out deep observations of eleven unassociated gamma-ray sources. Periodicity searches of these data have discovered two millisecond pulsars, PSR J1103-5403 (1FGL J1103.9-5355) and PSR J2241-5236 (1FGL J2241.9-5236), and a long period pulsar, PSR J1604-44 (1FGL J1604.7-4443). In addition we searched for but did not detect any radio pulsations from six gammaray pulsars discovered by the Fermi satellite to a level of - 0.04 mJy (for pulsars with a 10% duty cycle). Timing of the millisecond pulsar PSR J1103-5403 has shown that its position is 9' from the centroid of the gamma-ray source. Since these observations were carried out, independent evidence has shown that 1FGL J1103.9-5355 is associated with the flat spectrum radio source PKS 1101-536. It appears certain that the pulsar is not associated with the gamma-ray source, despite the seemingly low probability of a chance detection of a radio millisecond pulsar. We consider that PSR J1604-44 is a chance discovery of a weak, long period pulsar and is unlikely to be associated with 1FGL J1604.7-4443. PSR J2241-5236 has a spin period of 2.2 ms and orbits a very low mass companion with a 3.5 hour orbital period. The relatively high flux density and low dispersion measure of PSR J2241-5236 makes it an excellent candidate for high precision timing experiments. The gamma-rays of 1FGL J2241.9-5236 have a spectrum that is well modelled by a power law with exponential cutoff, and phasebinning with the radio ephemeris results in a multi-peaked gamma-ray pulse profile. Observations with Chandra have identified a coincident X-ray source within 0.1" of the position of the pulsar obtained by radio timing

We have used the Green Bank Telescope at 350MHz to search 50 faint, unidentified Fermi Gamma-ray sources for radio pulsations. So far, these searches have resulted in the discovery of 10 millisecond pulsars, which are plausible counterparts to these unidentified Fermi sources. Here we briefly describe this survey and the characteristics of the newly discovered MSPs.

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

We present an overview of precise pulsar timing using data from the Large Area Telescope (LAT) on Fermi. We describe the analysis techniques including a maximum likelihood method for determining pulse times of arrival from unbinned photon data. In addition to determining the spindown behavior of the pulsars and detecting glitches and timing noise, such timing analyses allow the precise determination of the pulsar position, thus enabling detailed multiwavelength follow up.

We searched for radio pulsars in 25 of the non-variable, unassociated sources in the Fermi LAT Bright Source List with the Green Bank Telescope at 820 MHz. We report the discovery of three radio and gamma-ray millisecond pulsars (MSPs) from a high Galactic latitude subset of these sources. All of the pulsars are in binary systems, which would have made them virtually impossible to detect in blind gamma-ray pulsation searches. They seem to be relatively normal, nearby (<=2 kpc) millisecond pulsars. These observations, in combination with the Fermi detection of gamma-rays from other known radio MSPs, imply that most, if not all, radio MSPs are efficient gamma-ray producers. The gamma-ray spectra of the pulsars are power-law in nature with exponential cutoffs at a few GeV, as has been found with most other pulsars. The MSPs have all been detected as X-ray point sources. Their soft X-ray luminosities of ~10^{30-31} erg/s are typical of the rare radio MSPs seen in X-rays.

This paper presents a search for radio transients at a frequency of 73.8 MHz (4 m wavelength) using the all-sky imaging capabilities of the Long Wavelength Demonstrator Array (LWDA). The LWDA was a 16-dipole phased array telescope, located on the site of the Very Large Array in New Mexico. The field of view of the individual dipoles was essentially the entire sky, and the number of dipoles was sufficiently small that a simple software correlator could be used to make all-sky images. From 2006 October to 2007 February, we conducted an all-sky transient search program, acquiring a total of 106 hr of data; the time sampling varied, being 5 minutes at the start of the program and improving to 2 minutes by the end of the program. We were able to detect solar flares, and in a special-purpose mode, radio reflections from ionized meteor trails during the 2006 Leonid meteor shower. We detected no transients originating outside of the solar system above a flux density limit of 500 Jy, equivalent to a limit of no more than about 10^{-2} events/yr/deg^2, having a pulse energy density >~ 1.5 x 10^{-20} J/m^2/Hz at 73.8 MHz for pulse widths of about 300 s. This event rate is comparable to that determined from previous all-sky transient searches, but at a lower frequency than most previous all-sky searches. We believe that the LWDA illustrates how an all-sky imaging mode could be a useful operational model for low-frequency instruments such as the Low Frequency Array, the Long Wavelength Array station, the low-frequency component of the Square Kilometre Array, and potentially the Lunar Radio Array.

Pulsars PSR J0248+6021 (rotation period P=217 ms and spin-down power Edot = 2.13E35 erg/s) and PSR J2240+5832 (P=140 ms, Edot = 2.12E35 erg/s) were discovered in 1997 with the Nancay radio telescope during a northern Galactic plane survey, using the Navy-Berkeley Pulsar Processor (NBPP) filter bank. GeV gamma-ray pulsations from both were discovered using the Fermi Large Area Telescope. Twelve years of radio and polarization data allow detailed investigations. The two pulsars resemble each other both in radio and in gamma-ray data. Both are rare in having a single gamma-ray pulse offset far from the radio peak. The high dispersion measure for PSR J0248+6021 (DM = 370 pc cm^-3) is most likely due to its being within the dense, giant HII region W5 in the Perseus arm at a distance of 2 kpc, not beyond the edge of the Galaxy as obtained from models of average electron distributions. Its high transverse velocity and the low magnetic field along the line-of-sight favor this small distance. Neither gamma-ray, X-ray, nor optical data yield evidence for a pulsar wind nebula surrounding PSR J0248+6021. The gamma-ray luminosity for PSR J0248+6021 is L_ gamma = (1.4 \pm 0.3)\times 10^34 erg/s. For PSR J2240+5832, we find either L_gamma = (7.9 \pm 5.2) \times 10^34 erg/s if the pulsar is in the Outer arm, or L_gamma = (2.2 \pm 1.7) \times 10^34 erg/s for the Perseus arm. These luminosities are consistent with an L_gamma ~ sqrt(Edot) rule. Comparison of the gamma-ray pulse profiles with model predictions, including the constraints obtained from radio polarization data, favor emission in the far magnetosphere. These two pulsars differ mainly in their inclination angles and acceleration gap widths, which in turn explains the observed differences in the gamma-ray peak widths.

Prompted by the Fermi LAT discovery of a radio-quiet gamma-ray pulsar inside the CTA 1 supernova remnant, we obtained a 130 ks XMM-Newton observation to assess the timing behavior of this pulsar. Exploiting both the unprecedented photon harvest and the contemporary Fermi LAT timing measurements, a 4.7 sigma single peak pulsation is detected, making PSR J0007+7303 the second example, after Geminga, of a radio-quiet gamma-ray pulsar also seen to pulsate in X-rays. Phase-resolved spectroscopy shows that the off-pulse portion of the light curve is dominated by a power-law, non-thermal spectrum, while the X-ray peak emission appears to be mainly of thermal origin, probably from a polar cap heated by magnetospheric return currents, pointing to a hot spot varying throughout the pulsar rotation.

In late 2008, the quasi-persistent neutron star X-ray transient and eclipsing binary EXO 0748-676 started a transition from outburst to quiescence, after it had been actively accreting for more than 24 years. In a previous work, we discussed Chandra and Swift observations obtained during the first five months after this transition. Here, we report on further X-ray observations of EXO 0748-676, extending the quiescent monitoring to 1.6 years. Chandra and XMM-Newton data reveal quiescent X-ray spectra composed of a soft, thermal component that is well-fitted by a neutron star atmosphere model. An additional hard powerlaw tail is detected that changes non-monotonically over time, contributing between 4 and 20 percent to the total unabsorbed 0.5-10 keV flux. The combined set of Chandra, XMM-Newton and Swift data reveals that the thermal bolometric luminosity fades from ~1E34 to 6E33 (D/7.4 kpc)^2 erg/s, whereas the inferred neutron star effective temperature decreases from ~124 to 109 eV. We interpret the observed decay as cooling of the neutron star crust and show that the fractional quiescent temperature change of EXO 0748-676 is markedly smaller than observed for three other neutron star X-ray binaries that underwent prolonged accretion outbursts.

LS I +61 303 is an exceptionally rare example of a high mass X-ray binary (HMXB) that also exhibits MeV-TeV emission, making it one of only a handful of "gamma-ray binaries". Here we present H-alpha spectra that show strong variability during the 26.5 day orbital period and over decadal time scales. We detect evidence of a spiral density wave in the Be circumstellar disk over part of the orbit. The H-alpha line profile also exhibits a dramatic emission burst shortly before apastron, observed as a redshifted shoulder in the line profile, as the compact source moves almost directly away from the observer. We investigate several possible origins for this red shoulder, including an accretion disk, mass transfer stream, and a compact pulsar wind nebula that forms via a shock between the Be star's wind and the relativistic pulsar wind.

We report on the analysis of Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) monitoring observations of the $\gamma$-ray binary system LS I +61 303, covering 35 full cycles of its orbital motion. This constitutes the largest continuous X-ray monitoring dataset analyzed to date for this source. Such an extended analysis allows us to report: a) the discovery of variability in the orbital profiles of the X- ray emission, b) the existence of a few (recent) short flares on top of the overall behavior typical of the source, which, given the PCA field-of-view, may or may not be associated with LS I +61 303, and c) the determination of the orbital periodicity using soft X-ray data alone.

The launch of the Fermi Gamma-ray Space Telescope has heralded a new era in the study of gamma-ray pulsars. The population of confirmed gamma-ray pulsars has gone from 6-7 to more than 60, and the superb sensitivity of the Large Area Telescope (LAT) on Fermi has allowed the detailed study of their spectra and light curves. Twenty-four of these pulsars were discovered in blind searches of the gamma-ray data, and twenty-one of these are, at present, radio quiet, despite deep radio follow-up observations. In addition, millisecond pulsars have been confirmed as a class of gamma-ray emitters, both individually and collectively in globular clusters. Recently, radio searches in the direction of LAT sources with no likely counterparts have been highly productive, leading to the discovery of a large number of new millisecond pulsars. Taken together, these discoveries promise a great improvement in the understanding of the gamma-ray emission properties and Galactic population of pulsars. We summarize some of the results stemming from these newly-detected pulsars and their timing and multi-wavelength follow-up observations.

The Advanced X-ray Timing Array (AXTAR) is a mission concept for X-ray timing of compact objects that combines very large collecting area, broadband spectral coverage, high time resolution, highly flexible scheduling, and an ability to respond promptly to time-critical targets of opportunity. It is optimized for submillisecond timing of bright Galactic X-ray sources in order to study phenomena at the natural time scales of neutron star surfaces and black hole event horizons, thus probing the physics of ultradense matter, strongly curved spacetimes, and intense magnetic fields. AXTAR's main instrument, the Large Area Timing Array (LATA) is a collimated instrument with 2-50 keV coverage and over 3 square meters effective area. The LATA is made up of an array of supermodules that house 2-mm thick silicon pixel detectors. AXTAR will provide a significant improvement in effective area (a factor of 7 at 4 keV and a factor of 36 at 30 keV) over the RXTE PCA. AXTAR will also carry a sensitive Sky Monitor (SM) that acts as a trigger for pointed observations of X-ray transients in addition to providing high duty cycle monitoring of the X-ray sky. We review the science goals and technical concept for AXTAR and present results from a preliminary mission design study.

We report the discovery of eight gamma-ray pulsars in blind frequency searches using the LAT, onboard the Fermi Gamma-ray Space Telescope. Five of the eight pulsars are young (tau_c<100 kyr), energetic (Edot>10^36 erg/s), and located within the Galactic plane (|b|<3 deg). The remaining three are older, less energetic, and located off the plane. Five pulsars are associated with sources included in the LAT bright gamma-ray source list, but only one, PSR J1413-6205, is clearly associated with an EGRET source. PSR J1023-5746 has the smallest characteristic age (tau_c=4.6 kyr) and is the most energetic (Edot=1.1E37 erg/s) of all gamma-ray pulsars discovered so far in blind searches. PSRs J1957+5033 and J2055+25 have the largest characteristic ages (tau_c~1 Myr) and are the least energetic (Edot~5E33 erg/s) of the newly-discovered pulsars. We present the timing models, light curves, and detailed spectral parameters of the new pulsars. We used recent XMM observations to identify the counterpart of PSR J2055+25 as XMMU J205549.4+253959. In addition, publicly available archival Chandra X-ray data allowed us to identify the likely counterpart of PSR J1023-5746 as a faint, highly absorbed source, CXOU J102302.8-574606. The large X-ray absorption indicates that this could be among the most distant gamma-ray pulsars detected so far. PSR J1023-5746 is positionally coincident with the TeV source HESS J1023-575, located near the young stellar cluster Westerlund 2, while PSR J1954+2836 is coincident with a 4.3 sigma excess reported by Milagro at a median energy of 35 TeV. Deep radio follow-up observations of the eight pulsars resulted in no detections of pulsations and upper limits comparable to the faintest known radio pulsars, indicating that these can be included among the growing population of radio-quiet pulsars in our Galaxy being uncovered by the LAT, and currently numbering more than 20.

We report detection of strong circularly polarized emission from the transient bursting source GCRT J1745-3009 based on new analysis of 325 MHz GMRT observations conducted on 28 September 2003. We place 8 Solar radius as the upper limit on the size of the emission region. The implied high brightness temperature required for an object beyond 1 pc and the high fraction of circular polarization firmly establish the emission as coherent. Electron cyclotron or plasma emission from a highly subsolar magnetically dominated dwarf located less than 4 kpc away could have given rise to the GCRT radio emission.

Sixteen pulsars have been discovered so far in blind searches of photons collected with the Large Area Telescope on the Fermi Gamma-ray Space Telescope. We here report the discovery of radio pulsations from two of them. PSR J1741-2054, with period P=413ms, was detected in archival Parkes telescope data and subsequently has been detected at the Green Bank Telescope (GBT). Its received flux varies greatly due to interstellar scintillation and it has a very small dispersion measure of DM=4.7pc/cc, implying a distance of ~0.4kpc and possibly the smallest luminosity of any known radio pulsar. At this distance, for isotropic emission, its gamma-ray luminosity above 0.1GeV corresponds to 25% of the spin-down luminosity of dE/dt=9.4e33erg/s. The gamma-ray profile occupies 1/3 of pulse phase and has three closely-spaced peaks with the first peak lagging the radio pulse by delta=0.29P. We have also identified a soft Swift source that is the likely X-ray counterpart. In many respects PSR J1741-2054 resembles the Geminga pulsar. The second source, PSR J2032+4127, was detected at the GBT. It has P=143ms, and its DM=115pc/cc suggests a distance of ~3.6kpc, but we consider it likely that it is located within the Cyg OB2 stellar association at half that distance. The radio emission is nearly 100% linearly polarized, and the main radio peak precedes by delta=0.15P the first of two narrow gamma-ray peaks that are separated by Delta=0.50P. Faint, diffuse X-ray emission in a Chandra image is possibly its pulsar wind nebula. PSR J2032+4127 likely accounts for the EGRET source 3EG J2033+4118, while its pulsar wind is responsible for the formerly unidentified HEGRA source TeV J2032+4130.

We report our complete database of X-ray eclipse timings of the low mass X-ray binary EXO0748-676 observed by the Rossi X-Ray Timing Explorer (RXTE) satellite. As of this writing we have accumulated 443 full X-ray eclipses, 392 of which have been observed with the Proportional Counter Array on RXTE. These include both observations where an eclipse was specifically targeted and those eclipses found in the RXTE data archive. Eclipse cycle count has been maintained since the discovery of the EXO0748-676 system in February 1985. We describe our observing and analysis techniques for each eclipse and describe improvements we have made since the last compilation by Wolff et al. (2002). The principal result of this paper is the database containing the timing results from a seven-parameter fit to the X-ray light curve for each observed eclipse along with the associated errors in the fitted parameters. Based on the standard O-C analysis, EXO0748-676 has undergone four distinct orbital period epochs since its discovery. In addition, EXO0748-676 shows small-scale events in the O-C curve that are likely due to short-lived changes in the secondary star.

LS I +61 303 and LS 5039 are two of only a handful of known high mass X-ray binaries (HMXBs) that exhibit very high energy emission in the MeV-TeV range, and these "gamma-ray binaries" are of renewed interest due to the recent launch of the Fermi Gamma-ray Space Telescope. Here we present new radial velocities of both systems based on recent red and blue optical spectra. Both systems have somewhat discrepant orbital solutions available in the literature, and our new measurements result in improved orbital elements and resolve the disagreements. The improved geometry of each orbit will aid in studies of the high energy emission region near each source.

The quasi-persistent neutron star X-ray transient and eclipsing binary EXO 0748-676 recently started the transition to quiescence following an accretion outburst that lasted more than 24 years. We report on two Chandra and twelve Swift observations performed within five months after the end of the outburst. The Chandra spectrum is composed of a soft, thermal component that fits to a neutron star atmosphere model with kT^inf~0.12 keV, joined by a hard powerlaw tail that contributes ~20% of the total 0.5-10 keV unabsorbed flux. The combined Chandra/Swift data set reveals a relatively hot and luminous quiescent system with a temperature of kT^inf~0.11-0.13 keV and a bolometric thermal luminosity of ~8.1E33-1.6E34 (d/7.4 kpc)^2 erg/s. We discuss our results in the context of cooling neutron star models.

The characteristic physical timescales near stellar-mass compact objects are measured in milliseconds. These timescales -- the free-fall time, the fastest stable orbital period, and stellar spin periods -- encode the fundamental physical properties of compact objects: mass, radius, and angular momentum. The characteristic temperature of matter in the vicinity of neutron stars is such that the principal electromagnetic window into their realms is the X-ray band. Because of these connections to the fundamental properties of neutron stars, X-ray timing studies remain today the most direct means of probing their structure and dynamics. While current X-ray observatories have revealed many relevant and fascinating phenomena, they lack the sensitivity to fully exploit them to uncover the fundamental properties of compact objects and their extreme physics. With this white paper, we summarize and highlight the science opportunities that will accompany an order-of-magnitude improvement in X-ray timing sensitivity, a goal attainable in the coming decade.

We report the detection of a new transient radio source, GCRT J1742-3001, located ~1 degree from the Galactic center. The source was detected ten times from late 2006 to 2007 May in our 235 MHz transient monitoring program with the Giant Metrewave Radio Telescope (GMRT). The radio emission brightened in about one month, reaching a peak observed flux density of ~100 mJy on 2007 January 28, and decaying to ~50 mJy by 2007 May when our last monitoring observation was made. Two additional faint, isolated 235 MHz detections were made in mid-2006, also with the GMRT. GCRT J1742-3001 is unresolved at each epoch, with typical resolutions of ~20 arcsec x 10 arcsec. No polarization information is available from the observations. Based on nondetections in observations obtained simultaneously at 610 MHz, we deduce that the spectrum of GCRT J1742-3001 is very steep, with a spectral index less than about -2. Follow-up radio observations in 2007 September at 330 MHz and 1.4 GHz, and in 2008 February at 235 MHz yielded no detections. No X-ray counterpart is detected in a serendipitous observation obtained with the X-ray telescope aboard the Swift satellite during the peak of the radio emission in early 2007. We consider the possibilities that GCRT J1742-3001 is either a new member of an existing class of radio transients, or is representative of a new class having no associated X-ray emission.

Multi-wavelength studies at radio, infrared, optical, X-ray, and TeV wavelengths have discovered probable counterparts to many Galactic sources of GeV emission detected by EGRET. These include pulsar wind nebulae, high mass X-ray binaries, and mixed morphology supernova remnants. Here we provide an overview of the observational properties of Galactic sources which emit across 19 orders of magnitude in energy. We also present new observations of several sources.

AXTAR is an X-ray observatory mission concept, currently under study in the U.S., that combines very large collecting area, broadband spectral coverage, high time resolution, highly flexible scheduling, and an ability to respond promptly to time-critical targets of opportunity. It is optimized for submillisecond timing of bright Galactic X-ray sources in order to study phenomena at the natural time scales of neutron star surfaces and black hole event horizons, thus probing the physics of ultradense matter, strongly curved spacetimes, and intense magnetic fields. AXTAR's main instrument is a collimated, thick Si pixel detector with 2-50 keV coverage and 8 square meters collecting area. For timing observations of accreting neutron stars and black holes, AXTAR provides at least an order of magnitude improvement in sensitivity over both RXTE and Constellation-X. AXTAR also carries a sensitive sky monitor that acts as a trigger for pointed observations of X-ray transients and also provides continuous monitoring of the X-ray sky with 20 times the sensitivity of the RXTE ASM. AXTAR builds on detector and electronics technology previously developed for other applications and thus combines high technical readiness and well understood cost.

The radio sky is poorly sampled for rapidly varying transients because of the narrow field-of-view of most imaging radio telescopes at cm and shorter wavelengths. The emergence of sensitive long wavelength observations with intrinsically larger fields-of-view are changing this situation, as partly illustrated by our ongoing meter-wavelength monitoring observations and archival studies of the Galactic Center. In this search, we discovered a transient, bursting, radio source in the direction of the Galactic Center, GCRT J1745-3009, with extremely unusual properties. Its flux and rapid variability imply a brightness temperature >10^12 K if it is at a distance >70 pc, implying that it is a coherent emitter. I will discuss the discovery of the source and the subsequent re-detections, as well as searches for counterparts at other wavelengths, and several proposed models.

We present an optical/near-infrared search for a counterpart to the perplexing radio transient GCRT J1745-3009, a source located ~1 degree from the Galactic Center. Motivated by some similarities to radio bursts from nearby ultracool dwarfs, and by a distance upper limit of 70 pc for the emission to not violate the 1e12 K brightness temperature limit for incoherent radiation, we searched for a nearby star at the position of GCRT J1745-3009. We found only a single marginal candidate, limiting the presence of any late-type star to >1 kpc (spectral types earlier than M9), >200 pc (spectral types L and T0-T4), and >100 pc (spectral types T4-T7), thus severely restricting the possible local counterparts to GCRT J1745-3009. We also exclude any white dwarf within 1 kpc or a supergiant star out to the distance of the Galactic Center as possible counterparts. This implies that GCRT J1745-3009 likely requires a coherent emission process, although whether or not it reflects a new class of sources is unclear.

We perform a study of cosmic evolution with an equation of state parameter $\omega(t)=\omega_0+\omega_1(t\dot H/H)$ by selecting a phenomenological $\Lambda$ model of the form, $\dot\Lambda\sim H^3$. This simple proposition explains both linearly expanding and inflationary Universes with a single set of equations. We notice that the inflation leads to a scaling in the equation of state parameter, $\omega(t)$, and hence in equation of state. In this approach, one of its two parameters have been pin pointed and the other have been delineated. It has been possible to show a connection between dark energy and Higgs-Boson.

We present a 7 yr timing study of the 2.5 ms X-ray pulsar SAX J1808.4-3658, an X-ray transient with a recurrence time of ~2 yr, using data from the Rossi X-ray Timing Explorer covering 4 transient outbursts (1998-2005). We verify that the 401 Hz pulsation traces the spin frequency fundamental and not a harmonic. Substantial pulse shape variability, both stochastic and systematic, was observed during each outburst. Analysis of the systematic pulse shape changes suggests that, as an outburst dims, the X-ray "hot spot" on the pulsar surface drifts longitudinally and a second hot spot may appear. The overall pulse shape variability limits the ability to measure spin frequency evolution within a given X-ray outburst (and calls previous nudot measurements of this source into question), with typical upper limits of |nudot| < 2.5x10^{-14} Hz/s (2 sigma). However, combining data from all the outbursts shows with high (6 sigma) significance that the pulsar is undergoing long-term spin down at a rate nudot = (-5.6+/-2.0)x10^{-16} Hz/s, with most of the spin evolution occurring during X-ray quiescence. We discuss the possible contributions of magnetic propeller torques, magnetic dipole radiation, and gravitational radiation to the measured spin down, setting an upper limit of B < 1.5x10^8 G for the pulsar's surface dipole magnetic field and and Q/I < 5x10^{-9} for the fractional mass quadrupole moment. We also measured an orbital period derivative of Pdot = (3.5+/-0.2)x10^{-12} s/s. This surprising large Pdot is reminiscent of the large and quasi-cyclic orbital period variation observed in the so-called "black widow" millisecond radio pulsars, supporting speculation that SAX J1808.4-3658 may turn on as a radio pulsar during quiescence. In an appendix we derive an improved (0.15 arcsec) source position from optical data.

We report evidence of magnetic activity associated with the secondary star in the EXO 0748-676 low mass X-ray binary system. An analysis of a sequence of five consecutive X-ray eclipses observed during December 2003 with the RXTE satellite brings out a feature occurring during ingress we interpret as the X-ray photoelectric absorption shadow, as seen by an observer at Earth, of a plasma structure suspended above the surface of the secondary star. The light curve feature consists of an initial drop in count rate to near zero (the absorption shadow) with a very short rebound to a significant fraction of the pre-ingress count rate and then a final plunge to totality over a total time scale of ~25 s. The ingress feature persists for at least 5 consecutive orbital periods (a total of ~19 hr), and possibly up to 5 days in our data. Our data also show significant post-egress dipping during this eclipse sequence, unusual for this source, indicating possible secondary star mass ejection during this episode.

In this paper we have investigated the occurrence rate of high energetic(E>10 MeV) solar electron flares measured by IMP-8 spacecraft of NASA for solar cycle 21 (June, 1976 to August, 1986) first time by three different methods to detect periodicities accurately. Power-spectrum analysis confirms a periodicity ~155 days which is in consistent with the result of Chowdhury and Ray (2006), that "Rieger periodicity" was operated throughout the cycle 21 and it is independent on the energy of the electron fluxes.

GCRT J1745-3009 is a transient bursting radio source located in the direction of the Galactic center. It was discovered in a 330 MHz VLA observation from 2002 September 30--October 1 and subsequently rediscovered in a 330 MHz GMRT observation from 2003 September 28 by Hyman et al. Here we report a new radio detection of the source in 330 MHz GMRT data taken on 2004 March 20. The observed properties of the single burst detected differ significantly from those measured previously, suggesting that GCRT J1745-3009 was detected in a new physical state. The 2004 flux density was ~0.05 Jy, ~10x weaker than the single 2003 burst and ~30x weaker than the five bursts detected in 2002. We derive a very steep spectral index, alpha = -13.5 +/- 3.0, across the bandpass, a new result previously not detectable due to limitations in the analysis of the 2002 and 2003 observations. Also, the burst was detected for only ~2 min., in contrast to the 10 min. duration observed in the earlier bursts. Due to sparse sampling, only the single burst was detected in 2004, as in the 2003 epoch, and we cannot rule out additional undetected bursts that may have occurred with the same ~77 min. periodicity observed in 2002 or with a different periodicity. Considering our total time on source throughout both our archival and active monitoring campaigns, we estimate the source exhibits detectable bursting activity ~7% of the time.

As part of a large-scale search for coherent pulsations from LMXBs in the RXTE archive, we have completed a detailed series of searches for coherent pulsations of 4U 1820-30 -- an ultracompact LMXB with a binary period of 11.4 min, located in the globular cluster NGC6624. The short binary period implies any coherent signal would be highly accelerated, so we used phase modulation searches, orbital-parameter-fitting coherent searches, and standard acceleration searches to give significant sensitivity to millisecond pulsations. We searched, in four energy bands and at a range of luminosities, a total of 34 archival RXTE observations, 32 of which had on-source integration times longer than 10 ks, and some of which were made consecutively which allowed us to combine them. We found no pulsations. Using our phase modulation search technique, which we ran on all 34 observations, we have been able to place the first stringent (95% confidence) pulsed fraction limits of <~0.8% for all realistic spin frequencies (i.e. <~2kHz) and likely companion masses (0.02Msun <= M_c <= 0.3Msun). Using our orbital-parameter-fitting coherent search, which we ran on only 11 selected observations, we have placed a pulsed fraction limit of <~0.3% for spin frequencies <~1.25kHz and companion masses M_ <= 0.106Msun. By contrast, all five LMXBs known to emit coherent pulsations have intrinsic pulsed fractions in the range 3% to 7% when pulsations are observed. Hence, our searches rule out pulsations with significantly lower pulsed fractions than those already observed.

GCRT J1745-3009 is a transient bursting radio source located in the direction of the Galactic center, discovered in 330 MHz VLA observations from 2002 September 30--October 1 by Hyman et al. We have searched for bursting activity from GCRT J1745-3009 in nearly all of the available 330 MHz VLA observations of the Galactic center since 1989 as well as in 2003 GMRT observations. We report a new radio detection of the source in 330 MHz GMRT data taken on 2003 September 28. A single ~0.5 Jy burst was detected, approximately 3x weaker than the five bursts detected in 2002. Due to the sparse sampling of the 2003 observation, only the decay portion of a single burst was detected. We present additional evidence indicating that this burst is an isolated one, but we cannot completely rule out additional undetected bursts that may have occured with the same ~77 min. periodicity observed in 2002 or with a different periodicity. Assuming the peak emission was detected, the decay time of the burst, ~2 min, is consistent with that determined for the 2002 bursts. Based on the total time for which we have observations, we estimate that the source has a duty cycle of roughly 10%.

We have observed an unusually strong X-ray burst as a part of our regular eclipse timing observations of the low mass binary system EXO0748-676. The burst peak flux was 5.2x10^-8 ergs cm^-2 s^-1, approximately five times the normal peak X-ray burst flux observed from this source by RXTE. Spectral fits to the data strongly suggest that photospheric radius expansion occurred during the burst. In this Letter we examine the properties of this X-ray burst, which is the first example of a radius expansion burst from EXO0748-676 observed by RXTE. We find no evidence for coherent burst oscillations. Assuming that the peak burst luminosity is the Eddington luminosity for a 1.4 solar mass neutron star we derive a distance to EXO0748-676 of 7.7 kpc for a helium-dominated burst photosphere and 5.9 kpc for a hydrogen-dominated burst photosphere.

Transient astronomical sources are typically powered by compact objects and usually signify highly explosive or dynamic events. While radio astronomy has an impressive record of obtaining high time resolution observations, usually it is achieved in quite narrow fields-of-view. Consequently, the dynamic radio sky is poorly sampled, in contrast to the situation in the X- and gamma-ray bands in which wide-field instruments routinely detect transient sources. Here we report a new transient source, GCRT J1745-3009, detected in 2002 during a moderately wide-field radio transient monitoring program of the Galactic center (GC) region at 0.33 GHz. The characteristics of its bursts are unlike those known for any other class of radio transient. If located in or near the GC, its brightness temperature (~10^16 K) and the implied energy density within GCRT J1745-3009 vastly exceeds that observed in most other classes of radio astronomical sources, and is consistent with coherent emission processes rarely observed. We conclude that GCRT J1745-3009 is the first member of a new class of radio transient sources, the first of possibly many new classes to be identified through current and upcoming radio surveys.

Orbital period changes are an important diagnostic for understanding low mass X-ray binary (LMXB) accretion-induced angular momentum exchange and overall system evolution. We present our most recent results for the eclipse timing of the LMXB EXO0748-676. Since its discovery in 1985 it has apparently undergone three distinct orbital period "epochs", each characterized by a different orbital period than the previous epoch. We outline the orbital period behavior for EXO0748-676 over the past 18 years and discuss the implications of this behavior in light of current theoretical ideas for LMXB evolution.

As part of a large-scale search for coherent pulsations from LMXBs in the RXTE archive, we have completed a detailed series of searches for coherent pulsations of 4U 1820-30 -- an ultracompact LMXB with a binary period of 11.4 minutes located in the globular cluster NGC6624. The small binary period leads to a very high acceleration, so we used phase modulation searches as well as acceleration searches to give significant sensitivity to millisecond pulsations. We searched a total of 34 archival RXTE observations, 32 of which had an on-source integration time longer than 10 ks, and some of which were made consecutively which allowed us to combine them. While we found no pulsations, we have been able to place the first stringent (95% confidence) pulsed fraction limits of < 0.8% for all realistic spin frequencies (i.e. < 1 kHz) and likely companion masses (0.02 Msun < Mc < 0.3 Msun). By contrast all five LMXBs known to emit coherent pulsations have intrinsic pulsed fractions in the range 3% to 7% when pulsations are observed.

Narayan and Heyl (2002) have developed a theoretical framework to convert suitable upper limits on type I X-ray bursts from accreting black hole candidates (BHCs) into evidence for an event horizon. However, no appropriate observational limit exists in the literature. In this paper we survey 2101.2 ks of data from the Unconventional Stellar Aspect (USA) X-ray timing experiment and 5142 ks of data from the Rossi X-ray Timing Explorer (RXTE) experiment to obtain a formal constraint of this type. 1122 ks of neutron star data yield a population averaged mean burst rate of 1.69 x 10-5 bursts s-1 while 6081 ks of BHC data yield a 95% confidence level upper limit of 4.9 x 10-7 bursts s-1. This is the first published limit of this type for Black Hole Candidates. Applying the theoretical framework of Narayan and Heyl (2002) we calculate regions of unstable luminosity where the neutron stars are expected to burst and the BHCs would be expected to burst if they had a surface. In this unstable luminosity region 464 ks of neutron star data yield an averaged mean burst rate of 4.1 x 10-5 bursts s-1 and 1512 ks of BHC data yield a 95% confidence level upper limit of 2.0 x 10-6 bursts s-1, and a limit of > 10\sigma that BHCs do not burst with a rate similar to the rate of neutron stars in these unstable regions. This gives further evidence that BHCs do not have surfaces unless there is some new physics occurring on their surface.

We present an analysis of long term X-ray monitoring observations of Circinus X-1 (Cir X-1) made with four different instruments: Vela 5B, Ariel V ASM, Ginga ASM, and RXTE ASM, over the course of more than 30 years. We use Lomb-Scargle periodograms to search for the ~16.5 day orbital period of Cir X-1 in each of these data sets and from this derive a new orbital ephemeris based solely on X-ray measurements, which we compare to the previous ephemerides obtained from radio observations. We also use the Phase Dispersion Minimization (PDM) technique, as well as FFT analysis, to verify the periods obtained from periodograms. we obtain dynamic periodograms (both Lomb-Scargle and PDM) of Cir X-1 during the RXTE era, showing the period evolution of Cir X-1, and also displaying some unexplained discrete jumps in the location of the peak power.

Observations of the 297 s X-ray pulsar 1E1145.1-6141 with the Rossi X-Ray Timing Explorer have revealed its 14.4 d eccentric orbit around its B supergiant companion. The best-fit orbital elements are: P_orb = 14.365(2) d, a_x sin i=99.4(18) light-s, e=0.20(3). No eclipses are detected, indicating that the binary inclination is less than 55 deg.

We present 7 eclipse timings of the low mass X-ray binary EXO0748-676 obtained with the USA experiment during 1999-2000 as well as 122 eclipse timings obtained with RXTE during 1996-2000. According to our analysis, the mean orbital period has increased by ~8 ms between the pre-RXTE era (1985-1990) and the RXTE/USA era (1996-2000). This corresponds to an orbital period derivative of P(orb)/(dP(orb)/dt)~2x10^7 years. However, neither a constant orbital period derivative nor any other simple ephemeris provides an acceptable fit to the data: individual timings of eclipse centers have residuals of up to 15 or more seconds away from our derived smooth ephemerides. When we consider all published eclipse timing data including those presented here, a model that includes observational measurement error, cumulative period jitter, and underlying period evolution is found to be consistent with the timing data. We discuss several physical mechanisms for LMXB orbital evolution in an effort to account for the change in orbital period and the observed intrinsic jitter in the mid-eclipse times.

This paper reports X-ray spectral observations of a relatively nearby z = 0.048 BL Lacertae (BL Lac) object 1ES1959+650, which is a potential TeV emitter. The observations include 31 short pointings made by the Unconventional Stellar Aspect (USA) Experiment on board the Advanced Research and Global Observation Satellite (ARGOS), and 17 pointings by the PCA on board the Rossi X-ray Timing Explorer (RXTE). Most of these observations were spaced by less than 1 day. \es was detected by the ARGOS USA detector in the range 1-16 keV, and by the PCA in the 2-16 keV range but at different times. During the closely spaced RXTE observations beginning on 2000 July 28, an ending of one flare and a start of another are visible, associated with spectral changes, where the photon index Gamma ranges between ~ 1.4 and 1.7, and the spectrum is harder when the source is brighter. This implies that 1ES1959 is an XBL-type blazar, with the X-ray emission likely to originate via the synchrotron process. The USA observations reveal another flare that peaked on 2000 November 14 and doubled the flux within a few days, again associated with spectral changes of the same form. The spectral variability correlated with the flux and timing characteristics of this object that are similar to those of other nearby BL Lacs, and suggest relativistic beaming with a Doppler factor delta > 1.6 and magnetic fields of the order of a few mG. We also suggest that the steady component of the X-ray emission -- present in this object as well as in other XBLs -- may be due to the large-scale relativistic jet (such as measured by Chandra in many radio-loud AGN), but pointing very closely to our line of sight.

We present a linear diffusion model for the evolution of the double-peaked outburst in the transient source XTEJ1118+480. The model treats the two outbursts as episodic mass deposition at the outer radius of the disk followed by evolution of disk structure according to a diffusion process. We demonstrate that light curves with fast-rise, exponential decay profile are a general consequence of the diffusion process. Deconvolution of the light curve proves to be feasible and gives an input function specifying mass deposition at the outer disk edge as well as the total mass of the disk, both as functions of time. The derived evolution of total disk mass can be correlated with the observed evolution of the ~0.1 Hz QPO in the source reported in Wood et al. (2000).

We report on timing and spectral observations of the 2000 outburst of XTE J1550-564 made by the Unconventional Stellar Aspect (USA) Experiment on board the Advanced Research and Global Observation Satellite (ARGOS). We observe a low-frequency quasi-periodic oscillation (LFQPO) with a centroid frequency that tends to increase with increasing flux and a fractional rms amplitude which is correlated with the hardness ratio. The evolution of the hardness ratio (4--16 keV/1--4 keV) with time and source flux is examined. The hardness-intensity diagram (HID) shows a cyclical movement in the clockwise direction and possibly indicates the presence of two independent accretion flows. We observe a relationship between the USA 4--16 keV count rate and radio observations and discuss this in the context of previously observed correlations between X-ray, radio, optical and IR data. We examine our results in the context of models invoking two accretion flows: a thin disk and a hot sub-Keplerian flow.

The USA experiment on ARGOS and RXTE have exensively observed the X-ray transient XTE J1118+480 during its recent outburst in 2000 April--June. We present detailed monitoring of the evolution of a low frequency QPO which drifts from 0.07 Hz to 0.15 Hz during the outburst. We examine possible correlations of the QPO frequency with the flux and spectral characteristics of the source, and compare this QPO to low frequency QPOs observed in other black hole candidates.

The USA Experiment is a new X-ray timing experiment with large collecting area and microsecond time resolution capable of conducting a broad program of studies of galactic X-ray binaries. USA is one of nine experiments aboard the Advanced Research and Global Observation Satellite which was launched February 23, 1999. USA is a collimated proportional counter X-ray telescope with about 1000 cm$^2$ of effective area per detector with two detectors sensitive to photons in the energy range 1-15 keV. A unique feature of USA is that photon events are time tagged by reference to an onboard GPS receiver allowing precise absolute time and location determination. We will present an overview of the USA instrument, capabilities, and scientific observing plan as well as the current status of the instrument.

The binary system, LSI+61 303, is unusual both because of the dramatic, periodic, radio outbursts, and because of its possible association with the 100 MeV gamma-ray source, 2CG135+01. We have performed simultaneous radio and Rossi X-ray Timing Explorer X-ray observations at eleven intervals over the 26.5 day orbit, and in addition searched for variability on timescales ranging from milliseconds to hours. We confirm the modulation of the X-ray emission on orbital timescales originally reported by Taylor et al. (1996), and in addition we find a significant offset between the peak of the X-ray and radio flux. We argue that based on these results, the most likely X-ray emission mechanism is inverse Compton scattering of stellar photons off of electrons accelerated at the shock boundary between the relativistic wind of a young pulsar and the Be star wind. In these observations we also detected 2 -- 150 keV flux from the nearby low-redshift quasar QSO~0241+622. Comparing these measurements to previous hard X-ray and gamma-ray observations of the region containing both LSI+61 303 and QSO~0241+622, it is clear that emission from the QSO dominates.

We present the results of a multiwavelength monitoring campaign targeting the gamma-ray source 2CG 135+1 in an attempt to confirm the association of this object with the radio/Be/X-ray binary system LSI +61o 303. The campaign included simultaneous radio, optical, infrared, and hard x-ray/gamma-ray observations carried out with a variety of instruments, covering (not continously) almost three binary cycles of LSI +61o 303 during the period April-July 1994. Three separate OSSE observations of the gamma-ray source were carried out, covering different phases of the radio lightcurve. Hard X-ray/gamma-ray emission was detected from the direction of 2CG 135+1 during the first of these OSSE observations. The signal to noise ratio of the OSSE observations was insufficient to establish a spectral or intensity correlation of the high-energy emission with simultaneous radio, optical and infrared emission of LSI +61o 303. We briefly discuss the theoretical implications of our observations.

LSI +61 303 is an exotic binary system consisting of a ~10 Msun B star and a compact object which is probably a neutron star. The system is associated with the interesting radio source GT0236+610 that exhibits bright radio outbursts with a period of 26.5 days. We report the results of continuous daily radio interferometric observations of GT0236+610 at 2.25 and 8.3 GHz from 1994 January to 1996 February. The observations cover 25 complete (and 3 partial) cycles with multiple observations each day. We detect substantial cycle-to-cycle variability of the radio emission characterized by a rapid onset of the radio flares followed by a more gradual decrease of the emission. We detect a systematic change of the radio spectral index alpha which typically becomes larger than zero at the onset of the radio outbursts. This behavior is suggestive of expansion of material initially optically thick to radio frequencies, indicating either that synchrotron or inverse Compton cooling are important or that the free-free optical depth to the source is rapidly changing. After two years of observations, we see only weak evidence for the proposed 4-year periodic modulation in the peak flux of the outbursts. We observe a secular trend in the outburst phases according the the best published ephemeris. This trend indicates either orbital period evolution, or a drift in outburst orbital phase in response to some other change in the system.

We report the results of a CGRO multi-instrument series of observations of the unidentified gamma-ray source 2CG 135+1 and of its possible counterpart, the peculiar radio source GT 0236+610 coincident with the Be star LSI+61 303. We monitored the time variable radio source GT 0236+610 during the CGRO observations. OSSE and COMPTEL data were obtained during the period May-June 1994, and BATSE data for the period April 1994-January 1995. We discuss the time variability of the gamma-ray emission and spectral properties of 2CG~135+1. Understanding the nature of 2CG 135+1 may be of crucial importance for the interpretation of a class of unidentified time variable gamma-ray sources in the Galaxy.