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

We study the properties of internal gravity waves (IGWs) detected in synthetic observations that are obtained from realistic numerical simulation of the solar atmosphere. We used four different simulations of the solar magneto-convection performed using the CO5BOLD code. A magnetic-field-free model and three magnetic models were simulated. The latter three models start with an initial vertical, homogeneous field of 10, 50, and 100 G magnetic flux density, representing different regions of the quiet solar surface. We used the NICOLE code to compute synthetic spectral maps from all the simulated models for the two magnetically insensitive neutral iron lines Fe I 5434 {\AA} and 5576 {\AA}. We find the signatures of the internal gravity waves in the synthetic spectra to be consistent with observations of the real Sun. The phase differences obtained using the spectral lines are significantly different from the phase differences in the simulation. The phase coherency between two atmospheric layers in the gravity wave regime is height dependent and is seen to decrease with the travel distance between the observed layers. We conclude that the energy flux of IGWs determined from the phase difference analysis may be overestimated by an order of magnitude. Spectral lines that are weak and less temperature sensitive may be better suited to detecting internal waves and accurately determining their energy flux in the solar atmosphere.

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

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

We present the results of three-dimensional numerical simulations to investigate the excitation of waves in the magnetic network of the Sun due to footpoint motions of a magnetic flux tube. We consider motions that typically mimic granular buffeting and vortex flows and implement them as driving motions at the base of the flux tube. The driving motions generates various MHD modes within the flux tube and acoustic waves in the ambient medium. The response of the upper atmosphere to the underlying photospheric motion and the role of the flux tube in channeling the waves is investigated. We compute the acoustic energy flux in the various wave modes across different boundary layers defined by the plasma and magnetic field parameters and examine the observational implications for chromospheric and coronal heating.

We present results of three-dimensional numerical simulations of magnetohydrodynamic (MHD) wave propagation in a solar magnetic flux tube. Our study aims at understanding the properties of a range of MHD wave modes generated by different photospheric motions. We consider two scenarios observed in the lower solar photosphere, namely, granular buffeting and vortex-like motion, among the simplest mechanism for the generation of waves within a strong, localized magnetic flux concentration. We show that granular buffeting is likely to generate stronger slow and fast magnetoacoustic waves as compared to swirly motions. Correspondingly, the energy flux transported differs as a result of the driving motions. We also demonstrate that the waves generated by granular buffeting are likely to manifest in stronger emission in the chromospheric network. We argue that different mechanisms of wave generation are active during the evolution of a magnetic element in the intergranular lane, resulting in temporally varying emission at chromospheric heights.

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

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

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