2025-11-22T19:25:16.707578

Selecting Clusters and Protoclusters via Stellar Mass Density: I. Method and tests on Mock HSC-SSP catalogs

Vicentin, Araya-Araya, Sodré et al.
We present an algorithm designed to identify galaxy (proto)clusters in wide-area photometric surveys by first selecting their dominant galaxy-i.e., the Brightest Cluster Galaxy (BCG) or protoBCG-through the local stellar mass density traced by massive galaxies. We focus on its application to the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Wide Survey to detect candidates up to $\rm z \sim 2$. In this work, we apply the method to mock galaxy catalogs that replicate the observational constraints of the HSC-SSP Wide Survey. We derive functions that describe the probability of a massive galaxy being the dominant galaxy in a structure as a function of its stellar mass density contrast within a given redshift interval. We show that galaxies with probabilities greater than 50\% yield a sample of BCGs/protoBCGs with $\gtrsim 65\%$ purity, where most of the contamination arises from galaxies in massive groups below our cluster threshold. Using the same threshold, the resulting (proto)cluster sample achieves 80\% purity and 50\% completeness for halos with $M_{\rm{halo}} \geq 10^{14} \ M_{\odot}$, reaching nearly 100\% completeness for $M_{\rm{halo}} \geq 10^{14.5} \ M_{\odot}$. We also assign probabilistic membership to surrounding galaxies based on stellar mass and distance to the dominant galaxy, from which we define the cluster richness as the number of galaxies more likely to be true members than contaminants. This allows us to derive a halo mass-richness relation. In a companion paper, we apply the algorithm to the HSC-SSP data and compare our catalog with others based on different cluster-finding techniques and X-ray detections.
academic

Sélection des Amas et Protoclusters via la Densité de Masse Stellaire : I. Méthode et tests sur les catalogues simulés HSC-SSP

Informations Fondamentales

  • ID de l'article : 2510.10735
  • Titre : Selecting Clusters and Protoclusters via Stellar Mass Density: I. Method and tests on Mock HSC-SSP catalogs
  • Auteurs : Marcelo C. Vicentin, Pablo Araya-Araya, Laerte Sodré Jr., Michael A. Strauss
  • Classification : astro-ph.CO astro-ph.GA
  • Date de publication : 14 octobre 2025 (version préliminaire)
  • Lien de l'article : https://arxiv.org/abs/2510.10735

Résumé

Cet article propose un nouvel algorithme pour identifier les amas (proto)galactiques via la densité de masse stellaire. La méthode sélectionne d'abord les galaxies dominantes (c'est-à-dire les galaxies les plus brillantes de l'amas BCG ou proto-BCG) en traçant la densité locale de masse stellaire via les galaxies massives. L'étude se concentre sur l'application de cette méthode au Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Wide Survey pour détecter des candidats dans la plage de décalage spectral z ~ 2. En appliquant la méthode à des catalogues de galaxies simulées contraints par les observations du HSC-SSP Wide Survey, les auteurs déduisent une fonction décrivant la probabilité qu'une galaxie massive devienne une galaxie dominante de structure dans un intervalle de décalage spectral donné. Les résultats montrent que les échantillons de BCGs/protoBCGs générés par des galaxies avec une probabilité supérieure à 50% présentent une pureté ≳65%, la plupart de la contamination provenant de galaxies dans des groupes massifs en dessous du seuil d'amas. En utilisant le même seuil, l'échantillon d'(proto)amas obtenu atteint une pureté de 80% et une complétude de 50% pour les halos avec Mhalo ≥ 10^14 M⊙, et une complétude proche de 100% pour Mhalo ≥ 10^14.5 M⊙.

Contexte et Motivation de la Recherche

Questions Scientifiques

Les amas de galaxies sont les plus grandes structures liées gravitationnellement de l'univers, traçant les régions les plus denses à l'échelle cosmique. Leur identification est cruciale pour comprendre l'évolution de la structure à grande échelle et les mécanismes d'évolution galactique. Cependant, l'identification des (proto)amas de galaxies à haut décalage spectral (z > 1) présente des défis majeurs :

  1. État d'évolution structurelle : Dans la plage 1 < z < 2, une fraction considérable de structures est encore en formation, appelées protoclusters de galaxies, et n'a pas encore atteint l'équilibre
  2. Limitations observationnelles : À haut décalage spectral, les échantillons spectroscopiques sont réduits et non uniformes, posant des défis pour l'estimation du décalage spectral photométrique et l'étalonnage des couleurs des galaxies de la séquence rouge
  3. Distribution spatiale : Les galaxies dans les protoclusters peuvent être dispersées sur plusieurs mégaparsecs comobiles

Limitations des Méthodes Existantes

Les algorithmes de détection d'amas de galaxies existants incluent principalement :

  • redMaPPer et CAMIRA : se concentrent sur la détection de la « séquence rouge » dans le diagramme couleur-magnitude des amas de galaxies
  • Autres méthodes : dépendent uniquement de la distribution des galaxies en séparation angulaire et décalage spectral photométrique

Ces méthodes réussissent avec un taux supérieur à 60% pour z < 1, mais leur efficacité diminue à haut décalage spectral, car la séquence rouge n'est souvent pas clairement définie.

Points Novateurs de la Recherche

La méthode proposée dans cet article présente les avantages suivants :

  1. Indépendance des hypothèses de couleur : ne dépend pas directement des hypothèses de couleur galactique
  2. Identification prioritaire de la galaxie dominante : identifie d'abord la galaxie dominante plutôt que les caractéristiques de la structure globale
  3. Applicabilité à haut décalage spectral : particulièrement adaptée aux régions à haut décalage spectral où les protoclusters sont plus courants

Contributions Principales

  1. Proposition d'un nouvel algorithme de détection d'amas de galaxies : nouvelle méthode basée sur le contraste de densité de masse stellaire pour identifier les galaxies dominantes
  2. Établissement d'un modèle probabiliste : déduction d'une fonction décrivant la probabilité qu'une galaxie devienne une galaxie dominante de structure
  3. Vérification de l'efficacité de l'algorithme : validation sur les données simulées PCcones, réalisant une pureté BCG/protoBCG supérieure à 65%
  4. Méthode de sélection des membres : fourniture d'une sélection probabiliste des membres basée sur la masse stellaire et la distance à la galaxie dominante
  5. Relation masse-abondance : établissement d'une relation masse-abondance du halo pour l'estimation de masse

Détails de la Méthode

Définition de la Tâche

Entrée : catalogue de galaxies d'un relevé photométrique (position, décalage spectral photométrique, masse stellaire) Sortie : candidats d'(proto)amas de galaxies avec galaxies dominantes, galaxies membres, abondance Contraintes : applicabilité aux limitations observationnelles du HSC-SSP Wide Survey, plage de détection 0.1 < z < 2

Architecture de l'Algorithme

1. Présélection des Candidats Galaxies Dominantes

Pour chaque candidat galaxie dominante massive présélectionné (i), on calcule la masse stellaire totale dans un volume cylindrique centré sur ce candidat :

  • Définition du cylindre : rayon r = 1 Mpc, hauteur correspondant à la distance comobile dans la tranche de décalage spectral
  • Tranche de décalage spectral : Δzi = zi - σz(1 + zi), zi + σz(1 + zi)

2. Calcul du Contraste de Densité

Le volume cylindrique est divisé en trois anneaux concentriques équidistants, avec application d'un facteur de pondération basé sur l'inverse de la distance radiale projetée :

ρ̂i = Σj=1³ [M★,i,j^tot / (π(rj² - rj-1²)dc(Δzi))] × (rj/cMpc)^(-w) / Σj=1³ (rj/cMpc)^(-w)

où w = 0.8 est le poids optimisé.

Le contraste de densité est défini comme :

δρi = (ρ̂i - ρ̄) / ρ̄

3. Modélisation Probabiliste

Construction d'un modèle probabiliste basé sur la distribution du contraste de densité, utilisant une fonction sigmoïde modifiée :

f(δρ; a,b,c,d) = a / {1 + exp[-b(δρ - c)]} + d

Algorithme de Sélection des Membres

  1. Sélection des caractéristiques : basée sur la masse stellaire photométrique et la distance à la BCG
  2. Calcul de probabilité : P(Member|M★,phot, ddominant, zdominant) vs P(Cont|M★,phot, ddominant, zdominant)
  3. Détermination de l'appartenance : les galaxies avec P(Member) > P(Cont) sont considérées comme membres
  4. Définition de l'abondance : nombre de galaxies satisfaisant les conditions

Configuration Expérimentale

Ensemble de Données

Cônes de lumière simulés PCcones :

  • Basés sur le modèle semi-analytique L-GALAXIES de la simulation Millennium
  • 10 cônes de lumière, chacun de 36 deg²
  • Simulation des contraintes observationnelles du HSC-SSP Wide Survey
  • Incluant les filtres grizY + données infrarouge W1/W2

Données de Vérification Observationnelle :

  • Candidats amas CAMIRA HSC-SSP wide
  • Catalogue Wen & Han 2021 HSC-SSP wide
  • Candidats amas redMaPPer

Métriques d'Évaluation

  • Complétude : proportion de galaxies dominantes réelles correctement identifiées
  • Pureté : proportion de galaxies dominantes réelles dans l'échantillon identifié
  • Taux de contamination : proportion de différents types de sources de contamination

Détails de Mise en Œuvre

  • Intervalles de décalage spectral : 6 intervalles de décalage spectral photométrique [0.1,0.45), [0.45,0.7), [0.7,1.05), [1.05,1.3), [1.3,1.5), [1.5,2)
  • Seuils de masse : log(M★,phot/M⊙) = 11, 11, 11, 10.5, 10.5, 10.5 correspondant à chaque intervalle de décalage spectral
  • Définition de structure : Mhalo ≥ 10^14 M⊙ comme seuil d'amas de galaxies

Résultats Expérimentaux

Résultats Principaux

Efficacité d'Identification BCG/protoBCG

  • Seuil Pdominant > 50% : pureté BCGs/protoBCGs ≥ 65%
  • Sources de contamination principales : galaxies dans des groupes massifs (en dessous du seuil d'amas)
  • Évolution avec le décalage spectral : proportion de sélection BCG diminue de 52% (bas z) à 33% (z>1.5)

Performance de Détection d'Amas de Galaxies

  • Pdominant ≥ 0.5 :
    • Mhalo ≥ 10^14 M⊙ : 80% pureté, 50% complétude
    • Mhalo ≥ 10^14.5 M⊙ : complétude proche de 100%
  • Pdominant ≥ 0.8 : pureté augmente à ~95%

Relation Masse-Abondance

Établissement de relations masse-abondance du halo pour différents intervalles de décalage spectral :

  • Utilisation de relation log-linéaire : log(Mhalo) = α × λ + β
  • Plage de pentes α : 0.022-0.053
  • Plage d'ordonnées à l'origine β : 13.140-13.769

Vérification de Cohérence

La comparaison avec les catalogues observationnels montre :

  • Propriétés BCG : magnitude en bande i, couleur r-i, masse stellaire cohérentes avec les observations
  • Distribution spatiale : dispersion de vitesse et profil radial cohérents avec redMaPPer (test KS p=0.245)

Travaux Connexes

Méthodes Traditionnelles

  1. Méthodes de séquence rouge : redMaPPer, CAMIRA basées sur la détection de séquence rouge de galaxies passives
  2. Méthodes de densité : méthodes FOF basées sur la distribution spatiale et le décalage spectral des galaxies
  3. Méthodes multi-longueurs d'onde : détection par rayons X, effet Sunyaev-Zeldovich

Innovations de cet Article

Par rapport aux méthodes traditionnelles, l'algorithme proposé :

  • Ne dépend pas de l'hypothèse de séquence rouge, applicable à haut décalage spectral
  • Identifie d'abord la galaxie dominante plutôt que la structure globale
  • Combine les informations de densité de masse stellaire pour améliorer la précision d'identification

Conclusions et Discussion

Conclusions Principales

  1. Validité de l'algorithme : le nouvel algorithme peut identifier efficacement les (proto)amas de galaxies dans la plage z < 2
  2. Indicateurs de performance : réalise une bonne performance dans un équilibre raisonnable pureté-complétude
  3. Applicabilité à haut décalage spectral : particulièrement adaptée aux régions à haut décalage spectral où les protoclusters dominent
  4. Sélection des membres : la méthode de sélection probabiliste des membres peut déterminer efficacement l'abondance de l'amas

Limitations

  1. Dépendance de la précision photométrique : la précision de l'estimation du décalage spectral photométrique et de la masse stellaire diminue à haut décalage spectral
  2. Couverture W1/W2 : la couverture réduite des bandes infrarouge à haut décalage spectral affecte les performances
  3. Effets de projection : les effets de projection à haut décalage spectral peuvent augmenter la contamination
  4. Seuil de masse : le choix du seuil 10^14 M⊙ peut affecter la définition des protoclusters

Directions Futures

  1. Application aux données réelles : appliquer l'algorithme aux données observationnelles réelles du HSC-SSP
  2. Combinaison multi-longueurs d'onde : combiner les informations multi-longueurs d'onde telles que les rayons X pour la vérification
  3. Généralisation de la méthode : adapter à d'autres projets de relevé (DES, LSST, etc.)
  4. Optimisation pour haut décalage spectral : optimiser davantage l'algorithme pour la région z > 1.5

Évaluation Approfondie

Points Forts

  1. Innovation méthodologique : propose une stratégie d'identification unique « priorité à la galaxie dominante »
  2. Expérimentation complète : vérification complète sur des données simulées détaillées
  3. Bonne performance : atteint un bon équilibre entre pureté et complétude
  4. Forte applicabilité : particulièrement adaptée à la détection de protoclusters à haut décalage spectral
  5. Rédaction claire : description détaillée de la méthode, conception expérimentale raisonnable

Insuffisances

  1. Dépendance aux simulations : les résultats sont principalement basés sur des données simulées, l'efficacité de l'application réelle reste à vérifier
  2. Optimisation des paramètres : le choix de plusieurs paramètres empiriques (comme w=0.8) manque de justification théorique
  3. Complexité computationnelle : l'analyse du coût computationnel du calcul de densité est insuffisante
  4. Erreurs systématiques : l'analyse de robustesse aux erreurs systématiques photométriques est limitée

Impact

  1. Valeur académique : contribution méthodologique nouvelle à la détection d'amas de galaxies
  2. Valeur pratique : susceptible de jouer un rôle important dans les prochains grands relevés
  3. Reproductibilité : description détaillée de la méthode, facilitant la reproduction et l'amélioration
  4. Avancement du domaine : fournit un outil puissant pour la recherche d'amas de galaxies à haut décalage spectral

Scénarios d'Application

  1. Grands relevés photométriques : HSC-SSP, DES, LSST, etc.
  2. Détection à haut décalage spectral : particulièrement adaptée à la recherche de protoclusters pour z > 1
  3. Observations spectroscopiques de suivi : sélection de cibles pour les instruments spectroscopiques tels que PFS
  4. Recherche cosmologique : évolution de la structure à grande échelle et études d'évolution galactique

Références

L'article cite de nombreuses recherches connexes, incluant principalement :

  • Algorithmes de détection d'amas de galaxies : Rykoff et al. (2014), Oguri (2014)
  • Relevé HSC-SSP : Aihara et al. (2018, 2022)
  • Travaux de simulation : Henriques et al. (2015), Springel (2005)
  • Recherche sur les protoclusters : Overzier (2016), Toshikawa et al. (2018)

Résumé : Ceci est un article de méthodologie en astrophysique de haute qualité proposant un algorithme innovant de détection d'amas de galaxies avec vérification complète. Cette méthode possède des avantages uniques dans la détection de protoclusters à haut décalage spectral, fournissant un outil précieux pour les prochains grands projets de relevé. Les contributions principales de l'article résident dans l'innovation méthodologique et la vérification systématique, apportant une contribution importante au développement de ce domaine.