Engineering Papers⌕ Search

Engineering topics

Carlsson, Mats

Publications and source records attributed to Carlsson, Mats.

At least 19 records

Modeling Mg II During Solar Flares. II. Nonequilibrium Effects

To extract the information that the Mg II NUV spectra (observed by the Interface Region Imaging Spectrograph) carry about the chromosphere during solar flares, and to validate models of energy transport via model–data comparison, forward modeling is required. The assumption of statistical equilibrium (SE) is typically used to obtain the atomic level populations from snapshots of flare atmospheres, due to computational necessity. However, it is possible that relying on SE could lead to spurious results. We compare solving the atomic level populations via SE versus a non-equilibrium (NEQ) time-dependent approach. This was achieved using flare simulations from RADYN alongside the minority species version MS_RADYN from which the time-dependent Mg II atomic level populations and radiation transfer were computed in complete frequency redistribution. The impacts on the emergent profiles, light curves, line ratios, and formation heights are discussed. In summary we note that NEQ effects during flares are typically important only in the initial stages and for a short period following the cessation of the energy injection. An analysis of the timescales of ionization equilibrium reveals that for most of the duration of the flare, when the temperatures and densities are sufficiently enhanced, the relaxation timescales are short (τrelax < 0.1 s), so that the equilibrium solution is an adequate approximation. These effects vary with the size of the flare, however. In weaker flares, effects can be more pronounced. We recommend that NEQ effects be considered when possible but that SE is sufficient at most stages of the flare.

Kerr, Graham S.↗

Modeling Mg II During Solar Flares. I. Partial Frequency Redistribution, Opacity, and Coronal Irradiation

The Interface Region Imaging Spectrograph has routinely observed the flaring Mg II near-ultraviolet (NUV) spectrum, offering excellent diagnostic potential and a window into the location of energy deposition. A number of studies have forward-modeled both the general properties of these lines and specific flare observations. Generally these have forward-modeled radiation via post-processing of snapshots from hydrodynamic flare simulations through radiation transfer codes. There has, however, not been a study of how the physics included in these radiation transport codes affects the solution. A baseline setup for forward-modeling Mg II in flares is presented and contrasted with approaches that add or remove complexity. It is shown for Mg II that (1) partial frequency distribution (PRD) is still required during flare simulations despite the increased densities; (2) using full angle-dependent PRD affects the solution but takes significantly longer to process a snapshot; (3) including Mg I in non-LTE (NLTE) results in negligible differences to the Mg II lines but does affect the NUV quasi-continuum; (4) only hydrogen and Mg II need to be included in NLTE; (5) ideally the nonequilibrium hydrogen populations, with nonthermal collisional rates, should be used rather than the statistical equilibrium populations; (6) an atom consisting of only the ground state, h and k upper levels, and continuum level is insufficient to model the resonance lines; and (7) irradiation from a hot, dense flaring transition region can affect the formation of Mg II. We discuss modifications to the RH code allowing straightforward inclusion of the transition region and coronal irradiation in flares.

Kerr, Graham S.↗

SI IV Resonance Line Emission during Solar Flares: Non-LTE, Non-equilibrium, Radiation Transfer Simulations

The Interface Region Imaging Spectrograph routinely observes the Si IV resonance lines. When analyzing quiescent observations of these lines, it has typically been assumed that they form under optically thin conditions. This is likely valid for the quiescent Sun, but this assumption has also been applied to the more extreme flaring scenario. We used 36 electron-beam-driven radiation hydrodynamic solar flare simulations, computed using the RADYN code, to probe the validity of this assumption. Using these simulated atmospheres, we solved the radiation transfer equations to obtain the non-LTE, nonequilibrium populations, line profiles, and opacities for a model silicon atom, including charge exchange processes. This was achieved using the "minority species" version of RADYN. The inclusion of charge exchange resulted in a substantial fraction of Si IV at cooler temperatures than those predicted by ionization equilibrium. All simulations with an injected energy flux F> 5 x 10(exp 10) erg cm(exp -2) s(exp -1) resulted in optical depth effects on the Si IV emission, with differences in both intensity and line shape compared to the optically thin calculation. Weaker flares (down to F ~ 10(exp 9) erg cm(exp -2) s(exp -1)) also resulted in Si IV emission forming under optically thick conditions, depending on the other beam parameters. When opacity was significant, the atmospheres generally had column masses in excess of 5 x 10-6 g cm-2 over the temperature range 40-100 kK, and the Si IV formation temperatures were between 30 and 60 kK. We urge caution when analyzing Si IV flare observations, or when computing synthetic emission without performing a full radiation transfer calculation.

Kerr, Graham S.↗

Two-dimensional Radiative Magnetohydrodynamic Simulations of Partial Ionization in the Chromosphere. II. Dynamics and Energetics of the Low Solar Atmosphere

We investigate the effects of interactions between ions and neutrals on the chromosphere and overlying corona using 2.5D radiative MHD simulations with the Bifrost code. We have extended the code capabilities implementing ion–neutral interaction effects using the generalized Ohm’s law, i.e., we include the Hall term and the ambipolar diffusion (Pedersen dissipation) in the induction equation. Our models span from the upper convection zone to the corona, with the photosphere, chromosphere, and transition region partially ionized. Our simulations reveal that the interactions between ionized particles and neutral particles have important consequences for the magnetothermodynamics of these modeled layers: (1) ambipolar diffusion increases the temperature in the chromosphere; (2) sporadically the horizontal magnetic field in the photosphere is diffused into the chromosphere, due to the large ambipolar diffusion; (3) ambipolar diffusion concentrates electrical currents, leading to more violent jets and reconnection processes, resulting in (3a) the formation of longer and faster spicules, (3b) heating of plasma during the spicule evolution, and (3c) decoupling of the plasma and magnetic field in spicules. Our results indicate that ambipolar diffusion is a critical ingredient for understanding the magnetothermodynamic properties in the chromosphere and transition region. The numerical simulations have been made publicly available, similar to previous Bifrost simulations. This will allow the community to study realistic numerical simulations with a wider range of magnetic field configurations and physics modules than previously possible.

79 ASTRONOMY AND ASTROPHYSICS↗

The Formation of IRIS Diagnostics. IX. The Formation of the C i 135.58 NM Line in the Solar Atmosphere

The C i 135.58 nm line is located in the wavelength range of NASA’s Interface Region Imaging Spectrograph ( IRIS ) small explorer mission. We study the formation and diagnostic potential of this line by means of non local-thermodynamic-equilibrium modeling, employing both 1D and 3D radiation-magnetohydrodynamic models. The C i/C ii ionization balance is strongly influenced by photoionization by Ly α emission. The emission in the C i 135.58 nm line is dominated by a recombination cascade and the line forming region is optically thick. The Doppler shift of the line correlates strongly with the vertical velocity in its line forming region, which is typically located at 1.5 Mm height. With IRIS , the C i 135.58 nm line is usually observed together with the O i 135.56 nm line, and from the Doppler shift of both lines, we obtain the velocity difference between the line forming regions of the two lines. From the ratio of the C i/O i line core intensity, we can determine the distance between the C i and the O i forming layers. Combined with the velocity difference, the velocity gradient at mid-chromospheric heights can be derived. The C i/O i total intensity line ratio is correlated with the inverse of the electron density in the mid-chromosphere. We conclude that the C i 135.58 nm line is an excellent probe of the middle chromosphere by itself, and together with the O i 135.56 nm line the two lines provide even more information, which complements other powerful chromospheric diagnostics of IRIS such as the Mg ii h and k lines and the C ii lines around 133.5 nm.

79 ASTRONOMY AND ASTROPHYSICS↗

RADYN Simulations of Non-thermal and Thermal Models of Ellerman Bombs

Ellerman bombs (EBs) are brightenings in the H α line wings that are believed to be caused by magnetic reconnection in the lower atmosphere. To study the response and evolution of the chromospheric line profiles, we perform radiative hydrodynamic simulations of EBs using both non-thermal and thermal models. Overall, these models can generate line profiles that are similar to observations. However, in non-thermal models we find dimming in the H α line wings and continuum when the heating begins, while for the thermal models dimming occurs only in the H α line core, and with a longer lifetime. This difference in line profiles can be used to determine whether an EB is dominated by non-thermal heating or thermal heating. In our simulations, if a higher heating rate is applied, then the H α line will be unrealistically strong and there are still no clear UV burst signatures.

79 ASTRONOMY AND ASTROPHYSICS↗

Hydrogen Balmer Line Broadening in Solar and Stellar Flares

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

79 ASTRONOMY AND ASTROPHYSICS↗

Hydrogen Balmer Line Broadening in Solar and Stellar Flares

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

Kowalski, Adam F.↗

The Atmospheric Response to High Nonthermal Electron Beam Fluxes in Solar Flares. I. Modeling the Brightest NUV Footpoints in the X1 Solar Flare of 2014 March 29

The 2014 March 29 X1 solar flare (SOL20140329T17:48) produced bright continuum emission in the far- and near-ultraviolet (NUV) and highly asymmetric chromospheric emission lines, providing long-sought constraints on the heating mechanisms of the lower atmosphere in solar flares. We analyze the continuum and emission line data from the Interface Region Imaging Spectrograph (IRIS) of the brightest flaring magnetic footpoints in this flare. We compare the NUV spectra of the brightest pixels to new radiative-hydrodynamic predictions calculated with the RADYN code using constraints on a nonthermal electron beam inferred from the collisional thick-target modeling of hard X-ray data from Reuven Ramaty High Energy Solar Spectroscopic Imager . We show that the atmospheric response to a high beam flux density satisfactorily achieves the observed continuum brightness in the NUV. The NUV continuum emission in this flare is consistent with hydrogen (Balmer) recombination radiation that originates from low optical depth in a dense chromospheric condensation and from the stationary beam-heated layers just below the condensation. A model producing two flaring regions (a condensation and stationary layers) in the lower atmosphere is also consistent with the asymmetric Fe ii chromospheric emission line profiles observed in the impulsive phase.

79 ASTRONOMY AND ASTROPHYSICS↗

The Atmospheric Response to High Nonthermal Electron Beam Fluxes in Solar Flares. I. Modeling the Brightest NUV Footpoints in the X1 Solar Flare of 2014 March 29

The 2014 March 29 X1 solar flare (SOL20140329T17:48) produced bright continuum emission in the far- and near-ultraviolet (NUV) and highly asymmetric chromospheric emission lines, providing long-sought constraints on the heating mechanisms of the lower atmosphere in solar flares. We analyze the continuum and emission line data from the Interface Region Imaging Spectrograph (IRIS) of the brightest flaring magnetic footpoints in this flare. We compare the NUV spectra of the brightest pixels to new radiative-hydrodynamic predictions calculated with the RADYN code using constraints on a nonthermal electron beam inferred from the collisional thick-target modeling of hard X-ray data from Reuven Ramaty High Energy Solar Spectroscopic Imager. We show that the atmospheric response to a high beam flux density satisfactorily achieves the observed continuum brightness in the NUV. The NUV continuum emission in this flare is consistent with hydrogen (Balmer) recombination radiation that originates from low optical depth in a dense chromospheric condensation and from the stationary beam-heated layers just below the condensation. A model producing two flaring regions (a condensation and stationary layers) in the lower atmosphere is also consistent with the asymmetric Fe II chromospheric emission line profiles observed in the impulsive phase.

Kowalski, Adam F.↗

ON THE MISALIGNMENT BETWEEN CHROMOSPHERIC FEATURES AND THE MAGNETIC FIELD ON THE SUN

Observations of the upper chromosphere show an enormous amount of intricate fine structure. Much of this comes in the form of linear features, which are most often assumed to be well aligned with the direction of the magnetic field in the low plasma β regime that is thought to dominate the upper chromosphere. We use advanced radiative magnetohydrodynamic simulations, including the effects of ion-neutral interactions (using the generalized Ohm’s law) in the partially ionized chromosphere, to show that the magnetic field is often not well aligned with chromospheric features. This occurs where the ambipolar diffusion is large, i.e., ions and neutral populations decouple as the ion-neutral collision frequency drops, allowing the field to slip through the neutral population; where currents perpendicular to the field are strong; and where thermodynamic timescales are longer than or similar to those of ambipolar diffusion. We find this often happens in dynamic spicule or fibril-like features at the top of the chromosphere. This has important consequences for field extrapolation methods, which increasingly use such upper chromospheric features to help constrain the chromospheric magnetic field: our results invalidate the underlying assumption that these features are aligned with the field. In addition, our results cast doubt on results from 1D hydrodynamic models, which assume that plasma remains on the same field lines. Finally, our simulations show that ambipolar diffusion significantly alters the amount of free energy available in the coronal part of our simulated volume, which is likely to have consequences for studies of flare initiation.

79 ASTRONOMY AND ASTROPHYSICS↗

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 ν ≈ 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.

79 ASTRONOMY AND ASTROPHYSICS↗

EMERGENCE OF GRANULAR-SIZED MAGNETIC BUBBLES THROUGH THE SOLAR ATMOSPHERE. III. THE PATH TO THE TRANSITION REGION

We study, for the first time, the ascent of granular-sized magnetic bubbles from the solar photosphere through the chromosphere into the transition region and above. Such events occurred in a flux emerging region in NOAA 11850 on 2013 September 25. During that time, the first co-observing campaign between the Swedish 1-m Solar Telescope (SST) and the Interface Region Imaging Spectrograph (IRIS) spacecraft was carried out. Simultaneous observations of the chromospheric H α 656.28 nm and Ca ii 854.2 nm lines, plus the photospheric Fe i 630.25 nm line, were made with the CRISP spectropolarimeter at the Spitzer Space Telescope ( SST ) reaching a spatial resolution of 0.″14. At the same time, IRIS was performing a four-step dense raster of the emerging flux region, taking slit jaw images at 133 (C ii, transition region), 140 (Si iv, transition region), 279.6 (Mg ii k, core, upper chromosphere), and 283.2 nm (Mg ii 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 Si iv 140.3 and the Mg ii 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. 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.

79 ASTRONOMY AND ASTROPHYSICS↗

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-equilibrium hydrogen ionization by performing a 2D radiation-magnetohydrodynamics simulation featuring non-equilibrium ionization of both hydrogen and helium. The simulation includes the effect of hydrogen Lyα 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 local thermodynamic equilibrium (LTE) ionization shows that non-equilibrium helium ionization leads to higher temperatures in wavefronts and lower temperatures in the gas between shocks. Assuming LTE ionization results in a thermostat-like behavior 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-equilibrium 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.

79 ASTRONOMY AND ASTROPHYSICS↗

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 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 redshifts; 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.

79 ASTRONOMY AND ASTROPHYSICS↗

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 (HD) of the atmospheric plasma. We combine the Stanford unified Fokker–Planck code that models particle acceleration and transport with the RADYN HD 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 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 ≪10{sup 10} erg s{sup −1} cm{sup −2}, can cause “explosive” chromospheric evaporation and drive stronger up- and downflows (and HD 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 Å, He ii 304 Å, Ca ii K 3934 Å, and Si iv 1393 Å. One interesting result is the unusual high temperature (up to a few times 10{sup 5} K) of the formation site of He ii 304 Å, which is expected owing 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 nonthermal electrons and thus the poorly understood particle acceleration mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

THE FORMATION OF IRIS DIAGNOSTICS. VII. THE FORMATION OF THE O i 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 study here the formation and diagnostic potential of this line by means of non-local thermodynamic equilibrium modeling employing both 1D semi-empirical and 3D radiation magnetohydrodynamic models. We study the basic formation mechanisms and derive a quintessential model atom that incorporates 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 and k lines and the C ii lines around 133.5 nm.

74 ATOMIC AND MOLECULAR PHYSICS↗