Dynamic modeling of coronal abundances during flares on M-dwarf stars

Solar atmospheric elemental abundances are now known to vary both in space and time. Dynamic modeling of these changes is therefore necessary to improve the accuracy of radiative hydrodynamic simulations. Recent studies have shown that including spatio-temporal variations in coronal abundances during solar flares leads to the formation of coronal condensations (rain), which are otherwise difficult to create in impulsively heated field aligned hydrodynamic flare models. These simulations start with a solar corona dominated by the first ionization potential (FIP) effect, and evaporate photospheric material into the post-flare loops. We here explore perhaps the most extreme non-solar starting condition for the coronal composition in these simulations: an initial corona dominated by the inverse FIP (iFIP) effect, such as is observed on active M-dwarf stars. We show that a flaring event in a corona enriched with high FIP elements leads to a solution similar to the solar case. Coronal rain is harder to form by this method during flares on M-dwarfs, however, if the corona is depleted of low FIP elements.

💡 Research Summary

This paper investigates how time‑dependent elemental abundances affect the dynamics of flare loops on active M‑dwarf stars, focusing on the formation of coronal rain. Recent solar studies have shown that allowing the coronal composition to evolve during a flare can trigger a localized increase in radiative losses at the loop apex, leading to catastrophic cooling and condensation (coronal rain). Those studies began with a solar‑like corona dominated by the First Ionization Potential (FIP) effect, where low‑FIP elements (e.g., Fe, Si, Mg) are enhanced relative to photospheric values. The authors ask whether the same mechanism operates under the most extreme opposite condition: an inverse FIP (iFIP) corona, as observed in many active M‑dwarfs, where high‑FIP elements are enriched or low‑FIP elements are depleted.

To answer this, they use the one‑dimensional hydrodynamic code HYDRAD, which solves the full set of mass, momentum, and energy equations along a magnetic flux tube, includes optically thick chromospheric radiation, and calculates optically thin radiative losses using the CHIANTI v10 database for the 15 most abundant elements. A novel feature is the implementation of a spatially and temporally varying abundance factor f(t,s) that obeys the advection equation ∂f/∂t + v ∂f/∂s = 0, allowing the composition to be advected with the bulk plasma flow. The authors adopt photospheric abundances from Asplund et al. (2009) as the baseline.

Two scenarios are simulated, each using the same flare heating parameters to isolate the effect of the initial coronal composition:

1. High‑FIP enrichment: The corona is initially enriched in high‑FIP elements (Ne, Ar, etc.) by a factor f_H = 4, while low‑FIP elements retain photospheric values.
2. Low‑FIP depletion: The corona is initially depleted in low‑FIP elements by a factor f_L = 0.25, while high‑FIP elements are at photospheric levels.

Both simulations employ a loop length of 150 Mm, appropriate for the larger magnetic structures typical of M‑dwarfs, and an electron‑beam heating profile identical to that used in previous solar studies: a constant energy flux of 2 × 10¹⁰ erg cm⁻² s⁻¹ applied for 10 s, with a low‑energy cutoff of 15 keV and a spectral index of 5. Adaptive mesh refinement ensures adequate resolution of the transition region.

Results – High‑FIP enrichment:
When the flare begins, chromospheric evaporation injects photospheric plasma into the loop. Because the initial coronal plasma is rich in high‑FIP elements, the abundance factor f_H remains high near the loop apex for several hundred seconds, creating a localized peak in the radiative loss function. This peak accelerates cooling at the apex, producing a pressure dip that draws additional material upward, raising the density locally. The runaway cooling culminates in a condensation at the loop top after roughly 5000 s (longer than the ~2000 s found in comparable solar simulations, primarily due to the longer loop). The condensation is identified as coronal rain. Thus, even with an iFIP corona dominated by high‑FIP enrichment, the mechanism that generates coronal rain in solar flares remains operative.

Results – Low‑FIP depletion:
In the second case, the initial low‑FIP depletion reduces the overall radiative loss rate to about 25 % of the photospheric value. As evaporation proceeds, the abundance factor f_L rises toward unity, but a localized dip in f_L persists at the apex, leading to a reduction—not an increase—in radiative losses there. Consequently, the apex cools more slowly, the pressure remains relatively uniform, and no condensation forms. The loop cools and drains in a more homogeneous fashion, mirroring simulations that assume fixed radiative losses.

Discussion and Implications:
The contrasting outcomes demonstrate that the sign of the abundance perturbation (enhancement vs. depletion) critically determines whether a flare loop will experience the runaway radiative cooling necessary for coronal rain. High‑FIP enrichment mimics the solar FIP case because the radiative loss curve steepens, whereas low‑FIP depletion flattens the curve and suppresses the instability. This sensitivity suggests that observations of coronal rain (or its absence) on M‑dwarfs could be used as diagnostics of the underlying coronal composition and, by extension, of the processes that generate the iFIP effect (e.g., ponderomotive forces, wave‑driven fractionation).

The authors acknowledge that their treatment of abundance evolution is simplified: the advection equation assumes perfect mixing with bulk flows and neglects diffusion, non‑equilibrium ionization, and magnetic field variations that could modify the transport of elements. Nevertheless, the approach captures the essential physics on the timescales relevant to flare evolution and is readily implementable in existing 1‑D codes. Future work should extend the model to three‑dimensional radiative MHD, incorporate more realistic fractionation mechanisms, and compare synthetic spectra with high‑resolution X‑ray/EUV observations from missions such as XMM‑Newton, XRISM, and upcoming facilities.

Conclusion:
The study confirms that an iFIP corona enriched in high‑FIP elements can still produce coronal rain during a flare, while a corona depleted of low‑FIP elements makes rain formation unlikely under the same heating conditions. These findings highlight the importance of incorporating dynamic abundance variations into flare models, especially for stars where the inverse FIP effect dominates, and provide a framework for interpreting future multi‑wavelength observations of M‑dwarf flares.