Astrometry: The Foundation for Observational Astronomy

Astronomy has seen unprecedented growth in the past century, due to the rise in multiwavelength observations. The foundation for multiwavelength astronomy is given by Astrometry; the science of position and motion determination of celestial bodies. We present a technique of determining equatorial coordinates of celestial bodies from their pixel coordinates. We also present the subsequent results of using this technique in achieving the initial few steps required for the multiwavelength studies of young open clusters

Comparative Study of Different Stellar Tracks and Isochrones

With better understanding of the working physics behind the formation and evolution of stars, it has become increasingly important for observers to use an appropriate set of stellar evolutionary tracks and isochrones with respect to the objects of interest. We present a comparative study of three widely used competent stellar models -- MIST, PARSEC and Siess. We analyze their input physics and the final tracks thus generated, especially focusing on the behaviour of the models during the pre-main sequence phase

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

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 figure

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

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 MNRA

Pulse Profile Modelling of Thermonuclear Burst Oscillations II: Handling variability

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

Atmospheric Effects on Neutron Star Parameter Constraints with NICER

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 parametrization 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. Accepted for publication in Ap

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

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-ray Pulse Simulation and Inference (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 not previously tested, 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 multimodal structure in the posterior surface. Noise particularly impacts our sensitivity to the analysis settings and widths of the posterior distributions. The multimodal 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

X-PSI: A Python package for neutron star X-ray pulse simulation and inference

International audienc

Atmospheric Effects on Neutron Star Parameter Constraints with NICER

International audienceWe 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 parametrization 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

Atmospheric Effects on Neutron Star Parameter Constraints with NICER

<p>Posterior sample files associated with the publication "Atmospheric Effects on Neutron Star Parameter Constraints with NICER" by Salmi et al. (2023; <a href="https://doi.org/10.48550/arXiv.2308.09319">arXiv:2308.09319</a>; <a href="https://doi.org/10.3847/1538-4357/acf49d">https://doi.org/10.3847/1538-4357/acf49d</a>).</p><p>Also included are: the data products; the numeric model files including the telescope calibration products; model modules in the Python language using the X-PSI framework; and Jupyter analysis notebooks.</p><p>Please refer to the README for detailed information.</p&gt