Publications

X-PSI has been applied in the following studies. In addition to giving an idea of the scientific applications, these may also be useful for performance benchmarking and planning of computational resource consumption.

If you have used X-PSI for a project and would like to link it here, please contact the X-PSI team and/or submit a pull-request on GitHub.

Early X-PSI development

X-PSI was initiated by Thomas Riley as part of his PhD work at the University of Amsterdam. Chapter 3 of his PhD thesis (Riley 2019) provides an extended technical description of v0.1 of X-PSI and contains results of some parameter recovery tests using synthetic data.

NICER papers

X-PSI has been used in several Pulse Profile Modeling analysis papers using NICER data of rotation-powered millisecond pulsars. These are typically published with a Zenodo repository containing data, model files, X-PSI scripts, samples and post-processing notebooks to enable full reproduction and independent analysis.

Vinciguerra et al. 2024 (ApJ, 961,62) An Updated Mass-Radius Analysis of the 2017-2018 NICER Data Set of PSR J0030+0451. See also the associated Zenodo repository.

Vinciguerra et al. 2023 (ApJ, 959, 55) X-PSI Parameter Recovery for Temperature Map Configurations Inspired by PSR J0030+0451. See also the associated Zenodo repository.

Salmi et al. 2023 (ApJ, 956, 138) Atmospheric Effects on Neutron Star Parameter Constraints with NICER. See also the associated Zenodo repository.

Salmi et al. 2022 (ApJ, 941, 150) The Radius of PSR J0740+6620 from NICER with NICER Background Estimates. See also the associated Zenodo repository.

Riley et al. 2021 (ApJL, 918, L27) A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy. See also the associated Zenodo repository.

Bogdanov et al. 2021 (ApJL, 914, L15) Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. III. Model Description and Verification of Parameter Estimation Codes.

Bogdanov et al. 2019 (ApJL, 887, L26) Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. III. Model Description and Verification of Parameter Estimation Codes.

Bilous et al. 2019 (ApJL, 887, L23) A NICER View of PSR J0030+0451: Evidence for a Global-scale Multipolar Magnetic Field.

Riley et al. 2019 (ApJL, 887, L21) A NICER View of PSR J0030+0451: Millisecond Pulsar Parameter Estimation. See also the associated Zenodo repository.

Other papers

Kini et al. 2024 (MNRAS, 527, 8118) Pulse profile modelling of thermonuclear burst oscillations - II. Handling variability. See also the associated Zenodo repository.

Kini et al. 2023 (MNRAS, 522, 3389) Pulse profile modelling of thermonuclear burst oscillations - I. The effect of neglecting variability. See also the associated Zenodo repository.

Settings

Below we give details of some settings for the Riley et al. 2019 paper. For later papers some changes were made - see papers for details.

Resource consumption: The calculations reported were performed on the Dutch national supercomputer Cartesius, mostly on the Broadwell nodes (i.e., using CPUs only, usually Broadwell architecture). A number of models were applied, with increasing complexity. Typically, ~1000+ cores were used for posterior sampling, when the parallelisation is weighted by time consumed. In total, the sampling problems in R19 required ~450,000 core hours.

Likelihood function settings: On compute nodes, the likelihood function evaluation times ranged from ~1 to ~3 seconds (single-threaded), depending on the model and point in parameter space.[1] The parallelisation was purely via MPI, with all but one process receiving likelihood function evaluation requests.

Resolution settings: Number of surface cells/elements per hot (closed) region:[2] 24x24; Number of phases (linearly spaced) = 100; Number of energies (logarithmically spaced) = 175; Number of rays (per parallel; linear in cosine of ray angle alpha):[3] 200

Interpolation settings: Steffen splines (GSL) were used everywhere apart from for the atmosphere, and for the atmosphere cubic polynomial interpolants were used in four dimensions.

Compilation: Intel compiler collection,[4] with the CORE-AVX2 instruction set, for X-PSI and dependencies (apart from numpy, which was centrally installed).

Footnotes

[1]

Variation between models on the same processor results from a combination of model complexity and resolution settings. Variation as a function of parameters also occurs due to issues like mesh construction on the surface subject to the resolution settings.

[2]

This is the approximate number of cells that have centres lying within a hot region. For models where a hot region has two temperature components these cells are split between the subregions according to area only such that the surface density of cells is commensurate—at least for subregions that are commensurate in shape. The fast precomputation mode was deactivated but if it were activated the distribution of cells would further be weighted by the approximate total number of counts generated by each subregion.

[3]

The rays were integrated for each likelihood function call instead of loading lookup tables from disk, or using an analytic treatment such as a high-order expansion (see, e.g., rayXpanda). Only primary images were included.

[4]

On Cartesius (and Lisa), one can simply execute module load intel/2017b to access these compilers from the toolchain.