We consider the formation of solar infrared (2-6 micron) rovibrational bands of carbon monoxide (CO) in CO5BOLD 3D convection models, with the aim to refine abundances of the heavy isotopes of carbon (13C) and oxygen (18O,17O), to compare with direct capture measurements of solar wind light ions by the Genesis Discovery Mission. We find that previous, mainly 1D, analyses were systematically biased toward lower isotopic ratios (e.g., R23= 12C/13C), suggesting an isotopically "heavy" Sun contrary to accepted fractionation processes thought to have operated in the primitive solar nebula. The new 3D ratios for 13C and 18O are: R23= 91.4 +/- 1.3 (Rsun= 89.2); and R68= 511 +/- 10 (Rsun= 499), where the uncertainties are 1 sigma and "optimistic." We also obtained R67= 2738 +/- 118 (Rsun= 2632), but we caution that the observed 12C17O features are extremely weak. The new solar ratios for the oxygen isotopes fall between the terrestrial values and those reported by Genesis (R68= 530, R6= 2798), although including both within 2 sigma error flags, and go in the direction favoring recent theories for the oxygen isotope composition of Ca-Al inclusions (CAI) in primitive meteorites. While not a major focus of this work, we derive an oxygen abundance of 603 +/- 9 ppm (relative to hydrogen; 8.78 on the logarithmic H= 12 scale). That the Sun likely is lighter than the Earth, isotopically speaking, removes the necessity to invoke exotic fractionation processes during the early construction of the inner solar system.

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.

Our aim is to determine the dominant processes and timescales for the ionisation equilibrium of calcium under solar chromospheric conditions. The study is based on numerical simulations with the RADYN code, which includes hydrodynamics, radiative transfer, and a detailed non-equilibrium treatment of hydrogen, calcium, and helium. The simulations are characterised by upwards propagating shock waves, which cause strong temperature fluctuations and variations of the ionisation degree of calcium. The passage of a hot shock front leads to a strong net ionisation of Ca II, rapidly followed by net recombination. The relaxation timescale of the Ca ionisation state is found to be of the order of a few seconds at the top of the photosphere and 10 to 30 s in the upper chromosphere. Generally, the timescales are significantly reduced in the wakes of hot shock fronts. The timescales can be reliably determined from a simple analysis of the eigenvalues of the transition rate matrix. The timescales are dominated by the radiative recombination from Ca III into the metastable Ca II energy levels of the 4d 2D term. These transitions depend strongly on the density of free electrons and therefore on the (non-equilibrium) ionisation degree of hydrogen, which is the main electron donor. The ionisation/recombination timescales derived here are too long for the assumption of an instantaneous ionisation equilibrium to be valid and, on the other hand, are not long enough to warrant an assumption of a constant ionisation fraction. Fortunately, the ionisation degree of Ca II remains small in the height range, where the cores of the H, K, and the infrared triplet lines are formed. We conclude that the difference due to a detailed treatment of Ca ionisation has only negligible impact on the modelling of spectral lines of Ca II and the plasma properties under the conditions in the quiet solar chromosphere.

We use joint observations obtained with the Hinode space observatory and the Interferometric Bidimensional Spectrometer (IBIS) installed at the DST of the NSO/SP to investigate the morphology and dynamics of (a) non-magnetic and (b) magnetic regions in the fluctosphere. In inter-network regions with no significant magnetic flux contributions above the detection limit of IBIS, we find intensity structures with similar characteristics as those seen in numerical simulations by Wedemeyer-B\"ohm (2008). The magnetic flux elements in the network are stable and seem to resemble the spatially extended counterparts to the underlying photospheric magnetic elements. We will explain some of the difficulties in deriving the magnetic field vector from observations of the fluctosphere.

We propose a method to derive the line-of-sight magnetic flux density from measurements in the chromospheric Ca II IR line at 854.2 nm. The method combines two well-understood techniques, the center-of-gravity and bisector methods, in a single hybrid technique. The technique is tested with magneto-static simulations of a flux tube. We apply the method to observations with the Interferometric Bidimensional Spectrometer (IBIS) installed at the Dunn Solar Telescope of the NSO/SP to investigate the morphology of the lower chromosphere, with focus on the chromospheric counterparts to the underlying photospheric magnetic flux elements.

Recent advances in observational performance and numerical simulations have revolutionised our understanding of the solar chromosphere. This concerns in particular the structure and dynamics on small spatial and temporal scales. As a result, the picture of the solar chromosphere changed from an idealised static and plane-parallel stratification to a complex compound of intermittent domains, which are dynamically coupled to the layers below and above. In this picture, the chromosphere in a stricter sense is associated with the typical fibrillar structure shaped by magnetic fields like it is known from images taken in the H alpha line core. In internetwork regions below this layer, there exists a domain with propagating shock waves and weak magnetic fields, which both probably interact with the overlying large scale field. The existence of such a sub-canopy domain certainly depends on the properties of the overlying field. Details of the structure of the lower atmosphere can therefore be expected to vary significantly from location to location. Here, high-resolution observations, which were obtained with the CRISP filter at the Swedish Solar Telescope, are used to derive qualitative constraints for the atmospheric structure of quiet Sun regions.

Recent progress in instrumentation enables solar observations with high resolution simultaneously in the spatial, temporal, and spectral domains. We use such high-resolution observations to study small-scale structures and dynamics in the chromosphere of the quiet Sun. We analyze time series of spectral scans through the Ca II 854.2nm spectral line obtained with the CRISP instrument at the Swedish 1-m Solar Telescope. The targets are quiet Sun regions inside coronal holes close to disc-centre. The line core maps exhibit relatively few fibrils compared to what is normally observed in quiet Sun regions outside coronal holes. The time series show a chaotic and dynamic scene that includes spatially confined "swirl" events. These events feature dark and bright rotating patches, which can consist of arcs, spiral arms, rings or ring fragments. The width of the fragments typically appears to be on the order of only 0.2", which is close to the effective spatial resolution. They exhibit Doppler shifts of -2 to -4 km/s but sometimes up to -7 km/s, indicating fast upflows. The diameter of a swirl is usually of the order of 2". At the location of these swirls, the line wing and wide-band maps show close groups of photospheric bright points that move with respect to each other. A likely explanation is that the relative motion of the bright points twists the associated magnetic field in the chromosphere above. Plasma or propagating waves may then spiral upwards guided by the magnetic flux structure, thereby producing the observed intensity signature of Doppler-shifted ring fragments.

The Interferometric Bidimensional Spectrometer (IBIS) installed at the Dunn Solar Telescope of the NSO/SP is used to investigate the morphology and dynamics of the lower chromosphere and the virtually non-magnetic fluctosphere below. The study addresses in particular the structure of magnetic elements that extend into these layers. We choose different quiet Sun regions in and outside coronal holes. In inter-network regions with no significant magnetic flux contributions above the detection limit of IBIS, we find intensity structures with the characteristics of a shock wave pattern. The magnetic flux elements in the network are long lived and seem to resemble the spatially extended counterparts to the underlying photospheric magnetic elements. We suggest a modification to common methods to derive the line-of-sight magnetic field strength and explain some of the difficulties in deriving the magnetic field vector from observations of the fluctosphere.

For many years, there seemed to be significant differences between the continuum intensity distributions derived from observations and simulations of the solar photosphere. In order to settle the discussion on these apparent discrepancies, we present a detailed comparison between simulations and seeing-free observations that takes into account the crucial influence of instrumental image degradation. We use a set of images of quiet Sun granulation taken in the blue, green and red continuum bands of the Broadband Filter Imager of the Solar Optical Telescope (SOT) onboard Hinode. The images are deconvolved with Point Spread Functions (PSF) that account for non-ideal contributions due to instrumental stray-light and imperfections. In addition, synthetic intensity images are degraded with the corresponding PSFs. The results are compared with respect to spatial power spectra, intensity histograms, and the centre-to-limb variation of the intensity contrast. The observational findings are well matched with corresponding synthetic observables from three-dimensional radiation (magneto-)hydrodynamic simulations. We conclude that the intensity contrast of the solar continuum intensity is higher than usually derived from ground-based observations and is well reproduced by modern numerical simulations. Properly accounting for image degradation effects is of crucial importance for comparisons between observations and numerical models. It finally settles the traditionally perceived conflict between observations and simulations.

We present a detailed comparison between simulations and seeing-free observations that takes into account the crucial influence of instrumental image degradation. We use images of quiet Sun granulation taken in the blue, green and red continuum bands of the Broadband Filter Imager of the Solar Optical Telescope (SOT) onboard Hinode. The images are deconvolved with Point Spread Functions (PSF) that account for non-ideal contributions due to instrumental stray-light and imperfections. In addition, synthetic intensity images are degraded with the corresponding PSFs. ... Removing the influence of the PSF unveils much broader intensity distributions with a secondary component that is otherwise only visible as an asymmetry between the darker and brighter than average part of the distribution. The contrast values increase to (26.7 +/- 1.3) %, (19.4 +/- 1.4) %, and (16.6 +/- 0.7) % for blue, green, and red continuum, respectively. The power spectral density of the images exhibits a pronounced peak at spatial scales characteristic for the granulation pattern and a steep decrease towards smaller scales. The observational findings like the absolute values and centre-to-limb variation of the intensity contrast, intensity histograms, and power spectral density are well matched with corresponding synthetic observables from three-dimensional radiation (magneto-)hydrodynamic simulations. We conclude that the intensity contrast of the solar continuum intensity is higher than usually derived from ground-based observations and is well reproduced by modern radiation (magneto-)hydrodynamic models. Properly accounting for image degradation effects is of crucial importance for comparisons between observations and numerical models.

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.

The atmosphere of the Sun is characterized by a complex interplay of competing physical processes: convection, radiation, conduction, and magnetic fields. The most obvious imprint of the solar convection and its overshooting in the low atmosphere is the granulation pattern. Beside this dominating scale there is a more or less smooth distribution of spatial scales, both towards smaller and larger scales, making the Sun essentially a multi-scale object. Convection and overshooting give the photosphere its face but also act as drivers for the layers above, namely the chromosphere and corona. The magnetic field configuration effectively couples the atmospheric layers on a multitude of spatial scales, for instance in the form of loops that are anchored in the convection zone and continue through the atmosphere up into the chromosphere and corona. The magnetic field is also an important structuring agent for the small, granulation-size scales, although (hydrodynamic) shock waves also play an important role -- especially in the internetwork atmosphere where mostly weak fields prevail. Based on recent results from observations and numerical simulations, we attempt to present a comprehensive picture of the atmosphere of the quiet Sun as a highly intermittent and dynamic system.

The combined PSF of the BFI and the SOT onboard the Hinode spacecraft is investigated. Observations of the Mercury transit from November 2006 and the solar eclipse(s) from 2007 are used to determine the PSFs of SOT for the blue, green, and red continuum channels of the BFI. For each channel large grids of theoretical point spread functions are calculated by convolution of the ideal diffraction-limited PSF and Voigt profiles. These PSFs are applied to artificial images of an eclipse and a Mercury transit. The comparison of the resulting artificial intensity profiles across the terminator and the corresponding observed profiles yields a quality measure for each case. The optimum PSF for each observed image is indicated by the best fit. The observed images of the Mercury transit and the eclipses exhibit a clear proportional relation between the residual intensity and the overall light level in the telescope. In addition there is a anisotropic stray-light contribution. ... BFI/SOT operate close to the diffraction limit and have only a rather small stray-light contribution. The FWHM of the PSF is broadened by only ~1% with respect to the diffraction-limited case, while the overall Strehl ratio is ~ 0.8. In view of the large variations -- best seen in the residual intensities of eclipse images -- and the dependence on the overall light level and position in the FOV, a range of PSFs should be considered instead of a single PSF per wavelength. The individual PSFs of that range allow then the determination of error margins for the quantity under investigation. Nevertheless the stray-light contributions are here found to be best matched with Voigt functions with the parameters sigma = 0."008 and gamma = 0."004, 0."005, and 0."006 for the blue, green, and red continuum channels, respectively.

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.

The chromosphere of the quiet Sun is a highly intermittent and dynamic phenomenon. Three-dimensional radiation (magneto-)hydrodynamic simulations exhibit a mesh-like pattern of hot shock fronts and cool expanding post-shock regions in the sub-canopy part of the inter-network. This domain might be called "fluctosphere". The pattern is produced by propagating shock waves, which are excited at the top of the convection zone and in the photospheric overshoot layer. New high-resolution observations reveal a ubiquitous small-scale pattern of bright structures and dark regions in-between. Although it qualitatively resembles the picture seen in models, more observations - e.g. with the future ALMA - are needed for thorough comparisons with present and future models. Quantitative comparisons demand for synthetic intensity maps and spectra for the three-dimensional (magneto-)hydrodynamic simulations. The necessary radiative transfer calculations, which have to take into account deviations from local thermodynamic equilibrium, are computationally very involved so that no reliable results have been produced so far. Until this task becomes feasible, we have to rely on careful qualitative comparisons of simulations and observations. Here we discuss what effects have to be considered for such a comparison. Nevertheless we are now on the verge of assembling a comprehensive picture of the solar chromosphere in inter-network regions as dynamic interplay of shock waves and structuring and guiding magnetic fields.

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 role of carbon monoxide as a cooling agent for the thermal structure of the mid-photospheric to low-chromospheric layers of the solar atmosphere in internetwork regions is investigated. The treatment of radiative cooling via spectral lines of carbon monoxide (CO) has been added to the radiation chemo-hydrodynamics code CO5BOLD. [...] The CO opacity indeed causes additional cooling at the fronts of propagating shock waves in the chromosphere. There, the time-dependent approach results in a higher CO number density compared to the equilibrium case and hence in a larger net radiative cooling rate. The average gas temperature stratification of the model atmosphere, however, is only reduced by roughly 100 K. Also the temperature fluctuations and the CO number density are only affected to small extent. A numerical experiment without dynamics shows that the CO cooling process works in principle and drives the atmosphere to a cool radiative equilibrium state. At chromospheric heights, the radiative relaxation of the atmosphere to a cool state takes several 1000 s. The CO cooling process thus would seem to be too slow compared to atmospheric dynamics to be responsible for the very cool temperature regions observed in the solar atmosphere. The hydrodynamical timescales in our solar atmosphere model are much too short to allow for the radiative relaxation to a cool state, thus suppressing the potential thermal instability due to carbon monoxide as a cooling agent. Apparently, the thermal structure and dynamics of the outer model atmosphere are instead determined primarily by shock waves.

The chromosphere in internetwork regions of the quiet Sun was regarded as a static and homogeneous layer for a long time. Thanks to advances in observations and numerical modelling, the wave nature of these atmospheric regions received increasing attention during the last decade. Recent three-dimensional radiation magnetohydrodynamic simulations with CO5BOLD feature the chromosphere of internetwork regions as a dynamic and intermittent phenomenon. It is a direct product of interacting waves that form a mesh-like pattern of hot shock fronts and cool post-shock regions. The waves are excited self-consistently at the top of the convection zone. In the middle chromosphere above an average height of 1000 km, plasma beta gets larger than one and magnetic fields become more important. The model chromosphere exhibits a magnetic field that is much more homogeneous than in the layers below and evolves much faster. That includes fast propagating (MHD) waves. Further improvements of the simulations like time-dependent hydrogen ionisation are currently in progress. This class of models is capable of explaining apparently contradicting diagnostics such as carbon monoxide and UV emission at the same time.

In this work we investigate the dynamic behavior of inter-network regions of the solar chromosphere. We observed the chromosphere of the quiet Sun using a narrow-band Lyot filter centered at the Ca II K 2v emission peak with a bandpass of 0.3A. We achieved a spatial resolution of on average 0.7" at a cadence of 10s. In the inter-network we find a mesh-like pattern that features bright grains at the vertices. The pattern has a typical spatial scale of 1.95" and a mean evolution time scale of 53s with a standard deviation of 10s. A comparison of our results with a recent three-dimensional radiation hydrodynamical model implies that the observed pattern is of chromospheric origin. The measured time scales are not compatible with those of reversed granulation in the photosphere although the appearance is similar. A direct comparison between network and inter-network structure shows that their typical time scales differ by at least a factor of two. The existence of a rapidly evolving small-scale pattern in the inter-network regions supports the picture of the lower chromosphere as a highly dynamical and intermittent phenomenon.