The hydrogen and helium accreted by X-ray bursting neutron stars is periodically consumed in runaway thermonuclear reactions that cause the entire surface to glow brightly in X-rays for a few seconds. With models of the emission, the mass and radius of the neutron star can be inferred from the observations. By simultaneously probing neutron star masses and radii, X-ray bursts are one of the strongest diagnostics of the nature of matter at extremely high densities. Accurate determinations of these parameters are difficult, however, due to the highly non-ideal nature of the atmospheres where X-ray bursts occur. Observations from X-ray telescopes such as RXTE and NuStar can potentially place strong constraints on nuclear matter once uncertainties in atmosphere models have been reduced. Here we discuss current progress on modeling atmospheres of X-ray bursting neutron stars and some of the challenges still to be overcome.

The heavy element ashes of rp-process hydrogen and helium burning in accreting neutron stars are compressed to high density where they freeze, forming the outer crust of the star. We calculate the chemical separation on freezing for a number of different nuclear mixtures resulting from a range of burning conditions for the rp-process. We confirm the generic result that light nuclei are preferentially retained in the liquid and heavy nuclei in the solid. This is in agreement with the previous study of a 17-component mixture of rp-process ashes by Horowitz et al. (2007), but extends that result to a much larger range of compositions. We also find an alternate phase separation regime for the lightest ash mixtures which does not demonstrate this generic behaviour. With a few exceptions, we find that chemical separation reduces the expected $Q_{\rm imp}$ in the outer crust compared to the initial rp-process ash, where $Q_{\rm imp}$ measures the mean-square dispersion in atomic number $Z$ of the nuclei in the mixture. We find that the fractional spread of $Z$ plays a role in setting the amount of chemical separation and is strongly correlated to the divergence between the two/three-component approximations and the full component model. The contrast in $Y_e$ between the initial rp-process ashes and the equilibrium liquid composition is similar to that assumed in earlier two-component models of compositionally driven convection, except for very light compositions which produce nearly negligible convective driving. We discuss the implications of these results for observations of accreting neutron stars.

We discuss the effect of chemical separation as matter freezes at the base of the ocean of an accreting neutron star, and the subsequent enrichment of the ocean in light elements and inward transport of heat through convective mixing. We extend the steady-state results of Medin & Cumming 2011 to transiently accreting neutron stars, by considering the time-dependent cases of heating during accretion outbursts and cooling during quiescence. Convective mixing is extremely efficient, flattening the composition profile in about one convective turnover time (weeks to months at the base of the ocean). During accretion outbursts, inward heat transport has only a small effect on the temperature profile in the outer layers until the ocean is strongly enriched in light elements, a process that takes hundreds of years to complete. During quiescence, however, inward heat transport rapidly cools the outer layers of the ocean while keeping the inner layers hot. We find that this leads to a sharp drop in surface emission at around a week followed by a gradual recovery as cooling becomes dominated by the crust. Such a dip should be observable in the light curves of these neutron star transients, if enough data is taken at a few days to a month after the end of accretion. If such a dip is definitively observed, it will provide strong constraints on the chemical composition of the ocean and outer crust.

We show that convection driven by chemical separation can significantly affect the cooling curves of accreting neutron stars after they go into quiescence. We calculate the thermal relaxation of the neutron star ocean and crust including the thermal and compositional fluxes due to convection. After the inward propagating cooling wave reaches the base of the neutron star ocean, the ocean begins to freeze, driving chemical separation. The resulting convection transports heat inward, giving much faster cooling of the surface layers than found assuming the ocean cools passively. The light curves including convection show a rapid drop in temperature weeks after outburst. Identifying this signature in observed cooling curves would constrain the temperature and composition of the ocean as well as offer a real time probe of the freezing of a classical multicomponent plasma.

X-ray observations of quiescent X-ray binaries have the potential to provide insight into the structure and the composition of neutron stars. EXO 0748-676 had been actively accreting for over 24 yr before its outburst ceased in late 2008. Subsequent X-ray monitoring revealed a gradual decay of the quiescent thermal emission that can be attributed to cooling of the accretion-heated neutron star crust. In this work, we report on new Chandra and Swift observations that extend the quiescent monitoring to ~5 yr post-outburst. We find that the neutron star temperature remained at ~117 eV between 2009 and 2011, but had decreased to ~110 eV in 2013. This suggests that the crust has not fully cooled yet, which is supported by the lower temperature of ~95 eV that was measured ~4 yr prior to the accretion phase in 1980. Comparing the data to thermal evolution simulations reveals that the apparent lack of cooling between 2009 and 2011 could possibly be a signature of convection driven by phase separation of light and heavy nuclei in the outer layers of the neutron star.

The neutron rich isotope 22Ne may be a significant impurity in carbon and oxygen white dwarfs and could impact how the stars freeze. We perform molecular dynamics simulations to determine the influence of 22Ne in carbon-oxygen-neon systems on liquid-solid phase equilibria. Both liquid and solid phases are present simultaneously in our simulation volumes. We identify liquid, solid, and interface regions in our simulations using a bond angle metric. In general we find good agreement for the composition of liquid and solid phases between our MD simulations and the semi analytic model of Medin and Cumming. The trace presence of a third component, neon, does not appear to strongly impact the chemical separation found previously for two component carbon and oxygen systems. This suggests that small amounts of 22Ne may not qualitatively change how the material in white dwarf stars freezes. However, we do find systematically lower melting temperatures (higher Gamma) in our MD simulations compared to the semi analytic model. This difference seems to grow with impurity parameter Q_imp and suggests a problem with simple corrections to the linear mixing rule for the free energy of multicomponent solid mixtures that is used in the semi analytic model.

We discuss the effect of chemical separation as matter freezes at the base of the ocean of an accreting neutron star, and argue that the retention of light elements in the liquid acts as a source of buoyancy that drives a slow but continual mixing of the ocean, enriching it substantially in light elements, and leading to a relatively uniform composition with depth. We first consider the timescales associated with different processes that can redistribute elements in the ocean, including convection, sedimentation, crystallization, and diffusion. We then calculate the steady state structure of the ocean of a neutron star for an illustrative model in which the accreted hydrogen and helium burns to produce a mixture of O and Se. Even though the H/He burning produces only 2% oxygen by mass, the steady state ocean has an oxygen abundance more than ten times larger, almost 40% by mass. Furthermore, we show that the convective motions transport heat inwards, with a flux of ~ 0.2 MeV per nucleon for an O-Se ocean, heating the ocean and steepening the outwards temperature gradient. The enrichment of light elements and heating of the ocean due to compositionally-driven convection likely have important implications for carbon ignition models of superbursts.

We develop a method for calculating the equilibrium properties of the liquid-solid phase transition in a classical, ideal, multi-component plasma. Our method is a semi-analytic calculation that relies on extending the accurate fitting formulae available for the one-, two-, and three-component plasmas to the case of a plasma with an arbitrary number of components. We compare our results to those of Horowitz, Berry, & Brown (Phys. Rev. E, 75, 066101, 2007), who use a molecular dynamics simulation to study the chemical properties of a 17-species mixture relevant to the ocean-crust boundary of an accreting neutron star, at the point where half the mixture has solidified. Given the same initial composition as Horowitz et al., we are able to reproduce to good accuracy both the liquid and solid compositions at the half-freezing point; we find abundances for most species within 10% of the simulation values. Our method allows the phase diagram of complex mixtures to be explored more thoroughly than possible with numerical simulations. We briefly discuss the implications for the nature of the liquid-solid boundary in accreting neutron stars.

We present numerical simulations of electron-positron pair cascades in the magnetospheres of magnetic neutron stars for a wide range of surface fields (B_p = 10^{12}--10^{15} G), rotation periods (0.1--10 s), and field geometries. This has been motivated by the discovery in recent years of a number of radio pulsars with inferred magnetic fields comparable to those of magnetars. Evolving the cascade generated by a primary electron or positron after it has been accelerated in the inner gap of the magnetosphere, we follow the spatial development of the cascade until the secondary photons and pairs leave the magnetosphere, and we obtain the pair multiplicity and the energy spectra of the cascade pairs and photons under various conditions. Going beyond previous works, which were restricted to weaker fields (B < a few x 10^{12} G), we have incorporated in our simulations detailed treatments of physical processes that are potentially important (especially in the high field regime) but were either neglected or crudely treated before, including photon splitting with the correct selection rules for photon polarization modes, one-photon pair production into low Landau levels for the e^+e^-, and resonant inverse Compton scattering from polar cap hot spots. We discuss the implications of our results for the radio pulsar death line and for the hard X-ray emission from magnetized neutron stars.

Recent observations of thermally emitting isolated neutron stars revealed spectral features that could be interpreted as radiative transitions of He in a magnetized neutron star atmosphere. We present Hartree-Fock calculations of the polarization-dependent photoionization cross sections of the He atom in strong magnetic fields ranging from 10^12 G to 10^14 G. Convenient fitting formulae for the cross sections are given as well as related oscillator strengths for various bound-bound transitions. The effects of finite nucleus mass on the radiative absorption cross sections are examined using perturbation theory.

For sufficiently strong magnetic fields and/or low temperatures, the neutron star surface may be in a condensed state with little gas or plasma above it. Such surface condensation can significantly affect the thermal emission from isolated neutron stars, and may lead to the formation of a charge-depleted acceleration zone ("vacuum gap") in the magnetosphere above the stellar polar cap. Using the latest results on the cohesive property of magnetic condensed matter, we quantitatively determine the conditions for surface condensation and vacuum gap formation in magnetic neutron stars. We find that condensation can occur if the thermal energy kT of the neutron star surface is less than about 8% of its cohesive energy Q_s, and that a vacuum gap can form if the neutron star's rotation axis and magnetic moment point in opposite directions and kT is less than about 4% of Q_s. Thus, vacuum gap accelerators may exist for some neutron stars. Motivated by this result, we also study the physics of pair cascades in the vacuum gap model for photon emission by accelerating electrons and positrons due to both curvature radiation and resonant/nonresonant inverse Compton scattering. Our calculations of the condition of cascade-induced vacuum breakdown and the related pulsar death line/boundary generalize previous works to the superstrong field regime. We find that inverse Compton scatterings do not produce a sufficient number of high energy photons in the gap and thus do not lead to pair cascades for most neutron star parameters. We discuss the implications of our results for the recent observations of neutron star thermal radiation as well as for the detection/non-detection of radio emission from high-B pulsars and magnetars.

The cohesive energy of condensed matter in strong magnetic fields is a fundamental quantity characterizing magnetized neutron star surfaces. The cohesive energy refers to the energy required to pull an atom out of the bulk condensed matter at zero pressure. Theoretical models of pulsar and magnetar magnetospheres depend on the cohesive properties of the surface matter in strong magnetic fields. For example, depending on the cohesive energy of the surface matter, an acceleration zone ("polar gap") above the polar cap of a pulsar may or may not form. Also, condensation of the neutron star surface, if it occurs, can significantly affect thermal emission from isolated neutron stars. We describe our calculations of the cohesive property of matter in strong magnetic fields, and discuss the implications of our results to the recent observations of neutron star surface emission as well as to the detection/non-detection of radio emission from magnetars.

We present new, ab initio calculations of the electronic structure of one-dimensional infinite chains and three-dimensional condensed matter in strong magnetic fields ranging from B=10^12 G to 2x10^15 G, appropriate for observed magnetic neutron stars. At these field strengths, the magnetic forces on the electrons dominate over the Coulomb forces, and to a good approximation the electrons are confined to the ground Landau level. Our calculations are based on the density functional theory, and use a local magnetic exchange-correlation function appropriate in the strong field regime. The band structures of electrons in different Landau orbitals are computed self-consistently. Numerical results of the ground-state energies and electron work functions are given for one-dimensional chains of H, He, C, and Fe. Fitting formulae for the B-dependence of the energies are also provided. For all the field strengths considered in this paper, hydrogen, helium, and carbon chains are found to be bound relative to individual atoms (although for B less than a few x 10^12 G, the relative binding between C and C_infinity is small). Iron chains are significantly bound for B>10^14 G and are weakly bound if at all at B<10^13 G. We also study the cohesive property of three-dimensional condensed matter of H, He, C, and Fe at zero pressure, constructed from interacting chains in a body-centered tetragonal lattice. Such three-dimensional condensed matter is found to be bound relative to individual atoms, with the cohesive energy increasing rapidly with increasing B.

We present new ab initio calculations of the electronic structure of various atoms and molecules in strong magnetic fields ranging from B=10^12 G to 2x10^15 G, appropriate for radio pulsars and magnetars. For these field strengths, the magnetic forces on the electrons dominate over the Coulomb forces, and to a good approximation the electrons are confined to the ground Landau level. Our calculations are based on the density functional theory, and use a local magnetic exchange-correlation function which is tested to be reliable in the strong field regime. Numerical results of the ground-state energies are given for H_N (up to N=10), He_N (up to N=8), C_N (up to N=5) and Fe_N (up to N=3), as well as for various ionized atoms. Fitting formulae for the B-dependence of the energies are also given. In general, as N increases, the binding energy per atom in a molecule, |E_N|/N, increases and approaches a constant value. For all the field strengths considered in this paper, hydrogen, helium, and carbon molecules are found to be bound relative to individual atoms (although for B less than a few x 10^12 G, the relative binding between C and C_2 is small). Iron molecules are not bound at B<10^13 G, but become energetically more favorable than individual atoms at larger field strengths.