Flare-generated acoustic oscillations in solar and stellar coronal loops

Long period longitudinal oscillations of a flaring coronal loop are studied numerically. In the recent work of Nakariakov et al. (2004) it has been shown that the time dependence of density and velocity in a flaring loop contain pronounced quasi-harmonic oscillations associated with the 2nd harmonic of a standing slow magnetoacoustic wave. In this work we investigate the physical nature of these oscillations in greater detail, namely, their spectrum (using the periodogram technique) and how heat positioning affects mode excitation. We found that excitation of such oscillations is practically independent of the location of the heat deposition in the loop. Because of the change of the background temperature and density, the phase shift between the density and velocity perturbations is not exactly a quarter of the period; it varies along the loop and is time dependent, especially in the case of one footpoint (asymmetric) heating

Three dimensional particle-in-cell simulation of particle acceleration by circularly polarised inertial Alfven waves in a transversely inhomogeneous plasma

The process of particle acceleration by left-hand, circularly polarised inertial Alfven waves (IAW) in a transversely inhomogeneous plasma is studied using 3D particle-in-cell simulation. A cylindrical tube with, transverse to the background magnetic field, inhomogeneity scale of the order of ion inertial length is considered on which IAWs with frequency $0.3 \omega_{ci}$ are launched that are allowed to develop three wavelength. As a result time-varying parallel electric fields are generated in the density gradient regions which accelerate electrons in the parallel to magnetic field direction. Driven perpendicular electric field of IAWs also heats ions in the transverse direction. Such numerical setup is relevant for solar flaring loops and earth auroral zone. This first, 3D, fully-kinetic simulation demonstrates electron acceleration efficiency in the density inhomogeneity regions, along the magnetic field, of the order of 45% and ion heating, in the transverse to the magnetic field direction, of 75%. The latter is a factor of two times higher than the previous 2.5D analogous study and is in accordance with solar flare particle acceleration observations. We find that the generated parallel electric field is localised in the density inhomogeneity region and rotates in the same direction and with the same angular frequency as the initially launched IAW. Our numerical simulations seem also to suggest that the "knee" often found in the solar flare electron spectra can alternatively be interpreted as the Landau damping (Cerenkov resonance effect) of IAWs due to the wave-particle interactions.Comment: Physics of Plasmas, in-press, September 2012 issue, final accepted versio

Radiative hydrodynamic modeling of the Bastille-Day flare (14 July, 2000). I, Numerical simulations

A 1D loop radiative hydrodynamic model that incorporates the effects of gravitational stratification, heat conduction, radiative losses, external heat input, presence of helium, and Braginskii viscosity is used to simulate elementary flare loops. The physical parameters for the input are taken from observations of the Bastille-Day flare of 2000 July 14. The present analysis shows that: a) the obtained maximum values of the electron density can be considerably higher (4.2 × 10 11 cm −3 or more) in the case of footpoint heating than in the case of apex heating (2.5 × 10 11 cm −3); b) the average cooling time after the flare peak takes less time in the case of footpoint heating than in the case of apex heating; c) the peak apex temperatures are significantly lower (by about 10 MK) for the case of footpoint heating than for apex heating (for the same average loop temperature of about 30 MK). This characteristic would allow to discriminate between different heating positioning; d) in both cases (of apex and footpoint heating), the maximum obtained apex temperature T max is practically independent of the heating duration σ t , but scales directly with the heating rate E H0 ; e) the maximum obtained densities at the loop apex, n max e, increase with the heating rate E H0 and heating duration σ t for both footpoint and apex heating. In Paper II we will use the outputs of these hydrodynamic simulations, which cover a wide range of the parameter space of heating rates and durations, as an input for forward-fitting of the multi-loop arcade of the Bastille-day flare

Transverse oscillations of systems of coronal loops

We study the collective kinklike normal modes of a system of several cylindrical loops using the T-matrix theory. Loops that have similar kink frequencies oscillate collectively with a frequency which is slightly different from that of the individual kink mode. On the other hand, if the kink frequency of a loop is different from that of the others, it oscillates individually with its own frequency. Since the individual kink frequency depends on the loop density but not on its radius for typical 1 MK coronal loops, a coupling between kink oscillations of neighboring loops take place when they have similar densities. The relevance of these results in the interpretation of the oscillations studied by \citet{schrijver2000} and \citet{verwichte2004}, in which transverse collective loop oscillations seem to be detected, is discussed. In the first case, two loops oscillating in antiphase are observed; interpreting this motion as a collective kink mode suggests that their densities are roughly equal. In the second case, there are almost three groups of tubes that oscillate with similar periods and therefore their dynamics can be collective, which again seems to indicate that the loops of each group share a similar density. All the other loops seem to oscillate individually and their densities can be different from the rest

Transverse oscillations of a multi-stranded loop

We investigate the transverse oscillations of a line-tied multi-stranded coronal loop composed of several parallel cylindrical strands. First, the collective fast normal modes of the loop are found with the T-matrix theory. There is a huge quantity of normal modes with very different frequencies and a complex structure of the associated magnetic pressure perturbation and velocity field. The modes can be classified as bottom, middle, and top according to their frequencies and spatial structure. Second, the temporal evolution of the velocity and magnetic pressure perturbation after an initial disturbance are analyzed. We find complex motions of the strands. The frequency analysis reveals that these motions are a combination of low and high frequency modes. The complexity of the strand motions produces a strong modulation of the whole tube movement. We conclude that the presumed internal fine structure of a loop influences its transverse oscillations and so its transverse dynamics cannot be properly described by those of an equivalent monolithic loop.Comment: Accepted in Ap

Electron versus Proton Timing Delays in Solar Flares

Both electrons and ions are accelerated in solar flares and carry nonthermal energy from the acceleration site to the chromospheric energy loss site, but the relative amount of energy carried by electrons versus ions is subject of debate. In this {\sl Letter} we test whether the observed energy-dependent timing delays of 20-200 keV HXR emission can be explained in terms of propagating electrons versus protons. For a typical flare, we show that the timing delays of fast (\lapprox 1 s) {\sl HXR pulses} is consistent with time-of-flight differences of directly precipitating electrons, while the timing delays of the {\sl smooth HXR} flux is consistent with collisional deflection times of trapped electrons. We show that these HXR timing delays cannot be explained either by $\le 1$ MeV protons (as proposed in a model by Simnett \& Haines 1990), because of their longer propagation and trapping times, or by $\approx 40$ MeV protons (which have the same velocity as $\approx 20$ keV electrons), because of their longer trapping times and the excessive fluxes required to generate the HXRs. Thus, the HXR timing results clearly rule out protons as the primary generators of $\ge 20$ keV HXR emission.Comment: 7 pages, TEX type, AASTeX macros, 1 Figure, to appear in Astrophysical Journal Letters, accepted 1996 July 2

Three-dimensional coronal slow modes: toward three-dimensional seismology

On 2008 January 10, the twin Solar Terrestrial Relations Observatory (STEREO) A and B spacecraft conducted a high time cadence study of the solar corona with the Extreme UltraViolet Imager (EUVI) instruments with the aim of investigating coronal dynamics. Observations of the three-dimensional propagation of waves within active region coronal loops and a measurement of the true coronal slow mode speed are obtained. Intensity oscillations with a period of approximately 12 minutes are observed to propagate outwards from the base of a loop system, consistent with the slow magnetoacoustic mode. A novel analysis technique is applied to measure the wave phase velocity in the observations of the A and B spacecraft. These stereoscopic observations are used to infer the three-dimensional velocity vector of the wave propagation, with an inclination of 37 +- 6 deg to the local normal and a magnitude of 132 +- 9 and 132 +- 11 km s-1, giving the first measurement of the true coronal longitudinal slow mode speed, and an inferred temperature of 0.84 +- 12 MK and 0.84 +- 15 MK

Coronal loop seismology using multiple transverse loop oscillation harmonics

Context. TRACE observations (23/11/1998 06:35:57−06:48:43 UT) in the 171 Å bandpass of an active region are studied. Coronal loop oscillations are observed after a violent disruption of the equilibrium. Aims. The oscillation properties are studied to give seismological estimates of physical quantities, such as the density scale height. Methods. A loop segment is traced during the oscillation, and the resulting time series is analysed for periodicities. Results. In the loop segment displacement, two periods are found: 435.6 ± 4.5 s and 242.7 ± 6.4 s, consistent with the periods of the fundamental and 2nd harmonic fast kink oscillation. The small uncertainties allow us to estimate the density scale height in the loop to be 109 Mm, which is about double the estimated hydrostatical value of 50 Mm. Because a loop segment is traced, the amplitude dependence along the loop is found for each of these oscillations. The obtained spatial information is used as a seismological tool to give details about the geometry of the observed loop