This paper identifies constraints on the growth of a small planetary core (0.3 M$_{\oplus}$) that accretes millimeter-sized pebbles from a gaseous disk. We construct time-dependent spherical envelope models that capture physical processes not included in existing global hydrodynamic simulations, including particle size evolution, dust transport, and realistic gas equations of state. We assume a low enough disk density that pebbles are marginally coupled to the gas and are trapped efficiently near the core Bondi radius. Pebbles then drift rapidly enough to experience erosion by sandblasting, mutual collisions, and sublimation of water ice. We find that pebble fragmentation is more efficient than dust re-sticking. Therefore the large pebble accretion rate $\dot{M}_p$ needed to build a core of mass > M$_{\oplus}$ leads to a high envelope metallicity and grain opacity. Above $\dot{M}_p$ $\sim 10^{-7}$ M$_{\oplus}$/yr, and without other luminosity sources, convective motions expand near the Bondi radius. The warm, dusty, and turbulent envelope buffers the inward drift of pebble debris: given a turbulent concentration factor $f_{turb} \geq$ 1 near the lower convective boundary, the core growth rate is limited to $1 \times 10^{-7} f_{turb}$ M$_{\oplus}$/yr and the e-folding time 3$f_{turb}^{-1}$ Myr. The remainder of the solid debris is expelled as highly processed silicates. Pebble ice never reaches the core, and the envelope contains comparable amounts of H2/He and metals. We interpret our results using simpler steady models and semi-analytical estimates. Future global simulations incorporating the processes modelled here are needed to understand the influence of rotation and vertical disk structure.

We evaluate several basic electrodynamic processes as modified by the presence of a very strong magnetic field, exceeding $B_{\rm Q} \equiv m^2/e = 4.4\times 10^{13}$ G. These results are needed to build models of dissipative phenomena outside magnetars and some other neutron stars. Differential and total cross sections and rates are presented for electron-photon scattering, the annihilation of an electron-positron pair into two photons, the inverse process of two-photon pair creation, and single-photon pair creation into the lowest Landau state. The relative importance of these interactions changes as the background magnetic field grows in strength. The particle phase space relevant for a given process may be restricted by single-photon pair creation, which also opens up efficient channels for pair multiplication, e.g. in combination with scattering. Our results are presented in the form of compact formulae that allow for relativistic electron (positron) motion, in the regime where Landau excitations can be neglected (corresponding to $10^3B_{\rm Q} \gg B \gg B_{\rm Q}$ for moderately relativistic motion along the magnetic field). Where a direct comparison is possible, our results are tested against earlier calculations, and a brief astrophysical context is provided. A companion paper considers electron-positron scattering, scattering of electrons and positrons by ions, and relativistic electron-ion bremsstrahlung.

This paper continues the approach of Kostenko \& Thompson to calculating quantum electrodynamic processes in the ultrastrong magnetic field near some neutron stars, such as magnetars or merging binary neutron stars. Here we consider electron-positron scattering, the Coulomb scattering of electrons and positrons off ions, and relativistic $e^\pm$-ion bremsstrahlung. The evaluation of differential and total cross sections simplifies considerably when the magnetic field lies in the range $10^3B_{\rm Q} \gg B \gg B_{\rm Q}$, where $B_{\rm Q} \equiv m^2/e = 4.4\times 10^{13}$ G. Then, relativistic motion of $e^\pm$ is possible even when restricted to the lowest Landau state. Accurate results for differential and total cross sections are obtained by truncating the sum over intermediate-state Landau levels and otherwise disregarding terms inversely proportional to the magnetic field, which are complicated enough to have inhibited previous attempts to calculate magnetic electron-positron scattering and relativistic bremsstrahlung. A quantitative account is made of the effects of Debye screening.

The aerosols present in the atmosphere of the Saturn s moon Titan are of particular planetary science interest and several spacecraft missions already allowed to gather spectroscopic data. Titan haze s analogs, so-called tholins, were produced on earth to push forward the comprehension of their formation and properties. In this study, this highly complex mixture was analyzed here for the first time by trapped ion mobility spectrometry coupled to ultra-high resolution mass spectrometry (FTICR MS). Electrospray ionization revealed the characteristic CHNx-class components, with CHN5-6 and DBE 6-7 most abundant. Deploying specialized visualization, enabled by accurate mass measurements and elemental composition assignments, the adapted Kendrick mass defect analysis highlights the C2H3N homolog series, whereas the nitrogen-modified van Krevelen diagram exhibits a clear trend towards H/C 1.5 and N/C 0.5. More interestingly, the representation of m/z versus collision cross section (CCS) allowed hypothesizing a ramified N-PAH structural motif. State-of-the-art IMS is currently not able to resolve the isomeric continuum of ultra-complex mixtures; thus peak parameter other than the CCS-value are explored. As such, analyzing the mobility peak width versus m/z shows a linear increase in isomeric diversity between m/z 170 and 350 and a near plateau in diversity at higher m/z for the N-PAH-like structure. Due to the high complexity of the sample, these structural insights are only to be revealed by TIMS-FTICR MS.

The electromagnetic signature of a point explosion near a Kerr black hole (BH) is evaluated. The first repetitions produced by gravitational lensing are not periodic in time; periodicity emerges only as the result of multiple circuits of the prograde and retrograde light rings and is accompanied by exponential dimming. Gravitational focusing creates a sequence of concentrated caustic features and biases the detection of a repeating source toward alignment of the BH spin with the plane of the sky. We consider the polarization pattern in the case of emission by the Lorentz upboosting and reflection of a magnetic field near the explosion site. Then the polarized fraction of the detected pulse approaches unity, and rays propagating near the equatorial plane maintain a consistent polarization direction. Near a slowly accreting supermassive BH (SMBH), additional repetitions are caused by reflection off annular fragments of an orbiting disk that has passed through an ionization instability. These results are applied to the repeating fast radio burst (FRB) source 121102, giving a concrete and predictive example of how FRB detectability may be biased by lensing. A gravitational lensing delay of 10-30 s, and reflection delay up to $\sim 10^4$ s, are found for emission near the innermost stable circular orbit of a $3\times 10^5\,M_\odot$ SMBH; these effects combine to produce interesting correlations between delay time and burst fluence. A similar repetitive pulse envelope could be seen in the gravitational wave signal produced by a collision between compact stars near a SMBH.

The internal rotation and magnetism of evolved massive stars are considered in response to i) the inward pumping of angular momentum through deep and slowly rotating convective layers; and ii) the winding up of a helical magnetic field in radiative layers. Field winding can transport angular momentum effectively even when the toroidal field is limited by kinking. Magnetic helicity is pumped into a growing radiative layer from an adjacent convective envelope (or core). The receding convective envelope that forms during the early accretion phase of a massive star is the dominant source of helicity in its core, yielding a $\sim 10^{13}$ G polar magnetic field in a collapsed neutron star (NS) remnant. Using MESA models of various masses, we find that the NS rotation varies significantly, from $P_{\rm NS} \sim 0.1-1$ s in a 13$\,M_\odot$ model to $P_{\rm NS} \sim 2$ ms in a $25\,M_\odot$ model with an extended core. Stronger inward pumping of angular momentum is found in more massive stars, due to the growing thickness of the convective shells that form during the later stages of thermonuclear burning. On the other hand, stars that lose enough mass to form blue supergiants in isolation end up as very slow rotators. The tidal spin-up of a 40$\,M_\odot$ star by a massive binary companion is found to dramatically increase the spin of the remnant black hole, allowing a rotationally supported torus to form during the collapse. The implications for post-collapse decay or amplification of the magnetic field are also considered.

This paper considers electromagnetic transients of a modest total energy (${\cal E} \gtrsim 10^{40}$ erg) and small initial size (${\cal R} \gtrsim 10^{-1}$ cm). They could be produced during collisions between relativistic field structures (e.g. macroscopic magnetic dipoles) that formed around, or before, cosmic electroweak symmetry breaking. The outflowing energy has a dominant electromagnetic component; a subdominant thermal component (temperature $> 1$ GeV) supplies inertia in the form of residual $e^\pm$. A thin shell forms that expands subluminally, attaining a Lorentz factor $\sim 10^{6-7}$ before decelerating. Drag is supplied by the reflection of an ambient magnetic field, and by deflection of ambient free electrons. Emission of low-frequency (GHz-THz) superluminal waves takes place through three channels: i) reflection of the ambient magnetic field; ii) direct linear conversion of the embedded magnetic field into a superluminal mode; and iii) excitation outside the shell by corrugation of its surface. The escaping electromagnetic pulse is very narrow (a few wavelengths) and so the width of the detected transient is dominated by propagation effects. GHz radio transients are emitted from i) the dark matter halos of galaxies and ii) the near-horizon regions of supermassive black holes that formed by direct gas collapse and now accrete slowly. Brighter and much narrower 0.01-1 THz pulses are predicted at a rate at least comparable to fast radio bursts, experiencing weaker scattering and absorption. The same explosions also accelerate protons up to $\sim 10^{19}$ eV and heavier nuclei up to $10^{20-21}$ eV.

Macroscopic magnetic dipoles are considered as cosmic dark matter. Permanent magnetism in relativistic field structures can involve some form of superconductivity, one example being current-carrying string loops (`springs') with vanishing net tension. We derive the cross section for free classical dipoles to collide, finding it depends weakly on orientation when mutual precession is rapid. The collision rate of `spring' loops with tension ${\cal T} \sim 10^{-8}c^4/G$ in galactic halos approaches the measured rate of fast radio bursts (FRBs) if the loops comprise most of the dark matter. A large superconducting dipole (LSD) with mass $\sim 10^{20}$ g and size $\sim 1$ mm will form a $\sim 100$ km magnetosphere moving through interstellar plasma. Although hydromagnetic drag is generally weak, it is strong enough to capture some LSDs into long-lived rings orbiting supermassive black holes (SMBHs) that form by the direct collapse of massive gas clouds. Repeated collisions near young SMBHs could dominate the global collision rate, thereby broadening the dipole mass spectrum. Colliding LSDs produce tiny, hot electromagnetic explosions. The accompanying paper shows that these explosions couple effectively to propagating low-frequency electromagnetic modes, with output peaking at 0.01-1 THz. We describe several constraints on, and predictions of, LSDs as cosmic dark matter. The shock formed by an infalling LSD triggers self-sustained thermonuclear burning in a C/O (ONeMg) white dwarf (WD) of mass $\gtrsim 1\,M_\odot$ ($1.3\,M_\odot$). The spark is generally located well off the center of the WD. The rate of LSD-induced explosions matches the observed rate of Type Ia supernovae.

This paper considers the yielding response of a neutron star crust to smooth, unbalanced Maxwell stresses imposed at the core-crust boundary, and the coupling of the dynamic crust to the external magnetic field. Stress buildup and yielding in a magnetar crust is a global phenomenon: an elastic distortion radiating from one plastically deforming zone is shown to dramatically increase the creep rate in distant zones. Runaway creep to dynamical rates is shown to be possible, being enhanced by in situ heating and suppressed by thermal conduction and shearing of an embedded magnetic field. A global and time-dependent model of elastic, plastic, magnetic, and thermal evolution is developed. Fault-like structures develop naturally, and a range of outburst timescales is observed. Transient events with time profiles similar to giant magnetar flares (millisecond rise, $\sim$ 0.1 s duration, and decaying power-law tails) result from runaway creep that starts in localized sub-km-sized patches and spreads across the crust. A one-dimensional model of stress relaxation in the vertically stratified crust shows that a modest increase in applied stress allows embedded magnetic shear to escape the star over $\sim$ 3-10 ms, dissipating greater energy if the exterior field is already sheared. Several such zones coupled to each other naturally yield a burst of duration $\sim$ 0.1 s, as is observed over a wide range of burst energies. The collective interaction of many plastic zones forces an overstability of global elastic modes of the crust, consistent with QPO activity extending over $\gtrsim$ 100 s. Giant flares probably involve sudden meltdown in localized zones, with high-frequency ($\gg$ 100 Hz) QPOs corresponding to standing Alfv\'en waves within these zones.

This paper studies the response of a thin accretion disk to an external radial magnetic field. Our focus is on protoplanetary disks (PPDs), which are exposed during their later evolution to an intense, magnetized wind from the central star. A radial magnetic field is mixed into a thin surface layer, is wound up by the disk shear, and is pushed downward by a combination of turbulent mixing and ambipolar and Ohmic drift. The toroidal field reaches much greater strengths than the seed vertical field that is usually invoked in PPD models, even becoming superthermal. Linear stability analysis indicates that the disk experiences the magnetorotational instability (MRI) at a higher magnetization than a vertically magnetized disk when both the effects of ambipolar and Hall drift are taken into account. Steady vertical profiles of density and magnetic field are obtained at several radii between 0.06 and 1 AU in response to a wind magnetic field $B_r \sim (10^{-4}$-$10^{-2})(r/{\rm AU})^{-2}$ G. Careful attention is given to the radial and vertical ionization structure resulting from irradiation by stellar X-rays. The disk is more strongly magnetized closer to the star, where it can support a higher rate of mass transfer. As a result, the inner $\sim 1$ AU of a PPD is found to evolve toward lower surface density. Mass transfer rates around $10^{-8}\,M_\odot$ yr$^{-1}$ are obtained under conservative assumptions about the MRI-generated stress. The evolution of the disk, and the implications for planet migration, are investigated in the accompanying paper.

We consider the inner $\sim$ AU of a protoplanetary disk (PPD), at a stage where angular momentum transport is driven by the mixing of a radial magnetic field into the disk from a T-Tauri wind. Because the radial profile of the imposed magnetic field is well constrained, a deterministic calculation of the disk mass flow becomes possible. The vertical disk profiles obtained in Paper I imply a stronger magnetization in the inner disk, faster accretion, and a secular depletion of the disk material. Inward transport of solids allows the disk to maintain a broad optical absorption layer even when the grain abundance becomes too small to suppress its ionization. Thus a PPD may show a strong middle-to-near infrared spectral excess even while its mass profile departs radically from the minimum-mass solar nebula. The disk surface density is buffered at $\sim 30$ g cm$^{-2}$: below this, X-rays trigger strong enough magnetorotational turbulence at the midplane to loft mm-cm sized particles high in the disk, followed by catastrophic fragmentation. A sharp density gradient bounds the inner depleted disk, and propagates outward to $\sim 1$-2 AU over a few Myr. Earth-mass planets migrate through the inner disk over a similar timescale, whereas the migration of Jupiters is limited by the supply of gas. Gas-mediated migration must stall outside 0.04 AU, where silicates are sublimated and the disk shifts to a much lower column. A transition disk emerges when the dust/gas ratio in the MRI-active layer falls below $X_d \sim 10^{-6}(a_d/\mu{\rm m})$, where $a_d$ is the grain size.

The magnetism and rotation of white dwarf (WD) stars are investigated in relation to a hydromagnetic dynamo operating in the progenitor during shell burning phases. The downward pumping of angular momentum in the convective envelope, in combination with the absorption of a planet or tidal spin-up from a binary companion, can trigger strong dynamo action near the core-envelope boundary. Several arguments point to the outer core as the source for a magnetic field in the WD remnant: the outer third of a $\sim 0.55\,M_\odot$ WD is processed during the shell burning phase(s) of the progenitor; the escape of magnetic helicity through the envelope mediates the growth of (compensating) helicity in the core, as is needed to maintain a stable magnetic field in the remnant; and the intense radiation flux at the core boundary facilitates magnetic buoyancy within a relatively thick tachocline layer. The helicity flux into the growing core is driven by a dynamical imbalance with a latitude-dependent rotational stress. The magnetic field deposited in an isolated massive WD is concentrated in an outer shell of mass $\lesssim 0.1\,M_\odot$ and can reach $\sim 10\,$MG. A buried toroidal field experiences moderate ohmic decay above an age $\sim 0.3$ Gyr, which may lead to growth or decay of the external magnetic field. The final WD spin period is related to a critical spin rate below which magnetic activity shuts off, and core and envelope decouple; it generally sits in the range of hours to days. WD periods ranging up to a year are possible if the envelope re-expands following a late thermal pulse.

The internal rotation of post-main sequence stars is investigated, in response to the convective pumping of angular momentum toward the stellar core, combined with a tight magnetic coupling between core and envelope. The spin evolution is calculated using model stars of initial mass 1, 1.5 and $5\,M_\odot$, taking into account mass loss on the giant branches. We also include the deposition of orbital angular momentum from a sub-stellar companion, as influenced by tidal drag along with the excitation of orbital eccentricity by a fluctuating gravitational quadrupole moment. A range of angular velocity profiles $\Omega(r)$ is considered in the envelope, extending from solid rotation to constant specific angular momentum. We focus on the back reaction of the Coriolis force, and the threshold for dynamo action in the inner envelope. Quantitative agreement with measurements of core rotation in subgiants and post-He core flash stars by Kepler is obtained with a two-layer angular velocity profile: uniform specific angular momentum where the Coriolis parameter ${\rm Co} \equiv \Omega \tau_{\rm con} \lesssim 1$ (here $\tau_{\rm con}$ is the convective time); and $\Omega(r) \propto r^{-1}$ where ${\rm Co} \gtrsim 1$. The inner profile is interpreted in terms of a balance between the Coriolis force and angular pressure gradients driven by radially extended convective plumes. Inward angular momentum pumping reduces the surface rotation of subgiants, and the need for a rejuvenated magnetic wind torque. The co-evolution of internal magnetic fields and rotation is considered in Kissin & Thompson, along with the breaking of the rotational coupling between core and envelope due to heavy mass loss.

Chondrules are millimeter-sized silicate spherules ubiquitous in primitive meteorites, but whose origin remains mysterious. One of the main proposed mechanisms for producing them is melting of solids in shock waves in the gaseous protoplanetary disk. However, evidence is mounting that chondrule-forming regions were enriched in solids well above solar abundances. Given the high velocities involved in shock models destructive collisions would be expected between differently sized grains after passage of the shock front as a result of differential drag. We investigate the probability and outcome of collisions of particles behind a 1D shock using analytic methods as well as a full integration of the coupled mass, momentum, energy and radiation equations. Destruction of protochondrules seems unavoidable for solid/gas ratios $\epsilon \gtrsim 0.1$, and possibly even for solar abundances because of "sandblasting" by finer dust. A flow with $\epsilon \gtrsim 10$ requires much smaller shock velocities ($\sim 2$ vs 8 km s$^{-1}$) in order to achieve chondrule-melting temperatures, and radiation trapping allows slow cooling of the shocked fragments. Initial destruction would still be extensive; although re-assembly of mm-sized particles would naturally occur by grain sticking afterward, the compositional heterogeneity of chondrules may be difficult to reproduce. We finally note that solids passing through small-scale bow shocks around few-km-sized planetesimals might experience partial melting and yet escape fragmentation.

Gamma-ray bursts (GRBs) show a dramatic pulse structure that requires bulk relativistic motion, but whose physical origin has remained murky. We focus on a hot, magnetized jet that is emitted by a black hole and interacts with a confining medium. Strongly relativistic expansion of the magnetic field, as limited by a corrugation instability, may commence only after it forms a thin shell. Then the observed $T_{90}$ burst duration is dominated by the curvature delay, and null periods arise from angular inhomogeneities, not the duty cycle of the engine. We associate the $O(1)$ s timescale observed in the pulse width distribution of long GRBs with the collapse of the central 2.5-3$M_\odot$ of a massive stellar core. A fraction of the baryons are shown to be embedded in the magnetized outflow by the hyper-Eddington radiation flux; they strongly disturb the magnetic field after the compactness drops below $\sim 4\times 10^3(Y_e/0.5)^{-1}$. The high-energy photons so created have a compressed pulse structure. Delayed breakout of magnetic field from heavier baryon shells is also a promising approach to X-ray flares. In the second part of the paper, we calculate the imprint of an expanding, scattering photosphere on pulse evolution. Two models for generating the high-energy spectral tail are contrasted: i) pair breakdown due to reheating of an optically thin pair plasma embedded in a thermal radiation field; and ii) continuous heating extending from large to small scattering depth. The second model is strongly inconsistent with the observed hard-to-soft evolution in GRB pulses. The first shows some quantitative differences if the emission is purely spherical, but we show that finite shell width, mild departures from spherical curvature, and latitudinal Lorentz factor gradients have interesting effects.

We consider delayed, volumetric heating in a magnetized outflow that has broken out of a confining medium and expanded to a high Lorentz factor ($\Gamma \sim 10^2-10^3$) and low optical depth to scattering ($\tau_{\rm T} \sim 10^{-3}-10^{-2}$). The energy flux at breakout is dominated by the magnetic field, with a modest contribution from quasi-thermal gamma rays whose spectrum was calculated in Paper I. We focus on the case of extreme baryon depletion in the magnetized material, but allow for a separate baryonic component that is entrained from a confining medium. Dissipation is driven by relativistic motion between these two components, which develops once the photon compactness drops below $ 4\times 10^3(Y_e/0.5)^{-1}$. We first calculate the acceleration of the magnetized component following breakout, showing that embedded MHD turbulence provides significant inertia, the neglect of which leads to unrealistically high estimates of flow Lorentz factor. After re-heating begins, the pair and photon distributions are evolved self-consistently using a one-zone kinetic code that incorporates an exact treatment of Compton scattering, pair production and annihilation, and Coulomb scattering. Heating leads to a surge in pair creation, and the scattering depth saturates at $\tau_{\rm T} \sim$ 1-4. The plasma maintains a very low ratio of particle to magnetic pressure, and can support strong anisotropy in the charged particle distribution, with cooling dominated by Compton scattering. High-energy power-law spectra with photon indices in the range observed in GRBs ($-3 < \beta < -3/2$) are obtained by varying the ratio of heat input to the seed energy in quasi-thermal photons. We contrast our results with those for continuous heating across an expanding photosphere, and show that the latter model produces soft-hard evolution that is inconsistent with observations of GRBs.

The evolution of a dilute electron-positron fireball is calculated in the regime of strong magnetization and very high compactness (l ~10^3-10^8). Heating is applied at a low effective temperature (< 25 keV), and the fireball is allowed to expand, so that the formation of a black-body spectral distribution is inhibited by pair annihilation. The diffusion equation for Compton scattering is coupled to a single-temperature pair gas and an exact (trans-relativistic) cyclo-synchrotron photon source. We find that the photon spectrum develops a quasi-thermal peak, with a power-law slope below it that is characteristic of gamma-ray bursts. The formation of a thermal high-frequency spectrum is checked using the full kinetic equations at l ~ 10^3. These results have several implications for the central engine of GRBs, and the mechanism of energy transport. 1. Baryon rest mass carries less than ~ 10^{-5} of the energy flux at jet breakout inside ~ 10^{12} cm from the engine, with most carried by the magnetic field. 2. This degree of baryon purity points to the presence of an event horizon in the engine, and neutrons play a negligible role in the prompt emission mechanism. 3. X-ray flashes are emitted by outflows carrying enough baryons that the photosphere is pair-depleted, which we show results in faster thermalization. 4. The relation between observed peak frequency and burst luminosity is bounded below by the observed Amati et al. relation if jet Lorentz factor ~ 1/(opening angle) at breakout. 5. Stellar models are used to demonstrate an inconsistency between the highest observed GRB energies, and a hydrodynamic nozzle: magnetic collimation is required. 6. The magnetized pair gas is dilute enough that high-frequency Alfven waves may become charge starved. Finally, we suggest that limitations on magnetic reconnection from plasma collisionality have been overestimated.

Relativistically expanding sources of X-rays and gamma-rays cover an enormous range of (central) compactness and Lorentz factor. The underlying physics is discussed, with an emphasis on how the dominant dissipative mode and the emergent spectrum depend on these parameters. Photons advected outward from high optical depth are a potentially important source of Compton seeds. Their characteristic energy is bounded below by ~1 MeV in pair-loaded outflows of relatively low compactness, and remains near ~1 MeV at very high compactness and low matter loading. This is compared with the characteristic energy of O(1) MeV observed in the rest frame spectra of many sources, including gamma-ray bursts, OSSE jet sources, MeV Blazars, and the intense initial 0.1 s pulse of the March 5 event. Additional topics discussed include the feedback of pair creation on electron heating and the formation of non-thermal spectra, their effectiveness at shielding the dissipative zone from ambient photons, direct Compton damping of irregularities in the outflow, the relative importance of various soft photon sources, and the softening of the emergent spectrum that results from heavy matter loading. The implications of this work for X-ray and optical afterglow from GRB's are briefly considered. Direct synchrotron emission behind the forward shock is inhibited by the extremely low energy density of the ambient magnetic field. Mildly relativistic ejecta off axis from the main gamma-ray emitting cone become optically thin to scattering on a timescale of ~1 day (E/10^{52} erg)^{1/2}, and can be a direct source of afterglow radiation.

We consider the interaction between radiation, matter and a magnetic field in a compact, relativistic jet. The entrained matter accelerates outward as the jet breaks out of a star or other confining medium. In some circumstances, such as gamma-ray bursts (GRBs), the magnetization of the jet is greatly reduced by an advected radiation field while the jet is optically thick to scattering. Where magnetic flux surfaces diverge rapidly, a strong outward Lorentz force develops and radiation and matter begin to decouple. The increase in magnetization is coupled to a rapid growth in Lorentz factor. We take two approaches to this problem. The first examines the flow outside the fast magnetosonic critical surface, and calculates the flow speed and the angular distribution of the radiation field over a range of scattering depths. The second considers the flow structure on both sides of the critical surface in the optically thin regime, using a relaxation method. In both approaches, we find how the terminal Lorentz factor, and radial profile of the outflow, depend on the radiation intensity and optical depth at breakout. The effect of bulk Compton scattering on the radiation spectrum is calculated by a Monte Carlo method, while neglecting the effects of internal dissipation. The peak of the scattered spectrum sits near the seed peak if radiation pressure dominates the acceleration, but is pushed to a higher frequency if the Lorentz force dominates, and especially if the seed photon cone is broadened by interaction with a slower component of the outflow.

The theory of cold, relativistic, magnetohydrodynamic outflows is generalized by the inclusion of an intense radiation source. In some contexts, such the breakout of a gamma-ray burst jet from a star, the outflow is heated to a high temperature at a large optical depth. Eventually it becomes transparent and is pushed to a higher Lorentz factor by a combination of the Lorentz force and radiation pressure. We obtain its profile, both inside and outside the fast magnetosonic critical point, when the poloidal magnetic field is radial and monopolar. Most of the energy flux is carried by the radiation field and the toroidal magnetic field that is wound up close to the rapidly rotating engine. Although the entrained matter carries little energy, it couples the radiation field to the magnetic field. Then the fast critical point is pushed inward from infinity and, above a critical radiation intensity, the outflow is accelerated mainly by radiation pressure. We identify a distinct observational signature of this hybrid outflow: a hardening of the radiation spectrum above the peak of the seed photon distribution, driven by bulk Compton scattering. The non-thermal spectrum -- obtained by a Monte Carlo method -- is most extended when the Lorentz force dominates the acceleration, and the seed photon beam is wider than the Lorentz cone of the MHD fluid. This effect is a generic feature of hot, magnetized outflows interacting with slower relativistic material. It may explain why some GRB spectra appear to peak at photon energies above the original Amati et al. scaling. A companion paper addresses the case of jet breakout, where diverging magnetic flux surfaces yield strong MHD acceleration over a wider range of Lorentz factor.

We present the discovery and timing solutions of five new pulsars by students involved in the Pulsar Search Collaboratory (PSC), a NSF-funded joint program between the National Radio Astronomy Observatory and West Virginia University designed to excite and engage high-school students in Science, Technology, Engineering, and Mathematics (STEM) and related fields. We encourage students to pursue STEM fields by apprenticing them within a professional scientific community doing cutting edge research, specifically by teaching them to search for pulsars. The students are analyzing 300 hours of drift-scan survey data taken with the Green Bank Telescope at 350 MHz. These data cover 2876 square degrees of the sky. Over the course of five years, more than 700 students have inspected diagnostic plots through a web-based graphical interface designed for this project. The five pulsars discovered in the data have spin periods ranging from 3.1 ms to 4.8 s. Among the new discoveries are - PSR J1926-1314, a long period, nulling pulsar; PSR J1821+0155, an isolated, partially recycled 33-ms pulsar; and PSR J1400-1438, a millisecond pulsar in a 9.5-day orbit whose companion is likely a white dwarf star.

We study the gravitational collapse of a magnetized neutron star using a novel numerical approach able to capture both the dynamics of the star and the behavior of the surrounding plasma. In this approach, a fully general relativistic magnetohydrodynamics implementation models the collapse of the star and provides appropriate boundary conditions to a force-free model which describes the stellar exterior. We validate this strategy by comparing with known results for the rotating monopole and aligned rotator solutions and then apply it to study both rotating and non-rotating stellar collapse scenarios, and contrast the behavior with what is obtained when employing the electrovacuum approximation outside the star. The non-rotating electrovacuum collapse is shown to agree qualitatively with a Newtonian model of the electromagnetic field outside a collapsing star. We illustrate and discuss a fundamental difference between the force-free and electrovacuum solutions, involving the appearance of large zones of electric-dominated field in the vacuum case. This provides a clear demonstration of how dissipative singularities appear generically in the non-linear time-evolution of force-free fluids. In both the rotating and non-rotating cases, our simulations indicate that the collapse induces a strong electromagnetic transient. In the case of sub-millisecond rotation, the magnetic field experiences strong winding and the transient carries much more energy. This result has important implications for models of gamma-ray bursts.

We investigate the effects of neutrino heating and alpha-particle recombination on the hydrodynamics of core-collapse supernovae. Our focus is on the non-linear dynamics of the shock wave that forms in the collapse, and the assembly of positive energy material below it. To this end, we perform time-dependent hydrodynamic simulations with FLASH2.5 in spherical and axial symmetry. These generalize our previous calculations by allowing for bulk neutrino heating and for nuclear statistical equilibrium between n, p and alpha. The heating rate is freely tunable, as is the starting radius of the shock relative to the recombination radius of alpha-particles. An explosion in spherical symmetry involves the excitation of an overstable mode, which may be viewed as the L=0 version of the `Standing Accretion Shock Instability'. In 2D simulations, non-spherical deformations of the shock are driven by plumes of material with positive Bernoulli parameter, which are concentrated well outside the zone of strong neutrino heating. The non-spherical modes of the shock reach a large amplitude only when the heating rate is also high enough to excite convection below the shock. The critical heating rate that causes an explosion depends sensitively on the initial position of the shock relative to the recombination radius. Weaker heating is required to drive an explosion in 2D than in 1D, but the difference also depends on the size of the shock. Forcing the infalling heavy nuclei to break up into n and p below the shock only causes a slight increase in the critical heating rate, except when the shock starts out at a large radius. This shows that heating by neutrinos (or some other mechanism) must play a significant role in pushing the shock far enough out that recombination heating takes over.

We examine the stability of a standing shock wave within a spherical accretion flow onto a gravitating star, in the context of core-collapse supernova explosions. Our focus is on the effect of nuclear dissociation below the shock on the linear growth, and non-linear saturation, of non-radial oscillations of the shocked fluid. We combine two-dimensional, time-dependent hydrodynamic simulations using FLASH2.5 with a solution to the linear eigenvalue problem, and demonstrate the consistency of the two approaches. Previous studies of this `Standing Accretion Shock Instability' (SASI) have focused either on zero-energy accretion flows without nuclear dissociation, or made use of a detailed finite-temperature nuclear equation of state and included strong neutrino heating. Our main goal in this and subsequent papers is to introduce equations of state of increasing complexity, in order to isolate the various competing effects. In this work we employ an ideal gas equation of state with a constant rate of nuclear dissociation below the shock, and do not include neutrino heating. We find that a negative Bernoulli parameter below the shock significantly lowers the real frequency, growth rate, and saturation amplitude of the SASI. A decrease in the adiabatic index has similar effects. The non-linear development of the instability is characterized by an expansion of the shock driven by turbulent kinetic energy at nearly constant internal energy. Our results also provide further insight into the instability mechanism: the rate of growth of a particular mode is fastest when the radial advection time from the shock to the accretor overlaps with the period of a standing lateral sound wave. The fastest-growing mode can therefore be modified by nuclear dissociation.

We consider the details of the QED processes that create electron-positron pairs in magnetic fields approaching and exceeding 10^{14} G. The formation of free and bound pairs is addressed, and the importance of positronium dissociation by thermal X-rays is noted. We calculate the collision cross section between an X-ray and a gamma ray, and point out a resonance in the cross section when the gamma ray is close to the threshold for pair conversion. We also discuss how the pair creation rate in the open-field circuit and the outer magnetosphere can be strongly enhanced by instabilities near the light cylinder. When the current has a strong fluctuating component, a cascade develops. We examine the details of particle heating, and show that a high rate of pair creation can be sustained close to the star, but only if the spin period is shorter than several seconds. The dissipation rate in this turbulent state can easily accommodate the observed radio output of the transient radio-emitting magnetars, and even their infrared emission. Finally, we outline how a very high rate of pair creation on the open magnetic field lines can help to stabilize a static twist in the closed magnetosphere and to regulate the loss of magnetic helicity by reconnection at the light cylinder.

We consider the voltage structure in the open-field circuit and outer magnetosphere of a magnetar. The standard polar-cap model for radio pulsars is modified significantly when the polar magnetic field exceeds 1.8x10^{14} G. Pairs are created by accelerated particles via resonant scattering of thermal X-rays, followed by the nearly instantaneous conversion of the scattered photon to a pair. A surface gap is then efficiently screened by e+- creation, which regulates the voltage in the inner part of the circuit to ~10^9 V. We also examine the electrostatic gap structure that can form when the magnetic field is somewhat weaker, and deduce a voltage 10-30 times larger over a range of surface temperatures. We examine carefully how the flow of charge back to the star above the gap depends on the magnitude of the current that is extracted from the surface of the star, on the curvature of the magnetic field lines, and on resonant drag. The rates of different channels of pair creation are determined self-consistently, including the non-resonant scattering of X-rays, and collisions between gamma rays and X-rays. We find that the electrostatic gap solution has too small a voltage to sustain the observed pulsed radio output of magnetars unless i) the magnetic axis is nearly aligned with the rotation axis and the light of sight; or ii) the gap is present on the closed as well as the open magnetic field lines. Several properties of the radio magnetars -- their rapid variability, broad pulses, and unusually hard radio spectra -- are consistent with a third possibility, that the current in the outer magnetosphere is strongly variable, and a very high rate of pair creation is sustained by a turbulent cascade.

We present an interpretation of the phenomenological relations between the spectral peak, isotropic luminosity and duration of long gamma ray bursts that have been discovered by Amati et al., Ghirlanda et al., Firmani et al., and Liang & Zhang. In our proposed model, a jet undergoes internal dissipation which prevents its bulk Lorentz factor from exceeding 1/theta (theta being the jet opening angle) until it escapes from the core of its progenitor star, at a radius of order 1e10 cm; dissipation may continue at larger radii. The dissipated radiation will be partially thermalized, and we identify its thermal peak (Doppler boosted by the outflow) with E_pk. The radiation comes, in effect, from within the jet photosphere. The non-thermal, high energy part of the GRB emission arises from Comptonization of this radiation by relativistic electrons and positrons outside the effective photosphere. This model can account naturally not only for the surprisingly small scatter in the various claimed correlations, but also for the normalization, as well as the slopes. It then has further implications for the jet energy, the limiting jet Lorentz factor, and the relation of the energy, opening angle and burst duration to the mass and radius of the stellar stellar progenitor. The observed relation between pulse width and photon frequency can be reproduced by inverse-Compton emission at ~ 1e14 cm from the engine, but there are significant constraints on the energy distribution and isotropy of the radiating particles.

Although the surface of a magnetar is a source of bright thermal X-rays, its spectrum contains substantial non-thermal components. The X-ray emission is pulsed, with pulsed fractions that can be as high as ~ 70%. Several properties of magnetars indicate the presence of persistent, static currents flowing across the stellar surface and closing within the magnetosphere. The charges supporting these currents supply a significant optical depth to resonant cyclotron scattering in the 1-100 keV band. Here we describe a Monte Carlo approach to calculating the redistribution of thermal seed photons in frequency and angle by multiple resonant scattering in the magnetosphere. The calculation includes the full angular dependence of the cyclotron scattering cross section, the relativistic Doppler effect due to the motion of the charges, and allows for an arbitrary particle velocity distribution and magnetic field geometry. We construct synthetic spectra and pulse profiles for arbitrary orientations of the spin axis, magnetic axis, and line of sight, using a self-similar, twisted dipole field geometry, and assuming that the seed photons are supplied by single-temperature black body emission from the stellar surface. Pulse profiles and 1-10 keV spectra typical of AXPs are easily produced by this model, with pulsed fractions of ~ 50%. However, this model cannot reproduce the hard, rising energy spectra that are observed from SGRs during periods of activity, without overproducing the thermal emission peak. This suggests that the 1-100 keV emission of SGRs has a common origin with the hard X-ray emission detected from some AXPs above ~20 keV.

Neutrino transport in spherically symmetric models of stellar core collapse and bounce has achieved a technically complete level, rewarded by the agreement among independent groups that a multi-dimensional treatment of the fluid-instabilities in the post-bounce phase is indispensable to model supernova explosions. While much effort is required to develop a reliable neutrino transport technique in axisymmetry, we explore neutrino physics approximations and parameterizations for an efficient three-dimensional simulation of the fluid-instabilities in the shock-heated matter that accumulates between the accretion shock and the protoneutron star. We demonstrate the reliability of a simple parameterization scheme in the collapse phase and extend our 3D magneto-hydrodynamical collapse simulations to a preliminary postbounce evolution. The growth of magnetic fields is investigated.

Persistent high-energy emission of magnetars is produced by a plasma corona around the neutron star, with total energy output of ~10^{36}erg/s. The corona forms as a result of occasional starquakes that twist the external magnetic field of the star and induce electric currents in the closed magnetosphere. Once twisted, the magnetosphere cannot untwist immediately because of its self-induction. The self-induction electric field lifts particles from the stellar surface, accelerates them, and initiates avalanches of pair creation in the magnetosphere. The created plasma corona maintains the electric current demanded by curl(B) and regulates the self-induction e.m.f. by screening. This corona persists in dynamic equilibrium: it is continually lost to the stellar surface on the light-crossing time and replenished with new particles. In essence, the twisted magnetosphere acts as an accelerator that converts the toroidal field energy to particle kinetic energy. The voltage along the magnetic field lines is maintained near threshold for ignition of pair production, in the regime of self-organized criticality. The voltage is found to be about 1 GeV which is in agreement with the observed dissipation rate \~10^{36}erg/s. The coronal particles impact the solid crust, knock out protons, and regulate the column density of the hydrostatic atmosphere of the star. The transition layer between the atmosphere and the corona is the likely source of the observed 100-keV emission from magnetars. The corona also emits curvature radiation up to 10^{14}Hz and can supply the observed IR-optical luminosity.

We develop a theoretical model that explains the formation of hot coronae around strongly magnetized neutron stars -- magnetars. The starquakes of a magnetar shear its external magnetic field, which becomes non-potential and is threaded by an electric current. Once twisted, the magnetosphere cannot untwist immediately because of its self-induction. The self-induction electric field lifts particles from the stellar surface, accelerates them, and initiates avalanches of pair creation in the magnetosphere. The created plasma corona maintains the electric current demanded by curl(B) and regulates the self-induction e.m.f. by screening. This corona persists in dynamic equilibrium: it is continually lost to the stellar surface on the light-crossing time of 10^{-4} s and replenished with new particles. In essence, the twisted magnetosphere acts as an accelerator that converts the toroidal field energy to particle kinetic energy. Using a direct numerical experiment, we show that the corona self-organizes quickly (on a millisecond timescale) into a quasi-steady state, with voltage ~1 GeV along the magnetic lines. The heating rate of the corona is ~10^{36} erg/s, in agreement with the observed persistent, high-energy output of magnetars. We deduce that a static twist that is suddenly implanted into the magnetosphere will decay on a timescale of 1-10 yrs. The particles accelerated in the corona impact the solid crust, knock out protons, and regulate the column density of the hydrostatic atmosphere of the star. The transition layer between the atmosphere and the corona is the likely source of the observed 100-keV emission from magnetars. The corona emits curvature radiation and can supply the observed IR-optical luminosity. (Abridged)

(Abridged) We consider the interaction of a relativistically-moving shell, composed of thermal photons, a reversing magnetic field and a small admixture of charged particles, with a dense Wolf-Rayet wind. A thin outer layer of Wolf-Rayet material is entrained by the jet head; it cools and becomes Rayleigh-Taylor unstable, thereby providing an additional source of inertia and variability. Pair creation in the wind material, and the associated pre-acceleration, defines a characteristic radiative compactness at the point where the reverse shock has completed its passage back through the shell. We argue that the prompt gamma-ray emission is triggered by this external braking, at an optical depth ~1 to electron scattering. Torsional waves, excited by the forced reconnection of the reversing magnetic field, carry a fluctuating current, and are damped at high frequencies by the electrostatic acceleration of electrons and positrons. We show that inverse Compton radiation by the accelerated charges is stronger than their synchrotron emission, and is beamed along the magnetic field. Thermal radiation that is advected out from the base of the jet cools the particles. The observed relation between peak energy and isotropic luminosity is reproduced if the blackbody seeds are generated in a relativistic jet core that is subject to Kelvin-Helmholtz instabilities with the Wolf-Rayet envelope. This relation is predicted to soften to E_peak ~ L_iso^{1/4} below an isotropic luminosity L_iso ~ 3x10^{50} ergs/s. The duration of spikes in the inverse-Compton emission is narrower at higher frequencies, in agreement with the observed relation. The transition from prompt gamma-ray emission to afterglow can be explained by the termination of the thermal X-ray seed and the onset of synchrotron-self-Compton emission.

The recently discovered soft gamma-ray emission from the anomalous X-ray pulsar 1E 1841-045 has a luminosity L_g ~ 10^{36} ergs/s. This luminosity exceeds the spindown power by three orders of magnitude and must be fed by an alternative source of energy such as an ultrastrong magnetic field. A gradual release of energy in the stellar magnetosphere is expected if it is twisted and a strong electric current is induced on the closed field lines. We examine two mechanisms of gamma-ray emission associated with the gradual dissipation of this current. (1) A thin surface layer of the star is heated by the downward beam of current-carrying charges, which excite Langmuir turbulence in the layer. As a result, it can reach a temperature kT ~ 100 keV and emit bremsstrahlung photons up to this characteristic energy. (2) The magnetosphere is also a source of soft gamma rays at a distance of ~100 km from the star, where the electron cyclotron energy is in the keV range. A large electric field develops in this region in response to the outward drag force felt by the current-carrying electrons from the flux of keV photons leaving the star. A seed positron injected in this region undergoes a runaway acceleration and upscatters keV photons above the threshold for pair creation. The created pairs emit a synchrotron spectrum consistent with the observed 20-100 keV emission. This spectrum is predicted to extend to higher energies and reach a peak at ~1 MeV.

The propagation and non-linear interactions of magnetohydrodynamic waves are considered in the force-free limit, where the inertia of the conducting matter which enforces the MHD condition E.B = 0 can be neglected in comparison with the inertia of the electromagnetic field. By extending the analysis beyond the WKB approximation, we are able to study the non-linearities induced by a gravitational field. We treat the perturbed electromagnetic field as a fluid of infinite conductivity. We calculate the scattering of a torsional (Alfven) wave by a gravitational potential, and demonstrate a nonlinear coupling with a compressive (fast) wave which is second order in the amplitude of the Alfven wave. In a cylindrically symmetric spacetime with slow rotation, the coupling is second order in g_{t\phi} and first order in the amplitude of the wave. We also give a fresh analysis of the non-linear interactions between compressive and torsional waves in Minkowski space, with a focus on the relative strengths of their three- and four-mode interactions. In contrast with non-relativistic magnetofluids, the effects of compression are always present. In the case of colliding fast waves, a net displacement of the field lines across (at least) one of the colliding wavepackets is shown to have a strong effect on the outgoing waveform, and to have a qualitatively different interpretation than was previously suggested for colliding Alfven waves. Finally, we show how spacetime curvature modifies the collision between two torsional waves, in both the weak- and strong-field regimes.

Soft-gamma-ray repeaters (SGRs) are galactic X-ray stars that emit numerous short-duration (about 0.1 s) bursts of hard X-rays during sporadic active periods. They are thought to be magnetars: strongly magnetized neutron stars with emissions powered by the dissipation of magnetic energy. Here we report the detection of a long (380 s) giant flare from SGR 1806-20, which was much more luminous than any previous transient event observed in our Galaxy. (In the first 0.2 s, the flare released as much energy as the Sun radiates in a quarter of a million years.) Its power can be explained by a catastrophic instability involving global crust failure and magnetic reconnection on a magnetar, with possible large-scale untwisting of magnetic field lines outside the star. From a great distance this event would appear to be a short-duration, hard-spectrum cosmic gamma-ray burst. At least a significant fraction of the mysterious short-duration gamma-ray bursts therefore may come from extragalactic magnetars.

(Abridged) In the binary radio pulsar system J0737-3039, the faster pulsar A is eclipsed once per orbit. We construct a simple geometric model which successfully reproduces the eclipse light curves, based on the idea that the radio pulses are attenuated by synchrotron absorption on the closed magnetic field lines of pulsar B. The model explains most of the properties of the eclipse: its asymmetric form, the nearly frequency-independent duration, and the modulation of the brightness of pulsar A at both once and twice the rotation frequency of pulsar B in different parts of the eclipse. This detailed agreement confirms the dipolar structure of the star's poloidal magnetic field. The model makes clear predictions for the degree of linear polarization of the transmitted radiation. The weak frequency dependence of the eclipse duration implies that the absorbing plasma is relativistic, with a density much larger than the corotation charge density. Such hot, dense plasma can be effectively stored in the outer magnetosphere, where cyclotron cooling is slow. The gradual loss of particles inward through the cooling radius is compensated by an upward flux driven by a fluctuating component of the current, and by the pumping of magnetic helicity on the closed field lines. The trapped particles are heated to relativistic energies by the damping of magnetospheric turbulence and, at a slower rate, by the absorption of the radio emission of the companion pulsar.

This article is a review of Soft Gamma Repeaters and Anomalous X-ray Pulsars. It contains a brief historical record of the emergence of these classes of neutron stars, a thorough overview of the observational data, a succinct summary of the magnetar model, and suggested directions for future research in this field.

All recent numerical simulations agree that stars in the main sequence mass range of 9-40 solar masses do not produce a prompt hydrodynamic ejection of the outer layers after core collapse and bounce. Rather they suggest that stellar core collapse and supernova explosion are dynamically distinct astrophysical events, separated by an unspectacular accretion phase of at least ~40 ms duration. As long as the neutrinospheres remain convectively stable, the explosion dynamics is determined by the neutrons, protons, electrons and neutrinos in the layer of impact-heated matter piling up on the protoneutron star. The crucial role of neutrino transport in this regime has been emphasized in many previous investigations. Here, we search for efficient means to address the role of magnetic fields and fluid instabilities in stellar core collapse and the postbounce phase.

We present the results obtained by INTEGRAL on the Soft-Gamma Ray Repeater SGR 1806-20. In particular we report on the temporal and spectral properties of the bursts detected during a moderately active period of the source in September and October 2003 and on the search for quiescent emission.

(abridged) An outburst of more than 80 individual bursts, similar to those seen from Soft Gamma Repeaters (SGRs), was detected from the Anomalous X-ray Pulsar (AXP) 1E 2259+586 in 2002 June. Coincident with this burst activity were gross changes in the pulsed flux, persistent flux, energy spectrum, pulse profile and spin down of the underlying X-ray source. We present RXTE and XMM-Newton observations of 1E 2259+586 that show the evolution of the aforementioned source parameters during and following this episode. Specifically, we observe an X-ray flux increase by more than an order of magnitude having two distinct components. The first component is linked to the burst activity and decays within ~2 days during which the energy spectrum is considerably harder than during the quiescent state of the source. The second component decays over the year following the glitch according to a power law in time with an exponent -0.22 +/- 0.01. The pulsed fraction decreased initially to ~15% RMS, and the pulse profile changed significantly during the outburst. A glitch was observed in 1E 2259+586 that preceded the observed burst activity. A fraction of the glitch (~19%) recovered, although the recovery was not purely exponential. An exponential rise of ~20% of the frequency jump with a time scale of ~14 days results in a significantly better fit to the data, however, contamination from a systematic drift in the phase of the pulse profile cannot be excluded. The long-term post-glitch spin-down rate decreased in magnitude relative to the pre-glitch value. A comparison with SGR outburst properties, a physical interpretation of the results, and implications on the number of magnetar candidates in our Galaxy are discussed.

We calculate the coupled thermal evolution and magnetic field decay in relativistic model neutron stars threaded by superstrong magnetic fields (B > 10^{15} G). Our main goal is to evaluate how such ``magnetars'' evolve with time and how field decay modifies the transitions to core superfluidity and cooling dominated by surface X-ray emission. Observations of a thermal X-ray spectral component and fast timing noise place strong constraints on the presence of a superfluid core. We find that the transition to core superfluidity can be significantly delayed by field decay in the age range ~ 10^3-10^5 yrs. The mechanism of Hall drift is related to the stability of the core magnetic field, and to currents flowing outward through the crust. The heating effect is enhanced if it is continuous rather than spasmodic. Condensation of a heavy element layer at the surface is shown to cause only modest changes in the outward conduction of heat.

We report on observations of SGR 1900+14 made with the Rossi X-ray Timing Explorer (RXTE) and BeppoSAX during the April 2001 burst activation of the source. Using these data, we measure the spindown torque on the star and confirm earlier findings that the torque and burst activity are not directly correlated. We compare the X-ray pulse profile to the gamma-ray profile during the April 18 intermediate flare and show that (i) their shapes are similar and (ii) the gamma-ray profile aligns closely in phase with the X-ray pulsations. The good phase alignment of the gamma-ray and X-ray profiles suggests that there was no rapid spindown following this flare of the magnitude inferred for the August 27 giant flare. We discuss how these observations further constrain magnetic field reconfiguration models for the large flares of SGRs.

The X-ray afterglow that is observed following large flares on magnetars can be accurately fit by simple and quantitative theoretical models: The long term afterglow, lasting of order weeks, can be understood as thermal radiation of a heated neutron star crust that is re-scattered in the magnetosphere. Short term afterglow is well fit by the cooling of a non-degenerate, pair-rich layer, which gradually shrinks and releases heat to a pair-free zone above it. Measurements of persistent optical and infrared emission directly probe long lived currents in the magnetosphere which are a likely source of collective plasma emission. The superstrong magnetic field plays an important role in generating these various emissions, and previous inference of its strength in magnetars is supported by the good fits with observation.

Transient X-ray emission, with an approximate t^{-0.7} decay, was observed from SGR 1900+14 over 40 days following the the giant flare of 27 Aug 1998. We calculate in detail the diffusion of heat to the surface of a neutron star through an intense 10^{14}-10^{15} G magnetic field, following the release of magnetic energy in its outer layers. We show that the power law index, the fraction of burst energy in the afterglow, and the return to persistent emission can all be understood if the star is composed of normal baryonic material.

(ABBREVIATED) We consider the structure of neutron star magnetospheres threaded by large-scale electrical currents, and the effect of resonant Compton scattering by the charge carriers (both electrons and ions) on the emergent X-ray spectra and pulse profiles. In the magnetar model for the SGRs and AXPs, these currents are maintained by magnetic stresses acting deep inside the star. We construct self-similar, force-free equilibria of the current-carrying magnetosphere with a power-law dependence of magnetic field on radius, B ~ r^(-2-p), and show that a large-scale twist softens the radial dependence to p < 1. The spindown torque acting on the star is thereby increased in comparison with a vacuum dipole. We comment on the strength of the surface magnetic field in the SGR and AXP sources, and the implications of this model for the narrow measured distribution of spin periods. A magnetosphere with a strong twist, B_\phi/B_\theta = O(1) at the equator, has an optical depth ~ 1 to resonant cyclotron scattering, independent of frequency (radius), surface magnetic field strength, or charge/mass ratio of the scattering charge. When electrons and ions supply the current, the stellar surface is also heated by the impacting charges at a rate comparable to the observed X-ray output of the SGR and AXP sources, if B_{dipole} ~ 10^{14} G. Redistribution of the emerging X-ray flux at the ion and electron cyclotron resonances will significantly modify the emerging pulse profile and, through the Doppler effect, generate a non-thermal tail to the X-ray spectrum. The sudden change in the pulse profile of SGR 1900+14 after the 27 August 1998 giant flare is related to an enhanced optical depth to electron cyclotron scattering, resulting from a sudden twist imparted to the external magnetic field.

After a couple of years of quiescence, the soft gamma repeater SGR 1900+14 suddenly reactivated on 18 April 2001, with the emission of a very intense, long and modulated flare, only second in intensity and duration to the 27 August 1998 giant flare. BeppoSAX caught the large flare with its Gamma Ray Burst Monitor and with one of the Wide Field Cameras. The Wide Field Cameras also detected X-ray bursting activity shortly before the giant flare. A target of opportunity observation was started only 8 hours after the large flare with the Narrow Field Instruments, composed of two 60-ks long pointings. These two observations show an X-ray afterglow of the persistent SGR 1900+14 source, decaying with time according to a power law of index -0.6.

In the framework of the magnetar model for the Soft Gamma Repeaters and Anomalous X-ray Pulsars, we consider the structure of neutron star magnetospheres threaded by large-scale electrical currents. We construct self-similar, force-free equilibria under the assumption of axisymmetry and a power law dependence of magnetic field on radius, ${\bf B} \propto r^{-(2+p)}$. A large-scale twist of the field lines softens the radial dependence to $p < 1$, thereby accelerating the spindown torque with respect to a vacuum dipole. A magnetosphere with a strong twist ($B_\phi/B_\theta = O(1)$ at the equator) has an optical depth $\sim 1$ to resonant cyclotron scattering, independent of frequency (radius), surface magnetic field strength, or the charge/mass ratio of the scattering charge. We investigate the effects of the resonant Compton scattering by the charge carriers (both electrons and ions) on the emergent X-ray spectra and pulse profiles.

(ABBREVIATED) The extraordinary 1998 August 27 giant flare places strong constraints on the physical properties of its source, SGR 1900+14. We make detailed comparisons of the published data with the magnetar model. The giant flare evolved through three stages, whose radiative mechanisms we address in turn. A triggering mechanism is proposed, whereby a helical distortion of the core magnetic field induces large-scale fracturing in the crust and a twisting deformation of the crust and exterior magnetic field. The envelope of the pulsating tail of the August 27 flare can be accurately fit, after ~40 s, by the contracting surface of a relativistically hot, but inhomogeneous, trapped fireball. We quantify the effects of direct neutrino-pair emission, thereby deducing a lower bound ~ 10^{32} G-cm^3 to the magnetic moment of the confining field. The radiative flux during the intermediate ~40 s of the burst appears to exceed the trapped fireball fit. The spectrum and lightcurve of this smooth tail are consistent with heating in an extended pair corona, possibly powered by continuing seismic activity in the star. We consider in detail the critical luminosity, below which a stable balance can be maintained between heating and radiative cooling in a confined, magnetized pair plasma; but above which the confined plasma runs away to local thermodynamic equilibrium. In the later pulsating tail, the best fit temperature equilibrates at a value which agrees well with the regulating effect of photon splitting. The remarkable four-peaked substructure within each 5.16-s pulse provides strong evidence for the presence of higher magnetic multipoles in SGR 1900+14. The corresponding collimation of the X-ray flux is related to radiative transport in a super-QED magnetic field.

I summarize recent developments in the magnetar model of the Soft Gamma Repeaters and Anomalous X-ray Pulsars, give a critical inventory of alternative models for the AXPs, and outline the improved diagnostics expected from present observational efforts.

We consider the amplification and transport of a magnetic field in the collapsed core of a massive star, including both the region between the neutrinosphere and the shock, and the central, opaque core. An analytical argument explains why rapid convective overturns persist within a newly formed neutron star for roughly 10 seconds ($> 10^3$ overturns), consistent with recent numerical models. A dynamical balance between turbulent and magnetic stresses within this convective layer corresponds to flux densities in excess of $10^{15}$G. Material accreting onto the core is heated by neutrinos and also becomes strongly convective. We compare the expected magnetic stresses in this convective `gain layer' with those deep inside the neutron core. Buoyant motions of magnetized fluid are greatly aided by the intense neutrino flux. We calculate the transport rate through a medium containing free neutrons protons, and electrons, in the limiting cases of degenerate or non-degenerate nucleons. Fields stronger than $\sim 10^{13}$ G are able to rise through the outer degenerate layers of the neutron core during the last stages of Kelvin-Helmholtz cooling (up to 10 seconds post-collapse), even though these layers have become stable to convection. We also find the equilibrium shape of a thin magnetic flux rope in the dense hydrostatic atmosphere of the neutron star, along with the critical separation of the footpoints above which the rope undergoes unlimited expansion against gravity. The implications of these results for pulsar magnetism are summarized, and applied to the case of late fallback over the first 1,000-10,000 s of the life of a neutron star

We study the temporal and coarse spectral properties of 268 bursts from SGR 1806-20 and 679 bursts from SGR 1900+14, all observed with the Rossi X-Ray Timing Explorer/Proportional Counter Array. Hardness ratios and temporal parameters, such as T_90 durations and tau_90 emission times are determined for these bursts. We find a lognormal distribution of burst durations, ranging over more than two orders of magnitude: T_90 ~ 10^-2 to >~ 1 s, with a peak at ~ 0.1 s. The burst light curves tend to be asymmetrical, with more than half of all events showing rise times t_r < 0.3 T_90. We find that there exists a correlation between the duration and fluence of bursts from both sources. We also find a significant anti-correlation between hardness ratio and fluence for SGR 1806-20 bursts and a marginal anti-correlation for SGR 1900+14 events. Finally, we discuss possible physical implications of these results within the framework of the magnetar model.

We report the detection of large flux changes in the persistent X-ray flux of SGR 1900+14 during its burst active episode in 1998. Most notably, we find a factor ~700 increase in the non-burst X-ray flux following the August 27th flare, which decayed in time as a power-law. Our measurements indicate that the pulse fraction remains constant throughout this decay. This suggests a global flux enhancement as a consequence of the August 27th flare rather than localized heating. While the persistent flux has since recovered to the pre-outburst level, the pulse profile has not. The pulse shape changed to a near sinusoidal profile within the tail of the August 27th flare (in gamma-rays) and this effect has persisted for more than 1.5 years (in X-rays). The results presented here suggest the magnetic field of the neutron star in SGR 1900+14 was significantly altered (perhaps globally) during the giant flare of August 27.

The giant flare of 1998 August 27 from SGR 1900+14 was extraordinary in many ways: it was the most intense flux of gamma rays ever detected from a source outside our solar system; it was longer than any previously detected burst from a soft gamma repeater (SGR) in our Galaxy by more than an order of magnitude; and it showed a remarkable four-peaked, periodic pattern in hard X-rays with the same rotation period that was found modulating soft X-rays from the star in quiescence. The event was detected by several gamma-ray experiments in space, including the Ulysses gamma-ray burst detector and the BeppoSAX Gamma Ray Burst Monitor. These instruments operate in different energy ranges, and comparisons of their measurements reveal complex patterns of spectral evolution as the intensity varies. In this paper, we present a joint analysis of the BeppoSAX and Ulysses data and discuss some implications of these results for the SGRs. We also present newly-analyzed Venera/SIGNE and ISEE-3 data on the 1979 March 5 giant flare from an SGR in the Large Magellanic Cloud (SGR 0526-66), and compare them with the August 27 event. Our results are consistent with the hypothesis that giant flares are due to catastrophic magnetic instabilities in highly magnetized neutron stars, or "magnetars".

I summarize the recent advances in our understanding of the Soft Gamma Repeaters: in particular their spin behavior, persistent emission and hyper-Eddington outbursts. The giant flares on 5 March 1979 and 27 August 1998 provide compelling physical evidence for magnetic fields stronger than 10 B_{QED} = 4.4 x 10^{14} G, consistent with the rapid spindown detected in two of these sources. The persistent X-ray emission and variable spindown of the 6-12 s Anomalous X-ray Pulsars are compared and contrasted with those of the SGRs, and the case made for a close connection between the two types of sources. Their collective properties point to the existence of {\it magnetars}: neutron stars in which a decaying magnetic field (rather than accretion or rotation) is the dominant source of energy for radiative and particle emissions. Observational tests of the magnetar model are outlined, along with current ideas about the trigger of SGR outbursts, new evidence for the trapped fireball model, and the influence of QED processes on X-ray spectra and lightcurves. A critical examination is made of coherent radio emission from bursting strong-field neutron stars. I conclude with an overview of the genetic connection between neutron star magnetism and the violent fluid motions in a collapsing supernova core.

The Soft Gamma Repeater SGR 1900+14 entered a remarkable phase of activity during the summer of 1998. This activity peaked on August 27, 1998 when a giant periodic gamma-ray flare resembling the famous March 5, 1979 event from SGR 0526-66 was recorded. Two days later (August 29), a strong, bright burst was detected with RXTE and BATSE. This event reveals several similarities to the giant flares of August 27 and March 5, and shows a number of unique features not previously seen in SGR bursts. Unlike typically short SGR bursts, this event features a 3.5 s burst peak that was preceded by an extended (~ 1 s) complex precursor, and followed by a long (~ 1000 s) periodic tail modulated at the 5.16 s stellar rotation period. Spectral analysis shows a striking distinction between the spectral behavior of the precursor, burst peak and extended tail. While the spectrum during the peak is uniform, a significant spectral evolution is detected in both the precursor and tail emissions. Temporal behavior shows a sharp rise (~ 9.8 ms) at the event onset and a rapid cutoff (~ 17 ms) at the end of the burst peak. The tail pulsations show a simple pulse profile consisting of one 5.16 s peak that did not evolve with time. The spectral and temporal signatures of this event imply that the precursor, main peak, and extended tail are produced by different physical mechanisms. We discuss these features and their implications in the context of the magnetar model. The energetics of the August 29 event, and its close proximity to the August 27 flare, suggest that it is an `aftershock' of the preceding giant flare. "P.S. This is an abbreviated version of the original abstract."

We have phase connected a sequence of RXTE PCA observations of SGR 1806-20 covering 178 days. We find a simple secular spin-down model does not adequately fit the data. The period derivative varies gradually during the observations between 8.1 and 11.7 * 10^-11 s/s (at its highest, ~40% larger than the long term trend), while the average burst rate as seen with BATSE drops throughout the time interval. The phase residuals give no compelling evidence for periodicity, but more closely resemble timing noise as seen in radio pulsars. The magnitude of the timing noise, however, is large relative to the noise level typically found in radio pulsars. Combining these results with the noise levels measured for some AXPs, we find all magnetar candidates have \Delta_8 values larger than those expected from a simple extrapolation of the correlation found in radio pulsars. We find that the timing noise in SGR 1806-20 is greater than or equal to the levels found in some accreting systems (e.g., Vela X-1, 4U 1538-52 and 4U 1626-67), but the spin-down of SGR 1806-20 has thus far maintained coherence over 6 years. Alternatively, an orbital model with a period P_orb = 733 days provides a statistically acceptable fit to the data. If the phase residuals are created by Doppler shifts from a gravitationally bound companion, then the allowed parameter space for the mass function (small) and orbital separation (large) rule out the possibility of accretion from the companion sufficient to power the persistent emission from the SGR.

We consider the effects of rapid pair creation by an intense pulse of gamma-rays propagating ahead of a relativistic shock. Side-scattered photons colliding with the main gamma-ray beam amplify the density of scattering charges. The acceleration rate of the pair-loaded medium is calculated, and its limiting bulk Lorentz factor related to the spectrum and compactness of the photon source. One obtains, as a result, a definite prediction for the relative inertia in baryons and pairs. The deceleration of a relativistic shock in the moving medium, and the resulting synchrotron emissivity, are compared with existing calculations for a static medium. The radiative efficiency is increased dramatically by pair loading. When the initial ambient density exceeds a critical value, the scattering depth traversed by the main gamma-ray pulse rises above unity, and the pulse is broadened. These considerations place significant constraints on burst progenitors: a pre-burst mass loss rate exceeding 10^{-5} M_\odot per year is difficult to reconcile with individual pulses narrower than 10 s, unless the radiative efficiency is low. An anisotropic gamma-ray flux (on an angular scale \Gamma^{-1} or larger) drives a large velocity shear that greatly increases the energy in the seed magnetic field forward of the propagating shock.

We present statistics of SGR 1806-20 bursts, combining 290 events detected with RXTE/PCA, 111 events detected with BATSE and 134 events detected with ICE. We find that the fluence distribution of bursts observed with each instrument are well described by power laws with indices 1.43, 1.76 and 1.67, respectively. The distribution of time intervals between successive bursts from SGR 1806-20 is described by a lognormal function with a peak at 103 s. There is no correlation between the burst intensity and either the waiting times till the next burst or the time elapsed since the previous burst. In all these statistical properties, SGR 1806-20 bursts resemble a self-organized critical system, similar to earthquakes and solar flares. Our results thus support the hypothesis that the energy source for SGR bursts is crustquakes due to the evolving, strong magnetic field of the neutron star, rather than any accretion or nuclear power.

The extraordinary giant flare of 1998 August 27 from SGR 1900+14 was the most intense event ever detected from this or any other cosmic source (even more intense than the famous March 5th 1979 event). It was longer than any previous burst from SGR1900+14 by more than one order of magnitude, and it displayed the same 5.16-s periodicity in hard X-rays that was detected in the low energy X-ray flux of its quiescent counterpart. The event was detected by several gamma-ray experiments in space, among them the Ulysses gamma-ray burst detector and the BeppoSAX Gamma Ray Burst Monitor. These instruments operate in different energy ranges, and a comparison of their data shows that the event emitted a strongly energy-dependent flux, and displayed strong spectral evolution during the outburst itself. Here we present a joint analysis of the BeppoSAX and Ulysses data, in order to identify the energy-dependent features of this event and understand some of the physical conditions in the environment of the neutron star which generated this flare.

We consider the radiative acceleration to relativistic bulk velocities of a cold, optically thin plasma which is exposed to an external source of gamma-rays. The flow is driven by radiative momentum input to the gas, the accelerating force being due to Compton scattering in the relativistic Klein-Nishina limit. The bulk Lorentz factor of the plasma, Gamma, derived as a function of distance from the radiating source, is compared with the corresponding result in the Thomson limit. Depending on the geometry and spectrum of the radiation field, we find that particles are accelerated to the asymptotic Lorentz factor at infinity much more rapidly in the relativistic regime; and the radiation drag is reduced as blueshifted, aberrated photons experience a decreased relativistic cross section and scatter preferentially in the forward direction. The random energy imparted to the plasma by gamma-rays can be converted into bulk motion if the hot particles execute many Larmor orbits before cooling. This `Compton afterburn' may be a supplementary source of momentum if energetic leptons are injected by pair creation, but can be neglected in the case of pure Klein-Nishina scattering. Compton drag by side-scattered radiation is shown to be more important in limiting the bulk Lorentz factor than the finite inertia of the accelerating medium. The processes discussed here may be relevant to a variety of astrophysical situations where luminous compact sources of hard X- and gamma-ray photons are observed, including active galactic nuclei, galactic black hole candidates, and gamma-ray bursts.

We study the statistics of soft gamma repeater (SGR) bursts, using a data base of 187 events detected with BATSE and 837 events detected with RXTE PCA, all from SGR 1900+14 during its 1998-1999 active phase. We find that the fluence or energy distribution of bursts is consistent with a power law of index 1.66, over 4 orders of magnitude. This scale-free distribution resembles the Gutenberg-Richter Law for earthquakes, and gives evidence for self-organized criticality in SGRs. The distribution of time intervals between successive bursts from SGR 1900+14 is consistent with a log-normal distribution. There is no correlation between burst intensity and the waiting times till the next burst, but there is some evidence for a correlation between burst intensity and the time elapsed since the previous burst. We also find a correlation between the duration and the energy of the bursts, but with significant scatter. In all these statistical properties, SGR bursts resemble earthquakes and solar flares more closely than they resemble any known accretion-powered or nuclear-powered phenomena. Thus our analysis lends support to the hypothesis that the energy source for SGR bursts is internal to the neutron star, and plausibly magnetic.

We present evidence for burst emission from SGR 1900+14 with a power-law high energy spectrum extending beyond 500 keV. Unlike previous detections of high energy photons during bursts from SGRs, these emissions are not associated with high-luminosity burst intervals. Not only is the emission hard, but the spectra are better fit by Band's GRB function rather than by the traditional optically-thin thermal bremsstrahlung model. We find that the spectral evolution within these hard events obeys a hardness/intensity anti-correlation. Temporally, these events are distinct from typical SGR burst emissions in that they are longer (~ 1 s) and have relatively smooth profiles. Despite a difference in peak luminosity of > 1E+11 between these bursts from SGR 1900+14 and cosmological GRBs, there are striking temporal and spectral similarities between the two kinds of bursts, aside from spectral evolution. We outline an interpretation of these events in the context of the magnetar model.

We consider the physical implications of the rapid spindown of Soft Gamma Repeater 1900+14, and of the apparent "braking glitch", \Delta P/P = l x 10^-4, that was concurrent with the Aug. 27th giant flare. A radiation-hydrodynamical outflow associated with the flare could impart the required torque, but only if the dipole magnetic field is stronger than ~ 10^14 G and the outflow lasts longer and/or is more energetic than the observed X-ray flare. A positive period increment is also a natural consequence of a gradual, plastic deformation of the neutron star crust by an intense magnetic field, which forces the neutron superfluid to rotate more slowly than the crust. Sudden unpinning of the neutron vortex lines during the August 27th event would then induce a glitch opposite in sign to those observed in young pulsars, but of a much larger magnitude as a result of the slower rotation. The change in the persistent X-ray lightcurve following the August 27 event is ascribed to continued particle heating in the active region of that outburst. The enhanced X-ray output can be powered by a steady current flowing through the magnetosphere, induced by the twisting motion of the crust. The long term rate of spindown appears to be accelerated with respect to a simple magnetic dipole torque. Accelerated spindown of a seismically-active magnetar will occur when its persistent output of Alfven waves and particles exceeds its spindown luminosity. We suggest that SGRs experience some episodes of relative inactivity, with diminished spindown rates, and that such inactive magnetars are observed as Anomalous X-ray Pulsars (AXPs). The rapid reappearence of persistent X-ray emission following August 27 flare gives evidence against accretion-powered models.

We have analyzed Rossi X-ray Timing Explorer Proportional Counter Array observations of the pulsed emission from SGR 1900+14 during September 1996, June - October 1998, and early 1999. Using these measurements and results reported elsewhere, we construct a period history of this source for 2.5 years. We find significant deviations from a steady spin-down trend during quiescence and the burst active interval. Burst and Transient Source Experiment observations of the burst emission are presented and correlations between the burst activity and spin-down rate of SGR 1900+14 are discussed. We find an 80 day interval during the summer of 1998 when the average spin-down rate is larger than the rate elsewhere by a factor ~ 2.3. This enhanced spin-down may be the result of a discontinuous spin-down event or ``braking glitch'' at the time of the giant flare on 27 August 1998. Furthermore, we find a large discrepancy between the pulsar period and average spin-down rate in X-rays as compared to radio observations for December 1998 and January 1999.

The soft-gamma repeater SGR 1900+14 became active again on June 1998 after a long period of quiescence; it remained at a low state of activity until August 1998, when it emitted a series of extraordinarily intense outbursts. We have observed the source with RXTE twice, during the onset of each active episode. We confirm the pulsations at the 5.16 s period reported earlier (Hurley et al. 1998b, Hurley et al. 1998 e) from SGR 1900+14. Here we report the detection of a secular spindown of the pulse period at an average rate of 1.1*10^{-10} s/s. In view of the strong similarities between SGRs, we attribute the spindown of SGR 1900+14 to magnetic dipole radiation, possibly accelerated by a quiescent flux, as in the case of SGR 1806-20 (Kouveliotou et al. 1998a). This allows an estimate of the pulsar dipolar magnetic field, which is 2-8*10^{14} G. Our results confirm that SGRs are magnetars.