Search

We preform hybrid simulations of super Alfvenic quasi-parallel shock, driven by a Coronal Mass Ejection (CME), propagating in the Outer Coronal or Solar Wind at distances of between 3 to 6 solar radii. The hybrid treatment of the problem enable the study of the shock propagation on the ion time scale, preserving ion kinetics and allowing for a self consistent treatment of the shock propagation and particle acceleration. The CME plasma drags the embedded magnetic field lines stretching from the sun, and propagates out into interplanetary space at a greater velocity than the in-situ solar wind, driving the shock, and producing very energetic particles. Our results show electromagnetic Alfven waves are generated at the shock front. The waves propagate upstream of the shock and are produced by the counter streaming ions of the solar wind plasma being reflected at the shock. A significant fraction of the particles are accelerated in two distinct phases: first, particles drift from the shock and are accelerated in the upstream region and, second, particles arriving at the shock get trapped, and are accelerated at the shock front. A fraction of the particles diffused back to the shock, which is consistent with the Fermi acceleration mechanism., Comment: Accepted for publication in Astrophysical Journal 18th June 2014

In this paper we present in-situ satellite data, theory and laboratory validation that show how small scale collisionless shocks and mini-magnetospheres can form on the electron inertial scale length. The resulting retardation and deflection of the solar wind ions could be responsible for the unusual "lunar swirl" patterns seen on the surface of the Moon., Comment: 5 pages, 5 figures

We explore the physics of shock evolution and particle acceleration in non-relativistic collisionless shocks using multidimensional hybrid simulations. We analyze a wide range of physical parameters relevant to the acceleration of cosmic rays (CRs) in astrophysical non-relativistic shock scenarios, such as in supernova remnant (SNR) shocks. We explore the evolution of the shock structure and particle acceleration efficiency as a function of Alfv\'enic Mach number and magnetic field inclination angle $\theta$. We show that there are fundamental differences between high and low Mach number shocks in terms of the electromagnetic turbulence generated in the pre-shock zone and downstream; dominant modes are resonant with the streaming CRs in the low Mach number regime, while both resonant and non-resonant modes are present for high Mach numbers. Energetic power law tails for ions in the downstream plasma can account for up to 15% of the incoming upstream flow energy, distributed over $\sim5%$ of the particles in a power law with slope $-2\pm0.2$ in energy. The energy conversion efficiency (for CRs) peaks at $\theta=15^\circ$ to $30^\circ$ and $M_A=6$, and decreases for higher Mach numbers, down to $\sim2%$ for $M_A=31$. Accelerated particles are produced by Diffusive Shock Acceleration (DSA) and by Shock Drift Acceleration (SDA) mechanisms, with the SDA contribution to the overall energy gain increasing with magnetic inclination. We also present a direct comparison between hybrid and fully kinetic particle-in-cell results at early times; the agreement between the two models justifies the use of hybrid simulations for longer-term shock evolution. In SNR shocks, particle acceleration will be significant for low Mach number quasi-parallel flows ($M_A < 30$, $\theta< 45$). This finding underscores the need for effective magnetic amplification mechanism in SNR shocks.

A non-resonant instability for the amplification of the interstellar magnetic field in young Supernova Remnant (SNR) shocks was predicted by Bell (2004), and is thought to be relevant for the acceleration of cosmic ray (CR) particles. For this instability, the CRs streaming ahead of SNR shock fronts drive electromagnetic waves with wavelengths much shorter than the typical CR Larmor radius, by inducing a current parallel to the background magnetic field. We explore the nonlinear regime of the non-resonant mode using Particle-in-Cell (PIC) hybrid simulations, with kinetic ions and fluid electrons, and analyze the saturation mechanism for realistic CR and background plasma parameters. In the linear regime, the observed growth rates and wavelengths match the theoretical predictions; the nonlinear stage of the instability shows a strong reaction of both the background plasma and the CR particles, with the saturation level of the magnetic field varying with the CR parameters. The simulations with CR-to-background density ratios of n_CR/n_b=10^(-5) reveal the highest magnetic field saturation levels, with energy also being transferred to the background plasma and to the perpendicular velocity components of the CR particles. The results show that amplification factors >10 for the magnetic field can be achieved, and suggest that this instability is important for the generation of magnetic field turbulence, and for the acceleration of CR particles.

Three-dimensional (3D) hybrid particle-in-cell (PIC) simulations, with kinetic ions and fluid electrons, of a plasma cloud expansion in the solar wind are presented, revealing the dynamics of the expansion, with shock formation, magnetic field compression, and the solar wind ions deflection around the plasma bubble. The similarities of this system with a magnetosphere are also pointed out.

Solar energetic ions are a known hazard to both spacecraft electronics and to manned space flights in interplanetary space missions that extend over a long period of time. A dipole-like magnetic field and a plasma source, forming a mini magnetosphere, are being tested in the laboratory as means of protection against such hazards. We investigate, via particle-in-cell hybrid simulations, using kinetic ions and fluid electrons, the characteristics of the mini magnetospheres. Our results, for parameters identical to the experimental conditions, reveal the formation of a mini-magnetosphere, whose features are scanned with respect to the plasma density, the plasma flow velocity, and the intensity of the dipole field. Comparisons with a simplified theoretical model reveal a good qualitative agreement and excellent quantitative agreement for higher plasma dynamic pressures and lower B-fields.

A massively parallel simulation code, called \textit{dHybrid}, has been developed to perform global scale studies of space plasma interactions. This code is based on an explicit hybrid model; the numerical stability and parallel scalability of the code are studied. A stabilization method for the explicit algorithm, for regions of near zero density, is proposed. Three-dimensional hybrid simulations of the interaction of the solar wind with unmagnetized artificial objects are presented, with a focus on the expansion of a plasma cloud into the solar wind, which creates a diamagnetic cavity and drives the Interplanetary Magnetic Field out of the expansion region. The dynamics of this system can provide insights into other similar scenarios, such as the interaction of the solar wind with unmagnetized planets., Comment: 15 pages, 6 figures, accepted for publication in Computer Physics Communications