We present the first high-resolution, simultaneous observations of the solar chromosphere in the optical and millimeter wavelength ranges, obtained with ALMA and the IBIS instrument at the Dunn Solar Telescope. In this paper we concentrate on the comparison between the brightness temperature observed in ALMA Band 3 (3 mm; 100 GHz) and the core width of the H$\alpha$ 656.3 nm line, previously identified as a possible diagnostic of the chromospheric temperature. We find that in the area of plage, network and fibrils covered by our FOV the two diagnostics are well correlated, with similar spatial structures observed in both. The strength of the correlation is remarkable, given that the source function of the mm-radiation obeys local thermodynamic equilibrium, while the H$\alpha$ line has a source function that deviates significantly from the local Planck function. The observed range of ALMA brightness temperatures is sensibly smaller than the temperature range that was previously invoked to explain the observed width variations in H$\alpha$. We employ analysis from forward modeling with the RH code to argue that the strong correlation between H$\alpha$ width and ALMA brightness temperature is caused by their shared dependence on the population number $n_2$ of the first excited level of hydrogen. This population number drives millimeter opacity through hydrogen ionization via the Balmer continuum, and H$\alpha$ width through a curve-of-growth-like opacity effect. Ultimately, the $n_2$ population is regulated by the enhancement or lack of downward Ly$\alpha$ flux, which coherently shifts the formation height of both diagnostics to regions with different temperature, respectively.

Observations from the \textit{Interface Region Imaging Spectrograph} (\textsl{IRIS}) often reveal significantly broadened and non-reversed profiles of the Mg II h, k and triplet lines at flare ribbons. To understand the formation of these optically thick Mg II lines, we perform plane parallel radiative hydrodynamics modeling with the RADYN code, and then recalculate the Mg II line profiles from RADYN atmosphere snapshots using the radiative transfer code RH. We find that the current RH code significantly underestimates the Mg II h \& k Stark widths. By implementing semi-classical perturbation approximation results of quadratic Stark broadening from the STARK-B database in the RH code, the Stark broadenings are found to be one order of magnitude larger than those calculated from the current RH code. However, the improved Stark widths are still too small, and another factor of 30 has to be multiplied to reproduce the significantly broadened lines and adjacent continuum seen in observations. Non-thermal electrons, magnetic fields, three-dimensional effects or electron density effect may account for this factor. Without modifying the RADYN atmosphere, we have also reproduced non-reversed Mg II h \& k profiles, which appear when the electron beam energy flux is decreasing. These profiles are formed at an electron density of $\sim 8\times10^{14}\ \mathrm{cm}^{-3}$ and a temperature of $\sim1.4\times10^4$ K, where the source function slightly deviates from the Planck function. Our investigation also demonstrates that at flare ribbons the triplet lines are formed in the upper chromosphere, close to the formation heights of the h \& k lines.

The inference of the underlying state of the plasma in the solar chromosphere remains extremely challenging because of the nonlocal character of the observed radiation and plasma conditions in this layer. Inversion methods allow us to derive a model atmosphere that can reproduce the observed spectra by undertaking several physical assumptions. The most advanced approaches involve a depth-stratified model atmosphere described by temperature, line-of-sight velocity, turbulent velocity, the three components of the magnetic field vector, and gas and electron pressure. The parameters of the radiative transfer equation are computed from a solid ground of physical principles. To apply these techniques to spectral lines that sample the chromosphere, NLTE effects must be included in the calculations. We developed a new inversion code STiC to study spectral lines that sample the upper chromosphere. The code is based the RH synthetis code, which we modified to make the inversions faster and more stable. For the first time, STiC facilitates the processing of lines from multiple atoms in non-LTE, also including partial redistribution effects. Furthermore, we include a regularization strategy that allows for model atmospheres with a complex stratification, without introducing artifacts in the reconstructed physical parameters, which are usually manifested in the form of oscillatory behavior. This approach takes steps toward a node-less inversion, in which the value of the physical parameters at each grid point can be considered a free parameter. In this paper we discuss the implementation of the aforementioned techniques, the description of the model atmosphere, and the optimizations that we applied to the code. We carry out some numerical experiments to show the performance of the code and the regularization techniques that we implemented. We made STiC publicly available to the community.

The next generation of solar observatories aim to understand the magnetism of the solar chromosphere. Therefore, it is crucial to understand the polarimetric signatures of chromospheric spectral lines. For this purpose, we here examine the suitability of the three Fraunhofer Mg I b1, b2, and b4 lines at 5183.6, 5172.7, and 5167.3 A, respectively. We start by describing a simplified atomic model of only 6 levels and 3 line transitions for computing the atomic populations of the 3p-4s (multiplet number 2) levels involved in the Mg I b line transitions assuming non-local thermodynamic conditions and considering only the Zeeman effect using the field-free approximation. We test this simplified atom against more complex ones finding that, although there are differences in the computed profiles, they are small compared with the advantages provided by the simple atom in terms of speed and robustness. After comparing the three Mg I lines, we conclude that the most capable one is the b2 line as b1 forms at similar heights and always show weaker polarization signals while b4 is severely blended with photospheric lines. We also compare Mg I b2 with the K I D1 and Ca II 8542 A lines finding that the former is sensitive to the atmospheric parameters at heights that are in between those covered by the latter two lines. This makes Mg I b2 an excellent candidate for future multi-line observations that aim to seamlessly infer the thermal and magnetic properties of different features in the lower solar atmosphere.

We characterize the K I D1 & D2 lines in order to determine whether they could complement the 850 nm window, containing the Ca II infrared triplet lines and several Zeeman sensitive photospheric lines, that was studied previously. We investigate the effect of partial redistribution on the intensity profiles, their sensitivity to changes in different atmospheric parameters, and the spatial distribution of Zeeman polarization signals employing a realistic magnetohydrodynamic simulation. The results show that these lines form in the upper photosphere at around 500 km and that they are sensitive to the line of sight velocity and magnetic field strength at heights where neither the photospheric lines nor the Ca II infrared lines are. However, at the same time, we found that their sensitivity to the temperature essentially comes from the photosphere. Then, we conclude that the K I lines provide a complement to the lines in the 850 nm window for the determination of atmospheric parameters in the upper photosphere, especially for the line of sight velocity and the magnetic field.

Solar chromospheric observations of sunspot umbrae offer an exceptional view of magneto-hydrodynamic wave phenomena. In recent years, a wealth of wave signatures related to propagating magneto-acoustic modes have been presented, which demonstrate complex spatial and temporal structuring of the wave components. Theoretical modelling has demonstrated how these ubiquitous waves are consistent with an m=0 slow magneto-acoustic mode, which are excited by trapped sub-photospheric acoustic (p-mode) waves. However, the spectrum of umbral waves is broad, suggesting that the observed signatures represent the superposition of numerous frequencies and/or modes. We apply Fourier filtering, in both spatial and temporal domains, to extract chromospheric umbral wave characteristics consistent with an m=1 slow magneto-acoustic mode. This identification has not been described before. Angular frequencies of 0.037 +/- 0.007 rad/s (2.1 +/- 0.4 deg/s), corresponding to a period approximately 170 s for the m=1 mode are uncovered for spatial wavenumbers in the range of 0.45<k<0.90 arcsec^-1 (5000-9000 km). Theoretical dispersion relations are solved, with corresponding eigenfunctions computed, which allows the density perturbations to be investigated and compared with our observations. Such magnetohydrodynamic modelling confirms our interpretation that the identified wave signatures are the first direct observations of an m=1 slow magneto-acoustic mode in the chromospheric umbra of a sunspot.

The broadening of the hydrogen lines during flares is thought to result from increased charge (electron, proton) density in the flare chromosphere. However, disagreements between theory and modeling prescriptions have precluded an accurate diagnostic of the degree of ionization and compression resulting from flare heating in the chromosphere. To resolve this issue, we have incorporated the unified theory of electric pressure broadening of the hydrogen lines into the non-LTE radiative transfer code RH. This broadening prescription produces a much more realistic spectrum of the quiescent, A0 star Vega compared to the analytic approximations used as a damping parameter in the Voigt profiles. We test recent radiative-hydrodynamic (RHD) simulations of the atmospheric response to high nonthermal electron beam fluxes with the new broadening prescription and find that the Balmer lines are over-broadened at the densest times in the simulations. Adding many simultaneously heated and cooling model loops as a "multithread" model improves the agreement with the observations. We revisit the three-component phenomenological flare model of the YZ CMi Megaflare using recent and new RHD models. The evolution of the broadening, line flux ratios, and continuum flux ratios are well-reproduced by a multithread model with high-flux nonthermal electron beam heating, an extended decay phase model, and a "hot spot" atmosphere heated by an ultrarelativistic electron beam with reasonable filling factors: 0.1%, 1%, and 0.1% of the visible stellar hemisphere, respectively. The new modeling motivates future work to understand the origin of the extended gradual phase emission.

New spatially scanned spectropolarimetry sunspot observations are made of photospheric atomic and molecular absorption lines near 4135nm. The relative splittings among several atomic lines are measured and shown to agree with values calculated with configuration interaction and intermediate coupling. Large splitting is seen in a line identified with Fe I at 4137nm, showing multiple Stokes V components and an unusual linear polarization. This line will be a sensitive probe of quiet Sun magnetic fields, with a magnetic sensitivity of 2.5 times larger than that of the well-known 1565nm Fe I line.

Transverse magnetohydrodynamic (MHD) waves have been shown to be ubiquitous in the solar atmosphere and can in principle carry sufficient energy to generate and maintain the Sun's million-degree outer atmosphere or corona. However, direct evidence of the dissipation process of these waves and subsequent heating has not yet been directly observed. Here we report on high spatial, temporal, and spectral resolution observations of a solar prominence that show a compelling signature of so-called resonant absorption, a long hypothesized mechanism to efficiently convert and dissipate transverse wave energy into heat. Aside from coherence in the transverse direction, our observations show telltale phase differences around 180 degrees between transverse motions in the plane-of-sky and line-of-sight velocities of the oscillating fine structures or threads, and also suggest significant heating from chromospheric to higher temperatures. Comparison with advanced numerical simulations support a scenario in which transverse oscillations trigger a Kelvin-Helmholtz instability (KHI) at the boundaries of oscillating threads via resonant absorption. This instability leads to numerous thin current sheets in which wave energy is dissipated and plasma is heated. Our results provide direct evidence for wave-related heating in action, one of the candidate coronal heating mechanisms.

Transverse magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere and may be responsible for generating the Sun's million-degree outer atmosphere. However, direct evidence of the dissipation process and heating from these waves remains elusive. Through advanced numerical simulations combined with appropriate forward modeling of a prominence flux tube, we provide the observational signatures of transverse MHD waves in prominence plasmas. We show that these signatures are characterized by thread-like substructure, strong transverse dynamical coherence, an out-of-phase difference between plane-of-the-sky motions and LOS velocities, and enhanced line broadening and heating around most of the flux tube. A complex combination between resonant absorption and Kelvin-Helmholtz instabilities (KHI) takes place in which the KHI extracts the energy from the resonant layer and dissipates it through vortices and current sheets, which rapidly degenerate into turbulence. An inward enlargement of the boundary is produced in which the turbulent flows conserve the characteristic dynamics from the resonance, therefore guaranteeing detectability of the resonance imprints. We show that the features described in the accompanying paper (Okamoto et al. 2015) through coordinated Hinode and IRIS observations match well the numerical results.

(abridged) The heating mechanism at high densities during M dwarf flares is poorly understood. Spectra of M dwarf flares in the optical and near-ultraviolet wavelength regimes have revealed three continuum components during the impulsive phase: 1) an energetically dominant blackbody component with a color temperature of T $\sim$ 10,000 K in the blue-optical, 2) a smaller amount of Balmer continuum emission in the near-ultraviolet at lambda $<$ 3646 Angstroms and 3) an apparent pseudo-continuum of blended high-order Balmer lines. These properties are not reproduced by models that employ a typical "solar-type" flare heating level in nonthermal electrons, and therefore our understanding of these spectra is limited to a phenomenological interpretation. We present a new 1D radiative-hydrodynamic model of an M dwarf flare from precipitating nonthermal electrons with a large energy flux of $10^{13}$ erg cm$^{-2}$ s$^{-1}$. The simulation produces bright continuum emission from a dense, hot chromospheric condensation. For the first time, the observed color temperature and Balmer jump ratio are produced self-consistently in a radiative-hydrodynamic flare model. We find that a T $\sim$ 10,000 K blackbody-like continuum component and a small Balmer jump ratio result from optically thick Balmer and Paschen recombination radiation, and thus the properties of the flux spectrum are caused by blue light escaping over a larger physical depth range compared to red and near-ultraviolet light. To model the near-ultraviolet pseudo-continuum previously attributed to overlapping Balmer lines, we include the extra Balmer continuum opacity from Landau-Zener transitions that result from merged, high order energy levels of hydrogen in a dense, partially ionized atmosphere. This reveals a new diagnostic of ambient charge density in the densest regions of the atmosphere that are heated during dMe and solar flares.

Spectroscopic observations of prominence eruptions associated with coronal mass ejections (CMEs), although relatively rare, can provide valuable plasma and 3D geometry diagnostics. We report the first observations by the Interface Region Imaging Spectrograph (IRIS) mission of a spectacular fast CME/prominence eruption associated with an equivalent X1.6 flare on 2014 May 9. The maximum plane-of-sky and Doppler velocities of the eruption are 1200 and 460 km/s, respectively. There are two eruption components separated by ~200 km/s in Doppler velocity: a primary, bright component and a secondary, faint component, suggesting a hollow, rather than solid, cone-shaped distribution of material. The eruption involves a left-handed helical structure undergoing counter-clockwise (viewed top-down) unwinding motion. There is a temporal evolution from upward eruption to downward fallback with less-than-free-fall speeds and decreasing nonthermal line widths. We find a wide range of Mg II k/h line intensity ratios (less than ~2 expected for optically-thin thermal emission): the lowest ever-reported median value of 1.17 found in the fallback material and a comparably high value of 1.63 in nearby coronal rain and intermediate values of 1.53 and 1.41 in the two eruption components. The fallback material exhibits a strong ($> 5 \sigma$) linear correlation between the k/h ratio and the Doppler velocity as well as the line intensity. We demonstrate that Doppler dimming of scattered chromospheric emission by the erupted material can potentially explain such characteristics.

The emergence of three-dimensional magneto-hydrodynamic (MHD) simulations of stellar atmospheres has sparked a need for efficient radiative transfer codes to calculate detailed synthetic spectra. We present RH 1.5D, a massively parallel code based on the RH code and capable of performing Zeeman polarised multi-level non-local thermodynamical equilibrium (NLTE) calculations with partial frequency redistribution for an arbitrary amount of chemical species. The code calculates spectra from 3D, 2D or 1D atmospheric models on a column-by-column basis (or 1.5D). While the 1.5D approximation breaks down in the cores of very strong lines in an inhomogeneous environment, it is nevertheless suitable for a large range of scenarios and allows for faster convergence with finer control over the iteration of each simulation column. The code scales well to at least tens of thousands of CPU cores, and is publicly available. In the present work we briefly describe its inner workings, strategies for convergence optimisation, its parallelism, and some possible applications.

We revisit an old question: what are the effects of observing stratified atmospheres on scales below a photon mean free path? The mean free path of photons emerging from the solar photosphere and chromosphere is near 100 km. Using current 1m-class telescopes, the mean free path is on the order of the angular resolution. But the Daniel K. Inoue Solar Telescope will have a diffraction limit of 0.020'' near the atmospheric cutoff at 310nm, corresponding to 14 km at the solar surface. Even a small amount of scattering in the source function leads to physical smearing due to this solar 'fog', with effects similar to a degradation of the telescope PSF. We discuss a unified picture that depends simply on the nature and amount of scattering in the source function. Scalings are derived from which the scattering in the solar atmosphere can be transcribed into an effective Strehl ratio, a quantity useful to observers. Observations in both permitted (e.g., Fe I 630.2 nm) and forbidden (Fe I 525.0 nm) lines will shed light on both instrumental performance as well as on small scale structures in the solar atmosphere.

Oxygen is the most abundant element on the Sun after Hydrogen and Helium. The intensity spectrum of resonance lines of neutral Oxygen namely O {\sc i} (1302, 1305 and 1306 \AA\,) has been studied in the literature for chromospheric diagnostics. In this paper we study the resonance scattering polarization in the O {\sc i} line at 1302 \AA\, using two-dimensional radiative transfer in a composite atmosphere constructed using a two-dimensional magneto-hydrodynamical snapshot in the photosphere and columns of the one-dimensional FALC atmosphere in the chromosphere. The methods developed by us recently in a series of papers to solve multi-dimensional polarized radiative transfer have been incorporated in our new code POLY2D which we use for our analysis. We find that multi-dimensional radiative transfer including XRD effects is important in reproducing the amplitude and shape of scattering polarization signals of the O {\sc i} line at 1302 \AA\,.

We present a new determination of the solar fluorine abundance together with abundance measurements of fluorine in two Galactic open clusters. We analyzed a sunspot spectrum, observed by L. Wallace and W. Livingston with the FTS at the McMath/Pierce Solar Telescope situated on Kitt Peak and spectra of four giants in the old cluster M 67 ($\sim$4.5 Gyr) and three giants in the young cluster NGC 6404 ($\sim$0.5 Gyr), obtained with the CRIRES spectrograph at VLT. Fluorine was measured through synthesis of the available HF lines. We adopted the recent set of experimental molecular parameters of HF delivered by the HITRAN database, and found a new solar fluorine abundance of $A(F) = 4.40\pm 0.25$, in good agreement with the M 67 average fluorine abundance of $A(F) = 4.49\pm 0.20$. The new solar abundance is in a very good agreement with the meteoritic value. The used modern spectrosynthesis tools, the agreement with the meteoritic value and with the results in open cluster M67, known to be a solar analogue, make our solar determination very robust. At the same time, the fluorine measurement in the above-mentioned open clusters is the first step in the understanding of its evolution during the last $\sim$10 Gyr in the Galactic disk. In order to develop this project, a larger sample of open clusters is required, so that it would allow us to trace the evolution of fluorine as a function of time and, in turn, to better understand its origin.

Measurements from the Solar Irradiance Monitor (SIM) onboard the SORCE mission indicate that solar spectral irradiance at Visible and IR wavelengths varies in counter phase with the solar activity cycle. The sign of these variations is not reproduced by most of the irradiance reconstruction techniques based on variations of surface magnetism employed so far, and it is not clear yet whether SIM calibration procedures need to be improved, or if instead new physical mechanisms must be invoked to explain such variations. We employ three-dimensional magneto hydrodynamic simulations of the solar photosphere to investigate the dependence of solar radiance in SIM Visible and IR spectral ranges on variations of the filling factor of surface magnetic fields. We find that the contribution of magnetic features to solar radiance is strongly dependent on the location on the disk of the features, being negative close to disk center and positive toward the limb. If features are homogeneously distributed over a region around the equator (activity belt) then their contribution to irradiance is positive with respect to the contribution of HD snapshots, but decreases with the increase of their magnetic flux for average magnetic flux larger than 50 G in at least two of the Visible and IR spectral bands monitored by SIM. Under the assumption that the 50 G snapshots are representative of quiet Sun regions we find thus that the Spectral Irradiance can be in counter-phase with the solar magnetic activity cycle.

The Interface Region Imaging Spectrograph (IRIS) small explorer spacecraft provides simultaneous spectra and images of the photosphere, chromosphere, transition region, and corona with 0.33-0.4 arcsec spatial resolution, 2 s temporal resolution and 1 km/s velocity resolution over a field-of-view of up to 175 arcsec x 175 arcsec. IRIS was launched into a Sun-synchronous orbit on 27 June 2013 using a Pegasus-XL rocket and consists of a 19-cm UV telescope that feeds a slit-based dual-bandpass imaging spectrograph. IRIS obtains spectra in passbands from 1332-1358, 1389-1407 and 2783-2834 Angstrom including bright spectral lines formed in the chromosphere (Mg II h 2803 Angstrom and Mg II k 2796 Angstrom) and transition region (C II 1334/1335 Angstrom and Si IV 1394/1403 Angstrom). Slit-jaw images in four different passbands (C II 1330, Si IV 1400, Mg II k 2796 and Mg II wing 2830 Angstrom) can be taken simultaneously with spectral rasters that sample regions up to 130 arcsec x 175 arcsec at a variety of spatial samplings (from 0.33 arcsec and up). IRIS is sensitive to emission from plasma at temperatures between 5000 K and 10 MK and will advance our understanding of the flow of mass and energy through an interface region, formed by the chromosphere and transition region, between the photosphere and corona. This highly structured and dynamic region not only acts as the conduit of all mass and energy feeding into the corona and solar wind, it also requires an order of magnitude more energy to heat than the corona and solar wind combined. The IRIS investigation includes a strong numerical modeling component based on advanced radiative-MHD codes to facilitate interpretation of observations of this complex region. Approximately eight Gbytes of data (after compression) are acquired by IRIS each day and made available for unrestricted use within a few days of the observation.

G-band bright points are small-sized features characterized by high photometric contrast. Theoretical investigations indicate that these features have associated magnetic field strengths between 1-2 kG. Results from observations instead lead to contradictory results, indicating magnetic fields of only kG strength in some and including hG strengths in others. In order to understand the differences between measurements reported in the literature, and to reconcile them with results from theory, we analyze the distribution of magnetic field strength of G-band bright features identified on synthetic images of the solar photosphere, and its sensitivity to observational and methodological effects. We investigate the dependence of magnetic field strength distributions of G-band bright points identified in 3D magnetohydrodynamic simulations on feature selection method, data sampling, alignment and spatial resolution. The distribution of magnetic field strength of G-band bright features shows two peaks, one at about 1.5 kG and one below 1 hG. The former corresponds to magnetic features,the second mostly to bright granules. Peaks at several hG are obtained only on spatially degraded or misalligned data. Simulations show that magnetic G-band bright points have typically associated field strengths of few kG. Field strengths in the hG range can result from observational effects, thus explaining the discrepancies presented in the literature. Our results also indicate that outcomes from spectro-polarimetric inversions with imposed unit filling-factor should be employed with great caution.

The Mg II h&k lines are the prime chromospheric diagnostics of NASA's Interface Region Imaging Spectrograph (IRIS). In the previous papers of this series we used a realistic three-dimensional radiative magnetohydrodynamics model to calculate the h&k lines in detail and investigated how their spectral features relate to the underlying atmosphere. In this work, we employ the same approach to investigate how the h&k diagnostics fare when taking into account the finite resolution of IRIS and different noise levels. In addition, we investigate the diagnostic potential of several other photospheric lines and near-continuum regions present in the near-ultraviolet (NUV) window of IRIS and study the formation of the NUV slit-jaw images. We find that the instrumental resolution of IRIS has a small effect on the quality of the h&k diagnostics; the relations between the spectral features and atmospheric properties are mostly unchanged. The peak separation is the most affected diagnostic, but mainly due to limitations of the simulation. The effects of noise start to be noticeable at a signal-to-noise ratio (S/N) of 20, but we show that with noise filtering one can obtain reliable diagnostics at least down to a S/N of 5. The many photospheric lines present in the NUV window provide velocity information for at least eight distinct photospheric heights. Using line-free regions in the h&k far wings we derive good estimates of photospheric temperature for at least three heights. Both of these diagnostics, in particular the latter, can be obtained even at S/Ns as low as 5.

The presence of photospheric magnetic reconnection has long been thought to give rise to short and impulsive events, such as Ellerman bombs (EBs) and Type II spicules. In this article, we combine high-resolution, high-cadence observations from the Interferometric BIdimensional Spectrometer (IBIS) and Rapid Oscillations in the Solar Atmosphere (ROSA) instruments at the Dunn Solar Telescope, National Solar Observatory, New Mexico with co-aligned Atmospheric Imaging Assembly (SDO/AIA) and Solar Optical Telescope (Hinode/SOT) data to observe small-scale events situated within an active region. These data are then compared with state-of-the-art numerical simulations of the lower atmosphere made using the MURaM code. It is found that brightenings, in both the observations and the simulations, of the wings of the H alpha line profile, interpreted as EBs, are often spatially correlated with increases in the intensity of the FeI 6302.5A line core. Bi-polar regions inferred from Hinode/SOT magnetic field data show evidence of flux cancellation associated, co-spatially, with these EBs, suggesting magnetic reconnection could be a driver of these high-energy events. Through the analysis of similar events in the simulated lower atmosphere, we are able to infer that line profiles analogous to the observations occur co-spatially with regions of strong opposite polarity magnetic flux. These observed events and their simulated counterparts are interpreted as evidence of photospheric magnetic reconnection at scales observable using current observational instrumentation.

NASA's Interface Region Imaging Spectrograph (IRIS) small explorer mission will study how the solar atmosphere is energized. IRIS contains an imaging spectrograph that covers the Mg II h&k lines as well as a slit-jaw imager centered at Mg II k. Understanding the observations requires forward modeling of Mg II h&k line formation from 3D radiation-MHD models. We compute the vertically emergent h&k intensity from a snapshot of a dynamic 3D radiation-MHD model of the solar atmosphere, and investigate which diagnostic information about the atmosphere is contained in the synthetic line profiles. We find that the Doppler shift of the central line depression correlates strongly with the vertical velocity at optical depth unity, which is typically located less than 200 km below the transition region (TR). By combining the Doppler shifts of the h and the k line we can retrieve the sign of the velocity gradient just below the TR. The intensity in the central line depression is anticorrelated with the formation height, especially in subfields of a few square Mm. This intensity could thus be used to measure the spatial variation of the height of the transition region. The intensity in the line-core emission peaks correlates with the temperature at its formation height, especially for strong emission peaks. The peaks can thus be exploited as a temperature diagnostic. The wavelength difference between the blue and red peaks provides a diagnostic of the velocity gradients in the upper chromosphere. The intensity ratio of the blue and red peaks correlates strongly with the average velocity in the upper chromosphere. We conclude that the Mg II h&k lines are excellent probes of the very upper chromosphere just below the transition region, a height regime that is impossible to probe with other spectral lines.

NASA's Interface Region Imaging Spectrograph (IRIS) space mission will study how the solar atmosphere is energized. IRIS contains an imaging spectrograph that covers the Mg II h&k lines as well as a slit-jaw imager centered at Mg II k. Understanding the observations will require forward modeling of Mg II h&k line formation from 3D radiation-MHD models. This paper is the first in a series where we undertake this forward modeling. We discuss the atomic physics pertinent to h&k line formation, present a quintessential model atom that can be used in radiative transfer computations and discuss the effect of partial redistribution (PRD) and 3D radiative transfer on the emergent line profiles. We conclude that Mg II h&k can be modeled accurately with a 4-level plus continuum Mg II model atom. Ideally radiative transfer computations should be done in 3D including PRD effects. In practice this is currently not possible. A reasonable compromise is to use 1D PRD computations to model the line profile up to and including the central emission peaks, and use 3D transfer assuming complete redistribution to model the central depression.

The contribution of the quiet Sun to solar irradiance variability, due either to changes of the solar effective temperature or to the presence of unresolved magnetic field, is still poorly understood. In this study we investigate spectral line diagnostics that are sensitive to both temperature variations and the presence of small scale unresolved magnetic features in these areas of the solar atmosphere. Specifically we study the dependence on the magnetic flux density of three parameters describing the shape of two magnetically sensitive FeI lines, at 630.2 nm and 617.3 nm, namely the line core intensity (IC), full width at half maximum (FWHM), and the equivalent width (EQW). To this aim we analyze observations of active region NOAA 11172, acquired with IBIS at the Dunn Solar Telescope, as well as results from numerical synthesis. Our results show that IC is sensitive to both temperature and magnetic flux density variations, FWHM is mostly affected by magnetic field changes, and EQW is mostly sensitive to temperature. Variations of a few percent of the measured line parameters are found in observational data that was spatially degraded to represent quiet-Sun, disk-centre, medium resolution observations. It is therefore possible to disentangle magnetic from pure thermodynamic effects by comparison of temporal variations of the EQW and the FWHM of either the two FeI lines.

Radiative transfer modeling of spectral lines including partial redistribution (PRD) effects requires the evaluation of the ratio of the emission to the absorption profile. This quantity requires a large amount of computational work if one employs the angle-dependent redistribution function, which prohibits its use in 3D radiative transfer computations with model atmospheres containing velocity fields. We aim to provide a method to compute the emission to absorption profile ratio that requires less computational work but retains the effect of angle-dependent scattering in the resulting line profiles. We present a method to compute the profile ratio that employs the angle-averaged redistribution function and wavelength transforms to and from the rest frame of the scattering particles. We compare the emergent line profiles of the \MgII\,k and \Lyalpha\ lines computed with angle-dependent PRD, angle-averaged PRD and our new method in two representative test atmospheres. The new method yields a good approximation of true angle-dependent profile ratio and the resulting emergent line profiles while keeping the computational speed and simplicity of angle-averaged PRD theory.

A recent paper states that ultraviolet backradiation from the solar transition region and upper chromosphere strongly affects the degree of ionization of minority stages at the top of the photosphere, i.e., in the temperature minimum of the one-dimensional static model atmospheres presented in that paper. We show that this claim is incompatible with bservations and we demonstrate that the pertinent ionization balances are instead dominated by outward photospheric radiation, as in older static models. We then analyze the formation of Halpha in the above model and show that it has significant backradiation across the opacity gap by which Halpha differs from other strong scatttering lines.

We have investigated the problem of sunspot magnetohydrostatic equilibrium with comprehensive IR sunspot magnetic field survey observations of the highly sensitive Fe I lines at 15650 \AA\ and nearby OH lines. We have found that some sunspots show isothermal increases in umbral magnetic field strength which cannot be explained by the simplified sunspot model with a single-component ideal gas atmosphere assumed in previous investigations. Large sunspots universally display non-linear increases in magnetic pressure over temperature, while small sunspots and pores display linear behavior. The formation of molecules provides a mechanism for isothermal concentration of the umbral magnetic field, and we propose that this may explain the observed rapid increase in umbral magnetic field strength relative to temperature. Existing multi-component sunspot atmospheric models predict that a significant amount of molecular hydrogen (H2) exists in the sunspot umbra. The formation of H2 can significantly alter the thermodynamic properties of the sunspot atmosphere and may play a significant role in sunspot evolution. In addition to the survey observations, we have performed detailed chemical equilibrium calculations with full consideration of radiative transfer effects to establish OH as a proxy for H2, and demonstrate that a significant population of H2 exists in the coolest regions of large sunspots.

We investigate the feasibility of representing a structured multi-dimensional stellar atmosphere with a single one-dimensional average stratification for the purpose of spectral diagnosis of the atmosphere's average spectrum. In particular we construct four different one-dimensional stratifications from a single snapshot of a magneto-hydrodynamic simulation of solar convection: one by averaging its properties over surfaces of constant height, and three different ones by averaging over surfaces of constant optical depth at 500 nm. Using these models we calculate continuum, and atomic and molecular line intensities and their center-to-limb variations. From analysis of the emerging spectra we identify three main reasons why these average representations are inadequate for accurate determination of stellar atmospheric properties through spectroscopic analysis. These reasons are: non-linearity in the Planck function with temperature, which raises the average emergent intensity of an inhomogeneous atmosphere above that of an average-property atmosphere, even if their temperature-optical depth stratification is identical; non-linearities in molecular formation with temperature and density, which raise the abundance of molecules of an inhomogeneous atmosphere over that in a one-dimensional model with the same average properties; the anisotropy of convective motions, which strongly affects the center-to-limb variation of line-core intensities. We argue therefore that a one-dimensional atmospheric model that reproduces the mean spectrum of an inhomogeneous atmosphere necessarily does not reflect the average physical properties of that atmosphere, and are therefore inherently unreliable.

Continuum intensity observations obtained with the Michelson Doppler Imager (MDI) on-board the SoHO mission provide long time series of filtergrams that are ideal for studying the evolution of large-scale phenomena in the solar atmosphere and their dependence on solar activity. These filtergrams, however, are not taken in a pure continuum spectral band, but are constructed from a proxy, namely a combination of filtergrams sampling the Ni I 676.8 nm line. We studied the sensitivity of this continuum proxy to the shape of the nickel line and to the degradation in the instrumental transmission profiles. We compared continuum intensity measurements in the nearby of nickel line with MDI proxy values in three sets of high resolution spectro-polarimetric data obtained with the Interferometric Bidimensional Spectrometer (IBIS), and in synthetic data, obtained from multi-dimensional simulations of magneto-convection and one-dimensional atmosphere models. We found that MDI continuum measurements require brightness corrections which depend on magnetic field strength, temperature and, to a smaller extent, plasma velocity. The correction ranges from 2% to 25% in sunspots, and is, on average, less than 2% for other features. The brightness correction also varies with position on the disk, with larger variations obtained for sunspots, and smaller variations obtained for quiet sun, faculae and micropores. Correction factors derived from observations agree with those deduced from the numerical simulations when observational effects are taken into account. Finally, we found that the investigated potential uncertainties in the transmission characteristics of MDI filters only slightly affect the brightness correction to proxy measurements.

We study the radiative emission of various types of solar features, such as quiet Sun, enhanced network, plage, and bright plage regions, identified on filtergrams taken in the Ca II K line. We analysed fulldisk images obtained with the PSPT, by using three interference filters that sample the Ca II K line with different bandpasses. We studied the dependence of the radiative emission of disk features on the filter bandpass. We also performed a NLTE spectral synthesis of the Ca II K line integrated over the bandpass of PSPT filters. The synthesis was carried out by utilizing both the PRD and CRD with the most recent set of semi empirical atmosphere models in the literature and some earlier atmosphere models. We measured the CLV of intensity values for various solar features identified on PSPT images and compared the results obtained with those derived from the synthesis. We find that CRD calculations derived using the most recent quiet Sun model, on average, reproduce the measured values of the quiet Sun regions slightly more accurately than PRD computations with the same model. This may reflect that the utilized atmospheric model was computed assuming CRD. Calculations with PRD on earlier quiet Sun model atmospheres reproduce measured quantities with a similar accuracy as to that achieved here by applying CRD to the recent model. We also find that the median contrast values measured for most of the identified bright features, disk positions, and filter widths are, on average, a factor 1.9 lower than those derived from PRD simulations performed using the recent bright feature models. The discrepancy between measured and modeled values decreases by 12% after taking into account straylight effects on PSPT images. PRD computations on either the most recent or the earlier atmosphere models of bright features reproduce measurements from plage and bright plage regions with a similar accuracy.

We investigated the radiative emission of different types of solar features in the spectral range of the Ca II K line. We analyzed full-disk 2k x 2k observations from the PSPT Precision Solar Photometric Telescope. The data were obtained by using three narrow-band interference filters that sample the Ca II K line with different pass bands. Two filters are centered in the line core, the other in the red wing of the line. We measured the intensity and contrast of various solar features, specifically quiet Sun (inter-network), network, enhanced network, plage, and bright plage (facula) regions. Moreover, we compared the results obtained with those derived from the numerical synthesis performed for the three PSPT filters with a widely used radiative code on a set of reference semi-empirical atmosphere models.

We discuss spectropolarimetric measurements of photospheric (Fe I 630.25 nm) and chromospheric (Ca II 854.21 nm) spectral lines. Our long-term goal is to diagnose properties of the magnetic field near the base of the corona. We compare ground-based two-dimensional spectropolarimetric measurements with (almost) simultaneous space-based slit spectropolarimetry. The ground-based observations were obtained May 20, 2008, with IBIS in spectropolarimetric mode, The space observations were obtained with the Spectro-Polarimeter aboard the HINODE satellite. The agreement between the near-simultaneous co-spatial IBIS and HINODE Stokes-V profiles at 630.25 nm is excellent, with V/I amplitudes compatible with to within 1 %. IBIS QU measurements are affected by residual crosstalk from V, arising from calibration inaccuracies, not from any inherent limitation of imaging spectroscopy. We use a PCA analysis to quantify the detected cross talk. Chromospheric magnetic fields are difficult to constrain by polarization of Ca II lines alone. However, we demonstrate that high cadence, high angular resolution monochromatic images of fibrils in Ca II and H-alpha, can be used to improve the magnetic field constraints, under conditions of high electrical conductivity. Such work is possible only with time series datasets from two-dimensional spectroscopic instruments under conditions of good seeing.

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

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

There are significant advantages in combining Hinode observations with ground-based instruments that can observe additional spectral diagnostics at higher data rates and with greater flexibility. However, ground-based observations, because of the random effects of weather and seeing as well as the complexities data analysis due to changing instrumental configurations, have traditionally been less efficient than satellite observations in producing useful datasets. Future large ground-based telescopes will need to find new ways to optimize both their operational efficiency and scientific output. We have begun experimenting with service-mode or queue-mode observations at the Dunn Solar Telescope using the Interferometric Bidimensional Spectrometer (IBIS) as part of joint Hinode campaigns. We describe our experiences and the advantages of such an observing mode for solar physics.

The structure and energy balance of the solar chromosphere remain poorly known. We have used the imaging spectrometer IBIS at the Dunn Solar Telescope to obtain fast-cadence, multi-wavelength profile sampling of Halpha and Ca II 854.2 nm over a sizable two-dimensional field of view encompassing quiet-Sun network. We provide a first inventory of how the quiet chromosphere appears in these two lines by comparing basic profile measurements in the form of image displays, temporal-average displays, time slices, and pixel-by-pixel correlations. We find that the two lines can be markedly dissimilar in their rendering of the chromosphere, but that, nevertheless, both show evidence of chromospheric heating, particularly in and around network: Halpha in its core width, Ca II 854.2 in its brightness. We discuss venues for improved modeling.

Filtergrams obtained in Ca II H, Ca II K and H-alpha are often employed as diagnostics of the solar chromosphere. However, the vastly disparate appearance between the typical filtergrams in these different lines calls into question the nature of what is actually being observed. We investigate the lack of obvious structures of magnetic origin such as fibrils and mottles in on-disk Ca II H and K images by directly comparing a temporal sequence of classical Ca II K filtergrams with a co-spatial and co-temporal sequence of spectrally resolved Ca II 854.2 images obtained with the Interferometric Bidimensional Spectrometer (IBIS), considering the effect of both the spectral and spatial smearing. We find that the lack of fine magnetic structuring in Ca II K filtergrams, even with the narrowest available filters, is due to observational effects. Signatures of fibrils remain however in the temporal evolution of the filtergrams, in particular with the evidence of magnetic shadows around the network elements. The Ca II K filtergrams do not appear, however, to properly reflect the high-frequency behavior of the chromosphere. Using the same analysis, we find no significant chromospheric signature in the Hinode/SOT Ca II H quiet-Sun filtergrams. The picture provided by H-alpha and Ca II 854.2, which show significant portions of the chromosphere dominated by magnetic structuring, appears to reflect the true and essential nature of the solar chromosphere. Data which do not resolve, spatially or spectrally, this aspect may misrepresent the behavior the chromosphere.

(Abridged) Aims: In this paper, we seek to establish the suitability of imaging spectroscopy performed in the Ca II 854.2 nm line as a means to investigate the solar chromosphere at high resolution. Methods: We utilize monochromatic images obtained with the Interferometric BIdimensional Spectrometer (IBIS) at multiple wavelengths within the Ca II 854.2 nm line and over several quiet areas. We analyze both the morphological properties derived from narrow-band monochromatic images and the average spectral properties of distinct solar features such as network points, internetwork areas and fibrils. Results: The spectral properties derived over quiet-Sun targets are in full agreement with earlier results obtained with fixed-slit spectrographic observations, highlighting the reliability of the spectral information obtained with IBIS. Furthermore, the very narrowband IBIS imaging reveals with much clarity the dual nature of the Ca II 854.2 nm line: its outer wings gradually sample the solar photosphere, while the core is a purely chromospheric indicator. The latter displays a wealth of fine structures including bright points, akin to the Ca II H2V and K2V grains, as well as fibrils originating from even the smallest magnetic elements. The fibrils occupy a large fraction of the observed field of view even in the quiet regions, and clearly outline atmospheric volumes with different dynamical properties, strongly dependent on the local magnetic topology. This highlights the fact that 1-D models stratified along the vertical direction can provide only a very limited representation of the actual chromospheric physics.

Two of the most attractive spectral windows for spectropolarimetric investigations of the physical properties of the plasma structures in the solar chromosphere and corona are the ones provided by the spectral lines of the He I 10830 A and 5876 A (or D3) multiplets, whose polarization signals are sensitive to the Hanle and Zeeman effects. However, in order to be able to carry out reliable diagnostics, it is crucial to have a good physical understanding of the sensitivity of the observed spectral line radiation to the various competing driving mechanisms. Here we report a series of off-the-limb non-LTE calculations of the He I D3 and 10830 A emission profiles, focusing our investigation on their sensitivity to the EUV coronal irradiation and the model atmosphere used in the calculations. We show in particular that the intensity ratio of the blue to the red components in the emission profiles of the He I 10830 A multiplet turns out to be a good candidate as a diagnostic tool for the coronal irradiance. Measurements of this observable as a function of the distance to the limb and its confrontation with radiative transfer modeling might give us valuable information on the physical properties of the solar atmosphere and on the amount of EUV radiation at relevant wavelengths penetrating the chromosphere from above.

We compare the rms contrast in observed speckle reconstructed G-band images with synthetic filtergrams computed from two magneto-hydrodynamic simulation snapshots. The observations consist of 103 bursts of 80 frames each taken at the Dunn Solar Telescope (DST), sampled at twice the diffraction limit of the telescope. The speckle reconstructions account for the performance of the Adaptive Optics (AO) system at the DST to supply reliable photometry. We find a considerable discrepancy in the observed rms contrast of 14.1% for the best reconstructed images, and the synthetic rms contrast of 21.5% in a simulation snapshot thought to be representative of the quiet Sun. The areas of features in the synthetic filtergrams that have positive or negative contrast beyond the minimum and maximum values in the reconstructed images have spatial scales that should be resolved. This leads us to conclude that there are fundamental differences in the rms G-band contrast between observed and computed filtergrams. On the basis of the substantially reduced granular contrast of 16.3% in the synthetic plage filtergram we speculate that the quiet-Sun may contain more weak magnetic field than previously thought.

We promote the use of narrow-band (0.05 -- 0.20 nm FWHM) imaging in the molecularvCN band head at 388.33 nm as an effective method for monitoring small-scale magnetic field elements because it renders them with exceptionally high contrast. We create synthetic narrow-band CN filtergrams from spectra computed from a three-dimensional snapshot of a magnetohydrodynamic simulation of the solar convection to illustrate the expected high contrast and explain its nature. In addition, we performed observations with the horizontal slit spectrograph at the Dunn Solar Tower at 388.3 nm to experimentally confirm the high bright-point contrast, and to characterize and optimize the transmission profile of a narrow-band (0.04 FWHM) Lyot filter, which was built by Lyot and tailored to the CN band at Sacramento Peak in the early 70's. The presented theoretical computations predict that bright-point contrast in narrow-band (0.04 FWHM) CN filtergrams is more than 3 times higher than in CN filtergrams taken with 1 nm FWHM wide filters, and in typical G-band filtergrams.

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

We present a comparative study of the intensity contrast in synthetic CH-band and violet CN-band filtergrams computed from a high-resolution simulation of solar magnetoconvection. The underlying simulation has an average vertical magnetic field of 250 G with kG fields concentrated in its intergranular lanes, and is representative of a plage region. To simulate filtergrams typically obtained in CH- and CN-band observations we computed spatially resolved spectra in both bands and integrated these spectra over 1 nm FWHM filter functions centred at 430.5 nm and 388.3 nm, respectively. We find that the average contrast of magnetic bright points in the simulated filtergrams is lower in the CN-band by a factor of 0.96. This result strongly contradicts earlier semi-empirical modeling and recent observations, which both etimated that the bright-point contrast in the CN-band is \emph{higher} by a factor of 1.4. We argue that the near equality of the bright-point contrast in the two bands in the present simulation is a natural consequence of the mechanism that causes magnetic flux elements to be particularly bright in the CN and CH filtergrams, namely the partial evacuation of these elements and the concomitant weakening of molecular spectral lines in the filter passbands. We find that the RMS intensity contrast in the whole field-of-view of the filtergrams is 20.5% in the G band and 22.0% in the CN band and conclude that this slight difference in contrast is caused by the shorter wavelength of the latter. Both the bright-point and RMS intensity contrast in the CN band are sensitive to the precise choice of the central wavelength of the filter.

Here we report on the first observational investigation of the Zeeman effect in the G band around 4305 A. Our spectropolarimetric observations of sunspots with the Zurich Imaging Polarimeter at the IRSOL facility confirm our previous theoretical prediction that the molecular Zeeman effect produces measurable circular polarization signatures in several CH lines that are not overlapped with atomic transitions. We also find both circular and linear polarization signals produced by atomic lines whose wavelengths lie in the G band spectral region. Together, such molecular and atomic lines are potentially important for empirical investigations on solar and stellar magnetism. For instance, a comparison between observed and calculated Stokes profiles suggests that the thermodynamical and/or magnetic properties of the photospheric regions of sunspot umbrae are horizontally structured with a component that might be associated with umbral dots.

We investigate the possibility of measuring magnetic field strength in G-band bright points through the analysis of Zeeman polarization in molecular CH lines. To this end we solve the equations of polarized radiative transfer in the G band through a standard plane-parallel model of the solar atmosphere with an imposed magnetic field, and through a more realistic snapshot from a simulation of solar magneto-convection. This region of the spectrum is crowded with many atomic and molecular lines. Nevertheless, we find several instances of isolated groups of CH lines that are predicted to produce a measurable Stokes V signal in the presence of magnetic fields. In part this is possible because the effective Land\'{e} factors of lines in the stronger main branch of the CH A$^{2}\Delta$--X$^{2}\Pi$ transition tend to zero rather quickly for increasing total angular momentum $J$, resulting in a Stokes $V$ spectrum of the G band that is less crowded than the corresponding Stokes $I$ spectrum. We indicate that, by contrast, the effective Land\'{e} factors of the $R$ and $P$ satellite sub-branches of this transition tend to $\pm 1$ for increasing $J$. However, these lines are in general considerably weaker, and do not contribute significantly to the polarization signal. In one wavelength location near 430.4 nm the overlap of several magnetically sensitive and non-sensitive CH lines is predicted to result in a single-lobed Stokes $V$ profile, raising the possibility of high spatial-resolution narrow-band polarimetric imaging. In the magneto-convection snapshot we find circular polarization signals of the order of 1% prompting us to conclude that measuring magnetic field strength in small-scale elements through the Zeeman effect in CH lines is a realistic prospect.