Merging massive black holes: the right place and the right time

The LIGO/Virgo detections of gravitational waves from merging black holes of $\simeq$ 30 solar mass suggest progenitor stars of low metallicity (Z/Z$_{\odot} \lesssim 0.3$). In this talk I will provide constrains on where the progenitors of GW150914 and GW170104 may have formed, based on advanced models of galaxy formation and evolution combined with binary population synthesis models. First I will combine estimates of galaxy properties (star-forming gas metallicity, star formation rate and merger rate) across cosmic time to predict the low redshift BBH merger rate as a function of present day host galaxy mass, formation redshift of the progenitor system and different progenitor metallicities. I will show that the signal is dominated by binaries formed at the peak of star formation in massive galaxies with and binaries formed recently in dwarf galaxies. Then, I will present what very high resolution hydrodynamic simulations of different galaxy types can learn us about their black hole populations.Comment: Proceeding of IAU Symposium 338 : "Gravitational Waves Astrophysics : Early results from GW searches and EM counterparts

Colliding wind binaries and gamma-ray binaries : relativistic version of the RAMSES code

Gamma-ray binaries are colliding wind binaries (CWB) composed of a massive star a non-accreting pulsar with a highly relativistic wind. Particle acceleration at the shocks results in emission going from extended radio emission to the gamma-ray band. The interaction region is expected to show common features with stellar CWB. Performing numerical simulations with the hydrodynamical code RAMSES, we focus on their structure and stability and find that the Kelvin-Helmholtz instability (KHI) can lead to important mixing between the winds and destroy the large scale spiral structure. To investigate the impact of the relativistic nature of the pulsar wind, we extend RAMSES to relativistic hydrodynamics (RHD). Preliminary simulations of the interaction between a pulsar wind and a stellar wind show important similarities with stellar colliding winds with small relativistic corrections.Comment: Proceeding of the 5th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma2012). arXiv admin note: text overlap with arXiv:1212.404

Merging massive black holes the right place and the right time

The LIGO/Virgo detections of gravitational waves from merging black holes of ≃ 30 solar mass suggest progenitor stars of low metallicity (Z/ZO ≤ 0.3). In this talk I will provide constrains on where the progenitors of GW150914 and GW170104 may have formed, based on advanced models of galaxy formation and evolution combined with binary population synthesis models. First I will combine estimates of galaxy properties (star-forming gas metallicity, star formation rate and merger rate) across cosmic time to predict the low redshift BBH merger rate as a function of present day host galaxy mass, formation redshift of the progenitor system and different progenitor metallicities. I will show that the signal is dominated by binaries formed at the peak of star formation in massive galaxies with and binaries formed recently in dwarf galaxies. Then, I will present what very high resolution hydrodynamic simulations of different galaxy types can learn us about their black hole populations

SHARP: A Spatially Higher-order, Relativistic Particle-in-Cell Code

Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to fifth order). We validate our algorithm against several test problems -- thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutions of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation.Comment: 26 pages, 19 figures, discussion about performance is added, published in Ap

Growth of beam-plasma instabilities in the presence of background inhomogeneity

We explore how inhomogeneity in the background plasma number density alters the growth of electrostatic unstable wavemodes of beam plasma systems. This is particularly interesting for blazar-driven beam-plasma instabilities, which may be suppressed by inhomogeneities in the intergalactic medium as was recently claimed in the literature. Using high resolution Particle-In-Cell simulations with the SHARP code, we show that the growth of the instability is local, i.e., regions with almost homogeneous background density will support the growth of the Langmuir waves, even when they are separated by strongly inhomogeneous regions, resulting in an overall slower growth of the instability. We also show that if the background density is continuously varying, the growth rate of the instability is lower; though in all cases, the system remains within the linear regime longer and the instability is not extinguished. In all cases, the beam loses approximately the same fraction of its initial kinetic energy in comparison to the uniform case at non-linear saturation. Thus, inhomogeneities in the intergalactic medium are unlikely to suppress the growth of blazar-driven beam-plasma instabilities.Comment: 10 pages, 6 figures, Accepted by ApJ, comments welcom

Importance of resolving the spectral support of beam-plasma instabilities in simulations

Many astrophysical plasmas are prone to beam-plasma instabilities. For relativistic and dilute beams, the {\it spectral} support of the beam-plasma instabilities is narrow, i.e., the linearly unstable modes that grow with rates comparable to the maximum growth rate occupy a narrow range of wave numbers. This places stringent requirements on the box-sizes when simulating the evolution of the instabilities. We identify the implied lower limits on the box size imposed by the longitudinal beam plasma instability, i.e., typically the most stringent condition required to correctly capture the linear evolution of the instabilities in multidimensional simulations. We find that sizes many orders of magnitude larger than the resonant wavelength are typically required. Using one-dimensional particle-in-cell simulations, we show that the failure to sufficiently resolve the spectral support of the longitudinal instability yields slower growth and lower levels of saturation, potentially leading to erroneous physical conclusion.Comment: 7 pages, 9 figures, accepted by Ap