Source author record

Caitlin M. Casey

Caitlin M. Casey appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

35 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' 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' infrared/millimeter SEDs. Here we aim to standardize the measurement of galaxies' 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' 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' 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' 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' 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 "needle in a haystack" 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' 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'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 "OIR-dark" 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' 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'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'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 ("FIR-$z$") 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'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$'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''), 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).

preprint2016arXiv

Large scale structure around a z=2.1 cluster

The most prodigious starburst galaxies are absent in massive galaxy clusters today, but their connection with large scale environments is less clear at $z\gtrsim2$. We present a search of large scale structure around a galaxy cluster core at $z=2.095$ using a set of spectroscopically confirmed galaxies. We find that both color-selected star-forming galaxies (SFGs) and dusty star-forming galaxies (DSFGs) show significant overdensities around the $z=2.095$ cluster. A total of 8 DSFGs (including 3 X-ray luminous active galactic nuclei, AGNs) and 34 SFGs are found within a 10 arcmin radius (corresponds to $\sim$15 cMpc at $z\sim2.1$) from the cluster center and within a redshift range of $Δz=0.02$, which leads to galaxy overdensities of $δ_{\rm DSFG}\sim12.3$ and $δ_{\rm SFG}\sim2.8$. The cluster core and the extended DSFG- and SFG-rich structure together demonstrate an active cluster formation phase, in which the cluster is accreting a significant amount of material from large scale structure while the more mature core may begin to virialize. Our finding of this DSFG-rich structure, along with a number of other protoclusters with excess DSFGs and AGNs found to date, suggest that the overdensities of these rare sources indeed trace significant mass overdensities. However, it remains puzzling how these intense star formers are triggered concurrently. Although an increased probability of galaxy interactions and/or enhanced gas supply can trigger the excess of DSFGs, our stacking analysis based on 850 $μ$m images and morphological analysis based on rest-frame optical imaging do not show such enhancements of merger fraction and gas content in this structure.

preprint2016arXiv

Spatially Resolved Spectroscopy of Sub-mm Galaxies at z~2

We present near-infrared integral-field spectroscopic observations targeting H$α$ in eight sub-millimeter galaxies (SMGs) at $z$=1.3-2.5 using VLT/SINFONI, obtaining significant detections for six of them. The star formation rates derived from the H$α$ emission are $\sim$100 M$_\odot$yr$^{-1}$, which account for only $\sim$ 20-30\% of the infrared-derived values, thus suggesting that these systems are very dusty. Two of these systems present [NII]/H$α$ ratios indicative of the presence of an Active Galactic Nucleus (AGN). We mapped the spatial distribution and kinematics of the star forming regions in these galaxies on kpc-scales. In general, the H$α$ morphologies tend to be highly irregular and/or clumpy, showing spatial extents of $\sim$3-11~kpc. We find evidence for significant spatial offsets, of $\sim$0.1-0.4$"$ or 1.2-3.4 kpc, between the H$α$ and the continuum emission in three of the sources. Performing a kinemetry analysis we conclude that the majority of the sample is not consistent with disk-like rotation-dominated kinematics. Instead, they tend to show irregular and/or clumpy and turbulent velocity and velocity dispersion fields. This can be interpreted as evidence for scenario in which these extreme star formation episodes are triggered by galaxy-galaxy interactions and major mergers. In contrast to recent results for SMGs, these sources appear to follow the same relations between gas and star forming rate densities as less luminous and/or normal star forming galaxies.

preprint2016arXiv

The Ubiquity of Coeval Starbursts in Massive Galaxy Cluster Progenitors

The Universe's largest galaxy clusters likely built the majority of their massive $>10^{11} M_\odot$ galaxies in simultaneous, short-lived bursts of activity well before virialization. This conclusion is reached from emerging datasets on $z>2$ proto-clusters and the characteristics of their member galaxies, in particular, rare starbursts and ultraluminous active galactic nuclei (AGN). The most challenging observational hurdle in identifying such structures is their very large volumes, $\sim10^{4}$ comoving Mpc$^3$ at $z>2$, subtending areas $\sim$half a degree on the sky. Thus the contrast afforded by an overabundance of very rare galaxies in comparison to the background can more easily distinguish overdense structures from the surrounding, normal density field. Five $2<z<3$ proto-clusters from the literature are discussed in detail and are found to contain up to 12 dusty starbursts or luminous AGN galaxies each, a phenomenon that is unlikely to occur by chance even in overdense environments. Measurements of gas depletion times suggest that they are indeed short-lived on $\sim$100 Myr timescales, and the probability of finding a structure containing more than 8 such systems is $\sim$0.2\%, unless their `triggering' is correlated on very large spatial scales, $\sim$10 Mpc across. The volume density of starburst-rich proto-clusters is found to be comparable to all $>10^{15} M_\odot$ galaxy clusters in the nearby Universe, a factor of five larger than expected in some simulations. Some tension yet exists between measurements of their volume density of starburst-rich proto-clusters and the expectation that they are generated via short-lived episodes. However, improved observations of proto-clusters over large regions of sky will certainly shed more light on the assembly of galaxy clusters, and whether or not they build their galaxies through episodic bursts as suggested here. [abridged]

preprint2015arXiv

A Turnover in the Galaxy Main Sequence of Star Formation at $M_{*} \sim 10^{10} M_{\odot}$ for Redshifts $z < 1.3$

The relationship between galaxy star formation rates (SFR) and stellar masses ($M_\ast$) is re-examined using a mass-selected sample of $\sim$62,000 star-forming galaxies at $z \le 1.3$ in the COSMOS 2-deg$^2$ field. Using new far-infrared photometry from $Herschel$-PACS and SPIRE and $Spitzer$-MIPS 24 $μ$m, along with derived infrared luminosities from the NRK method based on galaxies' locations in the restframe color-color diagram $(NUV - r)$ vs. $(r - K)$, we are able to more accurately determine total SFRs for our complete sample. At all redshifts, the relationship between median $SFR$ and $M_\ast$ follows a power-law at low stellar masses, and flattens to nearly constant SFR at high stellar masses. We describe a new parameterization that provides the best fit to the main sequence and characterizes the low mass power-law slope, turnover mass, and overall scaling. The turnover in the main sequence occurs at a characteristic mass of about $M_{0} \sim 10^{10} M_{\odot}$ at all redshifts. The low mass power-law slope ranges from 0.9-1.3 and the overall scaling rises in SFR as a function of $(1+z)^{4.12 \pm 0.10}$. A broken power-law fit below and above the turnover mass gives relationships of $SFR \propto M_{*}^{0.88 \pm 0.06}$ below the turnover mass and $SFR \propto M_{*}^{0.27 \pm 0.04}$ above the turnover mass. Galaxies more massive than $M_\ast \gtrsim 10^{10}\ M_{\rm \odot}$ have on average, a much lower specific star formation rate (sSFR) than would be expected by simply extrapolating the traditional linear fit to the main sequence found for less massive galaxies.

preprint2015arXiv

Kinematic classifications of local interacting galaxies: implications for the merger/disk classifications at high-z

The classification of galaxy mergers and isolated disks is key for understanding the relative importance of galaxy interactions and secular evolution during the assembly of galaxies. The kinematic properties of galaxies as traced by emission lines have been used to suggest the existence of a significant population of high-z star-forming galaxies consistent with isolated rotating disks. However, recent studies have cautioned that post-coalescence mergers may also display disk-like kinematics. To further investigate the robustness of merger/disk classifications based on kinematic properties, we carry out a systematic classification of 24 local (U)LIRGs spanning a range of galaxy morphologies: from isolated spiral galaxies, ongoing interacting systems, to fully merged remnants. We artificially redshift the WiFeS observations of these local (U)LIRGs to z=1.5 to make a realistic comparison with observations at high-z, and also to ensure that all galaxies have the same spatial sampling of ~900 pc. Using both kinemetry-based and visual classifications, we find that the reliability of kinematic classification shows a strong trend with the interaction stage of galaxies. Mergers with two nuclei and tidal tails have the most distinct kinematic properties compared to isolated disks, whereas a significant population of the interacting disks and merger remnants are indistinguishable from isolated disks. The high fraction of late-stage mergers showing disk-like kinematics reflects the complexity of the dynamics during galaxy interactions. However, the exact fractions of misidentified disks and mergers depend on the definition of kinematic asymmetries and the classification threshold when using kinemetry-based classifications. Our results suggest that additional indicators such as morphologies traced by stars or molecular gas are required to further constrain the merger/disk classifications at high-z.

preprint2015arXiv

Next Generation Very Large Array Memo No. 8 Science Working Group 3: Galaxy Assembly through Cosmic Time

The Next-Generation Very Large Array (ngVLA) will be critical for understanding how galaxies are built and evolve at the earliest epochs. The sensitivity and frequency coverage will allow for the detection of cold gas and dust in `normal' distant galaxies, including the low-J transitions of molecular gas tracers such as CO, HNC, and HCO+; synchrotron and free-free continuum emission; and even the exciting possibility of thermal dust emission at the highest (z~7) redshifts. In particular, by enabling the total molecular gas reservoirs to be traced to unprecedented sensitivities across a huge range of epochs simultaneously -- something no other radio or submillimeter facility will be capable of -- the detection of the crucial low-J transitions of CO in a diverse body of galaxies will be the cornerstone of ngVLA's contribution to high-redshift galaxy evolution science. The ultra-wide bandwidths will allow a complete sampling of radio SEDs, as well as the detection of emission lines necessary for spectroscopic confirmation of elusive dusty starbursts. The ngVLA will also deliver unique contributions to our understanding of cosmic magnetism and to science accessible through microwave polarimetry. Finally, the superb angular resolution will move the field beyond detection experiments and allow detailed studies of the morphology and dynamics of these systems, including dynamical modeling of disks/mergers, determining the properties of outflows, measuring black hole masses from gas disks, and resolving multiple AGN nuclei. We explore the contribution of a ngVLA to these areas and more, as well as synergies with current and upcoming facilities including ALMA, SKA, large single-dish submillimeter observatories, GMT/TMT, and JWST.

preprint2014arXiv

A comparison of the morphological properties between local and z~1 infrared luminous galaxies. Are local and high-z (U)LIRGs different?

Ultraluminous and luminous infrared galaxies (ULIRGs and LIRGs) are the most extreme star-forming galaxies in the universe, and dominate the total star formation rate density at z>1. In the local universe (z<0.3), the majority of ULIRGs and a significant portion of LIRGs are triggered by interactions between gas-rich spiral galaxies, yet it is unclear if this is still the case at high-z. To investigate the relative importance of galaxy interactions in infrared luminous galaxies, we carry out a comparison of optical morphological properties between local (U)LIRGs and (U)LIRGs at z=0.5-1.5 based on the same sample selection, morphology classification scheme, and optical morphology at similar rest-frame wavelengths. In addition, we quantify the systematics in comparing local and high-z datasets by constructing a redshifted dataset from local (U)LIRGs, in which its data quality mimics the high-z dataset. Based on the Gini-M20 classification scheme, we find that the fraction of interacting systems decreases by ~8% from local to z<~1, and it is consistent with the reduction between local and redshifted datasets (6(+14-6)%). Based on visual classifications, the merger fraction of local ULIRGs is found to be ~20% lower compared to published results, and the reduction due to redshifiting is 15(+10-8)%. Consequently, the differences of merger fractions between local and z<~1 (U)LIRGs is only ~17%. These results demonstrate that there is no strong evolution in the fraction of (U)LIRGs classified as mergers at least out to z~1. At z>1, the morphology types of ~30% of (U)LIRGs can not be determined due to their faintness in the F814W-band, and thus the merger fraction measured at z>1 suffers from large uncertainties.

preprint2014arXiv

Confirming the Quiescent Galaxy Population out to $z=3$: A Stacking Analysis of Mid-, Far-Infrared and Radio Data

We present stringent constraints on the average mid-, far-infrared and radio emissions of $\sim$14200 quiescent galaxies (QGs), identified out to $z=3$ in the COSMOS field via their rest-frame NUV$-$r and r$-$J colors, and with stellar masses $M_{\star}=10^{9.8-12.2} \,M_{\odot} $. Stacking in deep Spitzer (MIPS $24\,μ$m), Herschel (PACS and SPIRE), and VLA (1.4 GHz) maps reveals extremely low dust-obscured star formation rates for QGs (SFR $<0.1-3\,M_{\odot}$yr$^{-1}$ at $z \leqslant 2$ and $<6-18\,M_{\odot}$yr$^{-1}$ at $z > 2$), consistent with the low unobscured SFRs ($<0.01-1.2\,M_{\odot}$yr$^{-1}$) inferred from modeling their ultraviolet-to-near-infrared photometry. The average SFRs of QGs are $>10\times$ below those of star-forming galaxies (SFGs) within the $M_{\star}$- and $z$-ranges considered. The stacked 1.4 GHz signals (S/N $> 5$) are, if attributed solely to star formation, in excess of the total (obscured plus unobscured) SFR limits, suggestive of a widespread presence of low-luminosity active galactic nuclei (AGN) among QGs. Our results reaffirm the existence of a significant population QGs out to $z = 3$, thus corroborating the need for powerful quenching mechanism(s) to terminate star formation in galaxies at earlier epochs.

preprint2014arXiv

Dusty Star-Forming Galaxies at High Redshift

Far-infrared and submillimeter wavelength surveys have now established the important role of dusty, star-forming galaxies (DSFGs) in the assembly of stellar mass and the evolution of massive galaxies in the Universe. The brightest of these galaxies have infrared luminosities in excess of 10$^{13}$ L$_{\odot}$ with implied star-formation rates of thousands of solar masses per year. They represent the most intense starbursts in the Universe, yet many are completely optically obscured. Their easy detection at submm wavelengths is due to dust heated by ultraviolet radiation of newly forming stars. When summed up, all of the dusty, star-forming galaxies in the Universe produce an infrared radiation field that has an equal energy density as the direct starlight emission from all galaxies visible at ultraviolet and optical wavelengths. The bulk of this infrared extragalactic background light emanates from galaxies as diverse as gas-rich disks to mergers of intense starbursting galaxies. Major advances in far-infrared instrumentation in recent years, both space-based and ground-based, has led to the detection of nearly a million DSFGs, yet our understanding of the underlying astrophysics that govern the start and end of the dusty starburst phase is still in nascent stage. This review is aimed at summarizing the current status of DSFG studies, focusing especially on the detailed characterization of the best-understood subset (submillimeter galaxies, who were summarized in the last review of this field over a decade ago, Blain et al., 2002), but also the selection and characterization of more recently discovered DSFG populations. We review DSFG population statistics, their physical properties including dust, gas and stellar contents, their environments, and current theoretical models related to the formation and evolution of these galaxies.

preprint2013arXiv

A kinematic study of the Andromeda dwarf spheroidal system

We present a homogeneous kinematic analysis of red giant branch stars within 18 of the 28 Andromeda dwarf spheroidal (dSph) galaxies, obtained using the Keck I LRIS and Keck II DEIMOS spectrographs. Based on their g-i colors (taken with the CFHT MegaCam imager), physical positions on the sky, and radial velocities, we assign probabilities of dSph membership to each observed star. Using this information, the velocity dispersions, central masses and central densities of the dark matter halos are calculated for these objects, and compared with the properties of the Milky Way dSph population. We also measure the average metallicity ([Fe/H]) from the co-added spectra of member stars for each M31 dSph and find that they are consistent with the trend of decreasing [Fe/H] with luminosity observed in the Milky Way population. We find that three of our studied M31 dSphs appear as significant outliers in terms of their central velocity dispersion, And XIX, XXI and XXV, all of which have large half-light radii (>700 pc) and low velocity dispersions (sigma_v<5 km/s). In addition, And XXV has a mass-to-light ratio within its half-light radius of just [M/L]_{half}=10.3^{+7.0}_{-6.7}, making it consistent with a simple stellar system with no appreciable dark matter component within its 1 sigma uncertainties. We suggest that the structure of the dark matter halos of these outliers have been significantly altered by tides.

preprint2013arXiv

A SCUBA-2 850 micron survey of protoplanetary discs in the sigma Orionis cluster

We present the results from a large 850 micron survey of the sigma Orionis cluster using the SCUBA-2 camera on the James Clerk Maxwell Telescope. The 0.5-degree diameter circular region we surveyed contains 297 young stellar objects with an age estimated at about 3Myr. We detect 9 of these objects, 8 of which have infrared excesses from an inner disc. We also serendipitously detect 3 non-stellar sources at > 5sigma that are likely background submillimetre galaxies. The 9 detected stars have inferred disc masses ranging from 5 to about 17MJup, assuming similar dust properties as Taurus discs and an ISM gas-to-dust ratio of 100. There is a net positive signal toward the positions of the individually undetected infrared excess sources indicating a mean disc mass of 0.5 MJup . Stacking the emission toward those stars without infrared excesses constrains their mean disc mass to less than 0.3MJup, or an equivalent Earth mass in dust. The submillimetre luminosity distribution is significantly different from that in the younger Taurus region, indicating disc mass evolution as star forming regions age and the infrared excess fraction decreases. Submillimeter Array observations reveal CO emission toward 4 sources demonstrating that some, but probably not much, molecular gas remains in these relatively evolved discs. These observations provide new constraints on the dust and gas mass of protoplanetary discs during the giant planet building phase and provide a reference level for future studies of disc evolution.

preprint2013arXiv

Characterisation of SCUBA-2 450um and 850um-selected Galaxies in the COSMOS Field

We present deep 450um and 850um observations of a large, uniformly covered 394arcmin^2 area in the COSMOS field obtained with the SCUBA-2 instrument on the James Clerk Maxwell Telescope (JCMT). We achieve root-mean-square noise values of 4.13mJy at 450um and 0.80mJy at 850um. The differential and cumulative number counts are presented and compared to similar previous works. Individual point sources are identified at >3.6sigma significance, a threshold corresponding to a 3-5% sample contamination rate. We identify 78 sources at 450um and 99 at 850um, with flux densities S450=13-37mJy and S850=2-16mJy. Only 62-76% of 450um sources are 850um detected and 61-81% of 850um sources are 450um detected. The positional uncertainties at 450um are small (1-2.5") and therefore allow a precise identification of multiwavelength counterparts without reliance on detection at 24um or radio wavelengths; we find that only 44% of 450um-selected galaxies and 60% of 850um-sources have 24um or radio counterparts. 450um-selected galaxies peak at <z>=1.95+-0.19 and 850um=selected galaxies peak at <z>=2.16+-0.11. The two samples occupy similar parameter space in redshift and luminosity, while their median SED peak wavelengths differ by ~10-50um (translating to deltaTdust =8-12K, where 450um-selected galaxies are warmer). The similarities of the 450um and 850um populations, yet lack of direct overlap between them, suggests that submillimeter surveys conducted at any single far-infrared wavelength will be significantly incomplete (~>30%) at censusing infrared-luminous star formation at high-z.

preprint2013arXiv

Multi-Wavelength SEDs of Herschel Selected Galaxies in the COSMOS Field

We combine Herschel PACS and SPIRE maps of the full 2 deg^2 COSMOS field with existing multi-wavelength data to obtain template and model-independent optical-to-far-infrared spectral energy distributions (SEDs) for 4,218 Herschel-selected sources with log(L_IR/L_sun) = 9.4-13.6 and z = 0.02-3.54. Median SEDs are created by binning the optical to far-infrared (FIR) bands available in COSMOS as a function of infrared luminosity. Herschel probes rest-frame wavelengths where the bulk of the infrared radiation is emitted, allowing us to more accurately determine fundamental dust properties of our sample of infrared luminous galaxies. We find that the SED peak wavelength (lambda_peak) decreases and the dust mass (M_dust) increases with increasing total infrared luminosity (L_IR). In the lowest infrared luminosity galaxies (log(L_IR/L_sun) = 10.0-11.5), we see evidence of Polycyclic Aromatic Hydrocarbons (PAH) features (lambda=7-9 um), while in the highest infrared luminosity galaxies (L_IR > 10^12 L_sun) we see an increasing contribution of hot dust and/or power-law emission, consistent with the presence of heating from an active galactic nucleus (AGN). We study the relationship between stellar mass and star formation rate of our sample of infrared luminous galaxies and find no evidence that Herschel-selected galaxies follow the SFR/M_* "main sequence" as previously determined from studies of optically selected, star-forming galaxies. Finally, we compare the mid-infrared (MIR) to FIR properties of our infrared luminous galaxies using the previously defined diagnostic, IR8 = L_IR / L_8, and find that galaxies with L_IR > 10^11.3 L_sun tend to systematically lie above (~3-5) the IR8 "infrared main sequence", suggesting either suppressed PAH emission or an increasing contribution from AGN heating.

preprint2013arXiv

Studying Faint Ultra Hard X-ray Emission from AGN in GOALS LIRGs with Swift BAT

We present the first analysis of the all-sky Swift BAT ultra hard X-ray (14-195 keV) data for a targeted list of objects. We find the BAT data can be studied at 3x fainter limits than in previous blind detection catalogs based on prior knowledge of source positions and using smaller energy ranges for source detection. We determine the AGN fraction in 134 nearby (z<0.05) luminous infrared galaxies (LIRGS) from the GOALS sample. We find that LIRGs have a higher detection frequency than galaxies matched in stellar mass and redshift at 14-195 keV and 24-35 keV. In agreement with work at other wavelengths, the AGN detection fraction increases strongly at high IR luminosity with half of high luminosity LIRGs (50%, 6/12, log L_IR/L_sun>11.8) detected. The BAT AGN classification shows 97% (37/38) agreement with Chandra and XMM AGN classification using hardness ratios or detection of a iron K-alpha line. This confirms our statistical analysis and supports the use of the Swift BAT all-sky survey to study fainter populations of any category of sources in the ultra hard X-ray band. BAT AGN in LIRGs tend to show higher column densities with 40\pm9% showing 14-195 keV/2-10 keV hardness flux ratios suggestive of high or Compton-thick column densities (log N_H>24 cm^-2), compared to only 12\pm5% of non-LIRG BAT AGN. We also find that using specific energy ranges of the BAT detector can yield additional sources over total band detections with 24% (5/21) of detections in LIRGs at 24-35 keV not detected at 14-195 keV.

preprint2013arXiv

The masses of Local Group dwarf spheroidal galaxies: The death of the universal mass profile

We investigate the claim that all dwarf spheroidal galaxies (dSphs) reside within halos that share a common, universal mass profile as has been derived for dSphs of the Galaxy. By folding in kinematic information for 25 Andromeda dSphs, more than doubling the previous sample size, we find that a singular mass profile can not be found to fit all the observations well. Further, the best-fit dark matter density profile measured for solely the Milky Way dSphs is marginally discrepant (at just beyond the 1 sigma level) with that of the Andromeda dSphs, where a profile with lower maximum circular velocity, and hence mass, is preferred. The agreement is significantly better when three extreme Andromeda outliers, And XIX, XXI and XXV, all of which have large half-light radii (>600pc) and low velocity dispersions (sigma_v < 5km/s) are omitted from the sample. We argue that the unusual properties of these outliers are likely caused by tidal interactions with the host galaxy.

preprint2013arXiv

The role of galaxy interaction in the SFR-M relation: characterizing morphological properties of Herschel-selected galaxies at 0.2<z<1.5

Galaxy interactions/mergers have been shown to dominate the population of IR luminous galaxies (log(LIR)>11.6Lsun) in the local Universe (z<0.25). Recent studies based on the relation between galaxies' star formation rates and stellar mass (the SFR-M relation or the galaxy main sequence (MS)) have suggested that galaxy interaction/mergers may only become significant when galaxies fall well above the galaxy MS. Since the typical SFR at given M increases with redshift, the existence of galaxy MS implies that massive, IR-luminous galaxies at high-z may not necessarily be driven by galaxy interactions. We examine the role of galaxy interactions in the SFR-M relation by carrying out a morphological analysis of 2084 Herschel-selected galaxies at 0.2 < z < 1.5 in the COSMOS field. Herschel-PACS and -SPIRE observations covering the full 2-deg^2 COSMOS field provide one of the largest far-IR selected samples of high-redshift galaxies with well-determined redshifts to date, with sufficient sensitivity at z ~ 1, to sample objects lying on and above the galaxy MS. Using a detailed visual classification scheme, we show that the fraction of "disk galaxies" decreases and the fraction of "irregular" galaxies increases systematically with increasing LIR out to z ~ 1.5 and z ~ 1.0, respectively. At log(LIR) > 11.5 Lsun, >50% of the objects show evident features of strongly interacting/merger systems, where this percentage is similar to the studies of local IR-luminous galaxies. The fraction of interacting/merger systems also systematically increases with the deviation from the SFR-M relation, supporting the view that galaxies fall above the MS are more dominated by mergers than the MS galaxies. Meanwhile, we find that ~18% of massive IR-luminous MS galaxies are classified as interacting systems, where this population may not evolve through the evolutionary track predicted by a simple gas exhaustion model.

preprint2012arXiv

Faint Submillimter Galaxy Counts at 450 micron

We present the results of SCUBA2 observations at 450 micron and 850 micron of the field lensed by the massive cluster A370. With a total survey area > 100 arcmin2 and 1 sigma sensitivities of 3.92 and 0.82 mJy/beam at 450 and 850 micron respectively, we find a secure sample of 20 sources at 450 micron and 26 sources at 850 micron with a signal-to-noise ratio > 4. Using the latest lensing model of A370 and Monte Carlo simulations, we derive the number counts at both wavelengths. The 450 micron number counts probe a factor of four deeper than the counts recently obtained from the Herschel Space Telescope at similar wavelengths, and we estimate that ~47-61% of the 450 micron extragalactic background light (EBL) resolved into individual sources with 450 micron fluxes greater than 4.5 mJy. The faint 450 micron sources in the 4 sigma sample have positional accuracies of 3 arcseconds, while brighter sources (signal-to-noise > 6 sigma) are good to 1.4 arcseconds. Using the deep radio map (1 sigma ~ 6 uJy) we find that the percentage of submillimeter sources having secure radio counterparts is 85% for 450 micron sources with intrinsic fluxes > 6 mJy and 67% for 850 micron sources with intrinsic fluxes > 4 mJy. We also find that 67% of the > 4 sigma 450 micron sources are detected at 850 micron, while the recovery rate at 450 micron of > 4 sigma 850 micron sources is 54%. Combined with the source redshifts estimated using millimetric flux ratios, the recovered rate is consistent with the scenario where both 450 micron and 20 cm emission preferentially select lower redshift dusty sources, while 850 micron emission traces a higher fraction of dusty sources at higher redshifts. [Abridge]

preprint2012arXiv

Far-infrared Spectral Energy Distribution Fitting for Galaxies Near and Far

Spectral Energy Distribution (SED) fitting in the far-infrared (FIR) is greatly limited by a dearth of data and an excess of free parameters - from galaxies' dust composition, temperature, mass, orientation, opacity, to heating from AGN. This paper presents a simple FIR SED fitting technique joining a modified, single dust temperature greybody, representing the reprocessed starburst emission in the whole galaxy, to a mid-infrared powerlaw, which approximates hot-dust emission from AGN heating or clumpy, hot starbursting regions. This FIR SED can be used to measure infrared luminosities, dust temperatures and dust masses for both local and high-z galaxies with 3 to 10+ FIR photometric measurements. This fitting method is compared to infrared template SEDs in the literature using photometric data on 65 local luminous and ultraluminous infrared galaxies, (U)LIRGs. Despite relying only on 2-4 free parameters, the coupled greybody/powerlaw SED fitting described here produces better fits to photometric measurements than best-fit literature template SEDs (with residuals a factor of ~2 lower). A mean emissivity index of beta=1.60+-0.38 and mid-infrared powerlaw slope of alpha=2.0+-0.5 is measured; the former agrees with the widely presumed emissivity index of beta=1.5 and the latter is indicative of an optically-thin dust medium with a shallow radial density profile, ~r^-0.5. Adopting characteristic dust temperature as the inverse wavelength where the SED peaks, dust temperatures ~25-45K are measured for local (U)LIRGs, ~5-15K colder than previous estimates using only simple greybodies. This comparative study highlights the impact of SED fitting assumptions on the measurement of physical properties such as infrared luminosity (and thereby infrared-based star formation rate), dust temperature and dust mass, for both local and high-redshift galaxies. [abridged]

preprint2012arXiv

Spectral Energy Distributions of Local Luminous And Ultraluminous Infrared Galaxies

Luminous and ultraluminous infrared galaxies ((U)LIRGs) are the most extreme star forming galaxies in the universe. The local (U)LIRGs provide a unique opportunity to study their multi-wavelength properties in detail for comparison to their more numerous counterparts at high redshifts. We present common large aperture photometry at radio through X-ray wavelengths, and spectral energy distributions (SEDs) for a sample of 53 nearby LIRGs and 11 ULIRGs spanning log (LIR/Lsun) = 11.14-12.57 from the flux-limited Great Observatories All-sky LIRG Survey (GOALS). The SEDs for all objects are similar in that they show a broad, thermal stellar peak and a dominant FIR thermal dust peak, where nuLnu(60um) / nuLnu(V) increases from ~2-30 with increasing LIR. When normalized at IRAS-60um, the largest range in the luminosity ratio, R(lambda)=log[nuLnu(lambda)/nuLnu(60um)] observed over the full sample is seen in the Hard X-rays (HX=2-10 keV). A small range is found in the Radio (1.4GHz), where the mean ratio is largest. Total infrared luminosities, LIR(8-1000um), dust temperatures, and dust masses were computed from fitting thermal dust emission modified blackbodies to the mid-infrared (MIR) through submillimeter SEDs. The new results reflect an overall ~0.02 dex lower luminosity than the original IRAS values. Total stellar masses were computed by fitting stellar population synthesis models to the observed near-infrared (NIR) through ultraviolet (UV) SEDs. Mean stellar masses are found to be log(M/Msun) = 10.79+/-0.40. Star formation rates have been determined from the infrared (SFR_IR~45Msun/yr) and from the monochromatic UV luminosities (SFR_UV~1.3Msun/yr), respectively. Multiwavelength AGN indicators have be used to select putative AGN: about 60% of the ULIRGs would have been classified as an AGN by at least one of the selection criteria.

preprint2009arXiv

Submillimetre Detection of the z=2.83 Lyman-Break Galaxy, Westphal-MM8, and Implications for SCUBA2

We present confusion limited submillimetre (submm) observations with the SCUBA camera on the JCMT of the z=2.83 Lyman-break galaxy (LBG), Westphal-MM8, reaching an 850um sensitivity even greater than that achieved in the SCUBA map of the Hubble Deep Field region. The detection of Westphal-MM8 (S_{850um} = 1.98+-0.48mJy), along with the literature submm detections of lensed LBGs, suggest that the LBG population may contribute significantly to the source counts of submm selected galaxies in the 1-2mJy regime. Additionally, submm-luminous LBGs are a viable progenitor population for the recently discovered evolved galaxies at z~2-2.5. These observations represent an important baseline for SCUBA2 observations which will regularly map large regions of the sky to this depth.