Magnetic reconnection is thought to be the driver for many explosive phenomena in the universe. The energy release and particle acceleration during reconnection have been proposed as a mechanism for producing high-energy emissions and cosmic rays. We carry out two-and three-dimensional kinetic simulations to investigate relativistic magnetic reconnection and the associated particle acceleration. The simulations focus on electron-positron plasmas starting with a magnetically dominated, force-free current sheet (σ ≡ B 2 /(4πn e m e c 2 ) 1). For this limit, we demonstrate that relativistic reconnection is highly efficient at accelerating particles through a first-order Fermi process accomplished by the curvature drift of particles along the electric field induced by the relativistic flows. This mechanism gives rise to the formation of hard power-law spectra f ∝ (γ − 1) −p and approaches p = 1 for sufficiently large σ and system size. Eventually most of the available magnetic free energy is converted into nonthermal particle kinetic energy. An analytic model is presented to explain the key results and predict a general condition for the formation of power-law distributions. The development of reconnection in these regimes leads to relativistic inflow and outflow speeds and enhanced reconnection rates relative to non-relativistic regimes. In the three-dimensional simulation, the interplay between secondary kink and tearing instabilities leads to strong magnetic turbulence, but does not significantly change the energy conversion, reconnection rate, or particle acceleration. This study suggests that relativistic reconnection sites are strong sources of nonthermal particles, which may have important implications to a variety of high-energy astrophysical problems.

We suggest a new approach that could be used for modeling both the large-scale behavior of astrophysical jets and the magnetically dominated explosions in astrophysics. We describe a method for modeling the injection of magnetic fields and their subsequent evolution in a regime where the free energy is magnetically dominated. The injected magnetic fields, along with their associated currents, have both poloidal and toroidal components, and they are not force free. The dynamic expansion driven by the Lorentz force of the injected fields is studied using threedimensional ideal magnetohydrodynamic simulations. The generic behavior of magnetic field expansion, the interactions with the background medium, and the dependence on various parameters are investigated.

We investigate the observational signatures of super-Earths (i.e., planets withEarth-to-Neptunemass), which are the most common type of exoplanet discovered to date, in their natal disks of gas and dust. Combining two-fluid global hydrodynamics simulations with a radiative transfer code, we calculate the distributions of gas and of submillimeter-sized dust in a disk perturbed by a super-Earth, synthesizing images in near-infrared scattered light and the millimeter-wave thermal continuum for direct comparison with observations. In low-viscosity gas (  a -10 4 ), a super-Earth opens two annular gaps to either side of its orbit by the action of Lindblad torques. This double gap and its associated gas pressure gradients cause dust particles to be dragged by gas into three rings: one ring sandwiched between the two gaps, and two rings located at the gap edges farthest from the planet. Depending on thesystem parameters, additional rings may manifest for a single planet. A double gap located at tens of aufrom a host star in Taurus can be detected in the dust continuum by the Atacama Large Millimeter Array (ALMA) at an angular resolution of ∼  0. 03 after two hours of integration. Ring and gap features persist in a variety of background disk profiles, last for thousands of orbits, and change their relative positions and dimensions depending on the speed and direction of planet migration. Candidate double gaps have been observed by ALMA in systems such asHL Tau (D5 and D6) and TW Hya (at 37 and 43 au); we submit that each double gap is carved by one super-Earth in nearly inviscid gas.

Magnetic reconnection is a leading mechanism for dissipating magnetic energy and accelerating nonthermal particles in Poynting-flux dominated flows. In this letter, we investigate nonthermal particle acceleration during magnetic reconnection in a magnetically-dominated ion-electron plasma using fully kinetic simulations. For an ion-electron plasma with the total magnetization σ 0 = B 2 /(4πn(m i + m e )c 2 ), the magnetization for each species is σ i ∼ σ 0 and σ e ∼ (m i /m e )σ 0 , respectively. We have studied the magnetically dominated regime by varying σ e = 10 3 − 10 5 with initial ion and electron temperatures T i = T e = 5 − 20m e c 2 and mass ratio m i /m e = 1 − 1836. The results demonstrate that reconnection quickly establishes power-law energy distributions for both electrons and ions within several (2 − 3) light-crossing times. For the cases with periodic boundary conditions, the power-law index is 1 < s < 2 for both electrons and ions. The hard spectra limit the power-law energies for electrons and ions to be γ be ∼ σ e and γ bi ∼ σ i , respectively. The main acceleration mechanism is a Fermi-like acceleration through the drift motions of charged particles. When comparing the spectra for electrons and ions in momentum space, the spectral indices s p are identical as predicted in Fermi acceleration. We also find that the bulk flow can carry a significant amount of energy during the simulations. We discuss the implication of this study in the context of Poynting-flux dominated jets and pulsar winds especially the applications for explaining the nonthermal high-energy emissions.

While observations have suggested that power-law electron energy spectra are a common outcome of strong energy release during magnetic reconnection, e.g., in solar flares, kinetic simulations have not been able to provide definite evidence of power-laws in energy spectra of non-relativistic reconnection. By means of 3D large-scale fully kinetic simulations, we study the formation of power-law electron energy spectra in nonrelativistic low-β reconnection. We find that both the global spectrum integrated over the entire domain and local spectra within individual regions of the reconnection layer have power-law tails with a spectral index p ∼ 4 in the 3D simulation, which persist throughout the non-linear reconnection phase until saturation. In contrast, the spectrum in the 2D simulation rapidly evolves and quickly becomes soft. We show that 3D effects such as self-generated turbulence and chaotic magnetic field lines enable the transport of high-energy electrons across the reconnection layer and allow them to access several main acceleration regions. This leads to a sustained and nearly constant accel-