X-ray bursts and neutron star equation of state

Neutron stars (NSs) contain the most extreme forms of matter available in the Universe. They also serve as astrophysical laboratories to study physics under extreme conditions of strong gravity, ultrahigh densities, super-strong magnetic fields and high radiation densities. The densities in NS interiors can exceed the density of saturated nuclear matter (about 3 x 1014 g/cm3) by a large factor. One of the goals of modern physics is to understand the nature of the fundamental interactions. NSs serve as excellent natural laboratories to study strong interactions and to determine the most important properties of matter in atomic nuclei and NSs. The theory of quantum chromodynamics, which can, in principle, treat strong interactions, is computationally intractable for multi-nucleon systems at NS densities. Instead, physicists have developed empirical models of nucleonic interactions, which make conflicting predictions. Laboratory experiments and NS observations are vital to test these theories and drive the progress. NSs observations can be used to place constraints on strong interactions because the forces between the nuclear particles set the stiffness of NS matter. This is encoded in the equation of state (EoS), the relation between pressure and density. For a given NS mass M, the EoS sets its radius R via the stellar structure equations. By measuring the M-R relation, we can recover the EoS at supranuclear densities, which is of major importance to both fundamental physics and astrophysics. It is central to understanding NSs, supernovae, and compact object mergers including at least one NS. To distinguish among the models of strong interactions one needs to measure NS M and R to a precision of a few per cent for several NSs.

Our group has developed state-of-the-art neutron star atmosphere models and methods to measure NS parameters from the evolution of X-ray spectra observed during thermonuclear explosions at the NS surface known as X-ray bursts. For the first time, we have used the exact Compton scattering kernel to construct the most accurate atmosphere models for hot NSs (Suleimanov et al. 2012). We have proposed a novel “cooling tail” method to determine NS parameters. We have also discovered that the cooling tail properties strongly depend on the accretion rate and showed that only the low-accretion-rate, hard-state X-ray bursts can potentially be used for NS parameter determination (Kajava et al. 2014; Poutanen et al. 2014). Using the cooling tail method, we first found that the NS radius in a couple of known X-ray bursters is above 12.5 km (Suleimanov et al. 2011b; Poutanen et al. 2014), which supports stiff equation of state and is consistent with the existence of 2 solar mass pulsars. With the more advanced method of direct fitting atmosphere model to the data (Nättilä et al. 2017), we were able to the get the NS radius in 4U 1702-429 of R=12.4±0.4 km (see Fig. 1). The achived accuracy of the measurement of NS radius is much better than anything that was available before and even much better than is provided by the NICER mission, which is designed for that purpose.

Figure 1: Constraints on the mass-radius of the neutron star in low-mass X-ray binary 4U 1702-429 from direct fits with the atmosphere model to the bursts cooling tail spectra. From Nättilä J. et al. (2017).

Our recent advance to the field includes the development of the NS atmosphere models for rapidly rotating stars accounting for the oblatness of the NS shape and variations of the gravity along the NS surface (Suleimanov et al. 2020, see Fig. 2). Constraints on the neutron star equation of state can also be obtained using pulse profiles of rotation-powered and accreting millisecond pulsars (AMPs) and from X-ray polarization to be observed from AMPs by the Imaging X-ray Polarimeter Explorer (IXPE) to be launched in the end of 2021.

Figure 2: Images of a neutron star rotating at 700 Hz. Upper panels: blackbody case of local isotropic intensity. ower panels: electron-scattering dominated atmosphere. Panels from the left to the right correspond to the inclinations i= 0o, 45o, and 90o. The lines of constant latitudes (every 10o) and longitudes (every 15o) are shown in black. From Suleimanov, Poutanen, Werner (2020).

 

Selected publications: