Researcher profile

Caitlin M. Casey

Caitlin M. Casey contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2025arXiv

COSMOS2025: The COSMOS-Web galaxy catalog of photometry, morphology, redshifts, and physical parameters from JWST, HST, and ground-based imaging

We present COSMOS2025, the COSMOS-Web catalog of photometry, morphology, photometric redshifts and physical parameters for more than 700,000 galaxies in the Cosmic Evolution Survey (COSMOS) field. This catalog is based on our \textit{James Webb Space Telescope} 255\,h COSMOS-Web program, which provides deep near-infrared imaging in four NIRCam (F115W, F150W, F277W, F444W) and one MIRI (F770W) filter over the central $\sim 0.54 {\, \rm deg}^2$ ($\sim 0.2 {\, \rm deg}^2$ for MIRI) in COSMOS. These data are combined with ground- and space-based data to derive photometric measurements of NIRCam-detected sources using both fixed-aperture photometry (on the space-based bands) and a profile-fitting technique on all 37 bands spanning 0.3-8 micron. We provide morphology for all sources from complementary techniques including profile fitting and machine-learning classification. We derive photometric redshifts, physical parameters and non-parametric star formation histories from spectral energy distribution (SED) fitting. The catalog has been extensively validated against previous COSMOS catalogs and other surveys. Photometric redshift accuracy measured using spectroscopically confirmed galaxies out to $z\sim9$ reaches $σ_{\rm MAD} = 0.012$ at $m_{\rm F444W}<28$ and remains at $σ_{\rm MAD} \lesssim 0.03$ as a function of magnitude, color, and galaxy type. This represents a factor of $\sim 2$ improvement at 26 AB mag compared to COSMOS2020. The catalog is approximately 80\% complete at $\log(M_{\star}/{\rm M}_{\odot}) \sim 9$ at $z \sim 10$ and at $\log(M_{\star}/{\rm M}_{\odot}) \sim 7$ at $z \sim 0.2$, representing a gain of 1\,dex compared to COSMOS2020. COSMOS2025 represents the definitive COSMOS-Web catalog. It is provided with complete documentation, together with redshift probability distributions, and it is ready for scientific exploitation today.

preprint2022arXiv

No Redshift Evolution of Dust Temperatures from $0 < z < 2$

Some recent literature has claimed there to be an evolution in galaxies&#39; dust temperatures towards warmer (or colder) spectral energy distributions (SEDs) between low and high redshift. These conclusions are driven by both theoretical models and empirical measurement. Such claims sometimes contradict one another and are prone to biases in samples or SED fitting techniques. What has made direct comparisons difficult is that there is no uniform approach to fitting galaxies&#39; infrared/millimeter SEDs. Here we aim to standardize the measurement of galaxies&#39; dust temperatures with a python-based SED fitting procedure, MCIRSED. We draw on reference datasets observed by IRAS, Herschel, and Scuba-2 to test for redshift evolution out to $z\sim2$. We anchor our work to the L$_{IR}$-$λ_{peak}$ plane, where there is an empirically observed anti-correlation between IR luminosity and rest-frame peak wavelength (an observational proxy for luminosity-weighted dust temperature) such that $λ_{peak}=λ_{t}({L_{IR}}/{L_{t}})^η$ where $η=-0.09\pm0.01$, L$_{t}=10^{12}$ L$_{\odot}$, and $λ_{t}=92\pm2μ$m. We find no evidence for redshift evolution of galaxies&#39; temperatures, or $λ_{peak}$, at fixed L$_{IR}$ from $0<z<2$ with >99.99% confidence. Our finding does not preclude evolution in dust temperatures at fixed stellar mass, which is expected from a non-evolving L$_{IR}$-$λ_{peak}$ relation and a strongly evolving SFR-M$_\star$ relation. The breadth of dust temperatures at a given L$_{IR}$ is likely driven by variation in galaxies&#39; dust geometries and sizes and does not evolve. Testing for L$_{IR}$-$λ_{peak}$ evolution toward higher redshift ($z\sim5-6$) requires better sampling of galaxies&#39; dust SEDs near their peaks (observed $\sim$200-600$μ$m) with $<$1 mJy sensitivity. This poses a significant challenge to current instrumentation.

preprint2022arXiv

Probing cold gas in a massive, compact star-forming galaxy at z=6

Observations of low order CO transitions represent the most direct way to study galaxies&#39; cold molecular gas, the fuel of star formation. Here we present the first detection of CO(2-1) in a galaxy lying on the main-sequence of star-forming galaxies at z>6. Our target, G09-83808 at z=6.03, has a short depletion time-scale of T_dep~50Myr and a relatively low gas fraction of M_gas/M_star=0.30 that contrasts with those measured for lower redshift main-sequence galaxies. We conclude that this galaxy is undergoing a starburst episode with a high star formation efficiency that might be the result of gas compression within its compact rotating disk. Its starburst-like nature is further supported by its high star formation rate surface density, thus favoring the use of the Kennicutt-Schmidt relation as a more precise diagnostic diagram for starbursts. Without further significant gas accretion, this galaxy would become a compact, massive quiescent galaxy at z~5.5. In addition, we find that the calibration for estimating ISM masses from dust continuum emission satisfactorily reproduces the gas mass derived from the CO(2-1) transition (within a factor of ~2). This is in line with previous studies claiming a small redshift evolution in the gas-to-dust ratio of massive, metal-rich galaxies. In the absence of gravitational amplification, this detection would have required of order ~1000h of observing time. The detection of cold molecular gas in unlensed star-forming galaxies at high redshifts is thus prohibitive with current facilities and requires a ten-fold improvement in sensitivity, such as that envisaged for the ngVLA.

preprint2022arXiv

Searching Far and Long I: Pilot ALMA 2mm Follow-up of Bright Dusty Galaxies as a Redshift Filter

A complete census of dusty star-forming galaxies (DSFGs) at early epochs is necessary to constrain the obscured contribution to the cosmic star formation rate density (CSFRD), however DSFGs beyond $z \sim 4$ are both rare and hard to identify from photometric data alone due to degeneracies in submillimeter photometry with redshift. Here, we present a pilot study obtaining follow-up Atacama Large Millimeter Array (ALMA) $2\,$mm observations of a complete sample of 39 $850\,\rmμm$-bright dusty galaxies in the SSA22 field. Empirical modeling suggests $2\,$mm imaging of existing samples of DSFGs selected at $850\,\rmμm - 1\,$mm can quickly and easily isolate the &#34;needle in a haystack&#34; DSFGs that sit at $z>4$ or beyond. Combining archival submillimeter imaging with our measured ALMA $2\,$mm photometry ($1σ\sim 0.08\,$mJy$\,$beam$^{-1}$ rms), we characterize the galaxies&#39; IR SEDs and use them to constrain redshifts. With available redshift constraints fit via the combination of six submillimeter bands, we identify 6/39 high-$z$ candidates each with $>50\%$ likelihood to sit at $z > 4$, and find a positive correlation between redshift and $2\,$mm flux density. Specifically, our models suggest the addition of $2\,$mm to a moderately constrained IR SED will improve the accuracy of a millimeter-derived redshift from $Δz/(1+z) = 0.3$ to $Δz/(1+z) = 0.2$. Our IR SED characterizations provide evidence for relatively high emissivity spectral indices ($\langle β\rangle = 2.4\pm0.3$) in the sample. We measure that especially bright ($S_{850\rmμm}>5.55\,$mJy) DSFGs contribute $\sim10$% to the cosmic-averaged CSFRD from $2<z<5$, confirming findings from previous work with similar samples.

preprint2022arXiv

Snowmass2021 CMB-HD White Paper

CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10 h Mpc^(-1)), which probes dark matter particle properties. It will also allow 2.) measurements of the thermal and kinetic Sunyaev-Zel&#39;dovich effects on small scales to map the gas density and velocity, another probe of cosmic structure. In addition, CMB-HD would allow us to cross critical thresholds: 3.) ruling out or detecting any new, light (< 0.1 eV) particles that were in thermal equilibrium with known particles in the early Universe, 4.) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe, and 5.) ruling out or detecting inflationary magnetic fields. CMB-HD would also provide world-leading constraints on 6.) axion-like particles, 7.) cosmic birefringence, 8.) the sum of the neutrino masses, and 9.) the dark energy equation of state. The CMB-HD survey would be delivered in 7.5 years of observing 20,000 square degrees of sky, using two new 30-meter-class off-axis crossed Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

preprint2021arXiv

Characterization of Two 2mm-detected Optically-Obscured Dusty Star-Forming Galaxies

The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform a rigorous analysis of all available photometric data to fit spectral energy distributions and estimate redshifts before deriving physical properties based on our findings. Our results suggest the two galaxies, called MORA-5 and MORA-9, represent two extremes of the &#34;OIR-dark&#34; class of DSFGs. MORA-5 ($z_{\rm phot}=4.3^{+1.5}_{-1.3}$) is a significantly more active starburst with a star-formation rate of 830$^{+340}_{-190}$M$_\odot$yr$^{-1}$ compared to MORA-9 ($z_{\rm phot}=4.3^{+1.3}_{-1.0}$) whose star-formation rate is a modest 200$^{+250}_{-60}$M$_\odot$yr$^{-1}$. Based on the stellar masses (M$_{\star}\approx10^{10-11}$M$_\odot$), space density ($n\sim(5\pm2)\times10^{-6}$Mpc$^{-3}$, which incorporates two other spectroscopically confirmed OIR-dark DSFGs in the MORA sample at $z=4.6$ and $z=5.9$), and gas depletion timescales ($<1$Gyr) of these sources, we find evidence supporting the theory that OIR-dark DSFGs are the progenitors of recently discovered $3<z<4$ massive quiescent galaxies.

preprint2021arXiv

First HETDEX Spectroscopic Determinations of Ly$α$ and UV Luminosity Functions at $z=2-3$: Bridging a Gap Between Faint AGN and Bright Galaxies

We present Ly$α$ and ultraviolet-continuum (UV) luminosity functions (LFs) of galaxies and active galactic nuclei (AGN) at $z=2.0-3.5$ determined by the un-targetted optical spectroscopic survey of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We combine deep Subaru imaging with HETDEX spectra resulting in $11.4$ deg$^2$ of fiber-spectra sky coverage, obtaining $18320$ galaxies spectroscopically identified with Ly$α$ emission, $2126$ of which host type 1 AGN showing broad (FWHM$~>1000$ km s$^{-1}$) Ly$α$ emission lines. We derive the Ly$α$ (UV) LF over 2 orders of magnitude covering bright galaxies and AGN in $\log L_\mathrm{Lyα}/\mathrm{[erg~s^{-1}]}=43.3-45.5$ ($-27<M_\mathrm{UV}<-20$) by the $1/V_\mathrm{max}$ estimator. Our results reveal the bright-end hump of the Ly$α$ LF is composed of type 1 AGN. In conjunction with previous spectroscopic results at the faint end, we measure a slope of the best-fit Schechter function to be $α_\mathrm{Sch}=-1.70^{+0.13}_{-0.14}$, which indicates $α_\mathrm{Sch}$ steepens from $z=2-3$ towards high redshift. Our UV LF agrees well with previous AGN UV LFs, and extends to faint-AGN and bright-galaxy regimes. The number fraction of Ly$α$-emitting objects ($X_\mathrm{LAE}$) increases from $M_\mathrm{UV}^*\sim-21$ to bright magnitude due to the contribution of type 1 AGN, while previous studies claim that $X_\mathrm{Lyα}$ decreases from faint magnitude to $M_\mathrm{UV}^*$, suggesting a valley in the $X_\mathrm{Lyα}-$magnitude relation at $M_\mathrm{UV}^*$. Comparing our UV LF of type 1 AGN at $z=2-3$ with those at $z=0$, we find that the number density of faint ($M_\mathrm{UV}>-21$) type 1 AGN increases from $z\sim2$ to $z\sim0$ as opposed to the evolution of bright ($M_\mathrm{UV}<-21$) type 1 AGN, suggesting the AGN downsizing in the rest-frame UV luminosity.

preprint2021arXiv

Mapping Obscuration to Reionization with ALMA (MORA): 2mm Efficiently Selects the Highest-Redshift Obscured Galaxies

We present the characteristics of 2mm-selected sources from the largest Atacama Large Millimeter and submillimeter Array (ALMA) blank-field contiguous survey conducted to-date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184arcmin$^2$ at 2mm. Twelve of the thirteen detections above 5$σ$ are attributed to emission from galaxies, eleven of which are dominated by cold dust emission. These sources have a median redshift of $\langle z_{\rm 2mm}\rangle=3.6^{+0.4}_{-0.3}$ primarily based on optical/near-infrared (OIR) photometric redshifts with some spectroscopic redshifts, with 77$\pm$11% of sources at $z>3$ and 38$\pm$12% of sources at $z>4$. This implies that 2mm selection is an efficient method for identifying the highest redshift dusty star-forming galaxies (DSFGs). Lower redshift DSFGs ($z<3$) are far more numerous than those at $z>3$ yet likely to drop out at 2mm. MORA shows that DSFGs with star-formation rates in excess of 300M$_\odot$ yr$^{-1}$ and relative rarity of $\sim$10$^{-5}$ Mpc$^{-3}$ contribute $\sim$30% to the integrated star-formation rate density between $3<z<6$. The volume density of 2mm-selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at $z>2$. Analysis of MORA sources&#39; spectral energy distributions hint at steeper empirically-measured dust emissivity indices than typical literature studies, with $\langleβ\rangle=2.2^{+0.5}_{-0.4}$. The MORA survey represents an important step in taking census of obscured star-formation in the Universe&#39;s first few billion years, but larger area 2mm surveys are needed to more fully characterize this rare population and push to the detection of the Universe&#39;s first dusty galaxies.

preprint2020arXiv

Emergence of an Ultra-Red Ultra-Massive Galaxy Cluster Core at $z=4$

Recent simulations and observations of massive galaxy cluster evolution predict that the majority of stellar mass build up happens within cluster members by $z=2$, before cluster virialization. Protoclusters rich with dusty, star-forming galaxies (DSFGs) at $z>3$ are the favored candidate progenitors for these massive galaxy clusters at $z\sim0$. We present here the first study analyzing stellar emission along with cold dust and gas continuum emission in a spectroscopically confirmed $z=4.002$ protocluster core rich with DSFGs, the Distant Red Core (DRC). We combine new HST and Spitzer data with existing Gemini, Herschel, and ALMA observations to derive individual galaxy-level properties, and compare them to coeval field and other protocluster galaxies. All of the protocluster members are massive ($>10^{10}$ M$_\odot$), but not significantly more so than their coeval field counterparts. Within uncertainty, all are nearly indistinguishable from galaxies on the star-forming vs. stellar mass main-sequence relationship, and on the star formation efficiency plane. Assuming no future major influx of fresh gas, we estimate that these gaseous DSFGs will deplete their gas reservoirs in $\sim300$ Myr, becoming the massive quiescent ellipticals dominating cluster cores by $z\sim3$. Using various methodologies, we derive a total $z=4$ halo mass of $\sim10^{14}$ M$_\odot$, and estimate that the DRC will evolve to become an ultra-massive cluster core of mass $\gtrsim10^{15}$ M$_\odot$ by $z=0$.

preprint2020arXiv

Far-Infrared Photometric Redshifts: A New Approach to a Highly Uncertain Enterprise

I present a new approach at deriving far-infrared photometric redshifts for galaxies based on their reprocessed emission from dust at rest-frame far-infrared through millimeter wavelengths. Far-infrared photometric redshifts (&#34;FIR-$z$&#34;) have been used over the past decade to derive redshift constraints for highly obscured galaxies that lack photometry at other wavelengths like the optical/near-infrared. Most literature FIR-z fits are performed through $χ^2$minimization to a single galaxy&#39;s far-infrared template spectral energy distribution (SED). The use of a single galaxy template, or modest set of templates, can lead to an artificially low uncertainty estimate on FIR-$z$&#39;s because real galaxies display a wide range in intrinsic dust SEDs. I use the observed distribution of galaxy SEDs (for well-constrained samples across $0<z<5$) to motivate a new far-infrared through millimeter photometric redshift technique called MMpz. The MMpz algorithm asserts that galaxies are most likely drawn from the empirically observed relationship between rest-frame peak wavelength, $λ_{\rm peak}$, and total IR luminosity, L$_{\rm IR}$; the derived photometric redshift accounts for the measurement uncertainties and intrinsic variation in SEDs at the inferred L$_{\rm IR}$, as well as heating from the CMB at $z>5$. The MMpz algorithm has a precision of $σ_{Δz/(1+z)}\approx0.3-0.4$, similar to single-template fits, while providing a more accurate estimate of the FIR-$z$ uncertainty with reduced chi-squared of order $\mathcal{O}(χ^2_ν)=1$, compared to alternative far-infrared photometric redshift techniques (with $\mathcal{O}(χ^2_ν)\approx10-10^{3}$).

preprint2020arXiv

SuperCLASS -- I. The Super CLuster Assisted Shear Survey: Project overview and Data Release 1

The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers.

preprint2020arXiv

SuperCLASS -- II: Photometric Redshifts and Characteristics of Spatially-Resolved $μ$Jy Radio Sources

We present optical and near-infrared imaging covering a $\sim$1.53 deg$^2$ region in the Super-Cluster Assisted Shear Survey (SuperCLASS) field, which aims to make the first robust weak lensing measurement at radio wavelengths. We derive photometric redshifts for $\approx$176,000 sources down to $i^\prime_{\rm AB}\sim24$ and present photometric redshifts for 1.4 GHz $e$-MERLIN and VLA detected radio sources found in the central 0.26 deg$^{2}$. We compile an initial catalog of 149 radio sources brighter than S$_{1.4}>75$ $μ$Jy and find their photometric redshifts span $0<z_{\rm phot}<4$ with radio luminosities between $10^{21}-10^{25}$ W Hz$^{-1}$, with medians of $\langle z \rangle =0.55$ and $\langle L_{1.4}\rangle =1.9\times10^{23}$ W Hz$^{-1}$ respectively. We find 95% of the \uJy\ radio source sample (141/149) have SEDs best fit by star-forming templates while 5% (8/149) are better fit by AGN. Spectral indices are calculated for sources with radio observations from VLA and GMRT at 325 MHz, with an average spectral slope of $α=0.59\pm0.04$. Using the full photometric redshift catalog we construct a density map at the redshift of the known galaxy clusters, $z=0.20\pm0.08$. Four of the five clusters are prominently detected at $>7 σ$ in the density map and we confirm the photometric redshifts are consistent with previously measured spectra from a few galaxies at the cluster centers.

preprint2020arXiv

SuperCLASS -- III. Weak lensing from radio and optical observations in Data Release 1

We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1 data set). Using standard methodologies on the optical data only we make a significant (10 sigma) detection of the weak lensing signal (a shear power spectrum) due to the massive supercluster of galaxies in the targeted region. For the radio data we develop a new method to measure the shapes of galaxies from the interferometric data, and we construct a simulation pipeline to validate this method. We then apply this analysis to our radio observations, treating the e-MERLIN and VLA data independently. We achieve source densities of 0.5 per square arcmin in the VLA data and 0.06 per square arcmin in the e-MERLIN data, numbers which prove too small to allow a detection of a weak lensing signal in either the radio data alone or in cross-correlation with the optical data. Finally, we show preliminary results from a visibility-plane combination of the data from e-MERLIN and VLA which will be used for the forthcoming full SuperCLASS data release. This approach to data combination is expected to enhance both the number density of weak lensing sources available and the fidelity with which their shapes can be measured.

preprint2020arXiv

Three Dusty Star Forming Galaxies at $z\sim1.5$: Mergers and Disks on the Main Sequence

The main sequence of galaxies, a correlation between the star formation rates and stellar masses of galaxies, has been observed out to $z\sim4$. Galaxies within the scatter of the correlation are typically interpreted to be secularly evolving while galaxies with star formation rates elevated above the main sequence are interpreted to be undergoing interactions or to be Toomre-unstable disks with starbursting clumps. In this paper we investigate the recent merger histories of three dusty star forming galaxies, identified by their bright submillimeter emission at $z\sim1.5$. We analyze rest-frame optical and UV imaging, rest-frame optical emission line kinematics using slit spectra obtained with MOSFIRE on Keck I, and calculate Gini and M$_{20}$ statistics for each galaxy and conclude two are merger-driven while the third is an isolated disk galaxy. The disk galaxy lies $\sim$4$\times$ above the main sequence, one merger lies within the scatter of the main sequence, and one merger lies $\sim$4$\times$ below the main sequence. This hints that the location of a galaxy with respect to the main sequence may not be a useful discriminator of the recent star formation history of high-M$_{\star}$ galaxies at $z\sim1$.

preprint2019arXiv

Physical Characterization of an Unlensed Dusty Star-Forming Galaxy at $z=5.85$

We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9&#39;&#39;), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$μ$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $τ_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr).