Accretion disks around stellar-mass black holes (BHs) emit radiation peaking in the soft X-rays when the source is in the thermal state. The emerging photons are polarized and, for symmetry reasons, the polarization integrated over the source is expected to be either parallel or perpendicular to the (projected) disk symmetry axis, because of electron scattering in the disk. However, due to General Relativity effects photon polarization vectors will rotate with respect to their original orientation, by an amount depending on both the BH spin and the observer's inclination. Hence, X-ray polarization measurements may provide important information about strong gravity effects around these sources. Along with the spectral and polarization properties of radiation which reaches directly the observer once emitted from the disk, in this paper we also include the contribution of returning radiation, i.e. photons that are bent by the strong BH gravity to return again on the disk, where they scatter until eventually escaping to infinity. A comparison between our results and those obtained in previous works by different authors show an overall good agreement, despite the use of different code architectures. We finally consider the effects of absorption in the disk material by including more realistic albedo profiles for the disk surface. Our findings in this respect show that considering also the ionization state of the disk may deeply modify the behavior of polarization observables.

Magnetars are believed to host the strongest magnetic fields in the present universe ($B\gtrsim10^{14}$ G) and the study of their persistent emission in the X-ray band offers an unprecendented opportunity to gain insight into physical processes in the presence of ultra-strong magnetic fields. Up to now, most of our knowledge about magnetar sources came from spectral analysis, which allowed to test the resonant Compton scattering scenario and to probe the structure of the star magnetosphere. On the other hand, radiation emitted from magnetar surface is expected to be strongly polarized and its observed polarization pattern bears the imprint of both scatterings onto magnetospheric charges and QED effects as it propagates in the magnetized vacuum around the star. X-ray polarimeters scheduled to fly in the next years will finally allow to exploit the wealth of information stored in the polarization observables. Here we revisit the problem of assessing the spectro-polarimetric properties of magnetar persistent emission. At variance with previous investigations, proper account for more physical surface emission models is made by considering either a condensed surface or a magnetized atmosphere. Results are used to simulate polarimetric observations with the forthcoming Imaging X-ray Polarimetry Explorer (IXPE). We find that X-ray polarimetry will allow to detect QED vacuum effects for all the emission models we considered and to discriminate among them.

Since the birth of X-ray Astronomy, spectacular advances have been seen in the imaging, spectroscopic and timing studies of the hot and violent X-ray Universe, and further leaps forward are expected in the future. On the other hand, polarimetry is very much lagging behind: after the measurements of the Crab Nebula and Scorpius X-1, obtained by OSO-8 in the 70s, no more observations have been performed in the classical X-ray band, even if some interesting results have been obtained in hard X-rays and in soft gamma-rays. The NASA/ASI mission IXPE, scheduled for the launch in 2021, is going to provide for the first time imaging X-ray polarimetry in the 2-8 keV band thanks to its photoelectric polarimeter, coupled with ~25'' angular resolution X-ray mirrors. Its orders of magnitude improvement in sensitivity with respect to the OSO-8 Bragg polarimeter implies scientifically meaningful polarimetric measurements for at least the brightest specimens of most classes of X-ray sources. In 2027, the Chinese-led mission eXTP should also be launched. In addition to timing and spectroscopic instruments, eXTP will have on board photoelectric polarimeters very similar to those of IXPE, but with a total effective area 2-3 times larger. Building on IXPE results, eXTP will increase the number of sources for which significant polarimetric measurements could be obtained. However, further progresses, such as exploring a broader energy range, considering a larger effective area, improving the angular resolution, and performing wide-field polarization measurements, are needed to reach a mature phase for X-ray polarimetry. In the first part of this White Paper we will discuss a few scientific cases in which a next generation X-ray Polarimetry mission can provide significant advances. In the second part, a possible concept for a medium-class Next Generation X-ray Polarimetry (NGXP) mission will be sketched.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

We present the results of an XMM-Newton observation of the slowly rotating ($P = 3.4$ s), highly magnetized ($B \approx 3 \times 10^{13}$ G) radio pulsar PSR J0726-2612. A previous X-ray observation with the Chandra satellite showed that some of the properties of PSR J0726-2612 are similar to those of the X-ray Dim Isolated Neutron Stars (XDINSs), a small class of nearby slow pulsars characterized by purely thermal X-ray spectra and undetected in the radio band. We confirm the thermal nature of the X-ray emission of PSR J0726-2612, which can be fit by the sum of two blackbodies with temperatures $kT_1 = 0.074_{-0.011}^{+0.006}$ keV and $kT_2 = 0.14_{-0.02}^{+0.04}$ keV and emitting radii $R_1 = 10.4_{-2.8}^{+10.8}$ km and $R_2 = 0.5_{-0.3}^{+0.9}$ km, respectively (assuming a distance of 1 kpc). A broad absorption line modeled with a Gaussian profile centred at $0.39_{-0.03}^{+0.02}$ keV is required in the fit. The pulse profile of PSR J0726-2612 is characterized by two peaks with similar intensity separated by two unequal minima, a shape and pulsed fraction that cannot be reproduced without invoking magnetic beaming of the X-ray emission. The presence of a single radio pulse suggests that in PSR J0726-2612 the angles that the dipole axis and the line of sight make with the rotation axis, $\xi$ and $\chi$ respectively, are similar. This geometry differs from that of the two radio-silent XDINSs with a double peaked pulse profile similar to that of PSR J0726-2612, for which $\xi \sim 90^\circ$ and $\chi \sim 45^\circ$ have been recently estimated. These results strengthen the similarity between PSR J0726-2612 and the XDINSs and support the possibility that the lack of radio emission from the latter might simply be due to an unfavourable viewing geometry.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

PSR B0943+10 is a mode-switching radio pulsar characterized by two emission modes with different radio and X-ray properties. Previous studies, based on simple combinations of blackbody and power law models, showed that its X-ray flux can be decomposed in a pulsed thermal plus an unpulsed non-thermal components. However, if PSR B0943+10 is a nearly aligned rotator seen pole-on, as suggested by the radio data, it is difficult to reproduce the high observed pulsed fraction unless magnetic beaming is included. In this work we reanalyze all the available X-ray observations of PSR B0943+10 with simultaneous radio coverage, modeling its thermal emission with polar caps covered by a magnetized hydrogen atmosphere or with a condensed iron surface. The condensed surface model provides good fits to the spectra of both pulsar modes, but, similarly to the blackbody, it can not reproduce the observed pulse profiles, unless an additional power law with an ad hoc modulation is added. Instead, the pulse profiles and phase-resolved spectra are well described using the hydrogen atmosphere model to describe the polar cap emission, plus an unpulsed power law. For the X-ray brighter state (Q-mode) we obtain a best fit with a temperature kT~0.09 keV, an emitting radius R~260 m, a magnetic field consistent with the value of the dipole field of 4x10^12 G inferred from the timing parameters, and a small angle between the magnetic and spin axis, $\xi$=5. The corresponding parameters for the X-ray fainter state (B-mode) are kT~0.08 keV and R~170 m.

In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Sciences, the eXTP mission is expected to be launched in the mid 2020s.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

Isolated Neutron Stars are some of the most exciting stellar objects known to astronomers: they have the most extreme magnetic fields, with values up to $10^{15}$ G, and, with the exception of stellar-mass black holes, they are the most dense stars, with densities of $\approx 10^{14}$ g cm$^{-3}$. As such, they are perfect laboratories to test theories of electromagnetism and nuclear physics under conditions of magnetic field and density unattainable on Earth. In particular, the interaction of radiation with strong magnetic fields is the cause of the {\em vacuum birefringence}, an effect predicted by quantum electrodynamics in 1936 but that lacked an observational evidence until now. Here, we show how the study of the polarisation of the optical radiation from the surface of an isolated neutron star yielded such an observational evidence, opening exciting perspectives for similar studies at other wavelengths.

Recently, numerical calculations of the magnetic field evolution in neutron stars demonstrated the possible existence of a Hall attractor, a stage at which the evolution of the field driven by the Hall cascade ends. The existence of such a stage in neutron star magnetic evolution is very important, and can be potentially probed by observations. Here we discuss three types of objects which could have reached this stage. First, we briefly describe the evolution of normal radio pulsars with ages about a few hundred thousand years. Then we analyse in more detail observations of RX J1856.5-3754, one of the Magnificent Seven, focusing on the surface temperature distribution and comparing model predictions with the temperature map inferred from X-ray observations. Finally, we discuss the necessity of the Hall attractor stage to explain the hypothetical existence of accreting magnetars. We conclude that at the moment there is no direct confirmation of the Hall attractor stage in known sources. However, more detailed observations in the near future can demonstrate existence (or absence) of this stage of the crustal magnetic field evolution.

The recent detection of a large polarization degree in the optical emission of an isolated neutron star led to the suggestion that this has been the first evidence of vacuum polarization in a strong magnetic field, an effect predicted by quantum electrodynamics but never observed before. This claim was challanged in a paper by Capparelli, Maiani & Polosa (2017), according to whom a much higher polarization degree would be necessary to positively identify vacuum polarization. Here we show that their conclusions are biased by several inadequate assumptions and have no impact on the original claim.

Bursts and flares are among the distinctive observational manifestations of magnetars, isolated neutron stars endowed with an ultra-strong magnetic field ($B \approx 10^{14}$--$10^{15}$ G). It is believed that these events arise in a hot electron-positron plasma which remains trapped within the closed magnetic field lines. We developed a simple radiative transfer model to simulate magnetar flare emission in the case of a steady trapped fireball. After dividing the fireball surface in a number of plane-parallel slabs, the local spectral and polarization properties are obtained integrating the radiative transfer equations for the two normal modes. We assume that magnetic Thomson scattering is the dominant source of opacity, and neglect contributions from second-order radiative processes, although double-Compton scattering is accounted for in establishing local thermal equilibrium in the fireball atmospheric layers. The observed spectral and polarization properties as measured by a distant observer are obtained summing the contributions from the patches which are visible for a given viewing geometry by means of a ray-tracing code. The spectra we obtained in the 1-100 keV energy range are thermal and can be described in terms of the superposition of two blackbodies. The blackbody temperature and the emitting area ratio are in broad agreement with the available observations. The predicted linear polarization degree is in general greater than 80% over the entire energy range and should be easily detectable by new-generation X-ray polarimeters, like IXPE, XIPE and eXTP.

The "Magnificent Seven" (M7) are a group of radio-quiet Isolated Neutron Stars (INSs) discovered in the soft X-rays through their purely thermal surface emission. Owing to the large inferred magnetic fields ($B\approx 10^{13}$ G), radiation from these sources is expected to be substantially polarised, independently on the mechanism actually responsible for the thermal emission. A large observed polarisation degree is, however, expected only if quantum-electrodynamics (QED) polarisation effects are present in the magnetised vacuum around the star. The detection of a strongly linearly polarised signal would therefore provide the first observational evidence of QED effects in the strong-field regime. While polarisation measurements in the soft X-rays are not feasible yet, optical polarisation measurements are within reach also for quite faint targets, like the M7 which have optical counterparts with magnitudes $\approx 26$--$28$. Here, we report on the measurement of optical linear polarisation for the prototype, and brightest member, of the class, RX\, J1856.5$-$3754 ($V\sim 25.5$), the first ever for one of the M7, obtained with the Very Large Telescope. We measured a polarisation degree $\mathrm{P.D.} =16.43\% \pm5.26\%$ and a polarisation position angle $\mathrm{P.A.}=145\fdg39\pm9\fdg44$, computed east of the North Celestial Meridian. The $\mathrm{P.D.}$ that we derive is large enough to support the presence of vacuum birefringence, as predicted by QED.

The evolution of magnetic field in isolated neutron stars is one of the most important ingredients in the attempt to build a unified description of these objects. A prediction of field evolution models is the existence of an equilibrium configuration, in which the Hall cascade vanishes. Recent calculations have explored the field structure in this stage, called the Hall attractor. We use X-ray data of near-by, cooling neutron stars to probe this prediction, as these sources are surmised to be close to or at Hall attractor phase. We show that the source RX J1856.5-3754 might be closer to the attractor than other sources of its class. Our modelling indicates that the properties of surface thermal emission, assuming that the star is in the Hall attractor, are in contradiction with the spectral data of RX J1856.5-3754.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

The observed polarization properties of thermal radiation from isolated, cooling neutron stars depend on both the emission processes at the surface and the effects of the magnetized vacuum which surrounds the star. Here we investigate the polarized thermal emission from X-ray Dim Isolated Neutron Stars, taking RX J1856.5-3754 as a representative case. The physical conditions of the star outermost layers in these sources is still debated, and so we consider emission from a magnetized atmosphere and a condensed surface, accounting for the effects of vacuum polarization as the radiation propagates in the star magnetosphere. We have found that, for a significant range of viewing geometries, measurement of the phase-averaged polarization fraction and phase-averaged polarization angle at both optical and X-ray wavelengths allow us to determine whether this neutron star has an atmosphere or a condensed surface. Our results may therefore be relevant in view of future developments of soft X-ray polarimeters.

New-generation X-ray polarimeters currently under development promise to open a new window in the study of high-energy astrophysical sources. Among them, neutron stars appear particularly suited for polarization measurements. Radiation from the (cooling) surface of a neutron star is expected to exhibit a large intrinsic polarization degree due to the star strong magnetic field ($\approx 10^{12}-10^{15}$ G), which influences the plasma opacity in the outermost stellar layers. The polarization fraction and polarization angle as measured by an instrument, however, do not necessary coincide with the intrinsic ones derived from models of surface emission. This is due to the effects of quantum electrodynamics in the highly magnetized vacuum around the star (the vacuum polarization) coupled with the rotation of the Stokes parameters in the plane perpendicular to the line of sight induced by the non-uniform magnetic field. Here we revisit the problem and present an efficient method for computing the observed polarization fraction and polarization angle in the case of radiation coming from the entire surface of a neutron star, accounting for both vacuum polarization and geometrical effects due to the extended emitting region. Our approach is fairly general and is illustrated in the case of blackbody emission from a neutron star with either a dipolar or a (globally) twisted magnetic field.

The study of magnetars is of particular relevance since these objects are the only laboratories where the physics in ultra-strong magnetic fields can be directly tested. Until now, spectroscopic and timing measurements at X-ray energies in soft gamma-repeaters (SGRs) and anomalous X-ray pulsar (AXPs) have been the main source of information about the physical properties of a magnetar and of its magnetosphere. Spectral fitting in the ~ 0.5-10 keV range allowed to validate the "twisted magnetosphere" model, probing the structure of the external field and estimating the density and velocity of the magnetospheric currents. Spectroscopy alone, however, may fail in disambiguating the two key parameters governing magnetospheric scattering (the charge velocity and the twist angle) and is quite insensitive to the source geometry. X-ray polarimetry, on the other hand, can provide a quantum leap in the field by adding two extra observables, the linear polarization degree and the polarization angle. Using the bright AXP 1RXS J170849.0-400910 as a template, we show that phase-resolved polarimetric measurements can unambiguously determine the model parameters, even with a small X-ray polarimetry mission carrying modern photoelectric detectors and existing X-ray optics. We also show that polarimetric measurements can pinpoint vacuum polarization effects and thus provide an indirect evidence for ultra-strong magnetic fields.