Search

Recent constraints on neutron star mass and radius have advanced our understanding of the equation of state of cold dense matter. Some of them have been obtained by modeling the pulses of three millisecond X-ray pulsars observed by the Neutron Star Interior Composition Explorer (NICER). Here, we present a Bayesian parameter inference for a fourth pulsar, PSR J1231-1411, using the same technique with NICER and XMM-Newton data. When applying a broad mass-inclination prior from radio timing measurements and the emission region geometry model that can best explain the data, we find likely-converged results only when using a limited radius prior. If limiting the radius to be consistent with the previous observational constraints and equation of state analyses, we infer the radius to be $12.6 \pm 0.3$ km and the mass to be $1.04_{-0.03}^{+0.05}$ $M_\odot$, each reported as the posterior credible interval bounded by the $16\,\%$ and $84\,\%$ quantiles. If using an uninformative prior but limited between $10$ and $14$ km, we find otherwise similar results, but $R_{\mathrm{eq}} = 13.5_{-0.5}^{+0.3}$ km for the radius. In both cases, we find a non-antipodal hot region geometry where one emitting spot is at the equator or slightly above, surrounded by a large colder region, and where a non-circular hot region lies close to southern rotational pole. If using a wider radius prior, we only find solutions that fit the data significantly worse. We discuss the challenges in finding the better fitting solutions, possibly related to the weak interpulse feature in the pulse profile., Comment: 22 pages, 13 figures (2 of which are figure sets), 3 tables, accepted for publication in ApJ

Pulse profile modelling (PPM) is a technique for inferring mass, radius and hotspot properties of millisecond pulsars. PPM is now regularly used for analysis of rotation-powered millisecond pulsars (RMPs) with data from the Neutron Star Interior Composition ExploreR (NICER). Extending PPM to accreting millisecond pulsars (AMPs) is attractive, because they are a different source class featuring bright X-ray radiation from hotspots powered by accretion. In this paper, we present a modification of one of the PPM codes, X-PSI, so that it can be used for AMPs. In particular, we implement a model of an accretion disc and atmosphere model appropriate for the hotspots of AMPs, and improve the overall computational efficiency. We then test parameter recovery with synthetic NICER data in two scenarios with reasonable parameters for AMPs. We find in the first scenario, where the hotspot is large, that we are able to tightly and accurately constrain all parameters including mass and radius. In the second scenario, which is a high inclination system with a smaller hotspot, we find degeneracy between a subset of model parameters and a slight bias in the inferred mass and radius. This analysis of synthetic data lays the ground work for future analysis of AMPs with NICER data. Such an analysis could be complemented by future (joint) analysis of polarization data from the Imaging X-ray Polarimetry Explorer (IXPE)., Comment: 14 pages, 5 figures, 3 tables. Paper submitted to MNRAS

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 ($\chi$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 $\chi$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 $\chi$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., Comment: 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

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., Comment: 27 pages, 16 figures, + appendices. Accepted in ApJL. Data files and reproduction package available in Zenodo

Let $\delta$ be a derivation in a $K$-algebra $R$ and let ${\rm Aut_{\delta}}(R)$ be the isotropy group with respect to the natural conjugation action of $\rm Aut(R)$ of $K$-automorphisms on the set $\rm Der(R)$ of $K$-derivations: that is, the subgroup of automorphisms that commute with the derivation. We explore the characterization of ${\rm Aut_{\delta}}(R)$ for quantum Weyl algebras and we prove that in the case of the Jordanian plane it is always a finite cyclic group. Furthermore, we obtain an arithmetic characterization that allows us to trivially demonstrate that these algebras are not isomorphic.

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., Comment: 17 pages, 9 figures (2 of which are figure sets), 2 tables, accepted for publication in ApJ

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., Comment: Submitted to ApJ

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance., Comment: Accepted for publication in the A&A special issue`Euclid on Sky'

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., Comment: submitted to MNRAS. The Zenodo link will go public after peer review. Comments are welcome

Range filters allow checking whether a query range intersects a given set of keys with a chance of returning a false positive answer, thus generalising the functionality of Bloom filters from point to range queries. Existing practical range filters have addressed this problem heuristically, resulting in high false positive rates and query times when dealing with adversarial inputs, such as in the common scenario where queries are correlated with the keys. We introduce Grafite, a novel range filter that solves these issues with a simple design and clear theoretical guarantees that hold regardless of the input data and query distribution: given a fixed space budget of $B$ bits per key, the query time is $O(1)$, and the false positive probability is upper bounded by $\ell/2^{B-2}$, where $\ell$ is the query range size. Our experimental evaluation shows that Grafite is the only range filter to date to achieve robust and predictable false positive rates across all combinations of datasets, query workloads, and range sizes, while providing faster queries and construction times, and dominating all competitors in the case of correlated queries. As a further contribution, we introduce a very simple heuristic range filter whose performance on uncorrelated queries is very close to or better than the one achieved by the best heuristic range filters proposed in the literature so far., Comment: Accepted for publication in Proceedings of the ACM on Management of Data (SIGMOD 2024)

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., Comment: Accepted for publication in MNRAS. The Zenodo link is now public

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., Comment: 35 pages, 11 figures, 7 tables

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., Comment: 26 pages, 12 figures (2 of which are figure sets), 3 tables. Published in ApJ

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., Comment: 27 pages, 13 figures

A Monotone Minimal Perfect Hash Function (MMPHF) constructed on a set S of keys is a function that maps each key in S to its rank. On keys not in S, the function returns an arbitrary value. Applications range from databases, search engines, data encryption, to pattern-matching algorithms. In this paper, we describe LeMonHash, a new technique for constructing MMPHFs for integers. The core idea of LeMonHash is surprisingly simple and effective: we learn a monotone mapping from keys to their rank via an error-bounded piecewise linear model (the PGM-index), and then we solve the collisions that might arise among keys mapping to the same rank estimate by associating small integers with them in a retrieval data structure (BuRR). On synthetic random datasets, LeMonHash needs 34% less space than the next larger competitor, while achieving about 16 times faster queries. On real-world datasets, the space usage is very close to or much better than the best competitors, while achieving up to 19 times faster queries than the next larger competitor. As far as the construction of LeMonHash is concerned, we get an improvement by a factor of up to 2, compared to the competitor with the next best space usage. We also investigate the case of keys being variable-length strings, introducing the so-called LeMonHash-VL: it needs space within 13% of the best competitors while achieving up to 3 times faster queries than the next larger competitor.

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., Comment: Accepted for publication in MNRAS

Recent cosmic shear studies have shown that higher-order statistics (HOS) developed by independent teams now outperform standard two-point estimators in terms of statistical precision thanks to their sensitivity to the non-Gaussian features of large-scale structure. The aim of the Higher-Order Weak Lensing Statistics (HOWLS) project is to assess, compare, and combine the constraining power of ten different HOS on a common set of $Euclid$-like mocks, derived from N-body simulations. In this first paper of the HOWLS series, we computed the nontomographic ($\Omega_{\rm m}$, $\sigma_8$) Fisher information for the one-point probability distribution function, peak counts, Minkowski functionals, Betti numbers, persistent homology Betti numbers and heatmap, and scattering transform coefficients, and we compare them to the shear and convergence two-point correlation functions in the absence of any systematic bias. We also include forecasts for three implementations of higher-order moments, but these cannot be robustly interpreted as the Gaussian likelihood assumption breaks down for these statistics. Taken individually, we find that each HOS outperforms the two-point statistics by a factor of around two in the precision of the forecasts with some variations across statistics and cosmological parameters. When combining all the HOS, this increases to a $4.5$ times improvement, highlighting the immense potential of HOS for cosmic shear cosmological analyses with $Euclid$. The data used in this analysis are publicly released with the paper., Comment: 33 pages, 24 figures, main results in Fig. 19 & Table 5, version published in A&A

Constant-product market making functions were first introduced by Hayden Adams in 2017 to create Uniswap, a decentralised exchange on Ethereum. This enables users to exchange assets at any given rate. Some variations such as Balancer and Curve were later introduced. In this paper, we analyse the maths that rule this type of protocol. We show that splitting a trade in multiple smaller trades does not impact the final exchange rate. We also show that the protocol is safer and more profitable when no one recompounds their fees.

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., Comment: 29 pages, 18 figures (10 of which are figure sets), 1 animation, 10 tables. Published in ApJ

We perform a scaling and performance portability study of the particle-in-cell scheme for plasma physics applications through a set of mini-apps we name "Alpine", which can make use of exascale computing capabilities. The mini-apps are based on Independent Parallel Particle Layer, a framework that is designed around performance portable and dimension independent particles and fields. We benchmark the simulations with varying parameters such as grid resolutions ($512^3$ to $2048^3$) and number of simulation particles ($10^9$ to $10^{11}$) with the following mini-apps: weak and strong Landau damping, bump-on-tail and two-stream instabilities, and the dynamics of an electron bunch in a charge-neutral Penning trap. We show strong and weak scaling and analyze the performance of different components on several pre-exascale architectures such as Piz-Daint, Cori, Summit and Perlmutter. While the scaling and portability study helps identify the performance critical components of the particle-in-cell scheme in the current state-of-the-art computing architectures, the mini-apps by themselves can be used to develop new algorithms and optimize their high performance implementations targeting exascale architectures.

The initial Climate-Extended Risk Model (CERM) addresses the estimate of climate-related financial risk embedded within a bank loan portfolio, through a climatic extension of the Basel II IRB model. It uses a Gaussian copula model calibrated with non stationary macro-correlations in order to reflect the future evolution of climate-related financial risks. In this complementary article, we propose a stochastic forward-looking methodology to calibrate climate macro-correlation evolution from scientific climate data, for physical and transition efforts specifically. We assume a global physical and transition risk, likened to persistent greenhouse gas (GHG) concentration in the atmosphere. The economic risk is considered stationary and can therefore be calibrated with a backward-looking methodology. We present 4 key principles to model the GDP and we propose to model the economic, physical and transition effort factors with three interdependent stochastic processes allowing for a calibration with seven well defined parameters. These parameters can be calibrated using public data. This new approach means not only to evaluate climate risks without picking any specific scenario but also allows to fill the gap between current one year approach of regulatory and economic capital models and the necessarily long-term view of climate risks by designing a framework to evaluate the resulting credit loss on each step (typically yearly) of the transition path. This new approach could prove instrumental in the 2022 context of central banks weighing the pros and cons of a climate capital charge.