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

Context. Large-amplitude oscillations (LAOs) of the solar prominences are very spectacular but poorly understood phenomena. These motions have amplitudes larger than $10 {\mathrm{\, km \,s^{-1}}}$ and can be triggered by the external perturbations, e.g., Moreton or EIT waves. Aims. Our aim is to analyze the properties of large-amplitude oscillations using realistic prominence models and the triggering mechanism by external disturbances. Methods. We perform time-dependent numerical simulations of LAOs using a magnetic flux rope model with two values of the shear angle and the density contrast. We study the internal modes of the prominence using the horizontal and vertical triggering. In addition, we use the perturbation that arrives from outside in order to understand how such external disturbance can produce LAOs. Results. The period of longitudinal oscillations and its behavior with height show good agreement with the pendulum model. The period of transverse oscillations remains constant with height, suggesting a global normal mode. The transverse oscillations typically have shorter periods than the longitudinal oscillations. Conclusions. The periods of the longitudinal and transverse oscillations show only weak dependence on the shear angle of the magnetic structure and the prominence density contrast. The external disturbance perturbs the flux rope exciting oscillations of both polarizations. Their properties are a mixture of those excited by purely horizontal and vertical driving.

The European Solar Telescope (EST) is a research infrastructure for solar physics. It is planned to be an on-axis solar telescope with an aperture of 4 m and equipped with an innovative suite of spectro-polarimetric and imaging post-focus instrumentation. The EST project was initiated and is driven by EAST, the European Association for Solar Telescopes. EAST was founded in 2006 as an association of 14 European countries. Today, as of December 2019, EAST consists of 26 European research institutes from 18 European countries. The Preliminary Design Phase of EST was accomplished between 2008 and 2011. During this phase, in 2010, the first version of the EST Science Requirement Document (SRD) was published. After EST became a project on the ESFRI roadmap 2016, the preparatory phase started. The goal of the preparatory phase is to accomplish a final design for the telescope and the legal governance structure of EST. A major milestone on this path is to revisit and update the Science Requirement Document (SRD). The EST Science Advisory Group (SAG) has been constituted by EAST and the Board of the PRE-EST EU project in November 2017 and has been charged with the task of providing with a final statement on the science requirements for EST. Based on the conceptual design, the SRD update takes into account recent technical and scientific developments, to ensure that EST provides significant advancement beyond the current state-of-the-art. The present update of the EST SRD has been developed and discussed during a series of EST SAG meetings. The SRD develops the top-level science objectives of EST into individual science cases. Identifying critical science requirements is one of its main goals. Those requirements will define the capabilities of EST and the post-focus instrument suite. The technical requirements for the final design of EST will be derived from the SRD.

Aims. We aim to improve the understanding of the physical mechanisms behind the slower than free-fall motion and the two-stage evolution (an initial phase of acceleration followed by an almost constant velocity phase) detected in coronal rain events. Methods. Using the Mancha3D code, we solve the 2D ideal MHD equations. We represent the solar corona as an isothermal vertically stratified atmosphere with a uniform vertical magnetic field and the plasma condensation as a density enhancement described by a 2D Gaussian profile. We analyse the temporal evolution of the descending plasma and study its dependence on parameters such as density and magnetic field strength. Results. We confirm previous findings that the pressure gradient is the main force that opposes the action of gravity and slows down the blob descent and that larger densities require larger pressure gradients to reach the constant speed phase. We find that the shape of a condensation with a horizontal variation of density is distorted as it falls, due to the denser parts of the blob falling faster than the lighter ones. This is explained by the fact that the duration of the initial acceleration phase, and therefore the maximum falling speed attained by the plasma, increases with the ratio of blob to coronal density. We also find that the magnetic field plays a fundamental role in the evolution of the descending condensations. A strong enough magnetic field (greater than 10 G in our simulations) forces each plasma element to follow the path given by a particular field line, which allows to describe the evolution of each vertical slice of the blob in terms of 1D dynamics, without influence of the adjacent slices. In addition, under the typical conditions of the coronal rain events, the magnetic field prevents the development of the Kelvin-Helmholtz instability.

Waves and shocks traveling through the solar chromospheric plasma are influenced by its partial ionization and weak collisional coupling, and may become susceptible to multi-fluid effects, similar to interstellar shock waves. In this study, we consider fast magneto-acoustic shock wave formation and propagation in a stratified medium, that is permeated by a horizontal magnetic field, with properties similar to that of the solar chromosphere. The evolution of plasma and neutrals is modeled using a two-fluid code that evolves a set of coupled equations for two separate fluids. We observed that waves in neutrals and plasma, initially coupled at the upper photosphere, become uncoupled at higher heights in the chromosphere. This decoupling can be a consequence of either the characteristic spatial scale at the shock front, that becomes similar to the collisional scale, or the change in the relation between the wave frequency, ion cyclotron frequency, and the collisional frequency with height. The decoupling height is a sensitive function of the wave frequency, wave amplitude, and the magnetic field strength. We observed that decoupling causes damping of waves and an increase in the background temperature due to the frictional heating. The comparison between analytical and numerical results allows us to separate the role of the nonlinear effects from the linear ones on the decoupling and damping of waves.

Solar chromosphere consists of a partially ionized plasma, which makes modeling the solar chromosphere a particularly challenging numerical task. Here we numerically model chromospheric waves using a two-fluid approach with a newly developed numerical code. The code solves two-fluid equations of conservation of mass, momentum and energy, together with the induction equation, for the case of the purely hydrogen plasma with collisional coupling between the charged and neutral fluid components. The implementation of a semi-implicit algorithm allows us to overcome the numerical stability constraints due to the stiff collisional terms. We test the code against analytical solutions of acoustic and Alfv\'en wave propagation in uniform medium in several regimes of collisional coupling.The results of our simulations are consistent with the analytical estimates, and with other results described in the literature. In the limit of a large collisional frequency, the waves propagate with a common speed of a single fluid. In the other limit of a vanishingly small collisional frequency, the Alfv\'en waves propagate with an Alfv\'en speed of the charged fluid only, while the perturbation in neutral fluid is very small. The acoustic waves in these limits propagate with the sound speed corresponding to either the charges or the neutrals, while the perturbation in the other fluid component is very small. Otherwise, when the collision frequency is similar to the real part of the wave frequency, the interaction between charges and neutrals through momentum transfer collisions cause alterations of the waves frequencies and damping of the wave amplitudes.

Coupling between fast magneto-acoustic and Alfv\'en waves can be observe in fully ionized plasmas mediated by stratification and 3D geometrical effects. In Paper I, Cally & Khomenko (2015) have shown that in a weakly ionized plasma, such as the solar photosphere and chromosphere, the Hall current introduces a new coupling mechanism. The present study extends the results from Paper I to the case of warm plasma. We report on numerical experiments where mode transformation is studied using quasi-realistic stratification in thermodynamic parameters resembling the solar atmosphere. This redresses the limitation of the cold plasma approximation assumed in Paper I, in particular allowing the complete process of coupling between fast and slow magneto-acoustic modes and subsequent coupling of the fast mode to the Alfv\'en mode through the Hall current. Our results confirm the efficacy of the mechanism proposed in Paper I for the solar case. We observe that the efficiency of the transformation is a sensitive function of the angle between the wave propagation direction and the magnetic field, and of the wave frequency. The efficiency increases when the field direction and the wave direction are aligned for increasing wave frequencies. After scaling our results to typical solar values, the maximum amplitude of the transformed Alfv\'en waves, for a frequency of 1 Hz, corresponds to an energy flux (measured above the height of peak Hall coupling) of $\sim10^3$ $\rm W\,m^{-2}$, based on an amplitude of 500 $\rm m\,s^{-1}$ at $\beta=1$, which is sufficient to play a major role in both quiet and active region coronal heating.

Solar active regions show a wide variety of oscillatory phenomena. The presence of the magnetic field leads to the appearance of several wave modes, whose behavior is determined by the sunspot thermal and magnetic structure. We aim to study the relation between the umbral and penumbral waves observed at the high photosphere and the magnetic field topology of the sunspot. Observations of the sunspot in active region NOAA 12662 obtained with the GREGOR telescope (Observatorio del Teide, Spain) were acquired on 2017 June 17. The data set includes a temporal series in the Fe I 5435 \AA\ line obtained with the imaging spectrograph GREGOR Fabry-P\'erot Interferometer (GFPI) and a spectropolarimetric raster map acquired with the GREGOR Infrared Spectrograph (GRIS) in the 10830 \AA\ spectral region. The Doppler velocity deduced from the restored Fe I 5435 \AA\ line has been determined, and the magnetic field vector of the sunspot has been inferred from spectropolarimetric inversions of the Ca I 10839 \AA\ and the Si I 10827 \AA\ lines. A two-armed spiral wavefront has been identified in the evolution of the two-dimensional velocity maps from the Fe I 5435 \AA\ line. The wavefronts initially move counterclockwise in the interior of the umbra, and develop into radially outward propagating running penumbral waves when they reach the umbra-penumbra boundary. The horizontal propagation of the wavefronts approximately follows the direction of the magnetic field, which shows changes in the magnetic twist with height and horizontal position. The spiral wavefronts are interpreted as the visual pattern of slow magnetoacoustic waves which propagate upward along magnetic field lines. Their apparent horizontal propagation is due to their sequential arrival to different horizontal positions at the formation height of the Fe I 5435 \AA\ line, as given by the inclination and orientation of the magnetic field.

Over the last decades, realistic 3D radiative-MHD simulations have become the dominant theoretical tool for understanding the complex interactions between the plasma and the magnetic field on the Sun. Most of such simulations are based on approximations of magnetohydrodynamics, without directly considering the consequences of the very low degree of ionization of the solar plasma in the photosphere and bottom chromosphere. The presence of large amount of neutrals leads to a partial decoupling of the plasma and the magnetic field. As a consequence of that, a series of non-ideal effects (ambipolar diffusion, Hall effect and battery effect) arises. The ambipolar effect is the dominant one in the solar chromosphere. Here we report on the first three-dimensional realistic simulations of magneto-convection including ambipolar diffusion and battery effects. The simulations are done using the newly developed Mancha3D code. Our results reveal that ambipolar diffusion causes measurable effects on the amplitudes of waves excited by convection in the simulations, on the absorption of Poynting flux and heating and on the formation of chromospheric structures. We provide a low limit on the chromospheric temperature increase due to the ambipolar effect using the simulations with battery-excited dynamo fields.

The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere and low chromosphere. Under the single-fluid formalism the main non-ideal effects are described by a series of extra terms in the generalised induction equation and in the energy conservation equation. These effects are: Ohmic diffusion, ambipolar diffusion, the Hall effect, and the Biermann battery effect. From the point of view of the numerical solution of the single-fluid equations, when ambipolar diffusion or Hall effects dominate can introduce severe restrictions on the integration time step and can compromise the stability of the numerical scheme. In this paper we introduce two numerical schemes to overcome those limitations. The first of them is known as Super Time-Stepping (STS) and it is designed to overcome the limitations imposed when the ambipolar diffusion term is dominant. The second scheme is called the Hall Diffusion Scheme (HDS) and it is used when the Hall term becomes dominant. These two numerical techniques can be used together by applying Strang operator splitting. This paper describes the implementation of the STS and HDS schemes in the single-fluid code Mancha3D. The validation for each of these schemes is provided by comparing the analytical solution with the numerical one for a suite of numerical tests.

We investigate the properties of a sunspot light-bridge, focusing on the changes produced by the impact of a plasma blob ejected from a C-class flare. We observed a sunspot in active region NOAA 12544 using spectropolarimetric raster maps of the four Fe I lines around 15655 \AA\ with the GREGOR Infrared Spectrograph (GRIS), narrow-band intensity images sampling the Fe I 6173 \AA\ line with the GREGOR Fabry-P\'erot Interferometer (GFPI), and intensity broad band images in G-band and Ca II H band with the High-resolution Fast Imager (HiFI). All these instruments are located at the GREGOR telescope at the Observatorio del Teide, Tenerife, Spain. The data cover the time before, during, and after the flare event. The analysis is complemented with Atmospheric Imaging Assembly (AIA) and Helioseismic and Magnetic Imager (HMI) data from the Solar Dynamics Observatory (SDO). The physical parameters of the atmosphere at differents heights were inferred using spectral-line inversion techniques. We identify photospheric and chromospheric brightenings, heating events, and changes in the Stokes profiles associated to the flare eruption and the subsequent arrival of the plasma blob to the light bridge, after traveling along an active region loop. The measurements suggest that these phenomena are the result of reconnection events driven by the interaction of the plasma blob with the magnetic field topology of the light bridge.

Partially ionized plasmas are found across the Universe in many different astrophysical environments. They constitute an essential ingredient of the solar atmosphere, molecular clouds, planetary ionospheres and protoplanetary disks, among other environments, and display a richness of physical effects which are not present in fully ionized plasmas. This review provides an overview of the physics of partially ionized plasmas, including recent advances in different astrophysical areas in which partial ionization plays a fundamental role. We outline outstanding observational and theoretical questions and discuss possible directions for future progress.

The magnetic fields of the quiet Sun cover at any time more than 90\% of its surface and their magnetic energy budget is crucial to explain the thermal structure of the solar atmosphere. One of the possible origins of these fields is due to the action of local dynamo in the upper convection zone of the Sun. Existing simulations of the local solar dynamo require an initial seed field, and sufficiently high spatial resolution, in order to achieve the amplification of the seed field to the observed values in the quiet Sun. Here we report an alternative model of seeding based on the action of the Bierman battery effect. This effect generates a magnetic field due to the local imbalances in electron pressure in the partially ionized solar plasma. We show that the battery effect self-consistently creates from zero an initial seed field of a strength of the order of micro G, and together with dynamo amplification, allows the generation of quiet Sun magnetic fields of a similar strength to those from solar observations.

This work aims to understand the behavior of non-linear waves in the vicinity of a coronal null point. In previous works we have showed that high frequency waves are generated in such magnetic configuration. This paper studies those waves in detail in order to provide a plausible explanation of their generation. We demonstrate that slow magneto-acoustic shock waves generated in the chromosphere propagate through the null point and produce a train of secondary shocks that escape along the field lines. A particular combination of the shock wave speeds generates waves at a frequency of 80 mHz. We speculate that this frequency may be sensitive to the atmospheric parameters in the corona and therefore can be used to probe the structure of this solar layer.

Photospheric waves in sunspots moving radially outwards at speeds faster than the characteristic wave velocities have been recently detected. It has been suggested that they are the visual pattern of $p-$modes excited around 5 Mm beneath the sunspot. Using numerical simulations, we have performed a parametric study of the waves observed at the photosphere and higher layers produced by sources located at different depths. The observational measurements are consistent with waves driven between approximately 1 Mm and 5 Mm below the sunspot surface.

Solar photosphere and chromosphere are composed of weakly ionized plasma for which collisional coupling decreases with height. This implies a breakdown of some hypotheses underlying magnetohydrodynamics at low altitudes and gives rise to non-ideal MHD effects such as ambipolar diffusion, Hall effect, etc. Recently, there has been progress in understanding the role of these effects for the dynamics and energetics of the solar atmosphere. There are evidences that such phenomena as wave propagation and damping, magnetic reconnection, formation of stable magnetic field concentrations, magnetic flux emergence, etc. can be affected. This paper reviews the current state-of-the-art of multi-fluid MHD modeling of the coupled solar atmosphere.

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.

We provide estimates of the wave energy reaching the solar chromosphere and corona in a network-like magnetic field topology, including a coronal null point. The waves are excited by an instantaneous strong subphotospheric source and propagate through the subphotosphere, photosphere, chromosphere, transition region, and corona with the plasma beta and other atmospheric parameters varying by several orders of magnitude. We compare two regimes of the wave propagation: a linear and nonlinear regime. While the amount of energy reaching the corona is similar in both regimes, this energy is transmitted at different frequencies. In both cases the dominant periods of waves at each height strongly depend on the local magnetic field topology, but this distribution is only in accordance with observations in the nonlinear case.

We report a detection of differences in ion and neutral velocities in prominences using high resolution spectral data obtained in September 2012 at the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife). A time series of scans of a small portion of a solar prominence was obtained simultaneously with a high cadence using the lines of two elements with different ionization states, namely the CaII 8542 A and the HeI 10830 A. Displacements, widths and amplitudes of both lines were carefully compared to extract dynamical information about the plasma. Many dynamical features are detected, such as counterstreaming flows, jets and propagating waves. In all the cases we find very strong correlation between the parameters extracted from the lines of both elements, confirming that both trace the same plasma. Nevertheless, we also find short-lived transients where this correlation is lost. These transients are associated with the ion-neutral drift velocities of the order of several hundred m/s. The patches of non-zero drift velocity show coherence on time-distance diagrams.

This paper studies the dependence of the CaIIH line core brightness on the strength and inclination of photospheric magnetic field, and on the parameters of convective and wave motions in a facular region at the solar disc center. We use three simultaneous datasets obtained at the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife): (1) spectra of BaII 4554 A line registered with the instrument TESOS to measure the variations of intensity and velocity through the photosphere up to the temperature minimum; (2) spectropolarimetric data in FeI 1.56 $\mu$m lines (registered with the instrument TIP II) to measure photospheric magnetic fields; (3) filtergrams in CaIIH that give information about brightness fluctuations in the chromosphere. The results show that the CaIIH brightness in the facula strongly depends on the power of waves with periods in the 5-min range, that propagate upwards, and also on the phase shift between velocity oscillations at the bottom photosphere and around the temperature minimum height, measured from BaII line. The CaIIH brightness is maximum at locations where the phase shift between temperature and velocity oscillations lies within 0-100 degrees. There is an indirect influence of convective motions on the CaIIH brightness. Namely, the higher is the amplitude of convective velocities and the larger is the height where they change their direction of motion, the brighter is the facula. Altogether, our results lead to conclusions that facular regions appear bright not only because of the Wilson depression in magnetic structures, but also due to real heating.

In this paper, we show a "proof of concept" of the heating mechanism of the solar chromosphere due to wave dissipation caused by the effects of partial ionization. Numerical modeling of non-linear wave propagation in a magnetic flux tube, embedded in the solar atmosphere, is performed by solving a system of single-fluid quasi-MHD equations, which take into account the ambipolar term from the generalized Ohm's law. It is shown that perturbations caused by magnetic waves can be effectively dissipated due to ambipolar diffusion. The energy input by this mechanism is continuous and shown to be more efficient than dissipation of static currents, ultimately leading to chromospheric temperature increase in magnetic structures.

Large-amplitude longitudinal oscillations (LALOs) in prominences are spectacular manifestations of the solar activity. In such events nearby energetic disturbances induce periodic motions on filaments with displacements comparable to the size of the filaments themselves and with velocities larger than 20 km/s. The pendulum model, in which the gravity projected along a rigid magnetic field is the restoring force, was proposed to explain these events. However, it can be objected that in a realistic situation where the magnetic field reacts to the mass motion of the heavy prominence, the simplified pendulum model could be no longer valid. We have performed non-linear time-dependent numerical simulations of LALOs considering a dipped magnetic field line structure. In this work we demonstrate that for even relatively weak magnetic fields the pendulum model works very well. We therefore validate the pendulum model and show its robustness, with important implications for prominence seismology purposes. With this model it is possible to infer the geometry of the dipped field lines that support the prominence.

The photospheric temperature minimum in the Sun and solar-like stars is very weakly ionized, with ionization fraction $f$ as low as $10^{-4}$. In galactic star forming regions, $f$ can be $10^{-10}$ or lower. Under these circumstances, the Hall current can couple low frequency Alfv\'en and magneto\-acoustic waves via the dimensionless Hall parameter $\epsilon=\omega/\Omega_\text{i}f$, where $\omega$ is the wave frequency and $\Omega_\text{i}$ is the mean ion gyrofrequency. This is analysed in the context of a cold (zero-$\beta$) plasma, and in less detail for a warm plasma. It is found that Hall coupling preferentially occurs where the wave vector is nearly field-aligned. In these circumstances, Hall coupling in theory produces a continual oscillation between fast and Alfv\'en modes as the wave passes through the weakly ionized region. At low frequencies (mHz), characteristic of solar and stellar normal modes, $\epsilon$ is probably too small for more than a fraction of one oscillation to occur. On the other hand, the effect may be significant at the far higher frequencies (Hz) associated with magnetic reconnection events. In another context, characteristic parameters for star forming gas clouds suggest that $\mathcal{O}(1)$ or more full oscillations may occur in one cloud crossing. This mechanism is not expected to be effective in sunspots, due to their high ion gyrofrequencies and Alfv\'en speeds, since the net effect depends inversely on both and therefore inverse quadratically on field strength.

Recent high-resolution observations of sunspot oscillations using simultaneously operated ground- and space-based telescopes reveal the intrinsic connection between different layers of the solar atmosphere. However, it is not clear whether these oscillations are externally driven or generated in-situ. We address this question by using observations of propagating slow magneto-acoustic waves along a coronal fan loop system. In addition to the generally observed decreases in oscillation amplitudes with distance, the observed wave amplitudes are also found to be modulated with time, with similar variations observed throughout the propagation path of the wavetrain. Employing multi-wavelength and multi-instrument data we study the amplitude variations with time as the waves propagate through different layers of the solar atmosphere. By comparing the amplitude-modulation period in different layers, we find that slow magneto-acoustic waves observed in sunspots are externally driven by photospheric p-modes, which propagate upwards into the corona before becoming dissipated.

The amplitudes of the Evershed flow are measured using pairs of carefully selected FeI and FeII spectral lines located close in wavelength and registered simultaneously. A sunspot belonging to the NOAA 11582 group was scanned using the spectrograph of the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife). Velocities were extracted from intensity profiles using the lambda-meter technique. The formation heights of the observed spectral lines were calculated using semi-empirical models of a bright and dark penumbral filament taking into account the sunspot location at the limb. Our objective is to compare azimuthally averaged amplitudes of the Evershed flow extracted from neutral and ion lines. We find measurable differences in the radial component of the flow. All five pairs of lines show the same tendency, with a few hundred m/s larger amplitude of the flow measured from FeI lines compared to FeII lines. This tendency is preserved at all photospheric heights and radial distances in the penumbra. We discuss the possible origin of this effect.

Coronal waves are an important aspect of the dynamics of the plasma in the corona. Wavefront dislocations are topological features of most waves in nature and also of magnetohydrodynamic waves. Are there dislocations in coronal waves? The finding and explanation of dislocations may shed light on the nature and characteristics of the propagating waves, their interaction in the corona and in general on the plasma dynamics. We positively identify dislocations in coronal waves observed by the Coronal Multi-channel Polarimeter (CoMP) as singularities in the Doppler shifts of emission coronal lines. We study the possible singularities that can be expected in coronal waves and try to reproduce the observed dislocations in terms of localization and frequency of appearance. The observed dislocations can only be explained by the interference of a kink and a sausage wave modes propagating with different frequencies along the coronal magnetic field. In the plane transverse to the propagation, the cross-section of the oscillating plasma must be smaller than the spatial resolution, and the two waves result in net longitudinal and transverse velocity components that are mixed through projection onto the line of sight. Alfv\'en waves can be responsible of the kink mode, but a magnetoacoustic sausage mode is necessary in all cases. Higher (flute) modes are excluded. The kink mode has a pressure amplitude that is smaller than the pressure amplitude of the sausage mode, though its observed velocity is larger. This concentrates dislocations on the top of the loop. To explain dislocations, any model of coronal waves must include the simultaneous propagation and interference of kink and sausage wave modes of comparable but different frequencies, with a sausage wave amplitude much smaller than the kink one.

The temperature and density conditions in the magnetized photosphere and chromosphere of the Sun lead to a very small degree of atomic ionization. In addition, at particular height, the magnetic field may be strong enough to give rise to a cyclotron frequency larger than the collisional frequency for some species, while for others the opposite may happen. These circumstances influence the collective behavior of the particles and some of the hypotheses of magnetohydrodynamics may be relaxed, giving rise to non-ideal MHD effects. In this paper we discuss our recent developments in modeling non-ideal plasma effects derived from the presence of a large amount of neutrals in the solar photosphere and the chromosphere, as well as observational consequences of these effects.

The aim of this work is to study the energy transport by means of MHD waves propagating in quiet Sun magnetic topology from layers below the surface to the corona. Upward propagating waves find obstacles, such as the equipartition layer with plasma b=1 and the transition region, and get converted, reflected and refracted. Understanding the mechanisms by which MHD waves can reach the corona can give us information about the solar atmosphere and the magnetic structures. We carry out two-dimensional numerical simulations of wave propagation in a magnetic field structure that consists of two vertical flux tubes separated by an arcade shaped magnetic field. This configuration contains a null point in the corona, that significantly modifies the behaviour of the waves. We describe in detail the wave propagation through the atmosphere under different driving conditions. We also present the spatial distribution of the mean acoustic and magnetic energy fluxes and the spatial distribution of the dominant frequencies in the whole domain. We conclude that the energy reaches the corona preferably along vertical magnetic fields, inside the flux tubes, and it has an acoustic nature. Most of the magnetic energy keeps concentrated below the transition region due to the refraction of the magnetic waves and the continuous conversion of acoustic-like waves into fast magnetic waves in the equipartition layer located in the photosphere. However, part of the magnetic energy reaches the low corona when propagating in the region where the arcades are located, but waves are sent back downwards to the lower atmosphere at the null point surroundings. This phenomenon, together with the reflection and refraction of waves in the TR and the lower turning point, act as a re-feeding of the atmosphere. In the frequency distribution, we find that high frequency waves can reach the corona outside the vertical flux tubes.

Spectropolarimetric temporal series from Fe I $\lambda$ 6301.5 \AA\ and Ca II infrared triplet lines are obtained by applying the Stokes synthesis code NICOLE to a numerical simulation of wave propagation in a sunspot umbra from MANCHA code. The analysis of the phase difference between Doppler velocity and intensity core oscillations of the Fe I $\lambda$ 6301.5 \AA\ line reveals that variations in the intensity are produced by opacity fluctuations rather than intrinsic temperature oscillations, except for frequencies between 5 and 6.5 mHz. On the other hand, the photospheric magnetic field retrieved from the weak field approximation provides the intrinsic magnetic field oscillations associated to wave propagation. Our results suggest that this is due to the low magnetic field gradient of our sunspot model. The Stokes parameters of the chromospheric Ca II infrared triplet lines show striking variations as shock waves travel through the formation height of the lines, including emission self-reversals in the line core and highly abnormal Stokes V profiles. Magnetic field oscillations inferred from the Ca II infrared lines using the weak field approximation appear to be related with the magnetic field strength variation between the photosphere and the chromosphere.

We derive self-consistent formalism for the description of multi-component partially ionized solar plasma, by means of the coupled equations for the charged and neutral components for an arbitrary number of chemical species, and the radiation field. All approximations and assumptions are carefully considered. Generalized Ohm's law is derived for the single-fluid and two-fluid formalism. Our approach is analytical with some order-of-magnitude support calculations. After general equations are developed we particularize to some frequently considered cases as for the interaction of matter and radiation.

We study the Rayleigh-Taylor instability (RTI) at a prominence-corona transition region in a non-linear regime. Our aim is to understand how the presence of neutral atoms in the prominence plasma influences the instability growth rate, and the evolution of velocity, magnetic field vector and thermodynamic parameters of turbulent drops. We perform 2.5D numerical simulations of the instability initiated by a multi-mode perturbation at the corona-prominence interface using a single-fluid MHD approach including a generalized Ohm's law. The initial equilibrium configuration is purely hydrostatic and contains a homogeneous horizontal magnetic field forming an angle with the direction in which the plasma is perturbed. We analyze simulations with two different orientations of the magnetic field. For each field orientation we compare two simulations, one for the pure MHD case, and one including the ambipolar diffusion in the Ohm's law (AD case). Other than that, both simulations for each field orientation are identical. The numerical results in the initial stage of the instability are compared with the analytical linear calculations. We find that the configuration is always unstable in the AD case. The growth rate of the small-scale modes in the non-linear regime is up to 50% larger in the AD case than in the purely MHD case and the average velocities of flows are a few percent larger. Significant drift momenta are found at the interface between the coronal and the prominence material at all stages of the instability, produced by the faster downward motion of the neutral component with respect to the ionized component. The differences in temperature of the bubbles between the ideal and non-ideal case are also significant, reaching 30%. There is an asymmetry between large rising bubbles and small-scale down flowing fingers, favoring the detection of upward velocities in observations.

We study the modification of the classical criterion for the linear onset and growth rate of the Rayleigh-Taylor instability (RTI) in a partially ionized (PI) plasma in the one-fluid description, considering a generalized induction equation. The governing linear equations and appropriate boundary conditions, including gravitational terms, are derived and applied to the case of the RTI in a single interface between two partially ionized plasmas. The boundary conditions lead to an equation for the frequencies in which some of them have positive complex parts, marking the appearance of the RTI. We study the ambipolar term alone first, extending the result to the full induction equation later. We find that the configuration is always unstable because of the presence of a neutral species. In the classical stability regime the growth rate is small, since the collisions prevent the neutral fluid to fully develop the RTI. For parameters in the classical instability regime the growth rate is lowered, but for the considered theoretical values of the collision frequencies and diffusion coefficients for solar prominences the differences with the compressible MHD case are small. We conclude that PI modifies some aspects of the linear RTI instability, since it takes into account that neutrals do not feel the stabilizing effect of the magnetic field. For the set of parameters representative for solar prominences, our model gives the resulting timescale comparable with observed lifetimes of RTI plumes.

We study Rayleigh-Taylor instability (RTI) at the coronal-prominence boundary by means of 2.5D numerical simulations in a single-fluid MHD approach including a generalized Ohm's law. The initial configuration includes a homogeneous magnetic field forming an angle with the direction in which the plasma is perturbed. For each field inclination we compare two simulations, one for the pure MHD case, and one including the ambipolar diffusion in the Ohm's law, otherwise identical. We find that the configuration containing neutral atoms is always unstable. The growth rate of the small-scale modes in the non-linear regime is larger than in the purely MHD case.

We study the properties of waves in a facular region of moderate strength in the photosphere and chromosphere. Our aim is to analyse statistically the wave periods, power and phase relations as a function of the magnetic field strength and inclination. Our work is based on observations obtained at the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife) using two different instruments: the Triple Etalon SOlar Spectrometer (TESOS), in the BaII 4554 A line to measure velocity and intensity variations through the photosphere; and, simultaneously, the Tenerife Infrared Polarimeter (TIP-II), in the FeI 1.56 mm lines to the measure the Stokes parameters and magnetic field strength in the lower photosphere. Additionally, we use the simultaneous broad-band filtergrams in the CaIIH line to obtain information about intensity oscillations in the chromosphere. We find several clear trends in the oscillation behaviour: (i) the period of oscillation increases by 15-20 % with the magnetic field increasing from 500 to 1500 G; (ii) the temperature-velocity phase shifts show a strikingly different distribution in the facular region compared to the quiet region, a significant number of cases in the range from -180 to 180 degrees is detected in the facula. (iii) the most powerful chromospheric CaIIH intensity oscillations are observed at locations with strong magnetic fields (1.3-1.5 kG) inclined by 10-12 degrees, as a result of upward propagating waves with rather small phase speeds, and temperature-velocity phase shifts between 0 and 90 degrees; (iv) the power of the photospheric velocity oscillations from the \BaII\ line increases linearly with decreasing magnetic field inclination, reaching its maximum at strong field locations.

We describe the presence of wavefront dislocations in magnetohydrodynamic waves in stratified stellar atmospheres. Scalar dislocations such as edges and vortices can appear in Alfv\'en waves, as well as in general magneto-acoustic waves. We detect those dislocations in observations of magnetohydrodynamic waves in sunspots in the solar chromosphere. Through the measured charge of all the dislocations observed, we can give for the first time estimates of the modal contribution in the waves propagating along magnetic fields in solar sunspots.

Waves are observed at all layers of the solar atmosphere and the magnetic field plays a key role in their propagation. While deep down in the atmosphere the p-modes are almost entirely of acoustic nature, in the upper layers magnetic forces are dominating, leading to a large variety of new wave modes. Significant advances have been made recently in our understanding of the physics of waves interaction with magnetic structures, with the help of analytical theories, numerical simulations, as well as high-resolution observations. In this contribution, we review recent observational findings and current theoretical ideas in the field, with an emphasis on the following questions: (i) Peculiarities of the observed wave propagation in network, plage and facular regions; (ii) Role of the mode transformation and observational evidences of this process; (iii) Coupling of the photosphere, chromosphere, and above by means of waves propagating in magnetic structures.

We present results of the study of chromospheric and photospheric line-of-sight velocity fields in the young active region NOAA 11024. Multi-layer, multi-wavelength observational data were used for the analysis of the emerging flux in this active region. Spectropolarimetric observations were carried out with the telescope THEMIS on Tenerife (Canary Islands) on 4 July 2009. In addition, space-borne data from SOHO/MDI, STEREO and GOES were also considered. The combination of data from ground- and space-based telescopes allowed us to study the dynamics of the lower atmosphere of the active region with high spatial, spectral, and temporal resolutions. THEMIS spectra show strong temporal variations of the velocity in the chromosphere and photosphere for different activity features: two pores, active and quiet plage regions, and two surges. The range of variations of the chromospheric line-of-sight velocity at the heights of formation of the H-alpha core was extremely large. Both upward and downward motions were observed in these layers. In particular, a surge with upward velocities up to -73 km/s were detected. In the photosphere, predominantly upward motions were found, varying from -3.1 km/s upflows to 1.4 km/s downflows in different structures. The velocity variations at different levels in the lower atmosphere are compatible with magnetic flux emergence.

We investigate the impact on Fe abundance determination of including magnetic flux in series of 3D radiation-MHD simulations of solar convection which we used to synthesize spectral intensity profiles corresponding to disc centre. A differential approach is used to quantify the changes in theoretical equivalent width of a set of 28 iron spectral lines spanning a wide range in lambda, excitation potential, oscillator strength, Land\'e factor, and formation height. The lines were computed in LTE using the spectral synthesis code LILIA. We used input magnetoconvection snapshots covering 50 minutes of solar evolution and belonging to series having an average vertical magnetic flux density of 0, 50, 100 and 200 G. For the relevant calculations we used the Copenhagen Stagger code. The presence of magnetic fields causes both a direct (Zeeman-broadening) effect on spectral lines with non-zero Land\'e factor and an indirect effect on temperature-sensitive lines via a change in the photospheric T-tau stratification. The corresponding correction in the estimated atomic abundance ranges from a few hundredths of a dex up to |Delta log(Fe)| ~ 0.15 dex, depending on the spectral line and on the amount of average magnetic flux within the range of values we considered. The Zeeman-broadening effect gains relatively more importance in the IR. The largest modification to previous solar abundance determinations based on visible spectral lines is instead due to the indirect effect, i.e., the line-weakening caused by a warmer stratification on an optical depth scale. Our results indicate that the average solar iron abundance obtained when using magnetoconvection models can be 0.03-0.11 dex higher than when using the simpler HD convection approach. We demonstrate that accounting for magnetic flux is important in state-of-the-art solar photospheric abundance determinations based on 3D simulations.

In this paper, we study the properties of solar granulation in a facular region from the photosphere up to the lower chromosphere. Our aim is to investigate the dependence of granular structure on magnetic field strength. We use observations obtained at the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife) using two different instruments: Triple Etalon SOlar Spectrometer (TESOS), in the BaII 4554 A line to measure velocity and intensity variations along the photosphere; and, simultaneously, Tenerife Infrared Polarimeter (TIP-II), in the FeI 1.56 $\mu$m lines to the measure Stokes parameters and the magnetic field strength at the lower photosphere. We obtain that the convective velocities of granules in the facular area decrease with magnetic field while the convective velocities of intergranular lanes increase with the field strength. Similar to the quiet areas, there is a contrast and velocity sign reversal taking place in the middle photosphere. The reversal heights depend on the magnetic field strength and are, on average, about 100 km higher than in the quiet regions. The correlation between convective velocity and intensity decreases with magnetic field at the bottom photosphere, but increases in the upper photosphere. The contrast of intergranular lanes observed close to the disc center is almost independent of the magnetic field strength. The strong magnetic field of facular area seems to stabilize the convection and to promote more effective energy transfer in the upper layers of the solar atmosphere, since the convective elements reach larger heights.

We propose a mechanism for efficient heating of the solar chromosphere, based on non-ideal plasma effects. Three ingredients are needed for the work of this mechanism: (1) presence of neutral atoms; (2) presence of a non-potential magnetic field; (3) decrease of the collisional coupling of the plasma. Due to decrease of collisional coupling, a net relative motion appears between the neutral and ionized components, usually referred to as "ambipolar diffusion". This results in a significant enhancement of current dissipation as compared to the classical MHD case. We propose that the current dissipation in this situation is able to provide enough energy to heat the chromosphere by several kK on the time scale of minutes, or even seconds. In this paper, we show that this energy supply might be sufficient to balance the radiative energy losses of the chromosphere.

In this paper, we study the heating of the magnetized solar chromosphere induced by the large fraction of neutral atoms present in this layer. The presence of neutrals, together with the decrease with height of the collisional coupling, leads to deviations from the classical MHD behavior of the chromospheric plasma. A relative net motion appears between the neutral and ionized components, usually referred to as ambipolar diffusion. The dissipation of currents in the chromosphere is enhanced orders of magnitude due to the action of ambipolar diffusion, as compared to the standard ohmic diffusion. We propose that a significant amount of magnetic energy can be released to the chromosphere just by existing force-free 10--40 G magnetic fields there. As a consequence, we conclude that ambipolar diffusion is an important process that should be included in chromospheric heating models, as it has the potential to rapidly heat the chromosphere. We perform analytical estimations and numerical simulations to prove this idea.

We study the conversion of fast magneto-acoustic waves to Alfven waves by means of 2.5D numerical simulations in a sunspot-like magnetic configuration. A fast, essentially acoustic, wave of a given frequency and wave number is generated below the surface and propagates upward though the Alfven/acoustic equipartition layer where it splits into upgoing slow (acoustic) and fast (magnetic) waves. The fast wave quickly reflects off the steep Alfven speed gradient, but around and above this reflection height it partially converts to Alfven waves, depending on the local relative inclinations of the background magnetic field and the wavevector. To measure the efficiency of this conversion to Alfven waves we calculate acoustic and magnetic energy fluxes. The particular amplitude and phase relations between the magnetic field and velocity oscillations help us to demonstrate that the waves produced are indeed Alfven waves. We find that the conversion to Alfven waves is particularly important for strongly inclined fields like those existing in sunspot penumbrae. Equally important is the magnetic field orientation with respect to the vertical plane of wave propagation, which we refer to as "field azimuth". For field azimuth less than 90 degrees the generated Alfven waves continue upwards, but above 90 degrees downgoing Alfven waves are preferentially produced. This yields negative Alfven energy flux for azimuths between 90 and 180 degrees. Alfven energy fluxes may be comparable to or exceed acoustic fluxes, depending upon geometry, though computational exigencies limit their magnitude in our simulations.

Full Stokes filter-polarimeters are key instruments for investigating the rapid evolution of magnetic structures on the solar surface. To this end, the image quality is routinely improved using a-posteriori image reconstruction methods. We analyze the robustness of circular polarization asymmetries to phase-diversity image reconstruction techniques. We use snapshots of magneto-hydrodynamical simulations carried out with different initial conditions to synthesize spectra of the magnetically sensitive Fe I line at 5250.2 A. We degrade the synthetic profiles spatially and spectrally to simulate observations with the IMaX full Stokes filter-polarimeter. We also simulate the focused/defocused pairs of images used by the phase-diversity algorithm for reconstruction and the polarimetric modulation scheme. We assume that standard optimization methods are able to infer the projection of the wavefront on the Zernike polynomials with 10% precision. We also consider the less favorable case of 25% precision. We obtain reconstructed monochromatic modulated images that are later demodulated and compared with the original maps. Although asymmetries are often difficult to define in the quiet Sun due to the complexity of the Stokes V profiles, we show how asymmetries are degraded with spatial and spectral smearing. The results indicate that, although image reconstruction techniques reduce the spatial smearing, they can modify the asymmetries of the profiles, mainly caused by the appearance of spatially-correlated noise.

We aim at reproducing the height dependence of sunspot wave signatures obtained from spectropolarimetric observations through 3D MHD numerical simulations. A magneto-static sunspot model based on the properties of the observed sunspot is constructed and perturbed at the photosphere introducing the fluctuations measured with the \SiI\ $\lambda$ 10827 \AA\ line. The results of the simulations are compared with the oscillations observed simultaneously at different heights from the \HeI\ $\lambda$ 10830 \AA\ line, the \CaIIH\ core and the \FeI\ blends in the wings of the \CaIIH\ line. The simulations show a remarkable agreement with the observations. They reproduce the velocity maps and power spectra at the formation heights of the observed lines, as well as the phase and amplification spectra between several pair of lines. We find that the stronger shocks at the chromosphere are accompanied with a delay between the observed signal and the simulated one at the corresponding height, indicating that shocks shift the formation height of the chromospheric lines to higher layers. Since the simulated wave propagation matches very well the properties of the observed one, we are able to use the numerical calculations to quantify the energy contribution of the magneto-acoustic waves to the chromospheric heating in sunspots. Our findings indicate that the energy supplied by these waves is too low to balance the chromospheric radiative losses. The energy contained at the formation height of the lowermost \SiI\ $\lambda$ 10827 \AA\ line in the form of slow magneto-acoustic waves is already insufficient to heat the higher layers, and the acoustic energy which reaches the chromosphere is around 3-9 times lower than the required amount of energy. The contribution of the magnetic energy is even lower.

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.

We explore the effect of the magnetic field when using realistic three-dimensional convection experiments to determine solar element abundances. By carrying out magnetoconvection simulations with a radiation-hydro code (the Copenhagen stagger code) and through a-posteriori spectral synthesis of three Fe I lines, we obtain evidence that moderate amounts of mean magnetic flux cause a noticeable change in the derived equivalent widths compared with those for a non-magnetic case. The corresponding Fe abundance correction for a mean flux density of 200 G reaches up to ~0.1 dex in magnitude. These results are based on space- and time-averaged line profiles over a time span of 2.5 solar hours in the statistically stationary regime of the convection. The main factors causing the change in equivalent widths, namely the Zeeman broadening and the modification of the temperature stratification, act in different amounts and, for the iron lines considered here, in opposite directions; yet, the resulting $|\Delta\log\epsilon_{\odot}(Fe)|$ coincides within a factor two in all of them, even though the sign of the total abundance correction is different for the visible and infrared lines. We conclude that magnetic effects should be taken into account when discussing precise values of the solar and stellar abundances and that an extended study is warranted.

We study the propagation of waves from the photosphere to the chromosphere of sunspots. From time series of cospatial Ca II H (including its line blends) intensity spectra and polarimetric spectra of Si I 1082.7 nm and He I 1083.0 nm we retrieve the line-of-sight velocity at several heights. The analysis of the phase difference and amplification spectra shows standing waves for frequencies below 4 mHz and propagating waves for higher frequencies, and allows us to infer the temperature and height where the lines are formed. Using these observational data, we have constructed a model of sunspot, and we have introduced the velocity measured with the photospheric Si I 1082.7 nm line as a driver. The numerically propagated wave pattern fits reasonably well with the observed using the lines formed at higher layers, and the simulations reproduce many of the observed features. The observed waves are slow MHD waves propagating longitudinally along field lines.

We study the coupling of magneto-acoustic waves to Alven waves using 2.5D numerical simulations. In our experiment, a fast magnetoacoustic wave of a given frequency and wavenumber is generated below the surface. The magnetic field in the domain is assumed homogeneous and inclined. The efficiency of the conversion to Alfven waves near the layer of equal acoustic and Alfven speeds is measured calculating their energy flux. The particular amplitude and phase relations between the oscillations of magnetic field and velocity help us to demonstrate that the waves produced after the transformation and reaching upper atmosphere are indeed Alfven waves. We find that the conversion from fast magneto-acoustic waves to Alfven waves is particularly important for the inclination and azimuth angles of the magnetic field between 55 and 65 degrees, with the maximum shifted to larger inclinations for lower frequency waves. The maximum Alfven flux transmitted to the upper atmosphere is about 2-3 times lower than the corresponding acoustic flux.

We analyze the propagation of waves in sunspots from the photosphere to the chromosphere using time series of co-spatial Ca II H intensity spectra (including its line blends) and polarimetric spectra of Si I 10827 and the He I 10830 multiplet. From the Doppler shifts of these lines we retrieve the variation of the velocity along the line-of-sight at several heights. Phase spectra are used to obtain the relation between the oscillatory signals. Our analysis reveals standing waves at frequencies lower than 4 mHz and a continuous propagation of waves at higher frequencies, which steepen into shocks in the chromosphere when approaching the formation height of the Ca II H core. The observed non-linearities are weaker in Ca II H than in He I lines. Our analysis suggests that the Ca II H core forms at a lower height than the He I 10830 line: a time delay of about 20 s is measured between the Doppler signal detected at both wavelengths. We fit a model of linear slow magnetoacoustic wave propagation in a stratified atmosphere with radiative losses according to Newton's cooling law to the phase spectra and derive the difference in the formation height of the spectral lines. We show that the linear model describes well the wave propagation up to the formation height of Ca II H, where non-linearities start to become very important.

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.

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.

Waves observed in the photosphere and chromosphere of sunspots show complex dynamics and spatial patterns. The interpretation of high-resolution sunspot wave observations requires modeling of three-dimensional non-linear wave propagation and mode transformation in the sunspot upper layers in realistic spot model atmospheres. Here we present the first results of such modeling. We have developed a 3D non-linear numerical code specially designed to calculate the response of magnetic structures in equilibrium to an arbitrary perturbation. The code solves the 3D nonlinear MHD equations for perturbations; it is stabilized by hyper-diffusivity terms and is fully parallelized. The robustness of the code is demonstrated by a number of standard tests. We analyze several simulations of a sunspot perturbed by pulses of different periods at subphotospheric level, from short periods, introduced for academic purposes, to longer and realistic periods of three and five minutes. We present a detailed description of the three-dimensional mode transformation in a non-trivial sunspot-like magnetic field configuration, including the conversion between fast and slow magneto-acoustic waves and the Alfv\'en wave, by calculation of the wave energy fluxes. Our main findings are the following: (1) the conversion from acoustic to the Alfv\'en mode is only observed if the the driving pulse is located out of the sunspot axis, but this conversion is energetically inefficient; (2) as a consequence of the cut-off effects and refraction of the fast magneto-acoustic mode, the energy of the evanescent waves with periods around 5 minutes remains almost completely below the level beta=1; (3) waves with frequencies above the cut-off propagate field-aligned to the chromosphere and their power becomes dominating over that of evanescent 5-minute oscillations, in agreement with observations.

Three-dimensional numerical simulations of magnetoacoustic wave propagation are performed in a sunspot atmosphere with a computational domain covering from the photosphere to the chromosphere. The wave source, with properties resembling the solar spectrum, is located at different distances from the axis of the sunspot for each simulation. These results are compared with the theory of mode transformation and also with observational features. Simulations show that the dominant oscillation frequency in the chromosphere decreases with the radial distance from the sunspot axis. The energy flux of the different wave modes involved, including de Alfv\'en mode, is evaluated and discussed.

Time-resolved spectroscopic observations of rapidly oscillating Ap (roAp) stars show a complex picture of propagating magneto-acoustic pulsation waves, with amplitude and phase strongly changing as a function of atmospheric height. We have recently conducted numerical, non-linear MHD simulations to get an insight into the complex atmospheric dynamics of magnetic pulsators. Here we use the resulting time-dependent atmospheric structure and velocity field to predict line profile variations for roAp stars. These calculations use realistic atmospheric structure, account for vertical chemical stratification and treat the line formation in pulsating stellar atmosphere without relying on the simplistic single-layer approximation universally adopted for non-radial pulsators. The new theoretical calculations provide an essential tool for interpreting the puzzling complexity of the spectroscopic pulsations in roAp stars.

We present comparison of numerical simulations of propagation of MHD waves,excited by subphotospheric perturbations, in two different ("deep" and "shallow") magnetostatic models of the sunspots. The "deep" sunspot model distorts both the shape of the wavefront and its amplitude stronger than the "shallow" model. For both sunspot models, the surface gravity waves (f-mode) are affected by the sunspots stronger than the acoustic p-modes. The wave amplitude inside the sunspot depends on the photospheric strength of the magnetic field and the distance of the source from the sunspot axis. For the source located at 9 Mm from the center of the sunspot, the wave amplitude increases when the wavefront passes through the central part of the sunspot. For the source distance of 12 Mm, the wave amplitude inside the sunspot is always smaller than outside. For the same source distance from the sunspot center but for the models with different strength of the magnetic field, the wave amplitude inside the sunspot increases with the strength of the magnetic field. The simulations show that unlike the case of the uniform inclined background magnetic field, the p- and f-mode waves are not spatially separated inside the sunspot where the magnetic field is strongly non-uniform. These properties have to be taken into account for interpretation of observations of MHD waves traveling through sunspot regions.

We suggest an explanation for the high-frequency power excess surrounding active regions known as seismic halos. The idea is based on numerical simulations of magneto-acoustic waves propagation in sunspots. We propose that such an excess can be caused by the additional energy injected by fast mode waves refracted in the higher atmosphere due to the rapid increase of the Alfven speed. Our model qualitatively explains the magnitude of the halo and allows to make some predictions of its behavior that can be checked in future observations.

The aim of this paper is to analyse the validity of the Dopplergram and lambda-meter techniques for the Doppler diagnostics of solar photospheric velocities using the BaII 4554 A line. Both techniques are evaluated by means of NLTE radiative transfer calculations of the BaII 4554 A line in a three-dimensional hydrodynamical model of solar convection. We consider the cases of spatially unsmeared profiles and the profiles smeared to the resolution of ground-based observations. We find that: (i) Speckle-reconstructed Dopplergram velocities reproduce the ``true'' velocities well at heights around 300 km, except for intergranular lanes with strong downflows where the velocity can be overestimated. (ii) The lambda-meter velocities give a good representation of the ``true'' velocities through the whole photosphere, both under the original and reduced spatial resolutions. The velocities derived from the inner wing of smeared BaII 4554 A line profiles are more reliable than those for the outer wing. Only under high spatial resolution does the inner wing velocities calculated in intergranular regions give an underestimate (or even a sign reversal) compared with the model velocities. (iii) NLTE effects should be taken into account in modelling the BaII 4554 A line profiles. Such effects are more pronounced in intergranular regions. Our analysis supports the opinion that the Dopplergram technique applied to the BaII 4554 A line is a valuable tool for the Doppler diagnostics of the middle photosphere around 300 km. The \lambda-meter technique applied to this line gives us a good opportunity to ``trace'' the non-thermal motions along the whole photosphere up to the temperature minimum and lower chromosphere.

Rapidly oscillating Ap stars exhibit an astrophysically interesting combination of strong, dipolar-like magnetic fields and high-overtone p-mode pulsations similar to the Sun. Recent time-resolved spectroscopy of these stars unravelled a complex picture of propagating magneto-acoustic pulsation waves, with amplitude and phase strongly changing as a function of atmospheric height. To interpret these observations and gain a new insight into the atmospheric dynamics of roAp stars we have carried out 2-D time-dependent, non-linear magneto-hydrodynamical simulations of waves for a realistic atmospheric stratification of a cool Ap star. We explore a grid of simulations in a wide parameter space, treating oscillations of the velocity, magnetic field and thermodynamic quantities in a self-consistent manner. Our simulations foster a new understanding of the influence of the atmosphere and the magnetic field on the propagation and reflection properties of magneto-acoustic waves, formation of node surfaces, and relative variation of different quantities. Our simulations reproduce all main features of the observed pulsational behavior of roAp stars. We show, for the first time, that the overall depth dependence of the pulsations in roAp atmospheres is strongly influenced by the density inversion at the photospheric base.

The purpose of this paper is to characterize the statistical properties of solar granulation in the photosphere and low chromosphere up to 650 km. We use velocity and intensity variations obtained at different atmospheric heights from observations in BaII 4554 A. The observations were done during good seeing conditions at the VTT at the Observatorio del Teide on Tenerife. The line core forms rather high in the atmosphere and allows granulation properties to be studied at heights that have been not accessed before in similar studies. In addition, we analyze the synthetic profiles of the BaII 4554 A line by the same method computed taking NLTE effects into account in the 3D hydrodynamical model atmosphere. We suggest a 16-column model of solar granulation depending on the direction of motion and on the intensity contrast measured in the continuum and in the uppermost layer. We calculate the heights of intensity contrast sign reversal and velocity sign reversal. We show that both parameters depend strongly on the granulation velocity and intensity at the bottom photosphere. The larger the two parameters, the higher the reversal takes place in the atmosphere. On average, this happens at about 200-300 km. We suggest that this number also depends on the line depth of the spectral line used in observations. Despite the intensity and velocity reversal, about 40% of the column structure of granulation is preserved up to heights around 650 km.

Acoustic waves are one of the primary suspects besides magnetic fields for the chromospheric heating process to temperatures above radiative equilibrium (RE). We derived the mechanical wave energy as seen in line-core velocities to obtain a measure of mechanical energy flux with height for a comparison with the energy requirements in a semi-empirical atmosphere model. We analyzed a 1-hour time series and a large-area map of Ca II H spectra on the traces of propagating waves. We analyzed the velocity statistics of several spectral lines in the wing of Ca II H, and the line-core velocity of Ca II H. We converted the velocity amplitudes into volume and mass energy densities. For comparison, we used the increase of internal energy necessary to lift a RE atmosphere to the HSRA temperature stratification. We find that the velocity amplitude grows in agreement with linear wave theory and thus slower with height than predicted from energy conservation. The mechanical energy of the waves above around z~500 km is insufficient to maintain the chromospheric temperature rise in the semi-empirical HSRA model. The intensity variations of the Ca line core (z~1000 km) can be traced back to the velocity variations of the lowermost forming spectral line considered (z~ 250 km). The chromospheric intensity, and hence, (radiation) temperature variations are seen to be induced by passing waves originating in the photosphere.

We report 2D time-dependent non-linear magneto-hydrodynamical simulations of waves in the atmospheres of roAp stars. We explore a grid of simulations in a wide parameter space. The aim of our study is to understand the influence of the atmosphere and the magnetic field on the propagation and reflection properties of magneto-acoustic waves, formation of shocks and node layers.

We use 2D numerical simulations and eikonal approximation, to study properties of MHD waves traveling below the solar surface through the magnetic structure of sunspots. We consider a series of magnetostatic models of sunspots of different magnetic field strengths, from 10 Mm below the photosphere to the low chromosphere. The purpose of these studies is to quantify the effect of the magnetic field on local helioseismology measurements by modeling waves excited by sub-photospheric sources. Time-distance propagation diagrams and wave travel times are calculated for models of various field strength and compared to the non-magnetic case. The results clearly indicate that the observed time-distance helioseismology signals in sunspot regions correspond to fast MHD waves. The slow MHD waves form a distinctly different pattern in the time-distance diagram, which has not been detected in observations. The numerical results are in good agreement with the solution in the short-wavelength (eikonal) approximation, providing its validation. The frequency dependence of the travel times is in a good qualitative agreement with observations.

Waves connect all the layers of the Sun, from its interior to the upper atmosphere. It is becoming clear now the important role of magnetic field on the wave propagation. Magnetic field modifies propagation speed of waves, thus affecting the conclusions of helioseismological studies. It can change the direction of the wave propagation, help channeling them straight up to the corona, extending the dynamic and magnetic couplings between all the layers. Modern instruments provide measurements of complex patterns of oscillations in solar active regions and of tiny effects such as temporal oscillations of the magnetic field. The physics of oscillations in a variety of magnetic structures of the Sun is similar to that of pulsations of stars that posses strong magnetic fields, such as roAp stars. All these arguments point toward a need of systematic self-consistent modeling of waves in magnetic structures that is able to take into account the complexity of the magnetic field configurations. In this paper, we describe simulations of this kind, summarize our recent findings and bring together results from the theory and observations.

Magnetic field modifies the properties of waves in a complex way. Significant advances has been made recently in our understanding of the physics of waves in solar active regions with the help of analytical theories, numerical simulations, as well as hi-resolution observations. In this contribution we review the current ideas in the field, with the emphasis on theoretical models of waves in sunspots.

In order to understand the influence of magnetic fields on the propagation properties of waves, as derived from different local helioseismology techniques, forward modeling of waves is required. Such calculations need a model in magnetohydrostatic equilibrium as initial atmosphere to propagate oscillations through it. We provide a method to construct such a model in equilibrium for a wide range of parameters to be used for simulations of artificial helioseismologic data. The method combine the advantages of self-similar solutions and current-distributed models. A set of models is developed by numerical integration of magnetohydrostatic equations from the sub-photospheric to chromospheric layers.

We report two-dimensional MHD simulations which demonstrate that photospheric 5-min oscillations can leak into the chromosphere inside small-scale vertical magnetic flux tubes. The results of our numerical experiments are compatible with those inferred from simultaneous spectropolarimetric observations of the photosphere and chromosphere obtained with the Tenerife Infrared Polarimeter (TIP) at 10830 A. We conclude that the efficiency of energy exchange by radiation in the solar photosphere can lead to a significant reduction of the cut-off frequency and may allow for the propagation of the 5 minutes waves vertically into the chromosphere.

We present some results obtained from the synthesis of Stokes profiles in small-scale flux tubes with propagating MHD waves. To that aim, realistic flux tubes showing internal structure have been excited with 5 min period drivers, allowing non-linear waves to propagate inside the magnetic structure. The observational signatures of these waves in Stokes profiles of several spectral lines that are commonly used in spectropolarimetric measurements are discussed.

We present results of non-linear, 2D, numerical simulations of magneto-acoustic wave propagation in the photosphere and chromosphere of small-scale flux tubes with internal structure. Waves with realistic periods of three to five minutes are studied, after applying horizontal and vertical oscillatory perturbations to the equilibrium model. Spurious reflections of shock waves from the upper boundary are minimized thanks to a special boundary condition. This has allowed us to increase the duration of the simulations and to make it long enough to perform a statistical analysis of oscillations. The simulations show that deep horizontal motions of the flux tube generate a slow (magnetic) mode and a surface mode. These modes are efficiently transformed into a slow (acoustic) mode in the vA < cS atmosphere. The slow (acoustic) mode propagates vertically along the field lines, forms shocks and remains always within the flux tube. It might deposit effectively the energy of the driver into the chromosphere. When the driver oscillates with a high frequency, above the cut-off, non-linear wave propagation occurs with the same dominant driver period at all heights. At low frequencies, below the cut-off, the dominant period of oscillations changes with height from that of the driver in the photosphere to its first harmonic (half period) in the chromosphere. Depending on the period and on the type of the driver, different shock patterns are observed.

The reliability of quiet Sun magnetic field diagnostics based on the \ion{Fe}{1} lines at 6302 \AA has been questioned by recent work. We present here the results of a thorough study of high-resolution multi-line observations taken with the new spectro-polarimeter SPINOR, comprising the 5250 and 6302 \AA spectral domains. The observations were analyzed using several inversion algorithms, including Milne-Eddington, LTE with 1 and 2 components, and MISMA codes. We find that the line-ratio technique applied to the 5250 \AA lines is not sufficiently reliable to provide a direct magnetic diagnostic in the presence of thermal fluctuations and variable line broadening. In general, one needs to resort to inversion algorithms, ideally with realistic magneto-hydrodynamical constrains. When this is done, the 5250 \AA lines do not seem to provide any significant advantage over those at 6302 \AA . In fact, our results point towards a better performance with the latter (in the presence of turbulent line broadening). In any case, for very weak flux concentrations, neither spectral region alone provides sufficient constraints to fully disentangle the intrinsic field strengths. Instead, we advocate for a combined analysis of both spectral ranges, which yields a better determination of the quiet Sun magnetic properties. Finally, we propose the use of two other \ion{Fe}{1} lines (at 4122 and 9000 \AA) with identical line opacities that seem to work much better than the others.