Skip to main content

UMAP & HDBSCAN Clustering

Dataset

PHY-3500 · G1 · 43,019 spectra — UMAP / HDBSCAN results: run phy3500_umap_tsne_run_20260405T214233Z (timings, seeds) and its summary phy3500_umap_tsne_run_20260430T215407Z (HDBSCAN). Nothing on this page comes from the G2 runs.

UMAP projection​

UMAP is applied to the first 91 PCA components (95 % variance threshold) with: n_neighbors=15, min_dist=0.1, 200 epochs (t = 40.1 s on Ryzen 9 5950X).

The projection reveals a "continuous continent" structure where the Harvard spectral sequence (M → K → G → F → A) unfolds without supervision along the main axis.

Correlations:

  • ρ(UMAP axis 1, T_eff) = +0.464 — thermal gradient along first axis
  • The luminosity class (dwarfs vs giants) induces a visible topological bifurcation when coloured by log g

UMAP coloured by spectral type

UMAP projection coloured by LAMOST spectral type — Harvard sequence emerges without supervision

Negative control​

UMAP applied to data with randomly permuted columns → compact homogeneous cloud with no structure. The observed structure is of physical origin, not an algorithmic artefact.

Negative control

Permuted features → no structure. Physical origin confirmed.


HDBSCAN clustering — 20 groups without any label​

HDBSCAN (Hierarchical Density-Based Spatial Clustering of Applications with Noise) is applied to the UMAP embedding with min_cluster_size=75, min_samples=20.

Summary statistics​

Value
Total clusters20
Noise points2,643 (6.14 %)
Largest cluster (C11)16,716 stars — G-K main sequence
Second cluster (C13)10,373 stars — cool K dwarfs

Top-12 cluster populations​

ClusterN starsPhysical identification
C1116,716G-K main sequence
C1310,373Cool K dwarfs
C55,083—
C2915Subgiants / base RGB
C19900Subgiants / base RGB
C12851Subgiants / base RGB
C1654Subgiants / base RGB
C16620—
C17615—
C3818—
C9576—
C7474—

Physical identification on the HR diagram​

ClustersT_efflog gPhysical interpretation
C19, C1, C125000–5500 K≈ 3.0–3.5Subgiants / base RGB — strengthened Ca II, weakened Balmer → lower surface radiation pressure
C134500–5000 K≈ 4.5Cool K dwarfs — clearest physical identity
C115000–6000 K≈ 4.0–4.5G-K main sequence — dominant population
Original astrophysical result

Without any label or log g information, UMAP organises spectra such that HDBSCAN recovers the dwarf/subgiant distinction at fixed T_eff. Clusters C19, C1, and C12 are spectrally distinguished by strengthened Ca II lines and weakened Balmer lines — the signature of lower surface radiation pressure in giant atmospheres.

HDBSCAN clustersClusters on HR diagram
20 HDBSCAN clusters · 6.14 % noise · C11 = 16,716 starsClusters on HR diagram: C19/C1/C12 = subgiants · C13 = K dwarfs
Feature heatmap per clusterFeature profiles
SHAP-weighted feature heatmap by clusterSpectral feature profiles for each cluster

UMAP sensitivity analysis​

Hyperparameter sensitivity was assessed over: n_neighbors ∈ 50, min_dist ∈ 0.5.

Robust results:

  • The Harvard-sequence gradient along axis 1 is stable across all tested parameters
  • The dwarf/giant bifurcation appears consistently for n_neighbors ≥ 10
  • min_dist = 0.1 offers the best balance between local and global structure

UMAP sensitivity n_neighbors

UMAP sensitivity to n_neighbors parameter


See also​