• ### Improving timing sensitivity in the microhertz frequency regime: limits from PSR J1713$+$0747 on gravitational waves produced by super-massive black-hole binaries(1804.10571)

April 27, 2018 astro-ph.HE
We search for continuous gravitational waves (CGWs) produced by individual super-massive black-hole binaries (SMBHBs) in circular orbits using high-cadence timing observations of PSR J1713$+$0747. We observe this millisecond pulsar using the telescopes in the European Pulsar Timing Array (EPTA) with an average cadence of approximately 1.6 days over the period between April 2011 and July 2015, including an approximately daily average between February 2013 and April 2014. The high-cadence observations are used to improve the pulsar timing sensitivity across the GW frequency range of $0.008-5$ $\mu$Hz. We use two algorithms in the analysis, including a spectral fitting method and a Bayesian approach. For an independent comparison, we also use a previously published Bayesian algorithm. We find that the Bayesian approaches provide optimal results and the timing observations of the pulsar place a 95 per cent upper limit on the sky-averaged strain amplitude of CGWs to be $\lesssim3.5\times10^{-13}$ at a reference frequency of 1 $\mu$Hz. We also find a 95 per cent upper limit on the sky-averaged strain amplitude of low-frequency CGWs to be $\lesssim1.4\times10^{-14}$ at a reference frequency of 20~nHz.
• ### French SKA White Book - The French Community towards the Square Kilometre Array(1712.06950)

March 28, 2018 astro-ph.IM
The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.
• ### MeerTime - the MeerKAT Key Science Program on Pulsar Timing(1803.07424)

March 18, 2018 astro-ph.IM, astro-ph.HE
The MeerKAT telescope represents an outstanding opportunity for radio pulsar timing science with its unique combination of a large collecting area and aperture efficiency (effective area $\sim$7500 m$^2$), system temperature ($T<20$K), high slew speeds (1-2 deg/s), large bandwidths (770 MHz at 20cm wavelengths), southern hemisphere location (latitude $\sim -30^\circ$) and ability to form up to four sub-arrays. The MeerTime project is a five-year program on the MeerKAT array by an international consortium that will regularly time over 1000 radio pulsars to perform tests of relativistic gravity, search for the gravitational wave signature induced by supermassive black hole binaries in the timing residuals of millisecond pulsars, explore the interiors of neutron stars through a pulsar glitch monitoring programme, explore the origin and evolution of binary pulsars, monitor the swarms of pulsars that inhabit globular clusters and monitor radio magnetars. In addition to these primary programmes, over 1000 pulsars will have their arrival times monitored and the data made immediately public. The MeerTime pulsar backend comprises two server-class machines each of which possess four Graphics Processing Units. Up to four pulsars can be coherently dedispersed simultaneously up to dispersion measures of over 1000 pc cm$^{-3}$. All data will be provided in psrfits format. The MeerTime backend will be capable of producing coherently dedispersed filterbank data for timing multiple pulsars in the cores of globular clusters that is useful for pulsar searches of tied array beams. All MeerTime data will ultimately be made available for public use, and any published results will include the arrival times and profiles used in the results.
• ### Tests of Gravitational Symmetries with Pulsar Binary J1713+0747(1802.09206)

Feb. 26, 2018 astro-ph.HE
Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|\Delta|< 0.002$ for the UFF test and $-3\times10^{-20} < \hat{\alpha}_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hat{\alpha}_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests.
• ### PSR J1618-3921: a recycled pulsar in an eccentric orbit(1801.07123)

Jan. 22, 2018 astro-ph.SR, astro-ph.HE
PSR~J1618$-$3921 is an $11.99$-ms pulsar in a $22.7$-d orbit around a likely low-mass He white dwarf companion, discovered in a survey of the intermediate Galactic latitudes at 1400 MHz conducted with the Parkes radio telescope in the late 1990s. Although PSR~J1618$-$3921 was discovered more than 15 years ago, only limited information has been published about this pulsar which has a surprisingly large orbital eccentricity ($e \simeq 0.027$), considering its high spin frequency and likely small companion mass. The focus of this work is a precise measurement of the spin, astrometric and orbital characteristics of PSR J1618$-$3921. This was done with timing observations made at the Nan\c{c}ay Radio Telescope, from 2009 to 2017. We analyzed the timing data recorded at the Nan\c{c}ay Radio Telescope over several years to characterize the properties of PSR~J1618$-$3921. A rotation ephemeris for this pulsar was obtained by analyzing the times of arrival of the radio pulses at the telescope. We confirm the unusual eccentricity of PSR~J1618$-$3921 and discuss several hypotheses regarding its formation, in the context of other discoveries of recycled pulsars in eccentric orbits.
• ### Evidence for an intermediate-mass black hole in the globular cluster NGC 6624(1705.01612)

PSR B1820$-$30A is located in the globular cluster NGC 6624 and is the closest known pulsar to the centre of any globular cluster. We present more than 25 years of high-precision timing observations of this millisecond pulsar and obtain four rotational frequency time derivative measurements. Modelling these higher-order derivatives as being due to orbital motion, we find solutions which indicate that the pulsar is in either a low-eccentricity ($0.33\lesssim e\lesssim0.4$) smaller orbit with a low mass companion (such as a main sequence star, white dwarf, neutron star, or stellar mass black hole) or a high-eccentricity ($e\gtrsim0.9$) larger orbit with a massive companion. The cluster mass properties and the observed properties of 4U 1820$-$30 and the other pulsars in the cluster argue against the low-eccentricity possibility. The high-eccentricity solution reveals that the pulsar is most likely orbiting around an intermediate-mass black hole (IMBH) of mass $> 7,500$~M$_\odot$ located at the cluster centre. A gravitational model for the globular cluster, which includes such a central black hole (BH), predicts an acceleration that is commensurate with that measured for the pulsar. It further predicts that the model-dependent minimum mass of the IMBH is $\sim60,000$~M$_\odot$. Accounting for the associated contribution to the observed period derivative indicates that the $\gamma$-ray efficiency of the pulsar should be between 0.08 and 0.2. Our results suggest that other globular clusters may also contain central black holes and they may be revealed by the study of new pulsars found sufficiently close to their centres. Note that we found an erratum in Section 5 and thus, the $\sim$60,000~M$_\odot$ mass mentioned above has to be replaced by the correct model-dependent mass limit of $\sim$20,000~M$_\odot$. See the erratum appended.
• ### A millisecond pulsar in an extremely wide binary system(1604.00129)

April 1, 2016 astro-ph.SR, astro-ph.HE
We report on 22 yrs of radio timing observations of the millisecond pulsar J1024$-$0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869$-$0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, $\mathrm{[M/H]}=-1.0$, $T_\mathrm{eff}=4050\pm50$ K) and that its position, proper motion and distance are consistent with those of PSR J1024$-$0719. We conclude that PSR J1024$-$0719 and 2MASS J10243869$-$0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide ($P_\mathrm{b}>200$ yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of $384\pm45$ km s$^{-1}$ with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024$-$0719 in light of its inclusion in pulsar timing arrays.
• ### High-precision timing of 42 millisecond pulsars with the European Pulsar Timing Array(1602.08511)

We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with TempoNest yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semi-major axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz-Kelker bias) than the M2 and M3 models by Schnitzeler (2012). However, we measure an average uncertainty of 80\% (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600$-$3053 and J1918$-$0642, implying pulsar and companion masses $m_p=1.22_{-0.35}^{+0.5} \text{M}_{\odot}$, $m_c = 0.21_{-0.04}^{+0.06} \text{M}_{\odot }$ and $m_p=1.25_{-0.4}^{+0.6} \text{M}_{\odot}$, $m_c = 0.23_{-0.05}^{+0.07} \text{M}_{\odot }$, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012$+$5307 and J1909$-$3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600$-$3053 and J1909$-$3744.
• ### From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release(1602.05570)

We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.
• ### The International Pulsar Timing Array: First Data Release(1602.03640)

Feb. 11, 2016 astro-ph.IM
The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.
• ### The gamma-ray millisecond pulsar deathline, revisited - New velocity and distance measurements(1601.05987)

Jan. 22, 2016 astro-ph.HE
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.
• ### The noise properties of 42 millisecond pulsars from the European Pulsar Timing Array and their impact on gravitational wave searches(1510.09194)

Jan. 20, 2016 astro-ph.CO, astro-ph.IM
The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispersion and scattering effects, analysis errors and instrumental instabilities. We present the results from a noise analysis for 42 millisecond pulsars (MSPs) observed with the European Pulsar Timing Array. For characterising the low-frequency, stochastic and achromatic noise component, or "timing noise", we employ two methods, based on Bayesian and frequentist statistics. For 25 MSPs, we achieve statistically significant measurements of their timing noise parameters and find that the two methods give consistent results. For the remaining 17 MSPs, we place upper limits on the timing noise amplitude at the 95% confidence level. We additionally place an upper limit on the contribution to the pulsar noise budget from errors in the reference terrestrial time standards (below 1%), and we find evidence for a noise component which is present only in the data of one of the four used telescopes. Finally, we estimate that the timing noise of individual pulsars reduces the sensitivity of this data set to an isotropic, stochastic GW background by a factor of >9.1 and by a factor of >2.3 for continuous GWs from resolvable, inspiralling supermassive black-hole binaries with circular orbits.
• ### Limits on anisotropy in the nanohertz stochastic gravitational-wave background(1506.08817)

June 29, 2015 gr-qc, astro-ph.CO, astro-ph.HE
The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in prep.)]. Our analysis of the GWB in the $\sim 2 - 90$ nHz band shows consistency with isotropy, with the strain amplitude in $l>0$ spherical harmonic multipoles $\lesssim 40\%$ of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations.
• ### A Cosmic Census of Radio Pulsars with the SKA(1501.00056)

The Square Kilometre Array (SKA) will make ground breaking discoveries in pulsar science. In this chapter we outline the SKA surveys for new pulsars, as well as how we will perform the necessary follow-up timing observations. The SKA's wide field-of-view, high sensitivity, multi-beaming and sub-arraying capabilities, coupled with advanced pulsar search backends, will result in the discovery of a large population of pulsars. These will enable the SKA's pulsar science goals (tests of General Relativity with pulsar binary systems, investigating black hole theorems with pulsar-black hole binaries, and direct detection of gravitational waves in a pulsar timing array). Using SKA1-MID and SKA1-LOW we will survey the Milky Way to unprecedented depth, increasing the number of known pulsars by more than an order of magnitude. SKA2 will potentially find all the Galactic radio-emitting pulsars in the SKA sky which are beamed in our direction. This will give a clear picture of the birth properties of pulsars and of the gravitational potential, magnetic field structure and interstellar matter content of the Galaxy. Targeted searches will enable detection of exotic systems, such as the ~1000 pulsars we infer to be closely orbiting Sgr A*, the supermassive black hole in the Galactic Centre. In addition, the SKA's sensitivity will be sufficient to detect pulsars in local group galaxies. To derive the spin characteristics of the discoveries we will perform live searches, and use sub-arraying and dynamic scheduling to time pulsars as soon as they are discovered, while simultaneously continuing survey observations. The large projected number of discoveries suggests that we will uncover currently unknown rare systems that can be exploited to push the boundaries of our understanding of astrophysics and provide tools for testing physics, as has been done by the pulsar community in the past.
• ### Broadband Pulsations from PSR B1821-24: Implications for Emission Models and the Pulsar Population of M28(1310.1852)

Oct. 7, 2013 astro-ph.HE
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
• ### Six millisecond pulsars detected by the Fermi Large Area Telescope and the radio/gamma-ray connection of millisecond pulsars(1212.4360)

Jan. 29, 2013 astro-ph.HE
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.
• ### An Asteroid Belt Interpretation for the Timing Variations of the Millisecond Pulsar B1937+21(1301.6429)

Jan. 28, 2013 astro-ph.SR, astro-ph.EP
Pulsar timing observations have revealed companions to neutron stars that include other neutron stars, white dwarfs, main-sequence stars, and planets. We demonstrate that the correlated and apparently stochastic residual times of arrival from the millisecond pulsar B1937+21 are consistent with the signature of an asteroid belt having a total mass less than approximately 0.05 Earth masses. Unlike the solar system's asteroid belt, the best fit pulsar asteroid belt extends over a wide range of radii, consistent with the absence of any shepherding companions. We suggest that any pulsar that has undergone accretion-driven spin-up and subsequently evaporated its companion may harbor orbiting asteroid mass objects. The resulting timing variations may fundamentally limit the timing precision of some of the other millisecond pulsars. Observational tests of the asteroid belt model include identifying periodicities from individual asteroids, which are difficult; testing for statistical stationarity that become possible when observations are conducted over a longer observing span; and searching for reflected radio emission.
• ### Radio Searches of Fermi LAT Sources and Blind Search Pulsars: The Fermi Pulsar Search Consortium(1205.3089)

May 14, 2012 astro-ph.HE
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.
• ### Discovery of the millisecond pulsar PSR J2043+1711 in a Fermi source with the Nancay Radio Telescope(1202.1128)

Feb. 6, 2012 astro-ph.GA, astro-ph.HE
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.
• ### Pulsed Gamma Rays from the Original Millisecond and Black Widow Pulsars: a case for Caustic Radio Emission?(1110.1271)

Nov. 7, 2011 astro-ph.GA, astro-ph.HE
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.
• ### The INPOP10a planetary ephemeris and its applications in fundamental physics(1108.5546)

Aug. 29, 2011 gr-qc, astro-ph.EP
Compared to the previous INPOP versions, the INPOP10a planetary and lunar ephemeris has several improvements. For the planets of our solar system, no big change was brought in the dynamics but improvements were implemented in the fitting process, the data sets used in the fit and in the selection of fitted parameters. We report here the main characteristics of the planetary part of INPOP10a like the fit of the product of the Solar mass with the gravitational constant (GM$_{\odot}$) instead of the astronomical unit. Determinations of PPN parameters as well as adjustments of the Sun J2 and of asteroid masses are also presented. New advances of nodes and perihelia of planets were also estimated and are given here. As for INPOP08, INPOP10a provides to the user, positions and velocities of the planets, the moon, the rotation angles of the Earth and the Moon as well as TT-TDB chebychev polynomials at http://www.imcce.fr/inpop
• ### Detection of the Pulsar Wind Nebula HESS J1825-137 with the Fermi Large Area Telescope(1106.0184)

June 1, 2011 astro-ph.HE
We announce the discovery of 1 - 100 GeV gamma-ray emission from the archetypal TeV pulsar wind nebula HESS J1825-137 using 20 months of survey data from the Fermi Large Area Telescope (LAT). The gamma-ray emission detected by the LAT is significantly spatially extended, with a best-fit rms extension of sigma = 0.56{\deg} $\pm$ 0.07{\deg} for an assumed Gaussian model. The 1 - 100 GeV LAT spectrum of this source is well described by a power-law with a spectral index of 1.38 $\pm$ 0.12 $\pm$ 0.16 and an integral flux above 1 GeV of (6.50 $\pm$ 0.21 $\pm$ 3.90) x 10^{-9} cm^{-2} s^{-1}. The first errors represent the statistical errors on the fit parameters, while the second ones are the systematic uncertainties. Detailed morphological and spectral analyses bring new constraints on the energetics and magnetic field of the pulsar wind nebula system. The spatial extent and hard spectrum of the GeV emission are consistent with the picture of an inverse Compton origin of the GeV-TeV emission in a cooling-limited nebula powered by the pulsar PSR J1826-1334.
• ### Placing limits on the stochastic gravitational-wave background using European Pulsar Timing Array data(1103.0576)

March 2, 2011 gr-qc, astro-ph.CO, astro-ph.IM
Direct detection of low-frequency gravitational waves ($10^{-9} - 10^{-8}$ Hz) is the main goal of pulsar timing array (PTA) projects. One of the main targets for the PTAs is to measure the stochastic background of gravitational waves (GWB) whose characteristic strain is expected to approximately follow a power-law of the form $h_c(f)=A (f/\hbox{yr}^{-1})^{\alpha}$, where $f$ is the gravitational-wave frequency. In this paper we use the current data from the European PTA to determine an upper limit on the GWB amplitude $A$ as a function of the unknown spectral slope $\alpha$ with a Bayesian algorithm, by modelling the GWB as a random Gaussian process. For the case $\alpha=-2/3$, which is expected if the GWB is produced by supermassive black-hole binaries, we obtain a 95% confidence upper limit on $A$ of $6\times 10^{-15}$, which is 1.8 times lower than the 95% confidence GWB limit obtained by the Parkes PTA in 2006. Our approach to the data analysis incorporates the multi-telescope nature of the European PTA and thus can serve as a useful template for future intercontinental PTA collaborations.
• ### Evidence for gravitational quadrupole moment variations in the companion of PSR J2051-0827(1102.5646)

Feb. 28, 2011 astro-ph.GA
We have conducted radio timing observations of the eclipsing millisecond binary pulsar J2051-0827 with the European Pulsar Timing Array network of telescopes and the Parkes radio telescope, spanning over 13 years. The increased data span allows significant measurements of the orbital eccentricity, e = (6.2 {\pm} 1.3) {\times} 10^{-5} and composite proper motion, {\mu}_t = 7.3 {\pm} 0.4 mas/yr. Our timing observations have revealed secular variations of the projected semi-major axis of the pulsar orbit which are much more extreme than those previously published; and of the orbital period of the system. Investigation of the physical mechanisms producing such variations confirm that the variations of the semi-major axis are most probably caused by classical spin-orbit coupling in the binary system, while the variations in orbital period are most likely caused by tidal dissipation leading to changes in the gravitational quadrupole moment of the companion.
• ### Discovery of two millisecond pulsars in Fermi sources with the Nancay Radio Telescope(1102.4192)

Feb. 27, 2011 astro-ph.GA, astro-ph.HE
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.