Semi-analytical black hole growth using cosmological simulations
I am currently studying black hole growth in Milky-Way sized galaxies using dark-matter only
cosmological simulations. Here is a gif animation of the entire simulation, the circles show the halos
and the 'plus' symbols show the locations of the black hole particles. There are about 1400 black holes
in the box.
Hot Halo Gas in Galaxy Mergers (Thesis)
I work on high-resolution numerical simulations of colliding
galaxies with a baryonic halo. For this, I have first initialized a dark matter
NFW halo with an equilibrium gas fraction. The dark matter
position and velocities are obtained from a numerically calculated
distribution function for the joint dark matter and gas
distribution.
I studied the mergers of galaxies for a range of
mass-ratios, gas fractions and impact parameters (for a parabolic orbit); a total
of 36 simulations. Select results are directly shown in the following images.
The following figure shows a strong shock being formed
during the merger (Left panel) and the evolution of unbound + hot and unbound
gas fraction in the right panel. The figure is from an equal-mass merger with an impact
parameter of ~20 kpc.
Gif Animation
This is from a parabolic equal-mass merger with an impact parameter of ~110 kpc. The
contours show the projected X-ray luminosity while the arrows show the velocities
of a random fraction of the unbound particles (identified at the end of the simulation).
The following figure shows the shocks in the simulation.
The strongest shocks are created when the material flung out as the nuclei make a close pass interacts
with pre-heated material from the initial shock.
Gif Animation
Movie of fractional histogram of location of unbound particles in each of the galaxies. The numbers
in the legend indicate the fraction of unbound material within the range of the plot. This is from
the runs testing for the dependence of the results on artificial viscosity. Essentially there
were no differences for values of artificial viscosity between 0.1 and 1.0.
Movie
Can the IGM be fueled by mergers ? (Thesis)
I found that every merger resulted in some fraction of the halo gas getting unbound. Now,
mergers are quite frequent and play a critical role in the formation of structure in
the Universe. As such, if every mergers release some fraction of the gas, then the
cumulative amount could be significant. To estimate this, I use an semi-analytic
Extended Press-Schechter merger tree (using the public GALFORM code from the Durham
group,
paper) and track the gas fractions
of each halo accurately. All halos that appear for the first time are given the universal
gas fraction (all halos at the highest z, 6 is what I used qualify). Afterwards, I
track the unbound gas fractions consistently through the entire merger tree. Convolving
with the evolving number density of the halos, I get the fractional mass in the Universe
contributed by mergers of galaxies.
The 3:1 line represents the mass-fraction of the unbound gas when all mergers bigger than
3:1 release 10% of the gas in the incident halos. The 10:1 is the same but when considering
all mergers bigger than 10:1 release 10% of the gas. Since there are many more 10:1 mergers
in the Universe than 3:1, the 10:1 line is always greater than the 3:1. The shaded region
shows the 1-sigma deviation from the mean. Thus, ~35% of the gas can be released just
by mergers under somewhat conservative estimates.
Now, my simulations do not include any AGN feedback or starburst-driven galactic-scale winds.
Those will work in addition to the unbinding caused by the merger.
