Comparison of Observed Galaxy Properties with Semianalytic Model Predictions using Machine Learning

With current and upcoming experiments such as WFIRST, Euclid and LSST, we can observe up to billions of galaxies. While such surveys cannot obtain spectra for all observed galaxies, they produce galaxy magnitudes in color filters. This data set behaves like a high-dimensional nonlinear surface, an excellent target for machine learning. In this work, we use a lightcone of semianalytic galaxies tuned to match CANDELS observations from Lu et al. (2014) to train a set of neural networks on a set of galaxy physical properties. We add realistic photometric noise and use trained neural networks to predict stellar masses and average star formation rates on real CANDELS galaxies, comparing our predictions to SED fitting results. On semianalytic galaxies, we are nearly competitive with template-fitting methods, with biases of $0.01$ dex for stellar mass, $0.09$ dex for star formation rate, and $0.04$ dex for metallicity. For the observed CANDELS data, our results are consistent with template fits on the same data at $0.15$ dex bias in $M_{\rm star}$ and $0.61$ dex bias in star formation rate. Some of the bias is driven by SED-fitting limitations, rather than limitations on the training set, and some is intrinsic to the neural network method. Further errors are likely caused by differences in noise properties between the semianalytic catalogs and data. Our results show that galaxy physical properties can in principle be measured with neural networks at a competitive degree of accuracy and precision to template-fitting methods.Comment: 19 pages, 10 figures, 6 tables. Accepted for publication in Ap

Dependence of the IRX-β Dust Attenuation Relation on Metallicity and Environment

We use a sample of star-forming field and protocluster galaxies at z = 2.0–2.5 with Keck/MOSFIRE K-band spectra, a wealth of rest-frame ultraviolet (UV) photometry, and Spitzer/MIPS and Herschel/PACS observations, to dissect the relation between the ratio of infrared (IR) to UV luminosity (IRX) versus UV slope (β) as a function of gas-phase metallicity (12+log(O/H) ~ 8.2–8.7). We find no significant dependence of the IRX-β trend on environment. However, we find that at a given β, IRX is highly correlated with metallicity, and less correlated with mass, age, and specific star formation rate (sSFR). We conclude that, of the physical properties tested here, metallicity is the primary physical cause of the IRX-β scatter, and the IRX correlation with mass is presumably due to the mass dependence on metallicity. Our results indicate that the UV attenuation curve steepens with decreasing metallicity, and spans the full range of slope possibilities from a shallow Calzetti-type curve for galaxies with the highest metallicity in our sample (12+log(O/H) ~ 8.6) to a steep Small Magellanic Cloud (SMC)-like curve for those with 12+log(O/H) ~ 8.3. Using a Calzetti (SMC) curve for the low (high) metallicity galaxies can lead to up to a factor of 3 overestimation (underestimation) of the UV attenuation and obscured star formation rate. We speculate that this change is due to different properties of dust grains present in the interstellar medium of low- and high-metallicity galaxies

LATIS: The Stellar Mass-Metallicity Relation of Star-forming Galaxies at z∼2.5z\sim 2.5

We present the stellar mass - stellar metallicity relation for 3491 star-forming galaxies at $2 \lesssim z \lesssim 3$ using rest-frame far-ultraviolet (FUV) spectra from the Ly$\alpha$ Tomography IMACS Survey (LATIS). We fit stellar population synthesis models from the Binary Population And Spectral Synthesis code (BPASS v$2.2.1$) to medium resolution (R $\sim 1000$) and high signal-to-noise ($>30$ per 100 km/s over a wavelength range of 1221 - 1800 \r{A}) composite spectra of galaxies in bins of stellar mass to determine their stellar metallicity, primarily tracing $\rm Fe/H$. We find a strong correlation between stellar mass and stellar metallicity, with stellar metallicity monotonically increasing with stellar mass at low masses and flattening at high masses ($M_* \gtrsim 10^{10.3} M_\odot$). Additionally, we compare our stellar metallicity measurements with the gas-phase oxygen abundance of galaxies at similar redshift and estimate the average $\rm [\alpha/Fe] \sim 0.6$. Such high $\alpha$-enhancement indicates that high-redshift galaxies have not yet undergone significant iron enrichment through Type Ia supernovae. Moreover, we utilize an analytic chemical evolution model to constrain the mass loading parameter of galactic winds as a function of stellar mass. We find that as the stellar mass increases, the mass loading parameter decreases. The parameter then flattens or reaches a turning point at around $M_* \sim 10^{10.5} M_\odot$. Our findings may signal the onset of black hole-driven outflows at $z \sim 2.5$ for galaxies with $M_* \gtrsim 10^{10.5} M_\odot$.Comment: 21 pages, 17 figures, 3 tables, accepted for publication in Ap

Dependence of the IRX-β\beta dust attenuation relation on metallicity and environment

We use a sample of star-forming field and protocluster galaxies at z=2.0-2.5 with Keck/MOSFIRE K-band spectra, a wealth of rest-frame UV photometry, and Spitzer/MIPS and Herschel/PACS observations, to dissect the relation between the ratio of IR to UV luminosity (IRX) versus UV slope ($\beta$) as a function of gas-phase metallicity (12+log(O/H)~8.2-8.7). We find no significant dependence of the IRX-$\beta$ trend on environment. However, we find that at a given $\beta$, IRX is highly correlated with metallicity, and less correlated with mass, age, and sSFR. We conclude that, of the physical properties tested here, metallicity is the primary physical cause of the IRX-$\beta$ scatter, and the IRX correlation with mass is presumably due to the mass dependence on metallicity. Our results indicate that the UV attenuation curve steepens with decreasing metallicity, and spans the full range of slope possibilities from a shallow Calzetti-type curve for galaxies with the highest metallicity in our sample (12+log(O/H)~8.6) to a steep SMC-like curve for those with 12+log(O/H)~8.3. Using a Calzetti (SMC) curve for the low (high) metallicity galaxies can lead to up to a factor of 3 overestimation (underestimation) of the UV attenuation and obscured SFR. We speculate that this change is due to different properties of dust grains present in the ISM of low- and high-metallicity galaxies.Comment: Accepted for publication in ApJ

Bridging between the integrated and resolved main sequence of star formation

The position of galaxies on the stellar mass, star formation rate (SFR) plane with respect to the star-forming main sequence at each redshift is a convenient way to infer where the galaxy is in its evolution compared to the rest of the population. We use Hubble Space Telescope high-resolution images in the GOODS-S field from the the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) and fit multiwavelength lights in resolution elements of galaxies with stellar population synthesis models. We then construct resolved kpc-scale stellar mass, SFR surface density curves for galaxies at z ~ 1. Fitting these resolved main sequence curves with Schechter functions, we parameterize and explain the multiwavelength structure of galaxies with three variables: φ*, α, and M*. For quenched galaxies below the main sequence, we find an average high-mass slope (α) of the resolved main sequence curves to be ~−0.4. The scatter of this slope is higher among the lower mass star-forming galaxies and those above the main sequence compared to quenched galaxies, due to lack of an evolved bulge. Our findings agree well with an inside-out quenching of star formation. We find that the knee of the Schechter fits (M*) for galaxies below the main sequence occurs at lower stellar mass surface densities compared to star-forming galaxies, which hints at how far quenching has proceeded outward

Bridging between the integrated and resolved main sequence of star formation

The position of galaxies on the stellar mass, star formation rate (SFR) plane with respect to the star-forming main sequence at each redshift is a convenient way to infer where the galaxy is in its evolution compared to the rest of the population. We use Hubble Space Telescope high-resolution images in the GOODS-S field from the the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) and fit multiwavelength lights in resolution elements of galaxies with stellar population synthesis models. We then construct resolved kpc-scale stellar mass, SFR surface density curves for galaxies at z ~ 1. Fitting these resolved main sequence curves with Schechter functions, we parameterize and explain the multiwavelength structure of galaxies with three variables: φ*, α, and M*. For quenched galaxies below the main sequence, we find an average high-mass slope (α) of the resolved main sequence curves to be ~−0.4. The scatter of this slope is higher among the lower mass star-forming galaxies and those above the main sequence compared to quenched galaxies, due to lack of an evolved bulge. Our findings agree well with an inside-out quenching of star formation. We find that the knee of the Schechter fits (M*) for galaxies below the main sequence occurs at lower stellar mass surface densities compared to star-forming galaxies, which hints at how far quenching has proceeded outward

Exploring the Correlation between Hα\rm{H}\alpha-to-UV Ratio and Burstiness for Typical Star-forming Galaxies at z∼2z\sim2

The $\rm{H}\alpha$-to-UV luminosity ratio ($L(\rm H\alpha)/L(\rm UV)$) is often used to probe SFHs of star-forming galaxies and it is important to validate it against other proxies for burstiness. To address this issue, we present a statistical analysis of the resolved distribution of $\Sigma_{\rm{SFR}}$ as well as stellar age and their correlations with the globally measured $L(\rm H\alpha)/L(\rm UV)$ for a sample of 310 star-forming galaxies in two redshift bins of $1.37 < z < 1.70$ and $2.09 < z < 2.61$ observed by the MOSDEF survey. We use the multi-waveband CANDELS/3D-HST imaging of MOSDEF galaxies to construct $\Sigma_{\rm{SFR}}$ and stellar age maps. We analyze the composite rest-frame far-UV spectra of a subsample of MOSDEF targets obtained by the Keck/LRIS, which includes 124 star-forming galaxies (MOSDEF-LRIS) at redshifts $1.4 < z < 2.6$, to examine the average stellar population properties, and the strength of age-sensitive FUV spectral features in bins of $L(\rm H\alpha)/L(\rm UV)$. Our results show no significant evidence that individual galaxies with higher $L(\rm H\alpha)/L(\rm UV)$ are undergoing a burst of star formation based on the resolved distribution of $\Sigma_{\rm{SFR}}$ of individual star-forming galaxies. We segregate the sample into subsets with low and high $L(\rm H\alpha)/L(\rm UV)$. The high-$L(\rm H\alpha)/L(\rm UV)$ subset exhibits, on average, an age of $\log[\rm{Age/yr}]$ = 8.0, compared to $\log[\rm{Age/yr}]$ = 8.4 for the low-$L(\rm H\alpha)/L(\rm UV)$ galaxies, though the difference in age is significant at only the $2\sigma$ level. Furthermore, we find no variation in the strengths of Siiv$\lambda\lambda1393, 1402$ and Civ$\lambda\lambda1548, 1550$ P-Cygni features from massive stars between the two subsamples.Comment: 16 pages, 10 figures, published by the Monthly Notices of the Royal Astronomical Societ

The clustering of typical Lyαα emitters from z∼2.5−6z \sim 2.5 - 6: host halo masses depend on Lyαα and UV luminosities

We investigate the clustering and halo properties of $\sim 5000$ Ly$\alpha$-selected emission line galaxies (LAEs) from the Slicing COSMOS 4K (SC4K) and from archival NB497 imaging of SA22 split in 15 discrete redshift slices between $z \sim 2.5 - 6$. We measure clustering lengths of $r_0 \sim 3 - 6\ h^{-1}$ Mpc and typical halo masses of $\sim 10^{11}$ M$_\odot$ for our narrowband-selected LAEs with typical $L_{\rm{Ly}\alpha} \sim 10^{42 - 43}$ erg s$^{-1}$. The intermediate band-selected LAEs are observed to have $r_0 \sim 3.5 - 15\ h^{-1}$ Mpc with typical halo masses of $\sim 10^{11 - 12}$ M$_\odot$ and typical $L_{\rm{Ly}\alpha} \sim 10^{43 - 43.6}$ erg s$^{-1}$. We find a strong, redshift-independent correlation between halo mass and Ly$\alpha$ luminosity normalized by the characteristic Ly$\alpha$ luminosity, $L^\star(z)$. The faintest LAEs ($L \sim 0.1\ L^\star(z)$) typically identified by deep narrowband surveys are found in $10^{10}$ M$_\odot$ halos and the brightest LAEs ($L \sim 7\ L^\star(z)$) are found in $\sim 5 \times 10^{12}$ M$_\odot$ halos. A dependency on the rest-frame 1500 \AA~UV luminosity, M_\rm{UV}, is also observed where the halo masses increase from $10^{11}$ to $10^{13}$ M$_\odot$ for M_\rm{UV} \sim -19 to $-23.5$ mag. Halo mass is also observed to increase from $10^{9.8}$ to $10^{12.3}$ M$_\odot$ for dust-corrected UV star formation rates from $\sim 0.6$ to $10$ M$_\odot$ yr$^{-1}$ and continues to increase up to $10^{13.5}$ M$_\odot$ in halo mass, where the majority of those sources are AGN. All the trends we observe are found to be redshift-independent. Our results reveal that LAEs are the likely progenitors of a wide range of galaxies depending on their luminosity, from dwarf-like, to Milky Way-type, to bright cluster galaxies. LAEs therefore provide unique insight into the early formation and evolution of the galaxies we observe in the local Universe