Search

Understanding the effects of the lower solar atmosphere on the spectrum of standing kink oscillations of coronal loops, in both the decaying and decayless regime, is essential for developing more advanced tools for coronal seismology. We aim to reveal the effects of the chromosphere on the spatial profiles and frequencies of the standing kink modes, create synthetic emission maps to compare with observations, and study the results using spatial and temporal coronal seismology techniques. We excited transverse oscillations in a 3D straight flux tube using (a) a broadband footpoint driver, (b) a sinusoidal velocity pulse, and (c) an off-centre Gaussian velocity pulse, using the PLUTO code. The flux tube is gravitationally stratified, with footpoints embedded in chromospheric plasma. Using the FoMo code, we created synthetic observations of our data in the Fe IX 17.1 nm line and calculated the spectra with the Auto-NUWT code. We also numerically solved the generalised eigenvalue system for the 1D wave equation to determine the effects of the stratification on the kink modes of our system. The synthetic observations of the loops perturbed by the velocity pulses show a single dominant mode that our 1D analysis reveals to be the third harmonic of the system. For the broadband driver, the synthetic emission shows multiple frequency bands, associated with both the loop and the driver. Finally, using seismological techniques, we highlight the possibility of misidentifying the observed third, sixth, and ninth harmonics with the first, second, and third harmonics of the coronal part of longer loops. Unless more advanced techniques of spatial seismology are used with many data points from observations along the entire loop length, this misidentification can result in overestimating the mean magnetic field by a factor equal to the period ratio of the fundamental over the third harmonic., Comment: 14 pages, 9 figures, 1 appendix, accepted for publication in A&A

Recent observations have demonstrated the capability of mapping the solar coronal magnetic field using the technique of coronal seismology based on the ubiquitous propagating Alfvenic/kink waves through imaging spectroscopy. We established a magnetohydrodynamic (MHD) model of a gravitationally stratified open magnetic flux tube, exciting kink waves propagating upwards along the tube. Forward modeling was performed to synthesize the Fe XIII 1074.7 and 1079.8 nm spectral line profiles, which were then used to determine the wave phase speed, plasma density, and magnetic field with seismology method. A comparison between the seismologically inferred results and the corresponding input values verifies the reliability of the seismology method. In addition, we also identified some factors that could lead to errors during magnetic field measurements. Our results may serve as a valuable reference for current and future coronal magnetic field measurements based on observations of propagating kink waves., Comment: Accepted for publication in RAA

Context. In the coronal open-field regions, such as coronal holes, there are many transverse waves propagating along magnetic flux tubes, generally interpreted as kink waves. Previous studies have highlighted their potential in coronal heating, solar wind acceleration, and seismological diagnostics of various physical parameters. Aims. This study aims to investigate propagating kink waves, considering both vertical and horizontal density inhomogeneity, using three-dimensional magnetohydrodynamic (MHD) simulations. Methods. We establish a 3D MHD model of a gravitationally stratified open flux tube, incorporating a velocity driver at the lower boundary to excite propagating kink waves. Forward modelling is conducted to synthesise observational signatures of the Fe ix 17.1 nm line. Results. It is found that resonant absorption and density stratification both affect the wave amplitude. When diagnosing the relative density profile with velocity amplitude, resonant damping needs to be properly considered to avoid possible underestimation. In addition, unlike standing modes, propagating waves are believed to be Kelvin-Helmholtz stable. In the presence of vertical stratification, however, phase mixing of transverse motions around the tube boundary can still induce small scales, partially dissipating wave energy and leading to a temperature increase, especially at higher altitudes. Moreover, forward modeling is conducted to synthesise observational signatures, revealing the promising potential of future coronal imaging spectrometers such as MUSE in resolving these wave-induced signatures. Also, the synthesised intensity signals exhibit apparent periodic variations, offering a potential method to indirectly observe propagating kink waves with current EUV imagers., Comment: accepted for publication in A&A

Instabilities in oscillating loops are believed to be essential for dissipating the wave energy and heating the solar coronal plasma. Our aim is to study the development of the Kelvin-Helmholtz (KH) instability in an oscillating loop that is driven by random footpoint motions. Using the PLUTO code, we performed 3D simulations of a straight gravitationally stratified flux tube. The loop footpoints are embedded in chromospheric plasma, in the presence of thermal conduction and an artificially broadened transition region. Using drivers with a power-law spectrum, one with a red noise spectrum and one with the low-frequency part subtracted, we excited standing oscillations and the KH instability in our loops, after one-and-a-half periods of the oscillation. We see that our broadband drivers lead to fully deformed, turbulent loop cross-sections over the entire coronal part of the loop due to the spatially extended KH instability. The low RMS velocity of our driver without the low-frequency components supports the working hypothesis that the KH instability can easily manifest in oscillating coronal loops. We report for the first time in driven transverse oscillations of loops the apparent propagation of density perturbations due to the onset of the KH instability, from the apex towards the footpoints. Both drivers input sufficient energy to drive enthalpy and mass flux fluctuations along the loop, while also causing heating near the driven footpoint of the oscillating loop, which becomes more prominent when a low-frequency component is included in the velocity driver. Finally, our power-law driver with the low-frequency component provides a RMS input Poynting flux of the same order as the radiative losses of the quiet-Sun corona, giving us promising prospects for the contribution of decayless oscillations in coronal heating., Comment: 12 pages, 12 figures

Aims. We studied the manifestation of decayless oscillations in 3D simulations of coronal loops, driven by random motions. Methods. Using the PLUTO code, we ran magnetohydrodynamic (MHD) simulations of a straight gravitationally stratified flux tube, with its footpoints embedded in chromospheric plasma. We consider transverse waves drivers with a horizontally polarised red noise power-law spectrum. Results. Our broadband drivers lead to the excitation of standing waves with frequencies equal to the fundamental standing kink mode and its harmonics. These standing kink oscillations have non-decaying amplitudes, and spectra that depend on the characteristics of the loops, with the latter amplifying the resonant frequencies from the drivers. We thus report for the first time in 3D simulations the manifestation of decayless oscillations from broadband drivers. The spatial and temporal evolution of our oscillation spectra reveals the manifestation of a half harmonic, which exhibits half the frequency of the identified fundamental mode with a similar spatial profile. Our results suggest that this mode is related to the presence of the transition region in our model and could be interpreted as being the equivalent to the fundamental mode of standing sound waves driven on pipes closed at one end. Conclusions. The potential existence of this half harmonic has important implications for coronal seismology, since misinterpreting it for the fundamental mode of the system can lead to false estimations of the average kink speed profile along oscillating loops. Finally, its detection could potentially give us a tool for distinguishing between different excitation and driving mechanisms of decayless oscillations in observations., Comment: 6 pages, 5 figures

We model a coronal loop as a three-dimensional magnetic cylinder in a realistic solar atmosphere that extends from the chromosphere to the corona. Kink oscillations, believed ubiquitous in the solar corona, are launched in the loop. Heating is expected due to the dissipation of wave energy at small structures that develop from the Kelvin-Helmholtz instability induced by kink oscillations. Increases in temperature and internal energy can be observed in the coronal counterpart of the driven loop. With the presence of thermal conduction, chromospheric evaporation can also be seen. Although the volume averaged temperature and density changes seem slight ($\sim4\%$ relative to a non-driven loop), the enthalpy flow from the lower atmosphere redistributes the density and temperature in the vertical direction, thus enhancing the dissipation of wave energy in the corona. The efficient heating in the coronal counterpart of the loop can complement the thermal conductive losses shown in the current model and thus maintain the internal energy in the corona., Comment: Accepted for publication in ApJL

Context. Transverse oscillations are ubiquitously observed in the solar corona, both in coronal loops and open magnetic flux tubes. Numerical simulations suggest that their dissipation could heat coronal loops, counterbalancing radiative losses. These models rely on a continuous driver at the footpoint of the loops. However, analytical works predict that transverse waves are subject to a cut-off in the transition region. It is thus unclear whether they can reach the corona, and indeed heat coronal loops. Aims. Our aims are to determine how the cut-off of kink waves affects their propagation into the corona, and to characterize the variation of the cut-off frequency with altitude. Methods. Using 3D magnetohydrodynamic simulations, we modelled the propagation of kink waves in a magnetic flux tube, embedded in a realistic atmosphere with thermal conduction, that starts in the chromosphere and extends into the corona. We drove kink waves at four different frequencies, and determined whether they experienced a cut-off. We then calculated the altitude at which the waves were cut-off, and compared it to the prediction of several analytical models. Results. We show that kink waves indeed experience a cut-off in the transition region, and we identified the analytical model that gives the best predictions. In addition, we show that waves with periods shorter than approximately 500 s can still reach the corona by tunnelling through the transition region, with little to no attenuation of their amplitude. This means that such waves can still propagate from the footpoints of loop, and result in heating in the corona., Comment: Accepted for publication in A&A. 8 pages, 7 figures

Oscillatory reconnection is a relaxation process in magnetised plasma, with an inherent periodicity that is exclusively dependent on the properties of the background plasma. This study focuses on the seismological prospects of oscillatory reconnection in the solar corona. We perform three sets of parameter studies (for characteristic coronal values of the background magnetic field, density and temperature) using the PLUTO code to solve the fully compressive, resistive MHD equations for a 2D magnetic X-point. From each parameter study, we derive the period of the oscillatory reconnection. We find that this period is inversely proportional to the characteristic strength of the background magnetic field and the square root of the initial plasma temperature, while following a square root dependency upon the equilibrium plasma density. These results reveal an inverse proportionality between the magnitude of the Alfv\'en speed and the period, as well as the background sound speed and the period. Furthermore, we note that the addition of anisotropic thermal conduction only leads to a small increase in the mean value for the period. Finally, we establish an empirical formula that gives the value for the period in relation to the background magnetic field, density and temperature. This gives us a quantified relation for oscillatory reconnection, to be used as a plasma diagnostic in the solar corona, opening up the possibility of using oscillatory reconnection for coronal seismology., Comment: 16 pages, 10 figures, accepted for publication in ApJ

Oscillatory reconnection can manifest through the interaction between the ubiquitous MHD waves and omnipresent null points in the solar atmosphere and is characterized by an inherent periodicity. In the current study, we focus on the relationship between the period of oscillatory reconnection and the strength of the wave pulse initially perturbing the null point, in a hot coronal plasma. We use the PLUTO code to solve the fully compressive, resistive MHD equations for a 2D magnetic X-point. Using wave pulses with a wide range of amplitudes, we perform a parameter study to obtain values for the period, considering the presence and absence of anisotropic thermal conduction separately. In both cases, we find that the resulting period is independent of the strength of the initial perturbation. The addition of anisotropic thermal conduction only leads to an increase in the mean value for the period, in agreement with our previous study. We also consider a different type of initial driver and we obtain an oscillation period matching the independent trend previously mentioned. Thus, we report for the first time on the independence between the type and strength of the initializing wave pulse and the resulting period of oscillatory reconnection in a hot coronal plasma. This makes oscillatory reconnection a promising mechanism to be used within the context of coronal seismology., Comment: 19 pages, 14 figures, accepted for publication in ApJ

Oscillatory reconnection (a relaxation mechanism with periodic changes in connectivity) has been proposed as a potential physical mechanism underpinning several periodic phenomena in the solar atmosphere including, but not limited to, quasi-periodic pulsations (QPPs). Despite its importance, however, the mechanism has never been studied within a hot, coronal plasma. We investigate oscillatory reconnection in a one million Kelvin plasma by solving the fully-compressive, resistive MHD equations for a 2D magnetic X-point under coronal conditions using the PLUTO code. We report on the resulting oscillatory reconnection including its periodicity and decay rate. We observe a more complicated oscillating profile for the current density compared to that found for a cold plasma, due to mode-conversion at the equipartition layer. We also consider, for the first time, the effect of adding anisotropic thermal conduction to the oscillatory reconnection mechanism, and we find this simplifies the spectrum of the oscillation profile and increases the decay rate. Crucially, the addition of thermal conduction does not prevent the oscillatory reconnection mechanism from manifesting. Finally, we reveal a relationship between the equilibrium magnetic field strength, decay rate, and period of oscillatory reconnection, which opens the tantalising possibility of utilizing oscillatory reconnection as a seismological tool., Comment: 17 pages, 13 figures, accepted for publication in ApJ

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

Identifying the underlying mechanisms behind the excitation of transverse oscillations in coronal loops is essential for their role as diagnostic tools in coronal seismology and their potential use as wave heating mechanisms of the solar corona. In this paper, we explore the concept of these transverse oscillations being excited through a self-sustaining process, caused by Alfv\'{e}nic vortex shedding from strong background flows interacting with coronal loops. We show for the first time in 3D simulations that vortex shedding can generate transverse oscillations in coronal loops, in the direction perpendicular to the flow due to periodic "pushing" by the vortices. By plotting the power spectral density we identify the excited frequencies of these oscillations. We see that these frequencies are dependent both on the speed of the flow, as well as the characteristics of the oscillating loop. This, in addition to the fact that the background flow is constant and not periodic, makes us treat this as a self-oscillating process. Finally, the amplitudes of the excited oscillations are near constant in amplitude, and are comparable with the observations of decay-less oscillations. This makes the mechanism under consideration a possible interpretation of these undamped waves in coronal loops., Comment: 8 pages, 8 figures, accepted for publication in ApJL

In the quest to solve the long-standing coronal heating problem, it has been suggested half a century ago that coronal loops could be heated by waves. Despite the accumulating observational evidence of the possible importance of coronal waves, still no 3D MHD simulations exist that show significant heating by MHD waves. Here we report on the first 3D coronal loop model heating the plasma against radiative cooling. The coronal loop is driven at the footpoint by transverse oscillations and subsequently the induced Kelvin-Helmholtz instability deforms the loop cross-section and generates small-scale structures. Wave energy is transfered to smaller scales where it is dissipated, overcoming the internal energy losses by radiation. These results open up a new avenue to address the coronal heating problem., Comment: 17 pages, 8 figures, accepted for publication in ApJ

In recent years, the decay-less regime of standing transverse oscillations in coronal loops has been the topic of many observational and numerical studies, focusing on their physical characteristics, as well as their importance for coronal seismology and wave heating. However, no definitive answer has yet been given on the driving mechanism behind these oscillations, with most studies focusing on the use of periodic footpoint drivers as a means to excite them. In this paper, our goal is to explore the concept of these standing waves being self-sustained oscillations, driven by a constant background flow. To that end, we use the PLUTO code, to perform $3$D magnetohydrodynamic simulations of a gravitationally stratified straight flux tube in a coronal environment, in the presence of a weak flow around the loop, perpendicular to its axis. Once this flow is firmly set up, a transverse oscillation is initiated, dominated by the fundamental kink mode of a standing wave, while the existence of a second harmonic is revealed, with a frequency ratio to the fundamental mode near the observed ones in decay-less oscillations. The presence of vortex shedding is also established in our simulations, which is connected to the "slippery" interaction between the oscillator and its surrounding plasma. We thus present a proof-of-concept of a self-oscillation in a coronal loop, and we propose it as a mechanism that could interpret the observed decay-less transverse oscillations of coronal loops., Comment: 7 pages, 3 figures

Kelvin-Helmholtz {instability induced} turbulence is one promising mechanism by which loops in the solar corona can be heated by MHD waves. In this paper we present an analytical model of the dissipation rate of {Kelvin-Helmholtz instability induced} turbulence $\varepsilon_{\rm D}$, finding it scales as the wave amplitude ($d$) to the third power ($\varepsilon_{\rm D}\propto d^3$). Based on the concept of steady-state turbulence, we expect the turbulence heating throughout the volume of {the} loop to match the total energy injected through its footpoints. In situations where this holds, the wave amplitude has to vary as the cube-root of the injected energy. Comparing the analytic results with those of simulations shows that our analytic formulation captures the key aspects of the turbulent dissipation from the numerical work. Applying this model to the observed characteristics of decayless kink waves we predict that the amplitudes of these observed waves is insufficient to turbulently heat the solar corona., Comment: 6 pages, 1 figure, published in ApJL

It has been found that the Kelvin-Helmholtz instability (KHI) induced by both transverse and torsional oscillations in coronal loops can reinforce the effects of wave heating. In this study, we model a coronal loop as a system of individual strands, and we study wave heating effects by considering a combined transverse and torsional driver at the loop footpoint. We deposit the same energy into the multi-stranded loop and an equivalent monolithic loop, and then observe a faster increase in the internal energy and temperature in the multi-stranded model. Therefore, the multi-stranded model is more efficient in starting the heating process. Moreover, higher temperature is observed near the footpoint in the multi-stranded loop and near the apex in the monolithic loop. The apparent heating location in the multi-stranded loop agrees with the previous predictions and observations. Given the differences in the results from our multi-stranded loop and monolithic loop simulations, and given that coronal loops are suggested to be multi-stranded on both theoretical and observational grounds, our results suggest that the multi-strandedness of coronal loops needs to be incorporated in future wave-based heating mechanisms., Comment: 8 figures, accepted for publication in ApJ

We analyse transverse oscillations of a coronal loop excited by continuous monoperiodic motions of the loop footpoint at different frequencies in the presence of gravity. Using the MPI-AMRVAC code, we perform three-dimensional numerical magnetohydrodynamic simulations, considering the loop as a magnetic flux tube filled in with denser, hotter, and gravitationally stratified plasma. We show the resonant response of the loop to its external excitation and analyse the development of the Kelvin-Helmholtz instability at different heights. We also study the spatial distribution of plasma heating due to transverse oscillations along the loop. The positions of the maximum heating are in total agreement with those for the intensity of the Kelvin-Helmholtz instability, and correspond to the standing wave anti-nodes in the resonant cases. The initial temperature configuration and plasma mixing effect appear to play a significant role in plasma heating by transverse footpoint motions. In particular, the development of the Kelvin-Helmholtz instability in a hotter loop results in the enhancement of the mean plasma temperature in the domain., Comment: Published in ApJ