Vortex structures in the Sun's chromosphere are believed to channel energy between different layers of the solar atmosphere. We investigate the nature and dynamics of two small-scale quiet-Sun rotating structures in the chromosphere. We analyse two chromospheric structures that show clear rotational patterns in spectropolarimetric observations taken with the Interferometric Bidimensional Spectrometer (IBIS) at the Ca II 8542 \AA~ line. We present the detection of spectropolarimetric signals that manifest the magnetic nature of rotating structures in the chromosphere. Our observations show two long-lived structures of plasma that each rotate clockwise inside a 10 arcsec$^{2}$~ quiet-Sun region. Their circular polarization signals are 5-10 times above the noise level. Line-of-sight Doppler velocity and horizontal velocity maps from the observations reveal clear plasma flows at and around the two structures. An MHD simulation shows these two structures are plausibly magnetically connected. Wave analysis suggests that the observed rotational vortex pattern could be due to a combination of slow actual rotation and a faster azimuthal phase speed pattern of a magneto-acoustic mode. Our results imply that the vortex structures observed in the Sun's chromosphere are magnetic in nature and that they can be connected locally through the chromosphere.

We use observations of quiet Sun (QS) regions in the H$\alpha$ 6563 \AA, Ca~{\sc ii} 8542 \AA\, and Fe~{\sc i} 6302 \AA\ lines. We observe brightenings in the wings of the H$\alpha$ and Ca~{\sc ii} combined with observations of the interacting magnetic concentrations observed in the Stokes signals of Fe~{\sc i}. These brightenings are similar to Ellerman bombs (EBs), i.e. impulsive bursts in the wings of the Balmer lines which leave the line cores unaffected. Such enhancements suggest that these events have similar formation mechanisms to the classical EBs found in active regions, with the reduced intensity enhancements found in the QS regions due to a weaker feeding magnetic flux. The observations also show that the quiet Sun Ellerman bombs (QSEBs) are formed at a higher height in the upper photosphere than the photospheric continuum level. Using simulations, we investigate the formation mechanism associated with the events and suggest that these events are driven by the interaction of magnetic field-lines in the upper photospheric regions. The results of the simulation are in agreement with observations when comparing the light-curves, and in most cases we found that the peak in the Ca~{\sc ii} 8542 \AA\ wing occurred before the peak in H$\alpha$ wing. Moreover, in some cases, the line profiles observed in Ca~{\sc ii} are asymmetrical with a raised core profile. The source of heating in these events is shown by the MURaM simulations and is suggested to occur 430 km above the photosphere.

Two of the central problems in our understanding of the solar chromosphere are how the upper chromosphere is heated and what drives spicules. Estimates of the required chromospheric heating, based on radiative and conductive losses suggest a rate of ${\sim} 0.1 \mathrm{\:erg\:cm^{-3}\:s^{-1}}$ in the lower chromosphere dropping to ${\sim} 10^{-3} \mathrm{\:erg\:cm^{-3}\:s^{-1}}$ in the upper chromosphere (\citet{Avrett1981}). The chromosphere is also permeated by spicules, higher density plasma from the lower atmosphere propelled upwards at speeds of ${\sim} 10-20 \mathrm{\:km\:s^{-1}}$, for so called Type-I spicules (\citet{Pereira2012,Zhang2012}), reaching heights of ${\sim} 3000-5000 \mathrm{\:km}$ above the photosphere. A clearer understanding of chromospheric dynamics, its heating and the formation of spicules, is thus of central importance to solar atmospheric science. For over thirty years it has been proposed that photospheric driving of MHD waves may be responsible for both heating and spicule formation. This paper presents results from a high-resolution MHD treatment of photospheric driven Alfv\'en and kink waves propagating upwards into an expanding flux tube embedded in a model chromospheric atmosphere. We show that the ponderomotive coupling from Alfv\'en and kink waves into slow modes generates shocks which both heat the upper chromosphere and drive spicules. These simulations show that wave driving of the solar chromosphere can give a local heating rate which matches observations and drive spicules consistent with Type-I observations all within a single coherent model.

Physical processes which may lead to solar chromospheric heating are analyzed using high-resolution 1.5D non-ideal MHD modelling. We demonstrate that it is possible to heat the chromospheric plasma by direct resistive dissipation of high-frequency Alfv\'en waves through Pedersen resistivity. However this is unlikely to be sufficient to balance radiative and conductive losses unless unrealistic field strengths or photospheric velocities are used. The precise heating profile is determined by the input driving spectrum since in 1.5D there is no possibility of Alfv\'en wave turbulence. The inclusion of the Hall term does not affect the heating rates. If plasma compressibility is taken into account, shocks are produced through the ponderomotive coupling of Alfv\'en waves to slow modes and shock heating dominates the resistive dissipation. In 1.5D shock coalescence amplifies the effects of shocks and for compressible simulations with realistic driver spectra the heating rate exceeds that required to match radiative and conductive losses. Thus while the heating rates for these 1.5D simulations are an overestimate they do show that ponderomotive coupling of Alfv\'en waves to sound waves is more important in chromospheric heating than Pedersen dissipation through ion-neutral collisions.

Srivastava et al. (2010) have observed a highly twisted coronal loop, which was anchored in AR10960 during the period 04:43 UT-04:52 UT on 4 June 2007. The loop length and radius are approximately 80 Mm and 4 Mm, with a twist of 11.5 $\pi$. These observations are used as initial conditions in a three dimensional nonlinear magnetohydrodynamic simulation with parallel thermal conduction included. The initial unstable equilibrium evolves into the kink instability, from which synthetic observables are generated for various high-temperature filters of SDO/AIA. These observables include temporal and spatial averaging to account for the resolution and exposure times of SDO/AIA images. Using the simulation results, we describe the implications of coronal kink instability as observables in SDO/AIA filters.

The highly dynamical, complex nature of the solar atmosphere naturally implies the presence of waves in a topologically varied magnetic environment. Here, the interaction of waves with topological features such as null points is inevitable and potentially important for energetics. The low resistivity of the solar coronal plasma implies that non-MHD effects should be considered in studies of magnetic energy release in this environment. This paper investigates the role of the Hall term in the propagation and dissipation of waves, their interaction with 2D magnetic X-points and the nature of the resulting reconnection. A Lagrangian remap shock-capturing code (Lare2d) is used to study the evolution of an initial fast magnetoacoustic wave annulus for a range of values of the ion skin depth in resistive Hall MHD. A magnetic null-point finding algorithm is also used to locate and track the evolution of the multiple null-points that are formed in the system. Depending on the ratio of ion skin depth to system size, our model demonstrates that Hall effects can play a key role in the wave-null interaction. In particular, the initial fast-wave pulse now consists of whistler and ion-cyclotron components; the dispersive nature of the whistler wave leads to (i) earlier interaction with the null, (ii) the creation of multiple additional, transient nulls and, hence, an increased number of energy release sites. In the Hall regime, the relevant timescales (such as the onset of reconnection and the period of the oscillatory relaxation) of the system are reduced significantly, and the reconnection rate is enhanced.

It is known from numerical simulations that thermal conduction along magnetic field lines plays an important role in the evolution of the kink instability in coronal loops. This study presents the observational signatures of the kink instability in long coronal loops when parallel thermal conduction is included. The 3D nonlinear magnetohydrodynamic equations are solved numerically to simulate the evolution of a coronal loop that is initially in an unstable equilibrium. The loop has length 80 Mm, width 8 Mm and an initial maximum twist of Phi = 11.5 pi, where Phi is a function of the radius. The initial loop parameters are obtained from a highly twisted loop observed in the TRACE 171 A waveband. Synthetic observables are generated from the data. These observables include spatial and temporal averaging to account for the resolution and exposure times of TRACE images. Parallel thermal conduction reduces the maximum local temperature by up to an order of magnitude. This means that different spectral lines are formed and different internal loop structures are visible with or without the inclusion of thermal conduction. However, the response functions sample a broad range of temperatures. The result is that the inclusion of parallel thermal conductivity does not have as large an impact on observational signatures as the order of magnitude reduction in the maximum temperature would suggest; the net effect is a blurring of internal features of the loop structure.

For a magnetic flux tube, or indeed any flux, to emerge into the Solar corona from the convection zone it must pass through the partially ionised layers of the lower atmosphere: the photosphere and the chromosphere. In such regions the ion-neutral collisions lead to an increased resistivity for currents flowing across magnetic field lines. This Cowling resistivity can exceed the Spitzer resistivity by orders of magnitude and in 2.5D simulations has been shown to be sufficient to remove all cross field current from emerging flux. Here we extend this modelling into 3D. Once again it is found that the Cowling resistivity removes perpendicular current. However the presence of 3D structure prevents the simple comparison possible in 2.5D simulations. With a fully ionised atmosphere the flux emergence leads to an unphysically low temperature region in the overlying corona, lifting of chromospheric material and the subsequent onset of the Rayleigh-Taylor instability. Including neutrals removes the low temperature region, lifts less chromospheric matter and shows no signs of the Rayleigh-Taylor instability. Simulations of flux emergence therefore should include such a neutral layer in order to obtain the correct perpendicular current, remove the Rayleigh-Taylor instability and get the correct temperature profile. In situations when the temperature is not important, i.e. when no simulated spectral emission is required, a simple model for the neutral layer is demonstrated to adequately reproduce the results of fully consistent simulations.

Certain regions of the solar atmosphere are at sufficiently low temperatures to be only partially ionised. The lower chromosphere contains neutral atoms, the existence of which greatly increases the efficiency of the damping of waves due to collisional friction momentum transfer. More specifically the Cowling conductivity can be up to 12 orders of magnitude smaller than the Spitzer value, so that the main damping mechanism in this region is due to the collisions between neutrals and positive ions. Using values for the gas density and temperature as functions of height taken from the VAL C model of the quiet Sun, an estimate is made for the dependance of the Cowling conductivity on height and strength of magnetic field. Using both analytic and numerical approaches the passage of Alfven waves over a wide spectrum through this partially ionised region is investigated. Estimates of the efficiency of this region in the damping of Alfven waves are made and compared for both approaches. We find that Alfven waves with frequencies above 0.6Hz are completely damped and frequencies below 0.01 Hz unaffected.

In the recent work of Nakariakov et al. (2004), it has been shown that the time dependences of density and velocity in a flaring loop contain pronounced quasi-harmonic oscillations associated with the 2nd harmonic of a standing slow magnetoacoustic wave. That model used a symmetric heating function (heat deposition was strictly at the apex). This left outstanding questions: A) is the generation of the 2nd harmonic a consequence of the fact that the heating function was symmetric? B) Would the generation of these oscillations occur if we break symmetry? C) What is the spectrum of these oscillations? Is it consistent with a 2nd spatial harmonic? The present work (and partly Tsiklauri et al. (2004b)) attempts to answer these important outstanding questions. Namely, we investigate the physical nature of these oscillations in greater detail: we study their spectrum (using periodogram technique) and how heat positioning affects the mode excitation. We found that excitation of such oscillations is practically independent of 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. We also were able to model successfully SUMER oscillations observed in hot coronal loops.

Low-frequency longitudinal oscillations of a flaring coronal loop are studied numerically. In the recent work of Nakariakov et al., Astron. and Astrophys., 414, L25-L28 (2004) it has been shown that the time dependences of density and velocity in a flaring loop contain well-pronounced quasi-harmonic oscillations associated with a 2nd harmonic of a standing slow magnetoacoustic wave. In this work we investigate physical nature of these oscillations in greater detail, namely, we study their spectrum (using periodogram technique) and how does heat positioning affects the mode excitation. We found that excitation of such oscillations practically independent of the positioning 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 equal to the quarter of the period, it varies along the loop and is time dependent, especially in the case of one footpoint (asymmetric) heating. KEYWORDS -- Sun: flares -- Sun: oscillations -- Sun: Corona -- Stars: flare -- Stars: oscillations -- Stars: coronae

A 1-D 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 X 10^{11} cm^{-3} or more) in the case of footpoint heating than in the case of apex heating (2.5 X 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 allows us 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 sigma_{t}, but scales directly with the heating rate E_{H0}. (e) The maximum obtained densities at the loop apex, n_e^{max}, increase with the heating rate E_{H0} and heating duration sigma_{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. KEYWORDS: Sun: Flares -- Sun: Activity -- Sun: Corona

Evolution of a coronal loop in response to an impulsive energy release is numerically modelled. It is shown that the loop density evolution curves exhibit quasi-periodic perturbations with the periods given approximately by the ratio of the loop length to the average sound speed, associated with the second standing harmonics of an acoustic wave. The density perturbations have a maximum near the loop apex. The corresponding field-aligned flows have a node near the apex. We suggest that the quasi-periodic pulsations with periods in the range 10-300 s, frequently observed in flaring coronal loops in the radio, visible light and X-ray bands, may be produced by the second standing harmonic of the acoustic mode.

We investigate the interaction of a plane, linearly polarized, Alfvenic pulse with a one-dimensional, perpendicular to the magnetic field, plasma density inhomogeneity in the strongly nonlinear regime. Our numerical study of the full MHD equations shows that: (i) Plasma density inhomogeneity substantially enhances (by about a factor of 2) the generation of longitudinal compressive waves. (ii) Attained maximal values of the generated transverse compressive perturbations are insensitive to the strength of the plasma density inhomogeneity, plasma beta and the initial amplitude of the Alfven wave. Typically, they reach about 40% of the initial Alfven wave amplitude. (iii) Attained maximal values of the generated relative density perturbations are within 20-40% for 0.5 <= beta <= 2.0. They depend upon plasma beta strongly; and scale almost linearly with the initial Alfven wave amplitude.

A new mechanism for the enhanced generation of compressible fluctuations by Alfven waves is presented. A strongly nonlinear regime of Alfven wave phase-mixing is numerically simulated in a one-dimensionally inhomogeneous plasma of finite temperature. It is found that the inhomogeneity of the medium determines the efficiency of nonlinear excitation of magnetoacoustic waves. The level of the compressible fluctuations is found to be higher (up to the factor of two) in inhomogeneous regions. The amplitude of the generated magnetoacoustic wave can reach up to 30% of the source Alfven wave amplitude, and this value is practically independent of the Alfven wave amplitude and the steepness of Alfven speed profile. The highest amplitudes of compressible disturbances are reached in plasmas with beta of about 0.5. The further growth of the amplitude of compressible fluctuations is depressed by saturation.

We investigate the interaction of magnetohydrodynamic waves with plasma density inhomogeneities. Our numerical study of the full MHD equations shows that: (A) Plasma density inhomogeneities are a source of non-linear generation of transverse compressive waves by a plane Alfven wave, and substantially enhance (by about a factor of 2) the generation of the longitudinal compressive waves. (B) Attained maximal values of the generated transverse compressive waves are insensitive to the strength of the plasma density inhomogeneity and the initial amplitude of the Alfven wave. (C) Efficiency of the generation depends weakly upon the plasma beta parameter. The maximum generated amplitude of transverse compressive wave, up to 32% of the initial Alfven wave amplitude, is reached for about beta=0.5. The results obtained demonstrate that plasma inhomogeneities enhance the efficiency of the non-linear wave coupling.

The interaction of a weakly nonlinear Alfv\'enic pulse with an Alfv\'en speed inhomogeneity in the direction perpendicular to the magnetic field is investigated. Identically to the phase mixing experienced by a harmonic Alfv\'en wave, sharp transverse gradients are generated in the pulse by the inhomogeneity. In the initial stage of the evolution of an initially plane Alfv\'enic pulse, the transverse gradients efficiently generate transversely propagating fast magnetoacoustic waves. However, high resolution full MHD numerical simulations of the developed stage of the pulse evolution show that the generation saturates due to destructive wave interference. It is shown that the weakly non-linear description of the generated fast magnetoacoustic wave is well described by the driven wave equation proposed in Nakariakov, Roberts & Murawski (1997), and a simple numerical code (2D MacCromack), which solves it with minimal CPU resources produces identical results to those obtained from the full MHD code ({\it Lare2d}, Arber et al. 2001). A parametric study of the phenomenon is undertaken, showing that, contrary to one's expectations, steeper inhomogenities of the Alfv\'en speed do not produce higher saturation levels of the fast wave generation. There is a certain optimal gradient of the inhomogeneity that ensures the maximal efficiency of the fast wave generation.

We present numerical studies of the nonlinear, resistive magnetohydrodynamic (MHD) evolution of coronal loops. For these simulations we assume that the loops carry no net current, as might be expected if the loop had evolved due to vortex flows. Furthermore the initial equilibrium is taken to be a cylindrical flux tube with line-tied ends. For a given amount of twist in the magnetic field it is well known that once such a loop exceeds a critical length it becomes unstableto ideal MHD instabilities. The early evolution of these instabilities generates large current concentrations. Firstly we show that these current concentrations are consistent with the formation of a current sheet. Magnetic reconnection can only occur in the vicinity of these current concentrations and we therefore couple the resistivity to the local current density. This has the advantage of avoiding resistive diffusion in regions where it should be negligible. We demonstrate the importance of this procedure by comparison with simulations based on a uniform resistivity. From our numerical experiments we are able to estimate some observational signatures for unstable coronal loops. These signatures include: the timescale of the loop brightening; the temperature increase; the energy released and the predicted observable flow speeds. Finally we discuss to what extent these observational signatures are consistent with the properties of transient brightening loops.