Showing 1–50 of 136 results for author: Bogdanov, S

A NICER View of the Nearest and Brightest Millisecond Pulsar: PSR J0437$\unicode{x2013}$4715

Authors: Devarshi Choudhury, Tuomo Salmi, Serena Vinciguerra, Thomas E. Riley, Yves Kini, Anna L. Watts, Bas Dorsman, Slavko Bogdanov, Sebastien Guillot, Paul S. Ray, Daniel J. Reardon, Ronald A. Remillard, Anna V. Bilous, Daniela Huppenkothen, James M. Lattimer, Nathan Rutherford, Zaven Arzoumanian, Keith C. Gendreau, Sharon M. Morsink, Wynn C. G. Ho

Abstract: We report Bayesian inference of the mass, radius and hot X-ray emitting region properties - using data from the Neutron Star Interior Composition ExploreR (NICER) - for the brightest rotation-powered millisecond X-ray pulsar PSR J0437$\unicode{x2013}$4715. Our modeling is conditional on informative tight priors on mass, distance and binary inclination obtained from radio pulsar timing using the Pa… ▽ More We report Bayesian inference of the mass, radius and hot X-ray emitting region properties - using data from the Neutron Star Interior Composition ExploreR (NICER) - for the brightest rotation-powered millisecond X-ray pulsar PSR J0437$\unicode{x2013}$4715. Our modeling is conditional on informative tight priors on mass, distance and binary inclination obtained from radio pulsar timing using the Parkes Pulsar Timing Array (PPTA) (Reardon et al. 2024), and we use NICER background models to constrain the non-source background, cross-checking with data from XMM-Newton. We assume two distinct hot emitting regions, and various parameterized hot region geometries that are defined in terms of overlapping circles; while simplified, these capture many of the possibilities suggested by detailed modeling of return current heating. For the preferred model identified by our analysis we infer a mass of $M = 1.418 \pm 0.037$ M$_\odot$ (largely informed by the PPTA mass prior) and an equatorial radius of $R = 11.36^{+0.95}_{-0.63}$ km, each reported as the posterior credible interval bounded by the 16% and 84% quantiles. This radius favors softer dense matter equations of state and is highly consistent with constraints derived from gravitational wave measurements of neutron star binary mergers. The hot regions are inferred to be non-antipodal, and hence inconsistent with a pure centered dipole magnetic field. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 27 pages, 16 figures, + appendices. Accepted in ApJL. Data files and reproduction package available in Zenodo

The Magnificent Seven X-ray Isolated Neutron Stars Revisited. I. Improved Timing Solutions and Pulse Profile Analysis

Authors: Slavko Bogdanov, Wynn C. G. Ho

Abstract: We present the first systematic X-ray pulse timing analysis of the six members of the so-called "Magnificent Seven" nearby thermally-emitting isolated neutron stars (XINS) with detected pulsations. Using the extensive collection of archival XMM-Newton, Chandra, and NICER observations spanning over two decades, we obtain the first firm measurement of the spin-down rate for RX J2143.0+0654, while fo… ▽ More We present the first systematic X-ray pulse timing analysis of the six members of the so-called "Magnificent Seven" nearby thermally-emitting isolated neutron stars (XINS) with detected pulsations. Using the extensive collection of archival XMM-Newton, Chandra, and NICER observations spanning over two decades, we obtain the first firm measurement of the spin-down rate for RX J2143.0+0654, while for the rest we improve upon previously published spin ephemerides and extend them by up to an additional decade. Five of the XINS follow steady spin-down with no indication of major anomalies in their long-term timing behavior; the notable exception is RX J0720.4-3125, for which, in addition to confirming the previously identified glitch, we detect a second spin derivative. The high quality folded X-ray pulse profiles produced with the updated timing solutions exhibit diverse and complex morphologies, as well as striking energy dependence. These peculiarities cannot be readily explained by blackbody-like isotropic emission and simple hot spot configurations, hinting at the presence of complex multi-temperature surface heat distributions and highly anisotropic radiation patterns, such as may arise from a strongly magnetized atmospheric layer. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: 28 pages, 8 figures, 8 tables; accepted for publication in the Astrophysical Journal. The pulse timing solutions for the six neutron stars are included in the TeX source of this article in the form of Tempo-style parameter files

A More Precise Measurement of the Radius of PSR J0740+6620 Using Updated NICER Data

Authors: Alexander J. Dittmann, M. Coleman Miller, Frederick K. Lamb, Isiah Holt, Cecilia Chirenti, Michael T. Wolff, Slavko Bogdanov, Sebastien Guillot, Wynn C. G. Ho, Sharon M. Morsink, Zaven Arzoumanian, Keith C. Gendreau

Abstract: PSR J0740+6620 is the neutron star with the highest precisely determined mass, inferred from radio observations to be $2.08\pm0.07\,\rm M_\odot$. Measurements of its radius therefore hold promise to constrain the properties of the cold, catalyzed, high-density matter in neutron star cores. Previously, Miller et al. (2021) and Riley et al. (2021) reported measurements of the radius of PSR J0740+662… ▽ More PSR J0740+6620 is the neutron star with the highest precisely determined mass, inferred from radio observations to be $2.08\pm0.07\,\rm M_\odot$. Measurements of its radius therefore hold promise to constrain the properties of the cold, catalyzed, high-density matter in neutron star cores. Previously, Miller et al. (2021) and Riley et al. (2021) reported measurements of the radius of PSR J0740+6620 based on Neutron Star Interior Composition Explorer (NICER) observations accumulated through 17 April 2020, and an exploratory analysis utilizing NICER background estimates and a data set accumulated through 28 December 2021 was presented in Salmi et al. (2022). Here we report an updated radius measurement, derived by fitting models of X-ray emission from the neutron star surface to NICER data accumulated through 21 April 2022, totaling $\sim1.1$ Ms additional exposure compared to the data set analyzed in Miller et al. (2021) and Riley et al. (2021), and to data from X-ray Multi-Mirror (XMM-Newton) observations. We find that the equatorial circumferential radius of PSR J0740+6620 is $12.92_{-1.13}^{+2.09}$ km (68% credibility), a fractional uncertainty $\sim83\%$ the width of that reported in Miller et al. (2021), in line with statistical expectations given the additional data. If we were to require the radius to be less than 16 km, as was done in Salmi et al. (2024), then our 68% credible region would become $R=12.76^{+1.49}_{-1.02}$ km, which is close to the headline result of Salmi et al. (2024). Our updated measurements, along with other laboratory and astrophysical constraints, imply a slightly softer equation of state than that inferred from our previous measurements. △ Less

Submitted 30 June, 2024; v1 submitted 20 June, 2024; originally announced June 2024.

Comments: 18 pages, 8 figures, +appendices. Accepted in ApJ

The Radius of the High Mass Pulsar PSR J0740+6620 With 3.6 Years of NICER Data

Authors: Tuomo Salmi, Devarshi Choudhury, Yves Kini, Thomas E. Riley, Serena Vinciguerra, Anna L. Watts, Michael T. Wolff, Zaven Arzoumanian, Slavko Bogdanov, Deepto Chakrabarty, Keith Gendreau, Sebastien Guillot, Wynn C. G. Ho, Daniela Huppenkothen, Renee M. Ludlam, Sharon M. Morsink, Paul S. Ray

Abstract: We report an updated analysis of the radius, mass, and heated surface regions of the massive pulsar PSR J0740+6620 using NICER data from 2018 September 21 to 2022 April 21, a substantial increase in data set size compared to previous analyses. Using a tight mass prior from radio timing measurements and jointly modeling the new NICER data with XMM-Newton data, the inferred equatorial radius and gra… ▽ More We report an updated analysis of the radius, mass, and heated surface regions of the massive pulsar PSR J0740+6620 using NICER data from 2018 September 21 to 2022 April 21, a substantial increase in data set size compared to previous analyses. Using a tight mass prior from radio timing measurements and jointly modeling the new NICER data with XMM-Newton data, the inferred equatorial radius and gravitational mass are $12.49_{-0.88}^{+1.28}$ km and $2.073_{-0.069}^{+0.069}$ $M_\odot$ respectively, each reported as the posterior credible interval bounded by the $16\,\%$ and $84\,\%$ quantiles, with an estimated systematic error $\lesssim 0.1$ km. This result was obtained using the best computationally feasible sampler settings providing a strong radius lower limit but a slightly more uncertain radius upper limit. The inferred radius interval is also close to the $R=12.76_{-1.02}^{+1.49}$ km obtained by Dittmann et al. 2024, when they require the radius to be less than $16$ km as we do. The results continue to disfavor very soft equations of state for dense matter, with $R<11.15$ km for this high mass pulsar excluded at the $95\,\%$ probability. The results do not depend significantly on the assumed cross-calibration uncertainty between NICER and XMM-Newton. Using simulated data that resemble the actual observations, we also show that our pipeline is capable of recovering parameters for the inferred models reported in this paper. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 17 pages, 9 figures (2 of which are figure sets), 2 tables, accepted for publication in ApJ

Exploring Waveform Variations among Neutron Star Ray-tracing Codes for Complex Emission Geometries

Authors: Devarshi Choudhury, Anna L. Watts, Alexander J. Dittmann, M. Coleman Miller, Sharon M. Morsink, Tuomo Salmi, Serena Vinciguerra, Slavko Bogdanov, Sebastien Guillot, Michael T. Wolff, Zaven Arzoumanian

Abstract: Pulse Profile Modeling (PPM), the technique used to infer mass, radius and geometric parameters for rotation-powered millisecond pulsars using data from the Neutron Star Interior Composition Explorer (NICER), relies on relativistic ray-tracing of thermal X-ray photons from hot spots on the neutron star surface to the observer. To verify our ray-tracing codes we have in the past conducted cross-tes… ▽ More Pulse Profile Modeling (PPM), the technique used to infer mass, radius and geometric parameters for rotation-powered millisecond pulsars using data from the Neutron Star Interior Composition Explorer (NICER), relies on relativistic ray-tracing of thermal X-ray photons from hot spots on the neutron star surface to the observer. To verify our ray-tracing codes we have in the past conducted cross-tests for simple hot spot geometries, focusing primarily on the implementation of the space-time model. In this paper, we present verification for test problems that explore the more complex hot spot geometries that are now being employed in the NICER PPM analyses. We conclude that the accuracy of our computed waveforms is in general sufficiently high for analyses of current NICER data sets. We have however identified some extreme configurations where extra care may be needed. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: Submitted to ApJ

Pulse Profile Modelling of Thermonuclear Burst Oscillations III : Constraining the properties of XTE J1814-338

Authors: Yves Kini, Tuomo Salmi, Serena Vinciguerra, Anna L. Watts, Anna Bilous, Duncan K. Galloway, Emma van der Wateren, Guru Partap Khalsa, Slavko Bogdanov, Johannes Buchner, Valery Suleimanov

Abstract: Pulse profile modelling (PPM) is a comprehensive relativistic ray-tracing technique employed to determine the properties of neutron stars. In this study, we apply this technique to the Type I X-ray burster and accretion-powered millisecond pulsar XTE J1814-338, extracting its fundamental properties using PPM of its thermonuclear burst oscillations. Using data from its 2003 outburst, and a single u… ▽ More Pulse profile modelling (PPM) is a comprehensive relativistic ray-tracing technique employed to determine the properties of neutron stars. In this study, we apply this technique to the Type I X-ray burster and accretion-powered millisecond pulsar XTE J1814-338, extracting its fundamental properties using PPM of its thermonuclear burst oscillations. Using data from its 2003 outburst, and a single uniform temperature hot spot model, we infer XTE J1814-338 to be located at a distance of $7.2^{+0.3}_{-0.4}$ kpc, with a mass of $1.21^{+0.05}_{-0.05}$ M$_\odot$ and an equatorial radius of $7.0^{+0.4}_{-0.4}$ km. Our results also offer insight into the time evolution of the hot spot but point to some potential shortcomings of the single uniform temperature hot spot model. We explore the implications of this result, including what we can learn about thermonuclear burst oscillation mechanisms and the importance of modelling the accretion contribution to the emission during the burst. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: submitted to MNRAS. The Zenodo link will go public after peer review. Comments are welcome

Pulse Profile Modelling of Thermonuclear Burst Oscillations II: Handling variability

Authors: Yves Kini, Tuomo Salmi, Serena Vinciguerra, Anna L. Watts, Devarshi Choudhury, Slavko Bogdanov, Johannes Buchner, Zach Meisel, Valery Suleimanov

Abstract: Pulse profile modelling is a relativistic ray-tracing technique that can be used to infer masses, radii and geometric parameters of neutron stars. In a previous study, we looked at the performance of this technique when applied to thermonuclear burst oscillations from accreting neutron stars. That study showed that ignoring the variability associated with burst oscillation sources resulted in sign… ▽ More Pulse profile modelling is a relativistic ray-tracing technique that can be used to infer masses, radii and geometric parameters of neutron stars. In a previous study, we looked at the performance of this technique when applied to thermonuclear burst oscillations from accreting neutron stars. That study showed that ignoring the variability associated with burst oscillation sources resulted in significant biases in the inferred mass and radius, particularly for the high count rates that are nominally required to obtain meaningful constraints. In this follow-on study, we show that the bias can be mitigated by slicing the bursts into shorter segments where variability can be neglected, and jointly fitting the segments. Using this approach, the systematic uncertainties on the mass and radius are brought within the range of the statistical uncertainty. With about 10$^6$ source counts, this yields uncertainties of approximately 10% for both the mass and radius. However, this modelling strategy requires substantial computational resources. We also confirm that the posterior distributions of the mass and radius obtained from multiple bursts of the same source can be merged to produce outcomes comparable to that of a single burst with an equivalent total number of counts. △ Less

Submitted 22 November, 2023; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS. The Zenodo link is now public

An updated mass-radius analysis of the 2017-2018 NICER data set of PSR J0030+0451

Authors: Serena Vinciguerra, Tuomo Salmi, Anna L. Watts, Devarshi Choudhury, Thomas E. Riley, Paul S. Ray, Slavko Bogdanov, Yves Kini, Sebastien Guillot, Deepto Chakrabarty, Wynn C. G. Ho, Daniela Huppenkothen, Sharon M. Morsink, Zorawar Wadiasingh

Abstract: In 2019 the NICER collaboration published the first mass and radius inferred for PSR J0030+0451, thanks to NICER observations, and consequent constraints on the equation of state characterising dense matter. Two independent analyses found a mass of $\sim 1.3-1.4\,\mathrm{M_\odot}$ and a radius of $\sim 13\,$km. They also both found that the hot spots were all located on the same hemisphere, opposi… ▽ More In 2019 the NICER collaboration published the first mass and radius inferred for PSR J0030+0451, thanks to NICER observations, and consequent constraints on the equation of state characterising dense matter. Two independent analyses found a mass of $\sim 1.3-1.4\,\mathrm{M_\odot}$ and a radius of $\sim 13\,$km. They also both found that the hot spots were all located on the same hemisphere, opposite to the observer, and that at least one of them had a significantly elongated shape. Here we reanalyse, in greater detail, the same NICER data set, incorporating the effects of an updated NICER response matrix and using an upgraded analysis framework. We expand the adopted models and jointly analyse also XMM-Newton data, which enables us to better constrain the fraction of observed counts coming from PSR J0030+0451. Adopting the same models used in previous publications, we find consistent results, although with more stringent inference requirements. We also find a multi-modal structure in the posterior surface. This becomes crucial when XMM-Newton data is accounted for. Including the corresponding constraints disfavors the main solutions found previously, in favor of the new and more complex models. These have inferred masses and radii of $\sim [1.4 \mathrm{M_\odot}, 11.5$ km] and $\sim [1.7 \mathrm{M_\odot}, 14.5$ km], depending on the assumed model. They display configurations that do not require the two hot spots generating the observed X-rays to be on the same hemisphere, nor to show very elongated features, and point instead to the presence of temperature gradients and the need to account for them. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: 35 pages, 11 figures, 7 tables

Journal ref: ApJ, 961, 62 (2024)

Atmospheric Effects on Neutron Star Parameter Constraints with NICER

Authors: Tuomo Salmi, Serena Vinciguerra, Devarshi Choudhury, Anna L. Watts, Wynn C. G. Ho, Sebastien Guillot, Yves Kini, Bas Dorsman, Sharon M. Morsink, Slavko Bogdanov

Abstract: We present an analysis of the effects of uncertainties in the atmosphere models on the radius, mass, and other neutron star parameter constraints for the NICER observations of rotation-powered millisecond pulsars. To date, NICER has applied the X-ray pulse profile modeling technique to two millisecond-period pulsars: PSR J0030+0451 and the high-mass pulsar PSR J0740+6620. These studies have common… ▽ More We present an analysis of the effects of uncertainties in the atmosphere models on the radius, mass, and other neutron star parameter constraints for the NICER observations of rotation-powered millisecond pulsars. To date, NICER has applied the X-ray pulse profile modeling technique to two millisecond-period pulsars: PSR J0030+0451 and the high-mass pulsar PSR J0740+6620. These studies have commonly assumed a deep-heated, fully ionized hydrogen atmosphere model, although they have explored the effects of partial-ionization and helium composition in some cases. Here, we extend that exploration and also include new models with partially ionized carbon composition, externally heated hydrogen, and an empirical atmospheric beaming parameterization to explore deviations in the expected anisotropy of the emitted radiation. None of the studied atmosphere cases have any significant influence on the inferred radius of PSR J0740+6620, possibly due to its X-ray faintness, tighter external constraints, and/or viewing geometry. In the case of PSR J0030+0451, both the composition and ionization state could significantly alter the inferred radius. However, based on the evidence (prior predictive probability of the data), partially ionized hydrogen and carbon atmospheres are disfavored. The difference in the evidence for ionized hydrogen and helium atmospheres is too small to be decisive for most cases, but the inferred radius for helium models trends to larger sizes around or above 14-15 km. External heating or deviations in the beaming that are less than $5\,\%$ at emission angles smaller than 60 degrees, on the other hand, have no significant effect on the inferred radius. △ Less

Submitted 23 October, 2023; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: 26 pages, 12 figures (2 of which are figure sets), 3 tables. Published in ApJ

Journal ref: ApJ 956 138 (2023)

Green Bank Telescope Discovery of the Redback Binary Millisecond Pulsar PSR J0212+5321

Authors: Karen I. Perez, Slavko Bogdanov, Jules P. Halpern, Vishal Gajjar

Abstract: We report the discovery of a 2.11 ms binary millisecond pulsar during a targeted search of the redback optical candidate coincident with the $γ$-ray source 3FGL J0212.5+5320 using the Robert C. Byrd Green Bank Telescope (GBT) with the Breakthrough Listen backend at L-band. Over a seven month period, five pointings were made near inferior conjunction of the pulsar in its 20.9 hr orbit, resulting in… ▽ More We report the discovery of a 2.11 ms binary millisecond pulsar during a targeted search of the redback optical candidate coincident with the $γ$-ray source 3FGL J0212.5+5320 using the Robert C. Byrd Green Bank Telescope (GBT) with the Breakthrough Listen backend at L-band. Over a seven month period, five pointings were made near inferior conjunction of the pulsar in its 20.9 hr orbit, resulting in two detections, lasting 12 and 42 minutes. The pulsar dispersion measure (DM) of 25.7 pc cm$^{-3}$ corresponds to a distance of 1.15 kpc in the NE2001 Galactic electron density model, consistent with the Gaia parallax distance of $1.16\pm0.03$ kpc for the companion star. We suspect the pulsar experiences wide-orbit eclipses, similar to other redbacks, as well as scintillation and DM delays caused by its interaction with its companion and surroundings. Although the pulsar was only detected over $\approx3.7\%$ of the orbit, its measured acceleration is consistent with published binary parameters from optical radial velocity spectroscopy and light-curve modeling of the companion star, and it provides a more precise mass ratio and a projected semi-major axis for the pulsar orbit. We also obtained a refined optical photometric orbit ephemeris, and observed variability of the tidally distorted companion over 7 years. A hard X-ray light curve from NuSTAR shows expected orbit-modulated emission from the intrabinary shock. The pulsar parameters and photometric ephemeris greatly restrict the parameter space required to search for a coherent timing solution including pulsar spin-down rate, either using Fermi $γ$-rays, or further radio pulse detections. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: 20 pages, 8 figures, accepted to ApJ

Pulse Profile Modeling of Thermonuclear Burst Oscillations I: The Effect of Neglecting Variability

Authors: Yves Kini, Tuomo Salmi, Anna L. Watts, Serena Vinciguerra, Devarshi Choudhury, Siem Fenne, Slavko Bogdanov, Zach Meisel, Valery Suleimanov

Abstract: We study the effects of the time-variable properties of thermonuclear X-ray bursts on modeling their millisecond-period burst oscillations. We apply the pulse profile modeling technique that is being used in the analysis of rotation-powered millisecond pulsars by the Neutron Star Interior Composition Explorer (NICER) to infer masses, radii, and geometric parameters of neutron stars. By simulating… ▽ More We study the effects of the time-variable properties of thermonuclear X-ray bursts on modeling their millisecond-period burst oscillations. We apply the pulse profile modeling technique that is being used in the analysis of rotation-powered millisecond pulsars by the Neutron Star Interior Composition Explorer (NICER) to infer masses, radii, and geometric parameters of neutron stars. By simulating and analyzing a large set of models, we show that overlooking burst time-scale variability in temperatures and sizes of the hot emitting regions can result in substantial bias in the inferred mass and radius. To adequately infer neutron star properties, it is essential to develop a model for the time variable properties or invest a substantial amount of computational time in segmenting the data into non-varying pieces. We discuss prospects for constraints from proposed future X-ray telescopes. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: Accepted for publication in MNRAS

The neutron star population in M28: a joint Chandra/GBT look at pulsar paradise

Authors: Eda Vurgun, Manuel Linares, Scott Ransom, Alessandro Papitto, Slavko Bogdanov, Enrico Bozzo, Nanda Rea, Domingo Garcia-Senz, Paulo Freire, Ingrid Stairs

Abstract: We present the results of a deep study of the neutron star (NS) population in the globular cluster M28 (NGC 6626), using the full 330-ks 2002-2015 ACIS dataset from the Chandra X-ray Observatory and coordinated radio observations taken with the Green Bank Telescope (GBT) in 2015. We investigate the X-ray luminosity (Lx), spectrum, and orbital modulation of the 7 known compact binary millisecond pu… ▽ More We present the results of a deep study of the neutron star (NS) population in the globular cluster M28 (NGC 6626), using the full 330-ks 2002-2015 ACIS dataset from the Chandra X-ray Observatory and coordinated radio observations taken with the Green Bank Telescope (GBT) in 2015. We investigate the X-ray luminosity (Lx), spectrum, and orbital modulation of the 7 known compact binary millisecond pulsars (MSPs) in the cluster. We report two simultaneous detections of the redback PSR J1824-2452I (M28I) and its X-ray counterpart. We discover a double-peaked X-ray orbital flux modulation in M28I during its pulsar state, centered around pulsar inferior conjunction. We analyze the spectrum of the quiescent neutron star low-mass X-ray binary to constrain its mass and radius. Using both hydrogen and helium NS atmosphere models, we find a NS radius of R = 9.5-11.5 km and R = 13.5 - 16.7 km, respectively, for a neutron star mass of 1.4 Msun. We also search for long-term variability in the 46 brightest X-ray sources and report the discovery of six new variable low luminosity X-ray sources in M28. △ Less

Submitted 21 December, 2022; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 26 pages, Accepted for publication in The Astrophysical Journal

The Radius of PSR J0740+6620 from NICER with NICER Background Estimates

Authors: Tuomo Salmi, Serena Vinciguerra, Devarshi Choudhury, Thomas E. Riley, Anna L. Watts, Ronald A. Remillard, Paul S. Ray, Slavko Bogdanov, Sebastien Guillot, Zaven Arzoumanian, Cecilia Chirenti, Alexander J. Dittmann, Keith C. Gendreau, Wynn C. G. Ho, M. Coleman Miller, Sharon M. Morsink, Zorawar Wadiasingh, Michael T. Wolff

Abstract: We report a revised analysis for the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, studied previously with joint fits to NICER and XMM-Newton data by Riley et al. (2021) and Miller et al. (2021). We perform a similar Bayesian estimation for the pulse-profile model parameters, except that instead of fitting simultaneously the XMM-Newton data, we use the bes… ▽ More We report a revised analysis for the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, studied previously with joint fits to NICER and XMM-Newton data by Riley et al. (2021) and Miller et al. (2021). We perform a similar Bayesian estimation for the pulse-profile model parameters, except that instead of fitting simultaneously the XMM-Newton data, we use the best available NICER background estimates to constrain the number of photons detected from the source. This approach eliminates any potential issues in the cross-calibration between these two instruments, providing thus an independent check of the robustness of the analysis. The obtained neutron star parameter constraints are compatible with the already published results, with a slight dependence on how conservative the imposed background limits are. A tighter lower limit causes the inferred radius to increase, and a tighter upper limit causes it to decrease. We also extend the study of the inferred emission geometry to examine the degree of deviation from antipodality of the hot regions. We show that there is a significant offset to an antipodal spot configuration, mainly due to the non-half-cycle azimuthal separation of the two emitting spots. The offset angle from the antipode is inferred to be above 25 degrees with 84% probability. This seems to exclude a centered-dipolar magnetic field in PSR J0740+6620. △ Less

Submitted 13 January, 2023; v1 submitted 26 September, 2022; originally announced September 2022.

Comments: 29 pages, 18 figures (10 of which are figure sets), 1 animation, 10 tables. Published in ApJ

Journal ref: ApJ 941 150 (2022)

Snowmass 2021 Cosmic Frontier White Paper: The Dense Matter Equation of State and QCD Phase Transitions

Authors: Slavko Bogdanov, Emmanuel Fonseca, Rahul Kashyap, Aleksi Kurkela, James M. Lattimer, Jocelyn S. Read, Bangalore S. Sathyaprakash, H. Thankful Cromartie, Tim Dietrich, Arnab Dhani, Timothy Dolch, Tyler Gorda, Sebastien Guillot, Wynn C. G. Ho, Rachael Huxford, Frederick K. Lamb, Philippe Landry, Bradley W. Meyers, M. Coleman Miller, Joonas Nättilä, Risto Paatelainen, Chanda Prescod-Weinstein, Saga Säppi, Ingrid H. Stairs, Nikolaos Stergioulas , et al. (4 additional authors not shown)

Abstract: Our limited understanding of the physical properties of matter at ultra-high density, high proton/neutron number asymmetry, and low temperature is presently one of the major outstanding problems in physics. As matter in this extreme state is known to only exist stably in the cores of neutron stars (NSs), complementary measurements from electromagnetic and gravitational wave astrophysical observati… ▽ More Our limited understanding of the physical properties of matter at ultra-high density, high proton/neutron number asymmetry, and low temperature is presently one of the major outstanding problems in physics. As matter in this extreme state is known to only exist stably in the cores of neutron stars (NSs), complementary measurements from electromagnetic and gravitational wave astrophysical observations of NSs, combined with terrestrial laboratory constraints and further theoretical investigations, hold the promise to provide important insight into the properties of matter in a region of the quantum chromodynamics phase space that is otherwise inaccessible. This multidisciplinary endeavor imposes the following requirements for facilities and resources in the upcoming decade and beyond: * A next generation of gravitational wave detectors to uncover more double NS and neutron star-black hole mergers; * Sensitive radio telescopes to find the most massive and fastest spinning NSs; * Large-area, high-time-resolution and/or high angular resolution X-ray telescopes to constrain the NS mass-radius relation; * Suitable laboratory facilities for nuclear physics experiments to constrain the dense matter equation of state; * Funding resources for theoretical studies of matter in this regime; * The availability of modern large-scale high performance computing infrastructure. The same facilities and resources would also enable significant advances in other high-profile fields of inquiry in modern physics such as the nature of dark matter, alternative theories of gravity, nucleon superfluidity and superconductivity, as well as an array of astrophysics, including but not limited to stellar evolution, nucleosynthesis, and primordial black holes. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021) under Cosmic Frontier (CF07: Cosmic probes of fundamental physics); 30 pages, 8 figures

Luminous Optical and X-ray Flaring of the Putative Redback Millisecond Pulsar 1FGL J0523.5$-$2529

Authors: Jules P. Halpern, Karen I. Perez, Slavko Bogdanov

Abstract: Several redback and black widow millisecond pulsar binaries have episodes of flaring in X-rays and optical. We initially detected such behavior from the Fermi selected redback candidate 1FGL J0523.5$-$2529 during optical time-series monitoring. Triggered observations with the Neil Gehrels Swift Observatory over the next $\approx100$ days showed episodic flaring in X-rays with luminosity up to… ▽ More Several redback and black widow millisecond pulsar binaries have episodes of flaring in X-rays and optical. We initially detected such behavior from the Fermi selected redback candidate 1FGL J0523.5$-$2529 during optical time-series monitoring. Triggered observations with the Neil Gehrels Swift Observatory over the next $\approx100$ days showed episodic flaring in X-rays with luminosity up to $8\times10^{33}$ erg s$^{-1}$ ($\sim100$ times the minimum), and a comparable luminosity in the optical/UV, with similar power-law spectra of $f_ν\proptoν^{-0.7}$. These are the most luminous flares seen in any non-accreting "spider" pulsar system, which may be related to the large size of the companion through the fraction of the pulsar wind that it or its ablated wind intercepts. Simultaneously with an optical flare, we see Balmer-line and He I emission, not previously known in this object, which is evidence of a stellar wind that may also inhibit detection of radio pulsations. The quiescent optical light curves, while dominated by ellipsoidal modulation, show evidence of variable non-uniform temperature that could be due either to large starspots or asymmetric heating of the companion by the pulsar. This may explain a previous measurement of unusual non-zero orbital eccentricity as, alternatively, distortion of the radial-velocity curve by the surface temperature distribution of the large companion. △ Less

Submitted 17 July, 2022; originally announced July 2022.

Comments: 28 pages, 10 figures, to appear in The Astrophysical Journal

Timing six energetic rotation-powered X-ray pulsars, including the fast-spinning young PSR J0058-7218 and Big Glitcher PSR J0537-6910

Authors: Wynn C. G. Ho, Lucien Kuiper, Cristobal M. Espinoza, Sebastien Guillot, Paul S. Ray, D. A. Smith, Slavko Bogdanov, Danai Antonopoulou, Zaven Arzoumanian, Michal Bejger, Teruaki Enoto, Paolo Esposito, Alice K. Harding, Brynmor Haskell, Natalia Lewandowska, Chandreyee Maitra, Georgios Vasilopoulos

Abstract: Measuring a pulsar's rotational evolution is crucial to understanding the nature of the pulsar. Here we provide updated timing models for the rotational evolution of six pulsars, five of which are rotation phase-connected using primarily NICER X-ray data. For the newly-discovered fast energetic young pulsar, PSR J0058-7218, we increase the baseline of its timing model from 1.4 days to 8 months and… ▽ More Measuring a pulsar's rotational evolution is crucial to understanding the nature of the pulsar. Here we provide updated timing models for the rotational evolution of six pulsars, five of which are rotation phase-connected using primarily NICER X-ray data. For the newly-discovered fast energetic young pulsar, PSR J0058-7218, we increase the baseline of its timing model from 1.4 days to 8 months and not only measure more precisely its spin-down rate nudot = (-6.2324+/-0.0001)x10^-11 Hz s^-1 but also for the first time the second time derivative of spin rate nuddot = (4.2+/-0.2)x10^-21 Hz s^-2. For the fastest and most energetic young pulsar, PSR J0537-6910 (with 16 ms spin period), we detect 4 more glitches, for a total of 15 glitches over 4.5 years of NICER monitoring, and show that its spin-down behavior continues to set this pulsar apart from all others, including a long-term braking index n = -1.234+/-0.009 and interglitch braking indices that asymptote to <~ 7 for long times after a glitch. For PSR J1101-6101, we measure a much more accurate spin-down rate that agrees with a previous value measured without phase-connection. For PSR J1412+7922 (also known as Calvera), we extend the baseline of its timing model from our previous 1-year model to 4.4 years, and for PSR J1849-0001, we extend the baseline from 1.5 years to 4.7 years. We also present a long-term timing model of the energetic pulsar, PSR J1813-1749, by fitting previous radio and X-ray spin frequencies from 2009-2019 and new ones measured here using 2018 NuSTAR and 2021 Chandra data. △ Less

Submitted 29 July, 2022; v1 submitted 5 May, 2022; originally announced May 2022.

Comments: 18 pages, 17 figures; accepted for publication in ApJ; very minor edits and no changes to numbers

Journal ref: Astrophys. J. 939, 7 (2022)

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, A. Allocca, P. A. Altin, A. Amato , et al. (1636 additional authors not shown)

Abstract: Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational… ▽ More Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets. △ Less

Submitted 27 June, 2022; v1 submitted 21 December, 2021; originally announced December 2021.

Comments: 37 pages, 9 figures, submitted to ApJ

Report number: LIGO-P2100267

Journal ref: ApJ, 932, 133 (2022)

A Detection of Red Noise in PSR J1824$-$2452A and Projections for PSR B1937+21 using NICER X-ray Timing Data

Authors: Jeffrey S. Hazboun, Jack Crump, Andrea N. Lommen, Sergio Montano, Samantha J. H. Berry, Jesse Zeldes, Elizabeth Teng, Paul S. Ray, Matthew Kerr, Zaven Arzoumanian, Slavko Bogdanov, Julia Deneva, Natalia Lewandowska, Craig B. Markwardt, Scott Ransom, Teruaki Enoto, Kent S. Wood, Keith C. Gendreau, David A. Howe, Aditya Parthasarathy

Abstract: We have used X-ray data from the Neutron Star Interior Composition Explorer (NICER) to search for long time-scale, correlated variations ("red noise") in the pulse times of arrival from the millisecond pulsars PSR J1824$-$2452A and PSR B1937+21. These data more closely track intrinsic noise because X-rays are unaffected by the radio-frequency dependent propagation effects of the interstellar mediu… ▽ More We have used X-ray data from the Neutron Star Interior Composition Explorer (NICER) to search for long time-scale, correlated variations ("red noise") in the pulse times of arrival from the millisecond pulsars PSR J1824$-$2452A and PSR B1937+21. These data more closely track intrinsic noise because X-rays are unaffected by the radio-frequency dependent propagation effects of the interstellar medium. Our Bayesian search methodology yields strong evidence (natural log Bayes factor of $9.634 \pm 0.016$) for red noise in PSR J1824$-$2452A, but is inconclusive for PSR B1937+21. In the interest of future X-ray missions, we devise and implement a method to simulate longer and higher precision X-ray datasets to determine the timing baseline necessary to detect red noise. We find that the red noise in PSR B1937+21 can be reliably detected in a 5-year mission with a time-of-arrival (TOA) error of 2 microseconds and an observing cadence of 20 observations per month compared to the 5 microsecond TOA error and 11 observations per month that NICER currently achieves in PSR B1937+21. We investigate detecting red noise in PSR B1937+21 with other combinations of observing cadences and TOA errors. We also find that an injected stochastic gravitational wave background (GWB) with an amplitude of $A_{\rm GWB}=2\times10^{-15}$ and spectral index of $γ_{\rm GWB}=13/3$ can be detected in a pulsar with similar TOA precision to PSR B1937+21, but with no additional red noise, in a 10-year mission that observes the pulsar 15 times per month and has an average TOA error of 1 microsecond. △ Less

Submitted 3 December, 2021; originally announced December 2021.

Comments: 21 pages, 13 figures

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1672 additional authors not shown)

Abstract: We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both… ▽ More We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory. △ Less

Submitted 20 July, 2022; v1 submitted 25 November, 2021; originally announced November 2021.

Comments: 37 pages

Report number: LIGO-P2100049

Quasi-Simultaneous Radio/X-ray Observations of the Candidate Transitional Millisecond Pulsar 3FGL J1544.6-1125 During its Low-Luminosity Accretion-Disc State

Authors: Amruta Jaodand, Adam T. Deller, Nina Gusinskaia, Jason W. T. Hessels, James C. A. Miller-Jones, Anne M. Archibald, Slavko Bogdanov, Cees Bassa, Rudy Wijnands, Alessandro Patruno, Sotiris Sanidas

Abstract: 3FGL J1544.6-1125 is a candidate transitional millisecond pulsar (tMSP). Similar to the well-established tMSPs - PSR J1023+0038, IGR J18245-2452, and XSS J12270-4859 -- 3FGL J1544.6-1125 shows $γ$-ray emission and discrete X-ray "low" and "high" modes during its low-luminosity accretion-disk state. Coordinated radio/X-ray observations of PSR J1023+0038 in its current low-luminosity accretion-disk… ▽ More 3FGL J1544.6-1125 is a candidate transitional millisecond pulsar (tMSP). Similar to the well-established tMSPs - PSR J1023+0038, IGR J18245-2452, and XSS J12270-4859 -- 3FGL J1544.6-1125 shows $γ$-ray emission and discrete X-ray "low" and "high" modes during its low-luminosity accretion-disk state. Coordinated radio/X-ray observations of PSR J1023+0038 in its current low-luminosity accretion-disk state showed rapidly variable radio continuum emission-possibly originating from a compact, self-absorbed jet, the "propellering" of accretion material, and/or pulsar moding. 3FGL J1544.6-1125 is currently the only other (candidate) tMSP system in this state, and can be studied to see whether tMSPs are typically radio-loud compared to other neutron star binaries. In this work, we present a quasi-simultaneous Very Large Array and Swift radio/X-ray campaign on 3FGL J1544.6-1125. We detect 10 GHz radio emission varying in flux density from $47.7 \pm 6.0$ $μ$Jy down to $\sim$15 $μ$Jy (3$σ$ upper limit) at four epochs spanning three weeks. At the brightest epoch, the radio luminosity is $L_{5 GHz}$ $= (2.17 \pm 0.17) \times 10^{27}$ erg s$^{-1}$ for a quasi-simultaneous X-ray luminosity $L_{2-10 keV}$ $= (4.32 \pm 0.23) \times 10^{33}$ erg s$^{-1}$ (for an assumed distance of 3.8 kpc). These luminosities are close to those of PSR J1023+0038, and the results strengthen the case that 3FGL J1544.6-1125 is a tMSP showing similar phenomenology to PSR J1023+0038. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: Accepted for publication in the Astrophysical Journal

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

Authors: Thomas E. Riley, Anna L. Watts, Paul S. Ray, Slavko Bogdanov, Sebastien Guillot, Sharon M. Morsink, Anna V. Bilous, Zaven Arzoumanian, Devarshi Choudhury, Julia S. Deneva, Keith C. Gendreau, Alice K. Harding, Wynn C. G. Ho, James M. Lattimer, Michael Loewenstein, Renee M. Ludlam, Craig B. Markwardt, Takashi Okajima, Chanda Prescod-Weinstein, Ronald A. Remillard, Michael T. Wolff, Emmanuel Fonseca, H. Thankful Cromartie, Matthew Kerr, Timothy T. Pennucci , et al. (5 additional authors not shown)

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 (NICER XTI) event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint NANOGrav and CHIME/Pulsar wideban… ▽ More 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 (NICER XTI) event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint NANOGrav and CHIME/Pulsar wideband radio timing measurements of arXiv:2104.00880. We use XMM 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$_{\odot}$ 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. △ Less

Submitted 22 September, 2021; v1 submitted 14 May, 2021; originally announced May 2021.

Comments: 40 pages, 16 figures (3 of which are figure sets), 1 animation, 2 tables. ApJL accepted version. Software: https://github.com/ThomasEdwardRiley/xpsi. Zenodo: https://doi.org/10.5281/zenodo.4697624. (v1 = submitted version; v2 = accepted version; v3 = metadata edits)

Journal ref: The Astrophysical Journal Letters, Volume 918, Number 2, 2021

The Radius of PSR J0740+6620 from NICER and XMM-Newton Data

Authors: M. C. Miller, F. K. Lamb, A. J. Dittmann, S. Bogdanov, Z. Arzoumanian, K. C. Gendreau, S. Guillot, W. C. G. Ho, J. M. Lattimer, M. Loewenstein, S. M. Morsink, P. S. Ray, M. T. Wolff, C. L. Baker, T. Cazeau, S. Manthripragada, C. B. Markwardt, T. Okajima, S. Pollard, I. Cognard, H. T. Cromartie, E. Fonseca, L. Guillemot, M. Kerr, A. Parthasarathy , et al. (3 additional authors not shown)

Abstract: PSR J0740$+$6620 has a gravitational mass of $2.08\pm 0.07~M_\odot$, which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (N… ▽ More PSR J0740$+$6620 has a gravitational mass of $2.08\pm 0.07~M_\odot$, which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (NICER) and X-ray Multi-Mirror (XMM-Newton) X-ray observations. We find that the equatorial circumferential radius of PSR J0740$+$6620 is $13.7^{+2.6}_{-1.5}$ km (68%). We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030$+$0451, the masses of two other $\sim 2~M_\odot$ pulsars, and the tidal deformability constraints from two gravitational wave events, to three different frameworks for equation of state modeling, and find consistent results at $\sim 1.5-3$ times nuclear saturation density. For a given framework, when all measurements are included the radius of a $1.4~M_\odot$ neutron star is known to $\pm 4$% (68% credibility) and the radius of a $2.08~M_\odot$ neutron star is known to $\pm 5$%. The full radius range that spans the $\pm 1σ$ credible intervals of all the radius estimates in the three frameworks is $12.45\pm 0.65$ km for a $1.4~M_\odot$ neutron star and $12.35\pm 0.75$ km for a $2.08~M_\odot$ neutron star. △ Less

Submitted 14 May, 2021; originally announced May 2021.

Comments: 49 pages, 16 figures, submitted to The Astrophysical Journal Letters

NICER Detection of Thermal X-ray Pulsations from the Massive Millisecond Pulsars PSR J0740+6620 and PSR J1614-2230

Authors: M. T. Wolff, S. Guillot, S. Bogdanov, P. S. Ray, M. Kerr, Z. Arzoumanian, K. C. Gendreau, M. C. Miller, A. J. Dittmann, W. C. G. Ho, L. Guillemot, I. Cognard, G. Theureau, K. S. Wood

Abstract: We report the detection of X-ray pulsations from the rotation-powered millisecond-period pulsars PSR J0740+6620 and PSR J1614-2230, two of the most massive neutron stars known, using observations with the Neutron Star Interior Composition Explorer (NICER). We also analyze X-ray Multi-Mirror Mission (XMM-Newton) data for both pulsars to obtain their time-averaged fluxes and study their respective X… ▽ More We report the detection of X-ray pulsations from the rotation-powered millisecond-period pulsars PSR J0740+6620 and PSR J1614-2230, two of the most massive neutron stars known, using observations with the Neutron Star Interior Composition Explorer (NICER). We also analyze X-ray Multi-Mirror Mission (XMM-Newton) data for both pulsars to obtain their time-averaged fluxes and study their respective X-ray fields. PSR J0740+6620 exhibits a broad double-peaked profile with a separation of ~0.4 in phase. PSR J1614-2230, on the other hand, has a broad single-peak profile. The broad modulations with soft X-ray spectra of both pulsars are indicative of thermal radiation from one or more small regions of the stellar surface. We show the NICER detections of X-ray pulsations for both pulsars and also discuss the phase relationship to their radio pulsations. In the case of PSR J0740+6620, this paper documents the data reduction performed to obtain the pulsation detection and prepare for pulse profile modeling analysis. △ Less

Submitted 14 May, 2021; originally announced May 2021.

Comments: Submitted to The Astrophysical Journal Letters

Constraining the Neutron Star Mass--Radius Relation and Dense Matter Equation of State with NICER. III. Model Description and Verification of Parameter Estimation Codes

Authors: Slavko Bogdanov, Alexander J. Dittmann, Wynn C. G. Ho, Frederick K. Lamb, Simin Mahmoodifar, M. Coleman Miller, Sharon M. Morsink, Thomas E. Riley, Tod E. Strohmayer, Anna L. Watts, Devarshi Choudhury, Sebastien Guillot, Alice K. Harding, Paul S. Ray, Zorawar Wadiasingh, Michael T. Wolff, Craig B. Markwardt, Zaven Arzoumanian, Keith C. Gendreau

Abstract: We describe the X-ray pulse profile models we use, and how we use them, to analyze Neutron Star Interior Composition Explorer (NICER) observations of rotation-powered millisecond pulsars to obtain information about the mass-radius relation of neutron stars and the equation of state of the dense matter in their cores. Here we detail our modeling of the observed profile of PSR J0030+0451 that we ana… ▽ More We describe the X-ray pulse profile models we use, and how we use them, to analyze Neutron Star Interior Composition Explorer (NICER) observations of rotation-powered millisecond pulsars to obtain information about the mass-radius relation of neutron stars and the equation of state of the dense matter in their cores. Here we detail our modeling of the observed profile of PSR J0030+0451 that we analyzed in Miller et al. (2019) and Riley et al. (2019) and describe a cross-verification of computations of the pulse profiles of a star with R/M 3, in case stars this compact need to be considered in future analyses. We also present our early cross-verification efforts of the parameter estimation procedures used by Miller et al. (2019) and Riley et al. (2019) by analyzing two distinct synthetic data sets. Both codes yielded credible regions in the mass-radius plane that are statistically consistent with one another and both gave posterior distributions for model parameter values consistent with the values that were used to generate the data. We also summarize the additional tests of the parameter estimation procedure of Miller et al. (2019) that used synthetic pulse profiles and the NICER pulse profile of PSR J0030+0451. We then illustrate how the precision of mass and radius estimates depends on the pulsar's spin rate and the size of its hot spot by analyzing four different synthetic pulse profiles. Finally, we assess possible sources of systematic error in these estimates made using this technique, some of which may warrant further investigation. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: 26 pages, 5 figures; submitted for publication in the Astrophysical Journal Letters

Enhanced X-ray Emission Coinciding with Giant Radio Pulses from the Crab Pulsar

Authors: Teruaki Enoto, Toshio Terasawa, Shota Kisaka, Chin-Ping Hu, Sebastien Guillot, Natalia Lewandowska, Christian Malacaria, Paul S. Ray, Wynn C. G. Ho, Alice K. Harding, Takashi Okajima, Zaven Arzoumanian, Keith C. Gendreau, Zorawar Wadiasingh, Craig B. Markwardt, Yang Soong, Steve Kenyon, Slavko Bogdanov, Walid A. Majid, Tolga Güver, Gaurava K. Jaisawal, Rick Foster, Yasuhiro Murata, Hiroshi Takeuchi, Kazuhiro Takefuji , et al. (12 additional authors not shown)

Abstract: Giant radio pulses (GRPs) are sporadic bursts emitted by some pulsars, lasting a few microseconds. GRPs are hundreds to thousands of times brighter than regular pulses from these sources. The only GRP-associated emission outside radio wavelengths is from the Crab Pulsar, where optical emission is enhanced by a few percent during GRPs. We observed the Crab Pulsar simultaneously at X-ray and radio w… ▽ More Giant radio pulses (GRPs) are sporadic bursts emitted by some pulsars, lasting a few microseconds. GRPs are hundreds to thousands of times brighter than regular pulses from these sources. The only GRP-associated emission outside radio wavelengths is from the Crab Pulsar, where optical emission is enhanced by a few percent during GRPs. We observed the Crab Pulsar simultaneously at X-ray and radio wavelengths, finding enhancement of the X-ray emission by $3.8\pm0.7\%$ (a 5.4$σ$ detection) coinciding with GRPs. This implies that the total emitted energy from GRPs is tens to hundreds of times higher than previously known. We discuss the implications for the pulsar emission mechanism and extragalactic fast radio bursts. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 63 pages, 23 figures, 7 Tables. Published in Science

Discovery of UV millisecond pulsations and moding in the low mass X-ray binary state of transitional millisecond pulsar J1023+0038

Authors: Amruta D. Jaodand, Juan V. Hernández Santisteban, Anne M. Archibald, Jason W. T. Hessels, Slavko Bogdanov, Christian Knigge, Nathalie Degenaar, Adam T. Deller, Simone Scaringi, Alessandro Patruno

Abstract: PSR J1023+0038 is a rapidly-spinning neutron star with a low-mass-binary companion that switches between a radio pulsar and low-luminosity disk state. In 2013, it transitioned to its current disk state accompanied by brightening at all observed wavelengths. In this state, PSR J1023+0038 now shows optical and X-ray pulsations and abrupt X-ray luminosity switches between discrete 'low' and 'high' mo… ▽ More PSR J1023+0038 is a rapidly-spinning neutron star with a low-mass-binary companion that switches between a radio pulsar and low-luminosity disk state. In 2013, it transitioned to its current disk state accompanied by brightening at all observed wavelengths. In this state, PSR J1023+0038 now shows optical and X-ray pulsations and abrupt X-ray luminosity switches between discrete 'low' and 'high' modes. Continuum radio emission, denoting an outflow, is also present and brightens during the X-ray low modes. Here, we present a simultaneous optical, ultraviolet (UV) and X-ray campaign comprising Kepler ($400-800$ nm), Hubble Space Telescope ($180-280$ nm), XMM-Newton ($0.3-10$ keV) and NuSTAR ($3 - 79$ keV). We demonstrate that low and high luminosity modes in the UV band are strictly simultaneous with the X-ray modes and change the UV brightness by a factor of $\sim25$\% on top of a much brighter persistent UV component. We find strong evidence for UV pulsations (pulse fraction of $0.82\pm0.19$\%) in the high-mode, with a similar waveform as the X-ray pulsations making it the first known UV millisecond pulsar. Lastly, we find that the optical mode changes occur synchronously with the UV/X-ray mode changes, but optical modes are inverted compared to the higher frequencies. There appear to be two broad-band emission components: one from radio to near-infrared/optical that is brighter when the second component from optical to hard X-rays is dimmer (and vice-versa). We suggest that these components trace switches between accretion into the neutron star magnetosphere (high-energy high-mode) versus ejection of material (low-energy high-mode). Lastly, we propose that optical/UV/X-ray pulsations can arise from a shocked accretion flow channeled by the neutron star's magnetic field. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Comments: 18 pages, 10 figures

A Deep Chandra X-ray Observatory Study of the Millisecond Pulsar Population in the Globular Cluster Terzan 5

Authors: Slavko Bogdanov, Arash Bahramian, Craig O. Heinke, Scott M. Ransom

Abstract: We present an analysis of 745.6 ks of archival Chandra X-ray Observatory Advanced CCD Imaging Spectrometer data accumulated between 2000 and 2016 of the millisecond pulsar (MSP) population in the rich Galactic globular cluster Terzan 5. Eight of the 37 MSPs with precise positions are found to have plausible X-ray source matches. Despite the deep exposure, the remaining MSPs are either marginally d… ▽ More We present an analysis of 745.6 ks of archival Chandra X-ray Observatory Advanced CCD Imaging Spectrometer data accumulated between 2000 and 2016 of the millisecond pulsar (MSP) population in the rich Galactic globular cluster Terzan 5. Eight of the 37 MSPs with precise positions are found to have plausible X-ray source matches. Despite the deep exposure, the remaining MSPs are either marginally detected or have no obvious X-ray counterparts, which can be attributed to the typically soft thermal spectra of rotation-powered MSPs, which are strongly attenuated by the high intervening absorbing column (~$10^{22}$ cm$^{-2}$) towards the cluster, and in some instances severe source crowding/blending. For the "redback" MSP binaries, PSRs J1748-2446P and J1748-2446ad, and the "black widow" binary PSRs J1748-2446O, we find clear evidence for large-amplitude X-ray variability at the orbital period consistent with an intrabinary shock origin. The third redback MSP in the cluster, PSR J1748-2446A, shows large amplitude variations in flux on time scales of years, possibility due to state transitions or intense flaring episodes from the secondary star. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: 14 pages, 6 figures; submitted to the Astrophysical Journal

A NICER View of Spectral and Profile Evolution for Three X-ray Emitting Millisecond Pulsars

Authors: Dominick M. Rowan, Zaynab Ghazi, Lauren Lugo, Elizabeth Spano, Andrea Lommen, Alice Harding, Christo Venter, Renee Ludlam, Paul S. Ray, Matthew Kerr, Zaven Arzoumanian, Slavko Bogdanov, Julia Deneva, Sebastien Guillot, Natalia Lewandowska, Craig B. Markwardt, Scott Ransom, Teruaki Enoto, Kent S. Wood, Keith Gendreau

Abstract: 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 separatio… ▽ More 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σ$ 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. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: 16 pages, 13 figures, accepted for publication in ApJ

NICER X-ray Observations of Seven Nearby Rotation-Powered Millisecond Pulsars

Authors: S. Guillot, M. Kerr, P. S. Ray, S. Bogdanov, S. Ransom, J. S. Deneva, Z. Arzoumanian, P. Bult, D. Chakrabarty, K. C. Gendreau, W. C. G. Ho, G. K. Jaisawal, C. Malacaria, M. C. Miller, T. E. Strohmayer, M. T. Wolff, K. S. Wood, N. A. Webb, L. Guillemot, I. Cognard, G. Theureau

Abstract: 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σ$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sam… ▽ More 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σ$ 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 millisecond pulsars exhibit a significant phase offset between their X-ray and radio pulses. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 25 pages, 11 tables, 4 figures. In press in The Astrophysical Journal Letters

Journal ref: ApJL, 887, L27 (2019)

Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. II. Emission from Hot Spots on a Rapidly Rotating Neutron Star

Authors: Slavko Bogdanov, Frederick K. Lamb, Simin Mahmoodifar, M. Coleman Miller, Sharon M. Morsink, Thomas E. Riley, Tod E. Strohmayer, Anna L. Watts, Alexander J. Dittmann, Deepto Chakrabarty, Sebastien Guillot, Zaven Arzoumanian, Keith C. Gendreau

Abstract: We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verifica… ▽ More We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verification tests between our numerical codes for both the Schwarzschild + Doppler and Oblate Schwarzschild approximations, and compare both approximations against exact numerical calculations. We find superb agreement between the code outputs, as well as in comparison against a set of analytical and semi-analytical calculations, which combined with their speed, demonstrates that the codes are well-suited for large-scale statistical sampling applications. A set of verified, high-precision reference synthetic pulse profiles is provided to the community to facilitate testing of other independently developed codes. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 31 pages, 14 figures; accepted for publication in the Astrophysical Journal Letters

Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. I. The Millisecond Pulsar X-ray Data Set

Authors: Slavko Bogdanov, Sebastien Guillot, Paul S. Ray, Michael T. Wolff, Deepto Chakrabarty, Wynn C. G. Ho, Matthew Kerr, Frederick K. Lamb, Andrea Lommen, Renee M. Ludlam, Reilly Milburn, Sergio Montano, M. Coleman Miller, Michi Baubock, Feryal Ozel, Dimitrios Psaltis, Ronald A. Remillard, Thomas E. Riley, James F. Steiner, Tod E. Strohmayer, Anna L. Watts, Kent S. Wood, Jesse Zeldes, Teruaki Enoto, Takashi Okajima , et al. (5 additional authors not shown)

Abstract: 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 i… ▽ More 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. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 28 pages, 9 figures; accepted for publication in the Astrophysical Journal Letters

PSR J0030+0451 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter

Authors: M. C. Miller, F. K. Lamb, A. J. Dittmann, S. Bogdanov, Z. Arzoumanian, K. C. Gendreau, S. Guillot, A. K. Harding, W. C. G. Ho, J. M. Lattimer, R. M. Ludlam, S. Mahmoodifar, S. M. Morsink, P. S. Ray, T. E. Strohmayer, K. S. Wood, T. Enoto, R. Foster, T. Okajima, G. Prigozhin, Y. Soong

Abstract: 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 mas… ▽ More 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. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 49 pages, 16 figures, part of The Astrophysical Journal Letters focus issue on the Neutron Star Interior Composition Explorer

Journal ref: ApJL 2019, 887, L24

A NICER view of PSR J0030+0451: evidence for a global-scale multipolar magnetic field

Authors: Anna V. Bilous, Anna L. Watts, Alice K. Harding, Thomas E. Riley, Zaven Arzoumanian, Slavko Bogdanov, Keith C. Gendreau, Paul S. Ray, Sebastien Guillot, Wynn C. G. Ho, Deepto Chakrabarty

Abstract: 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 exa… ▽ More 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. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: Appears in ApJ Letters Focus Issue on NICER Constraints on the Dense Matter Equation of State

Journal ref: Bilous et al. 2019, ApJ Letters, 887, L23

A NICER view of PSR J0030+0451: Implications for the dense matter equation of state

Authors: G. Raaijmakers, T. E. Riley, A. L. Watts, S. K. Greif, S. M. Morsink, K. Hebeler, A. Schwenk, T. Hinderer, S. Nissanke, S. Guillot, Z. Arzoumanian, S. Bogdanov, D. Chakrabarty K. C. Gendreau, W. C. G. Ho, J. M. Lattimer, R. M. Ludlam, M. T. Wolff

Abstract: Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. U… ▽ More Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EOS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EOS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EOS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass-radius likelihood function derived by Riley et al. (2019) using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: Appears in ApJ Letters Focus Issue on NICER Constraints on the Dense Matter Equation of State, 17 pages, 5 figures

Journal ref: ApJ Letters, 887, L22 (2019)

A NICER View of PSR J0030+0451: Millisecond Pulsar Parameter Estimation

Authors: Thomas E. Riley, Anna L. Watts, Slavko Bogdanov, Paul S. Ray, Renee M. Ludlam, Sebastien Guillot, Zaven Arzoumanian, Charles L. Baker, Anna V. Bilous, Deepto Chakrabarty, Keith C. Gendreau, Alice K. Harding, Wynn C. G. Ho, James M. Lattimer, Sharon M. Morsink, Tod E. Strohmayer

Abstract: 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 p… ▽ More 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. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: Appears in ApJ Letters Focus Issue on NICER Constraints on the Dense Matter Equation of State; 76 pages, 24 figures, 7 tables, 8 figure sets (available in the online journal or from the authors)

Journal ref: ApJL, 887, L21 (2019)

The SUrvey for Pulsars and Extragalactic Radio Bursts IV: Discovery and polarimetry of a 12.1-second radio pulsar

Authors: V. Morello, E. F. Keane, T. Enoto, S. Guillot, W. C. G. Ho, A. Jameson, M. Kramer, B. W. Stappers, M. Bailes, E. D. Barr, S. Bhandari, M. Caleb, C. M. L. Flynn, F. Jankowski, S. Johnston, W. van Straten, Z. Arzoumanian, S. Bogdanov, K. C. Gendreau, C. Malacaria, P. S. Ray, R. A. Remillard

Abstract: 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 e… ▽ More 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. △ Less

Submitted 31 January, 2020; v1 submitted 9 October, 2019; originally announced October 2019.

Comments: Accepted for publication in MNRAS, 13 pages, 8 figures

Radio and X-ray monitoring of the accreting millisecond X-ray pulsar IGR J17591-2342 in outburst

Authors: N. V. Gusinskaia, T. D. Russell, J. W. T. Hessels, S. Bogdanov, N. Degenaar, A. T. Deller, J. van den Eijnden, A. D. Jaodand, J. C. A. Miller-Jones, R. Wijnands

Abstract: IGR J17591$-$2342 is a new accreting millisecond X-ray pulsar (AMXP) that was recently discovered in outburst in 2018. Early observations revealed that the source's radio emission is brighter than that of any other known neutron star low-mass X-ray binary (NS-LMXB) at comparable X-ray luminosity, and assuming its likely $\gtrsim 6$ kpc distance. It is comparably radio bright to black hole LMXBs at… ▽ More IGR J17591$-$2342 is a new accreting millisecond X-ray pulsar (AMXP) that was recently discovered in outburst in 2018. Early observations revealed that the source's radio emission is brighter than that of any other known neutron star low-mass X-ray binary (NS-LMXB) at comparable X-ray luminosity, and assuming its likely $\gtrsim 6$ kpc distance. It is comparably radio bright to black hole LMXBs at similar X-ray luminosities. In this work, we present the results of our extensive radio and X-ray monitoring campaign of the 2018 outburst of IGR J17591$-$2342. In total we collected 10 quasi-simultaneous radio (VLA, ATCA) and X-ray (Swift-XRT) observations, which make IGR J17591$-$2342 one of the best-sampled NS-LMXBs. We use these to fit a power-law correlation index $β= 0.37^{+0.42}_{-0.40}$ between observed radio and X-ray luminosities ( $L_\mathrm{R}\propto L_\mathrm{X}^β$). However, our monitoring revealed a large scatter in IGR J17591$-$2342's radio luminosity (at a similar X-ray luminosity, $L_\mathrm{X} \sim 10^{36}$ erg s$^{-1}$, and spectral state), with $L_\mathrm{R} \sim 4 \times 10^{29}$ erg s$^{-1}$ during the first three reported observations, and up to a factor of 4 lower $L_\mathrm{R}$ during later radio observations. Nonetheless, the average radio luminosity of IGR J17591$-$2342 is still one of the highest among NS-LMXBs, and we discuss possible reasons for the wide range of radio luminosities observed in such systems during outburst. We found no evidence for radio pulsations from IGR J17591$-$2342 in our Green Bank Telescope observations performed shortly after the source returned to quiescence. Nonetheless, we cannot rule out that IGR J17591$-$2342 becomes a radio millisecond pulsar during quiescence. △ Less

Submitted 6 December, 2019; v1 submitted 5 September, 2019; originally announced September 2019.

Comments: 12 pages, 3 figures, 2 tables, accepted for publication in MNRAS

A joint NICER and XMM-Newton view of the "Magnificent" thermally emitting X-ray Isolated Neutron Star RX J1605.3+3249

Authors: Christian Malacaria, Slavko Bogdanov, Wynn C. G. Ho, Teruaki Enoto, Paul S. Ray, Zaven Arzoumanian, Thoniel Cazeau, Keith C. Gendreau, Sebastien Guillot, Tolga Guver, Gaurava K. Jaisawal, Michael T. Wolff

Abstract: 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… ▽ More 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σ) 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. △ Less

Submitted 6 June, 2019; originally announced June 2019.

Comments: Accepted for publication in ApJ

Discovery of Soft X-ray Pulsations from PSR J1231-1411 using NICER

Authors: Paul S. Ray, Sebastien Guillot, Scott M. Ransom, Matthew Kerr, Slavko Bogdanov, Alice K. Harding, Michael T. Wolff, Christian Malacaria, Keith C. Gendreau, Zaven Arzoumanian, Craig Markwardt, Yang Soong, John P. Doty

Abstract: 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… ▽ More 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 $γ$ rays provides insight into the geometry of the system. △ Less

Submitted 28 May, 2019; originally announced May 2019.

Comments: ApJ Letters, in press

Journal ref: ApJL 878 L22 (2019)

Determining the Equation of State of Cold, Dense Matter with X-ray Observations of Neutron Stars

Authors: Slavko Bogdanov, Anna L. Watts, Deepto Chakrabarty, Zaven Arzoumanian, Sebastien Guillot, Keith C. Gendreau, Frederick K. Lamb, Thomas Maccarone, M. Coleman Miller, Feryal Ozel, Paul S. Ray, Colleen A. Wilson-Hodge

Abstract: 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.… ▽ More 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. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Submitted in response to the Astronomy and Astrophysics Decadal Survey (Astro2020) call for Science White Papers; 8 pages, 2 figures

STROBE-X: X-ray Timing and Spectroscopy on Dynamical Timescales from Microseconds to Years

Authors: Paul S. Ray, Zaven Arzoumanian, David Ballantyne, Enrico Bozzo, Soren Brandt, Laura Brenneman, Deepto Chakrabarty, Marc Christophersen, Alessandra DeRosa, Marco Feroci, Keith Gendreau, Adam Goldstein, Dieter Hartmann, Margarita Hernanz, Peter Jenke, Erin Kara, Tom Maccarone, Michael McDonald, Michael Nowak, Bernard Phlips, Ron Remillard, Abigail Stevens, John Tomsick, Anna Watts, Colleen Wilson-Hodge , et al. (134 additional authors not shown)

Abstract: 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 techniqu… ▽ More 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. △ Less

Submitted 8 March, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 50 pages, Probe class mission concept study report submitted to NASA for Astro2020 Decadal Survey

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015-2017 LIGO Data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, A. Ananyeva , et al. (1127 additional authors not shown)

Abstract: We present a search for gravitational waves from 221 pulsars with rotation frequencies $\gtrsim 10$ Hz. We use advanced LIGO data from its first and second observing runs spanning 2015-2017, which provides the highest-sensitivity gravitational-wave data so far obtained. In this search we target emission from both the $l = m = 2$ mass quadrupole mode, with a frequency at twice that of the pulsar's… ▽ More We present a search for gravitational waves from 221 pulsars with rotation frequencies $\gtrsim 10$ Hz. We use advanced LIGO data from its first and second observing runs spanning 2015-2017, which provides the highest-sensitivity gravitational-wave data so far obtained. In this search we target emission from both the $l = m = 2$ mass quadrupole mode, with a frequency at twice that of the pulsar's rotation, and from the $l = 2$, $m = 1$ mode, with a frequency at the pulsar rotation frequency. The search finds no evidence for gravitational-wave emission from any pulsar at either frequency. For the $l = m = 2$ mode search, we provide updated upper limits on the gravitational-wave amplitude, mass quadrupole moment, and fiducial ellipticity for 167 pulsars, and the first such limits for a further 55. For 20 young pulsars these results give limits that are below those inferred from the pulsars' spin-down. For the Crab and Vela pulsars our results constrain gravitational-wave emission to account for less than 0.017% and 0.18% of the spin-down luminosity, respectively. For the recycled millisecond pulsar J0711-6830 our limits are only a factor of 1.3 above the spin-down limit, assuming the canonical value of $10^{38}$ kg m$^2$ for the star's moment of inertia, and imply a gravitational-wave-derived upper limit on the star's ellipticity of $1.2\!\times\!10^{-8}$. We also place new limits on the emission amplitude at the rotation frequency of the pulsars. △ Less

Submitted 28 July, 2020; v1 submitted 22 February, 2019; originally announced February 2019.

Comments: Main paper: 39 pages, 9 figures, 2 tables. Accepted for publication in The Astrophysical Journal. Results from the tables are contained in the ancillary file Erratum: 11 pages, 4 figures, 1 table

Report number: LIGO-P1800344

Journal ref: The Astrophysical Journal 879 (2019) 10

High-Precision X-ray Timing of Three Millisecond Pulsars with NICER: Stability Estimates and Comparison with Radio

Authors: J. S. Deneva, P. S. Ray, A. Lommen, S. M. Ransom, S. Bogdanov, M. Kerr, K. S. Wood, Z. Arzoumanian, K. Black, J. Doty, K. C. Gendreau, S. Guillot, A. Harding, N. Lewandowska, C. Malacaria, C. B. Markwardt, S. Price, L. Winternitz, M. T. Wolff, L. Guillemot, I. Cognard, P. T. Baker, H. Blumer, P. R. Brook, H. T. Cromartie , et al. (23 additional authors not shown)

Abstract: 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 pul… ▽ More 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 $σ_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 $σ_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 $σ_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. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Comments: 30 pages, 15 figures; accepted to ApJ

Neutron Star Interior Composition Explorer X-ray Timing of the Radio and Gamma-ray Quiet Pulsars PSR J1412+7922 AND PSR J1849-0001

Authors: Slavko Bogdanov, Wynn C. G. Ho, Teruaki Enoto, Sebastien Guillot, Alice K. Harding, Gaurava K. Jaisawal, Christian Malacaria, Sridhar S. Manthripragada, Zaven Arzoumanian, Keith C. Gendreau

Abstract: We present new timing and spectral analyses of PSR J1412+7922 (Calvera) and PSR J1849-0001, which are only seen as pulsars in X-rays, based on observations conducted with the Neutron Star Interior Composition Explorer (NICER). We obtain updated and substantially improved pulse ephemerides compared to previous X-ray studies, as well as spectra that can be well-fit by simple blackbodies and/or a pow… ▽ More We present new timing and spectral analyses of PSR J1412+7922 (Calvera) and PSR J1849-0001, which are only seen as pulsars in X-rays, based on observations conducted with the Neutron Star Interior Composition Explorer (NICER). We obtain updated and substantially improved pulse ephemerides compared to previous X-ray studies, as well as spectra that can be well-fit by simple blackbodies and/or a power law. Our refined timing measurements enable deeper searches for pulsations at other wavelengths and sensitive targeted searches by LIGO/Virgo for continuous gravitational waves from these neutron stars. Using the sensitivity of LIGO's first observing run, we estimate constraints that a gravitational wave search of these pulsars would obtain on the size of their mass deformation and r-mode fluid oscillation. △ Less

Submitted 4 February, 2019; v1 submitted 31 January, 2019; originally announced February 2019.

Comments: 12 pages, 4 figures; submitted to the Astrophysical Journal

Dense matter with eXTP

Authors: Anna L. Watts, Wenfei Yu, Juri Poutanen, Shu Zhang, Sudip Bhattacharyya, Slavko Bogdanov, Long Ji, Alessandro Patruno, Thomas E. Riley, Pavel Bakala, Altan Baykal, Federico Bernardini, Ignazio Bombaci, Edward Brown, Yuri Cavecchi, Deepto Chakrabarty, Jérôme Chenevez, Nathalie Degenaar, Melania Del Santo, Tiziana Di Salvo, Victor Doroshenko, Maurizio Falanga, Robert D. Ferdman, Marco Feroci, Angelo F. Gambino , et al. (51 additional authors not shown)

Abstract: In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, b… ▽ More In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

Anti-glitches in the Ultraluminous Accreting Pulsar NGC 300 ULX-1 Observed with NICER

Authors: Paul S. Ray, Sebastien Guillot, Wynn C. G. Ho, Matthew Kerr, Teruaki Enoto, Keith C. Gendreau, Zaven Arzoumanian, Diego Altamirano, Slavko Bogdanov, Robert Campion, Deepto Chakrabarty, Gaurava K. Jaisawal, Robert Kozon, Christian Malacaria, Tod E. Strohmayer, Michael T. Wolff

Abstract: 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 $Δν= -23, -30,$ and $-43 \,μ$Hz (fractional amplitudes $Δν/ν= -4.4, -5.5,$ and $-7.7 \times 10^{-4}$). We deter… ▽ More 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 $Δν= -23, -30,$ and $-43 \,μ$Hz (fractional amplitudes $Δν/ν= -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. △ Less

Submitted 28 May, 2019; v1 submitted 22 November, 2018; originally announced November 2018.

Comments: ApJ, accepted

Variable Heating and Flaring of Three Redback Millisecond Pulsar Companions

Authors: Patricia B. Cho, Jules P. Halpern, Slavko Bogdanov

Abstract: We are monitoring established and putative redback millisecond pulsars (MSPs) in time-series photometry, repeatedly covering their 5-6 hr orbital light curves in $r^{\prime}$ or $R$. On timescales of months, PSR J1048+2339 and XMMU J083850.38$-$282756.8 exhibit similar variability of $\approx0.3$ mag on the heated side of the companion star. However, the heating light curve is rarely symmetric, su… ▽ More We are monitoring established and putative redback millisecond pulsars (MSPs) in time-series photometry, repeatedly covering their 5-6 hr orbital light curves in $r^{\prime}$ or $R$. On timescales of months, PSR J1048+2339 and XMMU J083850.38$-$282756.8 exhibit similar variability of $\approx0.3$ mag on the heated side of the companion star. However, the heating light curve is rarely symmetric, suggesting that the intrabinary shock generated by the pulsar wind is skewed in addition to being variable, or that changing magnetic fields intrinsic to the companion channel the pulsar wind. In addition to this variable heating, there are long-lived flaring states that increase the brightness by an additional 0.5 mag, with variability on $\approx10$ minute timescales. These flares also appear to originate on the heated side of the companion, while the "night"-side brightness remains relatively stable. Somewhat less active, PSR J1628$-$3205 has an optical light curve that is dominated by tidal distortion (ellipsoidal modulation), although it too shows evidence of variable and asymmetric heating due to shifting magnetic fields or migrating star spots. These effects frustrate any effort to derive system parameters such as inclination angle and Roche-lobe filling factor from optical light curves of redback MSPs. We also report on two Chandra X-ray observations of PSR J1048+2339 that show strong orbital modulation, possibly due to beaming along the intrabinary shock, and a third observation that is dominated by flaring. The peak flare luminosity in the 0.3-8 keV band is $\approx12\%$ of the pulsar's spin-down power, which may require magnetic reconnection. None of these three systems has yet shown a transition back to an accreting state. △ Less

Submitted 1 September, 2018; originally announced September 2018.

Comments: 11 pages, 6 figures, to appear in The Astrophysical Journal

NuSTAR and NICER reveal IGR J17591-2342 as a new accreting millisecond X-ray pulsar

Authors: A. Sanna, C. Ferrigno, P. S. Ray, L. Ducci, G. K. Jaisawal, T. Enoto, E. Bozzo, D. Altamirano, T. Di Salvo, T. E. Strohmayer, A. Papitto, A. Riggio, L. Burderi, P. M. Bult, S. Bogdanov, A. F. Gambino, A. Marino, R. Iaria, Z. Arzoumanian, D. Chakrabarty, K. C. Gendreau, S. Guillot, C. B. Markwardt, M. T. Wolff

Abstract: 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… ▽ More 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. △ Less

Submitted 30 August, 2018; originally announced August 2018.

Comments: 6 pages, 5 figures, 2 tables. Submitted to A&A

Journal ref: A&A 617, L8 (2018)

NICER Discovers the Ultracompact Orbit of the Accreting Millisecond Pulsar IGR J17062-6143

Authors: Tod E. Strohmayer, Zaven Arzoumanian, Slavko Bogdanov, Peter M. Bult, Deepto Chakrabarty, Teruaki Enoto, Keith C. Gendreau, Sebastien Guillot, Alice K. Harding, Wynn C. G. Ho, Jeroen Homan, Gaurava K. Jaisawal, Laurens Keek, Matthew Kerr, Simin Mahmoodifar, Craig B. Markwardt, Scott M. Ransom, Paul S. Ray, Ron Remillard, Michael T. Wolff

Abstract: 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… ▽ More 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). △ Less

Submitted 13 August, 2018; originally announced August 2018.

Comments: 19 pages, 7 figures, published in the Astrophysical Journal Letters

Journal ref: The Astrophysical Journal Letters, Volume 858, Issue 2, article id. L13, 7 pp. (2018)

The implementation of a Fast-Folding pipeline for long-period pulsar searching in the PALFA survey

Authors: E. Parent, V. M. Kaspi, S. M. Ransom, M. Krasteva, C. Patel, P. Scholz, A. Brazier, M. A. McLaughlin, M. Boyce, W. W. Zhu, Z. Pleunis, B. Allen, S. Bogdanov, K. Caballero, F. Camilo, R. Camuccio, S. Chatterjee, J. M. Cordes, F. Crawford, J. S. Deneva, R. Ferdman, P. C. C. Freire, J. W. T. Hessels, F. A. Jenet, B. Knispel , et al. (9 additional authors not shown)

Abstract: The Pulsar Arecibo L-Band Feed Array (PALFA) survey, the most sensitive blind search for radio pulsars yet conducted, is ongoing at the Arecibo Observatory in Puerto Rico. The vast majority of the 180 pulsars discovered by PALFA have spin periods shorter than 2 seconds. Pulsar surveys may miss long-period radio pulsars due to the summing of a finite number of harmonic components in conventional Fo… ▽ More The Pulsar Arecibo L-Band Feed Array (PALFA) survey, the most sensitive blind search for radio pulsars yet conducted, is ongoing at the Arecibo Observatory in Puerto Rico. The vast majority of the 180 pulsars discovered by PALFA have spin periods shorter than 2 seconds. Pulsar surveys may miss long-period radio pulsars due to the summing of a finite number of harmonic components in conventional Fourier analyses (typically $\sim$16), or due to the strong effect of red noise at low modulation frequencies. We address this reduction in sensitivity by using a time-domain search technique: the Fast-Folding Algorithm (FFA). We designed a program that implements a FFA-based search in the PALFA processing pipeline, and tested the efficiency of the algorithm by performing tests under both ideal, white noise conditions, as well as with real PALFA observational data. In the two scenarios, we show that the time-domain algorithm has the ability to outperform the FFT-based periodicity search implemented in the survey. We perform simulations to compare the previously reported PALFA sensitivity with that obtained using our new FFA implementation. These simulations show that for a pulsar having a pulse duty cycle of roughly 3%, the performance of our FFA pipeline exceeds that of our FFT pipeline for pulses with DM $\lesssim$ 40 pc cm$^{-3}$ and for periods as short as $\sim$500 ms, and that the survey sensitivity is improved by at least a factor of two for periods $\gtrsim$ 6 sec. Discoveries from the implementation of the algorithm in PALFA are also presented in this paper. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Comments: Accepted for publication in The Astrophysical Journal; 19 pages, 9 figures, 2 tables