We investigate the saturation of particle acceleration in relativistic reconnection using two-dimensional particle-in-cell simulations at various magnetizations \sigma. We find that the particle energy spectrum produced in reconnection quickly saturates as a hard power law that cuts off at \gamma~4\sigma, confirming previous work. Using particle tracing, we find that particle acceleration by the reconnection electric field in X-points determines the shape of the particle energy spectrum. By analyzing the current sheet structure, we show that physical cause of saturation is the spontaneous formation of secondary magnetic islands that can disrupt particle acceleration. By comparing the size of acceleration regions to the typical distance between disruptive islands, we show that the maximum Lorentz factor produced in reconnection is \gamma ~ 5 \sigma, which is very close to what we find in our particle energy spectra. We also show that the dynamic range in Lorentz factor of the power law spectrum in reconnection is < 40. The hardness of the power law combined with its narrow dynamic range implies that relativistic reconnection is capable of producing the hard narrowband flares observed in the Crab Nebula but has difficulty producing the softer broadband prompt GRB emission.

The maximum synchrotron burnoff limit of 160 MeV represents a fundamental limit to radiation resulting from electromagnetic particle acceleration in one-zone ideal plasmas. In magnetic reconnection, however, particle acceleration and radiation are decoupled because the electric field is larger than the magnetic field in the diffusion region. We carry out two-dimensional particle-in-cell simulations to determine the extent to which magnetic reconnection can produce synchrotron radiation above the burnoff limit. We use the test particle comparison (TPC) method to isolate the effects of cooling by comparing the trajectories and acceleration efficiencies of test particles incident on such a reconnection region with and without cooling them. We find that the cooled and uncooled particle trajectories are typically similar during acceleration in the reconnection region, and derive an effective limit on particle acceleration that is inversely proportional to the average magnetic field experienced by the particle during acceleration. Using the calculated distribution of this average magnetic field as a function of uncooled final particle energy, we find analytically that cooling does not affect power-law particle energy spectra except at energies far above the synchrotron burnoff limit. Finally, we compare fully cooled and uncooled simulations of reconnection, confirming that the synchrotron burnoff limit does not produce a cutoff in the particle energy spectrum. Our results indicate that the TPC method accurately predicts the effects of cooling on particle acceleration in relativistic reconnection, and that even far above the burnoff limit, the synchrotron energy of radiation produced in reconnection is not limited by cooling.

Relativistic reconnection has been invoked as a mechanism for particle acceleration in numerous astrophysical systems. According to idealised analytical models reconnection produces a bulk relativistic outflow emerging from the reconnection sites (X-points). The resulting radiation is therefore highly beamed. Using two-dimensional particle-in-cell (PIC) simulations, we investigate particle and radiation beaming, finding a very different picture. Instead of having a relativistic average bulk motion with isotropic electron velocity distribution in its rest frame, we find that the bulk motion of particles in X-points is similar to their Lorentz factor gamma, and the particles are beamed within about 5/gamma. On the way from the X-point to the magnetic islands, particles turn in the magnetic field, forming a fan confined to the current sheet. Once they reach the islands they isotropise after completing a full Larmor gyration and their radiation is not strongly beamed anymore. The radiation pattern at a given frequency depends on where the corresponding emitting electrons radiate their energy. Lower energy particles that cool slowly spend most of their time in the islands, and their radiation is not highly beamed. Only particles that quickly cool at the edge of the X-points generate a highly beamed fan-like radiation pattern. The radiation emerging from these fast cooling particles is above the burn-off limit (about 100 MeV in the overall rest frame of the reconnecting plasma.) This has significant implications for models of GRBs and AGNs that invoke beaming in that frame at much lower energies.

This review discusses the physics of magnetic reconnection, a process in which the magnetic field topology changes and magnetic energy is converted to kinetic energy, in pair plasmas in the relativistic regime. We focus on recent progress in the field driven by theory advances and the maturity of particle-in-cell codes. This work shows that fragmentation instabilities at the current sheet can play a critical role in setting the reconnection speed and affect the resulting particle acceleration, anisotropy, bulk flows, and radiation. Then, we discuss how this novel understanding of relativistic reconnection can be applied to high-energy astrophysical phenomena, with an emphasis on pulsars, pulsar wind nebulae, and active galactic nucleus jets.

The magnetorotational instability (MRI) is key physics in accretion disks and is widely considered to play some role in massive-star core collapse. Models of rotating massive stars naturally develop very strong shear at composition boundaries, a necessary condition for MRI instability, and the MRI is subject to triply-diffusive destabilizing effects in radiative regions. We have used the MESA stellar evolution code to compute magnetic effects due to the Spruit-Taylor mechanism and the MRI, separately and together, in a sample of massive star models. We find that the MRI can be active in the later stages of massive star evolution, leading to mixing effects that are not captured in models that neglect the MRI. The MRI and related magneto-rotational effects can move models of given ZAMS mass across "boundaries" from degenerate CO cores to degenerate O/Ne/Mg cores and from degenerate O/Ne/Mg cores to iron cores, thus affecting the final evolution and the physics of core collapse. The MRI acting alone can slow the rotation of the inner core in general agreement with the observed "initial" rotation rates of pulsars. The MRI analysis suggests that localized fields ~10^{12} G may exist at the boundary of the iron core. With both the ST and MRI mechanisms active in the 20 solar mass model, we find that the helium shell mixes entirely out into the envelope. Enhanced mixing could yield a population of yellow or even blue supergiant supernova progenitors that would not be standard SN IIP.

We calculate growth rates for nonaxisymmetric instabilities including the magnetorotational instability (MRI) throughout the Sun. We first derive a dispersion relation for nonaxisymmetric instability including the effects of shear, convective buoyancy, and three diffusivities (thermal conductivity, resistivity, and viscosity). We then use a solar model evolved with the stellar evolution code MESA and angular velocity profiles determined by Global Oscillations Network Group (GONG) helioseismology to determine the unstable modes present at each location in the Sun and the associated growth rates. The overall instability has unstable modes throughout the convection zone and also slightly below it at middle and high latitudes. It contains three classes of modes: large-scale hydrodynamic convective modes, large-scale hydrodynamic shear modes, and small-scale magnetohydrodynamic (MHD) shear modes, which may be properly called MRI modes. While large-scale convective modes are the most rapidly growing modes in most of the convective zone, MRI modes are important in both stably stratified and convectively unstable locations near the tachocline at colatitudes theta less than 53 degrees. Nonaxisymmetric MRI modes grow faster than the corresponding axisymmetric modes; for some poloidal magnetic fields, the nonaxisymmetric MRI growth rates are similar to the angular rotation frequency Omega, while axisymmetric modes are stabilized. We briefly discuss the saturation of the field produced by MRI modes, finding that the implied field at the base of the convective zone in the Sun is comparable to that derived based on dynamos active in the tachocline and that the saturation of field resulting from the MRI may be of importance even in the upper convection zone.

We investigate guide-field magnetic reconnection and particle acceleration in relativistic pair plasmas with three-dimensional particle-in-cell (PIC) simulations of a kinetic-scale current sheet in a periodic geometry at low magnetizations. The tearing instability is the dominant mode in the current sheet for all guide field strengths, while the linear kink mode is less important even without guide field. Oblique modes seem to be suppressed entirely. In its nonlinear evolution, the reconnection layer develops a network of interconnected and interacting magnetic flux ropes. As smaller flux ropes merge into larger ones, the reconnection layer evolves toward a three-dimensional, disordered state in which the resulting flux rope segments contain magnetic substructure on plasma skin depth scales. Embedded in the flux ropes, we detect spatially and temporally intermittent sites of dissipation reflected in peaks in the parallel electric field. Magnetic dissipation and particle acceleration persist until the end of the simulations, with simulations with higher magnetization and lower guide field strength exhibiting greater and faster energy conversion and particle energization. At the end of our largest simulation, the particle energy spectrum attains a tail extending to high Lorentz factors that is best modeled with a combination of two additional thermal components. We confirm that the primary energization mechanism is acceleration by the electric field in the X-line region. We discuss the implications of our results for macroscopic reconnection sites, and which of our results may be expected to hold in systems with higher magnetizations.

Supernovae (SNe) are thought to arise from two different physical processes. The cores of massive, short-lived stars undergo gravitational core collapse and typically eject a few solar masses during their explosion. These are thought to appear as as type Ib/c and II SNe, and are associated with young stellar populations. A type Ia SN is thought to arise from the thermonuclear detonation of a white dwarf star composed mainly of carbon and oxygen, whose mass approaches the Chandrasekhar limit. Such SNe are observed in both young and old stellar environments. Here we report our discovery of the faint type Ib SN 2005E in the halo of the nearby isolated galaxy, NGC 1032. The lack of any trace of recent star formation near the SN location (Fig. 1), and the very low derived ejected mass (~0.3 M_sun), argue strongly against a core-collapse origin for this event. Spectroscopic observations and the derived nucleosynthetic output show that the SN ejecta have high velocities and are dominated by helium-burning products, indicating that SN 2005E was neither a subluminous nor a regular SN Ia (Fig. 2). We have therefore found a new type of stellar explosion, arising from a low-mass, old stellar system, likely involving a binary with a primary white dwarf and a helium-rich secondary. The SN ejecta contain more calcium than observed in any known type of SN and likely additional large amounts of radioactive 44Ti. Such SNe may thus help resolve fundamental physical puzzles, extending from the composition of the primitive solar system and that of the oldest stars, to the Galactic production of positrons.

It has been observed that binary Kuiper Belt Objects (KBOs) exist contrary to theoretical expectations. Their creation presents problems to most current models. However, the inclusion of a third body (for example, one of the outer planets) may provide the conditions necessary for the formation of these objects. The presence of a third massive body not only helps to clear the primordial Kuiper Belt but can also result in long lived binary Kuiper belt objects. The gravitational interaction between the KBOs and the third body causes one of four effects; scattering into the Oort cloud, collisions with the growing protoplanets, formation of binary pairs, or creation of a single Kuiper belt object. Additionally, the initial location of the progenitors of the Kuiper belt objects also has a significant effect on binary formation.