Free Access
Issue
A&A
Volume 659, March 2022
Article Number A73
Number of page(s) 8
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/202141945
Published online 08 March 2022

© ESO 2022

1. Introduction

The outer atmosphere of the Sun, the solar corona, is a highly structured and dynamic, hot and rarefied plasma. The observed structures are generally thought to outline magnetic field lines, appearing as closed magnetic elements, known as coronal loops in active regions, and as plumes in expanding open field lines. Still, the observed fine field-aligned density filamentation of the coronal plasma (e.g., Raymond et al. 2014; Williams et al. 2020) remains a puzzle, being linked to the detailed mechanism of the coronal heating problem. Magnetohydrodynamic (MHD) waves are omnipresent and routinely observed in these structures, both standing and propagating (e.g., Nakariakov & Kolotkov 2020). These waves are of great scientific interest as they might play a role in coronal heating (e.g., Arregui 2015; Van Doorsselaere et al. 2020b). On the other hand, properties of the observed waves, such as wave speeds, periods, and so on, allow us to infer properties of the medium in which they propagate, a field known as coronal seismology (e.g., De Moortel & Nakariakov 2012; Arregui 2018). Oscillating displacements of coronal loops, which are interpreted as kink waves, were the first modes used for coronal seismology and have been intensively studied since. Kink modes are observed as impulsively triggered large amplitude standing waves or small amplitude decayless standing or ubiquitously propagating waves. The strong damping of observed large amplitude standing kink waves initially constituted a puzzle, as dissipative coefficients in the corona appear to be orders of magnitude smaller than required. Resonant absorption, that is the conversion of the global kink wave energy to localized Alfvén waves, was proposed as the main mechanism to explain the damping (e.g., Ruderman & Roberts 2002; Goossens et al. 2002). Recently, evidence on the amplitude-dependent damping of standing kink waves has emerged (Goddard & Nakariakov 2016; Nechaeva et al. 2019; Arregui 2021). Resonant absorption, which is a linear theory, cannot account for this amplitude dependence. Nonlinear mechanisms, such as the Kelvin-Helmholtz instability, or the initiation of a nonlinear azimuthal cascade have been proposed to explain this dependence (Van Doorsselaere et al. 2021). Nevertheless, resonant absorption continues to represent a viable mechanism to explain at least part of the observed damping. Analytical and numerical models building upon the original cylindrical model of Zaitsev & Stepanov (1975) and Edwin & Roberts (1983) have greatly improved our understanding of kink oscillations in coronal loops, including numerous effects related to loop geometry, structuring, and time-dependent backgrounds, among others. (see, e.g., Ruderman & Erdélyi 2009; Nakariakov et al. 2021, for reviews on the subject)

Another property of kink waves is their polarization. Most of the existing works focus on the linearly polarized kink wave, even though in analytical studies a circular polarization is assumed, as the mathematical treatment is more convenient in this case. One could argue that this is sufficient, since in the linear regime any state of polarization can be constructed from either linear or circularly polarized waves. However, some properties of the circularly polarized kink waves are not obvious, such as their resonant absorption. Additionally, the nonlinear evolution of circularly polarized kink waves might be different from their linearly polarized counterpart. In observations, standing kink waves are mostly categorized as “vertically” or “horizontally” polarized, meaning linear polarization with the displacement in or out of the plane of the oscillating loop, respectively. Nevertheless, using three-dimensional (3D) reconstructions from the STEREO EUVI/A and B spacecraft, Aschwanden (2009) found that the loop oscillation classified as a linearly polarized horizontal oscillation by Wang et al. (2008) using two-dimensional (2D) imaging, was in fact a circularly polarized oscillation. Moreover, Wang et al. (2008) found that out of the 14 loop oscillations studied, more than half of the cases were not classifiable as either horizontally or vertically polarized, which might be due to more complicated motions, for example, elliptical or circular polarization. Therefore, a significant proportion of observed kink oscillating loops could be elliptically or circularly polarized. Circularly polarized, propagating kink waves were also observed in chromospheric magnetic elements (Stangalini et al. 2017). Theoretically, a circular polarization of kink waves is expected to evolve in loops with twisted magnetic fields, even if the initial perturbation is linearly polarized (Ruderman & Terradas 2015). Although not transverse but a slow-type wave, it is worth mentioning that Jess et al. (2017) observed a circularly polarized m = 1 mode showing a spiral pattern in a sunspot.

In this paper, we study the circularly polarized kink oscillations, both standing and propagating, of a straight magnetic cylinder, for both linear and nonlinear perturbation amplitudes. The paper is organized as follows. In Sect. 2, we present the numerical method and describe the numerical models. In Sect. 3, we present the results, along with some discussion. Finally, in Sect. 4, we summarize the results and present some possibilities for future work.

2. Numerical method and model

We ran ideal 3D MHD simulations with MPI-AMRVAC1 (Porth et al. 2014; Xia et al. 2018). The second-order tvd solver was used with the woodward slope limiter, unless otherwise specified. The numerical setup at equilibrium consists of a straight cylindrical flux tube embedded in a uniform atmosphere. The background magnetic field B0 = 10 G is straight and uniform throughout the domain, and along the z-axis. The flux tube is defined by its three times higher density ρi, than that of the background, ρ0 ≈ 2.3 ⋅ 10−12 kg m−3. The inhomogeneous layer between the inside and outside of the flux tube is defined as a sinusoidal variation in density, with different widths l. Then, the density in the whole domain is defined as follows:

ρ ( r ) = { ρ i , r < a l / 2 , ρ i ( ρ i ρ 0 ) sin ( π 2 l ( r ( a l / 2 ) ) ) , a + l / 2 r a l / 2 . ρ 0 , r > a + l / 2 $$ \begin{aligned} \tiny \rho (r) = {\left\{ \begin{array}{ll}&\rho _i, r < a - l/2,\\&\rho _i - (\rho _i - \rho _0)\mathrm{sin} \left(\frac{\pi }{2 l}(r - (a-l/2)) \right), a + l/2 \le r \ge a - l/2. \\&\rho _0, r > a + l/2\\ \end{array}\right.} \end{aligned} $$(1)

Gravity and nonideal effects, such as optically thin radiative losses and thermal conduction, are not included in this model. On top of the background equilibrium, we imposed velocity and magnetic field perturbations of the kink wave solution (e.g., Zaitsev & Stepanov 1975; Edwin & Roberts 1983; Van Doorsselaere et al. 2020a) to study the circularly polarized standing kink wave. In Cartesian coordinates, the solution reads as follows:

v x = 1 ( 2 ) A 2 L π B 0 v x ( r ) sin ( π L z ) , b x = A 2 L π B 0 b x ( r ) cos ( π L z ) , v y = 1 ( 2 ) A 2 L π B 0 v y ( r ) sin ( π L z ) , b y = A 2 L π B 0 b y ( r ) cos ( π L z ) , v z = 0 , b z = 1 B 0 A b z ( r ) sin ( π L z ) , $$ \begin{aligned} \begin{array}{l} \mathbf v _x = -\frac{1}{\sqrt{(}2)} A \frac{2 L}{\pi B_0} { v}_x(r) \mathrm{sin} \left(\frac{\pi }{L} z\right) ,\\ \mathbf b _x = A \frac{2 L}{\pi B_0} b_x(r) \mathrm{cos} \left(\frac{\pi }{L} z\right),\\ \mathbf v _{ y} = \frac{1}{\sqrt{(}2)} A \frac{2 L}{\pi B_0} { v}_{ y}(r) \mathrm{sin} \left(\frac{\pi }{L} z\right),\\ \mathbf b _{ y} = A \frac{2 L}{\pi B_0} b_{ y}(r) \mathrm{cos} \left(\frac{\pi }{L} z\right),\\ \mathbf v _z = 0,\\ \mathbf b _z = \frac{1}{B_0} A b_z(r)\mathrm{sin} \left(\frac{\pi }{L} z\right) , \end{array} \end{aligned} $$(2)

where

v x ( r ) = { k i x 2 ( J 0 ( k i r ) J 2 ( k i r ) ) r 2 + 2 y 2 J 1 ( k i r ) r 3 2 J 1 ( k i a ) , r a , 2 y 2 K 1 ( k e r ) r 3 k e x 2 ( K 0 ( k e r ) + K 2 ( k e r ) ) r 2 2 K 1 ( k e a ) , r > a , v y ( r ) = { x y ( k i r ( J 0 ( k i r ) J 2 ( k i r ) ) 2 J 1 ( k i r ) ) 2 r 3 J 1 ( k i a ) , r a , x y ( 2 K 1 ( k e r ) + k e r ( K 0 ( k e r ) + K 2 ( k e r ) ) ) 2 r 3 K 1 ( k e a ) , r > a , b x ( r ) = { x y ( 2 J 1 ( k i r ) k i r ( J 0 ( k i r ) J 2 ( k i r ) ) ) 2 r 3 J 1 ( k i a ) , r a , x y ( 2 K 1 ( k e r ) + k e r ( K 0 ( k e r ) + K 2 ( k e r ) ) ) 2 r 3 K 1 ( k e a ) , r > a , b y ( r ) = { 2 x 2 J 1 ( k i r ) r 3 + k i y 2 ( J 0 ( k i r ) J 2 ( k i r ) ) r 2 2 J 1 ( k i a ) , r a , k e y 2 ( K 0 ( k e r ) + K 2 ( k e r ) ) r 2 2 x 2 K 1 ( k e r ) r 3 2 K 1 ( k e a ) , r > a , b z ( r ) = { x r J 1 ( k i r ) J 1 ( k i a ) , r a , x r K 1 ( k e r ) K 1 ( k e a ) , r > a , $$ \begin{aligned} \begin{array}{l} { v}_x(r) = {\left\{ \begin{array}{ll} \frac{\frac{k_i x^2 \left(J_0\left(k_i r \right)-J_2\left(k_i r \right)\right)}{r^2}+\frac{2 { y}^{\prime 2} J_1\left(k_i r \right)}{r^3}}{2 J_1(k_i a)}, r\le a,\\ \frac{\frac{2 { y}^{\prime 2} K_1\left(k_e r \right)}{r^3}-\frac{k_e x^2 \left(K_0\left(k_e r\right)+K_2\left(k_e r \right)\right)}{r^2}}{2 K_1(k_e a)}, r> a,\\ \end{array}\right.} \\ { v}_{ y}(r) = {\left\{ \begin{array}{ll} \frac{x { y}^{\prime } \left(k_i r \left(J_0\left(k_i r \right)-J_2\left(k_i r \right)\right)-2 J_1\left(k_i r \right)\right)}{2 r^3 J_1(k_i a)}, r \le a,\\ -\frac{x { y}^{\prime } \left(2 K_1\left(k_e r \right)+k_e r \left(K_0\left(k_e r \right)+K_2\left(k_e r \right)\right)\right)}{2 r^3 K_1(k_e a)}, r > a,\\ \end{array}\right.} \\ b_x(r) = {\left\{ \begin{array}{ll} \frac{x { y}^{\prime } \left(2 J_1\left(k_i r \right)- k_i r \left(J_0\left(k_i r \right)-J_2\left(k_i r \right)\right)\right)}{2 r^3 J_1(k_i a)}, r\le a,\\ \frac{x { y}^{\prime } \left(2 K_1\left(k_e r \right)+ k_e r \left(K_0\left(k_e r \right)+K_2\left(k_e r \right)\right)\right)}{2 r^3 K_1(k_e a)}, r > a,\\ \end{array}\right.} \\ b_{ y}(r) = {\left\{ \begin{array}{ll} -\frac{\frac{2 x^2 J_1\left(k_i r\right)}{r^3}+\frac{k_i { y}^{\prime 2} \left(J_0\left(k_i r \right)-J_2\left(k_i r \right)\right)}{r^2}}{2 J_1(k_i a)}, r \le a,\\ \frac{\frac{k_e { y}^{\prime 2} \left(K_0\left(k_e r \right)+K_2\left(k_e r \right)\right)}{r^2}-\frac{2 x^2 K_1\left(k_e r \right)}{r^3}}{2 K_1(k_e a)}, r > a,\\ \end{array}\right.} \\ b_z(r) = {\left\{ \begin{array}{ll} \frac{x}{r} \frac{J_1\left(k_i r \right)}{J_1(k_i a)}, r \le a,\\ \frac{x}{r} \frac{K_1\left(k_e r \right)}{K_1(k_e a)}, r > a,\\ \end{array}\right.} \end{array} \end{aligned} $$(3)

with r = x 2 + y 2 $ r = \sqrt{x^2 + {\mathit{y}^\prime}^2} $ and a = 1 Mm is the radius of the tube. We note that y′ = 0 denotes the center of the displaced flux tube along x = 0 and not the center of the flux tube in equilibrium. The initial displacement of the flux tube is required as a tube undergoing a circularly polarized kink wave does not cross the equilibrium position. Furthermore, A is the amplitude of the perturbation in units of pressure. Here, we use units in which the magnetic permeability μ = 1. The displacement is given by

y = y A 2 L 2 a π 2 B 0 2 sin ( π L z ) . $$ \begin{aligned} { y}^{\prime } = { y} - A \frac{2 L^2}{a \pi ^2 B_0^2}\mathrm{sin} \left(\frac{\pi }{L} z\right). \end{aligned} $$(4)

We note that while the displacement is in the y-direction, the initial velocity perturbation inside the flux tube is in the negative x-direction. The coefficients ki = ke ≈ 0.0112 are the fundamental radial wave numbers. Furthermore, Ji is the Bessel function of first kind and Ki is the modified Bessel function of second kind; L = 200 Mm is the length of the flux tube from foot point to foot point; A is the amplitude of the initial perturbation, in units of pressure; and vz was neglected due to a sufficiently low plasma beta, β ≈ 0.2 (Yuan & Van Doorsselaere 2016; Van Doorsselaere et al. 2020a).

For the study of the circularly polarized standing kink wave, the symmetrical properties of the standing mode were exploited in order to halve the computational costs. In this sense, the spatial extents are [ − 3.5a, 3.5a]2 × [0, L/2], where L/2 denotes the position of the antinode of the fundamental kink. The same spatial extent was used also for studying propagating waves. The boundary conditions were set “open” or continuous, zero divergence laterally, while the foot points were fixed by forcing the velocity components to be antisymmetric at the lower z-axis boundary. At the top z-axis boundary, the fundamental mode was mirrored by setting vz, bx, by antisymmetric, and the other variables symmetric. Simulations were run for both the standing and propagating mode. For propagating waves, at the bottom z-axis boundary, we employed a driver of circularly polarized kink waves, by setting the following:

v x ( t ) = A sin ( 2 π P t ) , v y ( t ) = A cos ( 2 π P t ) , $$ \begin{aligned} { v}_x(t) = A\, \mathrm{sin} \left(\frac{2 \pi }{P} t \right), \quad { v}_{ y}(t) = A\, \mathrm{cos} \left(\frac{2 \pi }{P} t \right), \end{aligned} $$(5)

where P ≈ 86 s is the wave period. We note that for propagating waves, the driven wave period is shorter than the fundamental standing mode period in order for the kink wave to evolve significantly by the time it reaches the top z-axis boundary. To simulate propagating waves, the top z-axis boundary was set to be continuous or having zero divergence.

The numerical domain consists of 2562 × 64 cells of uniform resolution, with fewer cells along the z-direction, along which the solution is expected to be smooth. We conducted a series of convergence studies with both lower and higher resolution, and qualitative differences can be observed for high-amplitude runs which show nonlinear small-scale generation, as expected. Stemming from the piecewise nature of the kink solution in Eq. (3) combined with the discretized numerical grid, a nonzero magnetic field divergence was generated at the flux tube boundary as an initial condition. Therefore, in order to eliminate the initial divergence and to maintain it at low levels, we employed the multigrid method of Teunissen & Keppens (2019), which uses a projection scheme to remove the divergence part of the magnetic field.

3. Results and discussion

The simulations evolve for ≈70 min, or ≈4.5 oscillation periods in the case of standing waves. We ran simulations for both standing and propagating waves, which are both circularly and linearly polarized, and also low (linear) and high (nonlinear) amplitude versions. Therefore, comparisons of damping times, oscillations periods, transverse-wave induced Kelvin-Helmholtz (TWIKH) growth rates, among others, can be made between the circularly and linearly polarized setups. First we present the results of low amplitude standing wave simulations. We set the amplitude to A = 3.17 μPa, which results in a maximal displacement of the loop of roughly 0.025a. Kink oscillations of coronal loops can be considered to evolve in the linear regime when the ratio of maximal displacement to the loop radius is much less than unity (Ruderman & Goossens 2014), which is satisfied in this case. The evolution of the circularly polarized standing kink wave is presented2 in Fig. 1.

thumbnail Fig. 1.

Kinetic energy and velocity vectors in the cross-sectional slice at the midpoint of the flux tube for the circularly polarized standing kink mode after approximately two fundamental kink mode periods of evolution. Axis units are 10 Mm. Kinetic energy is in code units. The size of the velocity vectors is proportional to their magnitude. An animation of this figure is available online.

In order to emphasize the dynamics in the inhomogeneous layer (where resonant absorption is taking place) for the simulation in Fig. 1, a thick inhomogeneous layer (l = a) combined with a diffusive, hll solver was employed. The flow of energy from the core region of the tube to the inhomogeneous layer, also called resonant absorption, is clearly visible in Fig. 1. For comparison, a simulation of the well-known case of linear polarization with the same parameters is shown in Fig. 2.

thumbnail Fig. 2.

Same as in Fig. 1, but for a linearly polarized standing kink mode. An animation of this figure is available online.

In these figures, resonant absorption can be identified as the kinetic energy growth in the inhomogeneous layer, around the resonant position, where the kink speed equals the local Alfvén speed. At the same time, the kinetic energy in the core region, which is associated with the kink wave, is decreasing. Resonant absorption is thus a mechanism to transfer energy from the global, large-scale transverse kink oscillation of the flux tube to local azimuthal flows around it. The differences in evolution between the two cases are obvious: while in the linearly polarized case the usual m  =  1 torsional Alfvén wave pattern is emerging, in the circularly polarized case the resonance leads to spiral patterns. The development of such a spiral appearance is a consequence of phase mixing in the inhomogeneous layer (Soler & Terradas 2015), as different parts of the kink wave propagate at different speeds within the inhomogeneous layer. Perturbations toward the outer edge of the flux tube overtake the inner edge ones, which in a circularly polarized case results in an increasingly wound up spiral over time. This is in analogy with the linearly polarized case, where over time it leads to a wavy velocity perturbation pattern around the resonant layer (e.g., Poedts et al. 1990; Soler & Terradas 2015). Although the dynamics of resonant absorption look different, we find that the resonant damping of the kink wave proceeds at the same rate, irrespective of the polarization, which just confirms the expectation from analytics. We verified this result for three different inhomogeneous layer thicknesses. As in previous studies (Magyar & Van Doorsselaere 2016), we find that for low amplitudes, theoretical damping rates due to resonant absorption (Ruderman & Roberts 2002; Goossens et al. 2002) constrain the measured damping well. The measured oscillation period of ≈16.2 min is within 5 percent of the analytically predicted one, and independent of the polarization state, as expected.

In order to determine whether circularly polarized kink waves are distinguishable from linearly polarized ones through Doppler shift measurements, we forward modeled the simulations, obtaining the emission from the data cubes corresponding to the AIA 171 Å channel, shown in Fig. 3.

thumbnail Fig. 3.

Time-distance plots of Doppler signal integrated through a slit along the y-axis at the midpoint of the flux tube, with the LOS the positive x-axis. Left: linearly polarized kink. Right: circularly polarized kink. Doppler shift is in code units. Length and time units are ≈30 km and 42 s, respectively.

It can be observed that despite resonant absorption appearing differently for a circularly polarized kink wave than for a linearly polarized one, this is not detectable in the Doppler signal. The only tell-tale sign of the circular polarization, albeit small due to the small displacement amplitude, is the alternating vertical displacement of the Doppler signals, which is absent in the linearly polarized case. We note, however, that while the Doppler signal is independent of the line-of-sight (LOS) for the circularly polarized kink, it varies for the linearly polarized one. Some additional discussion on the detectability of circularly polarized kink waves is presented in the next section.

The case of propagating circularly polarized kink waves in the linear regime naturally shares many similarities to the standing wave, which can be interpreted as the superposition of two counter-propagating waves. In Fig. 4 the cross-sectional kink wave dynamics for a flux tube with l = 0.25a is shown.

thumbnail Fig. 4.

Kinetic energy and velocity vectors in the cross-sectional slice at z = 50 Mm for the propagating circularly polarized kink wave. Axis units are 10 Mm. Kinetic energy is in code units. The size of the velocity vectors is proportional to their magnitude.

It is worth noting that in the propagating wave setup, the velocity perturbation that drives the kink wave at the bottom boundary of the numerical domain does not depend on position. That is, the whole boundary is driven. Therefore, while the driver excites propagating kink waves in the flux tube, in the homogeneous medium surrounding the flux tube there are also excited, propagating circularly polarized Alfvén waves. These waves propagate then in superposition through the numerical domain. The kink wave solution falls off exponentially outside the flux tube, leaving the Alfvén waves appearing dominant sufficiently far away from it. These can be seen as the velocity perturbations in the homogeneous region outside of the flux tube in Fig. 4 (readers can compare this to Fig. 1 for an estimation of the exponential scale length of the perturbation outside the flux tube). In a similar fashion as the standing wave, we compared the Doppler signals of propagating, linearly and circularly polarized kink waves. The result is very similar to the one in Fig. 3, albeit the time axis gets replaced by the z-axis, meaning that the difficulty in detecting kink wave polarization from Doppler signals is still valid for propagating waves.

In the following, the nonlinear evolution of the circularly polarized kink wave is investigated. The amplitude was set to A = 60.34 μPa, which led to a displacement of ≈0.6a, a value within the nonlinear regime. The inhomogeneous layer was set to be defined by only the diffusivity of the numerical solver, thus initially setting l = 0. This ensured that nonlinear processes at the edge of the flux tube developed rapidly (see, e.g., Terradas et al. 2018). In Fig. 5, the nonlinear dynamics are shown in the kink antinode cross section.

thumbnail Fig. 5.

Density and velocity vectors in the cross-sectional slice at the midpoint of the flux tube for the nonlinear setup after approximately two fundamental kink mode periods of evolution. Axis units are 10 Mm. Density is in code units. The size of the velocity vectors is proportional to their magnitude. An animation of this figure is available online.

As it is well known from simulations on the nonlinear evolution of the linearly polarized kink (e.g., Terradas et al. 2008; Antolin et al. 2014; Magyar et al. 2015, and many others), at the edge of the flux tube TWIKH-rolls develop. These roll-ups appear less symmetric than in a linearly polarized situation, with the asymmetry inducing deviations from circular polarization. This is illustrated in Fig. 6, where the trajectory of the flux tube midpoint is traced for the whole duration of the nonlinear simulation.

thumbnail Fig. 6.

Path line plot obtained by integrating the velocity map over the duration of the simulation in the tube cross section, showing the displacement of the flux tube midpoint over time. The color corresponds to the snapshot number, thus showing an advance in time, with “0” being the initial condition and “50” denoting the final snapshot. Axis units are 10 Mm.

We note that, initially, the trajectory is close to circular, albeit damping, leading to an inward spiral. Once TWIKH rolls develop, the kink oscillation transitions into elliptical polarization, with a more pronounced oscillation along the x-axis. Although both numerical (e.g., Terradas et al. 2018, see also above) and analytical (e.g., Barbulescu et al. 2019) works investigated the growth rate of the TWIKH rolls, these studies considered the linearly polarized case in which the shear flow is sinusoidal and it occurs at the same position with respect to the flux tube. In the case of a circularly polarized kink wave, at least in the idealized linear regime with a discontinuous boundary, the perturbation amplitude is constant, and it undergoes rotation with respect to the flux tube. Thus, from a Lagrangian perspective (i.e., following the motion of the flux tube during the oscillation), at a particular position along the flux tube boundary, the shear still undergoes a sinusoidal evolution. However, the difference compared to the linearly polarized case is that the shear is present all around the boundary within an oscillation period. This has implications on the growth rate of higher angular harmonics, specific to the development of TWIKH rolls. Moreover, the integrated energy over an oscillation period is two times higher for the same amplitude for a circularly polarized kink wave compared to a linearly polarized oscillation. Therefore, in order to test the effect of circular polarization on the growth of TWIKH rolls, the growth rate is compared to simulations with linear polarization of both the same amplitude and 2 $ \sqrt{2} $ times the amplitude. In Fig. 7, the power in high angular harmonics is shown for the linearly and circularly polarized simulations. The method of obtaining the power spectra of the tube perturbation is the following. Starting from the cross-sectional density plot in the kink antinode, a polar transformation was applied to the image. In the polar image of the cross section, an undisturbed loop appeared as a straight line at the constant radial distance r = a, while a kink wave of linear amplitude appeared as a sinusoidal curve fitting exactly one wavelength in the angular direction from 0 to 2π. Then, using smoothing and pattern identification methods, built in Mathematica3, the curve delimiting the tube boundary (r as a function of θ) was obtained, to which we applied the Fourier transform (see Fig. 8).

thumbnail Fig. 7.

Summed contribution of higher angular harmonics, summing over the Fourier components in a lower range, i.e., m = 3 5 | G ( m ) | 2 $ \sum\nolimits_{m=3}^{5} |G(m)|^2 $ (left) and higher range m = 6 15 | G ( m ) | 2 $ \sum\nolimits_{m=6}^{15} |G(m)|^2 $ (right), as a function of time, in units of ≈85.9 s (duration between snapshots).

thumbnail Fig. 8.

Steps for obtaining the boundary curve. Left: Polar transform of the density plot shown in Fig. 5, smoothed using the routine MedianFilter in Mathematica. Center: flux tube boundary identified from the polar transformed plot, using the routine MorphologicalPerimeter in Mathematica. Right: data points extracted from the flux tube boundary plot. Radius units are 10 Mm.

Smoothing removed the step-like appearance of the tube boundary due to the limited numerical resolution, and it helped in the detection of the boundary. We note that the boundary extraction was complicated by the smooth variation of the density and the generation of TWIKH rolls, which result in local density minima, and are recognized as small closed-curve patches. As the extracted boundary position as a function of the angle could not be multivalued at places where roll-ups appeared, the point with the largest radius value was chosen. This results in discontinuities appearing at these positions in the extracted curve, altering the relative contribution of highly developed roll-ups to the power spectra. Nevertheless, this technique still allows for the qualitative comparison of the growth rate of small scales for the two polarizations. For Fig. 7, we separated the higher amplitude harmonics into a lower (m ≈ 3 − 5) and upper range (m ≈ 6 − 15), which are summed up in these ranges, respectively. This is justified by the appearance of comparatively larger eddies in the circularly polarized case initially. The Fourier coefficients are defined as follows (Terradas et al. 2018):

G ( m ) = 1 N n = 0 N 1 g ( n ) exp i 2 π N n m , $$ \begin{aligned} G(m) = \frac{1}{N} \sum _{n=0}^{N-1} {g}(n) \exp ^{-i \frac{2 \pi }{N} n m}, \end{aligned} $$(6)

where N is the discrete total number of data points (n = 0, …, N − 1) extracted along the loop boundary, and g(n) is the extracted boundary curve (r as a function of θ). We define θn = 2πn/N for convenience. Then, the contribution of each mode m to the total signal is given by the inverse Fourier transform,

g ( θ n ) = m = 0 N 1 G ( m ) exp i m θ n . $$ \begin{aligned} {g}(\theta _n) = \sum _{m=0}^{N-1} G(m) \exp ^{i m \theta _n} . \end{aligned} $$(7)

The obtained spectra appear rapidly varying and noisy. In addition to the inherently noisy nature of the discrete Fourier transform (e.g., Vaughan 2005), the extraction process of the boundary curve is thought to account for some of this variability. Therefore, care should be taken when interpreting the spectral evolution presented in Fig. 7. However, it is clear that in the linearly polarized case, with the same amplitude as in the circularly polarized case, higher harmonics develop considerably slower; whereas for the linearly polarized case, with 2 $ \sqrt{2} $ times the amplitude, the spectral evolution is qualitatively similar to that for the circularly polarized case.

Not shown in Fig. 7, but visible4 in Fig. 5, is the coupling to the m = 2 fluting mode, as in the linearly polarized case (e.g., Magyar & Van Doorsselaere 2016), being especially prominent during the first oscillation period, manifesting as an area-conserving elliptical deformation of the circular cross section of the tube. A nonlinear propagating wave simulation was also carried out. In this case, the expected deformation and transition to turbulence of the unidirectionally propagating kink wave was observed (Magyar et al. 2017, 2019). The differences between the circularly polarized and the linearly polarized simulation are essentially the same as for the standing wave, in the sense that kink wave energy is a better indicator of a nonlinear growth rate than the velocity amplitude of the wave.

4. Conclusions

We have studied the circularly polarized kink wave in a cylindrical flux tube through numerical simulations. Both standing and propagating modes were studied for perturbation amplitudes, both in the linear and nonlinear regime. The most striking difference between the circularly and linearly polarized cases is the appearance of the wave perturbation once resonant absorption and phase mixing start to act. While in the linearly polarized case the resonant flow at the edge of the tube is in the direction of the transverse kink perturbation, in the circularly polarized case it is spiral-like and can be perpendicular to the kink perturbation. Although resonant absorption shows morphological differences for linear and circular polarization, no impact on the damping time or oscillation period was found, which is in agreement with analytical results. In the case of propagating waves, the same conclusions can be drawn on the differences between the circularly and linearly polarized kink waves, with a variation in time replaced by a variation over propagation distance for continuously driven waves.

In the nonlinear regime, for displacements comparable to the tube radius, a tube undergoing circularly polarized kink oscillations deforms and develops small scales or TWIKH rolls faster than a linearly polarized kink wave with the same amplitude. We find that the growth rate for nonlinear small-scale generation is governed by the energy of the oscillation, rather than the perturbation amplitude. Therefore, a growth rate comparable to the one in the circularly polarized kink was achieved with a linearly polarized kink wave with 2 $ \sqrt{2} $ times the former’s amplitude. This is also valid for nonlinear propagating kink waves.

An important point is the observability of the polarization state of kink waves in the corona. We find that despite the previously described differences between the appearance of resonant flows, there is no difference in the Doppler signal. The only tell-tale signature of circular polarization in the studied case is the displacement of alternating Doppler signals. However, we note that in the linear case, no displacement would be observed if only the LOS is parallel to the polarization direction. The Doppler signal is independent of LOS direction in the circularly polarized case. In the linearly polarized case, while the Doppler signal qualitatively remains the same for different LOS directions, the amplitude varies roughly with the cosine of the angle between LOS and polarization direction (see, e.g., Yuan & Van Doorsselaere 2016). In the near future, simultaneous observations of coronal loop oscillations by the Extreme Ultraviolet Imager on board the Solar Orbiter mission, and from SDO/AIA, could provide the opportunity to perform stereoscopic analysis with unprecedented resolution. With suitable vantage points (i.e., spacecraft separation angles ≤90°), the complete 3D structure of the loop may be deduced without the limiting assumption of semi-circularity. Then with sufficiently long-lived data, the polarization of the transverse oscillation with respect to the (modeled) plane of the loop may be calculated through comparison of the projected and apparent displacement in the two fields of view. The ubiquity of decayless kink oscillations in the solar corona implies that there will be at least some targets for stereoscopic analysis. Curiously, all known detections of decayless oscillations have been polarized in the horizontal direction, as with the large amplitude regime (Anfinogentov et al. 2015). Whether this is in fact an observational bias that is caused by the studies’ single field of view has important implications for the driver, and (apparent lack of) damping mechanisms of decayless oscillations. We note that even if all kink oscillations are initially horizontally polarized, the expansion or contraction of the hosting coronal loop (a common occurrence, e.g., Russell et al. 2015) may still provide the opportunity to observe transverse oscillations which are elliptically polarized. We suggest that observations looking for slowly expanding and contracting coronal loops exhibiting transverse kink oscillations, which were observed by both Solar Orbiter and SDO, may be used to substantiate the results presented here.

The caveats of this study are the simplifications used in modeling coronal loops. We neglect the curvature of the loops, opting for a straight flux tube. Nevertheless, Van Doorsselaere et al. (2004) demonstrated that loop curvature has a minimal effect on the kink eigenmodes and the differences between horizontal and vertical polarization are negligible, although Ruderman (2009) pointed out that this is only valid when the loop expansion (varying cross section along the loop) is neglected, as in the present study. We also note that Guo et al. (2020) have shown that any ellipticity of the loop cross section affects both the eigenfrequencies and the damping rates by resonant absorption for different polarizations. Gravity is also neglected, leading to a homogeneous density and Alfvén speed profile along the loop. The density scale height in the corona is around 50 Mm for a 1 MK plasma, leading to large variations in Alfvén speed within long loops. While the ratios of periods of different standing kink wave harmonics are affected by stratification (e.g., Andries et al. 2005), we do not expect a polarization-dependent effect. We have only studied two very specific polarization states, linear and circular. It seems improbable that various coronal mechanisms excite these specific polarization states of the kink oscillation exactly, with most oscillations probably in an elliptical polarization state. However, as elliptical polarization can be positioned between the two extremes of linear and circular polarization, our conclusions on the differences or similarities between these two states should still hold for the elliptical polarization state.

In the following, we discuss the potential applicability of this study. The polarization state of observed kink oscillations of coronal loops is often not determined, with a few exceptions noted in the Introduction. However, the polarization state could offer us important seismological information on the driving mechanism of kink oscillations, especially in the context of small amplitude decayless oscillations (e.g., Nisticò et al. 2013; Anfinogentov et al. 2015), which are thought to be continuously reinforced standing kink oscillations. In this sense, a statistical analysis on the polarization state of these oscillations could constrain some properties of the loop footpoint driver, such as variability, the amount of vorticity, and time correlation. In particular, some explanations of a decayless kink oscillation arising as a self-oscillatory process (Nakariakov et al. 2016) entail the flow inducing the oscillation is quasi-stationary, which would lead to linearly polarized oscillations.

Movies

Movie 1 associated with Fig. 1 (KE_circ) Access here

Movie 2 associated with Fig. 2 (KE_lin) Access here

Movie 3 associated with Fig. 5 (hamp_circ) Access here


2

In the online animated version.

4

In the online animated version.

Acknowledgments

N. M. was supported by a Newton International Fellowship of the Royal Society. T. D. and T. V. D. were supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 724326). T. V. D. was also supported by the C1 grant TRACEspace of Internal Funds KU Leuven. V. M. N. acknowledges the STFC consolidated grant ST/T000252/1, and the BK21 plus program through the National Research Foundation funded by the Ministry of Education of Korea.

References

  1. Andries, J., Arregui, I., & Goossens, M. 2005, ApJ, 624, L57 [NASA ADS] [CrossRef] [Google Scholar]
  2. Anfinogentov, S. A., Nakariakov, V. M., & Nisticò, G. 2015, A&A, 583, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Antolin, P., Yokoyama, T., & Van Doorsselaere, T. 2014, ApJ, 787, L22 [Google Scholar]
  4. Arregui, I. 2015, Philos, Trans. R. Soc. London Ser. A, 373, 20140261 [Google Scholar]
  5. Arregui, I. 2018, Adv. Space Res., 61, 655 [NASA ADS] [CrossRef] [Google Scholar]
  6. Arregui, I. 2021, ApJ, 915, L25 [NASA ADS] [CrossRef] [Google Scholar]
  7. Aschwanden, M. J. 2009, Space Sci. Rev., 149, 31 [CrossRef] [Google Scholar]
  8. Barbulescu, M., Ruderman, M. S., Van Doorsselaere, T., & Erdélyi, R. 2019, ApJ, 870, 108 [Google Scholar]
  9. De Moortel, I., & Nakariakov, V. M. 2012, Roy. Soc. London Philos. Trans. Ser. A, 370, 3193 [Google Scholar]
  10. Edwin, P. M., & Roberts, B. 1983, Sol. Phys., 88, 179 [Google Scholar]
  11. Goddard, C. R., & Nakariakov, V. M. 2016, A&A, 590, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Goossens, M., de Groof, A., & Andries, J. 2002, in SOLMAG 2002. Proceedings of the Magnetic Coupling of the Solar Atmosphere Euroconference, ed. H. Sawaya-Lacoste, ESA Spec. Pub., 505, 137 [NASA ADS] [Google Scholar]
  13. Guo, M., Li, B., & Van Doorsselaere, T. 2020, ApJ, 904, 116 [Google Scholar]
  14. Jess, D. B., Van Doorsselaere, T., Verth, G., et al. 2017, ApJ, 842, 59 [Google Scholar]
  15. Magyar, N., & Van Doorsselaere, T. 2016, A&A, 595, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Magyar, N., Van Doorsselaere, T., & Marcu, A. 2015, A&A, 582, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Magyar, N., Van Doorsselaere, T., & Gossens, M. 2017, Sci. Rep., 7, 14820 [NASA ADS] [CrossRef] [Google Scholar]
  18. Magyar, N., Van Doorsselaere, T., & Goossens, M. 2019, ApJ, 873, 56 [NASA ADS] [CrossRef] [Google Scholar]
  19. Nakariakov, V. M., & Kolotkov, D. Y. 2020, ARA&A, 58, 441 [Google Scholar]
  20. Nakariakov, V. M., Anfinogentov, S. A., Nisticò, G., & Lee, D. H. 2016, A&A, 591, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Nakariakov, V. M., Anfinogentov, S. A., Antolin, P., et al. 2021, Space Sci. Rev., 217, 73 [NASA ADS] [CrossRef] [Google Scholar]
  22. Nechaeva, A., Zimovets, I. V., Nakariakov, V. M., & Goddard, C. R. 2019, ApJS, 241, 31 [NASA ADS] [CrossRef] [Google Scholar]
  23. Nisticò, G., Nakariakov, V. M., & Verwichte, E. 2013, A&A, 552, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Poedts, S., Goossens, M., & Kerner, W. 1990, Comput. Phys. Commun., 59, 95 [NASA ADS] [CrossRef] [Google Scholar]
  25. Porth, O., Xia, C., Hendrix, T., Moschou, S. P., & Keppens, R. 2014, ApJS, 214, 4 [Google Scholar]
  26. Raymond, J. C., McCauley, P. I., Cranmer, S. R., & Downs, C. 2014, ApJ, 788, 152 [NASA ADS] [CrossRef] [Google Scholar]
  27. Ruderman, M. S. 2009, A&A, 506, 885 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Ruderman, M. S., & Erdélyi, R. 2009, Space Sci. Rev., 149, 199 [Google Scholar]
  29. Ruderman, M. S., & Goossens, M. 2014, Sol. Phys., 289, 1999 [NASA ADS] [CrossRef] [Google Scholar]
  30. Ruderman, M. S., & Roberts, B. 2002, ApJ, 577, 475 [Google Scholar]
  31. Ruderman, M. S., & Terradas, J. 2015, A&A, 580, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Russell, A. J. B., Simões, P. J. A., & Fletcher, L. 2015, A&A, 581, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Soler, R., & Terradas, J. 2015, ApJ, 803, 43 [Google Scholar]
  34. Stangalini, M., Giannattasio, F., Erdélyi, R., et al. 2017, ApJ, 840, 19 [NASA ADS] [CrossRef] [Google Scholar]
  35. Terradas, J., Andries, J., Goossens, M., et al. 2008, ApJ, 687, L115 [Google Scholar]
  36. Terradas, J., Magyar, N., & Van Doorsselaere, T. 2018, ApJ, 853, 35 [Google Scholar]
  37. Teunissen, J., & Keppens, R. 2019, Comput. Phys. Commun., 245, 106866 [Google Scholar]
  38. Van Doorsselaere, T., Debosscher, A., Andries, J., & Poedts, S. 2004, A&A, 424, 1065 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Van Doorsselaere, T., Srivastava, A. K., Antolin, P., et al. 2020a, Space Sci. Rev., 216, 140 [Google Scholar]
  40. Van Doorsselaere, T., Li, B., Goossens, M., Hnat, B., & Magyar, N. 2020b, ApJ, 899, 100 [NASA ADS] [CrossRef] [Google Scholar]
  41. Van Doorsselaere, T., Goossens, M., Magyar, N., Ruderman, M. S., & Ismayilli, R. 2021, ApJ, 910, 58 [NASA ADS] [CrossRef] [Google Scholar]
  42. Vaughan, S. 2005, A&A, 431, 391 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Wang, T. J., Solanki, S. K., & Selwa, M. 2008, A&A, 489, 1307 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Williams, T., Walsh, R. W., Winebarger, A. R., et al. 2020, ApJ, 892, 134 [Google Scholar]
  45. Xia, C., Teunissen, J., El Mellah, I., Chané, E., & Keppens, R. 2018, ApJS, 234, 30 [Google Scholar]
  46. Yuan, D., & Van Doorsselaere, T. 2016, ApJS, 223, 23 [Google Scholar]
  47. Zaitsev, V. V., & Stepanov, A. V. 1975, Issled. Geomagn. Aeron. Fiz. Solntsa, 37, 3 [NASA ADS] [Google Scholar]

All Figures

thumbnail Fig. 1.

Kinetic energy and velocity vectors in the cross-sectional slice at the midpoint of the flux tube for the circularly polarized standing kink mode after approximately two fundamental kink mode periods of evolution. Axis units are 10 Mm. Kinetic energy is in code units. The size of the velocity vectors is proportional to their magnitude. An animation of this figure is available online.

In the text
thumbnail Fig. 2.

Same as in Fig. 1, but for a linearly polarized standing kink mode. An animation of this figure is available online.

In the text
thumbnail Fig. 3.

Time-distance plots of Doppler signal integrated through a slit along the y-axis at the midpoint of the flux tube, with the LOS the positive x-axis. Left: linearly polarized kink. Right: circularly polarized kink. Doppler shift is in code units. Length and time units are ≈30 km and 42 s, respectively.

In the text
thumbnail Fig. 4.

Kinetic energy and velocity vectors in the cross-sectional slice at z = 50 Mm for the propagating circularly polarized kink wave. Axis units are 10 Mm. Kinetic energy is in code units. The size of the velocity vectors is proportional to their magnitude.

In the text
thumbnail Fig. 5.

Density and velocity vectors in the cross-sectional slice at the midpoint of the flux tube for the nonlinear setup after approximately two fundamental kink mode periods of evolution. Axis units are 10 Mm. Density is in code units. The size of the velocity vectors is proportional to their magnitude. An animation of this figure is available online.

In the text
thumbnail Fig. 6.

Path line plot obtained by integrating the velocity map over the duration of the simulation in the tube cross section, showing the displacement of the flux tube midpoint over time. The color corresponds to the snapshot number, thus showing an advance in time, with “0” being the initial condition and “50” denoting the final snapshot. Axis units are 10 Mm.

In the text
thumbnail Fig. 7.

Summed contribution of higher angular harmonics, summing over the Fourier components in a lower range, i.e., m = 3 5 | G ( m ) | 2 $ \sum\nolimits_{m=3}^{5} |G(m)|^2 $ (left) and higher range m = 6 15 | G ( m ) | 2 $ \sum\nolimits_{m=6}^{15} |G(m)|^2 $ (right), as a function of time, in units of ≈85.9 s (duration between snapshots).

In the text
thumbnail Fig. 8.

Steps for obtaining the boundary curve. Left: Polar transform of the density plot shown in Fig. 5, smoothed using the routine MedianFilter in Mathematica. Center: flux tube boundary identified from the polar transformed plot, using the routine MorphologicalPerimeter in Mathematica. Right: data points extracted from the flux tube boundary plot. Radius units are 10 Mm.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.