A key uncertainty in galaxy evolution is the physics regulating star formation, ranging from small-scale processes related to the life-cycle of molecular clouds within galaxies to large-scale processes such as gas accretion onto galaxies. We study the imprint of such processes on the time-variability of star formation with an analytical approach tracking the gas mass of galaxies ("regulator model"). Specifically, we quantify the strength of the fluctuation in the star-formation rate (SFR) on different timescales, i.e. the power spectral density (PSD) of the star-formation history, and connect it to gas inflow and the life-cycle of molecular clouds. We show that in the general case the PSD of the SFR has three breaks, corresponding to the correlation time of the inflow rate, the equilibrium timescale of the gas reservoir of the galaxy, and the average lifetime of individual molecular clouds. On long and intermediate timescales (relative to the dynamical timescale of the galaxy), the PSD is typically set by the variability of the inflow rate and the interplay between outflows and gas depletion. On short timescales, the PSD shows an additional component related to the life-cycle of molecular clouds, which can be described by a damped random walk with a power-law slope of $\beta\approx2$ at high frequencies with a break near the average cloud lifetime. We discuss star-formation "burstiness" in a wide range of galaxy regimes, study the evolution of galaxies about the main sequence ridgeline, and explore the applicability of our method for understanding the star-formation process on cloud-scale from galaxy-integrated measurements.

In the second work of this series, we analyse the connection between the availability of gas and the position of a region with respect to the spatially resolved main sequence (MS) relation. Following the procedure presented in Paper I we obtain 500pc scales estimates of stellar mass and star formation rate surface densities ($\Sigma_{\star}$ and $\Sigma_{\rm{SFR}}$). Our sample consists of five face-on, grand design spiral galaxies located on the MS. Thanks to HI 21cm and $^{12}$CO(2-1) maps, we connect the gas surface densities and gas fractions to the observed star formation properties of each region. We find that the spatially resolved MS ($\sigma=0.23$ dex) is the combination of two relations: the Kennicutt-Schmidt law ($\sigma=0.19$ dex) and the molecular gas MS (MGMS, $\sigma=0.22$ dex); $\Sigma_{\star}$, $\Sigma_{\rm{SFR}}$ and the surface density of the molecular gas, $\Sigma_{\rm{H_2}}$, define a 3D relation as proposed by \citet{2019ApJ...884L..33L}. We find that $\Sigma_{\rm{H_2}}$ steadily increases along the MS relation, varies little towards higher $\Sigma_{\rm{SFR}}$ at fixed stellar surface densities (not enough to sustain the change in SFR), and it is almost constant perpendicular to the relation. The surface density of neutral gas ($\Sigma_{\rm{HI}}$) is constant along the MS, and increases in its upper envelop. $\Sigma_{\rm{SFR}}$ can be expressed as a function of $\Sigma_{\star}$ and $\Sigma_{\rm{HI}}$, following the Equation: $\log\Sigma_{\rm{SFR}}$ = 0.97$\log\Sigma_{\star}$ + 1.99$\log\Sigma_{\rm{HI}}$ - 11.11. Finally, we show that f$_{\rm{gas}}$ increases significantly towards the starburst region in the $\log\Sigma_{\star}$ - $\log\Sigma_{\rm{SFR}}$ plane, accompanied by a slight increase in SFE.

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.

We present predictions for high redshift ($z=2-10$) galaxy populations based on the IllustrisTNG simulation suite and a full Monte Carlo dust radiative transfer post-processing. Specifically, we discuss the ${\rm H}_{\alpha}$ and ${\rm H}_{\beta}$ + $[\rm O \,III]$ luminosity functions up to $z=8$. The predicted ${\rm H}_{\beta}$ + $[\rm O \,III]$ luminosity functions are consistent with present observations at $z\lesssim 3$ with $\lesssim 0.1\,{\rm dex}$ differences in luminosities. However, the predicted ${\rm H}_{\alpha}$ luminosity function is $\sim 0.3\,{\rm dex}$ dimmer than the observed one at $z\simeq 2$. Furthermore, we explore continuum spectral indices, the Balmer break at $4000$\AA (D4000) and the UV continuum slope $\beta$. The median D4000 versus sSFR relation predicted at $z=2$ is in agreement with the local calibration despite a different distribution pattern of galaxies in this plane. In addition, we reproduce the observed $A_{\rm UV}$ versus $\beta$ relation and explore its dependence on galaxy stellar mass, providing an explanation for the observed complexity of this relation. We also find a deficiency in heavily attenuated, UV red galaxies in the simulations. Finally, we provide predictions for the dust attenuation curves of galaxies at $z=2-6$ and investigate their dependence on galaxy colors and stellar masses. The attenuation curves are steeper in galaxies at higher redshifts, with bluer colors, or with lower stellar masses. We attribute these predicted trends to dust geometry. Overall, our results are consistent with present observations of high redshift galaxies. Future JWST observations will further test these predictions.

The protagonists of cosmic reionization remain elusive. Faint star-forming galaxies are leading candidates because they are numerous and may have significant ionizing photon escape fractions ($f_{esc}$). Here we update this picture via an empirical model that successfully predicts latest observations (e.g., the drop in star-formation density at z>8). We generate an ionizing spectrum for each galaxy in our model and constrain $f_{esc}$ using latest measurements of the reionization timeline (e.g., Ly$\alpha$ damping of quasars and galaxies at z>7). Assuming a constant $f_{esc}$, we find $M_{UV}$<-13.5 galaxies need $f_{esc}=0.21^{+0.06}_{-0.04}$ to complete reionization. The inferred IGM neutral fraction is [0.9, 0.5, 0.1] at z=[8.2, 6.8, 6.2]$\pm$0.2, i.e., the bulk of reionization transpires in 300 Myrs. Inspired by the emergent sample of Lyman Continuum (LyC) leakers that overwhelmingly displays higher-than-average star-formation surface density ($\Sigma$), we propose a model relating $f_{esc}$ to $\Sigma$ and find $f_{esc}\propto\Sigma^{0.4\pm0.1}$. Since $\Sigma$ falls by ~2.5 dex between z=8 and z=0, our model explains the humble upper limits on $f_{esc}$ at lower redshifts and its required evolution to ~0.2 at z>6. Within this model, strikingly, <5% of galaxies with $M_{UV}$<-18 (the `oligarchs') account for >80% of the reionization budget. In fact, faint sources ($M_{UV}$>-16) must be relegated to a limited role to ensure high neutral fractions at z=7-8. Shallow faint-end slopes of the UV luminosity function ($\alpha$>-2) and/or $f_{esc}$ distributions skewed toward bright galaxies produce the required late and rapid reionization. We predict LyC leakers like COLA1 (z=6.6, $f_{esc}$~30%, $M_{UV}$=-21.5) become increasingly common towards z~6 and that the drivers of reionization do not lie hidden across the faint-end of the luminosity function, but are already known to us. (abridged)

We study the photometric properties and sizes of the reionization-epoch galaxies in high-resolution Simba cosmological hydrodynamical simulations with box sizes of $[25,50]~h^{-1}$Mpc. Assuming various attenuation laws, we compute photometry by extincting each star particle's spectrum using the line-of-sight gas metal column density. The predicted ultraviolet luminosity function (UVLF) generally agrees with observations at $z=6$, owing to a partial cancellation between the high metallicities of the simulated galaxies and lower dust-to-metal ratios. The simulated $z=8$ UVLF is low compared to observations, likely owing to excessive dust extinction. Simba predicts UV continuum slopes ($\beta$) in agreement with the $z=6$ observations, with the best agreement obtained using a Calzetti extinction law. Interestingly, the gas-phase mass-metallicity relation in Simba is higher at $z\sim 6$ than at $z\sim 2$, suggesting that rapid early enrichment (and dust growth) might be necessary to match the observed $\beta$. We find that $\beta$ is more sensitive to the dust extinction law than the UVLF. By generating mock James Webb Space Telescope (JWST) images and analysing in a manner similar to observations, we show that Simba's galaxy size-luminosity relation well reproduces the current $z=6$ Hubble observations. Unlike observations at lower redshifts, Simba predicts similar rest-UV and rest-optical sizes of $z=6$ galaxies, owing to weak age gradients and dust extinction in star-forming regions counteract each other to weaken the color gradients within galaxies. These predictions will be testable with JWST.

The NIRCam instrument on the upcoming James Webb Space Telescope (JWST) will offer an unprecedented view of the most distant galaxies. In preparation for future deep NIRCam extragalactic surveys, it is crucial to understand the color selection of high-redshift galaxies using the Lyman dropout technique. To that end, we have used the JAdes extraGalactic Ultradeep Artificial Realizations (JAGUAR) mock catalog to simulate a series of extragalactic surveys with realistic noise estimates. This enables us to explore different color selections and their impact on the number density of recovered high-redshift galaxies and lower-redshift interlopers. We explore how survey depth, detection signal-to-noise ratio, color selection method, detection filter choice, and the presence of the Ly$\alpha$ emission line affects the resulting dropout selected samples. We find that redder selection colors reduce the number of recovered high-redshift galaxies, but the overall accuracy of the final sample is higher. In addition, we find that methods that utilize two or three color cuts have higher accuracy because of their ability to select against low-redshift quiescent and faint dusty interloper galaxies. We also explore the near-IR colors of brown dwarfs and demonstrate that, while they are predicted to have low on-sky densities, they are most likely to be recovered in F090W dropout selection, but there are color cuts which help to mitigate this contamination. Overall, our results provide NIRCam selection methods to aid in the creation of large, pure samples of ultra high-redshift galaxies from photometry alone.

We have obtained data for 41 star forming galaxies in the MUSE Atlas of Disks (MAD) survey with VLT/MUSE. These data allow us, at high resolution of a few 100 pc, to extract ionized gas kinematics ($V, \sigma$) of the centers of nearby star forming galaxies spanning 3 dex in stellar mass. This paper outlines the methodology for measuring the ionized gas kinematics, which we will use in subsequent papers of this survey. We also show how the maps can be used to study the kinematics of diffuse ionized gas for galaxies of various inclinations and masses. Using two different methods to identify the diffuse ionized gas, we measure rotation velocities of this gas for a subsample of 6 galaxies. We find that the diffuse ionized gas rotates on average slower than the star forming gas with lags of 0-10 km/s while also having higher velocity dispersion. The magnitude of these lags is on average 5 km/s lower than observed velocity lags between ionized and molecular gas. Using Jeans models to interpret the lags in rotation velocity and the increase in velocity dispersion we show that most of the diffuse ionized gas kinematics are consistent with its emission originating from a somewhat thicker layer than the star forming gas, with a scale height that is lower than that of the stellar disk.

The redshift range z=4-6 marks a transition phase between primordial and mature galaxy formation in which galaxies considerably increase their stellar mass, metallicity, and dust content. The study of galaxies in this redshift range is therefore important to understand early galaxy formation and the fate of galaxies at later times. Here, we investigate the burstiness of the recent star-formation history (SFH) of 221 $z\sim4.5$ main-sequence galaxies at log(M) > 9.7 by comparing their ultra-violet (UV) continuum, H$\alpha$ luminosity, and H$\alpha$ equivalent-width (EW). The H$\alpha$ properties are derived from the Spitzer [3.6$\mu$m]-[4.5$\mu$m] broad-band color, thereby properly taking into account model and photometric uncertainties. We find a significant scatter between H$\alpha$ and UV-derived luminosities and star-formation rates (SFRs). About half of the galaxies show a significant excess in H$\alpha$ compared to expectations from a constant smooth SFH. We also find a tentative anti-correlation between H$\alpha$ EW and stellar mass, ranging from 1000$\r{A}$ at log(M) < 10 to below 100$\r{A}$ at log(M) > 11. Consulting models suggests that most $z\sim4.5$ galaxies had a burst of star-formation within the last 50 Myrs, increasing their SFRs by a factor of > 5. The most massive galaxies on the other hand might decrease their SFRs, and may be transitioning to a quiescent stage by z=4. We identify differential dust attenuation (f) between stars and nebular regions as the main contributor to the uncertainty. With local galaxies selected by increasing H$\alpha$ EW (reaching values similar to high-z galaxies), we predict that f approaches unity at $z>4$ consistent with the extrapolation of measurements out to z=2.

Modelling reionization often requires significant assumptions about the properties of ionizing sources. Here, we infer the total output of hydrogen-ionizing photons (the ionizing emissivity, $\dot{N}_\textrm{ion}$) at $z=4-14$ from current reionization constraints, being maximally agnostic to the properties of ionizing sources. We use a Bayesian analysis to fit for a non-parametric form of $\dot{N}_\textrm{ion}$, allowing us to flexibly explore the entire prior volume. We infer a declining $\dot{N}_\textrm{ion}$ with redshift at $z>6$, which can be used as a benchmark for reionization models. Model-independent reionization constraints from the CMB optical depth and Ly$\alpha$ and Ly$\beta$ forest dark pixel fraction produce $\dot{N}_\textrm{ion}$ evolution ($d\log_{10}\dot{N}_\textrm{ion}/dz|_{z=6\rightarrow8} = -0.31\pm0.35$ dex) consistent with the declining UV luminosity density of galaxies, assuming constant ionizing photon escape fraction and efficiency. Including measurements from Ly$\alpha$ damping of galaxies and quasars produces a more rapid decline: $d\log_{10}\dot{N}_\textrm{ion}/dz|_{z=6\rightarrow8} =-0.44\pm0.22$ dex, steeper than the declining galaxy luminosity density (if extrapolated beyond $M_\mathrm{UV} \lesssim -13$), and constrains the mid-point of reionization to $z = 6.93\pm0.14$.

The UVJ color-color diagram is a popular and efficient method to distinguish between quiescent and star-forming galaxies through their rest-frame $U-V$ vs. $V-J$ colors. Here we explore the information content of this color-color space using the Bayesian inference machine Prospector. We fit the same physical model to two datasets: (i) UVJ fluxes alone, and (ii) full UV-mid IR (MIR) broadband SEDs from the 3D-HST survey. Notably this model uses both nonparametric SFHs and a flexible dust attenuation curve, both of which have the potential to `break' the typical correlations observed in UVJ color-color space. Instead, these fits confirm observed trends between UVJ colors and observed galaxy properties, including specific star formation rate (sSFR), dust attenuation, stellar age, and stellar metallicity. They also demonstrate that UVJ colors do not, on their own, constrain stellar age or metallicity; the observed trends in the UVJ diagram are instead driven by galaxy scaling relationships and thus will evolve with cosmological time. We also show that UVJ colors 'saturate' below $\log(\mathrm{sSFR/yr}^{-1})\lesssim -10.5$, i.e. changing sSFR no longer produces substantial changes in UVJ colors. We show that far-UV and/or MIR fluxes continue to correlate with sSFR down to low sSFRs and can be used in color-color diagrams to efficiently target galaxies with much lower levels of ongoing star formation. We provide selection criteria in these new color-color spaces as a function of desired sample sSFR.

Using the IllustrisTNG simulations, we investigate the connection between galaxy morphology and star formation in central galaxies with stellar masses in the range $10^9-10^{11.5}~\mathrm{M}_{\odot}$. We quantify galaxy morphology by a kinematical decomposition of the stellar component into a spheroidal and a disc component (spheroid-to-total ratio, S/T) and by the concentration of the stellar mass density profile ($C_{82}$). S/T is correlated with stellar mass and star-formation activity, while $C_{82}$ correlates only with stellar mass. Overall, we find good agreement with observational estimates for both S/T and $C_{82}$. Low and high mass galaxies are dominated by random stellar motion, while only intermediate-mass galaxies ($M_{\star}\approx10^{10}-10^{10.5}~\mathrm{M}_{\odot}$) are dominated by ordered rotation. Whereas higher-mass galaxies are typical spheroids with high concentrations, lower-mass galaxies have low concentration, pointing to different formation channels. Although we find a correlation between S/T and star-formation activity, in the TNG model galaxies do not necessarily change their morphology when they transition through the green valley or when they cease their star formation, this depending on galaxy stellar mass and morphological estimator. Instead, the morphology (S/T and $C_{82}$) is generally set during the star-forming phase of galaxies. The apparent correlation between S/T and star formation arises because earlier-forming galaxies had, on average, a higher S/T at a given stellar mass. Furthermore, we show that mergers drive in-situ bulge formation in intermediate-mass galaxies and are responsible for the recent spheroidal mass assembly in the massive galaxies with $M_{\star}>10^{11}~\mathrm{M}_{\odot}$. In particular, these massive galaxies assemble about half of the spheroidal mass while star-forming and the other half through mergers while quiescent.

We have recently developed a post-processing framework to estimate the abundance of atomic and molecular hydrogen (HI and H2, respectively) in galaxies in large-volume cosmological simulations. Here we compare the HI and H2 content of IllustrisTNG galaxies to observations. We mostly restrict this comparison to $z \approx 0$ and consider six observational metrics: the overall abundance of HI and H2, their mass functions, gas fractions as a function of stellar mass, the correlation between H2 and star formation rate, the spatial distribution of gas, and the correlation between gas content and morphology. We find generally good agreement between simulations and observations, particularly for the gas fractions and the HI mass-size relation. The H2 mass correlates with star formation rate as expected, revealing an almost constant depletion time that evolves up to z = 2 as observed. However, we also discover a number of tensions with varying degrees of significance, including an overestimate of the total neutral gas abundance at z = 0 by about a factor of two and a possible excess of satellites with no or very little neutral gas. These conclusions are robust to the modelling of the HI/H2 transition. In terms of their neutral gas properties, the IllustrisTNG simulations represent an enormous improvement over the original Illustris run. All data used in this paper are publicly available as part of the IllustrisTNG data release.

The James Webb Space Telescop (JWST) promises to revolutionise our understanding of the early Universe, and contrasting its upcoming observations with predictions of the $\Lambda$CDM model requires detailed theoretical forecasts. Here, we exploit the large dynamic range of the IllustrisTNG simulation suite, TNG50, TNG100, and TNG300, to derive multi-band galaxy luminosity functions from $z=2$ to $z=10$. We put particular emphasis on the exploration of different dust attenuation models to determine galaxy luminosity functions for the rest-frame ultraviolet (UV), and apparent wide NIRCam bands. Our most detailed dust model is based on continuum Monte Carlo radiative transfer calculations employing observationally calibrated dust properties. This calibration results in constraints on the redshift evolution of the dust attenuation normalisation and dust-to-metal ratios yielding a stronger redshift evolution of the attenuation normalisation compared to most previous theoretical studies. Overall we find good agreement between the rest-frame UV luminosity functions and observational data for all redshifts, also beyond the regimes used for the dust-model calibrations. Furthermore, we also recover the observed high redshift ($z=4-6$) UV luminosity versus stellar mass relation, the H$\alpha$ versus star formation rate relation, and the H$\alpha$ luminosity function at $z=2$. The bright end ($M_{\rm UV}>-19.5$) cumulative galaxy number densities are consistent with observational data. For the F200W NIRCam band, we predict that JWST will detect $\sim 80$ ($\sim 200$) galaxies with a signal-to-noise ratio of $10$ ($\sim 5$) within the NIRCam field of view, $2.2\times2.2 \,{\rm arcmin}^{2}$, for a total exposure time of $10^5{\rm s}$ in the redshift range $z=8 \pm 0.5$. These numbers drop to $\sim 10$ ($\sim 40$) for an exposure time of $10^4{\rm s}$.

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

We investigate the relationship between star formation activity and outflow properties on kiloparsec scales in a sample of 28 star forming galaxies at $z\sim$ 2-2.6, using adaptive optics assisted integral field observations from SINFONI on the VLT. The narrow and broad components of the H$\alpha$ emission are used to simultaneously determine the local star formation rate surface density ($\Sigma_{\rm SFR}$), and the outflow velocity $v_{\rm out}$ and mass outflow rate $\dot{M}_{\rm out}$, respectively. We find clear evidence for faster outflows with larger mass loading factors at higher $\Sigma_{\rm SFR}$. The outflow velocities scale as $v_{\rm out}$ $\propto$ $\Sigma_{\rm SFR}^{0.34 \pm 0.10}$, which suggests that the outflows may be driven by a combination of mechanical energy released by supernova explosions and stellar winds, as well as radiation pressure acting on dust grains. The majority of the outflowing material does not have sufficient velocity to escape from the galaxy halos, but will likely be re-accreted and contribute to the chemical enrichment of the galaxies. In the highest $\Sigma_{\rm SFR}$ regions the outflow component contains an average of $\sim$45% of the H$\alpha$ flux, while in the lower $\Sigma_{\rm SFR}$ regions only $\sim$10% of the H$\alpha$ flux is associated with outflows. The mass loading factor, $\eta$ = $\dot{M}_{\rm out}$/SFR, is positively correlated with $\Sigma_{\rm SFR}$ but is relatively low even at the highest $\Sigma_{\rm SFR}$: $\eta \lesssim$ 0.5 $\times$ (380 cm$^{-3}$/n$_e$). This may be in tension with the $\eta$ $\gtrsim$ 1 required by cosmological simulations, unless a significant fraction of the outflowing mass is in other gas phases and has sufficient velocity to escape the galaxy halos.

We present a framework for modelling the star-formation histories of galaxies as a stochastic process. We define this stochastic process through a power spectrum density with a functional form of a broken power-law. Star-formation histories are correlated on short timescales, the strength of this correlation described by a power-law slope, $\alpha$, and they decorrelate to resemble white noise over a timescale that is proportional to the timescale of the break in the power spectrum density, $\tau_{\rm break}$. We use this framework to explore the properties of the stochastic process that, we assume, gives rise to the log-normal scatter about the relationship between star-formation rate and stellar mass, the so-called galaxy star-forming main sequence. Specifically, we show how the measurements of the normalisation and width ($\sigma_{\rm MS}$) of the main sequence, measured in several passbands that probe different timescales, give a constraint on the parameters of the underlying power spectrum density. We first derive these results analytically for a simplified case where we model observations by averaging over the recent star-formation history. We then run numerical simulations to find results for more realistic observational cases. As a proof of concept, we use observational estimates of the main sequence scatter at $z\sim0$ and $M_{\star}\approx10^{10}~M_{\odot}$ measured in H$\alpha$, UV+IR and the u-band, and show that combination of these point to $\tau_{\rm break}=178^{+104}_{-66}$ Myr, when assuming $\alpha=2$. This implies that star-formation histories of galaxies lose "memory" of their previous activity on a timescale of $\sim200$ Myr, highlighting the importance of baryonic effects that act over the dynamical timescales of galaxies.

We study the physical properties of the ionized gas in local disks using the sample of 38 nearby $\sim10^{8.5-11.2}$M$_\odot$ Star-Forming Main Sequence (SFMS) galaxies observed so far as part of the MUSE Atlas of Disks (MAD). Specifically, we use all strong emission lines in the MUSE wavelength range 4650-9300 \AA\ to investigate the resolved ionized gas properties on $\sim$100 pc scales. This spatial resolution enables us to disentangle HII regions from the Diffuse Ionized Gas (DIG) in the computation of gas metallicities and Star Formation Rates (SFRs) of star forming regions. The gas metallicities generally decrease with radius. The metallicity of the HII regions is on average $\sim$0.1 dex higher than that of the DIG, but the metallicity radial gradient in both components is similar. The mean metallicities within the inner galaxy cores correlate with the total stellar mass of the galaxies. On our <100 pc scales, we find two correlations previously reported at kpc scales: a spatially resolved Mass-Metallicity Relation (RMZR) and a spatially resolved SFMS (RSFMS). We find no secondary dependency of the RMZR with the SFR density. We find that both resolved relations have a local origin, as they do not depend on the total stellar mass. The observational results of this paper are consistent with the inside-out scenario for the growth of galactic disks.

Accreting supermassive black holes (SMBHs) can exhibit variable emission across the electromagnetic spectrum and over a broad range of timescales. The variability of active galactic nuclei (AGNs) in the ultraviolet and optical is usually at the few tens of per cent level over timescales of hours to weeks. Recently, rare, more dramatic changes to the emission from accreting SMBHs have been observed, including tidal disruption events, 'changing look' AGNs and other extreme variability objects. The physics behind the 're-ignition', enhancement and 'shut-down' of accretion onto SMBHs is not entirely understood. Here we present a rapid increase in ultraviolet-optical emission in the centre of a nearby galaxy, marking the onset of sudden increased accretion onto a SMBH. The optical spectrum of this flare, dubbed AT 2017bgt, exhibits a mix of emission features. Some are typical of luminous, unobscured AGNs, but others are likely driven by Bowen fluorescence - robustly linked here with high-velocity gas in the vicinity of the accreting SMBH. The spectral features and increased ultraviolet flux show little evolution over a period of at least 14 months. This disfavours the tidal disruption of a star as their origin, and instead suggests a longer-term event of intensified accretion. Together with two other recently reported events with similar properties, we define a new class of SMBH-related flares. This has important implications for the classification of different types of enhanced accretion onto SMBHs.

We have generated synthetic images of $\sim$27,000 galaxies from the IllustrisTNG and the original Illustris hydrodynamic cosmological simulations, designed to match Pan-STARRS observations of $\log_{10}(M_{\ast}/{\rm M}_{\odot}) \approx 9.8$-$11.3$ galaxies at $z \approx 0.05$. Most of our synthetic images were created with the SKIRT radiative transfer code, including the effects of dust attenuation and scattering, and performing the radiative transfer directly on the Voronoi mesh used by the simulations themselves. We have analysed both our synthetic and real Pan-STARRS images with the newly developed $\tt{statmorph}$ code, which calculates non-parametric morphological diagnostics -- including the Gini-$M_{20}$ and concentration-asymmetry-smoothness (CAS) statistics -- and performs two-dimensional S\'ersic fits. Overall, we find that the optical morphologies of IllustrisTNG galaxies are in good agreement with observations, and represent a substantial improvement compared to the original Illustris simulation. In particular, the locus of the Gini-$M_{20}$ diagram is consistent with that inferred from observations, while the median trends with stellar mass of all the morphological, size and shape parameters considered in this work lie within the $\sim$1$\sigma$ scatter of the observational trends. However, the IllustrisTNG model has some difficulty with more stringent tests, such as producing a strong morphology-colour relation. This results in a somewhat higher fraction of red discs and blue spheroids compared to observations. Similarly, the morphology-size relation is problematic: while observations show that discs tend to be larger than spheroids at a fixed stellar mass, such a trend is not present in IllustrisTNG.

We explore the connection between the UV luminosity functions (LFs) of high-$z$ galaxies and the distribution of stellar masses and star-formation histories (SFHs) in their host dark matter halos. We provide a baseline for a redshift-independent star-formation efficiency model to which observations and models can be compared. Our model assigns a star-formation rate (SFR) to each dark matter halo based on the growth rate of the halo and a redshift-independent star-formation efficiency. The dark matter halo accretion rate is obtained from a high-resolution $N$-body simulation in order to capture the stochasticity in accretion histories and to obtain spatial information for the distribution of galaxies. The halo mass dependence of the star-formation efficiency is calibrated at $z=4$ by requiring a match to the observed UV LF at this redshift. The model then correctly predicts the observed UV LF at $z=5-10$. We present predictions for the UV luminosity and stellar mass functions, \textit{JWST} number counts, and SFHs. In particular, we find a stellar-to-halo mass relation at $z=4-10$ that scales with halo mass at $M_{\rm h}<10^{11}~M_{\odot}$ as $M_{\star}\propto M_{\rm h}^2$, with a normalization that is higher than the relation inferred at $z=0$. The average SFRs increase as a function of time to $z=4$, although there is significant scatter around the average: about 6\% of the $z=4$ galaxies show no significant mass growth. Using these SFHs, we present redshift-dependent UV-to-SFR conversion factors, mass return fractions, and mass-to-light ratios for different intial mass functions and metallicities, finding that current estimates of the cosmic SFR density at $z\sim10$ may be overestimated by $\sim0.1-0.2~\mathrm{dex}$.

Star-forming galaxies have been found to follow a relatively tight relation between stellar mass ($M_{*}$) and star formation rate (SFR), dubbed the `star formation sequence'. A turnover in the sequence has been observed, where galaxies with $M_{*} < 10^{10} {\rm M}_{\odot}$ follow a steeper relation than their higher mass counterparts, suggesting that the low-mass slope is (nearly) linear. In this paper, we characterise the properties of the low-mass end of the star formation sequence between $7 \leq \log M_{*}[{\rm M}_{\odot}] \leq 10.5$ at redshift $0.11 < z < 0.91$. We use the deepest MUSE observations of the Hubble Ultra Deep Field and the Hubble Deep Field South to construct a sample of 179 star-forming galaxies with high signal-to-noise emission lines. Dust-corrected SFRs are determined from H$\beta$ $\lambda 4861$ and H$\alpha$ $\lambda 6563$. We model the star formation sequence with a Gaussian distribution around a hyperplane between $\log M_{*}$, $\log {\rm SFR}$, and $\log (1+z)$, to simultaneously constrain the slope, redshift evolution, and intrinsic scatter. We find a sub-linear slope for the low-mass regime where $\log {\rm SFR}[{\rm M}_{\odot}/{\rm yr}] = 0.83^{+0.07}_{-0.06} \log M_{*}[{\rm M}_{\odot}] + 1.74^{+0.66}_{-0.68} \log (1+z)$, increasing with redshift. We recover an intrinsic scatter in the relation of $\sigma_{\rm intr} = 0.44^{+0.05}_{-0.04}$ dex, larger than typically found at higher masses. As both hydrodynamical simulations and (semi-)analytical models typically favour a steeper slope in the low-mass regime, our results provide new constraints on the feedback processes which operate preferentially in low-mass halos.

We present the stellar mass profiles of 147 isolated quiescent galaxies in very low-density environments (i.e., void regions) in the local Universe ($0.01<z<0.06$) from the Sloan Digital Sky Survey. These galaxies have stellar masses between $ 9.8\lesssim \log(M_{\ast}/M_{\odot})\lesssim11.2$ and they represent $\sim15\%$ of the whole galaxy population in the void regions down to $M_{r} = -19$. We do not find any isolated quiescent galaxies with $\log(M_{\ast}/M_{\odot})\gtrsim11.2$. We compare the stellar mass profiles of these isolated quiescent galaxies with the profiles of stellar mass-matched samples of the quiescent galaxies in group and cluster environments. We find that, at fixed mass, quiescent galaxies in voids have similar central ($1$ kpc) mass densities ($\Sigma_1$) and central velocity dispersions ($\sigma_1$) compared to their counterparts in groups and clusters. We show that quiescent galaxies in voids have at most $10-25\%$ smaller half-mass (and half-light) sizes compared to quiescent galaxies in groups and clusters. We conclude that for the intermediate stellar mass range of $10^{10}-10^{11}M_{\odot}$ in the local Universe, environmental mechanisms have no significant additional effect on the mass profiles of the quiescent galaxies.

We derive two-dimensional dust attenuation maps at $\sim1~\mathrm{kpc}$ resolution from the UV continuum for ten galaxies on the $z\sim2$ Star-Forming Main Sequence (SFMS). Comparison with IR data shows that 9 out of 10 galaxies do not require further obscuration in addition to the UV-based correction, though our sample does not include the most heavily obscured, massive galaxies. The individual rest-frame $V$-band dust attenuation (A$_{\rm V}$) radial profiles scatter around an average profile that gently decreases from $\sim1.8$ mag in the center down to $\sim0.6$ mag at $\sim3-4$ half-mass radii. We use these maps to correct UV- and H$\alpha$-based star-formation rates (SFRs), which agree with each other. At masses $<10^{11}~M_{\rm sun}$, the dust-corrected specific SFR (sSFR) profiles are on average radially constant at a mass-doubling timescale of $\sim300~\mathrm{Myr}$, pointing at a synchronous growth of bulge and disk components. At masses $>10^{11}~M_{\rm sun}$, the sSFR profiles are typically centrally-suppressed by a factor of $\sim10$ relative to the galaxy outskirts. With total central obscuration disfavored, this indicates that at least a fraction of massive $z\sim2$ SFMS galaxies have started their inside-out star-formation quenching that will move them to the quenched sequence. In combination with other observations, galaxies above and below the ridge of the SFMS relation have respectively centrally-enhanced and centrally-suppressed sSFRs relative to their outskirts, supporting a picture where bulges are built due to gas `compaction' that leads to a high central SFR as galaxies move towards the upper envelope of SFMS.

We present an original phenomenological model to describe the evolution of galaxy number counts, morphologies, and spectral energy distributions across a wide range of redshifts (0.2<z<15) and stellar masses [Log10 M/Msun >6]. Our model follows observed mass and luminosity functions of both star-forming and quiescent galaxies, and reproduces the redshift evolution of colors, sizes, star-formation and chemical properties of the observed galaxy population. Unlike other existing approaches, our model includes a self-consistent treatment of stellar and photoionized gas emission and dust attenuation based on the BEAGLE tool. The mock galaxy catalogs generated with our new model can be used to simulate and optimize extragalactic surveys with future facilities such as the James Webb Space Telescope (JWST), and to enable critical assessments of analysis procedures, interpretation tools, and measurement systematics for both photometric and spectroscopic data. As a first application of this work, we make predictions for the upcoming JWST Advanced Deep Extragalactic Survey (JADES), a joint program of the JWST/NIRCam and NIRSpec Guaranteed Time Observations teams. We show that JADES will detect, with NIRCam imaging, thousands of galaxies at z>6, and tens at z>10 at m_AB<30 (5-sigma) within the 236 arcmin^2 of the survey. The JADES data will enable accurate constraints on the evolution of the UV luminosity function at z>8, and resolve the current debate about the rate of evolution of galaxies at z>8. Ready to use mock catalogs and software to generate new realizations are publicly available as the JAdes extraGalactic Ultradeep Artificial Realizations (JAGUAR) package.

We present the "SINS/zC-SINF AO survey" of 35 star-forming galaxies, the largest sample with deep adaptive optics-assisted (AO) near-infrared integral field spectroscopy at z~2. The observations, taken with SINFONI at the Very Large Telescope, resolve the Ha and [NII] line emission and kinematics on scales of ~1.5 kpc. In stellar mass, star formation rate, rest-optical colors and size, the AO sample is representative of its parent seeing-limited sample and probes the massive (M* ~ 2x10^9 - 3x10^11 Msun), actively star-forming (SFR ~ 10-600 Msun/yr) part of the z~2 galaxy population over a wide range in colors ((U-V)_rest ~ 0.15-1.5 mag) and half-light radii (R_e,H ~ 1-8.5 kpc). The sample overlaps largely with the "main sequence" of star-forming galaxies in the same redshift range to a similar K_AB = 23 magnitude limit; it has ~0.3 dex higher median specific SFR, ~0.1 mag bluer median (U-V)_rest color, and ~10% larger median rest-optical size. We describe the observations, data reduction, and extraction of basic flux and kinematic properties. With typically 3-4 times higher resolution and 4-5 times longer integrations (up to 23hr) than the seeing-limited datasets of the same objects, the AO data reveal much more detail in morphology and kinematics. The now complete AO observations confirm the majority of kinematically-classified disks and the typically elevated disk velocity dispersions previously reported based on subsets of the data. We derive typically flat or slightly negative radial [NII]/Ha gradients, with no significant trend with global galaxy properties, kinematic nature, or the presence of an AGN. Azimuthal variations in [NII]/Ha are seen in several sources and are associated with ionized gas outflows, and possible more metal-poor star-forming clumps or small companions. [Abridged]

We investigate the origin of the evolution of the population-averaged central stellar mass density ($\Sigma_1$) of quiescent galaxies (QGs) by probing the relation between stellar age and $\Sigma_1$ at $z\sim0$. We use the Zurich ENvironmental Study (ZENS), which is a survey of galaxy groups with a large fraction of satellite galaxies. QGs shape a narrow locus in the $\Sigma_1-M_{\star}$ plane, which we refer to as $\Sigma_1$ ridgeline. Colors of ($B-I$) and ($I-J$) are used to divide QGs into three age categories: young ($<2~\mathrm{Gyr}$), intermediate ($2-4~\mathrm{Gyr}$), and old ($>4~\mathrm{Gyr}$). At fixed stellar mass, old QGs on the $\Sigma_1$ ridgeline have higher $\Sigma_1$ than young QGs. This shows that galaxies landing on the $\Sigma_1$ ridgeline at later epochs arrive with lower $\Sigma_1$, which drives the zeropoint of the ridgeline down with time. We compare the present-day zeropoint of the oldest population at $z=0$ with the zeropoint of the quiescent population 4 Gyr back in time, at $z=0.37$. These zeropoints are identical, showing that the intrinsic evolution of individual galaxies after they arrive on the $\Sigma_1$ ridgeline must be negligible, or must evolve parallel to the ridgeline during this interval. The observed evolution of the global zeropoint of 0.07 dex over the last 4 Gyr is thus largely due to the continuous addition of newly quenched galaxies with lower $\Sigma_1$ at later times ("progenitor bias"). While these results refer to the satellite-rich ZENS sample as a whole, our work suggests a similar age-$\Sigma_1$ trend for central galaxies.

The rest-frame UV-optical (i.e., $NUV-B$) color is sensitive to both low-level recent star formation (specific star formation rate - sSFR) and dust. In this Letter, we extend our previous work on the origins of $NUV-B$ color gradients in star-forming galaxies (SFGs) at $z\sim1$ to those at $z\sim2$. We use a sample of 1335 large (semi-major axis radius $R_{\rm SMA}>0.''18$) SFGs with extended UV emission out to $2R_{\rm SMA}$ in the mass range $M_{\ast} = 10^{9}-10^{11}M_{\odot}$ at $1.5<z<2.8$ in the CANDELS/GOODS-S and UDS fields. We show that these SFGs generally have negative $NUV-B$ color gradients (redder centres), and their color gradients strongly increase with galaxy mass. We also show that the global rest-frame $FUV-NUV$ color is approximately linear with $A_{\rm V}$, which is derived by modeling the observed integrated FUV to NIR spectral energy distributions of the galaxies. Applying this integrated calibration to our spatially-resolved data, we find a negative dust gradient (more dust extinguished in the centers), which steadily becomes steeper with galaxy mass. We further find that the $NUV-B$ color gradients become nearly zero after correcting for dust gradients regardless of galaxy mass. This indicates that the sSFR gradients are negligible and dust reddening is likely the principal cause of negative UV-optical color gradients in these SFGs. Our findings support that the buildup of the stellar mass in SFGs at the Cosmic Noon is self-similar inside $2R_{\rm SMA}$.

Cosmological simulations suggest that most of the matter in the Universe is distributed along filaments connecting galaxies. Illuminated by the cosmic UV background (UVB), these structures are expected to glow in fluorescent Lyman alpha emission with a Surface Brightness (SB) that is well below current observational limits for individual detections. Here, we perform a stacking analysis of the deepest MUSE/VLT data using three-dimensional regions (subcubes) with orientations determined by the position of neighbouring Lyman alpha galaxies (LAEs) at 3<z<4. Our method should increase the probability of detecting filamentary Lyman alpha emission, provided that these structures are Lyman Limit Systems (LLSs). By stacking 390 oriented subcubes we reach a 2 sigma sensitivity level of SB ~ 0.44e-20 erg/s/cm^2/arcsec^2 in an aperture of 1 arcsec^2 x 6.25 Angstrom, which is three times below the expected fluorescent Lyman alpha signal from the Haardt-Madau 2012 (HM12) UVB at z~3.5. No detectable emission is found on intergalactic scales, implying that at least two thirds of our subcubes do not contain oriented LLSs for a HM12 UVB. On the other hand, significant emission is detected in the circum-galactic medium (CGM) of galaxies in the direction of the neighbours. The signal is stronger for galaxies with a larger number of neighbours and appears to be independent of any other galaxy properties such as luminosity, redshift and neighbour distance. We estimate that preferentially oriented satellite galaxies cannot contribute significantly to this signal, suggesting instead that gas densities in the CGM are typically larger in the direction of neighbouring galaxies on cosmological scales.

We use $>$9400 $\log(m/M_{\odot})>10$ quiescent and star-forming galaxies at $z\lesssim2$ in COSMOS/UltraVISTA to study the average size evolution of these systems, with focus on the rare, ultra-massive population at $\log(m/M_{\odot})>11.4$. The large 2-square degree survey area delivers a sample of $\sim400$ such ultra-massive systems. Accurate sizes are derived using a calibration based on high-resolution images from the Hubble Space Telescope. We find that, at these very high masses, the size evolution of star-forming and quiescent galaxies is almost indistinguishable in terms of normalization and power-law slope. We use this result to investigate possible pathways of quenching massive $m>M^*$ galaxies at $z<2$. We consistently model the size evolution of quiescent galaxies from the star-forming population by assuming different simple models for the suppression of star-formation. These models include an instantaneous and delayed quenching without altering the structure of galaxies and a central starburst followed by compaction. We find that instantaneous quenching reproduces well the observed mass-size relation of massive galaxies at $z>1$. Our starburst$+$compaction model followed by individual growth of the galaxies by minor mergers is preferred over other models without structural change for $\log(m/M_{\odot})>11.0$ galaxies at $z>0.5$. None of our models is able to meet the observations at $m>M^*$ and $z<1$ with out significant contribution of post-quenching growth of individual galaxies via mergers. We conclude that quenching is a fast process in galaxies with $ m \ge 10^{11} M_\odot$, and that major mergers likely play a major role in the final steps of their evolution.

We exploit the deep resolved Halpha kinematic data from the KMOS^3D and SINS/zC-SINF surveys to examine the largely unexplored outer disk kinematics of star-forming galaxies (SFGs) out to the peak of cosmic star formation. Our sample contains 101 SFGs representative of the more massive (9.3 < log(M*/Msun) < 11.5) main sequence population at 0.6<z<2.6. Through a novel stacking approach we are able to constrain a representative rotation curve extending out to ~4 effective radii. This average rotation curve exhibits a significant drop in rotation velocity beyond the turnover, with a slope of Delta(V)/Delta(R) = $-0.26^{+0.10}_{-0.09}$ in units of normalized coordinates V/V_max and R/R_turn. This result confirms that the fall-off seen previously in some individual galaxies is a common feature of our sample of high-z disks. We show that this outer fall-off strikingly deviates from the flat or mildly rising rotation curves of local spiral galaxies of similar masses. We furthermore compare our data with models including baryons and dark matter demonstrating that the falling stacked rotation curve can be explained by a high mass fraction of baryons relative to the total dark matter halo (m_d>~0.05) in combination with a sizeable level of pressure support in the outer disk. These findings are in agreement with recent studies demonstrating that star-forming disks at high redshift are strongly baryon dominated within the disk scale, and furthermore suggest that pressure gradients caused by large turbulent gas motions are present even in their outer disks. We demonstrate that these results are largely independent of our model assumptions such as the presence of a central stellar bulge, the effect of adiabatic contraction at fixed m_d, and variations in the concentration parameter.

In cold dark matter cosmology, the baryonic components of galaxies are thought to be mixed with and embedded in non-baryonic and non-relativistic dark matter, which dominates the total mass of the galaxy and its dark matter halo. In the local Universe, the mass of dark matter within a galactic disk increases with disk radius, becoming appreciable and then dominant in the outer, baryonic regions of the disks of star-forming galaxies. This results in rotation velocities of the visible matter within the disk that are constant or increasing with disk radius. Comparison between the dynamical mass and the sum of stellar and cold gas mass at the peak epoch of galaxy formation, inferred from ancillary data, suggest high baryon factions in the inner, star-forming regions of the disks. Although this implied baryon fraction may be larger than in the local Universe, the systematic uncertainties (stellar initial mass function, calibration of gas masses) render such comparisons inconclusive in terms of the mass of dark matter. Here we report rotation curves for the outer disks of six massive star-forming galaxies, and find that the rotation velocities are not constant, but decrease with radius. We propose that this trend arises because of two main factors: first, a large fraction of the massive, high-redshift galaxy population was strongly baryon dominated, with dark matter playing a smaller part than in the local Universe; and second, the large velocity dispersion in high-redshift disks introduces a substantial pressure term that leads to a decrease in rotation velocity with increasing radius. The effect of both factors appears to increase with redshift. Qualitatively, the observations suggest that baryons in the early Universe efficiently condensed at the centres of dark matter halos when gas fractions were high, and dark matter was less concentrated. [Abridged]

We study the history from $z\sim2$ to $z\sim0$ of the stellar mass assembly of quiescent and star-forming galaxies in a spatially resolved fashion. For this purpose we use multi-wavelength imaging data from the Hubble Space Telescope (HST) over the GOODS fields and the Sloan Digital Sky Survey (SDSS) for the local population. We present the radial stellar mass surface density profiles of galaxies with $M_{\ast}>10^{10} M_{\odot}$, corrected for mass-to-light ratio ($M_{\ast}/L$) variations, and derive the half-mass radius ($R_{m}$), central stellar mass surface density within 1 kpc ($\Sigma_{1}$) and surface density at $R_{m}$ ($\Sigma_{m}$) for star-forming and quiescent galaxies and study their evolution with redshift. At fixed stellar mass, the half-mass sizes of quiescent galaxies increase from $z\sim2$ to $z\sim0$ by a factor of $\sim3-5$, whereas the half-mass sizes of star-forming galaxies increase only slightly, by a factor of $\sim2$. The central densities $\Sigma_{1}$ of quiescent galaxies decline slightly (by a factor of $\lesssim1.7$) from $z\sim2$ to $z\sim0$, while for star-forming galaxies $\Sigma_{1}$ increases with time, at fixed mass. We show that the central density $\Sigma_{1}$ has a tighter correlation with specific star-formation rate (sSFR) than $\Sigma_{m}$ and for all masses and redshifts galaxies with higher central density are more prone to be quenched. Reaching a high central density ($\Sigma_{1} \gtrsim 10^{10} M_{\odot} \mathrm{kpc}^2$) seems to be a prerequisite for the cessation of star formation, though a causal link between high $\Sigma_{1}$ and quenching is difficult to prove and their correlation can have a different origin.

We analyze the role of first (leading) author gender on the number of citations that a paper receives, on the publishing frequency and on the self-citing tendency. We consider a complete sample of over 200,000 publications from 1950 to 2015 from five major astronomy journals. We determine the gender of the first author for over 70% of all publications. The fraction of papers which have a female first author has increased from less than 5% in the 1960s to about 25% today. We find that the increase of the fraction of papers authored by females is slowest in the most prestigious journals such as Science and Nature. Furthermore, female authors write 19$\pm$7% fewer papers in seven years following their first paper than their male colleagues. At all times papers with male first authors receive more citations than papers with female first authors. This difference has been decreasing with time and amounts to $\sim$6% measured over the last 30 years. To account for the fact that the properties of female and male first author papers differ intrinsically, we use a random forest algorithm to control for the non-gender specific properties of these papers which include seniority of the first author, number of references, total number of authors, year of publication, publication journal, field of study and region of the first author's institution. We show that papers authored by females receive 10.4$\pm$0.9% fewer citations than what would be expected if the papers with the same non-gender specific properties were written by the male authors. Finally, we also find that female authors in our sample tend to self-cite more, but that this effect disappears when controlled for non-gender specific variables.

Using the Sloan Digital Sky Survey, we adopt the sSFR-$\Sigma_{1kpc}$ diagram as a diagnostic tool to understand quenching in different environments. sSFR is the specific star formation rate, and $\Sigma_{1kpc}$ is the stellar surface density in the inner kpc. Although both the host halo mass and group-centric distance affect the satellite population, we find that these can be characterised by a single number, the quenched fraction, such that key features of the sSFR-$\Sigma_{1kpc}$ diagram vary smoothly with this proxy for the "environment". Particularly, the sSFR of star-forming galaxies decreases smoothly with this quenched fraction, the sSFR of satellites being 0.1 dex lower than in the field. Furthermore, $\Sigma_{1kpc}$ of the transition galaxies (i.e., the "green valley" or GV) decreases smoothly with the environment, by as much as 0.2 dex for $M_* = 10^{9.75-10} M_{\odot}$ from the field, and decreasing for satellites in larger halos and at smaller radial distances within same-mass halos. We interpret this shift as indicating the relative importance of today's field quenching track vs. the cluster quenching track. These environmental effects in the sSFR-$\Sigma_{1kpc}$ diagram are most significant in our lowest mass range ($9.75 < \log M_{*}/M_{\odot} < 10$). One feature that is shared between all environments is that at a given $M_{*}$ quenched galaxies have about 0.2-0.3 dex higher $\Sigma_{1kpc}$ than the star-forming population. These results indicate that either $\Sigma_{1kpc}$ increases (subsequent to satellite quenching), or $\Sigma_{1kpc}$ for individual galaxies remains unchanged, but the original $M_*$ or the time of quenching is significantly different from those now in the GV.

We investigate the origin of the evolution of the population-averaged size of quenched galaxies (QGs) through a spectroscopic analysis of their stellar ages. The two most favoured scenarios for this evolution are either the size growth of individual galaxies through a sequence of dry minor merger events, or the addition of larger, newly quenched galaxies to the pre-existing population (i.e., a progenitor bias effect). We use the 20k zCOSMOS-bright spectroscopic survey to select bona fide quiescent galaxies at 0.2<z<0.8. We stack their spectra in bins of redshift, stellar mass and size to compute stellar population parameters in these bins through fits to the rest-frame optical spectra and through Lick spectral indices. We confirm a change of behaviour in the size-age relation below and above the ~10^11 MSun stellar mass scale: In our 10.5 < log M*/MSun < 11 mass bin, over the entire redshift window, the stellar populations of the largest galaxies are systematically younger than those of the smaller counterparts, pointing at progenitor bias as the main driver of the observed average size evolution at sub-10^11 MSun masses. In contrast, at higher masses, there is no clear trend in age as a function of galaxy size, supporting a substantial role of dry mergers in increasing the sizes of these most massive QGs with cosmic time. Within the errors, the [alpha/Fe] abundance ratios of QGs are (i) above-solar over the entire redshift range of our analysis, hinting at universally short timescales for the buildup of the stellar populations of QGs, and (ii) similar at all masses and sizes, suggesting similar (short) timescales for the whole QG population and strengthening the role of mergers in the buildup of the most massive QGs in the Universe.

Using cosmological simulations, we address the interplay between structure and star formation in high-redshift galaxies via the evolution of surface density profiles. Our sample consists of 26 galaxies evolving in the redshift range $z=7-1$, spanning the stellar mass range $(0.2-6.4)\times 10^{10}M_\odot$ at $z=2$. We recover the main trends by stacking the profiles in accordance to their evolution phases. Following a wet compaction event that typically occurs when the stellar mass is $\sim10^{9.5}~M_{\odot}$ at $z\sim2-4$, the gas develops a cusp inside the effective radius, associated with a peak in star-formation rate (SFR). The SFR peak and the associated feedback, in the absence of further gas inflow to the centre, marks the onset of gas depletion from the central 1 kpc, leading to quenching of the central SFR. An extended, star-forming ring that forms by fresh gas during the central quenching process shows as a rising specific SFR (sSFR) profile, which is interpreted as inside-out quenching. Before quenching, the stellar density profile grows self-similarly, maintaining its log-log shape because the sSFR is similar at all radii. During the quenching process, the stellar density saturates to a constant value, especially in the inner 1 kpc. The stellar mass and SFR profiles deduced from observations show very similar shapes, consistent with the scenario of wet compaction leading to inside-out quenching and the subsequent saturation of a dense stellar core. We predict a cuspy gas profile during the blue nugget phase, and a gas-depleted core, sometimes surrounded by a ring, in the post-blue nugget phase.

Like many other late-type galaxies, the Milky Way contains a nuclear star cluster. In this work we obtain the basic properties of its dominant old stellar population. Firstly, we derive its structural properties by constructing a stellar surface density map of the central 1000" using extinction corrected star counts from VISTA, WFC3/IR and VLT/NACO data. We can describe the profile with a two-component model. The inner, slightly flattened (q=0.80+/-0.04) component is the nuclear cluster, while the outer component corresponds to the stellar component of the circumnuclear zone. We measure for the nuclear cluster a half-light radius of 178+/-51"~7+/-2 pc and a luminosity of M_Ks=-16.0+/-0.5. Secondly, we enlarge the range over which detailed dynamics are available from 1 pc to 4 pc. We obtain more than 10000 individual proper motions from NACO data, and more than 2500 radial velocities from VLT/SINFONI data. We determine the cluster mass by means of isotropic spherical Jeans modeling. We fix the distance to the Galactic Center and the mass of the supermassive black hole. We model the cluster either with a constant mass to light ratio or with a power law mass model. For the latter we obtain a slope of delta_M=1.18+/-0.06. Assuming spherical symmetry, we get a nuclear cluster mass within 100" of M_100"=6.09+/- 0.53_(fix_R_0) +/- 0.97_(R_0) * 10^6 M_sun for both modeling approaches. A model which includes the observed flattening gives a 47% larger mass, see Chatzopoulos et al. 2015. Our results slightly favor a core over a cusp in the mass profile. By minimizing the number of unbound stars within 8" in our sample we obtain a distance estimate of R$_0=8.53^{+0.21}_{-0.15}$ kpc, where an a priori relation between R_0 and SMBH mass from stellar orbits is used. Combining our mass and flux we obtain M/L=0.51+/-0.12 M_sun/L_(sun Ks). This is roughly consistent with a Chabrier IMF.

We study the relationship between stellar mass, star formation rate (SFR),ionization state, and gas-phase metallicity for a sample of 41 normal star-forming galaxies at $3 \lesssim z \lesssim 3.7$. The gas-phase oxygen abundance, ionization parameter, and electron density of ionized gas are derived from rest-frame optical strong emission lines measured on near-infrared spectra obtained with Keck/MOSFIRE. We remove the effect of these strong emission lines in the broad-band fluxes to compute stellar masses via spectral energy distribution fitting, while the SFR is derived from the dust-corrected ultraviolet luminosity. The ionization parameter is weakly correlated with the specific SFR, but otherwise the ionization parameter and electron density do not correlate with other global galaxy properties such as stellar mass, SFR, and metallicity. The mass-metallicity relation (MZR) at $z\simeq3.3$ shows lower metallicity by $\simeq 0.7$ dex than that at $z=0$ at the same stellar mass. Our sample shows an offset by $\simeq 0.3$ dex from the locally defined mass-metallicity-SFR relation, indicating that simply extrapolating such relation to higher redshift may predict an incorrect evolution of MZR. Furthermore, within the uncertainties we find no SFR-metallicity correlation, suggesting a less important role of SFR in controlling the metallicity at high redshift. We finally investigate the redshift evolution of the MZR by using the model by Lilly et al. (2013), finding that the observed evolution from $z=0$ to $z\simeq3.3$ can be accounted for by the model assuming a weak redshift evolution of the star formation efficiency.

Using cosmological simulations, we address the properties of high-redshift star-forming galaxies (SFGs) across their main sequence (MS) in the plane of star-formation rate (SFR) versus stellar mass. We relate them to the evolution of galaxies through phases of gas compaction, depletion, possible replenishment, and eventual quenching. We find that the high-SFR galaxies in the upper envelope of the MS are compact, with high gas fractions and short depletion times ("blue nuggets"), while the lower-SFR galaxies in the lower envelope have lower central gas densities, lower gas fractions and longer depletion times, consistent with observed gradients across the MS. Stellar-structure gradients are negligible. The SFGs oscillate about the MS ridge on timescales $\sim0.4~t_{\mathrm{Hubble}}$ ($\sim1$ Gyr at $z\sim3$). The propagation upwards is due to gas compaction, triggered, e.g., by mergers, counter-rotating streams, and/or violent disc instabilities. The downturn at the upper envelope is due to central gas depletion by peak star formation and outflows while inflow from the shrunken gas disc is suppressed. An upturn at the lower envelope can occur once the extended disc has been replenished by fresh gas and a new compaction can be triggered, namely as long as the replenishment time is shorter than the depletion time. The mechanisms of gas compaction, depletion and replenishment confine the SFGs to the narrow ($\pm0.3$ dex) MS. Full quenching occurs in massive haloes ($M_{\mathrm{vir}}>10^{11.5}~M_\odot$) and/or at low redshifts ($z<3$), where the replenishment time is long compared to the depletion time, explaining the observed bending down of the MS at the massive end.

We investigate the stellar population properties of a sample of 24 massive quenched galaxies at $1.25<z_\mathrm{spec}<2.09$ identified in the COSMOS field with our Subaru/MOIRCS near-IR spectroscopic observations. Tracing the stellar population properties as close to their major formation epoch as possible, we try to put constraints on the star formation history, post-quenching evolution, and possible progenitor star-forming populations for such massive quenched galaxies. By using a set of Lick absorption line indices on a rest-frame optical composite spectrum, the average age, metallicity [Z/H], and $\alpha$-to-iron element abundance ratio [$\alpha$/Fe] are derived as $\log(\mathrm{age}/\mathrm{Gyr})=0.04_{-0.08}^{+0.10}$, $\mathrm{[Z/H]}=0.24_{-0.14}^{+0.20}$, and $[\alpha/\mathrm{Fe}]=0.31_{-0.12}^{+0.12}$, respectively. If our sample of quenched galaxies at $\langle z \rangle = 1.6$ is evolved passively to $z=0$, their stellar population properties will align in excellent agreement with local counterparts at similar stellar velocity dispersions, which qualifies them as progenitors of local massive early-type galaxies. Redshift evolution of stellar population ages in quenched galaxies combined with low redshift measurements from the literature suggests a formation redshift of $z_\mathrm{f} \sim 2.3$ around which the bulk of stars in these galaxies have been formed. The measured [$\alpha$/Fe] value indicates a star formation timescale of $\lesssim 1$ Gyr, which can be translated into a specific star formation rate of $\simeq 1\,\mathrm{Gyr}^{-1}$ prior to quenching. Based on these findings, we discuss identifying possible progenitor star-forming galaxies at $z \simeq 2.3$. We identify normal star-forming galaxies, i.e, those on the star-forming main sequence, followed by a rapid quenching event, as likely precursors of the quenched galaxies at $\langle z \rangle = 1.6$ presented here.

Most present-day galaxies with stellar masses $\geq10^{11}$ solar masses show no ongoing star formation and are dense spheroids. Ten billion years ago, similarly massive galaxies were typically forming stars at rates of hundreds solar masses per year. It is debated how star formation ceased, on which timescales, and how this "quenching" relates to the emergence of dense spheroids. We measured stellar mass and star-formation rate surface density distributions in star-forming galaxies at redshift 2.2 with $\sim1$ kiloparsec resolution. We find that, in the most massive galaxies, star formation is quenched from the inside out, on timescales less than 1 billion years in the inner regions, up to a few billion years in the outer disks. These galaxies sustain high star-formation activity at large radii, while hosting fully grown and already quenched bulges in their cores.

We present the analysis of HST $J$- and $H$-band imaging for 29 galaxies on the star-forming main sequence at $z\sim2$, which have Adaptive Optics VLT SINFONI integral field spectroscopy from our SINS/zC-SINF program. The SINFONI H$\alpha$ data resolve the on-going star-formation and the ionized gas kinematics on scales of $1-2$ kpc; the near-IR images trace the galaxies' rest-frame optical morphologies and distributions of stellar mass in old stellar populations at a similar resolution. The global light profiles of most galaxies show disk-like properties well described by a single S\'ersic profile with $n\sim1$, with only $\sim15%$ requiring a high $n>3$ S\'ersic index, all more massive than $10^{10}M_\odot$. In bulge+disk fits, about $40%$ of galaxies have a measurable bulge component in the light profiles, with $\sim15%$ showing a substantial bulge-to-total ratio $B/T\ge0.3$. This is a lower limit to the frequency of $z\sim2$ massive galaxies with a developed bulge component in stellar mass because it could be hidden by dust and/or outshined by a thick actively star-forming disk component. The galaxies' rest-optical half-light radii range between $1-7$ kpc, with a median of 2.1 kpc, and lie slightly above the size-mass relation at these epochs reported in the literature. This is attributed to differences in sample selection and definitions of size and/or mass measurements. The $(u-g)_{rest}$ color gradient and scatter within individual $z\sim2$ massive galaxies with $\ge10^{11}M_\odot$ are as high as in $z=0$ low-mass, late-type galaxies, and are consistent with the high star-formation rates of massive $z\sim2$ galaxies being sustained at large galactocentric distances.

We use the Zurich ENvironmental Study (ZENS) database to investigate the environmental dependence of the merger fraction $\Gamma$ and merging galaxy properties in a sample of ~1300 group galaxies with $M>10^{9.2}M_\odot$ and 0.05<z<0.0585. In all galaxy mass bins investigated in our study, we find that $\Gamma$ decreases by a factor of ~2-3 in groups with halo masses $M_{HALO}>10^{13.5} M_\odot$ relative to less massive systems, indicating a suppression of merger activity in large potential wells. In the fiducial case of relaxed groups only, we measure a variation $\Delta\Gamma/\Delta \log (M_{HALO}) \sim - 0.07$ dex$^{-1}$, which is almost independent of galaxy mass and merger stage. At galaxy masses $>10^{10.2} M_\odot$, most mergers are dry accretions of quenched satellites onto quenched centrals, leading to a strong increase of $\Gamma$ with decreasing group-centric distance at these mass scales.Both satellite and central galaxies in these high mass mergers do not differ in color and structural properties from a control sample of nonmerging galaxies of equal mass and rank. At galaxy masses $<10^{10.2} M_\odot$, where we mostly probe satellite-satellite pairs and mergers between star-forming systems, close pairs (projected distance $<10-20$ kpc) show instead $\sim2\times$ enhanced (specific) star formation rates and $\sim1.5\times$ larger sizes than similar mass, nonmerging satellites. The increase in both size and SFR leads to similar surface star-formation densities in the merging and control-sample satellite populations.

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field.

In this paper we follow up on our previous detection of nuclear ionized outflows in the most massive (log(M*/Msun) >= 10.9) z~1-3 star-forming galaxies (Forster Schreiber et al.), by increasing the sample size by a factor of six (to 44 galaxies above log(M*/Msun) >= 10.9) from a combination of the SINS/zC-SINF, LUCI, GNIRS, and KMOS^3D spectroscopic surveys. We find a fairly sharp onset of the incidence of broad nuclear emission (FWHM in the Ha, [NII], and [SII] lines ~ 450-5300 km/s), with large [NII]/Ha ratios, above log(M*/Msun) ~ 10.9, with about two thirds of the galaxies in this mass range exhibiting this component. Broad nuclear components near and above the Schechter mass are similarly prevalent above and below the main sequence of star-forming galaxies, and at z~1 and ~2. The line ratios of the nuclear component are fit by excitation from active galactic nuclei (AGN), or by a combination of shocks and photoionization. The incidence of the most massive galaxies with broad nuclear components is at least as large as that of AGNs identified by X-ray, optical, infrared or radio indicators. The mass loading of the nuclear outflows is near unity. Our findings provide compelling evidence for powerful, high-duty cycle, AGN-driven outflows near the Schechter mass, and acting across the peak of cosmic galaxy formation.

We present the correlations between stellar mass, star formation rate (SFR) and [NII]/Ha flux ratio as indicator of gas-phase metallicity for a sample of 222 galaxies at 0.8 < z < 2.6 and log(M*/Msun)=9.0-11.5 from the LUCI, SINS/zC-SINF and KMOS3D surveys. This sample provides a unique analysis of the mass-metallicity relation (MZR) over an extended redshift range using consistent data analysis techniques and strong-line metallicity indicator. We find a constant slope at the low-mass end of the relation and can fully describe its redshift evolution through the evolution of the characteristic turnover mass where the relation begins to flatten at the asymptotic metallicity. At fixed mass and redshift, our data do not show a correlation between the [NII]/Ha ratio and SFR, which disagrees with the 0.2-0.3dex offset in [NII]/Ha predicted by the "fundamental relation" between stellar mass, SFR and metallicity discussed in recent literature. However, the overall evolution towards lower [NII]/Ha at earlier times does broadly agree with these predictions.

We report the detection of ubiquitous powerful nuclear outflows in massive (> 10^11 Msun) z~2 star-forming galaxies (SFGs), which are plausibly driven by an Active Galactic Nucleus (AGN). The sample consists of the eight most massive SFGs from our SINS/zC-SINF survey of galaxy kinematics with the imaging spectrometer SINFONI, six of which have sensitive high-resolution adaptive optics (AO) assisted observations. All of the objects are disks hosting a significant stellar bulge. The spectra in their central regions exhibit a broad component in Halpha and forbidden [NII] and [SII] line emission, with typical velocity FWHM ~ 1500 km/s, [NII]/Halpha ratio ~ 0.6, and intrinsic extent of 2 - 3 kpc. These properties are consistent with warm ionized gas outflows associated with Type 2 AGN, the presence of which is confirmed via independent diagnostics in half the galaxies. The data imply a median ionized gas mass outflow rate of ~ 60 Msun/yr and mass loading of ~ 3. At larger radii, a weaker broad component is detected but with lower FWHM ~ 485 km/s and [NII]/Halpha ~ 0.35, characteristic for star formation-driven outflows as found in the lower-mass SINS/zC-SINF galaxies. The high inferred mass outflow rates and frequent occurrence suggest the nuclear outflows efficiently expel gas out of the centers of the galaxies with high duty cycles, and may thus contribute to the process of star formation quenching in massive galaxies. Larger samples at high masses will be crucial to confirm the importance and energetics of the nuclear outflow phenomenon, and its connection to AGN activity and bulge growth.

Here we present the fundamental properties of the nuclear cluster of the Milky Way. First, we derive its structural properties by constructing a density map of the central 1000'' using extinction-corrected star counts. We can describe the data with a two-component model built from Sersic profiles. The inner nearly spherical component is the nuclear cluster. The outer, strongly flattened component can be identified with the stellar component of the circumnuclear zone. Second, we enlarge the radius inside which detailed dynamics are available from $1\,$pc to $4\,$pc. We use more than 10000 individual proper motions and more than 2700 radial velocities. We determine the cluster mass by means of isotropic spherical Jeans modeling. We get a nuclear cluster mass within 100'' of $M_{100''}=(6.11 \pm 0.52|_{\mathrm{fix} R_0}\pm 0.97|_{R_0}) \times 10^6 $ M$_{\odot}$, which corresponds to a total cluster mass of M$_{\mathrm{NC}}=(13.08 \pm 2.51|_{\mathrm{fix} R_0}\pm 2.08|_{R_0}) \times 10^6 $ M$_{\odot}$. By combination of our mass with the flux we calculate $M/L=0.50 \pm 0.12 M_{\odot}/L_{\odot,\mathrm{Ks}}$ for the central 100''. That is broadly consistent with a Chabrier IMF. With its mass and a luminosity of M$_{\mathrm{Ks}}=-15.30\pm0.26$ the nuclear cluster is a bright and massive specimen with a typical size.

As part of the SINS/zC-SINF surveys of high-z galaxy kinematics, we derive the radial distributions of H-alpha surface brightness, stellar mass surface density, and dynamical mass at ~2 kpc resolution in 19 z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. From these data we infer the radial distribution of the Toomre Q-parameter for these main-sequence star forming galaxies (SFGs), covering almost two decades of stellar mass (10^9.6 to 10^11.5 solar masses). In more than half of our SFGs, the H-alpha distributions cannot be fit by a centrally peaked distribution, such as an exponential, but are better described by a ring, or the combination of a ring and an exponential. At the same time the kinematic data indicate the presence of a mass distribution more centrally concentrated than a single exponential distribution for 5 of the 19 galaxies. The resulting Q-distributions are centrally peaked for all, and significantly exceed unity there for three quarters of the SFGs. The occurrence of H-alpha rings and of large nuclear Q-values is strongly correlated, and is more common for the more massive SFGs. While our sample is small and there remain substantial uncertainties and caveats, our observations are consistent with a scenario in which cloud fragmentation and global star formation are secularly suppressed in gas rich high-z disks from the inside out, as the central stellar mass density of the disks grows.

Gamma Ray Bursts (GRBs) and galaxies at high redshift represent complementary probes of the star formation history of the Universe. In fact, both the GRB rate and the galaxy luminosity density are connected to the underlying star formation. Here, we combine a star formation model for the evolution of the galaxy luminosity function from z=0 to z=10 with a metallicity-dependent efficiency for GRB formation to simultaneously predict the comoving GRB rate. Our model sheds light on the physical origin of the empirical relation often assumed between GRB rate and luminosity density-derived star formation rate: Rgrb(z) = \epsilon(z)*SFR_{obs}(z), with \epsilon(z) (1+z)^{1.2}. At z<4, \epsilon(z) is dominated by the effects of metallicity evolution in the GRB efficiency. Our best-fitting model only requires a moderate preference for low-metallicity, that is a GRB rate per unit stellar mass about four times higher for log(Z/Zsun)<-3 compared to log(Z/Zsun)>0. Models with total suppression of GRB formation at log(Z/Zsun)>0 are disfavored. At z>4, most of the star formation happens in low-metallicity hosts with nearly saturated efficiency of GRB production per unit stellar mass. However at the same epoch, galaxy surveys miss an increasing fraction of the predicted luminosity density because of flux limits, driving an accelerated evolution of \epsilon(z) compared to the empirical power-law fit from lower z. Our findings are consistent with the non-detections of GRB hosts in ultradeep imaging at z>5, and point toward current galaxy surveys at z>8 only observing the top 15-20 % of the total luminosity density.

We report the results of deep (4.6h) H band spectroscopy of the well studied z~12 H-band dropout galaxy candidate UDFj-39546284 with MOSFIRE on Keck-I. These data reach a sensitivity of 5-10 x 10^-19 erg s^-1 cm^-2 per 4.4 Angstrom resolution element between sky lines. Previous papers have argued this source could either be a large equivalent width line emitting galaxy at 2<z<3.5 or a luminous galaxy at z~12. We find a 2.2 sigma peak associated with a line candidate in deep Hubble-Space-Telescope Wide-Field-Camera 3 Infrared grism observations, but at a lower flux than what was expected. After considering several possibilities we conclude these data can not conclusively confirm or reject the previous line detection, and significantly deeper spectroscopic observations are required. We also search for low-redshift emission lines in ten other 7<z<10 z, Y, and J-dropout candidates in our mask and find no significant detections.

We present a model to understand the redshift evolution of the UV luminosity and stellar mass functions of Lyman Break Galaxies. Our approach is based on the assumption that the luminosity and stellar mass of a galaxy is related to its dark matter halo assembly and gas infall rate. Specifically, galaxies experience a burst of star formation at the halo assembly time, followed by a constant star formation rate, representing a secular star formation activity sustained by steady gas accretion. Star formation from steady gas accretion is the dominant contribution to the galaxy UV luminosity at all redshifts. The model is calibrated by constructing a galaxy luminosity versus halo mass relation at $z=4$ via abundance matching. After this luminosity calibration, the model naturally fits the $z=4$ stellar mass function, and correctly predicts the evolution of both luminosity and stellar mass functions from $z=0$ to $z=8$. While the details of star formation efficiency and feedback are hidden within our calibrated luminosity versus halo mass relation, our study highlights that the primary driver of galaxy evolution across cosmic time is the build-up of dark matter halos, without the need to invoke a redshift dependent efficiency in converting gas into stars.