Source author record

Andreas Reisenegger

Andreas Reisenegger appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2021arXiv

Evolution of random initial magnetic fields in stably stratified and barotropic stars

Long-lived magnetic fields are known to exist in upper main-sequence stars, white dwarfs, and neutron stars. In order to explore possible equilibrium configurations of the magnetic field inside these stars, we have performed 3D-magnetohydrodynamic simulations of the evolution of initially random magnetic fields in stably stratified and barotropic stars with an ideal-gas equation of state using the {\sc Pencil Code}, a high-order finite-difference code for compressible hydrodynamic flows in the presence of magnetic fields. In barotropic (isentropic) stars, we confirm previous results in the sense that all initial magnetic fields we tried decay away, unable to reach a stable equilibrium. In the case of stably stratified stars (with radially increasing specific entropy), initially random magnetic fields appear to always evolve to a stable equilibrium. However, the nature of this equilibrium depends on the dissipation mechanisms considered. If magnetic diffusivity (or hyperdiffusivity) is included, the final state is more axially symmetric and dominated by large wavelengths than the initial state, whereas this is not the case if only viscosity (or hyperviscosity) is present. In real stars, the main mechanism allowing them to relax to an equilibrium is likely to be phase mixing, which we argue is more closely mimicked by viscosity. Therefore, we conclude that, depending on its formation mechanism, the equilibrium magnetic field in these stars could in principle be very asymmetric.

preprint2020arXiv

A redshift database towards the Shapley Supercluster region

We present a database and velocity catalogue towards the region of the Shapley Supercluster based on 18,146 measured velocities for 10,719 galaxies in the approximately 300 square degree area between 12h 43mn 00s < R.A. < 14h 17mn 00s and -23° 30' 00" < Dec < -38° 30' 00". The data catalogue contains velocities from the literature found until 2015. It also includes 5,084 velocities, corresponding to 4,617 galaxies, observed by us at Las Campanas and CTIO observatories and not reported individually until now. Of these, 2,585 correspond to galaxies with no other previously published velocity measurement before 2015. Every galaxy in the velocity database has been identified with a galaxy extracted from the SuperCOSMOS photometric catalogues. We also provide a combined average velocity catalogue for all 10,719 galaxies with measured velocities, adopting the SuperCOSMOS positions as a homogeneous base. A general magnitude cut-off at R2=18.0 mag was adopted (with exceptions only for some of the new reported velocities). In general terms, we confirm the overall structure of the Shapley Supercluster, as found on earlier papers. However, the more extensive velocity data show finer structure, to be discussed in a future publication.

preprint2019arXiv

Neutron star radius measurement from the ultraviolet and soft X-ray thermal emission of PSR J0437-4715

We analyzed the thermal emission from the entire surface of the millisecond pulsar PSR J0437-4715 observed in the ultraviolet and soft X-ray bands. For this, we calculated non-magnetized, partially ionized atmosphere models of hydrogen, helium, and iron compositions and included plasma frequency effects that may affect the emergent spectrum. This is particularly true for the coldest atmospheres composed of iron (up to a few percent changes in the soft X-ray flux). Employing a Markov chain Monte Carlo method, we found that the spectral fits favour a hydrogen atmosphere, disfavour a helium composition and rule out iron atmosphere and blackbody models. By using a Gaussian prior on the dust extinction, based on the latest 3D map of Galactic dust, and accounting for the presence of hot polar caps found in previous work, we found that the hydrogen atmosphere model results in a well-constrained neutron star radius $R_\mathrm{NS} = 13.6^{+0.9}_{-0.8}~\mathrm{km}$ and bulk surface temperature $T^{\infty}_\mathrm{eff}=\left(2.3\pm0.1\right)\times10^5~\mathrm{K}$. This relatively large radius favours a stiff equation of state and disfavours a strange quark composition inside neutron stars.

preprint2019arXiv

On the weak magnetic field of millisecond pulsars: Does it decay before accretion?

Millisecond pulsars are old, fast spinning neutron stars thought to have evolved from classical pulsars in binary systems, where the rapid rotation is caused by the accretion of matter and angular momentum from their companion. During this transition between classical and millisecond pulsars, there is a magnetic field reduction of around 4 orders of magnitude, which is not well understood. According to the standard scenario, the magnetic field is reduced as a consequence of accretion, either through ohmic dissipation or through screening by the accreted matter. We explored an alternative hypothesis in which the magnetic field is reduced through ambipolar diffusion before the accretion. This is particularly effective during the long epoch in which the pulsar has cooled, but has not yet started accreting. This makes the final magnetic field dependent on the evolution time of the companion star and thus its initial mass. We use observed binary systems to constrain the time available for the magnetic field decay based on the current pulsar companion: a helium white dwarf, a carbon-oxygen white dwarf, or another neutron star. Based on a simplified model without baryon pairing, we show that the proposed process agrees with the general distribution of observed magnetic field strengths in binaries, but is not able to explain some mildly recycled pulsars where no significant decay appears to have occurred. We discuss the possibility of other formation channels for these systems and the conditions under which the magnetic field evolution would be set by the neutron star crust rather than the core.

preprint2015arXiv

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.

preprint2014arXiv

Rotochemical heating of millisecond and classical pulsars with anisotropic and density-dependent superfluid gap models

When a rotating neutron star loses angular momentum, the progressive reduction of the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium, inducing reactions that release heat (rotochemical heating). This effect has previously been studied by Fernández & Reisenegger (2005) for non-superfluid neutron stars and by Petrovich & Reisenegger (2010) for superfluid millisecond pulsars. Both studies found that pulsars reach a quasi-steady state in which the compression driving the matter out of beta equilibrium is balanced by the reactions trying to restore the equilibrium. We extend previous studies by considering the effect of density-dependence and anisotropy of the superfluid energy gaps, for the case in which the dominant reactions are the modified Urca processes, the protons are non-superconducting, and the neutron superfluidity is parametrized by models proposed in the literature. By comparing our predictions with the surface temperature of the millisecond pulsar PSR J0437-4715 and upper limits for twenty-one classical pulsars, we find the millisecond pulsar can be only explained by the models with the effectively largest energy gaps (type B models), the classical pulsars require with the gap models that vanish for some angle (type C) and two different envelope compositions. Thus, no single model for neutron superfluidity can simultaneously account for the thermal emission of all available observations of non-accreting neutron stars, possibly due to our neglect of proton superconductivity.

preprint2014arXiv

Stability of Hall equilibria in neutron star crusts

In the solid crusts of neutron stars, the advection of the magnetic field by the current-carrying electrons, an effect known as Hall drift, should play a very important role as the ions remain essentially fixed (as long as the solid does not break). Although Hall drift preserves the magnetic field energy, it has been argued that it may drive a turbulent cascade to scales at which Ohmic dissipation becomes effective, allowing a much faster decay in objects with very strong fields. On the other hand, it has been found that there are "Hall equilibria", i.e., field configurations that are unaffected by Hall drift. Here, we address the crucial question of the stability of these equilibria through axially symmetric (2D) numerical simulations of Hall drift and Ohmic diffusion, with the simplifying assumption of uniform electron density and conductivity. We demonstrate the 2D-stability of a purely poloidal equilibrium, for which Ohmic dissipation makes the field evolve towards an attractor state through adjacent stable configurations, around which damped oscillations occur. For this field, the decay scales with the Ohmic timescale. We also study the case of an unstable equilibrium consisting of both poloidal and toroidal field components that are confined within the crust. This field evolves into a stable configuration, which undergoes damped oscillations superimposed on a slow evolution towards an attractor, just as the purely poloidal one.

preprint2013arXiv

Magnetic fields of neutron stars

Neutron stars contain the strongest magnetic fields known in the Universe. In this paper, I discuss briefly how these magnetic fields are inferred from observations, as well as the evidence for their time-evolution. I show how these extremely strong fields are actually weak in terms of their effects on the stellar structure, as is also the case for magnetic stars on the upper main sequence and magnetic white dwarfs, which have similar total magnetic fluxes. I propose a scenario in which a stable hydromagnetic equilibrium (containing a poloidal and a toroidal field component) is established soon after the birth of the neutron star, aided by the strong compositional stratification of neutron star matter, and this state is slowly eroded by non-ideal magnetohydrodynamic processes such as beta decays and ambipolar diffusion in the core of the star and Hall drift and breaking of the solid in its crust. Over sufficiently long time scales, the fluid in the neutron star core will behave as if it were barotropic, because, depending on temperature and magnetic field strength, beta decays will keep adjusting the composition to the chemical equilibrium state, or ambipolar diffusion will decouple the charged component from the neutrons. Therefore, the still open question regarding stable hydromagnetic equilibria in barotropic fluids will become relevant for the evolution, at least for ``magnetar'' fields, too strong to be stabilized by the solid crust.

preprint2013arXiv

Magnetohydrodynamic equilibria in barotropic stars

Although barotropic matter does not constitute a realistic model for magnetic stars, it would be interesting to confirm a recent conjecture that states that magnetized stars with a barotropic equation of state would be dynamically unstable (Reisenegger 2009). In this work we construct a set of barotropic equilibria, which can eventually be tested using a stability criterion. A general description of the ideal MHD equations governing these equilibria is summarized, allowing for both poloidal and toroidal magnetic field components. A new finite-difference numerical code is developed in order to solve the so-called Grad-Shafranov equation describing the equilibrium of these configurations, and some properties of the equilibria obtained are briefly discussed.

preprint2013arXiv

Stability of magnetic fields in non-barotropic stars: an analytic treatment

Magnetic fields in upper main-sequence stars, white dwarfs, and neutron stars are known to persist for timescales comparable to their lifetimes. From a theoretical perspective this is problematic, as it can be shown that simple magnetic field configurations are always unstable. In non-barotropic stars, stable stratification allows for a much wider range of magnetic field structures than in barotropic stars, and helps stabilize them by making it harder to induce radial displacements. Recent simulations by Braithwaite and collaborators have shown that, in stably stratified stars, random initial magnetic fields evolve into nearly axisymmetric configurations with both poloidal and toroidal components, which then remain stable for some time. It is desirable to provide an analytic study of the stability of such fields. We write an explicit expression for a plausible equilibrium structure of an axially symmetric magnetic field with both poloidal and toroidal components of adjustable strengths, in a non-barotropic static fluid star, and study its stability using the energy principle. We construct a displacement field that should be a reasonable approximation to the most unstable mode of a toroidal field, and confirm Braithwaite's result that a given toroidal field can be stabilized by a poloidal field containing much less energy than the former. This is consistent with the speculation that the toroidal field is the main reservoir powering magnetar activity. The deformation of a neutron star caused by the hidden toroidal field can also cause emission of gravitational waves.

preprint2011arXiv

Gravitational wave emission from a magnetically deformed non-barotropic neutron star

A strong candidate for a source of gravitational waves is a highly magnetised, rapidly rotating neutron star (magnetar) deformed by internal magnetic stresses. We calculate the mass quadrupole moment by perturbing a zeroth-order hydrostatic equilibrium by an axisymmetric magnetic field with a \emph{linked poloidal-toroidal structure}. In this work, we do \emph{not} require the model star to obey a barotropic equation of state (as a realistic neutron star is not barotropic), allowing us to explore the hydromagnetic equilibria with fewer constraints. We derive the relation between the ratio of poloidal-to-total field energy $Λ$ and ellipticity $ε$ and briefly compare our results to those obtained using the barotropic assumption. Then, we present some examples of how our results can be applied to astrophysical contexts. First, we show how our formulae, in conjunction with current gravitational wave (non-)detections of the Crab pulsar and the Cassiopeia A central compact object (Cas A CCO), can be used to constrain the strength of the internal toroidal fields of those objects. We find that, for the Crab pulsar (whose canonical equatorial dipole field strength, inferred from spin down, is $4\times 10^8$ T) to emit detectable gravitational radiation, the neutron star must have a strong toroidal field component, with maximum internal toroidal field strength $B_{\mathrm{tm}}=7\times 10^{12}$ T; for gravitational waves to be detected from the Cas A CCO at 300 Hz, $B_{\mathrm{tm}}\sim 10^{13}$ T, whereas detection at 100 Hz would require $B_{\mathrm{tm}}\sim 10^{14}$ T. Using our results, we also show how the gravitational wave signal emitted by a magnetar immediately after its birth (assuming it is born rapidly rotating, with $Λ\lesssim 0.2$) makes such a newborn magnetar a stronger candidate for gravitational wave detection than, for example, an SGR giant flare.

preprint2011arXiv

Long-period thermal oscillations in superfluid millisecond pulsars

In previous papers, we have shown that, as the rotation of a neutron star slows down, it will be internally heated as a consequence of the progressively changing mix of particles (rotochemical heating). In previously studied cases (non-superfluid neutron stars or superfluid stars with only modified Urca reactions), this leads to a quasi-steady state in which the star radiates thermal photons for a long time, possibly accounting for the ultraviolet radiation observed from the millisecond pulsar J0437-4715. For the first time, we explore the phenomenology of rotochemical heating with direct Urca reactions and uniform and isotropic superfluid energy gaps of different sizes. We first do exploratory work by integrating the thermal and chemical evolution equations numerically for different energy gaps, which uncovers a rich phenomenology of stable and unstable solutions. To understand these, we perform a stability analysis around the quasi-steady state, identifying the characteristic times of growing, decaying, and oscillating solutions. For small gaps, the phenomenology is similar to the previously studied cases, in the sense that the solutions quickly converge to a quasi-steady state. For large gaps ($\gtrsim 0.05$ MeV), these solutions become unstable, leading to a limit-cycle behavior of periodicity $\sim 10^{6-7} yr, in which the star is hot ($T_s\gtrsim 10^5$ K) for a small fraction of the cycle ($\sim 5- 20 %$), and cold for a longer time.

preprint2011arXiv

Revisiting the Flowers-Ruderman instability of magnetic stars

In 1977, Flowers and Ruderman described a perturbation that destabilises a purely dipolar magnetic field in a fluid star. They considered the effect of cutting the star in half along a plane containing the symmetry axis and rotating each half by $90\degr$ in opposite directions, which would cause the energy of the magnetic field in the exterior of the star to be greatly reduced, just as it happens with a pair of aligned magnets. We formally solve for the energy of the external magnetic field and check that it decreases monotonously along the entire rotation. We also describe the instability using perturbation theory, and see that it happens due to the work done by the interaction of the magnetic field with surface currents. Finally, we consider the stabilising effect of adding a toroidal field by studying the potential energy perturbation when the rotation is not done along a sharp cut, but with a continuous displacement field that switches the direction of rotation across a region of small but finite width. Using these results, we estimate the relative strengths of the toroidal and poloidal field needed to make the star stable to this displacement and see that the energy of the toroidal field required for stabilisation is much smaller than the energy of the poloidal field. We also show that, contrary to a common argument, the Flowers-Ruderman instability cannot be applied many times in a row to reduce the external magnetic energy indefinitely.

preprint2010arXiv

Asymptotic, non-linear solutions for ambipolar diffusion in one dimension

We study the effect of the non-linear process of ambipolar diffusion (joint transport of magnetic flux and charged particles relative to neutral particles) on the long-term behavior of a non-uniform magnetic field in a one-dimensional geometry. Our main focus is the dissipation of magnetic energy inside neutron stars(particularly magnetars), but our results have a wider application, particularly to the interstellar medium and the loss of magnetic flux from collapsing molecular cloud cores. Our system is a weakly ionized plasma in which neutral and charged particles can be converted into each other through nuclear beta decays (or ionization-recombination processes). In the "weak-coupling" limit of infrequent inter-particle interactions, the evolution of the magnetic field is controlled by the beta decay rate and can be described by a non-linear partial integro-differential equation. In the opposite, "strong-coupling" regime, the evolution is controlled by the inter-particle collisions and can be modelled through a non-linear diffusion equation. We show numerically that, in both regimes, ambipolar diffusion tends to spread out the magnetic flux, but, contrary to the normal Ohmic diffusion, it produces sharp magnetic field gradients with associated current sheets around those regions where the magnetic field is weak.

preprint2010arXiv

Internal Heating of Old Neutron Stars: Contrasting Different Mechanisms

Context: The standard cooling models of neutron stars predict temperatures $T<10^{4}$ K for ages $t>10^{7}$ yr. However, the likely thermal emission detected from the millisecond pulsar J0437-4715, of spin-down age $t_s \sim 7\times10^9$ yr, implies a temperature $T\sim 10^5$ K. Thus, a heating mechanism needs to be added to the cooling models in order to obtain agreement between theory and observation. Aims: Several internal heating mechanisms could be operating in neutron stars, such as magnetic field decay, dark matter accretion, crust cracking, superfluid vortex creep, and non-equilibrium reactions ("rotochemical heating"). We study these mechanisms in order to establish which could be the dominant source of thermal emission from old pulsars. Methods: We show by simple estimates that magnetic field decay, dark matter accretion, and crust cracking mechanism are unlikely to have a significant effect on old neutron stars. The thermal evolution for the other mechanisms is computed using the code of Fernández and Reisenegger. Given the dependence of the heating mechanisms on the spin-down parameters, we study the thermal evolution for two types of pulsars: young, slowly rotating "classical" pulsars and old, fast rotating millisecond pulsars. Results: We find that magnetic field decay, dark matter accretion, and crust cracking do not produce detectable heating of old pulsars. Rotochemical heating and vortex creep can be important both for classical pulsars and millisecond pulsars. More restrictive upper limits on the surface temperatures of classical pulsars could rule out vortex creep as the main source of thermal emission. Rotochemical heating in classical pulsars is driven by the chemical imbalance built up during their early spin-down, and therefore strongly sensitive to their initial rotation period.

preprint2010arXiv

Internal Heating of Old Neutron Stars: Contrasting Different Mechanisms

The thermal emission detected from the millisecond pulsar J0437-4715 is not explained by standard cooling models of neutron stars without a heating mechanism. We investigated three heating mechanisms controlled by the rotational braking of the pulsar: breaking of the solid crust, superfluid vortex creep, and non-equilibrium reactions ('rotochemical heating'). We find that the crust cracking mechanism does not produce detectable heating. Given the dependence of the heating mechanisms on spin-down parameters, which leads to different temperatures for different pulsars, we study the thermal evolution for two types of pulsars: young, slowly rotating 'classical' pulsars and old, fast rotating millisecond pulsars (MSPs). We find that the rotochemical heating and vortex creep mechanism can be important both for classical pulsars and MSPs.

preprint2010arXiv

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.

preprint2010arXiv

Rotochemical heating in millisecond pulsars: modified Urca reactions with uniform Cooper pairing gaps

Context: 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 nonsuperfluid matter, finding that they reach a quasi-steady configuration in which the rate at which the spin-down modifies the equilibrium concentrations is the same at which neutrino reactions restore the equilibrium. Aims: We describe the thermal effects of Cooper pairing with spatially uniform energy gaps of neutrons Δ_n and protons Δ_p on the rotochemical heating in millisecond pulsars (MSPs) when only modified Urca reactions are allowed. By this, we may determine the amplitude of the superfluid energy gaps for the neutron and protons needed to produce different thermal evolution of MSPs. Results: We find that the chemical imbalances in the star grow up to the threshold value Δ_{thr}= min(Δ_n+ 3Δ_p, 3Δ_n+Δ_p), which is higher than the quasi-steady state achieved in absence of superfluidity. Therefore, the superfluid MSPs will take longer to reach the quasi-steady state than their nonsuperfluid counterparts, and they will have a higher a luminosity in this state, given by L_γ~ (1-4) 10^{32}Δ_{thr}/MeV \dot{P}_{-20}/P_{ms}^3 erg s^-1. We can explain the UV emission of the PSR J0437-4715 for 0.05 MeV<Δ_{thr}<0.45 MeV.

preprint2010arXiv

Rotochemical heating with a density-dependent superfluid energy gap in neutron stars

When a rotating neutron star loses angular momentum, the reduction of the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium, inducing reactions that release heat (rotochemical heating). This effect has previously been studied by Fernández and Reisenegger for neutron stars of non-superfluid matter and by Petrovich and Reisenegger for superfluid matter, finding that the system in both cases reaches a quasi-steady state, corresponding to a partial equilibration between compression, due to the loss of angular momentum, and reactions that try to restore the equilibrium. However, Petrovich and Reisenegger assumes a constant value of the superfluid energy gap, whereas theoretical models predict density-dependent gap amplitudes, and therefore gaps that depend on the location in the star. In this work, we try to discriminate between several proposed gap models, comparing predicted surface temperatures to the value measured for the nearest millisecond pulsar, J0437-4715.

preprint2009arXiv

Future Evolution of Bound Superclusters in an Accelerating Universe

The evolution of marginally bound supercluster-like objects in an accelerating LambdaCDM Universe is followed, by means of cosmological simulations, from the present time to an expansion factor a = 100. The objects are identified on the basis of the binding density criterion introduced by Dunner et al. (2006). superclusters are identified with the ones whose mass M > 10^15 M_sun/h, the most massive one with M ~ 8x10^15 M_sun/h, comparable to the Shapley supercluster. The spatial distribution of the superclusters remains essentially the same after the present epoch, reflecting the halting growth of the Cosmic Web as Lambda gets to dominate the expansion of the Universe. The same trend can be seen in the stagnation of the development of the mass function of virialized haloes and bound objects. The situation is considerably different when looking at the internal evolution, quantified in terms of their shape, compactness and density profile, and substructure in terms of their multiplicity function. We find a continuing evolution from a wide range of triaxial shapes at a = 1 to almost perfect spherical shapes at a = 100. We also find a systematic trend towards a higher concentration. Meanwhile, we see their substructure gradually disappearing, as the surrounding subclumps fall in and merge to form one coherent, virialized system.

preprint2009arXiv

Stable magnetic equilibria and their evolution in the upper main sequence, white dwarfs, and neutron stars

[abbreviated] Long-lived, large-scale magnetic field configurations exist in upper main sequence, white dwarf, and neutron stars. Externally, these fields have a strong dipolar component, while their internal structure and evolution are uncertain, but highly relevant for several problems in stellar and high-energy astrophysics. We discuss the main properties expected for the stable magnetic configurations in these stars from physical arguments, and how these properties may determine the modes of decay of these configurations. Stable magneto-hydrostatic equilibria are likely to exist in stars whenever the matter in their interior is stably stratified (not barotropic). These equilibria are not force-free and not required to satisfy the Grad-Shafranov equation, but they do involve both toroidal and poloidal field components. We argue that the main mode of decay for these configurations are processes that lift the constraints set by stable stratification, such as heat diffusion in main-sequence envelopes and white dwarfs, and beta decays or particle diffusion in neutron stars. In the former, heat diffusion is not fast enough to make these equilibria evolve over the stellar lifetime. In neutron stars, a strong enough field might decay by overcoming the compositional stratification through beta decays (at the highest field strengths) or through ambipolar diffusion (for somewhat weaker fields). These processes convert magnetic energy to thermal energy, and they occur at significant rates only once the latter is smaller than the former, and therefore they substantially delay the cooling of the neutron star, while slowly decreasing its magnetic energy.

preprint2005arXiv

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.

preprint1995arXiv

THE GUNN-PETERSON EFFECT FROM UNDERDENSE REGIONS IN A PHOTOIONIZED INTERGALACTIC MEDIUM

We use the Zel'dovich approximation and another analytical approximation to calculate the evolution under gravitational instability of the underdense regions of a photoionized intergalactic medium (IGM). We find that over most of the spectrum of a quasar, the optical depth to $\lya$ scattering originates from gas in underdense regions, or voids. This causes the ratio of the median Gunn-Peterson (GP) absorption to the value for a uniform medium containing all the baryons in the universe to be very small, and to decrease as gravitational collapse proceeds. We apply our calculations to the observations of the intensity distribution in a $z=4.11$ quasar by Webb and coworkers. We show that if $\lya$ clouds arise from gravitational collapse, their observations must be interpreted as the first detection of the (fluctuating) GP effect, with a median value $τ_{GP}\simeq 0.06$ at $z=4$. For typical low-density ($Ω\sim 0.4$) cosmological models, this is consistent with the predicted baryon density from primordial nucleosynthesis, and the intensity of the ionizing background derived from the proximity effect. From the numerical simulations of Cen \etal, such models also predict correctly the number of $\lya$ absorption lines observed.