We have developed a new Bayesian method to correct the flux densities of astronomical sources. The hybrid method combines a simulated likelihood to model survey selection together with an analytic ...source-count-based prior. The simulated likelihood captures the effect of complicated selection methods, such as multi-frequency filtering or imposed restrictions on recovered sample properties (e.g., color cuts). Simulations are also able to capture unanticipated sources of uncertainty. In this way, the method enables a broader application of Bayesian techniques. Use of an analytic prior allows variation of assumed source count models without re-simulating the likelihood. We present the method along with a detailed description of an application to real survey data from the Atacama Cosmology Telescope.
Abstract
We report on the serendipitous discovery of three transient millimeter-wave sources using data from the Atacama Cosmology Telescope. The first, detected at R.A
.
= 273.8138, decl. = −49.4628 ...at ∼50
σ
total, brightened from less than 5 mJy to at least 1100 mJy at 150 GHz with an unknown rise time shorter than 13 days, during which the increase from 250 mJy to 1100 mJy took only 8 minutes. Maximum flux was observed on 2019 November 8. The source’s spectral index in flux between 90–150 GHz was positive,
α
= 1.5 ± 0.2. The second, detected at R.A. = 105.1584, decl
.
= −11.2434 at ∼20
σ
total, brightened from less than 20 mJy to at least 300 mJy at 150 GHz with an unknown rise time shorter than 8 days. Maximum flux was observed on 2019 December 15. Its spectral index was also positive,
α
= 1.8 ± 0.2. The third, detected at R.A
.
= 301.9952, decl. = 16.1652 at ∼40
σ
total, brightened from less than 8 mJy to at least 300 mJy at 150 GHz over a day or less but decayed over a few days. Maximum flux was observed on 2018 September 11. Its spectrum was approximately flat, with a spectral index of
α
= −0.2 ± 0.1. None of the sources were polarized to the limits of these measurements. The two rising-spectrum sources are coincident in position with M and K stars, while the third is coincident with a G star.
Using a radio-quiet subsample of the Sloan Digital Sky Survey spectroscopic quasar catalogue, spanning redshifts 0.5–3.5, we derive the mean millimetre and far-infrared quasar spectral energy ...distributions (SEDs) via a stacking analysis of Atacama Cosmology Telescope and Herschel-Spectral and Photometric Imaging REceiver data. We constrain the form of the far-infrared emission and find 3σ–4σ evidence for the thermal Sunyaev–Zel'dovich (SZ) effect, characteristic of a hot ionized gas component with thermal energy (6.2 ± 1.7) × 1060 erg. This amount of thermal energy is greater than expected assuming only hot gas in virial equilibrium with the dark matter haloes of (1–5) × 1012 h
−1 M⊙ that these systems are expected to occupy, though the highest quasar mass estimates found in the literature could explain a large fraction of this energy. Our measurements are consistent with quasars depositing up to
$(14.5 \pm 3.3)\tau _8^{-1}$
per cent of their radiative energy into their circumgalactic environment if their typical period of quasar activity is τ8 × 108 yr. For high quasar host masses, ∼1013 h
−1 M⊙, this percentage will be reduced. Furthermore, the uncertainty on this percentage is only statistical and additional systematic uncertainties enter at the 40 per cent level. The SEDs are dust dominated in all bands and we consider various models for dust emission. While sufficiently complex dust models can obviate the SZ effect, the SZ interpretation remains favoured at the 3σ–4σ level for most models.
We present a catalog of 510 radio-loud active galactic nuclei (AGNs, primarily blazars) and 287 dusty star-forming galaxies (DSFGs) detected by the Atacama Cosmology Telescope at significance in ...frequency bands centered on 148 GHz (2 mm), 218 GHz (1.4 mm), and 277 GHz (1.1 mm), from a 480 deg2 strip centered at R.A. 00h on the celestial equator with additional 360 deg2 shallower auxiliary fields at other longitudes. The combination of the deepest available 218 GHz wide-field imaging, our 277 GHz data, and multiband filtering results in the most sensitive wide-field millimeter-wave DSFG selection to date, with rms noise level referenced to 218 GHz reaching below 2 mJy. We have developed new techniques to remove Galactic contamination (including evidence for CO (2−1) line emission) from the extragalactic catalog, yielding a catalog of 321 Galactic sources in addition to the extragalactic catalog. We employ a new flux debiasing method that accounts for the heterogeneous sample selection in the presence of Galactic cuts. We present the spectral properties and source counts of the AGNs and DSFGs. The DSFG spectra depart from the Rayleigh-Jeans regime of an optically thin modified blackbody between 218 and 277 GHz, consistent with optically thick emission or an additional cold dust component. For AGNs with 148 and 218 GHz flux density >50 mJy, we estimate the interyear rms fractional deviation in flux density due to source variability to be 40% with a 0.98 interband correlation coefficient. We provide source counts for AGNs in the range of 8-2870 mJy and for DSFGs in the range of 8-90 mJy. Our DSFG counts probe both the brighter, lensed population and the fainter, unlensed population. At 277 GHz we report the first measurements of source counts at these flux densities, finding an excess above most model count predictions. Finally, we present 30 of the brightest DSFGs that were selected for multifrequency study as candidate high-z lensed systems.
(ProQuest: ... denotes formulae and/or non-USASCII text omitted) We describe the measurement of the beam profiles and window functions for the Atacama Cosmology Telescope (ACT), which operated from ...2007 to 2010 with kilopixel bolometer arrays centered at 148, 218, and 277 GHz. Maps of Saturn are used to measure the beam shape in each array and for each season of observations. Radial profiles are transformed to Fourier space in a way that preserves the spatial correlations in the beam uncertainty to derive window functions relevant for angular power spectrum analysis. Several corrections are applied to the resulting beam transforms, including an empirical correction measured from the final cosmic microwave background (CMB) survey maps to account for the effects of mild pointing variation and alignment errors. Observations of Uranus made regularly throughout each observing season are used to measure the effects of atmospheric opacity and to monitor deviations in telescope focus over the season. Using the WMAP-based calibration of the ACT maps to the CMB blackbody, we obtain precise measurements of the brightness temperatures of the Uranus and Saturn disks at effective frequencies of 149 and 219 GHz. For Uranus we obtain thermodynamic brightness temperatures ... = 106.7 + or - 2.2 K and ... = 100.1+ or -3.1 K. For Saturn, we model the effects of the ring opacity and emission using a simple model and obtain resulting (unobscured) disk temperatures of ... = 137.3 + or - 3.2 K and ... = 137.3 + or - 4.7 K.
We present a comprehensive analysis of the rest-frame UV to near-IR spectral energy distributions (SEDs) and rest-frame optical spectra of four of the brightest gravitationally lensed galaxies in the ...literature: RCSGA 032727-132609 at z = 1.70, MS1512-cB58 at z = 2.73, SGAS J152745.1+065219 at z = 2.76, and SGAS J122651.3+215220 at z = 2.92. This includes new Spitzer imaging for RCSGA0327 as well as new spectra, near-IR imaging and Spitzer imaging for SGAS1527 and SGAS1226. Lensing magnifications of 3-4 mag allow a detailed study of the stellar populations and physical conditions. We compare star formation rates (SFRs) as measured from the SED fit, the H alpha and OII lambda3727 emission lines, and the UV+IR bolometric luminosity where 24 mu m photometry is available. The SFR estimate from the SED fit is consistently higher than the other indicators, which suggests that the Calzetti dust extinction law used in the SED fitting is too flat for young star-forming galaxies at z ~ 2. Our analysis finds similar stellar population parameters for all four lensed galaxies: stellar masses (3-7) x 10 super(9) M sub(middot in circle), young ages ~400 Myr, little dust content E(B - V) = 0.10-0.25, and SFRs around 20-100 M sub(middot in circle) yr super(-1). Compared to typical values for the galaxy population at z ~ 2, this suggests we are looking at newly formed, starbursting systems that have only recently started the buildup of stellar mass. These results constitute the first detailed, uniform analysis of a sample of the growing number of strongly lensed galaxies known at z ~ 2.
We report Northern Extended Millimeter Array CO(J = 3 − 2) observations of the dusty star-forming galaxy ACT-S J020941+001557 at z = 2.5528, which was detected as an unresolved source in the Atacama ...Cosmology Telescope (ACT) equatorial survey. Our spatially resolved spectral line data support the derivation of a gravitational lens model from 37 independent velocity channel maps using a pixel-based algorithm, from which we infer a velocity-dependent magnification factor 7-22 with a luminosity-weighted mean . The resulting source-plane reconstruction is consistent with a rotating disk, although other scenarios cannot be ruled out by our data. After correction for lensing, we derive a line luminosity , a cold gas mass , a dynamical mass , and a gas mass fraction . The line brightness temperature ratio of r3,1 1.6 relative to a Green Bank Telescope CO(J = 1 − 0) detection may be elevated by a combination of external heating of molecular clouds, differential lensing, and/or pointing errors.
We present the first weak-lensing-based scaling relation between galaxy cluster mass, M sub(WL), and integrated Compton parameter Y sub(sph). Observations of 18 galaxy clusters at z Asymptotically = ...to 0.2 were obtained with the Subaru 8.2 m telescope and the Sunyaev-Zel'dovich Array. The M sub(WL)-Y sub(sph) scaling relations, measured at Delta = 500, 1000, and 2500 rho sub(c), are consistent in slope and normalization with previous results derived under the assumption of hydrostatic equilibrium (HSE). We find an intrinsic scatter in M sub(WL) at fixed Y sub(sph) of 20%, larger than both previous measurements of M sub(HSE)-Y sub(sph) scatter as well as the scatter in true mass at fixed Y sub(sph) found in simulations. Moreover, the scatter in our lensing-based scaling relations is morphology dependent, with 30%-40% larger M sub(WL) for undisturbed compared to disturbed clusters at the same Y sub(sph) at r sub(500). Further examination suggests that the segregation may be explained by the inability of our spherical lens models to faithfully describe the three-dimensional structure of the clusters, in particular, the structure along the line of sight. We find that the ellipticity of the brightest cluster galaxy, a proxy for halo orientation, correlates well with the offset in mass from the mean scaling relation, which supports this picture. This provides empirical evidence that line-of-sight projection effects are an important systematic uncertainty in lensing-based scaling relations.
We conduct a statistical analysis of the radio source population in galaxy clusters as a function of redshift by matching radio sources from the Faint Images of the Radio Sky at Twenty-Centimeters ...catalog with 618 optically selected galaxy clusters from the first Red-Sequence Cluster Survey (RCS1). The total number of excess radio source components (above the background level) per cluster is 0.14 ? 0.02 for clusters with 0.35 < z < 0.65 and is 0.10 ? 0.02 for clusters with 0.65 < z < 0.95. The richest clusters in the sample have more radio sources than clusters with low or intermediate richness. When we divide our sample into bins according to cluster richness, we do not observe any significant difference (>1.5 Delta *s) in the number of radio sources per unit of cluster mass for the galaxy clusters with 0.35 < z < 0.65 as compared to the galaxy clusters with 0.65 < z < 0.95. Thus, the entire sample can be characterized by the number of (L 1.4 GHz > 4.1 X 1024 W Hz--1) radio sources per unit (1014 M ) mass, which we measure to be 0.031 ? 0.004. We further characterize the population of galaxy cluster-related radio sources through visual inspection of the RCS1 images, finding that although the radio activity of brightest cluster galaxies (BCGs) also does not strongly evolve between our high- and low-redshift samples, the lower-redshift, richest clusters are more likely to host radio-loud BCGs than the higher-redshift, richest clusters or poorer clusters at the 2 Delta *s level. We also find evidence that the number of radio-loud cluster members that are not BCGs per cluster per unit mass is higher at high redshift than at low redshift, at the 3 Delta *s level.