Neutron Star Surface Research Articles

MAXI J1957+032: a new accreting millisecond X-ray pulsar in an ultra-compact binary

ABSTRACT The detection of coherent X-ray pulsations at ∼314 Hz (3.2 ms) classifies MAXI J1957+032 as a fast-rotating, accreting neutron star. We present the temporal and spectral analysis performed using NICER observations collected during the latest outburst of the source. Doppler modulation of the X-ray pulsation revealed the ultra-compact nature of the binary system characterized by an orbital period of ∼1 h and a projected semimajor axis of 14 lt-ms. The neutron star binary mass function suggests a minimum donor mass of 1.7 × 10−2 M⊙, assuming a neutron star mass of 1.4 M⊙ and a binary inclination angle lower than 60 deg. This assumption is supported by the lack of eclipses or dips in the X-ray light curve of the source. We characterized the 0.5–10 keV energy spectrum of the source in outburst as the superposition of a relatively cold black-body-like thermal emission compatible with the emission from the neutron star surface and a Comptonization component with photon index consistent with a typical hard state. We did not find evidence for iron K α lines or reflection components.

Discovery of a 584.65 Hz Burst Oscillation in the Low-mass X-Ray Binary 4U 1730–22

Type-I X-ray burst oscillations are powered by thermonuclear energy released on the neutron star (NS) surface in low-mass X-ray binaries (LMXBs), where the burst oscillation frequencies are close to the NS spin rates. In this work, we report the detection of oscillation at 584.65 Hz during the cooling tail of type-I X-ray bursts observed from the accreting NS LMXB 4U 1730–22 on 2022 March 20, by the Neutron star Interior Composition Explorer telescope. The oscillation signal showed a strong Leahy power, P m ∼ 54.04, around 584.65 Hz, which has single-trial and multiple-trial confidence levels of 7.05σ and 4.73σ, respectively. The folded pulse profile of the oscillation in the 0.2–10 keV band showed a sinusoidal shape with the fractional rms amplitude of (8.0 ± 1.1)%. We found the oscillation frequency showed insignificant upward drifting, i.e., less than 0.3 Hz, during the cooling tail, similar to the behavior appearing in accreting millisecond X-ray pulsars (AMXP), and indicate the source could be an AMXP spinning at 1.71 ms.

Outflows and spectral evolution in the eclipsing AMXP SWIFT J1749.4–2807 with NICER, XMM-Newton, and NuSTAR

ABSTRACT The neutron star low-mass X-ray binary SWIFT J1749.4–2807 is the only known eclipsing accreting millisecond X-ray pulsar. In this manuscript, we perform a spectral characterization of the system throughout its 2021, 2-week-long outburst, analysing 11 NICER observations and quasi-simultaneous XMM-Newton and NuSTAR single observations at the outburst peak. The broad-band spectrum is well-modelled with a blackbody component with a temperature of ∼0.6 keV, most likely consistent with a hotspot on the neutron star surface, and a Comptonization spectrum with power-law index Γ ∼ 1.9, arising from a hot corona at ∼12 keV. No direct emission from the disc was found, possibly due to it being too cool. A high truncation radius for the disc, i.e. at ∼20–30 RG, was obtained from the analysis of the broadened profile of the Fe line in the reflection component. The significant detection of a blue-shifted Fe XXVI absorption line at ∼7 keV indicates weakly relativistic X-ray disc winds, which are typically absent in the hard state of X-ray binaries. By comparing the low flux observed during the outburst and the one expected in a conservative mass-transfer, we conclude that mass-transfer in the system is highly non-conservative, as also suggested by the wind detection. Finally, using the NICER spectra alone, we followed the system while it was fading to quiescence. During the outburst decay, as the spectral shape hardened, the hotspot on the neutron star surface cooled down and shrank, a trend which could be consistent with the pure power-law spectrum observed during quiescence.

Constraints on the magnetic field structure in accreting compact objects from aperiodic variability

ABSTRACTWe investigate the aperiodic variability for a relatively large sample of accreting neutron stars and intermediate polars, focusing on the properties of the characteristic break commonly observed in power spectra of accreting objects. In particular, we investigate the relation of the break frequency and the magnetic field strength, both of which are connected to the size of the magnetosphere. We find that for the majority of objects in our sample the measured break frequency values indeed agree with estimated inner radii of the accretion disc, which allows to use observed break frequencies to independently assess the magnetic field strength and structure in accreting compact objects. As a special case, we focus on Hercules X-1 which is a persistent, medium-luminosity X-ray pulsar accreting from its low-mass companion. In the literature, it has been suggested that the complex pulse profiles, the spin-up behaviour and the luminosity-correlation of the cyclotron energy seen in Her X-1 can be explained with a complex magnetic field structure of the neutron star. Here, we connect the measured break frequency to the magnetospheric radius and show that the magnetic field strength derived assuming a dipole configuration is nearly an order of magnitude smaller than the magnetic field strength corresponding to the cyclotron energy. Accordingly, this discrepancy can be explained with the magnetic field having strong multipole components. The multipolar structure would also increase the accreting area on the neutron star surface, explaining why the critical luminosity for accretion column formation is puzzlingly high in this source.

Principal-axis Analysis of the Eddington Tensor for the Early Post-bounce Phase of Rotational Core-collapse Supernovae

Using full Boltzmann neutrino transport, we performed 2D core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2 and 15.0 M ⊙, both rotational and nonrotational. We employed the results obtained in the early post-bounce phase (t ≲ 20 ms) to assess performance under rapid rotation of some closure relations commonly employed in the truncated moment method. We first made a comparison in 1D under spherical symmetry, though, of the Eddington factor p defined in the fluid rest frame (FR). We confirmed that the maximum entropy closure for the Fermionic distribution (MEFD) performs better than others near the proto–neutron star surface, where p < 1/3 occurs, but does not work well even in 1D when the phase-space occupancy satisfies e < 0.5 together with p < 1/3, the condition known to be not represented by MEFD. For the 2D models with the rapid rotation, we employed the principal-axis analysis of the Eddington tensor. We paid particular attention to the direction of the longest principal axis. We observed in FR that it is aligned neither with the radial direction nor with the neutrino flux in 2D, particularly so in convective and/or rapidly rotating regions, the fact not accommodated in the moment method. We repeated the same analysis in the laboratory frame and found again that the direction of the longest principal axis is not well reproduced by MEFD because the interpolation between the optically thick and thin limits is not very accurate in this frame.

A force explosion condition for spherically symmetric core-collapse supernovae

ABSTRACT Understanding which stars explode leaving behind neutron stars and which stars collapse forming black holes remains a fundamental astrophysical problem. We derive an analytic explosion condition for spherically symmetric core-collapse supernovae. The derivation starts with the exact governing equations, considers the balance of integral forces, includes the important dimensionless parameters, and includes an explicit set of self-consistent approximations. The force explosion condition is $\tilde{L}_\nu \tau _g - 0.06 \tilde{\kappa } \gt 0.38$, and only depends upon two dimensionless parameters. The first compares the neutrino power deposited in the gain region with the accretion power, $\tilde{L}_\nu \tau _g = L_{\nu } \tau _g R_{\rm NS}/ (G \dot{M} M_{\rm NS})$. The second, $\tilde{\kappa } = \kappa \dot{M} / \sqrt{G M_{\rm NS} R_{\rm NS}}$, parametrizes the neutrino optical depth in the accreted matter near the neutron star surface. Over the years, many have proposed approximate explosion conditions: the critical neutrino-luminosity, ante-sonic, and time-scale conditions. We are able to derive these other conditions from the force explosion condition, which unifies them all. Using numerical, steady-state and fully hydrodynamic solutions, we test the explosion condition. The success of these tests is promising in two ways. One, the force explosion condition helps to illuminate the underlying physics of explosions. Two, this condition may be a useful explosion diagnostic for more realistic, three-dimensional radiation hydrodynamic core-collapse simulations.

Common patterns in pulse profiles of high-mass X-ray binaries

Context. X-ray pulsars are binary systems which consist of a neutron star in orbit with a mass donor (companion). In these systems the neutron accretes matter from the companion star, which creates accretion columns or hot spots on the neutron star surface and gives rise to pulsations in the X-ray light curve. The pulse profiles carry information about the accretion and magnetic field geometry. Here we present a study and classification of energy-resolved pulse profiles of a sample of X-ray pulsars, focusing on high-mass X-ray binaries. Aims. Our goal is to perform a classification of X-ray pulsars based on their observed pulse profiles and look for correlations between this classification and their principle physical observables. The analysis pipeline is available online. Methods. We analysed the pulse profiles of a sample of X-ray pulsars using data obtained with the X-ray Multi-Mirror Mission (XMM-Newton) and the Nuclear Spectroscopic Telescope Array (NuSTAR). We fitted the energy-resolved pulse profiles with a Fourier series of up to five harmonics. We then used the energy evolution of the different Fourier components to classify the pulse profiles into groups. We investigated relationships between the pulse profile properties and other observables of the systems (e.g. orbital period, magnetic field strength, and luminosity) to study the extreme physics of these systems. Results. The sources were divided into three groups using a classification based on the shape, the dominance of the fitted Fourier harmonics, and their respective evolution with energy. We do not find a conclusive correlation between the pulse profile shapes or groups and other parameters of the systems. However, a weak trend was found when comparing our classification to the sources’ locations in the spin period-orbital period diagram. Further studies are required to confirm this trend. Conclusions. Despite the large variety of pulse profiles of the X-ray pulsars, we found that with our approach clear categories emerge which we use to classify their behaviour as a function of energy. As we do not find a clear relationship between our classification scheme and other parameters, like the luminosity, the magnetic field strength, or the orbital and spin periods, we conclude that X-ray pulse profiles are influenced by other hidden variables.

Explosion of a Minimum-Mass Neutron Star within Relativistic Hydrodynamics

The relativistic hydrodynamics equations are adapted for the spherically symmetric case and the Lagrangian form. They are used to model the explosive disruption of a minimum-mass neutron star: a key ingredient of the stripping model for short gamma-ray bursts. The shock breakout from the neutron star surface accompanied by the acceleration of matter to ultrarelativistic velocities is studied. A comparison with the results of previously published nonrelativistic calculations is made.

Systematic study for the surface properties of neutron stars

[Background] In our earlier work {\bf [Phys. Rev. C 104, 055804 (2021)]}, we studied the surface properties of a neutron star, assuming it as a huge finite nucleus containing protons, neutrons, electrons, and muons. For the first time, we reported these results of a neutron star for a few representative masses. In the present paper, we give a detailed study of these quantities to draw definite conclusions. [Method]To carry forward our earlier idea, the energy density functional of the momentum space of neutron star matter is converted to the coordinate space in a local density approximation. This functional is again used to derive the neutron star surface properties within the coherent density fluctuation model using the weight function obtained from the density profile of the neutron star using the recently developed G3 and widely used NL3 and IU-FSU parameter sets in the context of relativistic mean-field formalism. [Results] The systematic surface properties of the neutron star, such as incompressibility, symmetry energy, slope parameter, and curvature coefficient, is calculated. The volume and surface components of the total symmetry energy are decomposed with the help of the $\kappa$ factor obtained from the volume to surface ratio of the symmetry energies in the liquid drop limit of Danielewicz. The magnitude of the computed surface quantities increases with the neutron star's mass. [Conclusion] The incompressibility $K^{\rm star}$, symmetry energy $S^{\rm star}$, slope parameter $L_{\rm sym}^{\rm star}$ and curvature coefficient $K_{\rm sym}^{\rm star}$ of the neutron stars with different mass are analyzed and found to be model dependent. NL3 is the stiffest equation of state endue us with the higher magnitude of surface quantities as compared to the G3 and IU-FSU forces.

Trajectories and radiation of charged particles in the pulsar magnetosphere

ABSTRACT Trajectories and radiation of the accelerating electrons are studied in the pulsar magnetosphere approximated as the electromagnetic field of the Deutsch’s solutions. Because the electrons are accelerated rapidly to ultra-relativistic velocity near the neutron star surface, the electron velocity vector (and then its trajectory) is derived from the balance between Lorentz force and radiation reaction force, which makes the pitch angle between electron trajectories and magnetic field lines non-zero in most part of the magnetosphere. In such a case, the spectral energy distributions (SEDs) of synchro-curvature radiation for the accelerating electrons with a mono-energetic form are calculated. Our results indicate that: (i) the pitch angle is the function of electron position (r, θ, ϕ) in the open field line regions, and increases with increasing r and θ as well as increasing the inclination angle; (ii) the radius of curvature becomes large along the particle trajectory, and (iii) the SED appears a double peak structure depending on the emission position, where the synchrotron radiation plays an important role in X-ray band and curvature radiation mainly works in GeV band, which is only determined by parameters α and ζ.

Probing spectral and temporal evolution of the neutron star low-mass X-ray binary 4U 1724–30 with AstroSat

ABSTRACT We report the broad-band spectro-temporal study of the poorly studied accreting neutron star (NS) low mass X-ray binary (LMXB) 4U 1724–30 using data from Soft X-ray Telescope (SXT) and Large Area X-ray Proportional Counters (LAXPC) instruments on board AstroSat. The dim persistent LMXB source was observed with AstroSat over 4 epochs in 2017, all of which corresponded to a low-luminosity non-thermal emission dominated (hard/island) emission state with modest spectral evolution. All the X-ray broad-band spectra can be modelled by a combination of thermal emission from the NS boundary layer (BL) or NS surface and a non-thermal emission component possibly originating from the inverse Comptonization of the disc seed photons. We investigate the presence of frequency and energy-dependent variabilities to probe the origin of the disc/coronal fluctuations. We also report the detection of a Type-I X-ray burst displaying a photospheric radius expansion (PRE). During the burst, a hard X-ray shortage in the 30–80 keV energy band and the enhancement of the persistent emission reveal the burst feedback on the overall accretion process. Using the touch-down burst flux ∼4.25 × 10−8 erg s−1 cm−2, the distance of the source is estimated as ∼8.4 kpc.

Observational features of hot spots on neutron stars with a photon sphere

Abstract We describe and illustrate light deflection near neutron stars possessing a photon sphere motivated by modern confirmations of general relativity and its significance for relativistic astrophysics. In this paper, we model the spacetime outside the star with the Schwarzschild metric and summarize the properties of the photon orbits. More especially we illustrate the relativistic appearance of the hot spots on the neutron star surfaces over the celestial screen of an asymptotic observer. We visualize the relativistic images of the neutron star surface for various ratios of neutron star radius to Schwarzschild radius. We establish a noticeable increscent of the visual size of ultra-compact neutron stars with a photon sphere, which could distinguish stars with extremely strong surface gravity with respect to less compact neutron stars. We found an appearance of additional relativistic images of the hot spots leading to a higher magnification of the neutron star luminosity emission regions.

NICER observations of the evidence of Poynting-Robertson drag and disk reflection during type I X-ray bursts from 4U 1636–536

Type I X-ray bursts are the result of an unstable thermonuclear burning of accreting matter on the neutron star (NS) surface. The quick release of energetic X-ray photons during such bursts interacts with the surrounding accretion disk, which raises the accretion rate due to Poynting-Robertson drag and, thus, a fraction of the burst emission is reflected. We analyzed two photospheric radius expansion bursts in the NS low-mass X-ray binary 4U 1636–536 that took place in 2017, using data from Neutron star Interior Composition Explorer. The time-resolved burst spectra showed clear deviations from a blackbody model. The spectral fitting can be significantly improved by introducing either the enhanced persistent emission (thefamodel) or the reflection from the accretion disk (therelxillNSmodel). Thefamodel provides a higher blackbody temperature and higher burst flux compared with therelxillNSmodel. The peak fluxes of two bursts from thefamodel, 4.36 × 10−8erg cm−2s−1and 9.10 × 10−8erg cm−2s−1, are slightly higher than the Eddington limits of mixed hydrogen-helium and pure helium bursts from previous observations, respectively. When the disk reflections have been taken into account simultaneously, the peak fluxes are lower to match the preferred values. We find evidence to support the finding that both the Poynting-Robertson drag and disk reflection have been appeared during these two X-ray bursts. Moreover, the disk reflection may contribute ∼20−30% of the total burst emissions.

Angular dependence of coherent radio emission from magnetars with multipolar magnetic fields

ABSTRACT The recent detection of a fast radio burst (FRB) from a Galactic magnetar secured the fact that neutron stars (NSs) with superstrong magnetic fields are capable of producing these extremely bright coherent radio bursts. One of the leading mechanisms to explain the origin of such coherent radio emission is the curvature radiation process within the dipolar magnetic field structure. It has, however, already been demonstrated that magnetars likely have a more complex magnetic field topology. Here, we critically investigate curvature radio emission in the presence of inclined dipolar and quadrupolar (‘quadrudipolar’) magnetic fields and show that such field structures differ in their angular characteristics from a purely dipolar case. We analytically show that the shape of open field lines can be modified significantly depending on both the ratio of quadrupole to dipole field strength and their inclination angle at the NS surface. This creates multiple points along each magnetic field line that coincide with the observer’s line of sight, and may explain the complex spectral and temporal structure of the observed FRBs. We also find that in quadrudipole, the radio beam can take a wider angular range and the beamwidth can be wider than that in pure dipole. This may explain why the pulse width of the transient radio pulsation from magnetars is as large as that of ordinary radio pulsars.

X-ray spectral-timing variability of 1A 0535+262 during the 2020 giant outburst

Context. The Be/X-ray binary 1A 0535+262 underwent a giant X-ray outburst in November 2020, peaking at ∼1 × 1038 erg s−1 (1–100 keV, 1.8 kpc), the brightest outburst recorded for this source so far. The source was monitored over two orders of magnitude in luminosity with Insight-HXMT, which allowed us to probe the X-ray variability in an unprecedented range of accretion rates. Aims. Our goal is to search for patterns of correlated spectral and timing behavior that can be used to characterize the accretion states in hard X-ray transient pulsars. Methods. We have studied the evolution of the spectral continuum emission using hardness-intensity diagrams and the aperiodic variability of the source by analyzing power density spectra. We have used phenomenological models to fit the various broadband noise components. Results. The hardness-intensity diagram displays three distinct branches that can be identified with different accretion regimes. The characteristic frequency of the noise components correlates with the luminosity. Our observations cover the highest end of this correlation, at luminosities not previously sampled. We have found evidence for a flattening of the correlation at those high luminosities, which might indicate that the accretion disk reached the closest distance from the neutron star surface during the peak of the outburst. We also find evidence for hysteresis in the spectral and timing parameters: at the same luminosity level, the spectrum is harder and the characteristic noise frequency larger during the rise than during the decay of the outburst. Conclusions. As in black-hole binaries and low-mass X-ray binaries, the hardness-intensity diagram represents a useful diagnostic tool for defining the source state in an accreting pulsar. Our timing analysis confirms previous findings from spectral analyses of a hysteresis pattern of variability, where the spectral and timing parameters adopt different values at similar luminosity depending on whether the source is in the rising or decaying phase of the outburst.

Intensity and Polarization Characteristics of Extended Neutron Star Surface Regions

Abstract The surfaces of neutron stars are sources of strongly polarized soft X-rays due to the presence of strong magnetic fields. Radiative transfer mediated by electron scattering and free–free absorption is central to defining local surface anisotropy and polarization signatures. Scattering transport is strongly influenced by the complicated interplay between linear and circular polarizations. This complexity has been captured in a sophisticated magnetic Thomson scattering simulation we recently developed to model the outer layers of fully ionized atmospheres in such compact objects, heretofore focusing on case studies of localized surface regions. Yet, the interpretation of observed intensity pulse profiles and their efficacy in constraining key neutron star geometry parameters is critically dependent upon adding up emission from extended surface regions. In this paper, intensity, anisotropy, and polarization characteristics from such extended atmospheres, spanning considerable ranges of magnetic colatitudes, are determined using our transport simulation. These constitute a convolution of varied properties of Stokes parameter information at disparate surface locales with different magnetic field strengths and directions relative to the local zenith. Our analysis includes full general relativistic propagation of light from the surface to an observer at infinity. The array of pulse profiles for intensity and polarization presented highlights how powerful probes of stellar geometry are possible. Significant phase-resolved polarization degrees in the range of 10%–60% are realized when summing over a variety of surface field directions. These results provide an important background for observations to be acquired by NASA’s new Imaging X-ray Polarimetry Explorer X-ray polarimetry mission.

Quasi-periodic whispers from a transient ULX in M 101: signatures of a fast-spinning neutron star?

ABSTRACT We have studied the unusual time variability of an ultraluminous X-ray source in M 101, 4XMM J140314.2 + 541806 (henceforth, J1403), using Chandra and XMM-Newton data. Over the last two decades, J1403 has shown short-duration outbursts with an X-ray luminosity ∼1–3 × 1039 erg s−1, and longer intervals at luminosities ∼0.5–1 × 1038 erg s−1. The bimodal behaviour and fast outburst evolution (sometimes only a few days) are more consistent with an accretor/propeller scenario for a neutron star than with the canonical outburst cycles of stellar-mass black holes. If this scenario is correct, the luminosities in the accretor and propeller states suggest a fast spin (P ≈ 5 ms) and a low surface magnetic field (B ∼ 1010 G), despite our identification of J1403 as a high-mass X-ray binary. The most striking property of J1403 is the presence of strong ∼600-s quasi-periodic oscillations (QPOs), mostly around frequencies of ≈1.3–1.8 mHz, found at several epochs during the ultraluminous regime. We illustrate the properties of such QPOs, in particular their frequency and amplitude changes between and within observations, with a variety of techniques (Fast Fourier Transforms, Lomb–Scargle periodograms, weighted wavelet Z-transform analysis). The QPO frequency range &amp;lt;10 mHz is an almost unexplored regime in X-ray binaries and ultraluminous X-ray sources. We compare our findings with the (few) examples of very low frequency variability found in other accreting sources, and discuss possible explanations (Lense–Thirring precession of the inner flow or outflow; radiation pressure limit-cycle instability; marginally stable He burning on the neutron star surface).

Hybrid cooling of the Cassiopeia A neutron star

ABSTRACT The observed rapid cooling of the neutron star Cassiopeia A is usually interpreted as being caused by transitions of neutrons and protons in the star’s core from the normal state to the superfluid and superconducting state. However, this so-called ‘minimal’ cooling paradigm faces the problem of numerically simulating the observed anomalously fast drop in the neutron star surface temperature using theoretical neutrino energy losses from superfluid neutrons. As a solution to this problem, I propose a somewhat more complex cooling model, in which, in addition to superfluid neutrons, direct Urca processes from a very small central part of the neutron star core are also involved. Numerical simulations of the cooling trajectory in this scenario show excellent agreement with observations of the Cassiopeia A neutron star. The proposed cooling scenario unambiguously relates the used equation of state and the mass of the neutron star. For a neutron star constructed according to BSk25 equation of state, the most appropriate are the mass $M=1.62\, {\rm M}_{\odot }$ and the radius R = 12.36 km. If BSk24 equation of state is used, then the most suitable solution is $M=1.60\, {\rm M}_{\odot }$ and R = 12.55 km.

Influence of the Gravitational Darkening Effect on the Spectrum of a Hot, Rapidly Rotating Neutron Star

Abstract In this paper, we discuss the influence of the gravitational darkening effect on the emergent spectrum of a fast-rotating, flattened neutron star. Model atmosphere codes always calculate spectra of emergent intensities and fluxes emitted from the unit surface on the star in plane-parallel geometry. Here we took a step beyond that and calculated a small sample grid of theoretical spectra integrated over the distorted surface of a sample rotating neutron star seen by a distant observer at various inclination angles. We assumed parameters like two dimensionless angular velocities Ω ¯ 2 = 0.30 and 0.60, the effective temperature of a nonrotating star T eff = 2.20 × 107 K, the logarithm of the surface gravity of a spherical star log ( g ) = 14.40 (cgs), and inclination angles from i = 0° to i = 90° with step Δi = 10°. We assumed that the atmosphere consists of a mixture of hydrogen and helium with M H = 0.70 and M He = 0.30. At each point on the neutron star surface, we calculated true intensities for local values of parameters (T eff and log ( g ) ), and these monochromatic intensities are next integrated over the whole surface to obtain the emergent spectrum. In this paper, we compute for the first time theoretical spectra of the fast-rotating neutron star. Our work clearly shows that the gravitational darkening effect strongly influences the spectrum and should be included in realistic models of the atmospheres of rotating neutron stars.

Fitting XMM-Newton observations of the AXP 1RXS J170849.0−400910 with four magnetar surface emission models, and predictions for X-ray polarization observations with IXPE

Context. Phase-resolved spectral and spectropolarimetric X-ray observations of magnetars present us with the opportunity to test models of the origin of the X-ray emission from these objects, and to constrain the properties of the neutron star surface and atmosphere. Aims. Our first aim is to use archival XMM-Newton observations of the magnetar 1RXS J170849.0−400910 to ascertain how well four emission models describe the phase-resolved XMM-Newton energy spectra. Our second aim is to evaluate the scientific potential of future spectropolarimetric observations of 1RXS J170849.0−400910 with the Imaging X-ray Polarimetry Explorer (IXPE) scheduled for launch in late 2021. The most salient questions are whether IXPE is able to distinguish between the different emission models, and whether IXPE can unambiguously detect the signatures of quantum electrodynamics (QED) effects in strong magnetic fields. Methods. We used numerical radiation transport calculations for a large number of different system parameters to predict the X-ray flux and polarization energy spectra of the source 1RXS J170849.0−400910. Based on the numerical results, we developed a new model to fit phase-resolved and phase-averaged X-ray spectral (i.e., XMM-Newton and IXPE) and spectropolarimetric (IXPE) data. In order to test the sensitivity of IXPE to strong-field QED effects, we fit a simulated IXPE observation with two versions of the model, i.e., with and without QED effects accounted for. Results. The fixed-ions condensed surface model gives the best description of the phase-resolved XMM-Newton spectra, followed by the blackbody and free-ions condensed surface models. The magnetized atmosphere model gives a poor description of the data and seems to be largely excluded. Simulations show that the IXPE observations of sources such as 1RXS J170849.0−400910 will allow us to cleanly distinguish between high-polarization (blackbody, magnetized atmosphere) and low-polarization (condensed surface) models. If the blackbody or magnetized atmosphere models apply, IXPE can easily prove QED effects based on ∼200 ksec observations as studied here; longer IXPE observation times will be needed for a clear detection in the case of the condensed surface models. Conclusions. The XMM-Newton data have such a good signal-to-noise ratio that they reveal some limitations of the theoretical models. Notwithstanding this caveat, the fits clearly favor the fixed-ions condensed surface and blackbody models over the free-ions condensed surface and magnetized atmosphere models. The IXPE polarization information will greatly help us to figure out how to improve the models. The first detection of strong-field QED effects in the signal from astrophysical sources seems possible if an adequate amount of time is dedicated to the observations.