We present a Vlasov-DArwin numerical code (ViDA) specifically designed to address plasma physics problems, where small-scale high accuracy is requested even during the non linear regime to guarantee a clean description of the plasma dynamics at fine spatial scales. The algorithm provides a low-noise description of proton and electron kinetic dynamics, by splitting in time the multi-advection Vlasov equation in phase space. Maxwell equations for the electric and magnetic fields are reorganized according to Darwin approximation to remove light waves. Several numerical tests show that ViDA successfully reproduces the propagation of linear and nonlinear waves and captures the physics of magnetic reconnection. We also discuss preliminary tests of the parallelization algorithm efficiency, performed at CINECA on the Marconi-KNL cluster. ViDA will allow to run Eulerian simulations of a non-relativistic fully-kinetic collisionless plasma and it is expected to provide relevant insights on important problems of plasma astrophysics such as, for instance, the development of the turbulent cascade at electron scales and the structure and dynamics of electron-scale magnetic reconnection, such as the electron diffusion region.

Aims. We investigate the role played by a pre-supernova (SN) ambient magnetic field on the dynamics of the expanding remnant of SN 1987A and the origin and evolution of the radio emission from the remnant, in particular, during the interaction of the blast wave with the nebula surrounding the SN. Methods. We model the evolution of SN 1987A from the breakout of the shock wave at the stellar surface to the expansion of its remnant through the surrounding nebula by 3D MHD simulations. The model considers the radiative cooling, the deviations from equilibrium of ionization, the deviation from temperature-equilibration between electrons and ions, and a plausible configuration of the pre-SN ambient magnetic field. From the simulations, we synthesize both thermal X-ray and non-thermal radio emission and compare the model results with observations. Results. The presence of an ambient magnetic field reduces the erosion and fragmentation of the dense equatorial ring after the impact of the SN blast wave. As a result, the ring survives the passage of the blast, at least, during the time covered by the simulations (40 years). Our model is able to reproduce the morphology and lightcurves of SN 1987A in both X-ray and radio bands. The model reproduces the observed radio emission if the flux originating from the reverse shock is heavily suppressed. In this case, the radio emission originates mostly from the forward shock traveling through the H II region. Possible mechanisms for the suppression of emission from the reverse shock are investigated. We find that synchrotron self-absorption and free-free absorption have negligible effects on the emission during the interaction with the nebula. We suggest that the emission from the reverse shock at radio frequencies might be limited by highly magnetized ejecta.

We investigate the spectral properties of plasma turbulence from fluid to sub-ion scales by means of high-resolution three-dimensional (3D) numerical simulations performed with the hybrid particle-in-cell (HPIC) code CAMELIA. We produce extended turbulent spectra with well-defined power laws for the magnetic, ion bulk velocity, density, and electric fluctuations. The present results are in good agreement with previous two-dimensional (2D) HPIC simulations, especially in the kinetic range of scales, and reproduce several features observed in solar wind spectra. By providing scaling tests on many different architectures and convergence studies, we prove CAMELIA to represent a very efficient, accurate and reliable tool for investigating the develpoment of the turbulent cascade in the solar wind, being able to cover simultaneously several decades in wavenumber, also in 3D.

We perform MHD modeling of a single bright coronal loop to include the interaction with a non-uniform magnetic field. The field is stressed by random footpoint rotation in the central region and its energy is dissipated into heating by growing currents through anomalous magnetic diffusivity that switches on in the corona above a current density threshold. We model an entire single magnetic flux tube, in the solar atmosphere extending from the high-beta chromosphere to the low-beta corona through the steep transition region. The magnetic field expands from the chromosphere to the corona. The maximum resolution is ~30 km. We obtain an overall evolution typical of loop models and realistic loop emission in the EUV and X-ray bands. The plasma confined in the flux tube is heated to active region temperatures (~3 MK) after ~2/3 hr. Upflows from the chromosphere up to ~100 km/s fill the core of the flux tube to densities above 10^9 cm^-3. More heating is released in the low corona than the high corona and is finely structured both in space and time.

The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission measure. The cross sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the heating rate. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation in terms of steady heating models. We study the area response with a time-dependent 2D MHD loop model, including the description of the expanding magnetic field, coronal heating and losses by thermal conduction and radiation from optically thin plasma. We run a simulation for a loop 50 Mm long and quasi-statically heated to about 4 MK. We find that the area can change substantially with the quasi-steady heating rate, e.g. by ~40% at 0.5 MK as the loop temperature varies between 1 and 4 MK, and, therefore, affects the interpretation of DEM(T) curves.

One scenario proposed to explain the million degrees solar corona is a finely-stranded corona where each strand is heated by a rapid pulse. However, such fine structure has neither been resolved through direct imaging observations nor conclusively shown through indirect observations of extended superhot plasma. Recently it has been shown that the observed difference in appearance of cool and warm coronal loops (~1 MK, ~2-3 MK, respectively) -- warm loops appearing "fuzzier" than cool loops -- can be explained by models of loops composed of subarcsecond strands, which are impulsively heated up to ~10 MK. That work predicts that images of hot coronal loops (>~6 MK) should again show fine structure. Here we show that the predicted effect is indeed widely observed in an active region with the Solar Dynamics Observatory, thus supporting a scenario where impulsive heating of fine loop strands plays an important role in powering the active corona.

Coronal active regions are observed to get fuzzier and fuzzier (i.e. more and more confused and uniform) in harder and harder energy bands or lines. We explain this evidence as due to the fine multi-temperature structure of coronal loops. To this end, we model bundles of loops made of thin strands, each heated by short and intense heat pulses. For simplicity, we assume that the heat pulses are all equal and triggered only once in each strand at a random time. The pulse intensity and cadence are selected so as to have steady active region loops ($\sim 3$ MK), on the average. We compute the evolution of the confined heated plasma with a hydrodynamic loop model. We then compute the emission along each strand in several spectral lines, from cool ($\leq 1$ MK), to warm ($2-3$ MK) lines, detectable with Hinode/EIS, to hot X-ray lines. The strands are then put side-by-side to construct an active region loop bundle. We find that in the warm lines ($2-3$ MK) the loop emission fills all the available image surface. Therefore the emission appears quite uniform and it is difficult to resolve the single loops, while in the cool lines the loops are considerably more contrasted and the region is less fuzzy. The main reasons for this effect are that, during their evolution, i.e. pulse heating and slow cooling, each strand spends a relatively long time at temperatures around $2-3$ MK, and that it has a high emission measure during that phase, so the whole region appears more uniform or smudged. We make the prediction that the fuzziness should be reduced in the hot UV and X-ray lines.