Energy conversion via reconnecting current sheets is common in space and astrophysical plasmas. Frequently, current sheets disrupt at multiple reconnection sites, leading to the formation of plasmoid structures between sites, which might affect energy conversion. We present in situ evidence of the firehose instability in multiple reconnection in the Earth's magnetotail. The observed proton beams accelerated in the direction parallel to magnetic field and ion-scale fluctuations of whistler type imply the development of firehose instability between two active reconnection sites. The linear wave dispersion relation, estimated for the measured plasma parameters, indicates a positive growth rate of firehose-related electromagnetic fluctuations. Simulations of temporal evolution of the observed multiple reconnection by using a 2.5D implicit particle-in-cell code show that, as the plasmoid formed between two reconnection sites evolves, the plasma at its edge becomes anisotropic and overcomes the firehose marginal stability threshold, leading to the generation of magnetic field fluctuations. The combined results of observations and simulations suggest that the firehose instability, operating between reconnection sites, converts plasma kinetic energy into energy of magnetic field fluctuations, counteracting the conversion of magnetic energy into plasma energy occurring at reconnection sites. This suggests that magnetic energy conversion in multiple reconnection can be less efficient than in the case of the single-site reconnection.

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.

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.

We have revisited the sdB+K-star long-period binary PG 1018-047 based on 20 new high-resolution VLT/UVES spectra that provided regular coverage over a period of more than 26 months. We refine the period and establish that the orbit is significantly eccentric ($P=751.6\pm1.9\rm{d}$ and $e=0.049\pm0.008$). A simultaneous fit derived from the narrow metal lines visible in the spectrum of the sdB star and the metal lines in the red part of the spectrum that originate from the companion provides the mass ratio, $M_{\rm MS}/M_{\rm sdB} = 1.52 \pm 0.04$, for the system. From an NLTE model atmosphere analysis of the combined spectra, we find $T_{\rm eff}$=29900$\pm$330K, log$g$= 5.65$\pm$0.06dex and log($n_{\rm He}$/$n_{\rm H}$)=-3.98$\pm$0.16dex for the primary, consistent with a B-type hot subdwarf star. The spectral contribution of the companion is consistent with a K5V-type star. With the companion having a mass of only $\sim$0.7$\rm M_{\odot}$ this system lies close to the boundary below which stable Roche-lobe overflow (RLOF) cannot be supported. To model the evolution of such a system, we have extended earlier MESA models towards lower companion masses. We find that both phase-dependent mass loss during RLOF, when 30 to 40\% of the available mass is lost through the outer Lagrange point, and phase-dependent mass loss during RLOF in combination with a circumbinary disk of maximum $M_{\rm CB} = 0.001\rm M_{\odot}$ could have formed the PG1018-047 binary system.

We present a systematic attempt to study magnetic null points and the associated magnetic energy conversion in kinetic Particle-in-Cell simulations of various plasma configurations. We address three-dimensional simulations performed with the semi-implicit kinetic electromagnetic code iPic3D in different setups: variations of a Harris current sheet, dipolar and quadrupolar magnetospheres interacting with the solar wind; and a relaxing turbulent configuration with multiple null points. Spiral nulls are more likely created in space plasmas: in all our simulations except lunar magnetic anomaly and quadrupolar mini-magnetosphere the number of spiral nulls prevails over the number of radial nulls by a factor of 3-9. We show that often magnetic nulls do not indicate the regions of intensive energy dissipation. Energy dissipation events caused by topological bifurcations at radial nulls are rather rare and short-lived. The so-called X-lines formed by the radial nulls in the Harris current sheet and lunar magnetic anomaly simulations are rather stable and do not exhibit any energy dissipation. Energy dissipation is more powerful in the vicinity of spiral nulls enclosed by magnetic flux ropes with strong currents at their axes (their cross-sections resemble 2D magnetic islands). These null lines reminiscent of Z-pinches efficiently dissipate magnetic energy due to secondary instabilities such as the two-stream or kinking instability, accompanied by changes in magnetic topology. Current enhancements accompanied by spiral nulls may signal magnetic energy conversion sites in the observational data.

Long-orbital-period subdwarf B (sdB) stars with main-sequence companions are believed to be the product of stable Roche Lobe overflow (RLOF), a scenario challenged by recent observations. Here we represent the results of a systematic study of the orbital-period distribution of sdB binaries in this channel using detailed binary evolution calculations. We show that the observed orbital-period distribution of long-period sdB binaries can be well explained by this scenario. Furthermore, we find that, if the progenitors of the sdB stars have initial masses below the helium flash mass, the sdB binaries produced from stable RLOF follow a unique mass -- orbital period relation for a given metallicity $Z$; increasing the orbital period from $\sim 400$ to $\sim 1100$\,d corresponds to increasing the mass of the sdB star from $\sim 0.40$ to $\sim 0.49\,M_\odot$ for $Z=0.02$. We suggest that the longest sdB binaries (with orbital period $> 1100$\,d) could be the result of atmospheric RLOF. The mass -- orbital period relation can be tested observationally if the mass of the sdB star can be determined precisely, e.g.\ from asteroseismology. Using this relation, we revise the orbital period distribution of sdB binaries produced by the first stable RLOF channel for the best fitting model of Han et al (2003), and show that the orbital period has a peak around 830\,d.

About 50% of all known hot subdwarf B stars (sdBs) reside in close (short period) binaries, for which common envelope ejection is the most likely formation mechanism. However, Han et al. 2003 predict that the majority of sdBs should form through stable mass transfer leading to long period binaries. Determining orbital periods for these systems is challenging and while the orbital periods of ~100 short period systems have been measured, there are no periods measured above 30 days. As part of a large program to characterise the orbital periods of subdwarf B binaries and their formation history, we have found that PG1018-047 has an orbital period of 760\pm6 days, easily making it the longest period ever detected for a subdwarf B binary. Exploiting the Balmer lines of the subdwarf primary and the narrow absorption lines of the companion present in the spectra, we derive the radial velocity amplitudes of both stars, and estimate the mass ratio M_{MS}/M_{sdB} = 1.6\pm0.2. From the combination of visual and infrared photometry, the spectral type of the companion star is determined to be mid K.