We revisit the problem of how energy transfer through the turbulent cascade operates in compressible hydrodynamic turbulence. In general, there is no conservative compressible cascade since the kinetic and internal energy reservoirs can exchange energy through pressure dilatation. Moreover, statistically stationary turbulence at high Mach number can only be maintained in nearly isothermal gas, i.e., if excess heat produced by shock compression and kinetic energy dissipation is continuously removed from the system. We mimic this process by a linear cooling term in numerical simulations of turbulence driven by stochastic forcing. This allows us to investigate turbulence statistics for a broad range of Mach numbers. We compute the rate of change of kinetic and internal energy in wavenumber shells caused by advective, compressive, and pressure dilatation effects and constrain power-law fits to compressible turbulence energy spectra to a range of wavenumbers in which the total energy transfer is close to zero. The resulting scaling exponents are significantly affected by the forcing. Depending on the root mean square Mach number, we find a nearly constant advective component of the cross-scale flux of kinetic energy at intermediate wavenumbers for particular mixtures of solenoidal and compressive modes in the forcing. This suggests the existence of a natural, Mach number dependent mixture of forcing modes. Our findings also support an advection-dominated regime at high Mach numbers with specific scaling exponents (Burgers scaling for the pure velocity fluctuation $u$ and Kolmogorov scaling for the mass-weighted variable $v=\rho^{1/3}u$).

Magnetic fields are an elemental part of the interstellar medium in galaxies. However, their impact on gas dynamics and star formation in galaxies remains controversial. We use a suite of global magnetohydrodynamical simulations of isolated disc galaxies to study the influence of magnetic fields on the diffuse and dense gas in the discs. We find that the magnetic field acts in multiple ways. Stronger magnetised discs fragment earlier due to the shorter growth time of the Parker instability. Due to the Parker instability in the magnetised discs we also find cold ($T<50\,\mathrm{K}$) and dense ($n\sim10^3-10^4\,\mathrm{cm}^{-3}$) gas several hundred pc above/below the midplane without any form of stellar feedback. In addition, magnetic fields change the fragmentation pattern. While in the hydrodynamical case, the disc breaks up into ring-like structures, magnetised discs show the formation of filamentary entities that extent both in the azimuthal and radial direction. These kpc scale filaments become magnetically (super-)critical very quickly and allow for the rapid formation of massive giant molecular clouds. Our simulations suggest that major differences in the behaviour of star formation - due to a varying magnetisation - in galaxies could arise.

Observations reveal the presence of supermassive black holes (SMBH) as early as ~700 million years after the Big Bang. Their formation path is still subject to current debate. We explore the influence of magnetic fields, which are strongly amplified via the turbulent small-scale dynamo, on the formation of SMBH seeds within the direct collapse scenario. In this study, we perform for the first time cosmological magnetohydrodynamic large eddy simulations that employ a model for unresolved, compressible MHD turbulence. In total we perform 36 simulations for 9 haloes each with two different initial magnetic field strengths, and with and without employing the unresolved turbulence model. We make use of the adaptive mesh refinement approach to achieve an effective spatial resolution of less than one proper astronomical unit. We consider a regime where cooling is regulated by atomic hydrogen and the molecular hydrogen gets dissociated by a strong radiation field. Our main finding is that the majority of the gas properties in the haloes at the final output are predominantly determined by the run-away gravitational collapse. Turbulence is supersonic and super-Alfv\'enic in all cases, and magnetic fields are amplified to an approximately dynamically relevant regime. Finally, fragmentation during the collapse is intermittent and mass accretion rates range from 0.2-3 Msun/yr. This suggests that the presence of strongly amplified magnetic fields and turbulence provides additional pressure support on small scales and make the direct collapse a viable scenario for the formation of massive objects under the required ambient conditions.

This paper describes the Polarimetric and Helioseismic Imager on the Solar Orbiter mission (SO/PHI), the first magnetograph and helioseismology instrument to observe the Sun from outside the Sun-Earth line. It is the key instrument meant to address the top-level science question: How does the solar dynamo work and drive connections between the Sun and the heliosphere? SO/PHI will also play an important role in answering the other top-level science questions of Solar Orbiter, as well as hosting the potential of a rich return in further science. SO/PHI measures the Zeeman effect and the Doppler shift in the FeI 617.3nm spectral line. To this end, the instrument carries out narrow-band imaging spectro-polarimetry using a tunable LiNbO_3 Fabry-Perot etalon, while the polarisation modulation is done with liquid crystal variable retarders (LCVRs). The line and the nearby continuum are sampled at six wavelength points and the data are recorded by a 2kx2k CMOS detector. To save valuable telemetry, the raw data are reduced on board, including being inverted under the assumption of a Milne-Eddington atmosphere, although simpler reduction methods are also available on board. SO/PHI is composed of two telescopes; one, the Full Disc Telescope (FDT), covers the full solar disc at all phases of the orbit, while the other, the High Resolution Telescope (HRT), can resolve structures as small as 200km on the Sun at closest perihelion. The high heat load generated through proximity to the Sun is greatly reduced by the multilayer-coated entrance windows to the two telescopes that allow less than 4% of the total sunlight to enter the instrument, most of it in a narrow wavelength band around the chosen spectral line.

Convective motions in the solar atmosphere cause spectral lines to become asymmetric and shifted in wavelength. For photospheric lines, this differential Doppler shift varies from the solar disk center to the limb. Precise and comprehensive observations of the convective blueshift and its center-to-limb variation improve our understanding of the atmospheric hydrodynamics and ensuing line formation, and provide the basis to refine 3D models of the solar atmosphere. We performed systematical spectroscopic measurements of the convective blueshift of the quiet Sun with the Laser Absolute Reference Spectrograph (LARS) at the German Vacuum Tower Telescope. The spatial scanning of the solar disk covered four radial (meridional and equatorial) axes. The high-resolution spectra of 26 photospheric to chromospheric lines in the visible range were calibrated with a laser frequency comb to absolute wavelengths at the 1m/s accuracy. Applying ephemeris and reference corrections, the bisector analysis provided line asymmetries and Doppler shifts with an uncertainty of only few m/s. To allow for a comparison with other observations, we convolved the results to lower spectral resolutions. Typically, a blueshifted "C"-shaped curve at disk center transforms into a less blueshifted "\"-shape toward the solar limb. The comparison of all lines reveals the systematic dependence of the convective blueshift on the line depth. Synthetic models yield considerable deviations from the observed center-to-limb variation. The obtained Doppler shifts of the quiet Sun can serve as an absolute reference for other observations, the relative calibration of Dopplergrams, and the necessary refinement of atmospheric models. Based on this, the development of high-precision models of stellar surface convection will advance the detection of (potentially habitable) exoplanets by radial velocity measurements.

Granular convective motions reach into the lower solar atmosphere, typically causing photospheric spectral lines to exhibit a differential line shift. This Doppler shift to shorter wavelength is commonly known as convective blueshift. We performed systematic observations of the quiet Sun with the Laser Absolute Reference Spectrograph (LARS) at the German Vacuum Tower Telescope. The solar disk was scanned along the meridian and the equator, from the disk center toward the limb. The solar spectrum around 6173{\AA} was calibrated with a laser frequency comb on an absolute wavelength scale with an accuracy of a few m/s. We applied a bisector analysis on the spectral lines to reveal the changes of convective blueshift and line asymmetry at different heliocentric positions. Being a signature for convective motions, the bisector curve of Fe i 6173.3{\AA} describes a C-shape at disk center. When approaching the solar limb, the bisector transforms into a \-shape. The analysis of the time- and bisector-averaged line shifts yields three distinct results. Firstly, the center-to-limb variation of Doppler velocities measured with LARS reveals a significant discrepancy (up to 200m/s) to the full-disk Dopplergrams of the Helioseismic and Magnetic Imager (HMI). Secondly, we obtained a significant decrease of convective blueshift toward the solar limb. Thirdly, the line-of-sight effect of solar activity leads to a scatter of up to 100m/s at intermediate heliocentric positions. The accurate observation of the absolute convective blueshift with LARS allows the identification of systematic discrepancy with Doppler velocities measured by HMI. The center-to-limb variation of HMI suffers from an additional blueshift toward the limb that is incompatible with our results. LARS measurements can be taken as reference for the correction of systematic errors in the synoptic HMI Dopplergrams.

The solar magnetic field is responsible for all aspects of solar activity. Sunspots are the main manifestation of the ensuing solar activity. Combining high-resolution and synoptic observations has the ambition to provide a comprehensive description of the sunspot growth and decay processes. Active region NOAA 12396 emerged on 2015 August 3 and was observed three days later with the 1.5-meter GREGOR solar telescope on 2015 August 6. High-resolution spectropolarimetric data from the GREGOR Infrared Spectrograph (GRIS) are obtained in the photospheric Si I $\lambda$ 1082.7 nm and Ca I $\lambda$1083.9 nm lines, together with the chromospheric He I $\lambda$1083.0 nm triplet. These near-infrared spectropolarimetric observations were complemented by synoptic line-of-sight magnetograms and continuum images of the Helioseismic and Magnetic Imager (HMI) and EUV images of the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO).

Observations show that galaxies and their interstellar media are pervaded by strong magnetic fields with energies in the diffuse component being at least comparable to the thermal and even as large or larger than the turbulent energy. Such strong magnetic fields prevent the formation of stars because patches of the interstellar medium are magnetically subcritical. Here we present the results from global numerical simulations of strongly magnetised and self-gravitating galactic discs, which show that the buoyancy of the magnetic field due to the Parker instability leads at first to the formation of giant filamentary regions. These filamentary structures become gravitationally unstable and fragment into $\sim10^5 M_{\odot}$ clouds that attract kpc long, coherent filamentary flows that build them into GMCs. Our results thus provide a solution to the long-standing problem of how the transition from sub- to supercritical regions in the interstellar medium proceeds.

We present high-precision spectro-polarimetric data with high spatial resolution (0.4$''$) of the very quiet Sun at 1.56$\mu$m obtained with the GREGOR telescope to shed some light on this complex magnetism. Half of our observed quiet-Sun region is better explained by magnetic substructure within the resolution element. However, we cannot distinguish whether this substructure comes from gradients of the physical parameters along the line of sight or from horizontal gradients (across the surface). In these pixels, a model with two magnetic components is preferred, and we find two distinct magnetic field populations. The population with the larger filling factor has very weak ($\sim$150 G) horizontal fields similar to those obtained in previous works. We demonstrate that the field vector of this population is not constrained by the observations, given the spatial resolution and polarimetric accuracy of our data. The topology of the other component with the smaller filling factor is constrained by the observations for field strengths above 250 G: we infer hG fields with inclinations and azimuth values compatible with an isotropic distribution. The filling factors are typically below 30\%. We also find that the flux of the two polarities is not balanced. From the other half of the observed quiet-Sun area $\sim$50\% are two-lobed Stokes $V$ profiles, meaning that 23\% of the field of view can be adequately explained with a single constant magnetic field embedded in a non-magnetic atmosphere. The magnetic field vector and filling factor are reliable inferred in only 50\% based on the regular profiles. Therefore, 12\% of the field of view harbour hG fields with filling factors typically below 30\%. At our present spatial resolution, 70\% of the pixels apparently are non-magnetised.

In sunspot umbrae, convection is largely suppressed by the strong magnetic field. Previous measurements reported on negligible convective flows in umbral cores. Based on this, numerous studies have taken the umbra as zero reference to calculate Doppler velocities of the ambient active region. To clarify the amount of convective motion in the darkest part of umbrae, we directly measured Doppler velocities with an unprecedented accuracy and precision. We performed spectroscopic observations of sunspot umbrae with the Laser Absolute Reference Spectrograph (LARS) at the German Vacuum Tower Telescope. A laser frequency comb enabled the calibration of the high-resolution spectrograph and absolute wavelength positions. A thorough spectral calibration, including the measurement of the reference wavelength, yielded Doppler shifts of the spectral line Ti i 5713.9 {\AA} with an uncertainty of around 5 m s-1. The measured Doppler shifts are a composition of umbral convection and magneto-acoustic waves. For the analysis of convective shifts, we temporally average each sequence to reduce the superimposed wave signal. Compared to convective blueshifts of up to -350 m s-1 in the quiet Sun, sunspot umbrae yield a strongly reduced convective blueshifts around -30 m s-1. {W}e find that the velocity in a sunspot umbra correlates significantly with the magnetic field strength, but also with the umbral temperature defining the depth of the titanium line. The vertical upward motion decreases with increasing field strength. Extrapolating the linear approximation to zero magnetic field reproduces the measured quiet Sun blueshift. Simply taking the sunspot umbra as a zero velocity reference for the calculation of photospheric Dopplergrams can imply a systematic velocity error.

We observed an arch filament system (AFS) in a sunspot group with the GREGOR Infrared Spectrograph attached to the GREGOR solar telescope. The AFS was located between the leading sunspot of negative polarity and several pores of positive polarity forming the following part of the sunspot group. We recorded five spectro-polarimetric scans of this region. The spectral range included the spectral lines Si I 1082.7 nm, He I 1083.0 nm, and Ca I 1083.9 nm. In this work we concentrate on the silicon line which is formed in the upper photosphere. The line profiles are inverted with the code `Stokes Inversion based on Response functions' to obtain the magnetic field vector. The line-of-sight velocities are determined independently with a Fourier phase method. Maximum velocities are found close to the ends of AFS fibrils. These maximum values amount to 2.4 km/s next to the pores and to 4 km/s at the sunspot side. Between the following pores, we encounter an area of negative polarity that is decreasing during the five scans. We interpret this by new emerging positive flux in this area canceling out the negative flux. In summary, our findings confirm the scenario that rising magnetic flux tubes cause the AFS.

We investigate a region of the quiet solar photosphere exhibiting three-lobed Stokes V profiles in the Fe I spectral line at 15648 Angstroem. The data were acquired with the GRIS spectropolarimeter attached to the GREGOR telescope. We aim at investigating the thermal, kinematic and magnetic properties of the atmosphere responsible for these measured complex signals. The SIR inversion code is employed to retrieve the physical parameters of the lower solar photosphere from the observed polarization signals. We follow two different approaches. On the one hand, we consider that the multi-lobe circular polarization signals are only produced by the line-of-sight variation of the physical parameters. We therefore invert the data assuming a single atmospheric component that occupies the entire resolution element in the horizontal plane and where the physical parameters vary with optical depth (i.e., line-of-sight). On the other hand, we consider that the multi-lobe circular polarization signals are produced not by the optical depth variations of the physical parameters but instead by their horizontal variations. Here we invert the data assuming that the resolution element is occupied by two different atmospheric components where the kinematic and magnetic properties are constant along the line-of-sight. Both approaches reveal some common features about the topology responsible for the observed three-lobed Stokes V signals: both a strong (>1000 Gauss) and a very weak (< 10 Gauss) magnetic field with opposite polarities and harboring flows directed in opposite directions must co-exist (either vertically or horizontally interlaced) within the resolution element.

The solar convection manifests as granulation and intergranulation at the solar surface. In the photosphere, convective motions induce differential Doppler shifts to spectral lines. The observed convective blueshift varies across the solar disk. We focus on the impact of solar convection on the atmosphere and aim to resolve its velocity stratification in the photosphere. We performed high-resolution spectroscopic observations of the solar spectrum in the 6302\,\AA\ range with the Laser Absolute Reference Spectrograph (LARS) at the Vacuum Tower Telescope. A laser frequency comb enabled the calibration of the spectra to an absolute wavelength scale with an accuracy of ${\rm 1\,m\,s^{-1}}$. We systematically scanned the Quiet Sun from disk center to the limb at ten selected heliocentric positions. The analysis included 99 time sequences of up to 20\,min in length. By means of ephemeris and reference corrections, we translated wavelength shifts into absolute line-of-sight velocities. A bisector analysis on the line profiles yielded the shapes and convective shifts of seven photospheric lines. At disk center, the bisector profiles of the iron lines feature a pronounced C-shape with maximum convective blueshifts of up to ${\rm -450\,m\,s^{-1}}$ in the spectral line wings. Toward the solar limb, the bisectors change into a \textbackslash-shape with a saturation in the line core at a redshift of ${\rm +100\,m\,s^{-1}}$. The center-to-limb variation of the line core velocities shows a slight increase in blueshift when departing the disk center for larger heliocentric angles. This increase in blueshift is more pronounced for the magnetically less active meridian than for the equator. Toward the solar limb, the blueshift decreases and can turn into a redshift. Best spectroscopic measurements enabled the accurate determination of absolute convective shifts in the solar photosphere.

The observed magnetic field on the solar surface is characterized by a very complex spatial and temporal behavior. Although feature-tracking algorithms have allowed us to deepen our understanding of this behavior, subjectivity plays an important role in the identification and tracking of such features. In this paper, we continue studies Gorobets, A. Y., Borrero, J. M., & Berdyugina, S. 2016, ApJL, 825, L18 of the temporal stochasticity of the magnetic field on the solar surface without relying either on the concept of magnetic features or on subjective assumptions about their identification and interaction. We propose a data analysis method to quantify fluctuations of the line-of-sight magnetic field by means of reducing the temporal field's evolution to the regular Markov process. We build a representative model of fluctuations converging to the unique stationary (equilibrium) distribution in the long time limit with maximum entropy. We obtained different rates of convergence to the equilibrium at fixed noise cutoff for two sets of data. This indicates a strong influence of the data spatial resolution and mixing-polarity fluctuations on the relaxation process. The analysis is applied to observations of magnetic fields of the relatively quiet areas around an active region carried out during the second flight of the Sunrise/IMaX and quiet Sun areas at the disk center from the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory satellite.

Magnetic fields, compressibility and turbulence are important factors in many terrestrial and astrophysical processes. While energy dynamics, i.e. how energy is transferred within and between kinetic and magnetic reservoirs, has been previously studied in the context of incompressible magnetohydrodynamic (MHD) turbulence, we extend shell-to-shell energy transfer analysis to the compressible regime. We derive four new transfer functions specifically capturing compressibility effects in the kinetic and magnetic cascade, and capturing energy exchange via magnetic pressure. To illustrate their viability, we perform and analyze four simulations of driven isothermal MHD turbulence in the sub- and supersonic regime with two different codes. On the one hand, our analysis reveals robust characteristics across regime and numerical method. For example, energy transfer between individual scales is local and forward for both cascades with the magnetic cascade being stronger than the kinetic one. Magnetic tension and magnetic pressure related transfers are less local and weaker than the cascades. We find no evidence for significant nonlocal transfer. On the other hand, we show that certain functions, e.g., the compressive component of the magnetic energy cascade, exhibit a more complex behavior that varies both with regime and numerical method. Having established a basis for the analysis in the compressible regime, the method can now be applied to study a broader parameter space.

LARS is an Absolute Reference Spectrograph designed for ultra-precise solar observations. The high-resolution echelle spectrograph of the Vacuum Tower Telescope is supported by a state-of-the-art laser frequency comb to calibrate the solar spectrum on an absolute wavelength scale. In this article, we describe the scientific instrument and focus on the upgrades in the last two years to turn the prototype into a turn-key system. The pursued goal was to improve the short-term and long-term stability of the systems, and enable a user-friendly and more versatile operation of the instrument. The first upgrade involved the modernization of the frequency comb. The Fabry-Perot cavities were adjusted to filter to a repetition frequency of 8GHz. A technologically matured photonic crystal fiber was implemented for spectral broadening. The second, quite recent upgrade was performed on the optics feeding the sunlight into a single-mode fiber connected to the spectrograph. A motorized translation stage was deployed to allow the automated selection of three different fields-of-view with diameters of 1", 3", and 10" for the analysis of the solar spectrum. The successful upgrades allow for long-term observations of up to several hours per day with a stable spectral accuracy of 1 m/s limited by the spectrograph. Stable, user-friendly operation of the instrument is supported. The selection of the pre-aligned fiber to change the field of view can now be done within seconds. LARS offers the possibility to observe absolute wavelength positions of spectral lines and Doppler velocities in the solar atmosphere. First results demonstrate the capabilities of the instrument for solar science. The accurate measurement of the solar convection, p-modes, and atmospheric waves will enhance our knowledge of the solar atmosphere and its physical conditions to improve current atmospheric models.

Major mergers are considered to be a significant source of turbulence in clusters. We performed a numerical simulation of a major merger event using nested-grid initial conditions, adaptive mesh refinement, radiative cooling of primordial gas, and a homogeneous ultraviolet background. By calculating the microscopic viscosity on the basis of various theoretical assumptions and estimating the Kolmogorov length from the turbulent dissipation rate computed with a subgrid-scale model, we are able to demonstrate that most of the warm-hot intergalactic medium can sustain a fully turbulent state only if the magnetic suppression of the viscosity is considerable. Accepting this as premise, it turns out that ratios of turbulent and thermal quantities change only little in the course of the merger. This confirms the tight correlations between the mean thermal and non-thermal energy content for large samples of clusters in earlier studies, which can be interpreted as second self-similarity on top of the self-similarity for different halo masses. Another long-standing question is how and to which extent turbulence contributes to the support of the gas against gravity. From a global perspective, the ratio of turbulent and thermal pressures is significant for the clusters in our simulation. On the other hand, a local measure is provided by the compression rate, i.e. the growth rate of the divergence of the flow. Particularly for the intracluster medium, we find that the dominant contribution against gravity comes from thermal pressure, while compressible turbulence effectively counteracts the support. For this reason it appears to be too simplistic to consider turbulence merely as an effective enhancement of thermal energy.

Convective flows are known as the prime means of transporting magnetic fields on the solar surface. Thus, small magnetic structures are good tracers of the turbulent flows. We study the migration and dispersal of magnetic bright features (MBFs) in intergranular areas observed at high spatial resolution with Sunrise/IMaX. We describe the flux dispersal of individual MBFs as a diffusion process whose parameters are computed for various areas in the quiet Sun and the vicinity of active regions from seeing-free data. We find that magnetic concentrations are best described as random walkers close to network areas (diffusion index, gamma=1.0), travelers with constant speeds over a supergranule (gamma=1.9-2.0), and decelerating movers in the vicinity of flux emergence and/or within active regions (gamma=1.4-1.5). The three types of regions host MBFs with mean diffusion coefficients of 130 km^2/s, 80-90 km^2/s, and 25-70 km^2/s, respectively. The MBFs in these three types of regions are found to display a distinct kinematic behavior at a confidence level in excess of 95%.

We present observations of transverse oscillations in slender Ca II H fibrils (SCFs) in the lower solar chromosphere. We use a 1 hr long time series of high- (spatial and temporal-) resolution seeing-free observations in a 0.11 nm wide passband covering the line core of Ca II H 396.9 nm from the second flight of the Sunrise balloon-borne solar observatory. The entire field of view, spanning the polarity inversion line of an active region close to the solar disk center, is covered with bright, thin, and very dynamic fine structures. Our analysis reveals the prevalence of transverse waves in SCFs with median amplitudes and periods on the order of 2.4+-0.8 km/s and 83+-29 s, respectively (with standard deviations given as uncertainties). We find that the transverse waves often propagate along (parts of) the SCFs with median phase speeds of 9+-14 km/s. While the propagation is only in one direction along the axis in some of the SCFs, propagating waves in both directions, as well as standing waves are also observed. The transverse oscillations are likely Alfvenic and are thought to be representative of magnetohydrodynamic kink waves. The wave propagation suggests that the rapid high-frequency transverse waves, often produced in the lower photosphere, can penetrate into the chromosphere with an estimated energy flux of ~ 15 kW/m^2. Characteristics of these waves differ from those reported for other fibrillar structures, which, however, were observed mainly in the upper solar chromosphere.

A dense forest of slender bright fibrils near a small solar active region is seen in high-quality narrowband Ca II H images from the SuFI instrument onboard the Sunrise balloon-borne solar observatory. The orientation of these slender Ca II H fibrils (SCF) overlaps with the magnetic field configuration in the low solar chromosphere derived by magnetostatic extrapolation of the photospheric field observed with Sunrise/IMaX and SDO/HMI. In addition, many observed SCFs are qualitatively aligned with small-scale loops computed from a novel inversion approach based on best-fit numerical MHD simulation. Such loops are organized in canopy-like arches over quiet areas that differ in height depending on the field strength near their roots.

How and where are coronal loops rooted in the solar lower atmosphere? The details of the magnetic environment and its evolution at the footpoints of coronal loops are crucial to understanding the processes of mass and energy supply to the solar corona. To address the above question, we use high-resolution line-of-sight magnetic field data from the Imaging Magnetograph eXperiment instrument on the SUNRISE balloon-borne observatory and coronal observations from the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory of an emerging active region. We find that the coronal loops are often rooted at the locations with minor small-scale but persistent opposite-polarity magnetic elements very close to the larger dominant polarity. These opposite-polarity small-scale elements continually interact with the dominant polarity underlying the coronal loop through flux cancellation. At these locations we detect small inverse Y-shaped jets in chromospheric Ca II H images obtained from the SUNRISE Filter Imager during the flux cancellation. Our results indicate that magnetic flux cancellation and reconnection at the base of coronal loops due to mixed polarity fields might be a crucial feature for the supply of mass and energy into the corona.

Large eddy simulations (LES) are a powerful tool in understanding processes that are inaccessible by direct simulations due to their complexity, for example, in the highly turbulent regime. However, their accuracy and success depends on a proper subgrid-scale (SGS) model that accounts for the unresolved scales in the simulation. We evaluate the applicability of two traditional SGS models, namely the eddy-viscosity (EV) and the scale-similarity (SS) model, and one recently proposed nonlinear (NL) SGS model in the realm of compressible MHD turbulence. Using 209 simulations of decaying, supersonic (initial sonic Mach number of ~3) MHD turbulence with a shock-capturing scheme and varying resolution, SGS model and filter, we analyze the ensemble statistics of kinetic and magnetic energy spectra and structure functions. Furthermore, we compare the temporal evolution of lower and higher order statistical moments of the spatial distributions of kinetic and magnetic energy, vorticity, current density, and dilatation magnitudes. We find no statistical influence on the evolution of the flow by any model if grid-scale quantities are used to calculate SGS contributions. In addition, the SS models, which employ an explicit filter, have no impact in general. On the contrary, both EV and NL models change the statistics if an explicit filter is used. For example, they slightly increase the dissipation on the smallest scales. We demonstrate that the nonlinear model improves higher order statistics already with a small explicit filter, i.e. a three-point stencil. The results of e.g. the structure functions or the skewness and kurtosis of the current density distribution are closer to the ones obtained from simulations at higher resolution. We conclude that the nonlinear model with a small explicit filter is suitable for application in more complex scenarios when higher order statistics are important.

Magneto-static models may overcome some of the issues facing force-free magnetic field extrapolations. So far they have seen limited use and have faced problems when applied to quiet-Sun data. Here we present a first application to an active region. We use solar vector magnetic field measurements gathered by the IMaX polarimeter during the flight of the \sunrise{} balloon-borne solar observatory in June 2013 as boundary condition for a magneto-static model of the higher solar atmosphere above an active region. The IMaX data are embedded in active region vector magnetograms observed with SDO/HMI. This work continues our magneto-static extrapolation approach, which has been applied earlier ({\it Paper I}) to a quiet Sun region observed with \sunrise{} I. In an active region the signal-to-noise-ratio in the measured Stokes parameters is considerably higher than in the quiet Sun and consequently the IMaX measurements of the horizontal photospheric magnetic field allow us to specify the free parameters of the model in a special class of linear magneto-static equilibria. The high spatial resolution of IMaX (110-130 km, pixel size 40 km) enables us to model the non-force-free layer between the photosphere and the mid chromosphere vertically by about 50 grid points. In our approach we can incorporate some aspects of the mixed beta layer of photosphere and chromosphere, e.g., taking a finite Lorentz force into account, which was not possible with lower resolution photospheric measurements in the past. The linear model does not, however, permit to model intrinsic nonlinear structures like strongly localized electric currents.

We report the detection of oscillations in slender Ca II H fibrils (SCFs) from high-resolution observations acquired with the Sunrise balloon-borne solar observatory. The SCFs show obvious oscillations in their intensity, but also their width. The oscillatory behaviors are investigated at several positions along the axes of the SCFs. A large majority of fibrils show signs of oscillations in intensity. Their periods and phase speeds are analyzed using a wavelet analysis. The width and intensity perturbations have overlapping distributions of the wave period. The obtained distributions have median values of the period of $32\pm17$s and $36\pm25$s, respectively. We find that the fluctuations of both parameters propagate in the SCFs with speeds of ${11}^{+49}_{-11}$ km/s and ${15}^{+34}_{-15}$ km/s, respectively. Furthermore, the width and intensity oscillations have a strong tendency to be either in anti-phase, or, to a smaller extent, in phase. This suggests that the oscillations of both parameters are caused by the same wave mode and that the waves are likely propagating. Taking all the evidence together, the most likely wave mode to explain all measurements and criteria is the fast sausage mode.

A new generation of solar instruments provides improved spectral, spatial, and temporal resolution, thus facilitating a better understanding of dynamic processes on the Sun. High-resolution observations often reveal multiple-component spectral line profiles, e.g., in the near-infrared He I 10830 \AA\ triplet, which provides information about the chromospheric velocity and magnetic fine structure. We observed an emerging flux region, including two small pores and an arch filament system, on 2015 April 17 with the 'very fast spectroscopic mode' of the GREGOR Infrared Spectrograph (GRIS) situated at the 1.5-meter GREGOR solar telescope at Observatorio del Teide, Tenerife, Spain. We discuss this method of obtaining fast (one per minute) spectral scans of the solar surface and its potential to follow dynamic processes on the Sun. We demonstrate the performance of the 'very fast spectroscopic mode' by tracking chromospheric high-velocity features in the arch filament system.

The SUNRISE balloon-borne solar observatory, consisting of a 1~m aperture telescope that provided a stabilized image to a UV filter imager and an imaging vector polarimeter, carried out its second science flight in June 2013. It provided observations of parts of active regions at high spatial resolution, including the first high-resolution images in the Mg~{\sc ii}~k line. The obtained data are of very high quality, with the best UV images reaching the diffraction limit of the telescope at 3000~\AA\ after Multi-Frame Blind Deconvolution reconstruction accounting for phase-diversity information. Here a brief update is given of the instruments and the data reduction techniques, which includes an inversion of the polarimetric data. Mainly those aspects that evolved compared with the first flight are described. A tabular overview of the observations is given. In addition, an example time series of a part of the emerging active region NOAA AR~11768 observed relatively close to disk centre is described and discussed in some detail. The observations cover the pores in the trailing polarity of the active region, as well as the polarity inversion line where flux emergence was ongoing and a small flare-like brightening occurred in the course of the time series. The pores are found to contain magnetic field strengths ranging up to 2500~G and, while large pores are clearly darker and cooler than the quiet Sun in all layers of the photosphere, the temperature and brightness of small pores approach or even exceed those of the quiet Sun in the upper photosphere.

We use seeing-free high spatial resolution Ca II H data obtained by the SUNRISE observatory to determine properties of slender fibrils in the lower solar chromosphere. In this work we use intensity images taken with the SUFI instrument in the Ca II H line during the second scientific flight of the SUNRISE observatory to identify and track elongated bright structures. After the identification, we analyze theses structures in order to extract their morphological properties. We identify 598 slender Ca II H fibrils (SCFs) with an average width of around 180 km, a length between 500 km and 4000 km, an average lifetime of ~400 s, and an average curvature of 0.002 arcsec^-1. The maximum lifetime of the SCFs within our time series of 57 minutes is ~2000 s. We discuss similarities and differences of the SCFs with other small-scale, chromospheric structures such as spicules of type I and II, or Ca II K fibrils.

Aims. A recent debate on the decade-long stability of the interstellar He flow vector, and in particular the flow longitude, has prompted us to check for any variability in the interstellar H flow vector as observed by the SWAN instrument on board SOHO. Methods. We used a simple model-independent method to determine the interstellar H flow longitude, based on the parallax effects induced on the Lyman-{\alpha} intensity measured by SWAN following the satellite motion around the Sun. Results. Our results show that the interstellar H flow vector longitude does not vary significantly from an average value of 252.9$^{\circ}$ $\pm$ 1.4$^{\circ}$ throughout the 20-year span of the SWAN dataset, further strengthening the arguments for the stability of the interstellar gas flow.

We study the dynamics and topology of an emerging magnetic flux concentration using high spatial resolution spectropolarimetric data acquired with the Imaging Magnetograph eXperiment on board the Sunrise balloon-borne solar observatory. We obtain the full vector magnetic field and the line-of-sight (LOS) velocity through inversions of the Fe I line at 525.02 nm with the SPINOR code. The derived vector magnetic field is used to trace magnetic field lines. Two magnetic flux concentrations with different polarity and LOS velocities are found to be connected by a group of arch-shaped magnetic field lines. The positive polarity footpoint is weaker (1100 G) and displays an upflow, while the negative polarity footpoint is stronger (2200 G) and shows a downflow. This configuration is naturally interpreted as a siphon flow along an arched magnetic flux tube.

We present a new method of Stokes inversion of spectropolarimetric data and evaluate it by taking the example of a SUNRISE/IMaX observation. An archive of synthetic Stokes profiles is obtained by the spectral synthesis of state-of-the-art magnetohydrodynamics (MHD) simulations and a realistic degradation to the level of the observed data. The definition of a merit function allows the archive to be searched for the synthetic Stokes profiles that match the observed profiles best. In contrast to traditional Stokes inversion codes, which solve the Unno-Rachkovsky equations for the polarized radiative transfer numerically and fit the Stokes profiles iteratively, the new technique provides the full set of atmospheric parameters. This gives us the ability to start an MHD simulation that takes the inversion result as initial condition. After a relaxation process of half an hour solar time we obtain physically consistent MHD data sets with a target similar to the observation. The new MHD simulation is used to repeat the method in a second iteration, which further improves the match between observation and simulation, resulting in a factor of 2.2 lower mean $\chi^2$ value. One advantage of the new technique is that it provides the physical parameters on a geometrical height scale. It constitutes a first step towards inversions giving results consistent with the MHD equations.

Active regions are the most prominent manifestations of solar magnetic fields; their generation and dissipation are fundamental problems in solar physics. Light bridges are commonly present during sunspot decay, but a comprehensive picture of their role in the removal of photospheric magnetic field is still missing. We study the three dimensional configuration of a sunspot and in particular its light bridge during one of the last stages of its decay. We present the magnetic and thermodynamical stratification inferred from full Stokes inversions of the photospheric Si I 10827 \AA\ and Ca I 10839 \AA\ lines obtained with the GREGOR Infrared Spectrograph of the GREGOR telescope at Observatorio del Teide, Tenerife, Spain. The analysis is complemented by a study of continuum images covering the disk passage of the active region, which are provided by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. The sunspot shows a light bridge with penumbral continuum intensity that separates the central umbra from a smaller umbra. We find that in this region the magnetic field lines form a canopy with lower magnetic field strength in the inner part. The photospheric light bridge is dominated by gas pressure (high-$\beta$), as opposed to the surrounding umbra where the magnetic pressure is higher. A convective flow is observed in the light bridge. This flow is able to bend the magnetic field lines and to produce field reversals. The field lines close above the light bridge and become as vertical and strong as in the surrounding umbra. We conclude that it develops because of two highly magnetized regions which come closer during the sunspot evolution.

Context. Investigations of the magnetism of the quiet Sun are hindered by extremely weak polarization signals in Fraunhofer spectral lines. Photon noise, straylight, and the systematically different sensitivity of the Zeeman effect to longitudinal and transversal magnetic fields result in controversial results in terms of the strength and angular distribution of the magnetic field vector. Aims. The information content of Stokes measurements close to the diffraction limit of the 1.5 m GREGOR telescope is analyzed. We took the effects of spatial straylight and photon noise into account. Methods. Highly sensitive full Stokes measurements of a quiet-Sun region at disk center in the deep photospheric Fe I lines in the 1.56 {\mu}m region were obtained with the infrared spectropolarimeter GRIS at the GREGOR telescope. Noise statistics and Stokes V asymmetries were analyzed and compared to a similar data set of the Hinode spectropolarimeter (SOT/SP). Simple diagnostics based directly on the shape and strength of the profiles were applied to the GRIS data. We made use of the magnetic line ratio technique, which was tested against MHD simulations. Results. About 80% of the GRIS spectra of a very quiet solar region show polarimetric signals above a 3{\sigma} level. Area and amplitude asymmetries agree well with small-scale surface dynamo MHD simulations. The magnetic line ratio analysis reveals ubiquitous magnetic regions in the ten to hundred Gauss range with some concentrations of kilo-Gauss fields. Conclusions. The GRIS spectropolarimetric data at a spatial resolution of 0.40" are so far unique in the combination of high spatial resolution scans and high magnetic field sensitivity. Nevertheless, the unavoidable effect of spatial straylight and the resulting dilution of the weak Stokes profiles means that inversion techniques still bear a high risk of misinterpretating the data.

We study the relation between mesogranular flows, convectively driven sinks and magnetic fields using high spatial resolution spectropolarimetric data acquired with the Imaging Magnetograph eXperiment on board Sunrise. We obtain the horizontal velocity flow fields of two quiet-Sun regions (31.2 $\times$ 31.2 Mm$^{2}$) via local correlation tracking. Mesogranular lanes and the central position of sinks are identified using Lagrange tracers. We find $6.7\times10^{-2}$ sinks per Mm$^{2}$ in the two observed regions. The sinks are located at the mesogranular vertices and turn out to be associated with (1) horizontal velocity flows converging to a central point and (2) long-lived downdrafts. The spatial distribution of magnetic fields in the quiet Sun is also examined. The strongest magnetic fields are preferentially located at sinks. We find that 40 \% of the pixels with longitudinal component of the magnetic field stronger than 500 G are located in the close neighborhood of sinks. In contrast, the small-scale magnetic loops detected by Mart\'{i}nez Gonz\'{a}lez et al. in the same two observed areas do not show any preferential distribution at mesogranular scales. The study of individual examples reveals that sinks can play an important role in the evolution of quiet-Sun magnetic features.

In June 2013, the two scientific instruments onboard the second Sunrise mission witnessed, in detail, a small-scale magnetic flux emergence event as part of the birth of an active region. The Imaging Magnetograph Experiment (IMaX) recorded two small (~5 arcsec) emerging flux patches in the polarized filtergrams of a photospheric Fe I spectral line. Meanwhile, the Sunrise Filter Imager (SuFI) captured the highly dynamic chromospheric response to the magnetic fields pushing their way through the lower solar atmosphere. The serendipitous capture of this event offers a closer look at the inner workings of active region emergence sites. In particular, it reveals in meticulous detail how the rising magnetic fields interact with the granulation as they push through the Sun's surface, dragging photospheric plasma in their upward travel. The plasma that is burdening the rising field slides along the field lines, creating fast downflowing channels at the footpoints. The weight of this material anchors this field to the surface at semi-regular spatial intervals, shaping it in an undulatory fashion. Finally, magnetic reconnection enables the field to release itself from its photospheric anchors, allowing it to continue its voyage up to higher layers. This process releases energy that lights up the arch-filament systems and heats the surrounding chromosphere.

Spectropolarimetric observations from Sunrise II/IMaX obtained in June 2013 are used for a statistical analysis to determine the physical properties of moving magnetic features (MMFs) observed near a pore. MMFs of the same and opposite polarity with respect to the pore are found to stream from its border at an average speed of 1.3 km s$^{-1}$ and 1.2 km s$^{-1}$ respectively, with mainly same-polarity MMFs found further away from the pore. MMFs of both polarities are found to harbor rather weak, inclined magnetic fields. Opposite-polarity MMFs are blue-shifted, while same-polarity MMFs do not show any preference for up- or downflows. Most of the MMFs are found to be of sub-arcsecond size and carry a mean flux of $\sim$ 1.2$\times 10^{17}$ Mx.

Ellerman Bombs are signatures of magnetic reconnection, which is an important physical process in the solar atmosphere. How and where they occur is a subject of debate. In this paper we analyse Sunrise/IMaX data together with 3D MHD simulations that aim to reproduce the exact scenario proposed for the formation of these features. Although the observed event seems to be more dynamic and violent than the simulated one, simulations clearly confirm the basic scenario for the production of EBs. The simulations also reveal the full complexity of the underlying process. The simulated observations show that the Fe I 525.02 nm line gives no information on the height where reconnection takes place. It can only give clues about the heating in the aftermath of the reconnection. The information on the magnetic field vector and velocity at this spatial resolution is, however, extremely valuable because it shows what numerical models miss and how they can be improved.

Arch filament systems occur in active sunspot groups, where a fibril structure connects areas of opposite magnetic polarity, in contrast to active region filaments that follow the polarity inversion line. We used the GREGOR Infrared Spectrograph (GRIS) to obtain the full Stokes vector in the spectral lines Si I 1082.7 nm, He I 1083.0 nm, and Ca I 1083.9 nm. We focus on the near-infrared calcium line to investigate the photospheric magnetic field and velocities, and use the line core intensities and velocities of the helium line to study the chromospheric plasma. The individual fibrils of the arch filament system connect the sunspot with patches of magnetic polarity opposite to that of the spot. These patches do not necessarily coincide with pores, where the magnetic field is strongest. Instead, areas are preferred not far from the polarity inversion line. These areas exhibit photospheric downflows of moderate velocity, but significantly higher downflows of up to 30 km/s in the chromospheric helium line. Our findings can be explained with new emerging flux where the matter flows downward along the fieldlines of rising flux tubes, in agreement with earlier results.

The various mechanisms of magneto-convective energy transport determines the structure of sunspots and active regions. We characterise the appearance of light bridges and other fine-structure details and elaborate on their magneto-convective nature. We present speckle-reconstructed images taken with the broad-band imager at the 1.5 m GREGOR telescope in the 486nm and 589nm bands. We estimate the spatial resolution from the noise characteristics of the image bursts and obtain 0.08" at 589nm. We describe structure details in individual best images as well as the temporal evolution of selected features. We find branched dark lanes extending along thin (~1") light bridges in sunspots at various heliocentric angles. In thick (~2") light bridges the branches are disconnected from the central lane and have a Y shape with a bright grain toward the umbra. The images reveal that light bridges exist on varying intensity levels and that their small-scale features evolve on timescales of minutes. Faint light bridges show dark lanes outlined by the surrounding bright features. Dark lanes are very common and are also found in the boundary of pores. They have a characteristic width of 0.1" or smaller. Intergranular dark lanes of that width are seen in active region granulation. While central dark lanes in thin light bridges are elevated and associated with an density increase above upflows, the dark lane branches correspond to locations of downflows and are depressed relative to the adjacent bright plasma. Thick light bridges with central dark lanes show no projection effect. They have a flat elevated plateau that falls off steeply at the umbral boundary. There, Y-shaped filaments form as they do in the inner penumbra. This indicates the presence of inclined magnetic fields, meaning that the umbral magnetic field is wrapped around the convective light bridge.

The fine-structure of magnetic field of a sunspot penumbra in the upper chromosphere is to be explored and compared to that in the photosphere. High spatial resolution spectropolarimetric observations were recorded with the 1.5-meter GREGOR telescope using the GREGOR Infrared Spectrograph (GRIS). The observed spectral domain includes the upper chromospheric He I triplet at 1083.0 nm and the photospheric Si I 1082.7 nm and Ca I 1083.3 nm spectral lines. The upper chromospheric magnetic field is obtained by inverting the He I triplet assuming a Milne-Eddington type model atmosphere. A height dependent inversion was applied to the Si I 1082.7 nm and Ca I 1083.3 nm lines to obtain the photospheric magnetic field. We find that the inclination of the magnetic field shows variations in the azimuthal direction both in the photosphere, but also in the upper chromosphere. The chromospheric variations remarkably well coincide with the variations in the inclination of the photospheric field and resemble the well-known spine and inter-spine structure in the photospheric layers of penumbrae. The typical peak-to-peak variations in the inclination of the magnetic field in the upper chromosphere is found to be 10-15 degree, i.e., roughly half the variation in the photosphere. In contrast, the magnetic field strength of the observed penumbra does not show variations on small spatial scales in the upper chromosphere. Thanks to the high spatial resolution observations possible with the GREGOR telescope at 1.08 microns, we find that the prominent small-scale fluctuations in the magnetic field inclination, which are a salient part of the property of sunspot penumbral photospheres, also persist in the chromosphere, although at somewhat reduced amplitudes. Such a complex magnetic configuration may facilitate penumbral chromospheric dynamic phenomena, such as penumbral micro-jets or transient bright dots.

Context. A significant part of the penumbral magnetic field returns below the surface in the very deep photosphere. For lines in the visible, a large portion of this return field can only be detected indirectly by studying its imprints on strongly asymmetric and three-lobed Stokes V profiles. Infrared lines probe a narrow layer in the very deep photosphere, providing the possibility of directly measuring the orientation of magnetic fields close to the solar surface. Aims. We study the topology of the penumbral magnetic field in the lower photosphere, focusing on regions where it returns below the surface. Methods. We analyzed 71 spectropolarimetric datasets from Hinode and from the GREGOR infrared spectrograph. We inferred the quality and polarimetric accuracy of the infrared data after applying several reduction steps. Techniques of spectral inversion and forward synthesis were used to test the detection algorithm. We compared the morphology and the fractional penumbral area covered by reversed-polarity and three-lobed Stokes V profiles for sunspots at disk center. We determined the amount of reversed-polarity and three-lobed Stokes V profiles in visible and infrared data of sunspots at various heliocentric angles. From the results, we computed center-to-limb variation curves, which were interpreted in the context of existing penumbral models. Results. Observations in visible and near-infrared spectral lines yield a significant difference in the penumbral area covered by magnetic fields of opposite polarity. In the infrared, the number of reversed-polarity Stokes V profiles is smaller by a factor of two than in the visible. For three-lobed Stokes V profiles the numbers differ by up to an order of magnitude.

Some models for the topology of the magnetic field in sunspot penumbrae predict the existence of field-free or dynamically weak-field regions in the deep Photosphere. To confirm or rule out the existence of weak-field regions in the deepest photospheric layers of the penumbra. The magnetic field at $\log\tau_5=0$ is investigated by means of inversions of spectropolarimetric data of two different sunspots located very close to disk center with a spatial resolution of approximately 0.4-0.45 arcsec. The data have been recorded using the GRIS instrument attached to the 1.5-meters GREGOR solar telescope at El Teide observatory. It includes three Fe I lines around 1565 nm, whose sensitivity to the magnetic field peaks at half a pressure-scale-height deeper than the sensitivity of the widely used Fe I spectral line pair at 630 nm. Prior to the inversion, the data is corrected for the effects of scattered light using a deconvolution method with several point spread functions. At $\log\tau_5=0$ we find no evidence for the existence of regions with dynamically weak ($B<500$ Gauss) magnetic fields in sunspot penumbrae. This result is much more reliable than previous investigations done with Fe I lines at 630 nm. Moreover, the result is independent of the number of nodes employed in the inversion, and also independent of the point spread function used to deconvolve the data, and does not depend on the amount of straylight (i.e. wide-angle scattered light) considered.

Compressible magnetohydrodynamic (MHD) turbulence is ubiquitous in astrophysical phenomena ranging from the intergalactic to the stellar scales. In studying them, numerical simulations are nearly inescapable, due to the large degree of nonlinearity involved. However the dynamical ranges of these phenomena are much larger than what is computationally accessible. In large eddy simulations (LES), the resulting limited resolution effects are addressed explicitly by introducing to the equations of motion additional terms associated with the unresolved, subgrid-scale (SGS) dynamics. This renders the system unclosed. We derive a set of nonlinear structural closures for the ideal MHD LES equations with particular emphasis on the effects of compressibility. The closures are based on a gradient expansion of the finite-resolution operator (W.K. Yeo CUP 1993, ed. Galperin & Orszag) and require no assumptions about the nature of the flow or magnetic field. Thus the scope of their applicability ranges from the sub- to the hyper-sonic and -Alfvenic regimes. The closures support spectral energy cascades both up and down-scale, as well as direct transfer between kinetic and magnetic resolved and unresolved energy budgets. They implicitly take into account the local geometry, and in particular the anisotropy, of the flow. Their properties are $\textit{a priori}$ validated in an accompanying article (Grete et al. Phys. Plasmas, 2016) against alternative closures available in the literature with respect to a wide range of simulation data of homogeneous and isotropic turbulence.

Even though compressible plasma turbulence is encountered in many astrophysical phenomena, its effect is often not well understood. Furthermore, direct numerical simulations are typically not able to reach the extreme parameters of these processes. For this reason, large-eddy simulations (LES), which only simulate large and intermediate scales directly, are employed. The smallest, unresolved scales and the interactions between small and large scales are introduced by means of a subgrid-scale (SGS) model. We propose and verify a new set of nonlinear SGS closures for future application as an SGS model in LES of compressible magnetohydrodynamics (MHD). We use 15 simulations (without explicit SGS model) of forced, isotropic, homogeneous turbulence with varying sonic Mach number $\mathrm{M_s} = 0.2$ to $20$ as reference data for the most extensive \textit{a priori} tests performed so far in literature. In these tests we explicitly filter the reference data and compare the performance of the new closures against the most widely tested closures. These include eddy-viscosity and scale-similarity type closures with different normalizations. Performance indicators are correlations with the turbulent energy and cross-helicty flux, the average SGS dissipation, the topological structure and the ability to reproduce the correct magnitude and direction of the SGS vectors. We find that only the new nonlinear closures exhibit consistently high correlations (median value \textgreater$0.8$) with the data over the entire parameter space and outperform the other closures in all tests. Moreover, we show that these results are independent of resolution and chosen filter scale. Additionally, the new closures are effectively coefficient-free with a deviation of less than $20\%$.

The Visible Tunable Filter (VTF) is a narrowband tunable filter system for imaging spectroscopy and spectropolarimetry based. The instrument will be one of the first-light instruments of the Daniel K. Inouye Solar Telescope that is currently under construction on Maui (Hawaii). The VTF is being developed by the Kiepenheuer Institut fuer Sonnenphysik in Freiburg as a German contribution to the DKIST. We perform end-to-end simulations of spectropolarimetric observations with the VTF to verify the science requirements of the instrument. The instrument is simulated with two Etalons, and with a single Etalon. The clear aperture of the Etalons is 250 mm, corresponding to a field of view with a diameter of 60 arcsec in the sky (42,000 km on the Sun). To model the large-scale figure errors we employ low-order Zernike polynomials (power and spherical aberration) with amplitudes of 2.5 nm RMS. We use an ideal polarization modulator with equal modulation coefficients of 3-1/2 for the polarization modulation We synthesize Stokes profiles of two iron lines (630.15 nm and 630.25 nm) and for the 854.2 nm line of calcium, for a range of magnetic field values and for several inclination angles. We estimated the photon noise on the basis of the DKIST and VTF transmission values, the atmospheric transmission and the spectral flux from the Sun. For the Fe 630.25 nm line, we obtain a sensitivity of 20 G for the longitudinal component and for 150 G for the transverse component, in agreement with the science requirements for the VTF.

Aims. Combining high-resolution and synoptic observations aims to provide a comprehensive description of flux emergence at photospheric level and of the growth process that eventually leads to a mature active region. Methods. Small active region NOAA 12118 was observed on 2014 July 18 with the 1.5-meter GREGOR solar telescope on 2014 July 18. High-resolution time-series of blue continuum and G-band images acquired in the blue imaging channel (BIC) of the GREGOR Fabry-P\'erot Interferometer (GFPI) were complemented by LOS magnetograms and continuum images obtained with the HMI onboard the SDO. Horizontal proper motions and horizontal plasma velocities were computed with local correlation tracking (LCT) and the differential affine velocity estimator, respectively. Morphological image processing was employed to measure the photometric/magnetic area, magnetic flux, and the separation profile of the EFR during its evolution. Results. The computed growth rates for photometric area, magnetic area, and magnetic flux are about twice as high as the respective decay rates. The space-time diagram using HMI magnetograms of five days traces a leaf-like structure, which is determined by the initial separation of the two polarities, a rapid expansion phase, a time when the spread stalls, and a period when the region slowly shrinks again. The separation rate of 0.26 km\s is highest in the initial stage, and it decreases when the separation comes to a halt. Horizontal plasma velocities computed at four evolutionary stages indicate a changing pattern of inflows. In LCT maps we find persistent flow patterns such as outward motions in the outer part of the two major pores, a diverging feature near the trailing pore marking the site of upwelling plasma and flux emergence, and low velocities in the interior of pores. We detected many elongated rapidly expanding granules between the two major polarities.

A small flare ribbon above a sunspot umbra in active region 12205 was observed on November 7, 2014, at 12:00 UT in the blue imaging channel of the 1.5 m GREGOR telescope, using a 1 A Ca II H interference filter. Context observations from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO), the Solar Optical Telescope (SOT) onboard Hinode, and the Interface Region Imaging Spectrograph (IRIS) show that this ribbon is part of a larger one that extends through the neighboring positive polarities and also participates in several other flares within the active region. We reconstructed a time series of 140 seconds of Ca II H images by means of the multiframe blind deconvolution method, which resulted in spatial and temporal resolutions of 0.1 arcsec and 1 s. Light curves and horizontal velocities of small-scale bright knots in the observed flare ribbon were measured. Some knots are stationary, but three move along the ribbon with speeds of 7-11 km/s. Two of them move in the opposite direction and exhibit highly correlated intensity changes, which provides evidence of a slipping reconnection at small spatial scales.

Context. Comets are primitive objects, remnants of the volatile-rich planetesimals from which the solar system condensed. Knowing their structure and composition is thus crucial for the understanding of our origins. After the successful landing of Philae on the nucleus of 67P/Churyumov-Gerasimenko in November 2014, for the first time, the Rosetta mission provided the opportunity to measure the low frequency electrical properties of a cometary mantle with the permittivity probe SESAME-PP (Surface Electric Sounding and Acoustic Monitoring Experiment - Permittivity Probe). Aims. In this paper, we conduct an in-depth analysis of the data from active measurements collected by SESAME-PP at Abydos, which is the final landing site of Philae, to constrain the porosity and, to a lesser extent, the composition of the surface material down to a depth of about 1 meter. Methods. SESAME-PP observations on the surface are then analyzed by comparison with data acquired during the descent toward the nucleus and with numerical simulations that explore different possible attitudes and environments of Philae at Abydos using a method called the Capacity-Influence Matrix Method. Results. Reasonably assuming that the two receiving electrode channels have not drifted with respect to each other during the ten-year journey of the Rosetta probe to the comet, we constrain the dielectric constant of the first meter below the surface at Abydos to be > 2.45 p/m 0.20, which is consistent with a porosity < 50% if the dust phase is analogous to carbonaceous chondrites and < 75% in the case of less primitive ordinary chondrites. This indicates that the near surface of the nucleus of 67P/Churyumov-Gerasimenko nucleus is more compacted than its interior and suggests that it could consist of a sintered dust-ice layer.

The gas in galaxy clusters is heated by shock compression through accretion (outer shocks) and mergers (inner shocks). These processes additionally produce turbulence. To analyse the relation between the thermal and turbulent energies of the gas under the influence of non-adiabatic processes, we performed numerical simulations of cosmic structure formation in a box of 152 Mpc comoving size with radiative cooling, UV background, and a subgrid scale model for numerically unresolved turbulence. By smoothing the gas velocities with an adaptive Kalman filter, we are able to estimate bulk flows toward cluster cores. This enables us to infer the velocity dispersion associated with the turbulent fluctuation relative to the bulk flow. For halos with masses above $10^{13}\,M_\odot$, we find that the turbulent velocity dispersions averaged over the warm-hot intergalactic medium (WHIM) and the intracluster medium (ICM) are approximately given by powers of the mean gas temperatures with exponents around 0.5, corresponding to a roughly linear relation between turbulent and thermal energies and transonic Mach numbers. However, turbulence is only weakly correlated with the halo mass. Since the power-law relation is stiffer for the WHIM, the turbulent Mach number tends to increase with the mean temperature of the WHIM. This can be attributed to enhanced turbulence production relative to dissipation in particularly hot and turbulent clusters.

Accurate measurements of magnetic and velocity fields in and around solar active regions are key to unlocking the mysteries of the formation and the decay of sunspots. High spatial resolution image and spectral sequences with a high cadence obtained with the GREGOR solar telescope give us an opportunity to scrutinize 3-D flow fields with local correlation tracking and imaging spectroscopy. We present GREGOR early science data acquired in 2014 July - August with the GREGOR Fabry-P\'erot Interferometer and the Blue Imaging Channel. Time-series of blue continuum (? 450.6 nm) images of the small active region NOAA 12118 were restored with the speckle masking technique to derive horizontal proper motions and to track the evolution of morphological changes. In addition, high-resolution observations are discussed in the context of synoptic data from the Solar Dynamics Observatory.

The new generation of solar instruments provides better spectral, spatial, and temporal resolution for a better understanding of the physical processes that take place on the Sun. Multiple-component profiles are more commonly observed with these instruments. Particularly, the He I 10830 \r{A} triplet presents such peculiar spectral profiles, which give information on the velocity and magnetic fine structure of the upper chromosphere. The purpose of this investigation is to describe a technique to efficiently fit the two blended components of the He I 10830 \r{A} triplet, which are commonly observed when two atmospheric components are located within the same resolution element. The observations used in this study were taken on 2015 April 17 with the very fast spectroscopic mode of the GREGOR Infrared Spectrograph (GRIS) attached to the 1.5-meter GREGOR solar telescope, located at the Observatorio del Teide, Tenerife, Spain. We apply a double-Lorentzian fitting technique using Levenberg-Marquardt least-squares minimization. This technique is very simple and much faster than inversion codes. Line-of-sight Doppler velocities can be inferred for a whole map of pixels within just a few minutes. Our results show sub- and supersonic downflow velocities of up to 32 km/s for the fast component in the vicinity of footpoints of filamentary structures. The slow component presents velocities close to rest.

The magnetic field of a sunspot guides magnetohydrodynamic waves toward higher atmospheric layers. In the upper photosphere and lower chromosphere, wave modes with periods longer than the acoustic cut-off period become evanescent. The cut-off period essentially changes due to the atmospheric properties, e.g., increases for larger zenith inclinations of the magnetic field. In this work, we aim at introducing a novel technique of reconstructing the magnetic field inclination on the basis of the dominating wave periods in the sunspot chromosphere and upper photosphere. On 2013 August 21st, we observed an isolated, circular sunspot (NOAA11823) for 58 min in a purely spectroscopic multi-wavelength mode with the Interferometric Bidimensional Spectro-polarimeter (IBIS) at the Dunn Solar Telescope. By means of a wavelet power analysis, we retrieved the dominating wave periods and reconstructed the zenith inclinations in the chromosphere and upper photosphere. The results are in good agreement with the lower photospheric HMI magnetograms. The sunspot's magnetic field in the chromosphere inclines from almost vertical (0$^{\circ}$) in the umbra to around 60$^{\circ}$ in the outer penumbra. With increasing altitude in the sunspot atmosphere, the magnetic field of the penumbra becomes less inclined. We conclude that the reconstruction of the magnetic field topology on the basis of sunspot oscillations yields consistent and conclusive results. The technique opens up a new possibility to infer the magnetic field inclination in the solar chromosphere.

We report on the dynamical interaction of quiet-Sun magnetic fields and granular convection in the solar photosphere as seen by \textsc{Sunrise}. We use high spatial resolution (0\farcs 15--0\farcs 18) and temporal cadence (33 s) spectropolarimetric Imaging Magnetograph eXperiment data, together with simultaneous CN and Ca\,\textsc{ii}\,H filtergrams from \textsc{Sunrise} Filter Imager. We apply the SIR inversion code to the polarimetric data in order to infer the line of sight velocity and vector magnetic field in the photosphere. The analysis reveals bundles of individual flux tubes evolving as a single entity during the entire 23 minute data set. The group shares a common canopy in the upper photospheric layers, while the individual tubes continually intensify, fragment and merge in the same way that chains of bright points in photometric observations have been reported to do. The evolution of the tube cores are driven by the local granular convection flows. They intensify when they are "compressed" by surrounding granules and split when they are "squeezed" between two moving granules. The resulting fragments are usually later regrouped in intergranular lanes by the granular flows. The continual intensification, fragmentation and coalescence of flux results in magnetic field oscillations of the global entity. From the observations we conclude that the magnetic field oscillations first reported by \citet{2011ApJ...730L..37M} correspond to the forcing by granular motions and not to characteristic oscillatory modes of thin flux tubes.

We investigate the influence of different analytical parameterizations and fit functions for the local star formation rate in AMR simulations of an isolated disk galaxy with the Nyx code. Such parameterizations express the star formation efficiency as function of the local turbulent Mach number and virial parameter. By employing the method of adaptively refined large eddy simulations, we are able to evaluate these physical parameters from the numerically unresolved turbulent energy associated with the grid scale. We consider both single and multi free-fall variants of star formation laws proposed by Padoan & Nordlund, Hennebelle & Chabrier, and Krumholz & McKee, summarised and tested recently with numerical simulations by Federrath & Klessen. We find that the global star formation rate and the relation between the local star formation rate and the gas column density is reproduced in agreement with observational constraints by all multi free-fall models of star formation. Some models with obsolete calibration or a single free-fall time scale, however, result in an overly clumpy disk that does not resemble the structure of observed spirals.

We model driven, compressible, isothermal, turbulence with Mach numbers ranging from the subsonic ($\mathcal{M} \approx 0.65$) to the highly supersonic regime ($\mathcal{M}\approx 16 $). The forcing scheme consists both solenoidal (transverse) and compressive (longitudinal) modes in equal parts. We find a relation $\sigma_{s}^2 = \mathrm{b}\log{(1+\mathrm{b}^2\mathcal{M}^2)}$ between the Mach number and the standard deviation of the logarithmic density with $\mathrm{b} = 0.457 \pm 0.007$. The density spectra follow $\mathcal{D}(k,\,\mathcal{M}) \propto k^{\zeta(\mathcal{M})}$ with scaling exponents depending on the Mach number. We find $\zeta(\mathcal{M}) = \alpha \mathcal{M}^{\beta}$ with a coefficient $\alpha$ that varies slightly with resolution, whereas $\beta$ changes systematically. We extrapolate to the limit of infinite resolution and find $\alpha = -1.91 \pm 0.01,\, \beta =-0.30\pm 0.03$. The dependence of the scaling exponent on the Mach number implies a fractal dimension $D=2+0.96 \mathcal{M}^{-0.30}$. We determine how the scaling parameters depend on the wavenumber and find that the density spectra are slightly curved. This curvature gets more pronounced with increasing Mach number. We propose a physically motivated fitting formula $\mathcal{D}(k) = \mathcal{D}_0 k^{\zeta k^{\eta}}$ by using simple scaling arguments. The fit reproduces the spectral behaviour down to scales $k\approx 80$. The density spectrum follows a single power-law $\eta = -0.005 \pm 0.01$ in the low Mach number regime and the strongest curvature $\eta = -0.04 \pm 0.02$ for the highest Mach number. These values of $\eta$ represent a lower limit, as the curvature increases with resolution.

The solar photospheric abundance of oxygen is still a matter of debate. For about ten years some determinations have favoured a low oxygen abundance which is at variance with the value inferred by helioseismology. Among the oxygen abundance indicators, the forbidden line at 630nm has often been considered the most reliable even though it is blended with a NiI line. In Papers I and Paper II of this series we reported a discrepancy in the oxygen abundance derived from the 630nm and the subordinate [OI] line at 636nm in dwarf stars, including the Sun. Here we analyse several, in part new, solar observations of the the centre-to-limb variation of the spectral region including the blend at 630nm in order to separate the individual contributions of oxygen and nickel. We analyse intensity spectra observed at different limb angles in comparison with line formation computations performed on a CO5BOLD 3D hydrodynamical simulation of the solar atmosphere. The oxygen abundances obtained from the forbidden line at different limb angles are inconsistent if the commonly adopted nickel abundance of 6.25 is assumed in our local thermodynamic equilibrium computations. With a slightly lower nickel abundance, A(Ni)~6.1, we obtain consistent fits indicating an oxygen abundance of A(O)=8.73+/-0.05. At this value the discrepancy with the subordinate oxygen line remains. The derived value of the oxygen abundance supports the notion of a rather low oxygen abundance in the solar hotosphere. However, it is disconcerting that the forbidden oxygen lines at 630 and 636nm give noticeably different results, and that the nickel abundance derived here from the 630nm blend is lower than expected from other nickel lines.

Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct numerical simulations computationally intractable even for the simplest treatment -- magnetohydrodynamics (MHD). To overcome this problem one can use subgrid-scale (SGS) closures -- models for the influence of unresolved, subgrid-scales on the resolved ones. In this work we propose and validate a set of constant coefficient closures for the resolved, compressible, ideal MHD equations. The subgrid-scale energies are modeled by Smagorinsky-like equilibrium closures. The turbulent stresses and the electromotive force (EMF) are described by expressions that are nonlinear in terms of large scale velocity and magnetic field gradients. To verify the closures we conduct a priori tests over 137 simulation snapshots from two different codes with varying ratios of thermal to magnetic pressure ($\beta_\mathrm{p} = 0.25, 1, 2.5, 5, 25$) and sonic Mach numbers ($M_s = 2, 2.5, 4$). Furthermore, we make a comparison to traditional, phenomenological eddy-viscosity and $\alpha-\beta-\gamma$ closures. We find only mediocre performance of the kinetic eddy-viscosity and $\alpha-\beta-\gamma$ closures, and that the magnetic eddy-viscosity closure is poorly correlated with the simulation data. Moreover, three of five coefficients of the traditional closures exhibit a significant spread in values. In contrast, our new closures demonstrate consistently high correlation and constant coefficient values over time and and over the wide range of parameters tested. Important aspects in compressible MHD turbulence such as the bi-directional energy cascade, turbulent magnetic pressure and proper alignment of the EMF are well described by our new closures.

We investigate a new scheme for astronomical spectrograph calibration using the laser frequency comb at the Solar Vacuum Tower Telescope on Tenerife. Our concept is based upon a single-mode fiber channel, that simultaneously feeds the spectrograph with comb light and sunlight. This yields nearly perfect spatial mode matching between the two sources. In combination with the absolute calibration provided by the frequency comb, this method enables extremely robust and accurate spectroscopic measurements. The performance of this scheme is compared to a sequence of alternating comb and sunlight, and to absorption lines from Earth's atmosphere. We also show how the method can be used for radial-velocity detection by measuring the well-explored 5-minute oscillations averaged over the full solar disk. Our method is currently restricted to solar spectroscopy, but with further evolving fiber-injection techniques it could become an option even for faint astronomical targets.

The {\em Solar Orbiter} is the next solar physics mission of the European Space Agency, ESA, in collaboration with NASA, with a launch planned in 2018. The spacecraft is designed to approach the Sun to within 0.28\,AU at perihelion of a highly eccentric orbit. The proximity with the Sun will also allow its observation at uniformly high resolution at EUV and visible wavelengths. Such observations are central for learning more about the magnetic coupling of the solar atmosphere. At a later phase in the mission the spacecraft will leave the ecliptic and study the enigmatic poles of the Sun from a heliographic latitude of up to 33$^\circ$.

We compare two different codes for simulations of cosmological structure formation to investigate the sensitivity of hydrodynamical instabilities to numerics, in particular, the hydro solver and the application of adaptive mesh refinement (AMR). As a simple test problem, we consider an initially spherical gas cloud in a wind, which is an idealized model for the merger of a subcluster or galaxy with a big cluster. Based on an entropy criterion, we calculate the mass stripping from the subcluster as a function of time. Moreover, the turbulent velocity field is analyzed with a multi-scale filtering technique. We find remarkable differences between the commonly used PPM solver with directional splitting in the Enzo code and an unsplit variant of PPM in the Nyx code, which demonstrates that different codes can converge to systematically different solutions even when using uniform grids. For the test case of an unbound cloud, AMR simulations reproduce uniform-grid results for the mass stripping quite well, although the flow realizations can differ substantially. If the cloud is bound by a static gravitational potential, however, we find strong sensitivity to spurious fluctuations which are induced at the cutoff radius of the potential and amplified by the bow shock. This gives rise to substantial deviations between uniform-grid and AMR runs performed with Enzo, while the mass stripping in Nyx simulations of the subcluster is nearly independent of numerical resolution and AMR. Although many factors related to numerics are involved, our study indicates that unsplit solvers with advanced flux limiters help to reduce grid effects and to keep numerical noise under control, which is important for hydrodynamical instabilities and turbulent flows.

The evolution of the physical parameters of magnetic bright points (MBPs) located in the quiet Sun (mainly in the interwork) during their lifetime is studied. First we concentrate on the detailed description of the magnetic field evolution of three MBPs. This reveals that individual features follow different, generally complex, and rather dynamic scenarios of evolution. Next we apply statistical methods on roughly 200 observed MBP evolutionary tracks. MBPs are found to be formed by the strengthening of an equipartition field patch, which initially exhibits a moderate downflow. During the evolution, strong downdrafts with an average velocity of 2.4 km/s set in. These flows, taken together with the concurrent strengthening of the field, suggest that we are witnessing the occurrence of convective collapses in these features, although only 30% of them reach kG field strengths. This fraction might turn out to be larger when the new 4 m class solar telescopes are operational as observations of MBPs with current state of the art instrumentation could still be suffering from resolution limitations. Finally, when the bright point disappears (although the magnetic field often continues to exist) the magnetic field strength has dropped to the equipartition level and is generally somewhat weaker than at the beginning of the MBP's evolution. Noteworthy is that in about 10% of the cases we observe in the vicinity of the downflows small-scale strong (exceeding 2 km/s) intergranular upflows related spatially and temporally to these downflows.

We present a subgrid-scale model for the Multi-phase Interstellar medium, Star formation, and Turbulence (MIST) and explore its behaviour in high-resolution large-eddy simulations of isolated disc galaxies. MIST follows the evolution of a clumpy cold and a diffuse warm component of the gas within a volume element which exchange mass and energy via various cooling, heating and mixing processes. The star formation rate is dynamically computed from the state of the gas in the cold phase. An important feature of MIST is the treatment of unresolved turbulence in the two phases and its interaction with star formation and feedback by supernovae. This makes MIST a particularly suitable model for the interstellar medium in galaxy simulations. We carried out a suite of simulations varying fundamental parameters of our feedback implementation. Several observational properties of galactic star formation are reproduced in our simulations, such as an average star formation efficiency ~1%, a typical velocity dispersion around ~10 km/s in star-forming regions, and an almost linear relationship between the column densities of star formation and dense molecular gas.

Low-metallicity star formation poses a central problem of cosmology, as it determines the characteristic mass scale and distribution for the first and second generations of stars forming in our Universe. Here, we present a comprehensive investigation assessing the relative impact of metals and magnetic fields, which may both be present during low-metallicity star formation. We show that the presence of magnetic fields generated via the small-scale dynamo stabilises the protostellar disc and provides some degree of support against fragmentation. In the absence of magnetic fields, the fragmentation timescale in our model decreases by a factor of ~10 at the transition from Z=0 to Z>0, with subsequently only a weak dependence on metallicity. Similarly, the accretion timescale of the cluster is set by the large-scale dynamics rather than the local thermodynamics. In the presence of magnetic fields, the primordial disc can become completely stable, therefore forming only one central fragment. At Z>0, the number of fragments is somewhat reduced in the presence of magnetic fields, though the shape of the mass spectrum is not strongly affected in the limits of the statistical uncertainties. The fragmentation timescale, however, increases by roughly a factor of 3 in the presence of magnetic fields. Indeed, our results indicate comparable fragmentation timescales in primordial runs without magnetic fields and Z>0 runs with magnetic fields.

We update two kinds of results obtained with the SWAN instrument on board SOHO. First, we use H cell data recorded in 2001 and derive the H flow direction in the same way we performed the study at solar minimum. We compare with the Helium flow and doing so we correct for the coordinate system change between the Ulysses and SOHO mission. The deflection plane we obtain is compatible with the previous result within error bars, confirming the predominant role of the interstellar magnetic field. Secondly, we extend the derivation of solar wind ionization temporal evolution as a function of heliolatitude. The pattern for the present solar minimum is strikingly different from the previous minimum, with a much wider slow solar wind equatorial belt which persists until at least 2008. Comparing with synoptic LASCO/C2 electron densities we infer from a preliminary study that the acceleration of the high speed solar wind occurs at a higher altitude during this minimum, a change expansion models should be able to explain.

Isolated flux tubes are considered to be fundamental magnetic building blocks of the solar photosphere. Their formation is usually attributed to the concentration of magnetic field to kG strengths by the convective collapse mechanism. However, the small size of the magnetic elements in quiet-Sun areas has prevented this scenario from being studied in fully resolved structures. Here we report on the formation and subsequent evolution of one such photospheric magnetic flux tube, observed in the quiet Sun with unprecedented spatial resolution (0\farcs 15 - 0\farcs 18) and high temporal cadence (33 s). The observations were acquired by the Imaging Magnetograph Experiment (IMaX) aboard the \textsc{Sunrise} balloon-borne solar observatory. The equipartition field strength magnetic element is the result of the merging of several same polarity magnetic flux patches, including a footpoint of a previously emerged loop. The magnetic structure is then further intensified to kG field strengths by convective collapse. The fine structure found within the flux concentration reveals that the scenario is more complex than can be described by a thin flux tube model with bright points and downflow plumes being established near the edges of the kG magnetic feature. We also observe a daisy-like alignment of surrounding granules and a long-lived inflow towards the magnetic feature. After a subsequent weakening process, the field is again intensified to kG strengths. The area of the magnetic feature is seen to change in anti-phase with the field strength, while the brightness of the bright points and the speed of the downflows varies in phase. We also find a relation between the brightness of the bright point and the presence of upflows within it.

In this review, the methodology of large eddy simulations (LES) is introduced and applications in astrophysics are discussed. As theoretical framework, the scale decomposition of the dynamical equations for neutral fluids by means of spatial filtering is explained. For cosmological applications, the filtered equations in comoving coordinates are also presented. To obtain a closed set of equations that can be evolved in LES, several subgrid scale models for the interactions between numerically resolved and unresolved scales are discussed, in particular the subgrid scale turbulence energy equation model. It is then shown how model coefficients can be calculated, either by dynamical procedures or, a priori, from high-resolution data. For astrophysical applications, adaptive mesh refinement is often indispensable. It is shown that the subgrid scale turbulence energy model allows for a particularly elegant and physically well motivated way of preserving momentum and energy conservation in AMR simulations. Moreover, the notion of shear-improved models for inhomogeneous and non-stationary turbulence is introduced. Finally, applications of LES to turbulent combustion in thermonuclear supernovae, star formation and feedback in galaxies, and cosmological structure formation are reviewed.

We investigate turbulence generated by cosmological structure formation by means of large eddy simulations using adaptive mesh refinement. In contrast to the widely used implicit large eddy simulations, which resolve a limited range of length scales and treat the effect of turbulent velocity fluctuations below the grid scale solely by numerical dissipation, we apply a subgrid-scale model for the numerically unresolved fraction of the turbulence energy. For simulations with adaptive mesh refinement, we utilize a new methodology that allows us to adjust the scale-dependent energy variables in such a way that the sum of resolved and unresolved energies is globally conserved. We test our approach in simulations of randomly forced turbulence, a gravitationally bound cloud in a wind, and the Santa Barbara cluster. To treat inhomogeneous turbulence, we introduce an adaptive Kalman filtering technique that separates turbulent velocity fluctuations on resolved length scales from the non-turbulent bulk flow. From the magnitude of the fluctuating component and the subgrid-scale turbulence energy, a total turbulent velocity dispersion of several 100 km/s is obtained for the Santa Barbara cluster, while the low-density gas outside the accretion shocks is nearly devoid of turbulence. The energy flux through the turbulent cascade and the dissipation rate predicted by the subgrid-scale model correspond to dynamical time scales around 5 Gyr, independent of numerical resolution.

Observations of quasars at $\rm z> 6$ report the existence of a billion solar mass black holes. Comprehending their formation in such a short time scale is a matter of ongoing research. One of the most promising scenarios to assemble supermassive black holes is a monolithic collapse of protogalactic gas clouds in atomic cooling halos with $\rm T_{vir} \geq 10^{4} K$. In this article, we study the amplification and impact of magnetic fields during the formation of seed black holes in massive primordial halos. We perform high resolution cosmological magnetohydrodynamics simulations for four distinct halos and follow their collapse for a few free-fall times until the simulations reach a peak density of $\rm 7 \times 10^{-10} g/cm^{3}$. Our findings show that irrespective of the initial seed field, the magnetic field strength reaches a saturated state in the presence of strong accretion shocks. Under such conditions, the growth time becomes very short and amplification occurs rapidly within a small fraction of the free-fall time. We find that the presence of such strong magnetic fields provides additional support against gravity and helps in suppressing fragmentation. Massive clumps of a few hundred solar masses are formed at the end of our simulations and high accretion rates of $\rm 1 M_{\odot}/yr$ are observed. We expect that in the presence of such accretion rates, the clumps will grow to form supermassive stars of $\rm \sim 10^{5} M_{\odot}$. Overall, the role of the magnetic fields seems supportive for the formation of massive black holes.

We analyzed the size, intensity, and magnetic field strength of sunspot umbrae to compare the present cycle 24 with the previous one. We used data of the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory and selected all sunspots between May 2010 and October 2012, using one image per day. We created two subsets of this data with a manual tracking algorithm, both without duplication. One is containing each sunspot (910 umbrae within 488 spots) and was used to analyze the distribution of umbral areas, selected with an automated thresholding method. The other one contains 205 fully evolved sunspots. We find nonlinear relations between umbral minimum intensity and size and between maximum magnetic field strength and size. The field strength scales linear with the intensity and the umbral size scales roughly linear with the total magnetic flux, while the size and field strength level off with stronger flux. When separated in hemisphere and averaged temporally, the southern umbrae show a temporal increase in size and the northern umbrae stay constant. There is no temporal variation in the umbral mean intensity detectable. The probability density function of the umbral area in the ascending phase of the current solar cycle is similar to that of the last solar cycle. From our investigation of umbral area, magnetic field, magnetic flux and umbral intensity of the sunspots of the rising phase of cycle 24, we do not find a significant difference to the previous cycle, and hence no indication for a long-term decline of solar activity.

We present the first high-resolution solar images in the Mg II k 2796 \AA{} line. The images, taken through a 4.8 \AA{} broad interference filter, were obtained during the second science flight of SUNRISE in June 2013 by the SuFI instrument. The Mg II k images display structures that look qualitatively very similar to images taken in the core of Ca II H. The Mg II images exhibit reversed granulation (or shock waves) in the internetwork regions of the quiet Sun, at intensity contrasts that are similar to those found in Ca II H. Very prominent in Mg II are bright points, both in the quiet Sun and in plage regions, particularly near disk center. These are much brighter than at other wavelengths sampled at similar resolution. Furthermore, Mg II k images also show fibril structures associated with plage regions. Again, the fibrils are similar to those seen in Ca II H images, but tend to be more pronounced, particularly in weak plage.

Black holes of a billion solar masses are observed in the infant universe a few hundred million years after the Big Bang. The direct collapse of protogalactic gas clouds in primordial halos with $\rm T_{vir} \geq 10^{4} K$ provides the most promising way to assemble massive black holes. In this study, we aim to determine the characteristic mass scale of seed black holes and the time evolution of the accretion rates resulting from the direct collapse model. We explore the formation of supermassive black holes via cosmological large eddy simulations (LES) by employing sink particles and following their evolution for twenty thousand years after the formation of the first sink. As the resulting protostars were shown to have cool atmospheres in the presence of strong accretion, we assume here that UV feedback is negligible during this calculation. We confirm this result in a comparison run without sinks. Our findings show that black hole seeds with characteristic mass of $\rm 10^{5} M_{\odot}$ are formed in the presence of strong Lyman Werner flux which leads to an isothermal collapse. The characteristic mass is a about two times higher in LES compared to the implicit large eddy simulations (ILES). The accretion rates increase with time and reach a maximum value of 10 $\rm M_{\odot}/yr$ after $\rm 10^{4}$ years. Our results show that the direct collapse model is clearly feasible as it provides the expected mass of the seed black holes.

Population III stars forming in the infant universe at z=30 heralded the end of the cosmic dark ages. They are presumed to be assembled in so-called minihaloes with virial temperatures of a few thousand K where collapse is triggered by molecular hydrogen cooling. A central question concerns their final masses, and whether fragmentation occurs during their formation. While studies employing Lagrangian codes suggest fragmentation via a self-gravitating disk, recent high resolution simulations indicated that disk formation is suppressed. Here we report the first high-resolution large-eddy simulations performed with the Eulerian grid-based code Enzo following the evolution beyond the formation of the first peak, to investigate the accretion of the central massive clump and potential fragmentation. For a total of 3 halos, we see that a disk forms around the first clump. The central clump reaches $\sim10$ solar masses after 40 years, while subsequent accretion is expected at a rate of $10^{-2}$ solar masses per year. In one of these halos, additional clumps form as a result of fragmentation which proceeds at larger scales. We note that subgrid-scale turbulence yields relevant contributions to the stability of the protostellar disks. We conclude that the first protostar may reach masses up to $\rm 40-100 M_{\odot}$, which are only limited by the effect of radiative feedback.

Supermassive black holes with up to a $\rm 10^{9}~M_{\odot}$ dwell in the centers of present-day galaxies, and their presence has been confirmed at z $\geq$ 6. Their formation at such early epochs is still an enigma. Different pathways have been suggested to assemble supermassive black holes in the first billion years after the Big Bang. Direct collapse has emerged as a highly plausible scenario to form black holes as it provides seed masses of $\rm 10^{5}-10^{6}~M_{\odot}$. Gravitational collapse in atomic cooling haloes with virial temperatures T$_{vir} \geq 10^{4}$~K may lead to the formation of massive seed black holes in the presence of an intense background UV flux. Turbulence plays a central role in regulating accretion and transporting angular momentum. We present here the highest resolution cosmological large-eddy simulations to date which track the evolution of high-density regions on scales of $0.25$~AU beyond the formation of the first peak, and study the impact of subgrid-scale turbulence. The peak density reached in these simulations is $\rm 1.2 \times 10^{-8}~g~cm^{-3}$. Our findings show that while fragmentation occasionally occurs, it does not prevent the growth of a central massive object resulting from turbulent accretion and occasional mergers. The central object reaches $\rm \sim 1000~M_{\odot}$ within $4$ free-fall times, and we expect further growth up to $\rm 10^{6}~M_{\odot}$ through accretion in about 1 million years. The direct collapse model thus provides a viable pathway of forming high-mass black holes at early cosmic times.

While present standard model of cosmology yields no clear prediction for the initial magnetic field strength, efficient dynamo action may compensate for initially weak seed fields via rapid amplification. In particular, the small-scale dynamo is expected to exponentially amplify any weak magnetic field in the presence of turbulence. We explore whether this scenario is viable using cosmological magneto-hydrodynamics simulations modeling the formation of the first galaxies, which are expected to form in so-called atomic cooling halos with virial temperatures $\rm T_{vir} \geq 10^{4}$ K. As previous calculations have shown that a high Jeans resolution is needed to resolve turbulent structures and dynamo effects, our calculations employ resolutions of up to 128 cells per Jeans length. The presence of the dynamo can be clearly confirmed for resolutions of at least 64 cells per Jeans length, while saturation occurs at approximate equipartition with turbulent energy. As a result of the large Reynolds numbers in primordial galaxies, we expect saturation to occur at early stages, implying magnetic field strengths of \sim0.1 $\mu$G at densities of 10^4 cm^{-3}.

Comparisons of the integrated thermal pressure support of gas against its gravitational potential energy lead to critical mass scales for gravitational instability such as the Jeans and the Bonnor-Ebert masses, which play an important role in analysis of many physical systems, including the heuristics of numerical simulations. In a strict theoretical sense, however, neither the Jeans nor the Bonnor-Ebert mass are meaningful when applied locally to substructure in a self-gravitating turbulent medium. For this reason, we investigate the local support by thermal pressure, turbulence, and magnetic fields against gravitational compression through an approach that is independent of these concepts. At the centre of our approach is the dynamical equation for the divergence of the velocity field. We carry out a statistical analysis of the source terms of the local compression rate (the negative time derivative of the divergence) for simulations of forced self-gravitating turbulence in periodic boxes with zero, weak, and moderately strong mean magnetic fields (measured by the averages of the magnetic and thermal pressures). We also consider the amplification of the magnetic field energy by shear and by compression. Thereby, we are able to demonstrate that the support against gravity is dominated by thermal pressure fluctuations, although magnetic pressure also yields a significant contribution. The net effect of turbulence in the highly supersonic regime, however, is to enhance compression rather than supporting overdense gas even if the vorticity is very high. This is incommensurate with the support of the highly dynamical substructures in magneto-turbulent fluids being determined by local virial equilibria of volume energies without surface stresses.

In a previous work, we reported on the discovery of supersonic magnetic upflows on granular cells in data from the {\sc Sunrise}/IMaX instrument. In the present work we investigate the physical origin of these events employing data of the same instrument but with higher spectral sampling. By means of the inversion of Stokes profiles we are able to recover the physical parameters (temperature, magnetic field, line-of-sight velocity, etc) present in the solar photosphere at the time of these events. The inversion is performed in a Monte-Carlo-like fashion, that is, repeating it many times with different initializations and retaining only the best result. We find that many of the events are characterized by a reversal in the polarity of the magnetic field along the vertical direction in the photosphere, accompanied by an enhancement in the temperature and by supersonic line-of-sight velocities. In about half of the studied events, large blue-shifted and red-shifted line-of-sight velocities coexist above/below each other. These features can be explained in terms of magnetic reconnection, where the energy stored in the magnetic field is released in the form of kinetic and thermal energy when magnetic field lines of opposite polarities coalesce. However, the agreement with magnetic reconnection is not perfect and therefore, other possible physical mechanisms might also play a role.

(Abridged) The small-scale dynamo may play a substantial role in magnetizing the Universe under a large range of conditions, including subsonic turbulence at low Mach numbers, highly supersonic turbulence at high Mach numbers and a large range of magnetic Prandtl numbers Pm, i.e. the ratio of kinetic viscosity to magnetic resistivity. Low Mach numbers may in particular lead to the well-known, incompressible Kolmogorov turbulence, while for high Mach numbers, we are in the highly compressible regime, thus close to Burgers turbulence. In this study, we explore whether in this large range of conditions, a universal behavior can be expected. Our starting point are previous investigations in the kinematic regime. Here, analytic studies based on the Kazantsev model have shown that the behavior of the dynamo depends significantly on Pm and the type of turbulence, and numerical simulations indicate a strong dependence of the growth rate on the Mach number of the flow. Once the magnetic field saturates on the current amplification scale, backreactions occur and the growth is shifted to the next-larger scale. We employ a Fokker-Planck model to calculate the magnetic field amplification during the non-linear regime, and find a resulting power-law growth that depends on the type of turbulence invoked. For Kolmogorov turbulence, we confirm previous results suggesting a linear growth of magnetic energy. For more general turbulent spectra, where the turbulent velocity v_t scales with the characteristic length scale as u_\ell\propto \ell^{\vartheta}, we find that the magnetic energy grows as (t/T_{ed})^{2\vartheta/(1-\vartheta)}, with t the time-coordinate and T_{ed} the eddy-turnover time on the forcing scale of turbulence. For Burgers turbulence, \vartheta=1/2, a quadratic rather than linear growth may thus be expected, and a larger timescale until saturation is reached.

This paper is a numerical study of the condensation of the warm neutral medium (WNM) into cold neutral medium (CNM) structures under the effect of turbulence and thermal instability. Using low resolution simulations we explored the impact of the WNM initial density and properties of the turbulence (stirring in Fourier with a varying mix of solenoidal and compressive modes) on the cold gas formation. Two sets of initial conditions which match the observations were selected to produce high resolution simulations (1024^3) allowing to study in details the properties of the produced dense structures. For typical values of the density, pressure and velocity dispersion of the WNM in the solar neighborhood, the turbulent motions of the HI can not provoque the phase transition from WNM to CNM, whatever their amplitude and their distribution in solenoidal and compressive modes. On the other hand we show that a quasi-isothermal increase in WNM density of a factor of 2 to 4 is enough to induce the phase transition, leading to the transition of about 40 percent of the gas to the cold phase within 1 Myr. Given the observed properties of the HI in the local ISM, the WNM and individual CNM structures in the local ISM are sub or transsonic and their dynamics are tightly interwoven. The velocity field bears the evidence of subsonic turbulence with a 2D power spectrum following the Kolmogorov law as P(k) \sim k^{-8/3} while the density is highly contrasted with a singificantly shallower power spectrum, reminiscent of what is observed in the cold ISM. Supra-thermal line width observed for CNM might be the result of relative velocity between cold structures. Finally, the cold structures denser than 5 cm^{-3} reproduce well the laws M \sim L^{2.25-2.28} and sigma(v) \sim 0.5-0.8L^{1/3} generally observed in molecular clouds.

Atomic cooling haloes with virial temperatures $\rm T_{vir} \geq 10^{4}$ K are the most plausible sites for the formation of the first galaxies and the first intermediate mass black holes. It is therefore important to assess whether one can obtain robust results concerning their main properties from numerical simulations. A major uncertainty is the presence of turbulence, which is barely resolved in cosmological simulations. We explore the latter both by pursuing high-resolution simulations with up to 64 cells per Jeans length and by incorporating a subgrid-scale turbulence model to account for turbulent pressure and viscosity on unresolved scales. We find that the main physical quantities in the halo, in particular the density, temperature and energy density profile, are approximately converged. However, the morphologies in the central 500 AU change significantly with increasing resolution and appear considerably more turbulent. In a systematic comparison of three different haloes, we further found that the turbulence subgrid-scale model gives rise to more compact central structures, and decreases the amount of vorticity. Such compact morphologies may in particular favor the accretion onto the central object.

We explore the amplification of magnetic fields in the high-redshift Universe. For this purpose, we perform high-resolution cosmological simulations following the formation of primordial halos with \sim10^7 M_solar, revealing the presence of turbulent structures and complex morphologies at resolutions of at least 32 cells per Jeans length. Employing a turbulence subgrid-scale model, we quantify the amount of unresolved turbulence and show that the resulting turbulent viscosity has a significant impact on the gas morphology, suppressing the formation of low-mass clumps. We further demonstrate that such turbulence implies the efficient amplification of magnetic fields via the small-scale dynamo. We discuss the properties of the dynamo in the kinematic and non-linear regime, and explore the resulting magnetic field amplification during primordial star formation. We show that field strengths of \sim10^{-5} G can be expected at number densities of \sim5 cm^{-3}.

In this review, we look back upon the literature, which had the GREGOR solar telescope project as its subject including science cases, telescope subsystems, and post-focus instruments. The articles date back to the year 2000, when the initial concepts for a new solar telescope on Tenerife were first presented at scientific meetings. This comprehensive bibliography contains literature until the year 2012, i.e., the final stages of commissioning and science verification. Taking stock of the various publications in peer-reviewed journals and conference proceedings also provides the "historical" context for the reference articles in this special issue of Astronomische Nachrichten/Astronomical Notes.

We analyze the spectral asymmetry of Stokes V (circularly polarized) profiles of an individual network patch in the quiet Sun observed by Sunrise/IMaX. At a spatial resolution of 0.15"-0.18", the network elements contain substructure which is revealed by the spatial distribution of Stokes V asymmetries. The area asymmetry between the red and blue lobes of Stokes V increases from nearly zero at the core of the structure to values close to unity at its edges (one-lobed profiles). Such a distribution of the area asymmetry is consistent with magnetic fields expanding with height, i.e., an expanding magnetic canopy (which is required to fulfill pressure balance and flux conservation in the solar atmosphere). Inversion of the Stokes I and V profiles of the patch confirms this picture, revealing a decreasing field strength and increasing height of the canopy base from the core to the periphery of the network patch. However, the non-roundish shape of the structure and the presence of negative area and amplitude asymmetries reveal that the scenario is more complex than a canonical flux tube expanding with height surrounded by downflows.

We investigate the interplay between the ionization radiation from massive stars and the turbulence inside the surrounding molecular gas thanks to 3D numerical simulations. We used the 3D hydrodynamical code HERACLES to model an initial turbulent medium that is ionized and heated by an ionizing source. Three different simulations are performed with different mean Mach numbers (1, 2 and 4). A non-equilibrium model for the ionization and the associated thermal processes was used. This revealed to be crucial when turbulent ram pressure is of the same order as the ionized-gas pressure. The density structures initiated by the turbulence cause local curvatures of the dense shell formed by the ionization compression. When the curvature of the shell is sufficient, the shell collapse on itself to form a pillar while a smaller curvature leads to the formation of dense clumps that are accelerated with the shell and therefore remain in the shell during the simulation. When the turbulent ram pressure of the cold gas is sufficient to balance the ionized-gas pressure, some dense-gas bubbles have enough kinetic energy to penetrate inside the ionized medium, forming cometary globules. This suggests a direct relation in the observations between the presence of globules and the relative importance of the turbulence compared to the ionized-gas pressure. The probability density functions present a double peak structure when the turbulence is low relative to the ionized-gas pressure. This could be used in observations as an indication of the turbulence inside molecular clouds.

We study the variation in the magnetic field strength and the umbral intensity of sunspots during the declining phase of the solar cycle no.23 and in the beginning of cycle no.24. We analyze a sample of 183 sunspots observed from 1999 until 2011 with the Tenerife Infrared Polarimeter at the German Vacuum Tower Telescope. The magnetic field strength is derived from the Zeeman splitting of the Stokes-V signal in one near-infrared spectral line, either Fe I 1564.8 nm, Fe I 1089.6 nm, or Si I 1082.7 nm. This avoids the effects of the unpolarized stray light from the field-free quiet Sun surroundings. The minimum umbral continuum intensity and umbral area are also measured. We find that there is a systematic trend for sunspots in the late stage of the solar cycle no.23 to be weaker, i.e., to have a smaller maximum magnetic field strength than those at the start of the cycle. The decrease in the field strength with time of about 94 G/yr is well beyond the statistical fluctuations that would be expected because of the larger number of sunspots close to cycle maximum (14 G/yr). In the same time interval, the continuum intensity of the umbra increases with a rate of 1.3 (+- 0.4)% of Ic/yr, while the umbral area does not show any trend above the statistical variance. Sunspots in the new cycle no.24 show higher field strengths and lower continuum intensities than those at the end of cycle no.23, interrupting the trend. Sunspots have an intrinsically weaker field strength and brighter umbrae at the late stages of solar cycles compared to their initial stages, without any significant change in their area. The abrupt increase in field strength in sunspots of the new cycle suggests that the cyclic variations are dominating over any long-term trend that continues across cycles. We find a slight decrease in field strength and an increase in intensity as a long-term trend across the cycles.

2011 was a successful year for the GREGOR project. The telescope was finally completed in May with the installation of the 1.5-meter primary mirror. The installation of the first-light focal plane instruments was completed by the end of the year. At the same time, the preparations for the installation of the high-order adaptive optics were finished, its integration to the telescope is scheduled for early 2012. This paper describes the telescope and its instrumentation in their present first light configuration, and provides a brief overview of the science goals of GREGOR.

In this contribution we describe some recent observations of high-speed magnetized flows in the quiet Sun granulation. These observations were carried out with the Imaging Magnetograph eXperiment (IMaX) onboard the stratospheric balloon {\sc Sunrise}, and possess an unprecedented spatial resolution and temporal cadence. These flows were identified as highly shifted circular polarization (Stokes $V$) signals. We estimate the LOS velocity responsible for these shifts to be larger than 6 km s$^{-1}$, and therefore we refer to them as {\it supersonic magnetic flows}. The average lifetime of the detected events is 81.3 s and they occupy an average area of about 23\,000 km$^2$. Most of the events occur within granular cells and correspond therefore to upflows. However some others occur in intergranular lanes or bear no clear relation to the convective velocity pattern. We analyze a number of representative examples and discuss them in terms of magnetic loops, reconnection events, and convective collapse.

We report on magnetic field emergences covering significant areas of exploding granules. The balloon-borne mission SUNRISE provided high spatial and temporal resolution images of the solar photosphere. Continuum images, longitudinal and transverse magnetic field maps and Dopplergrams obtained by IMaX onboard SUNRISE are analyzed by Local Correlation Traking (LCT), divergence calculation and time slices, Stokes inversions and numerical simulations are also employed. We characterize two mesogranular-scale exploding granules where $\sim$ 10$^{18}$ Mx of magnetic flux emerges. The emergence of weak unipolar longitudinal fields ($\sim$100 G) start with a single visible magnetic polarity, occupying their respective granules' top and following the granular splitting. After a while, mixed polarities start appearing, concentrated in downflow lanes. The events last around 20 min. LCT analyses confirm mesogranular scale expansion, displaying a similar pattern for all the physical properties, and divergence centers match between all of them. We found a similar behaviour with the emergence events in a numerical MHD simulation. Granule expansion velocities are around 1 \kms while magnetic patches expand at 0.65 \kms. One of the analyzed events evidences the emergence of a loop-like structure. Advection of the emerging magnetic flux features is dominated by convective motion resulting from the exploding granule due to the magnetic field frozen in the granular plasma. Intensification of the magnetic field occurs in the intergranular lanes, probably because of being directed by the downflowing plasma.

We present a systematic study of the influence of different forcing types on the statistical properties of supersonic, isothermal turbulence in both the Lagrangian and Eulerian frameworks. We analyse a series of high-resolution, hydrodynamical grid simulations with Lagrangian tracer particles and examine the effects of solenoidal (divergence-free) and compressive (curl-free) forcing on structure functions, their scaling exponents, and the probability density functions of the gas density and velocity increments. Compressively driven simulations show a significantly larger density contrast, a more intermittent behaviour, and larger fractal dimension of the most dissipative structures at the same root mean square Mach number. We show that the absolute values of Lagrangian and Eulerian structure functions of all orders in the integral range are only a function of the root mean square Mach number, but independent of the forcing. With the assumption of a Gaussian distribution for the probability density function of the velocity increments on large scales, we derive a model that describes this behaviour.

We report on the photospheric evolution of an intermediate-scale (~4 Mm footpoint separation) magnetic bipole, from emergence to decay, observed in the quiet Sun at high spatial 0".3 and temporal (33 s) resolution. The observations were acquired by the IMaX imaging magnetograph during the first science flight of the Sunrise balloon-borne solar observatory. The bipole flux content is 6 x 10^17 Mx, representing a structure bridging the gap between granular scale bipoles and the smaller ephemeral regions. Footpoints separate at a speed of 3.5 km s-1 and reach a maximum distance of 4.5 Mm before the field dissolves. The evolution of the bipole is revealed to be very dynamic: we found a proper motion of the bipole axis and detected a change of the azimuth angle of 90{\deg} in 300 seconds. The overall morphology and behaviour are in agreement with previous analyses of bipolar structures emerging at granular scale, but we also found several similarities with emerging flux structures at larger scale. The flux growth rate is 2.6 x 15 Mx s-1, while the mean decay rate is one order of magnitude smaller. We describe in some detail the decay phase of the bipole footpoints which includes break up into smaller structures, interaction with pre-existing fields leading to cancellation but appears to be dominated by an as-yet unidentified diffusive process that removes most of the flux with an exponential flux decay curve. The diffusion constant (8 x 10^2 km^2 s-1) associated with this decay is similar to the values used to describe the large scale diffusion in flux transport models.

We present a semi-analytic model for the interstellar medium that considers local processes and structures of turbulent star-forming gas. A volume element of the interstellar medium is described as a multiphase system, comprising a cold and a warm gas phase in effective (thermal plus turbulent) pressure equilibrium, and a stellar component. Since turbulence has a critical impact on the shape of the gaseous phases, on the production of molecular hydrogen and on the formation of stars, the consistent treatment of turbulence energy -- the kinetic energy of unresolved motions -- is an important new feature of our model. Besides turbulence production by supernovae and by the cooling instability, we also take into account the forcing by large scale motions. We formulate a set of ordinary differential equations, which statistically describes star formation and the exchange between the different budgets of mass and energy in a region of the interstellar medium with given mean density, size, metallicity and external turbulence forcing. By exploring the behaviour of the solutions, we find equilibrium states, in which the star formation efficiencies are consistent with observations. Kennicutt-Schmidt-like relations naturally arise from the equilibrium solutions, while conventional star formation models in numerical simulations impose such relations with observed efficiency parameters as phenomenological calibrations. Beyond the semi-analytic approach, a potential application is a complete subgrid scale model of the unresolved multi-phase structure, star formation and turbulence in simulations of galaxies or in cosmological simulations. The formulation presented in this article combines various models focusing on particular processes and yet can be adopted to specific applications, depending on the range of resolved length scales.

FEARLESS (Fluid mEchanics with Adaptively Refined Large Eddy SimulationS) is a numerical scheme for modelling subgrid-scale turbulence in cosmological adaptive mesh refinement simulations. In this contribution, the main features of this tool will be outlined. We discuss the application of this method to cosmological simulations of the large-scale structure. The simulations show that the production of turbulence has a different redshift dependence in the intra-cluster medium and the warm-hot intergalactic medium, caused by the distinct stirring mechanisms (mergers and shock interactions) acting in them. Some properties of the non-thermal pressure support in the two baryon phases are also described.

The transit of Venus in 2004 offered the rare possibility to remotely sense a well-known planetary atmosphere using ground-based observations for absorption spectroscopy. Transmission spectra of Venus' atmosphere were obtained in the near infrared using the Vacuum Tower Telescope (VTT) in Tenerife. Since the instrument was designed to measure the very bright photosphere of the Sun, extracting Venus' atmosphere was challenging. CO_2 absorption lines could be identified in the upper Venus atmosphere. Moreover, the relative abundance of the three most abundant CO_2 isotopologues could be determined. The observations resolved Venus' limb, showing Doppler-shifted absorption lines that are probably caused by high-altitude winds. This paper illustrates the ability of ground-based measurements to examine atmospheric constituents of a terrestrial planet atmosphere which might be applied in future to terrestrial extrasolar planets.

We present observational evidence for oscillations of magnetic flux density in the quiet areas of the Sun. The majority of magnetic fields on the solar surface have strengths of the order of or lower than the equipartition field (300-500 G). This results in a myriad of magnetic fields whose evolution is largely determined by the turbulent plasma motions. When granules evolve they squash the magnetic field lines together or pull them apart. Here we report on the periodic deformation of the shapes of features in circular polarization observed at high resolution with Sunrise. In particular, we note that the area of patches with constant magnetic flux oscillates with time, which implies that the apparent magnetic field intensity oscillates in antiphase. The periods associated to this oscillatory pattern is compatible with the granular life-time and change abruptly, which suggests that these oscillations might not correspond to characteristic oscillatory modes of magnetic structures, but to the forcing by granular motions. In one particular case, we find three patches around the same granule oscillating in phase, which means that the spatial coherence of these oscillations can reach 1600 km. Interestingly, the same kind of oscillatory phenomenon is found also in the upper photosphere.

The injection and evolution of turbulence in the intergalactic medium is studied by means of mesh-based hydrodynamical simulations, including a subgrid scale (SGS) model for small-scale unresolved turbulence. The simulations show that the production of turbulence has a different redshift dependence in the intracluster medium (ICM) and the warm-hot intergalactic medium (WHIM). We show that turbulence in the ICM is produced chiefly by merger-induced shear flows, whereas the production in the WHIM is dominated by shock interactions. Secondly, the effect of dynamical pressure support on the gravitational contraction has been studied. This turbulent support is stronger in the WHIM gas at baryon overdensities 1 < delta < 100, and less relevant for the ICM. Although the relative mass fraction of the gas with large vorticity is considerable (52% in the ICM), we find that for only about 10% in mass this is dynamically relevant, namely not associated to an equally large thermal pressure support. According to this result, a significant non-thermal pressure support counteracting the gravitational contraction is a localised characteristic in the cosmic flow, rather than a widespread feature.

The final results of processing the data from the balloon-born experiment ATIC-2 (Antarctica, 2002-2003) for the energy spectra of protons and He, C, O, Ne, Mg, Si, and Fe nuclei, the spectrum of all particles, and the mean logarithm of atomic weight of primary cosmic rays as a function of energy are presented. The final results are based on improvement of the methods used earlier, in particular, considerably increased resolution of the charge spectrum. The preliminary conclusions on the significant difference in the spectra of protons and helium nuclei (the proton spectrum is steeper) and the non-power character of the spectra of protons and heavier nuclei (flattening of carbon spectrum at energies above 10 TeV) are confirmed. A complex structure of the energy dependence of the mean logarithm of atomic weight is found.

The relation of the solar surface magnetic field with mesogranular cells is studied using high spatial (~ 100 km) and temporal (~ 30 sec) resolution data obtained with the IMaX instrument aboard SUNRISE. First, mesogranular cells are identified using Lagrange tracers (corks) based on horizontal velocity fields obtained through Local Correlation Tracking. After ~ 20 min of integration, the tracers delineate a sharp mesogranular network with lanes of width below about 280 km. The preferential location of magnetic elements in mesogranular cells is tested quantitatively. Roughly 85% of pixels with magnetic field higher than 100 G are located in the near neighborhood of mesogranular lanes. Magnetic flux is therefore concentrated in mesogranular lanes rather than intergranular ones. Secondly, magnetic field extrapolations are performed to obtain field lines anchored in the observed flux elements. This analysis, therefore, is independent of the horizontal flows determined in the first part. A probability density function (PDF) is calculated for the distribution of distances between the footpoints of individual magnetic field lines. The PDF has an exponential shape at scales between 1 and 10 Mm, with a constant characteristic decay distance, indicating the absence of preferred convection scales in the mesogranular range. Our results support the view that mesogranulation is not an intrinsic convective scale (in the sense that it is not a primary energy-injection scale of solar convection), but also give quantitative confirmation that, nevertheless, the magnetic elements are preferentially found along mesogranular lanes.

Context: The thermal instability is one of the dynamical agents for turbulence in the diffuse interstellar medium, where both, turbulence and thermal instability interact in a highly non-linear manner. Aims: We study basic properties of turbulence in thermally bistable gas for variable simulation parameters. The resulting cold gas fractions can be applied as parameterisation in simulations on galactic scales. Methods: Turbulent flow is induced on large scales by means of compressive stochastic forcing in a periodic box. The compressible Euler equations with constant UV heating and a parameterised cooling function are solved on uniform grids. We investigate several values of the mean density of the gas and different magnitudes of the forcing. For comparison with other numerical studies, solenoidal forcing is applied as well. Results: After a transient phase, we observe that a state of statistically stationary turbulence is approached. Compressive forcing generally produces a two-phase medium, with a decreasing fraction of cold gas for increasing forcing strength. This behaviour can be explained on the basis of turbulent mixing. We also find indications for power-law tails of probability density functions of the gas density. Solenoidal forcing, on the other hand, appears to prevent the evolution into a two-phase-medium for certain parameter regions. Conclusions: The dynamics of thermally bistable turbulence shows a substantial sensitivity on the initial state and the forcing properties.

We formulate and implement the Euler equations with SGS dynamics and provide numerical tests of an SGS turbulence energy model that predicts the turbulent pressure of unresolved velocity fluctuations and the rate of dissipation for highly compressible turbulence. We test closures for the turbulence energy cascade by filtering data from high-resolution simulations of forced isothermal and adiabatic turbulence. Optimal properties and an excellent correlation are found for a linear combination of the eddy-viscosity closure that is employed in LES of weakly compressible turbulence and a term that is non-linear in the Jacobian matrix of the velocity. Using this mixed closure, the SGS turbulence energy model is validated in LES of turbulence with stochastic forcing. It is found that the SGS model satisfies several important requirements: 1. The mean SGS turbulence energy follows a power law for varying grid scale. 2. The root mean square (RMS) Mach number of the unresolved velocity fluctuations is proportional to the RMS Mach number of the resolved turbulence, independent of the forcing. 3. The rate of dissipation and the turbulence energy flux are constant. Moreover, we discuss difficulties with direct estimates of the turbulent pressure and the dissipation rate on the basis of resolved flow quantities that have recently been proposed. In combination with the energy injection by stellar feedback and other unresolved processes, the proposed SGS model is applicable to a variety of problems in computational astrophysics. Computing the SGS turbulence energy, the treatment of star formation and stellar feedback in galaxy simulations can be improved. Further, we expect that the turbulent pressure on the grid scale affects the stability of gas against gravitational collapse.

We study the energy flux carried by acoustic waves excited by convective motions at sub-photospheric levels. The analysis of high-resolution spectropolarimetric data taken with IMaX/Sunrise provides a total energy flux of ~ 6400--7700 Wm$^{-2}$ at a height of ~ 250 km in the 5.2-10 mHz range, i.e. at least twice the largest energy flux found in previous works. Our estimate lies within a factor of 2 of the energy flux needed to balance radiative losses from the chromosphere according to Anderson & Athay (1989) and revives interest in acoustic waves for transporting energy to the chromosphere. The acoustic flux is mainly found in the intergranular lanes but also in small rapidly-evolving granules and at the bright borders, forming dark dots and lanes of splitting granules.

We characterize the observational properties of the convectively driven vortex flows recently discovered on the quiet Sun, using magnetograms, Dopplergrams and images obtained with the 1-m balloon-borne Sunrise telescope. By visual inspection of time series, we find some 3.1e-3 vortices/(Mm^2 min), which is a factor of 1.7 larger than previous estimates. The mean duration of the individual events turns out to be 7.9 min, with a standard deviation of 3.2 min. In addition, we find several events appearing at the same locations along the duration of the time series (31.6 min). Such recurrent vortices show up in the proper motion flow field map averaged over the time series. The typical vertical vorticities are <= 6e-3 1/sec, which corresponds to a period of rotation of some 35 min. The vortices show a preferred counterclockwise sense of rotation, which we conjecture may have to do with the preferred vorticity impinged by the solar differential rotation.

Bright points (BPs) are manifestations of small magnetic elements in the solar photosphere. Their brightness contrast not only gives insight into the thermal state of the photosphere (and chromosphere) in magnetic elements, but also plays an important role in modulating the solar total and spectral irradiance. Here we report on simultaneous high-resolution imaging and spectropolarimetric observations of BPs using Sunrise balloon-borne observatory data of the quiet Sun at disk center. BP contrasts have been measured between 214 nm and 525 nm, including the first measurements at wavelengths below 388 nm. The histograms of the BP peak brightness show a clear trend toward broader contrast distributions and higher mean contrasts at shorter wavelengths. At 214 nm we observe a peak brightness of up to five times the mean quiet-Sun value, the highest BP contrast so far observed. All BPs are associated with a magnetic signal, although in a number of cases it is surprisingly weak. Most of the BPs show only weak downflows, the mean value being 240 m/s, but some display strong down- or upflows reaching a few km/s.

Using the IMaX instrument on-board the Sunrise stratospheric balloon-telescope we have detected extremely shifted polarization signals around the Fe I 5250.217 {\AA} spectral line within granules in the solar photosphere. We interpret the velocities associated with these events as corresponding to supersonic and magnetic upflows. In addition, they are also related to the appearance of opposite polarities and highly inclined magnetic fields. This suggests that they are produced by the reconnection of emerging magnetic loops through granular upflows. The events occupy an average area of 0.046 arcsec$^2$ and last for about 80 seconds, with larger events having longer lifetimes. These supersonic events occur at a rate of $1.3\times10^{-5}$ occurrences per second per arcsec$^{2}$.

Until today, the small size of magnetic elements in quiet Sun areas has required the application of indirect methods, such as the line-ratio technique or multi-component inversions, to infer their physical properties. A consistent match to the observed Stokes profiles could only be obtained by introducing a magnetic filling factor that specifies the fraction of the observed pixel filled with magnetic field. Here, we investigate the properties of a small magnetic patch in the quiet Sun observed with the IMaX magnetograph on board the balloon-borne telescope Sunrise with unprecedented spatial resolution and low instrumental stray light. We apply an inversion technique based on the numerical solution of the radiative transfer equation to retrieve the temperature stratification and the field strength in the magnetic patch. The observations can be well reproduced with a one-component, fully magnetized atmosphere with a field strength exceeding 1 kG and a significantly enhanced temperature in the mid- to upper photosphere with respect to its surroundings, consistent with semi-empirical flux tube models for plage regions. We therefore conclude that, within the framework of a simple atmospheric model, the IMaX measurements resolve the observed quiet-Sun flux tube.

Solar oscillations are expected to be excited by turbulent flows in the intergranular lanes near the solar surface. Time series recorded by the IMaX instrument aboard the {\sc Sunrise} observatory reveal solar oscillations at high resolution, which allow studying the properties of oscillations with short wavelengths. We analyze two times series with synchronous recordings of Doppler velocity and continuum intensity images with durations of 32\thinspace min and 23\thinspace min, resp., recorded close to the disk center of the Sun to study the propagation and excitation of solar acoustic oscillations. In the Doppler velocity data, both the standing acoustic waves and the short-lived, high-degree running waves are visible. The standing waves are visible as temporary enhancements of the amplitudes of the large-scale velocity field due to the stochastic superposition of the acoustic waves. We focus on the high-degree small-scale waves by suitable filtering in the Fourier domain. Investigating the propagation and excitation of $f$- and $p_1$-modes with wave numbers $k > 1.4$\thinspace 1/Mm we find that also exploding granules contribute to the excitation of solar $p$-modes in addition to the contribution of intergranular lanes.

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 investigate the fine structure of magnetic fields in the atmosphere of the quiet Sun. We use photospheric magnetic field measurements from {\sc Sunrise}/IMaX with unprecedented spatial resolution to extrapolate the photospheric magnetic field into higher layers of the solar atmosphere with the help of potential and force-free extrapolation techniques. We find that most magnetic loops which reach into the chromosphere or higher have one foot point in relatively strong magnetic field regions in the photosphere. $91%$ of the magnetic energy in the mid chromosphere (at a height of 1 Mm) is in field lines, whose stronger foot point has a strength of more than 300 G, i.e. above the equipartition field strength with convection. The loops reaching into the chromosphere and corona are also found to be asymmetric in the sense that the weaker foot point has a strength $B < 300$ G and is located in the internetwork. Such loops are expected to be strongly dynamic and have short lifetimes, as dictated by the properties of the internetwork fields.

This paper describes the wave-front correction system developed for the Sunrise balloon telescope, and provides information about its in-flight performance. For the correction of low-order aberrations, a Correlating Wave-Front Sensor (CWS) was used. It consisted of a six-element Shack-Hartmann wave-front sensor (WFS), a fast tip-tilt mirror for the compensation of image motion, and an active telescope secondary mirror for focus correction. The CWS delivered a stabilized image with a precision of 0.04 arcsec (rms), whenever the coarse pointing was better than 90 arcsec peak-to-peak. The automatic focus adjustment maintained a focus stability of 0.01 waves in the focal plane of the CWS. During the 5.5 day flight, good image quality and stability was achieved during 33 hours, containing 45 sequences that lasted between 10 and 45 minutes.

The first science flight of the balloon-borne \Sunrise telescope took place in June 2009 from ESRANGE (near Kiruna/Sweden) to Somerset Island in northern Canada. We describe the scientific aims and mission concept of the project and give an overview and a description of the various hardware components: the 1-m main telescope with its postfocus science instruments (the UV filter imager SuFI and the imaging vector magnetograph IMaX) and support instruments (image stabilizing and light distribution system ISLiD and correlating wavefront sensor CWS), the optomechanical support structure and the instrument mounting concept, the gondola structure and the power, pointing, and telemetry systems, and the general electronics architecture. We also explain the optimization of the structural and thermal design of the complete payload. The preparations for the science flight are described, including AIV and ground calibration of the instruments. The course of events during the science flight is outlined, up to the recovery activities. Finally, the in-flight performance of the instrumentation is briefly summarized.

We present high-resolution images of the Sun in the near ultraviolet spectral range between 214 nm and 397 nm as obtained from the first science flight of the 1-m Sunrise balloon-borne solar telescope. The quiet-Sun rms intensity contrasts found in this wavelength range are among the highest values ever obtained for quiet-Sun solar surface structures - up to 32.8% at a wavelength of 214 nm. We compare with theoretical intensity contrasts obtained from numerical magneto-hydrodynamic simulations. For 388 nm and 312 nm the observations agree well with the numerical simulations whereas at shorter wavelengths discrepancies between observed and simulated contrasts remain.

We describe the design of the Sunrise Filter Imager (SuFI) and the Image Stabilization and Light Distribution (ISLiD) unit onboard the Sunrise balloon borne solar observatory. This contribution provides the necessary information which is relevant to understand the instruments working principles, the relevant technical data, and the necessary information about calibration issues directly related to the science data.

The Imaging Magnetograph eXperiment (IMaX) is a spectropolarimeter built by four institutions in Spain that flew on board the Sunrise balloon-borne telesocope in June 2009 for almost six days over the Arctic Circle. As a polarimeter IMaX uses fast polarization modulation (based on the use of two liquid crystal retarders), real-time image accumulation, and dual beam polarimetry to reach polarization sensitivities of 0.1%. As a spectrograph, the instrument uses a LiNbO3 etalon in double pass and a narrow band pre-filter to achieve a spectral resolution of 85 mAA. IMaX uses the high Zeeman sensitive line of Fe I at 5250.2 AA and observes all four Stokes parameters at various points inside the spectral line. This allows vector magnetograms, Dopplergrams, and intensity frames to be produced that, after reconstruction, reach spatial resolutions in the 0.15-0.18 arcsec range over a 50x50 arcsec FOV. Time cadences vary between ten and 33 seconds, although the shortest one only includes longitudinal polarimetry. The spectral line is sampled in various ways depending on the applied observing mode, from just two points inside the line to 11 of them. All observing modes include one extra wavelength point in the nearby continuum. Gauss equivalent sensitivities are four Gauss for longitudinal fields and 80 Gauss for transverse fields per wavelength sample. The LOS velocities are estimated with statistical errors of the order of 5-40 m/s. The design, calibration and integration phases of the instrument, together with the implemented data reduction scheme are described in some detail.

The Sunrise balloon-borne solar observatory consists of a 1m aperture Gregory telescope, a UV filter imager, an imaging vector polarimeter, an image stabilization system and further infrastructure. The first science flight of Sunrise yielded high-quality data that reveal the structure, dynamics and evolution of solar convection, oscillations and magnetic fields at a resolution of around 100 km in the quiet Sun. After a brief description of instruments and data, first qualitative results are presented. In contrast to earlier observations, we clearly see granulation at 214 nm. Images in Ca II H display narrow, short-lived dark intergranular lanes between the bright edges of granules. The very small-scale, mixed-polarity internetwork fields are found to be highly dynamic. A significant increase in detectable magnetic flux is found after phase-diversity-related reconstruction of polarization maps, indicating that the polarities are mixed right down to the spatial resolution limit, and probably beyond.

We present the first observations of the transverse component of photospheric magnetic field acquired by the imaging magnetograph Sunrise/IMaX. Using an automated detection method, we obtain statistical properties of 4536 features with significant linear polarization signal. Their rate of occurrence is 1-2 orders of magnitude larger than values reported by previous studies. We show that these features have no characteristic size or lifetime. They appear preferentially at granule boundaries with most of them being caught in downflow lanes at some point in their evolution. Only a small percentage are entirely and constantly embedded in upflows (16%) or downflows (8%).

Based on IMaX/Sunrise data, we report on a previously undetected phenomenon in solar granulation. We show that in a very narrow region separating granules and intergranular lanes the spectral line width of the Fe I 5250.2 A line becomes extremely small. We offer an explanation of this observation with the help of magneto-convection simulations. These regions with extremely small line widths correspond to the places where the granular flows bend from mainly upflow in granules to downflow in intergranular lanes. We show that the resolution and image stability achieved by IMaX/Sunrise are important requisites to detect this interesting phenomenon.

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 study two limiting cases of turbulence forcing in numerical experiments: solenoidal (divergence-free) forcing, and compressive (curl-free) forcing, and compare our results to observations reported in the literature. We solve the equations of hydrodynamics on grids with up to 1024^3 cells for purely solenoidal and purely compressive forcing. Eleven lower-resolution models with mixtures of both forcings are also analysed. We find velocity dispersion--size relations consistent with observations and independent numerical simulations, irrespective of the type of forcing. However, compressive forcing yields stronger turbulent compression at the same RMS Mach number than solenoidal forcing, resulting in a three times larger standard deviation of volumetric and column density probability distributions (PDFs). We conclude that the strong dependence of the density PDF on the type of forcing must be taken into account in any theory using the PDF to predict properties of star formation. We supply a quantitative description of this dependence. We find that different observed regions show evidence of different mixtures of compressive and solenoidal forcing, with more compressive forcing occurring primarily in swept-up shells.

We update the SOHO/SWAN H Lyman-alpha brightness analysis to cover the 1996-2008 time interval. A forward model applied to the intensity maps provides the latitude and time dependence of the interstellar Hydrogen ionisation rate over more than a full solar cycle. The hydrogen ionisation, being almost entirely due to charge-exchange with solar wind ions, reflects closely the solar wind flux. Our results show that the solar wind latitudinal structure during the present solar minimum is strikingly different from the previous minimum, with a much wider slow solar wind equatorial belt which persists until at least the end of 2008. We compute absolute values of the in-ecliptic H ionisation rates using OMNI solar wind data and use them to calibrate our ionisation rates at all heliographic latitudes. We then compare the resulting fluxes with the synoptic LASCO/C2 electron densities at 6 solar radii. The two time-latitude patterns are strikingly similar over all the cycle. This comparison shows that densities at 6 solar radii can be used to infer the solar wind type close to its source, with high (resp. low) densities tracing the slow (resp. fast) solar wind, simply because the density reflects at which altitude occurs the acceleration. The comparison between the two minima suggests that the fast polar wind acceleration occurs at larger distance during the current minimum compared to the previous one. This difference, potentially linked to the magnetic field decrease or(and) the coronal temperature decrease should be reproduced by solar wind expansion models.

We investigate the influence of the turbulence forcing on the mass distributions of gravitationally unstable cores by postprocessing data from simulations of non-selfgravitating isothermal supersonic turbulence with varying resolution. In one set of simulations solenoidal forcing is applied, while the second set uses purely compressive forcing to excite turbulent motions. From the resulting density field, we compute the mass distribution of gravitationally unstable cores by means of a clump-finding algorithm. Using the time-averaged probability density functions of the mass density, semi-analytic mass distributions are calculated from analytical theories. We apply stability criteria that are based on the Bonnor-Ebert mass resulting from the thermal pressure and from the sum of thermal and turbulent pressure. Although there are uncertainties in the application of the clump-finding algorithm, we find systematic differences in the mass distributions obtained from solenoidal and compressive forcing. Compressive forcing produces a shallower slope in the high-mass power-law regime compared to solenoidal forcing. The mass distributions also depend on the Jeans length resulting from the choice of the mass in the computational box, which is freely scalable for non-selfgravitating isothermal turbulence. Provided that all cores are numerically resolved and most cores are small compared to the length scale of the forcing, the normalised core mass distributions are found to be close to the semi-analytic models. Especially for the high-mass tails, the Hennebelle-Chabrier theory implies that the additional support due to turbulent pressure is important.

The delayed detonation model describes the observational properties of the majority of type Ia supernovae very well. Using numerical data from a three-dimensional deflagration model for type Ia supernovae, the intermittency of the turbulent velocity field and its implications on the probability of a deflagration-to-detonation (DDT) transition are investigated. From structure functions of the turbulent velocity fluctuations, we determine intermittency parameters based on the log-normal and the log-Poisson models. On the other hand, the analysis of the turbulent velocity fluctuations in the vicinity of the flame front by Roepke suggests a much higher probability of large velocity fluctuations on the grid scale in comparison to the log-normal intermittency model. Following Pan et al., we computed probability density functions for a DDT for the different distributions. Assuming that a DDT can occur in the stirred flame regime, as proposed by Woosley et al., the log-normal model would imply a delayed detonation between 0.7 and 0.8 seconds after the beginning of the deflagration phase for the multi-spot ignition scenario used in the simulation. However, the probability drops to virtually zero if a DDT is further constrained by the requirement that the turbulent velocity fluctuations reach about 500 km/s. Under this condition, delayed detonations are only possible if the distribution of the velocity fluctuations is not log-normal. From our calculations follows that the distribution obtained by Roepke allow for multiple DDTs around 0.8 seconds after ignition at a transition density close to 1x10^7 g/cm^3.

Context: Turbulent deflagrations of Chandrasekhar mass White Dwarfs are commonly used to model Type Ia Supernova explosions. In this context, rapid rotation of the progenitor star is plausible but has so far been neglected. Aims: The aim of this work is to explore the influence of rapid rotation on the deflagration scenario. Methods: We use three dimensional hydrodynamical simulations to model turbulent deflagrations ignited within a variety of rapidly rotating CO WDs obeying rotation laws suggested by accretion studies. Results: We find that rotation has a significant impact on the explosion. The flame develops a strong anisotropy with a preferred direction towards the stellar poles, leaving great amounts of unburnt matter along the equatorial plane. Conclusions: The large amount of unburnt matter is contrary to observed spectral features of SNe Ia. Thus, rapid rotation of the progenitor star and the deflagration scenario are incompatible in order to explain SNe Ia.

Context: Superluminous type Ia supernovae (SNe Ia) may be explained by super-Chandrasekhar-mass explosions of rapidly rotating white dwarfs (WDs). In a preceding paper, we showed that the deflagration scenario applied to rapidly rotating WDs generates explosions that cannot explain the majority of SNe Ia. Aims: Rotation of the progenitor star allows super-Chandrasekhar-mass WDs to form that have a shallower density stratification. We use simple estimates of the production of intermediate and iron group elements in pure detonations of rapidly rotating WDs to assess their viability in explaining rare SNe Ia. Methods: We numerically construct WDs in hydrostatic equilibrium that rotate according to a variety of rotation laws. The explosion products are estimated by considering the density stratification and by evaluating the result of hydrodynamics simulations. Results: We show that a significant amount of intermediate mass elements is produced for theoretically motivated rotation laws, even for prompt detonations of WDs. Conclusions: Rapidly rotating WDs that detonate may provide an explanation of rare superluminous SNe Ia in terms of both burning species and explosion kinematics.

FEARLESS (Fluid mEchanics with Adaptively Refined Large Eddy SimulationS) is a new numerical scheme arising from the combined use of subgrid scale (SGS) model for turbulence at the unresolved length scales and adaptive mesh refinement (AMR) for resolving the large scales. This tool is especially suitable for the study of turbulent flows in strongly clumped media. In this contribution, the main features of FEARLESS are briefly outlined. We then summarize the main results of FEARLESS cosmological simulations of galaxy cluster evolution. In clusters, the production of turbulence is closely correlated with merger events; for minor mergers, we find that turbulent dissipation affects the cluster energy budget only locally. The level of entropy in the cluster core is enhanced in FEARLESS simulations, in accord with a better modeling of the unresolved flow, and with its feedback on the resolved mixing in the ICM.

We analyze the statistics and star formation rate obtained in high-resolution numerical experiments of forced supersonic turbulence, and compare with observations. We concentrate on a systematic comparison of solenoidal (divergence-free) and compressive (curl-free) forcing, which are two limiting cases of turbulence driving. Our results show that for the same RMS Mach number, compressive forcing produces a three times larger standard deviation of the density probability distribution. When self-gravity is included in the models, the star formation rate is more than one order of magnitude higher for compressive forcing than for solenoidal forcing.

Based on ideas by Woodward et al., a subgrid scale model that is applicable to highly compressible turbulence is presented. Applying the subgrid scale model in large eddy simulations of forced supersonic turbulence, the bottleneck effect is largely reduced and, thereby, approximate scaling laws can be obtained at relatively low numerical resolution. In agreement with previous results from PPM simulations without explicit subgrid scale model, it is found that the energy spectrum function for the velocity field with fractional density-weighing, rho^(1/3)u, varies substantially with the forcing, at least for the decade of wavenumbers next to the energy-containing range. Consequently, if universal scaling of compressible turbulence exists, it can be found on length scales much smaller than the forcing scale only.

Contradicting results have been reported in the literature with respect to the performance of the numerical techniques employed for the study of supersonic turbulence. We aim at characterising the performance of different particle-based and grid-based techniques on the modelling of decaying supersonic turbulence. Four different grid codes (ENZO, FLASH, TVD, ZEUS) and three different SPH codes (GADGET, PHANTOM, VINE) are compared. We additionally analysed two calculations denoted as PHANTOM A and PHANTOM B using two different implementations of artificial viscosity. Our analysis indicates that grid codes tend to be less dissipative than SPH codes, though details of the techniques used can make large differences in both cases. For example, the Morris & Monaghan viscosity implementation for SPH results in less dissipation (PHANTOM B and VINE versus GADGET and PHANTOM A). For grid codes, using a smaller diffusion parameter leads to less dissipation, but results in a larger bottleneck effect (our ENZO versus FLASH runs). As a general result, we find that by using a similar number of resolution elements N for each spatial direction means that all codes (both grid-based and particle-based) show encouraging similarity of all statistical quantities for isotropic supersonic turbulence on spatial scales k<N/32 (all scales resolved by more than 32 grid cells), while scales smaller than that are significantly affected by the specific implementation of the algorithm for solving the equations of hydrodynamics. At comparable numerical resolution, the SPH runs were on average about ten times more computationally intensive than the grid runs, although with variations of up to a factor of ten between the different SPH runs and between the different grid runs. (abridged)

We present a numerical scheme for modelling unresolved turbulence in cosmological adaptive mesh refinement codes. As a first application, we study the evolution of turbulence in the intra-cluster medium and in the core of a galaxy cluster. Simulations with and without subgrid scale model are compared in detail. Since the flow in the ICM is subsonic, the global turbulent energy contribution at the unresolved length scales is smaller than 1% of the internal energy. We find that the production of turbulence is closely correlated with merger events occurring in the cluster environment, and its dissipation locally affects the cluster energy budget. Because of this additional source of dissipation, the core temperature is larger and the density is smaller in the presence of subgrid scale turbulence than in the standard adiabatic run, resulting in a higher entropy core value.

Colliding flows are a commonly used scenario for the formation of molecular clouds in numerical simulations. Due to the thermal instability of the warm neutral medium, turbulence is produced by cooling. We carry out a two-dimensional numerical study of such colliding flows in order to test whether statistical properties inferred from adaptive mesh refinement (AMR) simulations are robust with respect to the applied refinement criteria. We compare probability density functions of various quantities as well as the clump statistics and fractal dimension of the density fields in AMR simulations to a static-grid simulation. The static grid with 2048^2 cells matches the resolution of the most refined subgrids in the AMR simulations. The density statistics is reproduced fairly well by AMR. Refinement criteria based on the cooling time or the turbulence intensity appear to be superior to the standard technique of refinement by overdensity. Nevertheless, substantial differences in the flow structure become apparent. In general, it is difficult to separate numerical effects from genuine physical processes in AMR simulations.

In a systematic study, we compare the density statistics in high resolution numerical experiments of supersonic isothermal turbulence, driven by the usually adopted solenoidal (divergence-free) forcing and by compressive (curl-free) forcing. We find that for the same rms Mach number, compressive forcing produces much stronger density enhancements and larger voids compared to solenoidal forcing. Consequently, the Fourier spectra of density fluctuations are significantly steeper. This result is confirmed using the Delta-variance analysis, which yields power-law exponents beta~3.4 for compressive forcing and beta~2.8 for solenoidal forcing. We obtain fractal dimension estimates from the density spectra and Delta-variance scaling, and by using the box counting, mass size and perimeter area methods applied to the volumetric data, projections and slices of our turbulent density fields. Our results suggest that compressive forcing yields fractal dimensions significantly smaller compared to solenoidal forcing. However, the actual values depend sensitively on the adopted method, with the most reliable estimates based on the Delta-variance, or equivalently, on Fourier spectra. Using these methods, we obtain D~2.3 for compressive and D~2.6 for solenoidal forcing, which is within the range of fractal dimension estimates inferred from observations (D~2.0-2.7). The velocity dispersion to size relations for both solenoidal and compressive forcing obtained from velocity spectra follow a power law with exponents in the range 0.4-0.5, in good agreement with previous studies.

We analyze the statistical properties of the turbulent velocity field in the deflagration model for Type Ia supernovae. In particular, we consider the question of whether turbulence is isotropic and consistent with the Kolmogorov theory at small length scales. Using numerical data from a high-resolution simulation of a thermonuclear supernova explosion, spectra of the turbulence energy and velocity structure functions are computed. We show that the turbulent velocity field is isotropic at small length scales and follows a scaling law that is consistent with the Kolmogorov theory until most of the nuclear fuel is burned. At length scales greater than a certain characteristic scale, turbulence becomes anisotropic. Here, the radial velocity fluctuations follow the scaling law of the Rayleigh-Taylor instability, whereas the angular component still obeys Kolmogorov scaling. In the late phase of the explosion, this characteristic scale drops below the numerical resolution of the simulation. The analysis confirms that a subgrid-scale model for the unresolved turbulence energy is required for the consistent calculation of the flame speed in deflagration models of Type Ia supernovae, and that the assumption of isotropy on these scales is appropriate.

The probability density function (PDF) of the gas density in turbulent supersonic flows is investigated with high-resolution numerical simulations. In a systematic study, we compare the density statistics of compressible turbulence driven by the usually adopted solenoidal forcing (divergence-free) and by compressive forcing (curl-free). Our results are in agreement with studies using solenoidal forcing. However, compressive forcing yields a significantly broader density distribution with standard deviation ~3 times larger at the same rms Mach number. The standard deviation-Mach number relation used in analytical models of star formation is reviewed and a modification of the existing expression is proposed, which takes into account the ratio of solenoidal and compressive modes of the turbulence forcing.

The statistical properties of turbulence are considered to be universal at sufficiently small length scales, i. e., independent of boundary conditions and large-scale forces acting on the fluid. Analyzing data from numerical simulations of supersonic turbulent flow driven by external forcing, we demonstrate that this is not generally true for the two-point velocity statistics of compressible turbulence. However, a reformulation of the refined similarity hypothesis in terms of the mass-weighted velocity rho^(1/3)v yields scaling laws that are almost insensitive to the forcing. The results imply that the most intermittent dissipative structures are shocks closely following the scaling of Burgers turbulence.

A direct measurement of the universe's expansion history could be made by observing in real time the evolution of the cosmological redshift of distant objects. However, this would require measurements of Doppler velocity drifts of about 1 centimeter per second per year, and astronomical spectrographs have not yet been calibrated to this tolerance. We demonstrate the first use of a laser frequency comb for wavelength calibration of an astronomical telescope. Even with a simple analysis, absolute calibration is achieved with an equivalent Doppler precision of approximately 9 meters per second at about 1.5 micrometers - beyond state-of-the-art accuracy. We show that tracking complex, time-varying systematic effects in the spectrograph and detector system is a particular advantage of laser frequency comb calibration. This technique promises an effective means for modeling and removal of such systematic effects to the accuracy required by future experiments to see direct evidence of the universe's putative acceleration.

We performed numerical simulations of supersonic isothermal turbulence driven by mostly compressive large-scale forcing, using both a static grid and adaptive mesh refinement with an effective resolution N=768^3. After a transient phase dominated by shocks, turbulence evolves into a steady state with an RMS Mach number about 2.5, in which cloud-like structures of over-dense gas are surrounded by highly rarefied gas. The index of the turbulence energy spectrum function beta = 2.0 in the shock-dominated phase. As the flow approaches statistical equilibrium, the spectrum flattens, with beta = 1.9. For the scaling exponent of the root mean square velocity fluctuation, we obtain gamma = 0.43 from the velocity structure functions of second order. These results are well within the range of observed scaling properties for the velocity dispersion in molecular clouds. Calculating structure functions of order p=1,...,5, we find for all scaling exponents significant deviations from the Kolmogorov-Burgers model proposed by Boldyrev. Our results are very well described by a general log-Poisson model with a higher degree of intermittency, which implies an influence of the forcing on the scaling properties. Contrary to previous numerical results for isothermal turbulence, we obtain a skewed probability density function of the mass density fluctuations that is not consistent with log-normal statistics and entails a substantially higher fraction of mass in the density peaks than implied by the Padoan-Nordlund relation between the variance of the density fluctuations and the Mach number. In conclusion, it seems necessary to account for the production mechanism of turbulence in the ISM.

The problem of the resolution of turbulent flows in adaptive mesh refinement (AMR) simulations is investigated by means of 3D hydrodynamical simulations in an idealised setup, representing a moving subcluster during a merger event. AMR simulations performed with the usual refinement criteria based on local gradients of selected variables do not properly resolve the production of turbulence downstream of the cluster. Therefore we apply novel AMR criteria which are optimised to follow the evolution of a turbulent flow. We demonstrate that these criteria provide a better resolution of the flow past the subcluster, allowing us to follow the onset of the shear instability, the evolution of the turbulent wake and the subsequent back-reaction on the subcluster core morphology. We discuss some implications for the modelling of cluster cold fronts.

We use 3-dimensional numerical simulations of self-gravitating compressible turbulent gas in combination with Lagrangian tracer particles to investigate the mixing process of molecular hydrogen (H2) in interstellar clouds. Tracer particles are used to represent shock-compressed dense gas, which is associated with H2. We deposit tracer particles in regions of density contrast in excess of ten times the mean density. Following their trajectories and using probability distribution functions, we find an upper limit for the mixing timescale of H2, which is of order 0.3 Myr. This is significantly smaller than the lifetime of molecular clouds, which demonstrates the importance of the turbulent mixing of H2 as a preliminary stage to star formation.

We study chromospheric emission to understand the temperature stratification in the solar chromosphere. We observed the intensity profile of the CaIIH line in a quiet Sun region close to the disk center at the German Vacuum Tower Telescope. We analyze over 10^5 line profiles from inter-network regions. For comparison with the observed profiles, we synthesize spectra for a variety of model atmospheres with a non local thermodynamic equilibrium (NLTE) radiative transfer code. A fraction of about 25% of the observed CaIIH line profiles do not show a measurable emission peak in H_{2v} and H_{2r} wavelength bands (reversal-free). All of the chosen model atmospheres with a temperature rise fail to reproduce such profiles. On the other hand, the synthetic calcium profile of a model atmosphere that has a monotonic decline of the temperature with height shows a reversal-free profile that has much lower intensities than any observed line profile. The observed reversal-free profiles indicate the existence of cool patches in the interior of chromospheric network cells, at least for short time intervals. Our finding is not only in conflict with a full-time hot chromosphere, but also with a very cool chromosphere as found in some dynamic simulations.

We analyze a 1-hour time series of Ca II H intensity spectra and polarimetric spectra around 630 nm to derive the signature of the chromospheric heating and to investigate its relation to magnetic fields. We derived several characteristic quantities of Ca II H to define the chromospheric atmosphere properties. We study the power of the Fourier transform at different wavelengths and their phase relations. We perform local thermodynamic equilibrium inversions of the data to obtain the magnetic field, once including the Ca spectra. We find that the emission in the Ca II H line core at locations without detectable photospheric polarization signal is due to waves that propagate from low forming continuum layers in the line wing up to the line core. The phase differences of intensity oscillations at different wavelengths indicate standing waves for v < 2 mHz and propagating waves for higher frequencies. The waves steepen into shocks in the chromosphere. On average, shocks are both preceded and followed by intensity reductions. In field-free regions, the profiles show emission about half of the time. The correlation between wavelengths and the decorrelation time is significantly higher in the presence of magnetic fields than for field-free areas. The average Ca II H profile in the presence of magnetic fields contains emission features symmetric to the line core and an asymmetric contribution, where mainly the blue H2V emission peak is increased. We find that acoustic waves steepening into shocks are responsible for the emission in the Ca II H line core for locations without photospheric magnetic fields. We suggest using wavelengths in the line wing of Ca II H, where LTE still applies, to compare theoretical heating models with observations.

In this contribution, I give an overview of the various approaches toward the numerical modelling of turbulence, particularly, in the interstellar medium. The discussion is placed in a physical context, i. e. computational problems are motivated from basic physical considerations. Presenting selected examples for solutions to these problems, I introduce the basic ideas of the most commonly used numerical methods.

We study the structure and evolution of the magnetic field of the quiet Sun by investigating weak spectro-polarimetric signals. To this end, we observed a quiet region close to the disk center with the German VTT in Tenerife, July 07, 2006. We recorded 38 scans of the same area. Each scan was eight arcsec wide and observed within about 100 seconds. We used POLIS to simultaneously observe Stokes profiles of the neutral iron lines at 630.15 and 630.25 nm, the Stokes-I profile of the Ca II H line at 396.8 nm, and a continuum speckle channel at 500 nm. We witness two examples of magnetic flux cancellation of small-scale opposite-polarity patches, followed by an enhanced chromospheric emission. In each case, the two opposite-polarity patches gradually became smaller and, within a few minutes, the smaller one completely disappeared. The larger patch also diminished significantly. We provide evidence for a cancellation scenario in the photosphere which leaves minor traces at the chromospheric level.

The idea of this work is to compare a new positive and entropy stable approximate Riemann solver by Francois Bouchut with a state-of the-art algorithm for astrophysical fluid dynamics. We implemented the new Riemann solver into an astrophysical PPM-code, the Prometheus code, and also made a version with a different, more theoretically grounded higher order algorithm than PPM. We present shock tube tests, two-dimensional instability tests and forced turbulence simulations in three dimensions. We find subtle differences between the codes in the shock tube tests, and in the statistics of the turbulence simulations. The new Riemann solver increases the computational speed without significant loss of accuracy.

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 ATIC balloon-borne experiment measures the energy spectra of elements from H to Fe in primary cosmic rays from about 100 GeV to 100 TeV. ATIC is comprised of a fully active bismuth germanate calorimeter, a carbon target with embedded scintillator hodoscopes, and a silicon matrix that is used as the main charge detector. The silicon matrix produces good charge resolution for protons and helium but only partial resolution for heavier nuclei. In the present paper, the charge resolution of ATIC was improved and backgrounds were reduced in the region from Be to Si by using the upper layer of the scintillator hodoscope as an additional charge detector. The flux ratios of nuclei B/C, C/O, N/O in the energy region from about 10 GeV/nucleon to 300 GeV/nucleon obtained from this high-resolution, high-quality charge spectra are presented, and compared with existing theoretical predictions.

A simulation of the thermonuclear explosion of a Chandrasekhar-mass C+O white dwarf, the most popular scenario of a type Ia supernova (SN Ia), is presented. The underlying modeling is pursued in a self-consistent way, treating the combustion wave as a turbulent deflagration using well tested methods developed for laboratory combustion and based on the concept of `large eddy simulations' (LES). Such consistency requires to capture the onset of the turbulent cascade on resolved scales. This is achieved by computing the dynamical evolution on a 1024$^3$ moving grid, which resulted in the best-resolved three-dimensional SN Ia simulation carried out thus far, reaching the limits of what can be done on present supercomputers. Consequently, the model has no free parameters other than the initial conditions at the onset of the explosion, and therefore it has considerable predictive power. Our main objective is to determine to which extent such a simulation can account for the observations of normal SNe Ia. Guided by previous simulations with less resolution and a less sophisticated flame model, initial conditions were chosen that yield a reasonably strong explosion and a sufficient amount of radioactive nickel for a bright display. We show that observables are indeed matched to a reasonable degree. In particular, good agreement is found with the light curves of normal SNe Ia. Moreover, the model reproduces the general features of the abundance stratification as inferred from the analysis of spectra. This indicates that it captures the main features of the explosion mechanism of SNe Ia. However, we also show that even a seemingly best-choice pure deflagration model has shortcomings that indicate the need for a different mode of nuclear burning at late times, perhaps the transition to a detonation at low density.

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.

GREGOR is a new open solar telescope with an aperture of 1.5m. It replaces the former 45-cm Gregory Coude telescope on the Canary island Tenerife. The optical concept is that of a double Gregory system. The main and the elliptical mirrors are made from a silicon-carbide material with high thermal conductivity. This is important to keep the mirrors on the ambient temperature avoiding local turbulence. GREGOR will be equipped with an adaptive optics system. The new telescope will be ready for operation in 2008. Post-focus instruments in the first stage will be a spectrograph for polarimetry in the near infrared and a 2-dimensional spectrometer based on Fabry-Perot interferometers for the visible.

We simultaneously observed the Stokes parameters of the photospheric iron line pair at 630.2 nm and the intensity profile of the chromospheric Ca II H line at 396.8 nm in a quiet Sun region at a heliocentric angle of 53 deg. We perform a statistical analysis of network and inter-network properties.The H-index is the integrated emission in a 0.1 nm band around the Ca core. We separate a non-magnetically, H_non, and a magnetically, H_mag, heated component from a non-heated component, H_co in the H-index. The average network and inter-network H-indices are equal to 12 and 10 pm, respectively. The emission in the network is correlated with the magnetic flux density, approaching a value of H 10 pm for vanishing flux. The inter-network magnetic field is dominated by weak field strengths with values down to 200 G and its mean absolute flux density is 11 Mx cm$^{-2}$. We find that a dominant fraction of the calcium emission caused by the heated atmosphere in the magnetic network, has non-magnetic origin (H_mag = 2pm, H_non = 3pm). Considering the effect of straylight, the contribution from an atmosphere with no temperature rise to the H-index (H_co = 6pm) is about half of the observed H-index in the inter-network. The H-index in the inter-network is not correlated to any property of the photospheric magnetic field, suggesting that magnetic flux concentrations have a negligible role in the chromospheric heating in this region. The height range of the thermal coupling between the photosphere and low/mid chromosphere increases in presence of magnetic field. In addition, we demonstrate that a poor signal-to-noise level in the Stokes profiles leads to a significant over-estimation of the magnetic field strength.

We present a statistical analysis of network and internetwork properties in the photosphere and the chromosphere. For the first time we simultaneously observed (a) the four Stokes parameters of the photospheric iron line pair at 630.2 nm and (b) the intensity profile of the Ca II H line at 396.8 nm. The vector magnetic field was inferred from the inversion of the iron lines. We aim at an understanding of the coupling between photospheric magnetic field and chromospheric emission.

In thermonuclear supernovae, intermediate mass elements are mostly produced by distributed burning provided that a deflagration to detonation transition does not set in. Apart from the two-dimensional study by Roepke & Hillebrandt (2005), very little attention has been payed so far to the correct treatment of this burning regime in numerical simulations. In this article, the physics of distributed burning is reviewed from the literature on terrestrial combustion and differences which arise from the very small Prandtl numbers encountered in degenerate matter are pointed out. Then it is shown that the level set method continues to be applicable beyond the flamelet regime as long as the width of the flame brush does not become smaller than the numerical cutoff length. Implementing this constraint with a simple parameterisation of the effect of turbulence onto the energy generation rate, the production of intermediate mass elements increases substantially compared to previous simulations, in which the burning process was stopped once the mass density dropped below 10^7 g/cm^3. Although these results depend on the chosen numerical resolution, an improvement of the constraints on the the total mass of burning products in the pure deflagration scenario can be achieved.

The results of ATIC-2 baloon experiment (2002-2003) for energy spectra of protons, He, C, O, Ne, Mg, Si, Fe, some groups of nuclei, and all-particle spectrum in primary cosmic rays are presented in energy region 50GeV--200TeV. The conclusion is that the spectra of protons and helium nuclei are essentially different (the spectrum of protons is steeper) and the spectra of protons and heavy nuclei have no-power form.

We investigate the morphological, dynamical, and evolutionary properties of the internetwork and network fine structure of the quiet sun at disk centre. The analysis is based on a $\sim$6 h time sequence of narrow-band filtergrams centred on the inner-wing \Ca II K$_{\rm 2v}$ reversal at 393.3 nm. The results for the internetwork are related to predictions derived from numerical simulations of the quiet sun. The average evolutionary time scale of the internetwork in our observations is 52 sec. Internetwork grains show a tendency to appear on a mesh-like pattern with a mean cell size of $\sim$4-5 arcsec. Based on this size and the spatial organisation of the mesh we speculate that this pattern is related to the existence of photospheric downdrafts as predicted by convection simulations. The image segmentation shows that typical sizes of both network and internetwork grains are in the order of 1.6 arcs.

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.

Interplanetary Lyman alpha line profiles are derived from the SWAN H cell data measurements. The measurements cover a 6-year period from solar minimum (1996) to after the solar maximum of 2001. This allows us to study the variations of the line profiles with solar activity. These line profiles were used to derive line shifts and line widths in the interplanetary medium for various angles of the LOS with the interstellar flow direction. The SWAN data results were then compared to an interplanetary background upwind spectrum obtained by STIS/HST in March 2001. We find that the LOS upwind velocity associated with the mean line shift of the IP \lya line varies from 25.7 km/s to 21.4 km/s from solar minimum to solar maximum. Most of this change is linked with variations in the radiation pressure. LOS kinetic temperatures derived from IP line widths do not vary monotonically with the upwind angle of the LOS. This is not compatible with calculations of IP line profiles based on hot model distributions of interplanetary hydrogen. We also find that the line profiles get narrower during solar maximum. The results obtained on the line widths (LOS temperature) show that the IP line is composed of two components scattered by two hydrogen populations with different bulk velocities and temperature. This is a clear signature of the heliospheric interface on the line profiles seen at 1 AU from the sun.

The dynamics of the explosive burning process is highly sensitive to the flame speed model in numerical simulations of type Ia supernovae. Based upon the hypothesis that the effective flame speed is determined by the unresolved turbulent velocity fluctuations, we employ a new subgrid scale model which includes a localised treatment of the energy transfer through the turbulence cascade in combination with semi-statistical closures for the dissipation and non-local transport of turbulence energy. In addition, subgrid scale buoyancy effects are included. In the limit of negligible energy transfer and transport, the dynamical model reduces to the Sharp-Wheeler relation. According to our findings, the Sharp-Wheeler relation is insuffcient to account for the complicated turbulent dynamics of flames in thermonuclear supernovae. The application of a co-moving grid technique enables us to achieve very high spatial resolution in the burning region. Turbulence is produced mostly at the flame surface and in the interior ash regions. Consequently, there is a pronounced anisotropy in the vicinity of the flame fronts. The localised subgrid scale model predicts significantly enhanced energy generation and less unburnt carbon and oxygen at low velocities compared to earlier simulations.

We present a one-equation subgrid scale model that evolves the turbulence energy corresponding to unresolved velocity fluctuations in large eddy simulations. The model is derived in the context of the Germano consistent decomposition of the hydrodynamical equations. The eddy-viscosity closure for the rate of energy transfer from resolved toward subgrid scales is localised by means of a dynamical procedure for the computation of the closure parameter. Therefore, the subgrid scale model applies to arbitrary flow geometry and evolution. For the treatment of microscopic viscous dissipation a semi-statistical approach is used, and the gradient-diffusion hypothesis is adopted for turbulent transport. A priori tests of the localised eddy-viscosity closure and the gradient-diffusion closure are made by analysing data from direct numerical simulations. As an a posteriori testing case, the large eddy simulation of thermonuclear combustion in forced isotropic turbulence is discussed. We intend the formulation of the subgrid scale model in this paper as a basis for more advanced applications in numerical simulations of complex astrophysical phenomena involving turbulence.

We apply an ad hoc model for dynamical ignition in three-dimensional numerical simulations of thermonuclear supernovae assuming pure deflagrations. The model makes use of the statistical description of temperature fluctuations in the pre-supernova core proposed by Wunsch & Woosley (2004). Randomness in time is implemented by means of a Poisson process. We are able to vary the explosion energy and nucleosynthesis depending on the free parameter of the model which controls the rapidity of the ignition process. However, beyond a certain threshold, the strength of the explosion saturates and the outcome appears to be robust with respect to number of ignitions. In the most energetic explosions, we find about 0.75 solar masses of iron group elements. Other than in simulations with simultaneous multi-spot ignition, the amount of unburned carbon and oxygen at radial velocities of a few 1000 km/s tends to be reduced for an ever increasing number of ignition events and, accordingly, more pronounced layering results.

We study the dynamics of thermonuclear flames propagating in fuel stirred by stochastic forcing. The fuel consists of carbon and oxygen in a state which is encountered in white dwarfs close to the Chandrasekhar limit. The level set method is applied to represent the flame fronts numerically. The computational domain for the numerical simulations is cubic, and periodic boundary conditions are imposed. The goal is the development of a suitable flame speed model for the small-scale dynamics of turbulent deflagration in thermonuclear supernovae. Because the burning process in a supernova explosion is transient and spatially inhomogeneous, the localised determination of subgrid scale closure parameters is essential. We formulate a semi-localised model based on the dynamical equation for the subgrid scale turbulence energy $k_{\mathrm{sgs}}$. The turbulent flame speed $s_{\mathrm{t}}$ is of the order $\sqrt{2k_{\mathrm{sgs}}}$. In particular, the subgrid scale model features a dynamic procedure for the calculation of the turbulent energy transfer from resolved toward subgrid scales, which has been successfully applied to combustion problems in engineering. The options of either including or suppressing inverse energy transfer in the turbulence production term are compared. In combination with the piece-wise parabolic method for the hydrodynamics, our results favour the latter option. Moreover, different choices for the constant of proportionality in the asymptotic flame speed relation, $s_{\mathrm{t}}\propto\sqrt{2k_{\mathrm{sgs}}}$, are investigated.

The piece-wise parabolic method (PPM) is applied to simulations of forced isotropic turbulence with Mach numbers $\sim 0.1... 1$. The equation of state is dominated by the Fermi pressure of an electron-degenerate fluid. The dissipation in these simulations is of purely numerical origin. For the dimensionless mean rate of dissipation, we find values in agreement with known results from mostly incompressible turbulence simulations. The calculation of a Smagorinsky length corresponding to the rate of numerical dissipation supports the notion of the PPM supplying an implicit subgrid scale model. In the turbulence energy spectra of various flow realisations, we find the so-called bottleneck phenomenon, i.e., a flattening of the spectrum function near the wavenumber of maximal dissipation. The shape of the bottleneck peak in the compensated spectrum functions is comparable to what is found in turbulence simulations with hyperviscosity. Although the bottleneck effect reduces the range of nearly inertial length scales considerably, we are able to estimate the value of the Kolmogorov constant. For steady turbulence with a balance between energy injection and dissipation, it appears that $C\approx 1.7$. However, a smaller value is found in the case of transonic turbulence with a large fraction of compressive components in the driving force. Moreover, we discuss length scales related to the dissipation, in particular, an effective numerical length scale $\Delta_{\mathrm{eff}}$, which can be regarded as the characteristic smoothing length of the implicit filter associated with the PPM.

We present numerical calculations of the parameters $C_{\nu}$, $C_{\epsilon}$ and $C_{\kappa}$ associated with the common closures for turbulence production, dissipation and diffusion. In the case of homogeneous and isotropic turbulence, these parameters are expected to be statistically scale-invariant within the inertial subrange. In order to scrutinise this conjecture, we utilised a generalisation of the Germano filtering formalism, which is applicable to compressible flows as well. The filtering of data obtained from three-dimensional direct numerical simulations of forced isotropic turbulence with Mach numbers in the range $\sim 0.1...1$ then yielded values of the closure parameters associated with different length scales. The results indicate that the closure parameters are nearly universal for subsonic or moderately transonic flows, although the resolution of $432^{3}$ grid cells in our simulations is not quite sufficient to clearly establish scale invariance. In addition, it was found that the customary assumption of a kinetic Prandtl number of about unity for the gradient-diffusion closure is flawed due to the misalignment between turbulent flux and the gradient of the turbulence energy. Nevertheless, sound correlation can be achieved if the flux magnitude rather than the flux vector is locally matched. This conclusion is particularly useful for the family of subgrid scale models based on the turbulence energy equation. Furthermore, the parameter of production $C_{\nu}$ was computed in the fashion of dynamical procedures. Thereby, superior agreement between modelled and explicitly evaluated turbulence stresses in comparison to the eddy-viscosity closure with constant $C_{\nu}$ was verified.

Type Ia supernovae have recently received considerable attention because it appears that they can be used as "standard candles" to measure cosmic distances out to billions of light years away from us. Observations of type Ia supernovae seem to indicate that we are living in a universe that started to accelerate its expansion when it was about half its present age. These conclusions rest primarily on phenomenological models which, however, lack proper theoretical understanding, mainly because the explosion process, initiated by thermonuclear fusion of carbon and oxygen into heavier elements, is difficult to simulate even on supercomputers. Here, we investigate a new way of modeling turbulent thermonuclear deflagration fronts in white dwarfs undergoing a type Ia supernova explosion. Our approach is based on a level set method which treats the front as a mathematical discontinuity and allows for full coupling between the front geometry and the flow field. New results of the method applied to the problem of type Ia supernovae are obtained. It is shown that in 2-D with high spatial resolution and a physically motivated subgrid scale model for the nuclear flames numerically "converged" results can be obtained, but for most initial conditions the stars do not explode. In contrast, simulations in 3-D, do give the desired explosions and many of their properties, such as the explosion energies, lightcurves and nucleosynthesis products, are in very good agreement with observed type Ia supernovae.

We present several new sets of grids of model stellar atmospheres computed with modified versions of the ATLAS9 code. Each individual set consists of several grids of models with different metallicities ranging from [M/H] = -2.0 to +1.0 dex. The grids range from 4000 to 10000 K in T_eff and from 2.0 to 5.0 dex in logg. The individual sets differ from each other and from previous ones essentially in the physics used for the treatment of the convective energy transport, in the higher vertical resolution of the atmospheres and in a finer grid in the (T_eff, logg) plane. These improvements enable the computation of derivatives of color indices accurate enough for pulsation mode identification. In addition, we show that the chosen vertical resolution is necessary and sufficient for the purpose of stellar interior modelling. To explain the physical differences between the model grids we provide a description of the currently available modifications of ATLAS9 according to their treatment of convection. Our critical analysis of the dependence of the atmospheric structure and observable quantities on convection treatment, vertical resolution and metallicity reveals that spectroscopic and photometric observations are best represented when using an inefficient convection treatment. This conclusion holds whatever convection formulation investigated here is used, i.e. MLT(alpha=0.5), CM and CGM are equivalent. We also find that changing the convection treatment can lead to a change in the effective temperature estimated from Stroemgren color indices from 200 to 400 K.

The detection of a new 16.4 micron emission feature in the ISO-SWS spectra of NGC 7023, M17, and the Orion Bar is reported. Previous laboratory experiments measured a mode near this wavelength in spectra of PAHs (Polycyclic Aromatic Hydrocarbons), and so we suggest the new interstellar 16.4 micron feature could be assigned to low-frequency vibrations of PAHs. The best carrier candidates seem to be PAH molecules containing pentagonal rings.