The problem of the extent of habitable zones in different kinds of galaxies is one of the outstanding challenges for contemporary astrobiology. In the present study, we investigate habitability in a large sample of simulated galaxies from the Illustris Project in order to at least roughly quantify the hospitality to life of different galactic types. The pioneering study of Dayal et al. (2015) is critically examined and some of its results are amended. In particular, we find a tentative evidence for a second mode of galactic habitability comprising metal-rich dwarfs similar to IC 225, LMC or M32. The role of the galactic environment and the observation selection effects is briefly discussed and prospects for further research on the topic outlined.

We follow trajectories of recoiling supermassive black holes (SMBHs) in analytical and numerical models of galaxy merger remnants with masses of $10^{11} \rm{M_{sun}}$ and $10^{12} \rm{M_{sun}}$. We construct various merger remnant galaxies in order to investigate how the central SMBH mass and the mass ratio of progenitor galaxies influence escape velocities of recoiling SMBHs. Our results show that static analytical models of major merger remnant galaxies overestimate the SMBHs escape velocities. During major mergers violent relaxation leads to the decrease of galaxy mass and lower potential at large remnant radii. This process is not depicted in static analytical potential but clearly seen in our numerical models. Thus, the evolving numerical model is a more realistic description of dynamical processes in galaxies with merging SMBHs. We find that SMBH escape velocities in numerical major merger remnant galaxies can be up to 25 per cent lower compared to those in analytical models. Consequently, SMBHs in numerical models generally reach greater galactocentric distances and spend more time on bound orbits outside of the galactic nuclei. Thus, numerical models predict a greater number of spatially-offset active galactic nuclei (AGNs).

Primordial stars are likely to be very massive >30 Msun, form in isolation, and will likely leave black holes as remnants in the centers of their host dark matter halos. We expect primordial stars to form in halos in the mass range 10^6-10^10 Msun. Some of these early black holes, formed at redshifts z>10, could be the seed black hole for a significant fraction of the supermassive black holes found in galaxies in the local universe. If the black hole descendants of the primordial stars exist, their mergers with nearby supermassive black holes may be a prime candidate for long wavelength gravitational wave detectors. We simulate formation and evolution of dark matter halos in LambdaCDM universe. We seed high-redshift dark matter halos with early black holes, and explore the merger history of the host halos and the implications of black hole's kick velocities arising from their coalescence. The central concentration of low mass early black holes in present day galaxies is reduced if they experience even moderate kicks of tens of km/s. Even such modest kicks allow the black holes to leave their parent halo, which consequently leads to dynamical friction being less effective on the low mass black holes that were ejected, compared to those still embedded in their parent halos. Therefore, merger rates with central supermassive black holes in the largest halos may be reduced by more than an order of magnitude. Using analytical and illustrative cosmological N-body simulations, we quantify the role of kicks on the merger rates of black holes formed from massive metal free stars with supermassive black holes in present day galaxies.

We have identified two channels for the formation of compact dwarf galaxies in the Illustris simulation by reconstructing mass and distance histories of candidates located in the vicinity of the simulation's most massive cluster galaxies. One channel is tidal stripping of Milky Way mass galaxies that form outside of clusters and eventually sink into them, spiraling in toward central massive objects. Second channel of formation is an in-situ formation (in reference to the parent cluster) of dwarf mass galaxies, with negligible evolution and limited change in stellar mass. We find 19 compact dwarf galaxies at the centers of 14 clusters, consistent with observations. 30% of them have external origin while 70% are formed in-situ.

Recent observations of AGN activity in massive galaxies (log Mstar / Msun > 10.4) show that: 1) at z < 1, AGN-hosting galaxies do not show enhanced merger signatures compared to normal galaxies, 2) also at z < 1, most AGNs are hosted by quiescent galaxies; and 3) at z > 1, percentage of AGNs in star forming galaxies increases and becomes comparable to AGN percentage in quiescent galaxies at z ~ 2. How can major mergers explain AGN activity in massive quiescent galaxies which have no merger features and no star formation to indicate recent galaxy merger? By matching merger events in a cosmological N-body simulation to the observed AGN incidence probability in the COSMOS survey, we show that major merger triggered AGN activity is consistent with the observations. By distinguishing between "peak" AGNs (recently merger triggered and hosted by star forming galaxies) and "faded" AGNs (merger triggered a long time ago and now residing in quiescent galaxies), we show that the AGN occupation fraction in star forming and quiescent galaxies simply follows the evolution of the galaxy merger rate. Since the galaxy merger rate drops dramatically at z < 1, the only AGNs left to be observed are the ones triggered by old mergers and are now in the declining phase of their nuclear activity, hosted by quiescent galaxies. As we go toward higher redshifts the galaxy merger rate increases and the percentages of "peak" AGNs and "faded" AGNs become comparable.

We present an isolated Milky Way-like simulation in GADGET2 N-body SPH code. The Galactic disk star formation rate (SFR) surface densities and stellar mass indicative of Solar neighbourhood are used as thresholds to model the distribution of stellar mass in life friendly environments. SFR and stellar component density are calculated averaging the GADGET2 particle properties on a 2D grid mapped on the Galactic plane. The peak values for possibly habitable stellar mass surface density move from $10$ to $15$ kpc cylindrical galactocentric distance in $10$ Gyr simulated time span. At $10$ Gyr the simulation results imply the following. Stellar particles which have spent almost all of their life time in habitable friendly conditions reside typically at $\sim16$ kpc from Galactic centre and are $\sim 3$ Gyr old. Stellar particles that have spent $\ge 90 \%$ of their $4-5$ Gyr long life time in habitable friendly conditions, are also predominantly found in the outskirts of the Galactic disk. Less then $1 \%$ of these particles can be found at a typical Solar system galactocentric distance of $8-10$ kpc. Our results imply that the evolution of an isolated spiral galaxy is likely to result in galactic civilizations emerging at the outskirts of the galactic disk around stellar hosts younger than the Sun.

In this work we model the observed evolution in comoving number density of Lyman-alpha blobs (LABs) as a function of redshift, and try to find which mechanism of emission is dominant in LAB. Our model calculates LAB emission both from cooling radiation from the intergalactic gas accreting onto galaxies and from star formation (SF). We have used dark matter (DM) cosmological simulation to which we applied empirical recipes for Ly$\alpha$ emission produced by cooling radiation and SF in every halo. In difference to the previous work, the simulated volume in the DM simulation is large enough to produce an average LABs number density. At a range of redshifts $z\sim 1-7$ we compare our results with the observed luminosity functions of LABs and LAEs. Our cooling radiation luminosities appeared to be too small to explain LAB luminosities at all redshifts. In contrast, for SF we obtained a good agreement with observed LFs at all redshifts studied. We also discuss uncertainties which could influence the obtained results, and how LAB LFs could be related to each other in fields with different density.

We make black hole (BH) merger trees from Millennium and Millennium-II simulations to find under what conditions 10^9Msun SMBH can form by redshift z=7. In order to exploit both: large box size in the Millennium simulation; and large mass resolution in the Millennium-II simulation, we develop a method to combine these two simulations together, and use the Millennium-II merger trees to predict the BH seeds to be used in the Millennium merger trees. We run multiple semi-analytical simulations where SMBHs grow through mergers and episodes of gas accretion triggered by major mergers. As a constraint, we use observed BH mass function at redshift z=6. We find that in the light of the recent observations of moderate super-Eddington accretion, low-mass seeds (100Msun) could be the progenitors of high-redshift SMBHs (z~7), as long as the accretion during the accretion episodes is moderately super-Eddington, where f_Edd=3.7 is the effective Eddington ratio averaged over 50 Myr.

We use detailed numerical simulations of the coupled chemical, thermal and dynamical evolution of the gas in a turbulent molecular cloud to study the usefulness of the [CI] 609 micron and 370 micron fine structure emission lines as tracers of cloud structure. Emission from these lines is observed throughout molecular clouds, and yet the question of what we can learn from them about the physics of the clouds remains largely unexplored. We show that the fact that [CI] emission is widespread within molecular clouds is a simple consequence of the fact that the clouds are dominated by turbulent motions. Turbulence creates large density inhomogeneities, allowing radiation to penetrate deeply into the clouds. As a result, [CI] emitting gas is found throughout the cloud, rather than being concentrated at the edges. We examine how well we can use [CI] emission to trace the structure of the cloud, and show that the integrated intensity of the 609 micron line traces column density accurately over a wide range of visual extinctions. For extinctions greater than a few, [CI] and 13CO both perform well, but [CI] becomes a superior tracer of column densities for visual extinctions A_V <= 3 We have also studied the distribution of [CI] excitation temperatures in the gas, and show that these are typically smaller than the kinetic temperature, indicating that most of the carbon atoms are not in local thermodynamic equilibrium. We discuss how best to determine T_ex from observations of the [CI] lines, and how to use these values to estimate the column density of neutral atomic carbon. We show that even in the best case, we tend to systematically underestimate the atomic carbon content of the gas. Our results suggest that observationally-derived estimates of the atomic carbon content of real GMCs could be in error by as much as a factor of two.

We study the influence of the choice of cooling function on the formation of molecular clouds in high-resolution three-dimensional simulations of converging flows. We directly compare the results obtained using the simple, parametrized cooling function introduced by Koyama & Inutsuka (2002) and used by a number of converging flow studies with the results of the detailed calculation of the non-equilibrium chemistry and thermal balance of the gas. We find that a number of the cloud properties, such as the mass and volume filling fractions of cold gas, are relatively insensitive to the choice of cooling function. On the other hand, the cloud morphology and the large-scale velocity distribution of the gas do strongly depend on the cooling function. We show that the differences that we see can largely be explained by differences in the way that Lyman-alpha cooling is treated in the two complementary approaches, and that a proper non-equilibrium treatment of the ionisation and recombination of the gas is necessary in order to model the high-temperature cooling correctly. We also investigate the properties of the dense clumps formed within the cloud. In agreement with previous models, we find that the majority of these clumps are not self-gravitating, suggesting that some form of large-scale collapse of the cloud may be required in order to produce gravitationally unstable clumps and hence stars. Overall, the physical properties of the dense clumps are similar in both simulations, suggesting that they do not depend strongly on the choice of cooling function. However, we do find a systematic difference of around 10K in the mean temperatures of the clumps produced by the two models.

We present a comparative study of recent works on merger-timescales with dynamical friction and find a strong contrast between idealized/isolated mergers (Boylan-Kolchin et al. 2008) and mergers from a cosmological volume (Jiang et al. 2008). Our study measures the duration of mergers in a cosmological N-body simulation of dark matter, with emphasis on higher redshifts (z < 10) and a lower mass range. In our analysis we consider and compare two merger definitions; tidal disruption and coalescence. We find that the merger-time formula proposed by Jiang et al. (2008) describes our results well and conclude that cosmologically motivated merger-time formulae provide a more versatile and statistically robust approximation for practical applications such as semi-analytic/hybrid models.

We present results from high-resolution three-dimensional simulations of the turbulent interstellar medium that study the influence of the nature of the turbulence on the formation of molecular hydrogen. We have examined both solenoidal (divergence-free) and compressive (curl-free) turbulent driving, and show that compressive driving leads to faster H2 formation, owing to the higher peak densities produced in the gas. The difference in the H2 formation rate can be as much as an order of magnitude at early times, but declines at later times as the highest density regions become fully molecular and stop contributing to the total H2 formation rate. We have also used our results to test a simple prescription suggested by Gnedin et al. (2009) for modeling the influence of unresolved density fluctuations on the H2 formation rate in large-scale simulations of the ISM. We find that this approach works well when the H2 fraction is small, but breaks down once the highest density gas becomes fully molecular.

We explore the growth of < 10^7 Msun black holes that reside at the centers of spiral and field dwarf galaxies in a Local Group type of environment. We use merger trees from a cosmological N-body simulation known as Via Lactea II (VL-2) as a framework to test two merger-driven semi-analytic recipes for black hole growth that include dynamical friction, tidal stripping, and gravitational wave recoil in over 20,000 merger tree realizations. First, we apply a Fundamental Plane limited (FPL) model to the growth of Sgr A*, which drives the central black hole to a maximum mass limited by the Black Hole Fundamental Plane after every merger. Next, we present a new model that allows for low-level Prolonged Gas Accretion (PGA) during the merger. We find that both models can generate a Sgr A* mass black hole. We predict a population of massive black holes in local field dwarf galaxies - if the VL-2 simulation is representative of the growth of the Local Group, we predict up to 35 massive black holes (< 10^6 Msun) in Local Group field dwarfs. We also predict that hundreds of < 10^5 Msun black holes fail to merge, and instead populate the Milky Way halo, with the most massive of them at roughly the virial radius. In addition, we find that there may be hundreds of massive black holes ejected from their hosts into the nearby intergalactic medium due to gravitational wave recoil. We discuss how the black hole population in the Local Group field dwarfs may help to constrain the growth mechanism for Sgr A*.

We calculate the gravitational wave signal from the growth of 10 million solar mass supermassive black holes (SMBH) from the remnants of Population III stars. The assembly of these lower mass black holes is particularly important because observing SMBHs in this mass range is one of the primary science goals for the Laser Interferometer Space Antenna (LISA), a planned NASA/ESA mission to detect gravitational waves. We use high resolution cosmological N-body simulations to track the merger history of the host dark matter halos, and model the growth of the SMBHs with a semi-analytic approach that combines dynamical friction, gas accretion, and feedback. We find that the most common source in the LISA band from our volume consists of mergers between intermediate mass black holes and SMBHs at redshifts less than 2. This type of high mass ratio merger has not been widely considered in the gravitational wave community; detection and characterization of this signal will likely require a different technique than is used for SMBH mergers or extreme mass ratio inspirals. We find that the event rate of this new LISA source depends on prescriptions for gas accretion onto the black hole as well as an accurate model of the dynamics on a galaxy scale; our best estimate yields about 40 sources with a signal-to-noise ratio greater than 30 occur within a volume like the Local Group during SMBH assembly -- extrapolated over the volume of the universe yields roughly 500 observed events over 10 years, although the accuracy of this rate is affected by cosmic variance.

(Abridged) We use high resolution cosmological N-body simulations to study the growth of intermediate to supermassive black holes from redshift 49 to zero. We track the growth of black holes from the seeds of population III stars to black holes in the range of 10^3 < M < 10^7 Msun -- not quasars, but rather IMBH to low-mass SMBHs. These lower mass black holes are the primary observable for the Laser Interferometer Space Antenna (LISA). The large-scale dynamics of the black holes are followed accurately within the simulation down to scales of 1 kpc; thereafter, we follow the merger analytically from the last dynamical friction phase to black hole coalescence. We find that the merger rate of these black holes is R~25 per year between 8 < z < 11 and R = 10 per year at z=3. Before the merger occurs the incoming IMBH may be observed with a next generation of X-ray telescopes as a ULX source with a rate of about ~ 3 - 7 per year for 1 < z < 5. We develop an analytic prescription that captures the most important black hole growth mechanisms: galaxy merger-driven gas accretion and black hole coalescence. Using this, we find that SMBH at the center of Milky Way type galaxy was in place with most of its mass by z = 4.7, and most of the growth was driven by gas accretion excited by major mergers. Hundreds of black holes have failed to coalesce with the SMBH by z=0, some with masses of 10000 Msun, orbiting within the dark matter halo with luminosities up to ~ 30000 Lsun. These X-ray sources can easily be observed with Chandra at ~ 100 kpc.

Understanding how seed black holes grow into intermediate and supermassive black holes (IMBHs and SMBHs, respectively) has important implications for the duty-cycle of active galactic nuclei (AGN), galaxy evolution, and gravitational wave astronomy. Most studies of the cosmological growth and merger history of black holes have used semianalytic models and have concentrated on SMBH growth in luminous galaxies. Using high resolution cosmological N-body simulations, we track the assembly of black holes over a large range of final masses -- from seed black holes to SMBHs -- over widely varying dynamical histories. We used the dynamics of dark matter halos to track the evolution of seed black holes in three different gas accretion scenarios. We have found that growth of Sagittarius A* - size SMBH reaches its maximum mass M_{SMBH}~10^6Msun at z~6 through early gaseous accretion episodes, after which it stays at near constant mass. At the same redshift, the duty-cycle of the host AGN ends, hence redshift z=6 marks the transition from an AGN to a starburst galaxy which eventually becomes the Milky Way. By tracking black hole growth as a function of time and mass, we estimate that the IMBH merger rate reaches a maximum of R_{max}=55 yr^-1 at z=11. From IMBH merger rates we calculate N_{ULX}=7 per Milky Way type galaxy per redshift in redshift range 2<z<6.

Primordial stars are likely to be very massive $\geq30\Msun$, form in isolation, and will likely leave black holes as remnants in the centers of their host dark matter halos in the mass range $10^{6}-10^{10}\Ms$. Such early black holes, at redshifts z$\gtsim10$, could be the seed black holes for the many supermassive black holes found in galaxies in the local universe. If they exist, their mergers with nearby supermassive black holes may be a prime signal for long wavelength gravitational wave detectors. We simulate formation of black holes in the center of high redshift dark matter halos and explore implications of initial natal kick velocities conjectured by some formation models. The central concentration of early black holes in present day galaxies is reduced if they are born even with moderate kicks of tens of km/s. The modest kicks allow the black holes to leave their parent halo, which consequently leads to dynamical friction being less effective on the lower mass black holes as compared to those still embedded in their parent halos. Therefore, merger rates may be reduced by more than an order of magnitude. Using analytical and illustrative cosmological N--body simulations we quantify the role of natal kicks of black holes formed from massive metal free stars on their merger rates with supermassive black holes in present day galaxies. Our results also apply to black holes ejected by the gravitational slingshot mechanism.

Based on a high resolution cosmological n-body simulation, we track the hierarchical growth of black holes in galaxy clusters from z=20 to z=0. We present a census of black holes as function of redshift and will determine their mass assembly history under a variety of assumptions regarding the importance of gas accretion in black hole growth, from early supercritical Eddington accretion to gas-poor hierarchical assembly. Following a galaxy merger, black holes are expected to form, inspiral and merge after strongly radiating energy via gravitational waves. For each binary black hole inspiral and merger, we determine the expected gravitational wave signal for the Laser Interferometer Space Antenna (LISA), and calculate the LISA event rate as a function of time. We will calculate the black hole mass assembly history for several black hole growth scenerios, so that we can explore tests to characterize each model observationally. In particular, we will study how well LISA observations will be able to distinguish between these very different assembly scenarios.