Showing 1–50 of 112 results for author: Watts, A L

Constraining the dense matter equation of state with new NICER mass-radius measurements and new chiral effective field theory inputs

Authors: Nathan Rutherford, Melissa Mendes, Isak Svensson, Achim Schwenk, Anna L. Watts, Kai Hebeler, Jonas Keller, Chanda Prescod-Weinstein, Devarshi Choudhury, Geert Raaijmakers, Tuomo Salmi, Patrick Timmerman, Serena Vinciguerra, Sebastien Guillot, James M. Lattimer

Abstract: Pulse profile modeling of X-ray data from NICER is now enabling precision inference of neutron star mass and radius. Combined with nuclear physics constraints from chiral effective field theory ($χ$EFT), and masses and tidal deformabilities inferred from gravitational wave detections of binary neutron star mergers, this has lead to a steady improvement in our understanding of the dense matter equa… ▽ More Pulse profile modeling of X-ray data from NICER is now enabling precision inference of neutron star mass and radius. Combined with nuclear physics constraints from chiral effective field theory ($χ$EFT), and masses and tidal deformabilities inferred from gravitational wave detections of binary neutron star mergers, this has lead to a steady improvement in our understanding of the dense matter equation of state (EOS). Here we consider the impact of several new results: the radius measurement for the 1.42$\,M_\odot$ pulsar PSR J0437$-$4715 presented by Choudhury et al. (2024), updates to the masses and radii of PSR J0740$+$6620 and PSR J0030$+$0451, and new $χ$EFT results for neutron star matter up to 1.5 times nuclear saturation density. Using two different high-density EOS extensions -- a piecewise-polytropic (PP) model and a model based on the speed of sound in a neutron star (CS) -- we find the radius of a 1.4$\,M_\odot$ (2.0$\,M_\odot$) neutron star to be constrained to the 95\% credible ranges $12.28^{+0.50}_{-0.76}\,$km ($12.33^{+0.70}_{-1.34}\,$km) for the PP model and $12.01^{+0.56}_{-0.75}\,$km ($11.55^{+0.94}_{-1.09}\,$km) for the CS model. The maximum neutron star mass is predicted to be $2.15^{+0.14}_{-0.16}\,$$M_\odot$ and $2.08^{+0.28}_{-0.16}\,$$M_\odot$ for the PP and CS model, respectively. We explore the sensitivity of our results to different orders and different densities up to which $χ$EFT is used, and show how the astrophysical observations provide constraints for the pressure at intermediate densities. Moreover, we investigate the difference $R_{2.0} - R_{1.4}$ of the radius of 2$\,M_\odot$ and 1.4$\,M_\odot$ neutron stars within our EOS inference. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 23 pages, 15 figures; accepted for publication in ApJ Letters. This paper was made to be as similar to the ApJ Letters version as possible

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 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)

X-PSI Parameter Recovery for Temperature Map Configurations Inspired by PSR J0030+0451

Authors: Serena Vinciguerra, Tuomo Salmi, Anna L. Watts, Devarshi Choudhury, Yves Kini, Thomas E. Riley

Abstract: In the last few years, the NICER collaboration has provided mass and radius inferences, via pulse profile modeling, for two pulsars: PSR J0030+0451 and PSR J0740+6620. Given the importance of these results for constraining the equation of state of dense nuclear matter, it is crucial to validate them and test their robustness. We therefore explore the reliability of these results and their sensitiv… ▽ More In the last few years, the NICER collaboration has provided mass and radius inferences, via pulse profile modeling, for two pulsars: PSR J0030+0451 and PSR J0740+6620. Given the importance of these results for constraining the equation of state of dense nuclear matter, it is crucial to validate them and test their robustness. We therefore explore the reliability of these results and their sensitivity to analysis settings and random processes, including noise, focusing on the specific case of PSR J0030+0451. We use X-PSI, one of the two main analysis pipelines currently employed by the NICER collaboration for mass and radius inferences. With synthetic data that mimic the PSR J0030+0451 NICER data set, we evaluate the recovery performances of X-PSI under conditions never tested before, including complex modeling of the thermally emitting neutron star surface. For the test cases explored, our results suggest that X-PSI is capable of recovering the true mass and radius within reasonable credible intervals. This work also reveals the main vulnerabilities of the analysis: a significant dependence on noise and the presence of multi-modal structure in the posterior surface. Noise particularly impacts our sensitivity to the analysis settings and widths of the posterior distributions. The multi-modal structure in the posterior suggests that biases could be present if the analysis is unable to exhaustively explore the parameter space. Convergence testing, to ensure an adequate coverage of the parameter space and a suitable representation of the posterior distribution, is one possible solution to these challenges. △ Less

Submitted 16 August, 2023; originally announced August 2023.

Comments: 27 pages, 13 figures

Journal ref: ApJ, 959, 55 (2023)

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

Constraining a relativistic mean field model using neutron star mass-radius measurements I: Nucleonic models

Authors: Chun Huang, Geert Raaijmakers, Anna L. Watts, Laura Tolos, Constança Providência

Abstract: Measurements of neutron star mass and radius or tidal deformability deliver unique insight into the equation of state (EOS) of cold dense matter. EOS inference is very often done using generalized parametric or non-parametric models which deliver no information on composition. In this paper we consider a microscopic nuclear EOS model based on a field theoretical approach. We show that current meas… ▽ More Measurements of neutron star mass and radius or tidal deformability deliver unique insight into the equation of state (EOS) of cold dense matter. EOS inference is very often done using generalized parametric or non-parametric models which deliver no information on composition. In this paper we consider a microscopic nuclear EOS model based on a field theoretical approach. We show that current measurements from NICER and gravitational wave observations constrain primarily the symmetric nuclear matter EOS. We then explore what could be delivered by measurements of mass and radius at the level anticipated for future large-area X-ray timing telescopes. These should be able to place very strong limits on the symmetric nuclear matter EOS, in addition to constraining the nuclear symmetry energy that determines the proton fraction inside the neutron star. △ Less

Submitted 22 March, 2024; v1 submitted 30 March, 2023; originally announced March 2023.

Comments: Accepted for publication in MNRAS

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

X-ray burst ignition location on the surface of accreting X-ray pulsars: Can bursts preferentially ignite at the hotspot?

Authors: A. J. Goodwin, A. Heger, F. R. N. Chambers, A. L. Watts, Y. Cavecchi

Abstract: Hotspots on the surface of accreting neutron stars have been directly observed via pulsations in the lightcurves of X-ray pulsars. They are thought to occur due to magnetic channelling of the accreted fuel to the neutron star magnetic poles. Some X-ray pulsars exhibit burst oscillations during Type I thermonuclear X-ray bursts which are thought to be caused by asymmetries in the burning. In rapidl… ▽ More Hotspots on the surface of accreting neutron stars have been directly observed via pulsations in the lightcurves of X-ray pulsars. They are thought to occur due to magnetic channelling of the accreted fuel to the neutron star magnetic poles. Some X-ray pulsars exhibit burst oscillations during Type I thermonuclear X-ray bursts which are thought to be caused by asymmetries in the burning. In rapidly rotating neutron stars, it has been shown that the lower gravity at the equator can lead to preferential ignition of X-ray bursts at this location. These models, however, do not include the effect of accretion hotspots at the neutron star surface. There are two accreting neutron star sources in which burst oscillations have been observed to track exactly the neutron star spin period. We analyse whether this could be due to the X-ray bursts igniting at the magnetic pole of the neutron star, because of heating in the accreted layers under the hotspot causing ignition conditions to be reached earlier. We investigate heat transport in the accreted layers using a 2D model and study the prevalence of heating down to the ignition depth of X-ray bursts for different hotspot temperatures and sizes. We perform calculations for accretion at the pole and at the equator, and infer that ignition could occur away from the equator at the magnetic pole for hotspots with temperatures greater than $1\times10^8$ K. However, current observations have not identified such high temperatures in accreting X-ray pulsars. △ Less

Submitted 8 June, 2021; originally announced June 2021.

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

Constraints on the dense matter equation of state and neutron star properties from NICER's mass-radius estimate of PSR J0740+6620 and multimessenger observations

Authors: G. Raaijmakers, S. K. Greif, K. Hebeler, T. Hinderer, S. Nissanke, A. Schwenk, T. E. Riley, A. L. Watts, J. M. Lattimer, W. C. G. Ho

Abstract: In recent years our understanding of the dense matter equation of state (EOS) of neutron stars has significantly improved by analyzing multimessenger data from radio/X-ray pulsars, gravitational wave events, and from nuclear physics constraints. Here we study the additional impact on the EOS from the jointly estimated mass and radius of PSR J0740+6620, presented in Riley et al. (2021) by analyzing… ▽ More In recent years our understanding of the dense matter equation of state (EOS) of neutron stars has significantly improved by analyzing multimessenger data from radio/X-ray pulsars, gravitational wave events, and from nuclear physics constraints. Here we study the additional impact on the EOS from the jointly estimated mass and radius of PSR J0740+6620, presented in Riley et al. (2021) by analyzing a combined dataset from X-ray telescopes NICER and XMM-Newton. We employ two different high-density EOS parameterizations: a piecewise-polytropic (PP) model and a model based on the speed of sound in a neutron star (CS). At nuclear densities these are connected to microscopic calculations of neutron matter based on chiral effective field theory interactions. In addition to the new NICER data for this heavy neutron star, we separately study constraints from the radio timing mass measurement of PSR J0740+6620, the gravitational wave events of binary neutron stars GW190425 and GW170817, and for the latter the associated kilonova AT2017gfo. By combining all these, and the NICER mass-radius estimate of PSR J0030+0451 we find the radius of a 1.4 solar mass neutron star to be constrained to the 95% credible ranges 12.33^{+0.76}_{-0.81} km (PP model) and 12.18^{+0.56}_{-0.79} km (CS model). In addition, we explore different chiral effective field theory calculations and show that the new NICER results provide tight constraints for the pressure of neutron star matter at around twice saturation density, which shows the power of these observations to constrain dense matter interactions at intermediate densities. △ Less

Submitted 13 June, 2021; v1 submitted 14 May, 2021; originally announced May 2021.

Comments: 17 pages, 8 figures; accepted for publication in the Astrophysical Journal Letters

Journal ref: The Astrophysical Journal Letters, Volume 918 (2021), L29

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

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

Search for Continuous Gravitational Waves from Scorpius X-1 in LIGO O2 Data

Authors: Yuanhao Zhang, Maria Alessandra Papa, Badri Krishnan, Anna L. Watts

Abstract: We present the results of a search in LIGO O2 public data for continuous gravitational waves from the neutron star in the low-mass X-ray binary Scorpius X-1. We search for signals with $\approx$ constant frequency in the range 40-180 Hz. Thanks to the efficiency of our search pipeline we can use a long coherence time and achieve unprecedented sensitivity, significantly improving on existing result… ▽ More We present the results of a search in LIGO O2 public data for continuous gravitational waves from the neutron star in the low-mass X-ray binary Scorpius X-1. We search for signals with $\approx$ constant frequency in the range 40-180 Hz. Thanks to the efficiency of our search pipeline we can use a long coherence time and achieve unprecedented sensitivity, significantly improving on existing results. This is the first search that has been able to probe gravitational wave amplitudes that could balance the accretion torque at the neutron star radius. Our search excludes emission at this level between 67.5 Hz and 131.5 Hz, for an inclination angle $44^\circ \pm 6^\circ$ derived from radio observations (Fomalont et al. 2001), and assuming that the spin axis is perpendicular to the orbital plane. If the torque arm is $\approx $ 26 km -- a conservative estimate of the \alfven\ radius -- our results are more constraining than the indirect limit across the band. This allows us to exclude certain mass-radius combinations and to place upper limits on the strength of the star's magnetic field. We also correct a mistake that appears in the literature in the equation that gives the gravitational wave amplitude at the torque balance (Abbott et al. 2017b, 2019a) and we re-interpret the associated latest LIGO/Virgo results in light of this. △ Less

Submitted 9 November, 2020; originally announced November 2020.

Comments: 12 pages, 4 figures

Deep model simulation of polar vortices in gas giant atmospheres

Authors: Ferran Garcia, Frank R. N. Chambers, Anna L. Watts

Abstract: The Cassini and Juno probes have revealed large coherent cyclonic vortices in the polar regions of Saturn and Jupiter, a dramatic contrast from the east-west banded jet structure seen at lower latitudes. Debate has centered on whether the jets are shallow, or extend to greater depths in the planetary envelope. Recent experiments and observations have demonstrated the relevance of deep convection m… ▽ More The Cassini and Juno probes have revealed large coherent cyclonic vortices in the polar regions of Saturn and Jupiter, a dramatic contrast from the east-west banded jet structure seen at lower latitudes. Debate has centered on whether the jets are shallow, or extend to greater depths in the planetary envelope. Recent experiments and observations have demonstrated the relevance of deep convection models to a successful explanation of jet structure and cyclonic coherent vortices away from the polar regions have been simulated recently including an additional stratified shallow layer. Here we present new convective models able to produce long-lived polar vortices. Using simulation parameters relevant for giant planet atmospheres we find flow regimes that are in agreement with geostrophic turbulence (GT) theory in rotating convection for the formation of large scale coherent structures via an upscale energy transfer fully three-dimensional. Our simulations generate polar characteristics qualitatively similar to those seen by Juno and Cassini: they match the structure of cyclonic vortices seen on Jupiter; or can account for the existence of a strong polar vortex extending downwards to lower latitudes with a marked spiral morphology and the hexagonal pattern seen on Saturn. Our findings indicate that these vortices can be generated deep in the planetary interior. A transition differentiating these two polar flows regimes is described, interpreted in terms of different force balances and compared with previous shallow atmospheric models which characterised polar vortex dynamics in giant planets. In addition, the heat transport properties are investigated confirming recent scaling laws obtained in the context of reduced models of GT. △ Less

Submitted 23 September, 2020; v1 submitted 26 June, 2020; originally announced June 2020.

Comments: 18 pages, 13 figures and 3 tables

MSC Class: 76-10; 76F35; 86-10

Waves in Thin Oceans on Oblate Neutron Stars

Authors: Bart F. A. van Baal, Frank R. N. Chambers, Anna L. Watts

Abstract: Waves in thin fluid layers are important in various stellar and planetary problems. Due to rapid rotation such systems will become oblate, with a latitudinal variation in the gravitational acceleration across the surface of the object. In the case of accreting neutron stars, rapid rotation could lead to a polar radius smaller than the equatorial radius by a factor $\sim 0.8$. We investigate how th… ▽ More Waves in thin fluid layers are important in various stellar and planetary problems. Due to rapid rotation such systems will become oblate, with a latitudinal variation in the gravitational acceleration across the surface of the object. In the case of accreting neutron stars, rapid rotation could lead to a polar radius smaller than the equatorial radius by a factor $\sim 0.8$. We investigate how the oblateness and a changing gravitational acceleration affect different hydrodynamic modes that exist in such fluid layers through analytic approximations and numerical calculations. The wave vectors of $g$-modes and Yanai modes increase for more oblate systems compared to spherical counterparts, although the impact of variations in the changing gravitational acceleration is effectively negligible. We find that for increased oblateness, Kelvin modes show less equatorial confinement and little change in their wave vector. For $r$-modes, we find that for more oblate systems the wave vector decreases. The exact manner of these changes for the $r$-modes depends on the model for the gravitational acceleration across the surface. △ Less

Submitted 11 June, 2020; originally announced June 2020.

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

The Multi-INstrument Burst ARchive (MINBAR)

Authors: D. K. Galloway, J. J. M. in 't Zand, J. Chenevez, H. Wörpel, L. Keek, L. Ootes, A. L. Watts, L. Gisler, C. Sanchez-Fernandez, E. Kuulkers

Abstract: We present the largest sample of type-I (thermonuclear) X-ray bursts yet assembled, comprising 7083 bursts from 85 bursting sources. The sample is drawn from observations with Xenon-filled proportional counters on the long-duration satellites RXTE, BeppoSAX and INTEGRAL, between 1996 February 8, and 2012 May 3. The burst sources were drawn from a comprehensive catalog of 115 burst sources, assembl… ▽ More We present the largest sample of type-I (thermonuclear) X-ray bursts yet assembled, comprising 7083 bursts from 85 bursting sources. The sample is drawn from observations with Xenon-filled proportional counters on the long-duration satellites RXTE, BeppoSAX and INTEGRAL, between 1996 February 8, and 2012 May 3. The burst sources were drawn from a comprehensive catalog of 115 burst sources, assembled from earlier catalogs and the literature. We carried out a consistent analysis for each burst lightcurve (normalised to the relative instrumental effective area), and provide measurements of rise time, peak intensity, burst timescale, and fluence. For bursts observed with the RXTE/PCA and BeppoSAX/WFC we also provide time-resolved spectroscopy, including estimates of bolometric peak flux and fluence, and spectral parameters at the peak of the burst. For 950 bursts observed with the PCA from sources with previously detected burst oscillations, we include an analysis of the high-time resolution data, providing information on the detectability and amplitude of the oscillations, as well as where in the burst they are found. We also present analysis of 118848 observations of the burst sources within the sample timeframe. We extracted 3-25 keV X-ray spectra from most observations, and (for observations meeting our signal-to-noise criterion), we provide measurements of the flux, spectral colours, and for selected sources, the position on the colour-colour diagram, for the best-fit spectral model. We present a description of the sample, a summary of the science investigations completed to date, and suggestions for further studies. △ Less

Submitted 12 August, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: 105 pages, 34 figures, 18 tables; machine-readable sample data tables available via https://doi.org/10.26180/5e4a697d9b8b6; web interface to sample data available via http://burst.sci.monash.edu. Updated observation numbers to reflect final data table versions, slight modifications to text to preserve table ordering. Now published by ApJS

Journal ref: The Astrophysical Journal Supplement Series, Volume 249, Number 2, 2020

Constraining the dense matter equation of state with joint analysis of NICER and LIGO/Virgo measurements

Authors: G. Raaijmakers, S. K. Greif, T. E. Riley, T. Hinderer, K. Hebeler, A. Schwenk, A. L. Watts, S. Nissanke, S. Guillot, J. M. Lattimer, R. M. Ludlam

Abstract: The NICER collaboration recently published a joint estimate of the mass and the radius of PSR J0030+0451, derived via X-ray pulse-profile modeling. Raaijmakers et al. (2019) explored the implications of this measurement for the dense matter equation of state (EOS) using two parameterizations of the high-density EOS: a piecewise-polytropic model, and a model based on the speed of sound in neutron s… ▽ More The NICER collaboration recently published a joint estimate of the mass and the radius of PSR J0030+0451, derived via X-ray pulse-profile modeling. Raaijmakers et al. (2019) explored the implications of this measurement for the dense matter equation of state (EOS) using two parameterizations of the high-density EOS: a piecewise-polytropic model, and a model based on the speed of sound in neutron stars. In this work we obtain further constraints on the EOS following this approach, but we also include information about the tidal deformability of neutron stars from the gravitational wave signal of the compact binary merger GW170817. We compare the constraints on the EOS to those set by the recent measurement of a 2.14 solar mass pulsar, included as a likelihood function approximated by a Gaussian, and find a small increase in information gain. To show the flexibility of our method, we also explore the possibility that GW170817 was a neutron star-black hole merger, which yields weaker constraints on the EOS. △ Less

Submitted 17 April, 2020; v1 submitted 23 December, 2019; originally announced December 2019.

Comments: 18 pages, 8 figures. Accepted for publication in ApJ Letters

Journal ref: The Astrophysical Journal Letters, 893, L21 (2020)

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

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)

Relativistic ocean $r$-modes during type-I X-ray bursts

Authors: Frank R. N. Chambers, Anna L. Watts

Abstract: Accreting neutron stars (NS) can exhibit high frequency modulations in their lightcurves during thermonuclear X-ray bursts, known as burst oscillations. These frequencies can be offset from the NS spin frequency by several Hz (where known independently) and can drift by 1-3 Hz. One plausible explanation is that a wave is present in the bursting ocean, the rotating frame frequency of which is the o… ▽ More Accreting neutron stars (NS) can exhibit high frequency modulations in their lightcurves during thermonuclear X-ray bursts, known as burst oscillations. These frequencies can be offset from the NS spin frequency by several Hz (where known independently) and can drift by 1-3 Hz. One plausible explanation is that a wave is present in the bursting ocean, the rotating frame frequency of which is the offset. The frequency of the wave should decrease (in the rotating frame) as the burst cools hence explaining the drift. A strong candidate is a buoyant $r$-mode. To date, models that calculated the frequency of this mode taking into account the radial structure neglected relativistic effects and predicted rotating frame frequencies of $\sim$ 4 Hz and frequency drifts of > 5 Hz; too large to be consistent with observations. We present a calculation that includes frame-dragging and gravitational redshift that reduces the rotating frame frequency by up to 30 % and frequency drift by up to 20 %. Updating previous models for the ocean cooling in the aftermath of the burst to a model more representative of detailed calculations of thermonuclear X-ray bursts reduces the frequency of the mode still further. This model, combined with relativistic effects, can reduce the rotating frequency of the mode to $\sim$ 2 Hz and frequency drift to $\sim$ 2 Hz, which is closer to the observed values. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Constraining the neutron star equation of state using Pulse Profile Modeling

Authors: Anna L. Watts

Abstract: One very promising technique for measuring the dense matter Equation of State exploits hotspots that form on the neutron star surface due to the pulsar mechanism, accretion streams, or during thermonuclear explosions in the neutron star ocean. This article explains how Pulse Profile Modeling of hotspots is being used by the Neutron Star Interior Composition Explorer (NICER), an X-ray telescope ins… ▽ More One very promising technique for measuring the dense matter Equation of State exploits hotspots that form on the neutron star surface due to the pulsar mechanism, accretion streams, or during thermonuclear explosions in the neutron star ocean. This article explains how Pulse Profile Modeling of hotspots is being used by the Neutron Star Interior Composition Explorer (NICER), an X-ray telescope installed on the International Space Station in 2017 - and why the technique is a mission driver for the next, larger-area generation of telescopes including the enhanced X-ray Timing and Polarimetry (eXTP) mission and the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X). △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: To appear in the AIP Conference Proceedings of the Xiamen-CUSTIPEN Workshop on the EOS of Dense Neutron-Rich Matter in the Era of Gravitational Wave Astronomy (January 3 - 7, 2019, Xiamen, China)

A search for the $835\,\text{Hz}$ superburst oscillation signal in the regular thermonuclear bursts of 4U 1636-536

Authors: Emma van der Wateren, Anna L. Watts, Laura S. Ootes

Abstract: Burst oscillations are brightness asymmetries that develop in the burning ocean during thermonuclear bursts on accreting neutron stars. They have been observed during H/He-triggered (Type I) bursts and Carbon-triggered superbursts. The mechanism responsible is not unknown, but the dominant burst oscillation frequency is typically within a few Hz of the spin frequency, where this is independently k… ▽ More Burst oscillations are brightness asymmetries that develop in the burning ocean during thermonuclear bursts on accreting neutron stars. They have been observed during H/He-triggered (Type I) bursts and Carbon-triggered superbursts. The mechanism responsible is not unknown, but the dominant burst oscillation frequency is typically within a few Hz of the spin frequency, where this is independently known. One of the best-studied burst oscillation sources, 4U 1636-536, has oscillations at $581\,\text{Hz}$ in both its regular Type I bursts and in one superburst. Recently however, Strohmayer \& Mahmoodifar reported the discovery of an additional signal at a higher frequency, $835\,\text{Hz}$, during the superburst. This higher frequency is consistent with the predictions for several types of global ocean mode, one of the possible burst oscillation mechanisms. If this is the case then the same physical mechanism may operate in the normal Type I bursts of this source. In this paper we report a stacked search for periodic signals in the regular Type I bursts: we found no significant signal at the higher frequency, with upper limits for the single trial root mean square (rms) fractional amplitude of 0.57(6)\%. Our analysis did however reveal that the dominant $581\,\text{Hz}$ burst oscillation signal is present at a weak level even in the sample of bursts where it cannot be detected in individual bursts. This indicates that any cutoff in the burst oscillation mechanism occurs below the detection threshold of existing X-ray telescopes. △ Less

Submitted 26 March, 2019; v1 submitted 24 March, 2019; originally announced March 2019.

Comments: 6 pages, 2 figures. Accepted for publication by ApJ

Journal ref: The Astrophysical Journal, 876:11 (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

A uniform search for thermonuclear burst oscillations in the RXTE legacy dataset

Authors: Anna V. Bilous, Anna L. Watts

Abstract: We describe a blind uniform search for thermonuclear burst oscillations (TBOs) in the majority of Type-I bursts observed by RXTE (2118 bursts from 57 neutron stars). We examined 2-2002 Hz power spectra from the Fourier transform in sliding 0.5-2 s windows, using fine-binned light curves in 2-60 keV energy range. The significance of the oscillation candidates was assessed by simulations which took… ▽ More We describe a blind uniform search for thermonuclear burst oscillations (TBOs) in the majority of Type-I bursts observed by RXTE (2118 bursts from 57 neutron stars). We examined 2-2002 Hz power spectra from the Fourier transform in sliding 0.5-2 s windows, using fine-binned light curves in 2-60 keV energy range. The significance of the oscillation candidates was assessed by simulations which took into account light curve variations, dead time and sliding time windows. Some of our sources exhibited multi-frequency variability below approximately 15 Hz that cannot be readily removed with light-curve modeling and may have an astrophysical (non-TBO) nature. Overall, we found that the number and strength of potential candidates depends strongly on the parameters of the search. We found candidates from all previously known RXTE TBO sources, with pulsations that had been detected at similar frequencies in multiple independent time windows, and discovered TBOs from SAX J1810.8-2658. We could not confirm most previously-reported tentative TBO detections or identify any obvious candidates just below the detection threshold at similar frequencies in multiple bursts. We computed fractional amplitudes of all TBO candidates and placed upper limits on non-detections. Finally, for a few sources we noted small excess of candidates with powers comparable to fainter TBOs, but appearing in single independent time bins at random frequencies. At least some of these candidates may be noise spikes that appear interesting due to selection effects. The potential presence of such candidates calls for extra caution if claiming single-bin TBO detections. △ Less

Submitted 10 June, 2019; v1 submitted 27 December, 2018; originally announced December 2018.

Comments: Replaced to match the accepted version. Accepted to ApJS after numerous minor corrections

Equation of state sensitivities when inferring neutron star and dense matter properties

Authors: S. K. Greif, G. Raaijmakers, K. Hebeler, A. Schwenk, A. L. Watts

Abstract: Understanding the dense matter equation of state at extreme conditions is an important open problem. Astrophysical observations of neutron stars promise to solve this, with NICER poised to make precision measurements of mass and radius for several stars using the waveform modelling technique. What has been less clear, however, is how these mass-radius measurements might translate into equation of… ▽ More Understanding the dense matter equation of state at extreme conditions is an important open problem. Astrophysical observations of neutron stars promise to solve this, with NICER poised to make precision measurements of mass and radius for several stars using the waveform modelling technique. What has been less clear, however, is how these mass-radius measurements might translate into equation of state constraints and what are the associated equation of state sensitivities. We use Bayesian inference to explore and contrast the constraints that would result from different choices for the equation of state parametrization; comparing the well-established piecewise polytropic parametrization to one based on physically motivated assumptions for the speed of sound in dense matter. We also compare the constraints resulting from Bayesian inference to those from simple compatibility cuts. We find that the choice of equation of state parametrization and particularly its prior assumptions can have a significant effect on the inferred global mass-radius relation and the equation of state constraints. Our results point to important sensitivities when inferring neutron star and dense matter properties. This applies also to inferences from gravitational wave observations. △ Less

Submitted 6 March, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: to appear in MNRAS

Observatory science with eXTP

Authors: Jean J. M. in 't Zand, Enrico Bozzo, Jinlu Qu, Xiang-Dong Li, Lorenzo Amati, Yang Chen, Immacolata Donnarumma, Victor Doroshenko, Stephen A. Drake, Margarita Hernanz, Peter A. Jenke, Thomas J. Maccarone, Simin Mahmoodifar, Domitilla de Martino, Alessandra De Rosa, Elena M. Rossi, Antonia Rowlinson, Gloria Sala, Giulia Stratta, Thomas M. Tauris, Joern Wilms, Xuefeng Wu, Ping Zhou, Iván Agudo, Diego Altamirano , et al. (159 additional authors not shown)

Abstract: In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to stu… ▽ More In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed 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)

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)

The enhanced X-ray Timing and Polarimetry mission - eXTP

Authors: ShuangNan Zhang, Andrea Santangelo, Marco Feroci, YuPeng Xu, FangJun Lu, Yong Chen, Hua Feng, Shu Zhang, Søren Brandt, Margarita Hernanz, Luca Baldini, Enrico Bozzo, Riccardo Campana, Alessandra De Rosa, YongWei Dong, Yuri Evangelista, Vladimir Karas, Norbert Meidinger, Aline Meuris, Kirpal Nandra, Teng Pan, Giovanni Pareschi, Piotr Orleanski, QiuShi Huang, Stephane Schanne , et al. (125 additional authors not shown)

Abstract: In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In ad… ▽ More In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment. △ Less

Submitted 10 December, 2018; originally announced December 2018.

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

Burning in the tail: implications for a burst oscillation model

Authors: Frank R. N. Chambers, Anna L. Watts, L. Keek, Yuri Cavecchi, F. Garcia

Abstract: Accreting neutron stars (NS) can exhibit high-frequency modulations, known as burst oscillations, in their lightcurves during thermonuclear X-ray bursts. Their frequencies can be offset from the spin frequency of the NS (known independently) by several Hz, and can drift by 1-3 Hz. One plausible explanation for this phenomenon is that a wave is present in the bursting ocean that decreases in freque… ▽ More Accreting neutron stars (NS) can exhibit high-frequency modulations, known as burst oscillations, in their lightcurves during thermonuclear X-ray bursts. Their frequencies can be offset from the spin frequency of the NS (known independently) by several Hz, and can drift by 1-3 Hz. One plausible explanation for this phenomenon is that a wave is present in the bursting ocean that decreases in frequency (in the rotating frame) as the burst cools. The strongest candidate is the buoyant $r$-mode; however, models for the burning ocean background used in previous studies over-predict frequency drifts by several Hz. Using new background models (which include shallow heating, and burning in the tail of the burst) the evolution of the buoyant $r$-mode is calculated. The resulting frequency drifts are smaller, in line with observations. This illustrates the importance of accounting for the detailed nuclear physics in these bursts. △ Less

Submitted 29 November, 2018; originally announced November 2018.

Comments: 7 pages, 2 figures

Thermal convection in rotating spherical shells: temperature-dependent internal heat generation using the example of triple-$α$ burning in neutron stars

Authors: F. Garcia, F. R. N Chambers, A. L. Watts

Abstract: We present an extensive study of Boussinesq thermal convection including a temperature-dependent internal heating source, based on numerical three-dimensional simulations. The temperature dependence mimics triple-$α$ nuclear reactions and the fluid geometry is a rotating spherical shell. These are key ingredients for the study of convective accreting neutron star oceans. A dimensionless parameter… ▽ More We present an extensive study of Boussinesq thermal convection including a temperature-dependent internal heating source, based on numerical three-dimensional simulations. The temperature dependence mimics triple-$α$ nuclear reactions and the fluid geometry is a rotating spherical shell. These are key ingredients for the study of convective accreting neutron star oceans. A dimensionless parameter ${\rm Ra}_n$, measuring the relevance of nuclear heating, is defined. We explore how flow characteristics change with increasing ${\rm Ra}_n$ and give an astrophysical motivation. The onset of convection is investigated with respect to this parameter and periodic, quasiperiodic, chaotic flows with coherent structures, and fully turbulent flows are exhibited as ${\rm Ra}_n$ is varied. Several regime transitions are identified and compared with previous results on differentially heated convection. Finally, we explore (tentatively) the potential applicability of our results to the evolution of thermonuclear bursts in accreting neutron star oceans. △ Less

Submitted 14 December, 2018; v1 submitted 13 July, 2018; originally announced July 2018.

Comments: 9 Figures, 6 Tables. Published in the Physical Review Fluids journal

MSC Class: 76E20; 76F06; 76U05

Journal ref: Phys. Rev. Fluids 3, 123501 (2018)

STROBE-X: A probe-class mission for X-ray spectroscopy and timing on timescales from microseconds to years

Authors: Paul S. Ray, Zaven Arzoumanian, Søren Brandt, Eric Burns, Deepto Chakrabarty, Marco Feroci, Keith C. Gendreau, Olivier Gevin, Margarita Hernanz, Peter Jenke, Steven Kenyon, José Luis Gálvez Thomas J. Maccarone, Takashi Okajima, Ronald A. Remillard, Stéphane Schanne, Chris Tenzer, Andrea Vacchi, Colleen A. Wilson-Hodge, Berend Winter, Silvia Zane, David R. Ballantyne, Enrico Bozzo, Laura W. Brenneman, Edward Cackett, Alessandra De Rosa , et al. (8 additional authors not shown)

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

Submitted 3 July, 2018; originally announced July 2018.

Comments: 20 pages, 10 figures, to appear in Proc. SPIE Astronomical Telescopes and Instrumentation 2018

Polar waves and chaotic flows in thin rotating spherical shells

Authors: F. Garcia, F. R. N. Chambers, A. L. Watts

Abstract: Convection in rotating spherical geometries is an important physical process in planetary and stellar systems. Using continuation methods at low Prandtl number, we find both strong equatorially asymmetric and symmetric polar nonlinear rotating waves in a model of thermal convection in thin rotating spherical shells with stress-free boundary conditions. For the symmetric waves convection is confine… ▽ More Convection in rotating spherical geometries is an important physical process in planetary and stellar systems. Using continuation methods at low Prandtl number, we find both strong equatorially asymmetric and symmetric polar nonlinear rotating waves in a model of thermal convection in thin rotating spherical shells with stress-free boundary conditions. For the symmetric waves convection is confined to high latitude in both hemispheres but is only restricted to one hemisphere close to the pole in the case of asymmetric waves. This is in contrast to what is previously known from studies in the field. These periodic flows, in which the pattern is rotating steadily in the azimuthal direction, develop a strong axisymmetric component very close to onset. Using stability analysis of periodic orbits the regions of stability are determined and the topology of the stable/unstable oscillatory flows bifurcated from the branches of rotating waves is described. By means of direct numerical simulations of these oscillatory chaotic flows, we show that these three-dimensional convective polar flows exhibit characteristics, such as force balance or mean physical properties, which are similar to flows occuring in planetary atmospheres. We show that these results may open a route to understanding unexplained features of gas giant atmospheres, in particular for the case of Jupiter. These include the observed equatorial asymmetry with a pronounced decrease at the equator (the so-called dimple), and the coherent vortices surrounding the poles recently observed by the Juno mission. △ Less

Submitted 6 July, 2019; v1 submitted 27 June, 2018; originally announced June 2018.

Comments: Published in Physical Review Fluids (2019). Contains 2 tables and 8 figures

MSC Class: 76E20 37G35 37G40 65P30 65P40

Journal ref: Phys. Rev. Fluids 4, 074802 (2019)

Millisecond oscillations in the bursting flux of SAX J1810.8-2609

Authors: Anna V Bilous, Anna L Watts, Duncan K Galloway, Jean J M in 't Zand

Abstract: SAX J1810.8-2609 is a faint X-ray transient, mostly known for its abnormally low quiescent thermal luminosity, which disagrees with standard cooling models. It is also one of a small sample of stars whose mass and radius have been estimated using spectral modeling of one of its thermonuclear bursts. Here we report the discovery of millisecond oscillations in a type I thermonuclear X-ray burst from… ▽ More SAX J1810.8-2609 is a faint X-ray transient, mostly known for its abnormally low quiescent thermal luminosity, which disagrees with standard cooling models. It is also one of a small sample of stars whose mass and radius have been estimated using spectral modeling of one of its thermonuclear bursts. Here we report the discovery of millisecond oscillations in a type I thermonuclear X-ray burst from SAX J1810.8-2609 observed by RXTE during the 2007 outburst. A strong signal (Leahy-normalized power of 71.5, 4.5e-9 chance of coincidence with a conservative estimate for the number of trials) was present at 531.8 Hz during the decay of one out of six bursts observed. Oscillations were detected for about 6 seconds, during which their frequency increased from 531.4 to 531.9 Hz in a manner similar to other burst oscillation sources. The millisecond oscillations discovered pinpoint the spin frequency of the neutron star, which is important for the spectral modeling, associated mass-radius inference, and understanding the evolutionary status and cooling behavior of the star. As of April 2018 the source is in outburst again, providing a fleeting opportunity to acquire new material for the burst oscillation searches. △ Less

Submitted 25 May, 2018; originally announced May 2018.

Comments: 4 pages, 2 figures, submitted to ApJ letters

A pitfall of piecewise-polytropic equation of state inference

Authors: G. Raaijmakers, T. E. Riley, A. L. Watts

Abstract: The only messenger radiation in the Universe which one can use to statistically probe the Equation of State (EOS) of cold dense matter is that originating from the near-field vicinities of compact stars. Constraining gravitational masses and equatorial radii of rotating compact stars is a major goal for current and future telescope missions, with a primary purpose of constraining the EOS. From a B… ▽ More The only messenger radiation in the Universe which one can use to statistically probe the Equation of State (EOS) of cold dense matter is that originating from the near-field vicinities of compact stars. Constraining gravitational masses and equatorial radii of rotating compact stars is a major goal for current and future telescope missions, with a primary purpose of constraining the EOS. From a Bayesian perspective it is necessary to carefully discuss prior definition; in this context a complicating issue is that in practice there exist pathologies in the general relativistic mapping between spaces of local (interior source matter) and global (exterior spacetime) parameters. In a companion paper, these issues were raised on a theoretical basis. In this study we reproduce a probability transformation procedure from the literature in order to map a joint posterior distribution of Schwarzschild gravitational masses and radii into a joint posterior distribution of EOS parameters. We demonstrate computationally that EOS parameter inferences are sensitive to the choice to define a prior on a joint space of these masses and radii, instead of on a joint space interior source matter parameters. We focus on the piecewise-polytropic EOS model, which is currently standard in the field of astrophysical dense matter study. We discuss the implications of this issue for the field. △ Less

Submitted 24 April, 2018; originally announced April 2018.

Comments: 16 pages, 9 figures. Accepted for publication in MNRAS

On parametrised cold dense matter equation of state inference

Authors: Thomas E. Riley, Geert Raaijmakers, Anna L. Watts

Abstract: Constraining the equation of state of cold dense matter in compact stars is a major science goal for observing programmes being conducted using X-ray, radio, and gravitational wave telescopes. We discuss Bayesian hierarchical inference of parametrised dense matter equations of state. In particular we generalise and examine two inference paradigms from the literature: (i) direct posterior equation… ▽ More Constraining the equation of state of cold dense matter in compact stars is a major science goal for observing programmes being conducted using X-ray, radio, and gravitational wave telescopes. We discuss Bayesian hierarchical inference of parametrised dense matter equations of state. In particular we generalise and examine two inference paradigms from the literature: (i) direct posterior equation of state parameter estimation, conditioned on observations of a set of rotating compact stars; and (ii) indirect parameter estimation, via transformation of an intermediary joint posterior distribution of exterior spacetime parameters (such as gravitational masses and coordinate equatorial radii). We conclude that the former paradigm is not only tractable for large-scale analyses, but is principled and flexible from a Bayesian perspective whilst the latter paradigm is not. The thematic problem of Bayesian prior definition emerges as the crux of the difference between these paradigms. The second paradigm should in general only be considered as an ill-defined approach to the problem of utilising archival posterior constraints on exterior spacetime parameters; we advocate for an alternative approach whereby such information is repurposed as an approximative likelihood function. We also discuss why conditioning on a piecewise-polytropic equation of state model - currently standard in the field of dense matter study - can easily violate conditions required for transformation of a probability density distribution between spaces of exterior (spacetime) and interior (source matter) parameters. △ Less

Submitted 24 April, 2018; originally announced April 2018.

Comments: Accepted for publication in MNRAS, 43 pages, 7 figures

Superburst oscillations: ocean and crustal modes excited by Carbon-triggered Type I X-ray bursts

Authors: Frank R. N. Chambers, Anna L. Watts, Yuri Cavecchi, F. Garcia, L. Keek

Abstract: Accreting neutron stars (NS) can exhibit high frequency modulations in their lightcurves during thermonuclear X-ray bursts, known as burst oscillations. The frequencies can be offset from the spin frequency of the NS by several Hz, and can drift by 1-3 Hz. One possible explanation is a mode in the bursting ocean, the frequency of which would decrease (in the rotating frame) as the burst cools, hen… ▽ More Accreting neutron stars (NS) can exhibit high frequency modulations in their lightcurves during thermonuclear X-ray bursts, known as burst oscillations. The frequencies can be offset from the spin frequency of the NS by several Hz, and can drift by 1-3 Hz. One possible explanation is a mode in the bursting ocean, the frequency of which would decrease (in the rotating frame) as the burst cools, hence explaining the drifts. Most burst oscillations have been observed during H/He triggered bursts, however there has been one observation of oscillations during a superburst; hours' long Type I X-ray bursts caused by unstable carbon burning deeper in the ocean. This paper calculates the frequency evolution of an oceanic r-mode during a superburst. The rotating frame frequency varies during the burst from 4-14 Hz, and is sensitive to the background parameters, in particular the temperature of the ocean and ignition depth. This calculation is compared to the superburst oscillations observed on 4U-1636-536. The predicted mode frequencies ($\sim$ 10 Hz) would require a spin frequency of $\sim$ 592 Hz to match observations; 6 Hz higher than the spin inferred from an oceanic r-mode model for the H/He triggered burst oscillations. This model also over-predicts the frequency drift during the superburst by 90 %. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Comments: Accepted for publication in MNRAS

The onset of low Prandtl number thermal convection in thin spherical shells

Authors: F. Garcia, F. R. N. Chambers, A. L. Watts

Abstract: This study considers the onset of stress-free Boussinesq thermal convection in rotating spherical shells with aspect ratio $η=r_i/r_o=0.9$ ($r_i$ and $r_o$ being the inner and outer radius), Prandtl numbers ${\rm Pr} \in[10^{-4},10^{-1}]$, and Taylor numbers ${\rm Ta}\in[10^{4},10^{12}]$. We are particularly interested in the form of the convective cell pattern that develops, and in its time scale… ▽ More This study considers the onset of stress-free Boussinesq thermal convection in rotating spherical shells with aspect ratio $η=r_i/r_o=0.9$ ($r_i$ and $r_o$ being the inner and outer radius), Prandtl numbers ${\rm Pr} \in[10^{-4},10^{-1}]$, and Taylor numbers ${\rm Ta}\in[10^{4},10^{12}]$. We are particularly interested in the form of the convective cell pattern that develops, and in its time scales, since this may have observational consequences. For a fixed ${\rm Ta}<10^{9}$ and by decreasing ${\rm Pr}$ from 0.1 to $10^{-4}$ a transition between spiralling columnar (SC) and equatorially-attached (EA) modes, and a transition between EA and equatorially antisymmetric or symmetric polar (AP/SP) weakly multicellular modes are found. The latter modes are preferred at very low ${\rm Pr}$. Surprisingly, for ${\rm Ta}>3\times 10^{9}$ the unicellular polar modes become also preferred at moderate ${\rm Pr}\sim10^{-2}$ because two new transition curves between EA and AP/SP and between AP/SP and SC modes are born at a triple-point bifurcation. The dependence on ${\rm Pr}$ and ${\rm Ta}$ of the transitions is studied to estimate the type of modes, and their critical parameters, preferred at different stellar regimes. △ Less

Submitted 15 May, 2018; v1 submitted 8 February, 2018; originally announced February 2018.

Comments: Accepted for publication in Physical Review Fluids. Contains 17 pages, 8 figures and 3 tables. Added brief erratum correcting values used for estimates of neutron star ocean viscosity

MSC Class: 76E20

Rotational phase dependence of magnetar bursts

Authors: C. Elenbaas, A. L. Watts, D. Huppenkothen

Abstract: The trigger for the short bursts observed in $γ$-rays from many magnetar sources remains unknown. One particular open question in this context is the localization of burst emission to a singular active region or a larger area across the neutron star. While several observational studies have attempted to investigate this question by looking at the phase dependence of burst properties, results have… ▽ More The trigger for the short bursts observed in $γ$-rays from many magnetar sources remains unknown. One particular open question in this context is the localization of burst emission to a singular active region or a larger area across the neutron star. While several observational studies have attempted to investigate this question by looking at the phase dependence of burst properties, results have been mixed. At the same time, it is not obvious a priori that bursts from a localized active region would actually give rise to a detectable phase-dependence, taking into account issues such as geometry, relativistic effects, and intrinsic burst properties such brightness and duration. In this paper, we build a simple theoretical model to investigate the circumstances under which the latter effects could affect detectability of a dependence of burst emission on rotational phase. We find that even for strongly phase-dependent emission, inferred burst properties may not show a rotational phase dependence depending on the geometry of the system and the observer. Furthermore, the observed properties of bursts with durations short as 10-20% of the spin period can vary strongly depending on the rotational phase at which the burst was emitted. We also show that detectability of a rotational phase dependence depends strongly on the minimum number of bursts observed, and find that existing burst samples may simply be too small to rule out a phase dependence. △ Less

Submitted 6 February, 2018; originally announced February 2018.

Comments: 15 pages, 17 figures

On the dependence of X-ray burst rate on accretion and spin rate

Authors: Yuri Cavecchi, Anna L. Watts, Duncan K. Galloway

Abstract: Nuclear burning and its dependence on the mass accretion rate are fundamental ingredients for describing the complicated observational phenomenology of neutron stars in binary systems. Motivated by high quality burst rate data emerging from large statistical studies, we report general calculations relating bursting rate to mass accretion rate and neutron star rotation frequency. In this first work… ▽ More Nuclear burning and its dependence on the mass accretion rate are fundamental ingredients for describing the complicated observational phenomenology of neutron stars in binary systems. Motivated by high quality burst rate data emerging from large statistical studies, we report general calculations relating bursting rate to mass accretion rate and neutron star rotation frequency. In this first work we neglect general relativistic effects and accretion topology, though we discuss where their inclusion should play a role. The relations we derive are suitable for different burning regimes and provide a direct link between parameters predicted by theory and what is to be expected in observations. We illustrate this for analytical relations of different unstable burning regimes that operate on the surface of an accreting neutron star. We also use the observed behaviour of burst rate to suggest new constraints on burning parameters. We are able to provide an explanation for the long standing problem of the observed decrease of burst rate with increasing mass accretion that follows naturally from these calculations: when accretion rate crosses a certain threshold, ignition moves away from its initially preferential site and this can cause a net reduction of the burst rate due to the effects of local conditions that set local differences in both burst rate and stabilization criteria. We show under which conditions this can happen even if locally the burst rate keeps increasing with accretion. △ Less

Submitted 12 November, 2017; originally announced November 2017.

Comments: Accepted for publication on ApJ

Thermonuclear burst oscillations and the dense matter equation of state

Authors: Anna L. Watts

Abstract: Matter in neutron star cores reaches extremely high densities, forming states of matter that cannot be generated in the laboratory. The Equation of State (EOS) of the matter links to macroscopic observables, such as mass M and radius R, via the stellar structure equations. A promising technique for measuring M and R exploits hotspots (burst oscillations) that form on the stellar surface when mater… ▽ More Matter in neutron star cores reaches extremely high densities, forming states of matter that cannot be generated in the laboratory. The Equation of State (EOS) of the matter links to macroscopic observables, such as mass M and radius R, via the stellar structure equations. A promising technique for measuring M and R exploits hotspots (burst oscillations) that form on the stellar surface when material accreted from a companion star undergoes a thermonuclear explosion. As the star rotates, the hotspot gives rise to a pulsation, and relativistic effects encode information about M and R into the pulse profile. However the burst oscillation mechanism remains unknown, introducing uncertainty when inferring the EOS. I review the progress that we are making towards cracking this long-standing problem, and establishing burst oscillations as a robust tool for measuring M and R. This is a major goal for future large area X-ray telescopes. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: 4 pages, 2 figures, to appear in the proceedings of "IAUS 337: Pulsar Astrophysics - The Next 50 Years" eds: P. Weltevrede, B.B.P. Perera, L. Levin Preston & S. Sanidas

Magnetar giant flare high-energy emission

Authors: C. Elenbaas, D. Huppenkothen, C. Omand, A. L. Watts, E. Bissaldi, I. Caiazzo, J. Heyl

Abstract: High energy ($> 250$ keV) emission has been detected persisting for several tens of seconds after the initial spike of magnetar giant flares. It has been conjectured that this emission might arise via inverse Compton scattering in a highly extended corona generated by super-Eddington outflows high up in the magnetosphere. In this paper we undertake a detailed examination of this model. We investig… ▽ More High energy ($> 250$ keV) emission has been detected persisting for several tens of seconds after the initial spike of magnetar giant flares. It has been conjectured that this emission might arise via inverse Compton scattering in a highly extended corona generated by super-Eddington outflows high up in the magnetosphere. In this paper we undertake a detailed examination of this model. We investigate the properties of the required scatterers, and whether the mechanism is consistent with the degree of pulsed emission observed in the tail of the giant flare. We conclude that the mechanism is consistent with current data, although the origin of the scattering population remains an open question. We propose an alternative picture in which the emission is closer to that star and is dominated by synchrotron radiation. The $RHESSI$ observations of the December 2004 flare modestly favor this latter picture. We assess the prospects for the Fermi Gamma-Ray Space Telescope to detect and characterize a similar high energy component in a future giant flare. Such a detection should help to resolve some of the outstanding issues. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: 20 pages, 14 figures

Chandra spectroscopy of Rapid Burster type-I X-ray bursts

Authors: J. J. M. in 't Zand, T. Bagnoli, C. D'Angelo, A. Patruno, D. K. Galloway, M. B. M. van der Klis, A. L. Watts, H. L. Marshall

Abstract: We observed the Rapid Burster with Chandra when it was in the 'banana' state that usually precedes the type-II X-ray bursting 'island' state for which the source is particularly known. We employed the High-Energy Transmission Grating Spectrometer in combination with the ACIS-S detector in continuous clocking mode. The observation yielded 20 thermonuclear type-I X-ray bursts emitted from the neutro… ▽ More We observed the Rapid Burster with Chandra when it was in the 'banana' state that usually precedes the type-II X-ray bursting 'island' state for which the source is particularly known. We employed the High-Energy Transmission Grating Spectrometer in combination with the ACIS-S detector in continuous clocking mode. The observation yielded 20 thermonuclear type-I X-ray bursts emitted from the neutron star surface with recurrence times between 0.9 and 1.2 hr, and an e-folding decay time scale of 1 min. We searched for narrow spectral features in the burst emission that could constrain the composition of the ashes of the nuclear burning and the compactness of the neutron star, but found none. The upper limit on the equivalent width of narrow absorption lines between 2 and 6 keV is between 5 and 20 eV (single trial 3 sigma confidence level) and on those of absorption edges between 150 and 400 eV. The latter numbers are comparable to the levels predicted by Weinberg, Bildsten & Schatz (2006) for Eddington-limited thermonuclear bursts. △ Less

Submitted 21 March, 2017; originally announced March 2017.

Comments: Submitted to Proceedings of 11th INTEGRAL conference (Amsterdam, October 2016), to be published in PoS, 10 pages, 5 figures