Search for Hα Emitters at z ∼ 7.8: A Constraint on the Hα-based Star Formation Rate Density

We search for Hα emitters at z ∼ 7.8 in four gravitationally lensed fields observed in the Hubble Frontier Fields program. We use the Lyman break method to select galaxies at the target redshift and perform photometry in the Spitzer/IRAC 5.8 μm band to detect Hα emission from the candidate galaxies. We find no significant detections of counterparts in the IRAC 5.8 μm band, and this gives a constraint on the Hα luminosity function (LF) at z ∼ 7.8. We compare the constraint with previous studies based on rest-frame UV and far-infrared observations using the correlation between the Hα luminosity and the star formation rate. Additionally, we convert the constraint on the Hα LF into an upper limit for the star formation rate density (SFRD) at this epoch assuming the shape of the LF. We examine two types of parameterization of the LF and obtain an upper limit for the SFRD of log 10 ( ρ SFR [ M ⊙ yr − 1 Mpc − 3 ] ) ≲ − 1.1 at z ∼ 7.8. With this constraint on the SFRD, we present an independent probe into the total star formation activity including dust-obscured and unobscured star formation at the Epoch of Reionization.

What is Important? Morphological Asymmetries are Useful Predictors of Star Formation Rates of Star-forming Galaxies in SDSS Stripe 82

The morphology and structure of galaxies reflect their star formation and assembly histories. We use the framework of mutual information (MI) to quantify the interdependence among several structural variables and to rank them according to their relevance for predicting the specific star formation rate (SSFR) by comparing the MI of the predictor variables with the SSFR and penalizing variables that are redundant. We apply this framework to study ∼3700 face-on star-forming galaxies (SFGs) with varying degrees of bulge dominance and central concentration and with stellar mass M ⋆ ≈ 109 M ⊙−5 × 1011 M ⊙ at redshift z = 0.02–0.12. We use the Sloan Digital Sky Survey (SDSS) Stripe 82 deep i-band imaging data, which improve measurements of asymmetry and bulge dominance indicators. We find that star-forming galaxies are a multiparameter family. In addition to M ⋆, asymmetry emerges as the most powerful predictor of SSFR residuals of SFGs, followed by bulge prominence/concentration. Star-forming galaxies with higher asymmetry and stronger bulges have higher SSFR at a given M ⋆. The asymmetry reflects both irregular spiral arms and lopsidedness in seemingly isolated SFGs and structural perturbations by galaxy interactions or mergers.

Does the reionization model influence the constraints on dark matter decay or annihilation?

If dark matter decay or annihilate, a large amount of energy and particles would be released into the cosmic plasma. Therefore, they could modify the thermal and ionization history of our universe, then leave footprints on the cosmic microwave background power spectra. In this paper, we take dark matter annihilation as an example and investigate whether different reionization models influence the constraints on dark matter annihilation. We consider the ionization history including both dark matter annihilation and star formation, then put constraints on DM annihilation. Combining the latest Planck data, BAO data, SNIa measurement, Q HII constraints from observations of quasars, as well as the star formation rate density from UV and IR data, the optical depth is τ = 0.0571+0.0005 -0.0006 at 68%C.L. and the upper limit of ϵ0 f d reads 2.7765 × 10-24 at 95%C.L.. By comparison, we also constrain dark matter annihilation in the instantaneous reionization model from the same data combination except the Q HII constraints and star formation rate density. We get τ = 0.0559+0.0069 -0.0076 at 68%C.L. and the upper limit of ϵ0 f d is 2.8468 × 10-24 at 95%C.L.. This indicates various reionization models have little influence (≲ 2.5%) on constraining parameters of dark matter decay or annihilation.

Extensive Lensing Survey of Optical and Near-infrared Dark Objects (El Sonido): HST H-faint Galaxies behind 101 Lensing Clusters

We present a Spitzer/IRAC survey of H-faint (H 160 ≳ 26.4, < 5σ) sources in 101 lensing cluster fields. Across a CANDELS/Wide-like survey area of ∼648 arcmin2 (effectively ∼221 arcmin2 in the source plane), we have securely discovered 53 sources in the IRAC Channel-2 band (CH2, 4.5 μm; median CH2 = 22.46 ± 0.11 AB mag) that lack robust HST/WFC3-IR F160W counterparts. The most remarkable source in our sample, namely ES-009 in the field of Abell 2813, is the brightest H-faint galaxy at 4.5 μm known so far (CH2 = 20.48 ± 0.03 AB mag). We show that the H-faint sources in our sample are massive (median M star = 1010.3±0.3 M ⊙), star-forming (median star formation rate = 100 − 40 + 60 M ⊙ yr−1), and dust-obscured (A V = 2.6 ± 0.3) galaxies around a median photometric redshift of z = 3.9 ± 0.4. The stellar continua of 14 H-faint galaxies can be resolved in the CH2 band, suggesting a median circularized effective radius (R e,circ; lensing corrected) of 1.9 ± 0.2 kpc and <1.5 kpc for the resolved and whole samples, respectively. This is consistent with the sizes of massive unobscured galaxies at z ∼ 4, indicating that H-faint galaxies represent the dusty tail of the distribution of a wider galaxy population. Comparing with the ALMA dust continuum sizes of similar galaxies reported previously, we conclude that the heavy dust obscuration in H-faint galaxies is related to the compactness of both stellar and dust continua (R e,circ ∼ 1 kpc). These H-faint galaxies make up 16 − 7 + 13 % of the galaxies in the stellar-mass range of 1010 − 1011.2 M ⊙ at z = 3 ∼ 5, contributing to 8 − 4 + 8 % of the cosmic star formation rate density in this epoch and likely tracing the early phase of massive galaxy formation.

sígame v3: Gas Fragmentation in Postprocessing of Cosmological Simulations for More Accurate Infrared Line Emission Modeling

We present an update to the framework called Simulator of Galaxy Millimeter/submillimeter Emission (sígame). sígame derives line emission in the far-infrared (FIR) for galaxies in particle-based cosmological hydrodynamics simulations by applying radiative transfer and physics recipes via a postprocessing step after completion of the simulation. In this version, a new technique is developed to model higher gas densities by parameterizing the probability distribution function (PDF) of the gas density in higher-resolution simulations run with the pseudo-Lagrangian, Voronoi mesh code arepo. The parameterized PDFs are used as a look-up table, and reach higher densities than in previous work. sígame v3 is tested on redshift z = 0 galaxies drawn from the simba cosmological simulation for eight FIR emission lines tracing vastly different phases of the interstellar medium. This version of sígame includes dust radiative transfer with Skirt and high-resolution photoionization models with Cloudy, the latter sampled according to the density PDF of the arepo simulations to augment the densities in the cosmological simulation. The quartile distributions of the predicted line luminosities overlap with the observed range for nearby galaxies of similar star formation rate (SFR) for all but two emission lines: [O i]63 and CO(3–2), which are overestimated by median factors of 1.3 and 1.0 dex, respectively, compared to the observed line–SFR relation of mixed-type galaxies. We attribute the remaining disagreement with observations to the lack of precise attenuation of the interstellar light on sub-grid scales (≲200 pc) and differences in sample selection.

Impact of relativistic jets on the star formation rate: A turbulence-regulated framework

We apply a turbulence-regulated model of star formation to calculate the star formation rate (SFR) of dense star-forming clouds in simulations of jet-ISM interactions. The method isolates individual clumps and accounts for the impact of virial parameter and Mach number of the clumps on the star formation activity. This improves upon other estimates of the SFR in simulations of jet–ISM interactions, which are often solely based on local gas density, neglecting the impact of turbulence. We apply this framework to the results of a suite of jet-ISM interaction simulations to study how the jet regulates the SFR both globally and on the scale of individual star-forming clouds. We find that the jet strongly affects the multi-phase ISM in the galaxy, inducing turbulence and increasing the velocity dispersion within the clouds. This causes a global reduction in the SFR compared to a simulation without a jet. The shocks driven into clouds by the jet also compress the gas to higher densities, resulting in local enhancements of the SFR. However, the velocity dispersion in such clouds is also comparably high, which results in a lower SFR than would be observed in galaxies with similar gas mass surface densities and without powerful radio jets. We thus show that both local negative and positive jet feedback can occur in a single system during a single jet event, and that the star-formation rate in the ISM varies in a complicated manner that depends on the strength of the jet-ISM coupling and the jet break-out time-scale.

A flexible subhalo abundance matching model for galaxy clustering in redshift space

ABSTRACT We develop an extension of subhalo abundance matching (SHAM) capable of accurately reproducing the real and redshift-space clustering of galaxies in a state-of-the-art hydrodynamical simulation. Our method uses a low-resolution gravity-only simulation and it includes orphan and tidal disruption prescriptions for satellite galaxies, and a flexible amount of galaxy assembly bias. Furthermore, it includes recipes for star formation rate (SFR) based on the dark matter accretion rate. We test the accuracy of our model against catalogues of stellar-mass- and SFR-selected galaxies in the TNG300 hydrodynamic simulation. By fitting a small number of free parameters, our extended SHAM reproduces the projected correlation function and redshift-space multipoles for number densities $10^{-3} - 10^{-2}\, h^{3}{\rm Mpc}^{-3}$, at z = 1 and z = 0, and for scales r ∈ [0.3 − 20]h−1Mpc. Simultaneously, the SHAM results also retrieve the correct halo occupation distribution, the level of galaxy assembly bias, and higher order statistics present in the TNG300 galaxy catalogues. As an application, we show that our model simultaneously fits the projected correlation function of the SDSS in three disjoint stellar mass bins, with an accuracy similar to that of TNG300 galaxies. This SHAM extension can be used to get accurate clustering prediction even when using low and moderate-resolution simulations.