A model for the formation of the active region corona driven by magnetic flux emergence

We present the first model that couples the formation of the corona of a solar active region to a model of the emergence of a sunspot pair. This allows us to study when, where, and why active region loops form, and how they evolve. We use a 3D radiation MHD simulation of the emergence of an active region through the upper convection zone and the photosphere as a lower boundary for a 3D MHD coronal model. The latter accounts for the braiding of the magnetic fieldlines, which induces currents in the corona heating up the plasma. We synthesize the coronal emission for a direct comparison to observations. Starting with a basically field-free atmosphere we follow the filling of the corona with magnetic field and plasma. Numerous individually identifiable hot coronal loops form, and reach temperatures well above 1 MK with densities comparable to observations. The footpoints of these loops are found where small patches of magnetic flux concentrations move into the sunspots. The loop formation is triggered by an increase of upwards-directed Poynting flux at their footpoints in the photosphere. In the synthesized EUV emission these loops develop within a few minutes. The first EUV loop appears as a thin tube, then rises and expands significantly in the horizontal direction. Later, the spatially inhomogeneous heat input leads to a fragmented system of multiple loops or strands in a growing envelope.Comment: 13 pages, 10 figures, accepted to publication in A&

Magnetic Jam in the Corona of the Sun

The outer solar atmosphere, the corona, contains plasma at temperatures of more than a million K, more than 100 times hotter that solar surface. How this gas is heated is a fundamental question tightly interwoven with the structure of the magnetic field in the upper atmosphere. Conducting numerical experiments based on magnetohydrodynamics we account for both the evolving three-dimensional structure of the atmosphere and the complex interaction of magnetic field and plasma. Together this defines the formation and evolution of coronal loops, the basic building block prominently seen in X-rays and extreme ultraviolet (EUV) images. The structures seen as coronal loops in the EUV can evolve quite differently from the magnetic field. While the magnetic field continuously expands as new magnetic flux emerges through the solar surface, the plasma gets heated on successively emerging fieldlines creating an EUV loop that remains roughly at the same place. For each snapshot the EUV images outline the magnetic field, but in contrast to the traditional view, the temporal evolution of the magnetic field and the EUV loops can be different. Through this we show that the thermal and the magnetic evolution in the outer atmosphere of a cool star has to be treated together, and cannot be simply separated as done mostly so far.Comment: Final version published online on 27 April 2015, Nature Physics 12 pages and 8 figure

A kpc-scale X-ray jet in the BL Lac source S5 2007+777

X-ray jets in AGN are commonly observed in FRII and FRI radio-galaxies, but rarely in BL Lacs, most probably due to their orientation close to the line of sight and the ensuing foreshortening effects. Only three BL Lacs are known so far to contain a kpc-scale X-ray jet. In this paper, we present the evidence for the existence of a fourth extended X-ray jet in the classical radio-selected source S5 2007+777, which for its hybrid FRI/II radio morphology has been classified as a HYMOR (HYbrid MOrphology Radio source). Our Chandra ACIS-S observations of this source revealed an X-ray counterpart to the 19"-long radio jet. Interestingly, the X-ray properties of the kpc-scale jet in S5 2007+777 are very similar to those observed in FRII jets. First, the X-ray morphology closely mirrors the radio one, with the X-rays being concentrated in the discrete radio knots. Second, the X-ray continuum of the jet/brightest knot is described by a very hard power law, with photon index Gamma_x~1, although the uncertainties are large. Third, the optical upper limit from archival HST data implies a concave radio-to-X-ray SED. If the X-ray emission is attributed to IC/CMB with equipartition, strong beaming (delta=13) is required, implying a very large scale (Mpc) jet. The beaming requirement can be somewhat relaxed assuming a magnetic field lower than equipartition. Alternatively, synchrotron emission from a second population of very high-energy electrons is viable. Comparison to other HYMOR jets detected with Chandra is discussed, as well as general implications for the origin of the FRI/II division.Comment: Accepted for publication in ApJ, 19 pages, 3 figure

A simplified model of the source channel of the Leksell Gamma Knife(R)^(R): testing multisource configurations with PENELOPE

A simplification of the source channel geometry of the Leksell Gamma Knife$^{\circledR}$, recently proposed by the authors and checked for a single source configuration (Al-Dweri et al 2004), has been used to calculate the dose distributions along the $x$, $y$ and $z$ axes in a water phantom with a diameter of 160~mm, for different configurations of the Gamma Knife including 201, 150 and 102 unplugged sources. The code PENELOPE (v. 2001) has been used to perform the Monte Carlo simulations. In addition, the output factors for the 14, 8 and 4~mm helmets have been calculated. The results found for the dose profiles show a qualitatively good agreement with previous ones obtained with EGS4 and PENELOPE (v. 2000) codes and with the predictions of GammaPlan$^{\circledR}$. The output factors obtained with our model agree within the statistical uncertainties with those calculated with the same Monte Carlo codes and with those measured with different techniques. Owing to the accuracy of the results obtained and to the reduction in the computational time with respect to full geometry simulations (larger than a factor 15), this simplified model opens the possibility to use Monte Carlo tools for planning purposes in the Gamma Knife$^{\circledR}$.Comment: 13 pages, 8 figures, 5 table

Jets from Sub-Parsec to Kiloparsec Scales: A Physical Connection

The Chandra discovery of bright X-ray emission from kpc-scale jets allows insight into the physical parameters of the jet flow at large scale. At the opposite extreme, extensive studies of the inner relativistic jets in Blazars with multiwavelength observations, yield comparable information on sub-parsec scales. In the framework of simple radiation models for the emission regions we compare the physical parameters of jets on these two very different scales in the only two well studied Blazars for which large-scale emission has been resolved by Chandra. Notably, we find that the relativistic Doppler factors and powers derived independently at the two scales are consistent, suggesting that the jet does not suffer severe deceleration or dissipation. Moreover the internal equipartition pressures in the inner jet and in the external X-ray bright knots scale inversely with the jet cross section as expected in the simple picture of a freely expanding jet in equipartition.Comment: 4 figures, accepted by Ap

The origin of the reversed granulation in the solar photosphere

We study the structure and reveal the physical nature of the reversed granulation pattern in the solar photosphere by means of 3-dimensional radiative hydrodynamics simulations. We used the MURaM code to obtain a realistic model of the near-surface layers of the convection zone and the photosphere. The pattern of horizontal temperature fluctuations at the base of the photosphere consists of relatively hot granular cells bounded by the cooler intergranular downflow network. With increasing height in the photosphere, the amplitude of the temperature fluctuations diminishes. At a height of z=130-140 km in the photosphere, the pattern of horizontal temperature fluctuations reverses so that granular regions become relatively cool compared to the intergranular network. Detailed analysis of the trajectories of fluid elements through the photosphere reveal that the motion of the fluid is non-adiabatic, owing to strong radiative cooling when approaching the surface of optical depth unity followed by reheating by the radiation field from below. The temperature structure of the photosphere results from the competition between expansion of rising fluid elements and radiative heating. The former acts to lower the temperature of the fluid whereas the latter acts to increase it towards the radiative equilibrium temperature with a net entropy gain. After the fluid overturns and descends towards the convection zone, radiative energy loss again decreases the entropy of the fluid. Radiative heating and cooling of fluid elements that penetrate into the photosphere and overturn do not occur in equal amounts. The imbalance in the cumulative heating and cooling of these fluid elements is responsible for the reversal of temperature fluctuations with respect to height in the photosphere

Thermal Diagnostics with the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory: A Validated Method for Differential Emission Measure Inversions

We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). The method yields positive definite DEM solutions by solving a linear program. This method has been validated against a diverse set of thermal models of varying complexity and realism. These include (1) idealized gaussian DEM distributions, (2) 3D models of NOAA Active Region 11158 comprising quasi-steady loop atmospheres in a non-linear force-free field, and (3) thermodynamic models from a fully-compressible, 3D MHD simulation of AR corona formation following magnetic flux emergence. We then present results from the application of the method to AIA observations of Active Region 11158, comparing the region's thermal structure on two successive solar rotations. Additionally, we show how the DEM inversion method can be adapted to simultaneously invert AIA and XRT data, and how supplementing AIA data with the latter improves the inversion result. The speed of the method allows for routine production of DEM maps, thus facilitating science studies that require tracking of the thermal structure of the solar corona in time and space.Comment: 21 pages, 18 figures, accepted for publication in Ap