explanation      blue bibcodes open ADS page with paths to full text
Author name code: martinez-sykora
ADS astronomy entries on 2022-09-14
author:"Martinez-Sykora, Juan"

---------------------------------------------------------
Title: A novel inversion method to determine the coronal magnetic
    field including the impact of bound-free absorption
Authors: Martinez-Sykora, Juan; Hansteen, Viggo H.; De Pontieu, Bart;
   Landi, Enrico
2022arXiv220813984M    Altcode:
  The magnetic field governs the corona; hence it is a crucial parameter
  to measure. Unfortunately, existing techniques for estimating its
  strength are limited by strong assumptions and limitations. These
  techniques include photospheric or chromospheric field extrapolation
  using potential or non-linear-force-free methods, estimates based on
  coronal seismology, or by direct observations via, e.g., the Cryo-NIRSP
  instrument on DKIST which will measure the coronal magnetic field,
  but only off the limb. Alternately, in this work we investigate a
  recently developed approach based on the magnetic-field-induced (MIT)
  transition of the \fex~257.261~Å. In order to examine this approach,
  we have synthesized several \fex\ lines from two 3D magnetohydrodynamic
  simulations, one modeling an emerging flux region and the second an
  established mature active region. In addition, we take bound-free
  absorption from neutral hydrogen and helium and singly ionised
  helium into account. The absorption from cool plasma that occurs at
  coronal heights has a significant impact on determining the magnetic
  field. We investigate in detail the challenges of using these \fex\
  lines to measure the field, considering their density and temperature
  dependence. We present a novel approach to deriving the magnetic field
  from the MIT using inversions of the differential emission measure as a
  function of the temperature, density, and magnetic field. This approach
  successfully estimates the magnetic field strength (up to \%18 relative
  error) in regions that do not suffer from significant absorption and
  that have relatively strong coronal magnetic fields ($>250$~G). This
  method allows the masking of regions where absorption is significant.

---------------------------------------------------------
Title: Detailed Description of the Collision Frequency in the Solar
    Atmosphere
Authors: Wargnier, Q. M.; Martínez-Sykora, J.; Hansteen, V. H.;
   De Pontieu, B.
2022ApJ...933..205W    Altcode:
  This work aims to provide an accurate description and calculations of
  collision frequencies in conditions relevant to the solar atmosphere. To
  do so, we focus on the detailed description of the collision
  frequency in the solar atmosphere based on a classical formalism with
  Chapman-Cowling collision integrals, as described by Zhdanov. These
  collision integrals allow linking the macroscopic transport fluxes
  of multifluid models to the kinetic scales involved in the Boltzmann
  equations. In this context, the collision frequencies are computed
  accurately while being consistent at the kinetic level. We calculate
  the collision frequencies based on this formalism and compare them with
  approaches commonly used in the literature for conditions typical of the
  solar atmosphere. To calculate the collision frequencies, we focus on
  the collision integral data provided by Bruno et al., which is based on
  a multicomponent hydrogen-helium mixture used for conditions typical for
  the atmosphere of Jupiter. We perform a comparison with the classical
  formalism of Vranjes & Krstic and Leake & Linton. We highlight
  the differences obtained in the distribution of the cross sections as
  functions of the temperature. Then, we quantify the disparities obtained
  in numerical simulations of a 2.5D solar atmosphere by calculating
  collision frequencies and ambipolar diffusion. This strategy allows
  us to validate and assess the accuracy of these collision frequencies
  for conditions typical of the solar atmosphere.

---------------------------------------------------------
Title: Modeling of small-scale phenomena
Authors: Martinez Sykora, Juan
2022cosp...44.2550M    Altcode:
  Over the last years, to mention a few small-scale phenomena, substantial
  progress has been made on modeling spicules, jets, surges, coronal rain,
  and the local dynamo. These advances are because models have reached
  smaller scales, but essential and complex physics is a key aspect of
  advancing our understanding of small-scale phenomena. To cite some of
  the relevant improvements: combination and interaction of even larger
  and smaller dynamic structures, ion-neutral interaction effects (e.g.,
  ambipolar diffusion and Hall term), or the ionization is treated in
  NEQ. In this review, I will present some of these advances in the
  numerical modeling of small-scale phenomena.

---------------------------------------------------------
Title: Multilfuid Alfven wave simulations to understand the chemical
    fractionation in the chromosphere and the role of the NEQ ionization.
Authors: Martinez Sykora, Juan
2022cosp...44.2575M    Altcode:
  We aim to study the role of Alfven waves on the first ionization
  potential (FIP) effects — the enrichment of low FIP elements in
  the outer solar atmosphere; this phenomenon is intimately tied to
  the physics of the chromosphere and the corona. For this study, we
  combine single fluid 2D radiative MHD models of the solar atmosphere
  using Bifrost, with a novel multi-fluids multi-species numerical
  code (Ebysus). With the former, we investigate the possible impact
  of non-equilibrium ionization within the region where the FIP may
  occur and its plasma properties. From the plasma properties from the
  Bifrost model, we initialize our multi-fluid models to investigate
  the fractionation and the role of the Ponderomotive force.

---------------------------------------------------------
Title: Probing the Physics of the Solar Atmosphere with the Multi-slit
    Solar Explorer (MUSE). II. Flares and Eruptions
Authors: Cheung, Mark C. M.; Martínez-Sykora, Juan; Testa, Paola;
   De Pontieu, Bart; Chintzoglou, Georgios; Rempel, Matthias; Polito,
   Vanessa; Kerr, Graham S.; Reeves, Katharine K.; Fletcher, Lyndsay; Jin,
   Meng; Nóbrega-Siverio, Daniel; Danilovic, Sanja; Antolin, Patrick;
   Allred, Joel; Hansteen, Viggo; Ugarte-Urra, Ignacio; DeLuca, Edward;
   Longcope, Dana; Takasao, Shinsuke; DeRosa, Marc L.; Boerner, Paul;
   Jaeggli, Sarah; Nitta, Nariaki V.; Daw, Adrian; Carlsson, Mats; Golub,
   Leon; The
2022ApJ...926...53C    Altcode: 2021arXiv210615591C
  Current state-of-the-art spectrographs cannot resolve the fundamental
  spatial (subarcseconds) and temporal (less than a few tens of
  seconds) scales 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 the Interface Region Imaging
  Spectrograph 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), subarcsecond-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 Extreme
  Ultraviolet Solar Telescope, 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., which focuses on investigating coronal heating with MUSE.

---------------------------------------------------------
Title: Probing the Physics of the Solar Atmosphere with the Multi-slit
    Solar Explorer (MUSE). I. Coronal Heating
Authors: De Pontieu, Bart; Testa, Paola; Martínez-Sykora, Juan;
   Antolin, Patrick; Karampelas, Konstantinos; Hansteen, Viggo; Rempel,
   Matthias; Cheung, Mark C. M.; Reale, Fabio; Danilovic, Sanja; Pagano,
   Paolo; Polito, Vanessa; De Moortel, Ineke; Nóbrega-Siverio, Daniel;
   Van Doorsselaere, Tom; Petralia, Antonino; Asgari-Targhi, Mahboubeh;
   Boerner, Paul; Carlsson, Mats; Chintzoglou, Georgios; Daw, Adrian;
   DeLuca, Edward; Golub, Leon; Matsumoto, Takuma; Ugarte-Urra, Ignacio;
   McIntosh, Scott W.; the MUSE Team
2022ApJ...926...52D    Altcode: 2021arXiv210615584D
  The Multi-slit Solar Explorer (MUSE) is a proposed mission composed of
  a multislit extreme ultraviolet (EUV) spectrograph (in three spectral
  bands around 171 Å, 284 Å, and 108 Å) and an EUV context imager (in
  two passbands around 195 Å and 304 Å). MUSE will provide unprecedented
  spectral and imaging diagnostics of the solar corona at high spatial
  (≤0.″5) and temporal resolution (down to ~0.5 s for sit-and-stare
  observations), thanks to its innovative multislit design. By obtaining
  spectra in four bright EUV lines (Fe IX 171 Å, Fe XV 284 Å, Fe XIX-Fe
  XXI 108 Å) covering a wide range of transition regions and coronal
  temperatures along 37 slits simultaneously, MUSE will, for the first
  time, "freeze" (at a cadence as short as 10 s) 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)
  to the large-scale (~170″ × 170″) atmospheric response. We use
  numerical modeling to showcase how MUSE will constrain the properties of
  the solar atmosphere on spatiotemporal scales (≤0.″5, ≤20 s) and
  the large field of view on which state-of-the-art models of the physical
  processes that drive coronal heating, flares, and coronal mass ejections
  (CMEs) make distinguishing and testable predictions. We describe the
  synergy between MUSE, the single-slit, high-resolution Solar-C EUVST
  spectrograph, and ground-based observatories (DKIST and others), and
  the critical role MUSE plays because of the multiscale nature of the
  physical processes involved. In this first paper, we focus on coronal
  heating mechanisms. An accompanying paper focuses on flares and CMEs.

---------------------------------------------------------
Title: Thermal Instability-Induced Fundamental Magnetic Field Strands
    in the Solar Corona
Authors: Antolin, Patrick; Martínez-Sykora, Juan; Şahin, Seray
2022ApJ...926L..29A    Altcode:
  Thermal instability is a fundamental process of astrophysical
  plasmas. It is expected to occur whenever the cooling is dominated
  by radiation and cannot be compensated for by heating. In this work,
  we conduct 2.5D radiation MHD simulations with the Bifrost code
  of an enhanced activity network in the solar atmosphere. Coronal
  loops are produced self-consistently, mainly through Joule heating,
  which is sufficiently stratified and symmetric to produce thermal
  nonequilibrium. During the cooling and driven by thermal instability,
  coronal rain is produced along the loops. Due to flux freezing,
  the catastrophic cooling process leading to a rain clump produces a
  local enhancement of the magnetic field, thereby generating a distinct
  magnetic strand within the loop up to a few Gauss stronger than the
  surrounding coronal field. These strands, which can be considered
  fundamental, are a few hundred kilometers in width, span most of
  the loop leg, and emit strongly in the UV and extreme UV (EUV),
  thereby establishing a link between the commonly seen rain strands
  in the visible spectrum with the observed EUV coronal strands at
  high resolution. The compression downstream leads to an increase in
  temperature that generates a plume-like structure, a strongly emitting
  spicule-like feature, and short-lived brightening in the UV during
  the rain impact, providing an explanation for similar phenomena seen
  with IRIS. Thermal instability and nonequilibrium can therefore be
  associated with localized and intermittent UV brightening in the
  transition region and chromosphere, as well as contribute to the
  characteristic filamentary morphology of the solar corona in the EUV.

---------------------------------------------------------
Title: Probing the Physics of the Solar Atmosphere with the Multi-slit
Solar Explorer (MUSE): II. Flares and Eruptions
Authors: Cheung, Chun Ming Mark; Martinez-Sykora, Juan; Testa, Paola;
   De Pontieu, Bart; Chintzoglou, Georgios; Rempel, Matthias; Polito,
   Vanessa; Kerr, Graham; Reeves, Katharine; Fletcher, Lyndsay; Jin,
   Meng; Nobrega, Daniel; Danilovic, Sanja; Antolin, Patrick; Allred,
   Joel; Hansteen, Viggo; Ugarte-Urra, Ignacio; DeLuca, Edward; Longcope,
   Dana; Takasao, Shinsuke; DeRosa, Marc; Boerner, Paul; Jaeggli, Sarah;
   Nitta, Nariaki; Daw, Adrian; Carlsson, Mats; Golub, Leon
2021AGUFMSH51A..08C    Altcode:
  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, also submitted to SH-17),
  which focuses on investigating coronal heating with MUSE.

---------------------------------------------------------
Title: Description of collisional frequencies for multifluid MHD
    models with Chapman-Cowling collision integrals
Authors: Wargnier, Q.; Martinez-Sykora, Juan; Hansteen, Viggo;
   Magin, Thierry
2021AGUFMSH45B2362W    Altcode:
  We focus on the detailed description of the collisional frequency in the
  solar atmosphere based on a classical formalism with Chapman-Cowling
  collision integrals, as described by Zhdanov (2002) in the context of
  the 13N-moment model derived with a method of Grad (Grad 1949). These
  collision integrals allow linking the macroscopic transport fluxes
  of multifluid models to the kinetic scales involved in the Boltzmann
  equations. In this context, the collisional frequencies are computed
  accurately while being consistent at the kinetic level. We calculate
  the collisional frequencies based on this formalism and compare them
  with approaches commonly used in the literature in solar atmosphere
  conditions. To calculate the collisional frequencies, we focus on
  the collision integrals data provided by Bruno et al. (2010), which
  is based on a multicomponent hydrogen-helium mixture used in Jupiter
  atmosphere conditions. We propose a comparison with the classical
  formalism of Vranjes & Krstic (2013) and Leake & Linton
  (2013). We compare it with the formalism used in the three approaches
  and highlight the differences obtained in the distribution of the
  cross sections as functions of the temperature. Then, we quantify
  the disparities obtained in postprocessed simulations of a 2.5D solar
  atmosphere with the Bifrost code (see Gudiksen et al. 2011). Finally,
  we assess the impact of the collisional frequency in a simulated
  2.5D solar atmosphere with a single-fluid radiative MHD model with
  ambipolar diffusion to consider ion-neutral interactions. Significant
  disparities in the cross sections have been obtained between these
  three formalisms. or instance, we note that Vranjes & Krstic 2013
  did no integrate the transport cross sections. We will describe the
  impact of these discrepancies from previous results and the importance
  of doing these calculations properly.

---------------------------------------------------------
Title: Multi-fluid Simulations of Small-scale Collisional Plasma
    Instabilities in the Solar Chromosphere
Authors: Evans, Samuel; Oppenheim, Meers; Martinez-Sykora, Juan;
   Dimant, Yakov; Xiao, Richard
2021AGUFMSH25A2073E    Altcode:
  The chromosphere may be the most complex region in the solar
  atmosphere. The neutral flows, metal ions, magnetic field structure,
  radiation, and non-local thermal equilibrium effects may all play
  an important role in heating the solar atmosphere from a few thousand
  Kelvin to over a million Kelvin. The chromosphere also spans temperature
  ranges that cause it to transition between predominantly neutral to
  predominantly ionized, and the ions to transition from demagnetized
  to magnetized. The flows, densities, and temperatures of the various
  species in this region, along with the electric and magnetic fields,
  create conditions which can trigger the multi-species thermal plus
  Farley-Buneman instability. This instability causes the plasma
  to develop waves that lead to turbulence and heating, which may
  help to explain the discrepancy between models and observations of
  heating in the chromosphere. In this work, we present simulations
  of this instability, using the multi-fluid multi-species (MFMS)
  code, Ebysus. These simulations model a small piece of the coldest
  regions of the chromosphere with a realistic, but externally imposed
  current. We analyze the resulting heating, and compare the simulation
  with results from a particle-in-cell (PIC) code. The ability to simulate
  this instability in a multi-fluid code should enable simulations with
  chromospheric parameters unobtainable by a PIC code. We expect the
  result of this study will be to determine the effects of this type of
  small-scale turbulence on heating and transport in the larger scale
  solar atmosphere.

---------------------------------------------------------
Title: Probing the physics of coronal heating with the Multi-slit
    Solar Explorer (MUSE)
Authors: De Pontieu, Bart; Testa, Paola; Martinez-Sykora, Juan;
   Cheung, Chun Ming Mark
2021AGUFMSH55B1836D    Altcode:
  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 171Å, 284Å, and 108Å)
  and an EUV context imager (in two narrow passbands around 195Å and
  304Å). MUSE will provide unprecedented spectral and imaging diagnostics
  of the solar corona at high spatial (~0.5 arcseconds), and temporal
  resolution (down to ~0.5 seconds) thanks to its innovative multi-slit
  design. By obtaining spectra in 4 bright EUV lines (Fe IX 171Å,
  Fe XV 284Å, Fe XIX-XXI 108Å) 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. We focus on how comparisons between MUSE observations and
  theoretical models will significantly further our understanding of
  coronal heating mechanisms. This is a companion paper to Cheung et
  al. (2021), also submitted to SH-17.

---------------------------------------------------------
Title: Evidence of the multi-thermal nature of spicular
    downflows. Impact on solar atmospheric heating
Authors: Bose, Souvik; Rouppe van der Voort, Luc; Joshi, Jayant;
   Henriques, Vasco M. J.; Nóbrega-Siverio, Daniel; Martínez-Sykora,
   Juan; De Pontieu, Bart
2021A&A...654A..51B    Altcode: 2021arXiv210802153B
  Context. Spectroscopic observations of the emission lines formed in the
  solar transition region commonly show persistent downflows on the order
  of 10−15 km s<SUP>−1</SUP>. The cause of such downflows, however, is
  still not fully clear and has remained a matter of debate. <BR /> Aims:
  We aim to understand the cause of such downflows by studying the coronal
  and transition region responses to the recently reported chromospheric
  downflowing rapid redshifted excursions (RREs) and their impact on the
  heating of the solar atmosphere. <BR /> Methods: We have used two sets
  of coordinated data from the Swedish 1 m Solar Telescope, the Interface
  Region Imaging Spectrograph, and the Solar Dynamics Observatory for
  analyzing the response of the downflowing RREs in the transition
  region and corona. To provide theoretical support, we use an already
  existing 2.5D magnetohydrodynamic simulation of spicules performed
  with the Bifrost code. <BR /> Results: We find ample occurrences of
  downflowing RREs and show several examples of their spatio-temporal
  evolution, sampling multiple wavelength channels ranging from the cooler
  chromospheric to the hotter coronal channels. These downflowing features
  are thought to be likely associated with the returning components of
  the previously heated spicular plasma. Furthermore, the transition
  region Doppler shifts associated with them are close to the average
  redshifts observed in this region, which further implies that these
  flows could (partly) be responsible for the persistent downflows
  observed in the transition region. We also propose two mechanisms -
  (i) a typical upflow followed by a downflow and (ii) downflows along a
  loop -from the perspective of a numerical simulation that could explain
  the ubiquitous occurrence of such downflows. A detailed comparison
  between the synthetic and observed spectral characteristics reveals a
  distinctive match and further suggests an impact on the heating of the
  solar atmosphere. <BR /> Conclusions: We present evidence that suggests
  that at least some of the downflowing RREs are the chromospheric
  counterparts of the transition region and lower coronal downflows. <P
  />Movies associated to Figs. 1-3, 8, and 10 are available at <A
  href="https://www.aanda.org/10.1051/0004-6361/202141404/olm">https://www.aanda.org</A>

---------------------------------------------------------
Title: Probing Uncertainties in Diagnostics of a Synthetic
    Chromosphere
Authors: Schmit, Don; Martinez-Sykora, Juan; Pereira, Tiago; Asensio
   Ramos, Andrés
2021ApJ...913...71S    Altcode:
  Effective spectroscopic diagnostics rely on the ability to convert a
  particular flux measurement into a physical interpretation. Knowledge
  of uncertainty is a central component of diagnostics. We present data
  from a simulated solar-like chromosphere, where we have addressed the
  question of whether degeneracy is a problem in mapping from a non-LTE
  chromospheric line profile to a particular vertical stratification
  of atmospheric properties along the line of sight. Our results
  indicate that stratification degeneracies do exist, at least in our
  simulated atmosphere. We quantify this effect through the creation
  of posterior densities for atmospheric properties based on the Mg
  II h line profile using the approximate Bayesian computation (ABC)
  technique. We find that the predictive power of the ABC temperature
  posterior systematically varies as a function of atmospheric column
  mass and ground-truth temperature. The ABC posteriors more effectively
  reproduce the spectral intensity in the Ca II 8542 Å line than they do
  temperature stratification, although residual error in the Ca II line
  core is common. Our results illustrate that some degeneracies should
  be alleviated through simultaneous analysis of multiple chromospheric
  lines.

---------------------------------------------------------
Title: A New View of the Solar Interface Region from the Interface
    Region Imaging Spectrograph (IRIS)
Authors: De Pontieu, Bart; Polito, Vanessa; Hansteen, Viggo; Testa,
   Paola; Reeves, Katharine K.; Antolin, Patrick; Nóbrega-Siverio,
   Daniel Elias; Kowalski, Adam F.; Martinez-Sykora, Juan; Carlsson,
   Mats; McIntosh, Scott W.; Liu, Wei; Daw, Adrian; Kankelborg, Charles C.
2021SoPh..296...84D    Altcode: 2021arXiv210316109D
  The Interface Region Imaging Spectrograph (IRIS) has been obtaining
  near- and far-ultraviolet images and spectra of the solar atmosphere
  since July 2013. IRIS is the highest resolution observatory to provide
  seamless coverage of spectra and images from the photosphere into the
  low corona. The unique combination of near- and far-ultraviolet spectra
  and images at sub-arcsecond resolution and high cadence allows the
  tracing of mass and energy through the critical interface between the
  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 waves, the acceleration of
  non-thermal particles, and various small-scale instabilities. IRIS has
  provided insights into a wide range of phenomena including the discovery
  of non-thermal particles in coronal nano-flares, the formation and
  impact of spicules and other jets, resonant absorption and dissipation
  of Alfvénic waves, energy release and jet-like dynamics associated
  with braiding of magnetic-field lines, the role of turbulence and the
  tearing-mode instability in reconnection, 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 various other mechanisms in triggering and driving CMEs. IRIS
  observations have also been used to elucidate the physical mechanisms
  driving the solar irradiance that impacts Earth's upper atmosphere,
  and the connections between solar and stellar physics. Advances in
  numerical modeling, inversion codes, and machine-learning techniques
  have played a key role. With the advent of exciting new instrumentation
  both on the ground, e.g. the Daniel K. Inouye Solar Telescope (DKIST)
  and the Atacama Large Millimeter/submillimeter Array (ALMA), and
  space-based, e.g. the Parker Solar Probe and the Solar Orbiter, we aim
  to review new insights based on IRIS observations or related modeling,
  and highlight some of the outstanding challenges.

---------------------------------------------------------
Title: Thermal instability-induced fundamental magnetic strands in
    coronal loops
Authors: Antolin, Patrick; Martinez-Sykora, Juan
2021cosp...43E.968A    Altcode:
  Thermal instability is a fundamental process of astrophysical
  plasmas. It is expected to occur whenever the cooling is dominated
  by radiation and cannot be compensated by heating. This mechanism has
  been invoked to explain structures at multiple scales in the Universe,
  from the filamentary structure of the ISM to the phenomenon of coronal
  rain in the solar corona. In this work we conduct 2.5-D Radiation MHD
  simulations with the Bifrost code of an enhanced activity network in
  the solar atmosphere. Coronal loops are produced self-consistently,
  mainly through Ohmic heating, which is stratified and of a high enough
  frequency as to produce thermal non-equilibrium. During the cooling
  and driven by thermal instability, coronal rain is produced along the
  loops. Due to flux freezing, the catastrophic cooling process leading
  to a rain clump produces a local enhancement of the magnetic field,
  thereby generating a distinct magnetic strand within the loop up to a
  few Gauss stronger than the ambient corona. The compression downstream
  leads to an increase in temperature that generates a strongly emitting
  spicule-like feature in the UV during the rain impact. The stronger
  magnetic field strength in the rarefied upstream region has a stronger
  Ohmic heating, leading to a filamentary coronal strand with enhanced
  EUV emission. Thermal instability and _x0005_non-equilibrium can
  therefore be associated with localised and intermittent UV brightening
  in the transition region and chromosphere, as well as contribute to
  the characteristic filamentary morphology of the solar corona in the
  EUV. An additional effect of a strand with enhanced magnetic field is to
  serve as a waveguide, which combined with the Ohmic heating can act as a
  seed to sustain the coronal loop and the thermal non-equilibrium cycle.

---------------------------------------------------------
Title: ALMA and IRIS Observations of the Solar
    Chromosphere. II. Structure and Dynamics of Chromospheric Plages
Authors: Chintzoglou, Georgios; De Pontieu, Bart; Martínez-Sykora,
   Juan; Hansteen, Viggo; de la Cruz Rodríguez, Jaime; Szydlarski,
   Mikolaj; Jafarzadeh, Shahin; Wedemeyer, Sven; Bastian, Timothy S.;
   Sainz Dalda, Alberto
2021ApJ...906...83C    Altcode: 2020arXiv201205970C
  We propose and employ a novel empirical method for determining
  chromospheric plage regions, which seems to better isolate a plage from
  its surrounding regions than other methods commonly used. We caution
  that isolating a plage from its immediate surroundings must be done
  with care in order to successfully mitigate statistical biases that,
  for instance, can impact quantitative comparisons between different
  chromospheric observables. Using this methodology, our analysis suggests
  that λ = 1.25 mm free-free emission in plage regions observed with
  the Atacama Large Millimeter/submillimeter Array (ALMA)/Band6 may
  not form in the low chromosphere as previously thought, but rather
  in the upper chromospheric parts of dynamic plage features (such as
  spicules and other bright structures), i.e., near geometric heights
  of transition-region temperatures. We investigate the high degree of
  similarity between chromospheric plage features observed in ALMA/Band6
  (at 1.25 mm wavelengths) and the Interface Region Imaging Spectrograph
  (IRIS)/Si IV at 1393 Å. We also show that IRIS/Mg II h and k are
  not as well correlated with ALMA/Band6 as was previously thought,
  and we discuss discrepancies with previous works. Lastly, we report
  indications of chromospheric heating due to propagating shocks supported
  by the ALMA/Band6 observations.

---------------------------------------------------------
Title: ALMA and IRIS Observations of the Solar Chromosphere. I. An
    On-disk Type II Spicule
Authors: Chintzoglou, Georgios; De Pontieu, Bart; Martínez-Sykora,
   Juan; Hansteen, Viggo; de la Cruz Rodríguez, Jaime; Szydlarski,
   Mikolaj; Jafarzadeh, Shahin; Wedemeyer, Sven; Bastian, Timothy S.;
   Sainz Dalda, Alberto
2021ApJ...906...82C    Altcode: 2020arXiv200512717C
  We present observations of the solar chromosphere obtained
  simultaneously with the Atacama Large Millimeter/submillimeter Array
  (ALMA) and the Interface Region Imaging Spectrograph. The observatories
  targeted a chromospheric plage region of which the spatial distribution
  (split between strongly and weakly magnetized regions) allowed the
  study of linear-like structures in isolation, free of contamination
  from background emission. Using these observations in conjunction with
  a radiative magnetohydrodynamic 2.5D model covering the upper convection
  zone all the way to the corona that considers nonequilibrium ionization
  effects, we report the detection of an on-disk chromospheric spicule
  with ALMA and confirm its multithermal nature.

---------------------------------------------------------
Title: Expected spectropolarimetric observables in the lower solar
    atmosphere from 3D radiative MHD models
Authors: Sainz Dalda, A.; Gosic, M.; Martinez-Sykora, J.
2020AGUFMSH0010019S    Altcode:
  No abstract at ADS

---------------------------------------------------------
Title: ALMA and IRIS Observations Highlighting the Dynamics and
    Structure of Chromospheric Plage
Authors: Chintzoglou, G.; De Pontieu, B.; Martinez-Sykora, J.;
   Hansteen, V. H.; de la Cruz Rodriguez, J.; Szydlarski, M.; Jafarzadeh,
   S.; Wedemeyer, S.; Bastian, T.; Sainz Dalda, A.
2020AGUFMSH0010009C    Altcode:
  We present observations of the solar chromosphere obtained
  simultaneously with the Atacama Large Millimeter/submillimeter Array
  (ALMA) and the Interface Region Imaging Spectrograph (IRIS). The
  observatories targeted a chromospheric plage region of which the spatial
  distribution (split between strongly and weakly magnetized regions)
  allowed the study of linear-like structures in isolation, free of
  contamination from background emission. Using these observations
  in conjunction with a radiative magnetohydrodynamic 2.5D model
  covering the upper convection zone all the way to the corona
  that considers non-equilibrium ionization effects, we report the
  detection of an on-disk chromospheric spicule with ALMA and confirm
  its multithermal nature. In addition, we discuss the strikingly high
  degree of similarity between chromospheric plage features observed
  in ALMA/Band6 and IRIS/\ion{Si}{4} (also reproduced in our model)
  suggesting that ALMA/Band6 does not observe in the low chromosphere as
  previously thought but rather observes the upper chromospheric parts
  of structures such as spicules and other bright structures above plage
  at geometric heights near transition region temperatures. We also show
  that IRIS/\ion{Mg}{2} is not as well correlated with ALMA/Band6 as was
  previously thought. For these comparisons, we propose and employ a novel
  empirical method for the determination of plage regions, which seems
  to better isolate plage from its surrounding regions as compared to
  other methods commonly used. We caution that isolating plage from its
  immediate surroundings must be done with care to mitigate statistical
  bias in quantitative comparisons between different chromospheric
  observables. Lastly, we report indications for chromospheric heating
  due to traveling shocks supported by the ALMA/Band6 observations.

---------------------------------------------------------
Title: On the velocity drift between ions in the solar atmosphere
Authors: Martínez-Sykora, J.; Szydlarski, M.; Hansteen, V. H.;
   De Pontieu, B.
2020AGUFMSH0010017M    Altcode:
  Very recent results that compare ALMA and IRIS observations with 2D
  radiative MHD including non-equilibrium ionization and ambipolar
  diffusion models of the type II spicules reveal that these models
  may underestimate the energy dissipated in the chromosphere. The
  solar atmosphere is composed of many species that are populated at
  different ionization and excitation levels. The upper chromosphere,
  transition region, and corona are nearly collisionless. Consequently,
  slippage between, for instance, ions and neutral particles, or
  interactions between separate species, may play an important role
  in the local momentum and energy balance. The interaction between
  species is missing in the 2D radiative MHD model. We have developed a
  3D multi-fluid and multi-species numerical code (Ebysus) to investigate
  such effects. Ebysus is capable of treating species (e.g., hydrogen,
  helium, etc) and fluids (neutrals, excited and ionized elements)
  separately. Treating different species as different fluids leads
  to drifts between different ions and an electric field that couple
  these motions. Different ionized species and momentum exchange can
  dissipate this velocity drift, i.e., convert wave kinetic energy into
  thermal energy. High-frequency Alfven waves, driven for instance by
  reconnection, thought to occur in the solar atmosphere, can drive such
  multi-ion velocity drifts.

---------------------------------------------------------
Title: Multi-Fluid Simulations of Collisional Plasma Instabilities
    in the Solar Chromosphere
Authors: Evans, S.; Martínez-Sykora, J.; Dimant, Y. S.; Oppenheim,
   M. M.
2020AGUFMSH0010016E    Altcode:
  The state-of-the-art instruments (IRIS-ALMA-DKIST) are challenging our
  understanding of the thermodynamics of the solar chromosphere. The
  chromosphere is the region between the solar surface and the
  million-degree corona. Chromospheric physics involves complex processes
  that are not well-described by the standard MHD assumptions, such as
  partially ionized regions, diffusive terms becoming non-negligible
  compared to collisional effects, and ion species which transition
  between magnetized and demagnetized. All these ingredients are necessary
  to develop thermal and Farley-Buneman instabilities in the solar
  chromosphere. By using the new multi-fluid multi-species (MFMS) code
  Ebysus, we work to simulate the effects of thermal and Farley-Buneman
  instabilities in the solar chromosphere. Ebysus models each ion species
  and ionization level as a separate fluid and tracks many fluids at
  once, making it well-suited to investigate these instabilities with
  a fluid-model. This paper will describe these instabilities as they
  appear in the solar chromosphere using Ebysus and making comparisons to
  PIC code simulations. <P />This work is supported by NSF Grant 1903416.

---------------------------------------------------------
Title: A Multi-Fluid Multi-Species (MFMS) numerical code for
    simulating the solar atmosphere
Authors: Wargnier, Q.; Martínez-Sykora, J.; Hansteen, V. H.;
   Szydlarski, M.; Evans, S.
2020AGUFMSH0370007W    Altcode:
  The solar atmosphere is characterized by multiple periods and sizes
  involving a large spectrum of temporal and spatial scales. It is
  regulated through complex interactions between different species and
  chemical reactions, amongst other physical processes . Because of
  this complexity, an accurate description of all of these multi-scale
  phenomena in the solar atmosphere is out of the reach of standard
  single-fluid MHD models (Hartlep et al. 2012). Furthermore, the
  enrichment of low first ionization potential elements in the outer
  layer of the solar atmosphere (the FIP effect) is not fully described
  by the current theoretical models (e.g. Laming et al. 2017), since
  they employ semi-empirical static atmospheres.

---------------------------------------------------------
Title: Chromospheric and TR diagnostics in a large scale numerical
simulation of flux emergence: Synthetic vs Real observables
Authors: Hansteen, V. H.; De Pontieu, B.; Testa, P.; Gosic, M.;
   Martinez-Sykora, J.
2020AGUFMSH0010021H    Altcode:
  Field stored just below or rising to the photosphere will break through
  the surface and enter the upper atmosphere once the gradient of the
  subphotospheric field strength becomes sufficiently large. Opposite
  polarity flux bundles will reconnect in the photosphere and above,
  to form steadily longer loops that expand into the outer solar
  atmosphere, forming the corona. Some of the emerging flux is likely
  due to a local dynamo, but also the direct emergence of large scale
  magnetic structures from below is important, even in the quiet Sun. A
  significant proportion of this field likely reaches the chromosphere
  and may leave imprint on chromospheric dynamics and energetics. Using
  large scale numerical models (72x72x60) Mm and the high resolution
  spectra and slit jaw images from IRIS, as well as photospheric data
  from Hinode/SOT, and SDO/HMI we study the interactions between the
  magnetic flux caught in the granular flow field and the chromosphere
  and chromospheric field above. We will compare synthetic observables
  of the photospheric Fe I 617.3 nm line, the chromospheric Mg II h&amp;k
  lines, and the transition region Si IV lines, with their observational
  counterparts. We will also generate synthetic ALMA band 3 images. The
  comparison of synthetic observational data will let us draw conclusions
  as to the validity of the numerical modeling and the importance of flux
  emergence for the dynamics and energetics of the outer solar atmosphere.

---------------------------------------------------------
Title: On the Velocity Drift between Ions in the Solar Atmosphere
Authors: Martínez-Sykora, Juan; Szydlarski, Mikolaj; Hansteen,
   Viggo H.; De Pontieu, Bart
2020ApJ...900..101M    Altcode: 2020arXiv200800069M
  The solar atmosphere is composed of many species, which are populated
  at different ionization and excitation levels. The upper chromosphere,
  transition region, and corona are nearly collisionless. Consequently,
  slippage between, for instance, ions and neutral particles, or
  interactions between separate species, may play important roles. We
  have developed a 3D MFMS numerical code (Ebysus) to investigate such
  effects. Ebysus is capable of treating species (e.g., hydrogen,
  helium, etc.) and fluids (neutrals, excited and ionized elements)
  separately, including nonequilibrium ionization, momentum exchange,
  radiation, thermal conduction, and other complex processes in the solar
  atmosphere. Treating different species as different fluids leads to
  drifts between different ions and an electric field that couples these
  motions. The coupling for two ionized fluids can lead to an anti-phase
  rotational motion between them. Different ionized species and momentum
  exchange can dissipate this velocity drift, i.e., convert wave kinetic
  energy into thermal energy. High-frequency Alfvén waves driven by,
  e.g., reconnection thought to occur in the solar atmosphere, can drive
  such multi-ion velocity drifts.

---------------------------------------------------------
Title: High-resolution observations of the solar photosphere,
    chromosphere, and transition region. A database of coordinated IRIS
    and SST observations
Authors: Rouppe van der Voort, L. H. M.; De Pontieu, B.; Carlsson,
   M.; de la Cruz Rodríguez, J.; Bose, S.; Chintzoglou, G.; Drews, A.;
   Froment, C.; Gošić, M.; Graham, D. R.; Hansteen, V. H.; Henriques,
   V. M. J.; Jafarzadeh, S.; Joshi, J.; Kleint, L.; Kohutova, P.;
   Leifsen, T.; Martínez-Sykora, J.; Nóbrega-Siverio, D.; Ortiz, A.;
   Pereira, T. M. D.; Popovas, A.; Quintero Noda, C.; Sainz Dalda, A.;
   Scharmer, G. B.; Schmit, D.; Scullion, E.; Skogsrud, H.; Szydlarski,
   M.; Timmons, R.; Vissers, G. J. M.; Woods, M. M.; Zacharias, P.
2020A&A...641A.146R    Altcode: 2020arXiv200514175R
  NASA's Interface Region Imaging Spectrograph (IRIS) provides
  high-resolution observations of the solar atmosphere through ultraviolet
  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 spectropolarimetric imaging
  in photospheric Fe I lines and spectrally resolved imaging in the
  chromospheric Ca II 8542 Å, Hα, 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, plages, the quiet Sun,
  and coronal holes.

---------------------------------------------------------
Title: Ambipolar diffusion in the Bifrost code
Authors: Nóbrega-Siverio, D.; Martínez-Sykora, J.; Moreno-Insertis,
   F.; Carlsson, M.
2020A&A...638A..79N    Altcode: 2020arXiv200411927N
  Context. Ambipolar diffusion is a physical mechanism related to the
  drift between charged and neutral particles in a partially ionized
  plasma that is key to 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. <BR /> Aims: 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. <BR /> Methods: 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 the main feature of the module, we have
  implemented the super time stepping (STS) technique, which allows an
  important acceleration of the calculations. We have also introduced
  hyperdiffusion terms to guarantee the stability of the code. <BR />
  Results: 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.

---------------------------------------------------------
Title: The Formation Height of Millimeter-wavelength Emission in
    the Solar Chromosphere
Authors: Martínez-Sykora, Juan; De Pontieu, Bart; de la Cruz
   Rodriguez, Jaime; Chintzoglou, Georgios
2020ApJ...891L...8M    Altcode: 2020arXiv200110645M
  In the past few years, the ALMA radio telescope has become available
  for solar observations. ALMA diagnostics of the solar atmosphere are of
  high interest because of the theoretically expected linear relationship
  between the brightness temperature at millimeter wavelengths and
  the local gas temperature in the solar atmosphere. Key for the
  interpretation of solar ALMA observations is understanding where in
  the solar atmosphere the ALMA emission originates. Recent theoretical
  studies have suggested that ALMA bands at 1.2 (band 6) and 3 mm
  (band 3) form in the middle and upper chromosphere at significantly
  different heights. We study the formation of ALMA diagnostics using
  a 2.5D radiative MHD model that includes the effects of ion-neutral
  interactions (ambipolar diffusion) and nonequilibrium ionization
  of hydrogen and helium. Our results suggest that in active regions
  and network regions, observations at both wavelengths most often
  originate from similar heights in the upper chromosphere, contrary to
  previous results. Nonequilibrium ionization increases the opacity in the
  chromosphere so that ALMA mostly observes spicules and fibrils along the
  canopy fields. We combine these modeling results with observations from
  IRIS, SDO, and ALMA to suggest a new interpretation for the recently
  reported "dark chromospheric holes," regions of very low temperatures
  in the chromosphere.

---------------------------------------------------------
Title: Ion-neutral Interactions and Nonequilibrium Ionization in
    the Solar Chromosphere
Authors: Martínez-Sykora, Juan; Leenaarts, Jorrit; De Pontieu,
   Bart; Nóbrega-Siverio, Daniel; Hansteen, Viggo H.; Carlsson, Mats;
   Szydlarski, Mikolaj
2020ApJ...889...95M    Altcode: 2019arXiv191206682M
  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 nonequilibrium ionization (NEI) of hydrogen and helium and
  ion-neutral interaction effects. We have performed a 2.5D radiative
  magnetohydrodynamic simulation using the Bifrost code. This model
  includes ion-neutral interaction effects by solving the generalized
  Ohm' s law (GOL) as well as NEI for hydrogen and helium. 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.

---------------------------------------------------------
Title: Nonequilibrium ionization and ambipolar diffusion in solar
    magnetic flux emergence processes
Authors: Nóbrega-Siverio, D.; Moreno-Insertis, F.; Martínez-Sykora,
   J.; Carlsson, M.; Szydlarski, M.
2020A&A...633A..66N    Altcode: 2019arXiv191201015N
  Context. Magnetic flux emergence from the solar interior has
  been shown to be a key mechanism for unleashing a wide variety of
  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. <BR /> Aims: 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. <BR /> Methods: Using the radiation-magnetohydrodynamic
  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 H<SUB>2</SUB> molecules, and the ambipolar
  diffusion term of the generalized Ohm's law. <BR /> Results: 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 H<SUB>2</SUB> molecules
  within the emerged region, thus mistakenly magnifying the exothermic
  contribution of the H<SUB>2</SUB> 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. <P
  />Movies associated to Figs. 2-5, 8, 9, and A.1 are available at <A
  href="https://www.aanda.org/10.1051/0004-6361/201936944/olm">https://www.aanda.org</A>

---------------------------------------------------------
Title: The Multi-slit Approach to Coronal Spectroscopy with the
    Multi-slit Solar Explorer (MUSE)
Authors: De Pontieu, Bart; Martínez-Sykora, Juan; Testa, Paola;
   Winebarger, Amy R.; Daw, Adrian; Hansteen, Viggo; Cheung, Mark C. M.;
   Antolin, Patrick
2020ApJ...888....3D    Altcode: 2019arXiv190908818D
  The Multi-slit Solar Explorer (MUSE) is a proposed mission aimed
  at understanding the physical mechanisms driving the heating of the
  solar corona and the eruptions that are at the foundation of space
  weather. MUSE contains two instruments, a multi-slit extreme ultraviolet
  (EUV) spectrograph and a context imager. It will simultaneously
  obtain EUV spectra (along 37 slits) and context images with the
  highest resolution in space (0.″33-0.″4) and time (1-4 s) ever
  achieved for the transition region (TR) and corona. The MUSE science
  investigation will exploit major advances in numerical modeling, and
  observe at the spatial and temporal scales on which competing models
  make testable and distinguishable predictions, thereby leading to a
  breakthrough in our understanding of coronal heating and the drivers
  of space weather. By obtaining spectra in four bright EUV lines (Fe
  IX 171 Å, Fe XV 284 Å, Fe XIX 108Å, Fe XXI 108 Å) covering a wide
  range of TR and coronal temperatures along 37 slits simultaneously,
  MUSE will be able to “freeze” the evolution of the dynamic
  coronal plasma. We describe MUSE’s multi-slit approach and show
  that the optimization of the design minimizes the impact of spectral
  lines from neighboring slits, generally allowing line parameters to
  be accurately determined. We also describe a Spectral Disambiguation
  Code to resolve multi-slit ambiguity in locations where secondary lines
  are bright. We use simulations of the corona and eruptions to perform
  validation tests and show that the multi-slit disambiguation approach
  allows accurate determination of MUSE observables in locations where
  significant multi-slit contamination occurs.

---------------------------------------------------------
Title: Ebysus: a multi-fluid and multi-species numerical code:
    on coupling between ionized species
Authors: Martínez-Sykora, J.; Szydlarski, M.; Hansteen, V. H.
2019AGUFMSH33D3412M    Altcode:
  The solar atmosphere is composed of many species which are populated at
  different ionized/excited levels. The upper chromosphere, transition
  region and corona are nearly collisionless. Consequently, ion-neutral
  interaction effects or interactions between species may play a role. We
  have developed a 3D multi-fluid and multi-species numerical code
  (Ebysus) to investigate this scenario. Ebysus is capable of treating
  species (e.g., hydrogen, helium etc) and fluids (neutrals, excited and
  ionized elements) separately including non-equilibrium ionization,
  momentum exchange, radiation, thermal conduction, and other complex
  processes in the solar atmosphere. The development of the Ebysus
  code started from the already existing and advanced 3D radiative MHD
  Bifrost code designed to model the upper convection zone and outer
  solar atmosphere. Treating different species as different fluids leads
  to the presence of an electric coupling due to drifts on different ion
  velocities. This coupling leads to anti-correlated high frequency waves
  between the different ionized species. Momentum exchange can dissipate
  this kinetic energy. Alfven waves and reconnection can produce such
  high frequency waves and thus result in reconnection instabilities,
  and in addition reduce the reconnection rate.

---------------------------------------------------------
Title: Unfolding Overlappogram Data: Preparing for the COOL-AID
    instrument on Hi-C FLARE
Authors: Winebarger, A. R.; De Pontieu, B.; Cheung, C. M. M.;
   Martinez-Sykora, J.; Hansteen, V. H.; Testa, P.; Golub, L.; Savage,
   S. L.; Samra, J.; Reeves, K.
2019AGUFMSH33A..06W    Altcode:
  During a solar flare, energy released in the corona streams to the solar
  chromosphere, where plasma is heated and then evaporated upward. The
  magnitude of these velocities and their evolution as a function of time
  can provide quantitative information on the magnitude of energy released
  and the method by which it is transported in a solar flare. Measuring
  these velocities, however, is quite challenging. Typically, they are
  measured with single slit spectrometers, where light passing through
  a long but narrow slit is dispersed and emission lines formed across
  a range of temperatures are observed. The main issue with using
  single slit spectrometers to make this measurement is that they are
  rarely pointed at the right place at the right time. Additionally,
  their fields of view are limited by narrow slit widths, and although
  rastering can effectively expand the field of view, it does so at the
  cost of time. This combination means that single slit spectrometers
  cannot adequately capture the evolution of the flare velocities. On
  the contrary, slitless spectrometers can make "overlappograms”,
  which provide both imaging and spectral information over a large field
  of view. However, spatial information from different spectral lines
  can overlap in the dispersion direction, making the data difficult
  to interpret. Furthermore, the spectral resolution of slitless
  spectrometers are limited and typically worse than single-slit
  spectrometers, since no line fitting (and hence sub-pixel sampling) is
  possible. <P />For the next generation of the High-resolution Coronal
  Imager (Hi-C) Rocket Experiment, which we are proposing to launch during
  a solar flare, we are including the COronal OverLapagram - Ancillary
  Imaging Diagnostics (COOL-AID) instrument. COOL-AID is a slitless
  spectrometer based on the COronal Spectrographic Imager in the EUV
  (COSIE) design, but with a narrow passband coating around 12.9 nm (the
  same passband as the primary Hi-C telescope), a spatial resolution of
  ~1"x2", and a velocity resolution of ~5 km/s. The goal of the COOL-AID
  instrument is to determine the velocity associated with the Fe XXI
  12.9 nm spectral line during a solar flare. In this talk, we will
  demonstrate the unfolding method developed by Cheung et al (2019) to
  determine the velocity information from a simulated COOL-AID data set.

---------------------------------------------------------
Title: A comprehensive three-dimensional radiative magnetohydrodynamic
    simulation of a solar flare
Authors: Cheung, M. C. M.; Rempel, M.; Chintzoglou, G.; Chen, F.;
   Testa, P.; Martínez-Sykora, J.; Sainz Dalda, A.; DeRosa, M. L.;
   Malanushenko, A.; Hansteen, V.; De Pontieu, B.; Carlsson, M.; Gudiksen,
   B.; McIntosh, S. W.
2019NatAs...3..160C    Altcode: 2018NatAs...3..160C
  Solar and stellar flares are the most intense emitters of X-rays and
  extreme ultraviolet radiation in planetary systems<SUP>1,2</SUP>. On
  the Sun, strong flares are usually found in newly emerging sunspot
  regions<SUP>3</SUP>. The emergence of these magnetic sunspot groups
  leads to the accumulation of magnetic energy in the corona. When
  the magnetic field undergoes abrupt relaxation, the energy released
  powers coronal mass ejections as well as heating plasma to temperatures
  beyond tens of millions of kelvins. While recent work has shed light
  on how magnetic energy and twist accumulate in the corona<SUP>4</SUP>
  and on how three-dimensional magnetic reconnection allows for rapid
  energy release<SUP>5,6</SUP>, a self-consistent model capturing how
  such magnetic changes translate into observable diagnostics has remained
  elusive. Here, we present a comprehensive radiative magnetohydrodynamics
  simulation of a solar flare capturing the process from emergence to
  eruption. The simulation has sufficient realism for the synthesis of
  remote sensing measurements to compare with observations at visible,
  ultraviolet and X-ray wavelengths. This unifying model allows us to
  explain a number of well-known features of solar flares<SUP>7</SUP>,
  including the time profile of the X-ray flux during flares, origin
  and temporal evolution of chromospheric evaporation and condensation,
  and sweeping of flare ribbons in the lower atmosphere. Furthermore,
  the model reproduces the apparent non-thermal shape of coronal X-ray
  spectra, which is the result of the superposition of multi-component
  super-hot plasmas<SUP>8</SUP> up to and beyond 100 million K.

---------------------------------------------------------
Title: Multi-component Decomposition of Astronomical Spectra by
    Compressed Sensing
Authors: Cheung, Mark C. M.; De Pontieu, Bart; Martínez-Sykora,
   Juan; Testa, Paola; Winebarger, Amy R.; Daw, Adrian; Hansteen, Viggo;
   Antolin, Patrick; Tarbell, Theodore D.; Wuelser, Jean-Pierre; Young,
   Peter; MUSE Team
2019ApJ...882...13C    Altcode: 2019arXiv190203890C
  The signal measured by an astronomical spectrometer may be due to
  radiation from a multi-component mixture of plasmas with a range of
  physical properties (e.g., temperature, Doppler velocity). Confusion
  between multiple components may be exacerbated if the spectrometer
  sensor is illuminated by overlapping spectra dispersed from different
  slits, with each slit being exposed to radiation from a different
  portion of an extended astrophysical object. We use a compressed sensing
  method to robustly retrieve the different components. This method can
  be adopted for a variety of spectrometer configurations, including
  single-slit, multi-slit (e.g., the proposed MUlti-slit Solar Explorer
  mission), and slot spectrometers (which produce overlappograms).

---------------------------------------------------------
Title: Radiative MHD Simulation of a Solar Flare
Authors: Cheung, Mark; Rempel, Matthias D.; Chintzoglou, Georgios;
   Chen, Feng; Testa, Paola; Martinez-Sykora, Juan; Sainz Dalda, Alberto;
   DeRosa, Marc L.; Malanushenko, Anna; Hansteen, Viggo; Carlsson, Mats;
   De Pontieu, Bart; Gudiksen, Boris; McIntosh, Scott W.
2019AAS...23431005C    Altcode:
  We present a radiative MHD simulation of a solar flare. The
  computational domain captures the near-surface layers of the convection
  zone and overlying atmosphere. Inspired by the observed evolution of
  NOAA Active Region (AR) 12017, a parasitic bipolar region is imposed
  to emerge in the vicinity of a pre-existing sunspot. The emergence of
  twisted magnetic flux generates shear flows that create a pre-existing
  flux rope underneath the canopy field of the sunspot. Following erosion
  of the overlying bootstrapping field, the flux rope erupts. Rapid
  release of magnetic energy results in multi-wavelength synthetic
  observables (including X-ray spectra, narrowband EUV images, Doppler
  shifts of EUV lines) that are consistent with flare observations. This
  works suggests the super-position of multi-thermal, superhot (up
  to 100 MK) plasma may be partially responsible for the apparent
  non-thermal shape of coronal X-ray sources in flares. Implications
  for remote sensing observations of other astrophysical objects is also
  discussed. This work is an important stepping stone toward high-fidelity
  data-driven MHD models.

---------------------------------------------------------
Title: On the Origin of the Magnetic Energy in the Quiet Solar
    Chromosphere
Authors: Martínez-Sykora, Juan; Hansteen, Viggo H.; Gudiksen, Boris;
   Carlsson, Mats; De Pontieu, Bart; Gošić, Milan
2019ApJ...878...40M    Altcode: 2019arXiv190404464M
  The presence of magnetic field is crucial in the transport of energy
  through the solar atmosphere. Recent ground-based and space-borne
  observations of the quiet Sun have revealed that magnetic field
  accumulates at photospheric heights, via a local dynamo or from
  small-scale flux emergence events. However, most of this small-scale
  magnetic field may not expand into the chromosphere due to the entropy
  drop with height at the photosphere. Here we present a study that uses
  a high-resolution 3D radiative MHD simulation of the solar atmosphere
  with non-gray and non-LTE radiative transfer and thermal conduction
  along the magnetic field to reveal that (1) the net magnetic flux
  from the simulated quiet photosphere is not sufficient to maintain a
  chromospheric magnetic field (on average), (2) processes in the lower
  chromosphere, in the region dominated by magnetoacoustic shocks,
  are able to convert kinetic energy into magnetic energy, (3) the
  magnetic energy in the chromosphere increases linearly in time until
  the rms of the magnetic field strength saturates at roughly 4-30 G
  (horizontal average) due to conversion from kinetic energy, (4) and
  that the magnetic features formed in the chromosphere are localized
  to this region.

---------------------------------------------------------
Title: Multi-component Decomposition of Astronomical Spectra by
    Compressed Sensing
Authors: Cheung, Mark; De Pontieu, Bart; Martinez-Sykora, Juan; Testa,
   Paola; Winebarger, Amy R.; Daw, Adrian N.; Hansteen, Viggo; Antolin,
   Patrick; Tarbell, Theodore D.; Wuelser, Jean-Pierre; Young, Peter R.
2019AAS...23411603C    Altcode:
  The signal measured by an astronomical spectrometer may be due to
  radiation from a multi-component mixture of plasmas with a range of
  physical properties (e.g. temperature, Doppler velocity). Confusion
  between multiple components may be exacerbated if the spectrometer
  sensor is illuminated by overlapping spectra dispersed from different
  slits, with each slit being exposed to radiation from a different
  portion of an extended astrophysical object. We use a compressed sensing
  method to robustly retrieve the different components. This method can
  be adopted for a variety of spectrometer configurations, including
  single-slit, multi-slit (e.g., the proposed MUlti-slit Solar Explorer
  mission; MUSE) and slot spectrometers (which produce overlappograms).

---------------------------------------------------------
Title: Impact of Type II Spicules in the Corona: Simulations and
    Synthetic Observables
Authors: Martínez-Sykora, Juan; De Pontieu, Bart; De Moortel, Ineke;
   Hansteen, Viggo H.; Carlsson, Mats
2018ApJ...860..116M    Altcode: 2018arXiv180506475M
  The role of type II spicules in the corona has been a much debated topic
  in recent years. This paper aims to shed light on the impact of type
  II spicules in the corona using novel 2.5D radiative MHD simulations,
  including ion-neutral interaction effects with the Bifrost code. We
  find that the formation of simulated type II spicules, driven by
  the release of magnetic tension, impacts the corona in various
  manners. Associated with the formation of spicules, the corona
  exhibits (1) magneto-acoustic shocks and flows, which supply mass
  to coronal loops, and (2) transversal magnetic waves and electric
  currents that propagate at Alfvén speeds. The transversal waves and
  electric currents, generated by the spicule’s driver and lasting
  for many minutes, are dissipated and heat the associated loop. These
  complex interactions in the corona can be connected with blueshifted
  secondary components in coronal spectral lines (red-blue asymmetries)
  observed with Hinode/EIS and SOHO/SUMER, as well as the EUV counterpart
  of type II spicules and propagating coronal disturbances observed with
  the 171 Å and 193 Å SDO/AIA channels.

---------------------------------------------------------
Title: Small-scale Magnetic Flux Emergence in the Quiet Sun
Authors: Moreno-Insertis, F.; Martinez-Sykora, J.; Hansteen, V. H.;
   Muñoz, D.
2018ApJ...859L..26M    Altcode: 2018arXiv180600489M
  Small bipolar magnetic features are observed to appear in the interior
  of individual granules in the quiet Sun, signaling the emergence of
  tiny magnetic loops from the solar interior. We study the origin
  of those features as part of the magnetoconvection process in the
  top layers of the convection zone. Two quiet-Sun magnetoconvection
  models, calculated with the radiation-magnetohydrodynamic (MHD)
  Bifrost code and with domain stretching from the top layers of the
  convection zone to the corona, are analyzed. Using 3D visualization
  as well as a posteriori spectral synthesis of Stokes parameters,
  we detect the repeated emergence of small magnetic elements in the
  interior of granules, as in the observations. Additionally, we identify
  the formation of organized horizontal magnetic sheets covering whole
  granules. Our approach is twofold, calculating statistical properties
  of the system, like joint probability density functions (JPDFs), and
  pursuing individual events via visualization tools. We conclude that
  the small magnetic loops surfacing within individual granules in the
  observations may originate from sites at or near the downflows in the
  granular and mesogranular levels, probably in the first 1 or 1.5 Mm
  below the surface. We also document the creation of granule-covering
  magnetic sheet-like structures through the sideways expansion of a
  small subphotospheric magnetic concentration picked up and pulled out
  of the interior by a nascent granule. The sheet-like structures that we
  found in the models may match the recent observations of Centeno et al.

---------------------------------------------------------
Title: Chromospheric Dynamics and Heating Processes Ion-Neutral
    Effects in the Solar Chromosphere and Type II Spicules
Authors: Haraldson Hansteen, Viggo; Martinez-Sykora, Juan
2018tess.conf40002H    Altcode:
  Three-dimensional (3D) Magnetohydrodynamic (MHD) "realistic" models
  of the solar chromosphere and lower corona have now progressed to
  the point where meaningful comparisons of synthetic observables
  and solar data can be made. These comparisons show that while much
  can be understood in the context of MHD, there are also important
  discrepancies. The solar chromosphere is largely neutral, this leads
  to a whole range of phenomena as neutrals slip with respect to charged
  particles, perhaps enough to resolve some of these issues. For example,
  in the lower solar atmosphere, the chromosphere is permeated by jets
  known as spicules, in which plasma is propelled at speeds of 50-150
  km/s into the corona. The origin of the spicules is poorly understood,
  although they are expected to play a role in heating the million-degree
  corona and are associated with Alfvén waves that help drive the solar
  wind. In this talk spicules are shown to occur when magnetic tension
  is amplified and transported upwards through interactions between ions
  and neutrals. The tension is impulsively released to drive flows, heat
  plasma and generate Alfvén waves. This talk aims to shed light on ion
  neutral effects and in particular the formation and impact of type
  II spicules in the chromosphere and in the corona, using novel 2.5D
  radiative MHD simulations including ion-neutral interaction effects
  with the Bifrost code.

---------------------------------------------------------
Title: Tracing non-vertical acoustic shock propagation in the
    chromosphere
Authors: Schmit, Don; Martinez-Sykora, Juan
2018tess.conf20441S    Altcode:
  We report on preliminary progress toward quantifying the energetic
  impact of acoustic shocks to the chromosphere in magnetic regions. We
  use the Bifrost simulation to track the evolution of shocks in
  both vertical and inclined field regions in 2D and 3D radiative-MHD
  simulations. It is well known that the telltale Ca II H&amp;K shock
  signature is near wing emission (grain) followed by a protracted
  blue-to-red Doppler shift in absorption (sawtooth) for a vertical flux
  tube. Given that the chromospheric magnetic field is highly structured
  and that the formation height of chromospheric lines vary, we build off
  that earlier work by identifying the signature of shocks in a broader
  environmental and spectral context. In particular, we focus on the Mg
  II h&amp;k lines which are observed by IRIS and the Ca II 8542A line
  which is observed by the Swedish Solar Telescope. We plan to use the
  diagnostics derived from the simulations to derive an occurrence rate
  and an energetic deposition rate for shocks in a plage region.

---------------------------------------------------------
Title: On the Importance of the Nonequilibrium Ionization of Si IV
    and O IV and the Line of Sight in Solar Surges
Authors: Nóbrega-Siverio, D.; Moreno-Insertis, F.; Martínez-Sykora,
   J.
2018ApJ...858....8N    Altcode: 2018arXiv180310251N
  Surges are ubiquitous cool ejections in the solar atmosphere that often
  appear associated with transient phenomena like UV bursts or coronal
  jets. Recent observations from the Interface Region Imaging Spectrograph
  show that surges, although traditionally related to chromospheric
  lines, can exhibit enhanced emission in Si IV with brighter spectral
  profiles than for the average transition region (TR). In this paper,
  we explain why surges are natural sites to show enhanced emissivity
  in TR lines. We performed 2.5D radiative-MHD numerical experiments
  using the Bifrost code including the nonequilibrium (NEQ) ionization of
  silicon and oxygen. A surge is obtained as a by-product of magnetic flux
  emergence; the TR enveloping the emerged domain is strongly affected
  by NEQ effects: assuming statistical equilibrium would produce an
  absence of Si IV and O IV ions in most of the region. Studying the
  properties of the surge plasma emitting in the Si IV λ1402.77 and O IV
  λ1401.16 lines, we find that (a) the timescales for the optically thin
  losses and heat conduction are very short, leading to departures from
  statistical equilibrium, and (b) the surge emits in Si IV more and has
  an emissivity ratio of Si IV to O IV larger than a standard TR. Using
  synthetic spectra, we conclude the importance of line-of-sight effects:
  given the involved geometry of the surge, the line of sight can cut the
  emitting layer at small angles and/or cross it multiple times, causing
  prominent, spatially intermittent brightenings in both Si IV and O IV.

---------------------------------------------------------
Title: Bridging the Gap: Capturing the Lyα Counterpart of a Type-II
    Spicule and Its Heating Evolution with VAULT2.0 and IRIS Observations
Authors: Chintzoglou, Georgios; De Pontieu, Bart; Martínez-Sykora,
   Juan; Pereira, Tiago M. D.; Vourlidas, Angelos; Tun Beltran, Samuel
2018ApJ...857...73C    Altcode: 2018arXiv180303405C
  We present results from an observing campaign in support of the
  VAULT2.0 sounding rocket launch on 2014 September 30. VAULT2.0 is a Lyα
  (1216 Å) spectroheliograph capable of providing spectroheliograms at
  high cadence. Lyα observations are highly complementary to the IRIS
  observations of the upper chromosphere and the low transition region
  (TR) but have previously been unavailable. The VAULT2.0 data provide new
  constraints on upper-chromospheric conditions for numerical models. The
  observing campaign was closely coordinated with the IRIS mission. Taking
  advantage of this simultaneous multi-wavelength coverage of target
  AR 12172 and by using state-of-the-art radiative-MHD simulations of
  spicules, we investigate in detail a type-II spicule associated with
  a fast (300 km s<SUP>-1</SUP>) network jet recorded in the campaign
  observations. Our analysis suggests that spicular material exists
  suspended high in the atmosphere but at lower temperatures (seen in
  Lyα) until it is heated and becomes visible in TR temperatures as a
  network jet. The heating begins lower in the spicule and propagates
  upwards as a rapidly propagating thermal front. The front is then
  observed as fast, plane-of-the-sky motion typical of a network jet,
  but contained inside the pre-existing spicule. This work supports
  the idea that the high speeds reported in network jets should not be
  taken as real mass upflows but only as apparent speeds of a rapidly
  propagating heating front along the pre-existing spicule.

---------------------------------------------------------
Title: Surges and Si IV Bursts in the Solar Atmosphere: Understanding
    IRIS and SST Observations through RMHD Experiments
Authors: Nóbrega-Siverio, D.; Martínez-Sykora, J.; Moreno-Insertis,
   F.; Rouppe van der Voort, L.
2017ApJ...850..153N    Altcode: 2017arXiv171008928N
  Surges often appear as a result of the emergence of magnetized
  plasma from the solar interior. Traditionally, they are observed
  in chromospheric lines such as Hα 6563 \mathringA and Ca II 8542
  \mathringA . However, whether there is a response to the surge
  appearance and evolution in the Si IV lines or, in fact, in many
  other transition region lines has not been studied. In this paper,
  we analyze a simultaneous episode of an Hα surge and a Si IV burst
  that occurred on 2016 September 03 in active region AR 12585. To that
  end, we use coordinated observations from the Interface Region Imaging
  Spectrograph and the Swedish 1-m Solar Telescope. For the first time,
  we report emission of Si IV within the surge, finding profiles that
  are brighter and broader than the average. Furthermore, the brightest
  Si IV patches within the domain of the surge are located mainly near
  its footpoints. To understand the relation between the surges and the
  emission in transition region lines like Si IV, we have carried out 2.5D
  radiative MHD (RMHD) experiments of magnetic flux emergence episodes
  using the Bifrost code and including the nonequilibrium ionization of
  silicon. Through spectral synthesis, we explain several features of
  the observations. We show that the presence of Si IV emission patches
  within the surge, their location near the surge footpoints and various
  observed spectral features are a natural consequence of the emergence of
  magnetized plasma from the interior to the atmosphere and the ensuing
  reconnection processes.

---------------------------------------------------------
Title: Observations and Modeling of Transition Region and Coronal
    Heating Associated with Spicules
Authors: De Pontieu, B.; Martinez-Sykora, J.; De Moortel, I.;
   Chintzoglou, G.; McIntosh, S. W.
2017AGUFMSH43A2793D    Altcode:
  Spicules have been proposed as significant contributorsto the coronal
  energy and mass balance. While previous observationshave provided
  a glimpse of short-lived transient brightenings in thecorona that
  are associated with spicules, these observations have beencontested
  and are the subject of a vigorous debate both on the modelingand
  the observational side so that it remains unclear whether plasmais
  heated to coronal temperatures in association with spicules. We use
  high-resolution observations of the chromosphere and transition region
  with the Interface Region Imaging Spectrograph (IRIS) and ofthe corona
  with the Atmospheric Imaging Assembly (AIA) onboard theSolar Dynamics
  Observatory (SDO) to show evidence of the formation of coronal
  structures as a result of spicular mass ejections andheating of
  plasma to transition region and coronaltemperatures. Our observations
  suggest that a significant fraction of the highly dynamic loop fan
  environment associated with plage regions may be the result of the
  formation of such new coronal strands, a process that previously had
  been interpreted as the propagation of transient propagating coronal
  disturbances (PCD)s. Our observationsare supported by 2.5D radiative
  MHD simulations that show heating tocoronal temperatures in association
  with spicules. Our results suggest that heating and strong flows play
  an important role in maintaining the substructure of loop fans, in
  addition to the waves that permeate this low coronal environment. Our
  models also matches observations ofTR counterparts of spicules and
  provides an elegant explanation forthe high apparent speeds of these
  "network jets".

---------------------------------------------------------
Title: Bridging the Gap: Capturing the Lyα Counterpart of a Type-II
    Spicule and its Heating Evolution with VAULT2.0 and IRIS Campaign
    Observations
Authors: Chintzoglou, G.; De Pontieu, B.; Martinez-Sykora, J.; Mendes
   Domingos Pereira, T.; Vourlidas, A.; Tun Beltran, S.
2017AGUFMSH43A2794C    Altcode:
  We present the analysis of data from the observing campaign in support
  to the VAULT2.0 sounding rocket launch on September 30, 2014. VAULT2.0
  is a Lyα (1216 Å) spectroheliograph capable of providing fast
  cadence spectroheliograms of high-spectral purity. High resolution
  Lyα observations are highly complementary with the IRIS observations
  of the upper chromosphere and the low transition region but have
  previously been unavailable. The VAULT2.0 data provide critical, new
  upper-chromospheric constraints for numerical models. The observing
  campaign was closely coordinated with the IRIS mission. Taking
  advantage of this simultaneous multi-wavelength coverage of target
  AR 12172 and by using state-of-the-art radiative-MHD simulations of
  spicules, we are able to perform a detailed investigation of a type-II
  spicule associated with a fast apparent network jet recorded in the
  campaign observations during the VAULT2.0 flight. Our unique analysis
  suggests that spicular material exists suspended in lower temperatures
  until it rapidly gets heated and becomes visible in transition-region
  temperatures as an apparent network jet.

---------------------------------------------------------
Title: Intermittent Reconnection and Plasmoids in UV Bursts in the
    Low Solar Atmosphere
Authors: Rouppe van der Voort, L.; De Pontieu, B.; Scharmer, G. B.;
   de la Cruz Rodríguez, J.; Martínez-Sykora, J.; Nóbrega-Siverio,
   D.; Guo, L. J.; Jafarzadeh, S.; Pereira, T. M. D.; Hansteen, V. H.;
   Carlsson, M.; Vissers, G.
2017ApJ...851L...6R    Altcode: 2017arXiv171104581R
  Magnetic reconnection is thought to drive a wide variety of dynamic
  phenomena in the solar atmosphere. Yet, the detailed physical mechanisms
  driving reconnection are difficult to discern in the remote sensing
  observations that are used to study the solar atmosphere. In this
  Letter, we exploit the high-resolution instruments Interface Region
  Imaging Spectrograph and the new CHROMIS Fabry-Pérot instrument at
  the Swedish 1-m Solar Telescope (SST) to identify the intermittency
  of magnetic reconnection and its association with the formation of
  plasmoids in so-called UV bursts in the low solar atmosphere. The Si IV
  1403 Å UV burst spectra from the transition region show evidence of
  highly broadened line profiles with often non-Gaussian and triangular
  shapes, in addition to signatures of bidirectional flows. Such profiles
  had previously been linked, in idealized numerical simulations, to
  magnetic reconnection driven by the plasmoid instability. Simultaneous
  CHROMIS images in the chromospheric Ca II K 3934 Å line now provide
  compelling evidence for the presence of plasmoids by revealing highly
  dynamic and rapidly moving brightenings that are smaller than 0.″2 and
  that evolve on timescales of the order of seconds. Our interpretation
  of the observations is supported by detailed comparisons with synthetic
  observables from advanced numerical simulations of magnetic reconnection
  and associated plasmoids in the chromosphere. Our results highlight
  how subarcsecond imaging spectroscopy sensitive to a wide range of
  temperatures combined with advanced numerical simulations that are
  realistic enough to compare with observations can directly reveal the
  small-scale physical processes that drive the wide range of phenomena
  in the solar atmosphere.

---------------------------------------------------------
Title: What Causes the High Apparent Speeds in Chromospheric and
    Transition Region Spicules on the Sun?
Authors: De Pontieu, Bart; Martínez-Sykora, Juan; Chintzoglou,
   Georgios
2017ApJ...849L...7D    Altcode: 2017arXiv171006803D
  Spicules are the most ubuiquitous type of jets in the solar
  atmosphere. The advent of high-resolution imaging and spectroscopy
  from the Interface Region Imaging Spectrograph (IRIS) and ground-based
  observatories has revealed the presence of very high apparent motions of
  order 100-300 km s<SUP>-1</SUP> in spicules, as measured in the plane of
  the sky. However, line of sight measurements of such high speeds have
  been difficult to obtain, with values deduced from Doppler shifts in
  spectral lines typically of order 30-70 km s<SUP>-1</SUP>. In this work,
  we resolve this long-standing discrepancy using recent 2.5D radiative
  MHD simulations. This simulation has revealed a novel driving mechanism
  for spicules in which ambipolar diffusion resulting from ion-neutral
  interactions plays a key role. In our simulation, we often see that
  the upward propagation of magnetic waves and electrical currents
  from the low chromosphere into already existing spicules can lead to
  rapid heating when the currents are rapidly dissipated by ambipolar
  diffusion. The combination of rapid heating and the propagation of these
  currents at Alfvénic speeds in excess of 100 km s<SUP>-1</SUP> leads
  to the very rapid apparent motions, and often wholesale appearance,
  of spicules at chromospheric and transition region temperatures. In
  our simulation, the observed fast apparent motions in such jets are
  actually a signature of a heating front, and much higher than the
  mass flows, which are of order 30-70 km s<SUP>-1</SUP>. Our results
  can explain the behavior of transition region “network jets” and
  the very high apparent speeds reported for some chromospheric spicules.

---------------------------------------------------------
Title: Two-dimensional Radiative Magnetohydrodynamic Simulations of
    Partial Ionization in the Chromosphere. II. Dynamics and Energetics
    of the Low Solar Atmosphere
Authors: Martínez-Sykora, Juan; De Pontieu, Bart; Carlsson, Mats;
   Hansteen, Viggo H.; Nóbrega-Siverio, Daniel; Gudiksen, Boris V.
2017ApJ...847...36M    Altcode: 2017arXiv170806781M
  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.

---------------------------------------------------------
Title: Impact of Type II Spicules into the Corona
Authors: Martinez-Sykora, Juan; De Pontieu, Bart; Carlsson, Mats;
   Hansteen, Viggo H.; Pereira, Tiago M. D.
2017SPD....4810403M    Altcode:
  In the lower solar atmosphere, the chromosphere is permeated by jets,
  in which plasma is propelled at speeds of 50-150 km/s into the Sun’s
  atmosphere or corona. Although these spicules may play a role in heating
  the million-degree corona and are associated with Alfvén waves that
  help drive the solar wind, their generation remains mysterious. We
  implemented in the radiative MHD Bifrost code the effects of partial
  ionization using the generalized Ohm’s law. This code also solves
  the full MHD equations with non-grey and non-LTE radiative transfer
  and thermal conduction along magnetic field lines. The ion-neutral
  collision frequency is computed using recent studies that improved the
  estimation of the cross sections under chromospheric conditions (Vranjes
  &amp; Krstic 2013). Self-consistently driven jets (spicules type II)
  in magnetohydrodynamic simulations occur ubiquitously when magnetic
  tension is confined and transported upwards through interactions
  between ions and neutrals, and impulsively released to drive flows,
  heat plasma, generate Alfvén waves, and may play an important role in
  maintaining the substructure of loop fans. This mechanism explains how
  spicular plasma can be heated to millions of degrees and how Alfvén
  waves are generated in the chromosphere.

---------------------------------------------------------
Title: Realistic radiative MHD simulation of a solar flare
Authors: Rempel, Matthias D.; Cheung, Mark; Chintzoglou, Georgios;
   Chen, Feng; Testa, Paola; Martinez-Sykora, Juan; Sainz Dalda, Alberto;
   DeRosa, Marc L.; Viktorovna Malanushenko, Anna; Hansteen, Viggo H.;
   De Pontieu, Bart; Carlsson, Mats; Gudiksen, Boris; McIntosh, Scott W.
2017SPD....4840001R    Altcode:
  We present a recently developed version of the MURaM radiative
  MHD code that includes coronal physics in terms of optically thin
  radiative loss and field aligned heat conduction. The code employs
  the "Boris correction" (semi-relativistic MHD with a reduced speed
  of light) and a hyperbolic treatment of heat conduction, which allow
  for efficient simulations of the photosphere/corona system by avoiding
  the severe time-step constraints arising from Alfven wave propagation
  and heat conduction. We demonstrate that this approach can be used
  even in dynamic phases such as a flare. We consider a setup in which
  a flare is triggered by flux emergence into a pre-existing bipolar
  active region. After the coronal energy release, efficient transport
  of energy along field lines leads to the formation of flare ribbons
  within seconds. In the flare ribbons we find downflows for temperatures
  lower than ~5 MK and upflows at higher temperatures. The resulting
  soft X-ray emission shows a fast rise and slow decay, reaching a peak
  corresponding to a mid C-class flare. The post reconnection energy
  release in the corona leads to average particle energies reaching 50
  keV (500 MK under the assumption of a thermal plasma). We show that
  hard X-ray emission from the corona computed under the assumption of
  thermal bremsstrahlung can produce a power-law spectrum due to the
  multi-thermal nature of the plasma. The electron energy flux into the
  flare ribbons (classic heat conduction with free streaming limit) is
  highly inhomogeneous and reaches peak values of about 3x10<SUP>11</SUP>
  erg/cm<SUP>2</SUP>/s in a small fraction of the ribbons, indicating
  regions that could potentially produce hard X-ray footpoint sources. We
  demonstrate that these findings are robust by comparing simulations
  computed with different values of the saturation heat flux as well as
  the "reduced speed of light".

---------------------------------------------------------
Title: Observations and Numerical Models of Solar Coronal Heating
    Associated with Spicules
Authors: De Pontieu, B.; De Moortel, I.; Martinez-Sykora, J.; McIntosh,
   S. W.
2017ApJ...845L..18D    Altcode: 2017arXiv171006790D
  Spicules have been proposed as significant contributors to the mass
  and energy balance of the corona. While previous observations have
  provided a glimpse of short-lived transient brightenings in the
  corona that are associated with spicules, these observations have
  been contested and are the subject of a vigorous debate both on the
  modeling and the observational side. Therefore, it remains unclear
  whether plasma is heated to coronal temperatures in association with
  spicules. We use high-resolution observations of the chromosphere and
  transition region (TR) with the Interface Region Imaging Spectrograph
  and of the corona with the Atmospheric Imaging Assembly on board
  the Solar Dynamics Observatory to show evidence of the formation
  of coronal structures associated with spicular mass ejections and
  heating of plasma to TR and coronal temperatures. Our observations
  suggest that a significant fraction of the highly dynamic loop fan
  environment associated with plage regions may be the result of the
  formation of such new coronal strands, a process that previously had
  been interpreted as the propagation of transient propagating coronal
  disturbances. Our observations are supported by 2.5D radiative MHD
  simulations that show heating to coronal temperatures in association
  with spicules. Our results suggest that heating and strong flows
  play an important role in maintaining the substructure of loop fans,
  in addition to the waves that permeate this low coronal environment.

---------------------------------------------------------
Title: On the generation of solar spicules and Alfvénic waves
Authors: Martínez-Sykora, J.; De Pontieu, B.; Hansteen, V. H.;
   Rouppe van der Voort, L.; Carlsson, M.; Pereira, T. M. D.
2017Sci...356.1269M    Altcode: 2017arXiv171007559M
  In the lower solar atmosphere, the chromosphere is permeated by jets
  known as spicules, in which plasma is propelled at speeds of 50 to
  150 kilometers per second into the corona. The origin of the spicules
  is poorly understood, although they are expected to play a role in
  heating the million-degree corona and are associated with Alfvénic
  waves that help drive the solar wind. We compare magnetohydrodynamic
  simulations of spicules with observations from the Interface Region
  Imaging Spectrograph and the Swedish 1-m Solar Telescope. Spicules
  are shown to occur when magnetic tension is amplified and transported
  upward through interactions between ions and neutrals or ambipolar
  diffusion. The tension is impulsively released to drive flows, heat
  plasma (through ambipolar diffusion), and generate Alfvénic waves.

---------------------------------------------------------
Title: On the Misalignment between Chromospheric Features and the
    Magnetic Field on the Sun
Authors: Martínez-Sykora, Juan; De Pontieu, Bart; Carlsson, Mats;
   Hansteen, Viggo
2016ApJ...831L...1M    Altcode: 2016arXiv160702551M
  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.

---------------------------------------------------------
Title: The Cool Surge Following Flux Emergence in a Radiation-MHD
    Experiment
Authors: Nóbrega-Siverio; D.; Moreno-Insertis, F.; Martínez-Sykora,
   J.
2016usc..confE..68N    Altcode:
  Cool and dense ejections, typically Hα surges, often appear alongside
  EUV or X-ray coronal jets as a result of the emergence of magnetized
  plasma from the solar interior. Idealized numerical experiments explain
  those ejections as being indirectly associated with the magnetic
  reconnection taking place between the emerging and preexisting
  systems. However, those experiments miss basic elements that can
  importantly affect the surge phenomenon. In this paper we study the
  cool surges using a realistic treatment of the radiation transfer and
  material plasma properties. To that end, the Bifrost code is used,
  which has advanced modules for the equation of state of the plasma,
  photospheric and chromospheric radiation transfer, heat conduction,
  and optically thin radiative cooling. We carry out a 2.5D experiment of
  the emergence of magnetized plasma through (meso) granular convection
  cells and the low atmosphere to the corona. Through detailed Lagrange
  tracing we study the formation and evolution of the cool ejection and,
  in particular, the role of the entropy sources; this allows us to
  discern families of evolutionary patterns for the plasma elements. In
  the launch phase, many elements suffer accelerations well in excess
  of gravity; when nearing the apex of their individual trajectories,
  instead, the plasma elements follow quasi-parabolic trajectories with
  accelerations close to the solar gravity . We show how the formation
  of the cool ejection is mediated by a wedge-like structure composed
  of two shocks, one of which leads to the detachment of the surge from
  the original emerged plasma dome.

---------------------------------------------------------
Title: The Cool Surge Following Flux Emergence in a Radiation-MHD
    Experiment
Authors: Nóbrega-Siverio, D.; Moreno-Insertis, F.; Martínez-Sykora,
   J.
2016ApJ...822...18N    Altcode: 2016arXiv160104074N
  Cool and dense ejections, typically Hα surges, often appear alongside
  EUV or X-ray coronal jets as a result of the emergence of magnetized
  plasma from the solar interior. Idealized numerical experiments explain
  those ejections as being indirectly associated with the magnetic
  reconnection taking place between the emerging and preexisting
  systems. However, those experiments miss basic elements that can
  importantly affect the surge phenomenon. In this paper we study the
  cool surges using a realistic treatment of the radiation transfer and
  material plasma properties. To that end, the Bifrost code is used,
  which has advanced modules for the equation of state of the plasma,
  photospheric and chromospheric radiation transfer, heat conduction,
  and optically thin radiative cooling. We carry out a 2.5D experiment of
  the emergence of magnetized plasma through (meso) granular convection
  cells and the low atmosphere to the corona. Through detailed Lagrange
  tracing we study the formation and evolution of the cool ejection and,
  in particular, the role of the entropy sources; this allows us to
  discern families of evolutionary patterns for the plasma elements. In
  the launch phase, many elements suffer accelerations well in excess
  of gravity; when nearing the apex of their individual trajectories,
  instead, the plasma elements follow quasi-parabolic trajectories with
  accelerations close to {g}<SUB>⊙ </SUB>. We show how the formation
  of the cool ejection is mediated by a wedge-like structure composed
  of two shocks, one of which leads to the detachment of the surge from
  the original emerged plasma dome.

---------------------------------------------------------
Title: Physics &amp; Diagnostics of the Drivers of Solar Eruptions
Authors: Cheung, Mark; Rempel, Matthias D.; Martinez-Sykora, Juan;
   Testa, Paola; Hansteen, Viggo H.; Viktorovna Malanushenko, Anna;
   Sainz Dalda, Alberto; DeRosa, Marc L.; De Pontieu, Bart; Carlsson,
   Mats; Chen, Feng; McIntosh, Scott W.; Gudiksen, Boris
2016SPD....47.0607C    Altcode:
  We provide an update on our NASA Heliophysics Grand Challenges Research
  (HGCR) project on the ‘Physics &amp; Diagnostics of the Drivers of
  Solar Eruptions’. This presentation will focus on results from a
  data-inspired, 3D radiative MHD model of a solar flare. The model
  flare results from the interaction of newly emerging flux with a
  pre-existing active region. Synthetic observables from the model
  reproduce observational features compatible with actual flares. These
  include signatures of coronal magnetic reconnection, chromospheric
  evaporation, EUV flare arcades, sweeping motion of flare ribbons
  and sunquakes.

---------------------------------------------------------
Title: Time Dependent Nonequilibrium Ionization of Transition Region
    Lines Observed with IRIS
Authors: Martínez-Sykora, Juan; De Pontieu, Bart; Hansteen, Viggo H.;
   Gudiksen, Boris
2016ApJ...817...46M    Altcode: 2015arXiv151200865M
  The properties of nonstatistical equilibrium ionization of silicon
  and oxygen ions are analyzed in this work. We focus on five solar
  targets (quiet Sun; coronal hole; plage; quiescent active region,
  AR; and flaring AR) as observed with the Interface Region Imaging
  Spectrograph (IRIS). IRIS is best suited for this work owing to the
  high cadence (up to 0.5 s), high spatial resolution (up to 0.″32),
  and high signal-to-noise ratios for O IV λ1401 and Si IV λ1402. We
  find that the observed intensity ratio between lines of three times
  ionized silicon and oxygen ions depends on their total intensity
  and that this correlation varies depending on the region observed
  (quiet Sun, coronal holes, plage, or active regions) and on the
  specific observational objects present (spicules, dynamic loops, jets,
  microflares, or umbra). In order to interpret the observations, we
  compare them with synthetic profiles taken from 2D self-consistent
  radiative MHD simulations of the solar atmosphere, where the
  statistical equilibrium or nonequilibrium treatment of silicon and
  oxygen is applied. These synthetic observations show vaguely similar
  correlations to those in the observations, I.e., between the intensity
  ratios and their intensities, but only in the nonequilibrium case do
  we find that (some of) the observations can be reproduced. We conclude
  that these lines are formed out of statistical equilibrium. We use
  our time-dependent nonequilibrium ionization simulations to describe
  the physical mechanisms behind these observed properties.

---------------------------------------------------------
Title: Impact of the Ion-Neutral Interaction Effects in the Solar
    Chromosphere
Authors: Martínez-Sykora, J.; De Pontieu, B.; Hansteen, V. H.;
   Carlsson, M.
2015AGUFMSH31B2411M    Altcode:
  The complexity of the chromosphere is due to various regime changes
  that take place across it. Consequently, the interpretation of
  chromospheric observations is a challenging task. It is thus crucial
  to combine these observations with advanced radiative-MHD numerical
  modeling. Because the photosphere, chromosphere and transition region
  are partially ionized, the interaction between ionized and neutral
  particles has important consequences on the magneto-thermodynamics
  of these regions. We implemented the effects of partial ionization
  using generalized Ohm's law in the Bifrost code (Gudiksen et al. 2011)
  which solves the full MHD equations with non-grey and non-LTE radiative
  transfer and thermal conduction along magnetic field lines. We perform
  2.5D simulations which combines large and small scales structures. This
  leads to a highly dynamic chromosphere with large variety of physical
  processes which have not been reproduced with smaller simulations. The
  implementation of partial ionization effects impact our modeled
  radiative-MHD atmosphere, such as producing chromospheric heating and
  diffusion of photospheric magnetic field into the upper-chromosphere. We
  will also focus on which observables of these processes can be revealed
  with chromospheric observations.

---------------------------------------------------------
Title: Multi-parametric Study of Rising 3D Buoyant Flux Tubes in an
    Adiabatic Stratification Using AMR
Authors: Martínez-Sykora, Juan; Moreno-Insertis, Fernando; Cheung,
   Mark C. M.
2015ApJ...814....2M    Altcode: 2015arXiv150701506M
  We study the buoyant rise of magnetic flux tubes embedded in
  an adiabatic stratification using two-and three-dimensional,
  magnetohydrodynamic simulations. We analyze the dependence of the tube
  evolution on the field line twist and on the curvature of the tube axis
  in different diffusion regimes. To be able to achieve a comparatively
  high spatial resolution we use the FLASH code, which has a built-in
  Adaptive Mesh Refinement (AMR) capability. Our 3D experiments reach
  Reynolds numbers that permit a reasonable comparison of the results
  with those of previous 2D simulations. When the experiments are run
  without AMR, hence with a comparatively large diffusivity, the amount
  of longitudinal magnetic flux retained inside the tube increases
  with the curvature of the tube axis. However, when a low-diffusion
  regime is reached by using the AMR algorithms, the magnetic twist is
  able to prevent the splitting of the magnetic loop into vortex tubes
  and the loop curvature does not play any significant role. We detect
  the generation of vorticity in the main body of the tube of opposite
  sign on the opposite sides of the apex. This is a consequence of the
  inhomogeneity of the azimuthal component of the field on the flux
  surfaces. The lift force associated with this global vorticity makes
  the flanks of the tube move away from their initial vertical plane in
  an antisymmetric fashion. The trajectories have an oscillatory motion
  superimposed, due to the shedding of vortex rolls to the wake, which
  creates a Von Karman street.

---------------------------------------------------------
Title: The role of partial ionization effects in the chromosphere
Authors: Martínez-Sykora, Juan; De Pontieu, Bart; Hansteen, Viggo;
   Carlsson, Mats
2015RSPTA.37340268M    Altcode: 2015arXiv150302723M
  The energy for the coronal heating must be provided from the
  convection zone. However, the amount and the method by which this
  energy is transferred into the corona depend on the properties of the
  lower atmosphere and the corona itself. We review: (i) how the energy
  could be built in the lower solar atmosphere, (ii) how this energy is
  transferred through the solar atmosphere, and (iii) how the energy is
  finally dissipated in the chromosphere and/or corona. Any mechanism of
  energy transport has to deal with the various physical processes in the
  lower atmosphere. We will focus on a physical process that seems to
  be highly important in the chromosphere and not deeply studied until
  recently: the ion-neutral interaction effects in the chromosphere. We
  review the relevance and the role of the partial ionization in the
  chromosphere and show that this process actually impacts considerably
  the outer solar atmosphere. We include analysis of our 2.5D radiative
  magnetohydrodynamic simulations with the Bifrost code (Gudiksen et
  al. 2011 Astron. Astrophys. 531, A154 (doi:10.1051/0004-6361/201116520))
  including the partial ionization effects on the chromosphere
  and corona and thermal conduction along magnetic field lines. The
  photosphere, chromosphere and transition region are partially ionized
  and the interaction between ionized particles and neutral particles
  has important consequences on the magneto-thermodynamics of these
  layers. The partial ionization effects are treated using generalized
  Ohm's law, i.e. we consider the Hall term and the ambipolar diffusion
  (Pedersen dissipation) in the induction equation. The interaction
  between the different species affects the modelled atmosphere as
  follows: (i) the ambipolar diffusion dissipates magnetic energy and
  increases the minimum temperature in the chromosphere and (ii) the
  upper chromosphere may get heated and expanded over a greater range
  of heights. These processes reveal appreciable differences between
  the modelled atmospheres of simulations with and without ion-neutral
  interaction effects.

---------------------------------------------------------
Title: Internetwork Chromospheric Bright Grains Observed With IRIS
    and SST
Authors: Martínez-Sykora, Juan; Rouppe van der Voort, Luc; Carlsson,
   Mats; De Pontieu, Bart; Pereira, Tiago M. D.; Boerner, Paul; Hurlburt,
   Neal; Kleint, Lucia; Lemen, James; Tarbell, Ted D.; Title, Alan;
   Wuelser, Jean-Pierre; Hansteen, Viggo H.; Golub, Leon; McKillop, Sean;
   Reeves, Kathy K.; Saar, Steven; Testa, Paola; Tian, Hui; Jaeggli,
   Sarah; Kankelborg, Charles
2015ApJ...803...44M    Altcode: 2015arXiv150203490M
  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 1403 Å) spectral imaging, and single-wavelength Stokes maps
  in Fe i 6302 Å at high spatial (0\buildrel{\prime\prime}\over{.}
  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 (SJIs) 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 Å SJIs and we conclude that the grain
  signal in these passbands comes mostly from the continuum. Although
  weak, the characteristic shock signatures of acoustic grains can often
  be detected in IRIS C ii spectra. For some grains, a 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.

---------------------------------------------------------
Title: Homologous Helical Jets: Observations By IRIS, SDO, and Hinode
    and Magnetic Modeling With Data-Driven Simulations
Authors: Cheung, Mark C. M.; De Pontieu, B.; Tarbell, T. D.; Fu, Y.;
   Tian, H.; Testa, P.; Reeves, K. K.; Martínez-Sykora, J.; Boerner,
   P.; Wülser, J. P.; Lemen, J.; Title, A. M.; Hurlburt, N.; Kleint,
   L.; Kankelborg, C.; Jaeggli, S.; Golub, L.; McKillop, S.; Saar, S.;
   Carlsson, M.; Hansteen, V.
2015ApJ...801...83C    Altcode: 2015arXiv150101593C
  We report on observations of recurrent jets by instruments on board
  the Interface Region Imaging Spectrograph, Solar Dynamics Observatory
  (SDO), and Hinode spacecraft. Over a 4 hr period on 2013 July 21,
  recurrent coronal jets were observed to emanate from NOAA Active Region
  11793. Far-ultraviolet spectra probing plasma at transition region
  temperatures show evidence of oppositely directed flows with components
  reaching Doppler velocities of ±100 km s<SUP>-1</SUP>. 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.

---------------------------------------------------------
Title: Why is Non-Thermal Line Broadening of Spectral Lines in the
    Lower Transition Region of the Sun Independent of Spatial Resolution?
Authors: De Pontieu, B.; McIntosh, S.; Martinez-Sykora, J.; Peter,
   H.; Pereira, T. M. D.
2015ApJ...799L..12D    Altcode: 2017arXiv171006807D
  Spectral observations of the solar transition region (TR) and
  corona show broadening of spectral lines beyond what is expected
  from thermal and instrumental broadening. The remaining non-thermal
  broadening is significant (5-30 km s<SUP>-1</SUP>) and correlated with
  intensity. Here we study spectra of the TR Si iv 1403 Å line obtained
  at high resolution with the Interface Region Imaging Spectrograph
  (IRIS). We find that the large improvement in spatial resolution
  (0.″33) of IRIS compared to previous spectrographs (2″) does
  not resolve the non-thermal line broadening which, in most regions,
  remains at pre-IRIS levels of about 20 km s<SUP>-1</SUP>. This
  invariance to spatial resolution indicates that the processes behind
  the broadening occur along the line-of-sight (LOS) and/or on spatial
  scales (perpendicular to the LOS) smaller than 250 km. Both effects
  appear to play a role. Comparison with IRIS chromospheric observations
  shows that, in regions where the LOS is more parallel to the field,
  magneto-acoustic shocks driven from below impact the TR and can lead
  to significant non-thermal line broadening. This scenario is supported
  by MHD simulations. While these do not show enough non-thermal line
  broadening, they do reproduce the long-known puzzling correlation
  between non-thermal line broadening and intensity. This correlation
  is caused by the shocks, but only if non-equilibrium ionization is
  taken into account. In regions where the LOS is more perpendicular
  to the field, the prevalence of small-scale twist is likely to play
  a significant role in explaining the invariance and correlation with
  intensity.

---------------------------------------------------------
Title: Observables of Ion-Neutral Interaction Effects in the Solar
    Chromosphere
Authors: Martínez-Sykora, J.; De Pontieu, B.; Hansteen, V. H.;
   Pereira, T. M. D.; Leenaarts, J.; Carlsson, M.
2014AGUFMSH51C4176M    Altcode:
  The chromosphere and transition region constitute the interface
  between the solar surface and the corona and modulate the flow of
  mass and energy into the upper atmosphere. IRIS was launched in 2013
  to study the chromosphere and transition region. The complexity of the
  chromosphere is due to various regime changes that take place across it,
  like: Hydrogen goes from predominantly neutral to predominantly ionized;
  the plasma behavior changes from collisional to collision-less; it goes
  from gas-pressure dominated to magnetically driven, etc. Consequently,
  the interpretation of chromospheric observations in general and those
  from IRIS, in particular, is a challenging task. It is thus crucial
  to combine IRIS observations with advanced radiative-MHD numerical
  modeling. Because the photosphere, chromosphere and transition region
  are partially ionized, the interaction between ionized and neutral
  particles has important consequences on the magneto-thermodynamics of
  these regions. We implemented the effects of partial ionization using
  generalized Ohm's law in the Bifrost code (Gudiksen et al. 2011) which
  solves the full MHD equations with non-grey and non-LTE radiative
  transfer and thermal conduction along magnetic field lines. The
  implementation of partial ionization effects impact our modeled
  radiative-MHD atmosphere, such as producing chromospheric heating and
  diffusion of photospheric magnetic field into the upper-chromosphere. We
  will focus on which observables of these processes can be revealed
  with IRIS.

---------------------------------------------------------
Title: Why Is Non-thermal Line Broadening of Lower Transition Region
    Lines Independent of Spatial Resolution?
Authors: De Pontieu, B.; Mcintosh, S. W.; Martínez-Sykora, J.; Peter,
   H.; Pereira, T. M. D.
2014AGUFMSH51C4175D    Altcode:
  Spectral observations of the solar transition region (TR) and
  corona typically show broadening of the spectral lines beyond what
  is expected from thermal and instrumental broadening. The remaining
  non-thermal broadening is significant (10-30 km/s), correlated with
  the intensity, and has been attributed to waves, macro and micro
  turbulence, nanoflares, etc... Here we study spectra of the low
  TR Si IV 1403 Angstrom line obtained at high spatial and spectral
  resolution with the Interface Region Imaging Spectrograph (IRIS). We
  find that the large improvement in spatial resolution (0.33 arcsec)
  of IRIS compared to previous spectrographs (2 arcsec) does not resolve
  the non-thermal line broadening which remains at pre-IRIS levels of
  20 km/s. This surprising invariance to spatial resolution indicates
  that the physical processes behind the non-thermal line broadening
  either occur along the line-of-sight (LOS) and/or on spatial scales
  (perpendicular to the LOS) smaller than 250 km. Both effects appear
  to play a role. Comparison with IRIS chromospheric observations
  shows that, in regions where the LOS is more parallel to the field,
  magneto-acoustic shocks driven from below impact the low TR leading to
  strong non-thermal line broadening from line-of-sight integration across
  the shock at the time of impact. This scenario is confirmed by advanced
  MHD simulations. In regions where the LOS is perpendicular to the field,
  the prevalence of small-scale twist is likely to play a significant
  role in explaining the invariance and the correlation with intensity.

---------------------------------------------------------
Title: Hot explosions in the cool atmosphere of the Sun
Authors: Peter, H.; Tian, H.; Curdt, W.; Schmit, D.; Innes, D.;
   De Pontieu, B.; Lemen, J.; Title, A.; Boerner, P.; Hurlburt, N.;
   Tarbell, T. D.; Wuelser, J. P.; Martínez-Sykora, Juan; Kleint,
   L.; Golub, L.; McKillop, S.; Reeves, K. K.; Saar, S.; Testa, P.;
   Kankelborg, C.; Jaeggli, S.; Carlsson, M.; Hansteen, V.
2014Sci...346C.315P    Altcode: 2014arXiv1410.5842P
  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 the Interface Region Imaging Spectrograph (IRIS) reveal
  that it is difficult to determine what is up and down, even in the cool
  6000-kelvin photosphere just above the solar surface: This region hosts
  pockets of hot plasma transiently heated to almost 100,000 kelvin. 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.

---------------------------------------------------------
Title: The unresolved fine structure resolved: IRIS observations of
    the solar transition region
Authors: Hansteen, V.; De Pontieu, B.; Carlsson, M.; Lemen, J.; Title,
   A.; Boerner, P.; Hurlburt, N.; Tarbell, T. D.; Wuelser, J. P.; Pereira,
   T. M. D.; De Luca, E. E.; Golub, L.; McKillop, S.; Reeves, K.; Saar,
   S.; Testa, P.; Tian, H.; Kankelborg, C.; Jaeggli, S.; Kleint, L.;
   Martínez-Sykora, J.
2014Sci...346E.315H    Altcode: 2014arXiv1412.3611H
  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 time
  scales 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.

---------------------------------------------------------
Title: Evidence of nonthermal particles in coronal loops heated
    impulsively by nanoflares
Authors: Testa, P.; De Pontieu, B.; Allred, J.; Carlsson, M.; Reale,
   F.; Daw, A.; Hansteen, V.; Martinez-Sykora, J.; Liu, W.; DeLuca, E. E.;
   Golub, L.; McKillop, S.; Reeves, K.; Saar, S.; Tian, H.; Lemen, J.;
   Title, A.; Boerner, P.; Hurlburt, N.; Tarbell, T. D.; Wuelser, J. P.;
   Kleint, L.; Kankelborg, C.; Jaeggli, S.
2014Sci...346B.315T    Altcode: 2014arXiv1410.6130T
  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 (~20 to
  60 seconds) of intensity and velocity on small spatial scales (≲500
  kilometers) at the footpoints of hot and dynamic coronal loops. The
  observations are consistent with numerical simulations of heating by
  beams of nonthermal electrons, which are generated in small impulsive
  (≲30 seconds) heating events called “coronal nanoflares.” The
  accelerated electrons deposit a sizable fraction of their energy
  (≲10<SUP>25 </SUP>erg) 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.

---------------------------------------------------------
Title: Prevalence of small-scale jets from the networks of the solar
    transition region and chromosphere
Authors: Tian, H.; DeLuca, E. E.; Cranmer, S. R.; De Pontieu, B.;
   Peter, H.; Martínez-Sykora, J.; Golub, L.; McKillop, S.; Reeves,
   K. K.; Miralles, M. P.; McCauley, P.; Saar, S.; Testa, P.; Weber,
   M.; Murphy, N.; Lemen, J.; Title, A.; Boerner, P.; Hurlburt, N.;
   Tarbell, T. D.; Wuelser, J. P.; Kleint, L.; Kankelborg, C.; Jaeggli,
   S.; Carlsson, M.; Hansteen, V.; McIntosh, S. W.
2014Sci...346A.315T    Altcode: 2014arXiv1410.6143T
  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 to 250 kilometers per second from the narrow
  bright network lanes of this interface region. These jets have lifetimes
  of 20 to 80 seconds and widths of ≤300 kilometers. They originate from
  small-scale bright regions, often preceded by footpoint brightenings
  and accompanied by transverse waves with amplitudes of ~20 kilometers
  per second. Many jets reach temperatures of at least ~10<SUP>5</SUP>
  kelvin 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.

---------------------------------------------------------
Title: On the prevalence of small-scale twist in the solar
    chromosphere and transition region
Authors: De Pontieu, B.; Rouppe van der Voort, L.; McIntosh, S. W.;
   Pereira, T. M. D.; Carlsson, M.; Hansteen, V.; Skogsrud, H.; Lemen,
   J.; Title, A.; Boerner, P.; Hurlburt, N.; Tarbell, T. D.; Wuelser,
   J. P.; De Luca, E. E.; Golub, L.; McKillop, S.; Reeves, K.; Saar,
   S.; Testa, P.; Tian, H.; Kankelborg, C.; Jaeggli, S.; Kleint, L.;
   Martinez-Sykora, J.
2014Sci...346D.315D    Altcode: 2014arXiv1410.6862D
  The solar chromosphere and transition region (TR) form an interface
  between the Sun’s surface and its hot outer atmosphere. There,
  most of the nonthermal energy that powers the solar atmosphere
  is transformed into heat, although the detailed mechanism remains
  elusive. High-resolution (0.33-arc second) observations with NASA’s
  Interface Region Imaging Spectrograph (IRIS) reveal a chromosphere
  and TR that are replete with twist or torsional motions on sub-arc
  second scales, occurring in active regions, quiet Sun regions, and
  coronal holes alike. We coordinated observations with the Swedish
  1-meter 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.

---------------------------------------------------------
Title: An Interface Region Imaging Spectrograph First View on Solar
    Spicules
Authors: Pereira, T. M. D.; De Pontieu, B.; Carlsson, M.; Hansteen,
   V.; Tarbell, T. D.; Lemen, J.; Title, A.; Boerner, P.; Hurlburt,
   N.; Wülser, J. P.; Martínez-Sykora, J.; Kleint, L.; Golub, L.;
   McKillop, S.; Reeves, K. K.; Saar, S.; Testa, P.; Tian, H.; Jaeggli,
   S.; Kankelborg, C.
2014ApJ...792L..15P    Altcode: 2014arXiv1407.6360P
  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.

---------------------------------------------------------
Title: The Interface Region Imaging Spectrograph (IRIS)
Authors: De Pontieu, B.; Title, A. M.; Lemen, J. R.; Kushner, G. D.;
   Akin, D. J.; Allard, B.; Berger, T.; Boerner, P.; Cheung, M.; Chou,
   C.; Drake, J. F.; Duncan, D. W.; Freeland, S.; Heyman, G. F.; Hoffman,
   C.; Hurlburt, N. E.; Lindgren, R. W.; Mathur, D.; Rehse, R.; Sabolish,
   D.; Seguin, R.; Schrijver, C. J.; Tarbell, T. D.; Wülser, J. -P.;
   Wolfson, C. J.; Yanari, C.; Mudge, J.; Nguyen-Phuc, N.; Timmons,
   R.; van Bezooijen, R.; Weingrod, I.; Brookner, R.; Butcher, G.;
   Dougherty, B.; Eder, J.; Knagenhjelm, V.; Larsen, S.; Mansir, D.;
   Phan, L.; Boyle, P.; Cheimets, P. N.; DeLuca, E. E.; Golub, L.;
   Gates, R.; Hertz, E.; McKillop, S.; Park, S.; Perry, T.; Podgorski,
   W. A.; Reeves, K.; Saar, S.; Testa, P.; Tian, H.; Weber, M.; Dunn, C.;
   Eccles, S.; Jaeggli, S. A.; Kankelborg, C. C.; Mashburn, K.; Pust, N.;
   Springer, L.; Carvalho, R.; Kleint, L.; Marmie, J.; Mazmanian, E.;
   Pereira, T. M. D.; Sawyer, S.; Strong, J.; Worden, S. P.; Carlsson,
   M.; Hansteen, V. H.; Leenaarts, J.; Wiesmann, M.; Aloise, J.; Chu,
   K. -C.; Bush, R. I.; Scherrer, P. H.; Brekke, P.; Martinez-Sykora,
   J.; Lites, B. W.; McIntosh, S. W.; Uitenbroek, H.; Okamoto, T. J.;
   Gummin, M. A.; Auker, G.; Jerram, P.; Pool, P.; Waltham, N.
2014SoPh..289.2733D    Altcode: 2014arXiv1401.2491D; 2014SoPh..tmp...25D
  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, two-second temporal resolution, and 1 km
  s<SUP>−1</SUP> velocity resolution over a field-of-view of up to
  175 arcsec × 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 Å, including bright spectral lines
  formed in the chromosphere (Mg II h 2803 Å and Mg II k 2796 Å) and
  transition region (C II 1334/1335 Å and Si IV 1394/1403 Å). Slit-jaw
  images in four different passbands (C II 1330, Si IV 1400, Mg II k
  2796, and Mg II wing 2830 Å) can be taken simultaneously with spectral
  rasters that sample regions up to 130 arcsec × 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.

---------------------------------------------------------
Title: Detection of Supersonic Downflows and Associated Heating
    Events in the Transition Region above Sunspots
Authors: Kleint, L.; Antolin, P.; Tian, H.; Judge, P.; Testa, P.;
   De Pontieu, B.; Martínez-Sykora, J.; Reeves, K. K.; Wuelser, J. P.;
   McKillop, S.; Saar, S.; Carlsson, M.; Boerner, P.; Hurlburt, N.; Lemen,
   J.; Tarbell, T. D.; Title, A.; Golub, L.; Hansteen, V.; Jaeggli, S.;
   Kankelborg, C.
2014ApJ...789L..42K    Altcode: 2014arXiv1406.6816K
  Interface Region Imaging Spectrograph data allow us to study the solar
  transition region (TR) with an unprecedented spatial resolution of
  0.”33. On 2013 August 30, we observed bursts of high Doppler shifts
  suggesting strong supersonic downflows of up to 200 km s<SUP>-1</SUP>
  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 the
  Atmospheric Imaging Assembly, 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.

---------------------------------------------------------
Title: IRIS observations of the transition region above sunspots:
    oscillations and moving penumbral dots
Authors: Tian, Hui; DeLuca, Ed; Weber, Mark A.; McKillop, Sean;
   Reeves, Kathy; Kleint, Lucia; Martinez-Sykora, Juan; De Pontieu,
   Bart; Carlsson, Mats
2014AAS...22431306T    Altcode:
  NASA's IRIS mission is providing high-cadence and high-resolution
  observations of the solar transition region and chromosphere. We
  present results from IRIS observation of the transition region above
  sunspots. The major findings can be summarized as following: (1) The C
  II and Mg II line profiles are almost Gaussian in the sunspot umbra and
  clearly exhibit a deep reversal at the line center in the plage region,
  suggesting a greatly reduced opacity in the sunspot atmosphere. (2)
  Strongly nonlinear sunspot oscillations can be clearly identified
  in not only the slit jaw images of 2796Å, 1400Å and 1330Å, but
  also in spectra of the bright Mg II, C II and Si IV lines. The Si
  iv oscillation lags those of C ii and Mg ii by 3 and 12 seconds,
  respectively. The temporal evolution of the line core is dominated by
  the following behavior: a rapid excursion to the blue side, accompanied
  by an intensity increase, followed by a linear decrease of the velocity
  to the red side. The maximum intensity slightly lags the maximum blue
  shift in Si iv , whereas the intensity enhancement slightly precedes the
  maximum blue shift in Mg ii . We find a positive correlation between
  the maximum velocity and deceleration. These results are consistent
  with numerical simulations of upward propagating magneto-acoustic
  shock waves. We also demonstrate that the strongly nonlinear line
  width oscillation, reported both previously and here, is spurious. (3)
  Many small-scale bright dots are present in the penumbral filaments and
  light bridges in SJI 1330Å and 1400Å images obtained in high-cadence
  observations. They are usually smaller than 1" and often just a couple
  of pixels wide. Some bright dots show apparent movement with a speed of
  20-60 km/s(either outward or inward). The lifetime of these penumbral
  dots is mostly less than 1 min. The most obvious feature of the Si IV
  profiles in the bright dots is the enhanced line width. Besides that,
  the profile looks normal and no obvious fast flows are detected. The
  bright dots in the light bridges even show oscillation patterns. It's
  not clear whether these oscillations are triggered by the umbral
  oscillations or not.

---------------------------------------------------------
Title: High-resolution Observations of the Shock Wave Behavior for
    Sunspot Oscillations with the Interface Region Imaging Spectrograph
Authors: Tian, H.; DeLuca, E.; Reeves, K. K.; McKillop, S.; De Pontieu,
   B.; Martínez-Sykora, J.; Carlsson, M.; Hansteen, V.; Kleint, L.;
   Cheung, M.; Golub, L.; Saar, S.; Testa, P.; Weber, M.; Lemen, J.;
   Title, A.; Boerner, P.; Hurlburt, N.; Tarbell, T. D.; Wuelser, J. P.;
   Kankelborg, C.; Jaeggli, S.; McIntosh, S. W.
2014ApJ...786..137T    Altcode: 2014arXiv1404.6291T
  We present the first results of sunspot oscillations from observations
  by the Interface Region Imaging Spectrograph. The strongly nonlinear
  oscillation is identified in both the slit-jaw images and the
  spectra of several emission lines formed in the transition region and
  chromosphere. We first apply a single Gaussian fit to the profiles of
  the Mg II 2796.35 Å, C II 1335.71 Å, and Si IV 1393.76 Å lines in the
  sunspot. The intensity change is ~30%. The Doppler shift oscillation
  reveals a sawtooth pattern with an amplitude of ~10 km s<SUP>-1</SUP>
  in Si IV. The Si IV oscillation lags those of C II and Mg II by ~3 and
  ~12 s, respectively. The line width suddenly increases as the Doppler
  shift changes from redshift to blueshift. However, we demonstrate
  that this increase is caused by the superposition of two emission
  components. We then perform detailed analysis of the line profiles at
  a few selected locations on the slit. The temporal evolution of the
  line core is dominated by the following behavior: a rapid excursion
  to the blue side, accompanied by an intensity increase, followed by a
  linear decrease of the velocity to the red side. The maximum intensity
  slightly lags the maximum blueshift in Si IV, whereas the intensity
  enhancement slightly precedes the maximum blueshift in Mg II. We find
  a positive correlation between the maximum velocity and deceleration,
  a result that is consistent with numerical simulations of upward
  propagating magnetoacoustic shock waves.

---------------------------------------------------------
Title: Small scale variability in quiet sun and coronal holes
Authors: Martinez-Sykora, Juan; De Pontieu, Bart
2014cosp...40E2020M    Altcode:
  IRIS (Interface Region Imaging Spectrograph) was launched in 2013 and
  revealed small-scale rapid brightening variations in large regions
  in the quiet sun and coronal holes. Their lifetime is smaller than
  a couple of minutes and the distance between them is of a granular
  scale. We are going to present a statistical study of their observed
  properties such as variability, lifetime, frequency using the images
  of the SJI, and velocities, and temperatures using the spectral data
  coming IRIS. The spectrograph reveals that most of this emission
  comes from the continuum and their properties reveals that most of
  these events result from acoustic shocks. We compare our observations
  with synthetic observables (using forward modeling with the RH and
  multi3D code) from recent numerical 3D radiative-MHD simulations using
  the Bifrost code (Gudiksen et al. 2011). Similar rapid brightening
  variability of chromospheric synthetic images has been also reproduced
  in our simulations with mainly unipolar field. We will describe their
  evolution, how they are driven and their thermodynamic properties from
  the simulations.

---------------------------------------------------------
Title: Measuring energy flux of magneto-acoustic wave in the magnetic
    elements by using IRIS
Authors: Kato, Yoshiaki; De Pontieu, Bart; Martinez-Sykora, Juan;
   Hansteen, Viggo; Pereira, Tiago; Leenaarts, Jorritt; Carlsson, Mats
2014cosp...40E1423K    Altcode:
  NASA's Interface Region Imaging Spectrograph (IRIS) has opened a new
  window to explore the chromospheric/coronal waves that potentially
  energize the solar atmosphere. By using an imaging spectrograph covering
  the Si IV and Mg II h&amp;k lines as well as a slit-jaw imager centered
  at Si IV and Mg II k onboard IRIS, we can determine the nature of
  propagating magneto-acoustic waves just below and in the transition
  region. In this study, we compute the vertically emergent intensity of
  the Si IV and Mg II h&amp;k lines from a time series of snapshots of
  a magnetic element in a two-dimensional Radiative MHD simulation from
  the Bifrost code. We investigate the synthetic line profiles to detect
  the slow magneto-acoustic body wave (slow mode) which becomes a slow
  shock at the lower chromosphere in the magnetic element. We find that
  the Doppler shift of the line core gives the velocity amplitude of the
  longitudinal magneto-acoustic body wave. The contribution function of
  the line core indicates that the formation of Mg II h&amp;k lines is
  associated with the propagating shocks and therefore the time evolution
  of the line core intensity represents the propagating shocks projected
  on the optical surface. We will report on measurement of the energy
  flux of slow modes in the magnetic elements by using IRIS observations.

---------------------------------------------------------
Title: Impact of the Partial Ionization in the solar atmosphere
    using 2.5D Radiative MHD Simulations
Authors: Martinez-Sykora, Juan; De Pontieu, Bart; Hansteen, Viggo;
   Carlsson, Mats
2014cosp...40E2019M    Altcode:
  The chromosphere/transition region constitute the interface between
  the solar surface and the corona and modulate the flow of mass and
  energy into the upper atmosphere. IRIS was launched in 2013 to study the
  chromosphere and transition region. The complexity of the chromosphere
  is due to various regime changes that take place across it, like:
  Hydrogen goes from predominantly neutral to predominantly ionized;
  the plasma behavior changes from collisional to collision-less; it goes
  from gas-pressure dominated to magnetically driven, etc. Consequently,
  the interpretation of chromospheric observations in general and those
  from IRIS, in particular, is a challenging task. It is thus crucial
  to combine IRIS observations with advanced radiative-MHD numerical
  modeling. Because the photosphere, chromosphere and transition region
  are partially ionized, the interaction between ionized and neutral
  particles has important consequences on the magneto-thermodynamics
  of these regions. We implemented the effects of partial ionization
  using generalized Ohm's law in the Bifrost code (Gudiksen et al. 2011)
  which includes full MHD equations with non-grey and non-LTE radiative
  transfer and thermal conduction along magnetic field lines. I will
  describe the importance and impact of taking into account partial
  ionization effects in the modeled radiative-MHD atmosphere, such as
  chromospheric heating, photospheric magnetic field diffused into the
  upper-chromosphere which expands into the upper atmosphere filling
  the corona with mass, magnetic flux, energy and current, etc.

---------------------------------------------------------
Title: A Detailed Comparison between the Observed and Synthesized
    Properties of a Simulated Type II Spicule
Authors: Martínez-Sykora, Juan; De Pontieu, Bart; Leenaarts, Jorrit;
   Pereira, Tiago M. D.; Carlsson, Mats; Hansteen, Viggo; Stern, Julie
   V.; Tian, Hui; McIntosh, Scott W.; Rouppe van der Voort, Luc
2013ApJ...771...66M    Altcode: 2013arXiv1305.2397M
  We have performed a three-dimensional 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's Solar Optical Telescope. 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 a wide
  variety of synthetic emission and absorption lines from the simulations
  including chromospheric (Ca II 8542 Å, Ca II H, and Hα) to transition
  region and coronal temperatures (10,000 K to several million K). We
  compare their synthetic intensities, line profiles, Doppler shifts,
  line widths, and asymmetries with observations from Hinode/SOT and
  EIS, SOHO/SUMER, the Swedish 1 m Solar Telescope, 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 line-of-sight superposition in the optically thin atmosphere
  requires the combination of Doppler shifts and spectral line asymmetry
  to determine the velocity in the jet. In our simulated type II spicule,
  the lifetime of the asymmetry of the transition region lines is shorter
  than that of the coronal lines. 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. The synthetic
  Ca II 8542 Å and Hα profiles therefore do not show signal resembling
  RBEs. 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.

---------------------------------------------------------
Title: Observing Coronal Nanoflares in Active Region Moss
Authors: Testa, Paola; De Pontieu, Bart; Martínez-Sykora, Juan;
   DeLuca, Ed; Hansteen, Viggo; Cirtain, Jonathan; Winebarger, Amy;
   Golub, Leon; Kobayashi, Ken; Korreck, Kelly; Kuzin, Sergey; Walsh,
   Robert; DeForest, Craig; Title, Alan; Weber, Mark
2013ApJ...770L...1T    Altcode: 2013arXiv1305.1687T
  The High-resolution Coronal Imager (Hi-C) has provided Fe XII 193Å
  images of the upper transition region moss at an unprecedented spatial
  (~0.”3-0.”4) and temporal (5.5 s) resolution. The Hi-C observations
  show in some moss regions variability on timescales down to ~15 s,
  significantly shorter than the minute-scale variability typically found
  in previous observations of moss, therefore challenging the conclusion
  of moss being heated in a mostly steady manner. These rapid variability
  moss regions are located at the footpoints of bright hot coronal
  loops observed by the Solar Dynamics Observatory/Atmospheric Imaging
  Assembly in the 94 Å channel, and by the Hinode/X-Ray Telescope. The
  configuration of these loops is highly dynamic, and suggestive of
  slipping reconnection. We interpret these events as signatures of
  heating events associated with reconnection occurring in the overlying
  hot coronal loops, i.e., coronal nanoflares. We estimate the order
  of magnitude of the energy in these events to be of at least a few
  10<SUP>23</SUP> erg, also supporting the nanoflare scenario. These
  Hi-C observations suggest that future observations at comparable
  high spatial and temporal resolution, with more extensive temperature
  coverage, are required to determine the exact characteristics of the
  heating mechanism(s).

---------------------------------------------------------
Title: Modeling small-scale flux emergence from the Convection Zone
    into the Corona
Authors: Martinez-Sykora, Juan
2013enss.confE..60M    Altcode:
  High resolution telescopes reveal small-scale flux in the photosphere
  and roughly 20% of these events seem to reach and impact the
  chromosphere. As a result of such flux emergence, reconnection
  with the ambient field or other processes that do not necessarily
  involve reconnection but nevertheless impact the chromosphere and
  lower corona may occur. I am going to present recent simulations
  that show small-scale flux emergence in a computational domain that
  captures the upper-convection zone, photosphere, chromosphere and lower
  corona. As we will see, small scale activity is strongly dependent on
  the physics that dominate in the various layers of the atmosphere,
  such as thermo-dynamics, radiative transfer in the photosphere and
  thermal conduction along field lines in the corona (we use for that
  Bifrost). In addition, small scale activity is also dependent on the
  ambient field which changes rapidly with height both in strength and
  topology through the different layers of the solar atmosphere. Some
  of these small-scale events erupts into the atmosphere destabilizes
  the pre-existing magnetic field and drives it to new configurations.

---------------------------------------------------------
Title: Small scale activity in the solar atmosphere.
Authors: Martinez-Sykora, J.
2012AGUFMSH53B..01M    Altcode:
  High resolution telescopes reveal a large variety of small-scale
  activity in the photosphere and chromosphere. These processes can be
  driven by convective motion and small scale flux emergence that in
  some cases can penetrate into the chromosphere and higher layers. As
  result of such flux emergence, reconnection with the ambient field,
  or other processes that do not necessarily involve reconnection but
  nevertheless impact the chromosphere and lower corona will occur. I
  will review the most recent simulations that describe the various
  small-scale processes that could impact the chromosphere and corona. As
  we will see small scale activity is strongly dependent on the physics
  that dominate in the various layers of the atmosphere. These processes
  include radiative transfer in the photosphere and chromosphere, partial
  ionization effects and time-dependent ionization in the chromosphere,
  and thermal conduction along field lines in the corona. In addition,
  small scale activity is also dependent on the ambient field which
  changes rapidly with height both in strength and topology through the
  different layers of the solar atmosphere.

---------------------------------------------------------
Title: Coupling of the chromosphere and corona: What physics is
    required?
Authors: Martinez-Sykora, J.
2012AGUFMSH31B..02M    Altcode:
  The chromosphere is the interface region between the solar surface
  and the corona. The observations indicate that the chromosphere is
  feeding the corona with mass, momentum, energy, dynamics and, of course,
  magnetic flux. It is therefore crucial to implement the most important
  physical processes active in these layers in order to build a coherent
  physical model of their coupling. These processes include, but are not
  limited to, radiative transfer, partial ionization, time-dependent
  Hydrogen ionization and thermal conduction along field lines. The
  Bifrost code takes these processes into account and in this talk I
  will discuss the impact and the importance of the partial ionization
  effects and time-dependent ionization. We find that ambipolar (Pedersen)
  dissipation is important in the chromosphere in these radiative-MHD
  simulations. Moreover, ambipolar dissipation is strongly dependent
  on electron density, and the ionization state must be calculated by
  taking the time-dependent Hydrogen ionization into account. As a result
  this effect impacts the chromosphere in terms of thermal properties,
  dynamics and magnetic evolution. Because of the coupling between the
  chromosphere and corona, we will show that these physical processes
  in the chromosphere also impact the corona.

---------------------------------------------------------
Title: Current status of self-consistent 3D radiative-MHD simulations
    of the solar atmosphere
Authors: Martinez-Sykora, J.
2012IAUSS...6E.105M    Altcode:
  In recent years, there has been major progress in the development of
  self-consistent models of the solar atmosphere. These simulations aim to
  capture the physics from the convection zone to the corona. However, the
  solar chromosphere, which is an interface layer between the photosphere
  and corona, has been difficult to model realistically. This is because
  the chromosphere is dominated by transitions: from optically thin to
  optically thick radiation, from continuum to energetically important
  atomic transitions, from gas pressure to magnetic pressure dominance. In
  addition, the chromosphere is where non-LTE and non-grey radiative
  transport are important, hydrogen ionization is time dependent,
  the plasma is partially ionized, etc. As a result, the chromosphere
  is highly dynamic and complex, and filled with dynamical features
  that may play a significant role in the mass and energy balance of
  the corona and solar wind. I will describe the physics included in
  current state-of-the-art numerical codes, and their importance for the
  thermodynamics of the solar atmosphere. I will focus on recent work that
  investigates the formation of spicule-like features, and the effects of
  partial ionization on the dynamics and energetics of the chromosphere.

---------------------------------------------------------
Title: Investigating the Reliability of Coronal Emission Measure
    Distribution Diagnostics using Three-dimensional Radiative
    Magnetohydrodynamic Simulations
Authors: Testa, Paola; De Pontieu, Bart; Martínez-Sykora, Juan;
   Hansteen, Viggo; Carlsson, Mats
2012ApJ...758...54T    Altcode: 2012arXiv1208.4286T
  Determining the temperature distribution of coronal plasmas can provide
  stringent constraints on coronal heating. Current observations with
  the Extreme ultraviolet Imaging Spectrograph (EIS) on board Hinode
  and the Atmospheric Imaging Assembly (AIA) on board 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 three-dimensional 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 the
  signal-to-noise ratio 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 and (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.

---------------------------------------------------------
Title: Two Types of Spicules “Observed” in 3D Realistic Models
Authors: Martínez-Sykora, J.
2012ASPC..454..133M    Altcode:
  Realistic numerical 3D models of the outer solar atmosphere show two
  different kind of spicule-like phenomena, as also observed on the solar
  limb. The numerical models are calculated using the Oslo Staggered Code
  (OSC) to solve the full MHD equations with non-grey and NLTE radiative
  transfer and thermal conduction along the magnetic field lines. The two
  types of spicules arise as a natural result of the dynamical evolution
  in the models. We discuss the different properties of these two types
  of spicules, their differences from observed spicules and what needs
  to be improved in the models.

---------------------------------------------------------
Title: Two-dimensional Radiative Magnetohydrodynamic Simulations of
    the Importance of Partial Ionization in the Chromosphere
Authors: Martínez-Sykora, Juan; De Pontieu, Bart; Hansteen, Viggo
2012ApJ...753..161M    Altcode: 2012arXiv1204.5991M
  The bulk of the solar chromosphere is weakly ionized and interactions
  between ionized particles and neutral particles likely have significant
  consequences for the thermodynamics of the chromospheric plasma. We
  investigate the importance of introducing neutral particles into the
  MHD equations using numerical 2.5D radiative MHD simulations obtained
  with the Bifrost code. The models span the solar atmosphere from the
  upper layers of the convection zone to the low corona, and solve the
  full MHD equations with non-gray and non-LTE radiative transfer, and
  thermal conduction along the magnetic field. The effects of partial
  ionization are implemented using the generalized Ohm's law, i.e.,
  we consider the effects of the Hall term and ambipolar diffusion
  in the induction equation. The approximations required in going
  from three fluids to the generalized Ohm's law are tested in our
  simulations. The Ohmic diffusion, Hall term, and ambipolar diffusion
  show strong variations in the chromosphere. These strong variations
  of the various magnetic diffusivities are absent or significantly
  underestimated when, as has been common for these types of studies,
  using the semi-empirical VAL-C model as a basis for estimates. In
  addition, we find that differences in estimating the magnitude
  of ambipolar diffusion arise depending on which method is used to
  calculate the ion-neutral collision frequency. These differences
  cause uncertainties in the different magnetic diffusivity terms. In
  the chromosphere, we find that the ambipolar diffusion is of the same
  order of magnitude or even larger than the numerical diffusion used
  to stabilize our code. As a consequence, ambipolar diffusion produces
  a strong impact on the modeled atmosphere. Perhaps more importantly,
  it suggests that at least in the chromospheric domain, self-consistent
  simulations of the solar atmosphere driven by magnetoconvection can
  accurately describe the impact of the dominant form of resistivity,
  i.e., ambipolar diffusion. This suggests that such simulations may
  be more realistic in their approach to the lower solar atmosphere
  (which directly drives the coronal volume) than previously assumed.

---------------------------------------------------------
Title: Using 3D MHD realistic simulations of the solar corona to
    test plasma diagnostics
Authors: Testa, P.; De Pontieu, B.; Martinez-Sykora, J.; Hansteen,
   V.; Carlsson, M.
2012decs.confE..27T    Altcode:
  We synthesize coronal images and spectra from advanced 3D MHD
  simulations obtained from the state-of-the art Bifrost code, and
  explore how well they reproduce coronal observations with SDO/AIA and
  Hinode/EIS. We apply standard diagnostic techniques (e.g., density, and
  temperature diagnostics) to the synthetic observations and investigate
  how accurately the derived physical information matches the plasma
  parameters of the model. We discuss the limitations of the diagnostics
  and their implications.

---------------------------------------------------------
Title: Importance of the partial ionization in the chromosphere
    using 2D radiative-MHD simulations
Authors: Martinez-Sykora, Juan; De Pontieu, Bart; Hansteen, Viggo H.
2012decs.confE..81M    Altcode:
  The bulk of the solar chromosphere is weakly ionized and interactions
  between ionized particles and neutral particles will have significant
  consequences for the thermodynamics of the chromospheric plasma. We
  investigate the importance of introducing neutral particles into the
  MHD equations using numerical 2.5D radiative MHD simulations obtained
  with the Bifrost code. The models span the solar atmosphere from
  upper layers of the convection zone to the low corona, and solve the
  full MHD equations with non-grey and non-LTE radiative transfer and
  thermal conduction along the magnetic field. The effects of partial
  ionization are implemented using the generalized Ohm's law, i.e.,
  we consider the effects of the Hall and ambipolar diffusion in the
  induction equation. The ohmic, Hall, and ambipolar diffusivities show
  variations of several orders of magnitude in the chromosphere. These
  strong variations of the various magnetic diffusivities are absent
  and significantly underestimated when using the semi-empirical VAL-C
  model as a basis for estimates. We find that in the chromosphere,
  the ambipolar diffusion is of the same order of magnitude or even
  larger than the numerical diffusion used to stabilize our code. As
  result of this, we can study the effects of it in the simulations. The
  ambipolar diffusion produces strong impact on the chromosphere changing
  the thermal properties, dynamics and magnetic field evolution.

---------------------------------------------------------
Title: Observation, inversion and numerical simulation of single-lobed
    Stokes V profiles in the quiet sun.
Authors: Sainz Dalda, A.; Martínez-Sykora, J.; Bellot Rubio, L.;
   Title, A.
2012decs.confE..89S    Altcode:
  We have studied characteristics and statistics of strong asymmetric
  profiles in Stokes V, i.e., single-lobed profiles, in quiet sun using
  Hinode/SOT. These profiles require the existence of a velocity gradient
  along the line-of-sight, possibly associated with gradients of magnetic
  field strength, inclination and/or azimuth. For a better understanding,
  observations, inversions and numerical simulations are compared. We
  focus our analysis of the observations on the statistical properties
  of the single-lobed Stokes V profiles and the results provided by the
  inversions using SIRJUMP, which is an LTE inversion code that can
  reproduce sharp discontinuities or jump in the magnetic field and
  line-of-sight velocity of the atmosphere model. In the quiet sun,
  magnetic field is continuously appearing and disappearing at small
  scales due to the convective motions and the input of new flux from
  deeper layers. From radiative MHD 3D simulations, using Bifrost code, we
  note that most of these small scale processes have stratifications with
  gradients of magnetic field strength, inclination and velocities. As
  result, those stratifications showing jumps in the magnetic field
  configuration are associated with the existence of single-lobe Stokes
  V profiles in the solar photosphere, as we previously assumed for the
  inversions. We show that most of these profiles come from emerging and
  disappearance magnetic flux in small scales in the simulations. Finally,
  we emphasize importance of the comparison between the synthetic profiles
  from the simulations with the observed ones and the atmospheres that
  produce them. This comparison will ultimately improve the realism of
  the simulations and quantify the emerging and disappearance flux in
  the quiet sun.

---------------------------------------------------------
Title: Study of Single-lobed Circular Polarization Profiles in the
    Quiet Sun
Authors: Sainz Dalda, A.; Martínez-Sykora, J.; Bellot Rubio, L.;
   Title, A.
2012ApJ...748...38S    Altcode: 2012arXiv1202.0593S
  The existence of asymmetries in the circular polarization (Stokes V)
  profiles emerging from the solar photosphere has been known since
  the 1970s. These profiles require the presence of a velocity gradient
  along the line of sight (LOS), possibly associated with gradients of
  magnetic field strength, inclination, and/or azimuth. We have focused
  our study on the Stokes V profiles showing extreme asymmetry in the
  form of only one lobe. Using Hinode spectropolarimetric measurements,
  we have performed a statistical study of the properties of these
  profiles in the quiet Sun. We show their spatial distribution, their
  main physical properties, how they are related with several physical
  observables, and their behavior with respect to their position on
  the solar disk. The single-lobed Stokes V profiles occupy roughly
  2% of the solar surface. For the first time, we have observed their
  temporal evolution and have retrieved the physical conditions of the
  atmospheres from which they emerged using an inversion code implementing
  discontinuities of the atmospheric parameters along the LOS. In
  addition, we use synthetic Stokes profiles from three-dimensional
  magnetoconvection simulations to complement the results of the
  inversion. The main features of the synthetic single-lobed profiles
  are in general agreement with the observed ones, lending support to
  the magnetic and dynamic topologies inferred from the inversion. The
  combination of all these different analyses suggests that most of the
  single-lobed Stokes V profiles are signals coming from the magnetic
  flux emergence and/or submergence processes taking place in small
  patches in the photosphere of the quiet Sun.

---------------------------------------------------------
Title: Potential for diagnostics with IRIS and Mg II lines
Authors: Pereira, Tiago M. D.; Carlsson, Mats; Leenaarts, Jorrit;
   Uitenbroek, Han; De Pontieu, Bart; Martinez-Sykora, Juan
2012decs.confE..13P    Altcode:
  The IRIS mission will open up a new window into the solar chromosphere
  and transition region. An important diagnostic that IRIS will bring
  is the Mg II H and K lines. Radiation from these lines is believed
  to be come from a wide range of formation depths, from the higher
  photosphere to the onset of the transition region. With a complex
  formation mechanism, Mg II H and K suffer from departures from LTE
  and partial redistribution (PRD). In this preliminary analysis we will
  look into the potential for diagnostics of Mg II H and K. Using a new
  parallel version of the RH code we synthesised Mg II H and K spectra
  from 3D rMHD simulations of the solar atmosphere. We will discuss
  the relevance of several approximations on the final observables,
  and will compare the Mg II H and K filtergrams with those of Ca II H,
  a robust chromospheric diagnostic line widely used with Hinode/SOT/BFI.

---------------------------------------------------------
Title: Forward Modeling of Emission in Solar Dynamics
    Observatory/Atmospheric Imaging Assembly Passbands from Dynamic
    Three-dimensional Simulations
Authors: Martínez-Sykora, Juan; De Pontieu, Bart; Testa, Paola;
   Hansteen, Viggo
2011ApJ...743...23M    Altcode: 2011arXiv1109.0704M
  It is typically assumed that emission in the passbands of the
  Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics
  Observatory (SDO) is dominated by single or several strong lines
  from ions that under equilibrium conditions are formed in a narrow
  range of temperatures. However, most SDO/AIA channels also contain
  contributions from lines of ions that have formation temperatures
  that are significantly different from the "dominant" ion(s). We
  investigate the importance of these lines by forward modeling the
  emission in the SDO/AIA channels with three-dimensional radiative MHD
  simulations of a model that spans the upper layer of the convection
  zone to the low corona. The model is highly dynamic. In addition,
  we pump a steadily increasing magnetic flux into the corona, in
  order to increase the coronal temperature through the dissipation
  of magnetic stresses. As a consequence, the model covers different
  ranges of coronal temperatures as time progresses. The model covers
  coronal temperatures that are representative of plasma conditions in
  coronal holes and quiet Sun. The 131, 171, and 304 Å AIA passbands
  are found to be the least influenced by the so-called non-dominant
  ions, and the emission observed in these channels comes mostly from
  plasma at temperatures near the formation temperature of the dominant
  ion(s). On the other hand, the other channels are strongly influenced
  by the non-dominant ions, and therefore significant emission in these
  channels comes from plasma at temperatures that are different from the
  "canonical" values. We have also studied the influence of non-dominant
  ions on the AIA passbands when different element abundances are assumed
  (photospheric and coronal), and when the effects of the electron
  density on the contribution function are taken into account.

---------------------------------------------------------
Title: Testing coronal plasma diagnostics using 3D MHD models of
    the solar atmosphere
Authors: Testa, P.; Martinez-Sykora, J.; Hansteen, V. H.; De Pontieu,
   B.; Carlsson, M.
2011AGUFMSH53C..06T    Altcode:
  We synthesize coronal images and spectra from advanced 3D radiative
  MHD simulations obtained from the state-of-the-art Bifrost code, and
  explore how well they reproduce coronal observations with SDO/AIA
  and Hinode/EIS and XRT. We apply standard diagnostic techniques
  (e.g., density, temperature, abundance diagnostics) to the synthetic
  observations and investigate how accurately the derived physical
  information matches the plasma parameters of the model. We discuss
  the limitations of the diagnostics and their implications.

---------------------------------------------------------
Title: Two Components of the Coronal Emission Revealed by
    Extreme-Ultraviolet Spectroscopic Observations
Authors: Tian, H.; Mcintosh, S. W.; De Pontieu, B.; Martinez-Sykora,
   J.; Wang, X.; Sechler, M.
2011AGUFMSH33A2027T    Altcode:
  Recent spectroscopic observations have revealed the ubiquitous presence
  of blueward asymmetries of emission lines formed in the solar corona
  and transition region. These asymmetries are most prominent in loop
  footpoint regions, where a clear correlation of the asymmetry with the
  Doppler shift and line width determined from the single-Gaussian fit
  is found. Such asymmetries suggest at least two emission components: a
  primary component accounting for the background emission and a secondary
  component associated with high-speed upflows. The latter has been
  proposed to play a vital role in the coronal heating process and there
  is no agreement on its properties. Here we slightly modify the initially
  developed technique of red-blue (RB) asymmetry analysis and apply it to
  both artificial spectra and spectra observed by the Extreme-ultraviolet
  Imaging Spectrometer on board Hinode, and demonstrate that the secondary
  component usually contributes a few percent of the total emission, has
  a velocity ranging from 50 to 150 km/s, and a Gaussian width comparable
  to that of the primary one in loop footpoint regions. The results of
  the RB asymmetry analysis are then used to guide a double-Gaussian fit
  and we find that the obtained properties of the secondary component
  are generally consistent with those obtained from the RB asymmetry
  analysis. Through a comparison of the location, relative intensity,
  and velocity distribution of the blueward secondary component with
  the properties of the upward propagating disturbances revealed in
  simultaneous images from the Atmospheric Imaging Assembly on board
  the Solar Dynamics Observatory, we find a clear association of the
  secondary component with the propagating disturbances.

---------------------------------------------------------
Title: Two Components of the Solar Coronal Emission Revealed by
    Extreme-ultraviolet Spectroscopic Observations
Authors: Tian, Hui; McIntosh, Scott W.; De Pontieu, Bart;
   Martínez-Sykora, Juan; Sechler, Marybeth; Wang, Xin
2011ApJ...738...18T    Altcode: 2011arXiv1106.1141T
  Recent spectroscopic observations have revealed the ubiquitous presence
  of blueward asymmetries of emission lines formed in the solar corona
  and transition region. These asymmetries are most prominent in loop
  footpoint regions, where a clear correlation of the asymmetry with the
  Doppler shift and line width determined from the single-Gaussian fit
  is found. Such asymmetries suggest at least two emission components: a
  primary component accounting for the background emission and a secondary
  component associated with high-speed upflows. The latter has been
  proposed to play a vital role in the coronal heating process and there
  is no agreement on its properties. Here we slightly modify the initially
  developed technique of red-blue (RB) asymmetry analysis and apply it to
  both artificial spectra and spectra observed by the Extreme-ultraviolet
  Imaging Spectrometer on board Hinode, and demonstrate that the secondary
  component usually contributes a few percent of the total emission,
  and has a velocity ranging from 50 to 150 km s<SUP>-1</SUP> and a
  Gaussian width comparable to that of the primary one in loop footpoint
  regions. The results of the RB asymmetry analysis are then used to
  guide a double-Gaussian fit and we find that the obtained properties of
  the secondary component are generally consistent with those obtained
  from the RB asymmetry analysis. Through a comparison of the location,
  relative intensity, and velocity distribution of the blueward secondary
  component with the properties of the upward propagating disturbances
  revealed in simultaneous images from the Atmospheric Imaging Assembly
  on board the Solar Dynamics Observatory, we find a clear association
  of the secondary component with the propagating disturbances.

---------------------------------------------------------
Title: On the Origin of the Type II Spicules: Dynamic
    Three-dimensional MHD Simulations
Authors: Martínez-Sykora, Juan; Hansteen, Viggo; Moreno-Insertis,
   Fernando
2011ApJ...736....9M    Altcode: 2010arXiv1011.4703M
  Recent high temporal and spatial resolution observations of the
  chromosphere have forced the definition of a new type of spicule, "type
  II's," that are characterized by rising rapidly, having short lives,
  and by fading away at the end of their lifetimes. Here, we report on
  features found in realistic three-dimensional simulations of the outer
  solar atmosphere that resemble the observed type II spicules. These
  features evolve naturally from the simulations as a consequence of
  the magnetohydrodynamical evolution of the model atmosphere. The
  simulations span from the upper layer of the convection zone to
  the lower corona and include the emergence of a horizontal magnetic
  flux. The state-of-art Oslo Staggered Code is used to solve the full
  MHD equations with non-gray and non-LTE radiative transfer and thermal
  conduction along the magnetic field lines. We describe in detail the
  physics involved in a process which we consider a possible candidate
  for the driver mechanism that produces type II spicules. The modeled
  spicule is composed of material rapidly ejected from the chromosphere
  that rises into the corona while being heated. Its source lies in
  a region with large field gradients and intense electric currents,
  which lead to a strong Lorentz force that squeezes the chromospheric
  material, resulting in a vertical pressure gradient that propels the
  spicule along the magnetic field, as well as Joule heating, which
  heats the jet material, forcing it to fade.

---------------------------------------------------------
Title: The stellar atmosphere simulation code Bifrost. Code
    description and validation
Authors: Gudiksen, B. V.; Carlsson, M.; Hansteen, V. H.; Hayek, W.;
   Leenaarts, J.; Martínez-Sykora, J.
2011A&A...531A.154G    Altcode: 2011arXiv1105.6306G
  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. <BR /> Aims:
  The design, implementation and validation of the massively parallel
  numerical code Bifrost for simulating stellar atmospheres from the
  convection zone to the corona. <BR /> 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.

---------------------------------------------------------
Title: What do Spectral Line Profile Asymmetries Tell us About the
    Solar Atmosphere?
Authors: Martínez-Sykora, Juan; De Pontieu, Bart; Hansteen, Viggo;
   McIntosh, Scott W.
2011ApJ...732...84M    Altcode:
  Recently, analysis of solar spectra obtained with the EUV Imaging
  Spectrograph (EIS) onboard the Hinode satellite has revealed the
  ubiquitous presence of asymmetries in transition region (TR) and coronal
  spectral line profiles. These asymmetries have been observed especially
  at the footpoints of coronal loops and have been associated with strong
  upflows that may play a significant role in providing the corona with
  hot plasma. Here, we perform a detailed study of the various processes
  that can lead to spectral line asymmetries, using both simple forward
  models and state-of-the-art three-dimensional radiative MHD simulations
  of the solar atmosphere using the Bifrost code. We describe a novel
  technique to determine the presence and properties of faint secondary
  components in the wings of spectral line profiles. This method is based
  on least-squares fitting of observed so-called R(ed)B(lue) asymmetry
  profiles with pre-calculated RB asymmetry profiles for a wide variety
  of secondary component properties. We illustrate how this method could
  be used to perform reliable double Gaussian fits that are not over- or
  under-constrained. We also find that spectral line asymmetries appear
  in TR and coronal lines that are synthesized from our three-dimensional
  MHD simulations. Our models show that the spectral asymmetries are a
  sensitive measure of the velocity gradient with height in the TR of
  coronal loops. The modeled TR shows a large gradient of velocity that
  increases with height: this occurs as a consequence of ubiquitous,
  episodic heating at low heights in the model atmosphere. We show
  that the contribution function of spectral lines as a function of
  temperature is critical for sensitivity to velocity gradients and thus
  line asymmetries: lines that are formed over a temperature range that
  includes most of the TR are the most sensitive. As a result, lines from
  lithium-like ions (e.g., O VI) are found to be the most sensitive to
  line asymmetries. We compare the simulated line profiles directly with
  line profiles observed in the quiet Sun with SOHO/SUMER and Hinode/EIS
  and find that the shape of the profiles is very similar. In addition,
  the simulated profiles with the strongest blueward asymmetry occur in
  footpoint regions of coronal loops, which is similar to what we observe
  with SUMER and EIS. There is however a significant discrepancy between
  the simulations and observations: the simulated RB asymmetries are
  an order of magnitude smaller than the observations. We discuss the
  possible reasons for this discrepancy. In summary, our analysis shows
  that observations of spectral line asymmetries can provide a powerful
  new diagnostic to help constrain coronal heating models.

---------------------------------------------------------
Title: The Origins of Hot Plasma in the Solar Corona
Authors: De Pontieu, B.; McIntosh, S. W.; Carlsson, M.; Hansteen,
   V. H.; Tarbell, T. D.; Boerner, P.; Martinez-Sykora, J.; Schrijver,
   C. J.; Title, A. M.
2011Sci...331...55D    Altcode:
  The Sun's outer atmosphere, or corona, is heated to millions of degrees,
  considerably hotter than its surface or photosphere. Explanations for
  this enigma typically invoke the deposition in the corona of nonthermal
  energy generated by magnetoconvection. However, the coronal heating
  mechanism remains unknown. We used observations from the Solar Dynamics
  Observatory and the Hinode solar physics mission to reveal a ubiquitous
  coronal mass supply in which chromospheric plasma in fountainlike jets
  or spicules is accelerated upward into the corona, with much of the
  plasma heated to temperatures between ~0.02 and 0.1 million kelvin (MK)
  and a small but sufficient fraction to temperatures above 1 MK. These
  observations provide constraints on the coronal heating mechanism(s)
  and highlight the importance of the interface region between photosphere
  and corona.

---------------------------------------------------------
Title: The role of the chromosphere in filling the corona with hot
    plasma (Invited)
Authors: de Pontieu, B.; McIntosh, S. W.; Carlsson, M.; Hansteen,
   V. H.; Tarbell, T. D.; Boerner, P.; Martinez-Sykora, J.; Schrijver,
   C. J.; Title, A. M.
2010AGUFMSH21C..03D    Altcode:
  We use coordinated observations from the Solar Dynamics Observatory
  (SDO), Hinode and the Swedish Solar Telescope (SST) to show how
  plasma is heated to coronal temperatures from its source in the
  chromosphere. Our observations reveal a ubiquitous mass supply
  for the solar corona in which chromospheric plasma is accelerated
  upward into the corona with much of the plasma heated to transition
  region temperatures, and a small, but significant fraction heated
  to temperatures in excess of 1 million K. Our observations show,
  for the first time, how chromospheric spicules, fountain-like jets
  that have long been considered potential candidates for coronal
  heating, are directly associated with heating of plasma to coronal
  temperatures. These results provide strong physical constraints on
  the mechanism(s) responsible for coronal heating and do not seem
  compatible with current models. The association with chromospheric
  spicules highlights the importance of the interface region between
  the photosphere and corona to gain a full understanding of the coronal
  heating problem.

---------------------------------------------------------
Title: Line profile asymmetries in the transition region: models
    and observations
Authors: Martinez-Sykora, J.; de Pontieu, B.; Hansteen, V. H.;
   McIntosh, S. W.
2010AGUFMSH31A1784M    Altcode:
  Asymmetries in spectral line profiles provide a wealth of
  information on the properties of the emitting plasma along the
  line-of-sight. Asymmetries can be produced by the superposition
  of profiles with different line-of-sight velocities and/or widths
  resulting from the variation of the velocity and/or temperature from
  emission sources along the line of sight. Spectral line asymmetries
  from synthetic transition region and coronal lines constructed
  from realistic 3D models appear similar to those observed with
  Hinode/EIS. The simulations span the upper layer of the convection zone
  to the lower corona and include horizontal magnetic flux emergence. We
  use the state of the art Bifrost code to solve the full MHD equations
  with non-grey and non-LTE radiative transfer and thermal conduction
  along the magnetic field line. Here, we perform a detailed study of
  the various physical, dynamical and observational processes that can
  lead to spectral line asymmetries at the transition region footpoints
  of loops in 3D radiative MHD simulations of the solar atmosphere and
  compare these with observations. Our models show that the spectral
  asymmetries are a sensitive measure of the velocity gradient with
  height in the transition region of coronal loops. In our models the
  TR shows a large gradient of velocity that increases with height:
  this occurs as a natural consequence of ubiquitous, episodic heating
  at low heights in the model atmosphere.

---------------------------------------------------------
Title: Forward modeling of emission in AIA passbands from advanced
    radiative MHD simulations
Authors: de Pontieu, B.; Martinez-Sykora, J.; Hansteen, V. H.
2010AGUFMSH11A1597D    Altcode:
  The emission from many of the passbands observed with the Atmospheric
  Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO)
  is dominated by single or several lines from ions that are formed in a
  narrow range of temperatures (under equilibrium conditions). However,
  most AIA passbands contain contributions from lines of ions that
  have formation temperatures that are significantly different from the
  dominant ion. We investigate the importance of these lines by forward
  modeling of the AIA passband emission from advanced radiative 3D MHD
  simulations calculated with the state of the art Bifrost code. We
  use simulations that span the upper layer of the convection zone to
  the low corona and solve the full magnetohydrodynamic equations with
  non-grey and non-LTE radiative transfer and thermal conduction along the
  magnetic field lines. We find that several of the AIA passbands often
  include significant contributions from plasma at different temperatures
  than the canonical temperature values. We describe under which solar
  conditions in the simulations these discrepancies can typically be
  expected to occur.

---------------------------------------------------------
Title: Comparison Of Observations And Advanced Numerical Simulations
    Of Type II Spicules
Authors: Martinez-Sykora, Juan; De Pontieu, B.; Hansteen, V.;
   Moreno-Insertis, F.
2010AAS...21640306M    Altcode: 2010BAAS...41..878M
  We have performed realistic 3D radiation MHD simulations of the
  solar atmosphere. These simulations show jet-like features that
  are similar to the type II spicules discovered with Hinode's Solar
  Optical Telescope. These type II spicules have been associated with
  so-called rapid blueshifted events (RBE's) on the solar disk, and with
  significant blueward asymmetries in transition region and coronal
  lines at the footpoints of coronal loops (discovered with Hinode's
  EIS). These observational results and their ubiquity suggest they may
  play a significant role in providing the corona with hot plasma. We
  will present a detailed comparison of the properties of the simulated
  jets, with those of type II spicules (observed with Hinode) and RBE's
  (with ground-based instruments). We will present analysis of a wide
  variety of synthetic emission lines from the simulations covering
  temperatures from 10,000 K to several million K, and compare their
  intensities, velocities, line widths and asymmetry with those of the
  observed phenomena. We will also show how the formation mechanism of
  these jets (reconnection at tangential discontinuities) complicates
  efforts to establish a firm link between observations of magnetic
  fields and of chromospheric flows, and suggests that magnetic field
  observations at chromospheric heights may be crucial to establish from
  observations how these jets are formed.

---------------------------------------------------------
Title: 2 types of spicules "observed" in 3D realistic models
Authors: Martínez-Sykora, Juan
2010arXiv1001.1256M    Altcode:
  Realistic numerical 3D models of the outer solar atmosphere show two
  different kind of spicule-like phenomena, as also observed on the solar
  limb. The numerical models are calculated using the 2 types of spicules
  "observed" in 3D realistic models Oslo Staggered Code (OSC) to solve
  the full MHD equations with non-grey and NLTE radiative transfer and
  thermal conduction along the magnetic field lines. The two types of
  spicules arise as a natural result of the dynamical evolution in the
  models. We discuss the different properties of these two types of
  spicules, their differences from observed spicules and what needs to
  be improved in the models.

---------------------------------------------------------
Title: Twisted Flux Tube Emergence from the Convection Zone to the
    Corona. II. Later States
Authors: Martínez-Sykora, Juan; Hansteen, Viggo; Carlsson, Mats
2009ApJ...702..129M    Altcode: 2009arXiv0906.5464M
  Three-dimensional numerical simulations of magnetic flux emergence
  are carried out in a computational domain spanning the upper layers
  of the convection zone to the lower corona. We use the Oslo Staggered
  Code to solve the full magnetohydrodynamic equations with non-gray
  and non-local thermodynamic equilibrium radiative transfer and thermal
  conduction along the magnetic field lines. In this paper, we concentrate
  on the later stages of the simulations and study the evolution of the
  structure of the rising flux in the upper chromosphere and corona, the
  interaction between the emerging flux and the weak coronal magnetic
  field initially present, and the associated dynamics. The flux tube
  injected at the bottom boundary rises to the photosphere where it
  largely remains. However, some parts of the flux tube become unstable
  and expand in patches into the upper chromosphere. The flux rapidly
  expands toward the corona, pushing the coronal and transition region
  material aside, and lifting and maintaining the transition region at
  heights greater than 5 Mm above the photosphere for extensive periods
  of time. The pre-existing magnetic field in the corona and transition
  region is perturbed by the incoming flux and reoriented by a series of
  high Joule heating events. Low-density structures form in the corona,
  while at later times a high-density filamentary structure appears in
  the lower part of the expanding flux. The dynamics of these and other
  structures is discussed. While Joule heating due to the expanding flux
  is episodic, it increases in relative strength as fresh magnetic field
  rises and becomes energetically important in the upper chromosphere and
  corona at later times. Chromospheric, transition region, and coronal
  lines are computed and their response to the perturbation caused by
  the expanding emerging flux is discussed.

---------------------------------------------------------
Title: Spicule-Like Structures Observed in Three-Dimensional Realistic
    Magnetohydrodynamic Simulations
Authors: Martínez-Sykora, Juan; Hansteen, Viggo; De Pontieu, Bart;
   Carlsson, Mats
2009ApJ...701.1569M    Altcode: 2009arXiv0906.4446M
  We analyze features that resemble type I spicules in two different
  three-dimensional 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 pre-existing ambient magnetic field
  strength and in the properties of the inserted flux tube. We use the
  Oslo Staggered Code to solve the full magnetohydrodynamic equations
  with nongray 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 versus
  time profiles, and are dominated by rapid upward motion at speeds
  of 10-30 km s<SUP>-1</SUP>, 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 pre-existing 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.

---------------------------------------------------------
Title: Flux emergence from the convection zone to the corona
Authors: Martínez-Sykora, Juan
2009PhDT.......417M    Altcode:
  No abstract at ADS

---------------------------------------------------------
Title: Twisted Flux Tube Emergence From the Convection Zone to
    the Corona
Authors: Martínez-Sykora, Juan; Hansteen, Viggo; Carlsson, Mats
2008ApJ...679..871M    Altcode: 2007arXiv0712.3854M
  Three-dimensional numerical simulations of a horizontal magnetic flux
  tube emergence with different twist are carried out in a computational
  domain spanning the upper layers of the convection zone to the lower
  corona. We use the Oslo Stagger Code to solve the full MHD equations
  with non-gray, non-LTE radiative transfer and thermal conduction along
  the magnetic lines. A magnetic flux tube is input at the bottom boundary
  into a weakly magnetized atmosphere. The photospheric and chromospheric
  response is described with magnetograms and synthetic continuum as
  well as Ca II H line images and velocity field distributions. In the
  photosphere the granular size increases when the flux tube approaches
  from below, as has been reported previously in the literature. In
  the convective overshoot region, some 200 km above the photosphere,
  adiabatic expansion produces cooling, darker regions with the structure
  of granulation cells. We also find evidence of collapsed granulation
  at the boundaries of the rising flux tube. Once the flux tube has
  crossed the photosphere, bright points related to concentrated magnetic
  field, vorticity, high vertical velocities, and heating by compressed
  material are found at heights up to 500 km above the photosphere. At
  greater heights, in the magnetized chromosphere, the rising flux tube
  produces a large, cool, magnetized bubble that tends to expel the
  usual chromospheric oscillations. In addition, the rising flux tube
  dramatically increases the chromospheric scale height, pushing the
  transition region and corona aside, such that the chromosphere extends
  up to 6 Mm above the photosphere. We find flux tube emergence through
  the photosphere to the lower corona to be a relatively slow process,
  taking of order 1 hr.