Galaxy formation and evolution

Abstract:

Understanding the processes that drive the formation of black holes (BHs) is a key topic in observational cosmology. While the observed $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ correlation in bulge-dominated galaxies is thought to be produced by major mergers, the existence of a $M_{\mathrm{BH}}$--$M_{\star}$ relation, across all galaxy morphological types, suggests that BHs may be largely built by secular processes. Recent evidence that bulge-less galaxies, which are unlikely to have had significant mergers, are offset from the $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ relation, but lie on the $M_{\mathrm{BH}}$--$M_{\star}$ relation, has strengthened this hypothesis. Nevertheless, the small size and heterogeneity of current datasets, coupled with the difficulty in measuring precise BH masses, makes it challenging to address this issue using empirical studies alone. Here, we use Horizon-AGN, a cosmological hydrodynamical simulation to probe the role of mergers in BH growth over cosmic time. We show that (1) as suggested by observations, simulated bulge-less galaxies lie offset from the main $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ relation, but on the $M_{\mathrm{BH}}$--$M_{\star}$ relation, (2) the positions of galaxies on the $M_{\mathrm{BH}}$--$M_{\star}$ relation are not affected by their merger histories and (3) only $\sim$35 per cent of the BH mass in today's massive galaxies is directly attributable to merging -- the majority ($\sim$65 per cent) of BH growth, therefore, takes place gradually, via secular processes, over cosmic time.

Abstract:

Observations find a median star formation efficiency per free-fall time in Milky Way Giant Molecular Clouds (GMCs) on the order of $\epsilon_{\rm ff}\sim 1\%$ and a four order of magnitude spread in values ($0.01\%-100\%$). The origin of the large range in $\epsilon_{\rm ff}$ is still debated and difficult to reproduce with analytical models. We track the formation, evolution and destruction of GMCs in a hydrodynamical simulation of a Milky Way-like galaxy and by deriving cloud properties in an observationally motivated way, measure the distribution of star formation efficiencies which are in excellent agreement with observations. We find no significant link between $\epsilon_{\rm ff}$ and any measured global property of GMCs (e.g. gas mass, velocity dispersion). Instead, a wide range of efficiencies exist in the entire parameter space. From the cloud evolutionary tracks, we find that each cloud follow a \emph{unique} evolutionary path which gives rise to wide diversity in all properties. We argue that it is this diversity in cloud properties, above all else, that results in scatter of $\epsilon_{\rm ff}$.

Abstract:

The extreme infrared (IR) luminosity of local luminous and ultra-luminous IR galaxies (U/LIRGs; 11 < log LIR /Lsun < 12 and log LIR /Lsun > 12, respectively) is mainly powered by star-formation processes triggered by mergers or interactions. While U/LIRGs are rare locally, at z > 1, they become more common, they dominate the star-formation rate (SFR) density, and a fraction of them are found to be normal disk galaxies. Therefore, there must be an evolution of the mechanism triggering these intense starbursts with redshift. To investigate this evolution, we present new optical SWIFT integral field spectroscopic H{\alpha}+[NII] observations of a sample of 9 intermediate-z (0.2 < z < 0.4) U/LIRG systems selected from Herschel 250{\mu}m observations. The main results are the following: (a) the ratios between the velocity dispersion and the rotation curve amplitude indicate that 10-25% (1-2 out of 8) might be compatible with being isolated disks while the remaining objects are interacting/merging systems; (b) the ratio between un-obscured and obscured SFR traced by H{\alpha} and LIR, respectively, is similar in both local and these intermediate-z U/LIRGs; and (c) the ratio between 250{\mu}m and the total IR luminosities of these intermediate-z U/LIRGs is higher than that of local U/LIRGs with the same LIR . This indicates a reduced dust temperature in these intermediate-z U/LIRGs. This, together with their already measured enhanced molecular gas content, suggests that the interstellar medium conditions are different in our sample of intermediate-z galaxies when compared to local U/LIRGs.

Abstract:

We investigate the evolution of triple, non-hierarchical, black hole (BH) systems making use of $2.9\times10^4$ 3-body simulations. Varying the mutual orbital inclination, the three BH masses and the inner and outer eccentricities, we show that retrograde, nearly planar configurations lead to a significant shrinkage of the inner binary. We find an universal trend of triple systems, that they tend to evolve toward prograde configurations, Moreover, we demonstrate that the orbital flip, driven by the torque exerted on the inner BH binary (BHB) by the outer BH, leads in general to tighter inner orbits. In some cases, the resulting BHB undergoes coalescence within a Hubble time, releasing gravitational waves (GWs). Frequently, the inner BHB merger occurs after a component swap between one of its components and the outer BH. The mass spectrum of the BHBs that underwent the component exchange differs significantly from the case in which the BHB merge without any swap. A large fraction of merging BHBs with initial separation $1$ AU enter the $10^{-3}-10^{-1}$ Hz frequency band with large eccentricities, thus representing potential LISA sources. Mergers originating from initially tighter BHB ($a\sim 0.01$ AU), instead, have a large probability to have eccentricities above 0.7 in the $1$ Hz band. We find that the mergers' mass distribution in this astrophysical channel maps the original BH binary spectrum. This might have interesting consequences in light of the growing population of BH mergers detected by LIGO.

Abstract:

Infrared and nebular lines provide some of our best probes of the physics regulating the properties of the interstellar medium (ISM) at high-redshift. However, interpreting the physical conditions of high-redshift galaxies directly from emission lines remains complicated due to inhomogeneities in temperature, density, metallicity, ionisation parameter, and spectral hardness. We present a new suite of cosmological, radiation-hydrodynamics simulations, each centred on a massive Lyman-break galaxy that resolves such properties in an inhomogeneous ISM. Many of the simulated systems exhibit transient but well defined gaseous disks that appear as velocity gradients in [CII]~158.6$\mu$m emission. Spatial and spectral offsets between [CII]~158.6$\mu$m and [OIII]~88.33$\mu$m are common, but not ubiquitous, as each line probes a different phase of the ISM. These systems fall on the local [CII]-SFR relation, consistent with newer observations that question previously observed [CII]~158.6$\mu$m deficits. Our galaxies are consistent with the nebular line properties of observed $z\sim2-3$ galaxies and reproduce offsets on the BPT and mass-excitation diagrams compared to local galaxies due to higher star formation rate (SFR), excitation, and specific-SFR, as well as harder spectra from young, metal-poor binaries. We predict that local calibrations between H$\alpha$ and [OII]~3727$\AA$ luminosity and galaxy SFR apply up to $z>10$, as do the local relations between certain strong line diagnostics (R23 and [OIII]~5007$\AA$/H$\beta$) and galaxy metallicity. Our new simulations are well suited to interpret the observations of line emission from current (ALMA and HST) and upcoming facilities (JWST and ngVLA).