We present the parallel particle-in-cell (PIC) code Apar-T and, more importantly, address the fundamental question of the relations between the PIC model, the Vlasov-Maxwell theory, and real plasmas. First, we present four validation tests: spectra from simulations of thermal plasmas, linear growth rates of the relativistic tearing instability and of the filamentation instability, and non-linear filamentation merging phase. For the filamentation instability we show that the effective growth rates measured on the total energy can differ by more than 50% from the linear cold predictions and from the fastest modes of the simulation. Second, we detail a new method for initial loading of Maxwell-J\"uttner particle distributions with relativistic bulk velocity and relativistic temperature, and explain why the traditional method with individual particle boosting fails. Third, we scrutinize the question of what description of physical plasmas is obtained by PIC models. These models rely on two building blocks: coarse-graining, i.e., grouping of the order of p~10^10 real particles into a single computer superparticle, and field storage on a grid with its subsequent finite superparticle size. We introduce the notion of coarse-graining dependent quantities, i.e., quantities depending on p. They derive from the PIC plasma parameter Lambda^{PIC}, which we show to scale as 1/p. We explore two implications. One is that PIC collision- and fluctuation-induced thermalization times are expected to scale with the number of superparticles per grid cell, and thus to be a factor p~10^10 smaller than in real plasmas. The other is that the level of electric field fluctuations scales as 1/Lambda^{PIC} ~ p. We provide a corresponding exact expression. Fourth, we compare the Vlasov-Maxwell theory, which describes a phase-space fluid with infinite Lambda, to the PIC model and its relatively small Lambda.

GNU Data Language (GDL) is an open-source interpreted language aimed at numerical data analysis and visualisation. It is a free implementation of the Interactive Data Language (IDL) widely used in Astronomy. GDL has a full syntax compatibility with IDL, and includes a large set of library routines targeting advanced matrix manipulation, plotting, time-series and image analysis, mapping, and data input/output including numerous scientific data formats. We will present the current status of the project, the key accomplishments, and the weaknesses - areas where contributions are welcome !

Using particle-in-cell (PIC) simulations, we study the evolution of a strongly magnetized plasma slab propagating into a finite density ambient medium. Like previous work, we find that the slab breaks into discrete magnetic pulses. The subsequent evolution is consistent with diamagnetic relativistic pulse acceleration of \cite{liangetal2003}. Unlike previous work, we use the actual electron to proton mass ratio and focus on understanding trapping vs. transmission of the ambient plasma by the pulses and on the particle acceleration spectra. We find that the accelerated electron distribution internal to the slab develops a double-power law. We predict that emission from reflected/trapped external electrons will peak after that of the internal electrons. We also find that the thin discrete pulses trap ambient electrons but allow protons to pass through, resulting in less drag on the pulse than in the case of trapping of both species. Poynting flux dominated scenarios have been proposed as the driver of relativistic outflows and particle acceleration in the most powerful astrophysical jets.

The isotropization process of a collisionless plasma with an electron temperature anisotropy along an external magnetic field ($T_\| ^e\gg T_\perp^e$, $\|$ and $\perp$ with respect to the background magnetic field) and isotropic protons is investigated using a particle-in-cell(PIC) code. Restricting wave growth mainly parallel to the external magnetic field, the isotropization mechanism is identified to be the Electron Firehose Instability (EFI). The free energy in the electrons is first transformed into left-hand circularly polarized transverse low-frequency waves by a non-resonant interaction. Fast electrons can then be scattered towards higher perpendicular velocities by gyroresonance, leading finally to a complete isotropization of the velocity distribution. During this phase of the instability, Langmuir waves are generated which may lead to the emission of radio waves. A large fraction of the protons is resonant with the left-hand polarized electromagnetic waves, creating a proton temperature anisotropy $T_\|^p < T_\perp^p$. The parameters of the simulated plasma are chosen compatible to solar flare conditions. The results indicate the significance of this mechanism in the particle acceleration context: The EFI limits the anisotropy of the electron velocity distribution, and thus provides the necessary condition for further acceleration. It enhances the pitch-angle of the electrons and heats the ions.

A new acousto-optic radio spectrometer has observed the 1 - 2 GHz radio emission of solar flares with unprecedented sensitivity. The number of detected decimeter type III bursts is greatly enhanced compared to observations by conventional spectrometers observing only one frequency at the time. The observations indicate a large number of electron beams propagating in dense plasmas. For the first time, we report weak, reversed drifting type III bursts at frequencies above simultaneous narrowband decimeter spikes. The type III bursts are reliable signatures of electron beams propagating downward in the corona, apparently away from the source of the spikes. The observations contradict the most popular spike model that places the spike sources at the footpoints of loops. Conspicuous also was an apparent bidirectional type U burst forming a fish-like pattern. It occurs simultaneously with an intense U-burst at 600-370 MHz observed in Tremsdorf. We suggest that it intermodulated with strong terrestrial interference (cellular phones) causing a spurious symmetric pattern in the spectrogram at 1.4 GHz. Symmetric features in the 1 - 2 GHz range, some already reported in the literature, therefore must be considered with utmost caution.

The minimum and the mean bandwidth of individual narrowband spikes in two events in decimetric radio waves is determined by means of multi-resolution analysis. Spikes of a few tens of millisecond duration occur at decimetric/microwave wavelength in the particle acceleration phase of solar flares. A first method determines the dominant spike bandwidth scale based on their scalegram, the mean squared wavelet coefficient at each frequency scale. This allows to measure the scale bandwidth independently of heuristic spike selection criteria, e.g. manual selection. The major drawback is a low resolution in the bandwidth. To overcome this uncertainty, a feature detection algorithm and a criterion for spike shape in the time-frequency plane is applied to locate the spikes. In that case, the bandwidth is measured by fitting an assumed spike profile into the denoised data. The smallest FWHM bandwidth of spikes was found at 0.17 % and 0.41 % of the center frequency in the two events. Knowing the shortest relevant bandwidth of spikes, the slope of the Fourier power spectrum of this two events was determined and no resemblance to a Kolmogorov spectrum detected. Additionally the correlation between spike peak flux and bandwidth was examined.