Search

The efficient release of magnetic energy in astrophysical plasmas, such as during solar flares, can in principle be achieved through magnetic diffusion, at a rate determined by the associated electric field. However, attempts at measuring electric fields in the solar atmosphere are scarce, and none exist for sites where the magnetic energy is presumably released. Here, we present observations of an energetic event using the National Science Foundation's Daniel K. Inouye Solar Telescope, where we detect the polarization signature of electric fields associated with magnetic diffusion. We measure the linear and circular polarization across the hydrogen H-epsilon Balmer line at 397 nm at the site of a brightening event in the solar chromosphere. Our spectro-polarimetric modeling demonstrates that the observed polarization signals can only be explained by the presence of electric fields, providing conclusive evidence of magnetic diffusion, and opening a new window for the quantitative study of this mechanism in space plasmas., Comment: 24 pages, 11 figures, published in Nature Communications

The distribution and evolution of the magnetic field at the solar poles through a solar cycle is an important parameter in understanding the solar dynamo. The accurate observations of the polar magnetic flux is very challenging from the ecliptic view, mainly due to (a) geometric foreshortening which limits the spatial resolution, and (b) the oblique view of predominantly vertical magnetic flux elements, which presents rather small line-of-sight component of the magnetic field towards the ecliptic. Due to these effects the polar magnetic flux is poorly measured. Depending upon the measurement technique, longitudinal versus full vector field measurement, where the latter is extremely sensitive to the SNR achieved and azimuth disamiguation problem, the polar magnetic flux measurements could be underestimated or overestimated. To estimate the extent of systematic errors in magetic flux measurements at the solar poles due to aforementioned projection effects we use MHD simulations of quiet sun network as a reference solar atmosphere. Using the numerical model of the solar atmosphere we simulate the observations from the ecliptic as well as from out-of-ecliptic vantage points, such as from a solar polar orbit at various heliographic latitudes. Using these simulated observations we make an assessment of the systematic errors in our measurements of the magnetic flux due to projection effects and the extent of under- or over estimation. We suggest that such effects could contribute to reported missing open magnetic flux in the heliosphere and that the multi-viewpoint observations from out-of-the-ecliptic plane together with innovative Compact Doppler Magnetographs provide the best bet for the future measurements., Comment: Presented as a virtual iPoster in AGU 2021 Fall meeting

Solar observations of carbon monoxide (CO) indicate the existence of lower-temperature gas in the lower solar chromosphere. We present an observation of pores, and quiet-Sun, and network magnetic field regions with CO 4.66 {\mu}m lines by the Cryogenic Infrared Spectrograph (CYRA) at Big Bear Solar Observatory. We used the strong CO lines at around 4.66 {\mu}m to understand the properties of the thermal structures of lower solar atmosphere in different solar features with various magnetic field strengths. AIA 1700 {\AA} images, HMI continuum images and magnetograms are also included in the observation. The data from 3D radiation magnetohydrodynamic (MHD) simulation with the Bifrost code are also employed for the first time to be compared with the observation. We used the RH code to synthesize the CO line profiles in the network regions. The CO 3-2 R14 line center intensity changes to be either enhanced or diminished with increasing magnetic field strength, which should be caused by different heating effects in magnetic flux tubes with different sizes. We find several "cold bubbles" in the CO 3-2 R14 line center intensity images, which can be classified into two types. One type is located in the quiet-Sun regions without magnetic fields. The other type, which has rarely been reported in the past, is near or surrounded by magnetic fields. Notably, some are located at the edge of the magnetic network. The two kinds of cold bubbles and the relationship between cold bubble intensities and network magnetic field strength are both reproduced by the 3D MHD simulation with the Bifrost and RH codes. The simulation also shows that there is a cold plasma blob near the network magnetic fields, causing the observed cold bubbles seen in the CO 3-2 R14 line center image. Our observation and simulation illustrate that the magnetic field plays a vital role in the generation of some CO cold bubbles., Comment: 8 pages, 8 figures. Accepted for publication in A&A

We present observations of a solar plage in the millimeter-continuum with the ALMA and in the Ca 8542 and Na 5896 spectral lines with the Interferometric BIdimensional Spectrometer (IBIS). Our goal is to compare the measurement of local gas temperatures provided by ALMA with the temperature diagnostics provided by non-LTE inversions using STIC. In performing these inversions, we find that using column mass as the reference height scale, rather than optical depth, provides more reliable atmospheric profiles above the temperature minimum and that the treatment of non- LTE hydrogen ionization brings the inferred chromospheric temperatures into better agreement with the ALMA measurements. The Band 3 brightness temperatures are higher but well correlated with the inversion-derived temperatures at the height of formation of the Ca 8542 line core. The Band 6 temperatures instead do not show good correlations with the temperatures at any specific layer in the inverted atmospheres. We then performed inversions that included the millimeter continuum intensities as an additional constraint. Incorporating Band 3 generally resulted in atmospheres showing a strong temperature rise in the upper atmosphere, while including Band 6 led to significant regions of anomalously low temperatures at chromospheric heights. This is consistent with the idea that the Band 6 emission can come from a range of heights. The poor constraints on the chromospheric electron density with existing inversion codes introduces difficulties in determining the height(s) of formation of the millimeter continuum as well as uncertainties in the temperatures derived from the spectral lines., Comment: 23 pages, 17 figures

Red-shifted components of chromospheric emission lines in the hard X-ray impulsive phase of solar flares have recently been studied through their 30 s evolution with the high resolution of IRIS. Radiative-hydrodynamic flare models show that these redshifts are generally reproduced by electron-beam generated chromospheric condensations. The models produce large ambient electron densities, and the pressure broadening of hydrogen Balmer series should be readily detected in observations. To accurately interpret upcoming spectral data of flares with the DKIST, we incorporate non-ideal, non-adiabatic line broadening profiles of hydrogen into the RADYN code. These improvements allow time-dependent predictions for the extreme Balmer line wing enhancements in solar flares. We study two chromospheric condensation models, which cover a range of electron beam fluxes ($1-5 \times 10^{11}$ erg s$^{-1}$ cm$^{-2}$) and ambient electron densities ($1 - 60 \times 10^{13}$ cm$^{-3}$) in the flare chromosphere. Both models produce broadening and redshift variations within 10 s of the onset of beam heating. In the chromospheric condensations, there is enhanced spectral broadening due to large optical depths at H$\alpha$, H$\beta$, and H$\gamma$, while the much lower optical depth of the Balmer series H12$-$H16 provides a translucent window into the smaller electron densities in the beam-heated layers below the condensation. The wavelength ranges of typical DKIST/ViSP spectra of solar flares will be sufficient to test the predictions of extreme hydrogen wing broadening and accurately constrain large densities in chromospheric condensations., Comment: 33 pages, 10 figures, accepted for publication in the Astrophysical Journal

The bright supergiant, Betelgeuse (Alpha Orionis, HD 39801) experienced a visual dimming during 2019 December and the first quarter of 2020 reaching an historic minimum 2020 February 7$-$13. During 2019 September-November, prior to the optical dimming event, the photosphere was expanding. At the same time, spatially resolved ultraviolet spectra using the Hubble Space Telescope/Space Telescope Imaging Spectrograph revealed a substantial increase in the ultraviolet spectrum and Mg II line emission from the chromosphere over the southern hemisphere of the star. Moreover, the temperature and electron density inferred from the spectrum and C II diagnostics also increased in this hemisphere. These changes happened prior to the Great Dimming Event. Variations in the Mg II k-line profiles suggest material moved outwards in response to the passage of a pulse or acoustic shock from 2019 September through 2019 November. It appears that this extraordinary outflow of material from the star, likely initiated by convective photospheric elements, was enhanced by the coincidence with the outward motions in this phase of the $\sim$400 day pulsation cycle. These ultraviolet observations appear to provide the connecting link between the known large convective cells in the photosphere and the mass ejection event that cooled to form the dust cloud in the southern hemisphere imaged in 2019 December, and led to the exceptional optical dimming of Betelgeuse in 2020 February., Comment: 11 pages, 8 figures, Astrophysical Journal, accepted

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., Comment: Accepted for publication in ApJ

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.

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

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., Comment: 31 pages, 13 figures, 2 tables, accepted for publication in the Astrophysical Journal

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., Comment: This is part 2 of a series of 2 papers. Part 1 corresponds to Okamoto et al. (2015, accepted). 36 Pages (single column), 10 figures. Accepted for publication in The Astrophysical Journal

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., Comment: 28 pages, 9 figures, accepted for publication in ApJ. Part II by Patrick Antolin et al. will appear soon

(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., Comment: 50 pages, 2 tables, 13 figures. Accepted for publication in the Solar Physics Topical Issue, "Solar and Stellar Flares". Version 2 (June 22, 2015): updated to include comments by Guest Editor. The final publication is available at Springer via http://dx.doi.org/10.1007/s11207-015-0708-x

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., Comment: 12 pages, 6 figures, accepted by ApJ (Feb 15, 2015)

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., Comment: 6 pages, 3 figures. A&A in press. Updated version reflects changes in latest proofs

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., Comment: Solar Physics, in press

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., Comment: submitted to The Astrophysical Journal

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., Comment: 16 pages, 13 figures. Accepted for publication in ApJ. Updated version with fixed typos in line list and language edits

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., Comment: Paper accepted by ApJ. Content differs from ApJ version for few sentencences

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., Comment: Accepted for publication in A&A

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., Comment: 17 pages, 19 figures, accepted for publication in ApJ

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., Comment: 27 pages, 9 figures

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., Comment: to appear in "The Second Hinode Science Meeting: Beyond Discovery - Toward Understanding", ASP Conference Series, Vol. 415, eds. B. Lites, M. Cheung, T. Magara, J. Mariska, and K. Reeves

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., Comment: 16 pages, 8 figures

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., Comment: 7 pages, 9 figures

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., Comment: 8 pages, 9 figures

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., Comment: 24 pages, 9 figures, submitted to ApJ

We would like the opportunity to thank NASA for supporting our efforts to construct tools to analyze the spectra of spatially inhomogeneous and temporally varying stellar atmospheres. This financial support has allowed us to a versatile radiative transfer code that can be used for many different applications. With this numerical code we have written a point-and-click analysis package written in IDL that can be used to look extensively at the generated output data. Below we describe the most recent results obtained with our transfer code and list papers that have appeared with these results. Although we have not been able to produce as many time-dependent calculations as we had hoped (mainly because of programmatic reasons; Sasselov took another position halfway through the grant), we believe we have

This grant funded an observational and theoretical program to study the structure and dynamics of the solar photosphere and low chromosphere, and the spectral signatures that result. The overall goal is to learn about mechanisms that cause heating of the overlying atmosphere, and produce variability of solar emission in spectral regions important for astrophysics and space physics. The program exploited two new ground-based observational capabilities: one using the Swedish Solar Telescope on La Palma for very high angular resolution observations of the photospheric intensity field (granulation) and proxies of the magnetic field (G-band images); and the other using the Near Infrared Magnetograph at the McMath-Pierce Solar Facility to map the spatial variation and dynamic behavior of the solar temperature minimum region using infrared CO lines. We have interpreted these data using a variety of theoretical and modelling approaches, some developed especially for this project. Previous annual reports cover the work done up to 31 May 1997. This final report summarizes our work for the entire period, including the period of no-cost extension from 1 June 1997 through September 30 1997. In Section 2 we discuss observations and modelling of the photospheric flowfields and their consequences for heating of the overlying atmosphere, and in Section 3 we discuss imaging spectroscopy of the CO lines at 4.67 mu.