Source author record

M. Carlsson

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

39 published item(s)

preprint2022arXiv

Chromospheric emission from nanoflare heating in RADYN simulations

Heating signatures from small-scale magnetic reconnection events in the solar atmosphere have proven to be difficult to detect through observations. Numerical models that reproduce flaring conditions are essential in the understanding of how nanoflares may act as a heating mechanism of the corona. We study the effects of non-thermal electrons in synthetic spectra from 1D hydrodynamic RADYN simulations of nanoflare heated loops to investigate the diagnostic potential of chromospheric emission from small-scale events. The Mg II h and k, Ca II H and K, Ca II 854.2 nm, H-alpha and H-beta chromospheric lines were synthesised from various RADYN models of coronal loops subject to electron beams of nanoflare energies. The contribution function to the line intensity was computed to better understand how the atmospheric response to the non-thermal electrons affects the formation of spectral lines and the detailed shape of their spectral profiles. The spectral line signatures arising from the electron beams highly depend on the density of the loop and the lower cutoff energy of the electrons. Low-energy (5 keV) electrons deposit their energy in the corona and transition region, producing strong plasma flows that cause both redshifts and blueshifts of the chromospheric spectra. Higher-energy (10 and 15 keV) electrons deposit their energy in the lower transition region and chromosphere, resulting in increased emission from local heating. Our results indicate that effects from small-scale events can be observed with ground-based telescopes, expanding the list of possible diagnostics for the presence and properties of nanoflares.

preprint2022arXiv

Quiet Sun Center to Limb Variation of the Linear Polarization Observed by CLASP2 Across the Mg II h & k Lines

The CLASP2 (Chromospheric LAyer SpectroPolarimeter 2) sounding rocket mission was launched on 2019 April 11. CLASP2 measured the four Stokes parameters of the Mg II h & k spectral region around 2800 Angstroms along a 200 arcsecond slit at three locations on the solar disk, achieving the first spatially and spectrally resolved observations of the solar polarization in this near ultraviolet region. The focus of the work presented here is the center-to-limb variation of the linear polarization across these resonance lines, which is produced by the scattering of anisotropic radiation in the solar atmosphere. The linear polarization signals of the Mg II h & k lines are sensitive to the magnetic field from the low to the upper chromosphere through the Hanle and magneto-optical effects. We compare the observations to theoretical predictions from radiative transfer calculations in unmagnetized semi-empirical models, arguing that magnetic fields and horizontal inhomogeneities are needed to explain the observed polarization signals and spatial variations. This comparison is an important step in both validating and refining our understanding of the physical origin of these polarization signatures, and also in paving the way toward future space telescopes for probing the magnetic fields of the solar upper atmosphere via ultraviolet spectropolarimetry.

preprint2021arXiv

A new view of the solar interface region from the Interface Region Imaging Spectrograph (IRIS)

The Interface Region Imaging Spectrograph (IRIS) has been obtaining near- and far-ultraviolet images and spectra of the solar atmosphere since July 2013. The unique combination of near and far-ultraviolet spectra and images at subarcsecond resolution and high cadence allows the tracing of mass and energy through the critical interface between the solar surface and the corona or solar wind. IRIS has enabled research into the fundamental physical processes thought to play a role in the low solar atmosphere such as ion-neutral interactions, magnetic reconnection, the generation, propagation, and dissipation of various types of waves, the acceleration of non-thermal particles, and various small-scale instabilities. These new findings have helped provide novel insights into a wide range of phenomena including the discovery of non-thermal particles in coronal nanoflares, the formation and impact of spicules and other jets, resonant absorption and dissipation of Alfvenic waves, energy release associated with braiding of magnetic field lines, the thermal instability in the chromosphere-corona mass and energy cycle, the contribution of waves, turbulence, and non-thermal particles in the energy deposition during flares and smaller-scale events such as UV bursts, and the role of flux ropes and other mechanisms in triggering CMEs. IRIS observations have also advanced studies of the connections between solar and stellar physics. Advances in numerical modeling, inversion codes, and machine learning techniques have played a key role in driving these new insights. With the advent of exciting new instrumentation both on the ground (e.g., DKIST, ALMA) and space-based (e.g., Parker Solar Probe, Solar Orbiter), we aim to review new insights based on IRIS observations or related modeling, and highlight some of the outstanding challenges that have been brought to the fore.

preprint2021arXiv

The formation and heating of chromospheric fibrils in a radiation-MHD simulation

Aims: We examine the movements of mass elements within dense fibrils using passive tracer particles (corks) in order to understand fibril creation and destruction processes. Methods: Simulated fibrils were selected at times when they were visible in an H$α$ image proxy. The corks were selected within fibril H$α$ formation regions. From this set, a cork was selected, and the field line passing through it was constructed. Other fibrilar corks close to this fieldline were also selected. Pathlines were constructed, revealing the locations of the mass elements forward and backward in time. The forces acting on these mass elements were analysed. Results: The main process of fibrilar loading in the simulation is different to the mass loading scenario in which waves steepen into shocks and push material upwards along the fieldlines from near their footpoints. Twisted low lying fieldlines were destabilised and then they untwisted, lifting the material trapped above their apexes via the Lorentz force. Subsequently, the majority of the mass drained down the fieldlines towards one or both footpoints under gravity. Material with large horizontal velocities could also elevated in rising fieldlines, creating somewhat parabolic motions, but material was not generally moving upward along a stationary magnetic fieldline during loading. Conclusions: The processes observed in the simulation are plausible additional scenarios. Criteria for observing such events are described. It is desirable that our simulations can also form more densely-packed fibrils from material fed from the base of field footpoints. Experimental parameters required to achieve this are discussed.

preprint2020arXiv

Acoustic-gravity wave propagation characteristics in 3D radiation hydrodynamic simulations of the solar atmosphere

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

preprint2020arXiv

Ambipolar diffusion in the Bifrost code

Ambipolar diffusion is a physical mechanism related to the drift between charged and neutral particles in a partially ionized plasma that is key in many different astrophysical systems. However, understanding its effects is challenging due to basic uncertainties concerning relevant microphysical aspects and the strong constraints it imposes on the numerical modeling. Our aim is to introduce a numerical tool that allows us to address complex problems involving ambipolar diffusion in which, additionally, departures from ionization equilibrium are important or high resolution is needed. The primary application of this tool is for solar atmosphere calculations, but the methods and results presented here may also have a potential impact on other astrophysical systems. We have developed a new module for the stellar atmosphere Bifrost code that improves its computational capabilities of the ambipolar diffusion term in the Generalized Ohm's Law. This module includes, among other things, collision terms adequate to processes in the coolest regions in the solar chromosphere. As a key feature of the module, we have implemented the Super Time-Stepping (STS) technique, that allows an important acceleration of the calculations. We have also introduced hyperdiffusion terms to guarantee the stability of the code. We show that to have an accurate value for the ambipolar diffusion coefficient in the solar atmosphere it is necessary to include as atomic elements in the equation of state not only hydrogen and helium but also the main electron donors like sodium, silicon and potassium. In addition, we establish a range of criteria to set up an automatic selection of the free parameters of the STS method that guarantees the best performance, optimizing the stability and speed for the ambipolar diffusion calculations. We validate the STS implementation by comparison with a self-similar analytical solution.

preprint2020arXiv

Capabilities of bisector analysis of the Si I 10827 A line for estimating line-of-sight velocities in the quiet Sun

We examine the capabilities of a fast and simple method to infer line-of-sight (LOS) velocities from observations of the photospheric Si I 10827 A line. This spectral line is routinely observed together with the chromospheric He I 10830 A triplet as it helps to constrain the atmospheric parameters. We study the accuracy of bisector analysis and a line core fit of Si I 10827 A. We employ synthetic profiles starting from the Bifrost enhanced network simulation. The profiles are computed solving the radiative transfer equation, including non-local thermodynamic equilibrium effects on the determination of the atomic level populations of Si I. We found a good correlation between the inferred velocities from bisectors taken at different line profile intensities and the original simulation velocity at given optical depths. This good correlation means that we can associate bisectors taken at different line-profile percentages with atmospheric layers that linearly increase as we scan lower spectral line intensities. We also determined that a fit to the line-core intensity is robust and reliable, providing information about atmospheric layers that are above those accessible through bisectors. Therefore, by combining both methods on the Si I 10827 A line, we can seamlessly trace the quiet-Sun LOS velocity stratification from the deep photosphere to higher layers until around $\log τ= -3.5$ in a fast and straightforward way. This method is ideal for generating quick-look reference images for future missions like the Daniel K. Inoue Solar Telescope and the European Solar Telescope, for example.

preprint2020arXiv

High-resolution observations of the solar photosphere, chromosphere and transition region. A database of coordinated IRIS and SST observations

NASA's Interface Region Imaging Spectrograph (IRIS) provides high resolution observations of the solar atmosphere through UV spectroscopy and imaging. Since the launch of IRIS in June 2013, we have conducted systematic observation campaigns in coordination with the Swedish 1-m Solar Telescope (SST) on La Palma. The SST provides complementary high-resolution observations of the photosphere and chromosphere. The SST observations include spectro-polarimetric imaging in photospheric Fe I lines and spectrally-resolved imaging in the chromospheric Ca II 8542 A, H-alpha, and Ca II K lines. We present a database of co-aligned IRIS and SST datasets that is open for analysis to the scientific community. The database covers a variety of targets including active regions, sunspots, plage, quiet Sun, and coronal holes.

preprint2020arXiv

Observational constraints on the origin of the elements II. 3D non-LTE formation of Ba ii lines in the solar atmosphere

Context. The pursuit of more realistic spectroscopic modelling and consistent abundances has led us to begin a new series of papers designed to improve current solar and stellar abundances of various atomic species. To achieve this, we have began updating the three-dimensional (3D) non-local thermodynamic equilibrium (non-LTE) radiative transfer code, Multi3D, and the equivalent one-dimensional (1D) non-LTE radiative transfer code, MULTI. Aims. We examine our improvements to these codes by redetermining the solar barium abundance. Barium was chosen for this test as it is an important diagnostic element of the s-process in the context of galactic chemical evolution. New Ba II + H collisional data for excitation and charge exchange reactions computed from first principles had recently become available and were included in the model atom. The atom also includes the effects of isotopic line shifts and hyperfine splitting. Method. A grid of 1D LTE barium lines were constructed with MULTI and fit to the four Ba II lines available to us in the optical region of the solar spectrum. Abundance corrections were then determined in 1D non-LTE, 3D LTE, and 3D non-LTE. A new 3D non-LTE solar barium abundance was computed from these corrections. Results. We present for the first time the full 3D non-LTE barium abundance of $A({\rm Ba})=2.27\pm0.02\pm0.01$, which was derived from four individual fully consistent barium lines. Errors here represent the systematic and random errors, respectively.

preprint2019arXiv

Nonequilibrium ionization and ambipolar diffusion in solar magnetic flux emergence processes

Magnetic flux emergence has been shown to be a key mechanism for unleashing a wide variety of solar phenomena. However, there are still open questions concerning the rise of the magnetized plasma through the atmosphere, mainly in the chromosphere, where the plasma departs from local thermodynamic equilibrium (LTE) and is partially ionized. We aim to investigate the impact of the nonequilibrium (NEQ) ionization and recombination and molecule formation of hydrogen, as well as ambipolar diffusion, on the dynamics and thermodynamics of the flux emergence process. Using the Bifrost code, we performed 2.5D numerical experiments of magnetic flux emergence from the convection zone up to the corona. The experiments include the NEQ ionization and recombination of atomic hydrogen, the NEQ formation and dissociation of H2 molecules, and the ambipolar diffusion term of the Generalized Ohm's Law. Our experiments show that the LTE assumption substantially underestimates the ionization fraction in most of the emerged region, leading to an artificial increase in the ambipolar diffusion and, therefore, in the heating and temperatures as compared to those found when taking the NEQ effects on the hydrogen ion population into account. We see that LTE also overestimates the number density of H2 molecules within the emerged region, thus mistakenly magnifying the exothermic contribution of the H2 molecule formation to the thermal energy during the flux emergence process. We find that the ambipolar diffusion does not significantly affect the amount of total unsigned emerged magnetic flux, but it is important in the shocks that cross the emerged region, heating the plasma on characteristic times ranging from 0.1 to 100 s. We also briefly discuss the importance of including elements heavier than hydrogen in the equation of state so as not to overestimate the role of ambipolar diffusion in the atmosphere.

preprint2016arXiv

ALMA Observations of the Sun in Cycle 4 and Beyond

This document was created by the Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) in preparation of the first regular observations of the Sun with the Atacama Large Millimeter/submillimeter Array (ALMA), which are anticipated to start in ALMA Cycle 4 in October 2016. The science cases presented here demonstrate that a large number of scientifically highly interesting observations could be made already with the still limited solar observing modes foreseen for Cycle 4 and that ALMA has the potential to make important contributions to answering long-standing scientific questions in solar physics. With the proposal deadline for ALMA Cycle 4 in April 2016 and the Commissioning and Science Verification campaign in December 2015 in sight, several of the SSALMON Expert Teams composed strategic documents in which they outlined potential solar observations that could be feasible given the anticipated technical capabilities in Cycle 4. These documents have been combined and supplemented with an analysis, resulting in recommendations for solar observing with ALMA in Cycle 4. In addition, the detailed science cases also demonstrate the scientific priorities of the solar physics community and which capabilities are wanted for the next observing cycles. The work on this White Paper effort was coordinated in close cooperation with the two international solar ALMA development studies led by T. Bastian (NRAO, USA) and R. Brajsa, (ESO). This document will be further updated until the beginning of Cycle 4 in October 2016. In particular, we plan to adjust the technical capabilities of the solar observing modes once finally decided and to further demonstrate the feasibility and scientific potential of the included science cases by means of numerical simulations of the solar atmosphere and corresponding simulated ALMA observations.

preprint2015arXiv

Halpha line profile asymmetries and the chromospheric flare velocity field

The asymmetries observed in the line profiles of solar flares can provide important diagnostics of the properties and dynamics of the flaring atmosphere. In this paper the evolution of the Halpha and Ca II 8542 Å lines are studied using high spatial, temporal and spectral resolution ground-based observations of an M1.1 flare obtained with the Swedish 1-m Solar Telescope. The temporal evolution of the Halpha line profiles from the flare kernel shows excess emission in the red wing (red asymmetry) before flare maximum, and excess in the blue wing (blue asymmetry) after maximum. However, the Ca II 8542 Å line does not follow the same pattern, showing only a weak red asymmetry during the flare. RADYN simulations are used to synthesise spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. We show that the red asymmetry observed in Halpha is not necessarily associated with plasma downflows, and the blue asymmetry may not be related to plasma upflows. Indeed, we conclude that the steep velocity gradients in the flaring chromosphere modifies the wavelength of the central reversal in the Halpha line profile. The shift in the wavelength of maximum opacity to shorter and longer wavelengths generates the red and blue asymmetries, respectively.

preprint2015arXiv

Homologous Helical Jets: Observations by IRIS, SDO and Hinode and Magnetic Modeling with Data-Driven Simulations

We report on observations of recurrent jets by instruments onboard the Interface Region Imaging Spectrograph (IRIS), Solar Dynamics Observatory (SDO) and Hinode spacecrafts. Over a 4-hour period on July 21st 2013, recurrent coronal jets were observed to emanate from NOAA Active Region 11793. FUV spectra probing plasma at transition region temperatures show evidence of oppositely directed flows with components reaching Doppler velocities of +/- 100 km/s. Raster Doppler maps using a Si IV transition region line show all four jets to have helical motion of the same sense. Simultaneous observations of the region by SDO and Hinode show that the jets emanate from a source region comprising a pore embedded in the interior of a supergranule. The parasitic pore has opposite polarity flux compared to the surrounding network field. This leads to a spine-fan magnetic topology in the coronal field that is amenable to jet formation. Time-dependent data-driven simulations are used to investigate the underlying drivers for the jets. These numerical experiments show that the emergence of current-carrying magnetic field in the vicinity of the pore supplies the magnetic twist needed for recurrent helical jet formation.

preprint2014arXiv

An Interface Region Imaging Spectrograph first view on Solar Spicules

Solar spicules have eluded modelers and observers for decades. Since the discovery of the more energetic type II, spicules have become a heated topic but their contribution to the energy balance of the low solar atmosphere remains unknown. Here we give a first glimpse of what quiet Sun spicules look like when observed with NASA's recently launched Interface Region Imaging Spectrograph (IRIS). Using IRIS spectra and filtergrams that sample the chromosphere and transition region we compare the properties and evolution of spicules as observed in a coordinated campaign with Hinode and the Atmospheric Imaging Assembly. Our IRIS observations allow us to follow the thermal evolution of type II spicules and finally confirm that the fading of Ca II H spicules appears to be caused by rapid heating to higher temperatures. The IRIS spicules do not fade but continue evolving, reaching higher and falling back down after 500-800 s. Ca II H type II spicules are thus the initial stages of violent and hotter events that mostly remain invisible in Ca II H filtergrams. These events have very different properties from type I spicules, which show lower velocities and no fading from chromospheric passbands. The IRIS spectra of spicules show the same signature as their proposed disk counterparts, reinforcing earlier work. Spectroheliograms from spectral rasters also confirm that quiet Sun spicules originate in bushes from the magnetic network. Our results suggest that type II spicules are indeed the site of vigorous heating (to at least transition region temperatures) along extensive parts of the upward moving spicular plasma.

preprint2014arXiv

Detection of supersonic downflows and associated heating events in the transition region above sunspots

IRIS data allow us to study the solar transition region (TR) with an unprecedented spatial resolution of 0.33 arcsec. On 2013 August 30, we observed bursts of high Doppler shifts suggesting strong supersonic downflows of up to 200 km/s and weaker, slightly slower upflows in the spectral lines Mg II h and k, C II 1336 Å, Si IV 1394 Å, and 1403 Å, that are correlated with brightenings in the slitjaw images (SJIs). The bursty behavior lasts throughout the 2 hr observation, with average burst durations of about 20 s. The locations of these short-lived events appear to be the umbral and penumbral footpoints of EUV loops. Fast apparent downflows are observed along these loops in the SJIs and in AIA, suggesting that the loops are thermally unstable. We interpret the observations as cool material falling from coronal heights, and especially coronal rain produced along the thermally unstable loops, which leads to an increase of intensity at the loop footpoints, probably indicating an increase of density and temperature in the TR. The rain speeds are on the higher end of previously reported speeds for this phenomenon, and possibly higher than the free-fall velocity along the loops. On other observing days, similar bright dots are sometimes aligned into ribbons, resembling small flare ribbons. These observations provide a first insight into small-scale heating events in sunspots in the TR.

preprint2014arXiv

Heating signatures in the disk counterparts of solar spicules in IRIS observations

We use coordinated observations with the Interface Region Imaging Spectrograph (IRIS) and the Swedish 1-m Solar Telescope (SST) to identify the disk counterpart of type II spicules in upper-chromospheric and transition region (TR) diagnostics. These disk counterparts were earlier identified through short-lived asymmetries in chromospheric spectral lines: rapid blue- or red-shifted excursions (RBEs or RREs). We find clear signatures of RBEs and RREs in Mg II h & k, often with excursions of the central h3 and k3 absorption features in concert with asymmetries in co-temporal and co-spatial H-alpha spectral profiles. We find spectral signatures for RBEs and RREs in C II 1335 and 1336 A and Si IV 1394 and 1403 A spectral lines and interpret this as a sign that type II spicules are heated to at least TR temperatures, supporting other recent work. These C II and Si IV spectral signals are weaker for a smaller network region than for more extended network regions in our data. A number of bright features around extended network regions observed in IRIS slit-jaw imagery SJI 1330 and 1400, recently identified as network jets, can be clearly connected to H-alpha RBEs and/or RREs in our coordinated data. We speculate that at least part of the diffuse halo around network regions in the IRIS SJI 1330 and 1400 images can be attributed to type II spicules with insufficient opacity in the C II and Si IV lines to stand out as single features in these passbands.

preprint2014arXiv

Hot Explosions in the Cool Atmosphere of the Sun

The solar atmosphere was traditionally represented with a simple one-dimensional model. Over the past few decades, this paradigm shifted for the chromosphere and corona that constitute the outer atmosphere, which is now considered a dynamic structured envelope. Recent observations by IRIS (Interface Region Imaging Spectrograph) reveal that it is difficult to determine what is up and down even in the cool 6000-K photosphere just above the solar surface: this region hosts pockets of hot plasma transiently heated to almost 100,000 K. The energy to heat and accelerate the plasma requires a considerable fraction of the energy from flares, the largest solar disruptions. These IRIS observations not only confirm that the photosphere is more complex than conventionally thought, but also provide insight into the energy conversion in the process of magnetic reconnection.

preprint2014arXiv

On the prevalence of small-scale twist in the solar chromosphere and transition region

The solar chromosphere and transition region (TR) form an interface between the Sun's surface and its hot outer atmosphere. Here most of the non-thermal energy that powers the solar atmosphere is transformed into heat, although the detailed mechanism remains elusive. High-resolution (0.33-arcsec) observations with NASA's Interface Region Imaging Spectrograph (IRIS) reveal a chromosphere and TR that are replete with twist or torsional motions on sub-arcsecond scales, occurring in active regions, quiet Sun regions, and coronal holes alike. We coordinated observations with the Swedish 1-m Solar Telescope (SST) to quantify these twisting motions and their association with rapid heating to at least TR temperatures. This view of the interface region provides insight into what heats the low solar atmosphere.

preprint2014arXiv

Prevalence of Small-scale Jets from the Networks of the Solar Transition Region and Chromosphere

As the interface between the Sun's photosphere and corona, the chromosphere and transition region play a key role in the formation and acceleration of the solar wind. Observations from the Interface Region Imaging Spectrograph reveal the prevalence of intermittent small-scale jets with speeds of 80-250 km/s from the narrow bright network lanes of this interface region. These jets have lifetimes of 20-80 seconds and widths of 300 km or less. They originate from small-scale bright regions, often preceded by footpoint brightenings and accompanied by transverse waves with ~20 km/s amplitudes. Many jets reach temperatures of at least ~100000 K and constitute an important element of the transition region structures. They are likely an intermittent but persistent source of mass and energy for the solar wind.

preprint2014arXiv

The Interface Region Imaging Spectrograph (IRIS)

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

preprint2014arXiv

The Unresolved Fine Structure Resolved - IRIS observations of the Solar Transition Region

The heating of the outer solar atmospheric layers, i.e., the transition region and corona, to high temperatures is a long standing problem in solar (and stellar) physics. Solutions have been hampered by an incomplete understanding of the magnetically controlled structure of these regions. The high spatial and temporal resolution observations with the Interface Region Imaging Spectrograph (IRIS) at the solar limb reveal a plethora of short, low lying loops or loop segments at transition-region temperatures that vary rapidly, on the timescales of minutes. We argue that the existence of these loops solves a long standing observational mystery. At the same time, based on comparison with numerical models, this detection sheds light on a critical piece of the coronal heating puzzle.

preprint2013arXiv

Heating of the magnetic chromosphere: observational constraints from Ca II 8542 spectra

The heating of the Sun's chromosphere remains poorly understood. While progress has been made on understanding what drives the quiet Sun internetwork chromosphere, chromospheric heating in strong magnetic field regions continues to present a difficult challenge, mostly because of a lack of observational constraints. We use high-resolution spectropolarimetric data from the Swedish 1-m Solar Telescope to identify the location and spatio-temporal properties of heating in the magnetic chromosphere. In particular, we report the existence of raised-core spectral line profiles in the Ca II 8542 line. These profiles are characterized by the absence of an absorption line core, showing a quasi-flat profile between +/- 0.5 Å, and are abundant close to magnetic bright-points and plage. Comparison with 3D MHD simulations indicates that such profiles occur when the line-of-sight goes through an "elevated temperature canopy" associated with the expansion with height of the magnetic field of flux concentrations. This temperature canopy in the simulations is caused by ohmic dissipation where there are strong magnetic field gradients. The raised-core profiles are thus indicators of locations of increased chromospheric heating. We characterize the location and temporal and spatial properties of such profiles in our observations, thus providing much stricter constraints on theoretical models of chromospheric heating mechanisms than before.

preprint2013arXiv

The formation of IRIS diagnostics I. A quintessential model atom of Mg II and general formation properties of the Mg II h&k lines

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

preprint2013arXiv

The formation of IRIS diagnostics II. The formation of the Mg II h&k lines in the solar atmosphere

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

preprint2012arXiv

A Sharp Look at Coronal Rain with Hinode/SOT and SST/CRISP

The tropical wisdom that when it is hot and dense we can expect rain might also apply to the Sun. Indeed, observations and numerical simulations have shown that strong heating at footpoints of loops, as is the case for active regions, puts their coronae out of thermal equilibrium, which can lead to a phenomenon known as catastrophic cooling. Following local pressure loss in the corona, hot plasma locally condenses in these loops and dramatically cools down to chromospheric temperatures. These blobs become bright in H-alpha and Ca II H in time scales of minutes, and their dynamics seem to be subject more to internal pressure changes in the loop rather than to gravity. They thus become trackers of the magnetic field, which results in the spectacular coronal rain that is observed falling down coronal loops. In this work we report on high resolution observations of coronal rain with the Solar Optical Telescope (SOT) on Hinode and CRISP at the Swedish Solar Telescope (SST). A statistical study is performed in which properties such as velocities and accelerations of coronal rain are derived. We show how this phenomenon can constitute a diagnostic tool for the internal physical conditions inside loops. Furthermore, we analyze transverse oscillations of strand-like condensations composing coronal rain falling in a loop, and discuss the possible nature of the wave. This points to the important role that coronal rain can play in the fields of coronal heating and coronal seismology.

preprint2012arXiv

Chromospheric magnetic fields. Observations, simulations and their interpretation

The magnetic field of the quiet-Sun chromosphere remains a mystery for solar physicists. The reduced number of chromospheric lines are intrinsically hard to model and only a few of them are magnetically sensitive. In this work, we use a 3D numerical simulation of the outer layers of the solar atmosphere, to asses the reliability of non-LTE inversions, in this case applied to the Ca II 8542 Å line. We show that NLTE inversions provide realistic estimates of physical quantities from synthetic observations.

preprint2012arXiv

On the Dynamics of the Fermi-Bose Model

We consider the exponential matrix representing the dynamics of the Fermi-Bose model in an undepleted bosonic field approximation. A recent application of this model is molecular dimers dissociating into its atomic compounds. The problem is solved in D spatial dimensions by dividing the system matrix into blocks with generalizations of Hankel matrices, here refered to as D-block-Hankel matrices. The method is practically useful for treating large systems, i.e. dense computational grids or higher spatial dimensions, either on a single standard computer or a cluster. In particular the results can be used for studies of three-dimensional physical systems of arbitrary geometry. We illustrate the generality of our approach by giving numerical results for the dynamics of Glauber type atomic pair correlation functions for a non-isotropic three-dimensional harmonically trapped molecular Bose-Einstein condensate.

preprint2012arXiv

Ubiquitous torsional motions in type II spicules

Spicules are long, thin, highly dynamic features that jut out ubiquitously from the solar limb. They dominate the interface between the chromosphere and corona and may provide significant mass and energy to the corona. We use high-quality observations with the Swedish 1-m Solar Telescope to establish that so-called type II spicules are characterized by the simultaneous action of three different types of motion: 1. field-aligned flows of order 50-100 km/s, 2. swaying motions of order 15-20 km/s, and 3. torsional motions of order 25-30 km/s. The first two modes have been studied in detail before, but not the torsional motions. Our analysis of many near-limb and off-limb spectra and narrow-band images yields strong evidence that most, if not all, type-II spicules undergo large torsional modulation and that these motions, like spicule swaying, represent Alfvenic waves propagating outward at several hundred km/s. The combined action of the different motions explains the similar morphology of spicule bushes in the outer red and blue wings of chromospheric lines, and needs to be taken into account when interpreting Doppler motions to derive estimates for field-aligned flows in spicules and determining the Alfvenic wave energy in the solar atmosphere. Our results also suggest that large torsional motion is an ingredient in the production of type-II spicules and that spicules play an important role in the transport of helicity through the solar atmosphere.

preprint2011arXiv

An exact lower energy bound for the infinite square well potential

We give a lower bound for the energy of a quantum particle in the infinite square well. We show that the bound is exact and identify the well-known element that fulfils the equality. Our approach is not directly dependent on the Schroedinger equation and illustrates an example where the wavefunction is obtained directly by energy minimization. The derivation presented can serve as an example of a variational method in an undergraduate level university course in quantum mechanics.

preprint2011arXiv

Quiet-Sun imaging asymmetries in NaI D1 compared with other strong Fraunhofer lines

Imaging spectroscopy of the solar atmosphere using the NaI D1 line yields marked asymmetry between the blue and red line wings: sampling a quiet-Sun area in the blue wing displays reversed granulation, whereas sampling in the red wing displays normal granulation. The MgI b2 line of comparable strength does not show this asymmetry, nor does the stronger CaII 8542 line. We demonstrate the phenomenon with near-simultaneous spectral images in NaI D1, MgI b2, and CaII 8542 from the Swedish 1-m Solar Telescope. We then explain it with line-formation insights from classical 1D modeling and with a 3D magnetohydrodynamical simulation combined with NLTE spectral line synthesis that permits detailed comparison with the observations in a common format. The cause of the imaging asymmetry is the combination of correlations between intensity and Dopplershift modulation in granular overshoot and the sensitivity to these of the steep profile flanks of the NaI D1 line. The MgI b2 line has similar core formation but much wider wings due to larger opacity buildup and damping in the photosphere. Both lines obtain marked core asymmetry from photospheric shocks in or near strong magnetic concentrations, less from higher-up internetwork shocks that produce similar asymmetry in the spatially averaged CaII 8542 profile.

preprint2011arXiv

Wave Propagation and Jet Formation in the Chromosphere

We present the results of numerical simulations of wave propagation and jet formation in solar atmosphere models with different magnetic field configurations. The presence in the chromosphere of waves with periods longer than the acoustic cutoff period has been ascribed to either strong inclined magnetic fields, or changes in the radiative relaxation time. Our simulations include a sophisticated treatment of radiative losses, as well as fields with different strengths and inclinations. Using Fourier and wavelet analysis techniques, we investigate the periodicity of the waves that travel through the chromosphere. We find that the velocity signal is dominated by waves with periods around 5 minutes in regions of strong, inclined field, including at the edges of strong flux tubes where the field expands, whereas 3-minute waves dominate in regions of weak or vertically oriented fields. Our results show that the field inclination is very important for long-period wave propagation, whereas variations in the radiative relaxation time have little effect. Furthermore, we find that atmospheric conditions can vary significantly on timescales of a few minutes, meaning that a Fourier analysis of wave propagation can be misleading. Wavelet techniques take variations with time into account and are more suitable analysis tools. Finally, we investigate the properties of jets formed by the propagating waves once they reach the transition region, and find systematic differences between the jets in inclined field regions and those in vertical field regions, in agreement with observations of dynamic fibrils.

preprint2010arXiv

High frequency waves in the solar atmosphere?

The present study addresses the following questions: How representative of the actual velocities in the solar atmosphere are the Doppler shifts of spectral lines? How reliable is the velocity signal derived from narrowband filtergrams? How well defined is the height of the measured Doppler signal? Why do phase difference spectra always pull to 0 degrees phase lag at high frequencies? Can we actually observe high frequency waves (P<= 70s)? What is the atmospheric MTF of high frequency waves? How reliably can we determine the energy flux of high frequency waves? We address these questions by comparing observations obtained with Hinode/NFI with results from two 3D numerical simulations (Oslo Stagger and CO5BOLD). Our results suggest that the observed high frequency Doppler velocity signal is caused by rapid height variations of the velocity response function in an atmosphere with strong velocity gradients and cannot be interpreted as evidence of propagating high frequency acoustic waves. Estimates of the energy flux of high frequency waves should be treated with caution, in particular those that apply atmospheric MTF corrections.

preprint2009arXiv

Neutral oxygen spectral line formation revisited with new collisional data: large departures from LTE at low metallicity

We present a detailed study of the non-Local Thermodynamic Equilibrium (non-LTE) formation of the high-excitation neutral oxygen 777 nm triplet in MARCS model atmospheres representative of late-type stars with spectral types F to K. We carried out the calculations using the statistical equilibrium code MULTI, including estimates of the impact on elemental abundance analysis. The atomic model employed includes, in particular, recent quantum-mechanical electron collision data. We confirm that the O I triplet lines form under non-LTE conditions in late-type stars, suffering negative abundance corrections with respect to LTE. At low metallicity, large line opacity stems from triplet-quintet intersystem electron collisions, a form of coupling previously not considered or seriously underestimated. The non-LTE effects become generally severe for models (both giants and dwarfs) with higher T_eff. Interestingly, in metal-poor turn-off stars, the negative non-LTE abundance corrections tend to rapidly become more severe towards lower metallicity. When neglecting H collisions, they amount to as much as ~ 0.9 dex and ~ 1.2 dex, respectively at [Fe/H]=-3 and [Fe/H]=-3.5. Even when such collisions are included, the LTE abundance remains a serious overestimate, correspondingly by ~ 0.5 dex and ~ 0.9 dex at such low metallicities. Although the poorly known inelastic hydrogen collisions thus remain an important uncertainty, the large metallicity-dependent non-LTE effects seem to point to a resulting "low" (compared to LTE) [O/Fe] in metal-poor halo stars.

preprint2009arXiv

On-disk counterparts of type II spicules in the Ca II 854.2 nm and Halpha lines

Recently a second type of spicules was discovered at the solar limb with the Solar Optical Telescope (SOT) onboard the Japanese Hinode spacecraft. These previously unrecognized type II spicules are thin chromospheric jets that are shorter-lived (10-60 s) and that show much higher apparent upward velocities (of order 50-100 km/s) than the classical spicules. Since they have been implicated in providing hot plasma to coronal loops, their formation, evolution and properties are important ingredients for a better understanding of the mass and energy balance of the low solar atmosphere. Here we report on the discovery of the disk counterparts of type II spicules using spectral imaging data in the Ca II 854.2 nm and Halpha lines with the CRisp Imaging SpectroPolarimeter (CRISP) at the Swedish Solar Telescope (SST) in La Palma. We find rapid blueward excursions in the line profiles of both chromospheric lines that correspond to thin, jet-like features that show apparent velocities of order 50 km/s. These blueward excursions seem to form a separate absorbing component with Doppler shifts of order 20 and 50 km/s for the Ca II 854.2 nm and Halpha line respectively. We show that the appearance, lifetimes, longitudinal and transverse velocities and occurrence rate of these rapid blue excursions on the disk are very similar to those of the type II spicules at the limb. A detailed study of the spectral line profiles in these events suggests that plasma is accelerated along the jet, and plasma is being heated throughout the short lifetime of the event.

preprint2009arXiv

Spicule-like structures observed in 3D realistic MHD simulations

We analyze features that resemble type i spicules in two different 3D numerical simulations in which we include horizontal magnetic flux emergence in a computational domain spanning the upper layers of the convection zone to the lower corona. The two simulations differ mainly in the preexisting ambient magnetic field strength and in the properties of the inserted flux tube. We use the Oslo Staggered Code (OSC) to solve the full MHD equations with non-grey and non-LTE radiative transfer and thermal conduction along the magnetic field lines. We find a multitude of features that show a spatiotemporal evolution that is similar to that observed in type i spicules, which are characterized by parabolic height vs. time profiles, and are dominated by rapid upward motion at speeds of 10-30 km/s, followed by downward motion at similar velocities. We measured the parameters of the parabolic profile of the spicules and find similar correlations between the parameters as those found in observations. The values for height (or length) and duration of the spicules found in the simulations are more limited in range than those in the observations. The spicules found in the simulation with higher preexisting ambient field have shorter length and smaller velocities. From the simulations, it appears that these kinds of spicules can, in principle, be driven by a variety of mechanisms that include p-modes, collapsing granules, magnetic energy release in the photosphere and lower chromosphere and convective buffeting of flux concentrations.

preprint2009arXiv

The quiet solar atmosphere observed and simulated in NaD1

The NaD1 line in the solar spectrum is sometimes attributed to the solar chromosphere. We study its formation in quiet-Sun network and internetwork. We first present high-resolution profile-resolved images taken in this line with the imaging spectrometer IBIS at the Dunn Solar Telescope and compare these to simultaneous chromospheric images taken in Ca 8542 and Halpha. We then model NaD1 formation by performing 3D NLTE profile synthesis for a snapshot from a 3D radiation-magnetohydrodynamics simulation. We find that most NaD1 brightness is not chromospheric but samples the magnetic concentrations that make up the quiet-Sun network in the photosphere, well below the height where they merge into chromospheric canopies, with aureoles from 3D resonance scattering. The line core is sensitive to magneto-acoustic shocks in and near magnetic concentrations, where shocks occur deeper than elsewhere, and may provide evidence of heating sited deep within magnetic concentrations.

preprint2009arXiv

Three-dimensional non-LTE radiative transfer computation of the Ca 8542 infrared line from a radiation-MHD simulation

Interpretation of imagery of the solar chromosphere in the widely used \CaIIIR infrared line is hampered by its complex, three-dimensional and non-LTE formation. Forward modelling is required to aid understanding. We use a 3D non-LTE radiative transfer code to compute synthetic \CaIIIR images from a radiation-MHD simulation of the solar atmosphere spanning from the convection zone to the corona. We compare the simulation with observations obtained with the CRISP filter at the Swedish 1--m Solar Telescope. We find that the simulation reproduces dark patches in the blue line wing caused by Doppler shifts, brightenings in the line core caused by upward-propagating shocks and thin dark elongated structures in the line core that form the interface between upward and downward gas motion in the chromosphere. The synthetic line core is narrower than the observed one, indicating that the sun exhibits both more vigorous large-scale dynamics as well as small scale motions that are not resolved within the simulation, presumably owing to a lack of spatial resolution.

preprint2003arXiv

Multi-level 3D non-LTE computations of lithium lines in the metal-poor halo stars HD140283 and HD84937

The lithium abundances in metal-poor halo stars are of importance for cosmology, galaxy evolution and stellar structure. In an attempt to study possible systematic errors in the derived Li abundances, the line formation of LiI lines has been investigated by means of realistic 3D hydrodynamical model atmospheres of halo stars and 3D non-LTE radiative transfer calculations. These are the first detailed 3D non-LTE computations reported employing a multi-level atomic model showing that such problems are now computationally tractable. The detailed computations reveal that the LiI population has a strong influence from the radiation field rather than the local gas temperature, indicating that the low derived Li abundances found by Asplund et al. (1999) are an artifact of their assumption of LTE. Relative to 3D LTE, the detailed calculations show pronounced over-ionization. In terms of abundances the 3D non-LTE values are within 0.05 dex of the 1D non-LTE results for the particular cases of HD140283 and HD84937, which is a consequence of the dominance of the radiation in determining the population density of LiI. Although 3D non-LTE can be expected to give results rather close (+/-0.1 dex) to 1D non-LTE for this reason, there may be systematic trends with metallicity and effective temperature.