Search

We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a ultraviolet (UV)-selected galaxy sample with the Atacama Large Millimeter Array (ALMA) spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at z = 6.4–7.7 with [C ii] 158 $\rm{\mu m}$ line scans. 16 sources exhibit dust detection, 15 of which are also spectroscopically confirmed through the [C ii] line. The infrared (IR) luminosities of the sample range from log LIR/L⊙ = 11.4 to 12.2. Using the UV luminosity function as a proxy to derive the effective volume for each of our target sources, we derive IRLF estimates, both for detections and for the full sample including IR luminosity upper limits. The resulting IRLFs are well reproduced by a Schechter function with the characteristic luminosity of $\log L_{*}/\mathrm{ L}_\odot =11.6^{+0.2}_{-0.1}$ . Our observational results are in broad agreement with the average of predicted IRLFs from simulations at z ∼ 7. Conversely, our IRLFs lie significantly below lower redshift estimates, suggesting a rapid evolution from z ∼ 4 to z ∼ 7, into the reionization epoch. The IR obscured contribution to the cosmic star formation rate density at z ∼ 7 amounts to $\mathrm{log(SFRD/{\rm M}_{\odot }\,yr^{-1}\,Mpc^{-3}) = -2.66^{+0.17}_{-0.14} }$ that is at least ∼10 per cent of UV-based estimates. We conclude that the presence of dust is already abundant in the EoR and discuss the possibility of unveiling larger samples of dusty galaxies with future ALMA and JWST observations.

We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rm\mu m$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a significant stellar mass of $M_{*}=8^{+4}_{-2}\times10^{9}$ M$_\odot$. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR $=199^{+101}_{-63}$ M$_\odot$ yr$^{-1}$. This is about four times the SFR estimated from an extrapolated main-sequence. We also infer a [CII]-based molecular gas mass of $M_{\rm H_2}=5.1^{+5.1}_{-2.6}\times10^{10}$ $M_\odot$, implying a molecular gas depletion time of $ t_{\rm depl, H_2}=0.3^{+0.3}_{-0.2}$ Gyr. We observe a [CII] velocity gradient consistent with disc rotation, but given the current resolution we cannot rule out a more complex velocity structure such as a merger. The spectrum exhibits excess [CII] emission at large positive velocities ($\sim500$ km s$^{-1}$), which we interpret as either a merging companion or an outflow. In the outflow scenario, we derive a lower limit of the mass outflow rate of 200 M$_\odot$ yr$^{-1}$, which is consistent with expectations for a star formation-driven outflow. Given its large stellar mass, SFR and molecular gas reservoir $\sim700$ Myr after the Big Bang, we explore the future evolution of REBELS-25. Considering a simple, conservative model assuming an exponentially declining star formation history, constant star formation efficiency, and no additional gas inflow, we find that REBELS-25 has the potential to evolve into a galaxy consistent with the properties of high-mass quiescent galaxies recently observed at $z\sim4$., Comment: Accepted for publication in MNRAS. 21 pages, 8 figures

We present the analysis of $\sim 100$pc-scale compact radio continuum sources detected in 63 local (Ultra) Luminous Infrared Galaxies (U/LIRGs; $L_{\rm IR} \ge 10^{11} L_\odot$), using FWHM $\lesssim 0''.1 - 0''.2$ resolution 15 and 33 GHz observations with the Karl G. Jansky Very Large Array. We identify a total of 133 compact radio sources with effective radii of 8 - 170pc, which are classified into four main categories -- "AGN" (AGN), "AGN/SBnuc" (AGN-starburst composite nucleus), "SBnuc" (starburst nucleus) and "SF" (star-forming clumps) -- based on ancillary datasets and the literature. We find that "AGN" and "AGN/SBnuc" more frequently occur in late-stage mergers and have up to 3 dex higher 33 GHz luminosities and surface densities compared with "SBnuc" and "SF", which may be attributed to extreme nuclear starburst and/or AGN activity in the former. Star formation rates (SFRs) and surface densities ($\Sigma_{\rm SFR}$) are measured for "SF" and "SBnuc" using both the total 33 GHz continuum emission (SFR $\sim 0.14 - 13$ M$_\odot$ yr$^{-1}$, $\Sigma_{\rm SFR} \sim 13 - 1600$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$) and the thermal free-free emission from HII regions (median SFR$_{\rm th} \sim 0.4$ M$_\odot$ yr$^{-1}$, $\Sigma_{\rm SFR_{th}} \sim 44$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$). These values are 1 - 2 dex higher than those measured for similar-sized clumps in nearby normal (non-U/LIRGs). The latter also have much flatter median 15 - 33 GHz spectral index ($\sim -0.08$) compared with "SBnuc" and "SF" ($\sim -0.46$), which may reflect higher non-thermal contribution from supernovae and/or ISM densities in local U/LIRGs that directly result from and/or lead to their extreme star-forming activities on 100\,pc scales., Comment: Accepted for publication in ApJ

We present the average [C ii] 158 μm emission line sizes of UV-bright star-forming galaxies at z ∼ 7. Our results are derived from a stacking analysis of [C ii] 158 μm emission lines and dust continua observed by the Atacama Large Millimeter/submillimeter Array (ALMA), taking advantage of the large program Reionization Era Bright Emission Line Survey. We find that the average [C ii] emission at z ∼ 7 has an effective radius r e of 2.2 ± 0.2 kpc. It is ≳2× larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for z ≲ 6 galaxies. Additionally, we compared the average [C ii] size with 4 < z < 6 galaxies observed by the ALMA Large Program to INvestigate [C ii] at Early times (ALPINE). By analyzing [C ii] sizes of 4 < z < 6 galaxies in two redshift bins, we find an average [C ii] size of r e = 2.2 ± 0.2 kpc and r e = 2.5 ± 0.2 kpc for z ∼ 5.5 and z ∼ 4.5 galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the [C ii] 158 μm emitting regions at redshift between z ∼ 7 and z ∼ 4. This finding suggests that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range.

The neutral atomic gas content of individual galaxies at large cosmological distances has until recently been difficult to measure due to the weakness of the hyperfine H i 21 cm transition. Here we estimate the H i gas mass of a sample of main-sequence star-forming galaxies at z ∼ 6.5–7.8 surveyed for [C ii] 158 μm emission as part of the Reionization Era Bright Emission Line Survey (REBELS), using a recent calibration of the [C ii]-to-H i conversion factor. We find that the H i gas mass excess in galaxies increases as a function of redshift, with an average of M Hi /M ⋆ ≈ 10, corresponding to H i gas mass fractions of f Hi = M Hi /(M ⋆ + M Hi ) = 90%, at z ≈ 7. Based on the [C ii] 158 μm luminosity function (LF) derived from the same sample of galaxies, we further place constraints on the cosmic H i gas mass density in galaxies (ρ Hi ) at this redshift, which we measure to be ρ H I = 7.1 − 3.0 + 6.4 × 10 6 M ⊙ Mpc − 3 . This estimate is substantially lower by a factor of ≈10 than that inferred from an extrapolation of damped Lyα absorber (DLA) measurements and largely depends on the exact [C ii] LF adopted. However, we find this decrease in ρ Hi to be consistent with recent simulations and argue that this apparent discrepancy is likely a consequence of the DLA sight lines predominantly probing the substantial fraction of H i gas in high-z galactic halos, whereas [C ii] traces the H i in the ISM associated with star formation. We make predictions for this buildup of neutral gas in galaxies as a function of redshift, showing that at z ≳ 5, only ≈10% of the cosmic H i gas content is confined in galaxies and associated with the star-forming ISM.

Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order to complete the census between LBGs and H-dropouts. Our criterion (H> 26.5 mag & [4.5] < 25 mag) combined with a de-blending technique is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies. In total, we identified 27 OFGs at z_phot > 3 (z_med=4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log($M_{\star}$/$M_{\odot}$) = 9.4-11.1. We find that up to 75% of the OFGs with log($M_{\star}$/$M_{\odot}$) = 9.5-10.5 were neglected by previous LBGs and H-dropout selection techniques. After performing stacking analyses, the OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than typical star-forming galaxies. Their SFR_tot (SFR_ IR+SFR_UV) is much larger than SFR_UVcorr (corrected for dust extinction), with SFR_tot/SFR_UVcorr = $8\pm1$, suggesting the presence of hidden dust regions in the OFGs that absorb all UV photons. The average dust size measured by a circular Gaussian model fit is R_e(1.13 mm)=1.01$\pm$0.05 kpc. We find that the cosmic SFRD at z>3 contributed by massive OFGs is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculate the combined contribution of OFGs and LBGs to the cosmic SFRD at z=4-5 to be 4 $\times$ 10$^{-2}$ $M_{\odot}$ yr$^{-1}$Mpc$^{-3}$, which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift., 26 pages, 14 figures, 8 tables. Accepted for publication in A&A

We present mid-infrared spectroscopic observations of the nucleus of the nearby Seyfert galaxy NGC 7469 taken with the MIRI instrument on the James Webb Space Telescope (JWST) as part of Directors Discretionary Time Early Release Science program 1328. The high-resolution nuclear spectrum contains 19 emission lines covering a wide range of ionization. The high-ionization lines show broad, blueshifted emission reaching velocities up to 1700 km s−1 and FWHM ranging from ∼500 to 1100 km s−1. The width of the broad emission and the broad-to-narrow line flux ratios correlate with ionization potential. The results suggest a decelerating, stratified, AGN-driven outflow emerging from the nucleus. The estimated mass outflow rate is 1–2 orders of magnitude larger than the current black hole accretion rate needed to power the AGN. Eight pure rotational H2 emission lines are detected with intrinsic widths ranging from FWHM ∼125 to 330 km s−1. We estimate a total mass of warm H2 gas of ∼1.2 × 107 M ⊙ in the central 100 pc. The PAH features are extremely weak in the nuclear spectrum, but a 6.2 μm PAH feature with an equivalent width of ∼0.07 μm and a flux of 2.7 × 10−17 W m−2 is detected. The spectrum is steeply rising in the mid-infrared, with a silicate strength of ∼0.02, significantly smaller than seen in most PG QSOs but comparable to other Seyfert 1s. These early MIRI mid-infrared IFU data highlight the power of JWST to probe the multiphase interstellar media surrounding actively accreting supermassive black holes.

We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard timescale $\tau_0= 30$ Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07-0.1\%$ and a UV-to-total star formation rate relation such that $log (\psi_{\rm UV}) = -0.05 ~[log (\psi)]^{2} + 0.86 ~log(\psi) -0.05$ (implying that 55-80\% of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9-10} M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at $z=7$. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94-1.1\%$) that are up to $18\times$ larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing $\tau_0$ to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor $ \sim 2$. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions., Comment: Accepted to MNRAS

ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature Td remains mostly unconstrained due to the few available FIR continuum data at redshift $z$ > 5. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [C $\scriptstyle \rm II$] 158-μm line and underlying dust continuum measurements, we derive Td in the continuum and [C $\scriptstyle \rm II$] detected $z$ ≈ 7 galaxies in the ALMA Large Project REBELS sample. We find 39 < Td < 58 K, and dust masses in the narrow range Md = (0.9-3.6) × 107 M⊙. These results allow us to extend for the first time the reported Td($z$) relation into the Epoch of Reionization. We produce a new physical model that explains the increasing Td($z$) trend with the decrease of gas depletion time, tdep = Mg/SFR, induced by the higher cosmological accretion rate at early times; this hypothesis yields Td ∝ (1 + $z$)0.4. The model also explains the observed Td scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, Td only depends on the gas column density (metallicity), $T_{\rm d} \propto N_{\rm H}^{1/6}$ (Td ∝ Z-1/6). REBELS galaxies are on average relatively transparent, with effective gas column densities around NH ≃ (0.03-1) × 1021 cm-2. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated Td ≫ 60 K, are significantly obscured, low-metallicity systems. In fact, Td is higher in metal-poor systems due to their smaller dust content, which for fixed LIR results in warmer temperatures.

We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $\beta$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158\mu$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have (28-90.5)% of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 < {\rm SFR}/ (M_\odot {\rm yr}^{-1}) < 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 M_\odot$ and $52 \pm 11$ K, respectively. In some galaxies dust is abundant (REBELS-14, $M'_d \approx 3.4 \times 10^7 M_\odot$), or hot (REBELS-18, $T'_d \approx 67$ K). The dust distribution is compact ($ 1 M_\odot$. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [CII] nicely follow the local $L_{\rm CII}-$SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of $0.11\, y_P^{-2}$ Gyr, where $y_P$ is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed ($\beta, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(\beta-\beta_{\rm int})$, where $\beta_{\rm int}$ is the intrinsic UV slope, exceeds $I_m^*\approx 1120$ for a MW curve. For these objects we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-$\beta$ relation., Comment: 22 pages, 7 figures, accepted by MNRAS. Comments welcome

Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived physical properties, such as star formation rates, gas fractions, depletion timescales, and dust temperatures for the galaxy sample built from the survey. There exists a subset of galaxies that exhibit starburst-like short depletion timescales, but they are located within the scatter of the so-called main sequence of SFGs. These are dubbed starbursts in the main sequence and display the most compact star formation and they are characterized by the shortest depletion timescales, lowest gas fractions, and highest dust temperatures of the galaxy sample, compared to typical SFGs at the same stellar mass and redshift. They are also very massive, accounting for $\sim 60\%$ of the most massive galaxies in the sample ($\log (M_{\rm{*}}/M_{\odot}) > 11.0$). We find trends between the areas of the ongoing star formation regions and the derived physical properties for the sample, unveiling the role of compact star formation as a physical driver of these properties. Starbursts in the main sequence appear to be the extreme cases of these trends. We discuss possible scenarios of galaxy evolution to explain the results drawn from our galaxy sample. Our findings suggest that the star formation rate is sustained in SFGs by gas and star formation compression, keeping them within the main sequence even when their gas fractions are low and they are presumably on the way to quiescence., Accepted by A&A. 26 pages, 13 figures

By compiling a comprehensive census of literature studies, we investigate the evolution of the Main Sequence (MS) of star-forming galaxies (SFGs) in the widest range of redshift ($0 < z < 6$) and stellar mass ($10^{8.5}-10^{11.5}$ $M_{\odot}$) ever probed. We convert all observations to a common calibration and find a remarkable consensus on the variation of the MS shape and normalization across cosmic time. The relation exhibits a curvature towards the high stellar masses at all redshifts. The best functional form is governed by two parameters: the evolution of the normalization and the turnover mass ($M_0(t)$), which both evolve as a power law of the Universe age. The turn-over mass determines the MS shape. It marginally evolves with time, making the MS slightly steeper towards $z\sim4-6$. At stellar masses below $M_0(t)$, SFGs have a constant specific SFR (sSFR), while above $M_0(t)$ the sSFR is suppressed. We find that the MS is dominated by central galaxies. This allows to turn $M_0(t)$ into the corresponding host halo mass. This evolves as the halo mass threshold between cold and hot accretion regimes, as predicted by the theory of accretion, where the central galaxy is fed or starved of cold gas supply, respectively. We, thus, argue that the progressive MS bending as a function of the Universe age is caused by the lower availability of cold gas in halos entering the hot accretion phase, in addition to black hole feedback. We also find qualitatively the same trend in the largest sample of star-forming galaxies provided by the IllustrisTNG simulation. Nevertheless, we still note large quantitative discrepancies with respect to observations, in particular at the high mass end. These can not be easily ascribed to biases or systematics in the observed SFRs and the derived MS., Accepted for publication by MNRAS, corrected table 2

Nuclear rings are excellent laboratories for studying intense star formation. We present results from a study of nuclear star-forming rings in five nearby normal galaxies from the Star Formation in Radio Survey (SFRS) and four local LIRGs from the Great Observatories All-sky LIRG Survey (GOALS) at sub-kpc resolutions using VLA high-frequency radio continuum observations. We find that nuclear ring star formation (NRSF) contributes 49 - 60\% of the total star formation of the LIRGs, compared to 7 - 40\% for the normal galaxies. We characterize a total of 58 individual star-forming regions in these rings, and find that with measured sizes of 10 - 200 pc, NRSF regions in the LIRGs have SFR and $\Sigma_\mathrm{SFR}$ up to 1.7 M$_\odot$yr$^{-1}$ and 402 M$_\odot$yr$^{-1}$kpc$^{-2}$, respectively, which are about 10 times higher than NRSF regions in the normal galaxies with similar sizes, and comparable to lensed high-$z$ star-forming regions. At $\sim 100 - 300$ pc scales, we estimate low contributions ($< 50\%$) of thermal free-free emission to total radio continuum emission at 33 GHz in the NRSF regions in the LIRGs, but large variations possibly exist at smaller physical scales. Finally, using archival sub-kpc resolution CO (J=1-0) data of nuclear rings in the normal galaxies and NGC 7469 (LIRG), we find a large scatter in gas depletion times at similar molecular gas surface densities, which tentatively points to a multi-modal star formation relation on sub-kpc scales., Comment: Accepted for publication in ApJ

We present the results of a Hubble Space Telescope WFC3 near-UV and ACS/WFC optical study into the star cluster populations of 10 luminous and ultra-luminous infrared galaxies (U/LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). Through integrated broadband photometry we have derived ages, masses, and extinctions for a total of 1027 star clusters in galaxies with $d_{L}, Comment: 70 pages, 21 figures, 5 tables, accepted for publication in ApJ

We present IRAM-30m Telescope $^{12}$CO and $^{13}$CO observations of a sample of 55 luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in the local universe. This sample is a subset of the Great Observatory All-Sky LIRG Survey (GOALS), for which we use ancillary multi-wavelength data to better understand their interstellar medium and star formation properties. Fifty-three (96%) of the galaxies are detected in $^{12}$CO, and 29 (52%) are also detected in $^{13}$CO above a 3$\sigma$ level. The median full width at zero intensity (FWZI) velocity of the CO line emission is 661km s$^{-1}$, and $\sim$54% of the galaxies show a multi-peak CO profile. Herschel photometric data is used to construct the far-IR spectral energy distribution of each galaxy, which are fit with a modified blackbody model that allows us to derive dust temperatures and masses, and infrared luminosities. We make the assumption that the gas-to-dust mass ratio of (U)LIRGs is comparable to local spiral galaxies with a similar stellar mass (i.e., gas/dust of mergers is comparable to their progenitors) to derive a CO-to-H$_2$ conversion factor of $\langle\alpha\rangle=1.8^{+1.3}_{-0.8}M_\odot$(K km s$^{-1}$pc$^{2}$)$^{-1}$; such a value is comparable to that derived for (U)LIRGs based on dynamical mass arguments. We derive gas depletion times of $400-600$Myr for the (U)LIRGs, compared to the 1.3Gyr for local spiral galaxies. Finally, we re-examine the relationship between the $^{12}$CO/$^{13}$CO ratio and dust temperature, confirming a transition to elevated ratios in warmer systems., Comment: 32 pages, 10 figures. Accepted for publication in Astronomy & Astrophysics

We present a search for [C II] line and dust continuum emission from optical dropout galaxies at z > 6 using ASPECS, our Atacama Large Millimeter submillimeter Array Spectroscopic Survey in the Hubble Ultra-deep Field (UDF). Our observations, which cover the frequency range of 212–272 GHz, encompass approximately the range of 6 4.5σ, two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60% of our candidates are expected to be spurious. For one of our blindly selected [C II] line candidates, we tentatively detect the CO(6-5) line in our parallel 3 mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all [C II] candidates results in a tentative detection with S_(1.2 mm) = 14 ± 5 μJy. This implies a dust-obscured star-formation rate (SFR) of (3 ± 1) M⊙ yr^(-1). We find that the two highest-SFR objects have candidate [C II] lines with luminosities that are consistent with the low-redshift L_([C II]) versus SFR relation. The other candidates have significantly higher [C II] luminosities than expected from their UV-based SFR. At the current sensitivity, it is unclear whether the majority of these sources are intrinsically bright [C II] emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for [C II] emitters at 6 < z < 8 that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.