Magnetic neutron star cooling and microphysics

This is the author accepted manuscript. The final version is available from EDP Sciences via the DOI in this record.We study the relative importance of several recent updates of microphysics input to the neutron star cooling theory and the effects brought about by superstrong magnetic fields of magnetars, including the effects of the Landau quantization in their crusts. We use a finite-difference code for simulation of neutron-star thermal evolution on timescales from hours to megayears with an updated microphysics input. The consideration of short timescales ($\lesssim1$ yr) is made possible by a treatment of the heat-blanketing envelope without the quasistationary approximation inherent to its treatment in traditional neutron-star cooling codes. For the strongly magnetized neutron stars, we take into account the effects of Landau quantization on thermodynamic functions and thermal conductivities. We simulate cooling of ordinary neutron stars and magnetars with non-accreted and accreted crusts and compare the results with observations. Suppression of radiative and conductive opacities in strongly quantizing magnetic fields and formation of a condensed radiating surface substantially enhance the photon luminosity at early ages, making the life of magnetars brighter but shorter. These effects together with the effect of strong proton superfluidity, which slows down the cooling of kiloyear-aged neutron stars, can explain thermal luminosities of about a half of magnetars without invoking heating mechanisms. Observed thermal luminosities of other magnetars are still higher than theoretical predictions, which implies heating, but the effects of quantizing magnetic fields and baryon superfluidity help to reduce the discrepancy.The work of A.P. was supported by the Russian Science Foundation (grant 14-12-00316)

Opacities and spectra of hydrogen atmospheres of moderately magnetized neutron stars

There is observational evidence that central compact objects (CCOs) in supernova remnants have moderately strong magnetic fields $B\sim10^{11}$ G. Meanwhile, available models of partially ionized hydrogen atmospheres of neutron stars with strong magnetic fields are restricted to $B\gtrsim10^{12}$ G. We extend the equation of state and radiative opacities, presented in previous papers for 10^{12}\mbox{ G}\lesssim B \lesssim 10^{15} G, to weaker fields. An equation of state and radiative opacities for a partially ionized hydrogen plasma are obtained at magnetic fields $B$, temperatures $T$, and densities $\rho$ typical for atmospheres of CCOs and other isolated neutron stars with moderately strong magnetic fields. The first- and second-order thermodynamic functions, monochromatic radiative opacities, and Rosseland mean opacities are calculated and tabulated, taking account of partial ionization, for 3\times10^{10}\mbox{ G}\lesssim B\lesssim 10^{12} G, $10^5$ K $\lesssim T\lesssim 10^7$ K, and a wide range of densities. Atmosphere models and spectra are calculated to verify the applicability of the results and to determine the range of magnetic fields and effective temperatures where the incomplete ionization of the hydrogen plasma is important.The work of AYP on calculation of the opacities and polarizabilities (Sect. 4) has been supported by the Russian Science Foundation (grant 14-12-00316). WCGH appreciates use of computer facilities at KIPAC

Ab initio based equation of state of dense water for planetary and exoplanetary modeling

This is the author accepted manuscript. The final version is available from EDP Sciences via the DOI in this record.As a first step toward a multi-phase equation of state for dense water, we develop a temperature-dependent equation of state for dense water covering the liquid and plasma regimes and extending to the super-ionic and gas regimes. This equation of state covers the complete range of conditions encountered in planetary modeling. We use first principles quantum molecular dynamics simulations and its Thomas-Fermi extension to reach the highest pressures encountered in giant planets several times the size of Jupiter. Using these results, as well as the data available at lower pressures, we obtain a parametrization of the Helmholtz free energy adjusted over this extended temperature and pressure domain. The parametrization ignores the entropy and density jumps at phase boundaries but we show that it is sufficiently accurate to model interior properties of most planets and exoplanets. We produce an equation of state given in analytical form that is readily usable in planetary modeling codes and dynamical simulations {\bf (a fortran implementation can be found at http://www.ioffe.ru/astro/H2O/)}. The EOS produced is valid for the entire density range relevant to planetary modeling, {\bf for densities where quantum effects for the ions can be neglected, and for temperatures below 50,000K. We use this equation of state to calculate the mass-radius relationship of exoplanets up to 5,000M_Earth, explore temperature effects in ocean and wet Earth-like planets, and quantify the influence of the water EOS for the core on the gravitational moments of Jupiter.s. Part of this work was supported by the SNR grant PLANETLAB 12-BS04-0015 and the Programme National de Planetologie (PNP) of CNRS-INSU co-funded by CNES. Funding and support from Paris Sciences et Lettres (PSL) university through the project origins and conditions for the emergence of life is also acknowledged. This work was performed using HPC resources from GENCI- TGCC (Grant 2017- A0030406113

Polarized emission from strongly magnetized sources

Anomalous X-ray pulsars (AXPs) and Soft gamma repeaters (SGRs) form together a single class of astrophysical sources, commonly associated to magnetars. New-generation X-ray polarimeters will play a key role in assessing the nature of these sources by directly probing the star magnetic field. In the highly magnetized environment radiation is expected to be strongly polarized and such a measure will be easily within reach of IXPE and eXTP. Polarization measurements will eventually confirm the presence of ultra-strong magnetic fields, probing the magnetar scenario. In this work we will discuss theoretical expectations for the polarization signature of AXPs and SGRs and present numerical simulations for the detector response of the polarimeters currently under construction. We will also show how these sources can be used to test vacuum birefringence, a QED effect predicted by Heisemberg and Euler in the Thirties and not experimentally verified as yet

Neutron star envelopes and thermal radiation from the magnetic surface

The thermal structure of neutron star envelopes is discussed with emphasis on analytic results. Recent progress on the effect of chemical constitution and high magnetic fields on the opacities and the thermal structure is further reviewed in view of the application to pulsar cooling and magnetars

Electron conduction opacities at the transition between moderate and strong degeneracy: Uncertainties and impact on stellar models

Electron conduction opacities are one of the main physics inputs for the calculation of low- and intermediate-mass stellar models, and a critical question is how to bridge calculations for moderate and strong degeneracy, which are necessarily performed adopting different methods. The density-temperature regime at the boundary between moderate and strong degeneracy is in fact crucial for modelling the helium cores of red giant branch stars and the hydrogen/helium envelopes of white dwarfs. Prompted by recently published new, improved calculations of electron thermal conductivities and opacities for moderate degeneracy, we study different, physically motivated prescriptions to bridge these new computations with well established results in the regime of strong degeneracy. We find that these different prescriptions have a sizable impact on the predicted He-core masses at the He-flash (up to 0.01$M_{\odot}$ for initial total masses far from the transition to non-degenerate He-cores, and up to $\sim 0.04M_{\odot}$ for masses around the transition), the tip of the red giant branch (up to $\sim$0.1~mag) and the zero age horizontal branch luminosities (up to 0.03~dex for masses far from the transition, and up to $\sim$0.2~dex around the transition), and white dwarf cooling times (up to 40-45\% at high luminosities, and up to $\sim$25\% at low luminosities). Current empirical constraints on the tip of the red giant branch and the zero age horizontal branch absolute magnitudes do not allow yet to definitely exclude any of these alternative options for the conductive opacities. Tests against observations of slowly-cooling faint WDs in old stellar populations will need to be performed to see whether they can set some more stringent constraints on how to bridge calculations of conductive opacities for moderate and strong degeneracy

Soft excess in the quiescent Be/X-ray pulsar RX J0812.4-3114

We report a 72 ks XMM-Newton observation of the Be/X-ray pulsar (BeXRP) RX J0812.4-3114 in quiescence (L-X approximate to 1.6 x 10(33) erg s(-1)). Intriguingly, we find a two-component spectrum, with a hard power-law (Gamma approximate to 1.5) and a soft blackbody-like excess below keV. The blackbody component is consistent in kT with a prior quiescent Chandra observation reported by Tsygankov et al. and has an inferred blackbody radius of;----10 km, consistent with emission from the entire neutron star (NS) surface. There is also mild evidence for an absorption line at approximate to and/or approximate to 1.4 keV. The hard component shows pulsations at P approximate to 31.908 s (pulsed fraction 0.84 +/- 0.10), agreeing with the pulse period seen previously in outbursts, but 110 pulsations were found in the soft excess (pulsed fraction less than or similar to 31 per cent). We conclude that the pulsed hard component suggests low-level accretion on to the NS poles, while the soft excess seems to originate from the entire NS surface. Wc speculate that, in quiescence, the source switches between a soft, thermal-dominated state (when the propeller effect is at work) and a relatively hard state with low-level accretion, and use the propeller cut-off to estimate the ma.gnetic field of the system to be less than or similar to 8.4 x 10(11)G. We compare the quiescent thermal L-X predicted by the standard deep crustal heating model to our observations and find that RX J0812.4-31.14 has a high thermal L-X, at or above the prediction for minimum cooling mechanisms. This suggests that RX J0812.4-3114 either contains a relatively low-mass NS with minimum cooling, or that the system may be young enough that the NS has not fully cooled from the supernova explosion

Cooling of Dark-Matter Admixed Neutron Stars with density-dependent Equation of State

We propose a dark-matter (DM) admixed density-dependent equation of state where the fermionic DM interacts with the nucleons via Higgs portal. Presence of DM can hardly influence the particle distribution inside neutron star (NS) but can significantly affect the structure as well as equation of state (EOS) of NS. Introduction of DM inside NS softens the equation of state. We explored the effect of variation of DM mass and DM Fermi momentum on the NS EOS. Moreover, DM-Higgs coupling is constrained using dark matter direct detection experiments. Then, we studied cooling of normal NSs using APR and DD2 EOSs and DM admixed NSs using dark-matter modified DD2 with varying DM mass and Fermi momentum. We have done our analysis by considering different NS masses. Also DM mass and DM Fermi momentum are varied for fixed NS mass and DM-Higgs coupling. We calculated the variations of luminosity and temperature of NS with time for all EOSs considered in our work and then compared our calculations with the observed astronomical cooling data of pulsars namely Cas A, RX J0822-43, 1E 1207-52, RX J0002+62, XMMU J17328, PSR B1706-44, Vela, PSR B2334+61, PSR B0656+14, Geminga, PSR B1055-52 and RX J0720.4-3125. It is found that APR EOS agrees well with the pulsar data for lighter and medium mass NSs but cooling is very fast for heavier NS. For DM admixed DD2 EOS, it is found that for all considered NS masses, all chosen DM masses and Fermi momenta agree well with the observational data of PSR B0656+14, Geminga, Vela, PSR B1706-44 and PSR B2334+61. Cooling becomes faster as compared to normal NSs in case of increasing DM mass and Fermi momenta. It is infered from the calculations that if low mass super cold NSs are observed in future that may support the fact that heavier WIMP can be present inside neutron stars.Comment: 24 Pages, 15 Figures and 2 Tables. Version accepted in The European Physical Journal

Reaction rates and transport in neutron stars

Understanding signals from neutron stars requires knowledge about the transport inside the star. We review the transport properties and the underlying reaction rates of dense hadronic and quark matter in the crust and the core of neutron stars and point out open problems and future directions.Comment: 74 pages; commissioned for the book "Physics and Astrophysics of Neutron Stars", NewCompStar COST Action MP1304; version 3: minor changes, references updated, overview graphic added in the introduction, improvements in Sec IV.A.

A variable absorption feature in the X-ray spectrum of a magnetar

Soft gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are slowly rotating, isolated neutron stars that sporadically undergo episodes of long-term flux enhancement (outbursts) generally accompanied by the emission of short bursts of hard X-rays. This behaviour can be understood in the magnetar model, according to which these sources are mainly powered by their own magnetic energy. This is supported by the fact that the magnetic fields inferred from several observed properties of AXPs and SGRs are greater than - or at the high end of the range of - those of radio pulsars. In the peculiar case of SGR 0418+5729, a weak dipole magnetic moment is derived from its timing parameters, whereas a strong field has been proposed to reside in the stellar interior and in multipole components on the surface. Here we show that the X-ray spectrum of SGR 0418+5729 has an absorption line, the properties of which depend strongly on the star's rotational phase. This line is interpreted as a proton cyclotron feature and its energy implies a magnetic field ranging from 2E14 gauss to more than 1E15 gauss.Comment: Nature, 500, 312 (including Supplementary Information