Search

We present an analysis of the physical and dynamical states of two sets of EAGLE zoom simulations of galaxy haloes, one at high redshift ($z=2-3$) and the other at low redshift ($z=0$), with masses of $\approx 10^{12} M_{\odot}$. Our focus is how the circumgalactic medium (CGM) of these $L^*$ star-forming galaxies change over the last 10 Gyr. We find that the high-$z$ CGM is almost equally divided between the "cool" ($T, 16 pages, 15 figures. Submitted to MNRAS. Comments welcome. (revision submitted on May 22, 2020 fixes Fig. 3. right panel that mistakenly displayed metal mass and now correctly displays total gas mass)

Early-type galaxies (ETGs) are of crucial importance to trace back the galaxy mass assembly across cosmic time, yet their formation and quenching remain remarkably elusive. The discoveries of massive, dead galaxies at ever-growing redshifts provided compelling evidence to push their formation up to redshift > 4–5 when the Universe was barely 1 Gyr old. In this talk I will present our results on the ages of a new sample of ETGs at z ∼ 3, built by exploiting HST WFC3/G141 rest-frame optical/near-UV grism spectroscopy to study the nature of 10 passive galaxy candidates at 2.5 < z < 3.5 in COSMOS.This work is part of a PhD project aimed at quantifying the parent space density of distant genuine passive galaxies. I will also discuss the importance of multi-wavelength data in clarifying the degree of contamination by dusty star-forming galaxies affecting the color selection.

This article presents an up-dated analysis of synthetic optical and UV emission lines of simulated galaxies over cosmic time. The strong emission lines are derived from self-consistently coupling novel spectral models accounting for nebular emission from young stars, AGN and Post-AGB stars to cosmological zoom-in as well as large-scale simulations. Investigating the evolution of optical line-ratios in the BPT diagrams, the simulations can successfully reproduce the observed trend of [OIII]/Hβ ratio increasing from low to high redshifts, due to evolving star formation rate and gas metallicity. Standard selection criteria in the BPT diagrams can appropriately distinguish the main ionising source(s) of galaxies at low redshifts, but they are less reliable for metal-poor galaxies, dominating the early Universe. To robustly classify the ionising radiation of such metal-poor galaxies, diagnostic diagrams based on luminosity ratios of UV lines are discussed. The novel interface between simulations and observations is potentially important for the interpretation of high-quality spectra of very distant galaxies to be gathered by next-generation telescopes, such as the James Webb Space Telescope.

Current observations suggest an accelerated evolution of the cosmic star formation rate density for 8 < z < 10, indicating that galaxy assembly experienced an extremely intense phase during the first ∼600Myr years of cosmic time. We performed a systematic search of ultrabright star-forming galaxies at z ∼ 8 over the COSMOS/UltraVISTA survey, identifying 16 candidate Lyman-break galaxies. The still large uncertainties on the associated volume density do not yet allow us to ascertain whether a different star-formation efficiency (SFE) existed at early cosmic epochs. Leveraging the deepest Spitzer/IRAC data available from the GREATS program over the CANDELS/GOODS fields, we also constructed stacked SEDs of sub- L* LBGs at z ∼ 8. We find extreme nebular line emission (EW0 (Hα) ∼ 1000Å), high specific star-formation rates (∼10/Gyr) and indication of an inverse Balmer break. These results point toward very young ages (

ISM comprises multiple components, including molecular, neutral, and ionized gas, and dust, which are related to each other mainly through star formation - some are fuel for star formation (molecular gas) while some are the products of it (ionized gas, dust). To fully understand the physics of star formation and its evolution throughout cosmic time, it is crucial to measure and observe different ISM components of galaxies out to high redshifts. I will review the current status of near-IR studies of galaxies during the peak of star formation activity (z~1-3). Using rest-frame optical emission lines, we measure dust, star formation, and gaseous properties of galaxies. JWST will advance such studies by probing lower luminosities and higher redshifts, owing to its significantly higher sensitivity. Incorporating ALMA observations of cold dust and molecular gas at z>1 will give us a nearly complete picture of the ISM in high-redshift galaxies over a large dynamic range in mass., To be published in the proceedings of the IAU Symposium 352 (12 pages, review article)

The next generation ground-based gravitational wave interferometers will possibly observe mergers of binary black holes (BBHs) and binary neutron stars (BNSs) to redshift $z\gtrsim{}10$ and $z\gtrsim{}2$, respectively. Here, we characterize the properties of merging BBHs, BNSs and neutron star-black hole binaries across cosmic time, by means of population-synthesis simulations combined with the Illustris cosmological simulation. We find that the mass of merging compact objects does not depend (or depends very mildly) on the merger redshift. Even the mass distribution of black holes depends only mildly on redshift, because BBHs originating from metal-poor progenitors ($Z\leq{}4\times{}10^{-3}$) dominate the entire population of merging BBHs across cosmic time. For a common-envelope efficiency $\alpha{}\ge{}3$, the main difference between the mass distribution of BBHs merging in the last Gyr and that of BBHs merging more than 11 Gyr ago is that there is an excess of heavy merging black holes ($20-35$ M$_\odot$) in the last Gyr. This excess is explained by the longer delay time of massive BBHs., Comment: 13 pages, 11 figures, 3 tables, MNRAS, in press. Fixed typo in Figure 9

Context. The properties of galaxies at redshift z > 6 hold the key to our understanding of the early stages of galaxy evolution and can potentially identify the sources of the ultraviolet radiation that give rise to the epoch of reionisation. The far-infrared cooling line of [C II] at 158 μm is known to be bright and correlate with the star formation rate (SFR) of low-redshift galaxies, and hence is also suggested to be an important tracer of star formation and interstellar medium properties for very high-redshift galaxies. Aims. With the aim to study the interstellar medium properties of gravitationally lensed galaxies at z > 6, we search for [C II] and thermal dust emission in a sample of 52 z ∼ 6 galaxies observed by the ALMA Lensing Cluster Survey. Methods. We perform our analysis using LINESTACKER, stacking both [C II] and continuum emission. The target sample is selected from multiple catalogues, and the sample galaxies have spectroscopic redshift or low-uncertainty photometric redshifts (σz Results. Our analyses find no detection of either [C II] or continuum. An upper limit on L[CII] is derived, implying that [C II] remains marginally consistent for low-SFR z > 6 galaxies but likely is under-luminous compared to the local L[CII]-SFR relationship. We discuss potential biases and possible physical effects that may be the cause of the non-detection. Further, the upper limit on the dust continuum implies that less than half of the star formation is obscured.

Post-starburst (or "E+A") galaxies trace the fastest and most dramatic processes in galaxy evolution. Recent work studying the evolution of galaxies through this phase have revealed insights on how galaxies undergo structural and stellar population changes as well as the role of various feedback mechanisms. In this review, I summarize recent work on identifying post-starburst galaxies; tracing the role of this phase through cosmic time; measuring stellar populations, on-going star formation, morphologies, kinematics, interstellar medium properties, and AGN activity; mechanisms to cause the recent starburst and its end; and the future evolution to quiescence (or not). The review concludes with a list of open questions and exciting possibilities for future facilities., 23 pages, 4 figures, Invited review for PASP, Accepted for publication in PASP

Fundamental differences in the radio properties of red quasars (QSOs), as compared to blue QSOs, have been recently discovered, positioning them as a potential key population in the evolution of galaxies and black holes across cosmic time. To elucidate their nature, we exploited a rich compilation of photometry and spectroscopic data to model their spectral energy distributions (SEDs) from the UV to the FIR and characterise their emission-line properties. Following a systematic comparison approach, we infer the AGN accretion, obscuration, and host galaxy properties in a sample of ~1800 QSOs at 0.21000 km/s) in red QSOs. We find that red QSOs that exhibit evidence for high-velocity winds present a stronger signature of the infrared excess, suggesting a causal connection between reddening and the presence of hot dust in QSO winds. We propose that dusty winds at nuclear scales are potentially the physical ingredient responsible for the colours in red QSOs, as well as a key parameter for the regulation of accretion material in the nucleus., Comment: Accepted for publication in Astronomy & Astrophysics. 22 pages, 16 figures. Abstract abridged

We present the first study of the large-scale clustering of post-starburst (PSB) galaxies in the high redshift Universe ($0.510^{14}$M$_{\odot}$), suggesting they are likely to be infalling satellite galaxies in cluster-like environments. High-mass PSB galaxies are more weakly clustered at low redshifts, but they reside in higher mass haloes with increasing look-back time, suggesting strong redshift-dependent halo downsizing. These key results are consistent with previous results suggesting that two main channels are responsible for the rapid quenching of galaxies. While high-redshift ($z>1$) galaxies appear to be quenched by secular feedback mechanisms, processes associated with dense environments are likely to be the key driver of rapid quenching in the low-redshift Universe ($z, Comment: 19 pages, 9 figures, 3 tables; accepted for publication in MNRAS. Comments/questions welcome: please send queries to aaron.wilkinson91@gmail.com

We analyse the metallicity histories of ~4,500 galaxies from the GAMA survey at z, 19 pages, 14 figures (inc. 2 animations), accepted for publication in MNRAS

Dwarf galaxies (M*, 22 pages, 19 figures, version accepted for publication in MNRAS. Some minor changes and additional appendices

While radio emission in quasars can be contributed to by a variety of processes (involving star forming regions, accretion disk coronas and winds, and jets), the powering of the radio loudest quasars must involve very strong jets, presumably launched by the Blandford-Znajek mechanism incorporating the magnetically arrested disk (MAD) scenario. We focus on the latter and investigate the dependence of their fraction on redshift. We also examine the dependence of the radio-loud fraction (RLF) on BH mass ($M_{\rm BH}$) and Eddington ratio ($��_{\rm Edd}$) while excluding the redshift bias by narrowing its range. In both these investigations we remove the bias associated with: (1) the diversity of source selection by constructing two well-defined, homogeneous samples of quasars (first within $0.7 \leq z < 1.9$, second within $0.5 \leq z < 0.7$); (2) a strong drop in the RLF of quasars at smaller BH masses by choosing those with BH masses larger than $10^{8.5} M_{\odot}$. We confirm some previous results showing the increase in the fraction of radio-loud quasars with cosmic time and that this trend can be even steeper if we account for the bias introduced by the dependence of the RLF on BH mass whereas the bias introduced by the dependence of the RLF on Eddington ratio is shown to be negligible. Assuming that quasar activities are triggered by galaxy mergers we argue that such an increase can result from the slower drop with cosmic time of mixed mergers than of wet mergers., 15 pages, 4 figures. Accepted for publication in ApJ

We present new stellar mass functions at z ∼ 6, z ∼ 7, z ∼ 8, z ∼ 9 and, for the first time, z ∼ 10, constructed from ∼800 Lyman-break galaxies previously identified over the eXtreme Deep Field and Hubble Ultra-Deep Field parallel fields and the five Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey fields. Our study is distinctive due to (1) the much deeper (∼200 hr) wide-area Spitzer/Infrared Array Camera (IRAC) imaging at 3.6 μm and 4.5 μm from the Great Observatories Origins Deep Survey Re-ionization Era Wide-area Treasury from Spitzer program (GREATS) and (2) consideration of z ∼ 6–10 sources over a 3× larger area than those of previous Hubble Space Telescope+Spitzer studies. The Spitzer/IRAC data enable ≥2σ rest-frame optical detections for an unprecedented 50% of galaxies down to a stellar mass limit of ∼ 10 8  ⊙ across all redshifts. Schechter fits to our volume densities suggest a combined evolution in the characteristic mass  * and normalization factor ϕ * between z ∼ 6 and z ∼ 8. The stellar mass density (SMD) increases by ∼1000× in the ∼500 Myr between z ∼ 10 and z ∼ 6, with indications of a steeper evolution between z ∼ 10 and z ∼ 8, similar to the previously reported trend of the star formation rate density. Strikingly, abundance matching to the Bolshoi–Planck simulation indicates halo mass densities evolving at approximately the same rate as the SMD between z ∼ 10 and z ∼ 4. Our results show that the stellar-to-halo mass ratios, a proxy for the star formation efficiency, do not change significantly over the huge stellar mass buildup occurred from z ∼ 10 to z ∼ 6, indicating that the assembly of stellar mass closely mirrors the buildup in halo mass in the first ∼1 Gyr of cosmic history. The James Webb Space Telescope is poised to extend these results into the “first galaxy” epoch at z ≳ 10.

We investigate the formation history of massive disk galaxies in hydro-dynamical simulation--the IllustrisTNG, to study why massive disk galaxies survive through cosmic time. 83 galaxies in the simulation are selected with M$_{*,z=0}$ $>8\times10^{10}$ M$_\odot$ and kinematic bulge-to-total ratio less than $0.3$. We find that 8.4 percent of these massive disk galaxies have quiet merger histories and preserve disk morphology since formed. 54.2 percent have a significant increase in bulge components in history, then become disks again till present time. The rest 37.3 percent experience prominent mergers but survive to remain disky. While mergers and even major mergers do not always turn disk galaxies into ellipticals, we study the relations between various properties of mergers and the morphology of merger remnants. We find a strong dependence of remnant morphology on the orbit type of major mergers. Specifically, major mergers with a spiral-in falling orbit mostly lead to disk-dominant remnants, and major mergers of head-on galaxy-galaxy collision mostly form ellipticals. This dependence of remnant morphology on orbit type is much stronger than the dependence on cold gas fraction or orbital configuration of merger system as previously studied., 12 pages, 1 table, 8 figures. Accepted for publication in MNRAS

Massive black holes (BHs) inhabiting galactic nuclei can be described by two parameters only, i.e. mass and spin, that change through cosmic time in response to accretion and merger events. While most numerical simulations accurately track the BH mass, spin evolution is rarely taken into account. In this work, we implement and validate a self-consistent sub-grid model for the evolution of the BH mass and spin via gas accretion in the hydrodynamics code GIZMO. The model assumes that accretion from resolved scales does not occur instantaneously, but is mediated by a sub-grid geometrically thin $\alpha$-disc. After validating our model semi-analytically, we test it in an idealized environment consisting of a circumnuclear disc, where gas accretion onto the accretion disc is consistently determined by GIZMO. In the absence of any accretion-related feedback, the spin evolution closely traces that observed in the semi-analytical models and depends on the free parameters of our implementation, such as the initial BH spin, angular momentum of the accretion disc, and the radius at which the gas inflow circularises. In GIZMO, we also couple our model with the biconical-outflow model presented in a companion paper, wherein the feedback axis is always aligned with the BH spin. In this last case, the evolution of the central BH differs significantly from the previous cases, since the feedback process modifies the gas dynamics and its inflow rates from resolved scales. Such an interaction cannot be modeled by simple semi-analytical models and should be treated using full N-body hydrodynamical simulations., Comment: 10 pages, 5 figures, accepted for publication in MNRAS. Updated with published version, minimal changes

The molecular gas content of normal galaxies at z>4 is poorly constrained, because the commonly used molecular gas tracers become hard to detect. We use the [CII]158um luminosity, recently proposed as a molecular gas tracer, to estimate the molecular gas content in a large sample of main-sequence star-forming galaxies at z=4.4-5.9, with a median stellar mass of 10^9.7 Msun, drawn from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. The agreement between molecular gas masses derived from [CII] luminosity, dynamical mass, and rest-frame 850um luminosity, extrapolated from the rest-frame 158um continuum, supports [CII] as a reliable tracer of molecular gas in our sample. We find a continuous decline of the molecular gas depletion timescale from z=0 to z=5.9, which reaches a mean value of (4.6+/-0.8)x10^8 yr at z~5.5, only a factor of 2-3 shorter than in present-day galaxies. This suggests a mild enhancement of star formation efficiency toward high redshifts, unless the molecular gas fraction significantly increases. Our estimates show that the rise in molecular gas fraction as reported previously, flattens off above z~3.7 to achieve a mean value of 63%+/-3 over z=4.4-5.9. This redshift evolution of the gas fraction is in line with the one of the specific star formation rate. We use multi-epoch abundance matching to follow the gas fraction evolution over cosmic time of progenitors of z=0 Milky Way-like galaxies in 10^13 Msun halos and of more massive z=0 galaxies in 10^14 Msun halos. Interestingly, the former progenitors show a monotonic decrease of the gas fraction with cosmic time, while the latter show a constant gas fraction from z=5.9 to z~2 and a decrease at z, 16 pages, 9 figures, 2 tables. Submitted to Astronomy & Astrophysics

We use the Fundamental Plane (FP) to measure the redshift evolution of the dynamical mass-to-light ratio ($M_{\mathrm{dyn}}/L$) and the dynamical-to-stellar mass ratio ($M_{\mathrm{dyn}}/M_*$). Although conventionally used to study the properties of early-type galaxies, we here obtain stellar kinematic measurements from the Large Early Galaxy Astrophysics Census (LEGA-C) Survey for a sample of $\sim1400$ massive ($\log( M_*/M_\odot) >10.5$) galaxies at $0.6, Comment: 7 pages, 4 figures; accepted for publication in ApJL

We use the IllustrisTNG simulations to show how the fractions of quenched galaxies vary across different environments and cosmic time, and to quantify the role AGN feedback and preprocessing play in quenching group and cluster satellites. At $z=0$, we select galaxies with $M_* = 10^{9-12} M_{\odot}$ residing within ($\leq R_{200c}$) groups and clusters of total host mass $M_{200c}=10^{13-15.2} M_{\odot}$. TNG predicts a quenched fraction of $\sim70-90\%$ (on average) for centrals and satellites $\gtrsim 10^{10.5} M_{\odot}$, regardless of host mass, cosmic time ($0\leq z\leq0.5$), clustercentric distance and time since infall in the $z=0$ host. Low-mass centrals ($\lesssim 10^{10} M_{\odot}$), instead, are rarely quenched unless they become members of groups ($10^{13-14} M_{\odot}$) or clusters ($\geq10^{14} M_{\odot}$), where the quenched fraction rises to $\sim80\%$. The fraction of low-mass passive galaxies is higher closer to the host center and for more massive hosts. The population of low-mass satellites accreted $\gtrsim$4-6 Gyr ago in massive hosts is almost entirely passive, thus suggesting an upper limit for the time needed for environmental quenching to occur. In fact, $\sim30\%$ of group and cluster satellites that are quenched at $z=0$ were already quenched before falling into their current host, and the bulk of them quenched as early as 4 to 10 billion years ago. For low-mass galaxies ($\lesssim10^{10-10.5}M_{\odot}$), this is due to preprocessing, whereby current satellites may have been members of other hosts, and hence have undergone environmental processes, before falling into their final host, this mechanism being more common and more effective for the purposes of quenching for satellites found today in more massive hosts. On the other hand, massive galaxies quench on their own and because of AGN feedback, regardless of whether they are centrals or satellites., 21 Pages, 11 figures, 2 tables, 1 appendix, published in MNRAS

Ultracompact massive galaxies (ucmgs), i.e., galaxies with stellar masses and effective radii , are very rare systems, in particular at low and intermediate redshifts. Their origin as well as their number density across cosmic time are still under scrutiny, especially because of the paucity of spectroscopically confirmed samples. We have started a systematic census of ucmg candidates within the ESO Kilo Degree Survey, together with a large spectroscopic follow-up campaign to build the largest possible sample of confirmed ucmgs. This is the third paper of the series and the second based on the spectroscopic follow-up program. Here, we present photometrical and structural parameters of 33 new candidates at redshifts and confirm 19 of them as ucmgs, based on their nominal spectroscopically inferred and . This corresponds to a success rate of , nicely consistent with our previous findings. The addition of these 19 newly confirmed objects allows us to fully assess the systematics on the system selection—and to finally reduce the number density uncertainties. Moreover, putting together the results from our current and past observational campaigns and some literature data, we build the largest sample of ucmgs ever collected, comprising 92 spectroscopically confirmed objects at . This number raises to 116, allowing for a 3σ tolerance on the and thresholds for the ucmg definition. For all these galaxies, we have estimated the velocity dispersion values at the effective radii, which have been used to derive a preliminary mass–velocity dispersion correlation.

We study the displacements between the centres of galaxies and their supermassive black holes (BHs) in the cosmological hydrodynamical simulation Horizon-AGN, and in a variety of observations from the literature. The BHs in Horizon-AGN feel a sub-grid dynamical friction force, sourced by the surrounding gas, which prevents recoiling BHs being ejected from the galaxy. We find that i) the fraction of spatially offset BHs increases with cosmic time, ii) BHs live on prograde orbits in the plane of the galaxy with an orbital radius that decays with time but stalls near $z=0$, and iii) the magnitudes of offsets from the galaxy centres are substantially larger in the simulation than in observations. We attribute the stalling of the infall and excessive offset magnitudes to the fact that dynamical friction from stars and dark matter is not modelled in the simulation, and hence provide a way to improve the black hole dynamics of future simulations., 19 pages, 15 figures, accepted for publication in MNRAS

The majority of galaxies with current star-formation rates (SFRs), SFRo >= 10^-3 Msun/yr, in the Local Cosmological Volume where observations should be reliable, have the property that their observed SFRo is larger than their average star formation rate. This is in tension with the evolution of galaxies described by delayed-tau models, according to which the opposite would be expected. The tension is apparent in that local galaxies imply the star formation timescale tau approx 6.7 Gyr, much longer than the 3.5-4.5 Gyr obtained using an empirically determined main sequence at several redshifts. Using models where the SFR is a power law in time of the form propto (t - t1)^eta for t1 = 1.8 Gyr (with no stars forming prior to t1) implies that eta = 0.18 +- 0.03. This suggested near-constancy of a galaxy's SFR over time raises non-trivial problems for the evolution and formation time of galaxies, but is broadly consistent with the observed decreasing main sequence with increasing age of the Universe., LaTeX, 8 pages, 5 figures, MNRAS, in press; the replacement is to fix an error in the metadata (paper remains unaffected)

Understanding the processes that trigger morphological transformation is central to understanding how and why the Universe transitions from being disc-dominated at early epochs to having the morphological mix that is observed today. We use Horizon-AGN, a cosmological hydrodynamical simulation, to perform a comprehensive study of the processes that drive morphological change in massive (M > 10^10 MSun) galaxies over cosmic time. We show that (1) essentially all the morphological evolution in galaxies that are spheroids at z=0 is driven by mergers with mass ratios greater than 1:10, (2) major mergers alone cannot produce today's spheroid population -- minor mergers are responsible for a third of all morphological transformation over cosmic time and are its dominant driver after z~1, (3) prograde mergers trigger milder morphological transformation than retrograde mergers -- while both types of events produce similar morphological changes at z>2, the average change due to retrograde mergers is around twice that due to their prograde counterparts at z~0, (4) remnant morphology depends strongly on the gas fraction of a merger, with gas-rich mergers routinely re-growing discs, and (5) at a given stellar mass, discs do not exhibit drastically different merger histories from spheroids -- disc survival in mergers is driven by acquisition of cold gas (via cosmological accretion and gas-rich interactions) and a preponderance of prograde mergers in their merger histories., Accepted for publication in MNRAS

We investigate the connection between environment and the different quenching channels that galaxies are prone to follow in the rest-frame NUVrK colour diagram, as identified by Moutard et al. (2016b). Namely, the fast quenching channel followed by $young$ low-mass galaxies and the slow quenching channel followed by $old$ high-mass ones. We make use of the >22 deg$^2$ covered the VIPERS Multi-Lambda Survey (VIPERS-MLS) to select a galaxy sample complete down to stellar masses of $M_* > 10^{9.4} M_\odot$ at $z < 0.65$ ($M_* > 10^{8.8} M_\odot$ at $z < 0.5$) and including 33,500 (43,000) quiescent galaxies properly selected at $0.2 < z < 0.65$, while being characterized by reliable photometric redshifts ($��_{��z/(1+z)} \leq 0.04$) that we use to measure galaxy local densities. We find that (1) the quiescence of low-mass [$M_* \leq 10^{9.7} M_\odot$] galaxies requires a strong increase of the local density, which confirms the lead role played by environment in their fast quenching and, therefore, confirms that the low-mass upturn observed in the stellar mass function of quiescent galaxies is due to $environmental$ $quenching$. We also observe that (2) the reservoir of low-mass star-forming galaxies located in very dense regions (prone to environmental quenching) has grown between $z \sim 0.6$ and $ z \sim 0.4$ whilst the share of low-mass quiescent galaxies (expected to being environmentally quenched) may have simultaneously increased, which would plead for a rising importance of $environmental$ $quenching$ with cosmic time, compared to $mass$ $quenching$. We finally discuss the composite picture of such environmental quenching of low-mass galaxies and, in particular, how this picture may be consistent with a $delayed$-$then$-$rapid$ quenching scenario., 14 pages, 5 figures, accepted for publication in MNRAS; published version

In this paper an explicit cosmological model which allows cosmological singularities are discussed in Kaluza–Klein space time. The generalized power-law and asymptotic expansions of the baro-tropic fluid index ω and equivalently the deceleration parameter q, in terms of cosmic time ‘t’ are considered. Finally, the strength of the found singularities is discussed.

The dearth of dwarf galaxies in the local universe is hard to reconcile with the large number of low mass haloes expected within the concordance $\Lambda$CDM paradigm. In this paper we perform a systematic evaluation of the uncertainties affecting the measurement of DM halo abundance using galaxy kinematics. Using a large sample of dwarf galaxies with spatially resolved kinematic data we derive a correction to obtain the observed abundance of galaxies as a function of their halo maximum circular velocity from the line-of-sight velocity function in the Local Volume. This estimate provides a direct means of comparing the predictions of theoretical models and simulations (including nonstandard cosmologies and novel galaxy formation physics) to the observational constraints. The new "galactic $V_{max}$" function is steeper than the line-of-sight velocity function but still shallower than the theoretical CDM expectation, showing that some unaccounted physical process is necessary to reduce the abundance of galaxies and/or drastically modify their density profiles compared to CDM haloes. Using this new galactic $V_{max}$ function, we investigate the viability of baryonic solutions such as feedback-powered outflows and photoevaporation of gas from an ionising radiation background. At the 3-$\sigma$ confidence level neither energetic feedback nor photoevaporation are effective enough to reconcile the disagreement. In the case of maximum baryonic effects, the theoretical estimate still deviates significantly from the observations for $V_{max} < 60$ km/s. CDM predicts at least 1.8 times more galaxies with $V_{max} = 50$ km/s and 2.5 times more than observed at $30$ km/s. Recent hydrodynamic simulations seem to resolve the discrepancy but disagree with the properties of observed galaxies with resolved kinematics. (abridged), Comment: 17 pages, 22 figures; major revisions include clarification of the method, expanded comparison with simulations with a new figure, analysis of uncertainties in model as well as pressure support corrections, and a new table with nomenclature

We explore the cosmic evolution of the bar length, strength, and light deficit around the bar for 379 barred galaxies at 0.2 < z $\leq$ 0.835 using F814W images from the COSMOS survey. Our sample covers galaxies with stellar mass 10.0 $\leq$ log(M*/Msun) $\leq$ 11.4 and various Hubble types. The bar length is strongly related to the galaxy mass, the disk scale length (h), R50, and R90, where the last two are the radii containing 50 and 90% of total stellar mass, respectively. Bar length remains almost constant, suggesting little or no evolution in bar length over the last 7 Gyrs. The normalized bar lengths (Rbar/h, Rbar/R50, and Rbar/R90) do not show any clear cosmic evolution. Also, the bar strength (A2 and Qb) and the light deficit around the bar reveal little or no cosmic evolution. The constancy of the normalized bar lengths over cosmic time implies that the evolution of bars and of disks is strongly linked over all times. We discuss our results in the framework of predictions from numerical simulations. We conclude there is no strong disagreement between our results and up-to-date simulations., 16 pages, 8 figures, Accepted for publication in ApJ

In Dou et al. (2021), we introduced the Fundamental Formation Relation (FFR), a tight relation between specific SFR (sSFR), H$_2$ star formation efficiency (SFE$_{\rm H_2}$), and the ratio of H$_2$ to stellar mass. Here we show that atomic gas HI does not follow a similar FFR as H$_2$. The relation between SFE$_{\rm HI}$ and sSFR shows significant scatter and strong systematic dependence on all of the key galaxy properties that we have explored. The dramatic difference between HI and H$_2$ indicates that different processes (e.g., quenching by different mechanisms) may have very different effects on the HI in different galaxies and hence produce different SFE$_{\rm HI}$-sSFR relations, while the SFE$_{\rm H_2}$-sSFR relation remains unaffected. The facts that SFE$_{\rm H_2}$-sSFR relation is independent of other key galaxy properties, and that sSFR is directly related to the cosmic time and acts as the cosmic clock, make it natural and very simple to study how different galaxy populations (with different properties and undergoing different processes) evolve on the same SFE$_{\rm H_2}$-sSFR $\sim t$ relation. In the gas regulator model (GRM), the evolution of a galaxy on the SFE$_{\rm H_2}$-sSFR($t$) relation is uniquely set by a single mass-loading parameter $\lambda_{\rm net,H_2}$. This simplicity allows us to accurately derive the H$_2$ supply and removal rates of the local galaxy populations with different stellar masses, from star-forming galaxies to the galaxies in the process of being quenched. This combination of FFR and GRM, together with the stellar metallicity requirement, provide a new powerful tool to study galaxy formation and evolution., Comment: 16 pages, 7 figures, accepted for publication in ApJ

We use high-resolution, multi-band imaging of ~16,500 galaxies in the CANDELS fields at 0 < z < 2.5 to study the evolution of color gradients and half-mass radii over cosmic time. We find that galaxy color gradients at fixed mass evolve rapidly between z~2.5 and z~1, but remain roughly constant below z~1. This result implies that the sizes of both star-forming and quiescent galaxies increase much more slowly than previous studies found using half-light radii. The half-mass radius evolution of quiescent galaxies is fully consistent with a model which uses observed minor merger rates to predict the increase in sizes due to the accretion of small galaxies. Progenitor bias may still contribute to the growth of quiescent galaxies, particularly if we assume a slower timescale for the minor merger growth model. The slower half-mass radius evolution of star-forming galaxies is in tension with cosmological simulations and semi-analytic galaxy models. Further detailed, consistent comparisons with simulations are required to place these results in context., Comment: 6 pages, 4 figures. Accepted to ApJL

We explore deep rest-frame UV to FIR data in the COSMOS field to measure the individual spectral energy distributions (SED) of the ~4000 SC4K (Sobral et al. 2018) Lyman-alpha (Lya) emitters (LAEs) at z~2-6. We find typical stellar masses of 10$^{9.3\pm0.6}$ M$_{\odot}$ and star formation rates (SFR) of SFR$_{SED}=4.4^{+10.5}_{-2.4}$ M$_{\odot}$/yr and SFR$_{Lya}=5.9^{+6.3}_{-2.6}$ M$_{\odot}$/yr, combined with very blue UV slopes of beta=-2.1$^{+0.5}_{-0.4}$, but with significant variations within the population. M$_{UV}$ and beta are correlated in a similar way to UV-selected sources, but LAEs are consistently bluer. This suggests that LAEs are the youngest and/or most dust-poor subset of the UV-selected population. We also study the Lya rest-frame equivalent width (EW$_0$) and find 45 "extreme" LAEs with EW$_0>240$ A (3 $��$), implying a low number density of $(7\pm1)\times10^{-7}$ Mpc$^{-3}$. Overall, we measure little to no evolution of the Lya EW$_0$ and scale length parameter ($w_0$) which are consistently high (EW$_0=140^{+280}_{-70}$ A, $w_0=129^{+11}_{-11}$ A) from z~6 to z~2 and below. However, $w_0$ is anti-correlated with M$_{UV}$ and stellar mass. Our results imply that sources selected as LAEs have a high Lya escape fraction (f$_{esc, Lya}$) irrespective of cosmic time, but f$_{esc, Lya}$ is still higher for UV-fainter and lower mass LAEs. The least massive LAEs ($, Accepted for publication in MNRAS. 23 pages, 11 figures, 3 tables + appendices. The full SC4K catalogue of LAEs with PSF photometry, stellar masses, SFRs and other properties for individual LAEs, is available at https://goo.gl/q9yfKo

Despite recent progress in understanding Lyman-alpha (Lya) emitters (LAEs), relatively little is known regarding their typical black hole activity across cosmic time. Here, we study the X-ray and radio properties of ~4000 LAEs at 2.23) and find an average luminosity of 10^{44.31+-0.01} erg/s and an average black hole accretion rate (BHAR) of 0.72+-0.01 Msun/yr, consistent with moderate to high accreting AGN. We detect 120 sources in deep radio data (radio AGN fraction of 3.2%+-0.3%). The global AGN fraction (8.6%+-0.4%) rises with Lya luminosity and declines with increasing redshift. For X-ray detected LAEs, Lya luminosities correlate with the BHARs, suggesting that Lya luminosity becomes an accretion rate indicator. Most LAEs (93.1+-0.6%) at 2, Comment: Accepted for publication in MNRAS. 21 pages, 15 figures, 2 tables + appendices. The full SC4K catalogue of LAEs with X-ray, radio, FIR and Lya measurements is available with the MNRAS refereed version, along with tables presenting the full stacking results in X-rays and radio

In the paper 'NIHAO IV: Core creation and destruction in dark matter density profiles across cosmic time' published in MNRAS main journal, Volume 456, Issue 4, p. 3542-3552, an error in table 2 came to our attention. It appears that there is a mismatch between the values on the second line of table 2 and the corresponding curve in fig. 5. Indeed the parameters given for the fit of the relation between α and Mh are incorrect. We realised that the published parameters correspond to a different fitting formula that was used internally during thiswork and that does not have the right asymptotic behavior when the mass goes to zero. We apologize for this mistake. (Table presented) The correct fitting parameters to reproduce the fitting curve on fig. 5 are given in Table Errata. These parameters were provided only to allow the reader to reproduce our work and therefore this error does not affect our conclusions.

We use the IllustrisTNG cosmological hydrodynamical simulation to study the evolution of star formation rate (SFR)-density relation over cosmic time. We construct several samples of galaxies at different redshifts from $z=2.0$ to $z=0.0$, which have the same comoving number density. The SFR of galaxies decreases with local density at $z=0.0$, but its dependence on local density becomes weaker with redshift. At $z\gtrsim1.0$, the SFR of galaxies increases with local density (reversal of the SFR-density relation), and its dependence becomes stronger with redshift. This change of SFR-density relation with redshift still remains even when fixing the stellar masses of galaxies. The dependence of SFR on the distance to a galaxy cluster also shows a change with redshift in a way similar to the case based on local density, but the reversal happens at a higher redshift, $z\sim1.5$, in clusters. On the other hand, the molecular gas fraction always decreases with local density regardless of redshift at $z=0.0-2.0$ even though the dependence becomes weaker when we fix the stellar mass. Our study demonstrates that the observed reversal of the SFR-density relation at $z\gtrsim1.0$ can be successfully reproduced in cosmological simulations. Our results are consistent with the idea that massive, star-forming galaxies are strongly clustered at high redshifts, forming larger structures. These galaxies then consume their gas faster than those in low-density regions through frequent interactions with other galaxies, ending up being quiescent in the local universe., 10 pages, 8 figures. To appear in MNRAS. Paper with high resolution figures is available at https://astro.kias.re.kr/~hshwang/Hwang_etal2019_TNG_Reversal_HR.pdf

We investigate the absorption features associated with a gas-rich dwarf galaxy using cosmological hydrodynamics simulations. Our goal is to explore whether the progenitors of the lowest mass dwarf galaxies known to harbor neutral hydrogen today (M_star~10^6 solar mass, M_halo=4x10^9 solar mass) could possibly be detected as Damped Lyman-alpha Absorbers (DLAs) over cosmic time. We trace the evolution of a single dwarf galaxy, pre-selected to contain DLAs, from the era of the first metal-free, so-called Population~III (Pop~III), stars, down to z=0, thus allowing us to study the metal enrichment history of DLAs associated with the simulated galaxy. We find that the progenitors of the simulated dwarf are expected to be seen for most of their evolution as DLAs that are contaminated by normal, Population~II, stars. The time period during which DLAs are only metal-enriched by Pop~III stars, on the other hand, is likely very brief, confined to high redshifts, z~6. The susceptibility of the dwarfs to the external UV radiation background allows them to preserve neutral gas only at the centre (a few ~100 pc). This results in a small probability that the simulated dwarf would be observed as a DLA. This study suggests that DLAs are unlikely to be hosted in the lowest mass dwarfs that can harbor neutral gas (M_halo~ 4x10^9 solar mass), below which neutral gas is unlikely to exist. However, this study does illustrate that, when detected, absorption lines provide a powerful method for probing ISM conditions inside the smallest dwarf galaxies at intermediate to high redshifts., Comment: 13 pages, 14 figures, Accepted for Publication in ApJ

The formation histories of globular clusters (GCs) are a key diagnostic for understanding their relation to the evolution of the Universe through cosmic time. We use the suite of 25 cosmological zoom-in simulations of present-day Milky Way-mass galaxies from the E-MOSAICS project to study the formation histories of stars, clusters, and GCs, and how these are affected by the environmental dependence of the cluster formation physics. We find that the median lookback time of GC formation in these galaxies is ${\sim}10.73~$Gyr ($z=2.1$), roughly $2.5~$Gyr earlier than that of the field stars (${\sim}8.34~$Gyr or $z=1.1$). The epoch of peak GC formation is mainly determined by the time evolution of the maximum cluster mass, which depends on the galactic environment and largely increases with the gas pressure. Different metallicity subpopulations of stars, clusters and GCs present overlapping formation histories, implying that star and cluster formation represent continuous processes. The metal-poor GCs ($-2.5, Comment: 15 pages, 6 figures; accepted by MNRAS (April 30, 2019)

Current galaxy observations suggest that a roughly linear correlation exists between the [CII] emission and the star formation rate, either as spatially-resolved or integrated quantities. Observationally, this correlation seems to be independent of metallicity, but the very large scatter does not allow to properly assess whether this is true. On the other hand, theoretical models tend to suggest a metallicity dependence of the correlation. In this study, we investigate the metallicity evolution of the correlation via a high-resolution zoom-in cosmological simulation of a dwarf galaxy employing state-of-the-art sub-grid modelling for gas cooling, star formation, and stellar feedback, and that self-consistently evolves the abundances of metal elements out of equilibrium. Our results suggest that the correlation should evolve with metallicity, in agreement with theoretical predictions, but also that this evolution can be hardly detected in observations, because of the large scatter. We also find that most of the [CII] emission is associated with neutral gas at low-intermediate densities, whereas the highest emissivity is produced by the densest regions around star-forming regions., 10 pages, 12 figures; 1 table; accepted for publication by MNRAS

Faint tidal features around galaxies record their merger and interaction histories over cosmic time. Due to their low surface brightnesses and complex morphologies, existing automated methods struggle to detect such features and most work to date has heavily relied on visual inspection. This presents a major obstacle to quantitative study of tidal debris features in large statistical samples, and hence the ability to be able to use these features to advance understanding of the galaxy population as a whole. This paper uses convolutional neural networks (CNNs) with dropout and augmentation to identify galaxies in the CFHTLS-Wide Survey that have faint tidal features. Evaluating the performance of the CNNs against previously-published expert visual classifications, we find that our method achieves high (76%) completeness and low (20%) contamination, and also performs considerably better than other automated methods recently applied in the literature. We argue that CNNs offer a promising approach to effective automatic identification of low surface brightness tidal debris features in and around galaxies. When applied to forthcoming deep wide-field imaging surveys (e.g. LSST, Euclid), CNNs have the potential to provide a several order-of-magnitude increase in the sample size of morphologically-perturbed galaxies and thereby facilitate a much-anticipated revolution in terms of quantitative low surface brightness science., 16 pages, 18 figures, accepted for publication in MNRAS

The light-cone (LC) effect causes the mean as well as the statistical properties of the redshifted 21-cm signal $T_{\rm b}(\hat{\bf n},\nu)$ to change with frequency $\nu$ (or cosmic time). Consequently, the statistical homogeneity (ergodicity) of the signal along the line of sight (LoS) direction is broken. This is a severe problem particularly during the Epoch of Reionization (EoR) when the mean neutral hydrogen fraction ($\bar{x}_{\rm HI}$) changes rapidly as the universe evolves. This will also pose complications for large bandwidth observations. These effects imply that the 3D power spectrum $P(k)$ fails to quantify the entire second-order statistics of the signal as it assumes the signal to be ergodic and periodic along the LoS. As a proper alternative to $P(k)$, we use the multi-frequency angular power spectrum (MAPS) ${\mathcal C}_{\ell}(\nu_1,\nu_2)$ which does not assume the signal to be ergodic and periodic along the LoS. Here, we study the prospects for measuring the EoR 21-cm MAPS using future observations with the upcoming SKA-Low. Ignoring any contribution from the foregrounds, we find that the EoR 21-cm MAPS can be measured at a confidence level $\ge 5\sigma$ at angular scales $\ell \sim 1300$ for total observation time $t_{\rm obs} \ge 128\,{\rm hrs}$ across $\sim 44\,{\rm MHz}$ observational bandwidth. We also quantitatively address the effects of foregrounds on MAPS detectability forecast by avoiding signal contained within the foreground wedge in $(k_\perp, k_\parallel)$ plane. These results are very relevant for the upcoming large bandwidth EoR experiments as previous predictions were all restricted to individually analyzing the signal over small frequency (or equivalently redshift) intervals., Comment: Published in Monthly Notices of the Royal Astronomical Society (MNRAS). Available at https://doi.org/10.1093/mnras/staa1026

We explore the host galaxies of compact-object binaries (black hole–black hole binaries, BHBs; neutron star–black hole binaries, NSBHs; double–neutron stars; DNSs) across cosmic time, by means of population-synthesis simulations combined with the Illustris cosmological simulation. At high redshift (z ≳ 4), the host galaxies of BHBs, NSBHs, and DNSs are very similar and are predominantly low-mass galaxies (stellar mass M < 1011 M⊙). If z ≳ 4, most compact objects form and merge in the same galaxy, with a short delay time. At low redshift (z ≤ 2), the host galaxy populations of DNSs differ significantly from the host galaxies of both BHBs and NSBHs. DNSs merging at low redshift tend to form and merge in the same galaxy, with relatively short delay time. The stellar mass of DNS hosts peaks around ∼1010–1011 M⊙. In contrast, BHBs and NSBHs merging at low redshift tend to form in rather small galaxies at high redshift and then to merge in larger galaxies with long delay times. This difference between DNSs and black hole binaries is a consequence of their profoundly different metallicity dependence.

Cosmological simulations almost invariably estimate the accretion of gas on to supermassive black holes using a Bondi-Hoyle-like prescription. Doing so ignores the effects of the angular momentum of the gas, which may prevent or significantly delay accreting material falling directly on to the black hole. We outline a black hole accretion rate prescription using a modified Bondi-Hoyle formulation that takes into account the angular momentum of the surrounding gas. Meaningful implementation of this modified Bondi-Hoyle formulation is only possible when the inner vorticity distribution is well resolved, which we achieve through the use of a super-Lagrangian refinement technique around black holes within our simulations. We then investigate the effects on black hole growth by performing simulations of isolated as well as merging disc galaxies using the moving-mesh code AREPO. We find that the gas angular momentum barrier can play an important role in limiting the growth of black holes, leading also to a several Gyr delay between the starburst and the quasar phase in major merger remnants. We stress, however, that the magnitude of this effect is highly sensitive to the thermodynamical state of the accreting gas and to the nature of the black hole feedback present., Comment: 16 pages, 8 figures, MNRAS accepted

Hydrodynamical shocks are a manifestation of the non-linearity of the Euler equations and play a fundamental role in cosmological gas dynamics. In this work, we identify and analyse shocks in the Illustris simulation, and contrast the results with those of non-radiative runs. We show that simulations with more comprehensive physical models of galaxy formation pose new challenges for shock finding algorithms due to radiative cooling and star-forming processes, prompting us to develop a number of methodology improvements. We find in Illustris a total shock surface area which is about 1.4 times larger at the present epoch compared to non-radiative runs, and an energy dissipation rate at shocks which is higher by a factor of around 7. Remarkably, shocks with Mach numbers above and below $\mathcal{M}\approx10$ contribute about equally to the total dissipation across cosmic time. This is in sharp contrast to non-radiative simulations, and we demonstrate that a large part of the difference arises due to strong black hole radio-mode feedback in Illustris. We also provide an overview of the large diversity of shock morphologies, which includes complex networks of halo-internal shocks, shocks on to cosmic sheets, feedback shocks due to black holes and galactic winds, as well as ubiquitous accretion shocks. In high redshift systems more massive than $10^{12}\,\mathrm{M}_\odot$ we discover the existence of a double accretion shock pattern in haloes. They are created when gas streams along filaments without being shocked at the outer accretion shock, but then forms a second, roughly spherical accretion shock further inside., Comment: 26 pages, 15 figures, 1 movie, published in mnras

We quantify the frequency of companions of low redshift ($0.013 < z < 0.0252$), dwarf galaxies ($2 \times 10^8$ M$_\odot 55'', projected separations r$_p < 150$ kpc and relative line-of-sight velocities $\Delta V_{\rm LOS} < 150$ km/s. While the mock catalogs predict a factor of 2 more isolated dwarfs than observed in SDSS, the mean number of observed companions per dwarf is $N_c \sim 0.04$, in good agreement with Illustris when accounting for SDSS sensitivity limits. Removing these limits in the mock catalogs predicts $N_c\sim 0.06$ for future surveys (LSST, DESI), which will be complete to M$_* = 2\times 10^8$ M$_\odot$. The 3D separations of mock dwarf multiples reveal a contamination fraction of $\sim$40% in observations from projection effects. Most isolated multiples are pairs; triples are rare and it is cosmologically improbable that bound groups of dwarfs with more than 3 members exist within the parameter range probed in this study. We find that $$1:4) steadily increases with decreasing Primary stellar mass, whereas the cosmological "Major Merger rate'' (per Gyr) has the opposite behaviour. We conclude that cosmological simulations can be reliably used to constrain the fraction of dwarf mergers across cosmic time., Comment: 21 pages, 16 Figures, submitted to MNRAS (Referee stage ; Comments are welcome)

One of the fundamental assumptions of chemical evolution models (CEMs) of the Milky Way (MW) and other spirals is that higher gas accretion rates are expected in the past, and in the inner regions of the Galaxy. This leads to the so-called `inside-out disc formation scenario'. Yet, these are probably the most unconstrained inputs of such models. In the present paper, we aim at investigating these main assumptions by studying how gas is accreted in four simulated MW-like galaxies assembled within the $\Lambda$CDM scenario. The galaxies were obtained using two different simulation techniques, cosmological setups and initial conditions. Two of them are MW candidates corresponding to the chemodynamical model of Minchev et al. (2013, 2014) (known as MCM) and the Local Group cosmological simulation of Nuza et al. (2014). We investigate vertical and radial gas accretion on to galaxy discs as a function of cosmic time and disc radius. We find that accretion in the MW-like galaxies seem to happen in two distinct phases, namely: an early, more violent period; followed by a subsequent, slowly declining phase. Our simulations seem to give support to the assumption that the amount of gas incorporated into the MW disc exponentially decreases with time, leading to current net accretion rates of $0.6-1\,$M$_\odot\,$yr$^{-1}$. In particular, accretion timescales on to the simulated thin-disc-like structures are within $\sim5-7\,$Gyr, consistent with expectations from CEMs. Moreover, our simulated MW discs are assembled from the inside-out with gas in the inner disc regions accreted in shorter timescales than in external ones, in qualitative agreement with CEMs of the Galaxy. However, this type of growth is not general to all galaxies and it is intimately linked to their particular merger and gas accretion history., Comment: Replaced to match published version. References and table added; 20 pages, 13 figures, 2 tables. Gas accretion timescale vs. galactocentric radius for simulated MW candidates shown in Table B1

Observationally, the fraction of blue satellite galaxies decreases steeply with host halo mass, and their radial distribution around central galaxies is significantly shallower in massive (M_* >10e11M_sun) than in Milky Way like systems. Theoretical models, based primarily on semi-analytical techniques, have had a long-standing problem with reproducing these trends, instead predicting too few blue satellites in general but also estimating a radial distribution that is too shallow, regardless of primary mass. In this Letter, we use the Illustris cosmological simulation to study the properties of satellite galaxies around isolated primaries. For the first time, we find good agreement between theory and observations. We identify the main source of this success relative to earlier work to be a consequence of the large gas contents of satellites at infall, a factor ~5-10 times larger than in semi-analytical models. Because of their relatively large gas reservoirs, satellites can continue to form stars long after infall, with a typical timescale for star-formation to be quenched ~2 Gyr in groups but more than ~5 Gyr for satellites around Milky Way like primaries. The gas contents we infer are consistent with z=0 observations of HI gas in galaxies, although we find large discrepancies among reported values in the literature. A testable prediction of our model is that the gas-to-stellar mass ratio of satellite progenitors should vary only weakly with cosmic time., 5 pages, 3 figures. Accepted for publication in MNRAS Letters

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., 12 pages, 9 figures, accepted for publication in MNRAS

Answers to the metal production of the Universe can be found in galaxy clusters, notably within their Intra-Cluster Medium (ICM). The X-ray Integral Field Unit (X-IFU) on board the next-generation European X-ray observatory Athena (2030s) will provide the necessary leap forward in spatially-resolved spectroscopy required to disentangle the intricate mechanisms responsible for this chemical enrichment. In this paper, we investigate the future capabilities of the X-IFU in probing the hot gas within galaxy clusters. From a test sample of four clusters extracted from cosmological hydrodynamical simulations, we present comprehensive synthetic observations of these clusters at different redshifts (up to z = 2) and within the scaled radius R500 performed using the instrument simulator SIXTE. Through 100 ks exposures, we demonstrate that the X-IFU will provide spatially-resolved mapping of the ICM physical properties with little to no biases (, 22 pages, 15 figures

We study the effects of the local environment on the molecular gas content of a large sample of log($M_{*}$/$M_{\odot}$) $\gtrsim$ 10 star-forming and starburst galaxies with specific star-formation rates (sSFRs) on and above the main sequence (MS) to $z$ $\sim$ 3.5. ALMA observations of the dust continuum in the COSMOS field are used to estimate molecular gas masses at $z$ $\approx$ 0.5-3.5. We also use a local universe sample from the ALFALFA HI survey after converting it into molecular masses. The molecular mass ($M_{ISM}$) scaling relation shows a dependence on $z$, $M_{*}$, and sSFR relative to the MS, but no dependence on environmental overdensity $\Delta$ ($M_{ISM}$ $\propto$ $\Delta^{0.03}$). Similarly, gas mass fraction (f$_{gas}$) and depletion timescale ($\tau$) show no environmental dependence to $z$ $\sim$ 3.5. At $\langle z\rangle$ $\sim$ 1.8, the average $\langle M_{ISM}\rangle$,$\langle$f$_{gas}\rangle$, and $\langle \tau \rangle$ in densest regions is (1.6$\pm$0.2)$\times$10$^{11}$ $M_{\odot}$, 55$\pm$2%, and 0.8$\pm$0.1 Gyr, respectively, similar to those in the lowest density bin. Independent of the environment, f$_{gas}$ decreases and $\tau$ increases with increasing cosmic time. Cosmic molecular mass density ($\rho$) in the lowest density bins peaks at $z$ $\sim$ 1-2, and this peak happens at $z$ $, Comment: ApJ published

Using a novel approach, we study the quenching and bursting of galaxies as a function of stellar mass ($M_{*}$), local environment ($\Sigma$), and specific star-formation rate (sSFR) using a large spectroscopic sample of $\sim$ 123,000 $GALEX$/SDSS and $\sim$ 420 $GALEX$/COSMOS/LEGA-C galaxies to $z$ $\sim$ 1. We show that out to $z$ $\sim$ 1 and at fixed sSFR and local density, on average, less massive galaxies are quenching, whereas more massive systems are bursting, with a quenching/bursting transition at log($M_{*}$/$M_{\odot}$) $\sim$ 10.5-11 and likely a short quenching/bursting timescale ($\lesssim$ 300 Myr). We find that much of the bursting of star-formation happens in massive (log($M_{*}$/$M_{\odot}$) $\gtrsim$ 11), high sSFR galaxies (log(sSFR/Gyr$^{-1}$) $\gtrsim$ -2), particularly those in the field (log($\Sigma$/Mpc$^{-2}$) $\lesssim$ 0; and among group galaxies, satellites more than centrals). Most of the quenching of star-formation happens in low-mass (log($M_{*}$/$M_{\odot}$) $\lesssim$ 9), low sSFR galaxies (log(sSFR/Gyr$^{-1}$) $\lesssim$ -2), in particular those located in dense environments (log($\Sigma$/Mpc$^{-2}$) $\gtrsim$ 1), indicating the combined effects of $M_{*}$ and $\Sigma$ in quenching/bursting of galaxies since $z$ $\sim$ 1. However, we find that stellar mass has stronger effects than environment on recent quenching/bursting of galaxies to $z$ $\sim$ 1. At any given $M_{*}$, sSFR, and environment, centrals are quenchier (quenching faster) than satellites in an average sense. We also find evidence for the strength of mass and environmental quenching being stronger at higher redshift. Our preliminary results have potential implications for the physics of quenching/bursting in galaxies across cosmic time., Comment: ApJ accepted

Recent discoveries of massive black holes (MBHs) in dwarf galaxies suggest that they may have a more common presence than once thought. Systematic searches are revealing more candidates, but this process could be accelerated by predictions from simulations. We perform a study of several high-resolution, cosmological, zoom-in simulations focusing on dwarf galaxies that host massive black holes at z = 0, with the aim of determining when the black holes are most observable. Larger dwarf galaxies are more likely to host MBHs than those of lower mass. About 50% of the MBHs in dwarfs are not centrally located, but rather are wandering within a few kpc of the galaxy center. The accretion luminosities of MBHs in dwarfs are low throughout cosmic time, rendering them extremely difficult to detect. However, the merger history of these MBHs is optimal for gravitational wave detection by LISA., Comment: replaced with accepted version