A program has been designed for accurately generating a potential magnetic field on a staggered grid by extrapolating the magnetic field normal to the photospheric surface. The code first calculates a magnetic potential using the Green's function method and then uses a finite differencing scheme to calculate the magnetic field from the potential. A new finite differencing formula was derived which accounts for grid staggering, it is shown that this formula gives a numerical approximation that is closest to the real potential field. It is also shown that extending the region over which normal photospheric field is specified can improve the accuracy of the potential field produced. The program is a FORTRAN 90 code that can be used to generate potential magnetic field inputs for Lare3d and other MHD solvers that use a staggered grid for magnetic field components. The program can be parallelised to run quickly over multiple computing cores. The code and supporting description are provided in the appendices and can also be found at https://github.com/calboo/Potential_Field_Calculator.

We studied the dynamics of all prominence tornadoes detected by the Solar Dynamics Observatory/Atmospheric Imaging Assembly from 2011 January 01 to December 31. In total, 361 events were identified during the whole year, but only 166 tornadoes were traced until the end of their lifetime. Out of 166 tornadoes, 80 (48%) triggered CMEs in hosting prominences, 83 (50%) caused failed coronal mass ejections (CMEs) or strong internal motion in the prominences, and only 3 (2%) finished their lifetimes without any observed activity. Therefore, almost all prominence tornadoes lead to the destabilization of their hosting prominences and half of them trigger CMEs. Consequently, prominence tornadoes may be used as precursors for CMEs and hence for space weather predictions.

In some laboratory and most astrophysical situations plasma wake-field acceleration of electrons is one dimensional, i.e. variation transverse to the beam's motion can be ignored. Thus, one dimensional (1D), particle-in-cell (PIC), fully electromagnetic simulations of electron plasma wake field acceleration are conducted in order to study the differences in electron plasma wake field acceleration in MeV versus GeV and linear versus blowout regimes. First, we show that caution needs to be taken when using fluid simulations, as PIC simulations prove that an approximation for an electron bunch not to evolve in time for few hundred plasma periods only applies when it is sufficiently relativistic. This conclusion is true irrespective of the plasma temperature. We find that in the linear regime and GeV energies, the accelerating electric field generated by the plasma wake is similar to the linear and MeV regime. However, because GeV energy driving bunch stays intact for much longer time, the final acceleration energies are much larger in the GeV energies case. In the GeV energy range and blowout regime the wake's accelerating electric field is much larger in amplitude compared to the linear case and also plasma wake geometrical size is much larger. Thus, the correct positioning of the trailing bunch is needed to achieve the efficient acceleration. For the considered case, optimally there should be approximately $(90-100) c/\omega_{pe}$ distance between trailing and driving electron bunches in the GeV blowout regime.

Three dimensional, particle-in-cell, fully electromagnetic simulations of electron plasma wake field acceleration applicable to solar atmosphere are presented. It is established that injecting driving and trailing electron bunches into solar coronal and chromospheric plasmas, results in electric fields ($-(20-5) \times 10^{6}$ V/m), leading to acceleration of the trailing bunch up to 52 MeV, starting from initial 36 MeV. The results provide one of potentially important mechanisms for the extreme energetic solar flare electrons, invoking plasma wake field acceleration.

Previous theoretical considerations of electron beam relaxation in inhomogeneous plasmas have indicated that the effects of the irregular solar wind may account for the poor agreement of homogeneous modelling with the observations. Quasi-linear theory and Hamiltonian models based on Zakharov's equations have indicated that when a level of density fluctuations is above a given threshold, density irregularities act to de-resonate the beam-plasma interaction, restricting Langmuir wave growth on the expense of beam energy. This work presents the first fully kinetic particle-in-cell (PIC) simulations of beam relaxation under the influence of density irregularities. We aim to independently determine the influence of background inhomogeneity on the beam-plasma system, and to test theoretical predictions and alternative models using a fully kinetic treatment. We carry out 1D PIC simulations of a bump-on-tail unstable electron beam in the presence of increasing levels of background inhomogeneity using the fully electromagnetic, relativistic EPOCH PIC code. We find that in the case of homogeneous background plasma density, Langmuir wave packets are generated at the resonant condition and then quasi-liear relaxation leads to a dynamic increase of wavenumbers generated. No electron acceleration is seen - unlike in the inhomogeneous experiments, all of which produce high-energy electrons. For the inhomogeneous experiments we also observe the generation of backwards propagating Langmuir waves, which is shown directly to be due to the refraction of the packets off the density gradients. Our fully kinetic PIC simulations broadly confirm the findings of quasi-linear theory and the Hamiltonian model based on Zakharov's equations. Strong density fluctuations modify properties of excited Langmuir waves altering their dispersion properties.

In the previous works harmonic, phase-mixed, Alfven wave dynamics was considered both in the kinetic and magnetohydrodynamic regimes. Up today only magnetohydrodynamic, phase-mixed, Gaussian Alfven pulses were investigated. In the present work we extend this into kinetic regime. Here phase-mixed, Gaussian Alfven pulses are studied, which are more appropriate for solar flares, than harmonic waves, as the flares are impulsive in nature. Collisionless, phase-mixed, dispersive, Gaussian Alfven pulse in transversely inhomogeneous plasma is investigated by particle-in-cell (PIC) simulations and by an analytical model. The pulse is in inertial regime with plasma beta less than electron-to-ion mass ratio and has a spatial width of 12 ion inertial length. The linear analytical model predicts that the pulse amplitude decrease is described by the linear Korteweg de Vries (KdV) equation. The numerical and analytical solution of the linear KdV equation produces the pulse amplitude decrease in time as $t^{-1}$. The latter scaling law is corroborated by full PIC simulations. It is shown that the pulse amplitude decrease is due to dispersive effects, while electron acceleration is due to Landau damping of the phase-mixed waves. The established amplitude decrease in time as $t^{-1}$ is different from the MHD scaling of $t^{-3/2}$. This can be attributed to the dispersive effects resulting in the different scaling compared to MHD, where the resistive effects cause the damping, in turn, enhanced by the inhomogeneity. Reducing background plasma temperature and increase in ion mass yields more efficient particle acceleration.

In previous simulations of collisionless 2D magnetic reconnection it was consistently found that the term in the generalised Ohm's law that breaks the frozen-in condition is the divergence of the electron pressure tensor's non-gyrotropic components. A fully relativistic particle-in-cell (PIC) code was used to model $X$-point collapse with a guide-field in two and three spatial dimensions. We show that in a 2D $X$-point collapse with a guide-field close to the strength of the in-plane field, the increased induced shear flows along the diffusion region lead to a new reconnection regime in which electron inertial terms play a dominant role at the $X$-point. This transition is marked by the emergence of a magnetic island - and hence a second reconnection site - as well as electron flow vortices moving along the current sheet. The reconnection electric field at the $X$-point is shown to exceed all lower guide-field cases for a brief period, indicating a strong burst in reconnection. By extending the simulation to three spatial dimensions it is shown that the locations of vortices along the current sheet (visualised by their $Q$-value) vary in the out-of-plane direction, producing tilted vortex tubes. The vortex tubes on opposite sides of the diffusion region are tilted in opposite directions, similarly to bifurcated current sheets in oblique tearing-mode reconnection. The tilt angles of vortex tubes were compared to a theoretical estimation and were found to be a good match. Particle velocity distribution functions for different guide-field runs, for 2.5D and 3D simulations, are analysed and compared.

The motivation for this study is to include the effect of plasma flow in Alfv\'en wave (AW) damping via phase mixing and to explore the observational implications. Our magnetohydrodynamic (MHD) simulations and analytical calculations show that, when a background flow is present, mathematical expressions for the AW damping via phase mixing are modified by the following substitution: $C_A^\prime(x) \to C_A^\prime(x)+V_0^\prime(x)$, where $C_A$ and $V_0$ are AW phase and the flow speeds, and the prime denotes a derivative in the direction across the background magnetic field. In uniform magnetic fields and over-dense plasma structures, where $C_A$ is smaller than in the surrounding plasma, the flow, which is confined to the structure and going in the same direction as the AW, reduces the effect of phase-mixing, because on the edges of the structure $C_A^\prime$ and $V_0^\prime$ have opposite signs. Thus, the wave damps by means of slower phase-mixing compared to the case without the flow. This is the result of the co-directional flow that reduces the wave front stretching in the transverse direction. We apply our findings to addressing the question why over-dense solar coronal open magnetic field structures (OMFS) are cooler than the background plasma. Observations show that the over-dense OMFS (e.g. solar coronal polar plumes) are cooler than surrounding plasma and that, in these structures, Doppler line-broadening is consistent with bulk plasma motions, such as AW. If over-dense solar coronal OMFS are heated by AW damping via phase-mixing, we show that, co-directional with AW, plasma flow in them reduces the phase-mixing induced-heating, thus providing an explanation of why they appear cooler than the background.

Aims. The simulation of three-wave interaction based plasma emission, thought to be the underlying mechanism for Type III solar radio bursts, is a challenging task requiring fully-kinetic, multi-dimensional models. This paper aims to resolve a contradiction in past attempts, whereby some studies indicate that no such processes occur. Methods. We self-consistently simulate three-waved based plasma emission through all stages by using 2D, fully kinetic, electromagnetic particle-in-cell simulations of relaxing electron beams using the EPOCH2D code. Results. Here we present the results of two simulations; Run 1 (nb/n0 = 0.0057, vb/{\Delta}vb = vb/Ve = 16) and Run 2 (nb/n0 = 0.05, vb/{\Delta}vb = vb/Ve = 8), which we find to permit and prohibit plasma emission respectively. We show that the possibility of plasma emission is contingent upon the frequency of the initial electrostatic waves generated by the bump-in-tail instability, and that these waves may be prohibited from participating in the necessary three-wave interactions due to frequency conservation requirements. In resolving this apparent contradiction through a comprehensive analysis, in this paper we present the first self-consistent demonstration of fundamental and harmonic plasma emission from a single-beam system via fully kinetic numerical simulation. We caution against simulating astrophysical radio bursts using unrealistically dense beams (a common approach which reduces run time), as the resulting non-Langmiur characteristics of the initial wave modes significantly suppresses emission. Comparison of our results also indicates that, contrary to the suggestions of previous authors, an alternative plasma emission mechanism based on two counter-propagating beams is unnecessary in an astrophysical context.

J. Graf von der Pahlen and D. Tsiklauri, Phys. Plas. 21, 060705 (2014), established that the generation of octupolar out-of-plane magnetic field structure in a stressed $X$-point collapse is due to ion currents. The field has a central region, comprising of the well-known qaudrupolar field (quadrupolar components), as well as four additional poles of reversed polarity closer to the corners of the domain (octupolar components). In this extended work, the dependence of the octupolar structure on domain size and ion mass variation is investigated. Simulations show that the strength and spatial structure of the generated octupolar magnetic field is independent of ion to electron mass ratio. Thus showing that ion currents play a significant role in out-of-plane magnetic structure generation in physically realistic scenarios. Simulations of different system sizes show that the width of the octupolar structure remains the same and has a spacial extent of the order of the ion inertial length. The width of the structure thus appears to be independent on boundary condition effects. The length of the octupolar structure however increases for greater domain sizes, prescribed by the external system size. This was found to be a consequence of the structure of the in-plane magnetic field in the outflow region halting the particle flow and thus terminating the in-plane currents that generate the out-of-plane field. The generation of octupolar magnetic field structure is also established in a tearing-mode reconnection scenario. The differences in the generation of the octupolar field and resulting qualitative differences between $X$-point collapse and tearing-mode are discussed.

The out-of-plane magnetic field, generated by fast magnetic reconnection, during collisionless, stressed $X$-point collapse, was studied with a kinetic, 2.5D, fully electromagnetic, relativistic particle-in-cell numerical code, using both closed (flux conserving) and open boundary conditions on a square grid. It was discovered that the well known quadrupolar structure in the out-of-plane magnetic field gains four additional regions of opposite magnetic polarity, emerging near the corners of the simulation box, moving towards the $X$-point. The emerging, outer, magnetic field structure has opposite polarity to the inner quadrupolar structure, leading to an overall octupolar structure. Using Ampere's law and integrating electron and ion currents, defined at grid cells, over the simulation domain, contributions to the out-of-plane magnetic field from electron and ion currents were determined. The emerging regions of opposite magnetic polarity were shown to be the result of ion currents. Magnetic octupolar structure is found to be a signature of $X$-point collapse, rather than tearing mode, and factors relating to potential discoveries in experimental scenarios or space-craft observations are discussed.

Previous studies [Malara et al ApJ, 533, 523 (2000)] considered small-amplitude Alfven wave (AW) packets in Arnold-Beltrami-Childress (ABC) magnetic field using WKB approximation. In this work linearly polarised Alfven wave dynamics in ABC magnetic field via direct 3D MHD numerical simulation is studied for the first time. Gaussian AW pulse with length-scale much shorter than ABC domain length and harmonic AW with wavelength equal to ABC domain length are studied for four different resistivities. While it is found that AWs dissipate quickly in the ABC field, surprisingly, AW perturbation energy increases in time. In the case of the harmonic AW perturbation energy growth is transient in time, attaining peaks in both velocity and magnetic perturbation energies within timescales much smaller than resistive time. In the case of the Gaussian AW pulse velocity perturbation energy growth is still transient in time, attaining a peak within few resistive times, while magnetic perturbation energy continues to grow. It is also shown that the total magnetic energy decreases in time and this is governed by the resistive evolution of the background ABC magnetic field rather than AW damping. On contrary, when background magnetic field is uniform, the total magnetic energy decrease is prescribed by AW damping, because there is no resistive evolution of the background. By considering runs with different amplitudes and by analysing perturbation spectra, possible dynamo action by AW perturbation-induced peristaltic flow and inverse cascade of magnetic energy have been excluded. Therefore, the perturbation energy growth is attributed to a new instability. The growth rate appears to be dependent on the value of the resistivity and spatial scale of the AW disturbance. Thus, when going beyond WKB approximation, AW damping, described by full MHD equations, does not guarantee decrease of perturbation energy.

Particle-in-cell code, EPOCH, is used for studying features of the wave component evident to propagate backwards from the front of the non-gyrotropic, relativistic beam of electrons injected in the Maxwellian, magnetised background plasma with decreasing density profile. According to recent findings presented in Tsiklauri (2011), Schmitz & Tsiklauri (2013) and Pechhacker & Tsiklauri (2012), in a 1.5-dimensional magnetised plasma system, the non-gyrotropic beam generates freely escaping electromagnetic radiation with properties similar to the Type-III solar radio bursts. In this study the backwards propagating wave component evident in the perpendicular components of the elecromagnetic field in such a system is presented for the first time. Background magnetic field strength in the system is varied in order to prove that the backwards propagating wave's frequency, prescribed by the whistler wave dispersion relation, is proportional to the specified magnetic field. Moreover, the identified whistlers are shown to be generated by the normal Doppler-shifted relativistic resonance. Large fraction of the energy of the perpendicular electromagnetic field components is found to be carried away by the whistler waves, while a small but sufficient fraction is going into L- and R- electromagnetic modes.

A possible solution to the unexplained high intensity hard x-ray (HXR) emission observable during solar flares was investigated via 3D fully relativistic, electromagnetic particle-in-cell (PIC) simulations with realistic ion to electron mass ratio. A beam of accelerated electrons was injected into a magnetised, Maxwellian, homogeneous and inhomogeneous background plasma. The electron distribution function was unstable to the beam-plasma instability and was shown to generate Langmuir waves, while relaxing to plateau formation. In order to estimate the role of the background density gradient on an unbound (infinite spatial extent) beam, three different scenarios were investigated: a) a uniform density background; b) a weak density gradient, n_R/n_L=3; c) a strong gradient case, n_R/n_L=10, where n_R and n_L denote background electron densities on the left and right edges of the simulation box respectively. The strong gradient case produced the largest fraction of electrons beyond 15 v_th. Further, two cases (uniform and strong gradient background) with spatially localized beam injections were performed aiming to show drifts of the generated Langmuir wave wavenumbers, as suggested in previous studies. For the strong gradient case, the Langmuir wave power is shown to drift to smaller wavenumbers, as found in previous quasi-linear simulations.

Works of D. Tsiklauri, T. Haruki, Phys. of Plasmas, 15, 102902 (2008) and D. Tsiklauri and T. Haruki, Phys. of Plasmas, 14, 112905, (2007) are extended by inclusion of the out-of-plane magnetic (guide) field. In particular, magnetic reconnection during collisionless, stressed $X$-point collapse for varying out-of-plane guide-fields is studied using a kinetic, 2.5D, fully electromagnetic, relativistic particle-in-cell numerical code. Cases for both open and closed boundary conditions are investigated, where magnetic flux and particles are lost and conserved respectively. It is found that reconnection rates and out-of-plane currents in the $X$-point increase more rapidly and peak sooner in the closed boundary case, but higher values are reached in the open boundary case. The normalized reconnection rate is fast: 0.10-0.25. In the open boundary case an increase of guide-field yields later onsets in the reconnection peak rates, while in the closed boundary case initial peak rates occur sooner but are suppressed. The reconnection current increases for low guide-fields but then decreases similarly. In the open boundary case, for guide-fields of the order of the in-plane magnetic field, the generation of electron vortices occurs. Possible causes of the vortex generation, based on the flow of particles in the diffusion region and localized plasma heating, are discussed. Before peak reconnection onset, oscillations in the out-of-plane electric field at the $X$-point are found, ranging in frequency from approximately 1 to 2 $\omega_{pe}$ and coinciding with oscillatory reconnection. These oscillations are found to be part of a larger wave pattern. Mapping the out-of-plane electric field over time and applying 2D Fourier transforms reveals that the waves predominantly correspond to the ordinary mode and may correspond to observable radio waves such as solar radio burst fine structure spikes.

Extensive particle-in-cell simulations of fast electron beams injected in a background magnetised plasma with a decreasing density profile were carried out. These simulations were intended to further shed light on a newly proposed mechanism for the generation of electromagnetic waves in type III solar radio bursts [D. Tsiklauri, Phys. Plasmas, 18, 052903 (2011)]. The numerical simulations were carried out using different density profiles and fast electron distribution functions. It is shown that electromagnetic L and R modes are excited by the transverse current, initially imposed on the system. In the course of the simulations no further interaction of the electron beam with the background plasma could be observed.

A beam of super-thermal, hot electrons was injected into maxwellian plasma with a density gradient along a magnetic field line. 1.5D particle-in-cell simulations were carried out which established that the EM emission is produced by the perpendicular component of the beam injection momentum. The beam has a positive slope in the distribution function in perpendicular momentum phase space, which is the characteristic feature of a cyclotron maser. The cyclotron maser in the overdense plasma generates emission at the electron cyclotron frequency. The frequencies of generated waves were too low to propagate away from the injection region, hence the wavelet transform shows a pulsating wave generation and decay process. The intensity pulsation frequency is twice the relativistic cyclotron frequency. Eventually, a stable wave packet formed and could mode couple on the density gradient to reach frequencies of the order of the plasma frequency, that allowed for propagation. The emitted wave is likely to be a z-mode wave. The total electromagnetic energy generated is of the order of 0.1% of the initial beam kinetic energy. The proposed mechanism is of relevance to solar type III radio bursts, as well as other situations, when the injected electron beam has a non-zero perpendicular momentum, e.g. magnetron.

1.5D Particle-In-Cell simulations of a hot, low density electron beam injected into magnetized, maxwellian plasma were used to further explore the alternative non-gyrotropic beam driven electromagnetic emission mechanism, first studied in Tsiklauri (2011). Variation of beam injection angle and background density gradient showed that the emission process is caused by the perpendicular component of the beam injection current, whereas the parallel component only produces Langmuir waves, which play no role in the generation of EM waves in our mechanism. Particular emphasis was put on the case, where the beam is injected perpendicularly to the background magnetic field, as this turned off any electrostatic wave generation along the field and left a purely electromagnetic signal in the perpendicular components. The simulations establish the following key findings: i) Initially waves at a few w_ce/gamma are excited, mode converted and emitted at w_pe ii) The emission intensity along the beam axis is proportional to the respective component of the kinetic energy of the beam; iii) The frequency of the escaping EM emission is independent of the injection angle; iv) A stronger background density gradient causes earlier emission; v) The beam electron distribution function in phase space shows harmonic oscillation in the perpendicular components at the relativistic gyrofrequency; vi) The requirement for cyclotron maser emission, df/dv_perp > 0, is fulfilled; vii) The degree of linear polarization of the emission is strongly dependent on the beam injection angle; viii) The generated electromagnetic emission is left-hand elliptically polarized as the pitch angle tends to 90 deg; ix) The generated electromagnetic energy is of the order of 0.1% of the initial beam kinetic energy.

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.

Dispersive Alfven waves (DAWs) offer, an alternative to magnetic reconnection, opportunity to accelerate solar flare particles. We study the effect of DAW polarisation, L-, R-, circular and elliptical, in different regimes inertial and kinetic on the efficiency of particle acceleration. We use 2.5D PIC simulations to study how particles are accelerated when DAW, triggered by a solar flare, propagates in transversely inhomogeneous plasma that mimics solar coronal loop. (i) In inertial regime, fraction of accelerated electrons (along the magnetic field), in density gradient regions is ~20% by the time when DAW develops 3 wavelengths and is increasing to ~30% by the time DAW develops 13 wavelengths. In all considered cases ions are heated in transverse to the magnetic field direction and fraction of the heated particles is ~35%. (ii) The case of R-circular, L- and R- elliptical polarisation DAWs, with the electric field in the non-ignorable transverse direction exceeding several times that of in the ignorable direction, produce more pronounced parallel electron beams and transverse ion beams in the ignorable direction. In the inertial regime such polarisations yield the fraction of accelerated electrons ~20%. In the kinetic regime this increases to ~35%. (iii) The parallel electric field that is generated in the density inhomogeneity regions is independent of m_i/m_e and exceeds the Dreicer value by 8 orders of magnitude. (iv) Electron beam velocity has the phase velocity of the DAW. Thus electron acceleration is via Landau damping of DAWs. For the Alfven speeds of 0.3c the considered mechanism can accelerate electrons to energies circa 20 keV. (v) The increase of mass ratio from m_i/m_e=16 to 73.44 increases the fraction of accelerated electrons from 20% to 30-35% (depending on DAW polarisation). For the mass ratio m_i/m_e=1836 the fraction of accelerated electrons would be >35%.

1.5D PIC, relativistic, fully electromagnetic (EM) simulations are used to model EM wave emission generation in the context of solar type III radio bursts. The model studies generation of EM waves by a super-thermal, hot beam of electrons injected into a plasma thread that contains uniform longitudinal magnetic field and a parabolic density gradient. In effect, a single magnetic line connecting Sun to earth is considered, for which several cases are studied. (i) We find that the physical system without a beam is stable and only low amplitude level EM drift waves (noise) are excited. (ii) The beam injection direction is controlled by setting either longitudinal or oblique electron initial drift speed, i.e. by setting the beam pitch angle. In the case of zero pitch angle, the beam excites only electrostatic, standing waves, oscillating at plasma frequency, in the beam injection spatial location, and only low level EM drift wave noise is also generated. (iii) In the case of oblique beam pitch angles, again electrostatic waves with same properties are excited. However, now the beam also generates EM waves with the properties commensurate to type III radio bursts. The latter is evidenced by the wavelet analysis of transverse electric field component, which shows that as the beam moves to the regions of lower density, frequency of the EM waves drops accordingly. (iv) When the density gradient is removed, electron beam with an oblique pitch angle still generates the EM radiation. However, in the latter case no frequency decrease is seen. Within the limitations of the model, the study presents the first attempt to produce simulated dynamical spectrum of type III radio bursts in fully kinetic plasma model. The latter is based on 1.5D non-zero pitch angle (non-gyrotropic) electron beam, that is an alternative to the plasma emission classical mechanism.

The heating of solar chromospheric inter-network regions by means of the absorption of electromagnetic (EM) waves that originate from the photospheric blackbody radiation is studied in the framework of a plasma slab model. The absorption is provided by the electron-neutral collisions in which electrons oscillate in the EM wave field and electron-neutral collisions damp the EM wave. Given the uncertain nature of the collision cross-section due to the plasma micro-turbulence, it is shown that for plausible physical parameters, the heating flux produced by the absorption of EM waves in the chromosphere is between $20 - 45$ % of the chromospheric radiative loss flux requirement. It is also established that there is an optimal value for the collision cross-section, $5 \times 10^{-18}$ m$^{2}$, that produces the maximal heating flux of 1990 W m$^{-2}$.

Using the Galaxy as an example, we study the effect of j x B force on the rotational curves of gas and plasma in galaxies. Acceptable model for the galactic magnetic field and plausible physical parameters are used to fit the flat rotational curve for gas and plasma based on the observed baryonic (visible) matter distribution and j x B force term in the static MHD equation of motion. We also study the effects of varied strength of the magnetic field, its pitch angle and length scale on the rotational curves. We show that j x B force does not play an important role on the plasma dynamics in the intermediate range of distances 6-12 kpc from the centre, whilst the effect is sizable for larger r (r > 15 kpc), where it is the most crucial.

1.5D Vlasov-Maxwell simulations are employed to model electromagnetic emission generation in a fully self-consistent plasma kinetic model for the first time in the solar physics context. The simulations mimic the plasma emission mechanism and Larmor drift instability in a plasma thread that connects the Sun to Earth with the spatial scales compressed appropriately. The effects of spatial density gradients on the generation of electromagnetic radiation are investigated. It is shown that 1.5D inhomogeneous plasma with a uniform background magnetic field directed transverse to the density gradient is aperiodically unstable to Larmor-drift instability. The latter results in a novel effect of generation of electromagnetic emission at plasma frequency. When density gradient is removed (i.e. when plasma becomes stable to Larmor-drift instability) and a $low$ density, super-thermal, hot beam is injected along the domain, in the direction perpendicular to the magnetic field, plasma emission mechanism generates non-escaping Langmuir type oscillations which in turn generate escaping electromagnetic radiation. It is found that in the spatial location where the beam is injected, the standing waves, oscillating at the plasma frequency, are excited. These can be used to interpret the horizontal strips observed in some dynamical spectra. Quasilinear theory predictions: (i) the electron free streaming and (ii) the beam long relaxation time, in accord with the analytic expressions, are corroborated via direct, fully-kinetic simulation. Finally, the interplay of Larmor-drift instability and plasma emission mechanism is studied by considering $dense$ electron beam in the Larmor-drift unstable (inhomogeneous) plasma. http://www.maths.qmul.ac.uk/~tsiklauri/movie1.mpg * http://www.maths.qmul.ac.uk/~tsiklauri/movie2.mpg * http://www.maths.qmul.ac.uk/~tsiklauri/movie3.mpg

The solar corona, anomalously hot outer atmosphere of the Sun, is traditionally described by magnetohydrodynamic, fluid-like approach. Here we review some recent developments when, instead, a full kinetic description is used. It is shown that some of the main unsolved problems of solar physics, such as coronal heating and solar flare particle acceleration can be viewed in a new light when the small-scale, kinetic plasma description methods are used.

Based on the first principles (i.e. (i) by balancing the magnetic field advection with the term containing electron pressure tensor non-gyrotropic components in the generalised Ohm's law; (ii) using the conservation of mass; and (iii) assuming that the weak magnetic field region width, where electron meandering motion supports electron pressure tensor off-diagonal (non-gyrotropic) components, is of the order of electron Larmor radius) a simple model of magnetic reconnection in a collisionless regime is formulated. The model is general, resembling its collisional Sweet-Parker analogue in that it is not specific to any initial configuration e.g. Harris type tearing unstable current sheet, X-point collapse or otherwise. In addition to its importance from the fundamental point of view, the collisionless reconnection model offers a much faster reconnection rate (M_c'less=(c/omega_pe)^2 /(r_L,e L)) than Sweet-Parker's classical one (M_sp=S^-1/2). The width of the diffusion region (current sheet) in the collisionless regime is found to be delta_c'less=(c/omega_pe)^2/r_L,e, which is independent of global reconnection scale L and is only prescribed by micro-physics (electron inertial length, c/omega_pe, and electron Larmor radius, r_L,e). Amongst other issues, the fastness of the reconnection rate alleviates e.g. the problem of interpretation of solar flares by means of reconnection, as for the typical solar coronal parameters the obtained collisionless reconnection time can be a few minutes, as opposed to Sweet-Parker's equivalent value of < a day. The new theoretical reconnection rate is compared to the MRX device experimental data by [Yamada et al., Phys. Plasmas, 13, 052119 (2006), Ji et al. GRL, 35, 13106 (2008)] and a good agreement is obtained.

Previous studies of phase mixing of ion cyclotron (IC), Alfv\'enic, waves in the collisionless regime have established the generation of parallel electric field and hence acceleration of electrons in the regions of transverse density inhomogeneity. However, outstanding issues were left open. Here we use 2.5D, relativistic, fully electromagnetic PIC (Particle-In-Cell) code and an analytic MHD (Magnetohydrodynamic) formulation, to establish the following points: (i) Using the generalised Ohm's law we find that the parallel electric field is supported mostly by the electron pressure tensor, with a smaller contribution from the electron inertia term. (ii) The generated parallel electric field and the fraction of accelerated electrons are independent of the IC wave frequency remaining at a level of six orders of magnitude larger than the Dreicer value and approximately 20% respectively. The generated parallel electric field and the fraction of accelerated electrons increase with the increase of IC wave amplitude. The generated parallel electric field seems to be independent of plasma beta, while the fraction of accelerated electrons strongly increases with the decrease of plasma beta (for plasma beta of 0.0001 the fraction of accelerated electrons can be as large as 47%). (iii) In the collisionless regime IC wave dissipation length (that is defined as the distance over which the wave damps) variation with the driving frequency shows a deviation from the analytical MHD result, which we attribute to a possible frequency dependence of the effective resistivity. (iv) Effective anomalous resistivity, inferred from our numerical simulations, is at least four orders of magnitude larger than the classical Spitzer value.

Recently, magnetic reconnection during collisionless, stressed, X-point collapse was studied using kinetic, 2.5D, fully electromagnetic, relativistic Particle-in-Cell numerical code [D. Tsiklauri and T. Haruki, Phys. Plasmas, 14, 112905 (2007)]. Here we finalise the investigation of this topic by addressing key outstanding physical questions. It has been established here that: (i) reconnection out-of-plane electric field at the magnetic null is generated by the electron pressure tensor off-diagonal terms, resembling to the case of tearing unstable Harris current sheet studied by the GEM reconnection challenge; (ii) For m_i / m_e >> 1 the time evolution of the reconnected flux is independent of ion-electron mass ratio; also, in the case of m_i / m_e = 1 we show that reconnection proceeds slowly as the Hall term is zero; when m_i / m_e >> 1 (i.e. the Hall term is non-zero) reconnection is fast and we conjecture that this is due to magnetic field being frozen into electron fluid, which moves significantly faster than ion fluid; (iii) within one Alfven time, somewhat less than half (40%) of the initial total (roughly magnetic) energy is converted into the kinetic energy of electrons, and somewhat more than half (60%) into kinetic energy of ions (similar to solar flare observations); (iv) in the strongly stressed X-point case, in about one Alfven time, a full isotropy in all three spatial directions of the velocity distribution is seen for super-thermal electrons (also commensurate to solar flare observations).

The propagation of shear Alfven wave packets in inhomogeneous magnetic fields, at the origin of their distortion, regardless of the occurrence of non-linear coupling, is considered. It is shown that the distortion mechanism can be regarded as mixing process and hence, standard "phase mixing" corresponds to the effect of an "Alfvenic" shear flow while enhanced dissipation at a magnetic X-point corresponds to mixing by an "Alfvenic" strain flow. The evolution of the global wave field is supposed to result from the dynamics of a superposition of wave packets and a kinetic equation for the wave energy is obtained following this eikonal (WKB) description. Since shear Alfven wave packets experience continuous shearing/straining while transported by an inhomogeneous Alfvenic flow $\mathbf{V}_{A}$, their mixing process, in physical space, is also a cascade of wave energy in k-space. The wave energy spectrum resulting from this linear mechanism of energy transfer is determined for the special case of waves propagating along chaotic magnetic field lines, the analog of a chaotic mixing process. The latter follows a $k^{-1}$ power-law, in the energy conserving range in $k$ space.

The generation of parallel electric fields by the propagation of ion cyclotron waves in the plasma with a transverse density inhomogeneity was studied. It was proven that the minimal model required to reproduce the previous kinetic simulation results of E_{||} generation [Tsiklauri et al 2005, Genot et al 2004] is the two-fluid, cold plasma approximation in the linear regime. By considering the numerical solutions it was also shown that the cause of E_{||} generation is the electron and ion flow separation induced by the transverse density inhomogeneity. We also investigate how E_{||} generation is affected by the mass ratio and found that amplitude attained by E_{||} decreases linearly as inverse of the mass ratio m_i/m_e. For realistic mass ratio of m_i/m_e=1836, such empirical scaling law, within a time corresponding to 3 periods of the driving ion cyclotron wave, is producing E_{||}=14 Vm^{-1} for solar coronal parameters. Increase in mass ratio does not have any effect on final parallel (magnetic field aligned) speed attained by electrons. However, parallel ion velocity decreases linearly with inverse of the mass ratio m_i/m_e. These results can be interpreted as following: (i) ion dynamics plays no role in the E_{||} generation; (ii) E_{||} \propto 1/m_i scaling is caused by the fact that omega_d = 0.3 omega_{ci} \propto 1/m_i is decreasing with the increase of ion mass, and hence the electron fluid can effectively "short-circuit" (recombine with) the slowly oscillating ions, hence producing smaller E_{||}.

Corrected analytical solutions describing the enhanced phase mixing of Alfven waves propagating in divergent stratified coronal structures are presented. These show that the enhanced phase mixing mechanism can dissipate Alfven waves at heights less than half that is predicted by the previous analytical solutions. The enhanced phase mixing of 0.1 Hz harmonic Alfven waves propagating in strongly divergent, H_b=5 Mm, stratified coronal structures, H_rho=50 Mm, can fulfill 100% of an active region heating requirement, by generating viscous heating fluxes of F_H~2100 J /m^2 /s. The Alfven waves in this configuration are fully dissipated within 20 Mm, which is six times lower than would occur as a result of standard phase mixing in uniform magnetic fields. This results in the heating scale height, s_H, being lowered by a factor of six, to less than half of an active regions density scale height. Using the corrected analytical solutions it was found that, for a given wave frequency, the generation of a heating scale height of s_H<=50 Mm, by enhanced phase mixing in strongly divergent magnetic fields, requires a shear viscosity eight orders of magnitude lower, than required by standard phase mixing in uniform magnetic fields. It was also found that the enhanced phase mixing of observable, 0.01 rads /s Alfven waves, in strongly divergent magnetic fields, H_b=5 Mm, can generate heating scale heights within a density scale height, H_rho=50 Mm, using classical Braginskii viscosity. It is therefore not necessary to invoke anomalous viscosity in corona, if phase mixing takes place in strongly divergent magnetic fields. This study shows that the importance of enhanced phase mixing as a mechanism for dissipating Alfven waves in the solar corona (a stratified and divergent medium), has been seriously underestimated.

Two cases of weakly and strongly stressed X-point collapse were considered. Here descriptors weakly and strongly refer to 20 % and 124 % unidirectional spatial compression of the X-point, respectively. In the weakly stressed case, the reconnection rate, defined as the out-of-plane electric field in the X-point (the magnetic null) normalised by the product of external magnetic field and Alfv\'en speeds, peaks at 0.11, with its average over 1.25 Alfv\'en times being 0.04. Electron energy distribution in the current sheet, at the high energy end of the spectrum, shows a power law distribution with the index varying in time, attaining a maximal value of -4.1 at the final simulation time step (1.25 Alfv\'en times). In the strongly stressed case, magnetic reconnection peak occurs 3.4 times faster and is more efficient. The peak reconnection rate now attains value 2.5, with the average reconnection rate over 1.25 Alfv\'en times being 0.5. The power law energy spectrum for the electrons in the current sheet attains now a steeper index of -5.5, a value close to the ones observed in the vicinity of X-type region in the Earth's magneto-tail. Within about one Alfv\'en time, 2% and 20% of the initial magnteic energy is converted into heat and accelerated particle energy in the case of weak and strong stress, respectively. In the both cases, during the peak of the reconnection, the quadruple out-of-plane magnetic field is generated, hinting possibly to the Hall regime of the reconnection. These results strongly suggest the importance of the collionless, stressed X-point collapse as a possible contributing factor to the solution of the solar coronal heating problem or more generally, as an efficient mechanism of converting magnetic energy into heat and super-thermal particle energy.

We study the generation of parallel electric fields by virtue of propagation of ion cyclotron waves (with frequency 0.3 w_{ci}) in the plasma with a transverse density inhomogeneity. Using two-fluid, cold plasma linearised equations, we show for the first time that E_|| generation can be understood by an analytic equation that couples E_|| to the transverse electric field of the driving ion cyclotron wave. In this simplified model, the generated E_|| amplitude e.g. for plausible solar coronal parameters attains values of 90 V/m for the mass ratio m_i/m_e=262, within a time corresponding to 3 periods of the driving ion cyclotron wave. By considering the numerical solutions we also show that the cause of E_|| generation is electron and ion flow separation (which is not the same as electrostatic charge separation) induced by the transverse density inhomogeneity. The model also correctly reproduces the previous kinetic results in that only electrons are accelerated (along the background magnetic field), while ions do not accelerate substantially. We also investigate how E_|| generation is affected by the mass ratio and found that amplitude attained by E_|| decreases linearly as inverse of the mass ratio m_i/m_e, i.e. E_|| \propto 1/m_i. This result contradicts to the earlier suggestion by Genot et al (1999, 2004) that the cause of E_|| generation is the polarisation drift of the driving wave, which scales as \propto m_i. Also, for realistic mass ratio of m_i/m_e=1836 our empirical scaling law is producing E_||=14 V/m (for solar coronal parameters).

Using Particle-In-Cell simulations i.e. in the kinetic plasma description the discovery of a new mechanism of parallel electric field generation was recently reported. Here we show that the electric field generation parallel to the uniform unperturbed magnetic field can be obtained in a much simpler framework using the ideal magnetohydrodynamics (MHD) description. In ideal MHD the electric field parallel to the uniform unperturbed magnetic field appears due to fast magnetosonic waves which are generated by the interaction of weakly non-linear Alfv\'en waves with the transverse density inhomogeneity. In the context of the coronal heating problem a new {\it two stage mechanism} of plasma heating is presented by putting emphasis, first, on the generation of parallel electric fields within an {\it ideal MHD} description directly, rather than focusing on the enhanced dissipation mechanisms of the Alfv\'en waves and, second, dissipation of these parallel electric fields via {\it kinetic} effects. It is shown that for a single Alfv\'en wave harmonic with frequency $\nu = 7$ Hz, and longitudinal wavelength $\lambda_A = 0.63$ Mm for a putative Alfv\'en speed of 4328 km s$^{-1}$, the generated parallel electric field could account for 10% of the necessary coronal heating requirement. We conjecture that wide spectrum (10$^{-4}-10^3$ Hz) Alfv\'en waves, based on the observationally constrained spectrum, could provide the necessary coronal heating requirement. By comparing MHD versus kinetic results we also show that there is a clear indication of the {\it anomalous resistivity} which is 100s of times greater than the classical Braginskii value.

Recently using Particle-In-Cell simulations i.e. in the kinetic plasma description Tsiklauri et al. and G\'enot et al. reported on a discovery of a new mechanism of parallel electric field generation, which results in electron acceleration. In this work we show that the parallel (to the uniform unperturbed magnetic field) electric field generation can be obtained in much simpler framework using ideal Magnetohydrodynamic (MHD) description, i.e. without resorting to complicated wave particle interaction effects such as ion polarisation drift and resulting space charge separation which seems to be an ultimate cause of the electron acceleration. Further, in the context of the coronal heating problem a new {\it two stage mechanism} of the plasma heating is presented by putting emphasis, first, on the generation of parallel electric fields within {\it ideal MHD} description directly, rather than focusing on the enhanced dissipation mechanisms of the Alfv\'en waves and, second, dissipation of these parallel electric fields via {\it kinetic} effects. It is shown that a single Alfv\'en wave harmonic with frequency ($\nu = 7$ Hz), (which has longitudinal wavelength $\lambda_A = 0.63$ Mm for putative Alfv\'en speed of 4328 km s$^{-1}$) the generated parallel electric field could account for the 10% of the necessary coronal heating requirement. We conjecture that wide spectrum (10$^{-4}-10^3$ Hz) Alfv\'en waves, based on observationally constrained spectrum, could provide necessary coronal heating requirement.

We solve numerically ideal, 2.5D, MHD equations in Cartesian coordinates, with a plasma beta of 0.0001 starting from the equilibrium that mimics a footpoint of a large curvature radius solar coronal loop or a polar region plume. On top of such an equilibrium, a purely Alfv\'enic, linearly polarised, plane wave is launched. In the context of the coronal heating problem a new two stage mechanism of plasma heating is presented by putting emphasis, first, on the generation of parallel electric fields within an ideal MHD description directly, rather than focusing on the enhanced dissipation mechanisms of the Alfv\'en waves and, second, dissipation of these parallel electric fields via {\it kinetic} effects. It is shown that a single Alfv\'en wave harmonic with frequency $\nu = 7$ Hz and longitudinal wavelength $\lambda_A = 0.63$ Mm, for a putative Alfv\'en speed of 4328 km s$^{-1}$, the generated parallel electric field could account for 10% of the necessary coronal heating requirement. We conjecture that wide spectrum (10$^{-4}-10^3$ Hz) Alfv\'en waves, based on the observationally constrained spectrum, could provide the necessary coronal heating requirement. The exact amount of energy that could be deposited by such waves through our mechanism of parallel electric field generation can only be calculated once a more complete parametric study is done. Thus, the "theoretical spectrum" of the energy stored in parallel electric fields versus frequency needs to be obtained.

We investigate a possibility of heating of the loops and other closed magnetic structures in active regions of the solar corona by the flow of solar wind (plus other flows that may be present) across the magnetic field lines (that are perpendicular to the flow), in a similar manner as a conventional MHD generator works. A simple model is formulated which allows to calculate the typical currents generated by this mechanism. This enables to make a judgement whether the loops in solar active regions can be heated by the proposed mechanism through the associated Ohmic dissipation. It is shown that even if the heating region width is of the order of a few proton Larmor radii (this effectively means that a loop is represented by {\it nano-threads}, by analogy with the nano-flare coronal heating mechanism, with width of a few proton Larmor radii) and the plasma flow speed is 50 km s$^{-1}$ then only about 1% of the heating requirement of coronal loops could be met under this conjecture. Other implications of the model, such as a predicted spatial heat distribution along the loop as well as application of the model to the other places e.g. chromosphere and corotating interaction region structures in the solar wind are discussed.

Particle-In-Cell (kinetic) simulations of shear Alfv\'en wave interaction with one dimensional, across the uniform magnetic field, density inhomogeneity (phase mixing) in collisionless plasma were performed for the first time. As a result a new electron acceleration mechanism is discovered. Progressive distortion of Alfv\'en wave front, due to the differences in local Alfv\'en speed, generates nearly parallel to the magnetic field electrostatic fields, which accelerate electrons via Landau damping. Surprisingly, amplitude decay law in the inhomogeneous regions, in the kinetic regime, is the same as in the MHD approximation described by Heyvaerts and Priest (1983).

In this work we used Particle-In-Cell simulations to study the interaction of circularly polarised Alfv\'en waves with one dimensional plasma density inhomogeneities transverse to the uniform magnetic field (phase mixing) in collisionless plasmas. In our preliminary work we reported discovery of a new electron acceleration mechanism, in which progressive distortion of the Alfv\'en wave front, due to the differences in local Alfv\'en speed, generates an oblique (nearly parallel to the magnetic field) electrostatic field. The latter accelerates electrons through the Landau resonance. Here we report a detailed study of this novel mechanism, including: (i) analysis of broadening of the ion distribution function due to the presence of Alfv\'en waves and (ii) the generation of compressive perturbations due to both weak non-linearity and plasma density inhomogeneity. The amplitude decay law in the inhomogeneous regions, in the kinetic regime, is demonstrated to be the same as in the MHD approximation described by Heyvaerts and Priest (1983).

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.

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

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

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.

Phase mixing of a three dimensional magnetohydrodynamic (MHD) pulse is studied in the compressive, three-dimensional (without an ignorable coordinate) regime. It is shown that the efficiency of decay of an Alfvenic part of a compressible MHD pulse is related linearly to the degree of non-planeness of the pulse in the homogeneous transverse direction. In the developed stage of phase mixing (for large times), coupling to its compressive part does not alter the power-law decay of an Alfvenic part of a compressible MHD pulse. The same applies to the dependence upon the resistivity of the Alfvenic part of the pulse. All this implies that the dynamics of Alfven waves can still be qualitatively understood in terms of the previous 2.5D models. Thus, the phase mixing remains a relevant paradigm for the coronal heating applications in the realistic 3D geometry and compressive plasma.

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.

Kink modes of solar coronal structures, perturbing the loop in the direction along the line-of-sight (LOS), can be observed as emission intensity disturbances propagating along the loop provided the angle between the LOS and the structure is not ninety degrees. The effect is based upon the change of the column depth of the loop (along the LOS) by the wave. The observed amplitude of the emission intensity variations can be larger than the actual amplitude of the wave by a factor of two and there is an optimal angle maximizing the observed amplitude. For other angles this effect can also attenuate the observed wave amplitude. The observed amplitude depends upon the ratio of the wave length of kink perturbations to the width of the structure and on the angle between the LOS and the axis of the structure. Sausage modes are always affected negatively from the observational point of view, as the observed amplitude is always less than the actual one. This effect should be taken into account in the interpretation of wave phenomena observed in the corona with space-borne and ground-based imaging telescopes.

Interaction of impulsively generated MHD waves with a one-dimensional plasma inhomogeneity, transverse to the magnetic field, is considered in the three-dimensional regime. Because of the transverse inhomogeneity, MHD fluctuations, even if they do not include initially any density perturbation, evolve toward states where the compressible components tend to become predominant. The propagating MHD pulse asymptotically reaches a quasi-steady state with the final levels of density perturbation weakly depending on the degree of non-planeness of the pulse in the homogeneous transverse direction and somewhat stronger depending on plasma beta. Our study demonstrates the necessity of incorporation of compressible and 3D effects in theory of Alfven wave phase mixing. However, as far as the dynamics of weakly non-plane Alfven waves is concerned it can still be qualitatively understood in terms of the previous 2.5D models.

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.

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.

A model interpreting variations of EUV brightness upward propagating in solar coronal loops as slow magnetoacoustic waves is developed. A loop is considered to have a non-zero plane inclination angle and offset of circular loop centre from the baseline. The model also incorporates effects of dissipation and gravitational stratification. A linear evolutionary equation is derived and applied to investigation of slow wave dynamics. Both the non-zero plane inclination and the offset decrease the growth of the wave amplitude due to stratification. It is shown that wide spectrum slow magnetoacoustic waves, consistent with currently available observations in the low frequency part of the spectrum, can provide the rate of heat deposition sufficient to heat the loop. In this scenario, the heat would be deposited near the loop footpoints which agrees with the current observational data.

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.

Linear stability of an isothermal, pressure-bounded, self-gravitating gas slab which is gravitationally coupled with the background weakly interacting massive particles (WIMPs) is investigated. Analytic dispersion relations describing such a configuration are derived. Two novel, distinct oscillatory modes are found. Astrophysical implications of the results are discussed.

It has been recently shown (Munyaneza, Tsiklauri and Viollier 1998) that the analysis of the orbits of the fast moving stars close to Sgr A$^{*}$ (Eckart and Genzel 1997) provides a valuable tool to probe the gravitational potential near the Galactic center. As an example, we present here the results on a calculation of possible orbits of the star S0-1 in both, the black hole and degenerate neutrino ball scenarios of the central mass, based on the recent measurements of stellar proper motions at the Galactic center by Ghez et al. 1998. Taking into account the error bars of their analysis, it is shown that within a few years time, the orbit of S0-1 may indeed reveal the nature of the supermassive compact dark object at the Galactic center.

We report here on a calculation of thermalization time-scale of the two temperature advection dominated accretion flow (ADAF) model. It is established that time required to equalize the electron and ion temperatures via electron-ion collisions in the ADAF with plausible physical parameters greatly exceeds age of the Universe, which corroborates validity one of the crucial assumptions of the ADAF model, namely the existence of a hot two temperature plasma. This work is motivated by the recent success (Mahadevan 1998a,b) of ADAF model in explaining the emitted spectrum of Sgr A*.

We report here on a calculation of possible orbits of the fast moving infrared source S1 which has been recently observed by Eckart and Genzel (1997) near the Galactic center. It is shown that tracking of the orbit of S1 or any other fast moving star near Sgr A^* offers a possibility of distinguishing between the supermassive black hole and extended object scenarios of Sgr A^*. In our calculations we assumed that the extended object at the Galactic center is a non-baryonic ball made of degenerate, self-gravitating heavy neutrino matter, as it has been recently proposed by Tsiklauri & Viollier (1998a,b).

It has been recently shown (Tsiklauri & Viollier, 1998a) that the matter concentration inferred from observed stellar motion at the galactic center (Eckart & Genzel, 1997, MNRAS, 284, 576 and Genzel et al., 1996, ApJ, 472, 153) is consistent with a supermassive object of $2.5 \times 10^6$ solar masses, composed of self-gravitating, degenerate heavy neutrinos. It has been furthermore suggested (Tsiklauri & Viollier, 1998a) that the neutrino ball scenario may have an advantage that it could possibly explain the so-called "blackness problem" of the galactic center. Here, we present a quantitative investigation of this statement, by calculating the emitted spectrum of Sgr A$^*$ in the framework of standard accretion disk theory.

Criterion of the Jeans instability of interstellar gas clouds which are gravitationally coupled with weakly interacting massive particles is revisited. It is established that presence of the dark matter always reduces the Jeans length, and in turn, Jeans mass of the interstellar gas clouds. Astrophysical implications of this effect are discussed.

It is shown that the matter concentration observed through stellar motion at the galactic center (Eckart & Genzel, 1997, MNRAS, 284, 576 and Genzel et al., 1996, ApJ, 472, 153) is consistent with a supermassive object of $2.5 \times 10^6$ solar masses composed of self-gravitating, degenerate heavy neutrinos, as an alternative to the black hole interpretation. According to the observational data, the lower bounds on possible neutrino masses are $m_\nu \geq 12.0$ keV$/c^2$ for $g=2$ or $m_\nu \geq 14.3$ keV$/c^2$ for $g=1$, where $g$ is the spin degeneracy factor. The advantage of this scenario is that it could naturally explain the low X-ray and gamma ray activity of Sgr A$^*$, i.e. the so called "blackness problem" of the galactic center.