The QED vacuum polarization in external monochromatic plane-wave electromagnetic fields is calculated with spatial and temporal variations of the external fields being taken into account. We develop a perturbation theory to calculate the induced electromagnetic current that appears in the Maxwell equations, based on Schwinger's proper-time method, and combine it with the so-called gradient expansion to handle the variation of external fields perturbatively. The crossed field, i.e., the long wavelength limit of the electromagnetic wave is first considered. The eigenmodes and the refractive indices as the eigenvalues associated with the eigenmodes are computed numerically for the probe photon propagating in some particular directions. In so doing, no limitation is imposed on the field strength and the photon energy unlike previous studies. It is shown that the real part of the refractive index becomes less than unity for strong fields, the phenomenon that has been known to occur for high-energy probe photons. We then evaluate numerically the lowest-order corrections to the crossed-field resulting from the field variations in space and time. It is demonstrated that the corrections occur mainly in the imaginary part of the refractive index.

In this paper, we systematically calculate the polarization in soft X-rays emitted from magnetized neutron stars, which are expected to be observed by the next-generation X-ray satellites. Magnetars are one of the targets for these observations. This is because thermal radiation is normally observed in the soft X-ray band, and it is thought to be linearly polarized because of different opacities for two polarization modes of photons in the magnetized atmosphere of neutron stars and the dielectric properties of the vacuum in strong magnetic fields. In their previous study, Taverna et al. illustrated how strong magnetic fields influence the behavior of the polarization observables for radiation propagating in vacuo without addressing a precise, physical emission model. In this paper, we pay attention to the conversion of photon polarization modes that can occur in the presence of an atmospheric layer above the neutron star surface, computing the polarization angle and fraction and systematically changing the magnetic field strength, radii of the emission region, temperature, mass, and radii of the neutron stars. We confirmed that if plasma is present, the effects of mode conversion cannot be neglected when the magnetic field is relatively weak, $B \sim 10^{13} \mathrm{G}$. Our results indicate that strongly magnetized ($B \gtrsim 10^{14} \mathrm{G}$) neutron stars are suitable to detect polarizations, but not-so-strongly magnetized ($B \sim 10^{13} \mathrm{G}$) neutron stars will be the ones to confirm the mode conversion.