The Neutron star Interior Composition Explorer observed several rotation-powered millisecond pulsars (MSPs) 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 (≥4.7 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 and almost 3 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 MSPs exhibit a significant phase offset between their X-ray and radio pulses.
Context.
The gravitational strong equivalence principle (SEP) is a cornerstone of the general theory of relativity (GR). Hence, testing the validity of SEP is of great importance when confronting GR, ...or its alternatives, with experimental data. Pulsars that are orbited by white dwarf companions provide an excellent laboratory, where the extreme difference in binding energy between neutron stars and white dwarfs allows for precision tests of the SEP via the technique of radio pulsar timing.
Aims.
To date, the best limit on the validity of SEP under strong-field conditions was obtained with a unique pulsar in a triple stellar system, PSR J0337+1715. We report here on an improvement of this test using an independent data set acquired over a period of 6 years with the Nançay radio telescope. The improvements arise from a uniformly sampled data set, a theoretical analysis, and a treatment that fixes some short-comings in the previously published results, leading to better precision and reliability of the test.
Methods.
In contrast to the previously published test, we use a different long-term timing data set, developed a new timing model and an independent numerical integration of the motion of the system, and determined the masses and orbital parameters with a different methodology that treats the parameter Δ, describing a possible strong-field SEP violation, identically to all other parameters.
Results.
We obtain a violation parameter Δ = ( + 0.5 ± 1.8) × 10
−6
at 95% confidence level, which is compatible with and improves upon the previous study by 30%. This result is statistics-limited and avoids limitation by systematics as previously encountered. We find evidence for red noise in the pulsar spin frequency, which is responsible for up to 10% of the reported uncertainty. We use the improved limit on SEP violation to place constraints on a class of well-studied scalar-tensor theories, in particular we find
ω
BD
> 140 000 for the Brans-Dicke parameter. The conservative limits presented here fully take into account current uncertainties in the equation for state of neutron-star matter.
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 yr. 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 semimajor 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. However, we measure an average uncertainty of 80 per cent (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_{\rm p}=1.22_{-0.35}^{+0.5}\ {\rm M}_{{\odot }}$
,
$m_{\rm c} = 0.21_{-0.04}^{+0.06}\ {\rm M}_{{\odot } }$
and
$m_{\rm p}=1.25_{-0.4}^{+0.6}\ {\rm M}_{{\odot }}$
,
$m_{\rm c} = 0.23_{-0.05}^{+0.07}\ {\rm 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.
Abstract
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 yr 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 that indicate the pulsar is in either a low-eccentricity (0.33 ≲ e ≲ 0.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 ≳ 0.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 >7500 M⊙ located at the cluster centre. A gravitational model for the globular cluster, which includes such a central 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 ∼60 000 M⊙. Accounting for the associated contribution to the observed period derivative indicates that the γ-ray efficiency of the pulsar should be between 0.08 and 0.2. Our results suggest that other globular clusters may also contain central BHs and they may be revealed by the study of new pulsars found sufficiently close to their centres.
Abstract
We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior ...Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be
12.39
−
0.98
+
1.30
km and
2.072
−
0.066
+
0.067
M
⊙
respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is
log
10
(
T
K
)
=
5.99
−
0.06
+
0.05
for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.
ABSTRACT We report on the results of a recent blind search survey for gamma-ray pulsars in Fermi Large Area Telescope (LAT) data being carried out on the distributed volunteer computing system, ...Einstein@Home. The survey has searched for pulsations in 118 unidentified pulsar-like sources, requiring about 10,000 years of CPU core time. In total, this survey has resulted in the discovery of 17 new gamma-ray pulsars, of which 13 are newly reported in this work, and an accompanying paper. These pulsars are all young, isolated pulsars with characteristic ages between 12 kyr and 2 Myr, and spin-down powers between 1034 and 4 × 1036 erg s−1. Two of these are the slowest spinning gamma-ray pulsars yet known. One pulsar experienced a very large glitch during the Fermi mission. In this, the first of two associated papers, we describe the search scheme used in this survey, and estimate the sensitivity of our search to pulsations in unidentified Fermi-LAT sources. One such estimate results in an upper limit of 57% for the fraction of pulsed emission from the gamma-ray source associated with the Cas A supernova remnant, constraining the pulsed gamma-ray photon flux that can be produced by the neutron star at its center. We also present the results of precise timing analyses for each of the newly detected pulsars.
Symmetries play a fundamental 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 universality of free fall (UFF) for strongly self-gravitating bodies; (3) the post-Newtonian parameter ˆα3 in gravitational Lorentz invariance. Based on the pulsar timing result of the combined data set from the North American
Nanohertz Gravitational Observatory and the European Pulsar Timing Array, we find G˙ /G = (−0.1 ± 0.9) × 10−12 yr−1, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain an improved test for a UFF violation by a strongly self-gravitating mass falling in the gravitational field of our Galaxy, with a limit of |Delta| < 0.002 (95 per cent C.L.). Finally, we derive an improved limit on the self-acceleration of a gravitationally bound rotating body, to a preferred reference frame in the Universe, with −3 × 10−20 < ˆα3 < 4 × 10−20 (95 per cent C.L.). These results are based on direct UFF and ˆα3 tests using pulsar binaries, and they overcome various limitations of previous tests of this kind.
We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest European Pulsar Timing Array ...(EPTA) data set, which consists of ultraprecise timing data on 41-ms pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95 per cent upper limit on the sky-averaged strain amplitude lies in the range ... This limit varies by a factor of five, depending on the assumed source position and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit - obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is ... These limits, the most stringent to date at f < 10...nHz, exclude the presence of sub-centiparsec binaries with chirp mass ... out to a distance of about 25 Mpc, and with ... out to a distance of about 1Gpc (...). We show that state-of-the-art SMBHB population models predict <1 per cent probability of detecting a CGW with the current EPTA data set, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years. (ProQuest: ... denotes formulae/symbols omitted.)
We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar data set spanning 18 yr of observations from the 2015 European Pulsar Timing Array data ...release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar system ephemeris errors, obtaining a robust 95 per cent upper limit on the dimensionless strain amplitude A of the background of A < 3.0 × 10−15 at a reference frequency of 1 yr−1 and a spectral index of 13/3, corresponding to a background from inspiralling supermassive black hole binaries, constraining the GW energy density to Ωgw(f)h
2 < 1.1 × 10−9 at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of ∼5 × 10−9 Hz. Finally, we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95 per cent upper limits on the string tension, Gμ/c
2, characterizing a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu–Goto field theory cosmic string network, we set a limit Gμ/c
2 < 1.3 × 10−7, identical to that set by the Planck Collaboration, when combining Planck and high-ℓ cosmic microwave background data from other experiments. For a stochastic relic background, we set a limit of
$\Omega ^\mathrm{relic}_\mathrm{gw}(f)h^2<1.2 \times 10^{-9}$
, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.
Context.
Accurate polarimetric calibration of the radio pulse profiles from pulsars is crucial for studying their radiation properties at these wavelengths. Additionally, inaccurate calibration can ...distort recorded pulse profiles, introducing noise in time of arrival (TOA) data and thus degrading pulsar timing analyses. One method for determining the full polarimetric response of a given telescope is to conduct observations of bright polarized pulsars over wide ranges of parallactic angles, to sample different orientations of their polarization angle and in turn determine the cross-couplings between polarization feeds.
Aims.
The Nançay decimetric Radio Telescope (NRT) is a 94 m equivalent meridian telescope, capable of tracking a given pulsar for approximately 1 h around transit. The NRT therefore cannot sample wide ranges of parallactic angles when observing a given pulsar, so until late 2019 the polarimetric calibration of 1.4 GHz pulsar observations with the NRT was rudimentary. We therefore aimed to develop a method for improving the calibration of NRT observations, overcoming the above-mentioned limitation. Ultimately, our goal was to improve the quality of NRT pulsar timing, with better calibrated pulsar pulse profiles.
Methods.
In November 2019, we began conducting regular observations of the bright and highly linearly polarized pulsar PSR J0742−2822, in a special observing mode in which the feed horn rotates by ~180° over the course of the 1 h observation, mimicking wide parallactic angle variations and in principle enabling us to determine the polarimetric response of the NRT at 1.4 GHz. In addition, we assessed the quality of the NRT timing of a selection of millisecond pulsars (MSPs), namely, J1730−2304, J1744−1134, and J1857+0953, with conventional TOAs extracted from total intensity pulse profiles, and TOAs extracted with the Matrix Template Matching (MTM) technique, designed to compensate for putative polarimetric calibration errors.
Results.
From the analysis of the rotating horn observations of PSR J0742−2822 we could determine the cross-couplings between the polarization feeds and also constrain the Stokes parameters of the noise diode signal, which prior to this work was erroneously assumed to be ideal and was used as the only reference source for the calibration of pulsar observations. The improved polarimetric response of the NRT as determined from these observations was applied to observations of a selection of MSPs with published polarimetric properties. We find that the new polarimetric profiles and polarization position angles are consistent with previous findings, unlike NRT polarimetric results obtained with the previously used method of calibration. The analysis of the timing data shows that the new calibration method improves the quality of the timing, and the MTM method proves very effective at reducing noise from imperfect calibration. For pulsars with sufficient degrees of polarization, the MTM method appears to be the preferred method of extracting TOAs from NRT observations.