Presupernova Evolution of Differentially Rotating Massive Stars Including Magnetic Fields

As a massive star evolves through multiple stages of nuclear burning on its way to becoming a supernova, a complex, differentially rotating structure is set up. Angular momentum is transported by a variety of classic instabilities, and also by magnetic torques from fields generated by the differential rotation. We present the first stellar evolution calculations to follow the evolution of rotating massive stars including, at least approximately, all these effects, magnetic and non-magnetic, from the zero-age main sequence until the onset of iron-core collapse. The evolution and action of the magnetic fields is as described by Spruit 2002 and a range of uncertain parameters is explored. In general, we find that magnetic torques decrease the final rotation rate of the collapsing iron core by about a factor of 30 to 50 when compared with the non-magnetic counterparts. Angular momentum in that part of the presupernova star destined to become a neutron star is an increasing function of main sequence mass. That is, pulsars derived from more massive stars will rotate faster and rotation will play a more dominant role in the star's explosion. The final angular momentum of the core is determined - to within a factor of two - by the time the star ignites carbon burning. For the lighter stars studied, around 15 solar masses, we predict pulsar periods at birth near 15 ms, though a factor of two range is easily tolerated by the uncertainties. Several mechanisms for additional braking in a young neutron star, especially by fall back, are also explored.Comment: 32 pages, 3 figures (8 eps files), submitted to Ap

Birth kicks as the origin of pulsar rotation

Radio pulsars are thought to born with spin periods of 0.02–0.5 s and space velocities of 100–1,000 kms^(-1), and they are inferred to have initial dipole magnetic fields of 10^(11)–10^(13) G. The average space velocity of their progenitor stars is less than 15 kms^(-1), which means that pulsars must receive a substantial ‘kick’ at birth. Here we propose that the birth characteristics of pulsars have a simple physical connection with each other. Magnetic fields maintained by differential rotation between the core and envelope of the progenitor would keep the whole star in a state of approximately uniform rotation until 10 years before the explosion. Such a slowly rotating core has 1,000 times less angular momentum than required to explain the rotation of pulsars. The specific physical process that ‘kicks’ the neutron star at birth has not been identified, but unless its force is exerted exactly head-on it will also cause the neutron star to rotate. We identify this process as the origin of the spin of pulsars. Such kicks may cause a correlation between the velocity and spin vectors of pulsars. We predict that many neutron stars are born with periods longer than 2 s, and never become radio pulsars

Presupernova Evolution of Rotating Massive Stars and the Rotation Rate of Pulsars

Rotation in massive stars has been studied on the main sequence and during helium burning for decades, but only recently have realistic numerical simulations followed the transport of angular momentum that occurs during more advanced stages of evolution. The results affect such interesting issues as whether rotation is important to the explosion mechanism, whether supernovae are strong sources of gravitational radiation, the star's nucleosynthesis, and the initial rotation rate of neutron stars and black holes. We find that when only hydrodynamic instabilities (shear, Eddington-Sweet, etc.) are included in the calculation, one obtains neutron stars spinning at close to critical rotation at their surface -- or even formally in excess of critical. When recent estimates of magnetic torques (Spruit 2002) are added, however, the evolved cores spin about an order of magnitude slower. This is still more angular momentum than observed in young pulsars, but too slow for the collapsar model for gamma-ray bursts.Comment: 10 pages, 2 figures, to appear in Proc. IAU 215 "Stellar Rotation

The kink-type instability of toroidal stellar magnetic fields with thermal diffusion

The stability of toroidal magnetic fields in rotating radiative stellar zones is studied for realistic values of both the Prandtl numbers. The two considered models for the magnetic geometry represent fields with odd and even symmetry with respect to the equator. In the linear theory in Boussinesq approximation the resulting complex eigenfrequency (including growth rate and drift rate) are calculated for a given radial wavenumber of a nonaxisymmetric perturbation with m=1. The ratio of the Alfven frequency, \Omega_A, to the rate of the basic rotation, \Omega, controls the eigenfrequency of the solution. For strong fields with \Omega_A > \Omega the solutions do not feel the thermal diffusion. The growth rate runs with \Omega_A and the drift rate is close to -\Omega so that the magnetic pattern will rest in the laboratory system. For weaker fields with \Omega_A < \Omega the growth rate strongly depends on the thermal conductivity. For fields with dipolar parity and for typical values of the heat conductivity the resulting very small growth rates are almost identical with those for vanishing gravity. For fields with dipolar symmetry the differential rotation of any stellar radiative zone (like the solar tachocline) is shown as basically stabilizing the instability independent of the sign of the shear. Finally, the current-driven kink-type instability of a toroidal background field is proposed as a model for the magnetism of Ap stars. The recent observation of a lower magnetic field treshold of about 300 Gauss for Ap stars is understood as corresponding to the minimum magnetic field producing the instability.Comment: 11 pages, 7 figures, acc. for publicatio

Gas-dynamic shock heating of post-flare loops due to retraction following localized, impulsive reconnection

We present a novel model in which shortening of a magnetic flux tube following localized, three-dimensional reconnection generates strong gas-dynamic shocks around its apex. The shortening releases magnetic energy by progressing away from the reconnection site at the Alfven speed. This launches inward flows along the field lines whose collision creates a pair of gas-dynamic shocks. The shocks raise both the mass density and temperature inside the newly shortened flux tube. Reconnecting field lines whose initial directions differ by more that 100 degrees can produce a concentrated knot of plasma hotter that 20 MK, consistent with observations. In spite of these high temperatures, the shocks convert less than 10% of the liberated magnetic energy into heat - the rest remains as kinetic energy of bulk motion. These gas-dynamic shocks arise only when the reconnection is impulsive and localized in all three dimensions; they are distinct from the slow magnetosonic shocks of the Petschek steady-state reconnection model

Magnetically-dominated jets inside collapsing stars as a model for gamma-ray bursts and supernova explosions

It has been suggested that magnetic fields play a dynamically-important role in core-collapse explosions of massive stars. In particular, they may be important in the collapsar scenario for gamma-ray bursts (GRB), where the central engine is a hyper-accreting black hole or a millisecond magnetar. The present paper is focussed on the magnetar scenario, with a specific emphasis on the interaction of the magnetar magnetosphere with the infalling stellar envelope. First, the ``Pulsar-in-a-Cavity'' problem is introduced as a paradigm for a magnetar inside a collapsing star. The basic set-up of this fundamental plasma-physics problem is described, outlining its main features, and simple estimates are derived for the evolution of the magnetic field. In the context of a collapsing star, it is proposed that, at first, the ram pressure of the infalling plasma acts to confine the magnetosphere, enabling a gradual build-up of the magnetic pressure. At some point, the growing magnetic pressure overtakes the (decreasing) ram pressure of the gas, resulting in a magnetically-driven explosion. The explosion should be highly anisotropic, as the hoop-stress of the toroidal field, confined by the surrounding stellar matter, collimates the magnetically-dominated outflow into two beamed magnetic-tower jets. This creates a clean narrow channel for the escape of energy from the central engine through the star, as required for GRBs. In addition, the delayed onset of the collimated-explosion phase can explain the production of large quantities of Nickel-56, as suggested by the GRB-Supernova connection. Finally, the prospects for numerical simulations of this scenario are discussed.Comment: Invited paper in the "Physics of Plasmas" (May 2007 special issue), based on an invited talk at the 48th Annual Meeting of the APS Division of Plasma Physics (Oct. 30 - Nov. 3, 2006, Philadelphia, PA); 24 pages, 7 figure

Numerical simulations of the Accretion-Ejection Instability in magnetised accretion disks

The Accretion-Ejection Instability (AEI) described by Tagger & Pellat (1999) is explored numerically using a global 2d model of the inner region of a magnetised accretion disk. The disk is initially currentless but threaded by a vertical magnetic field created by external currents, and frozen in the flow. In agreement with the theory a spiral instability, similar in many ways to those observed in self-gravitating disks, develops when the magnetic field is, within a factor of a few, at equipartition with the disk thermal pressure. Perturbations in the flow build up currents and create a perturbed magnetic field within the disk. The present non-linear simulations give good evidence that such an instability can occur in the inner region of accretion disks, and generate accretion of gas and vertical magnetic flux toward the central object, if the equilibrium radial profiles of density and magnetic flux exceed a critical threshold.Comment: single tar file with GIF figure

The Origin of Solar Activity in the Tachocline

Solar active regions, produced by the emergence of tubes of strong magnetic field in the photosphere, are restricted to within 35 degrees of the solar equator. The nature of the dynamo processes that create and renew these fields, and are therefore responsible for solar magnetic phenomena, are not well understood. We analyze the magneto-rotational stability of the solar tachocline for general field geometry. This thin region of strong radial and latitudinal differential rotation, between the radiative and convective zones, is unstable at latitudes above 37 degrees, yet is stable closer to the equator. We propose that small-scale magneto-rotational turbulence prevents coherent magnetic dynamo action in the tachocline except in the vicinity of the equator, thus explaining the latitudinal restriction of active regions. Tying the magnetic dynamo to the tachocline elucidates the physical conditions and processes relevant to solar magnetism.Comment: 10 pages, 1 figure, accepted for publication in ApJ