Search

The AGORA Cosmorun (arXiv:2106.09738) is a set of hydrodynamical cosmological zoom-in simulations carried out within the AGORA High-resolution Galaxy Simulations Comparison Project (arXiv:1308.2669,arXiv:1610.03066). These simulations show the formation and evolution of a Milky Way-sized galaxy using eight of the most widely used numerical codes in the community (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Gizmo, and Arepo). In this short report, we describe the public release of the raw output data from all of these simulations at z = 8, 7, 6, 5, 4, 3, 2 (plus at z=1, 0 when available), and several metadata files containing the halo centers, virial quantities, and merger trees. The data from even thinner timesteps will be released as soon as the upcoming collaboration papers (VII-IX) are submitted and accepted., Comment: ArXiv Data Release from the AGORA Collaboration. arXiv admin note: text overlap with arXiv:2001.04354

We address the formation of giant clumps in violently unstable gas-rich disc galaxies at cosmic noon. While these are commonly thought to originate from gravitational Toomre instability, cosmological simulations have indicated that clumps form even in regions where the Toomre $Q$ parameter is well above unity, which should be stable according to linear Toomre theory (Inoue et al., 2016). Examining one of these cosmological simulations, we find that it exhibits an excess in compressive modes of turbulence with converging motions. The energy in converging motions within proto-clump regions is $\sim 70\%$ of the total turbulent energy, compared to $\sim 17\%$ expected in equipartition. When averaged over the whole disc, $\sim 32\%$ of the turbulent energy is in converging motions, with a further $\sim 8\%$ in diverging motions. Thus, a total of $\sim 40\%$ of the turbulent energy is in compressive modes, with the rest in solenoidal modes, compared to the $(1/3)-(2/3)$ division expected in equipartition. By contrast, we find that in an isolated-disc simulation with similar properties, resembling high-$z$ star-forming galaxies, the energy in the different turbulence modes are in equipartition, both in proto-clump regions and over the whole disc. We conclude that the origin of the excessive converging motions in proto-clump regions is external to the disc, and propose several mechanisms that can induce them. This is an additional mechanism for clump formation, complementary to and possibly preceding gravitational instability., Comment: 6 pages, 4 figures, Submitted to MNRAS Letters Comments welcome!

We investigate the relation between galaxy structure and star formation rate (SFR) in a sample of $\sim2.9\times10^{4}$ central galaxies with $z<0.0674$ and axial ratios $b/a>0.5$. The star-forming main sequence (SFMS) shows a bend around the stellar mass of $M_\ast\leq{}M_c=2\times10^{10}{}M_{\odot}$. At $M_\ast\leq{}M_c$ the SFMS follows a power-law $\text{SFR}\propto{}M_\ast^{0.85}$, while at higher masses it flattens. $M_c$ corresponds to a dark matter halo mass of $M_\text{vir}\sim{}10^{11.8}M_{\odot}$ where virial shocks occurs. Some galaxy structure (e.g., half-light radius, $R_e$) exhibits a non-monotonic dependence across the SFMS at a fixed $M_\ast$. We find $\text{SFR}\propto{R_e^{-0.28}}$ at fixed $M_\ast$, consistent with the global Kennicutt-Schmidt (KS) law. This finding suggests that galaxy sizes contribute to the scatter of the SFMS. However, at $M_\ast>M_c$ the relationship between SFR and $R_e$ diminishes. Low-mass galaxies above the mean of the SFMS have smaller radii, exhibit compact and centrally concentrated profiles resembling green valley (GV) and quiescent galaxies at the same mass, and have higher $M_{\text{H}_2}/M_\text{HI}$. Conversely, those below the SFMS exhibit larger radii, lower densities, have no GV or quiescent counterparts at their mass and have lower $M_{\text{H}_2}/M_\text{HI}$. The above data suggest two pathways for quenching low-mass galaxies, $M_\ast\leq{}M_c$: a fast one that changes the morphology on the SFMS and a slow one that does not. Above $M_c$, galaxies below the SFMS resemble GV and quiescent galaxies structurally, implying that they undergo a structural transformation already within the SFMS. For these massive galaxies, CG are strongly bimodal, with SFMS galaxies exhibiting negative color gradients, suggesting most star formation occurs in their outskirts, maintaining them within the SFMS., Comment: 30 pages, 28 figures, submitted to MNRAS

We measure the clustering of low-redshift SDSS galaxies as a function of stellar mass ($10.0<\log(M_*/M_\odot)<11.5$) and specific star formation rate (sSFR) and compare the results to models of the galaxy--halo connection. We find that the auto-correlation functions of central galaxies exhibit little dependence on sSFR, with the well-known stronger clustering of quiescent galaxies mainly attributable to satellites. Because halo assembly history is known to affect distinct halo clustering, this result implies that there is little net correlation between halo assembly history and central galaxy sSFR. However, cross-correlations with satellites are stronger for quiescent centrals than star-forming centrals, consistent with quiescent centrals having more satellites in their haloes at fixed $M_*$, as found in SDSS group catalogues. We model the galaxy--halo connection in an $N$-body simulation by assigning sSFRs to central galaxies in three different ways. Two of the models depend on halo assembly history (being based on halo accretion rate or concentration), while the third is independent of halo assembly history (being based on peak halo circular velocity, $V_\text{peak}$, a proxy for halo mass). All three models replicate the observed auto-correlations of central galaxies, while only the $V_\text{peak}$ model reproduces the observed cross-correlations with satellites. This further suggests that the effects of halo assembly history may not be easily seen in auto-correlations of centrals and implies that a more complete understanding of central galaxy clustering may require more than auto-correlations of centrals alone. Additionally, the good agreement with the $V_\text{peak}$ model supports the idea that quiescent galaxies reside in more massive haloes than star-forming galaxies at fixed $M_*$., Comment: MNRAS accepted, 26 pages

In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at $z=4$ and 3, and before the last major merger, focusing on the formation of well-defined rotationally-supported disks, the mass-metallicity relation, the specific star formation rate, the gas metallicity gradients, and the non-axisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semi-analytic models at $z\sim$2. We see that almost all the hydro codes develop rotationally-supported structures at low redshifts. Most agree within 0.5 dex with the observed MZR at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the inter-code differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time-stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to $z=2$ and the auxiliary data will be made publicly available., Comment: Submitted to ApJ. Visit the AGORA Collaboration website (www.agorasimulations.org ) for more information

We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near $z=2$ for each code (hereafter ``$z\sim2$') at which the eight simulations are in the same stage in the target halo's merger history. After identifying the matched pairs of halos between the {\it CosmoRun} simulations and the DMO simulations, we discover that each {\it CosmoRun} halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at $z\sim2$, the number of satellite {\it galaxies} is significantly fewer than that of dark matter halos in all participating {\it AGORA} simulations, and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called ``missing satellite problem' is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in {\it AGORA} and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution ($\lesssim100$ proper pc at $z=2$). We also compare other properties such as the stellar mass$-$halo mass relation and the mass$-$metallicity relation. Our work highlights the value of comparison studies such as {\it AGORA}, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms., Comment: Accepted for publication to the ApJ, 19 pages, 10 figures. Visit the AGORA Collaboration website (www.agorasimulations.org) for more information. For summary video, please see https://www.youtube.com/watch?v=HShRJYvPe1k

We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the structure, composition, and phase dynamics of the CGM. We show properties such as metal distribution, ionization levels, and kinematics are effective tracers of the effects of the different code feedback and implementation methods, and as such they can be highly divergent between simulations. This is merely a fiducial set of models, against which we will in the future compare multiple feedback recipes for each code. Nevertheless, we find that the large parameter space these simulations establish can help disentangle the different variables that affect observable quantities in the CGM, e.g., showing that abundances for ions with higher ionization energy are more strongly determined by the simulation's metallicity, while abundances for ions with lower ionization energy are more strongly determined by the gas density and temperature., Comment: 32 pages, 17 figures. Published in ApJ

We present a novel method for identifying cosmic web filaments using the IllustrisTNG (TNG100) cosmological simulations and investigate the impact of filaments on galaxies. We compare the use of cosmic density field estimates from the Delaunay Tessellation Field Estimator (DTFE) and the Monte Carlo Physarum Machine (MCPM), which is inspired by the slime mold organism, in the DisPerSE structure identification framework. The MCPM-based reconstruction identifies filaments with higher fidelity, finding more low-prominence/diffuse filaments and better tracing the true underlying matter distribution than the DTFE-based reconstruction. Using our new filament catalogs, we find that most galaxies are located within 1.5-2.5 Mpc of a filamentary spine, with little change in the median specific star formation rate and the median galactic gas fraction with distance to the nearest filament. Instead, we introduce the filament line density, Sigma_fil(MCPM), as the total MCPM overdensity per unit length of a local filament segment, and find that this parameter is a superior predictor of galactic gas supply and quenching. Our results indicate that most galaxies are quenched and gas-poor near high-line density filaments at z<=1. At z=0, quenching in log(M*/Msun)>10.5 galaxies is mainly driven by mass, while lower-mass galaxies are significantly affected by the filament line density. In high-line density filaments, satellites are strongly quenched, whereas centrals have reduced star formation, but not gas fraction, at z<=0.5. We discuss the prospect of applying our new filament identification method to galaxy surveys with SDSS, DESI, Subaru PFS, etc. to elucidate the effect of large-scale structure on galaxy formation., Comment: Accepted for publication in ApJ. Data available at https://doi.org/10.5281/zenodo.10929036

The 3D geometry of high-redshift galaxies remains poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in JWST-CEERS observations with $\log M_*/M_{\odot}=9.0-10.5$ at $z=0.5-8.0$. We reproduce previous results from HST-CANDELS in a fraction of the computing time and constrain the mean ellipticity, triaxiality, size and covariances with samples as small as $\sim50$ galaxies. We find high 3D ellipticities for all mass-redshift bins suggesting oblate (disky) or prolate (elongated) geometries. We break that degeneracy by constraining the mean triaxiality to be $\sim1$ for $\log M_*/M_{\odot}=9.0-9.5$ dwarfs at $z>1$ (favoring the prolate scenario), with significantly lower triaxialities for higher masses and lower redshifts indicating the emergence of disks. The prolate population traces out a ``banana'' in the projected $b/a-\log a$ diagram with an excess of low $b/a$, large $\log a$ galaxies. The dwarf prolate fraction rises from $\sim25\%$ at $z=0.5-1.0$ to $\sim50-80\%$ at $z=3-8$. If these are disks, they cannot be axisymmetric but instead must be unusually oval (triaxial) unlike local circular disks. We simultaneously constrain the 3D size-mass relation and its dependence on 3D geometry. High-probability prolate and oblate candidates show remarkably similar S\'ersic indices ($n\sim1$), non-parametric morphological properties and specific star formation rates. Both tend to be visually classified as disks or irregular but edge-on oblate candidates show more dust attenuation. We discuss selection effects, follow-up prospects and theoretical implications., Comment: Accepted version to appear in ApJ, main body is 36 pages of which ~half are full-page figures

Giant star-forming clumps are a prominent feature of star-forming galaxies (SFGs) and contain important clues on galaxy formation and evolution. However, basic demographics of clumps and their host galaxies remain uncertain. Using the HST/WFC3 F275W images from the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (UVCANDELS), we detect and analyze giant star-forming clumps in galaxies at 0.5 $\leq$ z $\leq$ 1, connecting two epochs when clumps are common (at cosmic high-noon, z $\sim$ 2) and rare (in the local universe). We construct a clump sample whose rest-frame 1600 {\AA} luminosity is 3 times higher than the most luminous local HII regions (M$_{UV} \leq -$16 AB). In our sample, 35 $\pm$ 3$\%$ of low-mass galaxies (log[M$_{*}$/M$_{\odot}$] $<$ 10) are clumpy (i.e., containing at least one off-center clump). This fraction changes to 22 $\pm$ 3$\%$ and 22 $\pm$ 4$\%$ for intermediate (10 $\leq$ log[M$_{*}$/M$_{\odot}$] $\leq$ 10.5) and high-mass (log[M$_{*}$/M$_{\odot}$] $>$ 10.5) galaxies in agreement with previous studies. When compared to similar-mass non-clumpy SFGs, low- and intermediate-mass clumpy SFGs tend to have higher SFRs and bluer rest-frame U-V colors, while high-mass clumpy SFGs tend to be larger than non-clumpy SFGs. However, clumpy and non-clumpy SFGs have similar S\'ersic index, indicating a similar underlying density profile. Furthermore, we investigate how UV luminosity of star-forming regions correlates with the physical properties of host galaxies. On average, more luminous star-forming regions reside in more luminous, smaller, and/or higher-specific SFR galaxies and are found closer to their hosts' galactic center., Comment: 21 pages, 13 figures, accepted for publication in ApJ

We investigate how cosmic web structures affect galaxy quenching in the IllustrisTNG (TNG100) cosmological simulations by reconstructing the cosmic web within each snapshot using the DisPerSE framework. We measure the comoving distance from each galaxy with stellar mass $\log(M_{\ast}/\mathrm{M}_{\odot}) \geq 8$ to the nearest node ($d_{\mathrm{node}}$) and the nearest filament spine ($d_{\mathrm{fil}}$) to study the dependence of both median specific star formation rate () and median gas fraction (<$f_{\mathrm{gas}}$>) on these distances. We find that the of galaxies is only dependent on cosmic web environment at $z<2$, with the dependence increasing with time. At $z\leq0.5$, $8 \leq \log(M_{\ast}/\mathrm{M}_{\odot}) < 9$ galaxies are quenched at $d_{\mathrm{node}}\lesssim1$~Mpc, and have significantly-suppressed star formation at $d_{\mathrm{fil}}\lesssim1$~Mpc, trends driven mostly by satellite galaxies. At $z\leq1$, in contrast to the monotonic drop in of $\log(M_{\ast}/\mathrm{M}_{\odot}) <10$ galaxies with decreasing $d_{\mathrm{node}}$ and $d_{\mathrm{fil}}$, $\log(M_{\ast}/\mathrm{M}_{\odot}) \geq 10$ galaxies - both centrals and satellites - experience an upturn in at $d_{\mathrm{node}}\lesssim0.2$~Mpc. Much of this cosmic web dependence of star formation activity can be explained by an evolution in $$. Our results suggest that in the past $\sim$10 Gyr, low-mass satellites are quenched by rapid gas stripping in dense environments near nodes and gradual gas starvation in intermediate-density environments near filaments, while at earlier times cosmic web structures efficiently channeled cold gas into most galaxies. State-of-the-art ongoing spectroscopic surveys such as SDSS and DESI, as well as those planned with the Subaru Prime Focus Spectrograph, JWST and Roman, are required to test our predictions against observations., Comment: Accepted for publication in ApJ. 5 Figures

We utilize high-resolution cosmological simulations to reveal that high-redshift galaxies tend to undergo a robust `wet compaction' event when near a `golden' stellar mass of $\sim 10^{10} M_{\odot}$. This is a gaseous shrinkage to a compact star-forming phase, a `blue nugget' (BN), followed by central quenching of star formation to a compact passive stellar bulge, a `red nugget' (RN), and a buildup of an extended gaseous disc and ring. Such nuggets are observed at cosmic noon and seed today's early-type galaxies. The compaction is triggered by a drastic loss of angular momentum due to, e.g., wet mergers, counter-rotating cold streams, or violent disc instability. The BN phase marks drastic transitions in the galaxy structural, compositional and kinematic properties. The transitions are from star-forming to quenched inside-out, from diffuse to compact with an extended disc-ring and a stellar envelope, from dark matter to baryon central dominance, from prolate to oblate stellar shape, from pressure to rotation support, from low to high metallicity, and from supernova to AGN feedback. The central black hole growth, first suppressed by supernova feedback when below the golden mass, is boosted by the compaction, and the black hole keeps growing once the halo is massive enough to lock in the supernova ejecta., Comment: 33 pages, 26 figures in the main body (49 pages, 45 figures including appendix). Corresponds to accepted MNRAS version, minor revisions

The evolution of star-forming galaxies at high redshifts is very sensitive to the strength and nature of stellar feedback. Using two sets of cosmological, zoom-in simulations from the VELA suite, we compare the effects of two different models of feedback: with and without kinetic feedback from the expansion of supernovae shells and stellar winds. At a fixed halo mass and redshift, the stellar mass is reduced by a factor of 1-3 in the models with stronger feedback, so the stellar-mass-halo-mass relation is in better agreement with abundance matching results. On the other hand, the three-dimensional shape of low-mass galaxies is elongated along a major axis in both models. At a fixed stellar mass, Ms<10^10 Msun, galaxies are more elongated in the strong-feedback case. More massive, star-forming discs with high surface densities form giant clumps. However, the population of round, compact, old (age_c > 300 Myr), quenched, stellar (or gas-poor) clumps is absent in the model with strong feedback. On the other hand, giant star-forming clumps with intermediate ages (age_c= 100 - 300 Myr) can survive for several disc dynamical times, independently of feedback strength. The evolution through compaction followed by quenching in the plane of central surface density and specific star-formation rate is similar under the two feedback models., Comment: 14 pages, 12 figures, accepted version at MNRAS

Most studies of highly ionized plasmas have historically assumed ions are either in photoionization equilibrium, PIE, or collisional ionization equilibrium, CIE, sometimes including multiple phases with different relevant mechanisms. Simulation analysis packages, on the other hand, tend to use precomputed ion fraction tables which include both mechanisms, among others. Focusing on the low-density, high temperature phase space likely to be most relevant in the circumgalactic medium, in this work we show that most ions can be classified as 'PI' or 'CI' on an ion-by-ion basis. This means that for a cloud at a particular point in phase space, some ions will be created primarily by PI and others by CI, with other mechanisms playing only very minor roles. Specifically, we show that ions are generally CI if the thermal energy per particle is greater than $\sim6$\% of their ionization energy, and PI otherwise. We analyse the accuracy of this ansatz compared to usual PIE/CIE calculations, and show the surprisingly minor dependence of this conclusion on redshift and ionizing background., Comment: 12 pages, accepted by MNRAS

We study the survival versus disruption of the giant clumps in high-redshift disc galaxies, short-lived (S) versus long-lived (L) clumps and two L sub-types, via analytic modeling tested against simulations. We develop a criterion for clump survival, with or without their gas, based on a predictive survivability parameter $S$. It compares the energy sources by supernova feedback and gravitational contraction to the clump binding energy and losses by outflows and turbulence dissipation. The clump properties are derived from Toomre instability, approaching virial and Jeans equilibrium, and the supernova energy deposit is based on an up-to-date bubble analysis. For moderate feedback levels, we find that L clumps exist with circular velocities $\sim\!50\, km\, s^{-1}$ and masses $\geq\!10^8\, M_\odot$. They are likely in galaxies with circular velocities $\geq \!200\, km\, s^{-1}$, consistent at $z \sim 2$ with the favored stellar mass for discs, $\geq\!10^{9.3}\, M_\odot$. L clumps favor disc gas fractions $\geq\!0.3$, low-mass bulges and redshifts $z\!\sim\! 2$. The likelihood of L clumps is reduced if the feedback is more ejective, e.g., if the supernovae are optimally clustered, if radiative feedback is very strong, if the stellar initial mass function is top-heavy, or if the star-formation-rate efficiency is particularly high. A sub-type of L clumps (LS), which lose their gas in several free-fall times but retain bound stellar components, may be explained by a smaller contraction factor and stronger external gravitational effects, where clump mergers increase the SFR efficiency. The more massive L clumps (LL) retain most of their baryons for tens of free-fall times with a roughly constant star-formation rate., Comment: 26 pages, 14 figures

The stochastic order redshift technique (SORT) is a simple, efficient, and robust method to improve cosmological redshift measurements. The method relies upon having a small ($\sim$10 per cent) reference sample of high-quality redshifts. Within pencil-beam-like sub-volumes surrounding each galaxy, we use the precise dN/d$z$ distribution of the reference sample to recover new redshifts and assign them one-to-one to galaxies such that the original rank order of redshifts is preserved. Preserving the rank order is motivated by the fact that random variables drawn from Gaussian probability density functions with different means but equal standard deviations satisfy stochastic ordering. The process is repeated for sub-volumes surrounding each galaxy in the survey. This results in every galaxy with an uncertain redshift being assigned multiple "recovered" redshifts from which a new redshift estimate can be determined. An earlier paper applied SORT to a mock Sloan Digital Sky Survey at $z \lesssim$ 0.2 and accurately recovered the two-point correlation function on scales $\gtrsim$4 $h^{-1}$Mpc. In this paper, we test the performance of SORT in surveys spanning the redshift range 0.75$

We address the nature of the giant clumps in high-z galaxies that undergo Violent Disc Instability, attempting to distinguish between long-lived clumps that migrate inward and short-lived clumps that disrupt by feedback. We study the evolution of clumps as they migrate through the disc using an analytic model tested by simulations and confront theory with CANDELS observations. The clump ``bathtub" model, which considers gas and stellar gain and loss, is characterized by four parameters: the accretion efficiency, the star-formation-rate (SFR) efficiency, and the outflow mass-loading factors for gas and stars. The relevant timescales are all comparable to the migration time, two-three orbital times. A clump differs from a galaxy by the internal dependence of the accretion rate on the varying clump mass. The analytic solution, involving exponential growing and decaying modes, reveals a main evolution phase during the migration, where the SFR and gas mass are constant and the stellar mass is rising linearly with time. This makes the inverse of the specific SFR an observable proxy for clump age. Later, the masses and SFR approach an exponential growth with a constant specific SFR, but this phase is hypothetical as the clump disappears in the galaxy center. The model matches simulations with different, moderate feedback, both in isolated and cosmological settings. The observed clumps agree with our theoretical predictions, indicating that the massive clumps are long-lived and migrating. A non-trivial challenge is to model feedback that is non-disruptive in massive clumps but suppresses SFR to match the galactic stellar-to-halo mass ratio., Comment: 27 pages, 10 figures

We present a suite of high-resolution cosmological zoom-in simulations to $z=4$ of a $10^{12}\,{\rm M}_{\odot}$ halo at $z=0$, obtained using seven contemporary astrophysical simulation codes widely used in the numerical galaxy formation community. Physics prescriptions for gas cooling, heating, and star formation, are similar to the ones used in our previous {\it AGORA} disk comparison but now account for the effects of cosmological processes. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. After all the participating code groups successfully completed the calibration steps, we reach a suite of cosmological simulations with similar mass assembly histories down to $z=4$. With numerical accuracy that resolves the internal structure of a target halo, we find that the codes overall agree well with one another in e.g., gas and stellar properties, but also show differences in e.g., circumgalactic medium properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, the results of high-resolution cosmological zoom-in simulations can be robust and reproducible. New code groups are invited to join this comparison by generating equivalent models by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the simulations presented here will be in forthcoming reports from our Collaboration., Comment: Accepted for publication to the ApJ. The simulation snapshots used in this paper will be made publicly available in a user-friendly interface. The announcement of the data release will be posted on the AGORA Project webpage (www.agorasimulations.org). Early access to the data can be provided upon request to the authors, or to the project coordinator (sroca01@ucm.es, santacruzgalaxy@gmail.com)

Observed rotation curves in star-forming galaxies indicate a puzzling dearth of dark matter in extended flat cores within haloes of mass $\geq\! 10^{12}M_\odot$ at $z\!\sim\! 2$. This is not reproduced by current cosmological simulations, and supernova-driven outflows are not effective in such massive haloes. We address a hybrid scenario where post-compaction merging satellites heat up the dark-matter cusps by dynamical friction, allowing AGN-driven outflows to generate cores. Using analytic and semi-analytic models (SatGen), we estimate the dynamical-friction heating as a function of satellite compactness for a cosmological sequence of mergers. Cosmological simulations (VELA) demonstrate that satellites of initial virial masses $>\!10^{11.3}M_\odot$, that undergo wet compactions, become sufficiently compact for significant heating. Constituting a major fraction of the accretion onto haloes $\geq\!10^{12}M_\odot$, these satellites heat-up the cusps in half a virial time at $z\!\sim\! 2$. Using a model for outflow-driven core formation (CuspCore), we demonstrate that the heated dark-matter cusps develop extended cores in response to removal of half the gas mass, while the more compact stellar systems remain intact. The mergers keep the dark matter hot, while the gas supply, fresh and recycled, is sufficient for the AGN outflows. AGN indeed become effective in haloes $\geq\!10^{12}M_\odot$, where the black-hole growth is no longer suppressed by supernovae and its compaction-driven rapid growth is maintained by a hot CGM. For simulations to reproduce the dynamical-friction effects, they should resolve the compaction of the massive satellites and avoid artificial tidal disruption. AGN feedback could be boosted by clumpy black-hole accretion and clumpy response to AGN., Comment: 25+10 pages, 10+8 figures

Fine-grained estimation of galaxy merger stages from observations is a key problem useful for validation of our current theoretical understanding of galaxy formation. To this end, we demonstrate a CNN-based regression model that is able to predict, for the first time, using a single image, the merger stage relative to the first perigee passage with a median error of 38.3 million years (Myrs) over a period of 400 Myrs. This model uses no specific dynamical modeling and learns only from simulated merger events. We show that our model provides reasonable estimates on real observations, approximately matching prior estimates provided by detailed dynamical modeling. We provide a preliminary interpretability analysis of our models, and demonstrate first steps toward calibrated uncertainty estimation., Comment: Third Workshop on Machine Learning and the Physical Sciences (NeurIPS 2020), Vancouver, Canada

We present mock catalogs created to support the interpretation of the CANDELS survey. We extract halos along past lightcones from the Bolshoi Planck dissipationless N-body simulations and populate these halos with galaxies using two different independently developed semi-analytic models of galaxy formation and the empirical model UniverseMachine. Our mock catalogs have geometries that encompass the footprints of observations associated with the five CANDELS fields. In order to allow field-to-field variance to be explored, we have created eight realizations of each field. In this paper, we present comparisons with observable global galaxy properties, including counts in observed frame bands, luminosity functions, color-magnitude distributions and color-color distributions. We additionally present comparisons with physical galaxy parameters derived from SED fitting for the CANDELS observations, such as stellar masses and star formation rates. We find relatively good agreement between the model predictions and CANDELS observations for luminosity and stellar mass functions. We find poorer agreement for colors and star formation rate distributions. All of the mock lightcones as well as curated "theory friendly" versions of the observational CANDELS catalogs are made available through a web-based data hub., Comment: 22 pages, 16 figures. To appear in MNRAS

The diffuse extragalactic background light (EBL) is formed by ultraviolet (UV), optical, and infrared (IR) photons mainly produced by star formation processes over the history of the Universe, and contains essential information about galaxy evolution and cosmology. Here, we present a new determination of the evolving EBL spectral energy distribution using a novel approach purely based on galaxy data aiming to reduce current uncertainties on the higher redshifts and IR intensities. Our calculations use multiwavelength observations from the UV to the far-IR of a sample of approximately 150,000 galaxies detected up to $z\sim 6$ in the five fields of the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) from the Hubble Space Telescope. This is one of the most comprehensive and deepest multi-wavelength galaxy datasets ever obtained. These unprecedented resources allow us to derive the overall EBL evolution up to $z\sim 6$ and its uncertainties. Our results agree with cosmic observables estimated from galaxy surveys and $\gamma$-ray attenuation such as monochromatic luminosity densities, including those in the far-IR, and star formation rate densities, also at the highest redshits. Optical depths from our EBL approximation, which will be robust at high redshifts and for $\gamma$ rays up to tens of TeV, will be reported in a companion paper., Comment: 17 pages, 13 figures, 4 tables; accepted for publication in MNRAS on August 6, 2021. Online material available at https://www.ucm.es/blazars/ebl

We use deep learning to explore the nature of observed giant clumps in high-redshift disc galaxies, based on their identification and classification in cosmological simulations. Simulated clumps are detected using the 3D gas and stellar densities in the VELA zoom-in cosmological simulation suite, with $\sim \!\! 25\ \!\rm{pc}$ maximum resolution, targeting main sequence galaxies at $1\!<\!z\!<\!3$. The clumps are classified as long-lived clumps (LLCs) or short-lived clumps (SLCs) based on their longevity in the simulations. We then train neural networks to detect and classify the simulated clumps in mock, multi-color, dusty and noisy HST-like images. The clumps are detected using an encoder-decoder convolutional neural network (CNN), and are classified according to their longevity using a vanilla CNN. Tests using the simulations show our detector and classifier to be $\sim80\%$ complete and $\sim80\%$ pure for clumps more massive than $\sim10^{7.5}\rm{M_\odot}$. When applied to observed galaxies in the CANDELS/GOODS S+N fields, we find both types of clumps to appear in similar abundances in the simulations and the observations. LLCs are, on average, more massive than SLCs by $\sim 0.5\ \rm{dex}$, and they dominate the clump population above $M_{\rm c}\gtrsim 10^{7.6}\ \rm{M_\odot}$. LLCs tend to be found closer to the galactic centre, indicating clump migration to the centre or preferential formation at smaller radii. The LLCs are found to reside in high mass galaxies, indicating better clump survivability under supernova feedback there, due to clumps being more massive in these galaxies. We find the clump properties in the simulations and the observations to agree within a factor of $\sim\! 2$., Comment: 17 pages, 13 figures. Accepted for publication in MNRAS

The thermal evolution of rocky planets on geological timescales (Gyr) depends on the heat input from the long-lived radiogenic elements potassium, thorium, and uranium. Concentrations of the latter two in rocky planet mantles are likely to vary by up to an order of magnitude between different planetary systems because Th and U, like other heavy r-process elements, are produced by rare stellar processes. Here we discuss the effects of these variations on the thermal evolution of an Earth-size planet, using a 1D parameterized convection model. Assuming Th and U abundances consistent with geochemical models of the Bulk Silicate Earth based on chondritic meteorites, we find that Earth had just enough radiogenic heating to maintain a persistent dynamo. According to this model, Earth-like planets of stars with higher abundances of heavy r-process elements, indicated by the relative abundance of europium in their spectra, are likely to have lacked a dynamo for a significant fraction of their lifetimes, with potentially negative consequences for hosting a biosphere. Because the qualitative outcomes of our 1D model are strongly dependent on the treatment of viscosity, further investigations using fully 3D convection models are desirable.

We provide prescriptions to evaluate the dynamical mass ($M_{\rm dyn}$) of galaxies from kinematic measurements of stars or gas using analytic considerations and the VELA suite of cosmological zoom-in simulations at $z=1-5$. We find that Jeans or hydrostatic equilibrium is approximately valid for galaxies of stellar masses above $M_\star \!\sim\! 10^{9.5}M_\odot$ out to $5$ effective radii ($R_e$). When both measurements of the rotation velocity $v_\phi$ and of the radial velocity dispersion $\sigma_r$ are available, the dynamical mass $M_{\rm dyn} \!\simeq\! G^{-1} V_c^2 r$ can be evaluated from the Jeans equation $V_c^2= v_\phi^2 + \alpha \sigma_r^2$ assuming cylindrical symmetry and a constant, isotropic $\sigma_r$. For spheroids, $\alpha$ is inversely proportional to the S\'ersic index $n$ and $\alpha \simeq 2.5$ within $R_e$ for the simulated galaxies. The prediction for a self-gravitating exponential disc, $\alpha = 3.36(r/R_e)$, is invalid in the simulations, where the dominant spheroid causes a weaker gradient from $\alpha \!\simeq\! 1$ at $R_e$ to 4 at $5R_e$. The correction in $\alpha$ for the stars due to the gradient in $\sigma_r(r)$ is roughly balanced by the effect of the aspherical potential, while the effect of anisotropy is negligible. When only the effective projected velocity dispersion $\sigma_l$ is available, the dynamical mass can be evaluated as $M_{\rm dyn} = K G^{-1} R_e \sigma_l^2$, where the virial factor $K$ is derived from $\alpha$ given the inclination and $v_\phi/\sigma_r$. We find that the standard value $K=5$ is approximately valid only when averaged over inclinations and for compact and thick discs, as it ranges from 4.5 to above 10 between edge-on and face-on projections., Comment: 16 pages, 11 figures, published in MNRAS

Galaxy morphology and its evolution over the cosmic epoch hold important clues for understanding the regulation of star formation (SF). However, studying the relationship between morphology and SF has been hindered by the availability of consistent data at different redshifts. Our sample, combining CANDELS (0.8 < z < 2.5) and the GALEX-SDSS-WISE Legacy Catalog (GSWLC; z ~ 0), has physical parameters derived using consistent SED fitting with flexible dust attenuation laws. We adopt visual classifications from Kartaltepe et al. 2015 and expand them to z ~ 0 using SDSS images matching the physical resolution of CANDELS rest-frame optical images and deep FUV GALEX images matching the physical resolution of the CANDELS rest-frame FUV images. Our main finding is that disks with SF clumps at z ~ 0 make a similar fraction (~15%) of star-forming galaxies as at z ~ 2. The clumpy disk contribution to the SF budget peaks at z ~ 1, rather than z ~ 2, suggesting that the principal epoch of disk assembly continues to lower redshifts. Star-forming spheroids ("blue nuggets"), though less centrally concentrated than quenched spheroids, contribute significantly (~15%) to the SF budget at z ~ 1-2, suggesting that compaction precedes quenching. Among green valley and quiescent galaxies, the pure spheroid fraction drops since z ~ 1, whereas spheroids with disks (S0-like) become dominant. Mergers at or nearing coalescence are enhanced in SFR relative to the main sequence at all redshifts by a factor of ~2, but contribute $\lesssim$5% to the SF budget, with their contribution remaining small above the main sequence., Comment: 26 pages, 13 figures, 3 tables, accepted for publication in ApJ

We analyse the distribution and origin of OVI in the Circumgalactic Medium (CGM) of dark-matter haloes of $\sim 10^{12}$ M$_\odot$ at $z\sim1$ in the VELA cosmological zoom-in simulations. We find that the OVI in the inflowing cold streams is primarily photoionized, while in the bulk volume it is primarily collisionally ionized. The photoionized component dominates the observed column density at large impact parameters ($\gtrsim 0.3 R_{\rm vir}$), while the collisionally ionized component dominates closer in. We find that most of the collisional OVI, by mass, resides in the relatively thin boundaries of the photoionized streams. We discuss how the results are in agreement with analytic predictions of stream and boundary properties, and their compatibility with observations. This allows us to predict the profiles of OVI and other ions in future CGM observations and provides a toy model for interpreting them., Comment: 22 pages, 16 figures, submitted to MNRAS

An analytical approach is proposed to study the evolution of the star-forming galaxy (SFG) main sequence (MS) and the fraction of dust-obscured SF up to $z\sim4$. Far-ultraviolet (FUV) and infrared (IR) star formation rates, SFRs, are described as conditional probability functions of $M_{\ast}$. We convolve them with the galaxy stellar mass function (GSMF) of SFGs to derive the FUV and IR LFs. The 2 SF modes formalism is used to describe starburst galaxies. By fitting observed FUV and IR LFs, the parametrization of SFR$_{\rm FUV}-M_{\ast}$ and SFR$_{\rm IR}-M_{\ast}$ are constrained. Our derived SFR$_{\rm FUV+IR}-M_{\ast}$ reproduces the evolution of the MS as compared to other observational inferences. At any redshift, we find that the sSFR$_{\rm FUV+IR}-M_{\ast}$ relation for MS SFGs approaches to a power law at the high-mass end. At lower masses, it bends and eventually the slope sign changes from negative to positive at very low masses. At $z\sim0$, this change of sign is at $M_{\ast}\sim5\times10^{8}{\rm M}_{\odot}$ close to dust-obscured SF regime, $M_{\ast}\sim6\times10^{8}{\rm M}_{\odot}$. The slope sign change is related to the knee of the FUV LF. Our derived dust-obscured fractions agree with previous determinations at $0\leq z\leq2.5$. Dust-obscured fractions depend strongly on mass with almost no dependence with redshift at $z\gtrsim1.2$. At $z\lesssim0.75$ high-mass galaxies become more "transparent" compared to their high redshift counterparts. On the opposite, low- and intermediate-mass galaxies have become more obscured by dust. The joint evolution of the GSMF and the FUV and IR LFs is a promising approach to study mass growth and dust formation/destruction mechanisms., Comment: Submitted to ApJ, 20 pages, 12 figures. Comments are welcome. This is an improved version after the first Referee report. Results remained unchanged

Are there parallels between the furthest reaches of our universe, and the foundations of thought, awareness, perception, and emotion? What are the connections between the webs and structures that define both? What are the differences? "As Above As Below" was an exhibition that examined these questions. It consisted of six artworks, each of them the product of a collaboration that included at least one artist, astrophysicist, and neuroscientist. The installations explored new parallels between intergalactic and neuronal networks through media such as digital projection, virtual reality, and interactive multimedia, and served to illustrate diverse collaboration practices and ways to communicate across very different fields., Comment: In SciArt Magazine, https://www.sciartmagazine.com/collaboration-as-above-as-below.html ; exhibition video https://youtu.be/RM3gDjb3phU ; panel discussion https://youtu.be/qNhtyvPz8yw ; project website https://asaboveasbelow.com/

Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perform a stringent test of the fine morphological structure of galaxies coming from the Illustris and IllustrisTNG (TNG100 and TNG50) simulations against observations from a subsample of the Sloan Digital Sky Survey. Our framework is based on PixelCNN, an autoregressive model for image generation with an explicit likelihood. We adopt a strategy that combines the output of two PixelCNN networks in a metric that isolates the fine morphological details of galaxies from the sky background. We are able to \emph{quantitatively} identify the improvements of IllustrisTNG, particularly in the high-resolution TNG50 run, over the original Illustris. However, we find that the fine details of galaxy structure are still different between observed and simulated galaxies. This difference is driven by small, more spheroidal, and quenched galaxies which are globally less accurate regardless of resolution and which have experienced little improvement between the three simulations explored. We speculate that this disagreement, that is less severe for quenched disky galaxies, may stem from a still too coarse numerical resolution, which struggles to properly capture the inner, dense regions of quenched spheroidal galaxies., Comment: Accepted for publication in MNRAS. Comments welcome

LCDM cosmological models with Early Dark Energy (EDE) have been proposed to resolve tensions between the Hubble constant H0 = 100h km/s/Mpc measured locally, giving h ~ 0.73, and H0 deduced from Planck cosmic microwave background (CMB) and other early universe measurements plus LCDM, giving h ~ 0.67. EDE models do this by adding a scalar field that temporarily adds dark energy equal to about 10% of the cosmological energy density at the end of the radiation-dominated era at redshift z ~ 3500. Here we compare linear and nonlinear predictions of a Planck-normalized LCDM model including EDE giving h = 0.728 with those of standard Planck-normalized LCDM with h = 0.678. We find that nonlinear evolution reduces the differences between power spectra of fluctuations at low redshifts. As a result, at z = 0 the halo mass functions on galactic scales are nearly the same, with differences only 1-2%. However, the differences dramatically increase at high redshifts. The EDE model predicts 50% more massive clusters at z = 1 and twice more galaxy-mass halos at z = 4. Even greater increases in abundances of galaxy-mass halos at higher redshifts may make it easier to reionize the universe with EDE. Predicted galaxy abundances and clustering will soon be tested by JWST observations. Positions of baryonic acoustic oscillations (BAOs) and correlation functions differ by about 2% between the models -- an effect that is not washed out by nonlinearities. Both standard LCDM and the EDE model studied here agree well with presently available acoustic-scale observations, but DESI and Euclid measurements will provide stringent new tests., Comment: 13 pages, 12 figures, submitted to MNRAS

Using analytic modeling and simulations, we address the origin of an abundance of star-forming, clumpy, extended gas rings about massive central bodies in massive galaxies at $z \!<\! 4$. Rings form by high-angular-momentum streams and survive in galaxies of $M_{\rm star} \!>\! 10^{9.5-10} M_\odot$ where merger-driven spin flips and supernova feedback are ineffective. The rings survive after events of compaction to central nuggets. Ring longevity was unexpected based on inward mass transport driven by torques from violent disc instability. However, evaluating the torques from a tightly wound spiral structure, we find that the timescale for transport per orbital time is long and $\propto\! \delta_{\rm d}^{-3}$, with $\delta_{\rm d}$ the cold-to-total mass ratio interior to the ring. A long-lived ring forms when the ring transport is slower than its replenishment by accretion and the interior depletion by SFR, both valid for $\delta_{\rm d} \!<\! 0.3$. The central mass that lowers $\delta_{\rm d}$ is a compaction-driven bulge and/or dark matter, aided by the lower gas fraction at $z \!<\! 4$, provided that it is not too low. The ring is Toomre unstable for clump and star formation. The high-$z$ dynamic rings are not likely to arise form secular resonances or collisions. AGN feedback is not expected to affect the rings. Mock images of simulated rings through dust indicate qualitative consistency with observed rings about bulges in massive $z\!\sim\!0.5\!-\!3$ galaxies, in $H_{\alpha}$ and deep HST imaging. ALMA mock images indicate that $z\!\sim\!0.5\!-\!1$ rings should be detectable. We quote expected observable properties of rings and their central nuggets., Comment: 40 pages, 30 figures

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release., Comment: AGORA links: http://www.AGORAsimulations.org/ and http://sites.google.com/site/santacruzcomparisonproject/blogs/quicklinks/

We predict, analytically and by simulations, that gas discs tend to survive only in haloes above a threshold mass $\sim 2 \times 10^{11} M_\odot$ (stellar mass $\sim 10^9 M_\odot$), with only a weak redshift dependence. At lower masses, the disc spins typically flip in less than an orbital time due to mergers associated with a change in the pattern of the feeding cosmic-web streams. This threshold arises from the halo merger rate when accounting for the mass dependence of the ratio of galactic baryons and halo mass. Above the threshold, wet compactions lead to massive central nuggets that allow the longevity of extended clumpy gas rings. Supernova feedback has a major role in disrupting discs below the critical mass, by driving the stellar-to-halo mass ratio that affects the merger rate, by stirring up turbulence and suppressing high-angular-momentum gas supply, and by confining major compactions to the critical mass. Our predictions seem consistent with current observed fractions of gas discs, to be explored by future observations that will resolve galaxies below $10^9 M_\odot$ at high redshifts, e.g. by JWST., Comment: 23 pages, 11 figures

We use bulge-type classifications of 809 representative SDSS galaxies by Gadotti (2009) to classify a large sample of galaxies into real bulges (classical or elliptical) and pseudobulges using Random Forest. We use structural and stellar population predictors that can easily be measured without image decomposition. Multiple parameters such as the central mass density with 1 kpc, concentration index, S\'{e}rsic index and velocity dispersion do result in accurate bulge classifications when combined together. We classify $\sim 44,500$ face-on galaxies above stellar mass of 10$^{10}$ M$_\odot$ and redshift $ 0.02 < z < 0.07$ into real bulges or pseudobulges with $93 \pm 2$\% accuracy. We show that $\sim 75 - 90\%$ of AGNs identified by the optical line ratio diagnostic are hosted by real bulges. The pseudobulge fraction significantly decreases with AGN signature as the line ratios change from indicating pure star formation ($\sim 54 \pm 4$ \%), to composite of star formation and AGN ($\sim 18 \pm 3$\%), and to AGN-dominated galaxies ($\sim 5 \pm 3$\%). Using the dust-corrected [\ion{O}{3}] luminosity as an AGN accretion indicator, and the stellar mass and radius as proxies for a black hole mass, we find that AGNs in real bulges have lower Eddington ratios than AGNs in pseudobulges. Real bulges have a wide range of AGN and star formation activities, although most of them are weak AGNs. For both bulge-types, their Eddington ratios are correlated with specific star formation rates (SSFR). Real bulges have lower specific accretion rate but higher AGN fraction than pseudobulges do at similar SSFRs., Comment: 22 pages, 14 figures, accepted by ApJ

As computational resolution of modern cosmological simulations reach ever so close to resolving individual star-forming clumps in a galaxy, a need for "resolution-appropriate" physics for a galaxy-scale simulation has never been greater. To this end, we introduce a self-consistent numerical framework that includes explicit treatments of feedback from star-forming molecular clouds (SFMCs) and massive black holes (MBHs). In addition to the thermal supernovae feedback from SFMC particles, photoionizing radiation from both SFMCs and MBHs is tracked through full 3-dimensional ray tracing. A mechanical feedback channel from MBHs is also considered. Using our framework, we perform a state-of-the-art cosmological simulation of a quasar-host galaxy at z~7.5 for ~25 Myrs with all relevant galactic components such as dark matter, gas, SFMCs, and an embedded MBH seed of ~> 1e6 Ms. We find that feedback from SFMCs and an accreting MBH suppresses runaway star formation locally in the galactic core region. Newly included radiation feedback from SFMCs, combined with feedback from the MBH, helps the MBH grow faster by retaining gas that eventually accretes on to the MBH. Our experiment demonstrates that previously undiscussed types of interplay between gas, SFMCs, and a MBH may hold important clues about the growth and feedback of quasars and their host galaxies in the high-redshift Universe., Comment: 16 pages, 10 figures, Accepted for publication in the Astrophysical Journal, High-resolution version of this article also available at http://www.jihoonkim.org/index/research.html#mbh

This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR., Comment: 21 pages, 9 figures, ApJ accepted

Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws., Comment: 45 pages, 20 figures, 2 tables, accepted for publication in ApJ

This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $\Delta\Sigma_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $\Delta\Sigma_1$ to the established bulge-type indicator $\Delta\langle\mu_e\rangle$ from Gadotti (2009) and show that classifying by $\Delta\Sigma_1$ agrees well with $\Delta\langle\mu_e\rangle$. $\Delta\Sigma_1$ requires no bulge-disk decomposition and can be measured on SDSS images out to $z = 0.07$. Bulge types using it are mapped onto twenty different structural and stellar-population properties for 12,000 SDSS central galaxies with masses 10.0 < log $M_*$/$M_{\odot}$ < 10.4. New trends emerge from this large sample. Structural parameters show fairly linear log-log relations vs. $\Delta\Sigma_1$ and $\Delta\langle\mu_e\rangle$ with only moderate scatter, while stellar-population parameters show a highly non-linear "elbow" in which specific star-formation rate remains roughly flat with increasing central density and then falls rapidly at the elbow, where galaxies begin to quench. P-bulges occupy the low-density end of the horizontal arm of the elbow and are universally star-forming, while C-bulges occupy the elbow and the vertical branch and exhibit a wide range of star-formation rates at fixed density. The non-linear relation between central density and star-formation rate has been seen before, but this mapping onto bulge class is new. The wide range of star-formation rates in C-bulges helps to explain why bulge classifications using different parameters have sometimes disagreed in the past. The elbow-shaped relation between density and stellar indices suggests that central structure and stellar-populations evolve at different rates as galaxies begin to quench., Comment: 23 pages, 17 figures, published in MNRAS

Theory suggests that there are two primary modes of accretion through which dark matter halos acquire the gas to form and fuel galaxies, hot and cold mode accretion. In cold mode accretion, gas streams along cosmic web filaments to the center of the halo, allowing for the efficient delivery of star-forming fuel. Recently, two QSO-illuminated HI Lyman alpha (Ly{\alpha}) emitting objects were reported to have properties of cold, rotating structures (Martin et al. 2015, Martin et al. 2016). However, the spatial and spectral resolution available was insufficient to constrain radial flows associated with connecting filaments. With the Keck Cosmic Web Imager (KCWI) we now have eight times the spatial resolution, permitting the detection of these in-spiraling flows. In order to detect these inflows, we introduce a suite of models which incorporate zonal radial flows, demonstrate their performance on a numerical simulation that exhibits coldflow accretion, and show that they are an excellent match to KCWI velocity maps of two Ly{\alpha} emitters observed around high-redshift quasars. These Multi-Filament Inflow models kinematically isolate zones of radial inflow that correspond to extended filamentary emission. The derived gas flux and inflow path is sufficient to fuel the inferred central galaxy star formation rate and angular momentum. Thus, our kinematic emission maps provide strong evidence for the inflow of gas from the cosmic web building galaxies at the peak of star formation., Comment: Accepted for publication in Nature Astronomy

The role of major mergers in galaxy evolution remains a key open question. Existing empirical merger identification methods use non-parametric and subjective visual classifications which can pose systematic challenges to constraining merger histories. As a first step towards overcoming these challenges, we develop and share publicly a new Python-based software tool that identifies and extracts the flux-wise and area-wise significant contiguous regions from the model-subtracted "residual" images produced by popular parametric light-profile fitting tools (e.g., GALFIT). Using Hubble Space Telescope ($HST$) $H$-band single-S\'ersic residual images of $17$ CANDELS galaxies, we demonstrate the tool's ability to measure the surface brightness and improve the qualitative identification of a variety of common residual features (disk structures, spiral substructures, plausible tidal features, and strong gravitational arcs). We test our method on synthetic $HST$ observations of a $z\sim 1.5$ major merger from the VELA hydrodynamic simulations. We extract $H$-band residual features corresponding to the birth, growth, and fading of tidal features during different stages and viewing orientations at CANDELS depths and resolution. We find that the extracted features at shallow depths have noisy visual appearance and are susceptible to viewing angle effects. For a VELA $z\sim 3$ major merger, we find that James Webb Space Telescope NIRCam observations can probe high-redshift tidal features with considerable advantage over existing $HST$ capabilities. Further quantitative analysis of plausible tidal features extracted with our new software hold promise for the robust identification of hallmark merger signatures and corresponding improvements to merger rate constraints., Comment: 18 pages, 11 Figures, and accepted for publication in MNRAS. The data products discussed in this paper along with the software are publicly available at https://github.com/AgentM-GEG/residual_feature_extraction

We address a simple model where the Kennicutt-Schmidt (KS) relation between the macroscopic densities of star-formation rate (SFR, $\rho_{\rm sfr}$) and gas ($n$) in galactic discs emerges from self-regulation of the SFR via supernova feedback. It arises from the physics of supernova bubbles, insensitive to the microscopic SFR recipe and not explicitly dependent on gravity. The key is that the filling factor of SFR-suppressed supernova bubbles self-regulates to a constant, $f\sim 0.5$. Expressing the bubble fading radius and time in terms of $n$, the filling factor is $f \propto S\,n^{-s}$ with $s\sim 1.5$, where $S$ is the supernova rate density. A constant $f$ thus refers to $\rho_{\rm sfr} \propto n^{1.5}$, with a density-independent SFR efficiency per free-fall time $\sim 0.01$. The self-regulation to $f \sim 0.5$ and the convergence to a KS relation independent of the local SFR recipe are demonstrated in cosmological and isolated-galaxy simulations using different codes and recipes. In parallel, the spherical analysis of bubble evolution is generalized to clustered supernovae, analytically and via simulations, yielding $s \simeq 1.5 \pm 0.5$. An analysis of photo-ionized bubbles about pre-supernova stars yields a range of KS slopes but the KS relation is dominated by the supernova bubbles. Superbubble blowouts may lead to an alternative self-regulation by outflows and recycling. While the model is over-simplified, its simplicity and validity in the simulations may argue that it captures the origin of the KS relation.

Hubble Space Telescope observations show that low-mass ($M_*=10^9-10^{10}M_{\odot}$) galaxies at high redshift ($z=1.0-2.5$) tend to be elongated (prolate) rather than disky (oblate) or spheroidal. This is explained in zoom-in cosmological hydrodynamical simulations by the fact that these galaxies are forming in cosmic web filaments where accretion happens preferentially along the direction of elongation. We ask whether the elongated morphology of these galaxies allows them to be used as effective tracers of cosmic web filaments at high redshift via their intrinsic alignments. Using mock lightcones and spectroscopically-confirmed galaxy pairs from the CANDELS survey, we test two types of alignments: (1) between the galaxy major axis and the direction to nearby galaxies of any mass, and (2) between the major axes of nearby pairs of low-mass, likely prolate, galaxies. The mock lightcones predict strong signals in 3D real space, 3D redshift space, and 2D projected redshift space for both types of alignments (assuming prolate galaxy orientations are the same as those of their host prolate halos), but we do not detect significant alignment signals in CANDELS observations. However, we show that spectroscopic redshifts have been obtained for only a small fraction of highly elongated galaxies, and accounting for spectroscopic incompleteness and redshift errors significantly degrades the 2D mock signal. This may partly explain the alignment discrepancy and highlights one of several avenues for future work., Comment: Re-submitted to MNRAS after minor revisions

The majority of massive star-forming galaxies at $z\sim2$ have velocity gradients suggestive of rotation, in addition to large amounts of disordered motions. In this paper, we demonstrate that it is challenging to distinguish the regular rotation of a disk galaxy from the orbital motions of merging galaxies with seeing-limited data. However, the merger fractions at $z\sim2$ are likely too low for this to have a large effect on measurements of disk fractions. To determine how often mergers pass for disks, we look to galaxy formation simulations. We analyze $\sim$24000 synthetic images and kinematic maps of 31 high-resolution simulations of isolated galaxies and mergers at $z\sim2$. We determine if the synthetic observations pass criteria commonly used to identify disk galaxies, and whether the results are consistent with their intrinsic dynamical states. Galaxies that are intrinsically mergers pass the disk criteria for anywhere from 0 to 100$\%$ of sightlines. The exact percentage depends strongly on the specific disk criteria adopted, and weakly on the separation of the merging galaxies. Therefore, one cannot tell with certainty whether observations of an individual galaxy indicate a merger or a disk. To estimate the fraction of mergers passing as disks in current kinematics samples, we combine the probability that a merger will pass as a disk with theoretical merger fractions from a cosmological simulation. Taking the latter at face-value, the observed disk fractions are overestimated by small amounts: at most by $5\%$ at high stellar mass ($10^{10-11}$ M$_{\odot}$) and $15\%$ at low stellar mass ($10^{9-10}$ M$_{\odot}$)., Comment: 15 pages, 10 figures, accepted for publication in ApJ

We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component S\'ersic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxies against ~ 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc vs. bulge dominated) and the star formation activity (star-forming vs. quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our findings pose a challenge to models envisioning multiple channels for massive bulge growth, such as disc instabilities and mergers., Comment: 20 pages, 15 figures, Submitted to MNRAS

We model the projected b/a-log a distributions of CANDELS main sequence star-forming galaxies, where a (b) is the semi-major (semi-minor) axis of the galaxy images. We find that smaller-a galaxies are rounder at all stellar masses M and redshifts, so we include a when analyzing b/a distributions. Approximating intrinsic shapes of the galaxies as triaxial ellipsoids and assuming a multivariate normal distribution of galaxy size and two shape parameters, we construct their intrinsic shape and size distributions to obtain the fractions of prolate, oblate and spheroidal galaxies in each redshift and mass bin. We find that galaxies tend to be prolate at low m and high redshifts, and oblate at high M and low redshifts, qualitatively consistent with van der Wel et al. (2014), implying that galaxies tend to evolve from prolate to oblate. These results are consistent with the predictions from simulations (Ceverino et al. 2015, Tomassetti et al. 2016) that the transition from prolate to oblate is caused by a compaction event at a characteristic mass range, making the galaxy center baryon dominated. We give probabilities of a galaxy's being prolate, oblate or spheroidal as a function of its M, redshift, projected b/a and a, which can facilitate target selections of galaxies with specific shapes at hight redshifts. We also give predicted optical depths of galaxies, which are qualitatively consistent with the expected correlation that AV should be higher for edge-on disk galaxies in each log a slice at low redshift and high mass bins., Comment: 24 pages, 25 figures, submitted to MNRAS

We study the effects of the local environmental density and the cosmic web environment (filaments, walls, and voids) on key properties of dark matter halos using the Bolshoi-Planck LCDM cosmological simulation. The z = 0 simulation is analysed into filaments, walls, and voids using the SpineWeb method and also the VIDE package of tools, both of which use the watershed transform. The key halo properties that we study are the specific mass accretion rate, spin parameter, concentration, prolateness, scale factor of the last major merger, and scale factor when the halo had half of its z = 0 mass. For all these properties, we find that there is no discernible difference between the halo properties in filaments, walls, or voids when compared at the same environmental density. As a result, we conclude that environmental density is the core attribute that affects these properties. This conclusion is in line with recent findings that properties of galaxies in redshift surveys are independent of their cosmic web environment at the same environmental density at z ~ 0. We also find that the local web environment of the Milky Way and the Andromeda galaxies near the centre of a cosmic wall does not appear to have any effect on the properties of these galaxies' dark matter halos except for their orientation, although we find that it is rather rare to have such massive halos near the centre of a relatively small cosmic wall., Comment: 23 pages

The similarity between the distributions of spins for galaxies ($\lambda_{\rm g}$) and for dark-matter haloes ($\lambda_{\rm h}$), indicated both by simulations and observations, is naively interpreted as a one-to-one correlation between the spins of a galaxy and its host halo. This is used to predict galaxy sizes in semi-analytic models via $R_{\rm e}\simeq\lambda_{\rm h} R_{\rm v}$, with $R_{\rm e}$ the half-mass radius of the galaxy and $R_{\rm v}$ the halo radius. Utilizing two different suites of zoom-in cosmological simulations, we find that $\lambda_{\rm g}$ and $\lambda_{\rm h}$ are in fact only barely correlated, especially at $z\geq 1$. A general smearing of this correlation is expected based on the different spin histories, where the more recently accreted baryons through streams gain and then lose significant angular momentum compared to the gradually accumulated dark matter. Expecting the spins of baryons and dark matter to be correlated at accretion into $R_{\rm v}$, the null correlation at the end reflects an anti-correlation between $\lambda_{\rm g}/\lambda_{\rm h}$ and $\lambda_{\rm h}$, which can partly arise from mergers and a compact star-forming phase that many galaxies undergo. On the other hand, the halo and galaxy spin vectors tend to be aligned, with a median $\cos\theta=0.6$-0.7 between galaxy and halo, consistent with instreaming within a preferred plane. The galaxy spin is better correlated with the spin of the inner halo, but this largely reflects the effect of the baryons on the halo. Following the null spin correlation, $\lambda_{\rm h}$ is not a useful proxy for $R_{\rm e}$. While our simulations reproduce a general relation of the sort $R_{\rm e}=AR_{\rm vir}$, in agreement with observational estimates, the relation becomes tighter with $A=0.02(c/10)^{-0.7}$, where $c$ is the halo concentration, which in turn introduces a dependence on mass and redshift., Comment: 20 pages, 16 figures, submitted to MNRAS

Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0