Boundary layers in space and astrophysical plasmas are the location of complex dynamics where different mechanisms coexist and compete eventually leading to plasma mixing. In this work, we present fully kinetic Particle-In-Cell simulations of different boundary layers characterized by the following main ingredients: a velocity shear, a density gradient and a magnetic gradient localized at the same position. In particular, the presence of a density gradient drives the development of the lower hybrid drift instability (LHDI), which competes with the Kelvin-Helmholtz instability (KHI) in the development of the boundary layer. Depending on the density gradient, the LHDI can even dominate the dynamics of the layer. Because these two instabilities grow on different spatial and temporal scales, when the LHDI develops faster than the KHI an inverse cascade is generated, at least in 2D. This inverse cascade, starting at the LHDI kinetic scales, generates structures at scale lengths at which the KHI would typically develop. When that is the case, those structures can suppress the KHI itself because they significantly affect the underlying velocity shear gradient. We conclude that depending on the density gradient, the velocity jump and the width of the boundary layer, the LHDI in its nonlinear phase can become the primary instability for plasma mixing. These numerical simulations show that the LHDI is likely to be a dominant process at the magnetopause of Mercury. These results are expected to be of direct impact to the interpretation of the forthcoming BepiColombo observations.

A comet ionospheric model assuming the plasma to move radially outward with the same bulk speed as the neutral gas and not being subject to severe reduction through dissociative recombination has previously been tested in a series of case studies associated with the Rosetta mission at comet 67P/Churyumov-Gerasimenko. It has been found that at low activity and within several tens of km from the nucleus such models (which originally were developed for such conditions) generally work well in reproducing observed electron number densities, in particular when plasma production through both photoionization and electron-impact ionization is taken into account. Near perihelion, case studies have, on the contrary, showed that applying similar assumptions overestimates the observed electron number densities at the location of Rosetta. Here we compare ROSINA/COPS driven model results with RPC/MIP derived electron number densities for an extended time period (2015 November through 2016 March) during the post-perihelion phase with southern summer/spring. We observe a gradual transition from a state when the model grossly overestimates (by more than a factor of 10) the observations to being in reasonable agreement during 2016 March.

A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics.

Comets hold the key to the understanding of our solar system, its formation and its evolution, and to the fundamental plasma processes at work both in it and beyond it. A comet nucleus emits gas as it is heated by the sunlight. The gas forms the coma, where it is ionised, becomes a plasma and eventually interacts with the solar wind. Besides these neutral and ionised gases, the coma also contains dust grains, released from the comet nucleus. As a cometary atmosphere develops when the comet travels through the solar system, large-scale structures, such as the plasma boundaries, develop and disappear, while at planets such large-scale structures are only accessible in their fully grown, quasi-steady state. In situ measurements at comets enable us to learn both how such large-scale structures are formed or reformed and how small-scale processes in the plasma affect the formation and properties of these large scale structures. Furthermore, a comet goes through a wide range of parameter regimes during its life cycle, where either collisional processes, involving neutrals and charged particles, or collisionless processes are at play, and might even compete in complicated transitional regimes. Thus a comet presents a unique opportunity to study this parameter space, from an asteroid-like to a Mars- and Venus-like interaction. Fast flybys of comets have made many new discoveries, setting the stage for a multi-spacecraft mission to accompany a comet on its journey through the solar system. This white paper reviews the present-day knowledge of cometary plasmas, discusses the many questions that remain unanswered, and outlines a multi-spacecraft ESA mission to accompany a comet that will answer these questions by combining both multi-spacecraft observations and a rendezvous mission, and at the same time advance our understanding of fundamental plasma physics and its role in planetary systems.

This white paper proposes that AMBITION, a Comet Nucleus Sample Return mission, be a cornerstone of ESA's Voyage 2050 programme. We summarise some of the most important questions still open in cometary science after the successes of the Rosetta mission, many of which require sample analysis using techniques that are only possible in laboratories on Earth. We then summarise measurements, instrumentation and mission scenarios that can address these questions, with a recommendation that ESA select an ambitious cryogenic sample return mission. Rendezvous missions to Main Belt comets and Centaurs are compelling cases for M-class missions, expanding our knowledge by exploring new classes of comets. AMBITION would engage a wide community, drawing expertise from a vast range of disciplines within planetary science and astrophysics. With AMBITION, Europe will continue its leadership in the exploration of the most primitive Solar System bodies.

A major point of interest in cometary plasma physics has been the diamagnetic cavity, an unmagnetized region in the inner-most part of the coma. Here, we combine Langmuir and Mutual Impedance Probe measurements to investigate ion velocities and electron temperatures in the diamagnetic cavity of comet 67P, probed by the Rosetta spacecraft. We find ion velocities generally in the range 2-4 km/s, significantly above the expected neutral velocity $\lesssim$1~km/s, showing that the ions are (partially) decoupled from the neutrals, indicating that ion-neutral drag was not responsible for balancing the outside magnetic pressure. Observations of clear wake effects on one of the Langmuir probes showed that the ion flow was close to radial and supersonic, at least w.r.t. the perpendicular temperature, inside the cavity and possibly in the surrounding region as well. We observed spacecraft potentials $\lesssim$-5~V throughout the cavity, showing that a population of warm ($\sim$5~eV) electrons was present throughout the parts of the cavity reached by Rosetta. Also, a population of cold ($\lesssim0.1$~eV) electrons was consistently observed throughout the cavity, but less consistently in the surrounding region, suggesting that while Rosetta never entered a region of collisionally coupled electrons, such a region was possibly not far away during the cavity crossings.

The electron temperature of the plasma is one important aspect of the environment. Electrons created by photoionization or impact ionization of atmospheric gas have energies $\sim$10 eV. In an active comet coma, the gas density is high enough for rapid cooling of the electron gas to the neutral gas temperature (a few hundred kelvin). How cooling evolves in less active comets has not been studied before. Aims. We aim to investigate how electron cooling varied as comet 67P/Churyumov-Gerasimenko changed its activity by three orders of magnitude during the Rosetta mission. We used in situ data from the Rosetta plasma and neutral gas sensors. By combining Langmuir probe bias voltage sweeps and mutual impedance probe measurements, we determined at which time cold electrons formed at least 25\% of the total electron density. We compared the results to what is expected from simple models of electron cooling, using the observed neutral gas density as input. We demonstrate that the slope of the Langmuir probe sweep can be used as a proxy for the presence of cold electrons. We show statistics of cold electron observations over the two-year mission period. We find cold electrons at lower activity than expected by a simple model based on free radial expansion and continuous loss of electron energy. Cold electrons are seen mainly when the gas density indicates that an exobase may have formed. Collisional cooling of electrons following a radial outward path is not sufficient to explain the observations. We suggest that the ambipolar electric field keeps electrons in the inner coma for a much longer time, giving them time to dissipate energy by collisions with the neutrals. We conclude that better models are required to describe the plasma environment of comets. They need to include at least two populations of electrons and the ambipolar field.

We aim at analytically modelling the solar wind proton trajectories during their interaction with a partially ionised cometary atmosphere, not in terms of bulk properties of the flow but in terms of single particle dynamics. We first derive a generalised gyromotion, in which the electric field is reduced to its motional component. Steady-state is assumed, and simplified models of the cometary density and of the electron fluid are used to express the force experienced by individual solar wind protons during the interaction. A three-dimensional (3D) analytical expression of the gyration of two interacting plasma beams is obtained. Applying it to a comet case, the force on protons is always perpendicular to their velocity and has an amplitude proportional to 1/r^2. The solar wind deflection is obtained at any point in space. The resulting picture presents a caustic of intersecting trajectories, and a circular region is found that is completely free of particles. The particles do not lose any kinetic energy and this absence of deceleration, together with the solar wind deflection pattern and the presence of a solar wind ion cavity, is in good agreement with the general results of the Rosetta mission. The qualitative match between the model and the in situ data highlights how dominant the motional electric field is throughout most of the interaction region for the solar wind proton dynamics. The model provides a simple general kinetic description of how momentum is transferred between these two collisionless plasmas. It also shows the potential of this semi-analytical model for a systematic quantitative comparison to the data.

The first 1000 km of the ion tail of comet 67P/Churyumov-Gerasimenko were explored by the European Rosetta spacecraft, 2.7 au away from the Sun. We characterised the dynamics of both the solar wind and the cometary ions on the night-side of the comet's atmosphere. We analysed in situ ion and magnetic field measurements and compared the data to a semi-analytical model. The cometary ions are observed flowing close to radially away from the nucleus during the entire excursion. The solar wind is deflected by its interaction with the new-born cometary ions. Two concentric regions appear, an inner region dominated by the expanding cometary ions and an outer region dominated by the solar wind particles. The single night-side excursion operated by Rosetta revealed that the near radial flow of the cometary ions can be explained by the combined action of three different electric field components, resulting from the ion motion, the electron pressure gradients, and the magnetic field draping. The observed solar wind deflection is governed mostly by the motional electric field.

We present Rosetta RPC case study from four events at various radial distance, phase angle and local time from autumn 2015, just after perihelion of comet 67P/Churyumov-Gerasimenko. Pulse like (high amplitude, up to minutes in time) signatures are seen with several RPC instruments in the plasma density (LAP, MIP), ion energy and flux (ICA) as well as magnetic field intensity (MAG). Furthermore the cometocentric distance relative to the electron exobase is seen to be a good organizing parameter for the measured plasma variations. The closer Rosetta is to this boundary, the more pulses are measured. This is consistent with the pulses being filaments of plasma originating from the diamagnetic cavity boundary as predicted by simulations.

Strong electron cooling on the neutral gas in cometary comae has been predicted for a long time, but actual measurements of low electron temperature are scarce. We present in situ measurements of plasma density, electron temperature and spacecraft potential by the Rosetta Langmuir probe instrument, LAP. Data acquired within a few hundred km from the nucleus are dominated by a warm component with electron temperature typically 5--10 eV at all heliocentric distances covered (1.25 to 3.83 AU). A cold component, with temperature no higher than about 0.1 eV, appears in the data as short (few to few tens of seconds) pulses of high probe current, indicating local enhancement of plasma density as well as a decrease in electron temperature. These pulses first appeared around 3 AU and were seen for longer periods close to perihelion. The general pattern of pulse appearance follows that of neutral gas and plasma density. We have not identified any periods with only cold electrons present. The electron flux to Rosetta was always dominated by higher energies, driving the spacecraft potential to order -10 V. The warm (5--10 eV) electron population is interpreted as electrons retaining the energy they obtained when released in the ionisation process. The sometimes observed cold populations with electron temperatures below 0.1 eV verify collisional cooling in the coma. The cold electrons were only observed together with the warm population. The general appearance of the cold population appears to be consistent with a Haser-like model, implicitly supporting also the coupling of ions to the neutral gas. The expanding cold plasma is unstable, forming filaments that we observe as pulses.

We use measurements from the Rosetta plasma consortium (RPC) Langmuir probe (LAP) and mutual impedance probe (MIP) to study the spatial distribution of low-energy plasma in the near-nucleus coma of comet 67P/Churyumov-Gerasimenko. The spatial distribution is highly structured with the highest density in the summer hemisphere and above the region connecting the two main lobes of the comet, i.e. the neck region. There is a clear correlation with the neutral density and the plasma to neutral density ratio is found to be about 1-2x10^-6, at a cometocentric distance of 10 km and at 3.1 AU from the sun. A clear 6.2 h modulation of the plasma is seen as the neck is exposed twice per rotation. The electron density of the collisonless plasma within 260 km from the nucleus falls of with radial distance as about 1/r. The spatial structure indicates that local ionization of neutral gas is the dominant source of low-energy plasma around the comet.

Using a full kinetic, implicit particle-in-cell code, iPiC3D, we studied the properties of plasma kinetic turbulence, such as would be found at the interface between the solar wind and the Earth magnetosphere at low latitude during northwards periods. In this case, in the presence of a magnetic field B oriented mostly perpendicular to the velocity shear, turbulence is fed by the disruption of a Kelvin-Helmholtz vortex chain via secondary instabilities, vortex pairing and non-linear interactions. We found that the magnetic energy spectral cascade between ion and electron inertial scales, $d_i$ and $d_e$, is in agreement with satellite observations and other previous numerical simulations; however, in our case the spectrum ends with a peak beyond $d_e$ due to the occurrence of the lower hybrid drift instability. The electric energy spectrum is influenced by effects of secondary instabilities: anomalous resistivity, fed by the development of the lower hybrid drift instability, steepens the spectral decay and, depending on the alignment or anti-alignment of B and the shear vorticity, peaks due to ion-Bernstein waves may dominate the spectrum around $d_i$. These waves are generated by counter-streaming flow structures, through flux freezing also responsible for reconnection of the in-plane component of the magnetic field, which then generates electron pressure anisotropy and flattening of the field-aligned component of the electron distribution function.

The kinetic features of secondary magnetic reconnection in a single flux rope undergoing internal kink instability are studied by means of three-dimensional Particle-in-Cell simulations. Several signatures of secondary magnetic reconnection are identified in the plane perpendicular to the flux rope: a quadrupolar electron and ion density structure and a bipolar Hall magnetic field develop in proximity of the reconnection region. The most intense electric fields form perpendicularly to the local magnetic field, and a reconnection electric field is identified in the plane perpendicular to the flux rope. An electron current develops along the reconnection line in the opposite direction of the electron current supporting the flux rope magnetic field structure. Along the reconnection line, several bipolar structures of the electric field parallel to the magnetic field occur making the magnetic reconnection region turbulent. The reported signatures of secondary magnetic reconnection can help to localize magnetic reconnection events in space, astrophysical and fusion plasmas.

The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic, Hall-MHD, two-fluid, hybrid kinetic and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligable in the nonlinear regime. This study show that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the colisionless regime.

The dynamics of a plasmoid chain is studied with three dimensional Particle-in-Cell simulations. The evolution of the system with and without a uniform guide field, whose strength is 1/3 the asymptotic magnetic field, is investigated. The plasmoid chain forms by spontaneous magnetic reconnection: the tearing instability rapidly disrupts the initial current sheet generating several small-scale plasmoids, that rapidly grow in size coalescing and kinking. The plasmoid kink is mainly driven by the coalescence process. It is found that the presence of guide field strongly influences the evolution of the plasmoid chain. Without a guide field, a main reconnection site dominates and smaller reconnection regions are included in larger ones, leading to an hierarchical structure of the plasmoid-dominated current sheet. On the contrary in presence of a guide field, plasmoids have approximately the same size and the hierarchical structure does not emerge, a strong core magnetic field develops in the center of the plasmoid in the direction of the existing guide field, and bump-on-tail instability, leading to the formation of electron holes, is detected in proximity of the plasmoids.

Recent observation of large amplitude Langmuir waveforms during a Type III event in the solar wind have been interpreted as the signature of the electrostatic decay of beam-driven Langmuir waves. This mechanism is thought to be a first step to explain the generation of solar Type III radio emission. The threshold for this parametric instability in typical solar wind condition is investigated here by means of 1D-1V Vlasov-Poisson simulations. We show that the amplitude of the observed Langmuir beat-like waveforms is of the order of the effective threshold computed from the full kinetic simulations. The expected level of associated ion acoustic density fluctuations have also been computed for comparison with observations.

The electrostatic decay enables energy transfer from a finite amplitude Langmuir to a backscattered daughter Langmuir wave and ion acoustic density fluctuations. This mechanism is thought to be a first step for the generation of type III solar radio emissions at twice the plasma frequency. The electrostatic decay is here investigated through Vlasov-Poisson simulations by considering Langmuir localized wave packets in the case Te = Tp. Simulation results are found to be in good agreement with recently reported observations from the STEREO mission of the electrostatic decay of beam-driven Langmuir waves during a type III burst.

We investigate the behavior of the magnetic Prandtl number (ratio of microscopic viscosity to resistivity) for accretion sources. Generally this number is very small in standard accretion disk models, but can become larger than unity within $\sim 50$ Schwarzschild radii of the central mass. Recent numerical investigations suggest a marked dependence of the level of MHD turbulence on the value of the Prandtl number. Hence, black hole and neutron star accretors, i.e. compact X-ray sources, are affected. The astrophysical consequences of this could be significant, including a possible route to understanding the mysterious state changes that have long characterized these sources.