No Arabic abstract
High-mass microquasars consist of a massive star and a compact object, the latter producing jets that will interact with the stellar wind. The evolution of the jets, and ultimately their radiative outcome, could depend strongly on the inhomogeneity of the wind, which calls for a detailed study. The hydrodynamics of the interaction between a jet and a clumpy wind is studied, focusing on the global wind and single clump-jet interplay. We have performed, using the code textit{Ratpenat}, three-dimensional numerical simulations of a clumpy wind interacting with a mildly relativistic jet, and of individual clumps penetrating into a jet. For typical wind and jet velocities, filling factors of about > 0.1 are already enough for the wind to be considered as clumpy. An inhomogeneous wind makes the jet more unstable when crossing the system. Kinetic luminosities of the order 1.e37 erg/s allow the jet to reach the borders of a compact binary with an O star, as in the smooth wind case, although with a substantially higher degree of disruption. When able to enter into the jet, clumps are compressed and heated during a time of about their size divided by the sound speed in the shocked clump. Then, clumps quickly disrupt, mass-loading and slowing down the jet. We conclude that moderate wind clumpiness makes already a strong difference with the homogeneous wind case, enhancing jet disruption, mass-loading, bending, and likely energy dissipation in the form of emission. All this can have observational consequences at high-energies and also in the large scale radio jets.
In this paper we present steady-state RMHD simulations that include a mass-load term to study the process of jet deceleration. The mass-load mimics the injection of a proton-electron plasma from stellar winds within the host galaxy into initially pair plasma jets, with mean stellar mass-losses ranging from $10^{-14}$ to $10^{-9},{M_odot,yr^{-1}}$. The spatial jet evolution covers $sim 500,{rm pc}$ from jet injection in the grid at 10~pc from the jet nozzle. Our simulations use a relativistic gas equation of state and a pressure profile for the ambient medium. We compare these simulations with previous dynamical simulations of relativistic, non-magnetised jets. Our results show that toroidal magnetic fields can prevent fast jet expansion and the subsequent embedding of further stars via magnetic tension. In this sense, magnetic fields avoid a runaway deceleration process. Furthermore, when the mass-load is large enough to increase the jet density and produce fast, differential jet expansion, the conversion of magnetic energy flux into kinetic energy flux (i.e., magnetic acceleration), helps to delay the deceleration process with respect to non-magnetised jets. We conclude that the typical stellar population in elliptical galaxies cannot explain jet deceleration in classical FRI radio galaxies. However, we observe a significant change in the jet composition, thermodynamical parameters and energy dissipation along its evolution, even for moderate values of the mass-load.
The origin of Galactic cosmic rays remains a matter of debate, but supernova remnants are commonly considered to be the main place where high-energy cosmic rays are accelerated. Nevertheless, current models predict cosmic-ray spectra that do not match observations and the efficiency of the acceleration mechanism is still undetermined. On the other hand, the contribution of other kinds of sources to the Galactic cosmic-ray population is still unclear, and merits investigation. In this work we explore a novel mechanism through which microquasars might produce cosmic rays. In this scenario, microquasar jets generate relativistic neutrons, which escape and decay outside the system; protons and electrons, created when these neutrons decay, escape to the interstellar medium as cosmic rays. The most promising scenarios arise in extremely luminous systems ($L_mathrm{jet} sim 10^{40},mathrm{erg , s}^{-1}$), in which the fraction of jet power deposited in cosmic rays can reach $sim 0.001$. Slow jets ($Gamma lesssim 2$, where $Gamma$ is the bulk Lorentz factor) favour neutron production. The resulting cosmic-ray spectrum is similar for protons and electrons, which share the power in the ratio given by neutron decay. The spectrum peaks at roughly half the minimum energy of the relativistic protons in the jet; it is soft (spectral index $sim 3$) above this energy, and almost flat below. Values of spectral index steeper than $2$ are possible for cosmic rays in our model and these indeed agree with those required to explain the spectral signatures of Galactic cosmic rays, although only the most extreme microquasars provide power comparable to that of a typical supernova remnant. The mechanism explored in this work may provide stronger and softer cosmic-ray sources in the early Universe, and therefore contribute to the heating and reionisation of the intergalactic medium.
The properties of relativistic jets, their interaction with the ambient environment, and particle acceleration due to kinetic instabilities are studied self-consistently with Particle-in-Cell simulations. An important key issue is how a toroidal magnetic field affects the evolution of an e$^{pm}$ and an e$^{-}$ - p$^{+}$ jet, how kinetic instabilities such as the Weibel instability (WI), the mushroom instability (MI) and the kinetic Kelvin-Helmholtz instability (kKHI) are excited, and how such instabilities contribute to particle acceleration. We show that WI, MI and kKHI excited at the linear stage, generate a quasi-steady $x$-component of electric field which accelerates and decelerates electrons. In this work, we use a new jet injection scheme where an electric current is self-consistently generated at the jet orifice by the jet particles. We inject both e$^{pm}$ and e$^{-}$ - p$^{+}$ jets with a toroidal magnetic field (with a top-hat jet density profile) and for a sufficiently long time in order to examine the non-linear effects of the jet evolution. Despite the weakness of the initial magnetic field, we observe significant differences in the structure of the strong electromagnetic fields that are driven by the kinetic instabilities. We find that different jet compositions present different strongly excited instability modes. The magnetic field in the non-linear stage generated by different instabilities becomes dissipated and reorganized into a new topology. The 3-dimensional magnetic field topology indicates possible reconnection sites and the accelerated particles are significantly accelerated in the non-linear stage by the dissipation of the magnetic field and/or reconnection. This study will shed further light on the nature of astrophysical relativistic magnetized jet phenomena.
We investigate the interplay between jets from Active Galactic Nuclei (AGNs) and the surrounding InterStellar Medium (ISM) through full 3D, high resolution, Adaptive Mesh Refinement simulations performed with the FLASH code. We follow the jet- ISM system for several Myr in its transition from an early, compact source to an extended one including a large cocoon. During the jet evolution, we identify three major evolutionary stages and we find that, contrary to the prediction of popular theoretical models, none of the simulations shows a self-similar behavior. We also follow the evolution of the energy budget, and find that the fraction of input power deposited into the ISM (the AGN coupling constant) is of order of a few percent during the first few Myr. This is in broad agreement with galaxy formation models employing AGN feedback. However, we find that in these early stages, this energy is deposited only in a small fraction (< 1%) of the total ISM volume. Finally we demonstrate the relevance of backflows arising within the extended cocoon generated by a relativistic AGN jet within the ISM of its host galaxy, previously proposed as a mechanism for self-regulating the gas accretion onto the central object. These backflows tend later to be destabilized by the 3D dynamics, rather than by hydrodynamic (Kelvin- Helmholtz) instabilities. Yet, in the first few hundred thousand years, backflows may create a central accretion region of significant extent, and convey there as much as a few millions of solar masses.
We present Athena++ grid-based, hydrodynamic simulations of accretion onto Sagittarius A* via the stellar winds of the $sim 30$ Wolf-Rayet stars within the central parsec of the galactic center. These simulations span $sim$ 4 orders of magnitude in radius, reaching all the way down to 300 gravitational radii of the black hole, $sim 32$ times further in than in previous work. We reproduce reasonably well the diffuse thermal X-ray emission observed by Chandra in the central parsec. The resulting accretion flow at small radii is a superposition of two components: 1) a moderately unbound, sub-Keplerian, thick, pressure-supported disc that is at most (but not all) times aligned with the clockwise stellar disc, and 2) a bound, low-angular momentum inflow that proceeds primarily along the southern pole of the disc. We interpret this structure as a natural consequence of a few of the innermost stellar winds dominating accretion, which produces a flow with a broad distribution of angular momentum. Including the star S2 in the simulation has a negligible effect on the flow structure. Extrapolating our results from simulations with different inner radii, we find an accretion rate of $sim$ a few $times 10^{-8} M_odot$/yr at the horizon scale, consistent with constraints based on modeling the observed emission of Sgr A*. The flow structure found here can be used as more realistic initial conditions for horizon scale simulations of Sgr A*.