The goal of this work is to develop a simple analytic description of the emission properties (spectrum and polarization) of the condensed, strongly magnetized surface of neutron stars. We have improved the method of van Adelsberg et al. (2005) (arXiv:astro-ph/0406001) for calculating the spectral properties of condensed magnetized surfaces. Using the improved method, we calculate the reflectivity of an iron surface at magnetic field strengths B \sim (10^{12} - 10^{14}) G, with various inclinations of the magnetic field lines and radiation beam with respect to the surface and each other. We construct analytic expressions for the emissivity of this surface as functions of the photon energy, magnetic field strength, and the three angles that determine the geometry of the local problem. Using these expressions, we calculate X-ray spectra for neutron stars with condensed iron surfaces covered by thin partially ionized hydrogen atmospheres. We develop simple analytic descriptions of the intensity and polarization of radiation emitted or reflected by condensed iron surfaces of neutron stars with strong magnetic fields typical for isolated neutron stars. This description provides boundary conditions at the bottom of a thin atmosphere, which are more accurate than previously used approximations. The spectra calculated with this improvement show absorption features different from those in simplified models. The approach developed in this paper yields results that can facilitate modeling and interpretation of the X-ray spectra of isolated, strongly magnetized, thermally emitting neutron stars.

X-ray polarimetry, sometimes alone, and sometimes coupled to spectral and temporal variability measurements and to imaging, allows a wealth of physical phenomena in astrophysics to be studied. X-ray polarimetry investigates the acceleration process, for example, including those typical of magnetic reconnection in solar flares, but also emission in the strong magnetic fields of neutron stars and white dwarfs. It detects scattering in asymmetric structures such as accretion disks and columns, and in the so-called molecular torus and ionization cones. In addition, it allows fundamental physics in regimes of gravity and of magnetic field intensity not accessible to experiments on the Earth to be probed. Finally, models that describe fundamental interactions (e.g. quantum gravity and the extension of the Standard Model) can be tested. We describe in this paper the X-ray Imaging Polarimetry Explorer (XIPE), proposed in June 2012 to the first ESA call for a small mission with a launch in 2017 but not selected. XIPE is composed of two out of the three existing JET-X telescopes with two Gas Pixel Detectors (GPD) filled with a He-DME mixture at their focus and two additional GPDs filled with pressurized Ar-DME facing the sun. The Minimum Detectable Polarization is 14 % at 1 mCrab in 10E5 s (2-10 keV) and 0.6 % for an X10 class flare. The Half Energy Width, measured at PANTER X-ray test facility (MPE, Germany) with JET-X optics is 24 arcsec. XIPE takes advantage of a low-earth equatorial orbit with Malindi as down-link station and of a Mission Operation Center (MOC) at INPE (Brazil).

We present a general formalism for computing self-consistent, numerical solutions to the time-dependent radiative transfer equation in low velocity, multi-level ions undergoing radiative interactions. Recent studies of time-dependent radiative transfer have focused on radiation hydrodynamic and magnetohydrodynamic effects without lines, or have solved time-independent equations for the radiation field simultaneously with time-dependent equations for the state of the medium. In this paper, we provide a fully time-dependent numerical solution to the radiative transfer and atomic rate equations for a medium irradiated by an external source of photons. We use Accelerated Lambda Iteration to achieve convergence of the radiation field and atomic states. We perform calculations for a three-level atomic model that illustrates important time-dependent effects. We demonstrate that our method provides an efficient, accurate solution to the time-dependent radiative transfer problem. Finally, we characterize astrophysical scenarios in which we expect our solutions to be important.

Conversion of photons into axions under the presence of a strong magnetic field can dim the radiation from magnetized astrophysical objects. Here we perform a detailed calculation aimed at quantifying the signatures of photon-axion conversion in the spectra, light curves, and polarization of neutron stars (NSs). We take into account the energy and angle-dependence of the conversion probability and the surface thermal emission from NSs. The latter is computed from magnetized atmosphere models that include the effect of photon polarization mode conversion due to vacuum polarization. The resulting spectral models, inclusive of the general-relativistic effects of gravitational redshift and light deflection, allow us to make realistic predictions for the effects of photon to axion conversion on observed NS spectra, light curves, and polarization signals. We identify unique signatures of the conversion, such as an increase of the effective area of a hot spot as it rotates away from the observer line of sight. For a star emitting from the entire surface, the conversion produces apparent radii that are either larger or smaller (depending on axion mass and coupling strength) than the limits set by NS equations of state. For an emission region that is observed phase-on, photon-axion conversion results in an inversion of the plane of polarization with respect to the no-conversion case. While the quantitative details of the features that we identify depend on NS properties (magnetic field strength, temperature) and axion parameters, the spectral and polarization signatures induced by photon-axion conversion are distinctive enough to make NSs very interesting and promising probes of axion physics.

The X-ray spectra of some magnetized isolated neutron stars (NSs) show absorption features with equivalent widths (EWs) of 50 - 200 eV, whose nature is not yet well known. To explain the prominent absorption features in the soft X-ray spectra of the highly magnetized (B ~ 10^{14} G) X-ray dim isolated NSs (XDINSs), we theoretically investigate different NS local surface models, including naked condensed iron surfaces and partially ionized hydrogen model atmospheres, with semi-infinite and thin atmospheres above the condensed surface. We also developed a code for computing light curves and integral emergent spectra of magnetized neutron stars with various temperature and magnetic field distributions over the NS surface. We compare the general properties of the computed and observed light curves and integral spectra for XDINS RBS\,1223 and conclude that the observations can be explained by a thin hydrogen atmosphere above the condensed iron surface, while the presence of a strong toroidal magnetic field component on the XDINS surface is unlikely. We suggest that the harmonically spaced absorption features in the soft X-ray spectrum of the central compact object (CCO) 1E 1207.4-5209 (hereafter 1E 1207) correspond to peaks in the energy dependence of the free-free opacity in a quantizing magnetic field, known as quantum oscillations. To explore observable properties of these quantum oscillations, we calculate models of hydrogen NS atmospheres with B ~ 10^{10} - 10^{11} G (i.e., electron cyclotron energy E_{c,e} ~ 0.1 - 1 keV) and T_eff = 1 - 3 MK. Such conditions are thought to be typical for 1E 1207. We show that observable features at the electron cyclotron harmonics with EWs \approx 100 - 200 eV can arise due to these quantum oscillations.

In the X-ray spectra of most X-ray dim isolated neutron stars (XDINSs) absorption features with equivalent widths (EWs) of 50 -- 200 eV are observed. We theoretically investigate different models to explain absorption features and compare their properties with the observations. We consider various theoretical models for the magnetized neutron star surface: naked condensed iron surfaces and partially ionized hydrogen model atmospheres, including semi-infinite and thin atmospheres above a condensed surface. The properties of the absorption features (especially equivalent widths) and the angular distributions of the emergent radiation are described for all models. A code for computing light curves and integral emergent spectra of magnetized neutron stars is developed. We assume a dipole surface magnetic field distribution with a possible toroidal component and corresponding temperature distribution. A model with two uniform hot spots at the magnetic poles can also be employed. Light curves and spectra of highly magnetized neutron stars with parameters typical for XDINSs are computed using different surface temperature distributions and various local surface models. Spectra of magnetized model atmospheres are approximated by diluted blackbody spectra with one or two Gaussian lines having parameters, which allow us to describe the model absorption features. To explain the prominent absorption features in the soft X-ray spectra of XDINSs a thin atmosphere above the condensed surface can be invoked, whereas a strong toroidal magnetic field component on the XDINS surfaces can be excluded.

We review the polarization properties of X-ray emission from highly magnetized neutron stars, focusing on emission from the stellar surfaces. We discuss how x-ray polarization can be used to constrain neutron star magnetic field and emission geometry, and to probe strong-field quantum electrodynamics and possibly constrain the properties of axions.

Emission spectra from magnetars in the soft X-ray band likely contain a thermal component emerging directly from the neutron star surface. However, the lack of observed absorption-like features in quiescent spectra makes it difficult to directly constrain physical properties of the atmosphere. We argue that future X-ray polarization measurements represent a promising technique for directly constraining the magnetar magnetic field strength and geometry. We construct models of the observed polarization signal from a finite surface hotspot, using the latest NS atmosphere models for magnetic fields B = 4 x 10^13--5 x 10^14 G. Our calculations are strongly dependent on the NS magnetic field strength and geometry, and are more weakly dependent on the NS equation of state and atmosphere composition. We discuss how the complementary dependencies of phase-resolved spectroscopy and polarimetry might resolve degeneracies that currently hamper the determination of magnetar physical parameters using thermal models.

Balmer-dominated shocks in supernova remnants (SNRs) produce strong hydrogen lines with a two-component profile composed of a narrow contribution from cold upstream hydrogen atoms, and a broad contribution from hydrogen atoms that have undergone charge transfer reactions with hot protons. Observations of emission lines from edge-wise shocks in SNRs can constrain the gas velocity and collisionless electron heating at the shock front. Downstream hydrogen atoms engage in charge transfer, excitation and ionization reactions, defining an interaction region called the shock transition zone. The properties of hot hydrogen atoms produced by charge transfers (called broad neutrals) are critical for accurately calculating the structure and radiation from the shock transition zone. This paper is the third in a series describing the kinetic, fluid and emission properties of Balmer-dominated shocks, and is the first to properly treat the effect of broad neutral kinetics on shock transition zone structure. We use our models to extract shock parameters from observations of Balmer-dominated SNRs. We find that inferred shock velocities and electron temperatures are lower than those of previous calculations by <10% for v_s<1500 km/s, and by 10-30% for v_s>1500 km/s. This effect is primarily due to the fact that excitation by proton collisions and charge transfer to excited levels favor the high speed part of the neutral hydrogen velocity distribution. Our results have a strong dependence on the ratio of electron to proton temperatures, \beta=T_e/T_p, which allows us to construct a relation \beta(v_s) between the temperature ratio and shock velocity. We compare our calculations to previous results by Ghavamian et al. (2007).

We examine the structure of the post-shock region in supernova remnants (SNRs). The ``shock transition zone'' is set up by charge transfer and ionization events between atoms and ions, and has a width $\sim 10^{15}$ cm$^{-2}$ $n^{-1}_0$, where $n_0$ is the total pre-shock density (including both atoms and ions). For Balmer-dominated SNRs with shock velocity $v_s \gtrsim 1000$ km s$^{-1}$, the Rankine-Hugoniot conditions for ion velocity and temperature are obeyed instantly, leaving the full width at half-maximum (FWHM) of the broad H$\alpha$ line versus $v_s$ relation intact. However, the spatial variation in the post-shock densities is relevant to the problem of Ly$\alpha$ resonant scattering in young, core-collapse SNRs. Both two- (pre-shock atoms and ions) and three-component (pre-shock atoms, broad neutrals and ions) models are considered. We compute the spatial emissivities of the broad ($\xi_b$) and narrow ($\xi_n$) H$\alpha$ lines; a calculation of these emissivities in SN 1006 is in general agreement with the computed ones of Raymond et al. (2007). The (dimensionless) spatial shift, $\Theta_{\rm{shift}}$, between the centroids of $\xi_b$ and $\xi_n$ is unique for a given shock velocity and $f_{\rm{ion}}$, the pre-shock ion fraction. Measurements of $\Theta_{\rm{shift}}$ can be used to constrain $n_0$.

RX J1856.5-3754 is one of the brightest nearby isolated neutron stars, and considerable observational resources have been devoted to it. However, current models are unable to satisfactorily explain the data. We show that our latest models of a thin, magnetic, partially ionized hydrogen atmosphere on top of a condensed surface can fit the entire spectrum, from X-rays to optical, of RX J1856.5-3754, within the uncertainties. In our simplest model, the best-fit parameters are an interstellar column density N_H \approx 1x10^20 cm^-2 and an emitting area with R^infty \approx 17 km (assuming a distance to RX J1856.5-3754 of 140 pc), temperature T^infty \approx 4.3x10^5 K, gravitational redshift z_g \sim 0.22, atmospheric hydrogen column y_H \approx 1 g cm^-2, and magnetic field B \approx (3-4)x10^12 G; the values for the temperature and magnetic field indicate an effective average over the surface. We also calculate a more realistic model, which accounts for magnetic field and temperature variations over the neutron star surface as well as general relativistic effects, to determine pulsations; we find there exist viewing geometries that produce pulsations near the currently observed limits. The origin of the thin atmospheres required to fit the data is an important question, and we briefly discuss mechanisms for producing these atmospheres. Our model thus represents the most self-consistent picture to date for explaining all the observations of RX J1856.5-3754.

Theoretical modeling of surface emission from magnetized neutron stars (NSs) requires proper treatment of QED effects, in particular the effect of photon mode conversion due to the ``vacuum resonance'' between plasma and vacuum polarization. Previous NS atmosphere models incorporated this effect approximately, using transfer equations for the photon modes which are inadequate near the vacuum resonance, particularly for field strengths around $B_l\simeq 7\times 10^{13}$ G, where the vacuum resonance occurs near the photosphere. In this paper, we provide an accurate treatment of the QED-induced mode conversion effect in NS atmosphere models, employing both the modal radiative transfer equations, coupled with an accurate mode conversion probability at the resonance, and the full evolution equations for the photon Stokes parameters. In doing so, we are able to quantitatively calculate the effects of vacuum polarization on the atmosphere emission spectra, beam patterns, and polarizations for the entire range of field strengths, $B=10^{12}-10^{15}$ G. We find that for NSs with $B\ga 2 B_l$, vacuum polarization reduces the widths of spectral features, and softens the hard spectral tail typical of magnetized atmosphere models. For $B\la B_l/2$, vacuum polarization does not change the emission spectra, but can significantly affect the polarization signals. We show that vacuum polarization induces a unique energy-dependent linear polarization signature, and that circular polarization can be generated in the magnetospheres of rapidly rotating NSs. We discuss the implications of our results for recent observations of thermally emitting isolated NSs and magnetars, as well as the prospects for future spectral and polarization observations.

We study the equation of state, polarization and radiation properties for nonideal, strongly magnetized plasmas which compose outer envelopes of magnetic neutron stars. Detailed calculations are performed for partially ionized hydrogen atmospheres and for condensed hydrogen or iron surfaces of these stars. This is a companion paper to astro-ph/0511803

Recent observations show that the thermal X-ray spectra of many isolated neutron stars are featureless and in some cases (e.g., RX J1856.5-3754) well fit by a blackbody. Such a perfect blackbody spectrum is puzzling since radiative transport through typical neutron star atmospheres causes noticeable deviation from blackbody. Previous studies have shown that in a strong magnetic field, the outermost layer of the neutron star may be in a condensed solid or liquid form because of the greatly enhanced cohesive energy of the condensed matter. The critical temperature of condensation increases with the magnetic field strength, and can be as high as 10^6 K (for Fe surface at B \sim 10^{13} G or H surface at B \sim a few times 10^{14} G). Thus the thermal radiation can directly emerge from the degenerate metallic condensed surface, without going through a gaseous atmosphere. Here we calculate the emission properties (spectrum and polarization) of the condensed Fe and H surfaces of magnetic neutron stars in the regimes where such condensation may be possible. For a smooth condensed surface, the overall emission is reduced from the blackbody by less than a factor of 2. The spectrum exhibits modest deviation from blackbody across a wide energy range, and shows mild absorption features associated with the ion cyclotron frequency and the electron plasma frequency in the condensed matter. The roughness of the solid condensate (in the Fe case) tends to decrease the reflectivity of the surface, and make the emission spectrum even closer to blackbody. We discuss the implications of our results for observations of dim, isolated neutron stars and magnetars.