Free Access
Issue
A&A
Volume 637, May 2020
Article Number A75
Number of page(s) 13
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/201937083
Published online 19 May 2020

© ESO 2020

1. Introduction

Solar prominences (filaments when observed on the disk) are giant clouds of cold plasma suspended in the solar corona. The prominences are very dynamical structures subject to complex motions, including different types of oscillations. These oscillations have been classified into small- and large-amplitude oscillations (LAOs). The properties of LAOs differ greatly from small-amplitude oscillations. Velocities of LAOs are larger than 10 km s−1 and are often triggered by external perturbations such as Moreton and EIT waves (Eto et al. 2002; Gilbert et al. 2008; Liu et al. 2013) or extreme ultraviolet (EUV) waves (Liu et al. 2012; Shen et al. 2014a; Zhang & Ji 2018). In the case of LAOs, a large portion of the filament oscillates, reflecting global characteristics of the plasma and magnetic field structure. In contrast, small-amplitude oscillations are probably not excited externally, involving only a small part of the prominence (see review by Arregui et al. 2018). Large-amplitude oscillations can be classified according to the polarization of the oscillatory motions. Longitudinal LAOs are motions parallel to the magnetic field, whereas transverse LAOs are motions perpendicular to the magnetic field. In turn, transverse oscillations can be vertically or horizontally polarized. Detailed information on the prominence oscillations, their classification, and properties can be found in reviews by Tripathi et al. (2009) and, more recently, by Arregui et al. (2018).

Since the first observations of large-amplitude longitudinal oscillations (LALOs) by Jing et al. (2003, 2006), the number of the reported events has increased (e.g., Zhang et al. 2012; Li & Zhang 2012; Luna et al. 2014). These LALOs are characterized by a period of about one hour; the direction of these motions coincides with the direction of the magnetic field, forming an angle of approximately 25° with respect to the filament spine (Luna et al. 2014, 2018). This value of angle is in agreement with direct magnetic field measurements (e.g., Leroy et al. 1983, 1984; Casini et al. 2003). Several candidates have been proposed as a restoring force for LALOs. Luna & Karpen (2012) proposed the so-called pendulum model in which the main restoring force is gravity projected along the magnetic field, and the period depends only on the radius of curvature of the dipped magnetic field lines where the cool prominence plasma resides. Terradas et al. (2013), based on 2D magnetohydrodynamics (MHD) simulations, showed that the gas pressure gradient could make a significant contribution to the restoring force under some circumstances. Luna et al. (2012b) showed that the restoring force is a combination of the gravity projected along the field lines and the gas pressure gradient, but under realistic prominence conditions gravity is the dominant force. Numerical 2D and 3D MHD simulations (Luna et al. 2016; Zhou et al. 2018) have confirmed that gravity is the main restoring force in agreement with the pendulum model. However, the above numerical experiments are limited to a just single prominence model. It is desirable to investigate how the oscillatory parameters depend on the prominence properties by considering several prominence models.

The transverse oscillations of both vertical (Ramsey & Smith 1965, 1966; Eto et al. 2002; Okamoto et al. 2004; Gilbert et al. 2008; Zhang & Ji 2018) and horizontal polarizations (Isobe & Tripathi 2006; Hershaw et al. 2011) have also been observed. These oscillations usually have a short period between 10 min and 20 min. From a theoretical point of view, transverse oscillations were first studied in the works of Hyder (1966) and Kleczek & Kuperus (1969). The prominence was considered as a harmonic oscillator where the magnetic tension was the restoring force. Later, the properties of the transverse short-period oscillations were studied using 2D and 3D numerical simulations (Terradas et al. 2013; Luna et al. 2016; Zhou et al. 2018). Zhou et al. (2018) concluded that the magnetic tension is a dominant restoring force for both vertical and horizontal short-period oscillations. Terradas et al. (2013) and Luna et al. (2016) found that the vertical period remains constant with height, indicating that it is a global normal mode. Terradas et al. (2013) suggested that the vertical period has a simple dependence on the properties of the plasma and magnetic field. The period increases with the prominence density and decreases with the length of the field line. For prominence seismology, the relation between the vertical period and properties of the global magnetic structure can provide valuable information concerning the prominence structure.

As for the excitation mechanisms of LALOs, many observations reported LALO triggering by the disturbance close to the filament. In contrast, transverse LAOs are associated with distant flares and their corresponding EIT/Moreton waves (see review by Arregui et al. 2018). However, observations have also shown that different polarizations can be excited simultaneously by external disturbances. Gilbert et al. (2008) reported the event of the simultaneous longitudinal and transverse LAOs triggered by a Moreton wave. Using the Atmospheric Imaging Assembly (AIA) data from the Solar Dynamics Observatory (Lemen et al. 2012) and SMART Hα data (UeNo et al. 2004), Asai et al. (2012) found that the prominence oscillates in different directions simultaneously owing to external perturbation. Shen et al. (2014b) reported the excitation of the transverse oscillations of one filament and the longitudinal oscillations of another filament triggered by the same shock wave. The authors suggested that the motion type strongly depends on the propagation direction of the shock wave with respect to the prominence axis. Despite all these efforts, it is still unclear how the energetic EIT/Moreton disturbances propagate from a flare and produce the oscillations of different polarizations in a distant filament.

In this work, we aim to understand LAOs in a realistic prominence structure and how the properties of these LAOs depend on different parameters of the structure (i.e., density contrast and shear angle). Additionally, we study how an external disturbance excites different modes of oscillation of the structure.

This paper is organized as follows: in Sect. 2, an initial numerical configuration is described. In Sect. 3, we explain how the prominence was perturbed and how the plasma and magnetic field respond to the internal (horizontal and vertical) and external perturbations. In Sect. 5, the main results are summarized.

2. Numerical model

Our objective is to study the oscillatory properties of the prominence models triggered by internal and external perturbations. We numerically solve the equations of ideal MHD using the Mancha 3.0 code (Khomenko et al. 2008, 2014; Felipe et al. 2010; Khomenko & Collados 2012). We used the 2.5D approximation in which all the vectors have three spatial components, but the perturbation is only allowed to propagate in two dimensions. The governing equations and the corresponding source terms are described by Felipe et al. (2010). In this section, the numerical setup is outlined. We form the prominence structure using the cancelation process from a force-free arcade. Once the flux rope is formed, we load the cool prominence mass. The initial configuration and the processes of the flux rope formation, mass loading, and relaxation are explained in Sects. 2.12.3, respectively. The numerical domain consists of 800 × 600 grid points, which corresponds to the physical size of 192 Mm × 144 Mm (a spatial resolution of 0.24 Mm), except for one case in which we use a resolution that is four times higher. The Cartesian coordinate system is used with z- and x-axes defining the vertical and horizontal directions, respectively.

2.1. Initial equilibrium

The initial configuration consists of an isothermal atmosphere with a temperature of T = 2 MK gravitationally stratified along the vertical direction with a density at the bottom of the numerical domain of ρ = 1.31 × 10−12 kg m−3. We used an initial temperature larger than the typical value of the low corona of 1 MK. The reason for this selection is to have a large pressure scale height. This allowed us to have smoother stratification, avoiding huge Alfvén speeds at the top of the domain. The thermal conduction is not taken into account, but we included the radiative losses during the stage of the flux rope formation later on. We considered a monoatomic gas with an adiabatic constant, γ = 5/3. Using these values of the thermodynamical parameters, we obtained the maximum sound speed, c s = γ p 0 / ρ 0 = 166 $ c_{\mathrm{s}}=\sqrt{\gamma p_0/\rho_0}=166 $ km s−1, and density scale height, H = 47 Mm, in the initial atmosphere.

The initial force-free magnetic field is a periodic sheared arcade defined as

B x = 2 LB 0 π H b cos ( k ( x x 0 ) ) exp ( z z 0 H b ) , $$ \begin{aligned}&B_{x}=\frac{2{LB}_{0}}{\pi H_{\rm b}}\cos (k(x-x_{0}))\exp \left({-\frac{z-z_{0}}{H_{\rm b}}}\right), \end{aligned} $$(1)

B y = 1 ( 2 L π H b ) 2 B 0 cos ( k ( x x 0 ) ) exp ( z z 0 H b ) , $$ \begin{aligned}&B_{y}=\sqrt{1-\left(\frac{2L}{\pi H_{\rm b}}\right)^2}B_{0}\cos (k(x-x_{0}))\exp \left({-\frac{z-z_{0}}{H_{\rm b}}}\right), \end{aligned} $$(2)

B z = B 0 sin ( k ( x x 0 ) ) exp ( z z 0 H b ) , $$ \begin{aligned}&B_{z}=B_{0}\sin (k(x-x_{0}))\exp \left({-\frac{z-z_{0}}{H_{\rm b}}}\right), \end{aligned} $$(3)

where (x0, z0) = (0, 1) Mm are the coordinates of the center of the magnetic structure. The polarity inversion line (PIL), in this case, coincides with the line x = 0. The parameter L is related to the lateral extension k = 2π/L, and the constant A determines the magnetic scale height. In our experiments, we considered L = 48 Mm and B0 = 20 G.

It is useful to define an initial shear angle, θ, as

tan θ = B y B x = ( π H b 2 L ) 2 1 , $$ \begin{aligned} \tan \theta =\frac{B_{y}}{B_{x}}=\sqrt{\left(\frac{\pi H_{\rm b}}{2L}\right)^2-1}, \end{aligned} $$(4)

that is, the angle between the arcade field and the line perpendicular to the PIL. A large angle θ involves large magnetic shear. From Eq. (4), we see that there is a relation between the shear, θ, and the magnetic stratification, Hb. We considered the models with two initial shear angles, θ = 60°, 45°. The corresponding parameters of the structure have values of Hb = 61, 43 Mm, respectively. The initial magnetic field in the 2D projection is shown in Fig. 1a. In all the initial configurations considered the maximum Alfvén speed is v A = B / μ 0 ρ 0 = 1594 $ v_{\mathrm{A}}=B/\sqrt{\mu_0\rho_0}=1594 $ km s−1 and the plasma-β less than one everywhere in the domain.

thumbnail Fig. 1.

Magnetic field lines and density distribution in the central part of the computational domain during the prominence formation. Panel a: initial atmosphere and magnetic field. Panels b–d: processes of the convergence and the reconnection. Panels e, f: initial and final stages of the mass loading. Panel f: the dashed line indicates the magnetic field lines just before the mass loading process. The color bar scale indicates the logarithm of the density contrast.

2.2. Flux rope formation

Observations show that cool and dense prominence plasma is usually located in the dipped part of the curved magnetic field (see, e.g., reviews by Mackay et al. 2010; Vial & Engvold 2015) The dense prominence plasma is supported against gravity by the Lorentz magnetic force (Kippenhahn & Schlüter 1957; Kuperus & Raadu 1974). Several models have been proposed for the magnetic structure of prominences, but the flux rope structure is the most plausible. In this work, we consider the formation of the flux rope from the initial sheared arcade using the van Ballegooijen & Martens (1989) mechanism. This approach is similar to Kaneko & Yokoyama (2015). We only applied the converging motions at the bottom boundary in the direction toward the PIL. The horizontal velocity imposed at the base is written as

V x ( x , t ) = V 0 ( t ) sin [ 2 π ( x x 0 ) L ] exp [ ( x x 0 ) 2 2 σ 2 ] , $$ \begin{aligned} V_{x}(x, t)= -V_{0}(t)\sin {\left[\frac{2\pi (x-x_{0})}{L}\right]} \exp {\left[-\frac{(x-x_{0})^2}{2\sigma ^2}\right]}, \end{aligned} $$(5)

where parameter L = 189.6 Mm is the double width of the arcade, x0 = 0 Mm is the position of PIL, and σ = 13.6 Mm is related to the half-size of the converging region. We additionally set the Vy and Vz velocity components to zero at the bottom of the domain. We activated and disactivated the cancelation process smoothly with a function V0(t) defined as

V 0 ( t ) = V 00 { 0.5 [ erf ( t 2 λ t 1 λ ) + 1 ] 0.5 [ erf ( t 2 λ t 2 λ ) + 1 ] } , $$ \begin{aligned} V_0(t) =&V_{00} \left\{ 0.5\left[\mathrm{erf}\left(\frac{t-2\lambda -t_1}{\lambda }\right)+1\right] \right. \nonumber \\& \left.-0.5\left[\mathrm{erf}\left(\frac{t-2\lambda -t_2}{\lambda }\right)+1\right] \right\} , \end{aligned} $$(6)

where V00 = 3Vmax. The parameter Vmax = 38 km s−1 is the maximum converging velocity. Since we are interested in the final result of the flux rope formation, and not in the process, we used this high-velocity value to achieve the final stage quickly. The parameters t1 = 100 s and t2 = 3100 s determine the initiation and termination of the formation process. The parameter λ = 150 s determines the duration of the transition, and consequently, it is responsible for the smoothness of the cancelation process. The cancelation process produces a heating of the plasma of the flux rope. Consequently, the gas pressure inside of the rope and plasma-β also increase to values higher than desirable in a prominence (see, e.g., Hillier et al. 2012). This effect was already described by Kaneko & Yokoyama (2015). Following these authors, we include a radiative loss term in our energy equation (Eq. (3) from Felipe et al. 2010). We use Newton’s radiation law given by

Q rad = c v T 1 τ R , $$ \begin{aligned} Q_{\rm rad}=-c_{\rm v} \, \frac{T_1}{\tau _{\rm R}}, \end{aligned} $$(7)

where T1 is the perturbation of the temperature with respect to the initial 2 MK, τR is the radiative relaxation time, and cv is the specific heat at constant volume. This Qrad is only used during the flux rope formation process. We used a spatially constant and very small τR = 3 s to keep the temperature near 2 MK. At the end of the cancelation, the temperature is 2 MK, and then we deactivated the Qrad term. In Figs. 1a–d different stages of the formation are shown. Figure 1a shows the initial field of the magnetic arcade, described by Eqs. (1)–(3). Figure 1b shows how the footpoints of the field lines move toward the X = 0 axis following Eq. (5). When the field lines approach each other, they reconnect near the origin of the coordinates, forming a twisted magnetic structure (Fig. 1c). The flux rope slowly evolves through a series of stable equilibria and rises to the higher heights. The continuation of the converging process at the bottom could eventually lead the system to lose equilibrium and erupt. However, in the present work, we do not study the stability and eruption of the prominences. This problem can be addressed in future work. To avoid the eruption and to have a stable prominence, we stopped the convergence at t = 55 min. The resulting flux rope (see Fig. 1d) has its center at the height of 18.24 Mm. The magnetic field strength in the dipped part is about of 14 G. The characteristics of the flux rope are comparable to those observed in the solar prominences.

We used symmetric (for Vx, Vy, Bz) or antisymmetric (for Vz, Bx, By) boundary conditions at the top boundary. The gas pressure and density are fixed by assuming hydrostatic equilibrium at a constant temperature. The periodic boundary conditions are assumed at the left and right boundaries. At the bottom, we let the magnetic field evolve according to the velocity field. During the flux rope formation process, the magnetic field is computed from the induction equation using the converging velocity. Later on, the line-tied conditions are applied using zero-velocities. For the rest of the variables, zero-derivative conditions are assumed.

2.3. Mass loading and relaxation

In this work we do not consider the mechanisms of the prominence mass formation (see, e.g., Xia et al. 2011, 2012; Luna et al. 2012b; Xia & Keppens 2016). Instead of this, we placed the prominence mass artificially in the dip of the flux rope, adding a source term in the continuity equation, as in Terradas et al. (2013). Thereby, in a defined region of the domain, the density increases in comparison to the surrounding corona while the temperature decreases, keeping the gas pressure constant. The corresponding source for the continuity equation is determined as follows:

S ρ = χ ρ 0 t load exp ( ( x x p ) 4 σ x 4 ( z z p ) 4 σ z 4 ) , $$ \begin{aligned} S_{\rho }=\frac{\chi \rho _{0}}{t_{\rm load}}\exp \left({-\frac{(x-x_{\rm p})^4}{\sigma _x^4}-\frac{(z-z_{\rm p})^4}{\sigma _z^4}}\right), \end{aligned} $$(8)

where ρ0 and χ are the background density and density contrast, respectively. The time parameter tload = t4 − t3 = 500 s is associated with the loading rate, where t3 = 5000 s and t4 = 5500 s determine the initial and final stages of the mass loading process. The spatial parameters xp = 0 Mm, zp = 8.4 Mm are coordinates of the prominence center and σx = 3.6 Mm, σz = 2.4 Mm are related to the half-size of the mass distribution region. In this study we considered light and heavy prominences with χ = 100 and χ = 200.

The cancelation and the mass loading produce non-negligible vertical oscillations of the prominence structure. These oscillations are not desirable for our study. We defined a relaxation process that removes these oscillations. During the relaxation process, the velocities in the domain are decreased owing to artificial damping. At each iteration, the equations are updated with the new values of the velocities and total energy. It is necessary to let this system evolve and relax for a long time (4000 s in our case).

Figure 1e shows the density distribution at the beginning of the mass loading process. The flux rope center drops down to find a new equilibrium position. The dense mass modifies the magnetic field lines making them more elongated downward (see Fig. 1e). The magnetic forces increase as a result of lines stretching at the top and compression at the bottom. This way, magnetic forces compensate for the gravity increased due to the dense prominence mass. In the final stage of the relaxation, the prominence is close to equilibrium with velocities less than 2 km s−1. In the investigation of the prominence equilibrium in the 2.5D models, Hillier & van Ballegooijen (2013) showed that in the case of the absence of X-point, in the situation similar to ours, the prominence mass is supported by a combination of the magnetic tension force, magnetic pressure gradient, and gas pressure gradient.

In the final stage of the relaxation process, the prominence is located at a height of 7.4 Mm. The mass loading process produces a reconfiguration of the magnetic field and the plasma that produces a slight increase of the plasma-β in the dipped part of the flux rope. However, this increase is small, and the low-β situation is maintained.

3. Oscillations from internal perturbations

After flux rope formation, mass loading, and relaxation process, we obtained a prominence close to a static equilibrium at t = 150.0 min, as we described in Sect. 2. We first excited the system with internal perturbations to investigate the possible normal modes of different polarizations. This approach has already been used in several works (see, e.g., Terradas et al. 2013; Luna et al. 2016; Zhou et al. 2018) but for different configurations. We considered four prominence configurations corresponding to a combination of the two shear angles, θ = 45° ,60°, and the two maximum density contrasts, χ = 100, 200. Each system is perturbed in both vertical and horizontal directions. The perturbations can be produced by introducing a pulse in some of the quantities. For example, Terradas et al. (2013), Luna et al. (2016), and Zhou et al. (2018) imposed a perturbation exclusively in the velocity field and the rest of the magnitudes remain unchanged. These kinds of perturbations are unphysical because the rest of the magnitudes should also be modified to fulfill the governing equations. However, the system rapidly accommodates these perturbations in a few time steps. In some situations, if the initial perturbation is large, the numerical simulation crashes. We adopted an alternative approach to perturb the system. We incorporated a source term for each of the MHD equations in the Mancha 3.0 code (Felipe et al. 2010), similar to what we did for the mass loading (Sect. 2.3). The source term in the momentum equation represents an external force. This way, we produced self-consistent solutions. All the variables adjust to the source term following the MHD equations. The system is perturbed in either vertical or horizontal directions by the source terms in the z- or x-projection of the equation of motion, correspondingly. The source has the following form:

S m = ρ v pert t pert exp ( ( x x pert ) 4 σ x 4 ( z z pert ) 4 σ z 4 ) , $$ \begin{aligned} S_{\rm m}=\frac{\rho v_{\rm pert}}{t_{\rm pert}}\exp \left({-\frac{(x-x_{\rm pert})^4}{\sigma _x^4}-\frac{(z-z_{\rm pert})^4}{\sigma _z^4}}\right), \end{aligned} $$(9)

where tpert = 10 s is the duration of the initial perturbation, σx = σz = 12 Mm, xpert = 0 Mm and zpert = 7.4 Mm set the half-size and the x- and z-coordinates of the center of the perturbation, respectively. The maximum velocity of the perturbation, vpert = 20 km s−1, is in the range of the observed amplitudes.

3.1. Horizontal perturbation

Figure 2 shows the time sequence of the density and magnetic field evolution owing to the horizontal perturbation for the case (θ, χ) = (45° ,200). The three remaining cases are not shown because the temporal evolution is qualitatively similar. The solid lines represent the magnetic field, while the dashed lines indicate the magnetic field lines just before the perturbation at t = 150.0 min (at the end of the relaxation). The figure shows that the field lines change following the motion of the mass. However, the time-dependent changes in the shape of the magnetic field seem small. At the beginning, the prominence moves to the right and reaches its maximum displacement 3 min after the perturbation (Fig. 2a). In this initial stage, the prominence moves mainly as a whole (Figs. 2b and c). However, the bottom part is slightly delayed from the rest of the prominence. Later, the oscillations become less synchronized, the phase difference increases, resulting in a zigzag shape of the prominence (Figs. 2d and e). In the recent work by Luna et al. (2016), the zigzag appearance of LALOs is related to the differences in the period between the neighboring field lines. These phase shifts can make an impression of the vertical wave propagation. Such apparent wave motion has been called super-slow mode (Kaneko et al. 2015; Raes et al. 2017; Zapiór et al. 2019).

thumbnail Fig. 2.

Evolution of the density and magnetic field after the horizontal perturbation. The figure shows the dipped part of the flux rope. Panel a demonstrates the motion of the prominence just after the perturbation; panels b, c: lagging of the bottom part from the rest of the prominence body, and panels d, e demonstrate the zigzag shape of the motions. Panel f: state of the prominence in the final stage of simulation. The dashed lines indicate the magnetic field lines before the perturbation. The yellow lines represent the density isocontours.

From Fig. 2f, we see that the longitudinal oscillations are damped. This damping could be associated with the numerical dissipation. To test the convergence of our simulations and verify to which extent the damping is caused by the numerical reasons, we performed the simulations with a higher resolution (60 km). Figure 3 shows the velocity of the center of mass of the prominence projected along the magnetic field for the simulations with 240 km and 60 km resolution. This test shows that the damping in higher resolution simulation is weaker, indicating that the damping has a numerical contribution. It would be necessary to perform a dedicated study to understand this damping further by considering consecutively higher spatial resolutions. This can be a subject of future work. For both resolutions considered in this work, the periodicity is similar. We then conclude that the system has converged. An additional effect of the numerical diffusivity is a slight reduction of the prominence density with time. The high-resolution simulation shows a smaller reduction of the density. However, this density reduction does not seem to affect the prominence dynamics. Taking these results into account, for our purposes it is enough to consider a spatial resolution of 240 km.

thumbnail Fig. 3.

Longitudinal velocity variations at the position of the center of mass of the prominence for the simulations with different resolutions. The red and blue lines indicate the resolution of 240 km and 60 km, respectively.

As we see in Fig. 2, different parts of the structure oscillate differently, with different periods and phases. In addition, these motions are essentially nonlinear, and the plasma is a subject to important advection. The frozen-in condition applies in our fluid because we consider it to be perfectly conducting. Thus, it is necessary to advect every field line considered to catch the longitudinal and transverse motions of the plasma. We are interested in the cold plasma of the prominence, where the xz-projections of the field lines are more or less ellipses (Fig. 1f). We realized that the top part of these ellipses is approximately unperturbed. Thus, any field line starting at the top of the flux rope follows the plasma motion, simplifying the study of the oscillations. Following Luna et al. (2016), we calculated the longitudinal and transverse velocities at the positions of the center of mass of each field line. We selected 15 field lines of the flux rope with dips located at heights between 5 Mm and 10 Mm.

The temporal evolution of the longitudinal velocity for each of the selected field lines is shown in the left panel of Fig. 4. The velocity amplitude of these oscillations depends on the field line considered because of the shape of the triggering (Eq. (9)). The initial velocity is around 20 km s−1 for the field lines centered at the cold plasma around z = 7.5 Mm and is lower for the rest. For the field lines with z >  6 Mm, the oscillations show damping that depends on the height. However, for the field lines with z <  6 Mm, the amplitude is constant, or it slightly amplifies at the half time shown in the figure. From the edge of the flux rope to its center, the magnetic field lines become more twisted. Thus, we conclude that oscillations are strongly damped at places in which the lines have large curvature. In contrast to Luna et al. (2016), the phase differences between the lines are small. However, the combination of these small phase shifts and the different damping of the different layers produces the zigzag appearance.

thumbnail Fig. 4.

Temporal evolution of the longitudinal velocity (left panel) and vertical component of the transverse velocity (right panel) after the horizontal perturbation. The prominence model used has a shear angle of θ = 45° and a density contrast of χ = 200. Colored lines indicate velocity variations at the center of mass of different field lines. The colors denote the maximum initial density at the field lines according to the color bar scale. The right vertical axis indicates the height of the magnetic dips; the left vertical axis represents the scale of the maximum initial amplitude.

The horizontal perturbation excites not only the longitudinal oscillations but also the transverse oscillations, as we see in the right panel of Fig. 4. This figure shows that the initial vertical velocity is zero at all the height of the structure. However, immediately, the system starts to oscillate with the amplitudes of approximately 5 km s−1. Similarly to the longitudinal oscillations, for heights z >  6 Mm, the motions are clearly damped, whereas for z <  6 Mm, the plasma keeps oscillating.

We repeated the previous analysis for all four cases obtaining the temporal evolution at the different layers of the different prominence configurations. We analyzed the periods of oscillations at every field line considered as well as in the two polarization directions. The results for all four (θ, χ) configurations are plotted in Fig. 5. The red diamond symbols show the period of the longitudinal oscillations as a function of the height of the dips. In all four cases, we can see that the period tends to decrease with height. The difference between the minimum and the maximum value is 2 − 3 min (see Table 1). The longitudinal period slightly changes for different field lines. As a result, the prominence is shaped in the form of a zigzag after several cycles of oscillations (Fig. 2e). The vertical period due to the horizontal perturbation is plotted with the blue circle symbols in Fig. 5. The vertical oscillations have periods that are half of the longitudinal periods. The behavior with height is identical to that of the longitudinal oscillation. This result indicates that the vertical motion is the back-reaction to the longitudinal oscillations, as shown in Fig. 4. For example, when the prominence mass passes through the dip, it pushes the structure downward. This happens twice for each period explaining why the vertical oscillations have half of the period of the longitudinal oscillation.

thumbnail Fig. 5.

Oscillation periods as a function of height of the magnetic dips. Top panels: models with θ = 45° and bottom panels: models with θ = 60°. Left and right panels: periods of the heavy and light prominence models, respectively. Red diamonds: longitudinal period; blue circles: vertical period due to the horizontal perturbation; green triangles: vertical period due to the vertical perturbation; and solid line: period from the theoretical pendulum model.

Table 1.

Radii of curvature, the longitudinal and transverse vertical periods of the configurations with the different shear angle, θ, and the density contrast, χ.

We computed the radii of curvature averaged around the bottom of the dip to compare our results with the pendulum model. The magnetic field evolves as a consequence of the motion of the prominence plasma. Because of that, we used the time-averaged values of the radii of curvature, Rc. This allowed us to calculate the theoretical period of the longitudinal oscillations. We obtained the pendulum period shown as a solid line in Fig. 5 by using equation, P = 2 π R c / g $ P=2\pi \sqrt{R_{\mathrm{c}}/g} $, from Luna & Karpen (2012), where g = 273.94 m s−2 is the solar gravitational acceleration. For θ = 45° the periods obtained in the numerical experiment show a good agreement with the theoretical period for z <  7.5 Mm and larger discrepancy for larger heights. In contrast, for θ = 60° the agreement is much better in all the spatial domain. We conclude that, for this type of oscillation, the restoring force is the projected solar gravity. The remaining difference between the numerical experiment and the theory of about 7%−15% can be related to the local deformation of the magnetic field lines due to the motion of the heavy mass. Zhang et al. (2019) and Zhou et al. (2018) show using 2D and 3D simulations, respectively, that the pendulum model might lead to a large error in determining the radius of curvature using longitudinal periods when the gravity to Lorentz force ratio is close to unity. The authors define a dimensionless parameter to quantify this ratio as δ = 2 g L p / V A 2 $ {\delta=2gL_{\mathrm{p}}/V_{\mathrm{A}}^{2}} $ being Lp ≅ 10 Mm a typical dimension of the prominence and the local Alfvén speed, VA ≅ 100 km s−1, in the case of χ = 200. In our situation, this parameter is close to δ = 0.6, indicating that the deformation of the field lines by the gravity could have some moderate effect on the longitudinal oscillations.

A comparison between the different simulations in Fig. 5 shows that a change of the shear angle or the density contrast does not lead to a large difference in the longitudinal period. Comparing the heavy and light prominence models in Fig. 5, we can see that the longitudinal period is only 1 − 2 min longer in the case of the lower density contrast, χ = 100. Since we concluded that the projected solar gravity is the main restoring force, the period should essentially depend only on the radius of curvature. This means that the radii of curvature are slightly larger, where prominence is less massive. As mentioned before, the loaded cold mass deforms the magnetic field lines decreasing the radius of curvature (Fig. 1f). It is expected that heavier mass affects the shape of the field lines more, as a consequence, decreasing the longitudinal period.

The top and bottom panels of Fig. 5 compare oscillations for the models with θ = 45° and θ = 60°. The longitudinal period does not show a large difference in these two cases. The difference of 3 − 4 min can be related to the projection effect. The more sheared structure has a smaller radius of curvature in the 2D projection (Table 1). As a result, the longitudinal period is also slightly shorter.

In order to estimate the damping time, the signal of the longitudinal velocity at the position of the center of mass of the prominence is fitted by the decayed sinusoidal function v = v0et/τDsin(2πt/P + ϕ). We computed the damping time for the case shown in Fig. 4. For the lines at the bottom of the prominence, the damping time τD = 100 min, while at the top, the damping time is only 60 min. Thus, the damping time shows a decrease with height, in agreement with the visual impression from the left panel of Fig. 4. We repeated the analysis for each of the models and obtained a similar behavior of damping with height.

3.2. Vertical perturbation

We perturb the prominence along the vertical direction by applying the source term (Eq. (9)) to the z-component of the equation of motion. The temporal evolution of the prominence mass just after the vertical perturbation is shown in Fig. 6. Figure 6a shows the vertical upward motion of the prominence right after the perturbation. We see how the field follows the plasma motion due to the frozen-in condition and the dips flatten. In Fig. 6b, the direction of the motion has changed. The motion is downward and the magnetic field is bent in comparison to the initial configuration. In addition to the vertical motion, the prominence is a subject to the longitudinal motions. These longitudinal motions are due to the compression and rarefaction of the plasma. In Fig. 6b, we can see clearly how the prominence expands symmetrically with respect to the x = 0 axis along the field lines. In Figs. 6d and f, prominence shows symmetric longitudinal motions. This motion is more noticeable as the difference between the behavior at the bottom and at the upper prominence layer. In Fig. 6f, it is seen that the bottom part is more expanded than the upper part.

thumbnail Fig. 6.

Evolution of the density and magnetic field after the vertical perturbation. Panels b–f demonstrate the longitudinal motions at the bottom of the prominence due to compression and rarefaction. The dashed lines denote the magnetic field lines before the perturbation. The yellow lines represent the density isocontours.

The behavior of the different components of the velocity at several prominence layers is shown in Fig. 7. The right panel of Fig. 7 shows the vertical oscillations at the center of mass of each field line, as in Sect. 3.1. However, owing to the symmetry of the vertical oscillations, the longitudinal velocity is zero at the center of mass of each line. For this reason, we measured the longitudinal velocity outside of the center close to the position of the maximum velocity. Figure 7 (left) shows the temporal evolution of the longitudinal velocity at the selected field lines. We can see that the initial amplitude is zero, but the longitudinal oscillations start rapidly. Figure 6b shows the expansion of the prominence body when it moves downward. The amplitudes increase, reaching the value of 20 km s−1 at some lines. A very interesting feature of this numerical experiment that a vertical perturbation also produces an important horizontal oscillation of the cool plasma. It is unclear if, in the first stages of the oscillation, there is a mode coupling between vertical and horizontal polarization modes or, in contrast, it is just the transition phase of the establishment of modes, where the oscillation is composed of vertical and horizontal compression and rarefaction movements. Generally, the energy transfer could be related to the fast-to-slow mode conversion in the region where plasma-β is close to unity (Schunker & Cally 2006; Chen et al. 2016) or to the resonant absorption in the Alfvén continuum at the prominence-corona transition region (PCTR; Goossens et al. 2014; Antolin et al. 2015). Since we have plasma-β lower than unity everywhere in the numerical domain, the mechanism of the fast-to-slow mode conversion is not taking place. The resonant absorption could still be responsible for the energy transfer and the damping of the global prominence modes. If this mechanism is in action, the global mode of the prominence is transferred to the local modes in the magnetic surfaces of the PCTR. We studied the distribution of the kinetic energy with time and we have not detected any transfer of energy to magnetic surfaces of the structure. Unlike our case, the energy of transverse oscillations was shown to be transferred to a thin region at the PCTR in 3D MHD simulations by Terradas et al. (2016). The nonexistence of this phenomenon in our case could be associated with limitations of the 2.5D modeling. The resonant absorption in 3D flux rope structures will be the subject of a future study.

thumbnail Fig. 7.

Temporal evolution of the longitudinal velocity (left panel) and vertical component of the transverse velocity (right panel) after the vertical perturbation. The prominence model used has a shear angle of θ = 45° and a density contrast of χ = 200. The colored lines indicate velocity variations at the center of mass (right) or the certain distance from the center of mass (left) at each field line. The colors indicate the maximum initial density at the field line according to the color bar scale. The right vertical axis indicates the height of the magnetic dips; the left vertical axis represents the scale of the maximum initial amplitude.

After two minutes, the longitudinal motions seem to be established. However, the oscillations are very irregular with a complex pattern. These longitudinal motions are associated with compression and rarefaction, as shown in Figs. 6a–f, rather than related to the projected solar gravity. The periods of these longitudinal modes (see Table 1) are shorter than those of the pendulum-like motions of Sect. 3.1. Also, these periods are not identical or the double of the vertical periods indicating that the longitudinal modes are not related to the vertical modes. Altogether, this may indicate that the vertical perturbation also excites overtones of the structure with antisymmetric velocity profiles, as studied in Joarder & Roberts (1992), Oliver et al. (1993), and Luna et al. (2012a). The temporal evolution of the vertical velocity is shown in the right panel of Fig. 7. The vertical oscillations in the different layers of the structure show a very coherent motion with oscillations almost in phase at all heights. The periods are very uniform with values shown in Table 1 and in Fig. 5. Also, in all the field lines, the attenuation is moderate and very uniform along the vertical direction. As we explained previously, there is no resonant absorption or fast-to-slow mode conversion, and these can be discarded as possible damping mechanisms. However, this attenuation can be explained by numerical diffusivity or by the wave leakage. Figure 8 shows a time-distance diagram of the vertical velocity along a vertical cut at the center of the structure. The red and blue colors correspond to the upward and downward motions, respectively. The prominence is located between 5 Mm and 10 Mm. As we can see from this figure, strong vertical motion appears after the perturbation at heights above the prominence and it extends even higher. The vertical bands have a slight inclination. This inclination reveals that there is a propagation of the velocity perturbations in the form of waves. The inclination is consistent with a velocity of 1000 km s−1. This means that some energy leaves the system during this time owing to the wave leakage. Taking the attenuation into account shown in the right panel of Fig. 7, we conclude that this wave leakage contributes to the damping of the vertical oscillations. We calculated the damping time, as described in Sect. 3.1. The damping time is about of 60 min in almost all the lines and is slightly shorter at the bottom (about of 50 min).

thumbnail Fig. 8.

Time-distance diagram of the vertical component of the transverse velocity along z-axis at x = 0, in the model with the shear angle of θ = 45° and the density contrast of χ = 200. Black lines denote a region where the prominence mass is located.

In several works (Hyder 1966; Kleczek & Kuperus 1969; Luna et al. 2016; Zhou et al. 2018), the authors have concluded that the restoring force for the transverse oscillations has a magnetic origin. We calculated the magnetic forces and their relative variations as done by Zhou et al. (2018). The variations of the magnetic pressure gradient exceed the variations of the magnetic tension force for all considered configurations. This indicates that in our model, the magnetic pressure gradient is the main restoring force for the vertical oscillations.

We obtained the period of the vertical oscillations as in Sect. 3.1 (Table 1). The vertical period is plotted with green triangles as a function of the height of the dips in Fig. 5. The oscillatory period remains constant at different layers of the prominence for all models. This indicates that the vertical motion is associated with global normal mode (Terradas et al. 2013; Luna et al. 2016). From Fig. 5, we see that the transverse periods are much shorter than longitudinal periods, similar to as revealed by observations (see review Tripathi et al. 2009, Table 1). Same as the longitudinal period, the vertical period does not change much for different models with an average difference of around 1 min (see Table 1). Part of the variation of the vertical period for different density contrast can be due to the inertia effect. For the more massive prominence, with the density contrast χ = 200, the vertical period is slightly longer (Figs. 5a and c). Comparing the top and bottom panels of Fig. 5, we can see that the vertical period seems to be almost insensitive to the shear angle variation.

4. Oscillations from external perturbation

Usually, LAOs are triggered by external perturbations such as Moreton or EIT waves. It is crucial to understand how external perturbation interacts with the flux rope structure and produces oscillations of the cold prominence. We modeled this scenario using the perturbation that arrives from outside. In order to produce strong perturbation comparable to the realistic energetic events, we used the source term of the total energy equation instead of perturbation of pressure or velocity. In this case, we add energy to the system, all variables are adjusted following the MHD equations, and we have a self-consistent perturbation (see the discussion on Sect. 3.1). We placed the perturbation at some distance from the prominence center as follows:

S e = α t pert exp ( ( x x pert ) 2 σ x 2 ( z z pert ) 2 σ z 2 ) , $$ \begin{aligned} S_{\rm e}=\frac{\alpha }{t_{\rm pert}}\exp \left({-\frac{(x-x_{\rm pert})^2}{\sigma _x^2}-\frac{(z-z_{\rm pert})^2}{\sigma _z^2}}\right), \end{aligned} $$(10)

where the parameter α = 2 J m−3. In solar flares temperature and density can be up to 107 − 108 K and 1016 m−3 that corresponds to the parameter α from 2 J m−3 to 20 J m−3. With this parameter, we produce an initial corresponding velocity around of 360 km s−1. This value is typical for EIT waves that are usually associated with flares or coronal mass ejections (Moses et al. 1997; Thompson et al. 1998). In turn, tpert = 50 s, xpert = −45.6 Mm, zpert = 12 Mm, σx = σz = 12 Mm are the same parameters as in Eq. (9).

Our aim is to study how the wavefront propagates and interacts with the flux rope. The source of the wave should be placed at a certain distance from the flux rope. We consider a periodic system, virtually composed of an infinite number of arcades separated by a region of the vertical field. In our domain, we not only have the main arcade but also the neighboring arcades. If the perturbation were located far (in the area of the neighboring arcade), the wavefront would propagate along the field lines to the boundary, and most energy would remain trapped in the neighboring arcade with only a small part of the energy reaching the prominence. To avoid this behavior, we choose the initial position of the perturbation between the arcades, where the magnetic field is predominantly vertical. The temporal evolution of the system is shown in Fig. 9 (left panels). The right-hand panels show the same, but with a more detailed view around the cold prominence. Figures 9a and b show the initial distribution of the gas pressure. We can see high pressure inside the flux rope in comparison with the environment. In this region, pressure increases during the flux rope formation process. In order to visualize the propagation of the perturbation, we compute the running differences by subtracting two consecutive snapshots separated by 30 s (Figs. 9c–h). As shown in Figs. 9c and d, just after the source term activation, the external perturbation propagates mainly vertically, producing significant changes in the magnetic field. The field lines are displaced to the right with respect to the equilibrium structure (dashed lines). The external gas pressure perturbation is unable to cross the flux rope and does not reach the internal region of the prominence. However, the compression of the lines creates a gas pressure excess inside the flux rope near the cold prominence plasma (Figs. 9e and f). From Figs. 9e and f, we see that pressure perturbation inside the flux rope propagates crossing the cold prominence, but the prominence is not deformed or displaced by this perturbation. Figure 9g shows that the external wavefront continues propagating along the arcade field lines producing significant perturbations of the magnetic field. In Fig. 9h, we see that the internal gas pressure pulse has crossed the prominence, slightly deforming it. After this transient, an oscillation is established. From this temporal evolution, we see that the magnetic field is responsible for the movement of the prominence. The external perturbation produces an Alfvénic disturbance that displaces the flux rope, and then the cold prominence mass moves. In this sense, the motion of the prominence is secondary because the perturbation firstly perturbs the flux rope.

thumbnail Fig. 9.

Evolution of the magnetic field and gas pressure during the wavefront propagation. Panels a, b: initial gas pressure distribution; panels c–h: running difference of the two consecutive snapshots. Left-hand panels: same as the right-hand panels but focused at the location around the prominence. The dashed lines indicate the magnetic field lines before the perturbation. The yellow lines represent the density isocontours.

Chen et al. (2016) studied the interaction between the fast-mode shock wave and the arcade magnetic field structure using 2D MHD simulations of the periodic arcade magnetic field. The system did not include the prominence mass. The authors concluded that if a wave vector is parallel to the magnetic field and the wave propagates in the region where plasma-β is close to unity, the fast-to-slow mode conversion happens (see Chen et al. 2016, Fig. 7). In the present case, the shock wave propagates in the low-β environment and we do not observe the fast-to-slow mode conversion.

As in the previous sections, we obtain the longitudinal and vertical velocities at the field lines center of mass between z = 5 Mm and 10 Mm. From both panels of Fig. 10, we can see a complex dynamical behavior. In the longitudinal velocity, there are no clear signatures of oscillations in the first ten minutes after the perturbation. In the bottom part of the structure, there are strong peaks associated with the external disturbance that propagates along the circular magnetic field lines of the flux rope. After some time, it reaches the bottom part of the flux rope, producing these spikes. The oscillations have an initial amplitude of about 5 km s−1. The period of the longitudinal oscillation slightly decreases with height that resembles the behavior of the period with height shown in the left panel of Fig. 4. The vertical oscillations are established after the first five minutes of the perturbation (see Fig. 10 right panel). The amplitude of the vertical oscillations is around 20 km s−1. We can see that the signal reaches the maximum amplitude after several cycles of oscillations. The motions are coherent in the different prominence layers. The damping of the signal is similar to the situation shown in the right panel of Fig. 7. Therefore, as in the case of the vertical perturbation, this attenuation can be associated with the wave leakage (see Sect. 3.2). The periods of oscillations shown in Table 1 are in agreement with internal modes. From this numerical experiment, we conclude that external perturbation can excite both longitudinal and vertical oscillations. Even if the vertical oscillations are excited with larger amplitudes, they are also damped faster than the longitudinal modes. In this situation, the oscillation is mainly vertical just after the triggering, but the final phase contains mainly longitudinal oscillations.

thumbnail Fig. 10.

Temporal evolution of the longitudinal velocity (left panel) and vertical component of the transverse velocity (right panel) after the wavefront propagation. The prominence model used has a shear angle of θ = 45° and a density contrast of χ = 200. The colored lines indicate velocity variations at the center of mass of different field lines. The colors indicate the maximum initial density at the field lines according to the color bar scale. The right vertical axis indicates the height of the magnetic dips; the left vertical axis represents the scale of the maximum initial amplitude.

5. Summary and conclusions

To study the physics of the LAOs, we performed 2.5D numerical simulations using the mechanism of van Ballegooijen & Martens (1989) for the flux rope formation. Once the flux rope was formed, the cold mass is artificially loaded in the dipped part of the magnetic structure. After the relaxation process, we applied the horizontal and vertical perturbations at the prominence position and studied internal modes of different polarizations. We also considered the case of the external LAOs triggering using the disturbance placed out of the flux rope.

The horizontal perturbation triggers both longitudinal and transverse vertical oscillations. The period of the longitudinal oscillations, obtained numerically, shows a good agreement with the pendulum period. This confirms that, under the prominence conditions, gravity is the dominant restoring force for the longitudinal oscillations. The damping of the longitudinal motions is mostly related to the numerical dissipation. The vertical oscillations are associated with the back-reaction of the magnetic field to the longitudinal motions.

The vertical perturbation drives vertical motions together with the longitudinal component, associated with compression and rarefaction. The period of the vertical oscillations remains constant with height. This indicates that the vertical disturbance triggers global normal mode of the structure. The oscillations show attenuation. We found indications that the wave leakage could be responsible for the damping of the vertical oscillations.

We considered oscillations of different polarizations using models with the two values of shear angle and density contrast. We compared the periods of the different prominence configurations. The periods of longitudinal and vertical oscillations show a weak dependence on these parameters of the prominence structure. In the case of the longitudinal period, the small difference in the period is related to the variations of the radius of curvature. In contrast, the vertical periods are slightly different for the heavy and light prominences due to the inertia effect.

We investigated the triggering mechanism by external disturbance, choosing the location of the source at a certain distance from the flux rope. The external perturbation reaches the flux rope displacing and deforming the magnetic field lines. The flux rope displacement creates perturbation inside the prominence and, as a consequence, the prominence is also displaced from the equilibrium position. The signal analysis related to the different field lines has shown that the plasma motions have a complex character and external disturbance drives both longitudinal and transverse vertical LAOs.

Using the time-dependent numerical simulations, we demonstrated that the period of the longitudinal oscillations is in a good agreement with the pendulum model, and the period of the vertical oscillations only slightly depends on the prominence density and the shear angle of the magnetic field. We also proposed a model of the external LAOs excitation. This simulation shows that the external disturbance excites the normal prominence modes of the different polarizations.

In the future, it will be necessary to investigate the oscillatory properties, using 3D prominence models. Firstly, in a 3D model, the flux rope can be anchored at the footpoints, which allows the magnetic tension to make a significant contribution to the restoring force. Therefore, including the perpendicular direction can significantly change the oscillatory properties of this model. Secondly, 3D simulations can allow us to study external triggering in more detail. It is possible to reproduce the scenario when a wave arrives from different directions with respect to the prominence spine. For future work, it is essential to consider more realistic external disturbances such as an energy release due to the magnetic reconnection or the energetic wave associated with the prominence eruption.

Acknowledgments

V. Liakh acknowledges the support of the Instituto de Astrofísica de Canarias via an Astrophysicist Resident fellowship. M.L. acknowledges the support by the Spanish Ministry of Economy and Competitiveness (MINECO) through Severo Ochoa Program MINECO SEV-2015-0548 and the support from “Ministerio de Economia, Industria y Competitividad” for the Ramón y Cajal fellowship RYC2018-026129-I. E. K. thanks the support by the European Research Council through the Consolidator Grant ERC-2017-CoG-771310-PI2FA and by the Spanish Ministry of Economy, Industry and Competitiveness through the grant PGC2018-095832-B-I00 is acknowledged. V. Liakh, M. Luna, and E. Khomenko thankfully acknowledge the technical expertise and assistance provided by the Spanish Supercomputing Network (Red Española de Supercomputacíon), as well as the computer resources used: the LaPalma Supercomputer, located at the Instituto de Astrofísica de Canarias. V. L. also acknowledges the support from ISSI www.issibern.ch to team 413 and team 431 on “Large-Amplitude Oscillations as a Probe of Quiescent and Erupting Solar Prominences” led by M. Luna.

References

  1. Antolin, P., Okamoto, T. J., De Pontieu, B., et al. 2015, ApJ, 809, 72 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arregui, I., Oliver, R., & Ballester, J. L. 2018, Liv. Rev. Sol. Phys., 15, 3 [Google Scholar]
  3. Asai, A., Ishii, T. T., Isobe, H., et al. 2012, ApJ, 745, L18 [NASA ADS] [CrossRef] [Google Scholar]
  4. Casini, R., López Ariste, A., Tomczyk, S., & Lites, B. W. 2003, ApJ, 598, L67 [NASA ADS] [CrossRef] [Google Scholar]
  5. Chen, P. F., Fang, C., Chandra, R., & Srivastava, A. K. 2016, Sol. Phys., 291, 3195 [NASA ADS] [CrossRef] [Google Scholar]
  6. Eto, S., Isobe, H., Narukage, N., et al. 2002, PASJ, 54, 481 [NASA ADS] [CrossRef] [Google Scholar]
  7. Felipe, T., Khomenko, E., & Collados, M. 2010, ApJ, 719, 357 [Google Scholar]
  8. Gilbert, H. R., Daou, A. G., Young, D., Tripathi, D., & Alexander, D. 2008, ApJ, 685, 629 [NASA ADS] [CrossRef] [Google Scholar]
  9. Goossens, M., Soler, R., Terradas, J., Van Doorsselaere, T., & Verth, G. 2014, ApJ, 788, 9 [NASA ADS] [CrossRef] [Google Scholar]
  10. Hershaw, J., Foullon, C., Nakariakov, V. M., & Verwichte, E. 2011, A&A, 531, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Hillier, A., & van Ballegooijen, A. 2013, ApJ, 766, 126 [NASA ADS] [CrossRef] [Google Scholar]
  12. Hillier, A., Hillier, R., & Tripathi, D. 2012, ApJ, 761, 106 [NASA ADS] [CrossRef] [Google Scholar]
  13. Hyder, C. L. 1966, Z. Astrophys., 63, 78 [NASA ADS] [Google Scholar]
  14. Isobe, H., & Tripathi, D. 2006, A&A, 449, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Jing, J., Lee, J., Spirock, T. J., et al. 2003, ApJ, 584, L103 [NASA ADS] [CrossRef] [Google Scholar]
  16. Jing, J., Lee, J., Spirock, T. J., & Wang, H. 2006, Sol. Phys., 236, 97 [NASA ADS] [CrossRef] [Google Scholar]
  17. Joarder, P. S., & Roberts, B. 1992, A&A, 261, 625 [NASA ADS] [Google Scholar]
  18. Kaneko, T., & Yokoyama, T. 2015, ApJ, 806, 115 [NASA ADS] [CrossRef] [Google Scholar]
  19. Kaneko, T., Goossens, M., Soler, R., et al. 2015, ApJ, 812, 121 [NASA ADS] [CrossRef] [Google Scholar]
  20. Khomenko, E., & Collados, M. 2012, ApJ, 747, 87 [NASA ADS] [CrossRef] [Google Scholar]
  21. Khomenko, E., Collados, M., & Felipe, T. 2008, Sol. Phys., 251, 589 [NASA ADS] [CrossRef] [Google Scholar]
  22. Khomenko, E., Díaz, A., de Vicente, A., Collados, M., & Luna, M. 2014, A&A, 565, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Kippenhahn, R., & Schlüter, A. 1957, Z. Astrophys., 43, 36 [Google Scholar]
  24. Kleczek, J., & Kuperus, M. 1969, Sol. Phys., 6, 72 [NASA ADS] [CrossRef] [Google Scholar]
  25. Kuperus, M., & Raadu, M. A. 1974, A&A, 31, 189 [NASA ADS] [Google Scholar]
  26. Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 [NASA ADS] [CrossRef] [Google Scholar]
  27. Leroy, J. L., Bommier, V., & Sahal-Brechot, S. 1983, Sol. Phys., 83, 135 [Google Scholar]
  28. Leroy, J. L., Bommier, V., & Sahal-Brechot, S. 1984, A&A, 131, 33 [Google Scholar]
  29. Li, T., & Zhang, J. 2012, ApJ, 760, L10 [NASA ADS] [CrossRef] [Google Scholar]
  30. Liu, W., Ofman, L., Nitta, N. V., et al. 2012, ApJ, 753, 52 [NASA ADS] [CrossRef] [Google Scholar]
  31. Liu, R., Liu, C., Xu, Y., et al. 2013, ApJ, 773, 166 [Google Scholar]
  32. Luna, M., & Karpen, J. 2012, ApJ, 750, L1 [NASA ADS] [CrossRef] [Google Scholar]
  33. Luna, M., Karpen, J. T., & DeVore, C. R. 2012a, ApJ, 746, 30 [NASA ADS] [CrossRef] [Google Scholar]
  34. Luna, M., Díaz, A. J., & Karpen, J. 2012b, ApJ, 757, 98 [NASA ADS] [CrossRef] [Google Scholar]
  35. Luna, M., Knizhnik, K., Muglach, K., et al. 2014, ApJ, 785, 79 [NASA ADS] [CrossRef] [Google Scholar]
  36. Luna, M., Terradas, J., Khomenko, E., Collados, M., & de Vicente, A. 2016, ApJ, 817, 157 [NASA ADS] [CrossRef] [Google Scholar]
  37. Luna, M., Karpen, J., Ballester, J. L., et al. 2018, ApJS, 236, 35 [NASA ADS] [CrossRef] [Google Scholar]
  38. Mackay, D. H., Karpen, J. T., Ballester, J. L., Schmieder, B., & Aulanier, G. 2010, Space Sci. Rev., 151, 333 [NASA ADS] [CrossRef] [Google Scholar]
  39. Moses, D., Clette, F., Delaboudinière, J.-P., et al. 1997, Sol. Phys., 175, 571 [NASA ADS] [CrossRef] [Google Scholar]
  40. Okamoto, T. J., Nakai, H., Keiyama, A., et al. 2004, ApJ, 608, 1124 [NASA ADS] [CrossRef] [Google Scholar]
  41. Oliver, R., Ballester, J. L., Hood, A. W., & Priest, E. R. 1993, ApJ, 409, 809 [NASA ADS] [CrossRef] [Google Scholar]
  42. Raes, J. O., Van Doorsselaere, T., Baes, M., & Wright, A. N. 2017, A&A, 602, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Ramsey, H., & Smith, S. F. 1965, AJ, 70, 688 [NASA ADS] [CrossRef] [Google Scholar]
  44. Ramsey, H. E., & Smith, S. F. 1966, AJ, 71, 197 [NASA ADS] [CrossRef] [Google Scholar]
  45. Schunker, H., & Cally, P. S. 2006, MNRAS, 372, 551 [NASA ADS] [CrossRef] [Google Scholar]
  46. Shen, Y., Ichimoto, K., Ishii, T. T., et al. 2014a, ApJ, 786, 151 [NASA ADS] [CrossRef] [Google Scholar]
  47. Shen, Y., Liu, Y. D., Chen, P. F., & Ichimoto, K. 2014b, ApJ, 795, 130 [NASA ADS] [CrossRef] [Google Scholar]
  48. Terradas, J., Soler, R., Díaz, A. J., Oliver, R., & Ballester, J. L. 2013, ApJ, 778, 49 [NASA ADS] [CrossRef] [Google Scholar]
  49. Terradas, J., Soler, R., Luna, M., et al. 2016, ApJ, 820, 125 [NASA ADS] [CrossRef] [Google Scholar]
  50. Thompson, B. J., Plunkett, S. P., Gurman, J. B., et al. 1998, Geophys. Res. Lett., 25, 2465 [Google Scholar]
  51. Tripathi, D., Isobe, H., & Jain, R. 2009, Space Sci. Rev., 149, 283 [NASA ADS] [CrossRef] [Google Scholar]
  52. UeNo, S., Nagata, S.-I., Kitai, R., Kurokawa, H., & Ichimoto, K. 2004, Ground-based Instrumentation for Astronomy, eds. A. F. M. Moorwood, & M. Iye, Proc. SPIE, 5492, 958 [NASA ADS] [CrossRef] [Google Scholar]
  53. van Ballegooijen, A. A., & Martens, P. C. H. 1989, ApJ, 343, 971 [NASA ADS] [CrossRef] [Google Scholar]
  54. Vial, J. C., & Engvold, O. 2015, in Solar Prominences (Cham: Springer), Astrophys. Space Sci. Lib., 415 [Google Scholar]
  55. Xia, C., & Keppens, R. 2016, ApJ, 825, L29 [NASA ADS] [CrossRef] [Google Scholar]
  56. Xia, C., Chen, P. F., Keppens, R., & van Marle, A. J. 2011, ApJ, 737, 27 [NASA ADS] [CrossRef] [Google Scholar]
  57. Xia, C., Chen, P. F., & Keppens, R. 2012, ApJ, 748, L26 [Google Scholar]
  58. Zapiór, M., Schmieder, B., Mein, P., et al. 2019, A&A, 623, A144 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Zhang, Q. M., & Ji, H. S. 2018, ApJ, 860, 113 [NASA ADS] [CrossRef] [Google Scholar]
  60. Zhang, Q. M., Chen, P. F., Xia, C., & Keppens, R. 2012, A&A, 542, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Zhang, L. Y., Fang, C., & Chen, P. F. 2019, ApJ, 884, 74 [NASA ADS] [CrossRef] [Google Scholar]
  62. Zhou, Y.-H., Xia, C., Keppens, R., Fang, C., & Chen, P. F. 2018, ApJ, 856, 179 [Google Scholar]

All Tables

Table 1.

Radii of curvature, the longitudinal and transverse vertical periods of the configurations with the different shear angle, θ, and the density contrast, χ.

All Figures

thumbnail Fig. 1.

Magnetic field lines and density distribution in the central part of the computational domain during the prominence formation. Panel a: initial atmosphere and magnetic field. Panels b–d: processes of the convergence and the reconnection. Panels e, f: initial and final stages of the mass loading. Panel f: the dashed line indicates the magnetic field lines just before the mass loading process. The color bar scale indicates the logarithm of the density contrast.

In the text
thumbnail Fig. 2.

Evolution of the density and magnetic field after the horizontal perturbation. The figure shows the dipped part of the flux rope. Panel a demonstrates the motion of the prominence just after the perturbation; panels b, c: lagging of the bottom part from the rest of the prominence body, and panels d, e demonstrate the zigzag shape of the motions. Panel f: state of the prominence in the final stage of simulation. The dashed lines indicate the magnetic field lines before the perturbation. The yellow lines represent the density isocontours.

In the text
thumbnail Fig. 3.

Longitudinal velocity variations at the position of the center of mass of the prominence for the simulations with different resolutions. The red and blue lines indicate the resolution of 240 km and 60 km, respectively.

In the text
thumbnail Fig. 4.

Temporal evolution of the longitudinal velocity (left panel) and vertical component of the transverse velocity (right panel) after the horizontal perturbation. The prominence model used has a shear angle of θ = 45° and a density contrast of χ = 200. Colored lines indicate velocity variations at the center of mass of different field lines. The colors denote the maximum initial density at the field lines according to the color bar scale. The right vertical axis indicates the height of the magnetic dips; the left vertical axis represents the scale of the maximum initial amplitude.

In the text
thumbnail Fig. 5.

Oscillation periods as a function of height of the magnetic dips. Top panels: models with θ = 45° and bottom panels: models with θ = 60°. Left and right panels: periods of the heavy and light prominence models, respectively. Red diamonds: longitudinal period; blue circles: vertical period due to the horizontal perturbation; green triangles: vertical period due to the vertical perturbation; and solid line: period from the theoretical pendulum model.

In the text
thumbnail Fig. 6.

Evolution of the density and magnetic field after the vertical perturbation. Panels b–f demonstrate the longitudinal motions at the bottom of the prominence due to compression and rarefaction. The dashed lines denote the magnetic field lines before the perturbation. The yellow lines represent the density isocontours.

In the text
thumbnail Fig. 7.

Temporal evolution of the longitudinal velocity (left panel) and vertical component of the transverse velocity (right panel) after the vertical perturbation. The prominence model used has a shear angle of θ = 45° and a density contrast of χ = 200. The colored lines indicate velocity variations at the center of mass (right) or the certain distance from the center of mass (left) at each field line. The colors indicate the maximum initial density at the field line according to the color bar scale. The right vertical axis indicates the height of the magnetic dips; the left vertical axis represents the scale of the maximum initial amplitude.

In the text
thumbnail Fig. 8.

Time-distance diagram of the vertical component of the transverse velocity along z-axis at x = 0, in the model with the shear angle of θ = 45° and the density contrast of χ = 200. Black lines denote a region where the prominence mass is located.

In the text
thumbnail Fig. 9.

Evolution of the magnetic field and gas pressure during the wavefront propagation. Panels a, b: initial gas pressure distribution; panels c–h: running difference of the two consecutive snapshots. Left-hand panels: same as the right-hand panels but focused at the location around the prominence. The dashed lines indicate the magnetic field lines before the perturbation. The yellow lines represent the density isocontours.

In the text
thumbnail Fig. 10.

Temporal evolution of the longitudinal velocity (left panel) and vertical component of the transverse velocity (right panel) after the wavefront propagation. The prominence model used has a shear angle of θ = 45° and a density contrast of χ = 200. The colored lines indicate velocity variations at the center of mass of different field lines. The colors indicate the maximum initial density at the field lines according to the color bar scale. The right vertical axis indicates the height of the magnetic dips; the left vertical axis represents the scale of the maximum initial amplitude.

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.