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.
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.
We report on the results of a 4 year timing campaign of PSR J2222−0137, a 2.44 day binary pulsar with a massive white dwarf (WD) companion, with the Nançay, Effelsberg, and Lovell radio telescopes. ...Using the Shapiro delay for this system, we find a pulsar mass mp = 1.76 0.06 M and a WD mass mc = 1.293 0.025 M . We also measure the rate of advance of periastron for this system, which is marginally consistent with the general relativity prediction for these masses. The short lifetime of the massive WD progenitor star led to a rapid X-ray binary phase with little (< 10−2 M ) mass accretion onto the neutron star; hence, the current pulsar mass is, within uncertainties, its birth mass, which is the largest measured to date. We discuss the discrepancy with previous mass measurements for this system; we conclude that the measurements presented here are likely to be more accurate. Finally, we highlight the usefulness of this system for testing alternative theories of gravity by tightly constraining the presence of dipolar radiation. This is of particular importance for certain aspects of strong-field gravity, like spontaneous scalarization, since the mass of PSR J2222−0137 puts that system into a poorly tested parameter range.
Context.
Interstellar scintillation analysis of pulsars allows us to probe the small-scale distribution and inhomogeneities of the ionized interstellar medium. From the frequency scale of ...scintillation, one can estimate the geometric time delays from multipath propagation, a source of (typically) unmodeled, correlated noise in pulsar timing. Interstellar scintillation analysis of well-timed pulsars is useful to quantify the effects of time delays and may lead to improved timing precision, enhancing the probability of detecting gravitational waves.
Aims.
Our priority is to present the data set and the basic measurements of scintillation parameters of pulsars, employing long-term scintillation observations carried out from 2011 January to 2020 August by the European Pulsar Timing Array radio telescopes in the 21-cm and 11-cm bands. Additionally, we aim to identify future possible lines of study using this long-term scintillation data set.
Methods.
The autocorrelation function of dynamic spectra has been used to estimate the scintillation bandwidth
v
d
and scintillation timescale
τ
d
.
Results.
We present the long-term time series of
v
d
and
τ
d
for 13 pulsars. Sanity checks and comparisons indicate that the scintillation parameters of our work and previously published works are mostly consistent. For two pulsars, PSRs J1857+0943 and J1939+2134, we were able to obtain measurements of the
v
d
at both bands, which allowed us to derive the time series of frequency scaling indices with a mean and a standard deviation of 2.82 ± 1.95 and 3.18 ± 0.60, respectively. We found some interesting features which will be studied in more detail in subsequent papers in this series: (i) in the time series of PSR J1939+2134, where
v
d
and
τ
d
sharply decrease associated with a sharp increase in the dispersion measure; (ii) PSR J0613-0200 and PSR J0636+5126 show a strong annual variation in the time series of the
τ
d
; and (iii) PSR J1939+2134 shows a weak anticorrelation between the scintillation timescale and the dispersion in Westerbork Synthesis Radio Telescope data.
ABSTRACT
Propagation effects are one of the main sources of noise in high-precision pulsar timing. For pulsars below an ecliptic latitude of 5°, the ionized plasma in the solar wind can introduce ...dispersive delays of order $100\, \mu \mathrm{s}$ around solar conjunction at an observing frequency of 300 MHz. A common approach to mitigate this assumes a spherical solar wind with a time-constant amplitude. However, this has been shown to be insufficient to describe the solar wind. We present a linear, Gaussian-process piecewise Bayesian approach to fit a spherical solar wind of time-variable amplitude, which has been implemented in the pulsar software run_enterprise. Through simulations, we find that the current EPTA+InPTA data combination is not sensitive to such variations; however, solar wind variations will become important in the near future with the addition of new InPTA data and data collected with the low-frequency LOFAR telescope. We also compare our results for different high-precision timing data sets (EPTA+InPTA, PPTA, and LOFAR) of 3 ms pulsars (J0030+0451, J1022+1001, J2145−0450), and find that the solar-wind amplitudes are generally consistent for any individual pulsar, but they can vary from pulsar to pulsar. Finally, we compare our results with those of an independent method on the same LOFAR data of the three millisecond pulsars. We find that differences between the results of the two methods can be mainly attributed to the modelling of dispersion variations in the interstellar medium, rather than the solar wind modelling.
Annual variations of interstellar scintillation can be modelled to constrain parameters of the ionized interstellar medium. If a pulsar is in a binary system, then investigating the orbital ...parameters is possible through analysis of the orbital variation of scintillation. In observations carried out from 2011 to 2020 by the European Pulsar Timing Array radio telescopes, PSRs J0613–0200 and J0636+5128 show strong annual variations in their scintillation velocity, while the former additionally exhibits an orbital fluctuation. Bayesian theory and Markov-chain-Monte-Carlo methods are used to interpret these periodic variations. We assume a thin and anisotropic scattering screen model, and discuss the mildly and extremely anisotropic scattering cases. PSR J0613–0200 is best described by mildly anisotropic scattering, while PSR J0636+5128 exhibits extremely anisotropic scattering. We measure the distance, velocity, and degree of anisotropy of the scattering screen for our two pulsars, finding that scattering screen distances from Earth for PSRs J0613–0200 and J0636+5128 are 316
−20
+28
pc and 262
−38
+96
pc, respectively. The positions of these scattering screens are coincident with the shell of the Local Bubble towards both pulsars. These associations add to the growing evidence of the Local Bubble shell as a dominant region of scattering along many sightlines.
Timing stability of three black widow pulsars Bak Nielsen, Ann-Sofie; Janssen, Gemma H; Shaifullah, Golam ...
Monthly notices of the Royal Astronomical Society,
2020, Letnik:
494, Številka:
2
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
We study the timing stability of three black widow pulsars (BWPs), both in terms of their long-term spin evolution and their shorter term orbital stability. The erratic timing behaviour and ...radio eclipses of the first two BWP systems discovered (PSRs B1957+20 and J2051−0827) were assumed to be representative for this class of pulsars. With several new black widow systems added to this population in the last decade, there are now several systems known that do not show these typical orbital variations or radio eclipses. We present timing solutions using 7–8 yr of observations from four of the European Pulsar Timing Array telescopes for PSRs J0023+0923, J2214+3000, and J2234+0944, and confirm that two of these systems do not show any significant orbital variability over our observing time span, both in terms of secular or orbital parameters. The third pulsar PSR J0023+0923 shows orbital variability and we discuss the implications for the timing solution. Our results from the long-term timing of these pulsars provide several new or improved parameters compared to earlier works. We discuss our results regarding the stability of these pulsars, and the stability of the class of BWPs in general, in the context of the binary parameters, and discuss the potential of the Roche lobe filling factor of the companion star being an indicator for stability of these systems.
ABSTRACT
We describe how to implement the spectral kurtosis method of interference removal (zapping) on a digitized signal of averaged power values. Spectral kurtosis is a hypothesis test, analogous ...to the t-test, with a null hypothesis that the amplitudes from which power is formed belong to a ‘good’ distribution – typically Gaussian with zero mean – where power values are zapped if the hypothesis is rejected at a specified confidence level. We derive signal-to-noise ratios (SNRs) as a function of amount of zapping for folded radio pulsar observations consisting of a sum of signals from multiple telescopes in independent radio-frequency interference environments, comparing four methods to compensate for lost data with coherent (tied-array) and incoherent summation. For coherently summed amplitudes, scaling amplitudes from non-zapped telescopes achieves a higher SNR than replacing zapped amplitudes with artificial noise. For incoherently summed power values, the highest SNR is given by scaling power from non-zapped telescopes to maintain a constant mean. We use spectral kurtosis to clean a tied-array radio pulsar observation by the Large European Array for Pulsars: the signal from one telescope is zapped with time and frequency resolutions of $6.25\, \mathrm{ms}$ and $0.16\, \mathrm{MHz}$, removing interference, along with 0.27 per cent of ‘good’ data, giving an uncertainty of $0.25\, \mathrm{\mu \mathrm{ s}}$ in pulse time of arrival (TOA) for PSR J1022+1001. We use a single-telescope observation to demonstrate recovery of the pulse profile shape, with 0.6 per cent of data zapped and a reduction from 1.22 to $0.70\, \mathrm{\mu \mathrm{ s}}$ in TOA uncertainty.