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

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

MAUVE: A 6 kpc bipolar outflow launched from NGC 4383, one of the most HI-rich galaxies in the Virgo cluster

Authors: Adam B. Watts, Luca Cortese, Barbara Catinella, Amelia Fraser-McKelvie, Eric Emsellem, Lodovico Coccato, Jesse van de Sande, Toby H. Brown, Yago Ascasibar, Andrew Battisti, Alessandro Boselli, Timothy A. Davis, Brent Groves, Sabine Thater

Abstract: Stellar feedback-driven outflows are important regulators of the gas-star formation cycle. However, resolving outflow physics requires high resolution observations that can only be achieved in very nearby galaxies, making suitable targets rare. We present the first results from the new VLT/MUSE large program MAUVE (MUSE and ALMA Unveiling the Virgo Environment), which aims to understand the gas-st… ▽ More Stellar feedback-driven outflows are important regulators of the gas-star formation cycle. However, resolving outflow physics requires high resolution observations that can only be achieved in very nearby galaxies, making suitable targets rare. We present the first results from the new VLT/MUSE large program MAUVE (MUSE and ALMA Unveiling the Virgo Environment), which aims to understand the gas-star formation cycle within the context of the Virgo cluster environment. Outflows are a key part of this cycle, and we focus on the peculiar galaxy NGC 4383, which hosts a $\sim6\,$kpc bipolar outflow fuelled by one of Virgo's most HI-rich discs. The spectacular MUSE data reveal the clumpy structure and complex kinematics of the ionised gas in this M82-like outflow at 100 pc resolution. Using the ionised gas geometry and kinematics we constrain the opening half-angle to $θ=25-35^\circ$, while the average outflow velocity is $\sim210$ kms$^{-1}$. The emission line ratios reveal an ionisation structure where photoionisation is the dominant excitation process. The outflowing gas shows a marginally elevated gas-phase oxygen abundance compared to the disc but is lower than the central starburst, highlighting the contribution of mixing between the ejected and entrained gas. Making some assumptions about the outflow geometry, we estimate an integrated mass outflow-rate of $\sim1.8~$M$_\odot$yr$^{-1}$ and a corresponding mass-loading factor in the range 1.7-2.3. NGC 4383 is a useful addition to the few nearby examples of well-resolved outflows, and will provide a useful baseline for quantifying the role of outflows within the Virgo cluster. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 16 pages, 7 figures, 1 appendix. Accepted to MNRAS

VERTICO and IllustrisTNG: The spatially resolved effects of environment on galactic gas

Authors: Adam R. H. Stevens, Toby Brown, Benedikt Diemer, Annalisa Pillepich, Lars Hernquist, Dylan Nelson, Yannick M. Bahé, Alessandro Boselli, Timothy A. Davis, Pascal J. Elahi, Sara L. Ellison, María J. Jiménez-Donaire, Ian D. Roberts, Kristine Spekkens, Vicente Villanueva, Adam B. Watts, Christine D. Wilson, Nikki Zabel

Abstract: It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies w… ▽ More It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies with $m_* > \! 10^9\,{\rm M}_\odot$ and ${\rm SFR} \! > 0.05\,{\rm M}_\odot\,{\rm yr}^{-1}$. We select galaxies in the simulation that are analogous to those in the HERACLES and VERTICO surveys, and mock-observe them to the common specifications of the data. Although TNG50 does not quantitatively match the observed gas surface densities in the centers of galaxies, the simulation does qualitatively reproduce the trends of gas truncation and central density suppression seen in VERTICO in both HI and H$_2$. This result promises that modern cosmological hydrodynamic simulations can be used to reliably model the post-infall histories of cluster satellite galaxies. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 9 pages, 4 figures, accepted in ApJL

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

VERTICO VII: Environmental quenching caused by suppression of molecular gas content and star formation efficiency in Virgo Cluster galaxies

Authors: Toby Brown, Ian D. Roberts, Mallory Thorp, Sara L. Ellison, Nikki Zabel, Christine D. Wilson, Yannick M. Bahé, Dhruv Bisaria, Alberto D. Bolatto, Alessandro Boselli, Aeree Chung, Luca Cortese, Barbara Catinella, Timothy A. Davis, María J. Jiménez-Donaire, Claudia D. P. Lagos, Bumhyun Lee, Laura C. Parker, Rory Smith, Kristine Spekkens, Adam R. H. Stevens, Vicente Villanueva, Adam B. Watts

Abstract: We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with resp… ▽ More We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with respect to both HI-normal cluster and field galaxies (0.5 dex), suggesting that mechanisms regulating the global HI content are responsible for quenching local star formation. We demonstrate that the observed quenching in HI-poor galaxies is caused by environmental processes such as ram pressure stripping (RPS) simultaneously reducing molecular gas surface density and star formation efficiency (SFE), compared to regions in HI-normal systems (by 0.38 and 0.22 dex, respectively). We observe systematically elevated SFRs that are driven by increased molecular gas surface densities at fixed stellar mass surface density in the outskirts of early-stage RPS galaxies, while SFE remains unchanged with respect to the field sample. We quantify how RPS and starvation affect the star formation cycle of inner and outer galaxy discs as they are processed by the cluster. We show both are effective quenching mechanisms with the key difference being that RPS acts upon the galaxy outskirts while starvation regulates the star formation cycle throughout disc, including within the truncation radius. For both processes, the quenching is caused by a simultaneous reduction in molecular gas surface densities and SFE at fixed stellar mass surface density. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 17 pages, 1 table, 5 figures, accepted for publication in ApJ

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)

VERTICO VI: Cold-gas asymmetries in Virgo cluster galaxies

Authors: Ian D. Roberts, Toby Brown, Nikki Zabel, Christine D. Wilson, Aeree Chung, Laura C. Parker, Dhruv Bisaria, Alessandro Boselli, Barbara Catinella, Ryan Chown, Luca Cortese, Timothy A. Davis, Sara Ellison, Maria Jesus Jimenez-Donaire, Bumhyun Lee, Rory Smith, Kristine Spekkens, Adam R. H. Stevens, Mallory Thorp, Vincente Villanueva, Adam B. Watts, Charlotte Welker, Hyein Yoon

Abstract: We analyze cold-gas distributions in Virgo cluster galaxies using resolved CO(2-1) (tracing molecular hydrogen, H2) and HI observations from the Virgo Environment Traced In CO (VERTICO) and the VLA Imaging of Virgo in Atomic Gas (VIVA) surveys. From a theoretical perspective, it is expected that environmental processes in clusters will have a stronger influence on diffuse atomic gas compared to th… ▽ More We analyze cold-gas distributions in Virgo cluster galaxies using resolved CO(2-1) (tracing molecular hydrogen, H2) and HI observations from the Virgo Environment Traced In CO (VERTICO) and the VLA Imaging of Virgo in Atomic Gas (VIVA) surveys. From a theoretical perspective, it is expected that environmental processes in clusters will have a stronger influence on diffuse atomic gas compared to the relatively dense molecular gas component, and that these environmental perturbations can compress the cold interstellar medium in cluster galaxies leading to elevated star formation. In this work we observationally test these predictions for star-forming satellite galaxies within the Virgo cluster. We divide our Virgo galaxy sample into HI-normal, HI-tailed, and HI-truncated classes and show, unsurprisingly, that the HI-tailed galaxies have the largest quantitative HI asymmetries. We also compare to a control sample of non-cluster galaxies and find that Virgo galaxies, on average, have HI asymmetries that are 40 +/- 10 per cent larger than the control. There is less separation between control, HI-normal, HI-tailed, and HI-truncated galaxies in terms of H2 asymmetries, and on average, Virgo galaxies have H2 asymmetries that are only marginally (20 +/- 10 per cent) larger than the control sample. We find a weak correlation between HI and H2 asymmetries over our entire sample, but a stronger correlation for those specific galaxies being strongly impacted by environmental perturbations. Finally, we divide the discs of the HI-tailed Virgo galaxies into a leading half and trailing half according to the observed tail direction. We find evidence for excess molecular gas mass on the leading halves of the disc. This excess molecular gas on the leading half is accompanied by an excess in star formation rate such that the depletion time is, on average, unchanged. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 15 pages, 8 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 675, A78 (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

VERTICO V: The environmentally driven evolution of the inner cold gas discs of Virgo cluster galaxies

Authors: Adam B. Watts, Luca Cortese, Barbara Catinella, Toby Brown, Christine D. Wilson, Nikki Zabel, Ian D. Roberts, Timothy A. Davis, Mallory Thorp, Aeree Chung, Adam R. H. Stevens, Sara L. Ellison, Kristine Spekkens, Laura C. Parker, Yannick M. Bahé, Vicente Villanueva, María Jiménez-Donaire, Dhruv Bisaria, Alessandro Boselli, Alberto D. Bolatto, Bumhyun Lee

Abstract: The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virg… ▽ More The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virgo Environment traced in CO survey (VERTICO) and archival HI data, we study the HI and molecular gas within the optical discs of Virgo cluster galaxies on 1.2-kpc scales with spatially resolved scaling relations between stellar ($Σ_{\star}$), HI ($Σ_\mathrm{HI}$), and molecular gas ($Σ_\mathrm{mol}$) surface densities. Adopting HI deficiency as a measure of environmental impact, we find evidence that, in addition to removing the HI at large radii, the cluster processes also lower the average $Σ_\mathrm{HI}$ of the remaining gas even in the central 1.2 kpc. The impact on molecular gas is comparatively weaker than on the HI, and we show that the lower $Σ_\mathrm{mol}$ gas is removed first. In the most HI-deficient galaxies, however, we find evidence that environmental processes reduce the typical $Σ_\mathrm{mol}$ of the remaining gas by nearly a factor of 3. We find no evidence for environment-driven elevation of $Σ_\mathrm{HI}$ or $Σ_\mathrm{mol}$ in HI-deficient galaxies. Using the ratio of $Σ_\mathrm{mol}$-to-$Σ_\mathrm{HI}$ in individual regions, we show that changes in the ISM physical conditions, estimated using the total gas surface density and midplane hydrostatic pressure, cannot explain the observed reduction in molecular gas content. Instead, we suggest that direct stripping of the molecular gas is required to explain our results. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: Accepted for publication in PASA. Main text is 19 pages including 12 figures and 3 tables, plus 1 appendix. A 2.5 min, high-level summary can be found at https://youtu.be/7djMmVEpDVc

SAMI-HI: the connection between global asymmetry in the ionised and neutral atomic hydrogen gas in galaxies

Authors: Adam B. Watts, Luca Cortese, Barbara Catinella, Chris Power, Amelia Fraser-McKelvie, Julia J. Bryant, Scott M. Croom, Jesse van de Sande, Joss Bland-Hawthorn, Brent Groves

Abstract: Observations of the neutral atomic hydrogen (HI) gas in galaxies are predominantly spatially unresolved, in the form of a global HI spectral line. There has been substantial work on quantifying asymmetry in global HI spectra (`global HI asymmetry'), but due to being spatially unresolved, it remains unknown what physical regions of galaxies the asymmetry traces, and whether the other gas phases are… ▽ More Observations of the neutral atomic hydrogen (HI) gas in galaxies are predominantly spatially unresolved, in the form of a global HI spectral line. There has been substantial work on quantifying asymmetry in global HI spectra (`global HI asymmetry'), but due to being spatially unresolved, it remains unknown what physical regions of galaxies the asymmetry traces, and whether the other gas phases are affected. Using optical integral field spectrograph (IFS) observations from the Sydney AAO Multi-object IFS (SAMI) survey for which global HI spectra are also available (SAMI-HI), we study the connection between asymmetry in galaxies' ionised and neutral gas reservoirs to test if and how they can help us better understand the origin of global HI asymmetry. We reconstruct the global H$α$ spectral line from the IFS observations and find that, while some global H$α$ asymmetries can arise from disturbed ionised gas kinematics, the majority of asymmetric cases are driven by the distribution of H$α$-emitting gas. When compared to the HI, we find no evidence for a relationship between the global H$α$ and HI asymmetry. Further, a visual inspection reveals that cases where galaxies have qualitatively similar H$α$ and HI spectral profiles can be spurious, with the similarity originating from an irregular 2D H$α$ flux distribution. Our results highlight that comparisons between global H$α$ and HI asymmetry are not straightforward, and that many global HI asymmetries trace disturbances that do not significantly impact the central regions of galaxies. △ Less

Submitted 19 December, 2022; v1 submitted 9 December, 2022; originally announced December 2022.

Comments: 11 pages, 6 figures, 1 appendix, accepted for publication in MNRAS

SAMI-HI: The HI view of the H$α$ Tully-Fisher relation and data release

Authors: Barbara Catinella, Luca Cortese, Alfred L. Tiley, Steven Janowiecki, Adam B. Watts, Julia J. Bryant, Scott M. Croom, Francesco d'Eugenio, Jesse van de Sande, Joss Bland-Hawthorn, Amelia Fraser-McKelvie, Samuel N. Richards, Sarah M. Sweet, Daniel J. Pisano, Nickolas Pingel, Rebecca A. Koopmann, Dillion Cottrill, Meghan Hill

Abstract: We present SAMI-HI, a survey of the atomic hydrogen content of 296 galaxies with integral field spectroscopy available from the SAMI Galaxy Survey. The sample spans nearly 4 dex in stellar mass ($M_\star = 10^{7.4}-10^{11.1}~ \rm M_\odot$), redshift $z<0.06$, and includes new Arecibo observations of 153 galaxies, for which we release catalogues and HI spectra. We use these data to compare the rota… ▽ More We present SAMI-HI, a survey of the atomic hydrogen content of 296 galaxies with integral field spectroscopy available from the SAMI Galaxy Survey. The sample spans nearly 4 dex in stellar mass ($M_\star = 10^{7.4}-10^{11.1}~ \rm M_\odot$), redshift $z<0.06$, and includes new Arecibo observations of 153 galaxies, for which we release catalogues and HI spectra. We use these data to compare the rotational velocities obtained from optical and radio observations and to show how systematic differences affect the slope and scatter of the stellar-mass and baryonic Tully-Fisher relations. Specifically, we show that H$α$ rotational velocities measured in the inner parts of galaxies (1.3 effective radii in this work) systematically underestimate HI global measurements, with HI/H$α$ velocity ratios that increase at low stellar masses, where rotation curves are typically still rising and H$α$ measurements do not reach their plateau. As a result, the H$α$ stellar mass Tully-Fisher relation is steeper (when $M_\star$ is the independent variable) and has larger scatter than its HI counterpart. Interestingly, we confirm the presence of a small fraction of low-mass outliers of the H$α$ relation that are not present when HI velocity widths are used and are not explained by "aperture effects". These appear to be highly disturbed systems for which H$α$ widths do not provide a reliable estimate of the rotational velocity. Our analysis reaffirms the importance of taking into account differences in velocity definitions as well as tracers used when interpreting offsets from the Tully-Fisher relation, at both low and high redshifts and when comparing with simulations. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Comments: 18 pages, 12 figures. Accepted for publication in MNRAS

VERTICO IV: Environmental Effects on the Gas Distribution and Star Formation Efficiency of Virgo Cluster Spirals

Authors: Vicente Villanueva, Alberto D. Bolatto, Stuart Vogel, Tobias Brown, Christine D. Wilson, Nikki Zabel, Sara Ellison, Adam R. H. Stevens, Maria Jesus Jimenez Donaire, Kristine Spekkens, Mallory Thorp, Timothy A. Davis, Laura C. Parker, Ian D. Roberts, Dhruv Bisaria, Alessandro Boselli, Barbara Catinella, Aeree Chung, Luca Cortese, Bumhyun Lee, Adam Watts

Abstract: We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We determine their scale-lengths for the molecular and stellar components and find a roughly 3:5 ratio between them compared to $\sim$1:1 in field galaxies, indicating th… ▽ More We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We determine their scale-lengths for the molecular and stellar components and find a roughly 3:5 ratio between them compared to $\sim$1:1 in field galaxies, indicating that the CO emission is more centrally concentrated than the stars. We compute $R_{\rm mol}$ as a function of different physical quantities. While the spatially-resolved $R_{\rm mol}$ on average decreases with increasing radius, we find that the mean molecular-to-atomic gas ratio within the stellar effective radius $R_{\rm e}$, $R_{\rm mol}(r<R_{\rm e})$, shows a systematic increase with the level of H$_{\rm I}$, truncation and/or asymmetry (H$_{\rm I}$ perturbation). Analysis of the molecular- and the atomic-to-stellar mass ratios within $R_{\rm e}$, $R^{\rm mol}_{\star}(r<R_{\rm e})$ and $R^{\rm atom}_{\star}(r<R_{\rm e})$, shows that VERTICO galaxies have increasingly lower $R^{\rm atom}_{\star}(r<R_{\rm e})$ for larger levels of H$_{\rm I}$perturbation (compared to field galaxies matched in stellar mass), but no significant change in $R^{\rm mol}_{\star}(r<R_{\rm e})$. We also measure a clear systematic decrease of the SFE$_{\rm mol}$ within $R_{\rm e}$, SFE$_{\rm mol}(r<R_{\rm e})$, with increasingly perturbed H$_{\rm I}$. Therefore, compared to galaxies from the field, VERTICO galaxies are more compact in CO emission in relation to their stellar distribution, but increasingly perturbed atomic gas increases their $R_{\rm mol}$ and decreases the efficiency with which their molecular gas forms stars. (abridged) △ Less

Submitted 1 November, 2022; v1 submitted 11 October, 2022; originally announced October 2022.

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

Constraining bosonic asymmetric dark matter with neutron star mass-radius measurements

Authors: Nathan Rutherford, Geert Raaijmakers, Chanda Prescod-Weinstein, Anna Watts

Abstract: Neutron stars can accumulate asymmetric dark matter (ADM) in their interiors, which affects the neutron star's measurable properties and makes compact objects prime targets to search for ADM. In this work, we use Bayesian inference to explore potential neutron star mass-radius measurements, from current and future x-ray telescopes, to constrain the bosonic ADM parameters for the case where bosonic… ▽ More Neutron stars can accumulate asymmetric dark matter (ADM) in their interiors, which affects the neutron star's measurable properties and makes compact objects prime targets to search for ADM. In this work, we use Bayesian inference to explore potential neutron star mass-radius measurements, from current and future x-ray telescopes, to constrain the bosonic ADM parameters for the case where bosonic ADM has accumulated in the neutron star interior. We find that the current uncertainties in the baryonic equation of state do not allow for constraints on the ADM parameter space to be made. However, we also find that ADM cannot be excluded and the inclusion of bosonic ADM in neutron star cores relaxes the constraints on the baryonic equation of state space. If the baryonic equation of state were more tightly constrained independent of ADM, we find that statements about the ADM parameter space could be made. In particular, we find that the high bosonic ADM particle mass ($m_χ$) and low effective self-interaction strength ($g_χ/m_φ)$ regime is disfavored due to the observationally and theoretically motivated constraint that neutron stars must have at least a mass of $1 \, \mathrm{M_\odot}$. However, within the remaining parameter space, $m_χ$ and $g_χ/m_φ$ are individually unconstrained. On the other hand, the ADM mass-fraction, i.e., the fraction of ADM mass inside the neutron star, can be constrained by such neutron star measurements. △ Less

Submitted 2 June, 2023; v1 submitted 5 August, 2022; originally announced August 2022.

Comments: 19 pages, 7 figures. This paper was made to be as similar to the PRD version as possible

Journal ref: Phys. Rev. D 107, 103051 (2023)

Horizons: Nuclear Astrophysics in the 2020s and Beyond

Authors: H. Schatz, A. D. Becerril Reyes, A. Best, E. F. Brown, K. Chatziioannou, K. A. Chipps, C. M. Deibel, R. Ezzeddine, D. K. Galloway, C. J. Hansen, F. Herwig, A. P. Ji, M. Lugaro, Z. Meisel, D. Norman, J. S. Read, L. F. Roberts, A. Spyrou, I. Tews, F. X. Timmes, C. Travaglio, N. Vassh, C. Abia, P. Adsley, S. Agarwal , et al. (140 additional authors not shown)

Abstract: Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit… ▽ More Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 96 pages. Submitted to Journal of Physics G

Report number: LA-UR-22-23997

GRMHD Simulations of Accreting Neutron Stars with Non-Dipole Fields

Authors: Pushpita Das, Oliver Porth, Anna Watts

Abstract: NASA's NICER telescope has recently provided evidence for non-dipolar magnetic field structures in rotation-powered millisecond pulsars. These stars are assumed to have gone through a prolonged accretion spin-up phase, begging the question of what accretion flows onto stars with complex magnetic fields would look like. We present results from a suite of GRMHD simulations of accreting neutron stars… ▽ More NASA's NICER telescope has recently provided evidence for non-dipolar magnetic field structures in rotation-powered millisecond pulsars. These stars are assumed to have gone through a prolonged accretion spin-up phase, begging the question of what accretion flows onto stars with complex magnetic fields would look like. We present results from a suite of GRMHD simulations of accreting neutron stars for dipole, quadrupole, and quadrudipolar stellar field geometries. This is a first step towards simulating realistic hotspot shapes in a general relativistic framework to understand hotspot variability in accreting millisecond pulsars. We find that the location and size of the accretion columns resulting in hotspots changes significantly depending on initial stellar field strength and geometry. We also find that the strongest contributions to the stellar torque are from disk-connected fieldlines and the pulsar wind, leading to spin-down in almost all of the parameter regime explored here. We further analyze angular momentum transport in the accretion disk due to large scale magnetic stresses, turbulent stresses, wind- and compressible effects which we identify with convective motions. The disk collimates the initial open stellar flux forming jets. For dipoles, the disk-magnetosphere interaction can either enhance or reduce jet power compared to the isolated case. However for quadrupoles, the disk always leads to an enhanced net open flux making the jet power comparable to the dipolar case. We discuss our results in the context of observed neutron star jets and provide a viable mechanism to explain radio power both in the low- and high-magnetic field case. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: Resubmitted to MNRAS

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

On the relationship between gas content, star-formation, and global HI asymmetry of galaxies on the star-forming main-sequence

Authors: Adam B. Watts, Barbara Catinella, Luca Cortese, Chris Power, Sara L. Ellison

Abstract: Observations have revealed that disturbances in the cold neutral atomic hydrogen (HI) in galaxies are ubiquitous, but the reasons for these disturbances remain unclear. While some studies suggest that asymmetries in integrated HI spectra (global HI asymmetry) are higher in HI-rich systems, others claim that they are preferentially found in HI-poor galaxies. In this work, we utilise the ALFALFA and… ▽ More Observations have revealed that disturbances in the cold neutral atomic hydrogen (HI) in galaxies are ubiquitous, but the reasons for these disturbances remain unclear. While some studies suggest that asymmetries in integrated HI spectra (global HI asymmetry) are higher in HI-rich systems, others claim that they are preferentially found in HI-poor galaxies. In this work, we utilise the ALFALFA and xGASS surveys, plus a sample of post-merger galaxies, to clarify the link between global HI asymmetry and the gas properties of galaxies. Focusing on star-forming galaxies in ALFALFA, we find that elevated global HI asymmetry is not associated with a change in the HI content of a galaxy, and that only the galaxies with the highest global HI asymmetry show a small increase in specific star-formation rate (sSFR). However, we show that the lack of a trend with HI content is because ALFALFA misses the gas-poor tail of the star-forming main-sequence. Using xGASS to obtain a sample of star-forming galaxies that is representative in both sSFR and HI content, we find that global HI asymmetric galaxies are typically more gas-poor than symmetric ones at fixed stellar mass, with no change in sSFR. Our results highlight the complexity of the connection between galaxy properties and global HI asymmetry. This is further confirmed by the fact that even post-merger galaxies show both symmetric and asymmetric HI spectra, demonstrating that merger activity does not always lead to an asymmetric global HI spectrum. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Comments: 10 pages, 6 figures, accepted for publication in MNRAS

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

Global HI asymmetries in IllustrisTNG: a diversity of physical processes disturb the cold gas in galaxies

Authors: Adam B. Watts, Chris Power, Barbara Catinella, Luca Cortese, Adam R. H. Stevens

Abstract: Observations of the cold neutral atomic hydrogen (HI) in and around disc galaxies have revealed that spatial and kinematic asymmetries are commonplace, and are reflected in the global HI spectra. We use the TNG100 box from the IllustrisTNG suite of cosmological simulations to study the conditions under which these asymmetries may arise in current theoretical galaxy formation models. We find that m… ▽ More Observations of the cold neutral atomic hydrogen (HI) in and around disc galaxies have revealed that spatial and kinematic asymmetries are commonplace, and are reflected in the global HI spectra. We use the TNG100 box from the IllustrisTNG suite of cosmological simulations to study the conditions under which these asymmetries may arise in current theoretical galaxy formation models. We find that more than 50% of the sample has at least a 10% difference in integrated flux between the high- and low-velocity half of the spectrum, thus the typical TNG100 galaxy has an HI profile that is not fully symmetric. We find that satellite galaxies are a more asymmetric population than centrals, consistent with observational results. Using halo mass as a proxy for environment, this trend appears to be driven by the satellite population within the virial radius of haloes more massive than $10^{13} M_{\odot}$, typical of medium/large groups. We show that, while the excess of HI asymmetry in group satellites is likely driven by ram pressure, the bulk of the asymmetric HI profiles observed in TNG100 are driven by physical processes able to affect both the central and satellite populations. Our results highlight how asymmetries are not driven solely by environment, and multiple physical processes can produce the same asymmetric shape in global HI spectra. △ Less

Submitted 11 October, 2020; originally announced October 2020.

Comments: 16 pages, 13 figures, 1 appendix, accepted for publication in MRNRAS

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

xGASS: Robust quantification of asymmetries in global HI spectra and their relationship to environmental processes

Authors: Adam B. Watts, Barbara Catinella, Luca Cortese, Chris Power

Abstract: We present an analysis of asymmetries in global HI spectra from the extended GALEX Arecibo SDSS Survey (xGASS), a stellar mass-selected and gas fraction-limited survey which is representative of the HI properties of galaxies in the local Universe. We demonstrate that the asymmetry in a HI spectrum is strongly linked to its signal-to-noise meaning that, contrary to what was done in previous works,… ▽ More We present an analysis of asymmetries in global HI spectra from the extended GALEX Arecibo SDSS Survey (xGASS), a stellar mass-selected and gas fraction-limited survey which is representative of the HI properties of galaxies in the local Universe. We demonstrate that the asymmetry in a HI spectrum is strongly linked to its signal-to-noise meaning that, contrary to what was done in previous works, asymmetry distributions for different samples cannot be compared at face value. We develop a method to account for noise-induced asymmetry and find that the typical galaxy detected by xGASS exhibits higher asymmetry than what can be attributed to noise alone, with 37% of the sample showing asymmetry greater than 10% at an 80% confidence level. We find that asymmetric galaxies contain, on average, 29% less HI mass compared to their symmetric counterparts matched in both stellar mass and signal-to-noise. We also present clear evidence that satellite galaxies, as a population, exhibit more asymmetric HI spectra than centrals and that group central galaxies show a slightly higher rate of HI asymmetries compared to isolated centrals. All these results support a scenario in which environmental processes, in particular those responsible for gas removal, are the dominant driver of asymmetry in xGASS. △ Less

Submitted 12 January, 2020; originally announced January 2020.

Comments: 14 pages, 16 figures, accepted for publication in MNRAS

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.

The Initial Mass Function in the Extended Ultraviolet Disk of M83

Authors: Sarah M. Bruzzese, David A. Thilker, Gerhardt Meurer, Luciana Bianchi, Adam B. Watts, Annette M. N. Ferguson, Armando Gil de Paz, Barry F. Madore, D. Christopher Martin, R. Michael Rich

Abstract: Using Hubble Space Telescope ACS/WFC data we present the photometry and spatial distribution of resolved stellar populations of four fields within the extended ultraviolet disk (XUV disk) of M83. These observations show a clumpy distribution of main-sequence stars and a mostly smooth distribution of red giant branch stars. We constrain the upper-end of the initial mass function (IMF) in the outer… ▽ More Using Hubble Space Telescope ACS/WFC data we present the photometry and spatial distribution of resolved stellar populations of four fields within the extended ultraviolet disk (XUV disk) of M83. These observations show a clumpy distribution of main-sequence stars and a mostly smooth distribution of red giant branch stars. We constrain the upper-end of the initial mass function (IMF) in the outer disk using the detected population of main-sequence stars and an assumed constant star formation rate (SFR) over the last 300 Myr. By comparing the observed main-sequence luminosity function to simulations, we determine the best-fitting IMF to have a power law slope $α=-2.35 \pm 0.3$ and an upper-mass limit $\rm M_{u}=25_{-3}^{+17} \, M_\odot$. This IMF is consistent with the observed H$α$ emission, which we use to provide additional constraints on the IMF. We explore the influence of deviations from the constant SFR assumption, finding that our IMF conclusions are robust against all but strong recent variations in SFR, but these are excluded by causality arguments. These results, along with our similar studies of other nearby galaxies, indicate that some XUV disks are deficient in high-mass stars compared to a Kroupa IMF. There are over one hundred galaxies within 5 Mpc, many already observed with HST, thus allowing a more comprehensive investigation of the IMF, and how it varies, using the techniques developed here. △ Less

Submitted 11 November, 2019; v1 submitted 9 November, 2019; originally announced November 2019.

Comments: MNRAS accepted. 28 pages, 7 Tables, 19 Figures. See published article or contact third author for paper with figures at full resolution. (V2: small error in metadata corrected)

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

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

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

Abstract: We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniqu… ▽ More We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science. △ Less

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

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

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

Physics and astrophysics of strong magnetic field systems with eXTP

Authors: Andrea Santangelo, Silvia Zane, Hua Feng, RenXin Xu, Victor Doroshenko, Enrico Bozzo, Ilaria Caiazzo, Francesco Coti Zelati, Paolo Esposito, Denis González-Caniulef, Jeremy Heyl, Daniela Huppenkothen, Gianluca Israel, ZhaoSheng Li, Lin Lin, Roberto Mignani, Nanda Rea, Mauro Orlandini, Roberto Taverna, Hao Tong, Roberto Turolla, Cristina Baglio, Federico Bernardini, Niccoló Bucciantini, Marco Feroci , et al. (16 additional authors not shown)

Abstract: In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Conso… ▽ More In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Sciences, the eXTP mission is expected to be launched in the mid 2020s. △ Less

Submitted 11 December, 2018; originally announced December 2018.

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