Normalized to: Lamberts, A.
[1]
oai:arXiv.org:2006.13956 [pdf] - 2122111
CASTLE: performances and science cases
Lombardo, S.;
Prada, F.;
Hugot, E.;
Basa, S.;
Bautista, J. M.;
Boissier, S.;
Boselli, A.;
Bosma, A.;
Cuillandre, J. C.;
Duc, P. A.;
Ferrari, M.;
Grosso, N.;
Izzo, L.;
Joaquina, K.;
Junais;
Koda, J.;
Lamberts, A.;
Lemaitre, G. R.;
Longobardi, A.;
Martínez-Delgado, D.;
Muslimov, E.;
Ortiz, J. L.;
Perez, E.;
Porquet, D.;
Sicardy, B.;
Vola, P.
Submitted: 2020-06-24
We present here the Calar Alto Schmidt-Lemaitre Telescope (CASTLE) concept, a
technology demonstrator for curved detectors, that will be installed at the
Calar Alto Observatory (Spain). This telescope has a wide field of view
(2.36x1.56 deg^2) and a design, optimised to generate a Point Spread Function
with very low level wings and reduced ghost features, which makes it
considerably less susceptible to several systematic effects usually affecting
similar systems. These characteristics are particularly suited to study the low
surface brightness Universe. CASTLE will be able to reach surface brightness
orders of magnitude fainter than the sky background level and observe the
extremely extended and faint features around galaxies such as tidal features,
stellar halos, intra-cluster light, etc. CASTLE will also be used to search and
detect astrophysical transients such as gamma ray bursts (GRB), gravitational
wave optical counterparts, neutrino counterparts, etc. This will increase the
number of precisely localized GRBs from 20% to 60% (in the case of Fermi/GMB
GRBs).
[3]
oai:arXiv.org:2006.08695 [pdf] - 2115454
Revisiting the Impact of Dust Production from Carbon-Rich Wolf-Rayet
Binaries
Submitted: 2020-06-15
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).
[4]
oai:arXiv.org:2006.07382 [pdf] - 2114303
Stars stripped in binaries -- the living gravitational wave sources
Submitted: 2020-06-12
Binary interaction can cause stellar envelopes to be stripped, which
significantly reduces the radius of the star. The orbit of a binary composed of
a stripped star and a compact object can therefore be so tight that the
gravitational radiation the system produces reaches frequencies accessible to
the Laser Interferometer Space Antenna (LISA). Two such stripped stars in tight
orbits with white dwarfs are known so far (ZTF J2130+4420 and CD-30 11223), but
many more are expected to exist. These binaries provide important constraints
for binary evolution models and may be used as LISA verification sources. We
develop a Monte Carlo code that uses detailed evolutionary models to simulate
the Galactic population of stripped stars in tight orbits with either neutron
star or white dwarf companions. We predict 0-100 stripped star + white dwarf
binaries and 0-4 stripped star + neutron star binaries with SNR>5 after 10
years of observations with LISA. More than 90% of these binaries are expected
to show large radial velocity shifts of $\gtrsim 200$ km/s, which are
spectroscopically detectable. Photometric variability due to tidal deformation
of the stripped star is also expected and has been observed in ZTF J2130+4420
and CD-30 11223. In addition, the stripped star + neutron star binaries are
expected to be X-ray bright with $L_X \gtrsim 10^{33} - 10^{36}$ erg/s. Our
results show that stripped star binaries are promising multi-messenger sources
for the upcoming electromagnetic and gravitational wave facilities.
[5]
oai:arXiv.org:2005.04239 [pdf] - 2115328
GW190412 as a Third-Generation Black Hole Merger from a Super Star
Cluster
Rodriguez, Carl L.;
Kremer, Kyle;
Grudić, Michael Y.;
Hafen, Zachary;
Chatterjee, Sourav;
Fragione, Giacomo;
Lamberts, Astrid;
Martinez, Miguel A. S.;
Rasio, Frederic A.;
Weatherford, Newlin;
Ye, Claire S.
Submitted: 2020-05-08, last modified: 2020-05-27
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.
[7]
oai:arXiv.org:2003.02849 [pdf] - 2077023
The growth of the longitudinal beam-plasma instability in the presence
of an inhomogeneous background
Submitted: 2020-03-05
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.
[9]
oai:arXiv.org:1702.02586 [pdf] - 2048492
Constraints on the Intergalactic Magnetic Field from Bow Ties in the
Gamma-ray Sky
Submitted: 2017-02-08, last modified: 2020-02-13
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.
[11]
oai:arXiv.org:1912.09992 [pdf] - 2018333
Astro2020 APC White Paper: Theoretical Astrophysics 2020-2030
Kollmeier, Juna A.;
Anderson, Lauren;
Benson, Andrew;
Bogdanovic, Tamara;
Boylan-Kolchin, Michael;
Bullock, James S.;
Dave, Romeel;
Fraschetti, Federico;
Fuller, Jim;
Hopkins, Philip F.;
Kaplinghat, Manoj;
Kratter, Kaitlin;
Lamberts, Astrid;
Miller, M. Coleman;
Owen, James E.;
Phinney, E. Sterl;
Piro, Anthony L.;
Rix, Hans-Walter;
Robertson, Brant;
Wetzel, Andrew;
Wheeler, Coral;
Youdin, Andrew N.;
Zaldarriaga, Matias
Submitted: 2019-12-18
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.
[12]
oai:arXiv.org:1912.07627 [pdf] - 2065342
Finding binary black holes in the Milky Way with LISA
Submitted: 2019-12-16
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.
[13]
oai:arXiv.org:1908.11391 [pdf] - 1951581
Unveiling the Gravitational Universe at \mu-Hz Frequencies
Sesana, Alberto;
Korsakova, Natalia;
Sedda, Manuel Arca;
Baibhav, Vishal;
Barausse, Enrico;
Barke, Simon;
Berti, Emanuele;
Bonetti, Matteo;
Capelo, Pedro R.;
Caprini, Chiara;
Garcia-Bellido, Juan;
Haiman, Zoltan;
Jani, Karan;
Jennrich, Oliver;
Johansson, Peter;
Khan, Fazeel Mahmood;
Korol, Valeriya;
Lamberts, Astrid;
Lupi, Alessandro;
Mangiagli, Alberto;
Mayer, Lucio;
Nardini, Germano;
Pacucci, Fabio;
Petiteau, Antoine;
Raccanelli, Alvise;
Rajendran, Surjeet;
Regan, John;
Shao, Lijing;
Spallicci, Alessandro;
Tamanini, Nicola;
Volonteri, Marta;
Warburton, Niels;
Wong, Kaze;
Zumalacarregui, Miguel
Submitted: 2019-08-29
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.
[14]
oai:arXiv.org:1907.00014 [pdf] - 1979621
Predicting the LISA white dwarf binary population in the Milky Way with
cosmological simulations
Submitted: 2019-06-28
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.
[15]
oai:arXiv.org:1808.02959 [pdf] - 1790902
Missing Gamma-ray Halos and the Need for New Physics in the Gamma-ray
Sky
Submitted: 2018-08-08, last modified: 2018-10-26
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.
[16]
oai:arXiv.org:1708.09065 [pdf] - 1741992
From the Top Down and Back Up Again: Star Cluster Structure from
Hierarchical Star Formation
Submitted: 2017-08-29, last modified: 2018-08-30
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.
[17]
oai:arXiv.org:1801.03099 [pdf] - 1728847
Predicting the binary black hole population of the Milky Way with
cosmological simulations
Lamberts, Astrid;
Garrison-Kimmel, Shea;
Hopkins, Philip;
Quataert, Eliot;
Bullock, James;
Faucher-Giguère, Claude-André;
Wetzel, Andrew;
Keres, Dusan;
Drango, Kaliden;
Sanderson, Robyn
Submitted: 2018-01-09, last modified: 2018-07-20
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.
[18]
oai:arXiv.org:1806.08525 [pdf] - 1771573
The SPHERE view of Wolf-Rayet 104
Soulain, A.;
Millour, F.;
Lopez, B.;
Matter, A.;
Lagadec, E.;
Carbillet, M.;
Camera, A.;
Lamberts, A.;
Langlois, M;
Milli, J;
Avenhaus, H;
Magnard, Y;
Roux, A;
Moulin, T;
Carle, M;
Sevin, A;
Martinez, P;
Abe, L;
Ramos, J
Submitted: 2018-06-22
Context. WR104 is an emblematic dusty Wolf-Rayet star and the prototypical
member of a subgroup hosting spirals that are mainly observable with
high-angular resolution techniques. Previous aperture masking observations
showed that WR104 is likely an interacting binary star at the end of its life.
However, several aspects of the system are still unknown. This includes the
opening angle of the spiral, the dust formation locus, and the link between the
central binary star and a candidate companion star detected with the Hubble
Space Telescope (HST) at 1". Aims. Our aim was to directly image the dusty
spiral or "pinwheel" structure around WR104 for the first time and determine
its physical properties at large spatial scales. We also wanted to address the
characteristics of the candidate companion detected by the HST. Methods. For
this purpose, we used SPHERE and VISIR at the Very Large Telescope to
respectively image the system in the near-and mid-infrared. Both instruments
furnished an excellent view of the system at the highest angular resolution a
single, ground-based telescope can provide. Based on these direct images, we
then used analytical and radiative transfer models to determine several
physical properties of the system. Results. Employing a different technique
than previously used, our new images have allowed us to confirm the presence of
the dust pinwheel around the central star. We have also detected up to 5
revolutions of the spiral pattern of WR104 in the K-band for the first time.
The circumstellar dust extends up to 2 arcsec from the central binary star in
the N-band, corresponding to the past 20 years of mass loss. Moreover, we found
no clear evidence of a shadow of the first spiral coil onto the subsequent
ones, which likely points to a dusty environment less massive than inferred in
previous studies. We have also confirmed that the stellar candidate companion
previously detected by the HST is gravitationally bound to WR104 and herein
provide information about its nature and orbital elements.
[19]
oai:arXiv.org:1806.04143 [pdf] - 1886395
The Local Group on FIRE: Dwarf galaxy populations across a suite of
hydrodynamic simulations
Garrison-Kimmel, Shea;
Hopkins, Philip F.;
Wetzel, Andrew;
Bullock, James S.;
Boylan-Kolchin, Michael;
Keres, Dusan;
Faucher-Giguere, Claude-Andre;
El-Badry, Kareem;
Lamberts, Astrid;
Quataert, Eliot;
Sanderson, Robyn
Submitted: 2018-06-11
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.
[20]
oai:arXiv.org:1804.05071 [pdf] - 1692020
Growth of beam-plasma instabilities in the presence of background
inhomogeneity
Submitted: 2018-04-13
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.
[21]
oai:arXiv.org:1801.03121 [pdf] - 1851510
Merging massive black holes: the right place and the right time
Submitted: 2018-01-09
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.
[22]
oai:arXiv.org:1702.04362 [pdf] - 1608341
Relativistic simulations of long-lived reverse shocks in stratified
ejecta: the origin of flares in GRB afterglows
Submitted: 2017-02-14, last modified: 2017-12-05
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.
[23]
oai:arXiv.org:1704.00014 [pdf] - 1582282
Importance of resolving the spectral support of beam-plasma
instabilities in simulations
Submitted: 2017-03-31, last modified: 2017-10-04
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.
[24]
oai:arXiv.org:1702.04732 [pdf] - 1581514
SHARP: A Spatially Higher-order, Relativistic Particle-in-Cell Code
Submitted: 2017-02-15, last modified: 2017-06-11
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.
[25]
oai:arXiv.org:1703.05184 [pdf] - 1547841
Modelling Jets, Tori and Flares in Pulsar Wind Nebulae
Submitted: 2017-03-15
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.
[26]
oai:arXiv.org:1701.01124 [pdf] - 1570973
Numerical simulations and infrared spectro-interferometry reveal the
wind collision region in gamma2 Velorum
Lamberts, A.;
Millour, F.;
Liermann, A.;
Dessart, L.;
Driebe, T.;
Duvert, G.;
Finsterle, W.;
Girault, V.;
Massi, F.;
Petrov, R. G.;
Schmutz, W.;
Weigelt, G.;
Chesneau, O.
Submitted: 2017-01-04, last modified: 2017-03-09
Colliding stellar winds in massive binary systems have been studied through
their radio, optical lines and strong X-ray emission for decades. More
recently, near-infrared spectrointerferometric observations have become
available in a few systems, but isolating the contribution from the individual
stars and the wind collision region still remains a challenge. In this paper,
we study the colliding wind binary gamma2 Velorum and aim at identifying the
wind collision zone from infrared interferometric data, which provide unique
spatial information to determine the wind properties. Our analysis is based on
multi-epoch VLTI/AMBER data that allows us to separate the spectral components
of both stars. First, we determine the astrometric solution of the binary and
confirm previous distance measurements. We then analyse the spectra of the
individual stars, showing that the O star spectrum is peculiar within its
class. Then, we perform three-dimensional hydrodynamic simulations of the
system from which we extract model images, visibility curves and closure phases
which can be directly compared with the observed data. The hydrodynamic
simulations reveal the 3D spiral structure of the wind collision region, which
results in phase-dependent emission maps. Our model visibility curves and
closure phases provide a good match when the wind collision region accounts for
3 to 10 per cent of the total flux in the near infrared. The dialogue between
hydrodynamic simulations, radiative transfer models and observations allows us
to fully exploit the observations. Similar efforts will be crucial to study
circumstellar environments with the new generation of VLTI instruments like
GRAVITY and MATISSE.
[27]
oai:arXiv.org:1702.02585 [pdf] - 1602511
Bow Ties in the Sky II: Searching for Gamma-ray Halos in the Fermi Sky
Using Anisotropy
Submitted: 2017-02-08
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.
[28]
oai:arXiv.org:1610.02040 [pdf] - 1531930
The Linear Instability of Dilute Ultrarelativistic $e^{\pm}$ Pair Beams
Submitted: 2016-10-06
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.
[29]
oai:arXiv.org:1609.00387 [pdf] - 1528169
Bow Ties in the Sky I: The Angular Structure of Inverse Compton
Gamma-ray Halos in the Fermi Sky
Submitted: 2016-09-01
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.
[30]
oai:arXiv.org:1605.08783 [pdf] - 1475348
When and where did GW150914 form?
Submitted: 2016-05-27, last modified: 2016-07-28
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.
[31]
oai:arXiv.org:1502.07980 [pdf] - 1280798
Patchy blazar heating: diversifying the thermal history of the
intergalactic medium
Submitted: 2015-02-27, last modified: 2015-08-25
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.
[32]
oai:arXiv.org:1505.01026 [pdf] - 1269315
Modelling the high-energy emission from gamma-ray binaries using
numerical relativistic hydrodynamics
Submitted: 2015-05-05, last modified: 2015-05-07
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.
[33]
oai:arXiv.org:1410.3758 [pdf] - 882049
Gamma-ray binaries : a bridge between Be stars and high energy
astrophysics
Submitted: 2014-10-14
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.
[34]
oai:arXiv.org:1410.3797 [pdf] - 1222565
The Effect of Nonlinear Landau Damping on Ultrarelativistic Beam Plasma
Instabilities
Submitted: 2014-10-07, last modified: 2014-10-14
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.
[35]
oai:arXiv.org:1309.7629 [pdf] - 1179568
Simulating gamma-ray binaries with a relativistic extension of RAMSES
Submitted: 2013-09-29
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.
[36]
oai:arXiv.org:1212.4142 [pdf] - 1158558
Colliding wind binaries and gamma-ray binaries : relativistic version of
the RAMSES code
Submitted: 2012-12-17
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.
[37]
oai:arXiv.org:1212.4045 [pdf] - 603678
Modelling colliding wind binaries with RAMSES, extension to special
relativity
Submitted: 2012-12-17
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.
[38]
oai:arXiv.org:1202.2060 [pdf] - 1116515
Impact of orbital motion on the structure and stability of adiabatic
shocks in colliding wind binaries
Submitted: 2012-02-09
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.
[39]
oai:arXiv.org:1109.1434 [pdf] - 1083896
High resolution numerical simulations of unstable colliding stellar
winds
Submitted: 2011-09-07
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.
