Orbital decomposition of the nuclear regions in the early-type galaxy FCC 47: Unveiling the nuclear cluster origin

Nuclear star clusters (NSCs) are among the densest stellar systems in the Universe and often coexist with supermassive black holes (SMBHs) at galaxy centres. While SMBH formation histories are essentially lost, NSCs preserve evolutionary imprints through their stellar populations and stellar kinematics, reflecting the cumulative effects of mergers, accretion, and internal dynamical evolution. We aim to investigate the orbital structure of the unusually large NSC in FCC 47 (NGC 1336) by decomposing its stellar orbits into dynamically distinct components. We extract stellar kinematics, and in particular the line-of-sight velocity distributions (LOSVDs), from VLT/MUSE integral-field spectroscopy using the non-parametric Bayes-LOSVD approach, and apply triaxial Schwarzschild orbit-superposition modelling with the DYNAMITE software. We decompose the orbit library into hot, warm, cold, and counter-rotating components. We detect triple-peaked LOSVDs in the nucleus, indicating a complex orbital structure. The NSC forms a counter-rotating, kinematically decoupled component. A hot pressure-supported component, a warm counter-rotating structure and a counter-rotating cold disk in the centre suggest hierarchical assembly via early star cluster accretion and later in situ star formation. Our orbital decomposition of FCC 47 supports a hybrid formation scenario for this NSC. Dynamically distinct substructures reflect the interplay of accretion and in situ star formation during galaxy evolution.

💡 Research Summary

This paper presents a comprehensive dynamical study of the unusually massive nuclear star cluster (NSC) at the centre of the early‑type galaxy FCC 47 (NGC 1336) in the Fornax Cluster. The authors combine high‑resolution Hubble Space Telescope imaging with adaptive‑optics assisted VLT/MUSE integral‑field spectroscopy to resolve the NSC, whose effective radius (≈ 0.75″ or 66 pc) is comparable to the MUSE point‑spread function.

The first methodological innovation is the extraction of line‑of‑sight velocity distributions (LOSVDs) using the non‑parametric Bayes‑LOSVD code. By fitting the 4750–5450 Å spectral window with a PCA‑reduced set of MILES stellar templates, the authors obtain 828 Voronoi‑binned LOSVDs with a target signal‑to‑noise of 100 per bin. Crucially, the Bayes‑LOSVD approach preserves multi‑modal structures that would be smoothed out by traditional Gauss‑Hermite expansions. In the central few bins (within the NSC’s effective radius) the LOSVDs display a striking triple‑peaked shape, indicating the superposition of several kinematically distinct stellar components.

With these LOSVDs in hand, the team constructs triaxial Schwarzschild orbit‑superposition models using the DYNAMITE software. The underlying mass model is a Multi‑Gaussian Expansion (MGE) derived from HST imaging, scaled by a spatially varying stellar mass‑to‑light (M∗/L) map that incorporates age, metallicity, and a variable initial mass function (IMF) as measured from the same MUSE data. The orbit library contains on the order of 10⁴ orbits, each sampled across the full three‑dimensional phase space. Orbits are classified according to their angular momentum into four dynamical families: (i) hot, pressure‑supported orbits with low angular momentum; (ii) warm orbits with intermediate angular momentum; (iii) cold, high‑angular‑momentum orbits forming a rotating disc; and (iv) counter‑rotating orbits that rotate opposite to the main galaxy body.

The modelling reproduces both the observed LOSVDs and the previously measured supermassive black‑hole mass (M_BH ≈ 4.4 × 10⁷ M_⊙) and NSC mass (M_NSC ≈ 7.3 × 10⁸ M_⊙). The orbital decomposition reveals a complex, multi‑component structure in the nucleus:

1. A dominant hot, pressure‑supported component that is centrally concentrated and likely represents an old, spheroidal stellar population.
2. A warm, counter‑rotating component that retains the angular momentum of early accreted globular clusters or dwarf satellites.
3. A cold, counter‑rotating disc confined to the innermost region, indicative of in‑situ star formation from gas that settled with opposite spin.
4. An additional cold, co‑rotating disc associated with the main galaxy body, but with much lower mass contribution within the NSC radius.

The presence of both counter‑rotating warm and cold components, together with the hot spheroid, points to a hybrid formation scenario. The authors argue that the NSC’s mass and overall rotation could be explained by early, dissipationless accretion of massive globular clusters that delivered a net retrograde angular momentum. Subsequent gas inflow—perhaps triggered by a minor merger or tidal interaction—then formed stars in a counter‑rotating disc, creating the observed cold component. This two‑stage process naturally accounts for the triple‑peaked LOSVDs, the kinematically decoupled nature of the NSC (its rotation axis is nearly orthogonal to that of the host galaxy), and the coexistence of pressure‑supported and rotationally supported substructures.

The study is notable for being one of the first applications of non‑parametric LOSVD extraction combined with triaxial Schwarzschild modelling to an external galaxy’s NSC. It demonstrates that detailed orbital decomposition is feasible even when the NSC is only marginally resolved, provided that high‑quality integral‑field data and sophisticated kinematic extraction are employed. The authors suggest that extending this methodology to a larger sample of NSCs will enable statistical discrimination between pure in‑situ, pure cluster‑accretion, and hybrid formation pathways, and will shed light on the co‑evolution of NSCs and central supermassive black holes across different galaxy environments.