There has been tremendous progress in the degree of realism of three-dimensional radiation magneto-hydrodynamic simulations of the solar atmosphere in the past decades. Four of the most frequently used numerical codes are Bifrost, CO5BOLD, MANCHA3D, and MURaM. Here we test and compare the wave propagation characteristics in model runs from these four codes by measuring the dispersion relation of acoustic-gravity waves at various heights. We find considerable differences between the various models. The height dependence of wave power, in particular of high-frequency waves, varies by up to two orders of magnitude between the models, and the phase difference spectra of several models show unexpected features, including $\pm180^\circ$ phase jumps.

We study vortex dynamics in the solar atmosphere by employing and deriving the analytical evolution equations of two vortex identification criteria. The two criteria used are vorticity and the swirling strength. Vorticity can be biased in the presence of shear flows, but its dynamical equation is well known; the swirling strength is a more precise criterion for the identification of vortical flows, but its evolution equation is not known yet. Therefore, we explore the possibility of deriving a dynamical equation for the swirling strength. We then apply the two equations to analyze radiative MHD simulations of the solar atmosphere produced with the CO5BOLD code. We present a detailed review of the swirling strength criterion and the mathematical derivation of its evolution equation. This equation did not exist in the literature before and it constitutes a novel tool that is suitable for the analysis of a wide range of problems in (magneto-)hydrodynamics. By applying this equation to numerical models, we find that hydrodynamical and magnetic baroclinicities are the driving physical processes responsible for vortex generation in the convection zone and the photosphere. Higher up in the chromosphere, the magnetic terms alone dominate. Moreover, we find that the swirling strength is produced at small scales in a chaotic fashion, especially inside magnetic flux concentrations. The swirling strength represents an appropriate criterion for the identification of vortices in turbulent flows, such as those in the solar atmosphere. Moreover, its evolution equation, which is derived in this paper, is pivotal for obtaining precise information about the dynamics of these vortices and the physical mechanisms responsible for their production and evolution. Since this equation is available, the swirling strength is now the ideal quantity to study the dynamics of vortices in (magneto-)hydrodynamics.

In this second paper of the series on internal gravity waves (IGWs), we present a study of the generation and propagation of IGWs in a model solar atmosphere with diverse magnetic conditions. A magnetic field free, and three magnetic models that start with an initial, vertical, homogeneous field of 10 G, 50 G, and 100 G magnetic flux density, are simulated using the CO$^{\rm 5}$BOLD code. We find that the IGWs are generated in similar manner in all four models in spite of the differences in the magnetic environment. The mechanical energy carried by IGWs is significantly larger than that of the acoustic waves in the lower part of the atmosphere, making them an important component of the total wave energy budget. The mechanical energy flux (10$^{6}$ - 10$^{3}$ W m$^{-2}$) is few orders of magnitude larger than the Poynting flux (10$^{3}$ - 10$^{1}$ W m$^{-2}$). The Poynting fluxes show a downward component in the frequency range corresponding to the IGWs, which confirm that these waves do not propagate upwards in the atmosphere when the fields are predominantly vertical and strong. We conclude that, in the upper photosphere, the propagation properties of IGWs depend on the average magnetic field strength and therefore these waves can be potential candidate for magnetic field diagnostics of these layers. However, their subsequent coupling to Alfv\'{e}nic waves are unlikely in a magnetic environment permeated with predominantly vertical fields and therefore they may not directly or indirectly contribute to the heating of layers above plasma-$\beta$ less than 1.

The numerical computation of reliable and accurate Stokes profiles is of great relevance in solar physics. In the synthesis process, many actors play a relevant role: among them the formal solver, the discrete atmospheric model, and the spectral line. This paper tests the performances of different numerical schemes in the synthesis of polarized spectra for different spectral lines and atmospheric models. The hierarchy between formal solvers is enforced, stressing the peculiarities of high-order and low-order formal solvers. The density of grid points necessary for reaching a given accuracy requirement is quantitatively described for specific situations.

The discussion regarding the numerical integration of the polarized radiative transfer equation is still open and the comparison between the different numerical schemes proposed by different authors in the past is not fully clear. Aiming at facilitating the comprehension of the advantages and drawbacks of the different formal solvers, this work presents a reference paradigm for their characterization based on the concepts of order of accuracy, stability, and computational cost. Special attention is paid to understand the numerical methods belonging to the Diagonal Element Lambda Operator family, in an attempt to highlight their specificities.

When integrating the radiative transfer equation for polarized light, the necessity of high-order numerical methods is well known. In fact, well-performing high-order formal solvers enable higher accuracy and the use of coarser spatial grids. Aiming to provide a clear comparison between formal solvers, this work presents different high-order numerical schemes and applies the systematic analysis proposed by Janett et al. (2017), emphasizing their advantages and drawbacks in terms of order of accuracy, stability, and computational cost.

We characterize waves in small magnetic elements and investigate their propagation in the lower solar atmosphere from observations at high spatial and temporal resolution. We use the wavelet transform to analyze oscillations of both horizontal displacement and intensity in magnetic bright points found in the 300 nm and the Ca II H 396.8 nm passbands of the filter imager on board the Sunrise balloon-borne solar observatory. Phase differences between the oscillations at the two atmospheric layers corresponding to the two passbands reveal upward propagating waves at high frequencies (up to 30 mHz). Weak signatures of standing as well as downward propagating waves are also obtained. Both compressible and incompressible (kink) waves are found in the small-scale magnetic features. The two types of waves have different, though overlapping, period distributions. Two independent estimates give a height difference of approximately 450+-100 km between the two atmospheric layers sampled by the employed spectral bands. This value, together with the determined short travel times of the transverse and longitudinal waves provide us with phase speeds of 29+-2 km/s and 31+-2 km/s, respectively. We speculate that these phase speeds may not reflect the true propagation speeds of the waves. Thus, effects such as the refraction of fast longitudinal waves may contribute to an overestimate of the phase speed.

Observations of the solar atmosphere show that internal gravity waves are generated by overshooting convection, but are suppressed at locations of magnetic flux, which is thought to be the result of mode conversion into magneto-acoustic waves. Here, we present a study of the acoustic-gravity wave spectrum emerging from a realistic, self-consistent simulation of solar (magneto-)convection. A magnetic field free, hydrodynamic simulation and a magneto-hydrodynamic (MHD) simulation with an initial, vertical, homogeneous field of 50 G flux density were carried out and compared with each other to highlight the effect of magnetic fields on the internal gravity wave propagation in the Sun's atmosphere. We find that the internal gravity waves are absent or partially reflected back into the lower layers in the presence of magnetic fields and argue that the suppression is due to the coupling of internal gravity waves to slow magneto acoustic waves still within the high-beta region of the upper photosphere. The conversion to Alfven waves is highly unlikely in our model because there is no strongly inclined magnetic field present. We argue that the suppression of internal waves observed within magnetic flux-concentrations may also be due to non-linear breaking of internal waves due to vortex flows that are ubiquitously present in the upper photosphere and the chromosphere.

Context. Small-scale bright features in the photosphere of the Sun, such as faculae or G-band bright points, appear in connection with small-scale magnetic flux concentrations. Aims. Here we report on a new class of photospheric bright points that are free of magnetic fields. So far, these are visible in numerical simulations only. We explore conditions required for their observational detection. Methods. Numerical radiation (magneto-)hydrodynamic simulations of the near-surface layers of the Sun were carried out. The magnetic field-free simulations show tiny bright points, reminiscent of magnetic bright points, only smaller. A simple toy model for these non-magnetic bright points (nMBPs) was established that serves as a base for the development of an algorithm for their automatic detection. Basic physical properties of 357 detected nMBPs were extracted and statistically evaluated. We produced synthetic intensity maps that mimic observations with various solar telescopes to obtain hints on their detectability. Results. The nMBPs of the simulations show a mean bolometric intensity contrast with respect to their intergranular surroundings of approximately 20%, a size of 60-80 km, and the isosurface of optical depth unity is at their location depressed by 80-100 km. They are caused by swirling downdrafts that provide, by means of the centripetal force, the necessary pressure gradient for the formation of a funnel of reduced mass density that reaches from the subsurface layers into the photosphere. Similar, frequently occurring funnels that do not reach into the photosphere, do not produce bright points. Conclusions. Non-magnetic bright points are the observable manifestation of vertically extending vortices (vortex tubes) in the photosphere. The resolving power of 4-m-class telescopes, such as the DKIST, is needed for an unambiguous detection of them.

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.

We present the first radiation magnetohydrodynamics simulations of the atmosphere of white dwarf stars. We demonstrate that convective energy transfer is seriously impeded by magnetic fields when the plasma-beta parameter, the thermal to magnetic pressure ratio, becomes smaller than unity. The critical field strength that inhibits convection in the photosphere of white dwarfs is in the range B = 1-50 kG, which is much smaller than the typical 1-1000 MG field strengths observed in magnetic white dwarfs, implying that these objects have radiative atmospheres. We have then employed evolutionary models to study the cooling process of high-field magnetic white dwarfs, where convection is entirely suppressed during the full evolution (B > 10 MG). We find that the inhibition of convection has no effect on cooling rates until the effective temperature (Teff) reaches a value of around 5500 K. In this regime, the standard convective sequences start to deviate from the ones without convection owing to the convective coupling between the outer layers and the degenerate reservoir of thermal energy. Since no magnetic white dwarfs are currently known at the low temperatures where this coupling significantly changes the evolution, effects of magnetism on cooling rates are not expected to be observed. This result contrasts with a recent suggestion that magnetic white dwarfs with Teff < 10,000 K cool significantly slower than non-magnetic degenerates.

High-resolution observations with the Swedish 1-m Solar Telescope (SST) and the Solar Dynamics Observatory (SDO) reveal rotating magnetic field structures that extend from the solar surface into the chromosphere and the corona. These so-called magnetic tornadoes are primarily detected as rings or spirals of rotating plasma in the Ca II 854.2 nm line core (also known as chromospheric swirls). Detailed numerical simulations show that the observed chromospheric plasma motion is caused by the rotation of magnetic field structures, which again are driven by photospheric vortex flows at their footpoints. Under the right conditions, two vortex flow systems are stacked on top of each other. We refer to the lower vortex, which extends from the low photosphere into the convection zone, as intergranular vortex flow (IVF). Once a magnetic field structure is co-located with an IVF, the rotation is mediated into the upper atmospheric layers and an atmospheric vortex flow (AVF, or magnetic tornado) is generated. In contrast to the recent work by Shelyag et al., we demonstrate that particle trajectories in a simulated magnetic tornado indeed follow spirals and argue that the properties of the trajectories decisively depend on the location in the atmosphere and the strength of the magnetic field.

Chromospheric swirls are the observational signatures of rotating magnetic field structures in the solar atmosphere, also known as magnetic tornadoes. Swirls appear as dark rotating features in the core of the spectral line of singly ionized calcium at a wavelength of 854.2 nm. This signature can be very subtle and difficult to detect given the dynamic changes in the solar chromosphere. Important steps towards a systematic and objective detection method are the compilation and characterization of a statistically significant sample of observed and simulated chromospheric swirls. Here, we provide a more exact definition of the chromospheric swirl phenomenon and also present a first morphological classification of swirls with three types: (I) Ring, (II) Split, (III) Spiral. We also discuss the nature of the magnetic field structures connected to tornadoes and the influence of limited spatial resolution on the appearance of their photospheric footpoints.

We present first results from three-dimensional radiation magnetohydrodynamic simulations of M-type dwarf stars with CO5BOLD. The local models include the top of the convection zone, the photosphere, and the chromosphere. The results are illustrated for models with an effective temperature of 3240 K and a gravitational acceleration of log g = 4.5, which represent analogues of AD Leo. The models have different initial magnetic field strengths and field topologies. This first generation of models demonstrates that the atmospheres of M-dwarfs are highly dynamic and intermittent. Magnetic fields and propagating shock waves produce a complicated fine-structure, which is clearly visible in synthetic intensity maps in the core of the Ca II K spectral line and also at millimeter wavelengths. The dynamic small-scale pattern cannot be described by means of one-dimensional models, which has important implications for the construction of semi-empirical model atmospheres and thus for the interpretation of observations in general. Detailed three-dimensional numerical simulations are valuable in this respect. Furthermore, such models facilitate the analysis of small-scale processes, which cannot be observed on stars but nevertheless might be essential for understanding M-dwarf atmospheres and their activity. An example are so-called "magnetic tornadoes", which have recently been found on the Sun and are presented here in M-dwarf models for the first time.

We investigate the interaction of magneto-acoustic waves with magnetic network elements with the aim of finding possible signatures of the magnetic shadow phenomenon in the vicinity of network elements. We carried out three-dimensional numerical simulations of magneto-acoustic wave propagation in a model solar atmosphere that is threaded by a complexly structured magnetic field, resembling that of a typical magnetic network element and of internetwork regions. High-frequency waves of 10 mHz are excited at the bottom of the simulation domain. On their way through the upper convection zone and through the photosphere and the chromosphere they become perturbed, refracted, and converted into different mode types. We applied a standard Fourier analysis to produce oscillatory power-maps of the line-of-sight velocity. In the power maps of the upper photosphere and the lower chromosphere, we clearly see the magnetic shadow: a seam of suppressed power surrounding the magnetic network elements. We demonstrate that this shadow is linked to the mode conversion process and that power maps at these height levels show the signature of three different magneto-acoustic wave modes.

We review (i) observations and numerical simulations of vortical flows in the solar atmosphere and (ii) measurements of the horizontal magnetic field in quiet Sun regions. First, we discuss various manifestations of vortical flows and emphasize the role of magnetic fields in mediating swirling motion created near the solar surface to the higher layers of the photosphere and to the chromosphere. We reexamine existing simulation runs of solar surface magnetoconvection with regard to vortical flows and compare to previously obtained results. Second, we review contradictory results and problems associated with measuring the angular distribution of the magnetic field in quiet Sun regions. Furthermore, we review the Stokes-V-amplitude ratio method for the lines Fe I 630.15 and 630.25 nm. We come to the conclusion that the recently discovered two distinct populations in scatter plots of this ratio must not bee interpreted in terms of 'uncollapsed' and 'collapsed' fields but they stem from weak granular magnetic fields and weak canopy fields located at the boundaries between granules and the intergranular space. Based on new simulation runs, we reaffirm earlier findings of a predominance of the horizontal field components over the vertical one, particularly in the upper photosphere and at the base of the chromosphere.

The solar atmosphere is magnetically structured and highly dynamic. Owing to the dynamic nature of the regions in which the magnetic structures exist, waves can be excited in them. Numerical investigations of wave propagation in small-scale magnetic flux concentrations in the magnetic network on the Sun have shown that the nature of the excited modes depends on the value of plasma beta (the ratio of gas to magnetic pressure) where the driving motion occurs. Considering that these waves should give rise to observable characteristic signatures, we have attempted a study of synthesized emergent spectra from numerical simulations of magneto-acoustic wave propagation. We find that the signatures of wave propagation in a magnetic element can be detected when the spatial resolution is sufficiently high to clearly resolve it, enabling observations in different regions within the flux concentration. The possibility to probe various lines of sight around the flux concentration bears the potential to reveal different modes of the magnetohydrodynamic waves and mode conversion. We highlight the feasibility of using the Stokes-V asymmetries as a diagnostic tool to study the wave propagation within magnetic flux concentrations. These quantities can possibly be compared with existing and new observations in order to place constraints on different wave excitation mechanisms.

High-resolution images of the solar surface show a granulation pattern of hot rising and cooler downward-sinking material -- the top of the deep-reaching solar convection zone. Convection plays a role for the thermal structure of the solar interior and the dynamo acting there, for the stratification of the photosphere, where most of the visible light is emitted, as well as for the energy budget of the spectacular processes in the chromosphere and corona. Convective stellar atmospheres can be modeled by numerically solving the coupled equations of (magneto)hydrodynamics and non-local radiation transport in the presence of a gravity field. The CO5BOLD code described in this article is designed for so-called "realistic" simulations that take into account the detailed microphysics under the conditions in solar or stellar surface layers (equation-of-state and optical properties of the matter). These simulations indeed deserve the label "realistic" because they reproduce the various observables very well -- with only minor differences between different implementations. The agreement with observations has improved over time and the simulations are now well-established and have been performed for a number of stars. Still, severe challenges are encountered when it comes to extending these simulations to include ideally the entire star or substellar object: the strong stratification leads to completely different conditions in the interior, the photosphere, and the corona. Simulations have to account for a large range of spatial and time scales and also non-equilibrium processes. Last but not least, realistic simulations are based on detailed microphysics and depend on the quality of the input data. This article provides an overview of the physical problem and the numerical solution and the capabilities of CO5BOLD, illustrated with a number of applications.

We observe that intergranular jets, originating in the intergranular space surrounding individual granules, tend to be associated with granular fragmentation, in particular, with the formation and evolution of a bright granular lane (BGL) within individual granules. The BGLs have recently been identified as vortex tubes by Steiner et al. We further discover the development of a well-defined bright grain located between the BGL and the dark intergranular lane to which it is connected. Signatures of a BGL may reach the lower chromosphere and can be detected in off-band \ha images. Simulations also indicate that vortex tubes are frequently associated with small-scale magnetic fields. We speculate that the intergranular jets detected in the NST data may result from the interaction between the turbulent small-scale fields associated with the vortex tube and the larger-scale fields existing in the intergranular lanes. The intergranular jets are much smaller and weaker than all previously known jet-like events. At the same time, they appear much more numerous than the larger events, leading us to the speculation that the total energy release and mass transport by these tiny events may not be negligible in the energy and mass-flux balance near the temperature minimum atop the photosphere. The study is based on the photospheric TiO broadband (1.0 nm) filter data acquired with the 1.6 m New Solar Telescope (NST) operating at the Big Bear Solar Observatory. The data set also includes NST off-band \ha images collected through a Zeiss Lyot filter with a passband of 0.025 nm.

From radiation magnetohydrodynamic simulations of the solar atmosphere we find a new mechanism for the excitation of longitudinal slow modes within magnetic flux concentrations. We find that the convective downdrafts in the immediate surroundings of magnetic elements are responsible for the excitation of slow modes. The coupling between the external downdraft and the plasma motion internal to the flux concentration is mediated by the inertial forces of the downdraft that act on the magnetic flux concentration. These forces, in conjunction with the downward movement, pump the internal atmosphere in the downward direction, which entails a fast downdraft in the photospheric and chromospheric layers of the magnetic element. Subsequent to the transient pumping phase, the atmosphere rebounds, causing a slow mode traveling along the magnetic flux concentration in the upward direction. It develops into a shock wave in chromospheric heights, possibly capable of producing some kind of dynamic fibril. We propose an observational detection of this process.

We present two-dimensional simulations of wave propagation in a realistic, non-stationary model of the solar atmosphere. This model shows a granular velocity field and magnetic flux concentrations in the intergranular lanes similar to observed velocity and magnetic structures on the Sun and takes radiative transfer into account. We present three cases of magneto-acoustic wave propagation through the model atmosphere, where we focus on the interaction of different magneto-acoustic wave at the layer of similar sound and Alfv\'en speeds, which we call the equipartition layer. At this layer the acoustic and magnetic mode can exchange energy depending on the angle between the wave vector and the magnetic field vector. Our results show that above the equipartition layer and in all three cases the fast magnetic mode is refracted back into the solar atmosphere. Thus, the magnetic wave shows an evanescent behavior in the chromosphere. The acoustic mode, which travels along the magnetic field in the low plasma-$\beta$ regime, can be a direct consequence of an acoustic source within or outside the low-$\beta$ regime, or it can result from conversion of the magnetic mode, possibly from several such conversions when the wave travels across a series of equipartition layers.

We present simulations of magneto-acoustic wave propagation in a magnetic, plane-parallel stratified solar model atmosphere, employing the CO5BOLD-code. The tests are carried out for two models of the solar atmosphere, which are similar to the ones used by Cally (2007) and Schunker & Cally (2006). The two models differ only in the orientation of the magnetic field. A qualitative comparison shows good agreement between the numerical results and the results from ray theory. The tests are done in view of the application of the present numerical code for the computation of energy fluxes of propagating acoustic waves into a dynamically evolving magnetic solar atmosphere. For this, we consider waves with frequencies above the acoustic cut-off frequency.

We have investigated a time series of continuum intensity maps and corresponding Dopplergrams of granulation in a very quiet solar region at the disk center, recorded with the Imaging Magnetograph eXperiment (IMaX) on board the balloon-borne solar observatory Sunrise. We find that granules frequently show substructure in the form of lanes composed of a leading bright rim and a trailing dark edge, which move together from the boundary of a granule into the granule itself. We find strikingly similar events in synthesized intensity maps from an ab initio numerical simulation of solar surface convection. From cross sections through the computational domain of the simulation, we conclude that these `granular lanes' are the visible signature of (horizontally oriented) vortex tubes. The characteristic optical appearance of vortex tubes at the solar surface is explained. We propose that the observed vortex tubes may represent only the large-scale end of a hierarchy of vortex tubes existing near the solar surface.

We examine spectra of the Ca II H line, obtained under good seeing conditions with the VTT Echelle Spectrograph in June of 2007, and higher resolution data of the Ca II 8542 Angstrom line from Fabry-Perot instruments. The VTT targets were areas near disk center which included quiet Sun and some dispersed plage. The infrared data included quiet Sun and plage associated with small pores. Bright chromospheric network emission patches expand little with wavelength from line wing to line center, i.e. with increasing line opacity and height. We argue that this simple observation has implications for the force and energy balance of the chromosphere, since bright chromospheric network emission is traditionally associated with enhanced local mechanical heating which increases temperatures and pressures. Simple physical considerations then suggest that the network chromosphere may not be able to reach horizontal force balance with its surroundings, yet the network is a long-lived structure. We speculate on possible reasons for the observed behavior. By drawing attention to a potential conundrum, we hope to contribute to a better understanding of a long-standing unsolved problem: the heating of the chromospheric network.

We investigate wave propagation and energy transport in magnetic elements, which are representatives of small scale magnetic flux concentrations in the magnetic network on the Sun. This is a continuation of earlier work by Hasan et al. (2005). The new features in the present investigation include a quantitative evaluation of the energy transport in the various modes and for different field strengths, as well as the effect of the boundary-layer thickness on wave propagation. We carry out 2-D MHD numerical simulations of magnetic flux concentrations for strong and moderate magnetic fields. Waves are excited in the tube and ambient medium by a transverse impulsive motion of the lower boundary. The nature of the modes excited depends on the value of beta. Mode conversion occurs in the moderate field case when the fast mode crosses the beta=1 contour. In the strong field case the fast mode undergoes conversion from predominantly magnetic to predominantly acoustic when waves are leaking from the interior of the flux concentration to the ambient medium. We also estimate the energy fluxes in the acoustic and magnetic modes. The main conclusions of our work are twofold: firstly, for transverse, impulsive excitation, flux tubes/sheets with strong fields are more efficient than those with weak fields in providing acoustic flux to the chromosphere. However, there is insufficient energy in the acoustic flux to balance the chromospheric radiative losses in the network, even for the strong field case. Secondly, the acoustic emission from the interface between the flux concentration and the ambient medium decreases with the width of the boundary layer.

Three kinds of magnetic couplings in the quiet solar atmosphere are highlighted and discussed, all fundamentally connected to the Lorentz force. First the coupling of the convecting and overshooting fluid in the surface layers of the Sun with the magnetic field. Here, the plasma motion provides the dominant force, which shapes the magnetic field and drives the surface dynamo. Progress in the understanding of the horizontal magnetic field is summarized and discussed. Second, the coupling between acoustic waves and the magnetic field, in particular the phenomenon of wave conversion and wave refraction. It is described how measurements of wave travel times in the atmosphere can provide information about the topography of the wave conversion zone, i.e., the surface of equal Alfv\'en and sound speed. In quiet regions, this surface separates a highly dynamic magnetic field with fast moving magnetosonic waves and shocks around and above it from the more slowly evolving field of high-beta plasma below it. Third, the magnetic field also couples to the radiation field, which leads to radiative flux channeling and increased anisotropy in the radiation field. It is shown how faculae can be understood in terms of this effect. The article starts with an introduction to the magnetic field of the quiet Sun in the light of new results from the Hinode space observatory and with a brief survey of measurements of the turbulent magnetic field with the help of the Hanle effect.

Three-dimensional magnetohydrodynamic simulations of the surface layers of the Sun intrinsically produce a predominantly horizontal magnetic field in the photosphere. This is a robust result in the sense that it arises from simulations with largely different initial and boundary conditions for the magnetic field. While the disk-center synthetic circular and linear polarization signals agree with measurements from Hinode, their center-to-limb variation sensitively depends on the height variation of the horizontal and the vertical field component and they seem to be at variance with the observed behavior.

We present 2D local box simulations of near-surface radiative magneto-convection with prescribed magnetic flux, carried out with the MHD version of the CO5BOLD code for the Sun and a solar-like star with a metal-poor chemical composition (metal abundances reduced by a factor 100, [M/H]=-2). The resulting magneto-hydrodynamical models can be used to study the influence of the metallicity on the properties of magnetized stellar atmospheres. A preliminary analysis indicates that the horizontal magnetic field component tends to be significantly stronger in the optically thin layers of metal-poor stellar atmospheres.

Observations with the Hinode space observatory led to the discovery of predominantly horizontal magnetic fields in the photosphere of the quiet internetwork region. Here we investigate realistic numerical simulations of the surface layers of the Sun with respect to horizontal magnetic fields and compute the corresponding polarimetric response in the Fe I 630 nm line pair. We find a local maximum in the mean strength of the horizontal field component at a height of around 500 km in the photosphere, where it surpasses the vertical component by a factor of 2.0 or 5.6, depending on the initial and boundary conditions. From the synthesized Stokes profiles we derive a mean horizontal field component that is, respectively, 1.6 and 4.3 times stronger than the vertical component. This is a consequence of both the intrinsically stronger flux density of, and the larger area occupied by the horizontal fields. We find that convective overshooting expels horizontal fields to the upper photosphere, making the Poynting flux positive in the photosphere, while this quantity is negative in the convectively unstable layer below it.

We study the structure of the magnetic elements in network-cell interiors. A quiet Sun area close to the disc centre was observed with the spectro-polarimeter of the Solar Optical Telescope on board the Hinode space mission, which yielded the best spatial resolution ever achieved in polarimetric data of the Fe I 630 nm line pair. For comparison and interpretation, we synthesize a similar data set from a three-dimensional magneto-hydrodynamic simulation. We find several examples of magnetic elements, either roundish (tube) or elongated (sheet), which show a central area of negative Stokes-V area asymmetry framed or surrounded by a peripheral area with larger positive asymmetry. This pattern was predicted some eight years ago on the basis of numerical simulations. Here, we observationally confirm its existence for the first time. We gather convincing evidence that this pattern of Stokes-V area asymmetry is caused by the funnel-shaped boundary of magnetic elements that separates the flux concentration from the weak-field environment. We also conclude that this kind of magnetic element of the internetwork is accompanied by electric current sheets.

In the first part of these lecture notes, new high-resolution observations of small-scale magnetic flux concentrations are presented and compared to results from new three-dimensional magnetohydrodynamic simulations. Special attention is paid to the physics of faculae and to new three-dimensional radiation magnetohydrodynamic simulations of the integral layers from the top of the convection zone to the mid-chromosphere. The second part is devoted to a few basic properties of magnetic flux tubes, which can be considered to be an abstraction of the more complicated flux concentrations known from observations and numerical simulations. We treat electrical current sheets, the mechanical equilibrium condition at magnetic interfaces, the equations for constructing a magnetohydrostatic flux tube embedded in a gravitationally stratified atmosphere, the condition of radiative equilibrium, and the condition for interchange stability.

New high-resolution observations reveal that small-scale magnetic flux concentrations have a delicate substructure on a spatial scale of 0.1''. Its basic structure can be interpreted in terms of a magnetic flux sheet or tube that vertically extends through the ambient weak-field or field-free atmosphere with which it is in mechanical equilibrium. A more refined interpretation comes from new three-dimensional magnetohydrodynamic simulations that are capable of reproducing the corrugated shape of magnetic flux concentrations and their signature in the visible continuum. Furthermore it is shown that the characteristic asymmetric shape of the contrast profile of facular granules is an effect of radiative transfer across the rarefied atmosphere of the magnetic flux concentration. I also discuss three-dimensional radiation magnetohydrodynamic simulations of the integral layers from the top of the convection zone to the mid-chromosphere. They show a highly dynamic chromospheric magnetic field, marked by rapidly moving filaments of stronger than average magnetic field that form in the compression zone downstream and along propagating shock fronts. The simulations confirm the picture of flux concentrations that strongly expand through the photosphere into a more homogeneous, space filling chromospheric field. Future directions in the simulation of small-scale magnetic fields are indicated by a few examples of very recent work.

We study the structure of the photospheric magnetic field of the quiet Sun by investigating weak spectro-polarimetric signals. We took a sequence of Stokes spectra of the Fe I 630.15 nm and 630.25 nm lines in a region of quiet Sun near the disk center, using the POLIS spectro-polarimeter at the German VTT on Tenerife. The line cores of these two lines form at different heights in the atmosphere. The 3$\sigma$ noise level of the data is about 1.8 $\times 10^{-3} I_{c}$. We present co-temporal and co-spatial Stokes-$V$ profiles of the Fe I 630 nm line pair, where the two lines show opposite polarities in a single spectrum. We compute synthetic line profiles and reproduce these spectra with a two-component model atmosphere: a non-magnetic component and a magnetic component. The magnetic component consists of two magnetic layers with opposite polarity: the upper one moves upwards while the lower one moves downward. In-between, there is a region of enhanced temperature. The Stokes-$V$ line pair of opposite polarity in a single spectrum can be understood as a magnetic reconnection event in the solar photosphere. We demonstrate that such a scenario is realistic, but the solution may not be unique.

With numerical experiments we explore the feasibility of using high frequency waves for probing the magnetic fields in the photosphere and the chromosphere of the Sun. We track a plane-parallel, monochromatic wave that propagates through a non-stationary, realistic atmosphere, from the convection-zone through the photosphere into the magnetically dominated chromosphere, where it gets refracted and reflected. We compare the wave travel time between two fixed geometrical height levels in the atmosphere (representing the formation height of two spectral lines) with the topography of the surface of equal magnetic and thermal energy density (the magnetic canopy or beta=1 contour) and find good correspondence between the two. We conclude that high frequency waves indeed bear information on the topography of the `magnetic canopy'.

It is widely believed that the heating of the chromosphere in quiet-Sun internetwork regions is provided by dissipation of acoustic waves that are excited by the convective motions close to the top of the convection zone and in the photospheric overshoot layer. This view lately became challenged by observations suggesting that the acoustic energy flux into the chromosphere is too low, by a factor of at least ten. Based on a comparison of TRACE data with synthetic image sequences for a three-dimensional simulation extending from the top layers of the convection zone to the middle chromosphere, we come to the contradicting conclusion that the acoustic flux in the model provides sufficient energy for heating the solar chromosphere of internetwork regions. The role of a weak magnetic field and associated electric current sheets is also discussed.

The aim of this work is to identify the physical processes that occur in the network and contribute to its dynamics and heating. We model the network as consisting of individual flux tubes with a non-potential field structure that are located in intergranular lanes. With a typical horizontal size of 200 km at the base of the photosphere, they expand upward and merge with their neighbors at a height of about 600 km. Above a height of approximately 1000 km the magnetic field starts to become uniform. Waves are generated in this medium by means of motions at the lower boundary. We focus on transverse driving, which generates both fast and slow waves within a flux tube and acoustic waves at the interface of the tube and the field-free medium. The acoustic waves at the interface are due to compression of the gas on one side of the flux tube and expansion on the other. These waves travel upward along the two sides of the (2D) flux tube and enter it, where they become longitudinal waves. For impulsive excitation with a time constant of 120 s, we find that a dominant feature is the creation of vortical motions that propagate upward. We have identified an efficient mechanism for the generation of longitudinal waves and shock formation in the chromospheric part of the flux concentration. We examine some broad implications of our results.