Source author record

Astrid Lamberts

Astrid Lamberts appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

Great Balls of FIRE I: The formation of star clusters across cosmic time in a Milky Way-mass galaxy

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 $Σ_{\rm gas}$ and $Σ_{\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 $Σ_{\rm gas}$ and $Σ_{\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$.

preprint2020arXiv

GW190412 as a Third-Generation Black Hole Merger from a Super Star Cluster

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.

preprint2020arXiv

Revisiting the Impact of Dust Production from Carbon-Rich Wolf-Rayet Binaries

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, $χ_C$, between $0.002 - 40\%$. A large ($0.1 - 1.0$ $μ$m) dust grain size composition is favored for efficient dustars where $χ_C\gtrsim1\%$. Results for dustars with known orbital periods verify a power-law relation between $χ_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).

preprint2020arXiv

The growth of the longitudinal beam-plasma instability in the presence of an inhomogeneous background

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.

preprint2019arXiv

Finding binary black holes in the Milky Way with LISA

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.

preprint2016arXiv

Bow Ties in the Sky I: The Angular Structure of Inverse Compton Gamma-ray Halos in the Fermi Sky

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.

preprint2016arXiv

The Linear Instability of Dilute Ultrarelativistic $e^{\pm}$ Pair Beams

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{ü}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.

preprint2016arXiv

When and where did GW150914 form?

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.

preprint2015arXiv

Modelling the high-energy emission from gamma-ray binaries using numerical relativistic hydrodynamics

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 $σ$ and Lorentz factor $Γ_{p}$ of the pulsar wind. Our best compromise implies $σ\approx 1$ and $Γ_{p}\approx 5\times 10^{3}$, respectively higher and lower than the typical values in pulsar wind nebulae. The high value of $σ$ 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.

preprint2015arXiv

Patchy blazar heating: diversifying the thermal history of the intergalactic medium

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.

preprint2014arXiv

Gamma-ray binaries : a bridge between Be stars and high energy astrophysics

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.

preprint2014arXiv

The Effect of Nonlinear Landau Damping on Ultrarelativistic Beam Plasma Instabilities

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.

preprint2013arXiv

Simulating gamma-ray binaries with a relativistic extension of RAMSES

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.

preprint2012arXiv

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.

preprint2012arXiv

Impact of orbital motion on the structure and stability of adiabatic shocks in colliding wind binaries

The collision of winds from massive stars in binaries results in the formation of a double-shock structure with observed signatures from radio to X-rays. We study the structure and stability of the colliding wind region as it turns into a spiral due to orbital motion. We focus on adiabatic winds, where mixing between the two winds is expected to be restricted to the Kelvin-Helmholtz instability (KHI). Mixing of the Wolf-Rayet wind with hydrogen-rich material is important for dust formation in pinwheel nebulae such as WR 104, where the spiral structure has been resolved in infrared. We use the hydrodynamical code RAMSES with an adaptive grid. A wide range of binary systems with different wind velocities and mass loss rates are studied with 2D simulations. A specific 3D simulation is performed to model WR 104. Orbital motion leads to the formation of two distinct spiral arms where the KHI develops differently. We find that the spiral structure is destroyed when there is a large velocity gradient between the winds, unless the collimated wind is much faster. We argue that the KHI plays a major role in maintaining or not the structure. We discuss the consequences for various colliding wind binaries. When stable, there is no straightforward relationship between the spatial step of the spiral, the wind velocities, and the orbital period. Our 3D simulation of WR 104 indicates that the colder, well-mixed trailing arm has more favourable conditions for dust formation than the leading arm. The single-arm infrared spiral follows more closely the mixing map than the density map, suggesting the dust-to-gas ratio may vary between the leading and trailing density spirals. However, the density is much lower than what dust formation models require. Including radiative cooling would lead to higher densities, and also to thin shell instabilities whose impact on the large structure remains unknown.

preprint2012arXiv

Modelling colliding wind binaries with RAMSES, extension to special relativity

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.

preprint2011arXiv

High resolution numerical simulations of unstable colliding stellar winds

We investigate the hydrodynamics of the interaction of two supersonic winds in binary systems. The collision of the winds creates two shocks separated by a contact discontinuity. The overall structure depends on the momentum flux ratio eta of the winds. We use the code RAMSES with adaptive mesh refinement to study the shock structure up to smaller values of eta, higher spatial resolution and greater spatial scales than have been previously achieved. 2D and 3D simulations, neglecting orbital motion, are compared to widely-used analytic results and their applicability is discussed. In the adiabatic limit, velocity shear at the contact discontinuity triggers the Kelvin-Helmholtz instability. We quantify the amplitude of the resulting fluctuations and find that they can be significant even with a modest initial shear. Using an isothermal equation of state leads to the development of thin shell instabilities. The initial evolution and growth rates enables us to formally identify the non-linear thin shell instability (NTSI) close to the binary axis. Some analogue of the transverse acceleration instability is present further away. The NTSI produces large amplitude fluctuations and dominates the long-term behaviour. We point out the computational cost of properly following these instabilities. Our study provides a basic framework to which the results of more complex simulations, including additional physical effects, can be compared.