A new code for the Hall-driven magnetic evolution of neutron stars

Over the past decade, the numerical modeling of the magnetic field evolution in astrophysical scenarios has become an increasingly important field. In the crystallized crust of neutron stars the evolution of the magnetic field is governed by the Hall induction equation. In this equation the relative contribution of the two terms (Hall term and Ohmic dissipation) varies depending on the local conditions of temperature and magnetic field strength. This results in the transition from the purely parabolic character of the equations to the hyperbolic regime as the magnetic Reynolds number increases, which presents severe numerical problems. Up to now, most attempts to study this problem were based on spectral methods, but they failed in representing the transition to large magnetic Reynolds numbers. We present a new code based on upwind finite differences techniques that can handle situations with arbitrary low magnetic diffusivity and it is suitable for studying the formation of sharp current sheets during the evolution. The code is thoroughly tested in different limits and used to illustrate the evolution of the crustal magnetic field in a neutron star in some representative cases. Our code, coupled to cooling codes, can be used to perform long-term simulations of the magneto-thermal evolution of neutron stars.

Maximum mass of neutron stars

We determine the structure of neutron stars within a Brueckner-Hartree-Fock approach based on realistic nucleon-nucleon, nucleon-hyperon, and hyperon-hyperon interactions. Our results indicate rather low maximum masses below 1.4 solar masses. This feature is insensitive to the nucleonic part of the EOS due to a strong compensation mechanism caused by the appearance of hyperons and represents thus strong evidence for the presence of nonbaryonic "quark" matter in the interior of heavy stars.

Maximum mass of neutron stars

We determine the structure of neutron stars within a Brueckner-Hartree-Fock approach based on realistic nucleon-nucleon, nucleon-hyperon, and hyperon-hyperon interactions. Our results indicate rather low maximum masses below 1.4 solar masses. This feature is insensitive to the nucleonic part of the EOS due to a strong compensation mechanism caused by the appearance of hyperons and represents thus strong evidence for the presence of nonbaryonic "quark" matter in the interior of heavy stars.

On the mass distribution of neutron stars

The distribution of masses for neutron stars is analysed using the Bayesian statistical inference, evaluating the likelihood of the proposed Gaussian peaks by using 54 measured points obtained in a variety of systems. The results strongly suggest the existence of a bimodal distribution of the masses, with the first peak around 1.37 M(circle dot) and a much wider second peak at 1.73 M(circle dot). The results support earlier views related to the different evolutionary histories of the members for the first two peaks, which produces a natural separation (even if no attempt to `label` the systems has been made here). They also accommodate the recent findings of similar to M(circle dot) masses quite naturally. Finally, we explore the existence of a subgroup around 1.25 M(circle dot), finding weak, if any, evidence for it. This recently claimed low-mass subgroup, possibly related to the O-Mg-Ne core collapse events, has a monotonically decreasing likelihood and does not stand out clearly from the rest of the sample.

Force-free twisted magnetospheres of neutron stars

Context. The X-ray spectra observed in the persistent emission of magnetars are evidence for the existence of a magnetosphere. The high-energy part of the spectra is explained by resonant cyclotron upscattering of soft thermal photons in a twisted magnetosphere, which has motivated an increasing number of efforts to improve and generalize existing magnetosphere models. Aims. We want to build more general configurations of twisted, force-free magnetospheres as a first step to understanding the role played by the magnetic field geometry in the observed spectra. Methods. First we reviewed and extended previous analytical works to assess the viability and limitations of semi-analytical approaches. Second, we built a numerical code able to relax an initial configuration of a nonrotating magnetosphere to a force-free geometry, provided any arbitrary form of the magnetic field at the star surface. The numerical code is based on a finite-difference time-domain, divergence-free, and conservative scheme, based of the magneto-frictional method used in other scenarios. Results. We obtain new numerical configurations of twisted magnetospheres, with distributions of twist and currents that differ from previous analytical solutions. The range of global twist of the new family of solutions is similar to the existing semi-analytical models (up to some radians), but the achieved geometry may be quite different. Conclusions. The geometry of twisted, force-free magnetospheres shows a wider variety of possibilities than previously considered. This has implications for the observed spectra and opens the possibility of implementing alternative models in simulations of radiative transfer aiming at providing spectra to be compared with observations.

A highly resistive layer within the crust of X-ray pulsars limits their spin periods

The lack of isolated X-ray pulsars with spin periods longer than 12 s raises the question of where the population of evolved high-magnetic-field neutron stars has gone. Unlike canonical radiopulsars, X-ray pulsars are not subject to physical limits to the emission mechanism nor observational biases against the detection of sources with longer periods. Here we show that a highly resistive layer in the innermost part of the crust of neutron stars naturally limits the spin period to a maximum value of about 10–20 s. This highly resistive layer is expected if the inner crust is amorphous and heterogeneous in nuclear charge, possibly owing to the existence of a nuclear ‘pasta’ phase. Our findings suggest that the maximum period of isolated X-ray pulsars may be the first observational evidence for an amorphous inner crust, whose properties can be further constrained by future X-ray timing missions combined with more detailed models.

The role of nuclear physics in supernovae and the evolution of neutron stars, Neutrino Opacities, Equation of State, Transport Coefficients, and Dark Matter Production

Thesis (Ph.D.)--University of Washington, 2016-08

The saturated magnetic field of the neutron stars

Considering the ferromagnetic screening for the decay of the X-ray neutron star magnetic field in the binary accretion phase, the phase transition of ferromagnetic materials in the crust of neutron star induces the ferromagnetic screening saturation of the accreted crust, which results in the minimum surface magnetic field of the accreting neutron star, about 108 G, if the accreted matter has completely replaced the crust mass of the neutron star. The magnetic field evolution versus accreted mass is given as Bs ∝ ΔM-0.9, and the obtained magnetic field versus spin period relation is consistent with the distribution of the binary X-ray sources and recycled pulsars. The further thermal effect on the magnetic evolution is also studied.

The many lives of magnetized neutron stars

The magnetic field strength at birth is arguably one of the most important properties to determine the evolutionary path of a neutron star. Objects with very high fields, collectively known as magnetars, are characterized by high X-ray quiescent luminosities, occurrence of outbursts, and, for some of them, sporadic giant flares. While the magnetic field strength is believed to drive their collective behaviour, however, the diversity of their properties, and, especially, the observation of magnetar-like bursts from “low-field” pulsars, has been a theoretical puzzle. In this review, we discuss results of long-term simulations following the coupled evolution of the X-ray luminosity and the timing properties for a large, homogeneous sample of X-ray emitting isolated neutron stars, accounting for a range of initial magnetic field strengths, envelope compositions, and neutron star masses. In addition, by following the evolution of magnetic stresses within the neutron star crust, we can also relate the observed magnetar phenomenology to the physical properties of neutron stars, and in particular to their age and magnetic field strength and topology. The dichotomy of “high-B” field pulsars versus magnetars is naturally explained, and occasional outbursts from old, low B-field neutron stars are predicted. We conclude by speculating on the fate of old magnetars, and by presenting observational diagnostics of the neutron star crustal field topology.

Medium effects in the pion-pole mechanism [gamma gamma ->pi(0)->nu(R)(nu)over-bar(L)(nu(L)(nu)over-bar(R))] of neutron star cooling

Nuclear medium effects in the neutrino cooling of neutron stars through the reaction channel γγ→π0 →ν Rν̄L(νLν̄R) are incorporated. Throughout the paper we discuss different possibilities of right-handed neutrinos, massive left-handed neutrinos, and standard massless left-handed neutrinos (reaction is then allowed only with medium modified vertices). It is demonstrated that multiparticle effects suppress the rate of this reaction channel in the dense hadron matter by 6-7 orders of magnitude that does not allow to decrease existing experimental upper limit on the corresponding π0νν̄ coupling. Other possibilities of the manifestation of the given reaction channel in different physical situations, e.g., in the quark color superconducting cores of the most massive neutron stars, are also discussed. We demonstrate that in the color-flavor-locked superconducting phase for temperatures T≲ 0.1-10 MeV (depending on the effective pion mass and the decay width) the process is feasibly the most efficient neutrino cooling process, although the absolute value of the reaction rate is rather small.

Pulsar timing irregularities and the imprint of magnetic field evolution

Context. The rotational evolution of isolated neutron stars is dominated by the magnetic field anchored to the solid crust of the star. Assuming that the core field evolves on much longer timescales, the crustal field evolves mainly though Ohmic dissipation and the Hall drift, and it may be subject to relatively rapid changes with remarkable effects on the observed timing properties. Aims. We investigate whether changes of the magnetic field structure and strength during the star evolution may have observable consequences in the braking index n. This is the most sensitive quantity to reflect small variations of the timing properties that are caused by magnetic field rearrangements. Methods. We performed axisymmetric, long-term simulations of the magneto-thermal evolution of neutron stars with state-of-the-art microphysical inputs to calculate the evolution of the braking index. Relatively rapid magnetic field modifications can be expected only in the crust of neutron stars, where we focus our study. Results. We find that the effect of the magnetic field evolution on the braking index can be divided into three qualitatively different stages depending on the age and the internal temperature: a first stage that may be different for standard pulsars (with n ~ 3) or low field neutron stars that accreted fallback matter during the supernova explosion (systematically n < 3); in a second stage, the evolution is governed by almost pure Ohmic field decay, and a braking index n > 3 is expected; in the third stage, at late times, when the interior temperature has dropped to very low values, Hall oscillatory modes in the neutron star crust result in braking indices of a high absolute value and both positive and negative signs. Conclusions. Current magneto-thermal evolution models predict a large contribution to the timing noise and, in particular, to the braking index, from temporal variations of the magnetic field. Models with strong (≳ 1014 G) multipolar or toroidal components, even with a weak (~1012 G) dipolar field are consistent with the observed trend of the timing properties.

Proper motions of thermally emitting isolated neutron stars measured with Chandra

The remarkable astrometric capabilities of the Chandra Observatory offer the possibility to measure proper motions of X-ray sources with an unprecedented accuracy in this wavelength range. We recently completed a proper motion survey of three of the seven thermally emitting radio-quiet isolated neutron stars (INSs) discovered in the ROSAT all-sky survey. These INSs (RXJ0420.0-5022, RXJ0806.4-4123 and RXJ1308.6+2127) either lack an optical counterpart or have one so faint that ground based or space born optical observations push the current possibilities of the instrumentation to the limit. Pairs of ACIS observations were acquired 3 to 5 years apart to measure the displacement of the sources on the X-ray sky using as a reference the background of extragalactic or remote Galactic X-ray sources. We derive 2 sigma upper limits of 123 mas yr(-1) and 86 mas yr(-1) on the proper motion of RXJ0420.0-5022 and RXJ0806.4-4123, respectively. RXJ1308.6+2127 exhibits a very significant displacement (similar to 9 sigma) yielding mu = 220 +/- 25 mas yr(-1), the second fastest measured among all ROSAT-discovered INSs. The source is probably moving away rapidly from the Galactic plane at a speed which precludes any significant accretion of matter from the interstellar medium. Its transverse velocity of similar to 740 (d/700 pc) km s(-1) might be the largest of all ROSAT INSs and its corresponding spatial velocity lies among the fastest recorded for neutron stars. RXJ1308.6+2127 is thus a middle-aged (age similar to 1 My) high velocity cooling neutron star. We investigate its possible origin in nearby OB associations or from a field OB star. In most cases, the flight time from birth place appears significantly shorter than the characteristic age derived from spin down rate. Overall, the distribution in transverse velocity of the ROSAT INSs is not statistically different from that of normal radio pulsars.

Evolution of magnetic field and spin period in accreting neutron stars

Based on the accretion-induced magnetic field decay model, in which a frozen field and an incompressible fluid are assumed, we obtain the following results: (1) an analytic relation between the magnetic field and spin period, if the fastness parameter of the accretion disk is neglected (The evolutionary tracks of accreting neutron stars in the P-B diagram in our model are different from the equilibrium period lines when the influence of the fastness parameter is taken into account.); (2) the theoretical minimum spin period of an accreting neutron star is max(1.1ms (DeltaM/M(circle dot))(-1)R(6)(-5/14) I(45)(M/M(circle dot))(-1/2),1.1ms (M/M(circle dot))(-1/2) R(6)(17/14)), independent of the accretion rate (X-ray luminosity) but dependent on the total accretion mass, DeltaM; however, the minimum magnetic field depends on the accretion rate; (3) the magnetic field strength decreases faster with time than does the period.

Effects of vacuum structure on neutron stars

We study the equation of state for neutron matter using the Walecka model including quantum corrections for baryons and sigma mesons through a realignment of the vacuum. We next use this equation of state to calculate the radius, mass, and other properties of rotating neutron stars.