Galaxy formation and evolution

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).

Abstract:

The presence of spectroscopically confirmed Balmer breaks in galaxy spectral energy distributions (SEDs) at $z>9$ provides one of the best probes of the assembly history of the first generations of stars in our Universe. Recent observations of the gravitationally lensed source, MACS 1149_JD1 (JD1), indicate that significant amounts of star formation likely occurred at redshifts as high as $z\simeq15$. The inferred stellar mass, dust mass, and assembly history of JD1, or any other galaxy at these redshifts that exhibits a strong Balmer break, can provide a strong test of our best theoretical models from high-resolution cosmological simulations. In this work, we present the results from a cosmological radiation-hydrodynamics simulation of the region surrounding a massive Lyman-break galaxy. For two of our most massive systems, we show that dust preferentially resides in the vicinity of the young stars thereby increasing the strength of the measured Balmer break such that the simulated SEDs are consistent with the photometry of JD1 and two other $z>9$ systems (GN-z10-3 and GN-z9-1) that have proposed Balmer breaks at high redshift. We find strong variations in the shape and luminosity of the SEDs of galaxies with nearly identical stellar and halo masses, indicating the importance of morphology, assembly history, and dust distribution in making inferences on the properties of individual galaxies at high redshifts. Our results stress the importance that dust may play in modulating the observable properties of galaxies, even at the extreme redshifts of $z>9$.