Source author record

Juan Martinez-Sykora

Juan Martinez-Sykora appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

12 published item(s)

preprint2021arXiv

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

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

preprint2020arXiv

The Formation Height of Millimeter-wavelength Emission in the Solar Chromosphere

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 mm 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 non-equilibrium 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. Non-equilibrium ionization increases the opacity in the chromosphere so that ALMA mostly observe 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.

preprint2019arXiv

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

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

preprint2019arXiv

The Multi-slit Approach to Coronal Spectroscopy with the Multi-slit Solar Explorer (MUSE)

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 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 arcsec) and time (1-4 s) ever achieved for the transition region 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 4 bright EUV lines (Fe IX 171A, Fe XV 284A, Fe XIX-XXI 108A) covering a wide range of transition region 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.

preprint2015arXiv

Multi-Parametric Study of Rising 3D Buoyant Flux Tubes in an Adiabatic Stratification Using AMR

We study the buoyant rise of magnetic flux tubes embedded in an adiabatic stratification using two-and three-dimensional, MHD 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.

preprint2015arXiv

Time Dependent Non-Equilibrium Ionization of Transition Region Lines Observed with IRIS

The properties of non-statistical equilibrium ionization of silicon and oxygen ions are analyzed in this work. We focus on four 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 due to the high cadence (up to 0.5s), high spatial resolution (up to 0.32"), and high signal to noise ratios for O IV and Si IV. 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, micro-flares 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 non-equilibrium treatment silicon and oxygen is applied. These synthetic observations show vaguely similar correlations as in the observations, i.e. between the intensity ratios and their intensities, but only in the non-equilibrium 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 non-equilibrium ionization simulations to describe the physical mechanisms behind these observed properties.

preprint2014arXiv

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

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

preprint2013arXiv

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

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

preprint2013arXiv

Observing coronal nanoflares in active region moss

The High-resolution Coronal Imager (Hi-C) has provided Fe XII 193A images of the upper transition region moss at an unprecedented spatial (~0.3-0.4 arcsec) and temporal (5.5s) resolution. The Hi-C observations show in some moss regions variability on timescales down to ~15s, 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 SDO/AIA in the 94A channel, and by Hinode/XRT. 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^{23}rg, 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).

preprint2012arXiv

2D Radiative MHD Simulations of the Importance of Partial Ionization in the Chromosphere

The solar chromosphere is weakly ionized and interactions between ionized particles and neutral particles likely have significant consequences for the thermodynamics of the plasma. We investigate the importance of introducing neutral particles using numerical 2.5D radiative MHD simulations obtained with the Bifrost code. The models span from the upper layers of the convection zone to the low corona, and solve the full MHD equations with non-grey and NLTE radiative transfer, and thermal conduction. The effects of partial ionization are implemented using the generalized Ohm's law. The approximations required in going from three fluids to the generalized Ohm's law are tested in our simulations. The Ohmic diffusion, the 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 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 magneto-convection can accurately describe the impact of the dominant form of resistivity (ambipolar diffusion). This suggests that such simulations may be more realistic in their approach to the lower solar atmosphere than previously assumed.

preprint2012arXiv

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

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

preprint2011arXiv

Two components of the coronal emission revealed by EUV spectroscopic observations

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 EUV Imaging Spectrometer onboard 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-1 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 onboard the Solar Dynamics Observatory, we find a clear association of the secondary component with the propagating disturbances.