The solar corona as an active medium for magnetoacoustic waves

The presence and interplay of continuous cooling and heating processes maintaining the corona of the Sun at the observed one million K temperature were recently understood to have crucial effects on the dynamics and stability of magnetoacoustic waves. These essentially compressive waves perturb the coronal thermal equilibrium, leading to the phenomenon of a wave-induced thermal misbalance. Representing an additional natural mechanism for the exchange of energy between the plasma and the wave, thermal misbalance makes the corona an active medium for magnetoacoustic waves, so that the wave can not only lose but also gain energy from the coronal heating source (similarly to burning gases, lasers and masers). We review recent achievements in this newly emerging research field, focussing on the effects that slow-mode magnetoacoustic waves experience as a back-reaction of this perturbed coronal thermal equilibrium. The new effects include enhanced frequency-dependent damping or amplification of slow waves, and effective, not associated with the coronal plasma non-uniformity, dispersion. We also discuss the possibility to probe the unknown coronal heating function by observations of slow waves and linear theory of thermal instabilities. The manifold of the new properties that slow waves acquire from a thermodynamically active nature of the solar corona indicate a clear need for accounting for the effects of combined coronal heating/cooling processes not only for traditional problems of the formation and evolution of prominences and coronal rain, but also for an adequate modelling and interpretation of magnetohydrodynamic waves.


Introduction
The coronal heating problem remains one of the major puzzles in solar physics for almost eighty years since the pioneering works [1,2], which apparently makes it the longest-standing unsolved problem in all plasma astrophysics. Indeed, estimations show that 1 g of the coronal plasma with typical temperature of 1 MK radiates in the optically thin regime more than 10 11 erg s −1 , which would lead to the cooling of the corona in a few hours unless the radiative energy losses and parallel thermal conduction towards the chromosphere are re-supplied by some unknown yet heating mechanism. In reality, the observed lifetime of typical hot coronal plasma structures is much longer than the expected radiative cooling time (see e.g. [3,4], for comprehensive reviews of the properties of the corona and coronal loops), which indicates the Sun's corona has to be considered as a continuously cooling and heated medium, existing because of the everlasting competition and a delicate balance between these two processes.
The question of coronal heating is traditionally associated with the dynamics of magnetohydrodynamic (MHD) waves ubiquitously present in the corona (see e.g. [5], and also [6,7,8], for comprehensive reviews of specific coronal wave modes). The intrinsically filamentary nature of the coronal plasma and the associated with it generation of small spatial scales are considered as key ingredients for an effective dissipation of the wave energy to heat the plasma [9,10]. Despite an enormous effort in modelling and observational studies of coronal heating by MHD waves (see e.g. [11,12,13], for recent works) and sporadic indirect estimations that in some cases the wave energy could be sufficient (e.g. [14]), the results converge to the conclusion that realistic 3D MHD models accounting for cross-field coronal plasma inhomogeneities and observed wave amplitudes are not yet capable to reproduce the heating rate required to balance the colossal radiative losses of the corona [15].
In this review, we address the link between the coronal heating and MHD waves from the wave dynamics point of view, i.e. focussing on the effects that the coronal heating and cooling processes exert on the evolution of MHD waves. In other words, here we do not consider the waves as the heating agent, but discuss the implications of combined heating/cooling effects for the wave dynamics and a practical question of seismological diagnostics of thermal properties of the corona, including its enigmatic heating function. A vast majority of coronal wave modelling studies is carried out either in terms of ideal MHD, or with unrealistically high transport coefficients intrinsic for 3D numerical modelling. Thus, responding to the need to go beyond ideal MHD theory for an adequate description and interpretation of the wave processes in the corona, Ref. [16] considered the dynamics of fast magnetoacoustic kink waves in coronal loops undergoing rapid radiative cooling. In this approach, the hot plasma of the oscillating coronal loop was assumed to be not in a hydrostatic equilibrium, cooling down by intensive radiation which is not balanced by heating (see the schematic illustration in the left-hand panel of Fig. 1). Similar studies of the effects of radiative cooling on the damping profile of kink oscillations and on the behaviour of standing slow magnetoacoustic waves were carried out in more recent works [17,18] and [19], respectively. An interesting inherent feature of these models is that the coronal loop experiencing such strong radiation not compensated by heating should disappear from the observational waveband sensitive to the hot plasma emission (e.g. 94Å, 171Å, or 193Å of the Atmospheric Imaging Assembly onboard Solar Dynamics Observatory, SDO/AIA) in 10-20 min, i.e. after 2-3 cycles of oscillations. On the other hand, for example, kink oscillations of coronal loops are often observed to reside in long-lived loops for at least several clear cycles without signatures of the loop fading, or even in the decayless regime lasting for hours with more than ten oscillation cycles [20,21]. More details on the dissipation of kink oscillations including the effects of rapid radiative cooling are outlined in detail in the recent review [8].
Accounting for a constant heating term in the coronal energy balance allowed for considering the dynamics of fast and slow magnetoacoustic waves in long-lived coronal plasma structures (see e.g. [22,23,24,25]). In this scenario (see the middle panel in Fig. 1), the role of a constant heating term reduces to maintaining the initial thermal equilibrium of the coronal plasma, neither contributing to the wave dynamics nor being affected by it. The ratio of the wave period to the characteristic timescale of the optically thin radiative cooling was shown to determine the dynamic properties of the wave in this regime.
On the other hand, taking the link between Figure 1. A schematic illustration of three different scenarios for a small-amplitude compressive wave to evolve in a coronal plasma structure affected by local heating and/or cooling processes. Left: A wave in a coronal plasma structure continuously cooling by optically thin radiation and perhaps thermal conduction, with no local thermal equilibrium. Middle: A wave in a long-lived coronal plasma structure under thermal equilibrium with ρ 0 and T 0 , maintained by a constant heating process H 0 ≡ L 0 (ρ 0 , T 0 ). The wave-caused perturbations of the local plasma parameters affect the plasma cooling rate as L(ρ 0 + ρ 1 , T 0 + T 1 ) ≈ L 0 + L ρ ρ 1 + L T T 1 and do not affect the constant heating process. Right: A long-lived coronal plasma structure hosting a wave, where both heating and cooling processes are functions of the local plasma parameters, L(ρ, T ) and H(ρ, T, B), balancing each other in the equilibrium H 0 (ρ 0 , T 0 , B 0 ) = L 0 (ρ 0 , T 0 ). In this case, both the heating and cooling rates get modified by the wave-induced perturbations as , which leads to the phenomenon of wave-induced heating/cooling (thermal) misbalance.
the properties of coronal loops and the local plasma parameters into account, the coronal heating rate is often modelled as a function of the local plasma density and temperature (see e.g. [26,27,28,29]) and local value of the magnetic field strength [30]. Moreover, the dependences of the coronal heating and radiative cooling rates on the background plasma parameters are likely to be different. Hence, perturbations of the plasma parameters caused by essentially compressive (for example, slow magnetoacoustic) waves disturb not only the mechanical equilibrium of the coronal plasma, but also violate its thermal balance through the modification of both the heating and cooling rates, thus leading to the phenomenon of a wave-induced heating/cooling (thermal) misbalance (see the righthand panel of Fig. 1). Effects of the perturbed thermal equilibrium have been traditionally considered in the context of wave propagation and stability in the interstellar medium and molecular clouds (see e.g. [31,32] and more recent works [33,34]); electric discharges, lasers, and chemically-active media (e.g. [35]); and as mechanisms for formation and evolution of rapid condensations in the solar corona such as prominences and coronal rain (see e.g. [27,36,37] and a recent review [38]). In a series of more recent works, the phenomenon of local thermal misbalance was understood to strongly affect the dynamics and stability of essentially compressive slow magnetoacoustic waves in the solar corona too. Indeed, being abundantly present in the corona (see [39,40,7], for reviews), slow waves were shown to experience a back-reaction from the perturbed thermal equilibrium, either losing or gaining energy from the coronal heating source [41], which makes the continuously heated and cooling corona an active medium for magnetoacoustic waves (cf. burning gases, lasers and masers). Accounting for this intrinsic feature of the corona allowed for revealing a number of important properties of slow waves that were missing in previous models. These newly revealed properties include new mechanisms for enhanced and frequencydependent damping of slow coronal waves [42,43,44]; modification of the phase behaviour of propagating slow waves [45]; a new, not associated with the plasma non-uniformity mechanism, dispersion of slow waves [46,47]; formation of self-sustained nonlinear periodic wave structures and autosolitons (e.g. [48,49,50]); and a new tool for probing the unknown coronal heating function by slow waves [51].
In this review, we will address some of these recent findings, focussing on the description of slow coronal waves in terms of a thin flux tube model with thermal misbalance (Sec. 2), stability of slow and entropy waves in a thermodynamically active plasma of the solar corona and its implications for diagnostics of the coronal heating function (Sec. 3), and dispersion of slow waves, caused by thermal misbalance, and formation of quasi-periodic slow wave trains (Sec. 4). A digest of the results presented and future prospects are outlined in Sec. 5.

Thin flux tube model with thermal misbalance
The dynamics of long-wavelength axisymmetric magnetoacoustic (MA) waves in solar coronal loops can be described in terms of a so-called thin flux tube approximation [52]. This approach allows one to reduce the full set of MHD equations in cylindrical coordinates to a one-dimensional form, using the second-order Taylor expansion of the variables with respect to the radial coordinate and treating the ratio of the loop radius to the characteristic wavelength as a small parameter. Thus, assuming an untwisted and non-rotating flux tube stretched along the z-axis, the set of linearised equations describing the dynamics of MA waves can be written as In Eqs. (1)-(6), subscripts "0" and "1" indicate the unperturbed value and perturbation of the corresponding variable, respectively, ρ is the plasma density, T is the temperature, and P is the gas pressure. Also, u and B are the plasma velocity and magnetic field along the flux tube, v r is the radial derivative of the radial velocity taken at r = 0, A 0 is the flux tube cross-section area. We denote k B for the Boltzmann constant, C V for the specific heat capacity, and m for the mean particle mass (see also Eq. (13) for the set of typical coronal plasma parameters). The set of Eqs. (1)-(6) is similar to that used by Ref. [47]. The non-adiabatic processes of coronal heating H(ρ, T, B), optically thin radiative cooling L(ρ, T ), and field-aligned thermal conduction with the characteristic coefficient κ are accounted for on the RHS of energy Eq. (3). More specifically, we describe the interplay between the coronal heating and radiative cooling rates through the net heat-loss function Q(ρ, T, B) = L(ρ, T )−H(ρ, T, B), measured in W kg −1 (or erg g −1 s −1 ). In the equilibrium, it is equal to zero, Q(ρ 0 , T 0 , B 0 ) = 0, providing the local thermal balance. After a small-amplitude wave-induced perturbation, the coronal heat-loss function can be Taylorexpanded as Q(ρ, T, B) ≈ Q ρ ρ 1 +Q T T 1 +Q B B 1 , where the partial derivatives Q ρ = ∂Q/∂ρ, Q T = ∂Q/∂T , and Q B = ∂Q/∂B are evaluated at the initial equilibrium. The unperturbed loop is assumed isothermal along the axis, which is typical for the coronal part of solar atmospheric flux tubes.
It is important to stress that the heat-loss derivatives Q ρ , Q T , and Q B could be considered as effective transport coefficients of the coronal plasma, carrying the information about the coronal heating process. Indeed, using some a priori prescribed model for optically thin radiation per unit mass, L(ρ, T ) ∝ ρT α (for radiation per unit volume, it is ∝ ρ 2 T α ) and representing the unknown coronal heating function in a generic form H(ρ, T, B) ∝ ρ a T b B c with the power-law indices a, b, and c being free dimensionless parameters, those coefficients Q ρ , Q T , and Q B become As one can see from Eqs. (7)- (9), the coefficients Q ρ , Q T , and Q B and hence their effect on the wave dynamics depend on the unknown heating parameters a, b, and c and, on the other hand, are highly sensitive to the value and local gradients of the unperturbed radiative cooling function L 0 . In Fig. 2, we illustrate the dependence of the coronal radiative losses L 0 on the plasma temperature, as it is predicted by the Rosner-Tucker-Vaiana [26], Klimchuk-Raymond [53], and CHIANTI [54] models. Despite the similarity in a general tendency of the function L 0 to vary with temperature, the considered models have clearly different local gradients of L 0 . In this work, we use the most recent and therefore accurate version 10.  (10) with the following differential operators, Optically thin radiative losses per unit mass as function of the plasma temperature, obtained from the CHIANTI v. 10.1 atomic database [54] for the coronal abundance and number density n 0 = 10 9 cm −3 (the red solid line); Klimchuk-Raymond model ( [53], the blue dashed line); and Rosner-Tucker-Vaiana model ( [26], the green dotted line).
where c A = B 0 / √ 4πρ 0 and c s = γk B T 0 /m are the standard Alfvén and sound (with the adiabatic index γ = 5/3) speeds in the plasma.
Equation (10) is seen to be of the fifth-order in time, so that it describes the dynamics of two fast MA modes, two slow MA modes, and one entropy mode. However, as any change in the external (outside the flux tube) pressure is neglected on the RHS of Eq. (5) (see e.g. [55]), we focus on the dynamics of essentially compressive slow MA and entropy waves, which propagate/evolve inside the loop, always in a trapped regime. For the use of this model for an adequate description of fast MA waves, the interaction of the flux tube with the external medium should be taken into account (see e.g. [56]).
In the case of weak thermal conduction κ → 0 and Q B → 0 (no dependence of the coronal heating rate on the local magnetic field strength), Eq. (10) coincides with that derived in Ref. [47]. Likewise, neglecting the effects of finite loop width (A 0 k 2 → 0), Eq. (10) reduces to the model of slow MA waves affected by thermal misbalance and thermal conduction, considered by Ref. [43]. If one additionally assumes the loop's magnetic field is infinitely strong with c A c s , Eq. (10) yields a so-called infinite field approximation, in which slow MA waves do not perturb the magnetic field and its role effectively reduces to determining the propagation direction and 1D nature of the wave (see e.g. [42,46,57]).

Stability of acoustic and thermal modes and coronal heating function
In a highly magnetised plasma with c A c s typical for the solar corona, the dynamics of slow MA waves was shown to be insensitive to the dependence of the unknown coronal heating function on the local magnetic field strength (see Fig. 2 in Ref. [43]). This reduces the parametrisation of the coronal heating function through the local plasma parameters to H(ρ, T ) ∝ ρ a T b , with two unknowns a and b. In this low-β regime, Ref. [51] obtained the following conditions for thermal (entropy) and slow MA modes described by Eq. (10) to remain stable, respectively.
Here, τ c (k) = ρ 0 C V k −2 /κ is the wavelength-dependent characteristic time of the parallel thermal conduction, and the radiative loss function L(ρ, T ) is illustrated in Fig. 2. Equations (11)- (12) demonstrate that the term associated with the field-aligned thermal conduction is always positive so it tends to stabilise the perturbation and restore the initial equilibrium, while the terms associated with thermal misbalance (namely, the heating power-law indices a and b, and the local gradient of the radiative loss function ∂L 0 /∂T ) as well as their combinations could be both positive and negative. The latter implies that for different coronal plasma conditions and different coronal heating models, the phenomenon of thermal misbalance may lead to a number of different scenarios for the initial compressive perturbation to evolve in the solar corona, with the acoustic and thermal modes being stable (decaying) and/or unstable.
The regions of the heating power-indices a and b, corresponding to stable/unstable behaviour of the acoustic and thermal modes, prescribed by Eqs. (11)- (12), are demonstrated in Fig. 3 for the following combination of the plasma parameters common for many typical wave-hosting structures in the solar corona, Temperature, T 0 = 1 MK, Number density, n 0 = 10 9 cm −3 , Magnetic field, B 0 = 4 G (β ≈ 0.2), 40 G (β → 0), Parallel thermal conductivity, κ = 10 −11 T 5/2 0 W m −1 K −1 , Mean particle mass, m = 0.6 × 1.67 × 10 −27 kg, Adiabatic index, γ = 5/3, which gives the standard sound speed c s ≈ 152 km s −1 and the corresponding acoustic period P A = λ/c s ≈ 5.5 min, used as characteristic properties of slow waves in the ideal plasma for normalisation. The grey-shaded region in Fig. 3 shows values of the heating power-indices a and b, for which both the effect of thermal misbalance and parallel thermal conduction lead to the decay of thermal and acoustic modes. As such, the stability of the coronal plasma is wavelength-independent in this domain. In contrast, for the values of a and b outside the grey-shaded region in Fig. 3, thermal misbalance leads to an effective gain of energy by the perturbation from the coronal heating source, counteracting the damping by thermal conduction or even causing the amplification of thermal and acoustic modes. In this region, the stability of the coronal plasma depends on the wavelength of the perturbation. Thus, for some guessed heating model with a and b, there is a critical wavelength λ F below which thermal conduction is strong enough to suppress the instability of these harmonics, caused by thermal misbalance.
For longer-wavelength perturbations, thermal conduction is too weak to restrain the effect of thermal misbalance and thus the coronal plasma gets unstable to these perturbations. In particular, for typical coronal parameters (13), the RTV heating models [26] are seen to be stable to the perturbations shorter than ∼100 Mm. From Eq. (11), the critical wavelength at which thermal conduction becomes insufficient to sustain the instability of thermal (entropy) mode is which could be understood as a thermal Field's length generalised for a non-constant coronal heating rate (cf. Eq. (1) and its discussion in [38]). Likewise, using Eq. (12), one can obtain the wavelength at which the  (13) with β → 0, affected by thermal conduction and thermal misbalance with the unknown coronal heating function H ∝ ρ a T b (the dependence of the heating rate on the magnetic field has no effect on the wave dynamics in the considered low-β regime). For the heating models from the grey-shaded region, both thermal and acoustic modes are stable independently of the wavelength. Outside this region, the corona is stable to the perturbations shorter than the thermal and acoustic Field's lengths (14)- (15) shown by dashed and solid lines, respectively. The red triangles show the RTV heating models [26]. The green symbols show the heating models used for illustrations in Sec. 4.
acoustic mode grows, which can be referred to as an acoustic Field's length.

Non-waveguide dispersion of slow magnetoacoustic waves
Searching for the solution of Eq. (10) in a harmonic form ∝ e i(kz−ωt) with the cyclic frequency ω and the wavenumber k yields the following dispersion relation, where the coefficients are In terms of the considered model, the main sources of the wave dispersion described by Eq. (16) where ω R and ω I are real and imaginary parts of the complex cyclic frequency κ is the wavelengthdependent characteristic time of the field-aligned thermal conduction, and are the characteristic timescales of thermal misbalance, obtained through the temperature-derivatives of the coronal heat-loss function, taken at constant gas and magnetic pressures, respectively. For typical coronal plasma conditions, the timescales τ 1,2 were shown to be about several minutes, thus coinciding with periods and damping times of slow MA waves observed in the solar corona [51]. Solution (17)- (18) was obtained by Ref. [43] in the limit of weak non-adiabaticity, i.e. treating all non-adiabatic processes slow in comparison with the wave period, providing ωτ i 1 and ω R ω I . Similar solutions in the regime of strong nonadiabaticity with ωτ i 1 could also be found in Ref. [43].
In reality, slow MA waves are usually observed in the solar corona to damp rapidly, within a few oscillation cycles [39,40,7], which is equivalent to ω R ∼ ω I . Hence, strictly speaking, neither weak nor strong limits of non-adiabaticity are applicable. For ω R ∼ ω I , Refs. [46,47] demonstrated that the phenomenon of thermal misbalance leads to an effective dispersion of slow MA waves, manifested through the dependence of the wave propagation speed on the wave frequency. Thus, in Fig. 4, we illustrate the dependence of the slow MA wave phase speed V ph = ω R /k on the wave period P = 2π/ω R , obtained from the full solution of dispersion relation (16). More specifically, the dispersion curves shown in Fig. 4 are calculated for κ → 0, which allows for a direct comparison of the dispersive effects that slow MA waves experience from the coronal plasma waveguide and the process of wave-induced thermal misbalance. We adopt the set of coronal plasma parameters given by Eq. (13), the cooling rate obtained from the CHIANTI database (Fig. 2), and the heating rate in the form of a powerlaw function H ∝ ρ a T b B c . In order to satisfy stability conditions (11)-(12), we take the heating model with a = 2 and b = −3.5 (see Fig. 3), and consider three different dependences on the magnetic field, with c = 1, 0, and −1.
The left-hand panel of Fig. 4 shows the dependence of the phase speed V ph on the wave period P , obtained from Eq. (16) in the infinite magnetic field approximation of slow MA waves, with c A c s . In this regime, slow MA waves propagate strictly along the waveguide axis without perturbing its boundaries. In other words, the tube speed c T , which is the longperiod limit of the phase speed in the ideal plasma (without thermal misbalance), tends to the sound speed c s , making the effect of the geometrical dispersion on the dynamics of slow waves negligible. Thus, the dominant mechanism of slow-wave dispersion in this regime is thermal misbalance, due to which the phase speed varies from a short-period value c s to the long-period value c sQ , prescribed by the properties of the heating/cooling processes, plasma temperature, and density as It is worth mentioning that, in general, c sQ may be either greater or lower than c s . Also, the dependence of the heating rate on the magnetic field strength has no effect on the slow-wave phase speed in this low-β regime, which is consistent with Ref. [43]. The dependence of the phase speed on the slowwave period in the regime of finite plasma-β is shown in the right-hand panel of Fig. 4. In this case, the geometrical dispersion of slow waves appears. In the ideal plasma, it leads to the variation of the phase speed from c s at short periods to c T at long periods. However, for the considered value of β ≈ 0.2 (see Eq. (13)), c T differs from c s by several percent only, which is difficult to detect in observations. Hence, the waveguiding dispersion is usually thought to be weakly important for slow MA waves in the solar corona. However, accounting for the effect of thermal misbalance leads to stronger departure of the slow-wave phase speed from the standard c s towards a new value c TQ at longer periods, determined by the properties of heating/cooling processes and  (13), with and without thermal misbalance. The colour scheme is the same for both panels. The dispersion curves with thermal misbalance are calculated for three different heating models, in which the dependence on density and temperature is taken from the stability analysis (Fig. 3), and the dependences on the magnetic field are chosen arbitrarily for illustration. The long-period phase speed limits c sQ and c TQ are given by Eqs. (20) and (21). The wave period P M at which the misbalance-caused dispersion is most effective is given by Eq. (22). See also Ref. [47], for the case with H ∝ ρ 0.5 T −3.5 B 0 . equilibrium plasma parameters as For example, Ref. [47] demonstrated that the relative error in the seismological estimation of the coronal magnetic field by slow waves (see e.g. [58,59]), could reach up to 40% if one neglects the effect of misbalance and uses the standard tube speed c T instead of c TQ . In contrast to c sQ (20), the value of c TQ (21) is shown to be sensitive to the dependence of the coronal heating function on the magnetic field strength through the heat-loss derivative Q B (9). The characteristic wave period at which the effect of dispersion caused by thermal misbalance is maximum was determined by Ref. [46] in the infinite field approximation (β → 0 and A 0 k 2 → 0) through the maximum gradient of the phase speed, as where τ 1,2 are the timescales of thermal misbalance, given by Eq. (19). Moreover, the dependence of the heating rate on the magnetic field strength is seen to become important in the right-hand panel of Fig. 4 at wave periods comparable to or greater than P M . Typical values of the wave period P M are shown in Fig. 5, for broad intervals of coronal plasma densities and temperatures, and β → 0. In particular, in the quiescent corona and coronal holes with typical densities ≈(0.5-1)×10 9 cm −3 and temperatures around 1 MK, the wave period P M is seen to be from 20 min to 70 min, which coincides with observations of very long-period propagating compressive waves in coronal plumes and inter-plume regions [60,61]. The evolution of an essentially broadband velocity pulse of a slow-mode MA nature, affected by the above-described dispersion effects, is shown in Fig. 6, obtained from a numerical solution of Eq. (10) with A 0 k 2 → 0 and c A c s and for plasma parameters (13). Namely, the left-hand panel shows the regime with the heating model H ∝ ρ 2 T −3.5 , for which both processes of thermal conduction and heating/cooling misbalance lead to the damping of all slow-mode (and entropy) harmonics constituting the initial perturbation (see Fig. 3). In this case, the initially localised pulse is seen to rapidly decrease in the amplitude and disperse (broaden) in time and space, due to faster propagation and more effective damping of shorterperiod slow MA harmonics. In the right-hand panel of Fig. 6, we demonstrate the evolution of the same initially localised velocity pulse, using the heating model H ∝ ρ 2 T −1 , for which all entropy-mode harmonics decay, while for slow MA harmonics there is a critical wavelength λ F ≈ 170 Mm determined by Eq. (15), above which slow waves get overstable through the effective gain of energy from the coronal heating source (see Fig. 3). Thus, at the initial stage of the pulse evolution, its amplitude is seen to decrease by conductive damping of the harmonics shorter than λ F ≈ 170 Mm, followed by the amplification of longer-wavelength harmonics and formation of longperiod quasi-periodic patterns by the effect of thermal misbalance. The dominant period of these slowpropagating quasi-periodic wave trains is seen to be about P M (5), which is approximately 26 min for the combination of plasma parameters (13). More examples with κ → 0 can be found in Ref. [57].

Summary and prospects
We discussed the phenomenon of a wave-induced heating/cooling (thermal) misbalance, which occurs in a continuously heated and cooling plasma of the solar corona due to the violation of a local thermal balance by compressive magnetoacoustic waves. Acting as an additional natural mechanism for the exchange of energy between the plasma and the wave, the phenomenon of thermal misbalance makes the corona an active medium for magnetoacoustic waves. The presented results and physical effects that the wave experiences as a back-reaction of this perturbed thermal equilibrium could be summarised as • Both coronal heating and cooling processes are perturbed by any compressive perturbations (e.g. magnetoacoustic waves), causing a misbalance between heating and cooling rates. • The perturbed thermal equilibrium is important not only for traditional problems of the formation and evolution of prominences and coronal rain, but also for modelling and interpretation of MHD waves in the corona. • For typical coronal conditions, the characteristic timescales of thermal misbalance are demonstrated to be about observed oscillation periods of slowmode waves, i.e. about several minutes. • The back-reaction that the wave experiences from the phenomenon of thermal misbalance includes dispersion (not connected with the waveguiding or cutoff effects traditionally considered in the corona) and enhanced frequency-dependent damping or amplification (via the energy exchange between the plasma and the wave).
• The stability of the coronal plasma structures is sensitive to the heating model and the wavelength of compressive perturbations. Requiring a long-lived (stable) corona, one can use this effect for probing heating functions in various coronal structures. For example, for a specific set of coronal plasma parameters, the RTV heating models were shown to be unstable to the perturbations longer than ∼100 Mm. • The misbalance-caused dispersion of slow magnetoacoustic waves is most pronounced in the longerperiod part of the spectrum, leading to the modification of characteristic sound and tube speeds. In some regimes of misbalance, this may lead to the formation of long-period slow-propagating wave trains from initially aperiodic, broadband perturbations.
Below, we also list several potentially interesting and, in our opinion, promising avenues for future development of this research field • Given the importance of the thermal Field's length in the physics of solar prominences and coronal rain [38], its generalisation for a non-constant coronal heating rate would allow for a fresh look at this problem, from both modelling and observational points of view. Likewise, practical implications of the acoustic Field's length for the dynamics of slowmode waves in the corona are to be revealed.
• Search for theoretically predicted slow-propagating quasi-periodic compressive wave trains in observations, using Refs. [60,61] as a starting point.
• Development of a 2D theory of slow magnetoacoustic waves in coronal plasma loops with thermal misbalance, with a focus on the transverse fine structuring of the loop.
• Revealing the effects and implications of thermal misbalance on other MHD eigenmodes of coronal plasma structures, such as fast magnetoacosutic waves and torsional Alfvén waves. The ground for this has been recently seeded in, for example, Ref. [62], where nonlinear shear Alfvén waves in a plasma with thermal misbalance were considered.
• Development and applications of the theory of thermal misbalance, accounting for the effects of parallel non-uniformity of the coronal plasma, i.e. ρ(z), T (z), and B(z), and non-locality of the coronal heating function, H(∂ρ/∂z, ∂T /∂z, ∂B/∂z). An important question here is the possible modification of the acoustic cutoff frequency.
• Accounting for the link between the coronal heating function and global parameters of the loop, such as loop length and transit time, and parameters of the photospheric driver (see e.g. Table 1 in [3]). V max V max Figure 6. Evolution of a broadband velocity pulse, obtained from a numerical solution of Eq. (10) with A 0 k 2 → 0 and c A cs, plasma parameters (13) with β → 0, and heating models H ∝ ρ 2 T −3.5 (left) and H ∝ ρ 2 T −1 (right). Different colours show the same signal measured at different distances z from the site of initial perturbation. The vertical dotted lines show the expected locations of the pulse in the ideal plasma, i.e. in the absence of the slow wave dispersion caused by thermal misbalance. The time and z-coordinate are normalised to the acoustic period P A ≈ 5.5 min and wavelength λ = 50 Mm, determined for plasma parameters (13) in the ideal plasma. The signal's amplitude is normalised to the amplitude of the initial perturbation, Vmax.
• An interesting future study could also address the effect of thermal misbalance on the nonlinear cascade and shock wave dynamics in the corona. In particular, the results obtained in Ref. [63] could be re-considered through the prism of formation of quasi-periodic trains of magnetoacoustic shocks in the regime of wave amplification caused by the misbalance [50,64].