We present the first predictions for the $L_{\rm [CII]}$ - SFR relation and [CII] luminosity function (LF) in the Epoch of Reionization (EoR) based on cosmological hydrodynamics simulations using the SIMBA suite plus radiative transfer calculations via S\'IGAME. The sample consists of 11,137 galaxies covering halo mass $\log M_{\rm halo}\in$[9, 12.4] $M_\odot$, star formation rate SFR$\in$[0.01, 330] $M_\odot$ yr$^{-1}$, and metallicity $<Z_{\rm gas}>_{\rm SFR}\in$[0.1, 0.9] $Z_\odot$. The simulated $L_{\rm [CII]}$-SFR relation is consistent with the range observed, but with a spread of $\simeq$0.3 dex at the high end of SFR ($>$100 $M_\odot$ yr$^{-1}$) and $\simeq$0.6 dex at the lower end, and there is tension between our predictions and the values of $L_{\rm [CII]}$ above 10$^{8.5}$ $L_\odot$ observed in some galaxies reported in the literature. The scatter in the $L_{\rm [CII]}$-SFR relation is mostly driven by galaxy properties, such that at a given SFR, galaxies with higher molecular gas mass and metallicity have higher $L_{\rm [CII]}$. The [CII] LF predicted by SIMBA is consistent with the upper limits placed by the only existing untargeted flux-limited [CII] survey at the EoR (ASPECS) and those predicted by semi-analytic models. We compare our results with existing models and discuss differences responsible for the discrepant slopes in the $L_{\rm [CII]}$-SFR relatiion.

We study the properties of molecular-forming gas clumps (MGCs) at the epoch of reionization using cosmological zoom-in simulations. We identify MGCs in a z=6 prototypical galaxy ("Althaea") using an H2 density-based clump finder. We compare their mass, size, velocity dispersion, gas surface density, and virial parameter (alpha_vir) to observations. In Althaea, the typical MGC mass and size are Mgas=10^6.5 Msun and R=45-100 pc, which are comparable to those found in nearby spirals and starburst galaxies. MGCs are highly supersonic and supported by turbulence, with rms velocity dispersions of sigma_gas=20-100 km s^-1 and pressure of P/k_B=10^7.6 K cm^-3 (i.e., >1000x with respect to the Milky Way), similar to those found in nearby and z~2 gas-rich starburst galaxies. In addition, we perform stability analysis to understand the origin and dynamical properties of MGCs. We find that MGCs are globally stable in the main disk of Althaea. Densest regions where star formation is expected to take place in clouds and cores on even smaller scales instead have lower alpha_vir and Toomre-Q values. Detailed studies of the star-forming gas dynamics at the epoch of reionization thus require a spatial resolution of <40 pc (=0.01"), which is within reach with the Atacama Large (sub-)Millimeter Array and the Next Generation Very Large Array.

We report the detection of the far-infrared (FIR) fine-structure line of singly ionised nitrogen, \Nplusa, within the peak epoch of galaxy assembly, from a strongly lensed galaxy, hereafter ``The Red Radio Ring''; the RRR, at z = 2.55. We combine new observations of the ground-state and mid-J transitions of CO (J$_{\rm up} =$ 1,5,8), and the FIR spectral energy distribution (SED), to explore the multi-phase interstellar medium (ISM) properties of the RRR. All line profiles suggest that the HII regions, traced by \Nplusa, and the (diffuse and dense) molecular gas, traced by the CO, are co-spatial when averaged over kpc-sized regions. Using its mid-IR-to-millimetre (mm) SED, we derive a non-negligible dust attenuation of the \Nplusa line emission. Assuming a uniform dust screen approximation results a mean molecular gas column density $> 10^{24}$\, cm$^{-2}$, with a molecular gas-to-dust mass ratio of 100. It is clear that dust attenuation corrections should be accounted for when studying FIR fine-structure lines in such systems. The attenuation corrected ratio of $L_{\rm NII205} / L_{\rm IR(8-1000\mu m)} = 2.7 \times 10^{-4}$ is consistent with the dispersion of local and $z >$ 4 SFGs. We find that the lower-limit, \Nplusa -based star-formation rate (SFR) is less than the IR-derived SFR by a factor of four. Finally, the dust SED, CO line SED and $L_{\rm NII205}$ line-to-IR luminosity ratio of the RRR is consistent with a starburst-powered ISM.

ALMA 870$\mu$m continuum imaging has uncovered a population of blends of multiple dusty star-forming galaxies (DSFGs) in sources originally detected with the Herschel Space Observatory. However, their pairwise separations are much smaller that what is found by ALMA follow-up of other single-dish surveys or expected from theoretical simulations. Using ALMA and VLA, we have targeted three of these systems to confirm whether the multiple 870$\mu$m continuum sources lie at the same redshift, successfully detecting $^{12}$CO($J = 3$-2) and $^{12}$CO($J = 1$-0) lines and being able to confirm that in the three cases all the multiple DSFGs are likely physically associated within the same structure. Therefore, we report the discovery of two new gas-rich dusty protocluster cores (HELAISS02, $z = 2.171 \pm 0.004$; HXMM20, $z = 2.602 \pm 0.002$). The third target is located in the well known COSMOS overdensity at $z = 2.51$ (named CL J1001+0220 in the literature), for which we do not find any new secure CO(1-0) detection, although some of its members show only tentative detections and require further confirmation. From the gas, dust, and stellar properties of the two new protocluster cores, we find very large molecular gas fractions yet low stellar masses, pushing the sources above the main sequence, while not enhancing their star formation efficiency. We suggest that the sources might be newly formed galaxies migrating to the main sequence. The properties of the three systems compared to each other and to field galaxies may suggest a different evolutionary stage between systems.

We present far-IR photometry and infrared spectrum of the z=3.9114 quasar/starburst composite system APM 08279+5255 obtained using the Stratospheric Observatory for Infrared Astronomy (SOFIA)/HAWC+ and the Spitzer Space Telescope Infrared Spectrograph (IRS). We decompose the IR-to-radio spectral energy distribution (SED), sampled in 51 bands, using (i) a model comprised of two-temperature modified blackbodies (MBB) and radio power-laws and (ii) a semi-analytic model, which also accounts for emission from a clumpy torus. The latter is more realistic but requires a well-sampled SED, as possible here. In the former model, we find temperatures of T_warm = 296^17_15 K and T_cold = 110^3_3 K for the warm and cold dust components, respectively. This model suggests that the cold dust component dominates the FIR energy budget (66%) but contributes only 17% to the total IR luminosity. Based on the torus models, we infer an inclination angle of i=15^8_8 degree and the presence of silicate emission, in accordance with the Type-1 active galactic nucleus nature of APM 08279+5255. Accounting for the torus' contribution to the FIR luminosity, we find a lensing-corrected star formation rate of SFR=3075x(4/mu_L) Msun yr^-1. We find that the central quasar contributes 30% to the FIR luminosity but dominates the total IR luminosity (93%). The 30% correction is in contrast to the 90% reported in previous work. In addition, the IR luminosity inferred from the torus model is a factor of two higher. These differences highlight the importance of adopting physically motivated models to properly account for IR emission in high-z quasars, which is now possible with SOFIA/HAWC+.

We present CO(J= 1-0; 3-2; 5-4; 10-9) and 1.2-kpc resolution [CII] line observations of the dusty star-forming galaxy (SFG) HXMM05 -- carried out with the Karl G. Jansky Very Large Array, the Combined Array for Research in Millimeter-wave Astronomy, the Plateau de Bure Interferometer, and the Atacama Large Millimeter/submillimeter Array, measuring an unambiguous redshift of z = 2.9850+/-0.0009. We find that HXMM05 is a hyper-luminous infrared galaxy (LIR=(4+/-1)x10^13 Lsun) with a total molecular gas mass of (2.1+/-0.7)x10^11 (alpha_CO/0.8) Msun. The CO(J=1-0) and [CII] emission are extended over ~9 kpc in diameter, and the CO line FWHM exceeds 1100 km s^-1. The [CII] emission shows a monotonic velocity gradient consistent with a disk, with a maximum rotation velocity of v_c = 616+/-100 km s^-1 and a dynamical mass of (7.7+/-3.1)x10^11 Msun. We find a star formation rate (SFR) of 2900^750_-595 Msun yr^-1. HXMM05 is thus among the most intensely star-forming galaxies known at high redshift. Photo-dissociation region modeling suggests physical conditions similar to nearby SFGs, showing extended star formation, which is consistent with our finding that the gas and dust emission are co-spatial. Its molecular gas excitation resembles the local major merger Arp 220. The broad CO and [CII] lines and a pair of compact dust nuclei suggest the presence of a late-stage major merger at the center of the extended disk, again reminiscent of Arp 220. The observed gas kinematics and conditions together with the presence of a companion and the pair of nuclei suggest that HXMM05 is experiencing multiple mergers as a part of the evolution.

We report the detection of ADFS-27, a dusty, starbursting major merger at a redshift of z=5.655, using the Atacama Large Millimeter/submillimeter Array (ALMA). ADFS-27 was selected from Herschel/SPIRE and APEX/LABOCA data as an extremely red "870 micron riser" (i.e., S_250<S_350<S_500<S_870), demonstrating the utility of this technique to identify some of the highest-redshift dusty galaxies. A scan of the 3mm atmospheric window with ALMA yields detections of CO(5-4) and CO(6-5) emission, and a tentative detection of H2O(211-202) emission, which provides an unambiguous redshift measurement. The strength of the CO lines implies a large molecular gas reservoir with a mass of M_gas=2.5x10^11(alpha_CO/0.8)(0.39/r_51) Msun, sufficient to maintain its ~2400 Msun/yr starburst for at least ~100 Myr. The 870 micron dust continuum emission is resolved into two components, 1.8 and 2.1 kpc in diameter, separated by 9.0 kpc, with comparable dust luminosities, suggesting an ongoing major merger. The infrared luminosity of L_IR~=2.4x10^13Lsun implies that this system represents a binary hyper-luminous infrared galaxy, the most distant of its kind presently known. This also implies star formation rate surface densities of Sigma_SFR=730 and 750Msun/yr/kpc2, consistent with a binary "maximum starburst". The discovery of this rare system is consistent with a significantly higher space density than previously thought for the most luminous dusty starbursts within the first billion years of cosmic time, easing tensions regarding the space densities of z~6 quasars and massive quiescent galaxies at z>~3.

We present the results of combined deep Keck/NIRC2, HST/WFC3 near-infrared and Herschel far infrared observations of an extremely star forming dusty lensed galaxy identified from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS J133542.9+300401). The galaxy is gravitationally lensed by a massive WISE identified galaxy cluster at $z\sim1$. The lensed galaxy is spectroscopically confirmed at $z=2.685$ from detection of $\rm {CO (1 \rightarrow 0)}$ by GBT and from detection of $\rm {CO (3 \rightarrow 2)}$ obtained with CARMA. We use the combined spectroscopic and imaging observations to construct a detailed lens model of the background dusty star-forming galaxy (DSFG) which allows us to study the source plane properties of the target. The best-fit lens model provide magnification of $\mu_{\rm star}=2.10\pm0.11$ and $\mu_{\rm dust}=2.02\pm0.06$ for the stellar and dust components respectively. Multi-band data yields a magnification corrected star formation rate of $1900(\pm200)\,M_{\odot}{\rm yr^{-1}}$ and stellar mass of $6.8_{-2.7}^{+0.9}\times10^{11}\,M_{\odot}$ consistent with a main sequence of star formation at $z\sim2.6$. The CO observations yield a molecular gas mass of $8.3(\pm1.0)\times10^{10}\,M_{\odot}$, similar to the most massive star-forming galaxies, which together with the high star-formation efficiency are responsible for the intense observed star formation rates. The lensed DSFG has a very short gas depletion time scale of $\sim40$ Myr. The high stellar mass and small gas fractions observed indicate that the lensed DSFG likely has already formed most of its stellar mass and could be a progenitor of the most massive elliptical galaxies found in the local Universe.

The experimental efforts to detect the redshifted 21 cm signal from the Epoch of Reionization (EoR) are limited predominantly by the chromatic instrumental systematic effect. The delay spectrum methodology for 21 cm power spectrum measurements brought new attention to the critical impact of an antenna's chromaticity on the viability of making this measurement. This methodology established a straightforward relationship between time-domain response of an instrument and the power spectrum modes accessible to a 21 cm EoR experiment. We examine the performance of a prototype of the Hydrogen Epoch of Reionization Array (HERA) array element that is currently observing in Karoo desert, South Africa. We present a mathematical framework to derive the beam integrated frequency response of a HERA prototype element in reception from the return loss measurements between 100-200 MHz and determined the extent of additional foreground contamination in the delay space. The measurement reveals excess spectral structures in comparison to the simulation studies of the HERA element. Combined with the HERA data analysis pipeline that incorporates inverse covariance weighting in optimal quadratic estimation of power spectrum, we find that in spite of its departure from the simulated response, HERA prototype element satisfies the necessary criteria posed by the foreground attenuation limits and potentially can measure the power spectrum at spatial modes as low as $k_{\parallel} > 0.1h$~Mpc$^{-1}$. The work highlights a straightforward method for directly measuring an instrument response and assessing its impact on 21 cm EoR power spectrum measurements for future experiments that will use reflector-type antenna.

We report observations of CO(J=2-1) and CO(J=3-2) line emission towards the quadruply-lensed quasar RXS J1131-1231 at z = 0.654 obtained using the Plateau de Bure Interferometer (PdBI) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA). Our lens modeling shows that the asymmetry in the double-horned CO(J = 2-1) line profile is mainly a result of differential lensing, where the magnification factor varies from ~3 to ~9 across different kinematic components. The intrinsically symmetric line profile and a smooth source-plane velocity gradient suggest that the host galaxy is an extended rotating disk, with a CO size of R_CO~6 kpc and a dynamical mass of M_dyn~8x10^10 Msun. We also find a secondary CO-emitting source near RXS J1131-1231 whose location is consistent with the optically-faint companion reported in previous studies. The lensing-corrected molecular gas masses are M_gas = (1.4+/-0.3)x10^10 Msun and (2.0+/-0.1)x10^9 Msun for RXS J1131-1231 and the companion, respectively. We find a lensing-corrected stellar mass of M_* = (3+/-1)x10^10 Msun and a star formation rate of SFR_FIR = (120+/-63) Msun yr^-1 , corresponding to a specific SFR and star formation efficiency comparable to z~1 disk galaxies not hosting quasars. The implied gas mass fraction of ~18+/-4% is consistent with the previously-observed cosmic decline since z~2. We thus find no evidence for quenching of star formation in RXS J1131-1231. This agrees with our finding of an elevated MBH/Mbulge ratio of >0.27+0.11% -0.08 compared to the local value, suggesting that the bulk of its black hole mass is largely in place while its stellar bulge is still assembling.

We report the detection of CO(J=3-2) line emission in the strongly-lensed submillimeter galaxy (SMG) SMM J0939+8315 at z=2.221, using the Combined Array for Research in Millimeter-wave Astronomy. SMM J0939+8315 hosts a type-2 quasar, and is gravitationally lensed by the radio galaxy 3C220.3 and its companion galaxy at z=0.685. The 104 GHz continuum emission underlying the CO line is detected toward 3C220.3 with an integrated flux density of S_cont = 7.4 +/- 1.4 mJy. Using the CO(J=3-2) line intensity of I_(CO(3-2)) = (12.6 +/- 2.0) Jy km s^-1, we derive a lensing- and excitation-corrected CO line luminosity of L'(CO(3-2)) = (3.4 +/- 0.7) x 10^10 (10.1/mu_L) K km s^-1 pc^2 for the SMG, where mu_L is the lensing magnification factor inferred from our lens modeling. This translates to a molecular gas mass of M_gas = (2.7 +/- 0.6) x 10^10 (10.1/mu_L) Msun. Fitting spectral energy distribution models to the (sub)-millimeter data of this SMG yields a dust temperature of T = 63.1^{+1.1}_{-1.3} K, a dust mass of M_dust = (5.2 +/- 2.1) x 10^8 (10.1/mu_L) Msun, and a total infrared luminosity of L_IR = (9.1 +/- 1.2) x 10^12 (10.1/mu_L) Lsun. We find that the properties of the interstellar medium of SMM J0939+8315 overlap with both SMGs and type-2 quasars. Hence, SMM J0939+8315 may be transitioning from a star-bursting phase to an unobscured quasar phase as described by the "evolutionary link" model, according to which this system may represent an intermediate stage in the evolution of present-day galaxies at an earlier epoch.