Reproducible kk-means clustering in galaxy feature data from the GAMA survey

A fundamental bimodality of galaxies in the local Universe is apparent in many of the features used to describe them. Multiple sub-populations exist within this framework, each representing galaxies following distinct evolutionary pathways. Accurately identifying and characterising these sub-populations requires that a large number of galaxy features be analysed simultaneously. Future galaxy surveys such as LSST and Euclid will yield data volumes for which traditional approaches to galaxy classification will become unfeasible. To address this, we apply a robust $k$-means unsupervised clustering method to feature data derived from a sample of 7338 local-Universe galaxies selected from the Galaxy And Mass Assembly (GAMA) survey. This allows us to partition our sample into $k$ clusters without the need for training on pre-labelled data, facilitating a full census of our high dimensionality feature space and guarding against stochastic effects. We find that the local galaxy population natively splits into $2$, $3$, $5$ and a maximum of $6$ sub-populations, with each corresponding to a distinct ongoing evolutionary mechanism. Notably, the impact of the local environment appears strongly linked with the evolution of low-mass ($M_{*} < 10^{10}$ M$_{\odot}$) galaxies, with more massive systems appearing to evolve more passively from the blue cloud onto the red sequence. With a typical run time of $\sim3$ minutes per value of $k$ for our galaxy sample, we show how $k$-means unsupervised clustering is an ideal tool for future analysis of large extragalactic datasets, being scalable, adaptable, and providing crucial insight into the fundamental properties of the local galaxy population

What drives the evolution of gas kinematics in star-forming galaxies?

One important result from recent large integral field spectrograph (IFS) surveys is that the intrinsic velocity dispersion of galaxies traced by star-forming gas increases with redshift. Massive, rotation-dominated discs are already in place at z ∼ 2, but they are dynamically hotter than spiral galaxies in the local Universe. Although several plausible mechanisms for this elevated velocity dispersion (e.g. star formation feedback, elevated gas supply, or more frequent galaxy interactions) have been proposed, the fundamental driver of the velocity dispersion enhancement at high redshift remains unclear. We investigate the origin of this kinematic evolution using a suite of cosmological simulations from the FIRE (Feedback In Realistic Environments) project. Although IFS surveys generally cover a wider range of stellar masses than in these simulations, the simulated galaxies show trends between intrinsic velocity dispersion (σ intr ), SFR, and z in agreement with observations. In both observations and simulations, galaxies on the star-forming main sequence have median σ intr values that increase from z ∼ 0 to z ∼ 1–1.5, but this increasing trend is less evident at higher redshift. In the FIRE simulations, σ intr can vary significantly on time-scales of 100 Myr. These variations closely mirror the time evolution of the SFR and gas inflow rate (M gas ). By cross-correlating pairs of σ intr, M gas, and SFR, we show that increased gas inflow leads to subsequent enhanced star formation, and enhancements in σ intr tend to temporally coincide with increases in M gas and SFR

Stars made in outflows may populate the stellar halo of the Milky Way

We study stellar-halo formation using six Milky-Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that 5-40 per cent of the outer (50-300 kpc) stellar halo in each system consists of in-situ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by superbubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at z = 0. In the outer halo, outflow stars are more spatially homogeneous, metal-rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to ∼10 per cent of our kinematically identified halo stars were born in outflows; the fraction rises to as high as ∼40 per cent for the most metal-rich local halo stars ([Fe/H] &gt;-0.5). Such stars can be retrograde and create features similar to the recently discoveredMilkyWay 'Splash' in phase space.We conclude that theMilkyWay stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell &amp; Sandange, who suggested that halo stars formed in rapidly infalling gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme-velocity stars in the Milky Way

A Relationship Between Stellar Metallicity Gradients and Galaxy Age in Dwarf Galaxies

We explore the origin of stellar metallicity gradients in simulated and observed dwarf galaxies. We use FIRE-2 cosmological baryonic zoom-in simulations of 26 isolated galaxies as well as existing observational data for 10 Local Group dwarf galaxies. Our simulated galaxies have stellar masses between $10^{5.5}$ and 10^{8.6} \msun. Whilst gas-phase metallicty gradients are generally weak in our simulated galaxies, we find that stellar metallicity gradients are common, with central regions tending to be more metal-rich than the outer parts. The strength of the gradient is correlated with galaxy-wide median stellar age, such that galaxies with younger stellar populations have flatter gradients. Stellar metallicty gradients are set by two competing processes: (1) the steady "puffing" of old, metal-poor stars by feedback-driven potential fluctuations, and (2) the accretion of extended, metal-rich gas at late times, which fuels late-time metal-rich star formation. If recent star formation dominates, then extended, metal-rich star formation washes out pre-existing gradients from the "puffing" process. We use published results from ten Local Group dwarf galaxies to show that a similar relationship between age and stellar metallicity-gradient strength exists among real dwarfs. This suggests that observed stellar metallicity gradients may be driven largely by the baryon/feedback cycle rather than by external environmental effects.Comment: 14 pages, 13 figures, submitted to MNRA