Showing 1–50 of 98 results for author: Khomenko, E

Auction theory and demography

Authors: O. A. Malafeyev, I. E. Khomenko

Abstract: In economics, there are many ways to describe the interaction between a "seller" and a "buyer". The most common one, with which we interact almost every day, is selling for a fixed price. This option is perfect for selling a mass product, when we have a number of sellers and many buyers, and the price for the product varies depending on the conditions of the relationship between supply and demand.… ▽ More In economics, there are many ways to describe the interaction between a "seller" and a "buyer". The most common one, with which we interact almost every day, is selling for a fixed price. This option is perfect for selling a mass product, when we have a number of sellers and many buyers, and the price for the product varies depending on the conditions of the relationship between supply and demand. Another situation meets us already in markets, where a product can be either mass-produced or more unique, so this option is already closer to the object of our discussion.However, a one-on-one transaction is a much more unstable option, which is why it is also more difficult to model, since it is determined not so much by algorithms as by psychology and the difference in the bargaining ability of the two parties. An even closer example of an auction is price discrimination, when the price for the buyer is determined not only by supply and demand, but also by which group the buyer belongs to. But in this case, the product is not unique, and the final seller is the only one. Thus, we have identified the main auction criteria and their features of the "game". △ Less

Submitted 8 July, 2024; originally announced July 2024.

Hydrodynamic simulations of cool stellar atmospheres with \texttt{MANCHA}

Authors: Andrea Perdomo García, Nikola Vitas, Elena Khomenko, Manuel Collados

Abstract: Aims. Our aim is to test how different binning strategies previously studied in one-dimensional models perform in three-dimensional radiative hydrodynamic simulations of stellar atmospheres. Methods. Realistic box-in-a-star simulations of the near-surface convection and photosphere of three spectral types (G2V, K0V, and M2V) were run with the MANCHA code with grey opacity. After reaching the stati… ▽ More Aims. Our aim is to test how different binning strategies previously studied in one-dimensional models perform in three-dimensional radiative hydrodynamic simulations of stellar atmospheres. Methods. Realistic box-in-a-star simulations of the near-surface convection and photosphere of three spectral types (G2V, K0V, and M2V) were run with the MANCHA code with grey opacity. After reaching the stationary state, one snapshot of each of the three stellar simulations was used to compute the radiative energy exchange rate with grey opacity, opacity binned in four $τ$-bins, and opacity binned in 18 $\{ τ, λ\}$-bins. These rates were compared with the ones computed with opacity distribution functions. Then, stellar simulations were run with grey, four-bin, and 18-bin opacities to see the impact of the opacity setup on the mean stratification of the temperature and its gradient after time evolution. Results. The simulations of main sequence cool stars with the MANCHA code are consistent with those in the literature. For the three stars, the radiative energy exchange rates computed with 18 bins are remarkably close to the ones computed with the opacity distribution functions. The rates computed with four bins are similar to the rates computed with 18 bins, and present a significant improvement with respect to the rates computed with the Rosseland opacity, especially above the stellar surface. The Rosseland mean can reproduce the proper rates in sub-surface layers, but produces large errors for the atmospheric layers of the G2V and K0V stars. In the case of the M2V star, the Rosseland mean fails even in sub-surface layers, owing to the importance of the contribution from molecular lines in the opacity, underestimated by the harmonic mean. Similar conclusions are reached studying the mean stratification of the temperature and its gradient after time evolution. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: Accepted for publication in Astronomy and Astrophysics

The influence of thermal pressure gradients and ionization (im)balance on the ambipolar diffusion and charge-neutral drifts

Authors: M. M. Gómez-Míguez, D. Martínez-Gómez, E. Khomenko, N. Vitas

Abstract: Solar partially ionized plasma is frequently modeled using single-fluid (1F) or two-fluid (2F) approaches. In the 1F case, charge-neutral interactions are often described through ambipolar diffusion, while the 2F model fully considers charge-neutral drifts. Here, we expand the definition of the ambipolar diffusion coefficient to include inelastic collisions (ion/rec) in two cases: a VAL3C 1D model… ▽ More Solar partially ionized plasma is frequently modeled using single-fluid (1F) or two-fluid (2F) approaches. In the 1F case, charge-neutral interactions are often described through ambipolar diffusion, while the 2F model fully considers charge-neutral drifts. Here, we expand the definition of the ambipolar diffusion coefficient to include inelastic collisions (ion/rec) in two cases: a VAL3C 1D model and a 2F simulations of the Rayleigh-Taylor instability (RTI) in a solar prominence thread based on \cite{PopLukKho2021aa, PopLukKho2021ab}. On one side, we evaluate the relative importance of the inelastic contribution, compared to elastic and charge-exchange collisions. On the other side, we compare the contributions of ion/rec, thermal pressure, viscosity, and magnetic forces to the charge-neutral drift velocity of the turbulent flow of the RTI. Our analysis reveals that the contribution of inelastic collisions to the ambipolar diffusion coefficient is negligible across the chromosphere, allowing the classical definition of this coefficient to be safely used in 1F modeling. However, in the transition region, the contribution of inelastic collisions can become as significant as that of elastic collisions. Furthermore, we ascertain that the thermal pressure force predominantly influences the charge-neutral drifts in the RTI model, surpassing the impact of the magnetic force. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: Accepted

Mixing, heating and ion-neutral decoupling induced by Rayleigh-Taylor instability in prominence-corona transition regions

Authors: V. S. Lukin, E. Khomenko, B. Popescu Braileanu

Abstract: This study explores non-linear development of the magnetized Rayleigh-Taylor instability (RTI) in a prominence-corona transition region. Using a two-fluid model of a partially ionized plasma, we compare RTI simulations for several different magnetic field configurations. We follow prior descriptions of the numerical prominence model [Popescu Braileanu et al., 2021a,b, 2023] and explore the charged… ▽ More This study explores non-linear development of the magnetized Rayleigh-Taylor instability (RTI) in a prominence-corona transition region. Using a two-fluid model of a partially ionized plasma, we compare RTI simulations for several different magnetic field configurations. We follow prior descriptions of the numerical prominence model [Popescu Braileanu et al., 2021a,b, 2023] and explore the charged-neutral fluid coupling and plasma heating in a two-dimensional mixing layer for different magnetic field configurations. We also investigate how the shear in magnetic field surrounding a prominence may impact the release of the gravitational potential energy of the prominence material. We show that the flow decoupling is strongest in the plane normal to the direction of the magnetic field, where neutral pressure gradients drive ion-neutral drifts and frictional heating is balanced by adiabatic cooling of the expanding prominence material. We also show that magnetic field within the mixing plane can lead to faster non-linear release of the gravitational energy driving the RTI, while more efficiently heating the plasma via viscous dissipation of associated plasma flows. We relate the computational results to potential observables by highlighting how integrating over under-resolved two-fluid sub-structure may lead to misinterpretation of observational data. △ Less

Submitted 27 February, 2024; v1 submitted 6 January, 2024; originally announced January 2024.

Comments: accepted to Philosophical Transactions A

Journal ref: Phil. Trans. R. Soc. A. 382:20230417 (2024)

Mancha3D code: Multi-purpose Advanced Non-ideal MHD Code for High resolution simulations in Astrophysics

Authors: M. Modestov, E. Khomenko, N. Vitas, A. de Vicente, A. Navarro, P. A. Gonzalez-Morales, M. Collados, T. Felipe, D. Martinez-Gomez, P. Hunana, M. Luna, M. Koll Pistarini, B. Popescu Braileanu, A. Perdomo Garcia, V. Liakh, I. Santamaria, M. M. Gomez Miguez

Abstract: The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic processes in solar/stellar atmospheres. The code includes non-ideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high quality realistic simulations of magneto-convection, as well as idealized simulations of particular processes, suc… ▽ More The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic processes in solar/stellar atmospheres. The code includes non-ideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high quality realistic simulations of magneto-convection, as well as idealized simulations of particular processes, such as wave propagation, instabilities or energetic events. The paper summarizes the equations and methods used in the Mancha3D code. It also describes its numerical stability and parallel performance and efficiency. The code is based on a finite difference discretization and memory-saving Runge-Kutta (RK) scheme. It handles non-ideal effects through super-time stepping and Hall diffusion schemes, and takes into account thermal conduction by solving an additional hyperbolic equation for the heat flux. The code is easily configurable to perform different kinds of simulations. Several examples of the code usage are given. It is demonstrated that splitting variables into equilibrium and perturbation parts is essential for simulations of wave propagation in a static background. A perfectly matched layer (PML) boundary condition built into the code greatly facilitates a non-reflective open boundary implementation. Spatial filtering is an important numerical remedy to eliminate grid-size perturbations enhancing the code stability. Parallel performance analysis reveals that the code is strongly memory bound, which is a natural consequence of the numerical techniques used, such as split variables and PML boundary conditions. Both strong and weak scalings show adequate performance up till several thousands of CPUs. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Two fluid dynamics in solar prominences

Authors: S. J. González Manrique, E. Khomenko, M. Collados, C. Kuckein, T. Felipe, P. Gömöry

Abstract: Solar prominences contain a significant amount of neutral species. The dynamics of the ionised and neutral fluids composing the prominence plasma can be slightly different if the collisional coupling is not strong enough. Large-scale velocities can be quantified by Doppler effect. Small-scale velocities leave their imprint on the width of spectral lines. Here we use one spectral line of ionised an… ▽ More Solar prominences contain a significant amount of neutral species. The dynamics of the ionised and neutral fluids composing the prominence plasma can be slightly different if the collisional coupling is not strong enough. Large-scale velocities can be quantified by Doppler effect. Small-scale velocities leave their imprint on the width of spectral lines. Here we use one spectral line of ionised and two spectral lines of neutral elements to measure the resolved and unresolved velocities in a prominence with the aim to investigate the possible decoupling of the observed charged and neutral species. A prominence was observed with the German Vacuum Tower Telescope on June 17, 2017. Time series consisting of repeated 10-position scans were performed while recording simultaneously the intensity spectra of the Ca II 854.2 nm, Hα 656.28 nm, and HeD3 587.56 nm lines. The line-of-sight velocities and the Doppler width of the three spectral lines were determined at every spatial position and temporal moment. To make sure all spectral lines were sampling the same plasma volume, we applied selection criteria to identify locations with optically thin plasma. The velocities of the three spectral lines turned out to be very similar over this region, with the ionised Ca II showing velocity excursions systematically larger compared to those of the neutral lines of Hα and He I at some moments. The latter were found to be much closer to each other. The analysis of the Doppler widths indicated that the C aII line shows an excess of unresolved motions. The dynamics of the ionised and neutral plasma components in the observed prominence was very close one to the other. The differences found may indicate that a localised decoupling between ions and neutrals may appear at particular spatial locations or instants of time. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: Accepted for publication in Astronomy & Astrophysics, 12 pages, 12 figures

Opacity for realistic 3D MHD simulations of cool stellar atmospheres

Authors: A. Perdomo García, N. Vitas, E. Khomenko, M. Collados, C. Allende Prieto, I. Hubeny, Y. Osorio

Abstract: Context. Realistic 3D time-dependent simulations of stellar near-surface convection employ the opacity binning method for efficient and accurate computation of the radiative energy exchange. The method provides several orders of magnitude of speed-up, but its implementation includes a number of free parameters. Aims. Our aim is to evaluate the accuracy of the opacity binning method as a function o… ▽ More Context. Realistic 3D time-dependent simulations of stellar near-surface convection employ the opacity binning method for efficient and accurate computation of the radiative energy exchange. The method provides several orders of magnitude of speed-up, but its implementation includes a number of free parameters. Aims. Our aim is to evaluate the accuracy of the opacity binning method as a function of the choice of these free parameters. Methods. The monochromatic opacities computed with the SYNSPEC code are used to construct opacity distribution function (ODF) that is then verified through detailed comparison with the results of the ATLAS code. The opacity binning method is implemented with the SYNSPEC opacities for four representative cool main-sequence stellar spectral types (F3V, G2V, K0V, and M2V). Results. The ODFs from SYNSPEC and ATLAS show consistent results for the opacity and bolometric radiative energy exchange rate Q in case of the F, G, and K -- type stars. Significant differences, coming mainly from the molecular line lists, are found for the M -- type star. It is possible to optimise a small number of bins to reduce the deviation of the results coming from the opacity grouping with respect to the ODF for the F, G, and K -- type stars. In the case of the M -- type star, the inclusion of splitting in wavelength is needed in the grouping to get similar results, with a subsequent increase in computing time. In the limit of a large number of bins, the deviation for all the binning configurations tested saturates and the results do not converge to the ODF solution. Due to this saturation, the Q rate cannot be improved by increasing the number of bins to more than about 20 bins. The more effective strategy is to select the optimal location of fewer bins. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: 26 pages, 29 figures. Accepted for publication in section 15 (Numerical methods and codes) of Astronomy and Astrophysics

Journal ref: A&A 675, A160 (2023)

A study of the capabilities for inferring atmospheric information from high-spatial-resolution simulations

Authors: C. Quintero Noda, E. Khomenko, M. Collados, B. Ruiz Cobo, R. Gafeira, N. Vitas, M. Rempel, R. J. Campbell, A. Pastor Yabar, H. Uitenbroek, D. Orozco Suárez

Abstract: In this work, we study the accuracy that can be achieved when inferring the atmospheric information from realistic numerical magneto-hydrodynamic simulations that reproduce the spatial resolution we will obtain with future observations made by the 4m class telescopes DKIST and EST. We first study multiple inversion configurations using the SIR code and the Fe I transitions at 630 nm until we obtai… ▽ More In this work, we study the accuracy that can be achieved when inferring the atmospheric information from realistic numerical magneto-hydrodynamic simulations that reproduce the spatial resolution we will obtain with future observations made by the 4m class telescopes DKIST and EST. We first study multiple inversion configurations using the SIR code and the Fe I transitions at 630 nm until we obtain minor differences between the input and the inferred atmosphere in a wide range of heights. Also, we examine how the inversion accuracy depends on the noise level of the Stokes profiles. The results indicate that when the majority of the inverted pixels come from strongly magnetised areas, there are almost no restrictions in terms of the noise, obtaining good results for noise amplitudes up to 1$\times10^{-3}$ of $I_c$. At the same time, the situation is different for observations where the dominant magnetic structures are weak, and noise restraints are more demanding. Moreover, we find that the accuracy of the fits is almost the same as that obtained without noise when the noise levels are on the order of 1$\times10^{-4}$of $I_c$. We, therefore, advise aiming for noise values on the order of or lower than 5$\times10^{-4}$ of $I_c$ if observers seek reliable interpretations of the results for the magnetic field vector reliably. We expect those noise levels to be achievable by next-generation 4m class telescopes thanks to an optimised polarisation calibration and the large collecting area of the primary mirror. △ Less

Submitted 2 June, 2023; originally announced June 2023.

Comments: 14 pages, 13 figures

Improving the Understanding of Subsurface Structure and Dynamics of Solar Active Regions (A white paper submitted to the decadal survey for solar and space Physics (Heliophysics) -- SSPH 2024-2033)

Authors: S. C. Tripathy, K. Jain, D. Braun, P. Cally, M. Dikpati, T. Felipe, R. Jain, S. Kholikov, E. Khomenko, R. Komm, J. Leibacher, V. Martinez-Pillet, A. Pevtsov, S. P. Rajaguru, M. Roth, H. Uitenbroek, J. Zhao

Abstract: The goal of helioseismology is to provide accurate information about the Sun's interior from the observations of the wave field at its surface. In the last three decades, both global and local helioseismology studies have made significant advances and breakthroughs in solar physics. However, 3-d mapping of the structure and dynamics of sunspots and active regions below the surface has been a chall… ▽ More The goal of helioseismology is to provide accurate information about the Sun's interior from the observations of the wave field at its surface. In the last three decades, both global and local helioseismology studies have made significant advances and breakthroughs in solar physics. However, 3-d mapping of the structure and dynamics of sunspots and active regions below the surface has been a challenging task and are among the longest standing and intriguing puzzles of solar physics due to the complexity of the turbulent and dynamic nature of sunspots. Thus the key problems that need to be addressed during the next decade are: (i) Understanding the wave excitation mechanisms in the quiet Sun and magnetic regions, (ii) Characterizing the wave propagation and transformation in strong and inclined magnetic field regions and understanding the magnetic portals in the chromosphere, (iii) Improving helioseismology techniques and investigating the whole life cycle of active regions, from magnetic flux emergence to dissipation, and (iv) Detecting helioseismic signature of the magnetic flux of active regions before it becomes visible on the surface so as to provide warnings several days before the emergence. For a transformative progress on these problems require full disk, simultaneous Doppler and vector magnetic field measurements of the photosphere up to the chromosphere with a spatial resolution of about 2 arc-sec as well as large-scale radiative MHD simulations of the plasma dynamics from the sub-photosphere to the chromosphere. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: A White Paper Submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) -- SSPH 2024-2033

Numerical simulations of prominence oscillations triggered by external perturbations

Authors: Valeriia Liakh, Manuel Luna, Elena Khomenko

Abstract: Several energetic disturbances have been identified as triggers of the large-amplitude oscillations (LAOs) in prominences. However, the mechanisms for LAOs excitation are not well understood. We aim to study these mechanisms, performing time-dependent numerical simulations in 2.5D and 2D setups using magnetohydrodynamic (MHD) code MANCHA3D. Two types of disturbances are applied to excite prominenc… ▽ More Several energetic disturbances have been identified as triggers of the large-amplitude oscillations (LAOs) in prominences. However, the mechanisms for LAOs excitation are not well understood. We aim to study these mechanisms, performing time-dependent numerical simulations in 2.5D and 2D setups using magnetohydrodynamic (MHD) code MANCHA3D. Two types of disturbances are applied to excite prominence oscillations, such as a perturbation associated with an eruption and the waves caused by an artificial energy release. In the simulation with the eruption, we obtain that it does not produce LAOs in the prominence located in its vicinity. While the erupting flux rope rises, an elongated current sheet forms behind it, which becomes unstable and breaks into plasmoids. The downward-moving plasmoids cause perturbations in the velocity field by merging with the post-reconnection loops. This velocity perturbation propagates in the surroundings and perturbs the nearby prominence. The analysis of the oscillatory motions of the prominence plasma reveals the excitation of small-amplitude oscillations (SAOs), which are a mixture of longitudinal and vertical oscillations. In the simulation with a distant artificial perturbation, a fast-mode shock wave is produced, and it gradually reaches two flux rope prominences at different distances. This shock wave excites vertical LAOs and longitudinal SAOs with similar amplitudes, periods, and damping times in both prominences. Finally, in the experiment with the external triggering of LAOs in a dipped arcade prominence model, we find that, although the vector normal to the front of a fast-mode shock wave is parallel to the spine of the dipped arcade well before the contact, this wave does not excite longitudinal LAOs. When the wave front approaches the prominence, it pushes the dense plasma down, establishing vertical LAOs. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 18 pages, 18 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 673, A154 (2023)

Large ion-neutral drift velocities and plasma heating in partially ionized coronal rain blobs

Authors: David Martínez-Gómez, Ramón Oliver, Elena Khomenko, Manuel Collados

Abstract: In this paper we present a numerical study of the dynamics of partially ionized coronal rain blobs. We use a two-fluid model to perform a high-resolution 2D simulation that takes into account the collisional interaction between the charged and neutral particles contained in the plasma. We follow the evolution of a cold plasma condensation as it falls through an isothermal vertically stratified atm… ▽ More In this paper we present a numerical study of the dynamics of partially ionized coronal rain blobs. We use a two-fluid model to perform a high-resolution 2D simulation that takes into account the collisional interaction between the charged and neutral particles contained in the plasma. We follow the evolution of a cold plasma condensation as it falls through an isothermal vertically stratified atmosphere that represents the much hotter and lighter solar corona. We study the consequences of the different degrees of collisional coupling that are present in the system. On the one hand, we find that at the dense core of the blob there is a very strong coupling and the charged and neutral components of the plasma behave as a single fluid, with negligible drift velocities (of a few cm s^-1). On the other hand, at the edges of the blob the coupling is much weaker and larger drift velocities (of the order of 1 km s^-1) appear. In addition, frictional heating causes large increases of temperature at the transition layers between the blob and the corona. For the first time we show that such large drift velocities and temperature enhancements can develop as a consequence of ion-neutral decoupling associated to coronal rain dynamics. This can lead to enhanced emission coming from the plasma at the coronal rain-corona boundary, which possesses transition region temperature. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 11 pages, 7 figures, and 1 animation. Accepted for publication in The Astrophysical Journal Letters

Ambipolar Diffusion in the Lower Solar Atmosphere: MHD Simulations of a Sunspot

Authors: Conor D. MacBride, David B. Jess, Elena Khomenko, Samuel D. T. Grant

Abstract: Magnetohydrodynamic (MHD) simulations of the solar atmosphere are often performed under the assumption that the plasma is fully ionized. However, in the lower solar atmosphere a reduced temperature often results in only the partial ionization of the plasma. The interaction between the decoupled neutral and ionized components of such a partially ionized plasma produces ambipolar diffusion. To inves… ▽ More Magnetohydrodynamic (MHD) simulations of the solar atmosphere are often performed under the assumption that the plasma is fully ionized. However, in the lower solar atmosphere a reduced temperature often results in only the partial ionization of the plasma. The interaction between the decoupled neutral and ionized components of such a partially ionized plasma produces ambipolar diffusion. To investigate the role of ambipolar diffusion in propagating wave characteristics in the photosphere and chromosphere, we employ the Mancha3D numerical code to model magnetoacoustic waves propagating through the atmosphere immediately above the umbra of a sunspot. We solve the non-ideal MHD equations for data-driven perturbations to the magnetostatic equilibrium and the effect of ambipolar diffusion is investigated by varying the simulation to include additional terms in the MHD equations that account for this process. Analyzing the energy spectral densities for simulations with/without ambipolar diffusion, we find evidence to suggest that ambipolar diffusion plays a pivotal role in wave characteristics in the weakly ionized low density regions, hence maximizing the local ambipolar diffusion coefficient. As a result, we propose that ambipolar diffusion is an important mechanism that requires careful consideration into whether it should be included in simulations, and whether it should be utilized in the analysis and interpretation of particular observations of the lower solar atmosphere. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 13 pages, 10 figures. Accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal, 938(2), 154, 2022

Observational and numerical characterization of a recurrent arc-shaped front propagating along a coronal fan

Authors: M. V. Sieyra, S. Krishna Prasad, G. Stenborg, E. Khomenko, T. Van Doorsselaere, A. Costa, A. Esquivel, J. M. Riedl

Abstract: Recurrent, arc-shaped intensity disturbances were detected by EUV channels in an active region. The fronts were observed to propagate along a coronal loop bundle rooted in a small area within a sunspot umbra. Previous works have linked these intensity disturbances to slow magnetoacoustic waves that propagate from the lower atmosphere to the corona along the magnetic field. The slow magnetoacoustic… ▽ More Recurrent, arc-shaped intensity disturbances were detected by EUV channels in an active region. The fronts were observed to propagate along a coronal loop bundle rooted in a small area within a sunspot umbra. Previous works have linked these intensity disturbances to slow magnetoacoustic waves that propagate from the lower atmosphere to the corona along the magnetic field. The slow magnetoacoustic waves propagate at the local cusp speed. However, the measured propagation speeds from the intensity images are usually smaller as they are subject to projection effects due to the inclination of the magnetic field with respect to the line-of-sight. Here, we aim to understand the effect of projection by comparing observed speeds with those from a numerical model. Using multi-wavelength data we determine the periods present in the observations at different heights of the solar atmosphere through Fourier analysis. We calculate the plane-of-sky speeds along one of the loops from the cross-correlation time lags obtained as a function of distance along the loop. We perform a 2D ideal MHD simulation of an active region embedded in a stratified atmosphere. We drive slow waves from the photosphere with a 3 minutes periodicity. Synthetic time-distance maps are generated from the forward-modelled intensities in coronal wavelengths and the projected propagation speeds are calculated. The intensity disturbances show a dominant period between [2-3] minutes at different heights of the atmosphere. The apparent propagation speeds calculated for coronal channels exhibit an accelerated pattern with values increasing from 40 to 120 km/s as the distance along the loop rises. The propagation speeds obtained from the synthetic time-distance maps also exhibit accelerated profiles within a similar range of speeds. We conclude that the accelerated propagation in our observations is due to the projection effect. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 12 pages, 12 figures

Journal ref: A&A 667, A21 (2022)

The European Solar Telescope

Authors: C. Quintero Noda, R. Schlichenmaier, L. R. Bellot Rubio, M. G. Löfdahl, E. Khomenko, J. Jurcak, J. Leenaarts, C. Kuckein, S. J. González Manrique, S. Gunar, C. J. Nelson, J. de la Cruz Rodríguez, K. Tziotziou, G. Tsiropoula, G. Aulanier, M. Collados, the EST team

Abstract: The European Solar Telescope (EST) is a project aimed at studying the magnetic connectivity of the solar atmosphere, from the deep photosphere to the upper chromosphere. Its design combines the knowledge and expertise gathered by the European solar physics community during the construction and operation of state-of-the-art solar telescopes operating in visible and near-infrared wavelengths: the Sw… ▽ More The European Solar Telescope (EST) is a project aimed at studying the magnetic connectivity of the solar atmosphere, from the deep photosphere to the upper chromosphere. Its design combines the knowledge and expertise gathered by the European solar physics community during the construction and operation of state-of-the-art solar telescopes operating in visible and near-infrared wavelengths: the Swedish 1m Solar Telescope (SST), the German Vacuum Tower Telescope (VTT) and GREGOR, the French Télescope Héliographique pour l'Étude du Magnétisme et des Instabilités Solaires (THÉMIS), and the Dutch Open Telescope (DOT). With its 4.2 m primary mirror and an open configuration, EST will become the most powerful European ground-based facility to study the Sun in the coming decades in the visible and near-infrared bands. EST uses the most innovative technological advances: the first adaptive secondary mirror ever used in a solar telescope, a complex multi-conjugate adaptive optics with deformable mirrors that form part of the optical design in a natural way, a polarimetrically compensated telescope design that eliminates the complex temporal variation and wavelength dependence of the telescope Mueller matrix, and an instrument suite containing several (etalon-based) tunable imaging spectropolarimeters and several integral field unit spectropolarimeters. This publication summarises some fundamental science questions that can be addressed with the telescope, together with a complete description of its major subsystems. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: 38 pages, 24 figures

Journal ref: A&A 666, A21 (2022)

Modeling the thermal conduction in the solar atmosphere with the code MANCHA3D

Authors: Anamaría Navarro, E. Khomenko, M. Modestov, N. Vitas

Abstract: Thermal conductivity is one of the important mechanisms of heat transfer in the solar corona. In the limit of strongly magnetized plasma, it is typically modeled by Spitzer's expression where the heat flux is aligned with the magnetic field. This paper describes the implementation of the heat conduction into the code MANCHA3D with an aim of extending single-fluid MHD simulations from the upper con… ▽ More Thermal conductivity is one of the important mechanisms of heat transfer in the solar corona. In the limit of strongly magnetized plasma, it is typically modeled by Spitzer's expression where the heat flux is aligned with the magnetic field. This paper describes the implementation of the heat conduction into the code MANCHA3D with an aim of extending single-fluid MHD simulations from the upper convection zone into the solar corona. Two different schemes to model heat conduction are implemented: (1) a standard scheme where a parabolic term is added to the energy equation, and (2) a scheme where the hyperbolic heat flux equation is solved. The first scheme limits the time step due to the explicit integration of a parabolic term, which makes the simulations computationally expensive. The second scheme solves the limitations on the time step by artificially limiting the heat conduction speed to computationally manageable values. The validation of both schemes is carried out with standard tests in one, two, and three spatial dimensions. Furthermore, we implement the model for heat flux derived by Braginskii (1965) in its most general form, when the expression for the heat flux depends on the ratio of the collisional to cyclotron frequencies of the plasma, and, therefore on the magnetic field strength. Additionally, our implementation takes into account the heat conduction in parallel, perpendicular, and transverse directions, and provides the contributions from ions and electrons separately. The model also transitions smoothly between field-aligned conductivity and isotropic conductivity for regions with a low or null magnetic field. Finally, we present a two-dimensional test for heat conduction using realistic values of the solar atmosphere where we prove the robustness of the two schemes implemented. △ Less

Submitted 20 May, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: 11 pages, 8 figures

Journal ref: A&A 663, A96 (2022)

Generalized Fluid Models of the Braginskii Type

Authors: P. Hunana, T. Passot, E. Khomenko, D. Martinez-Gomez, M. Collados, A. Tenerani, G. P. Zank, Y. Maneva, M. L. Goldstein, G. M. Webb

Abstract: Several generalizations of the well-known fluid model of Braginskii (Rev. of Plasma Phys., 1965) are considered. We use the Landau collisional operator and the moment method of Grad. We focus on the 21-moment model that is analogous to the Braginskii model, and we also consider a 22-moment model. Both models are formulated for general multi-species plasmas with arbitrary masses and temperatures, w… ▽ More Several generalizations of the well-known fluid model of Braginskii (Rev. of Plasma Phys., 1965) are considered. We use the Landau collisional operator and the moment method of Grad. We focus on the 21-moment model that is analogous to the Braginskii model, and we also consider a 22-moment model. Both models are formulated for general multi-species plasmas with arbitrary masses and temperatures, where all the fluid moments are described by their evolution equations. The 21-moment model contains two "heat flux vectors" (3rd and 5th-order moments) and two "viscosity-tensors" (2nd and 4th-order moments). The Braginskii model is then obtained as a particular case of a one ion-electron plasma with similar temperatures, with de-coupled heat fluxes and viscosity-tensors expressed in a quasi-static approximation. We provide all the numerical values of the Braginskii model in a fully analytic form (together with the 4th and 5th-order moments). For multi-species plasmas, the model makes calculation of transport coefficients straightforward. Formulation in fluid moments (instead of Hermite moments) is also suitable for implementation into existing numerical codes. It is emphasized that it is the quasi-static approximation which makes some Braginskii coefficients divergent in a weakly-collisional regime. Importantly, we show that the heat fluxes and viscosity-tensors are coupled even in the linear approximation, and that the fully contracted (scalar) perturbations of the 4th-order moment, which are accounted for in the 22-moment model, modify the energy exchange rates. We also provide several Appendices, which can be useful as a guide for deriving the Braginskii model with the moment method of Grad. △ Less

Submitted 2 June, 2022; v1 submitted 27 January, 2022; originally announced January 2022.

Comments: Updated version. Title change from "Braginskii-type" to "Braginskii Type" to match the published version (accepted to ApJ Supplements, 29 Jan 2022). Bookmarks included, TOC is now at the beginning, numbering in Appendix I is fixed, few typos corrected

Magnetic field amplification and structure formation by the Rayleigh-Taylor instability

Authors: B. Popescu Braileanu, V. S. Lukin, E. Khomenko

Abstract: We report on results of high resolution two fluid non-linear simulations of the Rayleigh Taylor Instability (RTI) at the interface between a solar prominence and the corona. These follow results reported earlier by Popescu Braileanu et al. (2021a,b) on linear and early non-linear RTI dynamics in this environment. The simulations use a two fluid model that includes collisions between neutrals and c… ▽ More We report on results of high resolution two fluid non-linear simulations of the Rayleigh Taylor Instability (RTI) at the interface between a solar prominence and the corona. These follow results reported earlier by Popescu Braileanu et al. (2021a,b) on linear and early non-linear RTI dynamics in this environment. The simulations use a two fluid model that includes collisions between neutrals and charges, including ionization/recombination, energy and momentum transfer, and frictional heating. High resolution 2.5D magnetized RTI simulations with the magnetic field dominantly normal to and slightly sheared with respect to the prominence plane demonstrate that in a fully developed state of RTI a large fraction of the gravitational energy of a prominence thread can be converted into quasi-turbulent energy of the magnetic field. RTI magnetic energy generation is further accompanied by magnetic and plasma density structure formation, including dynamic formation, break-up, and merging of current sheets and plasmoid sub-structures. The simulations show the role of flow decoupling and ionization/recombination reactions between the neutrals and charges on the structure formation in magnetized RTI. We provide a careful examination of sources and form of numerical dissipation of the evolving magnetic field structures. △ Less

Submitted 22 December, 2022; v1 submitted 24 December, 2021; originally announced December 2021.

Journal ref: A&A 670, A31 (2023)

Transverse kink oscillations of inhomogeneous prominence threads: numerical analysis and H$α$ forward modelling

Authors: David Martínez-Gómez, Roberto Soler, Jaume Terradas, Elena Khomenko

Abstract: Prominence threads are very long and thin flux tubes which are partially filled with cold plasma. Observations have shown that transverse oscillations are frequent in these solar structures. The observations are usually interpreted as the fundamental kink mode, while the detection of the first harmonic remains elusive. Here, we aim to study how the density inhomogeneity in the longitudinal and rad… ▽ More Prominence threads are very long and thin flux tubes which are partially filled with cold plasma. Observations have shown that transverse oscillations are frequent in these solar structures. The observations are usually interpreted as the fundamental kink mode, while the detection of the first harmonic remains elusive. Here, we aim to study how the density inhomogeneity in the longitudinal and radial directions modify the periods and damping times of kink oscillations, and how this effect would be reflected in observations. We solve the ideal magnetohydrodynamics equations through two different methods: a) performing 3D numerical simulations, and b) solving a 2D generalised eigenvalue problem. We study the dependence of the periods, damping times and amplitudes of transverse kink oscillations on the ratio between the densities at the centre and at the ends of the tube, and on the average density. We apply forward modelling on our 3D simulations to compute synthetic H$α$ profiles. We confirm that the ratio of the period of the fundamental oscillation mode to the period of the first harmonic increases as the ratio of the central density to the footpoint density is increased or as the averaged density of the tube is decreased. We find that the damping times due to resonant absorption decrease as the central to footpoint density ratio increases. Contrary to the case of longitudinally homogeneous tubes, we find that the damping time to period ratio also increases as the density ratio is increased or the average density is reduced. We present snapshots and time-distance diagrams of the emission in the H$α$ line. The results presented here have implications for the field of prominence seismology. While the H$α$ emission can be used to detect the fundamental mode, the first harmonic is barely detectable in H$α$. This may explain the lack of detections of the first harmonic. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: 17 pages, 16 figures, 4 movies. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 658, A106 (2022)

Large-amplitude longitudinal oscillations in solar prominences simulated with different resolutions

Authors: Valeriia Liakh, Manuel Luna, Elena Khomenko

Abstract: Large-amplitude longitudinal oscillations (LALOs) in solar prominences have been widely studied in the last decades. However, their damping and amplification mechanisms are not well understood. In this study, we investigate the attenuation and amplification of LALOs using high-resolution numerical simulations with progressively increasing spatial resolutions. We performed time-dependent numerical… ▽ More Large-amplitude longitudinal oscillations (LALOs) in solar prominences have been widely studied in the last decades. However, their damping and amplification mechanisms are not well understood. In this study, we investigate the attenuation and amplification of LALOs using high-resolution numerical simulations with progressively increasing spatial resolutions. We performed time-dependent numerical simulations of LALOs using the 2D magnetic configuration that contains a dipped region. After the prominence mass loading in the magnetic dips, we triggered LALOs by perturbing the prominence mass along the magnetic field. We performed the experiments with four values of spatial resolution. In the simulations with the highest resolution, the period shows a good agreement with the pendulum model. The convergence experiment revealed that the damping time saturates at the bottom prominence region with improving the resolution, indicating the existence of a physical reason for the damping of oscillations. At the prominence top, the oscillations are amplified during the first minutes and then are slowly attenuated. The characteristic time suggests more significant amplification in the experiments with the highest spatial resolution. The analysis revealed that the energy exchange between the bottom and top prominence regions is responsible for the attenuation and amplification of LALOs. The high-resolution experiments are crucial for the study of the periods and the damping mechanism of LALOs. The period agrees with the pendulum model only when using high enough spatial resolution. The results suggest that numerical diffusion in simulations with insufficient spatial resolution can hide important physical mechanisms, such as amplification of oscillations. △ Less

Submitted 2 August, 2021; originally announced August 2021.

Comments: 12 pages, 13 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 654, A145 (2021)

Simulations of the Biermann battery mechanism in two-fluid partially ionised plasmas

Authors: David Martínez-Gómez, Beatrice Popescu Braileanu, Elena Khomenko, Peter Hunana

Abstract: In the absence of an initial seed, the Biermann battery term of a non-ideal induction equation acts as a source that generates weak magnetic fields. These fields are then amplified via a dynamo mechanism. The Kelvin-Helmholtz instability is a fluid phenomenon that takes place in many astrophysical scenarios and can trigger the action of the Biermann battery and dynamo processes. We aim to investig… ▽ More In the absence of an initial seed, the Biermann battery term of a non-ideal induction equation acts as a source that generates weak magnetic fields. These fields are then amplified via a dynamo mechanism. The Kelvin-Helmholtz instability is a fluid phenomenon that takes place in many astrophysical scenarios and can trigger the action of the Biermann battery and dynamo processes. We aim to investigate the effect that the ionisation degree of the plasma and the interaction between the charged and neutral species have on the generation and amplification of magnetic fields during the different stages of the instability. We use the two-fluid model implemented in the numerical code Mancha-2F. We perform 2D simulations starting from a configuration with no initial magnetic field and which is unstable due to a velocity shear. We vary the ionisation degree of the plasma and we analyse the role that the different collisional terms included in the equations of the model play on the evolution of the instability and the generation of magnetic field. We find that when no collisional coupling is considered between the two fluids, the effect of the Biermann battery mechanism does not depend on the ionisation degree. However, when elastic collisions are taken into account, the generation of magnetic field is increased as the ionisation degree is reduced. This behaviour is slightly enhanced if the process of charge-exchange is also considered. We also find a dependence on the total density of the plasma related to the dependence on the coupling degree between the two fluids. As the total density is increased, the results from the two-fluid model converge to the predictions of single-fluid models. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: 20 pages, 13 figures, 2 movies. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 650, A123 (2021)

Chromospheric Heating by MHD Waves and Instabilities

Authors: A. K. Srivastava, J. L. Ballester, P. S. Cally, M. Carlsson, M. Goossens, D. B. Jess, E. Khomenko, M. Mathioudakis, K. Murawski, T. V. Zaqarashvili

Abstract: The importance of the chromosphere in the mass and energy transport within the solar atmosphere is now widely recognised. This review discusses the physics of magnetohydrodynamic (MHD) waves and instabilities in large-scale chromospheric structures as well as in magnetic flux tubes. We highlight a number of key observational aspects that have helped our understanding of the role of the solar chrom… ▽ More The importance of the chromosphere in the mass and energy transport within the solar atmosphere is now widely recognised. This review discusses the physics of magnetohydrodynamic (MHD) waves and instabilities in large-scale chromospheric structures as well as in magnetic flux tubes. We highlight a number of key observational aspects that have helped our understanding of the role of the solar chromosphere in various dynamic processes and wave phenomena, and the heating scenario of the solar chromosphere is also discussed. The review focuses on the physics of waves and invokes the basics of plasma instabilities in the context of this important layer of the solar atmosphere. Potential implications, future trends and outstanding questions are also delineated. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: This is an invited review article in the special issue "Solar and Heliospheric Plasma Structures: Waves, Turbulence, and Dissipation" of JGR - Space Physics (57 pages, 9 figures, Paper #2020JA029097R) with Editors: Leon Ofman and Shreekrishna Tripathi

Two-fluid simulations of Rayleigh-Taylor instability in a magnetized solar prominence thread II. Effects of collisionality

Authors: B. Popescu Braileanu, V. S. Lukin, E. Khomenko, A. de Vicente

Abstract: In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid m… ▽ More In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid model takes into account neutral viscosity, thermal conductivity, and collisional interaction between neutrals and charges: ionization/recombination, energy and momentum transfer, and frictional heating. In this paper II, the sensitivity of the RTI dynamics to collisional effects for different magnetic field configurations supporting the prominence thread is explored. This is done by artificially varying, or eliminating, effects of both elastic and inelastic collisions by modifying the model equations. We find that ionization and recombination reactions between ionized and neutral fluids, if in equilibrium prior to the onset of the instability, do not substantially impact the development of the primary RTI. However, such reactions can impact development of secondary structures during mixing of the cold prominence and hotter surrounding coronal material. We find that collisionality within and between ionized and neutral particle populations play an important role in both linear and non-linear development of RTI, with ion-neutral collision frequency as the primary determining factor in development or damping of small scale structures. We also observe that degree and signatures of flow decoupling between ion and neutral fluids can depend both on the inter-particle collisionality and the magnetic field configuration of the prominence thread. △ Less

Submitted 8 April, 2021; v1 submitted 29 January, 2021; originally announced January 2021.

Journal ref: A&A 650, A181 (2021)

Influence of ambipolar and Hall effects on vorticity in 3D simulations of magneto-convection

Authors: E. Khomenko, M. Collados, N. Vitas, P. A. Gonzalez-Morales

Abstract: This paper presents the results of the analysis of 3D simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation-runs are compared: one including both ambipolar diffusion and Hall effect; one including only ambipolar diffusion; and one without any of these two effects. The magnetic field is amplified from initial field to… ▽ More This paper presents the results of the analysis of 3D simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation-runs are compared: one including both ambipolar diffusion and Hall effect; one including only ambipolar diffusion; and one without any of these two effects. The magnetic field is amplified from initial field to saturation level by the action of turbulent local dynamo. In each of these cases, we study 2 hours of simulated solar time after the local dynamo reaches the saturation regime. We analyze the power spectra of vorticity, of magnetic field fluctuations and of the different components of the magnetic Poynting flux responsible for the transport of vertical or horizontal perturbations. Our preliminary results show that the ambipolar diffusion produces a strong reduction of vorticity in the upper chromospheric layers and that it dissipates the vortical perturbations converting them into thermal energy. The Hall effect acts in the opposite way, strongly enhancing the vorticity. When the Hall effect is included, the magnetic field in the simulations becomes, on average, more vertical and long-lived flux tube-like structures are produced. We trace a single magnetic structure to study its evolution pattern and the magnetic field intensification, and their possible relation to the Hall effect. △ Less

Submitted 21 September, 2020; originally announced September 2020.

Comments: to appear in Philosophical Transactions of the Royal Society A, issue "High resolution wave dynamics in the lower solar atmosphere"

Joint action of Hall and ambipolar effects in 3D magneto-convection simulations of the quiet Sun. I. Dissipation and generation of waves

Authors: P. A. Gonzalez-Morales, E. Khomenko, N. Vitas, M. Collados

Abstract: The partial ionization of the solar plasma causes several nonideal effects such as the ambipolar diffusion, the Hall effect, and the Biermann battery effect. Here we report on the first three-dimensional realistic simulations of solar local dynamo where all three effects were taken into account. The simulations started with a snapshot of already saturated battery-seeded dynamo, where two new serie… ▽ More The partial ionization of the solar plasma causes several nonideal effects such as the ambipolar diffusion, the Hall effect, and the Biermann battery effect. Here we report on the first three-dimensional realistic simulations of solar local dynamo where all three effects were taken into account. The simulations started with a snapshot of already saturated battery-seeded dynamo, where two new series were developed: one with solely ambipolar diffusion and another one also taking into account the Hall term in the generalized Ohm's law. The simulations were then run for about 4 hours of solar time to reach the stationary regime and improve the statistics. In parallel, a purely MHD dynamo simulation was also run for the same amount of time. The simulations are compared in a statistical way. The results show that, with the inclusion of the ambipolar diffusion, the amplitudes of the incompressible perturbations related to Alfven waves are reduced, and the Poynting flux is absorbed, with a frequency dependence. The Hall effect causes the opposite action: significant excess of incompressible perturbations is generated and an excess of the Poynting flux is observed in the chromospheric layers. The model with ambipolar diffusion shows, on average, sharper current sheets and slightly more abundant fast magneto-acoustic shocks in the chromosphere. The model with the Hall effect has higher temperatures at the lower chromosphere and stronger and more vertical magnetic field concentrations all over the chromosphere. The study of high-frequency waves reveals that significant power of incompressible perturbations is associated with areas with intense and more vertical magnetic fields and larger temperatures. We find a positive correlation between the magnitude of the ambipolar heating and the temperature increase at the same location after a characteristic time of 10^2 sec. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Comments: To appear in Astronomy & Astrophysics

Critical Science Plan for the Daniel K. Inouye Solar Telescope (DKIST)

Authors: Mark P. Rast, Nazaret Bello González, Luis Bellot Rubio, Wenda Cao, Gianna Cauzzi, Edward DeLuca, Bart De Pontieu, Lyndsay Fletcher, Sarah E. Gibson, Philip G. Judge, Yukio Katsukawa, Maria D. Kazachenko, Elena Khomenko, Enrico Landi, Valentin Martínez Pillet, Gordon J. D. Petrie, Jiong Qiu, Laurel A. Rachmeler, Matthias Rempel, Wolfgang Schmidt, Eamon Scullion, Xudong Sun, Brian T. Welsch, Vincenzo Andretta, Patrick Antolin , et al. (62 additional authors not shown)

Abstract: The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With… ▽ More The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With this Critical Science Plan (CSP) we attempt to anticipate some of what those capabilities will enable, providing a snapshot of some of the scientific pursuits that the Daniel K. Inouye Solar Telescope hopes to engage as start-of-operations nears. The work builds on the combined contributions of the DKIST Science Working Group (SWG) and CSP Community members, who generously shared their experiences, plans, knowledge and dreams. Discussion is primarily focused on those issues to which DKIST will uniquely contribute. △ Less

Submitted 20 August, 2020; v1 submitted 18 August, 2020; originally announced August 2020.

Two-fluid simulations of Rayleigh-Taylor instability in a magnetized solar prominence thread. I. Effects of prominence magnetization and mass loading

Authors: B. Popescu Braileanu, V. S. Lukin, E. Khomenko, A. de Vicente

Abstract: Solar prominences are formed by partially ionized plasma with inter-particle collision frequencies generally warranting magnetohydrodynamic treatment. In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. We perform 2.5D two-fluid (charges - neutrals) sim… ▽ More Solar prominences are formed by partially ionized plasma with inter-particle collision frequencies generally warranting magnetohydrodynamic treatment. In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. We perform 2.5D two-fluid (charges - neutrals) simulations of the Rayleigh-Taylor instability (RTI) at a smoothly changing interface between a solar prominence thread and the corona. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid model takes into account viscosity, thermal conductivity, and collisional interaction between neutrals and charges: ionization/recombination, energy and momentum transfer, and frictional heating. In this paper I, the sensitivity of the RTI dynamics to the prominence equilibrium configuration, including the impact of the magnetic field strength and shear supporting the prominence thread, and the amount of prominence mass-loading is explored. We show that, at small scales, a realistically smooth prominence-corona interface leads to qualitatively different linear RTI evolution than that expected for a discontinuous interface, while magnetic field shear has the stabilizing effect of reducing the growth rate or eliminating the instability. In the non-linear phase, we observe that in the presence of field shear the development of the instability leads to formation of coherent and interacting 2.5D magnetic structures, which, in turn, can lead to substantial plasma flow across magnetic field lines and associated decoupling of the fluid velocities of charges and neutrals. △ Less

Submitted 15 December, 2020; v1 submitted 31 July, 2020; originally announced July 2020.

Journal ref: A&A 646, A93 (2021)

Accurately constraining velocity information from spectral imaging observations using machine learning techniques

Authors: Conor D. MacBride, David B. Jess, Samuel D. T. Grant, Elena Khomenko, Peter H. Keys, Marco Stangalini

Abstract: Determining accurate plasma Doppler (line-of-sight) velocities from spectroscopic measurements is a challenging endeavour, especially when weak chromospheric absorption lines are often rapidly evolving and, hence, contain multiple spectral components in their constituent line profiles. Here, we present a novel method that employs machine learning techniques to identify the underlying components pr… ▽ More Determining accurate plasma Doppler (line-of-sight) velocities from spectroscopic measurements is a challenging endeavour, especially when weak chromospheric absorption lines are often rapidly evolving and, hence, contain multiple spectral components in their constituent line profiles. Here, we present a novel method that employs machine learning techniques to identify the underlying components present within observed spectral lines, before subsequently constraining the constituent profiles through single or multiple Voigt fits. Our method allows active and quiescent components present in spectra to be identified and isolated for subsequent study. Lastly, we employ a Ca II 8542 Å spectral imaging dataset as a proof-of-concept study to benchmark the suitability of our code for extracting two-component atmospheric profiles that are commonly present in sunspot chromospheres. Minimisation tests are employed to validate the reliability of the results, achieving median reduced $χ^2$ values equal to 1.03 between the observed and synthesised umbral line profiles. △ Less

Submitted 4 January, 2021; v1 submitted 15 July, 2020; originally announced July 2020.

Comments: 23 pages, 8 figures. Improved formatting of abstract and references

Journal ref: Philosophical Transactions of the Royal Society A 379, 20200171 (2020)

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

Authors: B. Fleck, M. Carlsson, E. Khomenko, M. Rempel, O. Steiner, G. Vigeesh

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

Submitted 11 July, 2020; originally announced July 2020.

Comments: 19 pages, 15 figures

Numerical simulations of large-amplitude oscillations in flux-rope solar prominences

Authors: Valeriia Liakh, Manuel Luna, Elena Khomenko

Abstract: Context. Large-amplitude oscillations (LAOs) of the solar prominences are very spectacular but poorly understood phenomena. These motions have amplitudes larger than $10 {\mathrm{\, km \,s^{-1}}}$ and can be triggered by the external perturbations, e.g., Moreton or EIT waves. Aims. Our aim is to analyze the properties of large-amplitude oscillations using realistic prominence models and the trig… ▽ More Context. Large-amplitude oscillations (LAOs) of the solar prominences are very spectacular but poorly understood phenomena. These motions have amplitudes larger than $10 {\mathrm{\, km \,s^{-1}}}$ and can be triggered by the external perturbations, e.g., Moreton or EIT waves. Aims. Our aim is to analyze the properties of large-amplitude oscillations using realistic prominence models and the triggering mechanism by external disturbances. Methods. We perform time-dependent numerical simulations of LAOs using a magnetic flux rope model with two values of the shear angle and the density contrast. We study the internal modes of the prominence using the horizontal and vertical triggering. In addition, we use the perturbation that arrives from outside in order to understand how such external disturbance can produce LAOs. Results. The period of longitudinal oscillations and its behavior with height show good agreement with the pendulum model. The period of transverse oscillations remains constant with height, suggesting a global normal mode. The transverse oscillations typically have shorter periods than the longitudinal oscillations. Conclusions. The periods of the longitudinal and transverse oscillations show only weak dependence on the shear angle of the magnetic structure and the prominence density contrast. The external disturbance perturbs the flux rope exciting oscillations of both polarizations. Their properties are a mixture of those excited by purely horizontal and vertical driving. △ Less

Submitted 9 March, 2020; originally announced March 2020.

Comments: 15 pages, 10 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 637, A75 (2020)

Science Requirement Document (SRD) for the European Solar Telescope (EST) (2nd edition, December 2019)

Authors: R. Schlichenmaier, L. R. Bellot Rubio, M. Collados, R. Erdelyi, A. Feller, L. Fletcher, J. Jurcak, E. Khomenko, J. Leenaarts, S. Matthews, L. Belluzzi, M. Carlsson, K. Dalmasse, S. Danilovic, P. Gömöry, C. Kuckein, R. Manso Sainz, M. Martinez Gonzalez, M. Mathioudakis, A. Ortiz, T. L. Riethmüller, L. Rouppe van der Voort, P. J. A. Simoes, J. Trujillo Bueno, D. Utz , et al. (1 additional authors not shown)

Abstract: The European Solar Telescope (EST) is a research infrastructure for solar physics. It is planned to be an on-axis solar telescope with an aperture of 4 m and equipped with an innovative suite of spectro-polarimetric and imaging post-focus instrumentation. The EST project was initiated and is driven by EAST, the European Association for Solar Telescopes. EAST was founded in 2006 as an association o… ▽ More The European Solar Telescope (EST) is a research infrastructure for solar physics. It is planned to be an on-axis solar telescope with an aperture of 4 m and equipped with an innovative suite of spectro-polarimetric and imaging post-focus instrumentation. The EST project was initiated and is driven by EAST, the European Association for Solar Telescopes. EAST was founded in 2006 as an association of 14 European countries. Today, as of December 2019, EAST consists of 26 European research institutes from 18 European countries. The Preliminary Design Phase of EST was accomplished between 2008 and 2011. During this phase, in 2010, the first version of the EST Science Requirement Document (SRD) was published. After EST became a project on the ESFRI roadmap 2016, the preparatory phase started. The goal of the preparatory phase is to accomplish a final design for the telescope and the legal governance structure of EST. A major milestone on this path is to revisit and update the Science Requirement Document (SRD). The EST Science Advisory Group (SAG) has been constituted by EAST and the Board of the PRE-EST EU project in November 2017 and has been charged with the task of providing with a final statement on the science requirements for EST. Based on the conceptual design, the SRD update takes into account recent technical and scientific developments, to ensure that EST provides significant advancement beyond the current state-of-the-art. The present update of the EST SRD has been developed and discussed during a series of EST SAG meetings. The SRD develops the top-level science objectives of EST into individual science cases. Identifying critical science requirements is one of its main goals. Those requirements will define the capabilities of EST and the post-focus instrument suite. The technical requirements for the final design of EST will be derived from the SRD. △ Less

Submitted 18 December, 2019; originally announced December 2019.

Comments: 2nd edition, December 2019, 138 pages

Two-dimensional simulations of coronal rain dynamics. I. Model with vertical magnetic field and an unbounded atmosphere

Authors: David Martínez-Gómez, Ramón Oliver, Elena Khomenko, Manuel Collados

Abstract: Aims. We aim to improve the understanding of the physical mechanisms behind the slower than free-fall motion and the two-stage evolution (an initial phase of acceleration followed by an almost constant velocity phase) detected in coronal rain events. Methods. Using the Mancha3D code, we solve the 2D ideal MHD equations. We represent the solar corona as an isothermal vertically stratified atmosph… ▽ More Aims. We aim to improve the understanding of the physical mechanisms behind the slower than free-fall motion and the two-stage evolution (an initial phase of acceleration followed by an almost constant velocity phase) detected in coronal rain events. Methods. Using the Mancha3D code, we solve the 2D ideal MHD equations. We represent the solar corona as an isothermal vertically stratified atmosphere with a uniform vertical magnetic field and the plasma condensation as a density enhancement described by a 2D Gaussian profile. We analyse the temporal evolution of the descending plasma and study its dependence on parameters such as density and magnetic field strength. Results. We confirm previous findings that the pressure gradient is the main force that opposes the action of gravity and slows down the blob descent and that larger densities require larger pressure gradients to reach the constant speed phase. We find that the shape of a condensation with a horizontal variation of density is distorted as it falls, due to the denser parts of the blob falling faster than the lighter ones. This is explained by the fact that the duration of the initial acceleration phase, and therefore the maximum falling speed attained by the plasma, increases with the ratio of blob to coronal density. We also find that the magnetic field plays a fundamental role in the evolution of the descending condensations. A strong enough magnetic field (greater than 10 G in our simulations) forces each plasma element to follow the path given by a particular field line, which allows to describe the evolution of each vertical slice of the blob in terms of 1D dynamics, without influence of the adjacent slices. In addition, under the typical conditions of the coronal rain events, the magnetic field prevents the development of the Kelvin-Helmholtz instability. △ Less

Submitted 15 November, 2019; originally announced November 2019.

Comments: 14 pages, 11 figures; submitted to Astronomy & Astrophysics

Journal ref: A&A 634, A36 (2020)

Two-fluid simulations of waves in the solar chromosphere II. Propagation and damping of fast magneto-acoustic waves and shocks

Authors: B. Popescu Braileanu, V. S. Lukin, E. Khomenko, A. de Vicente

Abstract: Waves and shocks traveling through the solar chromospheric plasma are influenced by its partial ionization and weak collisional coupling, and may become susceptible to multi-fluid effects, similar to interstellar shock waves. In this study, we consider fast magneto-acoustic shock wave formation and propagation in a stratified medium, that is permeated by a horizontal magnetic field, with propertie… ▽ More Waves and shocks traveling through the solar chromospheric plasma are influenced by its partial ionization and weak collisional coupling, and may become susceptible to multi-fluid effects, similar to interstellar shock waves. In this study, we consider fast magneto-acoustic shock wave formation and propagation in a stratified medium, that is permeated by a horizontal magnetic field, with properties similar to that of the solar chromosphere. The evolution of plasma and neutrals is modeled using a two-fluid code that evolves a set of coupled equations for two separate fluids. We observed that waves in neutrals and plasma, initially coupled at the upper photosphere, become uncoupled at higher heights in the chromosphere. This decoupling can be a consequence of either the characteristic spatial scale at the shock front, that becomes similar to the collisional scale, or the change in the relation between the wave frequency, ion cyclotron frequency, and the collisional frequency with height. The decoupling height is a sensitive function of the wave frequency, wave amplitude, and the magnetic field strength. We observed that decoupling causes damping of waves and an increase in the background temperature due to the frictional heating. The comparison between analytical and numerical results allows us to separate the role of the nonlinear effects from the linear ones on the decoupling and damping of waves. △ Less

Submitted 14 August, 2019; originally announced August 2019.

Journal ref: A&A 630, A79 (2019)

Two-fluid simulations of waves in the solar chromosphere I: numerical code verification

Authors: B. Popescu Braileanu, V. S. Lukin, E. Khomenko, A. de Vicente

Abstract: Solar chromosphere consists of a partially ionized plasma, which makes modeling the solar chromosphere a particularly challenging numerical task. Here we numerically model chromospheric waves using a two-fluid approach with a newly developed numerical code. The code solves two-fluid equations of conservation of mass, momentum and energy, together with the induction equation, for the case of the pu… ▽ More Solar chromosphere consists of a partially ionized plasma, which makes modeling the solar chromosphere a particularly challenging numerical task. Here we numerically model chromospheric waves using a two-fluid approach with a newly developed numerical code. The code solves two-fluid equations of conservation of mass, momentum and energy, together with the induction equation, for the case of the purely hydrogen plasma with collisional coupling between the charged and neutral fluid components. The implementation of a semi-implicit algorithm allows us to overcome the numerical stability constraints due to the stiff collisional terms. We test the code against analytical solutions of acoustic and Alfvén wave propagation in uniform medium in several regimes of collisional coupling.The results of our simulations are consistent with the analytical estimates, and with other results described in the literature. In the limit of a large collisional frequency, the waves propagate with a common speed of a single fluid. In the other limit of a vanishingly small collisional frequency, the Alfvén waves propagate with an Alfvén speed of the charged fluid only, while the perturbation in neutral fluid is very small. The acoustic waves in these limits propagate with the sound speed corresponding to either the charges or the neutrals, while the perturbation in the other fluid component is very small. Otherwise, when the collision frequency is similar to the real part of the wave frequency, the interaction between charges and neutrals through momentum transfer collisions cause alterations of the waves frequencies and damping of the wave amplitudes. △ Less

Submitted 9 May, 2019; originally announced May 2019.

Journal ref: A&A 627, A25 (2019)

An introductory guide to fluid models with anisotropic temperatures Part 2 -- Kinetic theory, Padé approximants and Landau fluid closures

Authors: P. Hunana, A. Tenerani, G. P. Zank, M. L. Goldstein, G. M. Webb, E. Khomenko, M. Collados, P. S. Cally, L. Adhikari, M. Velli

Abstract: In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do… ▽ More In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do not converge to the kinetic description, even in the long-wavelength low-frequency limit. We start with a brief introduction to kinetic theory, where we discuss in detail the plasma dispersion function $Z(ζ)$, and the associated plasma response function $R(ζ)=1+ζZ(ζ)=-Z'(ζ)/2$. We then consider a 1D (electrostatic) geometry and make a significant effort to map all possible Landau fluid closures that can be constructed at the 4th-order moment level. These closures for parallel moments have general validity from the largest astrophysical scales down to the Debye length, and we verify their validity by considering examples of the (proton and electron) Landau damping of the ion-acoustic mode, and the electron Landau damping of the Langmuir mode. We proceed by considering 1D closures at higher-order moments than the 4th-order, and as was concluded in Part 1, this is not possible without Landau fluid closures. We show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing the convergence of the fluid and collisionless kinetic descriptions. We then consider a 3D (electromagnetic) geometry in the gyrotropic (long-wavelength low-frequency) limit and map all closures that are available at the 4th-order moment level. In the Appendix A, we provide comprehensive tables with Padé approximants of $R(ζ)$ up to the 8th-pole order, with many given in an analytic form. △ Less

Submitted 7 November, 2019; v1 submitted 27 January, 2019; originally announced January 2019.

Comments: Improved version, now accepted to JPP Lecture Notes. Some parts were shortened and some parts were expanded. The text now contains Conclusions

An introductory guide to fluid models with anisotropic temperatures Part 1 -- CGL description and collisionless fluid hierarchy

Authors: P. Hunana, A. Tenerani, G. P. Zank, E. Khomenko, M. L. Goldstein, G. M. Webb, P. S. Cally, M. Collados, M. Velli, L. Adhikari

Abstract: We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modeling turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the inte… ▽ More We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modeling turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the CGL fluid model in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius (FLR) corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the 4th-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the "normal" closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the "normal" closure, and beyond the 4th-order moment, Landau fluid closures are required. △ Less

Submitted 27 October, 2019; v1 submitted 27 January, 2019; originally announced January 2019.

Comments: Improved version, accepted to JPP Lecture Notes. Some parts were shortened and some parts were expanded. The text now contains Conclusions

Fast-to-Alfvén mode conversion mediated by Hall current. II Application to the solar atmosphere

Authors: P. A. González-Morales, E. Khomenko, P. S. Cally

Abstract: Coupling between fast magneto-acoustic and Alfvén waves can be observe in fully ionized plasmas mediated by stratification and 3D geometrical effects. In Paper I, Cally & Khomenko (2015) have shown that in a weakly ionized plasma, such as the solar photosphere and chromosphere, the Hall current introduces a new coupling mechanism. The present study extends the results from Paper I to the case of w… ▽ More Coupling between fast magneto-acoustic and Alfvén waves can be observe in fully ionized plasmas mediated by stratification and 3D geometrical effects. In Paper I, Cally & Khomenko (2015) have shown that in a weakly ionized plasma, such as the solar photosphere and chromosphere, the Hall current introduces a new coupling mechanism. The present study extends the results from Paper I to the case of warm plasma. We report on numerical experiments where mode transformation is studied using quasi-realistic stratification in thermodynamic parameters resembling the solar atmosphere. This redresses the limitation of the cold plasma approximation assumed in Paper I, in particular allowing the complete process of coupling between fast and slow magneto-acoustic modes and subsequent coupling of the fast mode to the Alfvén mode through the Hall current. Our results confirm the efficacy of the mechanism proposed in Paper I for the solar case. We observe that the efficiency of the transformation is a sensitive function of the angle between the wave propagation direction and the magnetic field, and of the wave frequency. The efficiency increases when the field direction and the wave direction are aligned for increasing wave frequencies. After scaling our results to typical solar values, the maximum amplitude of the transformed Alfvén waves, for a frequency of 1 Hz, corresponds to an energy flux (measured above the height of peak Hall coupling) of $\sim10^3$ $\rm W\,m^{-2}$, based on an amplitude of 500 $\rm m\,s^{-1}$ at $β=1$, which is sufficient to play a major role in both quiet and active region coronal heating. △ Less

Submitted 15 November, 2018; originally announced November 2018.

Spiral-shaped wavefronts in a sunspot umbra

Authors: T. Felipe, C. Kuckein, E. Khomenko, I. Thaler

Abstract: Solar active regions show a wide variety of oscillatory phenomena. The presence of the magnetic field leads to the appearance of several wave modes, whose behavior is determined by the sunspot thermal and magnetic structure. We aim to study the relation between the umbral and penumbral waves observed at the high photosphere and the magnetic field topology of the sunspot. Observations of the sunspo… ▽ More Solar active regions show a wide variety of oscillatory phenomena. The presence of the magnetic field leads to the appearance of several wave modes, whose behavior is determined by the sunspot thermal and magnetic structure. We aim to study the relation between the umbral and penumbral waves observed at the high photosphere and the magnetic field topology of the sunspot. Observations of the sunspot in active region NOAA 12662 obtained with the GREGOR telescope (Observatorio del Teide, Spain) were acquired on 2017 June 17. The data set includes a temporal series in the Fe I 5435 Å line obtained with the imaging spectrograph GREGOR Fabry-Pérot Interferometer (GFPI) and a spectropolarimetric raster map acquired with the GREGOR Infrared Spectrograph (GRIS) in the 10830 Å spectral region. The Doppler velocity deduced from the restored Fe I 5435 Å line has been determined, and the magnetic field vector of the sunspot has been inferred from spectropolarimetric inversions of the Ca I 10839 Å and the Si I 10827 Å lines. A two-armed spiral wavefront has been identified in the evolution of the two-dimensional velocity maps from the Fe I 5435 Å line. The wavefronts initially move counterclockwise in the interior of the umbra, and develop into radially outward propagating running penumbral waves when they reach the umbra-penumbra boundary. The horizontal propagation of the wavefronts approximately follows the direction of the magnetic field, which shows changes in the magnetic twist with height and horizontal position. The spiral wavefronts are interpreted as the visual pattern of slow magnetoacoustic waves which propagate upward along magnetic field lines. Their apparent horizontal propagation is due to their sequential arrival to different horizontal positions at the formation height of the Fe I 5435 Å line, as given by the inclination and orientation of the magnetic field. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Comments: Accepted for publication in A&A

Journal ref: A&A 621, A43 (2019)

Three-dimensional simulations of solar magneto-convection including effects of partial ionization

Authors: E. Khomenko, N. Vitas, M. Collados, A. de Vicente

Abstract: Over the last decades, realistic 3D radiative-MHD simulations have become the dominant theoretical tool for understanding the complex interactions between the plasma and the magnetic field on the Sun. Most of such simulations are based on approximations of magnetohydrodynamics, without directly considering the consequences of the very low degree of ionization of the solar plasma in the photosphere… ▽ More Over the last decades, realistic 3D radiative-MHD simulations have become the dominant theoretical tool for understanding the complex interactions between the plasma and the magnetic field on the Sun. Most of such simulations are based on approximations of magnetohydrodynamics, without directly considering the consequences of the very low degree of ionization of the solar plasma in the photosphere and bottom chromosphere. The presence of large amount of neutrals leads to a partial decoupling of the plasma and the magnetic field. As a consequence of that, a series of non-ideal effects (ambipolar diffusion, Hall effect and battery effect) arises. The ambipolar effect is the dominant one in the solar chromosphere. Here we report on the first three-dimensional realistic simulations of magneto-convection including ambipolar diffusion and battery effects. The simulations are done using the newly developed Mancha3D code. Our results reveal that ambipolar diffusion causes measurable effects on the amplitudes of waves excited by convection in the simulations, on the absorption of Poynting flux and heating and on the formation of chromospheric structures. We provide a low limit on the chromospheric temperature increase due to the ambipolar effect using the simulations with battery-excited dynamo fields. △ Less

Submitted 3 July, 2018; originally announced July 2018.

Comments: To appear in Astronomy & Astrophysics

Journal ref: A&A 618, A87 (2018)

MHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma

Authors: P. A. González-Morales, E. Khomenko, T. P. Downes, A. de Vicente

Abstract: The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere and low chromosphere. Under the single-fluid formalism the main non-ideal effects are described by a series of extra terms in the generalised induction equation… ▽ More The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere and low chromosphere. Under the single-fluid formalism the main non-ideal effects are described by a series of extra terms in the generalised induction equation and in the energy conservation equation. These effects are: Ohmic diffusion, ambipolar diffusion, the Hall effect, and the Biermann battery effect. From the point of view of the numerical solution of the single-fluid equations, when ambipolar diffusion or Hall effects dominate can introduce severe restrictions on the integration time step and can compromise the stability of the numerical scheme. In this paper we introduce two numerical schemes to overcome those limitations. The first of them is known as Super Time-Stepping (STS) and it is designed to overcome the limitations imposed when the ambipolar diffusion term is dominant. The second scheme is called the Hall Diffusion Scheme (HDS) and it is used when the Hall term becomes dominant. These two numerical techniques can be used together by applying Strang operator splitting. This paper describes the implementation of the STS and HDS schemes in the single-fluid code Mancha3D. The validation for each of these schemes is provided by comparing the analytical solution with the numerical one for a suite of numerical tests. △ Less

Submitted 13 March, 2018; originally announced March 2018.

Journal ref: A&A 615, A67 (2018)

Signatures of the impact of flare ejected plasma on the photosphere of a sunspot light-bridge

Authors: T. Felipe, M. Collados, E. Khomenko, S. P. Rajaguru, M. Franz, C. Kuckein, A. Asensio Ramos

Abstract: We investigate the properties of a sunspot light-bridge, focusing on the changes produced by the impact of a plasma blob ejected from a C-class flare. We observed a sunspot in active region NOAA 12544 using spectropolarimetric raster maps of the four Fe I lines around 15655 Å with the GREGOR Infrared Spectrograph (GRIS), narrow-band intensity images sampling the Fe I 6173 Å line with the GREGOR Fa… ▽ More We investigate the properties of a sunspot light-bridge, focusing on the changes produced by the impact of a plasma blob ejected from a C-class flare. We observed a sunspot in active region NOAA 12544 using spectropolarimetric raster maps of the four Fe I lines around 15655 Å with the GREGOR Infrared Spectrograph (GRIS), narrow-band intensity images sampling the Fe I 6173 Å line with the GREGOR Fabry-Pérot Interferometer (GFPI), and intensity broad band images in G-band and Ca II H band with the High-resolution Fast Imager (HiFI). All these instruments are located at the GREGOR telescope at the Observatorio del Teide, Tenerife, Spain. The data cover the time before, during, and after the flare event. The analysis is complemented with Atmospheric Imaging Assembly (AIA) and Helioseismic and Magnetic Imager (HMI) data from the Solar Dynamics Observatory (SDO). The physical parameters of the atmosphere at differents heights were inferred using spectral-line inversion techniques. We identify photospheric and chromospheric brightenings, heating events, and changes in the Stokes profiles associated to the flare eruption and the subsequent arrival of the plasma blob to the light bridge, after traveling along an active region loop. The measurements suggest that these phenomena are the result of reconnection events driven by the interaction of the plasma blob with the magnetic field topology of the light bridge. △ Less

Submitted 21 August, 2017; originally announced August 2017.

Comments: Accepted for publication in A&A

Journal ref: A&A 608, A97 (2017)

Partially Ionized Plasmas in Astrophysics

Authors: Jose Luis Ballester, Igor Alexeev, Manuel Collados, Turlough Downes, Robert F. Pfaff, Holly Gilbert, Maxim Khodachenko, Elena Khomenko, Ildar F. Shaikhislamov, Roberto Soler, Enrique Vazquez-Semadeni, Teimuraz Zaqarashvili

Abstract: Partially ionized plasmas are found across the Universe in many different astrophysical environments. They constitute an essential ingredient of the solar atmosphere, molecular clouds, planetary ionospheres and protoplanetary disks, among other environments, and display a richness of physical effects which are not present in fully ionized plasmas. This review provides an overview of the physics of… ▽ More Partially ionized plasmas are found across the Universe in many different astrophysical environments. They constitute an essential ingredient of the solar atmosphere, molecular clouds, planetary ionospheres and protoplanetary disks, among other environments, and display a richness of physical effects which are not present in fully ionized plasmas. This review provides an overview of the physics of partially ionized plasmas, including recent advances in different astrophysical areas in which partial ionization plays a fundamental role. We outline outstanding observational and theoretical questions and discuss possible directions for future progress. △ Less

Submitted 26 July, 2017; v1 submitted 25 July, 2017; originally announced July 2017.

Comments: 163 pages, submitted to Space Science Reviews

Numerical simulations of quiet Sun magnetic fields seeded by Biermann battery

Authors: E. Khomenko, N. Vitas, M. Collados, A. de Vicente

Abstract: The magnetic fields of the quiet Sun cover at any time more than 90\% of its surface and their magnetic energy budget is crucial to explain the thermal structure of the solar atmosphere. One of the possible origins of these fields is due to the action of local dynamo in the upper convection zone of the Sun. Existing simulations of the local solar dynamo require an initial seed field, and sufficien… ▽ More The magnetic fields of the quiet Sun cover at any time more than 90\% of its surface and their magnetic energy budget is crucial to explain the thermal structure of the solar atmosphere. One of the possible origins of these fields is due to the action of local dynamo in the upper convection zone of the Sun. Existing simulations of the local solar dynamo require an initial seed field, and sufficiently high spatial resolution, in order to achieve the amplification of the seed field to the observed values in the quiet Sun. Here we report an alternative model of seeding based on the action of the Bierman battery effect. This effect generates a magnetic field due to the local imbalances in electron pressure in the partially ionized solar plasma. We show that the battery effect self-consistently creates from zero an initial seed field of a strength of the order of micro G, and together with dynamo amplification, allows the generation of quiet Sun magnetic fields of a similar strength to those from solar observations. △ Less

Submitted 19 June, 2017; originally announced June 2017.

Comments: To appear in Astronomy & Astrophysics

Journal ref: A&A 604, A66 (2017)

High frequency waves in the corona due to null points

Authors: I. C. Santamaria, E. Khomenko, M. Collados, A. de Vicente

Abstract: This work aims to understand the behavior of non-linear waves in the vicinity of a coronal null point. In previous works we have showed that high frequency waves are generated in such magnetic configuration. This paper studies those waves in detail in order to provide a plausible explanation of their generation. We demonstrate that slow magneto-acoustic shock waves generated in the chromosphere pr… ▽ More This work aims to understand the behavior of non-linear waves in the vicinity of a coronal null point. In previous works we have showed that high frequency waves are generated in such magnetic configuration. This paper studies those waves in detail in order to provide a plausible explanation of their generation. We demonstrate that slow magneto-acoustic shock waves generated in the chromosphere propagate through the null point and produce a train of secondary shocks that escape along the field lines. A particular combination of the shock wave speeds generates waves at a frequency of 80 mHz. We speculate that this frequency may be sensitive to the atmospheric parameters in the corona and therefore can be used to probe the structure of this solar layer. △ Less

Submitted 21 April, 2017; originally announced April 2017.

Journal ref: A&A 602, A43 (2017)

Dependence of sunspot photospheric waves on the depth of the source of solar $p-$modes

Authors: T. Felipe, E. Khomenko

Abstract: Photospheric waves in sunspots moving radially outwards at speeds faster than the characteristic wave velocities have been recently detected. It has been suggested that they are the visual pattern of $p-$modes excited around 5 Mm beneath the sunspot. Using numerical simulations, we have performed a parametric study of the waves observed at the photosphere and higher layers produced by sources loca… ▽ More Photospheric waves in sunspots moving radially outwards at speeds faster than the characteristic wave velocities have been recently detected. It has been suggested that they are the visual pattern of $p-$modes excited around 5 Mm beneath the sunspot. Using numerical simulations, we have performed a parametric study of the waves observed at the photosphere and higher layers produced by sources located at different depths. The observational measurements are consistent with waves driven between approximately 1 Mm and 5 Mm below the sunspot surface. △ Less

Submitted 3 February, 2017; originally announced February 2017.

Comments: Accepted for publication in A&A

Journal ref: A&A 599, L2 (2017)

On the effects of ion-neutral interactions in solar plasmas

Authors: Elena Khomenko

Abstract: Solar photosphere and chromosphere are composed of weakly ionized plasma for which collisional coupling decreases with height. This implies a breakdown of some hypotheses underlying magnetohydrodynamics at low altitudes and gives rise to non-ideal MHD effects such as ambipolar diffusion, Hall effect, etc. Recently, there has been progress in understanding the role of these effects for the dynamics… ▽ More Solar photosphere and chromosphere are composed of weakly ionized plasma for which collisional coupling decreases with height. This implies a breakdown of some hypotheses underlying magnetohydrodynamics at low altitudes and gives rise to non-ideal MHD effects such as ambipolar diffusion, Hall effect, etc. Recently, there has been progress in understanding the role of these effects for the dynamics and energetics of the solar atmosphere. There are evidences that such phenomena as wave propagation and damping, magnetic reconnection, formation of stable magnetic field concentrations, magnetic flux emergence, etc. can be affected. This paper reviews the current state-of-the-art of multi-fluid MHD modeling of the coupled solar atmosphere. △ Less

Submitted 18 November, 2016; originally announced November 2016.

Comments: "Plasma Physics and Controlled Fusion", Volume 59, Number 1, 2016; Special issue featuring the invited talks from the 43rd EPS Conference on Plasma Physics, Leuven, 4-8 July 2016

Three-dimensional structure of a sunspot light bridge

Authors: T. Felipe, M. Collados, E. Khomenko, C. Kuckein, A. Asensio Ramos, H. Balthasar, T. Berkefeld, C. Denker, A. Feller, M. Franz, A. Hofmann, C. Kiess, A. Lagg, H. Nicklas, D. Orozco Suárez, A. Pastor Yabar, R. Rezaei, R. Schlichenmaier, D. Schmidt, W. Schmidt, M. Sigwarth, M. Sobotka, S. K. Solanki, D. Soltau, J. Staude , et al. (4 additional authors not shown)

Abstract: Active regions are the most prominent manifestations of solar magnetic fields; their generation and dissipation are fundamental problems in solar physics. Light bridges are commonly present during sunspot decay, but a comprehensive picture of their role in the removal of photospheric magnetic field is still missing. We study the three dimensional configuration of a sunspot and in particular its li… ▽ More Active regions are the most prominent manifestations of solar magnetic fields; their generation and dissipation are fundamental problems in solar physics. Light bridges are commonly present during sunspot decay, but a comprehensive picture of their role in the removal of photospheric magnetic field is still missing. We study the three dimensional configuration of a sunspot and in particular its light bridge during one of the last stages of its decay. We present the magnetic and thermodynamical stratification inferred from full Stokes inversions of the photospheric Si I 10827 Å and Ca I 10839 Å lines obtained with the GREGOR Infrared Spectrograph of the GREGOR telescope at Observatorio del Teide, Tenerife, Spain. The analysis is complemented by a study of continuum images covering the disk passage of the active region, which are provided by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. The sunspot shows a light bridge with penumbral continuum intensity that separates the central umbra from a smaller umbra. We find that in this region the magnetic field lines form a canopy with lower magnetic field strength in the inner part. The photospheric light bridge is dominated by gas pressure (high-$β$), as opposed to the surrounding umbra where the magnetic pressure is higher. A convective flow is observed in the light bridge. This flow is able to bend the magnetic field lines and to produce field reversals. The field lines close above the light bridge and become as vertical and strong as in the surrounding umbra. We conclude that it develops because of two highly magnetized regions which come closer during the sunspot evolution. △ Less

Submitted 15 November, 2016; originally announced November 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 596, A59 (2016)

Simulated interaction of MHD shock waves with a complex network-like region

Authors: Irantzu C. Santamaria, Elena Khomenko, Manuel Collados, Angel de Vicente

Abstract: We provide estimates of the wave energy reaching the solar chromosphere and corona in a network-like magnetic field topology, including a coronal null point. The waves are excited by an instantaneous strong subphotospheric source and propagate through the subphotosphere, photosphere, chromosphere, transition region, and corona with the plasma beta and other atmospheric parameters varying by severa… ▽ More We provide estimates of the wave energy reaching the solar chromosphere and corona in a network-like magnetic field topology, including a coronal null point. The waves are excited by an instantaneous strong subphotospheric source and propagate through the subphotosphere, photosphere, chromosphere, transition region, and corona with the plasma beta and other atmospheric parameters varying by several orders of magnitude. We compare two regimes of the wave propagation: a linear and nonlinear regime. While the amount of energy reaching the corona is similar in both regimes, this energy is transmitted at different frequencies. In both cases the dominant periods of waves at each height strongly depend on the local magnetic field topology, but this distribution is only in accordance with observations in the nonlinear case. △ Less

Submitted 29 April, 2016; originally announced April 2016.

Comments: 4 pages

Journal ref: A&A 590, L3 (2016)

Observational detection of drift velocity between ionized and neutral species in solar prominences

Authors: Elena Khomenko, Manuel Collados, Antonio J. Diaz

Abstract: We report a detection of differences in ion and neutral velocities in prominences using high resolution spectral data obtained in September 2012 at the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife). A time series of scans of a small portion of a solar prominence was obtained simultaneously with a high cadence using the lines of two elements with different ionization states, name… ▽ More We report a detection of differences in ion and neutral velocities in prominences using high resolution spectral data obtained in September 2012 at the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife). A time series of scans of a small portion of a solar prominence was obtained simultaneously with a high cadence using the lines of two elements with different ionization states, namely the CaII 8542 A and the HeI 10830 A. Displacements, widths and amplitudes of both lines were carefully compared to extract dynamical information about the plasma. Many dynamical features are detected, such as counterstreaming flows, jets and propagating waves. In all the cases we find very strong correlation between the parameters extracted from the lines of both elements, confirming that both trace the same plasma. Nevertheless, we also find short-lived transients where this correlation is lost. These transients are associated with the ion-neutral drift velocities of the order of several hundred m/s. The patches of non-zero drift velocity show coherence on time-distance diagrams. △ Less

Submitted 5 April, 2016; originally announced April 2016.

Comments: To appear in The Astrophysical Journal

Heating of the partially ionized solar chromosphere by waves in magnetic structures

Authors: S. Shelyag, E. Khomenko, A. de Vicente, D. Przybylski

Abstract: In this paper, we show a "proof of concept" of the heating mechanism of the solar chromosphere due to wave dissipation caused by the effects of partial ionization. Numerical modeling of non-linear wave propagation in a magnetic flux tube, embedded in the solar atmosphere, is performed by solving a system of single-fluid quasi-MHD equations, which take into account the ambipolar term from the gener… ▽ More In this paper, we show a "proof of concept" of the heating mechanism of the solar chromosphere due to wave dissipation caused by the effects of partial ionization. Numerical modeling of non-linear wave propagation in a magnetic flux tube, embedded in the solar atmosphere, is performed by solving a system of single-fluid quasi-MHD equations, which take into account the ambipolar term from the generalized Ohm's law. It is shown that perturbations caused by magnetic waves can be effectively dissipated due to ambipolar diffusion. The energy input by this mechanism is continuous and shown to be more efficient than dissipation of static currents, ultimately leading to chromospheric temperature increase in magnetic structures. △ Less

Submitted 10 February, 2016; originally announced February 2016.

Comments: 7 pages, 5 figures

On the origin of facular brightness

Authors: R. Kostik, E. Khomenko

Abstract: This paper studies the dependence of the CaIIH line core brightness on the strength and inclination of photospheric magnetic field, and on the parameters of convective and wave motions in a facular region at the solar disc center. We use three simultaneous datasets obtained at the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife): (1) spectra of BaII 4554 A line registered with the… ▽ More This paper studies the dependence of the CaIIH line core brightness on the strength and inclination of photospheric magnetic field, and on the parameters of convective and wave motions in a facular region at the solar disc center. We use three simultaneous datasets obtained at the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife): (1) spectra of BaII 4554 A line registered with the instrument TESOS to measure the variations of intensity and velocity through the photosphere up to the temperature minimum; (2) spectropolarimetric data in FeI 1.56 $μ$m lines (registered with the instrument TIP II) to measure photospheric magnetic fields; (3) filtergrams in CaIIH that give information about brightness fluctuations in the chromosphere. The results show that the CaIIH brightness in the facula strongly depends on the power of waves with periods in the 5-min range, that propagate upwards, and also on the phase shift between velocity oscillations at the bottom photosphere and around the temperature minimum height, measured from BaII line. The CaIIH brightness is maximum at locations where the phase shift between temperature and velocity oscillations lies within 0-100 degrees. There is an indirect influence of convective motions on the CaIIH brightness. Namely, the higher is the amplitude of convective velocities and the larger is the height where they change their direction of motion, the brighter is the facula. Altogether, our results lead to conclusions that facular regions appear bright not only because of the Wilson depression in magnetic structures, but also due to real heating. △ Less

Submitted 10 February, 2016; originally announced February 2016.

Comments: to appear in Astronomy & Astrophysics

Journal ref: A&A 589, A6 (2016)