No Arabic abstract
We use hydrodynamical simulations to study the color transformations induced by star formation and active galactic nuclei (AGN) during major mergers of spiral galaxies. Our modeling accounts for radiative cooling, star formation, and supernova feedback. Moreover, we include a treatment of accretion onto supermassive black holes embedded in the nuclei of the merging galaxies. We assume that a small fraction of the bolometric luminosity of an accreting black hole couples thermally to surrounding gas, providing a feedback mechanism that regulates its growth. The encounter and coalescence of the galaxies triggers nuclear gas inflow which fuels both a powerful starburst and strong black hole accretion. Comparing simulations with and without black holes, we show that AGN feedback can quench star formation and accretion on a short timescale,particularly in large galaxies where the black holes can drive powerful winds once they become sufficiently massive. The color evolution of the remnant differs markedly between mergers with and without central black holes. Without AGN, gas-rich mergers lead to ellipticals which remain blue owing to residual star formation, even after more than 7 Gyrs have elapsed. In contrast, mergers with black holes produce ellipticals that redden much faster, an effect that is more pronounced in massive remnants where a nearly complete termination of star formation occurs, allowing them to redden to u-r ~ 2.3 in less than one Gyr. AGN feedback may thus be required to explain the population of extremely red massive early type-galaxies, and it appears to be an important driver in generating the observed bimodal color distribution of galaxies in the Local Universe.
This study explored the GALEX ultraviolet (UV) properties of optical red sequence galaxies in 4 rich Abell clusters at z leq 0.1. In particular, we tried to find a hint of merger-induced recent star formation (RSF) in red sequence galaxies. Using the NUV - r colors of the galaxies, RSF fractions were derived based on various criteria for post-merger galaxies and normal galaxies. Following k-correction, about 36% of the post-merger galaxies were classified as RSF galaxies with a conservative criterion (NUV - r leq 5), and that number was doubled (~ 72%) when using a generous criterion (NUV - r leq 5.4). The trend was the same when we restricted the sample to galaxies within 0.5xR_{200}. Post-merger galaxies with strong UV emission showed more violent, asymmetric features in the deep optical images. The RSF fractions did not show any trend along the clustocentric distance within R_{200}. We performed a Dressler-Shectman test to check whether the RSF galaxies had any correlation with the sub-structures in the galaxy clusters. Within R_{200} of each cluster, the RSF galaxies did not appear to be preferentially related to the clusters sub-structures. Our results suggested that only 30% of RSF red sequence galaxies show morphological hints of recent galaxy mergers. This implies that internal processes (e.g., stellar mass-loss or hot gas cooling) for the supply of cold gas to early-type galaxies may play a significant role in the residual star formation of early-type galaxies at a recent epoch.
Based on numerical simulations of gas-rich galaxy mergers, we discuss a model in which quasar activity is tied to the self-regulated growth of supermassive black holes in galaxies. Nuclear inflow of gas attending a galaxy collision triggers a starburst and feeds black hole growth, but for most of the duration of the starburst, the black hole is heavily obscured by surrounding gas and dust which limits the visibility of the quasar, especially at optical and UV wavelengths. Eventually, feedback energy from accretion heats the gas and expels it in a powerful wind, leaving a dead quasar. Between buried and dead phases there is a window during which the galaxy would be seen as a luminous quasar. Because the black hole mass, radiative output, and distribution of obscuring gas and dust all evolve strongly with time, the duration of this phase of observable quasar activity depends on both the waveband and imposed luminosity threshold. We determine the observed and intrinsic lifetimes as a function of luminosity and frequency, and calculate observable lifetimes ~10 Myr for bright quasars in the optical B-band, in good agreement with empirical estimates and much smaller than the black hole growth timescales ~100 Myr, naturally producing a substantial population of buried quasars. However, observed and intrinsic energy outputs converge in the IR and hard X-ray bands as attenuation becomes weaker and chances of observation greatly increase. We obtain the distribution of column densities along sightlines in which the quasar is seen above a given luminosity, and find that our result agrees remarkably well with observed estimates of the column density distribution from the SDSS for appropriate luminosity thresholds. (Abridged)
We examine the pairing process of supermassive black holes (SMBHs) down to scales of 20-100 pc using a set of N-body/SPH simulations of binary mergers of disk galaxies with mass ratios of 1:4 and 1:10. Our numerical experiments are designed to represent merger events occurring at various cosmic epochs. The initial conditions of the encounters are consistent with the LambdaCDM paradigm of structure formation, and the simulations include the effects of radiative cooling, star formation, and supernovae feedback. We find that the pairing of SMBHs depends sensitively on the amount of baryonic mass preserved in the center of the companion galaxies during the last phases of the merger. In particular, due to the combination of gasdynamics and star formation, we find that a pair of SMBHs can form in 1:10 minor mergers provided that galaxies are relatively gas-rich (gas fractions of 30% of the disk mass) and that the mergers occur at relatively high redshift (z~3), when dynamical friction timescales are shorter. Since 1:10 mergers are most common events during the assembly of galaxies, and mergers are more frequent at high redshift when galaxies are also more gas-rich, our results have positive implications for future gravitational wave experiments such as the Laser Interferometer Space Antenna.
We present novel 3D multi-scale SPH simulations of gas-rich galaxy mergers between the most massive galaxies at $z sim 8 - 10$, designed to scrutinize the direct collapse formation scenario for massive black hole seeds proposed in citet{mayer+10}. The simulations achieve a resolution of 0.1 pc, and include both metallicity-dependent optically-thin cooling and a model for thermal balance at high optical depth. We consider different formulations of the SPH hydrodynamical equations, including thermal and metal diffusion. When the two merging galaxy cores collide, gas infall produces a compact, optically thick nuclear disk with densities exceeding $10^{-10}$ g cm$^3$. The disk rapidly accretes higher angular momentum gas from its surroundings reaching $sim 5$ pc and a mass of $gtrsim 10^9$ $M_{odot}$ in only a few $10^4$ yr. Outside $gtrsim 2$ pc it fragments into massive clumps. Instead, supersonic turbulence prevents fragmentation in the inner parsec region, which remains warm ($sim 3000-6000$ K) and develops strong non-axisymmetric modes that cause prominent radial gas inflows ($> 10^4$ $M_{odot}$ yr$^{-1}$), forming an ultra-dense massive disky core. Angular momentum transport by non-axisymmetric modes should continue below our spatial resolution limit, quickly turning the disky core into a supermassive protostar which can collapse directly into a massive black hole of mass $10^8-10^9$ $M_{odot}$ via the relativistic radial instability. Such a cold direct collapse explains naturally the early emergence of high-z QSOs. Its telltale signature would be a burst of gravitational waves in the frequency range $10^{-4} - 10^{-1}$ Hz, possibly detectable by the planned eLISA interferometer.
The coalescence of massive black hole binaries (BHBs) in galactic mergers is the primary source of gravitational waves (GWs) at low frequencies. Current estimates of GW detection rates for the Laser Interferometer Space Antenna and the Pulsar Timing Array vary by three orders of magnitude. To understand this variation, we simulate the merger of equal-mass, eccentric, galaxy pairs with central massive black holes and shallow inner density cusps. We model the formation and hardening of a central BHB using the Fast Multiple Method as a force solver, which features a $O(N)$ scaling with the number $N$ of particles and obtains results equivalent to direct-summation simulations. At $N sim 5times 10^5$, typical for contemporary studies, the eccentricity of the BHBs can vary significantly for different random realisations of the same initial condition, resulting in a substantial variation of the merger timescale. This scatter owes to the stochasticity of stellar encounters with the BHB and decreases with increasing $N$. We estimate that $N sim 10^7$ within the stellar half-light radius suffices to reduce the scatter in the merger timescale to $sim 10$%. Our results suggest that at least some of the uncertainty in low-frequency GW rates owes to insufficient numerical resolution.