ABSTRACT
We demonstrate the use of Haralick features for the automated classification of radio galaxies. The set of thirteen Haralick features represent an extremely compact non-parametric ...representation of image texture, and are calculated directly from imagery using the Grey Level Co-occurrence Matrix (GLCM). The GLCM is an encoding of the relationship between the intensity of neighbouring pixels in an image. Using 10 000 sources detected in the first data release of the LOFAR Two-metre Sky Survey (LoTSS), we demonstrate that Haralick features are highly efficient, rotationally invariant descriptors of radio galaxy morphology. After calculating Haralick features for LoTSS sources, we employ the fast density-based hierarchical clustering algorithm hdbscan to group radio sources into a sequence of morphological classes, illustrating a simple methodology to classify and label new, unseen galaxies in large samples. By adopting a ‘soft’ clustering approach, we can assign each galaxy a probability of belonging to a given cluster, allowing for more flexibility in the selection of galaxies according to combinations of morphological characteristics and for easily identifying outliers: those objects with a low probability of belonging to any cluster in the Haralick space. Although our demonstration focuses on radio galaxies, Haralick features can be calculated for any image, making this approach also relevant to large optical imaging galaxy surveys.
ABSTRACT
We present AstroVaDEr (Astronomical Variational Deep Embedder), a variational autoencoder designed to perform unsupervised clustering and synthetic image generation using astronomical ...imaging catalogues. The model is a convolutional neural network that learns to embed images into a low-dimensional latent space, and simultaneously optimizes a Gaussian Mixture Model (GMM) on the embedded vectors to cluster the training data. By utilizing variational inference, we are able to use the learned GMM as a statistical prior on the latent space to facilitate random sampling and generation of synthetic images. We demonstrate AstroVaDEr’s capabilities by training it on grey-scaled gri images from the Sloan Digital Sky Survey, using a sample of galaxies that are classified by Galaxy Zoo 2. An unsupervised clustering model is found that separates galaxies based on learned morphological features such as axial ratio, surface brightness profile, orientation, and the presence of companions. We use the learned mixture model to generate synthetic images of galaxies based on the morphological profiles of the Gaussian components. AstroVaDEr succeeds in producing a morphological classification scheme from unlabelled data, but unexpectedly places high importance on the presence of companion objects – demonstrating the importance of human interpretation. The network is scalable and flexible, allowing for larger data sets to be classified, or different kinds of imaging data. We also demonstrate the generative properties of the model, which allow for realistic synthetic images of galaxies to be sampled from the learned classification scheme. These can be used to create synthetic image catalogues or to perform image processing tasks such as deblending.
ABSTRACT
Matching the number counts of high-z submillimetre-selected galaxies (SMGs) has been a long-standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to ...model the submm emission from galaxies in the simba cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850 μm-selected sources. We find good agreement with the shape of the integrated 850 μm luminosity function, and the normalization is within 0.25 dex at >3 mJy, unprecedented for a fully cosmological hydrodynamic simulation, along with good agreement in the redshift distribution of bright SMGs. The agreement is driven primarily by simba’s good match to infrared measures of the star formation rate (SFR) function between z = 2 and 4 at high SFRs. Also important is the self-consistent on-the-fly dust model in simba, which predicts, on average, higher dust masses (by up to a factor of 2.5) compared to using a fixed dust-to-metals ratio of 0.3. We construct a light-cone to investigate the effect of far-field blending, and find that 52 per cent of sources are blends of multiple components, which makes a small contribution to the normalization of the bright end of the number counts. We provide new fits to the 850 μm luminosity as a function of SFR and dust mass. Our results demonstrate that solutions to the discrepancy between submm counts in simulations and observations, such as a top-heavy initial mass function, are unnecessary, and that submillimetre-bright phases are a natural consequence of massive galaxy evolution.
ABSTRACT
We demonstrate how galaxy morphologies can be represented by weighted sums of ‘eigengalaxies’ and how eigengalaxies can be used in a probabilistic framework to enable principled and ...simplified approaches in a variety of applications. Eigengalaxies can be derived from a Principal Component Analysis (PCA) of sets of single- or multiband images. They encode the image space equivalent of basis vectors that can be combined to describe the structural properties of large samples of galaxies in a massively reduced manner. As an illustration, we show how a sample of 10243 galaxies in the Hubble Space Telescope CANDELS survey can be represented by just 12 eigengalaxies. We show in some detail how this image space may be derived and tested. We also describe a probabilistic extension to PCA (PPCA) which enables the eigengalaxy framework to assign probabilities to galaxies. We present four practical applications of the probabilistic eigengalaxy framework that are particularly relevant for the next generation of large imaging surveys: we (i) show how low likelihood galaxies make for natural candidates for outlier detection; (ii) demonstrate how missing data can be predicted; (iii) show how a similarity search can be performed on exemplars; (iv) demonstrate how unsupervised clustering of objects can be implemented.
Abstract
Obscuration in quasars may arise from steep viewing angles along the dusty torus, or instead may represent a distinct phase of supermassive black hole growth. We test these scenarios by ...probing the host dark matter halo environments of ∼1.4 million Wide-field Infrared Survey Explorer-selected obscured and unobscured quasars at 〈
z
〉 = 1.4 using angular clustering measurements as well as cross-correlation measurements of quasar positions with the gravitational lensing of the cosmic microwave background. We interpret these signals within a halo occupation distribution framework to conclude that obscured systems reside in more massive effective halos (∼10
12.9
h
−1
M
⊙
) than their unobscured counterparts (∼10
12.6
h
−1
M
⊙
), though we do not detect a difference in the satellite fraction. We find excellent agreement between the clustering and lensing analyses and show that this implies the observed difference is robust to uncertainties in the obscured quasar redshift distribution, highlighting the power of combining angular clustering and weak lensing measurements. This finding appears in tension with models that ascribe obscuration exclusively to orientation of the dusty torus along the line of sight, and instead may be consistent with the notion that some obscured quasars are attenuated by galaxy-scale or circumnuclear material during an evolutionary phase.
Abstract
We have used the Submillimeter Array (SMA) at 860 $\mu$m to observe the brightest sources in the Submillimeter Common User Bolometer Array-2 (SCUBA-2) Cosmology Legacy Survey (S2CLS). The ...goal of this survey is to exploit the large field of the S2CLS along with the resolution and sensitivity of the SMA to construct a large sample of these rare sources and to study their statistical properties. We have targeted 70 of the brightest single-dish SCUBA-2 850 $\mu$m sources down to S850 ≈ 8 mJy, achieving an average synthesized beam of 2.4 arcsec and an average rms of σ860 = 1.5 mJy beam−1 in our primary beam-corrected maps. We searched our SMA maps for 4σ peaks, corresponding to S860 ≳ 6 mJy sources, and detected 62, galaxies, including three pairs. We include in our study 35 archival observations, bringing our sample size to 105 bright single-dish submillimetre sources with interferometric follow-up. We compute the cumulative and differential number counts, finding them to overlap with previous single-dish survey number counts within the uncertainties, although our cumulative number count is systematically lower than the parent S2CLS cumulative number count by 14 ± 6 per cent between 11 and 15 mJy. We estimate the probability that a ≳10 mJy single-dish submillimetre source resolves into two or more galaxies with similar flux densities to be less than 15 per cent. Assuming the remaining 85 per cent of the targets are ultraluminous starburst galaxies between z = 2 and 3, we find a likely volume density of ≳400 M⊙ yr−1 sources to be ${\sim }\,3^{+0.7}_{-0.6}\,{\times }\,10^{-7}$ Mpc−3. We show that the descendants of these galaxies could be ≳4 × 1011 M⊙ local quiescent galaxies, and that about 10 per cent of their total stellar mass would have formed during these short bursts of star formation.
Radio emission from radio-quiet quasars may be due to star formation in the quasar host galaxy, to a jet launched by the supermassive black hole, or to relativistic particles accelerated in a ...wide-angle radiatively driven outflow. In this paper, we examine whether radio emission from radio-quiet quasars is a byproduct of star formation in their hosts. To this end, we use infrared spectroscopy and photometry from Spitzer and Herschel to estimate or place upper limits on star formation rates in hosts of ∼300 obscured and unobscured quasars at z < 1. We find that low-ionization forbidden emission lines such as Ne ii and Ne iii are likely dominated by quasar ionization and do not provide reliable star formation diagnostics in quasar hosts, while polycyclic aromatic hydrocarbon (PAH) emission features may be suppressed due to the destruction of PAH molecules by the quasar radiation field. While the bolometric luminosities of our sources are dominated by the quasars, the 160 μm fluxes are likely dominated by star formation, but they too should be used with caution. We estimate median star formation rates to be 6–29 M⊙ yr−1, with obscured quasars at the high end of this range. This star formation rate is insufficient to explain the observed radio emission from quasars by an order of magnitude, with log (L
radio, obs/L
radio, SF) = 0.6–1.3 depending on quasar type and star formation estimator. Although radio-quiet quasars in our sample lie close to the 8–1000 μm infrared/radio correlation characteristic of the star-forming galaxies, both their infrared emission and their radio emission are dominated by the quasar activity, not by the host galaxy.
This paper presents new deep and wide narrow-band surveys undertaken with United Kingdom Infrared Telescope (UKIRT), Subaru and the Very Large Telescope (VLT), a unique combined effort to select ...large, robust samples of Hα star-forming galaxies at z = 0.40, 0.84, 1.47 and 2.23 (corresponding to look-back times of 4.2, 7.0, 9.2 and 10.6 Gyr) in a uniform manner over ∼2 deg2 in the Cosmological Evolution Survey and Ultra Deep Survey fields. The deep multi-epoch Hα surveys reach a matched 3σ flux limit of 3 M yr−1 out to z = 2.2 for the first time, while the wide area and the coverage over two independent fields allow us to greatly overcome cosmic variance and assemble by far the largest samples of Hα emitters. Catalogues are presented for a total of 1742, 637, 515 and 807 Hα emitters, robustly selected at z = 0.40, 0.84, 1.47 and 2.23, respectively, and used to determine the Hα luminosity function and its evolution. The faint-end slope of the Hα luminosity function is found to be α = −1.60 ± 0.08 over z = 0-2.23, showing no significant evolution. The characteristic luminosity of star-forming galaxies, L*Hα, evolves significantly as log L*Hα(z) = 0.45z + log L*
z = 0. This is the first time Hα has been used to trace star formation activity with a single homogeneous survey at z = 0.4-2.23. Overall, the evolution seen with Hα is in good agreement with the evolution seen using inhomogeneous compilations of other tracers of star formation, such as far-infrared and ultraviolet, jointly pointing towards the bulk of the evolution in the last 11 Gyr being driven by a statistically similar star-forming population across cosmic time, but with a strong luminosity increase from z ∼ 0 to ∼2.2. Our uniform analysis allows us to derive the Hα star formation history (SFRH) of the Universe, showing a clear rise up to z ∼ 2.2, for which the simple parametrization log10ρSFR = −2.1(1 + z)−1 is valid over 80 per cent of the age of the Universe. The results reveal that both the shape and normalization of the Hα SFRH are consistent with the measurements of the stellar mass density growth, confirming that our Hα SFRH is tracing the bulk of the formation of stars in the Universe for z < 2.23. The star formation activity over the last ∼11 Gyr is responsible for producing ∼95 per cent of the total stellar mass density observed locally, with half of that being assembled in 2 Gyr between z = 1.2 and 2.2, and the other half in 8 Gyr (since z < 1.2). If the star formation rate density continues to decline with time in the same way as seen in the past ∼11 Gyr, then the stellar mass density of the Universe will reach a maximum which is only 5 per cent higher than the present-day value.