Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Galaxy clusters”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

Formation of S0s in extreme environments Ⅲ: the role of environment in the formation pathways

It is well established that there are at least two main channels to form lenticular (or S0) galaxies. The first, which we name ‘faded spiral’ scenario, includes quenching events that led to consumption or removal of gas from a spiral progenitor. The second, which we call ‘merger’ scenario, includes merger-like events and interactions between galaxies. Each scenario leaves characteristic signatures in the newly formed lenticular galaxy. However, the conditions that trigger one mechanism over another are still unknown. This paper is the third of a series aimed at understanding the role of the environment in the formation of lenticular galaxies. In this study, we combine the kinematics, morphology, and properties of the stellar populations of 329 S0s from the SAMI and MaNGA surveys in order to highlight the role of the environment in the process. We divide the S0s into two classes (A and B) according to their global properties, that we can associate to the products of a faded spiral scenario (class A) or a merger scenario (class B). Here we then study how the various classes are distributed within different environments. Our study reveals that the ‘faded spiral’ pathway is the most efficient channel to produce S0s, and it becomes more efficient as the mass of the group or cluster or local density of galaxies increase. The merger pathway is also a viable channel, and its efficiency becomes higher with decreasing local density or environment mass.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterization of DESI fiber assignment incompleteness effect on 2-point clustering and mitigation methods for DR1 analysis

We present an in-depth analysis of the fiber assignment incompleteness in the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1). This incompleteness is caused by the restricted mobility of the robotic fiber positioner in the DESI focal plane, which limits the number of galaxies that can be observed at the same time, especially at small angular separations. As a result, the observed clustering amplitude is suppressed in a scale-dependent manner, which, if not addressed, can severely impact the inference of cosmological parameters. We discuss the methods adopted for simulating fiber assignment on mocks and data. In particular, we introduce the fast fiber assignment (FFA) emulator, which was employed to obtain the power spectrum covariance adopted for the DR1 full-shape analysis. We present the mitigation techniques, organised in two classes: measurement stage and model stage. We then use high fidelity mocks as a reference to quantify both the accuracy of the FFA emulator and the effectiveness of the different measurement-stage mitigation techniques. This complements the studies conducted in a parallel paper for the model-stage techniques, namely the θ-cut approach. We find that pairwise inverse probability (PIP) weights with angular upweighting recover the “true” clustering in all the cases considered, in both Fourier and configuration space. Notably, we present the first ever power spectrum measurement with PIP weights from real data.

cosmological simulations↗

The cross-correlation of galaxies in absorption with the Lyman α forest

ABSTRACT We present the first clustering measurement of Strong Blended Lyman α (SBLA) absorption systems by measuring their cross-correlation with the Lyman α forest. SBLAs are a new population of absorbers detected within the Lyman α forest. We find a bias of 2.329 ± 0.057, consistent with that of Damped Lyman α absorbers (DLAs). For DLAs, we recover a bias of 2.331 ± 0.057 larger than previously reported. We also find a redshift space distortion parameter β = 0.417 ± 0.010, also consistent with the recovered value for DLAs (β = 0.416 ± 0.010). This is consistent with SBLA and DLA systems tracing different portions of the circumgalactic medium of a broadly common population of galaxies. Given these common clustering properties, we combined them to perform a cross-correlation of galaxies in absorption with the Lyα forest. We find that the BAO scale uncertainty of this new measurement is 1.75 × that of Lyα autocorrelation and 1.6 × that of the quasar cross-correlation with the Lyα forest. We note that the current preferred metal contamination model for fitting the correlation functions with respect to the Lyα forest is not realistic enough for SBLA systems, likely due to their status as high redshift precision sites of high metal enrichment. Mock spectra including SBLA systems and their associated metal absorption are required to understand this sample fully. We conclude that SBLAs have the potential to complement the standard Lyα cosmological analyses in future surveys.

Astronomy & Astrophysics↗

Using Gravitational Lensing to Probe for Candidate Quintessential Galaxies in the Epoch of Reionization

Abstract Understanding the properties of the first generation of galaxies is an ongoing challenge in observational astrophysics. James Webb Space Telescope deep-field observations have led to the identification of large numbers of galaxies within the Epoch of Reionization (EoR); however, it remains observationally challenging to study the detailed properties of the most numerous, faint high- z galaxies without strong gravitational lensing. To this end, we present five m AB ≃ 24 galaxies, all projected behind the cores of massive clusters where high lensing magnifications significantly (≳1.5 mag) boost their observed brightness. These galaxies are each characterized by a drop-off in their spectra between the near-IR and optical wavelengths, suggesting the presence of a Lyman break at z > 6. The bright sources enable precise dropout color constraints. We combine the dropout data with follow-up imaging to constrain continuum slopes redward of their breaks. Using the open-source spectral energy distribution modeling software Prospector , we characterize two of these galaxies as probable high- z ( z ∼ 6.5–7) and three as probable low- z ( z ∼ 2) sources, demonstrating that for very bright ( m AB < 24) dropout candidates even limited photometric data have good power to distinguish between high- z galaxies and red or dusty low- z galaxies. The two candidate high- z dropouts are comparably bright to the brightest galaxies identified at such redshifts in deep fields, but due to lensing magnifications are intrinsically fainter by at least 1.5 mag. If confirmed, these galaxies represent exciting laboratories for performing detailed follow-up studies of galaxy properties within the EoR.

Roberson, Joshua [University of Cincinnati] (ORCID↗

Preliminary clustering properties of the DESI BGS bright targets using DR9 Legacy Imaging Surveys

ABSTRACT We characterize the selection cuts and clustering properties of a magnitude-limited sample of bright galaxies that is part of the Bright Galaxy Survey (BGS) of the Dark Energy Spectroscopic Instrument (DESI) using the ninth data release of the Legacy Imaging Surveys (DR9). We describe changes in the DR9 selection compared to the DR8 one and we also compare the DR9 selection in three distinct regions: BASS/MzLS in the north Galactic Cap (NGC), DECaLS in the NGC, and DECaLS in the south Galactic Cap (SGC). We investigate the systematics associated with the selection and assess its completeness by matching the BGS targets with the Galaxy and Mass Assembly (GAMA) survey. We measure the angular clustering for the overall bright sample (rmag ≤ 19.5) and as function of apparent magnitude and colour. This enables to determine the clustering strength r0 and slope γ by fitting a power-law model that can be used to generate accurate mock catalogues for this tracer. We use a counts-in-cells technique to explore higher order statistics and cross-correlations with external spectroscopic data sets in order to check the evolution of the clustering with redshift and the redshift distribution of the BGS targets using clustering redshifts. While this work validates the properties of the BGS bright targets, the final target selection pipeline and clustering properties of the entire DESI BGS will be fully characterized and validated with the spectroscopic data of Survey Validation.

79 ASTRONOMY AND ASTROPHYSICS↗

The DESI DR1 Peculiar Velocity Survey: Mock Catalog

We describe the production of the official set of mock catalogs for the Dark Energy Spectroscopic Instrument Peculiar Velocity Survey (DESI-PV) Data Release 1 (DR1). Our mock catalogs reproduce the Bright Galaxy Survey number density and clustering at low redshift $(z<0.1)$ and the DESI PV samples of Fundamental plane and Tully-Fisher distances, from which we derive peculiar velocities. We carefully match mock and data properties and we mimic measurements of distance indicators and peculiar velocities, which follow the same statistical properties as real data. Mock samples of type-Ia supernovae also complement the other two distance indicators. Our 675 available mock realizations were used consistently by our three different methodologies described in our companion papers that measure the growth rate of structure $fσ_8$ with DESI PV DR1. Those mocks allow us to perform precise tests of clustering models and uncertainty estimation to an unprecedented level of accuracy, and compute correlations between methodologies. The consensus value for the DESI DR1 PV growth rate measurement is $fσ_8 = 0.450 \pm 0.055$. This sample of mock catalogs represents the largest and most realistic set for cosmological measurements with peculiar velocities to date.

Bautista, J. [Marseille, CPPM]↗

Weak lensing mass-richness relation of redMaPPer clusters in LSST DESC DC2 simulations

Cluster scaling relations are key ingredients in cluster abundance-based cosmological studies. In optical cluster cosmology, where clusters are detected through their richness, cluster-weak gravitational lensing has proven to be a powerful tool to constrain the cluster mass-richness relation. This work is conducted as part of the Dark Energy Science Collaboration (DESC), which aims to analyze the Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory, starting in 2026. Cluster properties inferred from weak lensing, such as mass, suffer from several sources of bias. In this paper, we aim to test the impact of modeling choices and observational systematics in cluster lensing on the inference of the mass-richness relation. We constrained the mass-richness relation of 3600 clusters detected by the redMaPPer algorithm in the cosmoDC2 extragalactic mock catalog of the LSST DESC DC2 simulation, covering 440 deg 2 , using number count measurements and either stacked weak lensing profiles or mean cluster masses in several intervals of richness (20 ≤ λ ≤ 200) and redshift (0.2 ≤ z ≤ 1). We provide the first constraints on the redMaPPer cluster mass-richness relation detected in cosmoDC2. We find that for an LSST-like source galaxy density, our constraints are robust to changes in the concentration-mass relation, as well as the dark matter density profile modeling choices, when source redshifts and shapes are perfectly known. We find that photometric redshift uncertainties can introduce bias at the 1σ level, which could be mitigated by an overall correction factor fitted jointly with the scaling parameters. We find that including positive shear-richness covariance in the fit shifts the results by up to 0.5σ. Our constraints also offer a fair comparison to a fiducial mass-richness relation, obtained from matching cosmoDC2 halo masses to redMaPPer-detected cluster richness results.

galaxy clusters↗

Signatures of large-scale cold fronts in the optically-selected merging cluster HSC J085024+001536

Abstract We represent a joint X-ray, weak-lensing, and optical analysis of the optically-selected merging cluster, HSC J085024+001536, from the Subaru HSC-SSP survey. Both the member galaxy density and the weak-lensing mass map show that the cluster is composed of south-east and north-west components. The two-dimensional weak-lensing analysis shows that the south-east component is the main cluster, and the sub-cluster and main cluster mass ratio is $0.32^{+0.75}_{-0.23}$. The north-west sub-cluster is offset by ∼700 kpc from the main cluster center, and their relative line-of-sight velocity is ∼1300 km s−1 from spectroscopic redshifts of member galaxies. The X-ray emission is concentrated around the main cluster, while the gas mass fraction within a sphere of 1′ radius of the sub-cluster is only $f_{\mathrm{gas}}=4.0^{+2.3}_{-3.3}\%$, indicating that the sub-cluster gas was stripped by ram pressure. An X-ray residual image shows three arc-like excess patterns, of which two are symmetrically located at ∼550 kpc from the X-ray morphological center, and the other is close to the X-ray core. The excess close to the sub-cluster has a cold-front feature where dense-cold gas and thin-hot gas contact. The two outer excesses are tangentially elongated about ∼450–650 kpc, suggesting that the cluster is merged with a non-zero impact parameter. Overall features revealed by the multi-wavelength datasets indicate that the cluster is at the second impact or later. Since the optically-defined merger catalog is unbiased for merger boost of the intracluster medium, X-ray follow-up observations will pave the way to understand merger physics at various phases.

Tanaka, Keigo↗

Galaxy Alignments with Surrounding Structure in the Sloan Digital Sky Survey

Using data from the Sloan Digital Sky Survey Legacy Survey, we study the alignment of luminous galaxies with spectroscopic data with the surrounding larger-scale structure as defined by galaxies with only photometric data. We find that galaxies from the red sequence have a statistically significant tendency for their apparent long axes to align parallel to the projected surrounding structure. Red galaxies more luminous than the median of our sample (M r < –21.78) have a mean alignment angle $\langle$Φ$\rangle$ < 45°, indicating preferred parallel alignment, at a significance level >4.5σ on projected scales 0.1 Mpc < r p ≤ 7.5 Mpc. Fainter red galaxies have $\langle$Φ$\rangle$ < 45° at a significance level >4.3σ at scales 1 Mpc < r p < 3 Mpc. At a projected scale r p = 3.0 Mpc, the mean alignment angle decreases steadily with increasing luminosity for red galaxies with M r ≲ –22.5, reaching $\langle$Φ$\rangle$ = 40$^{°}_{.}$49 ± 0$^{°}_{.}$56 for the most luminous 1% (M r ~ –23.57). Galaxies from the blue sequence show no statistically significant tendency for their axes to align with larger-scale structure, regardless of galaxy luminosity. Galaxies in higher-density regions do not show a statistically significant difference in the mean alignment angle from galaxies in lower-density regions; this holds true for the faint blue, luminous blue, faint red, and luminous red subsets.

79 ASTRONOMY AND ASTROPHYSICS↗

Deep transfer learning for star cluster classification: I. application to the PHANGS– HST survey

ABSTRACT We present the results of a proof-of-concept experiment that demonstrates that deep learning can successfully be used for production-scale classification of compact star clusters detected in Hubble Space Telescope(HST) ultraviolet-optical imaging of nearby spiral galaxies ($D\lesssim 20\, \textrm{Mpc}$) in the Physics at High Angular Resolution in Nearby GalaxieS (PHANGS)–HST survey. Given the relatively small nature of existing, human-labelled star cluster samples, we transfer the knowledge of state-of-the-art neural network models for real-object recognition to classify star clusters candidates into four morphological classes. We perform a series of experiments to determine the dependence of classification performance on neural network architecture (ResNet18 and VGG19-BN), training data sets curated by either a single expert or three astronomers, and the size of the images used for training. We find that the overall classification accuracies are not significantly affected by these choices. The networks are used to classify star cluster candidates in the PHANGS–HST galaxy NGC 1559, which was not included in the training samples. The resulting prediction accuracies are 70 per cent, 40 per cent, 40–50 per cent, and 50–70 per cent for class 1, 2, 3 star clusters, and class 4 non-clusters, respectively. This performance is competitive with consistency achieved in previously published human and automated quantitative classification of star cluster candidate samples (70–80 per cent, 40–50 per cent, 40–50 per cent, and 60–70 per cent). The methods introduced herein lay the foundations to automate classification for star clusters at scale, and exhibit the need to prepare a standardized data set of human-labelled star cluster classifications, agreed upon by a full range of experts in the field, to further improve the performance of the networks introduced in this study.

Wei, Wei↗

Clustering redshift distribution calibration of weak lensing surveys using the DESI-DR1 spectroscopic dataset

We estimate the source redshift distribution of current weak lensing surveys by applying the clustering-based redshift calibration technique, using the galaxy redshift sample provided by the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1). We cross-correlate the Bright Galaxy Survey (BGS), Luminous Red Galaxies (LRGs) and Emission Line Galaxies (ELGs) from DESI, within the redshift range $0.1 < z < 1.6$, with overlapping tomographic source samples from the Dark Energy Survey (DES), Kilo-Degree Survey (KiDS), and Hyper Suprime-Cam (HSC) survey. Using realistic mock catalogues, we test the stability of the clustering-redshift signal to fitting scale, reference-sample choice, and the evolution of source galaxy bias, and we explicitly model and marginalise over magnification contributions, which become non-negligible at $z \gtrsim 1$ due to the depth of the DESI ELG sample. We then compare the resulting bias-weighted redshift distributions to those calibrated using self-organising map (SOM) techniques, finding agreement within uncertainties for all surveys and tomographic bins. Our results demonstrate that clustering redshifts enabled by DESI's unprecedented spectroscopic sample provides a robust, complementary, and independent constraint capable of reducing one of the dominant systematic uncertainties in weak lensing cosmology.

Ruggeri, R. [Queensland U.]↗

Milky way bulge extra-tidal star survey: MWBest

The Milky Way Bulge Extra-Tidal Star Survey (MWBest) is a spectroscopic initiative aimed at identifying stars that have been tidally stripped from globular clusters in the inner regions of the Galaxy. As part of the Visiting Faculty Program at Fermilab, we have developed theoretical tools and simulations to compare against observed signatures of extra-tidal stars surrounding both bulge globular clusters (GCs) and low-luminosity clusters. By combining the spatial distributions of simulated stripped stars with new MWBest observations, we are gaining insight into the gravitational potential of the inner Galaxy and the role that bulge GCs may play in enriching the surrounding field star population. We further explore the hypothesis that some low-luminosity bulge GCs may in fact be dense clumps of debris, stripped from more massive GCs as they migrate inward due to dynamical friction.

Kunder, Andrea [Unlisted, US, WA]↗

The DESI One-Percent survey: exploring the Halo Occupation Distribution of Emission Line Galaxies with AbacusSummit simulations

The One-Percent survey of the Dark Energy Spectroscopic Instrument collected ~ 270k emission line galaxies (ELGs) at 0.8 < z < 1.6. The high completeness of the sample allowed the clustering to be measured down to scales never probed before, 0.04 Mpc/h in rp for the projected 2-point correlation function (2PCF) and 0.17 Mpc/h in galaxy pair separation s for the 2PCF monopole and quadrupole. The most striking feature of the measurements is a strong signal at the smallest scales, below 0.2 Mpc/h in rp and 1 Mpc/h in s. We analyse these data in the halo occupation distribution framework. We consider different distributions for central galaxies, a standard power law for satellites with no condition on the presence of a central galaxy and explore several extensions of these models. For all considered models, the mean halo mass of the sample is found to be log10 $\langle$Mh$\rangle$ ~ 11.9. We obtain a satellite mean occupation function which agrees with physically motivated ELG models only if we introduce central-satellite conformity, meaning that the satellite occupation is conditioned by the presence of central galaxies of the same type. To achieve in addition a good modelling of the clustering between 0.1 and 1 Mpc/h in rp, we allow for ELG positioning outside of the halo virial radius and find 0.5% of ELGs residing in the outskirts of halos. Furthermore, the satellite velocity dispersion inside halos is found to be ~ 30% larger than that of the halo dark matter particles. These are the main findings of our work. Here we investigate assembly bias as a function of halo concentration, local density or local density anisotropies and observe no significant change in our results. We split the data sample in two redshift bins and report no significant evolution with redshift. Lastly, changing the cosmology in the modelling impacts only slightly our results.

79 ASTRONOMY AND ASTROPHYSICS↗

The One-hundred-deg 2 DECam Imaging in Narrowbands (ODIN): Survey Design and Science Goals

We describe the survey design and science goals for One-hundred-deg 2 DECam Imaging in Narrowbands (ODIN), a NOIRLab survey using the Dark Energy Camera (DECam) to obtain deep (AB ~ 25.7) narrowband images over an unprecedented area of sky. The three custom-built narrowband filters, N419, N501, and N673, have central wavelengths of 419, 501, and 673 nm and respective FWHM of 7.5, 7.6, and 10.0 nm, corresponding to Lyα at z = 2.4, 3.1, and 4.5 and cosmic times of 2.8, 2.1, and 1.4 Gyr, respectively. When combined with even deeper, public broadband data from the Hyper Suprime-Cam, DECam, and in the future, the Legacy Survey of Space and Time, the ODIN narrowband images will enable the selection of over 100,000 Lyα-emitting (LAE) galaxies at these epochs. ODIN-selected LAEs will identify protoclusters as galaxy overdensities, and the deep narrowband images enable detection of highly extended Lyα blobs (LABs). Primary science goals include measuring the clustering strength and dark matter halo connection of LAEs, LABs, and protoclusters, and their respective relationship to filaments in the cosmic web. The three epochs allow for the redshift evolution of these properties to be determined during the period known as Cosmic Noon, where star formation was at its peak. The narrowband filter wavelengths are designed to enable interloper rejection and further scientific studies by revealing [O II] and [O III] at z = 0.34, Lyα and He II 1640 at z = 3.1, and Lyman continuum plus Lyα at z = 4.5. Ancillary science includes similar studies of the lower-redshift emission-line galaxy samples and investigations of nearby starforming galaxies resolved into numerous [O III] and [S II] emitting regions.

79 ASTRONOMY AND ASTROPHYSICS↗

Environmental Quenching of Low-surface-brightness Galaxies Near Hosts from Large Magellanic Cloud to Milky Way Mass Scales

Low-surface-brightness galaxies (LSBGs) are excellent probes of quenching and other environmental processes near massive galaxies. We study an extensive sample of LSBGs near massive hosts in the local universe that are distributed across a diverse range of environments. The LSBGs with surface-brightness ${\mu }_{\mathrm{eff},{g}}\gt 24.2\,\mathrm{mag}\,{\mathrm{arcsec}}^{-2}$ are drawn from the Dark Energy Survey Year 3 catalog while the hosts with masses $9.0\lt \mathrm{log}({{ \mathcal M }}_{\star }/{M}_{\odot })\lt 11.0$ comparable to the Milky Way and the Large Magellanic Cloud are selected from the z0MGS sample. We study the projected radial density profiles of LSBGs as a function of their color and surface brightness around hosts in both the rich Fornax–Eridanus cluster environment and the low-density field. We detect an overdensity with respect to the background density, out to 2.5 times the virial radius for both hosts in the cluster environment and the isolated field galaxies. When the LSBG sample is split by g − i color or surface brightness μ eff, g , we find the LSBGs closer to their hosts are significantly redder and brighter, like their high-surface-brightness counterparts. The LSBGs form a clear “red sequence” in both the cluster and isolated environments that is visible beyond the virial radius of the hosts. This suggests preprocessing of infalling LSBGs and a quenched backsplash population around both host samples. More so, the relative prominence of the “blue cloud” feature implies that preprocessing is ongoing near the isolated hosts compared to the cluster environment where the LSBGs are already well processed.

79 ASTRONOMY AND ASTROPHYSICS↗

A comparison between ShapeFit compression and Full-Modelling method with PyBird for DESI 2024 and beyond

DESI aims to provide one of the tightest constraints on cosmological parameters by analysing the clustering of more than thirty million galaxies. However, obtaining such constraints requires special care in validating the methodology and efforts to reduce the computational time required through data compression and emulation techniques. In this work, we perform a rigorous validation of the PyBird power spectrum modelling code with both a traditional emulated Full-Modelling approach and the model-independent ShapeFit compression approach. By using cubic box simulations that accurately reproduce the clustering and precision of the DESI survey, we find that the cosmological constraints from ShapeFit and Full-Modelling are consistent with each other at the ∼ 0.5σ level for the ΛCDM model. Both ShapeFit and Full-Modelling are also consistent with the true ΛCDM simulation cosmology down to a scale of k max = 0.20 hMpc -1 even after including the hexadecapole. For extended models such as the wCDM and the oCDM models, we find that including the hexadecapole can significantly improve the constraints and reduce the modelling errors with the same k max . While their discrepancies between the constraints from ShapeFit and Full-Modelling are more significant than ΛCDM, they remain consistent within 0.7σ. Lastly, we also show that the constraints on cosmological parameters with the correlation function evaluated from PyBird down to s min = 30h -1 Mpc are unbiased and consistent with the constraints from the power spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing the galaxy–halo connection with total satellite luminosity

ABSTRACT We demonstrate how the total luminosity in satellite galaxies is a powerful probe of dark matter haloes around central galaxies. The method cross-correlates central galaxies in spectroscopic galaxy samples with fainter galaxies detected in photometric surveys. Using models, we show that the total galaxy luminosity, Lsat, scales linearly with host halo mass, making Lsat an excellent proxy for Mh. Lsat is also sensitive to the formation time of the halo. We demonstrate that probes of galaxy large-scale environment can break this degeneracy. Although this is an indirect probe of the halo, it yields a high signal-to-noise ratio measurement for galaxies expected to occupy haloes at <1012 M⊙, where other methods suffer from larger errors. In this paper, we focus on observational and theoretical systematics in the Lsat method. We test the robustness of our method of finding central galaxies and our methods of estimating the number of background galaxies. We implement this method on galaxies in the Sloan Digital Sky Survey (SDSS) data, with satellites identified in fainter imaging data. We find excellent agreement between our theoretical predictions and the observational measurements. Finally, we compare our Lsat measurements to weak lensing estimates of Mh for red and blue subsamples. In the stellar mass range where the measurements overlap, we find consistent results, where red galaxies live in larger haloes. However, the Lsat approach allows us to probe significantly lower mass galaxies. At these masses, the Lsat values are equivalent. This example shows the potential of Lsat as a probe of dark haloes.

79 ASTRONOMY AND ASTROPHYSICS↗

Beyond the hubble sequence – exploring galaxy morphology with unsupervised machine learning

We explore unsupervised machine learning for galaxy morphology analyses using a combination of feature extraction with a vector-quantized variational autoencoder (VQ-VAE) and hierarchical clustering (HC). We propose a new methodology that includes: (1) consideration of the clustering performance simultaneously when learning features from images; (2) allowing for various distance thresholds within the HC algorithm; (3) using the galaxy orientation to determine the number of clusters. This set-up provides 27 clusters created with this unsupervised learning that we show are well separated based on galaxy shape and structure (e.g. Sérsic index, concentration, asymmetry, Gini coefficient). These resulting clusters also correlate well with physical properties such as the colour–magnitude diagram, and span the range of scaling relations such as mass versus size amongst the different machine-defined clusters. When we merge these multiple clusters into two large preliminary clusters to provide a binary classification, an accuracy of $\sim 87{{\ \rm per\ cent}}$ is reached using an imbalanced data set, matching real galaxy distributions, which includes 22.7 per cent early-type galaxies and 77.3 per cent late-type galaxies. Comparing the given clusters with classic Hubble types (ellipticals, lenticulars, early spirals, late spirals, and irregulars), we show that there is an intrinsic vagueness in visual classification systems, in particular galaxies with transitional features such as lenticulars and early spirals. Based on this, the main result in this work is not how well our unsupervised method matches visual classifications and physical properties, but that the method provides an independent classification that may be more physically meaningful than any visually based ones.

79 ASTRONOMY AND ASTROPHYSICS↗