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.

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.

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 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.