Source author record

Mats Carlsson

Mats Carlsson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph astro-ph.IM astro-ph.HE Programming Languages

Catalog footprint

What is connected

63works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

Chromosphere of the quiet sun -- I. Shock and current-sheet dynamics and heating

The solar chromosphere is a crucial interface between the solar interior and its interplanetary environment, regulating how energy is locally deposited into heat and transported into the upper atmospheric layers. Despite significant progress, the dominant processes responsible for chromospheric heating remain debated, particularly under quiet-Sun (QS) conditions. We aim to disentangle and quantify the respective roles of shocks and current sheets (CS) in QS chromospheric modeling. We use a simulation performed with the radiation-magnetohydrodynamics code Bifrost. In order to identify shocks and CS events across space and time, we develop and apply physics-based criteria, allowing us to describe their dynamics and evaluate their contributions to both dissipative (viscous and ohmic) and mechanical (including compressive work) heating. Shocks are found to dominate the energy deposition in the lower chromosphere (up to $59\%$ of the mechanical heating), while CS become the primary contributor in the upper chromosphere, as both plasma $β$ and Mach number $Ma$ drop. Overall, $66\%$ of the mechanical chromospheric heating is powered by the combined action of shocks and CS. These results support a multi-process view of the chromospheric heating in the QS, dominated by shocks, CS, and non-steep gradient dynamics. In addition to viscous and ohmic dissipation, compressive heating can play a major role locally in the model, particularly in chromospheric shock structures, where it offsets non-reversibly cooling from expansion and radiation, and therefore constitutes a key heating contribution to consider in the energy budget. This study further highlights the need for next-generation observations to resolve the intermittent and small-scale nature of chromospheric dynamics, in order to bring new constraints on the coupling between the different layers of the solar atmosphere.

preprint2022arXiv

An approximate recipe of chromospheric radiative losses for solar flares

Radiative losses in the chromosphere are very important in the energy balance. There have been efforts to make simple lookup tables for chromospheric radiative losses in the quiet Sun. During solar flares, the atmospheric conditions are quite different, and the currently available recipe of Gan & Fang (1990) is constructed from semi-empirical models. It remains to be evaluated how these recipes work in flare conditions. We aim to construct an approximate recipe of chromospheric radiative losses for solar flares. We follow the method of Carlsson & Leenaarts (2012) to tabulate the optically thin radiative loss, escape probability, and ionization fraction, while using a grid of flare models from radiative hydrodynamic simulations as our dataset. We provide new lookup tables to calculate chromospheric radiative losses for flares. Compared with previous recipes, our recipe provides a better approximation to the detailed radiative losses for flares.

preprint2022arXiv

Effects of spatial resolution on inferences of atmospheric quantities from simulations

Small scale processes are thought to be important for the dynamics of the solar atmosphere. While numerical resolution fundamentally limits their inclusion in MHD simulations, real observations at the same nominal resolution should still contain imprints of sub-resolution effects. This means that the synthetic observables from a simulation of given resolution might not be directly comparable to real observables at the same resolution. It is thus of interest to investigate how inferences based on synthetic spectra from simulations with different numerical resolutions compare, and whether these differences persist after the spectra have been spatially degraded to a common resolution. We aim to compare synthetic spectra obtained from realistic 3D radiative magnetohydrodynamic (rMHD) simulations run at three different numerical resolutions from the same initial atmosphere, using very simple methods for inferring line-of-sight velocities and magnetic fields. Additionally we examine how the differing spatial resolution impacts the results retrieved from the STiC inversion code. We find that while the simple inferences for all three simulations reveal the same large-scale tendencies, the higher resolutions yield more fine-grained structures and more extreme line-of-sight velocities/magnetic fields in concentrated spots even after spatial smearing. We also see indications that the imprints of sub-resolution effects on the degraded spectra result in systematic errors in the inversions, and that these errors increase with the amount of sub-resolution effects included. Fortunately, however, we find that including successively more sub-resolution yields smaller additional effects; i.e. there is a clear trend of diminishing importance for progressively finer sub-resolution effects.

preprint2022arXiv

The Atmospheric Response to High Nonthermal Electron Beam Fluxes in Solar Flares. II. Hydrogen Broadening Predictions for Solar Flare Observations with the Daniel K. Inouye Solar Telescope

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$α$, H$β$, and H$γ$, 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.

preprint2021arXiv

Mapping Solar Magnetic Fields from the Photosphere to the Base of the Corona

Routine ultraviolet imaging of the Sun's upper atmosphere shows the spectacular manifestation of solar activity; yet we remain blind to its main driver, the magnetic field. Here we report unprecedented spectropolarimetric observations of an active region plage and its surrounding enhanced network, showing circular polarization in ultraviolet (Mg II $h$ & $k$ and Mn I) and visible (Fe I) lines. We infer the longitudinal magnetic field from the photosphere to the very upper chromosphere. At the top of the plage chromosphere the field strengths reach more than 300 gauss, strongly correlated with the Mg II $k$ line core intensity and the electron pressure. This unique mapping shows how the magnetic field couples the different atmospheric layers and reveals the magnetic origin of the heating in the plage chromosphere.

preprint2021arXiv

Probing the physics of the solar atmosphere with the Multi-slit Solar Explorer (MUSE): I. Coronal Heating

The Multi-slit Solar Explorer (MUSE) is a proposed NASA MIDEX mission, currently in Phase A, composed of a multi-slit EUV spectrograph (in three narrow spectral bands centered around 171A, 284A, and 108A) and an EUV context imager (in two narrow passbands around 195A and 304A). MUSE will provide unprecedented spectral and imaging diagnostics of the solar corona at high spatial (<0.5 arcsec), and temporal resolution (down to ~0.5s) thanks to its innovative multi-slit design. By obtaining spectra in 4 bright EUV lines (Fe IX 171A , Fe XV 284A, Fe XIX-Fe XXI 108A) covering a wide range of transition region and coronal temperatures along 37 slits simultaneously, MUSE will for the first time be able to "freeze" (at a cadence as short as 10 seconds) with a spectroscopic raster the evolution of the dynamic coronal plasma over a wide range of scales: from the spatial scales on which energy is released (~0.5 arcsec) to the large-scale often active-region size (170 arcsec x 170 arcsec) atmospheric response. We use advanced numerical modeling to showcase how MUSE will constrain the properties of the solar atmosphere on the spatio-temporal scales (~0.5 arcsec, ~20 seconds) and large field-of-view on which various state-of-the-art models of the physical processes that drive coronal heating, solar flares and coronal mass ejections (CMEs) make distinguishing and testable predictions. We describe how the synergy between MUSE, the single-slit, high-resolution Solar-C EUVST spectrograph, and ground-based observatories (DKIST and others) can address how the solar atmosphere is energized, and the critical role MUSE plays because of the multi-scale nature of the physical processes involved. In this first paper, we focus on how comparisons between MUSE observations and theoretical models will significantly further our understanding of coronal heating mechanisms.

preprint2021arXiv

Probing the Physics of the Solar Atmosphere with the Multi-slit Solar Explorer (MUSE): II. Flares and Eruptions

Current state-of-the-art spectrographs cannot resolve the fundamental spatial (sub-arcseconds) and temporal scales (less than a few tens of seconds) of the coronal dynamics of solar flares and eruptive phenomena. The highest resolution coronal data to date are based on imaging, which is blind to many of the processes that drive coronal energetics and dynamics. As shown by IRIS for the low solar atmosphere, we need high-resolution spectroscopic measurements with simultaneous imaging to understand the dominant processes. In this paper: (1) we introduce the Multi-slit Solar Explorer (MUSE), a spaceborne observatory to fill this observational gap by providing high-cadence (<20 s), sub-arcsecond resolution spectroscopic rasters over an active region size of the solar transition region and corona; (2) using advanced numerical models, we demonstrate the unique diagnostic capabilities of MUSE for exploring solar coronal dynamics, and for constraining and discriminating models of solar flares and eruptions; (3) we discuss the key contributions MUSE would make in addressing the science objectives of the Next Generation Solar Physics Mission (NGSPM), and how MUSE, the high-throughput EUV Solar Telescope (EUVST) and the Daniel K Inouye Solar Telescope (and other ground-based observatories) can operate as a distributed implementation of the NGSPM. This is a companion paper to De Pontieu et al. (2021; arXiv:2106.15584), which focuses on investigating coronal heating with MUSE.

preprint2020arXiv

Characterisation of shock wave signatures at millimetre wavelengths from Bifrost simulations

Observations at millimetre wavelengths provide a valuable tool to study the small scale dynamics in the solar chromosphere. We evaluate the physical conditions of the atmosphere in the presence of a propagating shock wave and link that to the observable signatures in mm-wavelength radiation, providing valuable insights into the underlying physics of mm-wavelength observations. A realistic numerical simulation from the 3D radiative Magnetohydrodynamic (MHD) code Bifrost is used to interpret changes in the atmosphere caused by shock wave propagation. High-cadence (1 s) time series of brightness temperature (T$_\text{b}$) maps are calculated with the Advanced Radiative Transfer (ART) code at the wavelengths $1.309$ mm and $1.204$ mm, which represents opposite sides of spectral band~$6$ of the Atacama Large Millimeter/submillimeter Array (ALMA). An example of shock wave propagation is presented. The brightness temperatures show a strong shock wave signature with large variation in formation height between $\sim0.7$ to $1.4$ Mm. The results demonstrate that millimetre brightness temperatures efficiently track upwardly propagating shock waves in the middle chromosphere. In addition, we show that the gradient of the brightness temperature between wavelengths within ALMA band $6$ can potentially be utilised as a diagnostics tool in understanding the small-scale dynamics at the sampled layers.

preprint2020arXiv

Chromospheric Bubbles in Solar Flares

We analyze a grid of radiative hydrodynamic simulations of solar flares to study the energy balance and response of the atmosphere to nonthermal electron beam heating. The appearance of chromospheric bubbles is one of the most notable features that we find in the simulations. These pockets of chromospheric plasma get trapped between the transition region and the lower atmosphere as it is superheated by the particle beam. The chromospheric bubbles are seen in the synthetic spectra, appearing as an additional component to Balmer line profiles with high Doppler velocities as high as 200 km/s. Their signatures are also visible in the wings of Ca II 8542A line profiles. These bubbles of chromospheric plasma are driven upward by a wave front that is induced by the shock of energy deposition, and require a specific heating rate and atmospheric location to manifest.

preprint2020arXiv

On the formation of Lyman $β$ and the O I 1027 and 1028 Å spectral lines

We study the capabilities of Lyman $β$ and the O I 1027 and 1028 Å spectral lines for understanding the properties of the chromosphere and transition region. The oxygen transitions are located in the wing of Lyman $β$ that is a candidate spectral line for the solar missions Solar Orbiter/SPICE and Solar-C (EUVST). We examine general spectroscopic properties of the three transitions in the quiet Sun by synthesizing them assuming non-local thermal equilibrium taking into account partial redistribution effects. We estimate the heights where the spectral lines are sensitive to the physical parameters computing the response functions to temperature and velocity using a 1D semi-empirical atmospheric model. We also synthesize the intensity spectrum using the 3D enhanced network simulation computed with the BIFROST code. The results indicate that Lyman $β$ is sensitive to the temperature from the middle chromosphere to the transition region while it is mainly sensitive to the line-of-sight velocity at the latter atmospheric layers, around 2000 km above the optical surface. The O I lines form lower in the middle chromosphere, being sensitive to the LOS velocities at lower heights than those covered by Lyman $β$. The spatial distribution of intensity signals computed with the BIFROST atmosphere, as well as the inferred velocities from the line core Doppler shift confirm the previous results. Therefore, these results indicate that the spectral window at 1025 Å contains several spectral lines that complement each other to seamlessly trace the thermal structure and gas dynamics from the middle chromosphere to the lower transition region.

preprint2020arXiv

The Sun at millimeter wavelengths I. Introduction to ALMA Band 3 observations

We present an initial study of one of the first ALMA Band 3 observations of the Sun with the aim to characterise the diagnostic potential of brightness temperatures measured with ALMA on the Sun. The observation covers 48min at a cadence of 2s targeting a Quiet Sun region at disk-centre. Corresponding time series of brightness temperature maps are constructed with the first version of the Solar ALMA Pipeline (SoAP) and compared to simultaneous SDO observations. The angular resolution of the observations is set by the synthesized beam (1.4x2.1as). The ALMA maps exhibit network patches, internetwork regions and also elongated thin features that are connected to large-scale magnetic loops as confirmed by a comparison with SDO maps. The ALMA Band 3 maps correlate best with the SDO/AIA 171, 131 and 304 channels in that they exhibit network features and, although very weak in the ALMA maps, imprints of large-scale loops. A group of compact magnetic loops is very clearly visible in ALMA Band 3. The brightness temperatures in the loop tops reach values of about 8000-9000K and in extreme moments up to 10 000K. ALMA Band 3 interferometric observations from early observing cycles already reveal temperature differences in the solar chromosphere. The weak imprint of magnetic loops and the correlation with the 171, 131, and 304 SDO channels suggests though that the radiation mapped in ALMA Band 3 might have contributions from a larger range of atmospheric heights than previously assumed but the exact formation height of Band 3 needs to be investigated in more detail. The absolute brightness temperature scale as set by Total Power measurements remains less certain and must be improved in the future. Despite these complications and the limited angular resolution, ALMA Band 3 observations have large potential for quantitative studies of the small-scale structure and dynamics of the solar chromosphere.

preprint2019arXiv

Cosmic-Ray Interactions in the Solar Atmosphere

High-energy particles enter the solar atmosphere from Galactic or solar coronal sources, producing an "albedo'' source from the quiet Sun, now observable across a wide range of photon energies. The interaction of high-energy particles in a stellar atmosphere depends essentially upon the joint variation of the magnetic field and the gas, which heretofore has been characterized parametrically as P ~ B^alpha with P the gas pressure and B the magnitude of the magnetic field. We re-examine that parametrization by using a self-consistent 3D MHD model (Bifrost) and show that this relationship tends to P ~ B^{2.9+-0.1} based on the visible portions of the sample of open-field flux tubes in such a model, but with large variations from point to point. This scatter corresponds to the strong meandering of the open-field flux tubes in the lower atmosphere, which will have a strong effect on the prediction of the emission anisotropy (limb brightening). The simulations show that much of the open flux in coronal holes originates in weak-field regions within the granular pattern of the convective motions seen in the simulations.

preprint2019arXiv

Ion-neutral interactions and non-equilibrium ionization in the solar chromosphere

The thermal structure of the chromosphere is regulated through a complex interaction of various heating processes, radiative cooling, and the ionization degree of the plasma. Here we study the impact on the thermal properties of the chromosphere when including the combined action of non-equilibrium ionization (NEI) of hydrogen and helium and ion-neutral interaction effects. We have performed a 2.5D radiative magnetohydrodynamic simulation including ion-neutral interaction effects by solving the generalized Ohm's law (GOL) as well as NEI for hydrogen and helium using the Bifrost code. The GOL equation includes ambipolar diffusion and the Hall term. We compare this simulation with another simulation that computes the ionization in local thermodynamic equilibrium (LTE) including ion-neutral interaction effects. Our numerical models reveal substantial thermal differences in magneto-acoustic shocks, the wake behind the shocks, spicules, low-lying magnetic loops, and the transition region. In particular, we find that heating through ambipolar diffusion in shock wakes is substantially less efficient, while in the shock fronts themselves it is more efficient, under NEI conditions than when assuming LTE.

preprint2016arXiv

Chromospheric and Coronal Wave Generation in a Magnetic Flux Sheath

Using radiation magnetohydrodynamic simulations of the solar atmospheric layers from the upper convection zone to the lower corona, we investigate the self-consistent excitation of slow magneto-acoustic body waves (slow modes) in a magnetic flux concentration. We find that the convective downdrafts in the close surroundings of a two-dimensional flux slab "pump" the plasma inside it in the downward direction. This action produces a downflow inside the flux slab, which encompasses ever higher layers, causing an upwardly propagating rarefaction wave. The slow mode, excited by the adiabatic compression of the downflow near the optical surface, travels along the magnetic field in the upward direction at the tube speed. It develops into a shock wave at chromospheric heights, where it dissipates, lifts the transition region, and produces an offspring in the form of a compressive wave that propagates further into the corona. In the wake of downflows and propagating shock waves, the atmosphere inside the flux slab in the chromosphere and higher tends to oscillate with a period of $ν\approx 4$~mHz. We conclude that this process of "magnetic pumping" is a most plausible mechanism for the direct generation of longitudinal chromospheric and coronal compressive waves within magnetic flux concentrations, and it may provide an important heat source in the chromosphere. It may also be responsible for certain types of dynamic fibrils.

preprint2016arXiv

Emergence of granular-sized magnetic bubbles through the solar atmosphere. III. The path to the transition region

We study the ascent of granular-sized magnetic bubbles from the solar photosphere through the chromosphere into the transition region and above, for the first time. Such events occurred in a flux emerging region in NOAA 11850 on September 25, 2013. During that time, the first co-observing campaign between the Swedish 1-m Solar Telescope and the IRIS spacecraft was carried out. Simultaneous observations of the chromospheric H$α$ 656.28 nm and \ion{Ca}{2} 854.2 nm lines, plus the photospheric \ion{Fe}{1} 630.25 nm line, were made with the CRISP spectropolarimeter at the SST reaching a spatial resolution of 0."14. At the same time, IRIS was performing a four-step dense raster of the said emerging flux region, taking slit-jaw images at 133 (C~{\sc ii}, transition region), 140 (\ion{Si}{4}, transition region), 279.6 (\ion{Mg}{2} k, core, upper chromosphere), and 283.2 nm (\ion{Mg}{2} k, wing, photosphere). Spectroscopy of several lines was performed by the IRIS spectrograph in the far and near ultraviolet, of which we have used the \ion{Si}{4} 140.3 and the \ion{Mg}{2} k 279.6 nm lines. Coronal images from the Atmospheric Imaging Assembly of the Solar Dynamics Observatory were used to investigate the possible coronal signatures of the flux emergence events. The photospheric and chromospheric properties of small-scale emerging magnetic bubbles have been described in detail in Ortiz et al. (2014). Here we are able to follow such structures up to the transition region. We describe the properties, including temporal delays, of the observed flux emergence in all layers. We believe this may be an important mechanism of transporting energy and magnetic flux from subsurface layers to the transition region and corona.

preprint2016arXiv

Modeling Repeatedly Flaring $δ$ Sunspots

Active regions (AR) appearing on the surface of the Sun are classified into $α$, $β$, $γ$, and $δ$ by the rules of the Mount Wilson Observatory, California on the basis of their topological complexity. Amongst these, the $δ$-sunspots are known to be super-active and produce the most X-ray flares. Here, we present results from a simulation of the Sun by mimicking the upper layers and the corona, but starting at a more primitive stage than any earlier treatment. We find that this initial state consisting of only a thin sub-photospheric magnetic sheet breaks into multiple flux-tubes which evolve into a colliding-merging system of spots of opposite polarity upon surface emergence, similar to those often seen on the Sun. The simulation goes on to produce many exotic $δ$-sunspot associated phenomena: repeated flaring in the range of typical solar flare energy release and ejective helical flux ropes with embedded cool-dense plasma filaments resembling solar coronal mass ejections.

preprint2016arXiv

The Appearance of Spicules in High Resolution Observations of Ca II H and H-alpha

Solar spicules are chromospheric fibrils that appear everywhere on the Sun, yet their origin is not understood. Using high resolution observations of spicules obtained with the Swedish 1-m Solar Telescope, we aim to understand how spicules appear in filtergrams and Dopplergrams, how they compare in Ca II H and H-alpha, and what can make them appear and disappear. We find thatspicules display a rich and detailed spatial structure, and show a distribution of transverse velocities that when aligned with the line of sight can make them appear at different H-alpha wing positions. They become more abundant at positions closer to the line core, reflecting a distribution of Doppler shifts and widths. In H-alpha width maps they stand out as bright features both on disk and off limb, reflecting their large Doppler motions and possibly higher temperatures than in the typical H-alpha formation region. Spicule lifetimes measured from narrowband images at only a few positions will be an underestimate because Doppler shifts can make them disappear prematurely from such images; for such cases, width maps are a more robust tool. In H-alpha and Ca II H filtergrams, off limb spicules essentially have the same properties, appearance, and evolution. We find that the sudden appearance of spicules can be explained by Doppler shifts from their transverse motions, and does not require other convoluted explanations.

preprint2016arXiv

The cause of spatial structure in solar He I 1083 nm multiplet images

Context. The He i 1083 nm is a powerful diagnostic for inferring properties of the upper solar chromosphere, in particular for the magnetic field. The basic formation of the line in one-dimensional models is well understood, but the influence of the complex 3D structure of the chromosphere and corona has however never been investigated. This structure must play an essential role because images taken in He i 1083 nm show structures with widths down to 100 km. Aims. To understand the effect of the three-dimensional temperature and density structure in the solar atmosphere on the formation of the He i 1083 nm line. Methods. We solve the non-LTE radiative transfer problem assuming statistical equilibrium for a simple 9-level helium atom that nevertheless captures all essential physics. As a model atmosphere we use a snapshot from a 3D radiation-MHD simulation computed with the Bifrost code. Ionising radiation from the corona is self-consistently taken into account. Results. The emergent intensity in the He i 1083 nm is set by the source function and the opacity in the upper chromosphere. The former is dominated by scattering of photospheric radiation and does not vary much with spatial location. The latter is determined by the photonionisation rate in the He i ground state continuum, as well as the electron density in the chromosphere. The spatial variation of the flux of ionising radiation is caused by the spatially-structured emissivity of the ionising photons from material at T = 100 kK in the transition region. The hotter coronal material produces more ionising photons, but the resulting radiation field is smooth and does not lead to small-scale variation of the UV flux. The corrugation of the transition region further increases the spatial variation of the amount of UV radiation in the chromosphere.

preprint2015arXiv

A publicly available simulation of an enhanced network region of the Sun

Context. The solar chromosphere is the interface between the solar surface and the solar corona. Modelling of this region is difficult because it represents the transition from optically thick to thin radiation escape, from gas-pressure domination to magnetic-pressure domination, from a neutral to an ionised state, from MHD to plasma physics, and from near-equilibrium (LTE) to non-equilibrium conditions. Aims. Our aim is to provide the community with realistic simulations of the magnetic solar outer atmosphere. This will enable detailed comparison of existing and upcoming observations with synthetic observables from the simulations, thereby elucidating the complex interactions of magnetic fields and plasma that are crucial for our understanding of the dynamic outer atmosphere. Methods. We used the radiation magnetohydrodynamics code Bifrost to perform simulations of a computational volume with a magnetic field topology similar to an enhanced network area on the Sun. Results. The full simulation cubes are made available online. The general properties of the simulation are discussed, and limitations are discussed.

preprint2015arXiv

A Unified Computational Model for Solar and Stellar Flares

We present a unified computational framework which can be used to describe impulsive flares on the Sun and on dMe stars. The models assume that the flare impulsive phase is caused by a beam of charged particles that is accelerated in the corona and propagates downward depositing energy and momentum along the way. This rapidly heats the lower stellar atmosphere causing it to explosively expand and dramatically brighten. Our models consist of flux tubes that extend from the sub-photosphere into the corona. We simulate how flare-accelerated charged particles propagate down one-dimensional flux tubes and heat the stellar atmosphere using the Fokker-Planck kinetic theory. Detailed radiative transfer is included so that model predictions can be directly compared with observations. The flux of flare-accelerated particles drives return currents which additionally heat the stellar atmosphere. These effects are also included in our models. We examine the impact of the flare-accelerated particle beams on model solar and dMe stellar atmospheres and perform parameter studies varying the injected particle energy spectra. We find the atmospheric response is strongly dependent on the accelerated particle cutoff energy and spectral index.

preprint2015arXiv

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares. II. Inclusion of Radiative Transfer with RADYN

Solar flares involve complex processes that are coupled and span a wide range of temporal, spatial, and energy scales. Modeling such processes self-consistently has been a challenge in the past. Here we present results from simulations that couple particle kinetics with hydrodynamics of the atmospheric plasma. We combine the Stanford unified Fokker-Planck code that models particle acceleration and transport with the RADYN hydrodynamic code that models the atmospheric response to collisional heating by accelerated electrons through detailed radiative transfer calculations. We perform simulations using two different electron spectra, one an {\it ad hoc} power law and the other predicted by the model of stochastic acceleration by turbulence or plasma waves. Surprisingly, the later model, even with energy flux $\ll 10^{10}$ erg s$^{-1}$ cm$^{-2}$, can cause "explosive" chromospheric evaporation and drive stronger up- and downflows (and hydrodynamic shocks). This is partly because our acceleration model, like many others, produces a spectrum consisting of a quasi-thermal component plus a power-law tail. We synthesize emission line profiles covering different heights in the lower atmosphere, including H$α$ 6563 Å, HeII 304 Å, CaII K 3934 Å and SiIV 1393 Å. One interesting result is the unusual high temperature (up to a few $10^5$ K) of the formation site of HeII 304 Å, which is expected due to photoionization-recombination under flare conditions, compared to those in the quiet Sun dominated by collisional excitation. When compared with observations, our results can constrain the properties of non-thermal electrons and thus the poorly understood particle acceleration mechanism.

preprint2015arXiv

Ellerman bombs at high resolution III. Simultaneous observations with IRIS and SST

Ellerman bombs are transient brightenings of the extended wings of the solar Balmer lines in emerging active regions. We describe their properties in the ultraviolet lines sampled by the Interface Region Imaging Spectrograph (IRIS), using simultaneous imaging spectroscopy in H$α$ with the Swedish 1-m Solar Telescope (SST) and ultraviolet images from the Solar Dynamics Observatory for Ellerman bomb detection and identification. We select multiple co-observed Ellerman bombs for detailed analysis. The IRIS spectra strengthen the view that Ellerman bombs mark reconnection between bipolar kilogauss fluxtubes with the reconnection and the resulting bi-directional jet located within the solar photosphere and shielded by overlying chromospheric fibrils in the cores of strong lines. The spectra suggest that the reconnecting photospheric gas underneath is heated sufficiently to momentarily reach stages of ionization normally assigned to the transition region and the corona. We also analyze similar outburst phenomena that we classify as small flaring arch filaments and ascribe to higher-located reconnection. They have different morphology and produce hot arches in million-Kelvin diagnostics.

preprint2015arXiv

Evidence for a Transition Region response to penumbral microjets in sunspots

Penumbral microjets are short-lived, fine-structured and bright jets that are generally observed in chromospheric imaging of the penumbra of sunspots. Here we investigate their potential transition region signature, by combining observations with the Swedish 1-m Solar Telescope (SST) in the Ca II H and Ca II 8542Å lines with ultraviolet imaging and spectroscopy obtained with the Interface Region Imaging Spectrograph (IRIS), which includes the C II 1334/1335Å, Si IV 1394/1403Å and Mg II h & k 2803/2796Å lines. We find a clear corresponding signal in the IRIS Mg II k, C II and Si IV slit-jaw images, typically offset spatially from the Ca II signature in the direction along the jets: from base to top, the penumbral microjets are predominantly visible in Ca II, Mg II k and C II/Si IV, suggesting progressive heating to transition region temperatures along the jet extent. Hence, these results support the suggestion from earlier studies that penumbral microjets may heat to transition region temperatures.

preprint2015arXiv

First High-resolution Spectroscopic Observations of an Erupting Prominence Within a Coronal Mass Ejection by the Interface Region Imaging Spectrograph (IRIS)

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 σ$) 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.

preprint2015arXiv

Internetwork chromospheric bright grains observed with IRIS

The Interface Region Imaging Spectrograph (IRIS) reveals small-scale rapid brightenings in the form of bright grains all over coronal holes and the quiet sun. These bright grains are seen with the IRIS 1330 Å, 1400 Å and 2796 Å slit-jaw filters. We combine coordinated observations with IRIS and from the ground with the Swedish 1-m Solar Telescope (SST) which allows us to have chromospheric (Ca II 8542 Å, Ca II H 3968 Å, Hα, and Mg II k 2796 Å), and transition region (C II 1334 Å, Si IV 1402) spectral imaging, and single-wavelength Stokes maps in Fe I 6302 Åat high spatial (0.33"), temporal and spectral resolution. We conclude that the IRIS slit-jaw grains are the counterpart of so-called acoustic grains, i.e., resulting from chromospheric acoustic waves in a non-magnetic environment. We compare slit-jaw images with spectra from the IRIS spectrograph. We conclude that the grain intensity in the 2796 Å slit-jaw filter comes from both the Mg II k core and wings. The signal in the C II and Si IV lines is too weak to explain the presence of grains in the 1300 and 1400 Å slit-jaw images and we conclude that the grain signal in these passbands comes mostly from the continuum. Even though weak, the characteristic shock signatures of acoustic grains can often be detected in IRIS C II spectra. For some grains, spectral signature can be found in IRIS Si IV. This suggests that upward propagating acoustic waves sometimes reach all the way up to the transition region.

preprint2015arXiv

New Insights into White-Light Flare Emission from Radiative-Hydrodynamic Modeling of a Chromospheric Condensation

(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.

preprint2015arXiv

Non-equilibrium helium ionization in an MHD simulation of the solar atmosphere

The ionization state of the gas in the dynamic solar chromosphere can depart strongly from the instantaneous statistical equilibrium commonly assumed in numerical modeling. We improve on earlier simulations of the solar atmosphere that only included non-equilbrium hydrogen ionization by performing a 2D radiation-magneto-hydrodynamics simulation featuring non-equilibrium ionization of both hydrogen and helium. The simulation includes the effect of hydrogen Lyman-$α$ and the EUV radiation from the corona on the ionization and heating of the atmosphere. Details on code implementation are given. We obtain helium ion fractions that are far from their equilibrium values. Comparison with models with LTE ionization shows that non-equilibrium helium ionization leads to higher temperatures in wave fronts and lower temperatures in the gas between shocks. Assuming LTE ionization results in a thermostat-like behaviour with matter accumulating around the temperatures where the LTE ionization fractions change rapidly. Comparison of DEM curves computed from our models shows that non-equilibrium ionization leads to more radiating material in the temperature range 11-18 kK compared to models with LTE helium ionization. We conclude that non-equilbrium helium ionization is important for the dynamics and thermal structure of the upper chromosphere and transition region. It might also help resolve the problem that intensities of chromospheric lines computed from current models are smaller than those observed.

preprint2015arXiv

Numerical Simulations of Coronal Heating through Footpoint Braiding

Advanced 3D radiative MHD simulations now reproduce many properties of the outer solar atmosphere. When including a domain from the convection zone into the corona, a hot chromosphere and corona are self-consistently maintained. Here we study two realistic models, with different simulated area, magnetic field strength and topology, and numerical resolution. These are compared in order to characterize the heating in the 3D-MHD simulations which self-consistently maintains the structure of the atmosphere. We analyze the heating at both large and small scales and find that heating is episodic and highly structured in space, but occurs along loop shaped structures, and moves along with the magnetic field. On large scales we find that the heating per particle is maximal near the transition region and that widely distributed opposite-polarity field in the photosphere leads to a greater heating scale height in the corona. On smaller scales, heating is concentrated in current sheets, the thicknesses of which are set by the numerical resolution. Some current sheets fragment in time, this process occurring more readily in the higher-resolution model leading to spatially highly intermittent heating. The large scale heating structures are found to fade in less than about five minutes, while the smaller, local, heating shows time scales of the order of 2 minutes in one model and 1 minutes in the other, higher-resolution, model.

preprint2015arXiv

Observed Variability of the Solar Mg II h Spectral Line

The Mg II h&k doublet are two of the primary spectral lines observed by the Sun-pointing Interface Region Imaging Spectrograph (IRIS). These lines are tracers of the magnetic and thermal environment that spans from the photosphere to the upper chromosphere. We use a double gaussian model to fit the Mg II h profile for a full-Sun mosaic dataset taken 24-Aug-2014. We use the ensemble of high-quality profile fits to conduct a statistical study on the variability of the line profile as it relates the magnetic structure, dynamics, and center-to-limb viewing angle. The average internetwork profile contains a deeply reversed core and is weakly asymmetric at h2. In the internetwork, we find a strong correlation between h3 wavelength and profile asymmetry as well h1 width and h2 width. The average reversal depth of the h3 core is inversely related to the magnetic field. Plage and sunspots exhibit many profiles which do not contain a reversal. These profiles also occur infrequently in the internetwork. We see indications of magnetically aligned structures in plage and network in statistics associated with the line core, but these structures are not clear or extended in the internetwork. The center-to-limb variations are compared with predictions of semi-empirical model atmospheres. We measure a pronounced limb darkening in the line core which is not predicted by the model. The aim of this work is to provide a comprehensive measurement baseline and preliminary analysis on the observed structure and formation of the Mg II profiles observed by IRIS.

preprint2015arXiv

On fibrils and field lines: The nature of H$α$ fibrils in the solar chromosphere

Observations of the solar chromosphere in the line-core of the \Halpha\ line show dark elongated structures called fibrils that show swaying motion. We performed a 3D radiation-MHD simulation of a network region, and computed synthetic \Halpha\ images from this simulation to investigate the relation between fibrils and the magnetic field lines in the chromosphere. The periods, amplitudes and phase-speeds of the synthetic fibrils are consistent with those observed. We analyse the relation between the synthetic fibrils and the field lines threading through them, and find that some fibrils trace out the same field line along the fibril's length, but there are also fibrils that sample different field lines at different locations along their length. Fibrils sample the same field lines on a time scale of $\sim200$~s. This is shorter than their own lifetime. We analysed the evolution of the atmosphere along a number of field lines that thread through fibrils and find that they carry slow-mode waves that load mass into the field line, as well as transverse waves that propagate with the Alfvén speed. Transverse waves propagating in opposite directions cause an interference pattern with complex apparent phase speeds. The relationship between fibrils and field lines is complex. It is governed by constant migration and swaying of the field lines, their mass loading by slow modes and subsequent draining, and their actual visibility in \Halpha. Field lines are visible where they lie close to the optical depth unity surface. The location of the latter is governed by the height at which the column mass in the chromosphere reaches a certain value. We conclude that using the swaying motion of fibrils as a tracer of chromospheric transverse oscillations must be done with caution.

preprint2015arXiv

Probing the Sun with ALMA: observations and simulations

ALMA will open a new chapter in the study of the Sun by providing a leap in spatial resolution and sensitivity compared to currently available mm wave- length observations. In preparation of ALMA, we have carried out a large number of observational tests and state-of-the-art radiation MHD simulations. Here we review the best available observations of the Sun at millimeter wavelengths.Using state of the art radiation MHD simulations of the solar atmosphere we demonstrate the huge potential of ALMA observations for uncovering the nature of the solar chromosphere. We show that ALMA will not only provide a reliable probe of the thermal structure and dynamics of the chromosphere, it will also open up a powerful new diagnostic of magnetic field at chromospheric heights, a fundamentally important, but so far poorly known parameter.

preprint2015arXiv

The formation of IRIS diagnostics V. A quintessential model atom of C II and general formation properties of the C II lines at 133.5 nm

The 133.5 nm lines are important observables for the NASA/SMEX mission Interface Region Imaging Spectrograph (IRIS). To make 3D non-LTE radiative transfer computationally feasible it is crucial to have a model atom with as few levels as possible while retaining the main physical processes. We here develop such a model atom and we study the general formation properties of the C II lines. We find that a nine-level model atom of C I-C III with the transitions treated assuming complete frequency redistribution (CRD) suffices to describe the 133.5 nm lines. 3D scattering effects are important for the intensity in the core of the line. The lines are formed in the optically thick regime. The core intensity is formed in layers where the temperature is about 10kK at the base of the transition region. The lines are 1.2-4 times wider than the atomic absorption profile due to the formation in the optically thick regime. The smaller opacity broadening happens for single peak intensity profiles where the chromospheric temperature is low with a steep source function increase into the transition region, the larger broadening happens when there is a temperature increase from the photosphere to the low chromosphere leading to a local source function maximum and a double peak intensity profile with a central reversal. Assuming optically thin formation with the standard coronal approximation leads to several errors: Neglecting photoionization severly underestimates the amount of C II at temperatures below 16kK, erroneously shifts the formation from 10kK to 25kK and leads to too low intensities.

preprint2015arXiv

The formation of IRIS diagnostics VI. The Diagnostic Potential of the C II Lines at 133.5 nm in the Solar Atmosphere

We use 3D radiation magnetohydrodynamic models to investigate how the thermodynamic quantities in the simulation are encoded in observable quantities, thus exploring the diagnostic potential of the 133.5 nm lines. We find that the line core intensity is correlated with the temperature at the formation height but the correlation is rather weak, especially when the lines are strong. The line core Doppler shift is a good measure of the line-of-sight velocity at the formation height. The line width is both dependent on the width of the absorption profile (thermal and non-thermal width) and an opacity broadening factor of 1.2-4 due to the optically thick line formation with a larger broadening for double peak profiles. The 133.5 nm lines can be formed both higher and lower than the core of the Mg II k line depending on the amount of plasma in the 14-50 kK temperature range. More plasma in this temperature range gives a higher 133.5 nm formation height relative to the Mg II k line core. The synthetic line profiles have been compared with IRIS observations. The derived parameters from the simulated line profiles cover the parameter range seen in observations but on average the synthetic profiles are too narrow. We interpret this discrepancy as a combination of a lack of plasma at chromospheric temperatures in the simulation box and too small non-thermal velocities. The large differences in the distribution of properties between the synthetic profiles and the observed ones show that the 133.5 nm lines are powerful diagnostics of the upper chromosphere and lower transition region.

preprint2015arXiv

The formation of IRIS diagnostics VIII. IRIS observations in the C II 133.5 nm multiplet

The C II 133.5 nm multiplet has been observed by NASA's Interface Region Imaging Spectrograph (IRIS) in unprecedented spatial resolution. The aims of this work are to characterize these new observations of the C II lines, place them in context with previous work, and to identify any additional value the C II lines bring when compared with other spectral lines. We make use of wide, long exposure IRIS rasters covering the quiet Sun and an active region. Line properties such as velocity shift and width are extracted from individual spectra and analyzed. The lines have a variety of shapes (mostly single-peak or double-peak), are strongest in active regions and weaker in the quiet Sun. The ratio between the 133.4 nm and 133.5 nm components is always less than 1.8, indicating that their radiation is optically thick in all locations. Maps of the C II line widths are a powerful new diagnostic of chromospheric structures, and their line shifts are a robust velocity diagnostic. Compared with earlier quiet Sun observations, we find similar absolute intensities and mean line widths, but smaller red shifts; this difference can perhaps be attributed to differences in spectral resolution and spatial coverage. The C II intensity maps are somewhat similar to those of transition region lines, but also share some features with chromospheric maps such as those from the Mg II k line, indicating that they are formed between the upper chromosphere and transition region. C II intensity, width, and velocity maps can therefore be used to gather additional information about the upper chromosphere.

preprint2015arXiv

The Formation of IRIS diagnostics. IV. The Mg II triplet lines as a new diagnostic for lower chromospheric heating

A triplet of subordinate lines of Mg II exists in the region around the h&k lines. In solar spectra these lines are seen mostly in absorption, but in some cases can become emission lines. The aim of this work is to study the formation of this triplet, and investigate any diagnostic value they can bring. Using 3D radiative magnetohydrodynamic simulations of quiet Sun and flaring flux emergence, we synthesize spectra and investigate how spectral features respond to the underlying atmosphere. We find that emission in the lines is rare and is typically caused by a steep temperature increase in the lower chromosphere (above 1500 K, with electron densities above 10$^{18}$ m$^{-3}$). In both simulations the lines are sensitive to temperature increases taking place at column masses >= 5e-4 g cm$^{-2}$. Additional information can also be inferred from the peak-to-wing ratio and shape of the line profiles. Using observations from NASA's Interface Region Imaging Spectrograph we find both absorption and emission line profiles with similar shapes to the synthetic spectra, which suggests that these lines represent a useful diagnostic that complements the MgII h&k lines.

preprint2015arXiv

The Formation of IRIS Diagnostics. VII. The Formation of the OI 135.56 nm Line in the Solar Atmosphere

The O I 135.56 nm line is covered by NASA's Interface Region Imaging Spectrograph (IRIS) small explorer mission which studies how the solar atmosphere is energized. We here study the formation and diagnostic potential of this line by means of non-LTE modelling employing both 1D semi-empirical and 3D radiation-Magneto Hydrodynamic (RMHD) models. We study the basic formation mechanisms and derive a quintessential model atom that incorporates the essential atomic physics for the formation of the O I 135.56 nm line. This atomic model has 16 levels and describes recombination cascades through highly excited levels by effective recombination rates. The ionization balance O I/O II is set by the hydrogen ionization balance through charge exchange reactions. The emission in the O I 135.56 nm line is dominated by a recombination cascade and the line is optically thin. The Doppler shift of the maximum emission correlates strongly with the vertical velocity in its line forming region, which is typically located at 1.0 - 1.5 Mm height. The total intensity of the line emission is correlated with the square of the electron density. Since the O I 135.56 nm line is optically thin, the width of the emission line is a very good diagnostic of non-thermal velocities. We conclude that the O I 135.56 nm line is an excellent probe of the middle chromosphere, and compliments other powerful chromospheric diagnostics of IRIS such as the Mg II h & k lines and the C II lines around 133.5 nm.

preprint2015arXiv

Three-dimensional radiative transfer simulations of the scattering polarization of the hydrogen Ly$α$ line in a MHD model of the chromosphere-corona transition region

Probing the magnetism of the upper solar chromosphere requires measuring and modeling the scattering polarization produced by anisotropic radiation pumping in UV spectral lines. Here we apply PORTA (a novel radiative transfer code) to investigate the hydrogen Ly$α$ line in a 3D model of the solar atmosphere resulting from a state of the art MHD simulation. At full spatial resolution the linear polarization signals are very significant all over the solar disk, with a large fraction of the field of view showing line-center amplitudes well above the 1% level. Via the Hanle effect the line-center polarization signals are sensitive to the magnetic field of the model's transition region, even when its mean field strength is only 15 G. The breaking of the axial symmetry of the radiation field produces significant forward-scattering polarization in Ly$α$, without the need of an inclined magnetic field. Interestingly, the Hanle effect tends to decrease such forward-scattering polarization signals in most of the points of the field of view. When the spatial resolution is degraded, the line-center polarization of Ly$α$ drops below the 1% level, reaching values similar to those previously found in 1D semi-empirical models (i.e., up to about 0.5%). The center to limb variation of the spatially-averaged polarization signals is qualitatively similar to that found in 1D models, with the largest line-center amplitudes at $μ=\cosθ\approx 0.4$ ($θ$ being the heliocentric angle). These results are important, both for designing the needed space-based instrumentation and for a reliable interpretation of future observations of the Ly$α$ polarization.

preprint2015arXiv

What do iris observations of Mg II k tell us about the solar plage chromosphere?

We analyze observations from the Interface Region Imaging Spectrograph of the Mg II k line, the Mg II UV subordinate lines, and the O I 135.6 nm line to better understand the solar plage chromosphere. We also make comparisons with observations from the Swedish 1 m Solar Telescope of the Hα line, the Ca II 8542 line, and Solar Dynamics Observatory/Atmospheric Imaging Assembly observations of the coronal 19.3 nm line. To understand the observed Mg II profiles, we compare these observations to the results of numerical experiments. The single-peaked or flat-topped Mg II k profiles found in plage imply a transition region at a high column mass and a hot and dense chromosphere of about 6500 K. This scenario is supported by the observed large-scale correlation between moss brightness and filled-in profiles with very little or absent self-reversal. The large wing width found in plage also implies a hot and dense chromosphere with a steep chromospheric temperature rise. The absence of emission in the Mg II subordinate lines constrain the chromospheric temperature and the height of the temperature rise while the width of the O I 135.6 nm line sets a limit to the non-thermal velocities to around 7 km/s.

preprint2014arXiv

Detailed and simplified non-equilibrium helium ionization in the solar atmosphere

Helium ionization plays an important role in the energy balance of the upper chromosphere and transition region. Helium spectral lines are also often used as diagnostics of these regions. We carry out 1D radiation-hydrodynamics simulations of the solar atmosphere and find that the helium ionization is mostly set by photoionization and direct collisional ionization, counteracted by radiative recombination cascades. By introducing an additional recombination rate mimicking the recombination cascades, we construct a simplified 3 level helium model atom consisting of only the ground states. This model atom is suitable for modeling non-equilibrium helium ionization in 3D numerical models. We perform a brief investigation of the formation of the He I 10830 and He II 304 spectral lines. Both lines show non-equilibrium features that are not recovered with statistical equilibrium models, and caution should therefore be exercised when such models are used as a basis in the interpretation of observations.

preprint2014arXiv

Evidence of Non-Thermal Particles in Coronal Loops Heated Impulsively by Nanoflares

The physical processes causing energy exchange between the Sun's hot corona and its cool lower atmosphere remain poorly understood. The chromosphere and transition region (TR) form an interface region between the surface and the corona that is highly sensitive to the coronal heating mechanism. High resolution observations with the Interface Region Imaging Spectrograph (IRIS) reveal rapid variability (about 20 to 60 seconds) of intensity and velocity on small spatial scales at the footpoints of hot dynamic coronal loops. The observations are consistent with numerical simulations of heating by beams of non-thermal electrons, which are generated in small impulsive heating events called "coronal nanoflares". The accelerated electrons deposit a sizable fraction of their energy in the chromosphere and TR. Our analysis provides tight constraints on the properties of such electron beams and new diagnostics for their presence in the nonflaring corona.

preprint2014arXiv

The effect of isotopic splitting on the bisector and inversions of the solar Ca II 854.2 nm line

The Ca II 854.2 nm spectral line is a common diagnostic of the solar chromosphere. The average line profile shows an asymmetric core, and its bisector shows a characteristic inverse-C shape. The line actually consists of six components with slightly different wavelengths depending on the isotope of calcium. This isotopic splitting of the line has been taken into account in studies of non-solar stars, but never for the Sun. We performed non-LTE radiative transfer computations from three models of the solar atmosphere and show that the asymmetric line-core and inverse C-shape of the bisector of the 854.2 nm line can be explained by isotopic splitting. We confirm this finding by analysing observations and showing that the line asymmetry is present irrespective of conditions in the solar atmosphere. Finally, we show that inversions based on the Ca II 854.2 nm line should take the isotopic splitting into account, otherwise the inferred atmospheres will contain erroneous velocity gradients and temperatures.

preprint2013arXiv

A Detailed Comparison Between The Observed and Synthesized Properties of a Simulated Type ii Spicule

We performed a 3D radiative MHD simulation of the solar atmosphere. This simulation shows a jet-like feature that shows similarities to the type II spicules observed for the first time with Hinode. Rapid Blueshifted Events (RBEs) on the solar disk are associated with these spicules. Observational results suggest they may contribute significantly in supplying the corona with hot plasma. We perform a detailed comparison of the properties of the simulated jet with those of type II spicules (observed with Hinode) and RBEs (with ground-based instruments). We analyze variety of synthetic emission and absorption lines from the simulations including chromospheric Ca II and Ha to TR and coronal temperatures (10E4 to several 10E6K). We compare their synthetic intensities, line profiles, Doppler shifts, line widths and asymmetries with observations from Hinode/SOT and EIS, SOHO/SUMER, SST and SDO/AIA. Many properties of the synthetic observables resemble the observations, and we describe in detail the physical processes that lead to these observables. Detailed analysis of the synthetic observables provides insight into how observations should be analyzed to derive information about physical variables in such a dynamic event. For example, we find that LOS superposition in the optically thin atmosphere requires the combination of Doppler shifts and spectral line asymmetry to determine the velocity in the jet. Other properties differ from the observations, especially in the chromospheric lines. The mass density of the part of the spicule with a chromospheric temperature is too low to produce significant opacity in chromospheric lines. These and other discrepancies are described in detail, and we discuss which mechanisms and physical processes may need to be included in the MHD simulations to mimic the thermodynamic processes of the chromosphere and corona, in particular to reproduce type II spicules.

preprint2013arXiv

The formation of IRIS diagnostics. III. Near-ultraviolet Spectra and Images

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.

preprint2012arXiv

Approximations for radiative cooling and heating in the solar chromosphere

Context. The radiative energy balance in the solar chromosphere is dominated by strong spectral lines that are formed out of LTE. It is computationally prohibitive to solve the full equations of radiative transfer and statistical equilibrium in 3D time dependent MHD simulations. Aims. To find simple recipes to compute the radiative energy balance in the dominant lines under solar chromospheric conditions. Methods. We use detailed calculations in time-dependent and 2D MHD snapshots to derive empirical formulae for the radiative cooling and heating. Results. The radiative cooling in neutral hydrogen lines and the Lyman continuum, the H and K and intrared triplet lines of singly ionized calcium and the h and k lines of singly ionized magnesium can be written as a product of an optically thin emission (dependent on temperature), an escape probability (dependent on column mass) and an ionization fraction (dependent on temperature). In the cool pockets of the chromosphere the same transitions contribute to the heating of the gas and similar formulae can be derived for these processes. We finally derive a simple recipe for the radiative heating of the chromosphere from incoming coronal radiation. We compare our recipes with the detailed results and comment on the accuracy and applicability of the recipes.

preprint2012arXiv

Investigating the reliability of coronal emission measure distribution diagnostics using 3D radiative MHD simulations

Determining the temperature distribution of coronal plasmas can provide stringent constraints on coronal heating. Current observations with the Extreme ultraviolet Imaging Spectrograph onboard Hinode and the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory provide diagnostics of the emission measure distribution (EMD) of the coronal plasma. Here we test the reliability of temperature diagnostics using 3D radiative MHD simulations. We produce synthetic observables from the models, and apply the Monte Carlo Markov chain EMD diagnostic. By comparing the derived EMDs with the "true" distributions from the model we assess the limitations of the diagnostics, as a function of the plasma parameters and of the signal-to-noise of the data. We find that EMDs derived from EIS synthetic data reproduce some general characteristics of the true distributions, but usually show differences from the true EMDs that are much larger than the estimated uncertainties suggest, especially when structures with significantly different density overlap along the line-of-sight. When using AIA synthetic data the derived EMDs reproduce the true EMDs much less accurately, especially for broad EMDs. The differences between the two instruments are due to the: (1) smaller number of constraints provided by AIA data, (2) broad temperature response function of the AIA channels which provide looser constraints to the temperature distribution. Our results suggest that EMDs derived from current observatories may often show significant discrepancies from the true EMDs, rendering their interpretation fraught with uncertainty. These inherent limitations to the method should be carefully considered when using these distributions to constrain coronal heating.

preprint2012arXiv

Non-local thermodynamic equilibrium inversions from a 3D MHD chromospheric model

The structure of the solar chromosphere is believed to be governed by magnetic fields, even in quiet-Sun regions that have a relatively weak photospheric field. During the past decade inversion methods have emerged as powerful tools for analyzing the chromosphere of active regions. The applicability of inversions to infer the stratification of the physical conditions in a dynamic 3D solar chromosphere has not yet been studied in detail. This study aims to establish the diagnostic capabilities of non-local thermodynamical equilibrium (NLTE) inversion techniques of Stokes profiles induced by the Zeeman effect in the Ca II 8542 line. We computed the Ca II atomic level populations in a snapshot from a 3D radiation-MHD simulation of the quiet solar atmosphere in non-LTE using the 3D radiative transfer code Multi3d. These populations were used to compute synthetic full-Stokes profiles in the Ca II 8542 line using 1.5D radiative transfer and the inversion code Nicole. The profiles were then spectrally degraded to account for finite filter width and Gaussian noise was added to account for finite photon flux. These profiles were inverted using Nicole and the results were compared with the original model atmosphere. Our NLTE inversions applied to quiet-Sun synthetic observations provide reasonably good estimates of the chromospheric magnetic field, line-of-sight velocities and somewhat less accurate, but still very useful, estimates of the temperature. Three dimensional scattering of photons cause cool pockets in the chromosphere to be invisible in the line profile and consequently they are also not recovered by the inversions. To successfully detect Stokes linear polarization in this quiet snapshot, a noise level below 10^{-3.5} is necessary.

preprint2012arXiv

Quantifying Spicules

Understanding the dynamic solar chromosphere is fundamental in solar physics. Spicules are an important feature of the chromosphere, connecting the photosphere to the corona, potentially mediating the transfer of energy and mass. The aim of this work is to study the properties of spicules over different regions of the sun. Our goal is to investigate if there is more than one type of spicules, and how spicules behave in the quiet sun, coronal holes, and active regions. We make use of high-cadence and high-spatial resolution Ca II H observations taken by Hinode/SOT. Making use of a semi-automated detection algorithm, we self-consistently track and measure the properties of 519 spicules over different regions. We find clear evidence of two types of spicules. Type I spicules show a rise and fall and have typical lifetimes of 150-400 s and maximum ascending velocities of 15-40 km/s, while type II spicules have shorter lifetimes of 50-150 s, faster velocities of 30-110 km/s, and are not seen to fall down, but rather fade at around their maximum length. Type II spicules are the most common, seen in quiet sun and coronal holes. Type I spicules are seen mostly in active regions. There are regional differences between quiet sun and coronal hole spicules, likely attributable to the different field configurations. The properties of type II spicules are consistent with published results of Rapid Blueshifted Events (RBEs), supporting the hypothesis that RBEs are their disk counterparts. For type I spicules we find the relations between their properties to be consistent with a magnetoacoustic shock wave driver, and with dynamic fibrils as their disk counterpart. The driver of type II spicules remains unclear from limb observations.

preprint2012arXiv

The formation of the Halpha line in the solar chromosphere

We use state-of-the-art radiation-MHD simulations and 3D non-LTE radiative transfer computations to investigate \Halpha\ line formation in the solar chromosphere and apply the results of this investigation to develop the potential of \Halpha\ as diagnostic of the chromosphere. We show that one can accurately model \Halpha\ line formation assuming statistical equilibrium and complete frequency redistribution provided the computation of the model atmosphere included non-equilibrium ionization of hydrogen, and the Lyman-$α$ and Lyman-$β$ line profiles are described by Doppler profiles. We find that 3D radiative transfer is essential in modeling hydrogen lines due to the low photon destruction probability in \Halpha. The \Halpha\ opacity in the upper chromosphere is mainly sensitive to the mass density and only weakly sensitive to temperature. We find that the \Halpha\ line-core intensity is correlated with the average formation height: the larger the average formation height, the lower the intensity. The line-core width is a measure of the gas temperature in the line-forming region. The fibril-like dark structures seen in \Halpha\ line-core images computed from our model atmosphere are tracing magnetic field lines. These structures are caused by field-aligned ridges of enhanced chromospheric mass density that raise their average formation height, and therefore makes them appear dark against their deeper-formed surroundings. We compare with observations, and find that the simulated line-core widths are very similar to the observed ones, without the need for additional microturbulence.

preprint2012arXiv

The Hanle effect of Ly$α$ in an MHD model of the Solar Transition Region

In order to understand the heating of the solar corona it is crucial to obtain empirical information on the magnetic field in its lower boundary (the transition region). To this end, we need to measure and model the linear polarization produced by scattering processes in strong UV lines, such as the hydrogen Ly$α$ line. The interpretation of the observed Stokes profiles will require taking into account that the outer solar atmosphere is highly structured and dynamic, and that the height of the transition region may well vary from one place in the atmosphere to another. Here we report on the Ly$α$ scattering polarization signals we have calculated in a realistic model of an enhanced network region, resulting from a state-of-the-art radiation MHD simulation. This model is characterized by spatially complex variations of the physical quantities at transition region heights. The results of our investigation lead us to emphasize that scattering processes in the upper solar chromosphere should indeed produce measurable linear polarization in Ly$α$. More importantly, we show that via the Hanle effect the model's magnetic field produces significant changes in the emergent $Q/I$ and $U/I$ profiles. Therefore, we argue that by measuring the polarization signals produced by scattering processes and the Hanle effect in Ly$α$ and contrasting them with those computed in increasingly realistic atmospheric models, we should be able to decipher the magnetic, thermal and dynamic structure of the upper chromosphere and transition region of the Sun.

preprint2011arXiv

Non-equilibrium calcium ionisation in the solar atmosphere

Our aim is to determine the dominant processes and timescales for the ionisation equilibrium of calcium under solar chromospheric conditions. The study is based on numerical simulations with the RADYN code, which includes hydrodynamics, radiative transfer, and a detailed non-equilibrium treatment of hydrogen, calcium, and helium. The simulations are characterised by upwards propagating shock waves, which cause strong temperature fluctuations and variations of the ionisation degree of calcium. The passage of a hot shock front leads to a strong net ionisation of Ca II, rapidly followed by net recombination. The relaxation timescale of the Ca ionisation state is found to be of the order of a few seconds at the top of the photosphere and 10 to 30 s in the upper chromosphere. Generally, the timescales are significantly reduced in the wakes of hot shock fronts. The timescales can be reliably determined from a simple analysis of the eigenvalues of the transition rate matrix. The timescales are dominated by the radiative recombination from Ca III into the metastable Ca II energy levels of the 4d 2D term. These transitions depend strongly on the density of free electrons and therefore on the (non-equilibrium) ionisation degree of hydrogen, which is the main electron donor. The ionisation/recombination timescales derived here are too long for the assumption of an instantaneous ionisation equilibrium to be valid and, on the other hand, are not long enough to warrant an assumption of a constant ionisation fraction. Fortunately, the ionisation degree of Ca II remains small in the height range, where the cores of the H, K, and the infrared triplet lines are formed. We conclude that the difference due to a detailed treatment of Ca ionisation has only negligible impact on the modelling of spectral lines of Ca II and the plasma properties under the conditions in the quiet solar chromosphere.

preprint2011arXiv

On the minimum temperature of the quiet solar chromosphere

aims: We aim to provide an estimate of the minimum temperature of the quiet solar chromosphere. methods: We perform a 2D radiation-MHD simulation spanning the upper convection zone to the lower corona. The simulation includes non-LTE radiative transfer and {an equation-of-state that includes non-equilibrium ionization of hydrogen and non-equilibrium H_2 molecule formation}. We analyze the reliability of the various assumptions made in our model in order to assess the realism of the simulation. results: Our simulation contains pockets of cool gas with down to 1660 K from 1 Mm up to 3.2 Mm height. It overestimates the radiative heating, and contains non-physical heating below 1660 K. Therefore we conclude that cool pockets in the quiet solar chromosphere might have even lower temperatures than in the simulation, provided that there exist areas in the chromosphere without significant magnetic heating. We suggest off-limb molecular spectroscopy to look for such cool pockets and 3D simulations including a local dynamo and a magnetic carpet to investigate Joule heating in the quiet chromosphere.

preprint2011arXiv

Solar velocity references from 3D HD photospheric models

The measurement of Doppler velocities in spectroscopic solar observations requires a reference for the local frame of rest. The rotational and radial velocities of the Earth and the rotation of the Sun introduce velocity offsets in the observations. Normally, good references for velocities are missing (e.g. telluric lines), especially in filter-based spectropolarimetric observations. We determine an absolute reference for line-of-sight velocities measured from solar observations for any heliocentric angle, calibrating the convective line shift of spatially-averaged profiles on quiet sun from a 3D hydrodynamical simulation. This method works whenever there is quiet sun in the field-of-view, and it has the advantage of being relatively insensitive to uncertainties in the atomic data. We carry out radiative transfer computations in LTE for selected C I and Fe I lines, whereas the Ca II infrared lines are synthesized in non-LTE. The degree of asymmetry and the line shifts, however, show a clear dependence on the heliocentric angle and the properties of the lines. The profiles at μ= 1 are compared with observed profiles to prove their reliability, and they are tested against errors induced by the LTE calculations, inaccuracies in the atomic data and the 3D simulation. Theoretical quiet-sun profiles of lines commonly used by solar observers are provided to the community. Those can be used as absolute references for line-of-sight velocities. The limb effect is produced by the projection of the 3D atmosphere along the line of sight. We estimate the precision of the disk-center line shifts to be approximately 50 m s^-1, but the off-center profiles remain to be tested against observations.

preprint2011arXiv

The stellar atmosphere simulation code Bifrost

Context: Numerical simulations of stellar convection and photospheres have been developed to the point where detailed shapes of observed spectral lines can be explained. Stellar atmospheres are very complex, and very different physical regimes are present in the convection zone, photosphere, chromosphere, transition region and corona. To understand the details of the atmosphere it is necessary to simulate the whole atmosphere since the different layers interact strongly. These physical regimes are very diverse and it takes a highly efficient massively parallel numerical code to solve the associated equations. Aims: The design, implementation and validation of the massively parallel numerical code Bifrost for simulating stellar atmospheres from the convection zone to the corona. Methods: The code is subjected to a number of validation tests, among them the Sod shock tube test, the Orzag-Tang colliding shock test, boundary condition tests and tests of how the code treats magnetic field advection, chromospheric radiation, radiative transfer in an isothermal scattering atmosphere, hydrogen ionization and thermal conduction. Results: Bifrost completes the tests with good results and shows near linear efficiency scaling to thousands of computing cores.

preprint2010arXiv

Chromospheric heating and structure as determined from high resolution 3D simulations

We have performed 3D radiation MHD simulations extending from the convection zone to the corona covering a box 16 Mm$^3$ at 32 km spatial resolution. The simulations show very fine structure in the chromosphere with acoustic shocks interacting with the magnetic field. Magnetic flux concentrations have a temperature lower than the surroundings in the photosphere but higher in the low chromosphere. The heating is there mostly through ohmic dissipation preferentially at the edges of the flux concentrations. The magnetic field is often wound up around the flux concentrations. When acoustic waves travel up along the field this topology leads to swirling motions seen in chromospheric diagnostic lines such as the calcium infrared triplet.

preprint2010arXiv

On red shifs in the transition region and corona

We present evidence that transition region red-shifts are naturally produced in episodically heated models where the average volumetric heating scale height lies between that of the chromospheric pressure scale height of 200 km and the coronal scale height of 50 Mm. In order to do so we present results from 3d MHD models spanning the upper convection zone up to the corona, 15 Mm above the photosphere. Transition region and coronal heating in these models is due both the stressing of the magnetic field by photospheric and convection `zone dynamics, but also in some models by the injection of emerging magnetic flux.

preprint2010arXiv

On the solar chromosphere observed at the limb with Hinode

Broad-band images in the Ca II H line, from the BFI instrument on the Hinode spacecraft, show emission from spicules emerging from and visible right down to the observed limb. Surprisingly, little absorption of spicule light is seen along their lengths. We present formal solutions to the transfer equation for given (ad-hoc) source functions, including a stratified chromosphere from which spicules emanate. The model parameters are broadly compatible with earlier studies of spicules. The visibility of Ca II spicules down to the limb in Hinode data seems to require that spicule emission be Doppler shifted relative to the stratified atmosphere, either by supersonic turbulent or organized spicular motion. The non-spicule component of the chromosphere is almost invisible in the broad band BFI data, but we predict that it will be clearly visible in high spectral resolution data. Broad band Ca II H limb images give the false impression that the chromosphere is dominated by spicules. Our analysis serves as a reminder that the absence of a signature can be as significant as its presence.

preprint2010arXiv

SICStus Prolog -- the first 25 years

SICStus Prolog has evolved for nearly 25 years. This is an appropriate point in time for revisiting the main language and design decisions, and try to distill some lessons. SICStus Prolog was conceived in a context of multiple, conflicting Prolog dialect camps and a fledgling standardization effort. We reflect on the impact of this effort and role model implementations on our development. After summarizing the development history, we give a guided tour of the system anatomy, exposing some designs that were not published before. We give an overview of our new interactive development environment, and describe a sample of key applications. Finally, we try to identify key good and not so good design decisions.

preprint2009arXiv

The Temperature Diagnostic Properties of the MG I 457.1 NM Line

We analyze the important formation processes for the Mg I 457.1 nm line. This line is an intercom- bination line and the source function is close to the local thermodynamic equilibrium (LTE) value. The strong coupling to the local temperature and the relatively high population of the lower level (the ground state of Mg I) makes this line an ideal candidate for temperature diagnostics in the lower chromosphere/temperature minimum region. Linking the temperature probed to an absolute phys- ical height is non trivial because of Non-LTE ionization. We analyze the Non-LTE effects and find that photo-ionization from the lower energy levels together with cascading collisional recombination dominate the ionization balance. Taking properly into account the line-blanketing in the UV is essen- tial for obtaining the right photoionization rates. The identification of the main Non-LTE effects in the line allows us to construct a quintessential model atom, ideal for computationally demanding tasks, e.g. full 3D and/or time-dependent radiative transfer. Furthermore we analyze the diagnostic potential to temperature of this line in solar-like atmospheres, by synthesizing the line from a series of parametrized atmospheric models. These models have been constructed with fixed effective temper- ature, but with a variable heat term in the energy equation to obtain a chromospheric temperature rise at different heights. We conclude that the line has a significant potential in the diagnostics of the lower chromosphere temperature structure, especially for cooler atmospheres, such as sunspots.

preprint2007arXiv

A Tale Of Two Spicules: The Impact of Spicules on the Magnetic Chromosphere

We use high-resolution observations of the Sun in Ca II H 3968 A from the Solar Optical Telescope on Hinode to show that there are at least two types of spicules that dominate the structure of the magnetic solar chromosphere. Both types are tied to the relentless magnetoconvective driving in the photosphere, but have very different dynamic properties. ``Type-I'' spicules are driven by shock waves that form when global oscillations and convective flows leak into the upper atmosphere along magnetic field lines on 3-7 minute timescales. ``Type-II'' spicules are much more dynamic: they form rapidly (in ~10s), are very thin (<200km wide), have lifetimes of 10-150s (at any one height) and seem to be rapidly heated to (at least) transition region temperatures, sending material through the chromosphere at speeds of order 50-150 km/s. The properties of Type II spicules suggest a formation process that is a consequence of magnetic reconnection, typically in the vicinity of magnetic flux concentrations in plage and network. Both types of spicules are observed to carry Alfven waves with significant amplitudes of order 20 km/s.

preprint2007arXiv

Can High Frequency Acoustic Waves Heat the Quiet Sun Chromosphere?

We use Hinode/SOT Ca II H-line and blue continuum broadband observations to study the presence and power of high frequency acoustic waves at high spatial resolution. We find that there is no dominant power at small spatial scales; the integrated power using the full resolution of Hinode (0.05'' pixels, 0.16'' resolution) is larger than the power in the data degraded to 0.5'' pixels (TRACE pixel size) by only a factor of 1.2. At 20 mHz the ratio is 1.6. Combining this result with the estimates of the acoustic flux based on TRACE data of Fossum & Carlsson (2006), we conclude that the total energy flux in acoustic waves of frequency 5-40 mHz entering the internetwork chromosphere of the quiet Sun is less than 800 W m$^{-2}$, inadequate to balance the radiative losses in a static chromosphere by a factor of five.

preprint2007arXiv

On connecting the dynamics of the chromosphere and transition region with Hinode SOT and EIS

We use coordinated Hinode SOT/EIS observations that include high-resolution magnetograms, chromospheric and TR imaging and TR/coronal spectra in a first test to study how the dynamics of the TR are driven by the highly dynamic photospheric magnetic fields and the ubiquitous chromospheric waves. Initial analysis shows that these connections are quite subtle and require a combination of techniques including magnetic field extrapolations, frequency-filtered time-series and comparisons with synthetic chromospheric and TR images from advanced 3D numerical simulations. As a first result, we find signatures of magnetic flux emergence as well as 3 and 5 mHz wave power above regions of enhanced photospheric magnetic field in both chromospheric, transition region and coronal emission.

preprint2005arXiv

Non-LTE neutral carbon spectral line formation in late-type stars

We present non-Local Thermodynamic Equilibrium (non-LTE) calculations for neutral carbon spectral line formation, carried out for a grid of model atmospheres covering the range of late-type stars. The results of our detailed calculations suggest that the carbon non-LTE corrections in these stars are higher than usually adopted, remaining substantial even at low metallicity. For the most metal-poor stars in the sample of Akerman et al. (2004), the non-LTE abundance corrections are of the order of -0.35...-0.45 dex (when neglecting H collisions). Applying our results to those observations, the apparent [C/O] upturn seen in their LTE analysis is no longer present, thus revealing no need to invoke contributions from Pop. III stars to the carbon nucleosynthesis.

preprint2002arXiv

Dynamic Hydrogen Ionization

We investigate the ionization of hydrogen in a dynamic Solar atmosphere. We show that the time scale for ionization/recombination can be estimated from the eigenvalues of a modified rate matrix where the optically thick Lyman transitions that are in detailed balance have been excluded. We find that the time scale for ionization/recombination is dominated by the slow collisional leakage from the ground state to the first excited state. Throughout the chromosphere the time scale is long ($10^3$-$10^5$ s), except in shocks where the increased temperature and density shorten the time scale for ionization/recombination, especially in the upper chromosphere. Because the relaxation time scale is much longer than dynamic time scales, hydrogen ionization does not have time to reach its equilibrium value and its fluctuations are much smaller than the variation of its statistical equilibrium value appropriate for the instantaneous conditions. The ionization state tends to represent the higher temperature of the shocks, and the mean electron density is up to a factor of six higher than the electron density calculated in statistical equilibrium from the mean atmosphere. The simulations show that a static picture and a dynamic picture of the chromosphere are fundamentally different and that time variations are crucial for our understanding of the chromosphere itself and the spectral features formed there.

Mats Carlsson

What is connected

Connect this record

See the researcher in context

Building this map preview

63 published item(s)

Chromosphere of the quiet sun -- I. Shock and current-sheet dynamics and heating

An approximate recipe of chromospheric radiative losses for solar flares

Effects of spatial resolution on inferences of atmospheric quantities from simulations

The Atmospheric Response to High Nonthermal Electron Beam Fluxes in Solar Flares. II. Hydrogen Broadening Predictions for Solar Flare Observations with the Daniel K. Inouye Solar Telescope

Mapping Solar Magnetic Fields from the Photosphere to the Base of the Corona

Probing the physics of the solar atmosphere with the Multi-slit Solar Explorer (MUSE): I. Coronal Heating

Probing the Physics of the Solar Atmosphere with the Multi-slit Solar Explorer (MUSE): II. Flares and Eruptions

Characterisation of shock wave signatures at millimetre wavelengths from Bifrost simulations

Chromospheric Bubbles in Solar Flares

On the formation of Lyman $β$ and the O I 1027 and 1028 Å spectral lines

The Sun at millimeter wavelengths I. Introduction to ALMA Band 3 observations

Cosmic-Ray Interactions in the Solar Atmosphere

Ion-neutral interactions and non-equilibrium ionization in the solar chromosphere

Chromospheric and Coronal Wave Generation in a Magnetic Flux Sheath

Emergence of granular-sized magnetic bubbles through the solar atmosphere. III. The path to the transition region

Modeling Repeatedly Flaring $δ$ Sunspots

The Appearance of Spicules in High Resolution Observations of Ca II H and H-alpha

The cause of spatial structure in solar He I 1083 nm multiplet images

A publicly available simulation of an enhanced network region of the Sun

A Unified Computational Model for Solar and Stellar Flares

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares. II. Inclusion of Radiative Transfer with RADYN

Ellerman bombs at high resolution III. Simultaneous observations with IRIS and SST

Evidence for a Transition Region response to penumbral microjets in sunspots

First High-resolution Spectroscopic Observations of an Erupting Prominence Within a Coronal Mass Ejection by the Interface Region Imaging Spectrograph (IRIS)

Internetwork chromospheric bright grains observed with IRIS

New Insights into White-Light Flare Emission from Radiative-Hydrodynamic Modeling of a Chromospheric Condensation

Non-equilibrium helium ionization in an MHD simulation of the solar atmosphere

Numerical Simulations of Coronal Heating through Footpoint Braiding

Observed Variability of the Solar Mg II h Spectral Line

On fibrils and field lines: The nature of H$α$ fibrils in the solar chromosphere

Probing the Sun with ALMA: observations and simulations

The formation of IRIS diagnostics V. A quintessential model atom of C II and general formation properties of the C II lines at 133.5 nm

The formation of IRIS diagnostics VI. The Diagnostic Potential of the C II Lines at 133.5 nm in the Solar Atmosphere

The formation of IRIS diagnostics VIII. IRIS observations in the C II 133.5 nm multiplet

The Formation of IRIS diagnostics. IV. The Mg II triplet lines as a new diagnostic for lower chromospheric heating

The Formation of IRIS Diagnostics. VII. The Formation of the OI 135.56 nm Line in the Solar Atmosphere

Three-dimensional radiative transfer simulations of the scattering polarization of the hydrogen Ly$α$ line in a MHD model of the chromosphere-corona transition region

What do iris observations of Mg II k tell us about the solar plage chromosphere?

Detailed and simplified non-equilibrium helium ionization in the solar atmosphere

Evidence of Non-Thermal Particles in Coronal Loops Heated Impulsively by Nanoflares

The effect of isotopic splitting on the bisector and inversions of the solar Ca II 854.2 nm line

A Detailed Comparison Between The Observed and Synthesized Properties of a Simulated Type ii Spicule

The formation of IRIS diagnostics. III. Near-ultraviolet Spectra and Images

Approximations for radiative cooling and heating in the solar chromosphere

Investigating the reliability of coronal emission measure distribution diagnostics using 3D radiative MHD simulations

Non-local thermodynamic equilibrium inversions from a 3D MHD chromospheric model

Quantifying Spicules

The formation of the Halpha line in the solar chromosphere

The Hanle effect of Ly$α$ in an MHD model of the Solar Transition Region

Non-equilibrium calcium ionisation in the solar atmosphere

On the minimum temperature of the quiet solar chromosphere

Solar velocity references from 3D HD photospheric models

The stellar atmosphere simulation code Bifrost

Chromospheric heating and structure as determined from high resolution 3D simulations

On red shifs in the transition region and corona

On the solar chromosphere observed at the limb with Hinode

SICStus Prolog -- the first 25 years

The Temperature Diagnostic Properties of the MG I 457.1 NM Line

A Tale Of Two Spicules: The Impact of Spicules on the Magnetic Chromosphere

Can High Frequency Acoustic Waves Heat the Quiet Sun Chromosphere?

On connecting the dynamics of the chromosphere and transition region with Hinode SOT and EIS

Non-LTE neutral carbon spectral line formation in late-type stars

Dynamic Hydrogen Ionization