Magnetic field evolution in neutron stars

Neutron stars contain persistent, ordered magnetic fields that are the strongest known in the Universe. However, their magnetic fluxes are similar to those in magnetic A and B stars and white dwarfs, suggesting that flux conservation during gravitational collapse may play an important role in establishing the field, although it might also be modified substantially by early convection, differential rotation, and magnetic instabilities. The equilibrium field configuration, established within hours (at most) of the formation of the star, is likely to be roughly axisymmetric, involving both poloidal and toroidal components. The stable stratification of the neutron star matter (due to its radial composition gradient) probably plays a crucial role in holding this magnetic structure inside the star. The field can evolve on long time scales by processes that overcome the stable stratification, such as weak interactions changing the relative abundances and ambipolar diffusion of charged particles with respect to neutrons. These processes become more effective for stronger magnetic fields, thus naturally explaining the magnetic energy dissipation expected in magnetars, at the same time as the longer-lived, weaker fields in classical and millisecond pulsars.Comment: To appear in Astronomische Nachrichten (Astronomical Notes) as part of the Proceedings of the 5th Potsdam Thinkshop, "Meridional Circulation, Differential Rotation, Solar and Stellar Activity", held 2007 June 24-29. 5 pages, no figure

Rotochemical Heating in Millisecond Pulsars. Formalism and Non-superfluid case

Rotochemical heating originates in a departure from beta equilibrium due to spin-down compression in a rotating neutron star. The main consequence is that the star eventually arrives at a quasi-equilibrium state, in which the thermal photon luminosity depends only on the current value of the spin-down power, which is directly measurable. Only in millisecond pulsars the spin-down power remains high long enough for this state to be reached with a substantial luminosity. We report an extensive study of the effect of this heating mechanism on the thermal evolution of millisecond pulsars, developing a general formalism in the slow-rotation approximation of general relativity that takes the spatial structure of the star fully into account, and using a sample of realistic equations of state to solve the non-superfluid case numerically. We show that nearly all observed millisecond pulsars are very likely to be in the quasi-equilibrium state. Our predicted quasi-equilibrium temperatures for PSR J0437-4715 are only 20% lower than inferred from observations. Accounting for superfluidity should increase the predicted value.Comment: 34 pages, 8 figures, AASTeX. Accepted for publication in Ap

Rotochemical heating in millisecond pulsars with Cooper pairing

When a rotating neutron star loses angular momentum, the reduction in the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium that enhance the neutrino emissivity and produce thermal energy. This mechanism is named rotochemical heating and has previously been studied for neutron stars of non-superfluid matter, finding that they reach a quasi-steady state in which the rate that the spin-down modifies the equilibrium concentrations is the same to that of the neutrino reactions restoring the equilibrium. On the other hand, the neutron star interior is believed to contain superfluid nucleons, which affect the thermal evolution of the star by suppressing the neutrino reactions and the specific heat, and opening new Cooper pairing reactions. In this work we describe the thermal effects of Cooper pairing with spatially uniform energy gaps of neutrons and protons on rotochemical heating in millisecond pulsars (MSPs) when only modified Urca reactions are allowed. We find that the chemical imbalances grow up to a value close to the energy gaps, which is higher than the one of the nonsuperfluid case. Therefore, the surface temperatures predicted with Cooper pairing are higher and explain the recent measurement of MSP J0437-4715.Comment: VIII Symposium in Nuclear Physics and Applications: Nuclear and Particle astrophysics. Appearing in the American Institute of Physics (AIP) conference proceeding

Order-of-magnitude physics of neutron stars

We use basic physics and simple mathematics accessible to advanced undergraduate students to estimate the main properties of neutron stars. We set the stage and introduce relevant concepts by discussing the properties of "everyday" matter on Earth, degenerate Fermi gases, white dwarfs, and scaling relations of stellar properties with polytropic equations of state. Then, we discuss various physical ingredients relevant for neutron stars and how they can be combined in order to obtain a couple of different simple estimates of their maximum mass, beyond which they would collapse, turning into black holes. Finally, we use the basic structural parameters of neutron stars to briefly discuss their rotational and electromagnetic properties.Comment: 13 pages, 3 figures, accepted for publication in European Physical Journal

Constraining a possible time-variation of the gravitational constant through "gravitochemical heating" of neutron stars

A hypothetical time-variation of the gravitational constant $G$ would make neutron stars expand or contract, so the matter in their interiors would depart from beta equilibrium. This induces non-equilibrium weak reactions, which release energy that is invested partly in neutrino emission and partly in internal heating. Eventually, the star arrives at a stationary state in which the temperature remains nearly constant, as the forcing through the change of $G$ is balanced by the ongoing reactions. Using the surface temperature of the nearest millisecond pulsar (PSR J0437$-$4715) inferred from ultraviolet observations and results from theoretical modelling of the thermal evolution, we estimate two upper limits for this variation: (1) $|\dot G/G| < 2 \times 10^{-10}\mathrm{yr}^{-1},$ if the fast, "direct Urca" reactions are allowed, and (2) $|\dot G/G|<4\times 10^{-12}\mathrm{yr}^{-1},$ considering only the slower, "modified Urca" reactions. The latter is among the most restrictive upper limits obtained by other methods.Comment: IAU 2009 JD9 conference proceedings. MmSAIt, vol.80, in press. Paolo Molaro & Elisabeth Vangioni, eds. - 4 pages, 2 figure

Magnetic field evolution and equilibrium configurations in neutron star cores: the effect of ambipolar diffusion

As another step towards understanding the long-term evolution of the magnetic field in neutron stars, we provide the first simulations of ambipolar diffusion in a spherical star. Restricting ourselves to axial symmetry, we consider a charged-particle fluid of protons and electrons carrying the magnetic flux through a motionless, uniform background of neutrons that exerts a collisional drag force on the former. We also ignore the possible impact of beta decays, proton superconductivity, and neutron superfluidity. All initial magnetic field configurations considered are found to evolve on the analytically expected time-scales towards "barotropic equilibria" satisfying the "Grad-Shafranov equation", in which the magnetic force is balanced by the degeneracy pressure gradient, so ambipolar diffusion is choked. These equilibria are so-called "twisted torus" configurations, which include poloidal and toroidal components, the latter restricted to the toroidal volumes in which the poloidal field lines close inside the star. In axial symmetry, they appear to be stable, although they are likely to undergo non-axially symmetric instabilities.Comment: MNRAS, accepte