Researcher profile

Sandro Tacchella

Sandro Tacchella contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
32works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2026arXiv

Abundant Population of Broad H$α$ Emitters in the GOODS-N Field Revealed by CONGRESS, FRESCO, and JADES

We present a spectroscopic search for broad H$α$ emitters at z$\approx$3.7-6.5 in the GOODS-N field, utilizing JWST/NIRCam slitless spectroscopy from FRESCO and CONGRESS, complemented by JADES imaging. We identify 19 broad H$α$ emitters with FWHM$>$1000 km/s at z$\approx$4-5.5, including 9 new sources. The broad H$α$ luminosity function (LF) derived from our sample is consistent with those of other JWST-selected broad-line AGN reported in the literature. The black hole masses and AGN bolometric luminosities, inferred from the broad H$α$ components, indicate that most sources are accreting at ~10% of the Eddington limit. We derive their host stellar masses via SED fitting and find higher $M_{BH}/M_{*}$ ratios relative to the local $M_{BH}-M_{*}$ relations, consistent with previous studies. We find that 42% of the sample do not satisfy the widely-used color selection criteria for Little Red Dots (LRDs), with the majority of these sources lacking the characteristic steep rest-optical red slope, indicating that the LRD selection is highly incomplete when selecting AGN galaxies. A comparison of the average SEDs between our sample and LRDs selected in the same field reveals that the steep red slopes observed in some LRDs are likely due to line-boosting effects as previously suggested. Furthermore, we find that 68% of color-selected LRDs with H$α$ detections in the NIRCam/Grism spectra do not exhibit broad-line features. While the limited sensitivity of the grism spectra may hinder the detection of broad-line components in faint sources, our findings still highlight the enigmatic nature of the LRD population.

preprint2026arXiv

Accelerating Redshift-Conditioned Galaxy Image Synthesis with One-step Generative Modeling

Understanding galaxy morphology evolution across cosmic time requires models that can generate realistic galaxy populations conditioned on redshift. In this work, we study efficient redshift-conditioned generative modeling for astrophysical image synthesis using diffusion models and pixel-MeanFlow. We first review the connections between score-based diffusion models, Flow Matching, one-step generative models, and modern diffusion samplers. We then evaluate DDPM, DDIM, DEIS-AB2, DPM++2M, and one-step pixel-MeanFlow on the GalaxiesML-64 dataset using morphology-based metrics, including ellipticity, semi-major axis, Sérsic index, and isophotal area. Our results show a clear accuracy-efficiency trade-off: standard DDPM sampling achieves the best distributional fidelity but requires high computational cost, while second-order samplers substantially improve efficiency over DDIM. Pixel-MeanFlow enables single-step generation and achieves competitive performance on several morphology statistics, though it remains weaker than many-step DDPM for fine-grained structure. Our results demonstrate that one-step generative models can recover key galaxy morphology statistics at orders-of-magnitude lower computational cost, opening a path toward efficient conditional simulators for large cosmological surveys and simulation-based scientific inference.

preprint2026arXiv

JADES: the chemical enrichment pattern of distant galaxies -- $α$ enhancement, silicon depletion, and iron enhancement

We present gas-phase abundances of carbon (C), $α$-elements (O, Ne, Si, and Ar) and iron (Fe) obtained from stacked spectra of high-$z$ star-forming galaxies with the deep Near Infrared Spectrograph medium-resolution data from the James Webb Space Telescope Advanced Deep Extragalactic Survey. Our 564 sources at $z=4$--7 have a median stellar mass of $\log(M_{*}/M_{\odot})=8.46$ and a median star-formation rate of $\log(\mathrm{SFR}/M_{\odot}\,\mathrm{yr^{-1}})=0.30$, placing them close to the star-formation main sequence. We find that the stacked spectrum of all our 564 sources has relatively low [C/O]$=-0.70$, moderate [Ne/O]$=-0.09$, and low [Ar/O]$=-0.28$ values at a low gas-phase metallicity of $12+\log(\mathrm{O/H})=7.71$ ($Z\sim 0.1~Z_\odot$), suggesting dominant yields of core-collapse supernovae evolved from massive stars. The detection of a weak SiIII] emission line in our stacked spectrum provides a silicon-to-oxygen abundance ratio of [Si/O]$=-0.63$, which is lower than that of stars in the Milky Way disc and lower than expected by chemical evolution models, suggesting silicon depletion onto dust grains. Likewise, this Si/O value is lower than that we newly derive for two individual $z>6$ galaxies (GN-z11 and RXCJ2248) with negligible dust attenuation. By performing spectral stacking in bins of $M_{*}$, SFR, specific SFR (sSFR), and ultra-violet (UV) continuum slope $β_{\mathrm{UV}}$, we identify [FeIII] line detections in the high-sSFR bin and the blue-$β_{\mathrm{UV}}$ bin, both of which exhibit supersolar Fe/O ratios, while their C/O, Ar/O, and Si/O ratios are comparable to those of the all-sources stack. Our findings support a chemically young gas composition with rapid dust depletion in the general population of high-$z$ star-forming galaxies, while raising the possibility of anomalous, selective Fe/O enhancement at the very early epoch of star formation.

preprint2026arXiv

On the origins of oxygen: ALMA and JWST characterise the multi-phase, metal-enriched, star-bursting medium within a 'normal' $z > 11$ galaxy

The unexpectedly high abundance of galaxies at $z > 11$ revealed by JWST has sparked a debate on the nature of early galaxies and the physical mechanisms regulating their formation. The Atacama Large Millimeter/submillimeter Array (ALMA) has begun to provide vital insights on their gas and dust content, but so far only for extreme &#39;blue monsters&#39;. Here we present new, deep ALMA observations of JADES-GS-z11-0, a more typical (sub-$L^*$) $z > 11$ galaxy that bridges the discovery space of JWST and the Hubble Space Telescope. These data confirm the presence of the [O III] 88 $μ$m line at $4.5σ$ significance, precisely at the redshift of several faint emission lines previously seen with JWST/NIRSpec, while the underlying dust continuum remains undetected ($F_ν< 9.0 \, \mathrm{μJy}$), implying an obscured star formation rate (SFR) of $\text{SFR}_\text{IR} \lesssim 6 \, \mathrm{M_\odot \, yr^{-1}}$ and dust mass of $M_\text{dust} \lesssim 1.0 \times 10^{6} \, \mathrm{M_\odot}$ (all $3σ$). The accurate ALMA redshift of $z_\text{[O III]} = 11.1221 \pm 0.0006$ ($\gtrsim \! 5\times$ refined over NIRSpec) helps confirm that redshifts measured purely from the Lyman-$α$ break, even spectroscopically, should properly take into account the effects of potential damped Lyman-$α$ absorption (DLA) systems to avoid systematic overestimates of up to $Δz \approx 0.5$. The [O III] 88 $μ$m luminosity of $L_\text{[O III]} = (1.1 \pm 0.3) \times 10^{8} \, \mathrm{L_\odot}$, meanwhile, agrees well with the scaling relation for local metal-poor dwarfs given the SFR measured by NIRCam, NIRSpec, and MIRI. The spatially resolved MIRI and ALMA emission also underscores that JADES-GS-z11-0 is likely to consist of two low-mass components that are undergoing strong bursts of star formation yet are already pre-enriched in oxygen (~20-30% solar), only 400 Myr after the Big Bang.

preprint2025arXiv

The THESAN-ZOOM project: Star-formation efficiencies in high-redshift galaxies

Recent JWST observations hint at unexpectedly intense cosmic star-formation in the early Universe, often attributed to enhanced star-formation efficiencies (SFEs). Here, we analyze the SFE in THESAN-ZOOM, a novel zoom-in radiation-hydrodynamic simulation campaign of high-redshift ($z \gtrsim 3$) galaxies employing a state-of-the-art galaxy formation model resolving the multiphase interstellar medium (ISM). The halo-scale SFE ($ε^{\ast}_{\rm halo}$) - the fraction of baryons accreted by a halo that are converted to stars - follows a double power-law dependence on halo mass, with a mild redshift evolution above $M_{\rm halo} \gtrsim 10^{9.5}\,{\rm M}_{\odot}$. The power-law slope is roughly $1/3$ at large halo masses, consistent with expectations when gas outflows are momentum-driven. At lower masses, the slope is roughly $2/3$ and is more aligned with the energy-driven outflow scenario. $ε^{\ast}_{\rm halo}$ is a factor of $2-3$ larger than commonly assumed in empirical galaxy-formation models at $M_{\rm halo} \lesssim 10^{11}\,{\rm M}_{\odot}$. On galactic (kpc) scales, the Kennicutt-Schmidt (KS) relation of neutral gas is universal in THESAN-ZOOM, following $Σ_{\rm SFR} \propto Σ_{\rm gas}^2$, indicative of a turbulent energy balance in the ISM maintained by stellar feedback. The rise of $ε^{\ast}_{\rm halo}$ with halo mass can be traced primarily to increasing gas surface densities in massive galaxies, while the underlying KS relation and neutral, star-forming gas fraction remain unchanged. Although the increase in $ε^{\ast}_{\rm halo}$ with redshift is relatively modest, it is sufficient to explain the large observed number density of UV-bright galaxies at $z \gtrsim 12$. However, reproducing the brightest sources at $M_{\rm UV} \lesssim -21$ may require extrapolating the SFE beyond the halo mass range directly covered by THESAN-ZOOM.

preprint2024arXiv

Inside the bubble: exploring the environments of reionisation-era Lyman-$α$ emitting galaxies with JADES and FRESCO

We present a study of the environments of 17 Lyman-$α$ (Ly$α$) emitting galaxies (LAEs) in the reionisation era ($5.8 < z < 8$) identified by JWST/NIRSpec as part of the JWST Advanced Deep Extragalactic Survey (JADES). Unless situated in sufficiently (re)ionised regions, Ly$α$ emission from these galaxies would be strongly absorbed by neutral gas in the intergalactic medium (IGM). We conservatively estimate sizes of the ionised regions required to reconcile the relatively low Ly$α$ velocity offsets ($Δv_\text{Ly$α$}<300\,\mathrm{km\,s^{-1}}$) with moderately high Ly$α$ escape fractions ($f_\mathrm{esc,\,Lyα}>5\%$) observed in our sample of LAEs, suggesting the presence of ionised hydrogen along the line of sight towards at least eight out of 17 LAEs. We find minimum physical `bubble&#39; sizes of the order of $R_\text{ion}\sim0.1$-$1\,\mathrm{pMpc}$ are required in a patchy reionisation scenario where ionised bubbles containing the LAEs are embedded in a fully neutral IGM. Around half of the LAEs in our sample are found to coincide with large-scale galaxy overdensities seen in FRESCO at $z \sim 5.8$-$5.9$ and $z\sim7.3$, suggesting Ly$α$ transmission is strongly enhanced in such overdense regions, and underlining the importance of LAEs as tracers of the first large-scale ionised bubbles. Considering only spectroscopically confirmed galaxies, we find our sample of UV-faint LAEs ($M_\text{UV}\gtrsim-20\,\mathrm{mag}$) and their direct neighbours are generally not able to produce the required ionised regions based on the Ly$α$ transmission properties, suggesting lower-luminosity sources likely play an important role in carving out these bubbles. These observations demonstrate the combined power of JWST multi-object and slitless spectroscopy in acquiring a unique view of the early Universe during cosmic reionisation via the most distant LAEs.

preprint2023arXiv

Molecular Gas Reservoirs in Massive Quiescent Galaxies at $\mathrm{z\sim0.7}$ Linked to Late Time Star Formation

We explore how the presence of detectable molecular gas depends on the inferred star formation histories (SFHs) in 8 massive, quiescent galaxies at $\mathrm{z\sim0.7}$. Half of the sample have clear detections of molecular gas, traced by CO(2-1). We find that the molecular gas content is unrelated to the rate of star formation decline prior to the most recent 1 Gyr, suggesting that the gas reservoirs are not leftover from their primary star formation epoch. However, the recent SFHs of CO-detected galaxies demonstrate evidence for secondary bursts of star formation in their last Gyr. The fraction of stellar mass formed in these secondary bursts ranges from $\mathrm{f_{burst}\approx0.3-6\%}$, and ended between $\mathrm{t_{end\mbox{-}burst}\approx0-330~Myr}$ ago. The CO-detected galaxies form a higher fraction of mass in the last Gyr ($\mathrm{f_{M_{1Gyr}}=2.6\pm1.8\%}$) compared to the CO-undetected galaxies ($\mathrm{f_{M_{1Gyr}}=0.2\pm0.1\%}$). The galaxies with gas reservoirs have enhanced late-time star formation, highlighting this as a contributing factor to the observed heterogeneity in the gas reservoirs in high-redshift quiescent galaxies. We find that the amount of gas and star formation driven by these secondary bursts are inconsistent with that expected from dry minor mergers, and instead are likely driven by recently-accreted gas i.e., gas-rich minor mergers. This conclusion would not have been made based on $\mathrm{SFR_{UV+IR}}$ measurements alone, highlighting the power of detailed SFH modeling in the interpretation of gas reservoirs. Larger samples are needed to understand the frequency of low-level rejuvenation among quiescent galaxies at intermediate redshifts, and to what extent this drives the diversity of molecular gas reservoirs.

preprint2022arXiv

(Re)Solving Reionization with Lyα: How Bright Lyα Emitters account for the $z\approx2-8$ Cosmic Ionizing Background

The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here we introduce an emissivity framework based on Ly$α$ emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between the production and escape of Ly$α$ and LyC photons. Using constraints on the escape fractions of bright LAEs ($L_{\rm{Lyα}}>0.2 L^{*}$) at $z\approx2$ obtained from resolved Ly$α$ profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at $z\approx2-6$, which is non-trivial given the strong evolution in both the star-formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between $z\approx6-9$ under plausible assumptions. Within this framework, the $>10\times$ rise in the UV population-averaged $f_{\rm{esc}}$ between $z\approx3-7$ naturally arises due to the same phenomena that drive the growing Ly$α$ emitter fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters ($f_{\rm{esc}}=50\%$, $ξ_{\rm{ion}}=10^{25.9}$ Hz erg$^{-1}$ for half the bright LAEs), a highly ionizing minority of galaxies with $M_{\rm UV}<-17$ accounts for the entire ionizing budget from star-forming galaxies. Rapid flashes of LyC from such rare galaxies produce a &#34;disco&#34; ionizing background. We conclude proposing tests to further develop our suggested Ly$α$-anchored formalism.

preprint2022arXiv

Dissecting Nearby Galaxies with piXedfit: II. Spatially Resolved Scaling Relations Among Stars, Dust, and Gas

We study spatially resolved scaling relations among stars, dust, and gas in ten nearby spiral galaxies. In a preceding paper Abdurro&#39;uf et al. (2022), we have derived spatially resolved properties of the stellar population and dust by panchromatic spectral energy distribution (SED) fitting using piXedfit. Now, we investigate resolved star formation ($Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$) and dust scaling relations. While the relations with all sub-galactic regions of the galaxies are reasonably tight ($σ\lesssim 0.3$ dex), we find that most of the scaling relations exhibit galaxy-to-galaxy variations in normalization and shape. Only two relations of $Σ_{\rm dust}$--$Σ_{\rm gas}$ and $Σ_{\rm dust}$--$Σ_{\rm H_{2}}$ do not show noticeable galaxy-to-galaxy variations among our sample galaxies. We further investigate correlations among the scaling relations. We find significant correlations among the normalization of the $Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$ relations, which suggest that galaxies with higher levels of resolved $\text{H}_{2}$ fraction ($f_{\rm H_{2}}$) tend to have higher levels of resolved star formation efficiency (SFE) and specific star formation rate (sSFR). We also observe that galaxies with higher levels of resolved dust-to-stellar mass ratios tend to have higher levels of resolved sSFR, SFE, and $f_{\rm H_{2}}$. Moreover, we find that galaxies with higher global sSFR and less compact morphology tend to have higher levels of the resolved sSFR, SFE, and $f_{\rm H_{2}}$, which can explain the variations in the normalization of the $Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$ relationships. Overall, we observe indications of the contributions of both global and local factors in governing the star formation process in galaxies.

preprint2022arXiv

Fast, Slow, Early, Late: Quenching Massive Galaxies at z~0.8

We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at $\langle z_{\rm obs} \rangle=0.8$ with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector, we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including non-parametric star-formation histories, emission lines, variable dust attenuation law, and dust and AGN emission) together with an uncertainty and outlier model. We show that both spectroscopy and photometry are needed to break the dust-age-metallicity degeneracy. We find a large diversity of star-formation histories: although the most massive ($M_{\star}>2\times10^{11}~M_{\odot}$) galaxies formed the earliest (formation redshift of $z_{\rm f}\approx5-10$ with a short star-formation timescale of $τ_{\rm SF}\lesssim1~\mathrm{Gyr}$), lower-mass galaxies have a wide range of formation redshifts, leading to only a weak trend of $z_{\rm f}$ with $M_{\star}$. Interestingly, several low-mass galaxies with have formation redshifts of $z_{\rm f}\approx5-8$. Star-forming galaxies evolve about the star-forming main sequence, crossing the ridgeline several times in their past. Quiescent galaxies show a wide range and continuous distribution of quenching timescales ($τ_{\rm quench}\approx0-5~\mathrm{Gyr}$) with a median of $\langleτ_{\rm quench}\rangle=1.0_{-0.9}^{+0.8}~\mathrm{Gyr}$ and of quenching epochs of $z_{\rm quench}\approx0.8-5.0$ ($\langle z_{\rm quench}\rangle=1.3_{-0.4}^{+0.7}$). This large diversity of quenching timescales and epochs points toward a combination of internal and external quenching mechanisms. In our sample, rejuvenation and &#34;late bloomers&#34; are uncommon. In summary, our analysis supports the &#34;grow & quench&#34; framework and is consistent with a wide and continuously-populated diversity of quenching timescales.

preprint2022arXiv

H-alpha emission in local galaxies: star formation, time variability and the diffuse ionized gas

The nebular recombination line H$α$ is widely used as a star-formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H$α$ radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H$α$ emission, and its connection to the underlying gas and star formation properties. The H$α$ and H$β$ radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H$α$ emission from collisional excitation amounts to $f_{\rm col}\sim5-10\%$, only weakly dependent on radius and vertical height, and that scattering boosts the H$α$ luminosity by $\sim40\%$. The dust correction via the Balmer decrement works well (intrinsic H$α$ emission recoverable within $25\%$), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H$α$-SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $f_{\rm abs}\approx28\%$ and $f_{\rm He}\approx9\%$, respectively. Together with an escape fraction of $f_{\rm esc}\approx6\%$, this reduces the available budget for hydrogen line emission by nearly half ($f_{\rm H}\approx57\%$). We discuss the impact of the diffuse ionized gas, showing - among other things - that the extraplanar H$α$ emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope.

preprint2022arXiv

Identifying Galaxy Mergers in Simulated CEERS NIRCam Images using Random Forests

Identifying merging galaxies is an important - but difficult - step in galaxy evolution studies. We present random forest classifications of galaxy mergers from simulated JWST images based on various standard morphological parameters. We describe (a) constructing the simulated images from IllustrisTNG and the Santa Cruz SAM, and modifying them to mimic future CEERS observations as well as nearly noiseless observations, (b) measuring morphological parameters from these images, and (c) constructing and training the random forests using the merger history information for the simulated galaxies available from IllustrisTNG. The random forests correctly classify $\sim60\%$ of non-merging and merging galaxies across $0.5 < z < 4.0$. Rest-frame asymmetry parameters appear more important for lower redshift merger classifications, while rest-frame bulge and clump parameters appear more important for higher redshift classifications. Adjusting the classification probability threshold does not improve the performance of the forests. Finally, the shape and slope of the resulting merger fraction and merger rate derived from the random forest classifications match with theoretical Illustris predictions, but are underestimated by a factor of $\sim 0.5$.

preprint2022arXiv

No Evolution in the Half-mass Radius of Milky Way-type Galaxies over the Last 10 Gyr

The Milky Way (MW) galaxy is in focus, thanks to new observational data. Here we shed new light on the MW&#39;s past by studying the structural evolution of MW progenitors, which we identify from extragalactic surveys. Specifically, we constrain the stellar-mass growth history (SMGH) of the MW with two methods: ($i$) direct measurement of the MW&#39;s star formation history, and ($ii$) assuming the MW is a typical star-forming galaxy that remains on the star-forming main sequence. We select MW progenitors based on these two SMGHs at $z=0.2-2.0$ from the CANDELS/3D-HST data. We estimate the structural parameters (including half-mass radius $r_{50}$ and Sérsic index) from the stellar-mass profiles. Our key finding is that the progenitors of the MW galaxy grow self-similarly on spatially resolved scales with roughly a constant half-mass radius ($\sim2-3$ kpc) over the past 10 Gyr, while their stellar masses increase by about 1 dex, implying little-to-no inside-out growth. We discover that the radius containing $20\%$ of the stellar mass ($r_{20}$) decreases by $60\%$ between redshifts of $z=2.0$ and $z=0.7$, while the central stellar-mass density ($Σ_1$) increases by a factor of 1.3 dex over the same time and the Sérsic index changes as $n\propto(1+z)^{-1.41\pm0.19}$. This is consistent with an early ($z>1$) formation of a thick disk, followed by the formation of a bar that led to an increase in the mass in the core. The formation and evolution of the thin disk had only little impact on the overall half-mass size. We also show that the constant-size evolution of the MW progenitors challenges semiempirical approaches and numerical simulations.

preprint2022arXiv

On the Stellar Populations of Galaxies at z=9-11: The Growth of Metals and Stellar Mass at Early Times

We present a detailed stellar population analysis of 11 bright ($H<26.6$) galaxies at $z=9-11$ (three spectroscopically confirmed) to constrain the chemical enrichment and growth of stellar mass of early galaxies. We use the flexible Bayesian spectral energy distribution (SED) fitting code Prospector with a range of star-formation histories (SFHs), a flexible dust attenuation law and a self-consistent modeling of emission lines. This approach allows us to assess how different priors affect our results, and how well we can break degeneracies between dust attenuation, stellar ages, metallicity and emission lines using data which probe only the rest-frame ultraviolet to optical wavelengths. We measure a median observed ultraviolet spectral slope $β=-1.87_{-0.43}^{+0.35}$ for relatively massive star-forming galaxies ($9<\log(M_{\star}/M_{\odot})<10$), consistent with no change from $z=4$ to $z=9-10$ at these stellar masses, implying rapid enrichment. Our SED-fitting results are consistent with a star-forming main sequence with sub-linear slope ($0.7\pm0.2$) and specific star-formation rates of $3-10~\mathrm{Gyr}^{-1}$. However, the stellar ages and SFHs are less well constrained. Using different SFH priors, we cannot distinguish between median mass-weighted ages of $\sim50-150$ Myr, which corresponds to 50\% formation redshifts of $z_{50}\sim10-12$ at $z\sim9$ and is of the order of the dynamical timescales of these systems. Importantly, the models with different SFH priors are able to fit the data equally well. We conclude that the current observational data cannot tightly constrain the mass-buildup timescales of these $z=9-11$ galaxies, with our results consistent with SFHs implying both a shallow and steep increase of the cosmic SFR density with time at $z>10$.

preprint2021arXiv

A Census of the Bright z=8.5-11 Universe with the Hubble and Spitzer Space Telescopes in the CANDELS Fields

We present the results from a new search for candidate galaxies at z ~ 8.5-11 discovered over the 850 arcmin^2 area probed by the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). We use a photometric redshift selection including both Hubble and Spitzer Space Telescope photometry to robustly identify galaxies in this epoch at F160W < 26.6. We use a detailed vetting procedure, including screening for persistence, stellar contamination, inclusion of ground-based imaging, and followup space-based imaging to build a robust sample of 11 candidate galaxies, three presented here for the first time. The inclusion of Spitzer/IRAC photometry in the selection process reduces contamination, and yields more robust redshift estimates than Hubble alone. We constrain the evolution of the rest-frame ultraviolet luminosity function via a new method of calculating the observed number densities without choosing a prior magnitude bin size. We find that the abundance at our brightest probed luminosities (M_UV=-22.3) is consistent with predictions from simulations which assume that galaxies in this epoch have gas depletion times at least as short as those in nearby starburst galaxies. Due to large Poisson and cosmic variance uncertainties we cannot conclusively rule out either a smooth evolution of the luminosity function continued from z=4-8, or an accelerate decline at z > 8. We calculate that the presence of seven galaxies in a single field (EGS) is an outlier at the 2-sigma significance level, implying the discovery of a significant overdensity. These scenarios will be imminently testable to high confidence within the first year of observations of the James Webb Space Telescope.

preprint2021arXiv

Dissecting Nearby Galaxies with piXedfit: I. Spatially Resolved Properties of Stars, Dust, and Gas as Revealed by Panchromatic SED Fitting

We study spatially resolved properties (on spatial scales of $\sim 1-2$ kpc out to at least $3$ effective radii) of the stars, dust, and gas in ten nearby spiral galaxies. The properties of the stellar population and dust are derived by fitting the spatially resolved spectral energy distribution (SED) with more than 20 photometric bands ranging from far-ultraviolet to far-infrared. Our newly developed software piXedfit performs point spread function matching of images, pixel binning, and models the stellar light, dust attenuation, dust emission, and emission from a dusty torus heated by an active galactic nucleus simultaneously through the energy balance approach. With this self-consistent analysis, we present the spatially resolved version of the IRX--$β$ relation, finding that it is consistent with the relationship from the integrated photometry. We show that the old stellar populations contribute to the dust heating, which causes an overestimation of star formation rate (SFR) derived from the total ultraviolet and infrared luminosities on kpc scales. With archival high-resolution maps of atomic and molecular gas, we study the radial variation of the properties of the stellar populations (including stellar mass, age, metallicity, and SFR), dust (including dust mass, dust temperature, and abundance of polycyclic aromatic hydrocarbon), gas, as well as dust-to-stellar mass and dust-to-gas mass ratios. We observe a depletion of molecular gas mass fraction in the central region of the majority of the galaxies, suggesting that the lack of available fuel is an important factor in suppressing the specific SFR at the center.

preprint2021arXiv

Evidence from Disrupted Halo Dwarfs that $r$-process Enrichment via Neutron Star Mergers is Delayed by $\gtrsim500$ Myrs

The astrophysical origins of $r$-process elements remain elusive. Neutron star mergers (NSMs) and special classes of core-collapse supernovae (rCCSNe) are leading candidates. Due to these channels&#39; distinct characteristic timescales (rCCSNe: prompt, NSMs: delayed), measuring $r$-process enrichment in galaxies of similar mass, but differing star-formation durations might prove informative. Two recently discovered disrupted dwarfs in the Milky Way&#39;s stellar halo, Kraken and \textit{Gaia}-Sausage Enceladus (GSE), afford precisely this opportunity: both have $M_{\star}\approx10^{8}M_{\rm{\odot}}$, but differing star-formation durations of ${\approx}2$ Gyrs and ${\approx}3.6$ Gyrs. Here we present $R\approx50,000$ Magellan/MIKE spectroscopy for 31 stars from these systems, detecting the $r$-process element Eu in all stars. Stars from both systems have similar [Mg/H]$\approx-1$, but Kraken has a median [Eu/Mg]$\approx-0.1$ while GSE has an elevated [Eu/Mg]$\approx0.2$. With simple models we argue NSM enrichment must be delayed by $500-1000$ Myrs to produce this difference. rCCSNe must also contribute, especially at early epochs, otherwise stars formed during the delay period would be Eu-free. In this picture, rCCSNe account for $\approx50\%$ of the Eu in Kraken, $\approx25\%$ in GSE, and $\approx15\%$ in dwarfs with extended star-formation durations like Sagittarius. The inferred delay time for NSM enrichment is $10-100\times$ longer than merger delay times from stellar population synthesis -- this is not necessarily surprising because the enrichment delay includes time taken for NSM ejecta to be incorporated into subsequent generations of stars. For example, this may be due to natal kicks that result in $r$-enriched material deposited far from star-forming gas, which then takes $\approx10^{8}-10^{9}$ years to cool in these galaxies.

preprint2021arXiv

On the formation of massive quiescent galaxies with diverse morphologies in the TNG50 simulation

Observations have shown that the star-formation activity and the morphology of galaxies are closely related, but the underlying physical connection is not well understood. Using the TNG50 simulation, we explore the quenching and the morphological evolution of the 102 massive quiescent galaxies in the mass range of $10.5<\log(M_{\rm stellar}/M_{\odot})<11.5$ selected at $z=0$. The morphology of galaxies is quantified based on their kinematics, and we measure the quenching timescale of individual galaxies directly from star formation history. We show that galaxies tend to be quenched more rapidly if they: (i) are satellites in massive halos, (ii) have lower star-forming gas fractions, or (iii) inject a larger amount of black hole kinetic feedback energy. By following the global evolutionary pathways, we conclude that quiescent discs are mainly disc galaxies that are recently and slowly quenched. Approximately half of the quiescent ellipticals at $z=0$ are rapidly quenched at higher redshifts while still disc-like. While being quiescent, they gradually become more elliptical mostly by disc heating, yet these ellipticals still retain some degree of rotation. The other half of quiescent ellipticals with the most random motion-dominated kinematics build up large spheroidal components before quenching primarily by mergers, or in some cases, misaligned gas accretion. However, the mergers that contribute to morphological transformation do not immediately quench galaxies in many cases. In summary, we find that quenching and morphological transformation are decoupled. We conclude that the TNG black hole feedback -- in combination with the stochastic merger history of galaxies -- leads to a large diversity of quenching timescales and a rich morphological landscape.

preprint2021arXiv

The Diverse Molecular Gas Content of Massive Galaxies Undergoing Quenching at z~1

We present a detailed study of the molecular gas content and stellar population properties of three massive galaxies at 1 < z < 1.3 that are in different stages of quenching. The galaxies were selected to have a quiescent optical/near-infrared spectral energy distribution and a relatively bright emission at 24 micron, and show remarkably diverse properties. CO emission from each of the three galaxies is detected in deep NOEMA observations, allowing us to derive molecular gas fractions Mgas/Mstar of 13-23%. We also reconstruct the star formation histories by fitting models to the observed photometry and optical spectroscopy, finding evidence for recent rejuvenation in one object, slow quenching in another, and rapid quenching in the third system. To better constrain the quenching mechanism we explore the depletion times for our sample and other similar samples at z~0.7 from the literature. We find that the depletion times are highly dependent on the method adopted to measure the star formation rate: using the UV+IR luminosity we obtain depletion times about 6 times shorter than those derived using dust-corrected [OII] emission. When adopting the star formation rates from spectral fitting, which are arguably more robust, we find that recently quenched galaxies and star-forming galaxies have similar depletion times, while older quiescent systems have longer depletion times. These results offer new, important constraints for physical models of galaxy quenching.

preprint2021arXiv

The galaxy-halo connection of emission-line galaxies in IllustrisTNG

We employ the hydrodynamical simulation IllustrisTNG-300-1 to explore the halo occupation distribution (HOD) and environmental dependence of luminous star-forming emission-line galaxies (ELGs) at $z \sim 1$. Such galaxies are key targets for current and upcoming cosmological surveys. We select model galaxies through cuts in colour-colour space allowing for a direct comparison with the Extended Baryon Oscillation Spectroscopic Survey and the Dark Energy Spectroscopic Instrument (DESI) surveys and then compare them with galaxies selected based on specific star-formation rate (sSFR) and stellar mass. We demonstrate that the ELG populations are twice more likely to reside in lower-density regions (sheets) compared with the mass-selected populations and twice less likely to occupy the densest regions of the cosmic web (knots). We also show that the colour-selected and sSFR-selected ELGs exhibit very similar occupation and clustering statistics, finding that the agreement is best for lower redshifts. In contrast with the mass-selected sample, the occupation of haloes by a central ELG peaks at $\sim$20\%. We furthermore explore the dependence of the HOD and the auto-correlation on environment, noticing that at fixed halo mass, galaxies in high-density regions cluster about 10 times more strongly than low-density ones. This result suggests that we should model carefully the galaxy-halo relation and implement assembly bias effects into our models (estimated at $\sim$4\% of the clustering of the DESI colour-selected sample at $z = 0.8$). Finally, we apply a simple mock recipe to recover the clustering on large scales ($r \gtrsim 1 \ {\rm Mpc}/h$) to within 1\% by augmenting the HOD model with an environment dependence, demonstrating the power of adopting flexible population models.

preprint2020arXiv

A panchromatic spatially resolved analysis of nearby galaxies -- II. The main sequence - gas relation at sub-kpc scale in grand-design spirals

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 ($Σ_{\star}$ and $Σ_{\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 ($σ=0.23$ dex) is the combination of two relations: the Kennicutt-Schmidt law ($σ=0.19$ dex) and the molecular gas MS (MGMS, $σ=0.22$ dex); $Σ_{\star}$, $Σ_{\rm{SFR}}$ and the surface density of the molecular gas, $Σ_{\rm{H_2}}$, define a 3D relation as proposed by \citet{2019ApJ...884L..33L}. We find that $Σ_{\rm{H_2}}$ steadily increases along the MS relation, varies little towards higher $Σ_{\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 ($Σ_{\rm{HI}}$) is constant along the MS, and increases in its upper envelop. $Σ_{\rm{SFR}}$ can be expressed as a function of $Σ_{\star}$ and $Σ_{\rm{HI}}$, following the Equation: $\logΣ_{\rm{SFR}}$ = 0.97$\logΣ_{\star}$ + 1.99$\logΣ_{\rm{HI}}$ - 11.11. Finally, we show that f$_{\rm{gas}}$ increases significantly towards the starburst region in the $\logΣ_{\star}$ - $\logΣ_{\rm{SFR}}$ plane, accompanied by a slight increase in SFE.

preprint2020arXiv

Dust, gas, and metal content in star-forming galaxies at $z\sim3.3$ revealed with ALMA and Near-IR spectroscopy

We conducted sub-millimeter observations with the Atacama Large Millimeter/sub-millimeter Array (ALMA) of star-forming galaxies at $z\sim3.3$, whose gas-phase metallicities have been previously measured. We investigate the dust and gas contents of the galaxies at $z\sim3.3$ and study how galaxies are interacting with their circumgalactic/intergalactic medium at this epoch by probing their gas mass fractions and gas-phase metallicities. Single-band dust continuum emission tracing dust mass and the relation between the gas-phase metallicity and gas-to-dust mass ratio are used to estimate the gas masses. The estimated gas mass fractions and depletion timescales are $f_{\rm gas}=$ 0.20-0.75 and $t_{\rm dep}=$ 0.09-1.55 Gyr, respectively. Although the galaxies appear to tightly distribute around the star-forming main sequence at $z\sim3.3$, both quantities show a wider spread at a fixed stellar mass than expected from the scaling relation, suggesting a large diversity of fundamental gas properties among star-forming galaxies apparently on the main sequence. Comparing gas mass fraction and gas-phase metallicity between the star-forming galaxies at $z\sim3.3$ and at lower redshifts, star-forming galaxies at $z\sim3.3$ appear to be more metal-poor than local galaxies with similar gas mass fractions. Using the gas regulator model to interpret this offset, we find that it can be explained by a higher mass-loading factor, suggesting that the mass-loading factor in outflows increases at earlier cosmic times.

preprint2020arXiv

Galaxy Sizes Since $z=2$ from the Perspective of Stellar Mass Distribution within Galaxies

How stellar mass assembles within galaxies is still an open question. We present measurements of the stellar mass distribution on kpc-scale for $\sim5500$ galaxies with stellar masses above $\log(M_{\ast}/M_{\odot})\geqslant9.8$ up to the redshift $2.0$. We create stellar mass maps from Hubble Space Telescope observations by means of the pixel-by-pixel SED fitting method. These maps are used to derive radii encompassing $20\%$, $50\%$, and $80\%$ ($r_{20}$, $r_{50}$ and $r_{80}$) of the total stellar mass from the best-fit Sérsic models. The reliability and limitations of the structural parameter measurements are checked extensively using a large sample ($\sim3000$) of simulated galaxies. The size-mass relations and redshift evolution of $r_{20}$, $r_{50}$ and $r_{80}$ are explored for star-forming and quiescent galaxies. At fixed mass, the star-forming galaxies do not show significant changes in their $r_{20}$, $r_{50}$ and $r_{80}$ sizes, indicating self-similar growth. Only above the pivot stellar mass of $\log(M_{\ast}/M_{\odot})\simeq10.5$, $r_{80}$ evolves as $r_{80}\propto(1+z)^{-0.85\pm0.20}$, indicating that mass builds up in the outskirts of these systems (inside-out growth). The Sérsic values also increase for the massive star-forming galaxies towards late cosmic time. Massive quiescent galaxies show stronger size evolution at all radii, in particular the $r_{20}$ sizes. For these massive galaxies, Sérsic values remain almost constant since at least $z\sim1.3$, indicating that the strong size evolution is related to the changes in the outer parts of these galaxies. We make all the structural parameters publicly available.

preprint2020arXiv

High redshift JWST predictions from IllustrisTNG: II. Galaxy line and continuum spectral indices and dust attenuation curves

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}_α$ and ${\rm H}_β$ + $[\rm O \,III]$ luminosity functions up to $z=8$. The predicted ${\rm H}_β$ + $[\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}_α$ 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$Å(D4000) and the UV continuum slope $β$. 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 $β$ 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.

preprint2020arXiv

Photometric properties of reionization-epoch galaxies in the Simba simulations

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&#39;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 ($β$) 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 $β$. We find that $β$ 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&#39;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.

preprint2020arXiv

Quenching as a Contest between Galaxy Halos and their Central Black Holes

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&#39;&#39; 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.

preprint2020arXiv

Rapid Reionization by the Oligarchs: The Case for Massive, UV-Bright, Star-Forming Galaxies with High Escape Fractions

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$α$ 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 ($Σ$), we propose a model relating $f_{esc}$ to $Σ$ and find $f_{esc}\proptoΣ^{0.4\pm0.1}$. Since $Σ$ 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&#39;) 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 ($α$>-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)

preprint2020arXiv

Simulating JWST/NIRCam Color Selection of High-Redshift Galaxies

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$α$ 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.

preprint2020arXiv

Stochastic modelling of star-formation histories II: star-formation variability from molecular clouds and gas inflow

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 (&#34;regulator model&#34;). 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 $β\approx2$ at high frequencies with a break near the average cloud lifetime. We discuss star-formation &#34;burstiness&#34; 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.

preprint2020arXiv

The Diversity and Variability of Star Formation Histories in Models of Galaxy Evolution

Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIRE-2, g14, and Marvel/Justice League), semi-analytic models (Santa Cruz SAM) and empirical models (UniverseMachine), and quantify the variability of these SFHs on different timescales using the power spectral density (PSD) formalism. We find that the PSDs are well described by broken power-laws, and variability on long timescales ($\gtrsim1$ Gyr) accounts for most of the power in galaxy SFHs. Most hydrodynamical models show increased variability on shorter timescales ($\lesssim300$ Myr) with decreasing stellar mass. Quenching can induce $\sim0.4-1$ dex of additional power on timescales $>1$ Gyr. The dark matter accretion histories of galaxies have remarkably self-similar PSDs and are coherent with the in-situ star formation on timescales $>3$ Gyr. There is considerable diversity among the different models in their (i) power due to SFR variability at a given timescale, (ii) amount of correlation with adjacent timescales (PSD slope), (iii) evolution of median PSDs with stellar mass, and (iv) presence and locations of breaks in the PSDs. The PSD framework is a useful space to study the SFHs of galaxies since model predictions vary widely. Observational constraints in this space will help constrain the relative strengths of the physical processes responsible for this variability.

preprint2019arXiv

High redshift JWST predictions from IllustrisTNG: Dust modelling and galaxy luminosity functions

The James Webb Space Telescop (JWST) promises to revolutionise our understanding of the early Universe, and contrasting its upcoming observations with predictions of the $Λ$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$α$ versus star formation rate relation, and the H$α$ 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}$.

preprint2019arXiv

The MUSE Atlas of Disks (MAD): Ionized gas kinematic maps and an application to Diffuse Ionized Gas

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, σ$) 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.