Collisionless magnetic reconnection is a prime candidate to account for flare-like or steady emission, outflow launching, or plasma heating, in a variety of high-energy astrophysical objects, including ones with relativistic ion-electron plasmas. But the fate of the initial magnetic energy in a reconnection event remains poorly known. What are the amounts assigned to kinetic energy, the ion and electron distribution, and the hardness of the particle distributions? We explored these questions with 2D particle-in-cell simulations of ion-electron plasmas. We find that 45 to 75% of the total initial magnetic energy ends up in kinetic energy, this fraction increasing with the inflow magnetization. Depending on the guide field strength, ions get from 30% to 60% of the total kinetic energy. Particles can be separated into two populations that mix only weakly: (i) particles initially in the current sheet heated by its initial tearing and subsequent contraction of the islands, and (ii) particles from the background plasma that primarily gain energy via the reconnection electric field when passing near the X-point. Particles of (ii) tend to form a power law with an index p = −dlog n(γ)/dlog γ that depends mostly on the inflow Alfvén speed V A and magnetization σ s of species s. For electrons p = 5 to 1.2 for increasing σ e . The highest particle Lorentz factor for ions or electrons increases roughly linearly with time for all the relativistic simulations. This is faster, and the spectra can be harder, than for collisionless shock acceleration. We discuss applications to microquasar and AGN coronae, to extragalactic jets, and to radio lobes. We point out situations where effects, such as Compton drag or pair creation, are important.

Magnetic reconnection is a leading mechanism for magnetic energy conversion and high-energy non-thermal particle production in a variety of high-energy astrophysical objects, including ones with relativistic ion-electron plasmas (e.g., microquasars or AGNs), a regime where first principle studies are scarce. We present 2D particle-in-cell (PIC) simulations of low β ion-electron plasmas under relativistic conditions, i.e., with inflow magnetic energy exceeding the plasma restmass energy. We identify outstanding properties: (i) For relativistic inflow magnetizations (here 10 ≤ σ e ≤ 360), the reconnection outflows are dominated by thermal agitation instead of bulk kinetic energy. (ii) At high inflow electron magnetization (σ e ≥ 80), the reconnection electric field is sustained more by bulk inertia than by thermal inertia. It challenges the thermal-inertia paradigm and its implications. (iii) The inflows feature sharp transitions at the entrance of the diffusion zones. These are not shocks but results from particle ballistic motions, all bouncing at the same location, provided that the thermal velocity in the inflow is far lower than the inflow E × B bulk velocity. (iv) Island centers are magnetically isolated from the rest of the flow and can present a density depletion at their center. (v) The reconnection rates are slightly higher than in non-relativistic studies. They are best normalized by the inflow relativistic Alfvén speed projected in the outflow direction, which then leads to rates in a close range (0.14-0.25), thus allowing for an easy estimation of the reconnection electric field.

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 nonlinear 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. We link these discrepancies to the superparticle number per cell and to the level of field fluctuations. Second, we detail a new method for initial loading of Maxwell-Jüttner particle distributions with relativistic bulk velocity and relativistic temperature, and explain why the traditional method with individual particle boosting fails. The formulation of the relativistic Harris equilibrium is generalized to arbitrary temperature and mass ratios. Both are required for the tearing instability setup. Third, we turn to the key point of this paper and scrutinize the question of what description of (weakly coupled) 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 Λ PIC , which we show to behave as Λ PIC ∝ 1/p. We explore two important 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, a fact that we confirm with simulations. The other is that the level of electric field fluctuations scales as 1/Λ PIC ∝ p. We provide a corresponding exact expression, taking into account the finite superparticle size. We confirm both expectations with simulations. Fourth, we compare the Vlasov-Maxwell theory, often used for code benchmarking, to the PIC model. The former describes a phase-space fluid with Λ = +∞ and no correlations, while the PIC plasma features a small Λ and a high level of correlations when compared to a real plasma. These differences have to be kept in mind when interpreting and validating PIC results against the Vlasov-Maxwell theory and when modeling real physical plasmas.

Using a new parallel algorithm implemented within the VisIt framework, we analysed large cosmological grid simulations to study the properties of baryons in filaments. The procedure allows us to build large catalogues with up to ∼ 3 · 10 4 filaments per simulated volume and to investigate the properties of cosmic filaments for very large volumes at high resolution (up to 300 3 Mpc 3 simulated with 2048 3 cells). We determined scaling relations for the mass, volume, length and temperature of filaments and compared them to those of galaxy clusters. The longest filaments have a total length of about 200 Mpc with a mass of several 10 15 M . We also investigated the effects of different gas physics. Radiative cooling significantly modifies the thermal properties of the warm-hot-intergalactic medium of filaments, mainly by lowering their mean temperature via line cooling. On the other hand, powerful feedback from active galactic nuclei in surrounding halos can heat up the gas in filaments. The impact of shock-accelerated cosmic rays from diffusive shock acceleration on filaments is small and the ratio of between cosmic ray and gas pressure within filaments is of the order of ∼ 10 − 20 percent.

Aims. Observations show nebulae around some massive stars but not around others. If observed, their chemical composition is far from homogeneous. Our goal is to put these observational features into the context of the evolution of massive stars and their circumstellar medium (CSM) and, more generally, to quantify the role of massive stars for the chemical and dynamical evolution of the ISM. Methods. Using the A-MAZE code, we perform 2d-axisymmetric hydrodynamical simulations of the evolution of the CSM, shaped by stellar winds, for a whole grid of massive stellar models from 15 to 120 M and following the stellar evolution from the zero-age main-sequence to the time of supernova explosion. In addition to the usual quantities, we also follow five chemical species: H, He, C, N, and O. Results. We show how various quantities evolve as a function of time: size of the bubble, position of the wind termination shock, chemical composition of the bubble, etc. The chemical composition of the bubble changes considerably compared to the initial composition, particularly during the red-supergiant (RSG) and Wolf-Rayet (WR) phases. In some extreme cases, the inner region of the bubble can be completely depleted in hydrogen and nitrogen, and is mainly composed of carbon, helium, and oxygen. We argue why the bubble typically expands at a lower rate than predicted by self-similarity theory. In particular, the size of the bubble is very sensitive to the density of the ISM, decreasing by a factor of ∼2.5 for each additional dex in ISM density. The bubble size also decreases with the metallicity of the central star, because low-metallicity stars have weaker winds. Our models qualitatively fit the observations of WR ejecta nebulae.