The hot loop structures in the solar corona can be well modeled by three dimensional magnetohydrodynamic simulations, where the corona is heated by field line braiding driven at the photosphere. To be able to reproduce the emission comparable to observations, one has to use realistic values for the Spitzer heat conductivity, which puts a large constraint on the time step of these simulations and make them therefore computationally expensive. Here, we present a non-Fourier description of the heat flux evolution, which allow us to speed up the simulations significantly. Together with the semi-relativistic Boris correction, we are able to limit the time step constraint of the Alfv\'en speed and speed up the simulations even further. We discuss the implementation of these two methods to the Pencil Code and present their implications on the time step, and the temperature structures, the ohmic heating rate and the emission in simulations of the solar corona. Using a non-Fourier description of the heat flux evolution together with the Boris correction, we can increase the time step of the simulation significantly without moving far away from the reference solution. However, for values of the Alfv\'en speed limit of 3 000 km/s and below, the simulation moves away from the reference solution und produces much higher temperatures and much structures with stronger emission.

We study the magnetic field and current structure associated with a coronal loop. Through this we investigate to what extent the assumptions of a force-free magnetic field break down and where they might be justified. We analyse a 3D MHD model of the solar corona in an emerging active region with the focus on the structure of the forming coronal loops. The lower boundary of this simulation is taken from a model of an emerging active region. As a consequence of the emerging magnetic flux and the horizontal motions at the surface a coronal loop forms self-consistently. We investigate the current density along magnetic field inside (and outside) this loop and study the magnetic and plasma properties in and around it. We find that the total current along the loop changes its sign from being antiparallel to parallel to the magnetic field. This is caused by the inclination of the loop together with the footpoint motion. Around the loop the currents form a complex non-force-free helical structure. This is directly related to a bipolar current structure at the loop footpoints at the base of the corona and a local reduction of the background magnetic field (i.e. outside the loop) caused by the plasma flow into and along the loop. The locally reduced magnetic pressure in the loop allows the loop to sustain a higher density, which is crucial for the emission in extreme UV. The acting of the flow on the magnetic field hosting the loop turns out to be also responsible for the observed squashing of the loop. The complex magnetic field and current system surrounding it can be modeled only in 3D MHD models where the magnetic field has to balance the plasma pressure. A 1D coronal loop model or a force-free extrapolation can not capture the current system and the complex interaction of the plasma and the magnetic field in the coronal loop, despite the fact that the loop is under low-$\beta$ conditions.

Aims. The goal is to employ a 3D magnetohydrodynamics (MHD) model including spectral synthesis to model the corona in an observed solar active region. This will allow us to judge the merits of the coronal heating mechanism built into the 3D model. Methods. Photospheric observations of the magnetic field and horizontal velocities in an active region are used to drive our coronal simulation from the bottom. The currents induced by this heat the corona through Ohmic dissipation. Heat conduction redistributes the energy that is lost in the end through optically thin radiation. Based on the MHD model, we synthesized profiles of coronal emission lines which can be directly compared to actual coronal observations of the very same active region. Results. In the synthesized model data we find hot coronal loops which host siphon flows or which expand and lose mass through draining. These synthesized loops are at the same location as and show similar dynamics in terms of Doppler shifts to the observed structures. This match is shown through a comparison with Hinode data as well as with 3D stereoscopic reconstructions of data from STEREO. Conclusions. The considerable match to the actual observations shows that the field-line braiding mechanism leading to the energy input in our corona provides the proper distribution of heat input in space and time. From this we conclude that in an active region the field-line braiding is the dominant heating process, at least at the spatial scales available to current observations.

Context. The structure and heating of coronal loops are investigated since decades. Established scaling laws relate fundamental quantities like the loop apex temperature, pressure, length, and the coronal heating. Aims. We test such scaling laws against a large-scale 3D MHD model of the Solar corona, which became feasible with nowadays high-performance computing. Methods. We drive an active region simulation a with photospheric observations and found strong similarities to the observed coronal loops in X-rays and EUV wavelength. A 3D reconstruction of stereoscopic observations showed that our model loops have a realistic spatial structure. We compare scaling laws to our model data extracted along an ensemble of field lines. Finally, we fit a new scaling law that represents well hot loops and also cooler structures, which was not possible before only based on observations. Results. Our model data gives some support for scaling laws that were established for hot and EUV-emissive coronal loops. For the RTV scaling law we find an offset to our model data, which can be explained by 1D considerations of a static loop with a constant heat input and conduction. With a fit to our model data we set up a new scaling law for the coronal heat input along magnetic field lines. Conclusions. RTV-like scaling laws were fitted to hot loops and therefore do not predict well the coronal heat input for cooler structures that are hardly observable. The basic differences between 1D and self-consistent 3D modeling account for deviations between our and earlier scaling laws. We also conclude that a heating mechanism by MHD-turbulent dissipation within a braided flux tube would heat the corona stronger than consistent with our model corona.

Context. We have conducted a 3D MHD simulation of the solar corona above an active region (AR) in full scale and high resolution, which shows coronal loops, and plasma flows within them, similar to observations. Aims. We want to find the connection between the photospheric energy input by field-line braiding with the coronal energy conversion by Ohmic dissipation of induced currents. Methods. To this end we compare the coronal energy input and dissipation within our simulation domain above different fields of view, e.g. for a small loops system in the AR core. We also choose an ensemble of field lines to compare, e.g., the magnetic energy input to the heating per particle along these field lines. Results. We find an enhanced Ohmic dissipation of currents in the corona above areas that also have enhanced upwards-directed Poynting flux. These regions coincide with the regions where hot coronal loops within the AR core are observed. The coronal density plays a role in estimating the coronal temperature due to the generated heat input. A minimum flux density of about 200 Gauss is needed in the photosphere to heat a field line to coronal temperatures of about 1 MK. Conclusions. This suggests that the field-line braiding mechanism provides the coronal energy input and that the Ohmic dissipation of induced currents dominates the coronal heating mechanism.

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.

We conducted a high-resolution numerical simulation of the solar corona above a stable active region. The aim is to test the field-line braiding mechanism for a sufficient coronal energy input. We also check the applicability of scaling laws for coronal loop properties like the temperature and density. Our 3D-MHD model is driven from below by Hinode observations of the photosphere, in particular a high-cadence time series of line-of-sight magnetograms and horizontal velocities derived from the magnetograms. This driving applies stress to the magnetic field and thereby delivers magnetic energy into the corona, where currents are induced that heat the coronal plasma by Ohmic dissipation. We compute synthetic coronal emission that we directly compare to coronal observations of the same active region taken by Hinode. In the model, coronal loops form at the same places as they are found in coronal observations. Even the shapes of the synthetic loops in 3D space match those found from a stereoscopic reconstruction based on STEREO spacecraft data. Some loops turn out to be slightly over-dense in the model, as expected from observations. This shows that the spatial and temporal distribution of the Ohmic heating produces the structure and dynamics of a coronal loops system close to what is found in observations.

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.

We will use new data from the High-resolution Coronal Imager (Hi-C) with unprecedented spatial resolution of the solar corona to investigate the structure of coronal loops down to 0.2 arcsec. During a rocket flight Hi-C provided images of the solar corona in a wavelength band around 193 A that is dominated by emission from Fe XII showing plasma at temperatures around 1.5 MK. We analyze part of the Hi-C field-of-view to study the smallest coronal loops observed so far and search for the a possible sub-structuring of larger loops. We find tiny 1.5 MK loop-like structures that we interpret as miniature coronal loops. These have length of the coronal segment above the chromosphere of only about 1 Mm and a thickness of less than 200 km. They could be interpreted as the coronal signature of small flux tubes breaking through the photosphere with a footpoint distance corresponding to the diameter of a cell of granulation. We find loops that are longer than 50 Mm to have a diameter of about 2 arcsec or 1.5 Mm, consistent with previous observations. However, Hi-C really resolves these loops with some 20 pixels across the loop. Even at this greatly improved spatial resolution the large loops seem to have no visible sub-structure. Instead they show a smooth variation in cross-section. The fact that the large coronal loops do not show a sub-structure at the spatial scale of 0.1 arcsec per pixel implies that either the densities and temperatures are smoothly varying across these loops or poses an upper limit on the diameter of strands the loops might be composed of. We estimate that strands that compose the 2 arcsec thick loop would have to be thinner than 15 km. The miniature loops we find for the first time pose a challenge to be properly understood in terms of modeling.

We investigate the difference in the spatial distribution of the energy input for parametrizations of different mechanisms to heat the corona of the Sun and possible impacts on the coronal emission. We use a 3D MHD model of a solar active region as a reference and compare the Ohmic-type heating in this model to parametrizations for alternating current (AC) and direct current (DC) heating models, in particular, we use Alfven wave and MHD turbulence heating. We extract the quantities needed for these two parametrizations from the reference model and investigate the spatial distribution of the heat input in all three cases, globally and along individual field lines. To study differences in the resulting coronal emission we employ 1D loop models with a prescribed heat input based on the heating rate we extracted along a bundle of field lines. On average, all heating implementations show a roughly drop of the heating rate with height. This also holds for individual field lines. While all mechanism show a concentration of the energy input towards the low parts of the atmosphere, for individual field lines the concentration towards the footpoints is much stronger for the DC mechanisms than for the Alfven wave AC case. In contrast, the AC model gives a stronger concentration of the emission towards the footpoints. This is because the more homogeneous distribution of the energy input leads to higher coronal temperatures and a more extended transition region. The significant difference in the concentration of the heat input towards the foot points for the AC and DC mechanisms, and the pointed difference in the spatial distribution of the coronal emission for these cases shows that the two mechanisms should be discriminable by observations. Before drawing final conclusions, these parametrizations should be implemented in new 3D models in a more self-consistent way.

Context. We investigate the statistics of the spatial and temporal distribution of the coronal heating in a three-dimensional magneto- hydrodynamical (3D MHD) model. The model describes the temporal evolution of the corona above an observed active region. The model is driven by photospheric granular motions which braid the magnetic field lines. This induces currents and their dissipation heats the plasma. We evaluate the transient heating as subsequent heating events and analyze their statistics. The results are then interpreted in the context of observed flare statistics and coronal heating mechanisms. Methods. To conduct the numerical experiment we use a high order finite difference code which solves the partial differential equations for the conservation of mass, the momentum and energy balance, and the induction equation. The energy balance includes the Spitzer heat conduction and the optical thin radiative loss in the corona. Results. The temporal and spatial distribution of the Ohmic heating in the 3D MHD model follow a power law and can therefore be explained by system in a self-organized critical state. The slopes of the power law are similar to the results based on flare observations. We find that the corona is heated foot point dominated and the coronal heating is dominated by events called nanoflares with energies on the order of 1017 J.

The corona of the Sun is dominated by emission from loop-like structures. When observed in X-ray or extreme ultraviolet emission, these million K hot coronal loops show a more or less constant cross section. In this study we show how the interplay of heating, radiative cooling, and heat conduction in an expanding magnetic structure can explain the observed constant cross section. We employ a three-dimensional magnetohydrodynamics (3D MHD) model of the corona. The heating of the coronal plasma is the result of braiding of the magnetic field lines through footpoint motions and subsequent dissipation of the induced currents. From the model we synthesize the coronal emission, which is directly comparable to observations from, e.g., the Atmospheric Imaging Assembly on the Solar Dynamics Observatory (AIA/SDO). We find that the synthesized observation of a coronal loop seen in the 3D data cube does match actually observed loops in count rate and that the cross section is roughly constant, as observed. The magnetic field in the loop is expanding and the plasma density is concentrated in this expanding loop; however, the temperature is not constant perpendicular to the plasma loop. The higher temperature in the upper outer parts of the loop is so high that this part of the loop is outside the contribution function of the respective emission line(s). In effect, the upper part of the plasma loop is not bright and thus the loop actually seen in coronal emission appears to have a constant width. From this we can conclude that the underlying field-line-braiding heating mechanism provides the proper spatial and temporal distribution of the energy input into the corona --- at least on the observable scales.

Condensations in the more than 10^6 K hot corona of the Sun are commonly observed in the extreme ultraviolet (EUV). While their contribution to the total solar EUV radiation is still a matter of debate, these condensations certainly provide a valuable tool for studying the dynamic response of the corona to the heating processes. We investigate different distributions of energy input in time and space to investigate which process is most relevant for understanding these coronal condensations. For a comparison to observations we synthesize EUV emission from a time-dependent, one-dimensional model for coronal loops, where we employ two heating scenarios: simply shutting down the heating and a model where the heating is very concentrated at the loop footpoints, while keeping the total heat input constant. The heating off/on model does not lead to significant EUV count rates that one observes with SDO/AIA. In contrast, the concentration of the heating near the footpoints leads to thermal non-equilibrium near the loop top resulting in the well-known catastrophic cooling. This process gives a good match to observations of coronal condensations. This shows that the corona needs a steady supply of energy to support the coronal plasma, even during coronal condensations. Otherwise the corona would drain very fast, too fast to even form a condensation.

We investigate the processes that lead to the formation, ejection and fall of a confined plasma ejection that was observed in a numerical experiment of the solar corona. By quantifying physical parameters such as mass, velocity, and orientation of the plasma ejection relative to the magnetic field, we provide a description of the nature of this particular phenomenon. The time-dependent three-dimensional magnetohydrodynamic (3D MHD) equations are solved in a box extending from the chromosphere to the lower corona. The plasma is heated by currents that are induced through field line braiding as a consequence of photospheric motions. Spectra of optically thin emission lines in the extreme ultraviolet range are synthesized, and magnetic field lines are traced over time. Following strong heating just above the chromosphere, the pressure rapidly increases, leading to a hydrodynamic explosion above the upper chromosphere in the low transition region. The explosion drives the plasma, which needs to follow the magnetic field lines. The ejection is then moving more or less ballistically along the loop-like field lines and eventually drops down onto the surface of the Sun. The speed of the ejection is in the range of the sound speed, well below the Alfven velocity. The plasma ejection is basically a hydrodynamic phenomenon, whereas the rise of the heating rate is of magnetic nature. The granular motions in the photosphere lead (by chance) to a strong braiding of the magnetic field lines at the location of the explosion that in turn is causing strong currents which are dissipated. Future studies need to determine if this process is a ubiquitous phenomenon on the Sun on small scales. Data from the Atmospheric Imaging Assembly on the Solar Dynamics Observatory (AIA/SDO) might provide the relevant information.

The origin of solar transition region redshifts is not completely understood. Current research is addressing this issue by investigating three-dimensional magneto-hydrodynamic models that extend from the photosphere to the corona. By studying the average properties of emission line profiles synthesized from the simulation runs and comparing them to observations with present-day instrumentation, we investigate the origin of mass flows in the solar transition region and corona. Doppler shifts were determined from the emission line profiles of various extreme-ultraviolet emission lines formed in the range of $T=10^4-10^6$ K. Plasma velocities and mass flows were investigated for their contribution to the observed Doppler shifts in the model. In particular, the temporal evolution of plasma flows along the magnetic field lines was analyzed. Comparing observed vs. modeled Doppler shifts shows a good correlation in the temperature range $\log(T$/[K])=4.5-5.7, which is the basis of our search for the origin of the line shifts. The vertical velocity obtained when weighting the velocity by the density squared is shown to be almost identical to the corresponding Doppler shift. Therefore, a direct comparison between Doppler shifts and the model parameters is allowed. A simple interpretation of Doppler shifts in terms of mass flux leads to overestimating the mass flux. Upflows in the model appear in the form of cool pockets of gas that heat up slowly as they rise. Their low temperature means that these pockets are not observed as blueshifts in the transition region and coronal lines. For a set of magnetic field lines, two different flow phases could be identified. The coronal part of the field line is intermittently connected to subjacent layers of either strong or weak heating, leading either to mass flows into the loop or to the draining of the loop.

We investigate the spatial and temporal evolution of the heating of the corona of a cool star such as our Sun in a three-dimensional magneto-hydrodynamic (3D MHD) model. We solve the 3D MHD problem numerically in a box representing part of the (solar) corona. The energy balance includes Spitzer heat conduction along the magnetic field and optically thin radiative losses. The self-consistent heating mechanism is based on the braiding of magnetic field lines rooted in the convective photosphere. Magnetic stress induced by photospheric motions leads to currents in the atmosphere which heat the corona through Ohmic dissipation. While the horizontally averaged quantities, such as heating rate, temperature or density, are relatively constant in time, the simulated corona is highly variable and dynamic, on average reaching temperatures and densities as found in observations. The strongest heating per particle is found in the transition region from the chromosphere to the corona. The heating is concentrated in current sheets roughly aligned with the magnetic field and is transient in time and space. This supports the idea that numerous small heating events heat the corona, often referred to a nanoflares.

We study extreme-ultraviolet emission line spectra derived from three-dimensional magnetohydrodynamic models of structures in the corona. In order to investigate the effects of increased magnetic activity at photospheric levels in a numerical experiment, a much higher magnetic flux density is applied at photospheric levels as compared to the Sun. Thus, we can expect our results to highlight the differences between the Sun and more active, but still solar-like stars. We discuss signatures seen in extreme-ultraviolet emission lines synthesized from these models and compare them to signatures found in the spatial distribution and temporal evolution of Doppler shifts in lines formed in the transition region and corona. This is of major interest to test the quality of the underlying magnetohydrodynamic model to heat the corona, i.e. currents in the corona driven by photospheric motions (flux braiding).