Search

Using integrated spectra for two gravitationally lensed galaxies from the JWST TEMPLATES Early Release Science program, we analyze faint auroral lines, which provide direct measurements of the gas-phase chemical abundance. For the brighter galaxy, SGAS1723$+$34 ($z = 1.3293$), we detect the [OIII]$\lambda4363$, [SIII]$\lambda6312$, and [OII]$\lambda\lambda$7320,7330 auroral emission lines, and set an upper limit for the [NII]$\lambda5755$ line. For the second galaxy, SGAS1226$+$21 ($z = 2.925$), we do not detect any auroral lines, and report upper limits. With these measurements and upper limits, we constrain the electron temperatures in different ionization zones within both of these galaxies. For SGAS1723$+$34, where auroral lines are detected, we calculate direct oxygen and nitrogen abundances, finding an N/O ratio consistent with observations of nearby ($z\sim 0$) galaxies. These observations highlight the potent combination of JWST and gravitational lensing to measure faint emission lines in individual distant galaxies and to directly study the chemical abundance patterns in those galaxies., Comment: Submitted to the Astrophysical Journal's Focus Issue on the TEMPLATES JWST Early Release Science Program. 15 pages, 6 figures. Comments welcome

We describe a custom outlier rejection algorithm for JWST/NIRSpec integral field spectroscopy. This method uses a layered sigma clipping approach that adapts clipping thresholds based upon the spatial profile of the science target. We find that this algorithm produces a robust outlier rejection while simultaneously preserving the signal of the science target. Originally developed as a response to unsatisfactory initial performance of the jwst pipeline outlier detection step, this method works either as a standalone solution, or as a supplement to the current pipeline software. Comparing leftover (i.e., not flagged) artifacts with the current pipeline's outlier detection step, we find that our method results in one fifth as many residual artifacts as the jwst pipeline. However, we find a combination of both methods removes nearly all artifacts -- an approach that takes advantage of both our algorithm's robust outlier rejection and the pipeline's use of individual dithers. This combined approach is what the TEMPLATES Early Release Science team has converged upon for our NIRSpec observations. Finally, we publicly release the code and Jupyter notebooks for the custom outlier rejection algorithm., Comment: 10 pages, including 5 figures. Submitted to the Publications of the Astronomical Society of the Pacific (PASP). Code associated with paper released at https://github.com/aibhleog/baryon-sweep

This paper gives an overview of TEMPLATES, a JWST Early Release Science program that targeted four extremely bright, gravitationally lensed galaxies: two extremely dusty, two with low attenuation, as templates for galaxy evolution studies with JWST. TEMPLATES obtains a common set of spectral diagnostics for these 1.3 < z < 4.2 galaxies, in particular H alpha, Paschen alpha, and the rest-frame optical and near-infrared continua. In addition, two of the four targets have JWST coverage of [O III] 5007 Angstrom and H beta; the other two targets have have JWST coverage of PAH 3.3 micron and complementary ALMA data covering the [C II] 158 micron emission line. The science goals of TEMPLATES are to demonstrate attenuation-robust diagnostics of star formation, map the distribution of star formation, compare the young and old stellar populations, and measure the physical conditions of star formation and their spatial variation across the galaxies. In addition, TEMPLATES has technical goals to establish best practices for the Integral Field Units (IFU) within the NIRSpec and MIRI instruments, both in terms of observing strategy and in terms of data reduction. The paper describes TEMPLATES's observing program, scientific and technical goals, data reduction methods, and deliverables, including high-level data products and data reduction cookbooks., Comment: 28 pages, including 8 figures and 3 tables. Submitted to The Astrophysical Journal's Focus Issue on the TEMPLATES JWST Early Release Science Program, https://iopscience.iop.org/collections/apj-231204-01-templates

We present JWST and ALMA results for the lensing system SPT0418-47, which includes a strongly-lensed, dusty star-forming galaxy at redshift z=4.225 and an associated multiply-imaged companion. JWST NIRCam and MIRI imaging observations presented in this paper were acquired as part of the Early Release Science program Targeting Extremely Magnified Panchromatic Lensed Arcs and Their Extended Star Formation (TEMPLATES). This data set provides robust, mutiwavelength detection of stellar light in both the main (SPT0418A) and companion (SPT0418B) galaxies, while the ALMA detection of [C II] emission confirms that SPT0418B lies at the same redshift as SPT0418A. From a source plane reconstruction, we infer that the projected physical separation of the two galaxies is $4.42\pm 0.05$ kpc. We derive total magnifications of $\mu=29.5\pm1.2$ and $\mu=4.2\pm 0.9$ for SPT0418A and SPT0418B, respectively. We use both CIGALE and PROSPECTOR to derive stellar masses. The stellar mass ratio of SPT0418A and SPT0418B is approximately 4 to 1 ($4.5\pm 1.0$ for CIGALE and $4.2^{+1.9}_{-1.6}$ for PROSPECTOR). We also see evidence of extended structure associated with SPT0418A in the lensing reconstruction that is suggestive of a tidal feature. Interestingly, the star formation rates and stellar masses of both galaxies are consistent with the main sequence of star-forming galaxies at this epoch, indicating that this ongoing interaction has not noticeably elevated the star formation levels., Comment: Submitted to AAS Journals,14 pages, 7 figures

We present the first spatially resolved maps of gas-phase metallicity for two dust-obscured star-forming galaxies (DSFGs) at $z\sim$ 4, from the JWST TEMPLATES Early Release Science program, derived from NIRSpec integral field unit spectroscopy of the H$\alpha$ and [NII] emission lines. Empirical optical line calibrations are used to determine that the sources are globally enriched to near-solar levels. While one source shows elevated [NII]/H$\alpha$ ratios and broad H$\alpha$ emission consistent with the presence of an AGN in a $\gtrsim$1kpc region, we argue that both systems have already undergone significant metal enrichment as a result of their extremely high star formation rates. Utilizing ALMA rest-frame 380$\mu$m continuum and [CI]($^3$P$_2$-$^3$P$_1$) line maps we compare the spatial variation of the metallicity and gas-to-dust ratio in the two galaxies, finding the two properties to be anticorrelated on highly resolved spatial scales, consistent with various literature studies of $z\sim$ 0 galaxies. The data are indicative of the enormous potential of JWST to probe the enrichment of the interstellar medium on $\sim$kpc scales in extremely dust-obscured systems at $z\sim$ 4 and beyond., Comment: 14 pages, 6 figures, accepted for publication in ApJ

Dust grains absorb half of the radiation emitted by stars throughout the history of the universe, re-emitting this energy at infrared wavelengths. Polycyclic aromatic hydrocarbons (PAHs) are large organic molecules that trace millimeter-size dust grains and regulate the cooling of the interstellar gas within galaxies. Observations of PAH features in very distant galaxies have been difficult due to the limited sensitivity and wavelength coverage of previous infrared telescopes. Here we present JWST observations that detect the 3.3um PAH feature in a galaxy observed less than 1.5 billion years after the Big Bang. The high equivalent width of the PAH feature indicates that star formation, rather than black hole accretion, dominates the infrared emission throughout the galaxy. The light from PAH molecules, large dust grains, and stars and hot dust are spatially distinct from one another, leading to order-of-magnitude variations in the PAH equivalent width and the ratio of PAH to total infrared luminosity across the galaxy. The spatial variations we observe suggest either a physical offset between the PAHs and large dust grains or wide variations in the local ultraviolet radiation field. Our observations demonstrate that differences in the emission from PAH molecules and large dust grains are a complex result of localized processes within early galaxies., Comment: Published in Nature 5 June 2023 at https://www.nature.com/articles/s41586-023-05998-6. MIRI MRS reduction notebook is available at https://github.com/jwst-templates

We have observed the z=4.3 protocluster SPT2349-56 with ATCA with the aim of detecting radio-loud active galactic nuclei (AGN) amongst the ~30 submillimeter galaxies identified in the structure. We detect the central complex of SMGs at 2.2\,GHz with a luminosity of L_2.2=(4.42pm0.56)x10^{25} W/Hz. The ASKAP also detects the source at 888 MHz, constraining the radio spectral index to alpha=-1.6pm0.3, consistent with ATCA non-detections at 5.5 and 9GHz, and implying L_1.4(rest)=(2.4pm0.3)x10^{26}W/Hz. This radio luminosity is about 100 times higher than expected from star formation, assuming the usual FIR-radio correlation, which is a clear indication of an AGN driven by a forming brightest cluster galaxy (BCG). None of the SMGs in SPT2349-56 show signs of AGN in any other diagnostics available to us (notably 12CO out to J=16, OH163um, CII/IR, and optical spectra), highlighting the radio continuum as a powerful probe of obscured AGN in high-z protoclusters. No other significant radio detections are found amongst the cluster members, consistent with the FIR-radio correlation. We compare these results to field samples of radio sources and SMGs, along with the 22 SPT-SMG gravitational lenses also observed in the ATCA program, as well as powerful radio galaxies at high redshifts. Our results allow us to better understand the effects of this gas-rich, overdense environment on early supermassive black hole (SMBH) growth and cluster feedback. We estimate that (3.3pm0.7)x10^{38} W of power are injected into the growing ICM by the radio-loud AGN, whose energy over 100Myr is comparable to the binding energy of the gas mass of the central halo. The AGN power is also comparable to the instantaneous energy injection from supernova feedback from the 23 catalogued SMGs in the core region of 120kpc projected radius. The SPT2349-56 radio-loud AGN may be providing strong feedback on a nascent ICM., Comment: 31 pages, submitted to ApJ, Dec17,2022

The SPT0311-58 system at z=6.900 is an extremely massive structure within the reionization epoch, and offers a chance to understand the formation of galaxies in an extreme peak in the primordial density field. We present 70mas Atacama Large Millimeter/submillimeter Array observations of the dust continuum and CII 158um emission in the central pair of galaxies and reach physical resolution ~100-350pc, among the most detailed views of any reionization-era system to date. The observations resolve the source into at least a dozen kiloparsec-size clumps. The global kinematics and high turbulent velocity dispersion within the galaxies present a striking contrast to recent claims of dynamically cold thin-disk kinematics in some dusty galaxies just 800Myr later at z~4. We speculate that both gravitational interactions and fragmentation from massive parent disks have likely played a role in the overall dynamics and formation of clumps in the system. Each clump individually is comparable in mass to other 6

We present a sample of 706, $z < 1.5$ active galactic nuclei (AGNs) selected from optical photometric variability in three of the Dark Energy Survey (DES) deep fields (E2, C3, and X3) over an area of 4.64 deg$^2$. We construct light curves using difference imaging aperture photometry for resolved sources and non-difference imaging PSF photometry for unresolved sources, respectively, and characterize the variability significance. Our DES light curves have a mean cadence of 7 days, a 6 year baseline, and a single-epoch imaging depth of up to $g \sim 24.5$. Using spectral energy distribution (SED) fitting, we find 26 out of total 706 variable galaxies are consistent with dwarf galaxies with a reliable stellar mass estimate ($M_{\ast}<10^{9.5}\ M_\odot$; median photometric redshift of 0.9). We were able to constrain rapid characteristic variability timescales ($\sim$ weeks) using the DES light curves in 15 dwarf AGN candidates (a subset of our variable AGN candidates) at a median photometric redshift of 0.4. This rapid variability is consistent with their low black hole masses. We confirm the low-mass AGN nature of one source with a high S/N optical spectrum. We publish our catalog, optical light curves, and supplementary data, such as X-ray properties and optical spectra, when available. We measure a variable AGN fraction versus stellar mass and compare to results from a forward model. This work demonstrates the feasibility of optical variability to identify AGNs with lower black hole masses in deep fields, which may be more "pristine" analogs of supermassive black hole seeds., Comment: 21 pages, 22 figures incl. 3 appendices; accepted for publication in MNRAS

We present Northern Extended Millimeter Array (NOEMA) CO(2-1) maps of the z=0.4418 cluster-central QSO IRAS 09104+4109, which trace ~4.5x10^10 Msol of molecular gas in and around the galaxy. As in many low redshift cool core clusters, the molecular gas is located in a series of clumps extending along the old radio jets and lobes. It has a relatively low velocity dispersion (336 [+39,-35] km/s FWHM) and shows no velocity gradients indicative of outflow or infall. Roughly half the gas is located in a central clump on the northeast side of the galaxy, overlapping a bright ionized gas filament and a spur of excess X-ray emission, suggesting that this is a location of rapid cooling. The molecular gas is unusually extended, out to ~55 kpc radius, comparable to the scale of the filamentary nebula in the Perseus cluster, despite the much higher redshift of this system. The extent falls within the thermal instability radius of the intracluster medium (ICM), with t_cool/t_ff<25 and t_cool}/t_eddy~1 within ~70 kpc. Continuum measurements at 159.9 GHz from NOEMA and 850 micron from the JCMT SCUBA-2 show excess far infrared emission, which we interpret as free-free emission arising from the ongoing starburst. These observations suggest that ICM cooling is not strongly affected by the buried QSO, and that cooling from the ICM can build gas reservoirs sufficient to fuel quasar-mode activity and drive the reorientation of the central AGN., Comment: 16 pages, 13 figures, 2 tables, accepted for publication in MNRAS, v2 updates references and corrects minor typographical errors identified at the proof stage

The protocluster SPT2349$-$56 at $z\,{=}\,4.3$ contains one of the most actively star-forming cores known, yet constraints on the total stellar mass of this system are highly uncertain. We have therefore carried out deep optical and infrared observations of this system, probing rest-frame ultraviolet to infrared wavelengths. Using the positions of the spectroscopically-confirmed protocluster members, we identify counterparts and perform detailed source deblending, allowing us to fit spectral energy distributions in order to estimate stellar masses. We show that the galaxies in SPT2349$-$56 have stellar masses proportional to their high star-formation rates, consistent with other protocluster galaxies and field submillimetre galaxies (SMGs) around redshift 4. The galaxies in SPT2349$-$56 have on average lower molecular gas-to-stellar mass fractions and depletion timescales than field SMGs, although with considerable scatter. We construct the stellar-mass function for SPT2349$-$56 and compare it to the stellar-mass function of $z\,{=}\,1$ galaxy clusters, finding consistent shapes between the two. We measure rest-frame galaxy ultraviolet half-light radii from our HST-F160W imaging, finding that on average the galaxies in our sample are similar in size to typical star-forming galaxies at these redshifts. However, the brightest HST-detected galaxy in our sample, found near the luminosity-weighted centre of the protocluster core, remains unresolved at this wavelength. Hydrodynamical simulations predict that the core galaxies will quickly merge into a brightest cluster galaxy, thus our observations provide a direct view of the early formation mechanisms of this class of object., Comment: Submitted to MNRAS Sept. 3, 2021 Accepted by MNRAS Dec, 4, 2021

SPT0311-58 is the most massive infrared luminous system discovered so far during the Epoch of Reionization (EoR). In this paper, we present a detailed analysis of the molecular interstellar medium at z = 6.9, through high-resolution observations of the CO(6-5), CO(7-6), CO(10-9), [CI](2-1), and p-H2O(211-202) lines and dust continuum emission with the Atacama Large Millimeter/submillimeter Array (ALMA). The system consists of a pair of intensely star-forming gravitationally lensed galaxies (labelled West and East). The intrinsic far-infrared luminosity is (16 $\pm$ 4)$\times\rm 10^{12} \ \rm L_{\odot}$ in West and (27 $\pm$ 4)$\times\rm 10^{11} \ \rm L_{\odot}$ in East. We model the dust, CO, and [CI] using non-local thermodynamic equilibrium radiative transfer models and estimate the intrinsic gas mass to be (5.4 $\pm$ 3.4)$\times\rm 10^{11} \ \rm M_{\odot}$ in West and (3.1 $\pm$ 2.7)$\times\rm 10^{10} \ \rm M_{\odot}$ in East. We find that the CO spectral line energy distribution in West and East are typical of high-redshift sub-millimeter galaxies (SMGs). The CO-to-H2 conversion factor ($\alpha_{CO}$) and the gas depletion time scales estimated from the model are consistent with the high-redshift SMGs in the literature within the uncertainties. We find no evidence of evolution of depletion time with redshift in SMGs at z > 3. This is the most detailed study of molecular gas content of a galaxy in the EoR to-date, with the most distant detection of H2O in a galaxy without any evidence for active galactic nuclei in the literature., Comment: Accepted for publication in The Astrophysical Journal. 34 pages, 17 figures

Galactic outflows of molecular gas are a common occurrence in galaxies and may represent a mechanism by which galaxies self-regulate their growth, redistributing gas that could otherwise have formed stars. We previously presented the first survey of molecular outflows at z > 4 towards a sample of massive, dusty galaxies. Here we characterize the physical properties of the molecular outflows discovered in our survey. Using low-redshift outflows as a training set, we find agreement at the factor-of-two level between several outflow rate estimates. We find molecular outflow rates 150-800Msun/yr and infer mass loading factors just below unity. Among the high-redshift sources, the molecular mass loading factor shows no strong correlations with any other measured quantity. The outflow energetics are consistent with expectations for momentum-driven winds with star formation as the driving source, with no need for energy-conserving phases. There is no evidence for AGN activity in our sample, and while we cannot rule out deeply-buried AGN, their presence is not required to explain the outflow energetics, in contrast to nearby obscured galaxies with fast outflows. The fraction of the outflowing gas that will escape into the circumgalactic medium (CGM), though highly uncertain, may be as high as 50%. This nevertheless constitutes only a small fraction of the total cool CGM mass based on a comparison to z~2-3 quasar absorption line studies, but could represent >~10% of the CGM metal mass. Our survey offers the first statistical characterization of molecular outflow properties in the very early universe., Comment: ApJ accepted. 25 pages, 16 figures. Data and tables from Papers I and II available at https://github.com/spt-smg/publicdata

Massive galaxy-scale outflows of gas are one of the most commonly-invoked mechanisms to regulate the growth and evolution of galaxies throughout the universe. While the gas in outflows spans a large range of temperatures and densities, the cold molecular phase is of particular interest because molecular outflows may be capable of suppressing star formation in galaxies by removing the star-forming gas. We have conducted the first survey of molecular outflows at z > 4, targeting 11 strongly-lensed dusty, star-forming galaxies (DSFGs) with high-resolution Atacama Large Millimeter Array (ALMA) observations of OH 119um absorption as an outflow tracer. In this first paper, we give an overview of the survey, focusing on the detection rate and structure of molecular outflows. We find unambiguous evidence for outflows in 8/11 (73%) galaxies, more than tripling the number known at z > 4. This implies that molecular winds in z > 4 DSFGs must have both a near-unity occurrence rate and large opening angles to be detectable in absorption. Lensing reconstructions reveal that 500pc-scale clumpy structures in the outflows are common. The individual clumps are not directly resolved, but from optical depth arguments we expect that future observations will require 50-200pc spatial resolution to do so. We do not detect high-velocity [CII] wings in any of the sources with clear OH outflows, indicating that [CII] is not a reliable tracer of molecular outflows. Our results represent a first step toward characterizing molecular outflows at z > 4 at the population level, demonstrating that large-scale outflows are ubiquitous among early massive, dusty galaxies., Comment: ApJ accepted. 28 pages, 12 figures + appendix. Data and tables from Papers I and II available at https://github.com/spt-smg/publicdata

We present APEX-LABOCA 870 micron observations of the fields surrounding the nine brightest, high-redshift, unlensed objects discovered in the South Pole Telescope's (SPT) 2500 square degrees survey. Initially seen as point sources by SPT's 1-arcmin beam, the 19-arcsec resolution of our new data enables us to deblend these objects and search for submillimetre (submm) sources in the surrounding fields. We find a total of 98 sources above a threshold of 3.7 sigma in the observed area of 1300 square arcminutes, where the bright central cores resolve into multiple components. After applying a radial cut to our LABOCA sources to achieve uniform sensitivity and angular size across each of the nine fields, we compute the cumulative and differential number counts and compare them to estimates of the background, finding a significant overdensity of approximately 10 at 14 mJy. The large overdensities of bright submm sources surrounding these fields suggest that they could be candidate protoclusters undergoing massive star-formation events. Photometric and spectroscopic redshifts of the unlensed central objects range from 3 to 7, implying a volume density of star-forming protoclusters of approximately 0.1 per giga-parsec cube. If the surrounding submm sources in these fields are at the same redshifts as the central objects, then the total star-formation rates of these candidate protoclusters reach 10,000 solar masses per year, making them much more active at these redshifts than what has been seen so far in both simulations and observations.

We report the identification of a low-mass AGN, DES J0218$-$0430, in a redshift $z = 0.823$ galaxy in the Dark Energy Survey (DES) Supernova field. We select DES J0218$-$0430 as an AGN candidate by characterizing its long-term optical variability alone based on DES optical broad-band light curves spanning over 6 years. An archival optical spectrum from the fourth phase of the Sloan Digital Sky Survey shows both broad Mg II and broad H$\beta$ lines, confirming its nature as a broad-line AGN. Archival XMM-Newton X-ray observations suggest an intrinsic hard X-ray luminosity of $L_{{\rm 2-12\,keV}}\sim7.6\pm0.4\times10^{43}$ erg s$^{-1}$, which exceeds those of the most X-ray luminous starburst galaxies, in support of an AGN driving the optical variability. Based on the broad H$\beta$ from SDSS spectrum, we estimate a virial BH mass of $M_{\bullet}\approx10^{6.43}$-$10^{6.72}M_{\odot}$ (with the error denoting 1$\sigma$ statistical uncertainties only), consistent with the estimation from OzDES, making it the lowest mass AGN with redshift $>$ 0.4 detected in optical. We estimate the host galaxy stellar mass to be $M_{\ast}\sim10^{10.5\pm0.3}M_{\odot}$ based on modeling the multi-wavelength spectral energy distribution. DES J0218$-$0430 extends the $M_{\bullet}$-$M_{\ast}$ relation observed in luminous AGNs at $z\sim1$ to masses lower than being probed by previous work. Our work demonstrates the feasibility of using optical variability to identify low-mass AGNs at higher redshift in deeper synoptic surveys with direct implications for the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory., Comment: 13 pages, 8 figures, 1 table, accepted to MNRAS

We present an extensive ALMA spectroscopic follow-up programme of the $z\,{=}\,4.3$ structure SPT2349$-$56, one of the most actively star-forming proto-cluster cores known, to identify additional members using their [C{\sc ii}] 158\,$\mu$m and \mbox{CO(4--3)} lines. In addition to robustly detecting the 14 previously published galaxies in this structure, we identify a further 15 associated galaxies at $z\,{=}\,4.3$, resolving 55$\,{\pm}\,$5\,per cent of the 870-$\mu$m flux density at 0.5\,arcsec resolution compared to 21\,arcsec single-dish data. These galaxies are distributed into a central core containing 23 galaxies extending out to 300\,kpc in diameter, and a northern extension, offset from the core by 400\,kpc, containing three galaxies. We discovered three additional galaxies in a red {\it Herschel\/}-SPIRE source 1.5\,Mpc from the main structure, suggesting the existence of many other sources at the same redshift as SPT2349$-$56 that are not yet detected in the limited coverage of our data. An analysis of the velocity distribution of the central galaxies indicates that this region may be virialized with a mass of (9$\pm$5)$\,{\times}\,$10$^{12}$\,M$_{\odot}$, while the two offset galaxy groups are about 30 and 60\,per cent less massive and show significant velocity offsets from the central group. We calculate the [C{\sc ii}] and far-infrared number counts, and find evidence for a break in the [C{\sc ii}] luminosity function. We estimate the average SFR density within the region of SPT2349$-$56 containing single-dish emission (a proper diametre of 720\,kpc), assuming spherical symmetry, to be roughly 4$\,{\times}\,10^4$\,M$_{\odot}$\,yr$^{-1}$\,Mpc$^{-3}$; this may be an order of magnitude greater than the most extreme examples seen in simulations., Comment: Submitted to MNRAS Feb. 25, 2020. Accepted Apr. 29, 2020

We present a study of six far-infrared fine structure lines in the z=4.225 lensed dusty star-forming galaxy SPT0418-47 to probe the physical conditions of its InterStellar Medium (ISM). In particular, we report Atacama Pathfinder EXperiment (APEX) detections of the [OI]145um and [OIII]88um lines and Atacama Compact Array (ACA) detections of the [NII]122 and 205um lines. The [OI]145um / [CII]158um line ratio is ~5x higher compared to the average of local galaxies. We interpret this as evidence that the ISM is dominated by photo-dissociation regions with high gas densities. The line ratios, and in particular those of [OIII]88um and [NII]122um imply that the ISM in SPT0418-47 is already chemically enriched close to solar metallicity. While the strong gravitational amplification was required to detect these lines with APEX, larger samples can be observed with the Atacama Large Millimeter/submillimeter Array (ALMA), and should allow to determine if the observed dense, solar metallicity ISM is common among these highly star-forming galaxies., Comment: Accepted for publication in Astronomy & Astrophysics. Revised to match language-edited version and adding 1 missing citation

Water ($\rm H_{2}O$), one of the most ubiquitous molecules in the universe, has bright millimeter-wave emission lines easily observed at high-redshift with the current generation of instruments. The low excitation transition of $\rm H_{2}O$, p$-$$\rm H_{2}O$(202 $-$ 111) ($\nu_{rest}$ = 987.927 GHz) is known to trace the far-infrared (FIR) radiation field independent of the presence of active galactic nuclei (AGN) over many orders-of-magnitude in FIR luminosity (L$_{\rm FIR}$). This indicates that this transition arises mainly due to star formation. In this paper, we present spatially ($\sim$0.5 arcsec corresponding to $\sim$1 kiloparsec) and spectrally resolved ($\sim$100 kms$^{-1}$) observations of p$-$$\rm H_{2}O$(202 $-$ 111) in a sample of four strong gravitationally lensed high-redshift galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA). In addition to increasing the sample of luminous ($ > $ $10^{12}$L$_{\odot}$) galaxies observed with $\rm H_{2}O$, this paper examines the L$_{\rm H_{2}O}$/L$_{\rm FIR}$ relation on resolved scales for the first time at high-redshift. We find that L$_{\rm H_{2}O}$ is correlated with L$_{\rm FIR}$ on both global and resolved kiloparsec scales within the galaxy in starbursts and AGN with average L$_{\rm H_{2}O}$/L$_{\rm FIR}$ =$2.76^{+2.15}_{-1.21}\times10^{-5}$. We find that the scatter in the observed L$_{\rm H_{2}O}$/L$_{\rm FIR}$ relation does not obviously correlate with the effective temperature of the dust spectral energy distribution (SED) or the molecular gas surface density. This is a first step in developing p$-$$\rm H_{2}O$(202 $-$ 111) as a resolved star formation rate (SFR) calibrator.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of high-J CO lines ($J_\mathrm{up}=6$, 7, 8) and associated dust continuum towards five strongly lensed, dusty, star-forming galaxies (DSFGs) at redshift $z = 2.7$-5.7. These galaxies, discovered in the South Pole Telescope survey, are observed at $0.2''$-$0.4''$ resolution with ALMA. Our high-resolution imaging coupled with the lensing magnification provides a measurement of the structure and kinematics of molecular gas in the background galaxies with spatial resolutions down to kiloparsec scales. We derive visibility-based lens models for each galaxy, accurately reproducing observations of four of the galaxies. Of these four targets, three show clear velocity gradients, of which two are likely rotating disks. We find that the reconstructed region of CO emission is less concentrated than the region emitting dust continuum even for the moderate-excitation CO lines, similar to what has been seen in the literature for lower-excitation transitions. We find that the lensing magnification of a given source can vary by 20-50% across the line profile, between the continuum and line, and between different CO transitions. We apply Large Velocity Gradient (LVG) modeling using apparent and intrinsic line ratios between lower-J and high-J CO lines. Ignoring these magnification variations can bias the estimate of physical properties of interstellar medium of the galaxies. The magnitude of the bias varies from galaxy to galaxy and is not necessarily predictable without high resolution observations., Comment: 19 pages, 12 figures, accepted by ApJ

SPT0346-52 is one of the most most luminous and intensely star-forming galaxies in the universe, with L_FIR > 10^13 L_sol and Sigma_SFR ~ 4200 M_sol yr^-1 kpc^-2. In this paper, we present ~0.15'' ALMA observations of the [CII]158micron emission line in this z=5.7 dusty star-forming galaxy. We use a pixellated lensing reconstruction code to spatially and kinematically resolve the source-plane [CII] and rest-frame 158 micron dust continuum structure at ~700 pc (~0.12'') resolution. We discuss the [CII] deficit with a pixellated study of the L_[CII]/L_FIR ratio in the source plane. We find that individual pixels within the galaxy follow the same trend found using unresolved observations of other galaxies, indicating that the deficit arises on scales <700 pc. The lensing reconstruction reveals two spatially and kinematically separated components (~1 kpc and ~500 km s^-1 apart) connected by a bridge of gas. Both components are found to be globally unstable, with Toomre Q instability parameters << 1 everywhere. We argue that SPT0346-52 is undergoing a major merger, which is likely driving the intense and compact star formation., Comment: Accepted for publication in ApJ