Observational and theoretical evidence suggests that beams of accelerated particles are produced in flaring events of all sizes in the solar atmosphere, from X-class flares to nanoflares. Current models of such particles in flaring loops assume an isolated 1D atmosphere. A more realistic environment for modelling accelerated particles can be provided by 3D radiative magnetohydrodynamics codes. Here we present a simple model for particle acceleration and propagation in the context of a 3D simulation of the quiet solar atmosphere spanning from convection zone to corona. We then examine the additional transport of energy introduced by the particle beams. Locations of particle acceleration associated with magnetic reconnection are identified by detecting changes in magnetic topology. At each location, the parameters of the accelerated particle distribution are estimated from local conditions. The particle distributions are then propagated along the magnetic field, and the energy deposition due to Coulomb collisions with the ambient plasma is computed. We find that particle beams originate in extended acceleration regions distributed across the corona. Upon reaching the transition region, they converge and produce strands of intense heating penetrating the chromosphere. Within these strands, beam heating consistently dominates conductive heating below the bottom of the transition region. This indicates that particle beams will qualitatively alter the energy transport even outside of active regions.

The presence of magnetic field is crucial in the transport of energy through the solar atmosphere. Recent ground-based and space-borne observations of the quiet Sun have revealed that magnetic field accumulates at photospheric heights, via a local dynamo or from small-scale flux emergence events. However, most of this small-scale magnetic field may not expand into the chromosphere due to the entropy drop with height at the photosphere. Here we present a study that uses a high resolution 3D radiative MHD simulation of the solar atmosphere with non-grey and non-LTE radiative transfer and thermal conduction along the magnetic field to reveal that: 1) the net magnetic flux from the simulated quiet photosphere is not sufficient to maintain a chromospheric magnetic field (on average), 2) processes in the lower chromosphere, in the region dominated by magneto-acoustic shocks, are able to convert kinetic energy into magnetic energy, 3) the magnetic energy in the chromosphere increases linearly in time until the r.m.s. of the magnetic field strength saturates at roughly 4 to 30 G (horizontal average) due to conversion from kinetic energy, 4) and that the magnetic features formed in the chromosphere are localized to this region.

Context. We introduce a model for including accelerated particles in pure magnetohydrodynamics (MHD) simulations of the solar atmosphere. Aims. We show that the method is viable and produces results that enhance the realism of MHD simulations of the solar atmosphere. Methods. The acceleration of high-energy electrons in solar flares is an accepted fact, but is not included in the most advanced 3D simulations of the solar atmosphere. The effect of the acceleration is not known, and here we introduce a simple method to account for the ability of the accelerated electrons to move energy from the reconnection sites and into the dense transition zone and chromosphere. Results. The method was only run for a short time and with low reconnection energies, but this showed that the reconnection process itself changes, and that there is a clear effect on the observables at the impact sites of the accelerated electrons. Further work will investigate the effect on the reconnection sites and the impact sites in detail.

One candidate-model for heating the solar corona is magnetic reconnection that embodies Ohmic dissipation of current sheets. When numerous small-scaled magnetic reconnection occur, then it is possible to heat the corona. Due to the limitations of current instrumentation, nanoflares cannot be resolved. But their importance is evaluated via statistics by finding the power-law index of the energy distribution. This method is however biased due to technical and physical reasons. We aim to overcome limitations imposed by observations and statistical analysis. This way, we will identify, and study the small scale impulsive events. We employ a three-dimensional magnetohydrodynamic (3D-MHD) simulation using the \bifrost code. We also employ a new technique to identify the evolution of 3D Joule heating events in the corona. Then, we derive parameters describing the heating events in these locations. We report on the identification of heating events. We obtain the distribution of duration, released energy, and volume. We also find weak power-law correlation between these parameters. In addition, we extract information about geometrical parameters of 2D slices of 3D events, and about the evolution of resolved Joule heating compared to the total Joule heating and the magnetic energy in the corona. Even though the energy power index is less than 2, when classifying the energy release into three categories regarding with respect to the energy release (pico-, nano-, and micro-events), we find that nano-events release $82 \ \%$ of the resolved energy. This fraction corresponds to an energy flux larger than the one needed to heat the corona. It seems that the most popular population is the one containing short-lived, with small spatial extend events.

The measured average velocities in solar and stellar spectral lines formed at transition region temperatures have been difficult to interpret. However, realistic three-dimensional radiation magnetohydrodynamics (3D rMHD) models of the solar atmosphere are able to reproduce the observed dominant line shifts and may thus hold the key to resolve these issues. Our new 3D rMHD simulations aim to shed light on how mass flows between the chromosphere and corona and on how the coronal mass is maintained. Passive tracer particles, so-called corks, allow the tracking of parcels of plasma over time and thus the study of changes in plasma temperature and velocity not only locally, but also in a co-moving frame. By following the trajectories of the corks, we can investigate mass and energy flows and understand the composition of the observed velocities. Our findings show that most of the transition region mass is cooling. The preponderance of transition region redshifts in the model can be explained by the higher percentage of downflowing mass in the lower and middle transition region. The average upflows in the upper transition region can be explained by a combination of both stronger upflows than downflows and a higher percentage of upflowing mass. The most common combination at lower and middle transition region temperatures are corks that are cooling and traveling downward. For these corks, a strong correlation between the pressure gradient along the magnetic field line and the velocity along the magnetic field line has been observed, indicating a formation mechanism that is related to downward propagating pressure disturbances. Corks at upper transition region temperatures are subject to a rather slow and highly variable but continuous heating process.

We investigate the effects of interactions between ions and neutrals on the chromosphere and overlying corona using 2.5D radiative MHD simulations with the Bifrost code. We have extended the code capabilities implementing ion-neutral interaction effects using the Generalized Ohm's Law, i.e., we include the Hall term and the ambipolar diffusion (Pedersen dissipation) in the induction equation. Our models span from the upper convection zone to the corona, with the photosphere, chromosphere and transition region partially ionized. Our simulations reveal that the interactions between ionized particles and neutral particles have important consequences for the magneto-thermodynamics of these modeled layers: 1) ambipolar diffusion increases the temperature in the chromosphere; 2) sporadically the horizontal magnetic field in the photosphere is diffused into the chromosphere due to the large ambipolar diffusion; 3) ambipolar diffusion concentrates electrical currents leading to more violent jets and reconnection processes, resulting in 3a) the formation of longer and faster spicules, 3b) heating of plasma during the spicule evolution, and 3c) decoupling of the plasma and magnetic field in spicules. Our results indicate that ambipolar diffusion is a critical ingredient for understanding the magneto-thermo-dynamic properties in the chromosphere and transition region. The numerical simulations have been made publicly available, similar to previous Bifrost simulations. This will allow the community to study realistic numerical simulations with a wider range of magnetic field configurations and physics modules than previously possible.

The solar coronal heating problem is an open question since 1939. One proposed model for the transport and release of mechanical energy generated in the sub-phorospheric layers and photosphere is the nanoflare model that incorporates Ohmic heating which releases a part of the energy stored in the magnetic field via magnetic reconnection. The problem with the verification of this model is that we cannot resolve observationally small scale events. Histograms of observable characteristics of flares, show powerlaw behavior, for both energy release rate, size and total energy. Depending on the powerlaw index of the energy release, nanoflares might be an important candidate for coronal heating; we seek to find that index. In this paper, we employ a numerical 3D-MHD simulation produced by the numerical code Bifrost, and a new technique to identify the 3D heating events at a specific instant. The quantity we explore is the Joule heating, which is explicitly correlated with the magnetic reconnection because depends on the curl of the magnetic field. We are able to identify 4136 events in a volume $24 \times 24 \times 9.5 \ \textrm{Mm}^3$ (i.e. $768 \times 786 \times 331$ grid cells) of a specific snapshot. We find a powerlaw slope of the released energy per second, and two powerlaw slopes of the identified volume. The identified energy events do not represent all the released energy, but of the identified events, the total energy of the largest events dominate the energy release. Most of the energy release happens in the lower corona, while heating drops with height. We find that with a specific identification method that large events can be resolved into smaller ones, but at the expense of the total identified energy releases. The energy release which cannot be identified as an event favours a low energy release mechanism.

Using radiation magnetohydrodynamic simulations of the solar atmospheric layers from the upper convection zone to the lower corona, we investigate the self-consistent excitation of slow magneto-acoustic body waves (slow modes) in a magnetic flux concentration. We find that the convective downdrafts in the close surroundings of a two-dimensional flux slab "pump" the plasma inside it in the downward direction. This action produces a downflow inside the flux slab, which encompasses ever higher layers, causing an upwardly propagating rarefaction wave. The slow mode, excited by the adiabatic compression of the downflow near the optical surface, travels along the magnetic field in the upward direction at the tube speed. It develops into a shock wave at chromospheric heights, where it dissipates, lifts the transition region, and produces an offspring in the form of a compressive wave that propagates further into the corona. In the wake of downflows and propagating shock waves, the atmosphere inside the flux slab in the chromosphere and higher tends to oscillate with a period of $\nu\approx 4$~mHz. We conclude that this process of "magnetic pumping" is a most plausible mechanism for the direct generation of longitudinal chromospheric and coronal compressive waves within magnetic flux concentrations, and it may provide an important heat source in the chromosphere. It may also be responsible for certain types of dynamic fibrils.

The properties of non-statistical equilibrium ionization of silicon and oxygen ions are analyzed in this work. We focus on four solar targets (quiet sun, coronal hole, plage, quiescent active region, AR, and flaring AR) as observed with the Interface Region Imaging Spectrograph (IRIS). IRIS is best suited for this work due to the high cadence (up to 0.5s), high spatial resolution (up to 0.32"), and high signal to noise ratios for O IV and Si IV. We find that the observed intensity ratio between lines of three times ionized silicon and oxygen ions depends on their total intensity and that this correlation varies depending on the region observed (quiet sun, coronal holes, plage or active regions) and on the specific observational objects present (spicules, dynamic loops, jets, micro-flares or umbra). In order to interpret the observations, we compare them with synthetic profiles taken from 2D self-consistent radiative MHD simulations of the solar atmosphere, where the statistical equilibrium or non-equilibrium treatment silicon and oxygen is applied. These synthetic observations show vaguely similar correlations as in the observations, i.e. between the intensity ratios and their intensities, but only in the non-equilibrium case do we find that (some of) the observations can be reproduced. We conclude that these lines are formed out of statistical equilibrium. We use our time-dependent non-equilibrium ionization simulations to describe the physical mechanisms behind these observed properties.

Context. The solar chromosphere is the interface between the solar surface and the solar corona. Modelling of this region is difficult because it represents the transition from optically thick to thin radiation escape, from gas-pressure domination to magnetic-pressure domination, from a neutral to an ionised state, from MHD to plasma physics, and from near-equilibrium (LTE) to non-equilibrium conditions. Aims. Our aim is to provide the community with realistic simulations of the magnetic solar outer atmosphere. This will enable detailed comparison of existing and upcoming observations with synthetic observables from the simulations, thereby elucidating the complex interactions of magnetic fields and plasma that are crucial for our understanding of the dynamic outer atmosphere. Methods. We used the radiation magnetohydrodynamics code Bifrost to perform simulations of a computational volume with a magnetic field topology similar to an enhanced network area on the Sun. Results. The full simulation cubes are made available online. The general properties of the simulation are discussed, and limitations are discussed.

The structure and dynamics of the outer solar atmosphere are reviewed with emphasis on the role played by the magnetic field. Contemporary observations that focus on high resolution imaging over a range of temperatures, as well as UV, EUV and hard X-ray spectroscopy, demonstrate the presence of a vast range of temporal and spatial scales, mass motions, and particle energies present. By focussing on recent developments in the chromosphere, corona and solar wind, it is shown that small scale processes, in particular magnetic reconnection, play a central role in determining the large-scale structure and properties of all regions. This coupling of scales is central to understanding the atmosphere, yet poses formidable challenges for theoretical models.

Context: Numerical simulations of stellar convection and photospheres have been developed to the point where detailed shapes of observed spectral lines can be explained. Stellar atmospheres are very complex, and very different physical regimes are present in the convection zone, photosphere, chromosphere, transition region and corona. To understand the details of the atmosphere it is necessary to simulate the whole atmosphere since the different layers interact strongly. These physical regimes are very diverse and it takes a highly efficient massively parallel numerical code to solve the associated equations. Aims: The design, implementation and validation of the massively parallel numerical code Bifrost for simulating stellar atmospheres from the convection zone to the corona. Methods: The code is subjected to a number of validation tests, among them the Sod shock tube test, the Orzag-Tang colliding shock test, boundary condition tests and tests of how the code treats magnetic field advection, chromospheric radiation, radiative transfer in an isothermal scattering atmosphere, hydrogen ionization and thermal conduction. Results: Bifrost completes the tests with good results and shows near linear efficiency scaling to thousands of computing cores.

aims: We aim to provide an estimate of the minimum temperature of the quiet solar chromosphere. methods: We perform a 2D radiation-MHD simulation spanning the upper convection zone to the lower corona. The simulation includes non-LTE radiative transfer and {an equation-of-state that includes non-equilibrium ionization of hydrogen and non-equilibrium H_2 molecule formation}. We analyze the reliability of the various assumptions made in our model in order to assess the realism of the simulation. results: Our simulation contains pockets of cool gas with down to 1660 K from 1 Mm up to 3.2 Mm height. It overestimates the radiative heating, and contains non-physical heating below 1660 K. Therefore we conclude that cool pockets in the quiet solar chromosphere might have even lower temperatures than in the simulation, provided that there exist areas in the chromosphere without significant magnetic heating. We suggest off-limb molecular spectroscopy to look for such cool pockets and 3D simulations including a local dynamo and a magnetic carpet to investigate Joule heating in the quiet chromosphere.

Three-dimensional (3D) radiative hydrodynamic model atmospheres of metal-poor late-type stars are characterized by cooler upper photospheric layers than their 1D counterparts. This property of 3D models can dramatically affect elemental abundances derived from temperature-sensitive spectral lines. We investigate whether the cool surface temperatures predicted by metal-poor 3D models can be ascribed to the approximated treatment of scattering in the radiative transfer. We use the Bifrost code to test three different ways to handle scattering in 3D model atmospheres of metal-poor stars. First, we solve self-consistently the radiative transfer equation for a source function with a coherent scattering term. Second, we solve the radiative transfer equation for a Planckian source function, neglecting the contribution of continuum scattering to extinction in the optically thin layers; this has been the default mode in previous models of ours. Third, we treat scattering as pure absorption everywhere, which is the standard case in CO5BOLD models. We find that the second approach produces temperature structures with cool upper photospheric layers very similar to the correct coherent scattering solution. In contrast, treating scattering as pure absorption leads to significantly hotter and shallower temperature stratifications. The main differences in temperature structure between our published models and those generated with the CO5BOLD code can be traced to the different treatments of scattering. Neglecting the contribution of continuum scattering to extinction in optically thin layers provides a good approximation to the full radiative transfer solution for metal-poor stars. Our results demonstrate that the cool temperature stratifications predicted for metal-poor late-type stellar atmospheres by previous models of ours are not an artifact of the approximated treatment of scattering.

We present the implementation of a radiative transfer solver with coherent scattering in the new BIFROST code for radiative magneto-hydrodynamical (MHD) simulations of stellar surface convection. The code is fully parallelized using MPI domain decomposition, which allows for large grid sizes and improved resolution of hydrodynamical structures. We apply the code to simulate the surface granulation in a solar-type star, ignoring magnetic fields, and investigate the importance of coherent scattering for the atmospheric structure. A scattering term is added to the radiative transfer equation, requiring an iterative computation of the radiation field. We use a short-characteristics-based Gauss-Seidel acceleration scheme to compute radiative flux divergences for the energy equation. The effects of coherent scattering are tested by comparing the temperature stratification of three 3D time-dependent hydrodynamical atmosphere models of a solar-type star: without scattering, with continuum scattering only, and with both continuum and line scattering. We show that continuum scattering does not have a significant impact on the photospheric temperature structure for a star like the Sun. Including scattering in line-blanketing, however, leads to a decrease of temperatures by about 350\,K below log tau < -4. The effect is opposite to that of 1D hydrostatic models in radiative equilibrium, where scattering reduces the cooling effect of strong LTE lines in the higher layers of the photosphere. Coherent line scattering also changes the temperature distribution in the high atmosphere, where we observe stronger fluctuations compared to a treatment of lines as true absorbers.

We have performed 3D radiation MHD simulations extending from the convection zone to the corona covering a box 16 Mm$^3$ at 32 km spatial resolution. The simulations show very fine structure in the chromosphere with acoustic shocks interacting with the magnetic field. Magnetic flux concentrations have a temperature lower than the surroundings in the photosphere but higher in the low chromosphere. The heating is there mostly through ohmic dissipation preferentially at the edges of the flux concentrations. The magnetic field is often wound up around the flux concentrations. When acoustic waves travel up along the field this topology leads to swirling motions seen in chromospheric diagnostic lines such as the calcium infrared triplet.

We perform a number of numerical simulations of the solar corona with the aim to understand how it responds to different conditions in the photosphere. By changing parameters which govern the motion of the plasma at the photosphere we study the behavior of the corona, in particular, the effects on the current density generated. An MHD code is used to run simulations, using a 20x20x20 Mm^3 box with time spans ranging from one hundred to several hundreds of minutes. All the experiments show a fast initial increase of the current density, followed by a stabilization around an asymptotic value which depends on the photospheric conditions. These asymptotic average current densities as well as the turn-over points are discussed.

A major goal in solar physics has during the last five decades been to find how energy flux generated in the solar convection zone is transported and dissipated in the outer solar layers. Progress in this field has been slow and painstaking. However, advances in computer hardware and numerical methods, vastly increased observational capabilities and growing physical insight seem finally to be leading towards understanding. Here we present exploratory numerical MHD models that span the entire solar atmosphere from the upper convection zone to the lower corona. These models include non-grey, non-LTE radiative transport in the photosphere and chromosphere, optically thin radiative losses as well as magnetic field-aligned heat conduction in the transition region and corona.

A forward model is described in which we synthesize spectra from an ab-initio 3D MHD simulation of an outer stellar atmosphere, where the coronal heating is based on braiding of magnetic flux due to photospheric footpoint motions. We discuss the validity of assumptions such as ionization equilibrium and investigate the applicability of diagnostics like the differential emission measure inversion. We find that the general appearance of the synthesized corona is similar to the solar corona and that, on a statistical basis, integral quantities such as average Doppler shifts or differential emission measures are reproduced remarkably well. The persistent redshifts in the transition region, which have puzzled theorists since their discovery, are explained by this model as caused by the flows induced by the heating through braiding of magnetic flux. While the model corona is only slowly evolving in intensity, as is observed, the amount of structure and variability in Doppler shift is very large. This emphasizes the need for fast coronal spectroscopy, as the dynamical response of the corona to the heating process manifests itself in a comparably slow evolving coronal intensity but rapid changes in Doppler shift.

Cool stars like our Sun are surrounded by a million degree hot outer atmosphere, the corona. Since more than 60 years the physical nature of the processes heating the corona to temperatures well in excess of those on the stellar surface remain puzzling. Recent progress in observational techniques and numerical modeling now opens a new window to approach this problem. We present the first coronal emission line spectra synthesized from three-dimensional numerical models describing the evolution of the dynamics and energetics as well as of the magnetic field in the corona. In these models the corona is heated through motions on the stellar surface that lead to a braiding of magnetic field lines inducing currents which are finally dissipated. These forward models enable us to synthesize observed properties like (average) emission line Doppler shifts or emission measures in the outer atmosphere, which until now have not been understood theoretically, even though many suggestions have been made in the past. As our model passes these observational tests, we conclude that the flux braiding mechanism is a prime candidate for being the dominant heating process of the magnetically closed corona of the Sun and solar-like stars.

Magnetic activity in M-dwarfs present enigmatic questions: On the one hand they have higher field strengths and larger filling factors than the magnetic field on the Sun, on the other hand, they are fully convective and their atmospheres are more neutral, hence they do not have an undershoot layer for magnetic flux storage and as we show here, cannot have small-scale dynamo action in their photospheres either. We present a discussion of these facts and propose a new numerical model to investigate M-dwarf magnetism.

Data from recent numerical simulations of the solar corona and transition region are analysed and the magnetic field connection between the low corona and the photosphere is found to be close to that of a potential field. The fieldline to fieldline displacements follow a power law distribution with typical displacements of just a few Mm. Three loops visible in emulated Transition Region And Coronal Explorer (TRACE) filters are analysed in detail and found to have significantly different heating rates and distributions thereof, one of them showing a small scale heating event. The dynamical structure is complicated even though all the loops are visible in a single filter along most of their lengths. None of the loops are static, but are in the process of evolving into loops with very different characteristics. Differential Emission Measure (DEM) curves along one of the loops illustrate that DEM curves have to be treated carefully if physical characteristics are to be extracted.

We present an ab initio approach to the solar coronal heating problem by modelling a small part of the solar corona in a computational box using a 3D MHD code including realistic physics. The observed solar granular velocity pattern and its amplitude and vorticity power spectra, as reproduced by a weighted Voronoi tessellation method, are used as a boundary condition that generates a Poynting flux in the presence of a magnetic field. The initial magnetic field is a potential extrapolation of a SOHO/MDI high resolution magnetogram, and a standard stratified atmosphere is used as a thermal initial condition. Except for the chromospheric temperature structure, which is kept fixed, the initial conditions are quickly forgotten because the included Spitzer conductivity and radiative cooling function have typical timescales much shorter than the time span of the simulation. After a short initial start up period, the magnetic field is able to dissipate 3-4 10^6 ergs cm^{-2} s^{-1} in a highly intermittent corona, maintaining an average temperature of $\sim 10^6$ K, at coronal density values for which emulated images of the Transition Region And Coronal Explorer(TRACE) 171 and 195 pass bands reproduce observed photon count rates.

Numerical simulations of the multi-phase interstellar medium have been carried out, using a 3D, nonlinear, magnetohydrodynamic, shearing-box model, with random motions driven by supernova explosions. These calculations incorporate the effects of magnetic fields and rotation in 3D; these play important dynamical roles in the galaxy, but are neglected in many other simulations. The supernovae driving the motions are not arbitrarily imposed, but occur where gas accumulates into cold, dense clouds; their implementation uses a physically motivated model for the evolution of such clouds. The process is self-regulating, and produces mean supernova rates as part of the solution. Simulations with differing mean density show a power law relation between the supernova rate and density, with exponent 1.7; this value is within the range suggested from observations (taking star formation rate as a proxy for supernova rate). The global structure of the supernova driven medium is strongly affected by the presence of magnetic fields; e.g. for one solution the filling factor of hot gas is found to vary from 0.19 (with no field) to 0.12 (with initial mid-plane field B_0=6 muG).

Statistical parameters of the ISM driven by thermal energy injections from supernova explosions have been obtained from 3D, nonlinear, magnetohydrodynamic, shearing-box simulations for spiral arm and interarm regions. The density scale height obtained for the interarm regions is 50% larger than within the spiral arms because of the higher gas temperature. The filling factor of the hot gas is also significantly larger between the arms and depends sensitively on magnetic field strength.

The heating of the solar corona and the puzzle of the slender high reaching magnetic loops seen in observations from the Transition Region And Coronal Explorer(TRACE) has been investigated through 3D numerical simulations, and found to be caused by the well observed plasma flows in the photosphere displacing the footpoints of magnetic loops in a nearly potential configuration. It is found that even the small convective displacements cause magnetic dissipation sufficient to heat the corona to temperatures of the order of a million Kelvin. The heating is intermittent in both space and time - at any one height and time it spans several orders of magnitude, and localized heating causes transonic flows along field lines, which explains the observed non-hydrostatic equilibrium of loops that are bright in emission measure.

We present 3-d numerical magneto-hydrodynamic simulations of a buoyant, twisted magnetic flux rope embedded in a stratified, solar-like model convection zone. The flux rope is given an initial twist such that it neither kinks nor fragments during its ascent. Moreover, its magnetic energy content with respect to convection is chosen so that the flux rope retains its basic geometry while being deflected from a purely vertical ascent by convective flows. The simulations show that magnetic flux is advected away from the core of the flux rope as it interacts with the convection. The results thus support the idea that the amount of toroidal flux stored at or near the bottom of the solar convection zone may currently be underestimated.