Search

We combined data from the Sloan Digital Sky Survey (SDSS) and the Arecibo Legacy Fast ALFA Survey (ALFALFA) to establish the HI mass vs. stellar mass and halo mass scaling relations using an abundance matching method that is free of the Malmquist bias. To enable abundance matching, a cross-match between the SDSS DR7 galaxy group sample and the ALFALFA HI sources provides a catalog of 16,520 HI-galaxy pairs within 14,270 galaxy groups (halos). By applying the observational completeness reductions for both optical and HI observations, we used the remaining 8,180 ALFALFA matched sources to construct the model constraints. Taking into account the dependence of HI mass on both the galaxy and group properties, we establish two sets of scaling relations: one with a combination of stellar mass, $({g-r})$ color and halo mass, and the other with stellar mass, specific star-formation rate ($\rm sSFR$), and halo mass. We demonstrate that our models can reproduce the HI mass component as both a stellar and halo mass. Additional tests showed that the conditional HI mass distributions as a function of the cosmic web type and the satellite fractions were well recovered., Comment: arXiv admin note: text overlap with arXiv:2008.09804

In this paper, we set out to construct a set of reference mock galaxy redshift surveys (MGRSs) for the future Chinese Space-station Survey Telescope (CSST) observation, where subsequent survey selection effects can be added and evaluated. This set of MGRSs is generated using the dark matter subhalos extracted from a high-resolution Jiutian $N$-body simulation of the standard $\Lambda$CDM cosmogony with $\Omega_m=0.3111$, $\Omega_{\Lambda}=0.6889$, and $\sigma_8=0.8102$. The simulation has a boxsize of $1~h^{-1} {\rm Gpc}$, and consists of $6144^3$ particles with mass resolution $3.723 \times 10^{8} h^{-1} M_\odot $. In order to take into account the effect of redshift evolution, we first use all 128 snapshots in the Jiutian simulation to generate a light-cone halo/subhalo catalog. Next, galaxy luminosities are assigned to the main and subhalo populations using the subhalo abundance matching (SHAM) method with the DESI $z$-band luminosity functions at different redshifts. Multi-band photometries, as well as images, are then assigned to each mock galaxy using a 3-dimensional parameter space nearest neighbor sampling of the DESI LS observational galaxies and groups. Finally, the CSST and DESI LS survey geometry and magnitude limit cuts are applied to generate the required MGRSs. As we have checked, this set of MGRSs can generally reproduce the observed galaxy luminosity/mass functions within 0.1 dex for galaxies with $L > 10^8 L_\odot$ (or $M_* > 10^{8.5} M_\odot$) and within 1-$\sigma$ level for galaxies with $L < 10^8L_\odot$ (or $M_* < 10^{8.5} M_\odot$). Together with the CSST slitless spectra and redshifts for our DESI LS seed galaxies that are under construction, we will set out to test various slitless observational selection effects in subsequent probes., Comment: 13 pages, 9 figures, accepted for publication in MNRAS

Dust attenuation in star-forming galaxies (SFGs), as parameterized by the infrared excess (IRX $\equiv L_{\rm IR}/L_{\rm UV}$), is found to be tightly correlated with star formation rate (SFR), metallicity and galaxy size, following a universal IRX relation up to $z=3$. This scaling relation can provide a fundamental constraint for theoretical models to reconcile galaxy star formation, chemical enrichment, and structural evolution across cosmic time. We attempt to reproduce the universal IRX relation over $0.1\leq z\leq 2.5$ using the EAGLE hydrodynamical simulations and examine sensitive parameters in determining galaxy dust attenuation. Our findings show that while the predicted universal IRX relation from EAGLE approximately aligns with observations at $z\leq 0.5$, noticeable disparities arise at different stellar masses and higher redshifts. Specifically, we investigate how modifying various galaxy parameters can affect the predicted universal IRX relation in comparison to the observed data. We demonstrate that the simulated gas-phase metallicity is the critical quantity for the shape of the predicted universal IRX relation. We find that the influence of the infrared luminosity and infrared excess is less important while galaxy size has virtually no significant effect. Overall, the EAGLE simulations are not able to replicate some of the observed characteristics between IRX and galaxy parameters of SFGs, emphasizing the need for further investigation and testing for our current state-of-the-art theoretical models., Comment: 19 pages, 15 figures, accepted for publication in MNRAS

In this investigation, we leverage the combination of Dark Energy Spectroscopic Instrument Legacy imaging Surveys Data Release 9 (DESI LS DR9), Survey Validation 3 (SV3), and Year 1 (Y1) data sets to estimate the conditional luminosity and stellar mass functions (CLFs & CSMFs) of galaxies across various halo mass bins and redshift ranges. To support our analysis, we utilize a realistic DESI Mock Galaxy Redshift Survey (MGRS) generated from a high-resolution Jiutian simulation. An extended halo-based group finder is applied to both MGRS catalogs and DESI observation. By comparing the r and z-band luminosity functions (LFs) and stellar mass functions (SMFs) derived using both photometric and spectroscopic data, we quantified the impact of photometric redshift (photo-z) errors on the galaxy LFs and SMFs, especially in the low redshift bin at low luminosity/mass end. By conducting prior evaluations of the group finder using MGRS, we successfully obtain a set of CLF and CSMF measurements from observational data. We find that at low redshift the faint end slopes of CLFs and CSMFs below $10^{9}h^{-2}L_{\odot}$ (or $h^{-2}M_{\odot}$) evince a compelling concordance with the subhalo mass functions. After correcting the cosmic variance effect of our local Universe following arXiv:1809.00523, the faint end slopes of the LFs/SMFs turn out to be also in good agreement with the slope of the halo mass function., Comment: 28 pages, 13 figures, Accepted for publication in ApJ

Star-forming galaxies (SFGs) adhere to a surprisingly tight scaling relation of dust attenuation parameterized by the infrared excess (IRX=$L_{\rm IR}/L_{\rm UV}$), being jointly determined by the star formation rate (SFR), galaxy size ($R_{\rm e}$), metallicity ($Z$/Z$_\odot$) and axial ratio ($b/a$). We examine how these galaxy parameters determine the effective dust attenuation and give rise to the universal IRX relation, utilizing a simple two-component star-dust geometry model in which dust in the dense and diffuse interstellar medium (ISM) follows exponential mass density profiles, connected with but not necessarily identical to the stellar mass profiles. Meanwhile, empirical relations are adopted to link galaxy properties, including the gas--star formation relation, the dust-to-stellar size relation, as well as the dust-to-gas ratio versus metallicity relation. By fitting a large sample of local SFGs with the model, we obtain the best-fitting model parameters as a function of metallicity, showing that the two-component geometry model is able to successfully reproduce the dependence of IRX on SFR, $R_{\rm e}$, $b/a$ at given $Z$/Z$_\odot$, as well as the dependence of power-law indices on metallicity. Moreover, we also retrieve constraints on the model geometry parameters, including the optical depth of birth clouds (BCs), BC-to-total dust mass fraction, BC covering factor of UV-emitting stars, and star-to-total dust disc radius ratio, which all evolve with galaxy metallicity. Finally, a consistent picture of how the star-dust geometry in SFGs evolves with galaxy metallicity is discussed., Comment: 20 pages, 10 figures, published in MNRAS (2024, Volume 528, Issue 1, pp.658-675); A PHTHON package IRX_TAU_TOT is available at https://github.com/LvZF/irx_tau_tot/ to calculate the total dust optical depth of a galaxy with given metallicity and best-fitting geometry parameters

We present a tentative constraint on cosmological parameters $\Omega_m$ and $\sigma_8$ from a joint analysis of galaxy clustering and galaxy-galaxy lensing from DESI Legacy Imaging Surveys Data Release 9 (DR9), covering approximately 10000 square degrees and spanning the redshift range of 0.1 to 0.9. To study the dependence of cosmological parameters on lens redshift, we divide lens galaxies into seven approximately volume-limited samples, each with an equal width in photometric redshift. To retrieve the intrinsic projected correlation function $w_{\rm p}(r_{\rm p})$ from the lens samples, we employ a novel method to account for redshift uncertainties. Additionally, we measured the galaxy-galaxy lensing signal $\Delta\Sigma(r_{\rm p})$ for each lens sample, using source galaxies selected from the shear catalog by applying our \texttt{Fourier\_Quad} pipeline to DR9 images. We model these observables within the flat $\Lambda$CDM framework, employing the minimal bias model. To ensure the reliability of the minimal bias model, we apply conservative scale cuts: $r_{\rm p} > 8$ and $12 ~h^{-1}{\rm Mpc}$, for $w_{\rm p}(r_{\rm p})$ and $\Delta\Sigma(r_{\rm p})$, respectively. Our findings suggest a mild tendency that $S_8 \equiv \sigma_8 \sqrt{\Omega_m/0.3} $ increases with lens redshift, although this trend is only marginally significant. When we combine low redshift samples, the value of $S_8$ is determined to be $0.84 \pm 0.02$, consistent with the Planck results but significantly higher than the 3$\times$ 2pt analysis by 2-5$\sigma$. Despite the fact that further refinements in measurements and modeling could improve the accuracy of our results, the consistency with standard values demonstrates the potential of our method for more precise and accurate cosmology in the future., Comment: slightly different with the published version

The Kilo Degree Survey (KiDS) is currently the only sky survey providing optical ($ugri$) plus near-infrared (NIR, $ZYHJK_S$) seeing matched photometry over an area larger than 1000 $\rm deg^2$. This is obtained by incorporating the NIR data from the VISTA Kilo Degree Infrared Galaxy (VIKING) survey, covering the same KiDS footprint. As such, the KiDS multi-wavelength photometry represents a unique dataset to test the ability of stellar population models to return robust photometric stellar mass ($M_*$) and star-formation rate (SFR) estimates. Here we use a spectroscopic sample of galaxies for which we possess $u g r i Z Y J H K_s$ ``gaussianized'' magnitudes from KiDS data release 4. We fit the spectral energy distribution from the 9-band photometry using: 1) three different popular libraries of stellar {population} templates, 2) single burst, simple and delayed exponential star-formation history models, and 3) a wide range of priors on age and metallicity. As template fitting codes we use two popular softwares: LePhare and CIGALE. We investigate the variance of the stellar masses and the star-formation rates from the different combinations of templates, star formation recipes and codes to assess the stability of these estimates and define some ``robust'' median quantities to be included in the upcoming KiDS data releases. As a science validation test, we derive the mass function, the star formation rate function, and the SFR-$M_*$ relation for a low-redshift ($z<0.5$) sample of galaxies, that result in excellent agreement with previous literature data. The final catalog, containing $\sim290\,000$ galaxies with redshift $0.01

We present a physical model of the Mass Accretion Histories (MAH) of haloes in concordance with the {\it observed} cosmic star formation rate density (CSFRD). We model the MAHs of dark matter haloes using a Gamma ($\Gamma$) functional form: $M_h(T) = \frac{M_0}{f_{0}} \, \times \frac{\gamma(\alpha_h, ~\beta_h \times (T-Th))}{\Gamma(\alpha_h)}$, where $M_0$ is the halo mass at present time, $T$ is time, $\alpha_h$ and $\beta_h$ are parameters we explore, $f_{0}$ is the percentage of the mass of the halo at z = 0 with respect to the final mass of the halo achieved at $T = \infty$. We use the MAHs of haloes obtained from cosmological simulations and analytical models to constrain our model. $f_{0}$ can be described by a power-law ($f_{0} = 1- c \times M_{0}^{d}$). Haloes with small masses have already on average attained most of their final masses. The average $$ of haloes in the Universe is $ > 0.95$ pointing to the direction that the cosmic MAH/CSFRD is saturated at our era. The average $<\beta_{h}>$ parameter (the depletion rate of the available dark matter for halo growth) is related to the dynamical timescales of haloes. The $\alpha$ parameter is a power-law index of $M_{0}$ and represents the early growth a halo experiences before the expansion of the Universe starts to slow it down. Finally, $T_{h}$ (the time that marks the co-evolution/growth of galaxies and haloes after the Big Bang) is found to be 150-300 million years., Comment: 18 pages, 8 figures, Accepted at MNRAS

We use spatially-resolved spectroscopy of a distant giant gravitational arc to test orientation effects on MgII absorption equivalent width (EW) and covering fraction (kappa) in the circumgalactic medium of a foreground star-forming galaxy (G1) at z~0.77. Forty-two spatially-binned arc positions uniformly sample impact parameters (D) to G1 between 10 and 30 kpc and azimuthal angles alpha between 30 and 90 degrees (minor axis). We find an EW-D anti-correlation, akin to that observed statistically in quasar absorber studies, and an apparent correlation of both EW and kappa with alpha, revealing a non-isotropic gas distribution. In line with our previous results on MgII kinematics suggesting the presence of outflows in G1, at minimum a simple 3-D static double-cone model (to represent the trace of bipolar outflows) is required to recreate the EW spatial distribution. The D and alpha values probed by the arc cannot confirm the presence of a disc, but the data highly disfavor a disc alone. Our results support the interpretation that the EW-alpha correlation observed statistically using other extant probes is partly shaped by bipolar metal-rich winds., Comment: Accepted in Monthly Notices of the Royal Astronomical Society

Based on a large group/cluster catalog recently constructed from the DESI Legacy Imaging Surveys DR9 using an extended halo-based group finder, we measure and model the group-galaxy weak lensing signals for groups/clusters in a few redshift bins within redshift range $0.1 \leqslant z<0.6$. Here, the background shear signals are obtained based on the DECaLS survey shape catalog derived with the \textsc{Fourier\_Quad} method. We divide the lens samples into 5 equispaced redshift bins and 7 mass bins, which allow us to probe the redshift and mass dependence of the lensing signals and hence the resulting halo properties. In addition to these sample selections, we have also checked the signals around different group centers, e.g., brightest central galaxy (BCG), luminosity weighted center and number weighted center. We use a lensing model that includes off-centering to describe the lensing signals we measure for all mass and redshift bins. The results demonstrate that our model predictions for the halo masses, bias and concentrations are stable and self-consistent among different samples for different group centers. Taking advantage of the very large and complete sample of groups/clusters, as well as the reliable estimation of their halo masses, we provide measurements of the cumulative halo mass functions up to redshift $z=0.6$, with a mass precision at $0.03\sim0.09$ dex., Comment: revised version submitted to ApJ

Using the extended halo-based group finder developed by Yang et al. (2021), which is able to deal with galaxies via spectroscopic and photometric redshifts simultaneously, we construct galaxy group and candidate protocluster catalogs in a wide redshift range ($0 < z < 6$) from the joint CFHT Large Area $U$-band Deep Survey (CLAUDS) and Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) deep data set. Based on a selection of 5,607,052 galaxies with $i$-band magnitude $m_{i} < 26$ and a sky coverage of $34.41\ {\rm deg}^2$, we identify a total of 2,232,134 groups, within which 402,947 groups have at least three member galaxies. We have visually checked and discussed the general properties of those richest groups at redshift $z>2.0$. By checking the galaxy number distributions within a $5-7\ h^{-1}\mathrm{Mpc}$ projected separation and a redshift difference $\Delta z \le 0.1$ around those richest groups at redshift $z>2$, we identified a list of 761, 343 and 43 protocluster candidates in the redshift bins $2\leq z<3$, $3\leq z<4$ and $z \geq 4$, respectively. In general, these catalogs of galaxy groups and protocluster candidates will provide useful environmental information in probing galaxy evolution along the cosmic time., Comment: 25 pages, 16 figures, 2 tables, accepted for publication in ApJ

We use weak lensing observations to make the first measurement of the characteristic depletion radius, one of the three radii that characterize the region where matter is being depleted by growing haloes. The lenses are taken from the halo catalog produced by the extended halo-based group/cluster finder applied to DESI Legacy Imaging Surveys DR9, while the sources are extracted from the DECaLS DR8 imaging data with the Fourier_Quad pipeline. We study halo masses $12 < \log ( M_{\rm grp} ~[{\rm M_{\odot}}/h] ) \leq 15.3$ within redshifts $0.2 \leq z \leq 0.3$. The virial and splashback radii are also measured and used to test the original findings on the depletion region. When binning haloes by mass, we find consistency between most of our measurements and predictions from the CosmicGrowth simulation, with exceptions to the lowest mass bins. The characteristic depletion radius is found to be roughly $2.5$ times the virial radius and $1.7 - 3$ times the splashback radius, in line with an approximately universal outer density profile, and the average enclosed density within the characteristic depletion radius is found to be roughly $29$ times the mean matter density of the Universe in our sample. When binning haloes by both mass and a proxy for halo concentration, we do not detect a significant variation of the depletion radius with concentration, on which the simulation prediction is also sensitive to the choice of concentration proxy. We also confirm that the measured splashback radius varies with concentration differently from simulation predictions., Comment: 17 pages, 13 figures, accepted: MNRAS May 03, 2022

A debate is emerging regarding the recent inconsistent results of different studies for the Cosmic Star Formation Rate Density (CSFRD) at high-z. We employ UV and IR datasets to investigate the star formation rate function (SFRF) at ${\rm z \sim 0-9}$. We find that the SFRFs derived from the dust corrected ${\rm UV}$ (${\rm UV_{corr}}$) data contradict those from IR on some key issues since they are described by different distributions (Schechter vs double-power law), imply different physics for galaxy formation (${\rm UV_{corr}}$ data suggest a SFR limit/strong mechanism that diminish the number density of high star forming systems with respect IR) and compare differently with the stellar mass density evolution obtained from SED fitting (${\rm UV_{corr}}$ is in agreement, while IR in tension up to 0.5 dex). However, both tracers agree on a constant CSFRD evolution at ${\rm z \sim 1-4}$ and point to a plateau instead of a peak. In addition, using both indicators we demonstrate that the evolution of the {\it observed} CSFRD can be described by only {\bf 2} parameters and a function that has the form of a Gamma distribution (${\bf \Gamma(a,bt)}$). In contrast to previous parameterizations used in the literature our framework connects the parameters to physical properties like the star formation rate depletion time and cosmic baryonic gas density. The build up of stellar mass occurs in $\Gamma(a,bt)$ distributed steps and is the result of gas consumption up to the limit that there is no eligible gas for SF at t = ${\rm \infty}$, resulting to a final cosmic stellar mass density of $\sim 0.5 \times 10^9 \, {\rm \frac{M_{\odot}}{Mpc^3}}$., Comment: 21 pages, 3 Figures, Accepted to ApJ

We present gravitational-arc tomography of the cool-warm enriched circumgalactic medium (CGM) of an isolated galaxy (``G1'') at $z \approx 0.77$. Combining VLT/MUSE adaptive-optics and Magellan/MagE echelle spectroscopy we obtain partially-resolved kinematics of MgII in absorption and [OII] in emission. The unique arc configuration allows us to probe 42 spatially independent arc positions transverse to G1, plus 4 positions in front of it. The transverse positions cover G1's minor and major axes at impact parameters of $\approx 10-30$ kpc and $\approx 60$ kpc, respectively. We observe a direct kinematic connection between the cool-warm enriched CGM (traced by MgII) and the interstellar medium (traced by [OII]). This provides strong evidence for the existence of an extended disc that co-rotates with the galaxy out to tens of kiloparsecs. The MgII velocity dispersion ($\sigma \approx 30-100$ km s$^{-1}$, depending on position) is of the same order as the modeled galaxy rotational velocity ($v_{\rm rot} \approx 80$ km s$^{-1}$), providing evidence for the presence of a turbulent and pressure-supported CGM component. We regard the absorption to be modulated by a galactic-scale outflow, as it offers a natural scenario for the observed line-of-sight dispersion and asymmetric profiles observed against both the arcs and the galaxy. An extended enriched co-rotating disc together with the signatures of a galactic outflow, are telltale signs of metal recycling in the $z\sim 1$ CGM., Comment: MNRAS accepted. Moderate changes after addressing referee feedback, including a correction in the systemic redshift; all qualitative results and conclusions remain the same

We extend the halo-based group finder developed by \citet[][]{Yang2005a} to use data {\it simultaneously} with either photometric or spectroscopic redshifts. A mock galaxy redshift survey constructed from a high-resolution N-body simulation is used to evaluate the performance of this extended group finder. For galaxies with magnitude ${\rm z\le 21}$ and redshift $0

We present a comprehensive analysis of supersymmetry anomalies in the free and massless Wess-Zumino (WZ) model in perturbation theory. At the classical level the model possesses ${\mathcal N}=1$ superconformal symmetry, which is partially broken by quantum anomalies. The form of the anomalies and the part of the symmetry they break depend on the multiplet of conserved currents used. It was previously shown that the R-symmetry anomaly of the conformal current multiplet induces an anomaly in Q-supersymmetry, which appears first in 4-point functions. Here we confirm this result by an explicit 1-loop computation using a supersymmetric Pauli-Villars regulator. The conformal current multiplet does not exist in the regulated theory because the regulator breaks conformal invariance, R-symmetry and S-supersymmetry explicitly. The minimal massive multiplet is the Ferrara-Zumino (FZ) one and the supersymmetry preserved by the regulator is a specific field dependent combination of Q- and S- supersymmetry of the conformal multiplet. While this supersymmetry is non anomalous, conformal invariance, R-symmetry and the original Q- and S-supersymmetries are explicitly broken by finite contact terms, both in the regulated and renormalized theories. A conformal current multiplet does exist for the renormalized theory and may be obtained from the FZ multiplet by a set of finite local counterterms that eliminate the explicit symmetry breaking, thus restoring superconformal invariance up to anomalies. However, this necessarily renders both Q- and S-supersymmetries anomalous, as is manifest starting at 4-point functions of conformal multiplet currents. The paper contains a detailed discussion of a number of issues and subtleties related to Ward identities that may be useful in a wider context., Comment: 37 pages + 45 pages appendices v2: minor improvements and amendments; v3: minor additions, version to appear in JHEP

We present the eddington bias corrected Specific Star Formation Rate Function (sSFRF) at different stellar mass scales from a sub-sample of the Sloan Digital Sky Survey Data Release DR7 (SDSS), which is considered complete both in terms of stellar mass (${\rm M_{\star}}$) and star formation rate (SFR). The above enable us to study qualitatively and quantitatively quenching, the distribution of passive/star-forming galaxies and perform comparisons with the predictions from state-of-the-art cosmological models, within the same ${\rm M_{\star}}$ and SFR limits. We find that at the low mass end (${\rm M_{\star}} = 10^{9.5} - 10^{10} \, {\rm M_{\odot}}$) the sSFRF is mostly dominated by star-forming objects. However, moving to the two more massive bins (${\rm M_{\star}} = 10^{10} - 10^{10.5} \, {\rm M_{\odot}}$ and ${\rm M_{\star}} = 10^{10.5} - 10^{11} \, {\rm M_{\odot}}$) a bi-modality with two peaks emerges. One peak represents the star-forming population, while the other describes a rising passive population. The bi-modal form of the sSFRFs is not reproduced by a range of cosmological simulations (e.g. Illustris, EAGLE, Mufasa, IllustrisTNG) which instead generate mostly the star-forming population, while a bi-modality emerges in others (e.g. L-Galaxies, Shark, Simba). Our findings reflect the need for the employed quenching schemes in state-of-the-art models to be reconsidered, involving prescriptions that allow "quenched galaxies" to retain a small level of SF activity (sSFR $=$ ${\rm 10^{-11} {\rm yr^{-1}}}$-${\rm 10^{-12} {\rm yr^{-1}}}$) and generate an adequate passive population/bi-modality even at intermediate masses (${\rm M_{\star}} = 10^{10} - 10^{10.5} \, {\rm M_{\odot}}$)., Comment: 13 pages, 4 Figures, Accepted to MNRAS

We combine data from the Sloan Digital Sky Survey (SDSS) and the Arecibo Legacy Fast ALFA Survey (ALFALFA) to establish an empirical model for the HI gas content within dark matter halos. A cross-match between our SDSS DR7 galaxy group sample and the ALFALFA HI sources provides a catalog of 16,520 HI-galaxy pairs within 14,270 galaxy groups (halos). Using these matched pairs, we model the HI gas mass distributions within halos using two components: 1) {\it in situ} galaxy relations that involve the HI masses, colors $({\rm g-r})$ and stellar masses 2) an {\it ex situ} dependence of the HI mass on the halo mass/environment. We find that if we solely use galaxy associated scaling relations to predict the HI gas distribution (solely component 1), the number of HI detections is significantly over-predicted with respect the ALFALFA observations. We introduce a concept for the survival of the HI masses/members within halos of different masses labelled as the `efficiency' factor, in order to describe the probability that a halo has in retaining its HI detections. Taking the above consideration into account we construct a `halo based HI mass model' which does not only predict the HI masses of galaxies, but also yields similar number, stellar, halo mass and satellite fraction distributions to the HI detections retrieved from observational data., Comment: 13 pages, 11 figures

The star formation rate function (SFRF) and specific star formation rate function (sSFRF) from the observation are impacted by the Eddington bias, due to the uncertainties on the estimated SFR. We develop a novel method to correct the Eddington bias and obtained the intrinsic SFRF and sSFRF from the Sloan Digital Sky Survey Data Release 7. The intrinsic SFRF is in good agreement with measurements from previous data in the literature that relied on UV SFRs but its high star-forming end is slightly lower than those IR and radio tracers. We demonstrate that the intrinsic sSFRF from SDSS has a bi-modal form with the one peak found at ${\rm sSFR \sim 10^{-9.7} yr^{-1}}$ representing the star-forming objects while the other peak is found at ${\rm sSFR \sim 10^{-12} yr^{-1}}$ representing the quenched population. Furthermore, we compare our observations with the predictions from the IllustrisTNG and Illustris simulations and show that the ``TNG'' model performs much better than its predecessor. However, we show that the simulated SFRF and cosmic star formation density (CSFRD) of TNG simulations are highly dependent on resolution, reflecting the limitations of the model and today state-of-the-art simulations. We demonstrate that the bi-modal, two peaked sSFRF implied by the SDSS observations does not appear in TNG regardless of the adopted box-size or resolution. This tension reflects the need for inclusion of an additional efficient quenching mechanism to the TNG model., Comment: 19 pages, 10 figures, accepted for publication in RAA

We present infrared luminosity functions and dust mass functions for the EAGLE cosmological simulation, based on synthetic multi-wavelength observations generated with the SKIRT radiative transfer code. In the local Universe, we reproduce the observed infrared luminosity and dust mass functions very well. Some minor discrepancies are encountered, mainly in the high luminosity regime, where the EAGLE-SKIRT luminosity functions mildly but systematically underestimate the observed ones. The agreement between the EAGLE-SKIRT infrared luminosity functions and the observed ones gradually worsens with increasing lookback time. Fitting modified Schechter functions to the EAGLE-SKIRT luminosity and dust mass functions at different redshifts up to $z=1$, we find that the evolution is compatible with pure luminosity/mass evolution. The evolution is relatively mild: within this redshift range, we find an evolution of $L_{\star,250}\propto(1+z)^{1.68}$, $L_{\star,\text{TIR}}\propto(1+z)^{2.51}$ and $M_{\star,\text{dust}}\propto(1+z)^{0.83}$ for the characteristic luminosity/mass. For the luminosity/mass density we find $\varepsilon_{250}\propto(1+z)^{1.62}$, $\varepsilon_{\text{TIR}}\propto(1+z)^{2.35}$ and $\rho_{\text{dust}}\propto(1+z)^{0.80}$, respectively. The mild evolution of the dust mass density is in relatively good agreement with observations, but the slow evolution of the infrared luminosity underestimates the observed luminosity evolution significantly. We argue that these differences can be attributed to increasing limitations in the radiative transfer treatment due to increasingly poorer resolution, combined with a slower than observed evolution of the SFR density in the EAGLE simulation and the lack of AGN emission in our EAGLE-SKIRT post-processing recipe., Comment: 13 pages, 6 pages, accepted for publication in MNRAS

We compare the spectral energy distributions (SEDs) and inferred physical properties for simulated and observed galaxies at low redshift. We exploit UV-submillimetre mock fluxes of ~7000 z=0 galaxies from the EAGLE suite of cosmological simulations, derived using the radiative transfer code SKIRT. We compare these to ~800 observed galaxies in the UV-submillimetre range, from the DustPedia sample of nearby galaxies. To derive global properties, we apply the SED fitting code CIGALE consistently to both data sets, using the same set of ~80 million models. The results of this comparison reveal overall agreement between the simulations and observations, both in the SEDs and in the derived physical properties, with a number of discrepancies. The optical and far-infrared regimes, and the scaling relations based upon the global emission, diffuse dust and stellar mass, show high levels of agreement. However, the mid-infrared fluxes of the EAGLE galaxies are overestimated while the far-UV domain is not attenuated enough, compared to the observations. We attribute these discrepancies to a combination of galaxy population differences between the samples, and limitations in the subgrid treatment of star-forming regions in the EAGLE-SKIRT post-processing recipe. Our findings show the importance of detailed radiative transfer calculations and consistent comparison, and provide suggestions for improved numerical models., Comment: 17 pages, 14 figures, accepted for publication in MNRAS

There is a severe tension between the observed star formation rate (SFR) - stellar mass (${\rm M}_{\star}$) relations reported by different authors at $z = 1-4$. In addition, the observations have not been successfully reproduced by state-of-the-art cosmological simulations which tend to predict a factor of 2-4 smaller SFRs at a fixed ${\rm M}_{\star}$. We examine the evolution of the SFR$-{\rm M}_{\star}$ relation of $z = 1-4 $ galaxies using the SKIRT simulated spectral energy distributions of galaxies sampled from the EAGLE simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR$-{\rm M}_{\star}$ relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 ${\rm \mu m}$, MIPS - 24, 70 and 160 ${\rm \mu m}$ or SPIRE - 250, 350, 500 ${\rm \mu m}$) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by -0.5 dex at $z \sim 4$ but overpredict the measurements by 0.3 dex at $z \sim 1$. Relations relying on dust-corrected rest-frame UV luminosities, are flatter since they overpredict/underpredict SFRs for low/high star forming objects and yield deviations from the intrinsic relation from 0.10 dex to -0.13 dex at $z \sim 4$. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs., Comment: 15 pages, 6 Figures, 4 tables, accepted for publication in MNRAS, Figures 1, 2, 3, 4 summarize results

We present spatially-resolved echelle spectroscopy of an intervening MgII-FeII-MgI absorption-line system detected at $z_{\rm abs}=0.73379$ toward the giant gravitational arc PSZ1 G311.65-18.48. The absorbing gas is associated to an inclined disk-like star-forming galaxy, whose major axis is aligned with the two arc-segments reported here. We probe in absorption the galaxy's extended disk continuously, at $\approx 3$ kpc sampling, from its inner region out to $15\times$ the optical radius. We detect strong ($W_0^{2796}>0.3$ \r{A}) coherent absorption along $13$ independent positions at impact parameters $D=0$--$29$ kpc on one side of the galaxy, and no absorption at $D=28$--$57$ kpc on the opposite side (all de-lensed distances at $z_{\rm abs}$). We show that: (1) the gas distribution is anisotropic; (2) $W_0^{2796}$, $W_0^{2600}$, $W_0^{2852}$, and the ratio $W_0^{2600}\!/W_0^{2796}$, all anti-correlate with $D$; (3) the $W_0^{2796}$-$D$ relation is not cuspy and exhibits significantly less scatter than the quasar-absorber statistics; (4) the absorbing gas is co-rotating with the galaxy out to $D \lesssim 20$ kpc, resembling a `flat' rotation curve, but at $D\gtrsim 20$ kpc velocities decline below the expectations from a 3D disk-model extrapolated from the nebular [OII] emission. These signatures constitute unambiguous evidence for rotating extra-planar diffuse gas, possibly also undergoing enriched accretion at its edge. Arguably, we are witnessing some of the long-sought processes of the baryon cycle in a single distant galaxy expected to be representative of such phenomena., Comment: Accepted for publication in MNRAS

The scatter (${\rm\sigma_{\text{sSFR}}}$) of the specific star formation rates (sSFRs) of galaxies is a measure of the diversity in their star formation histories (SFHs) at a given mass. In this paper we employ the EAGLE simulations to study the dependence of the ${\rm \sigma_{\text{sSFR}}}$ of galaxies on stellar mass (${\rm M_{\star}}$) through the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation in $ {\rm z \sim 0-4}$. We find that the relation evolves with time, with the dispersion depending on both stellar mass and redshift. The models point to an evolving U-shape form for the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation with the scatter being minimal at a characteristic mass $M^{\star}$ of ${\rm 10^{9.5}}$ ${\rm M_{\odot}}$ and increasing both at lower and higher masses. This implication is that the diversity of SFHs increases towards both at the low- and high-mass ends. We find that active galactic nuclei feedback is important for increasing the ${\rm \sigma_{\text{sSFR}}}$ for high mass objects. On the other hand, we suggest that SNe feedback increases the ${\rm \sigma_{\text{sSFR}}}$ of galaxies at the low-mass end. We also find that excluding galaxies that have experienced recent mergers does not significantly affect the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation. Furthermore, we employ the combination of the EAGLE simulations with the radiative transfer code SKIRT to evaluate the effect of SFR/stellar mass diagnostics in the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation and find that the ${\rm SFR/M_{\star}}$ methodologies (e.g. SED fitting, UV+IR, UV+IRX-$\beta$) widely used in the literature to obtain intrinsic properties of galaxies have a large effect on the derived shape and normalization of the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation., Comment: 23 pages, 10 figures, Accepted to ApJ

We study the shape of the gas-phase mass-metallicity relation (MZR) of a combined sample of present-day dwarf and high-mass star-forming galaxies using IZI, a Bayesian formalism for measuring chemical abundances presented in Blanc et al. 2015. We observe a characteristic stellar mass scale at $M_* \simeq 10^{9.5}$M$_{\odot}$, above which the ISM undergoes a sharp increase in its level of chemical enrichment. In the $10^{6}-10^{9.5}$M$_{\odot}$ range the MZR follows a shallow power-law ($Z\propto M^{\alpha}_*$) with slope $\alpha=0.14\pm0.08$. At approaching $M_* \simeq 10^{9.5}$M$_{\odot}$ the MZR steepens significantly, showing a slope of $\alpha=0.37\pm0.08$ in the $10^{9.5}-10^{10.5}$M$_{\odot}$ range, and a flattening towards a constant metallicity at higher stellar masses. This behavior is qualitatively different from results in the literature that show a single power-law MZR towards the low mass end. We thoroughly explore systematic uncertainties in our measurement, and show that the shape of the MZR is not induced by sample selection, aperture effects, a changing N/O abundance, the adopted methodology used to construct the MZR, secondary dependencies on star formation activity, nor diffuse ionized gas (DIG) contamination, but rather on differences in the method used to measure abundances. High resolution hydrodynamical simulations can qualitatively reproduce our result, and suggest a transition in the ability of galaxies to retain their metals for stellar masses above this threshold. The MZR characteristic mass scale also coincides with a transition in the scale height and clumpiness of cold gas disks, and a typical gas fraction below which the efficiency of star formation feedback for driving outflows is expected to decrease sharply., Comment: 24 pages, 11 figures, 4 tables, accepted for publication in ApJ

The distribution of gas-phase abundances in galaxy disks encodes the history of nucleosynthesis and transport through the interstellar medium (ISM) over cosmic time. Multi-object and high resolution integral-field spectroscopy have started to measure these distributions across hundreds of HII regions individually resolved at $\lesssim 100$ pc scales in a handful of objects, but in the coming decade these studies will expand to larger samples of galaxies. This will allow us to understand the role of feedback and turbulence in driving the mixing and diffusion of metals in the ISM, and statistically assess the role of galaxy environment and disk dynamics in modifying how mixing proceeds. Detailed searches for over- and under-enriched regions can address to what extent star formation is triggered by previous generations of star formation and by pristine and recycled gas flows. Local galaxies, for which these detailed measurements will be possible within the next decade, will inform the interpretation of integrated measurements at high-z, where very different dynamical gas-rich environments are found in early disk galaxies. Currently, progress in the field is severely hampered by the 0.2-0.3 dex level systematic uncertainties plaguing nebular abundance diagnostics. Improving our detailed understanding of ionized nebulae at $<$20 pc scales will help us find a solution to this problem, which will prove key to the study of metal enrichment and mixing across the galaxy population in the next decade., Comment: Submitted to the Astro2020 call for science white papers

We show that supersymmetry is anomalous in ${\cal N}=1$ superconformal quantum field theories (SCFTs) with an anomalous R-symmetry. This anomaly was originally found in holographic SCFTs at strong coupling. Here we show that this anomaly is present in general and demonstrate it for the massless superconformal Wess-Zumino model via a one loop computation. The anomaly appears first in four-point functions of two supercurrents with either two R-currents or with an R-current and an energy momentum tensor. In fact, the Wess-Zumino consistency conditions together with the standard R-symmetry anomaly imply the existence of the anomaly. We outline the implications of this anomaly., Comment: 5 pages; v2 - contextual discussion and references added; v3 - final PRL version, with comments added

Recently a number of studies have proposed that the dispersion along the star formation rate - stellar mass relation ($\sigma_{\mathrm{sSFR}}$-M$_{*}$) is indicative of variations in star-formation history (SFH) driven by feedback processes. They found a 'U'-shaped dispersion and attribute the increased scatter at low and high stellar masses to stellar and active galactic nuclei feed-back respectively. However, measuring $\sigma_{\mathrm{sSFR}}$ and the shape of the $\sigma_{\mathrm{sSFR}}$-M$_{*}$ relation is problematic and can vary dramatically depending on the sample selected, chosen separation of passive/star-forming systems, and method of deriving star-formation rates ($i.e.$ H$\alpha$ emission vs spectral energy distribution fitting). As such, any astrophysical conclusions drawn from measurements of $\sigma_{\mathrm{sSFR}}$ must consider these dependencies. Here we use the Galaxy And Mass Assembly survey to explore how $\sigma_{\mathrm{sSFR}}$ varies with SFR indicator for a variety of selections for disc-like `main sequence' star-forming galaxies including colour, star-formation rate, visual morphology, bulge-to-total mass ratio, S\'{e}rsic index and mixture modelling. We find that irrespective of sample selection and/or SFR indicator, the dispersion along the sSFR-M$_{*}$ relation does follow a 'U'-shaped distribution. This suggests that the shape is physical and not an artefact of sample selection or method. We then compare the $\sigma_{\mathrm{sSFR}}$-M$_{*}$ relation to state-of-the-art hydrodynamical and semi-analytic models and find good agreement with our observed results. Finally, we find that for group satellites this 'U'-shaped distribution is not observed due to additional high scatter populations at intermediate stellar masses., Comment: 21 pages, 15 figures, Accepted MNRAS

We investigate the evolution of the galaxy Star Formation Rate Function (SFRF) and Cosmic Star Formation Rate Density (CSFRD) of $z\sim 0-8 $ galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. In addition, we present a compilation of UV, IR and H$\alpha$ SFRFs and compare these with the predictions from the EAGLE suite of cosmological hydrodynamic simulations. We find that the constraints implied by different indicators are inconsistent with each other for the highest star-forming objects at z < 2, a problem that is possibly related to selection biases and the uncertainties of dust attenuation effects. EAGLE's feedback parameters were calibrated to reproduce realistic galaxy sizes and stellar masses at z = 0.1. In this work we test if and why those choices yield realistic Star Formation Rates (SFRs) for $z \sim 0-8$ as well. We demonstrate that SNe feedback plays a major role at setting the abundance of galaxies at all star-forming regimes, especially at high redshifts. On the contrary, Active Galactic Nuclei (AGN) feedback becomes more prominent at lower redshifts and is a major mechanism that affects only the highest star-forming systems. Furthermore, we find that galaxies with SFR $\sim 1-10 \, {\rm M_{\odot} \, yr^{-1}}$ dominate the CSFRD at redshifts z < 5, while rare high star-forming galaxies (SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$) contribute significantly only briefly around the peak era ($z \sim 2$) and then are quenched by AGN feedback. In the absence of this prescription objects with SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$ would dominate the CSFRD, while the cosmic budget of star formation would be extremely high. Finally, we demonstrate that the majority of the cosmic star formation occurs in relatively rare high mass halos ($ {\rm M_{Halo}} \sim 10^{11-13} \, {\rm M_{\odot}}$) even at the earliest epochs., Comment: 21 pages, 12 figures, accepted for publication in MNRAS

An interesting transient has been detected in one of our three Dark Energy Camera deep fields. Observations of these deep fields take advantage of the high red sensitivity of DECam on the Cerro Tololo Interamerican Observatory Blanco telescope. The survey includes the Y band with rest wavelength 1430{\AA} at z = 6. Survey fields (the Prime field 0555-6130, the 16hr field 1600-75 and the SUDSS New Southern Field) are deeper in Y than other infrared surveys. They are circumpolar, allowing all night to be used efficiently, exploiting the moon tolerance of 1 micron observations to minimize conflict with the Dark Energy Survey. As an i-band dropout (meaning that the flux decrement shortward of Lyman alpha is in the i bandpass), the transient we report here is a supernova candidate with z ~ 6, with a luminosity comparable to the brightest known current epoch superluminous supernova (i.e., ~ 2 x 10^11 solar luminosities)., Comment: Reference added

Superluminous supernovae are beginning to be discovered at redshifts as early as the epoch of reionization. A number of candidate mechanisms is reviewed, together with the discovery programs., Comment: for publication in PASA

We investigate the evolution of the galaxy Star Formation Rate Function (SFRF) and Cosmic Star Formation Rate Density (CSFRD) of $z\sim 1-4 $ galaxies, using cosmological Smoothed Particle Hydrodynamic (SPH) simulations and a compilation of UV, IR and H$\alpha$ observations. These tracers represent different populations of galaxies with the IR light being a probe of objects with high star formation rates and dust contents, while UV and H$\alpha$ observations provide a census of low star formation galaxies where mild obscuration occurs. We compare the above SFRFs with the results of SPH simulations run with the code {\small{P-GADGET3(XXL)}}. We focus on the role of feedback from Active Galactic Nuclei (AGN) and supernovae in form of galactic winds. The AGN feedback prescription that we use decreases the simulated CSFRD at $z < 3$ but is not sufficient to reproduce the observed evolution at higher redshifts. We explore different wind models and find that the key factor for reproducing the evolution of the observed SFRF and CSFRD at $z \sim1-4$ is the presence of a feedback prescription that is prominent at high redshifts ($z \ge 4$) and becomes less efficient with time. We show that variable galactic winds which are efficient at decreasing the SFRs of low mass objects are quite successful in reproducing the observables., Comment: 20 pages, 7 figures, 10 tables, accepted for publication in MNRAS

The relation between the Star Formation Rate (SFR) and stellar mass (${\rm M}_{\star}$) of galaxies represents a fundamental constraint on galaxy formation and has been studied extensively both in observations and cosmological simulations. However, the observed amplitude has not been successfully reproduced in simulations, indicating either that the halo accretion history and baryonic physics are poorly modeled or that observations contain biases. We examine the evolution of the SFR$-{\rm M}_{\star}$ relation of $z\sim1-4 $ galaxies and display the inconsistency between observed relations that are obtained using different techniques. We employ cosmological hydrodynamic simulations and compare these with a range of observed SFR$-{\rm M}_{\star}$ relations. We find that numerical results are consistent with observations that use Spectral Energy Distribution (SED) techniques to estimate star formation rates and dust corrections. On the contrary, simulations are not able to reproduce results that were obtained by combining only UV and IR luminosities. These imply SFRs at a fixed stellar mass that are larger almost by a factor of 5 than those of SED measurements for $z \sim1.5-4$. Furthermore, we find remarkable agreement between the numerical results from various authors who have employed different cosmological codes and run simulations with different resolutions. This is interesting for two reasons. A) simulations can produce realistic populations of galaxies within representative cosmological volumes even at relatively modest resolutions. B) It is likely that current numerical codes that rely on similar subgrid multiphase Inter-Stellar Medium (ISM) models and are tuned to reproduce statistical properties of galaxies, produce similar results for the SFR$-{\rm M}_{\star}$ relation by construction, regardless of resolution, box size and, to some extent, the adopted feedback prescriptions., Comment: 16 pages, 6 figures, 3 tables, accepted for publication in PASA

We investigate the evolution of the star formation rate-stellar mass relation (SFR-M*) and Galaxy Stellar Mass Function (GSMF) of z ~ 4-7 galaxies, using cosmological simulations run with the smoothed particle hydrodynamics code P-GADGET3(XXL). We explore the effects of different feedback prescriptions (supernova driven galactic winds and AGN feedback), initial stellar mass functions and metal cooling. We show that our fiducial model, with strong energy-driven winds and early AGN feedback, is able to reproduce the observed stellar mass function obtained from Lyman-break selected samples of star forming galaxies at redshift 6 < z < 7. At z ~ 4, observed estimates of the GSMF vary according to how the sample was selected. Our simulations are more consistent with recent results from K-selected samples, which provide a better proxy of stellar masses and are more complete at the high mass end of the distribution. We find that in some cases simulated and observed SFR-M* relations are in tension, and this can lead to numerical predictions for the GSMF in excess of the GSMF observed. By combining the simulated SFR(M*) relationship with the observed star formation rate function at a given redshift, we argue that this disagreement may be the result of the uncertainty in the SFR-M* (Luv-M*) conversion. Our simulations predict a population of faint galaxies not seen by current observations., Comment: 23 Pages, 13 figures, modified to match accepted version to MNRAS