Search

With the growing number of detections of binary black hole mergers, we are beginning to probe structure in the distribution of masses. A recent study by Schneider et al. proposes that isolated binary evolution of stripped stars naturally gives rise to the peaks at chirp masses $\sim 8 M_\odot$, $14 M_\odot$ in the chirp mass distribution and explains the dearth of black holes between $\approx 10-12 M_\odot$ in chirp mass. The gap in chirp mass results from an apparent gap in the component mass distribution between $m_1, m_2 \approx 10-15 M_\odot$ and the specific pairing of these black holes. This component mass gap results from the variation in core compactness of the progenitor, where a drop in compactness of Carbon-Oxygen core mass will no longer form black holes from core collapse. We develop a population model motivated by this scenario to probe the structure of the component mass distribution of binary black holes consisting of two populations: 1) two peak components to represent black holes formed in the compactness peaks, and 2) a powerlaw component to account for any polluting events, these are binaries that may have formed from different channels (e.g. dynamical). We perform hierarchical Bayesian inference to analyse the events from the third gravitational-wave transient catalogue (GWTC-3) with this model. We find that there is a preference for the lower mass peak to drop off sharply at $\sim 11 M_\odot$ and the upper mass peak to turn on at $\sim 13 M_\odot$, in line with predictions from Schneider et al. However, there is no clear evidence for a gap. We also find mild support for the two populations to have different spin distributions. In addition to these population results, we highlight observed events of interest that differ from the expected population distribution of compact objects formed from stripped stars., Comment: 11 pages, 8 figures

We present results from a systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, using data from the Palomar Gattini IR (PGIR) survey. R Coronae Borealis stars are dusty, erratic variable stars presumably formed from the merger of a He-core and a CO-core white dwarf (WD). PGIR is a 30 cm $J$-band telescope with a 25 deg$^{2}$ camera that surveys 18000 deg$^{2}$ of the northern sky ($\delta>-28^{o}$) at a cadence of 2 days. Using PGIR J-band lightcurves for $\sim$60 million stars together with mid-IR colors from WISE, we selected a sample of 530 candidate RCB stars. We obtained near-IR spectra for these candidates and identified 53 RCB stars in our sample. Accounting for our selection criteria, we find that there are a total of $\approx350^{+150}_{-100}$ RCB stars in the Milky Way. Assuming typical RCB lifetimes, this corresponds to an RCB formation rate of 0.8 - 5 $\times$ 10$^{-3}$ yr$^{-1}$, consistent with observational and theoretical estimates of the He-CO WD merger rate. We searched for quasi-periodic pulsations in the PGIR lightcurves of RCB stars and present pulsation periods for 16 RCB stars. We also examined high-cadenced TESS lightcurves for RCB and the chemically similar, but dustless hydrogen-deficient carbon (dLHdC) stars. We find that dLHdC stars show variations on timescales shorter than RCB stars, suggesting that they may have lower masses than RCB stars. Finally, we identified 3 new spectroscopically confirmed and 12 candidate Galactic DY Per type stars - believed to be colder cousins of RCB stars - doubling the sample of Galactic DY Per type stars., Comment: accepted for publication in PASP

Galactic white dwarf binaries (WDBs) and black hole binaries (BHBs) will be gravitational wave (GW) sources for LISA. Their detection will provide insights into binary evolution and the evolution of our Galaxy through cosmic history. Here, we make predictions of the expected WDB and BHB population within our Galaxy. We combine predictions of the compact remnant binary populations expected by stellar evolution by using the detailed Binary Population and Spectral Synthesis code (BPASS) with a Milky Way analogue galaxy model from the Feedback In Realistic Environment (FIRE) simulations. We use \textsc{PhenomA} and \textsc{LEGWORK} to simulate LISA observations. Both packages make similar predictions that on average four Galactic BHBs and 673 Galactic WDBs above the signal-to-noise ratio (SNR) threshold of 7 after a four-year mission. We compare these predictions to earlier results using the Binary Star Evolution (BSE) code with the same FIRE model galaxy. We find that BPASS predicts a few more LISA observable Galactic BHBs and a twentieth of the Galactic WDBs. The differences are due to the different physical assumptions that have gone into the binary evolution calculations. These results indicate that the expected population of compact binaries that LISA will detect depends very sensitively on the binary population synthesis models used and thus observations of the LISA population will provide tight constraints on our modelling of binary stars. Finally, from our synthetic populations we have created mock LISA signals that can be used to test and refine data processing methods of the eventual LISA observations.

Current gravitational wave (GW) detection pipelines for compact binary coalescence based on matched-filtering have reported over 90 confident detections during the first three observing runs of the LIGO-Virgo-KAGRA (LVK) detector network. Decreasing the latency of detection, in particular for future detectors anticipated to have high detection rates, remains an ongoing effort. In this paper, we develop and test a sparse dictionary learning (SDL) algorithm for the rapid detection of GWs. We evaluate the algorithms biases and estimate its GW detection rate for an astrophysical population of binary black holes. The SDL algorithm is assessed using both, simulated data injected into the proposed A+ detector sensitivity and real data containing confident detections from the third LVK observing run. We find that our SDL algorithm can reconstruct a single binary black hole signal in less than 1 s. This suggests that SDL could be regarded as a promising approach for rapid, efficient GW detection in future observing runs of ground-based detectors.

In this work, we investigate the impact of tidal torques and mass transfer on the population of double white dwarfs (DWDs) that will be observed with LISA. Starting from a distribution of DWDs at formation predicted by numerical simulations, we use a semi-analytical model to evolve DWDs under different hypotheses for the efficiency of tidal coupling and the birth spins of white dwarfs. We then estimate the stochastic foreground and the population of resolvable binaries for LISA in each scenario. Our predicted DWD binary distribution can differ substantially from the distribution expected if only gravitational waves (GWs) are considered. If white dwarfs spin slowly, then we predict an excess of systems around a few mHz, due to binaries that outspiral after the onset of mass transfer. This excess of systems leads to differences in the confusion noise, which are most pronounced for strong tidal coupling. In that case, we find a significantly higher number of resolvable binaries than in the GW-only scenario. If instead white dwarfs spin rapidly and tidal coupling is weak, then we find no excess around a few mHz, and the confusion noise due to DWDs is very small. In that scenario, we also predict a subpopulation of outspiralling binaries below 0.1 mHz. Using the Fisher matrix approximation, we estimate the uncertainty on the GW-frequency derivative of resolvable systems. We estimate that, even for non-accreting systems, the mismodelling error due to neglect of effects other than GWs is larger than the statistical uncertainty, and thus this neglect would lead to biased estimates for mass and distance. Our results highlight the need for flexible tools in LISA data analysis. Because our semi-analytical model hinges upon a simplistic approach to determining the stability of mass accretion it will be important to deepen our comprehension of stability in mass-transferring DWD binaries., Comment: 11 pages, 7 figures, with appendix 12 and 9 respectively

The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase., Comment: 155 pages, with executive summary and table of contents

After the nearly hundred gravitational-wave detections reported by the LIGO-Virgo-KAGRA Collaboration, the question of the cosmological origin of merging binary black holes (BBHs) remains open. The two main formation channels generally considered are from isolated field binaries or via dynamical assembly in dense star clusters. Here, we focus on understanding the dynamical formation of merging BBHs within massive clusters in galaxies of different masses. To this end, we apply a new framework to consistently model the formation and evolution of massive star clusters in zoom-in cosmological simulations of galaxies. Each simulation, taken from the FIRE project, provides a realistic star formation environment with a unique star formation history and hosts realistic giant molecular clouds that constitute the birthplace of star clusters. Combined with the code for star cluster evolution CMC, we are able to produce populations of dynamically formed merging BBHs across cosmic time in different environments. As the most massive star clusters preferentially form in dense massive clouds of gas, we find that, despite their low metallicities favourable to the creation of black holes, low-mass galaxies contain few massive clusters and therefore have a limited contribution to the global production of dynamically formed merging BBHs. Furthermore, we find that massive clusters can host hierarchical BBH mergers with clear identifiable physical properties. Looking at the evolution of the BBH merger rate in different galaxies, we find strong correlations between BBH mergers and the most extreme episodes of star formation. Finally, we discuss the implications for future LIGO-Virgo-KAGRA gravitational wave observations., Comment: 14 pages, 9 figures, 3 tables

We investigate the Stochastic Gravitational Wave Background (SGWB) produced by merging binary black holes (BBHs) and binary neutron stars (BNSs) in the frequency ranges of LIGO/Virgo/Kagra and LISA. We develop three analytical models, that are calibrated to the measured local merger rates, and complement them with three population synthesis models based on the COSMIC code. We discuss the uncertainties, focusing on the impact of the BBH mass distribution, the effect of the metallicity of the progenitor stars and the time delay distribution between star formation and compact binary merger. We also explore the effect of uncertainties in binary stellar evolution on the background. For BBHs, our analytical models predict $\Omega_{GW}$ in the range $[4.10^{-10}-1.10^{-9}]$ (25 Hz) and $[1.10^{-12}-4.10^{-12}]$ (3 mHz), and between $[2.10^{-10}-2.10^{-9}]$ (25 Hz) and $[7.10^{-13}- 7.10^{-12}]$ (3 mHz) for our population synthesis models. This background is unlikely to be detected during the LIGO/Virgo/Kagra O4 run, but could be detectable with LISA. We predict about 10 BBH and no BNS mergers that could be individually detectable by LISA for a period of observation of 4 years. Our study provides new insights into the population of compact binaries and the main sources of uncertainty in the astrophysical SGWB., Comment: 10 pages, 8 figures, accepted 2023 September 18 to MNRAS

With nearly a hundred gravitational wave detections, the origin of black hole mergers has become a key question. Here, we focus on understanding the typical galactic environment in which binary black hole mergers arise. To this end, we synthesize progenitors of binary black hole mergers as a function of the redshift of progenitor formation, present-day formation galaxy mass, and progenitor stellar metallicity for $240$ star formation and binary evolution models. We provide guidelines to infer the formation galaxy properties and time of formation, highlighting the interplay between the star formation rate and the efficiency of forming merging binary black holes from binary stars, both of which strongly depend on metallicity. We find that across models, over 50% of BBH mergers have a progenitor metallicity of a few tenths of Solar metallicity, however, inferring formation galaxy properties strongly depends on both the binary evolution model and global metallicity evolution. The numerous, low-mass black holes ($\mathrm{\lesssim 15\,M_{\odot}}$) trace the bulk of the star formation in galaxies heavier than the Milky Way ($M_\mathrm{Gal}$ $\mathrm{\gtrsim 10^{10.5}\,M_{\odot}}$). In contrast, heavier BBH mergers typically stem from larger black holes forming in lower metallicity dwarf galaxies ($M_\mathrm{Gal}$ $\mathrm{\lesssim 10^{9}\,M_{\odot}}$). We find that the progenitors of detectable binary black holes tend to arise from dwarf galaxies at a lower formation redshift ($\lesssim \, 1$). We also produce a posterior probability of the progenitor environment for any detected gravitational wave signal. For the massive GW150914 merger, we show that it likely came from a very low metallicity ($Z$ $\mathrm{\lesssim}\,0.025\,\mathrm{Z_{\odot}}$) environment., Comment: 17 pages, 15 figures, 1 table, Accpeted for publication by MNRAS

Massive colliding-wind binaries that host a Wolf-Rayet (WR) star present a potentially important source of dust and chemical enrichment in the interstellar medium (ISM). However, the chemical composition and survival of dust formed from such systems is not well understood. The carbon-rich WR (WC) binary WR~140 presents an ideal astrophysical laboratory for investigating these questions given its well-defined orbital period and predictable dust-formation episodes every 7.93 years around periastron passage. We present observations from our Early Release Science program (ERS1349) with the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) Medium-Resolution Spectrometer (MRS) and Imager that reveal the spectral and spatial signatures of nested circumstellar dust shells around WR~140. MIRI MRS spectroscopy of the second dust shell and Imager detections of over 17 shells formed throughout the past $\gtrsim130$ years confirm the survival of carbonaceous dust grains from WR~140 that are likely carriers of "unidentified infrared" (UIR)-band features at 6.4 and 7.7 $\mu$m. The observations indicate that dust-forming WC binaries can enrich the ISM with organic compounds and carbonaceous dust., Comment: Published in Nature Astronomy on Oct 12, 2022; 21 pages, 5 figures, 2 tables

With its last observing run, the LIGO, Virgo, and KAGRA collaboration has detected almost one hundred gravitational waves from compact binary coalescences. A common approach to studying the population properties of the observed binaries is to use phenomenological models to describe the spin, mass, and redshift distributions. More recently, with the aim of providing a clearer link to astrophysical processes forming the observed compact binaries coalescences, several authors have proposed to employ synthetic catalogs for population studies. In this paper, we review how to employ and interpret synthetic binary catalogs for gravitational-wave progenitors studies. We describe how to build multi-channel merger rates and describe their associated probabilities focusing on stellar progenitor properties. We introduce a method to quantify the match between the phenomenological reconstruction of merger rates with synthetic catalogs. We detail the implementation of synthetic catalogs for multi-channel hierarchical Bayesian inference, highlighting computational aspects and issues related to hyper-prior choice. We find that when inferring stellar progenitors' properties from gravitational-wave observations, the relative efficiency in compact objects production should be taken into account. Finally, by simulating binary black hole detections with LIGO and Virgo sensitivity expected for the O4 observing run, we present two case studies related to the inference of the common envelope efficiency and progenitor metallicity of the binary black holes. We finally discuss how progenitors' properties can be linked to binary black hole properties., Comment: 13 pages, 12 figures, submitted to MNRAS

The current generation of galaxy simulations can resolve individual giant molecular clouds, the progenitors of dense star clusters. But the evolutionary fate of these young massive clusters, and whether they can become the old globular clusters (GCs) observed in many galaxies, is determined by a complex interplay of internal dynamical processes and external galactic effects. We present the first star-by-star $N$-body models of massive ($N\sim10^5-10^7$) star clusters formed in a FIRE-2 MHD simulation of a Milky Way-mass galaxy, with the relevant initial conditions and tidal forces extracted from the cosmological simulation. We select 895 ($\sim 30\%$) of the YMCs with $ > 6\times10^4M_{\odot}$ from Grudi\'c et al.~2022 and integrate them to $z=0$ using the Cluster Monte Carlo Code, \texttt{CMC}. This procedure predicts a MW-like system with 148 GCs, predominantly formed during the early, bursty mode of star formation. Our GCs are younger, less massive, and more core-collapsed than clusters in the Milky Way or M31. This results from the assembly history and age-metallicity relationship of the host galaxy: younger clusters are preferentially born in stronger tidal fields and initially retain fewer stellar-mass black holes, causing them to lose mass faster and reach core collapse sooner than older GCs. Our results suggest that the masses and core/half-light radii of GCs are shaped not only by internal dynamical processes, but also by the specific evolutionary history of their host galaxies. These results emphasize that $N$-body studies with realistic stellar physics are crucial to understanding the evolution and present-day properties of GC systems., Comment: 24 pages, 16 figures, matches version accepted by MNRAS

The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.

The properties of young star clusters formed within a galaxy are thought to vary in different interstellar medium (ISM) conditions, but the details of this mapping from galactic to cluster scales are poorly understood due to the large dynamic range involved in galaxy and star cluster formation. We introduce a new method for modeling cluster formation in galaxy simulations: mapping giant molecular clouds (GMCs) formed self-consistently in a FIRE-2 MHD galaxy simulation onto a cluster population according to a GMC-scale cluster formation model calibrated to higher-resolution simulations, obtaining detailed properties of the galaxy's star clusters in mass, metallicity, space, and time. We find $\sim 10\%$ of all stars formed in the galaxy originate in gravitationally-bound clusters overall, and this fraction increases in regions with elevated $\Sigma_{\rm gas}$ and $\Sigma_{\rm SFR}$, because such regions host denser GMCs with higher star formation efficiency. These quantities vary systematically over the history of the galaxy, driving variations in cluster formation. The mass function of bound clusters varies -- no single Schechter-like or power-law distribution applies at all times. In the most extreme episodes, clusters as massive as $7\times 10^6 M_\odot$ form in massive, dense clouds with high star formation efficiency. The initial mass-radius relation of young star clusters is consistent with an environmentally-dependent 3D density that increases with $\Sigma_{\rm gas}$ and $\Sigma_{\rm SFR}$. The model does not reproduce the age and metallicity statistics of old ($>11\rm Gyr$) globular clusters found in the Milky Way, possibly because it forms stars more slowly at $z>3$., Comment: Accepted to MNRAS

In its observation band, the Laser Interferometer Space Antenna (LISA) will simultaneously observe stochastic gravitational-wave background (SGWB) signals of different origins; orbitally modulated waveforms from galactic white dwarf binaries, a binary black hole produced background, and possibly a cosmologically produced SGWB. We simulate the emission of gravitational waves from galactic white dwarf binaries based on the Lamberts \cite{Lamberts} distributions and determine a complex waveform from the galactic foreground. We generate the modulated galactic signal detected by LISA due to its orbital motion, and present a data analysis strategy to address it. The Fisher Information and Markov Chain Monte Carlo methods give an estimate of the LISA noise and parameters for the different signal sources. We simultaneously estimate the galactic foreground, the astrophysical and cosmological backgrounds, and estimate detection limits for the future LISA observation of the SGWB in the spectral domain with the 3 LISA channels $ A $, $ E $ and $ T $. In the context of the expected astrophysical background and a galactic foreground, a cosmological background energy density of about $ \Omega_{GW,Cosmo} \approx 8 \times 10^{-13} $ could be detected by LISA with our spectral separation strategy., Comment: Contribution to the 2021 Gravitation session of the 55th Rencontres de Moriond

Within its observational band the Laser Interferometer Space Antenna, LISA, will simultaneously observe orbital modulated waveforms from Galactic white dwarf binaries, a binary black hole produced gravitational-wave background, and potentially a cosmologically created stochastic gravitational-wave background (SGWB). The overwhelming majority of stars end their lives as white dwarfs, making them very numerous in the Milky Way. We simulate Galactic white dwarf binary gravitational-wave emission based on distributions from various mock catalogs and determine a complex waveform from the Galactic foreground with $3.5 \times 10^{7}$ binaries. We describe the effects from the Galactic binary distribution population across mass, position within the Galaxy, core type, and orbital frequency distribution. We generate the modulated Galactic white dwarf signal detected by \textit{LISA} due to its orbital motion, and present a data analysis strategy to address it. The Fisher Information and Markov Chain Monte Carlo methods give an estimation of the \textit{LISA} noise and the parameters for the different signal classes. We estimate the detectable limits for the future LISA observation of the SGWB in the spectral domain with the 3 \textit{LISA} channels $A$, $E$, and $T$. We simultaneously estimate the Galactic foreground, the astrophysical and cosmological backgrounds. Assuming the expected astrophysical background and a Galactic foreground, a cosmological background energy density of around $\Omega_{GW,Cosmo} \approx 8 \times 10^{-13}$ could be detected by LISA. LISA will either detect a cosmologically produced SGWB, or set a limit that will have important consequences.

We present a dust spectral energy distribution (SED) and binary stellar population analysis revisiting the dust production rates (DPRs) in the winds of carbon-rich Wolf-Rayet (WC) binaries and their impact on galactic dust budgets. DustEM SED models of 19 Galactic WC ``dustars" reveal DPRs of $\dot{M}_d\sim10^{-10}-10^{-6}$ M$_\odot$ yr$^{-1}$ and carbon dust condensation fractions, $\chi_C$, between $0.002 - 40\%$. A large ($0.1 - 1.0$ $\mu$m) dust grain size composition is favored for efficient dustars where $\chi_C\gtrsim1\%$. Results for dustars with known orbital periods verify a power-law relation between $\chi_C$, orbital period, WC mass-loss rate, and wind velocity consistent with predictions from theoretical models of dust formation in colliding-wind binaries. We incorporated dust production into Binary Population and Spectral Synthesis (BPASS) models to analyze dust production rates from WC dustars, asymptotic giant branch stars (AGBs), red supergiants (RSGs), and core-collapse supernovae (SNe). BPASS models assuming constant star formation (SF) and a co-eval $10^6$ M$_\odot$ stellar population were performed at low, Large Magellanic Cloud (LMC)-like, and solar metallicities (Z = 0.001, 0.008, and 0.020). Both constant SF and co-eval models show that SNe are net dust destroyers at all metallicities. Constant SF models at LMC-like metallicities show that AGB stars slightly outproduce WC binaries and RSGs by factors of $2-3$, whereas at solar metallicites WC binaries are the dominant source of dust for $\sim60$ Myr until the onset of AGBs, which match the dust input of WC binaries. Co-eval population models show that for "bursty" SF, AGB stars dominate dust production at late times ($t\gtrsim 70$ Myr)., Comment: 34 pages, 12 figures, 7 tables, Accepted for publication in ApJ

We explore the possibility that GW190412, a binary black hole merger with a non-equal-mass ratio and significantly spinning primary, was formed through repeated black hole mergers in a dense super star cluster. Using a combination of semi-analytic prescriptions for the remnant spin and recoil kick of black hole mergers, we show that the mass ratio and spin of GW190412 are consistent with a binary black hole whose primary component has undergone two successive mergers from a population of $\sim 10M_{\odot}$ black holes in a high-metallicity environment. We then explore the production of GW190412-like analogs in the CMC Cluster Catalog, a grid of 148 $N$-body star cluster models, as well as a new model, behemoth, with nearly $10^7$ particles and initial conditions taken from a cosmological MHD simulation of galaxy formation. We show that the production of binaries with GW190412-like masses and spins is dominated by massive super star clusters with high metallicities and large central escape speeds. While many are observed in the local universe, our results suggest that a careful treatment of these massive clusters, many of which may have been disrupted before the present day, is necessary to characterize the production of unique gravitational-wave events produced through dynamics., Comment: 9 pages, 4 figures, matches version accepted by ApJL

We study the longitudinal stability of beam-plasma systems in the presence of a density inhomogeneity in the background plasma. Previous works have focused on the non-relativistic regime where hydrodynamical models are used to evolve pre-existing Langmuir waves within inhomogeneous background plasmas. Here, for the first time we study the problem with kinetic equations in a fully-relativistic way. We do not assume the existence of Langmuir waves, and we focus on the rate and the mechanism by which waves are excited in such systems from an initial perturbation. We derive the structure of the unstable modes and compute an analytical approximation for their growth rates. Our computation is limited to dilute and cold beams, and shows an excellent agreement with particle-in-cell simulations performed using the SHARP code. We show that, due to such an inhomogeneity, the virulent beam-plasma instabilities in the intergalactic medium are not suppressed but their counterparts in the solar wind can be suppressed as evidenced by propagating type-III solar radio bursts., Comment: 23 pages, 5 figures, accepted for publication in the Journal of Plasma Physics; comments are welcome

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work., Comment: 7 pages, Submitted to the Astro2020 call for APC white papers

We determine the main properties of the Galactic binary black hole (BBH) population detectable by LISA and strategies to distinguish them from the much more numerous white dwarf binaries. We simulate BBH populations based on cosmological simulations of Milky Way-like galaxies and binary evolution models. We then determine their gravitational wave emission as observed by LISA and build mock catalogs. According to our model LISA will detect $\approx4(6)$ binary black holes assuming 4(10) years of operations. Those figures grow to $\approx6(9)$ when models are re-normalized to the inferred LIGO/Virgo merger rates. About 40\%(70\%) of the sources will have a good enough chirp mass measurement to separate them from the much lighter white dwarf and neutron star binaries. Most of the remaining sources should be identifiable by their lack of electromagnetic counterpart within $\approx100$ pc. These results are robust with respect to the current uncertainties of the BBH merger rate as measured by LIGO/Virgo as well as the global mass spectrum of the binaries. We determine there is a 94 per cent chance that LISA finds at least one of these systems, which will allow us to pinpoint the conditions where they were formed and possibly find unique electromagnetic signatures., Comment: 6 pages, 5 figures, submitted to MNRAS Letters

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe., Comment: 28 pages, 8 figures; White Paper submitted to ESA's Voyage 2050 call for papers on behalf of the LISA Consortium 2050 Task Force

Short-period (P<1 hour) white dwarf binaries will be the most numerous sources for the space-based gravitational wave detector LISA. Based on thousands of resolved systems, we will be able to constrain binary evolution and provide a new map of the Milky Way and its surroundings. Here we predict the main properties of populations of different types of detached white dwarf binaries detected by LISA. For the first time, we combine a high-resolution cosmological simulation of a Milky Way-mass galaxy (from the FIRE project) with a binary population synthesis model for low and intermediate mass stars. Our model therefore provides a cosmologically realistic star formation and metallicity history for the galaxy and naturally produces its different components such as the thin and thick disk, the bulge, the stellar halo, and satellite galaxies and streams. With the simulation, we show how different galactic components contribute differently to the gravitational wave signal, due to their typical age and distance distributions. We find that the dominant LISA sources will be He-He systems and He-CO systems with important contributions from the thick disk and bulge but also a few systems in the stellar halo. The resulting sky map of the sources is different from previous models, with important consequences for the searches for electromagnetic counterparts and data analysis. We also emphasize that much of the science-enabling information regarding white dwarf binaries, such as the chirp mass and the sky localisation, becomes increasingly rich with long observations, including an extended mission up to 8 years., Comment: 16 pages, 15 figures, submitted to MNRAS. Comments are welcome

An intergalactic magnetic field stronger than $3\times10^{-13}$~G would explain the lack of a bright, extended degree-scale, GeV-energy inverse Compton component in the gamma-ray spectra of TeV-blazars. A robustly predicted consequence of the presence of such a field is the existence of degree-scale GeV-energy gamma-ray halos -- gamma-ray bow ties -- about TeV-bright active galactic nuclei, corresponding to more than half of all radio galaxies. However, the emitting regions of these halos are confined to and aligned with the direction of the relativistic jets associated with gamma-ray sources. Based on the orientation of radio jets, we align and stack corresponding degree-scale gamma-ray images of isolated Fanaroff-Riley class I and II objects and exclude the existence of these halos at overwhelming confidence, limiting the intergalactic field strength to $<10^{-15}$~G for large-scale fields and progressively larger in the diffusive regime when the correlation length of the field becomes small in comparison to 1 Mpc. When combined with prior limits on the strength of the intergalactic magnetic field, this excludes a purely magnetic explanation for the absence of halos. Thus, it requires the existence of novel physical processes that preempt the creation of halos, e.g., the presence of beam-plasma instabilities in the intergalactic medium or a drastic cutoff of the very high energy spectrum of these sources., Comment: Accepted to ApJ, 20 pages, 15 figures

We present a new set of high-resolution hydrodynamic cosmological zoom-in simulations that apply the Feedback In Realistic Environments (FIRE) physics to both Local Group (LG)-like and isolated Milky Way (MW)-like volumes (ten host systems in total with baryonic particle mass $\simeq 3,500-7,000\,M_\odot$). We study the stellar mass functions, circular velocity or mass profiles, and velocity dispersions of the dwarf galaxy populations. The simulations reproduce the stellar mass function and central densities of MW satellite dwarfs for $M_\ast \geq10^{5.5}\,M_\odot$ and predict the existence of $\sim3$ unidentified galaxies with $M_\ast\sim10^5\,M_\odot$ within $300$ kpc of the MW. Overall, we find no evidence for the classical missing satellites or too-big-to-fail (TBTF) problems for satellite galaxies in our sample. Among the satellites, TBTF is resolved primarily by subhalo disruption and overall mass loss; central density profiles of subhalos are of secondary importance. For non-satellite galaxies, our LG-like simulations predict as many as $\sim10$ as-of-yet unseen galaxies at distances $0.3-1$ Mpc from both hosts, with $M_\ast\simeq10^{5-6}\,M_\odot$ (in halos with $V_\mathrm{max}\sim 20~\mathrm{km\,s}^{-1}$), albeit with large halo-to-halo variance. None of our simulations produces a compact, baryon-dominated, high-density dwarf elliptical-type galaxy (with $V_\mathrm{circ} \gtrsim 35~\mathrm{km\,s}^{-1}$ at $r<1\,$kpc), of which six may appear in the LG (but none in the MW). It may therefore remain a challenge to reproduce the full diversity of the dwarf population, including both the highest and lowest density systems., Comment: Main text: 14 pages, 6 figures, 1 table. Appendices: 3 figures, 1 table. Submitted to MNRAS

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 welcome

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"

Binary black holes are the primary endpoint of massive stellar evolution. Their properties provide a unique opportunity to constrain binary evolution, which is still poorly understood. In this paper, we predict the inventory of binary black holes and their merger products in/around the Milky Way, and detail their main properties. We present the first combination of a high-resolution cosmological simulation of a Milky Way-mass galaxy with a binary population synthesis model. The hydrodynamic simulation, taken from the FIRE project, provides a cosmologically realistic star formation history for the galaxy and its stellar halo and satellites. We apply a metallicity-dependent evolutionary model to the star particles to produce individual binary black holes. We find that a million binary black holes have merged in the model Milky Way, and 3 million binaries are still present, with an average mass of 28 Msun per binary. Because the black hole progenitors are biased towards low metallicity stars, half reside in the stellar halo and satellites and 40 per cent of the binaries were formed outside the main galaxy. This trend increases with the masses of the black holes. The numbers and mass distribution of the merged systems is compatible with the LIGO/Virgo detections. Observations of these black holes will be challenging, both with electromagnetic methods and LISA. We find that a cosmologically realistic star formation history, with self-consistent metal enrichment and Galactic accretion history, are key ingredients for determining binary black hole rates that can be compared with observations to constrain massive binary evolution., Comment: accepted in MNRAS, data available at https://fire.northwestern.edu/data/

Young massive star clusters spanning $\sim 10^4 - 10^8 M_\odot$ in mass have been observed to have similar surface brightness profiles. Recent hydrodynamical simulations of star cluster formation have also produced star clusters with this structure. We argue analytically that this type of mass distribution arises naturally in the relaxation from a hierarchically-clustered distribution of stars into a monolithic star cluster through hierarchical merging. We show that arbitrary initial profiles will tend to converge to a universal profile under hierarchical merging, owing to phase-space mixing obeying certain conservation constraints. We perform $N$-body simulations of a pairwise merger of model star clusters and find that mergers readily produce the shallow surface brightness profiles observed in young massive clusters. Finally, we simulate the relaxation of a hierarchically-clustered mass distribution constructed from an idealized fragmentation model. Assuming only power-law spatial and kinematic scaling relations, these numerical experiments are able to reproduce the surface density profiles of observed young massive star clusters. Thus we provide physical motivation for the structure of young massive clusters within the paradigm of hierarchical star formation. This has important implications for the structure of nascent globular clusters., Comment: 16 pages, 10 figures

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 ApJ

In this contribution we review the recent progress in the modeling of Pulsar Wind Nebulae (PWN). We start with a brief overview of the relevant physical processes in the magnetosphere, the wind-zone and the inflated nebula bubble. Radiative signatures and particle transport processes obtained from 3D simulations of PWN are discussed in the context of optical and X-ray observations. We then proceed to consider particle acceleration in PWN and elaborate on what can be learned about the particle acceleration from the dynamical structures called "wisps" observed in the Crab nebula. We also discuss recent observational and theoretical results of gamma-ray flares and the inner knot of the Crab nebula, which had been proposed as the emission site of the flares. We extend the discussion to GeV flares from binary systems in which the pulsar wind interacts with the stellar wind from a companion star. The chapter concludes with a discussion of solved and unsolved problems posed by PWN., Comment: To appear in "Jets and Winds in Pulsar Wind Nebulae, Gamma-ray Bursts and Blazars: Physics of Extreme Energy Release" of the Space Science Reviews series. The final publication is available at Springer via http://dx.doi.org/10.1007/s11214-017-0344-x

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 ApJ

The X-ray light curves of the early afterglow phase from gamma-ray bursts present a puzzling variability, including flares. The origin of these flares is still debated, and often associated with a late activity of the central engine. We discuss an alternative scenario where the central engine remains short-lived and flares are produced by the propagation of a long-lived reverse shock in a stratified ejecta. Here we focus on the hydrodynamics of the shock interactions. We perform one-dimensional ultrarelativistic hydrodynamic simulations with different initial internal structure in the gamma-ray burst ejecta. We use them to extract bolometric light curves and compare with a previous study based on a simplified ballistic model. We find a good agreement between both approaches, with similar slopes and variability in the light curves, but identify several weaknesses in the ballistic model: the density is underestimated in the shocked regions, and more importantly, late shock reflections are not captured. With accurate dynamics provided by our hydrodynamic simulations, we confirm that internal shocks in the ejecta lead to the formation of dense shells. The interaction of the long-lived reverse shock with a dense shell then produces a fast and intense increase of the dissipated power. Assuming that the emission is due to the synchrotron radiation from shock-accelerated electrons, and that the external forward shock is radiatively inefficient, we find that this results in a bright flare in the X-ray lightcurve, with arrival times, shapes, and duration in agreement with the observed properties of X-ray flares in GRB afterglows., Comment: Accepted for publication in MNRAS

Many-degree-scale gamma-ray halos are expected to surround extragalactic high-energy gamma ray sources. These arise from the inverse Compton emission of an intergalactic population of relativistic electron/positron pairs generated by the annihilation of >100 GeV gamma rays on the extragalactic background light. These are typically anisotropic due to the jetted structure from which they originate or the presence of intergalactic magnetic fields. Here we propose a novel method for detecting these inverse-Compton gamma-ray halos based upon this anisotropic structure. Specifically, we show that by stacking suitably defined angular power spectra instead of images it is possible to robustly detect gamma-ray halos with existing Fermi Large Area Telescope (LAT) observations for a broad class of intergalactic magnetic fields. Importantly, these are largely insensitive to systematic uncertainties within the LAT instrumental response or associated with contaminating astronomical sources., Comment: Submitted to ApJ, 22 pages, 27 figures, 9 pages in appendix

Pair creation on the cosmic infrared background and subsequent inverse-Compton scattering on the CMB potentially reprocesses the TeV emission of blazars into faint GeV halos with structures sensitive to intergalactic magnetic fields (IGMF). Previous work has shown that these halos are then highly-anisotropic and extended. If the coherence length of the IGMF is greater than the inverse-Compton cooling length of the pairs, then the orientation of the gamma-ray halo will be correlated with the direction of the magnetic field which is unknown and expected to change for each source. In order to constructively add each source we then use angular power spectra which are insensitive to the jet orientation. By looking at known GeV blazars detected by Fermi, we exclude the existence of an IGMF with coherence lengths $>100$ Mpc at greater than $3.9\sigma$ with current-day strengths in the range $10^{-16}$ to $10^{-15}$ G, and at 2$\sigma$ from $10^{-17}$ to $10^{-14}$ G. This provides a direct measurement of the non-existence of gamma-ray halos, providing an important check on previous results., Comment: 12 pages, 14 figures, 1 appendix. Accepted to ApJ

The annihilation of TeV photons from extragalactic TeV sources and the extragalactic background light produces ultrarelativistic $e^{\pm}$ beams, which are subject to powerful plasma instabilities that sap their kinetic energy. Here we study the linear phase of the plasma instabilities that these pair beams drive. To this end, we calculate the linear growth rate of the beam plasma and oblique instability in the electrostatic approximation in both the reactive and kinetic regimes, assuming a Maxwell-J{\"u}ttner distribution for the pair beam. We reproduce the well-known reactive and kinetic growth rates for both the beam plasma and oblique mode. We demonstrate for the oblique instability that there is a broad spectrum of unstable modes that grow at the maximum rate for a wide range of beam temperatures and wave vector orientations relative to the beam. We also delineate the conditions for applicability for the reactive and kinetic regimes and find that the beam plasma mode transitions to the reactive regime at a lower Lorentz factor than the oblique mode due to a combination of their different scalings and the anisotropy of the velocity dispersions. Applying these results to the ultrarelativistic $e^{\pm}$ beams from TeV blazars, we confirm that these beams are unstable to both, the kinetic oblique mode and the reactive beam-plasma mode. These results are important in understanding how powerful plasma instabilities may sap the energy of the ultrarelativistic $e^{\pm}$ beams as they propagate through intergalactic space., Comment: 13 pages, 4 figure, accepted to ApJ

Extended inverse Compton halos are generally anticipated around extragalactic sources of gamma rays with energies above 100 GeV. These result from inverse Compton scattered cosmic microwave background photons by a population of high-energy electron/positron pairs produced by the annihilation of the high-energy gamma rays on the infrared background. Despite the observed attenuation of the high-energy gamma rays, the halo emission has yet to be directly detected. Here, we demonstrate that in most cases these halos are expected to be highly anisotropic, distributing the up-scattered gamma rays along axes defined either by the radio jets of the sources or oriented perpendicular to a global magnetic field. We present a pedagogical derivation of the angular structure in the inverse Compton halo and provide an analytic formalism that facilitates the generation of mock images. We discuss exploiting this fact for the purpose of detecting gamma-ray halos in a set of companion papers., Comment: Accepted by ApJ, 18 pages, 17 figures, 13 pages in appendix

The recent LIGO detection of gravitational waves (GW150914), likely originating from the merger of two $\sim 30 M_\odot$ black holes suggests progenitor stars of low metallicity ($[Z/Z_\odot] \lesssim 0.3$), constraining when and where the progenitor of GW150914 may have formed. We combine estimates of galaxy properties (metallicity, star formation rate and merger rate) across cosmic time to predict the low redshift black hole - black hole merger rate as a function of present day host galaxy mass, $M_\mathrm{gal}$, and the formation redshift of the progenitor system $z_\mathrm{form}$ for different progenitor metallicities $Z_\mathrm{c}$. At $Z_\mathrm{c}=0.1 Z_\odot$, the signal is dominated by binaries in massive galaxies with $z_\mathrm{form}\simeq 2$, with a small contribution from binaries formed around $z_\mathrm{form}\simeq 0.5$ in dwarf galaxies. For $Z_\mathrm{c}=0.01Z_\odot$, fast mergers are possible and very recent star formation in dwarfs likely dominates. Additional gravitational wave detections from merging massive black holes will provide constraints on the mass-metallicity relation and massive star formation at high redshifts., Comment: 5 pages, 3 figures. accepted in MNRAS Letters, comments welcome

Detailed modeling of the high-energy emission from gamma-ray binaries has been propounded as a path to pulsar wind physics. Fulfilling this ambition requires a coherent model of the flow and its emission in the region where the pulsar wind interacts with the stellar wind of its companion. We developed a code that follows the evolution and emission of electrons in the shocked pulsar wind based on inputs from a relativistic hydrodynamical simulation. The code is used to model the well-documented spectral energy distribution and orbital modulations from LS 5039. The pulsar wind is fully confined by a bow shock and a back shock. The particles are distributed into a narrow Maxwellian, emitting mostly GeV photons, and a power law radiating very efficiently over a broad energy range from X-rays to TeV gamma rays. Most of the emission arises from the apex of the bow shock. Doppler boosting shapes the X-ray and VHE lightcurves, constraining the system inclination to $i\approx 35^{\rm o}$. There is a tension between the hard VHE spectrum and the level of X-ray to MeV emission, which requires differing magnetic field intensities that are hard to achieve with a constant magnetisation $\sigma$ and Lorentz factor $\Gamma_{p}$ of the pulsar wind. Our best compromise implies $\sigma\approx 1$ and $\Gamma_{p}\approx 5\times 10^{3}$, respectively higher and lower than the typical values in pulsar wind nebulae. The high value of $\sigma$ derived here, where the wind is confined close to the pulsar, supports the classical picture that has pulsar winds highly magnetised at launch. However, such magnetisations will require further investigations to be based on relativistic MHD simulations., Comment: 18 pages, 12 figures, under review with A&A (appendix with movies could not be uploaded to ArXiV)

TeV-blazars potentially heat the intergalactic medium (IGM) as their gamma rays interact with photons of the extragalactic background light to produce electron-positron pairs, which lose their kinetic energy to the surrounding medium through plasma instabilities. This results in a heating mechanism that is only weakly sensitive to the local density, and therefore approximately spatially uniform, naturally producing an inverted temperature-density relation in underdense regions. In this paper we go beyond the approximation of uniform heating and quantify the heating rate fluctuations due to the clustered distribution of blazars and how this impacts on the thermal history of the IGM. We analytically compute a filtering function that relates the heating rate fluctuations to the underlying dark matter density field. We implement it in the cosmological code GADGET-3 and perform large scale simulations to determine the impact of inhomogeneous heating. We show that, because of blazar clustering, blazar heating is inhomogeneous for z>= 2. At high redshift, the temperature-density relation shows an important scatter and presents a low temperature envelope of unheated regions, in particular at low densities and within voids. However, the median temperature of the IGM is close to that in the uniform case, albeit slightly lower at low redshift. We find that blazar heating is more complex than initially assumed and that the temperature-density relation is not unique. Our analytic model for the heating rate fluctuations couples well with large scale simulations and provides a cost-effective alternative to subgrid models.

Advances in X-ray and gamma-ray astronomy have opened a new window on our universe and revealed a wide variety of binaries composed of a compact object and a Be star. In Be X-ray binaries, a neutron star accretes the Be disk and truncates it through tidal interactions. Such systems have important X-ray outbursts, some related to the disk structure. In other systems, strong gamma ray emission is observed. In gamma-ray binaries, the neutron star is not accreting but driving a highly relativistic wind. The wind collision region presents similarities to colliding wind binaries composed of massive stars. The high energy emission is coming from particles being accelerated at the relativistic shock. I will review the physics of X-ray and gamma-ray binaries, focusing particularly on the recent developments on gamma-ray binaries. I will describe physical mechanisms such as relativistic hydrodynamics, tidal forces and non thermal emission. I will highlight how high energy astrophysics can shed a new light on Be star physics and vice-versa. A video of the talk can be found at http://activebstars.iag.usp.br/index.php/talk-conference-recordings/bestars-2014-in-london-ontario/session-6/video/astrid-lamberts-interacting-binaries-be-stars-and-high-energy-astrophysics., Comment: Proceeding of an invited review at the "Bright emissaries" conference held in London, Canada, Aug 10-13, 2014

Very-high energy gamma-rays from extragalactic sources pair-produce off of the extragalactic background light, yielding an electron-positron pair beam. This pair beam is unstable to various plasma instabilities, especially the "oblique" instability, which can be the dominant cooling mechanism for the beam. However, recently, it has been claimed that nonlinear Landau damping renders it physically irrelevant by reducing the effective damping rate to a low level. Here, we show with numerical calculations that the effective damping rate is $8\times 10^{-4}$ of the growth rate of the linear instability, which is sufficient for the "oblique" instability to be the dominant cooling mechanism of these pair beams. In particular, we show that previous estimates of this rate ignored the exponential cutoff in the scattering amplitude at large wavenumber and assumed that the damping of scattered waves entirely depends on collisions, ignoring collisionless processes. We find that the total wave energy eventually grows to approximate equipartition with the beam by increasingly depositing energy into long wavelength modes. As we have not included the effect of nonlinear wave-wave interactions on these long wavelength modes, this scenario represents the "worst-case" scenario for the oblique instability. As it continues to drain energy from the beam at a faster rate than other processes, we conclude that the "oblique" instability is sufficiently strong to make it the physically dominant cooling mechanism for high-energy pair beams in the intergalactic medium., Comment: 6 pages, 3 figures, accepted to ApJ

Gamma-ray binaries are composed of a massive star and a rotation-powered pulsar with a highly relativistic wind. The collision between the winds from both objects creates a shock structure where particles are accelerated, resulting in the observed high energy emission. We study the impact of special relativity on the structure and stability of the colliding wind region and highlight the differences with colliding winds from massive stars. We focus on evolution with increasing values of the Lorentz factor of the pulsar wind, keeping in mind that current simulations are unable to reach the expected values of the pulsar wind Lorentz factors by orders of magnitude. We use high resolution numerical simulations with a relativistic extension to the hydrodynamics code RAMSES we have developed. Using 2D simulations, we focus on the region close to the binary, neglecting orbital motion. We use different values of the Lorentz factor of the pulsar wind, up to 16. We find analytic scaling relations between stellar wind collisions and gamma-ray binaries. They provide the position of the contact discontinuity. The position of the shocks strongly depends on the Lorentz factor, the relativistic wind is more collimated than expected based on non-relativistic simulations. Beyond a certain distance, the shocked flow is accelerated to its initial velocity and follows adiabatic expansion. We provide guidance for extrapolation towards more realistic values of the Lorentz factor of the pulsar wind. We extended the adaptive mesh refinement code RAMSES to relativistic hydrodynamics. This code is suited for the study of astrophysical objects such as pulsar wind nebulae, gamma-ray bursts or relativistic jets and will be part of the next public release of RAMSES. Using this code we performed simulations of gamma-ray binaries, highlighting the limits and possibilities of current hydrodynamic models of such systems., Comment: Accepted for publication in Astronomy and Astrophysics

We present high resolution simulations with RAMSES of supersonic colliding stellar winds. The collision results in a double shock structure which is subject to different instabilities. The Kelvin-Helmholtz instability (KHI) introduces some mixing and variability. For isothermal winds, the Non-linear Thin Shell Instability violently affects the interaction region. Properly modelling these instabilities requires a high enough resolution and an adapted numerical method, especially when one of the winds strongly dominates the other one. At large scale, orbital motion is expected to turn the shocked zone into a spiral but we find that in some configurations the KHI may disrupt the spiral. A colliding wind structure is also expected in gamma-ray binaries composed of a massive star and a young pulsar which emits a highly relativistic wind. Numerical simulations are necessary to understand the geometry of such systems and should take into account the relativistic nature of the pulsar wind. We implemented a second order Godunov method to solve the equations of relativistic hydrodynamics with RAMSES, including the possibility for Adaptive Mesh Refinement. After a brief overview of our numerical implementation, we will present preliminary simulation that contrast the structure of the flow in gamma-ray and stellar binaries., Comment: Proceeding of Astronum 2012 Conference

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.4045