Reconnection current sheet structure in a turbulent medium

In the presence of turbulence, magnetic field lines lose their dynamical identity and particles entrained on field lines diffuse through space at a rate determined by the amplitude of the turbulence. In previous work (Lazarian and Vishniac, 1999; Kowal et al., 2009; Eyink et al., 2011) we showed that this leads to reconnection speeds which are independent of resistivity. In particular, in Kowal et al. (2009) we showed that numerical simulations were consistent with the predictions of this model. Here we examine the structure of the current sheet in simulations of turbulent reconnection. Laminar flows consistent with the Sweet-Parker reconnection model produce very thin and well ordered currents sheets. On the other hand, the simulations of Kowal et al. (2009) show a strongly disordered state even for relatively low levels of turbulence. Comparing data cubes with and without reconnection, we find that large scale field reversals are the cumulative effect of many individual eddies, each of which has magnetic properties which are not very different from turbulent eddies in a homogeneous background. This implies that the properties of stationary and homogeneous MHD turbulence are a reasonable guide to understanding turbulence during large scale magnetic reconnection events. In addition, dissipation and high energy particle acceleration during reconnection events take place over a macroscopic volume, rather than being confined to a narrow zone whose properties depend on microscopic transport coefficients.

Reconnection studies under different types of turbulence driving

We study a model of fast magnetic reconnection in the presence of weak turbulence proposed by Lazarian and Vishniac (1999) using three-dimensional direct numerical simulations. The model has been already successfully tested in Kowal et al. (2009) confirming the dependencies of the reconnection speed V-rec on the turbulence injection power P-inj and the injection scale l(inj) expressed by a constraint V-rec similar to P(inj)(1/2)l(inj)(3/4)and no observed dependency on Ohmic resistivity. In Kowal et al. (2009), in order to drive turbulence, we injected velocity fluctuations in Fourier space with frequencies concentrated around k(inj) = 1/l(inj), as described in Alvelius (1999). In this paper, we extend our previous studies by comparing fast magnetic reconnection under different mechanisms of turbulence injection by introducing a new way of turbulence driving. The new method injects velocity or magnetic eddies with a specified amplitude and scale in random locations directly in real space. We provide exact relations between the eddy parameters and turbulent power and injection scale. We performed simulations with new forcing in order to study turbulent power and injection scale dependencies. The results show no discrepancy between models with two different methods of turbulence driving exposing the same scalings in both cases. This is in agreement with the Lazarian and Vishniac (1999) predictions. In addition, we performed a series of models with varying viscosity nu. Although Lazarian and Vishniac (1999) do not provide any prediction for this dependence, we report a weak relation between the reconnection speed with viscosity, V-rec similar to nu(-1/4).

Magnetic flux transport by turbulent reconnection in astrophysical flows

The role of magnetohydrodynamics (MHD) turbulence in astrophysical environments is still highly debated. An important question that permeates this debate is the transport of magnetic flux. This is particularly important, for instance, in the context of star formation. When clouds collapse gravitationally to form stars, there must be some magnetic flux transport. Otherwise, the newborn stars would have magnetic fields several orders of magnitude larger than the observed ones. Also, the magnetic flux that is dragged in the late stages of the formation of a star can remove all the rotational support from the accretion disc that grows around the protostar. The efficiency of the mechanism that is often invoked to allow transport of magnetic fields at different stages of star formation, namely ambipolar diffusion, has recently been put in check. We discuss here an alternative mechanism for magnetic flux transport which is based on turbulent fast magnetic reconnection. We review recent results from three-dimensional MHD numerical simulations that indicate that this mechanism is very efficient in decoupling and transporting magnetic flux from the inner denser regions to the outskirts of collapsing clouds at different stages of star formation. We discuss this mechanism also in the context of dynamo processes and speculate that it can play a role both in solar dynamo and in accretion disc dynamo processes.

THE ROLE OF TURBULENT MAGNETIC RECONNECTION IN THE FORMATION OF ROTATIONALLY SUPPORTED PROTOSTELLAR DISKS

The formation of protostellar disks out of molecular cloud cores is still not fully understood. Under ideal MHD conditions, the removal of angular momentum from the disk progenitor by the typically embedded magnetic field may prevent the formation of a rotationally supported disk during the main protostellar accretion phase of low-mass stars. This has been known as the magnetic braking problem and the most investigated mechanism to alleviate this problem and help remove the excess of magnetic flux during the star formation process, the so-called ambipolar diffusion (AD), has been shown to be not sufficient to weaken the magnetic braking at least at this stage of the disk formation. In this work, motivated by recent progress in the understanding of magnetic reconnection in turbulent environments, we appeal to the diffusion of magnetic field mediated by magnetic reconnection as an alternative mechanism for removing magnetic flux. We investigate numerically this mechanism during the later phases of the protostellar disk formation and show its high efficiency. By means of fully three-dimensional MHD simulations, we show that the diffusivity arising from turbulent magnetic reconnection is able to transport magnetic flux to the outskirts of the disk progenitor at timescales compatible with the collapse, allowing the formation of a rotationally supported disk around the protostar of dimensions similar to 100 AU, with a nearly Keplerian profile in the early accretion phase. Since MHD turbulence is expected to be present in protostellar disks, this is a natural mechanism for removing magnetic flux excess and allowing the formation of these disks. This mechanism dismisses the necessity of postulating a hypothetical increase of the ohmic resistivity as discussed in the literature. Together with our earlier work which showed that magnetic flux removal from molecular cloud cores is very efficient, this work calls for reconsidering the relative role of AD in the processes of star and planet formation.

Particle Acceleration in Turbulence and Weakly Stochastic Reconnection

In this Letter we analyze the energy distribution evolution of test particles injected in three dimensional (3D) magnetohydrodynamic (MHD) simulations of different magnetic reconnection configurations. When considering a single Sweet-Parker topology, the particles accelerate predominantly through a first-order Fermi process, as predicted in [3] and demonstrated numerically in [8]. When turbulence is included within the current sheet, the acceleration rate is highly enhanced, because reconnection becomes fast and independent of resistivity [4,11] and allows the formation of a thick volume filled with multiple simultaneously reconnecting magnetic fluxes. Charged particles trapped within this volume suffer several head-on scatterings with the contracting magnetic fluctuations, which significantly increase the acceleration rate and results in a first-order Fermi process. For comparison, we also tested acceleration in MHD turbulence, where particles suffer collisions with approaching and receding magnetic irregularities, resulting in a reduced acceleration rate. We argue that the dominant acceleration mechanism approaches a second order Fermi process in this case.

MHD numerical simulations of colliding winds in massive binary systems - I. Thermal versus non-thermal radio emission

In the past few decades detailed observations of radio and X-ray emission from massive binary systems revealed a whole new physics present in such systems. Both thermal and non-thermal components of this emission indicate that most of the radiation at these bands originates in shocks. O and B-type stars and WolfRayet (WR) stars present supersonic and massive winds that, when colliding, emit largely due to the freefree radiation. The non-thermal radio and X-ray emissions are due to synchrotron and inverse Compton processes, respectively. In this case, magnetic fields are expected to play an important role in the emission distribution. In the past few years the modelling of the freefree and synchrotron emissions from massive binary systems have been based on purely hydrodynamical simulations, and ad hoc assumptions regarding the distribution of magnetic energy and the field geometry. In this work we provide the first full magnetohydrodynamic numerical simulations of windwind collision in massive binary systems. We study the freefree emission characterizing its dependence on the stellar and orbital parameters. We also study self-consistently the evolution of the magnetic field at the shock region, obtaining also the synchrotron energy distribution integrated along different lines of sight. We show that the magnetic field in the shocks is larger than that obtained when the proportionality between B and the plasma density is assumed. Also, we show that the role of the synchrotron emission relative to the total radio emission has been underestimated.

Non-linear magnetohydrodynamic simulations of density evolution in Tore Supra sawtoothing plasmas

The plasma density evolution in sawtooth regime on the Tore Supra tokamak is analyzed. The density is measured using fast-sweeping X-mode reflectometry which allows tomographic reconstructions. There is evidence that density is governed by the perpendicular electric flows, while temperature evolution is dominated by parallel diffusion. Postcursor oscillations sometimes lead to the formation of a density plateau, which is explained in terms of convection cells associated with the kink mode. A crescent-shaped density structure located inside q = 1 is often visible just after the crash and indicates that some part of the density withstands the crash. 3D full MHD nonlinear simulations with the code XTOR-2F recover this structure and show that it arises from the perpendicular flows emerging from the reconnection layer. The proportion of density reinjected inside the q = 1 surface is determined, and the implications in terms of helium ash transport are discussed. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4766893]

Secondary nontwist phenomena in area-preserving maps

Phenomena as reconnection scenarios, periodic-orbit collisions, and primary shearless tori have been recognized as features of nontwist maps. Recently, these phenomena and secondary shearless tori were analytically predicted for generic maps in the neighborhood of the tripling bifurcation of an elliptic fixed point. In this paper, we apply a numerical procedure to find internal rotation number profiles that highlight the creation of periodic orbits within islands of stability by a saddle-center bifurcation that emerges out a secondary shearless torus. In addition to the analytical predictions, our numerical procedure applied to the twist and nontwist standard maps reveals that the atypical secondary shearless torus occurs not only near a tripling bifurcation of the fixed point but also near a quadrupling bifurcation. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4750040]

Flares and variability from Sagittarius A(star): five nights of simultaneous multi-wavelength observations

Aims. We report on simultaneous observations and modeling of mid-infrared (MIR), near-infrared (NIR), and submillimeter (sub-mm) emission of the source Sgr A * associated with the supermassive black hole at the center of our Galaxy. Our goal was to monitor the activity of Sgr A* at different wavelengths in order to constrain the emitting processes and gain insight into the nature of the close environment of Sgr A*. Methods. We used the MIR instrument VISIR in the BURST imaging mode, the adaptive optics assisted NIR camera NACO, and the sub-mm antenna APEX to monitor Sgr A* over several nights in July 2007. Results. The observations reveal remarkable variability in the NIR and sub-mm during the five nights of observation. No source was detected in the MIR, but we derived the lowest upper limit for a flare at 8.59 mu m (22.4 mJy with A(8.59 mu m) = 1.6 +/- 0.5). This observational constraint makes us discard the observed NIR emission as coming from a thermal component emitting at sub-mm frequencies. Moreover, comparison of the sub-mm and NIR variability shows that the highest NIR fluxes (flares) are coincident with the lowest sub-mm levels of our five-night campaign involving three flares. We explain this behavior by a loss of electrons to the system and/or by a decrease in the magnetic field, as might conceivably occur in scenarios involving fast outflows and/or magnetic reconnection.

Motor recovery and synaptic preservation after ventral root avulsion and repair with a fibrin sealant derived from snake venom

Background: Ventral root avulsion is an experimental model of proximal axonal injury at the central/peripheral nervous system interface that results in paralysis and poor clinical outcome after restorative surgery. Root reimplantation may decrease neuronal degeneration in such cases. We describe the use of a snake venom-derived fibrin sealant during surgical reconnection of avulsed roots at the spinal cord surface. The present work investigates the effects of this fibrin sealant on functional recovery, neuronal survival, synaptic plasticity, and glial reaction in the spinal motoneuron microenvironment after ventral root reimplantation. Methodology/Principal Findings: Female Lewis rats (7 weeks old) were subjected to VRA and root replantation. The animals were divided into two groups: 1) avulsion only and 2) replanted roots with fibrin sealant derived from snake venom. Post-surgical motor performance was evaluated using the CatWalk system twice a week for 12 weeks. The rats were sacrificed 12 weeks after surgery, and their lumbar intumescences were processed for motoneuron counting and immunohistochemistry (GFAP, Iba-1 and synaptophysin antisera). Array based qRT-PCR was used to evaluate gene regulation of several neurotrophic factors and receptors as well as inflammatory related molecules. The results indicated that the root reimplantation with fibrin sealant enhanced motor recovery, preserved the synaptic covering of the motoneurons and improved neuronal survival. The replanted group did not show significant changes in microglial response compared to VRA-only. However, the astroglial reaction was significantly reduced in this group. Conclusions/Significance: In conclusion, the present data suggest that the repair of avulsed roots with snake venom fibrin glue at the exact point of detachment results in neuroprotection and preservation of the synaptic network at the microenvironment of the lesioned motoneurons. Also such procedure reduced the astroglial reaction and increased mRNA levels to neurotrophins and anti-inflammatory cytokines that may in turn, contribute to improving recovery of motor function.