We present a detailed chemical abundance analysis of the three brightest known stars from each of the Turranburra and Willka Yaku stellar streams using high-resolution Magellan/Magellan Inamori Kyocera Echelle spectra. Abundances for 27 elements, ranging from carbon to dysprosium, were derived. Our results support the original classification that Turranburra, with a low average metallicity of [Fe/H] = −2.45 ± 0.07, likely originates from a dwarf galaxy progenitor. Willka Yaku has a low average metallicity of [Fe/H] = −2.35 ± 0.03 with a small scatter in the abundances, consistent with a globular cluster progenitor as suggested by previous studies. Both streams exhibit mild enhancements in neutron-capture elements, with averages of [Eu II /Fe] = 0.47 ± 0.09 for Turranburra and 0.44 ± 0.05 for Willka Yaku, consistent with enrichment from an r -process event. A similar enrichment is observed in other stellar streams, and we further discuss this signature as it relates to the potential enrichment histories of these two streams.

The Fornax dwarf spheroidal (dSph) galaxy is one of only two Local Group dSphs that host a population of globular star clusters (GCs), the present-day properties of which have often been used to investigate the nature of the dark matter. An additional overdensity of stars called Fornax 6, lying at a projected distance of ∼0.3 kpc from the galaxy centre, was recently identified as a sixth GC residing in Fornax. However, this conclusion was based on shallow, low-resolution photometric observations and a limited number of spectra. Here, we reinvestigate the nature of Fornax 6 by analysing MUSE/Wide-Field-Mode observations alongside deep GMOS-S imagery. Using a sample of 132 spectra of red giant and horizontal branch stars, we confirm the Fornax 6 overdensity as a distinct chemo-dynamical component with respect to the surrounding field population. Specifically, we identify 43 likely members associated with the cluster, from which we measure [Fe/H] =-0.61±0.03dex, v_ los = 50.9^+0.8_-0.7 km s^-1, σ_v = 3.2^+1.4_-1.5 km s^-1. The main-sequence turn-off, observed here for the first time, strongly suggests that these stars are coeval and well described by a ∼3 Gyr old isochrone. The cluster is characterised by an irregular morphology, a large half-light radius (R_h=8.8^+1.1_-1.3pc), a small flattening (e=0.14^+0.12_-0.10), and low luminosity (M_V = -5.0 ± 0.4). Our improved characterisation of Fornax 6 supports its classification as a genuine low-mass cluster likely undergoing tidal disruption, making it the youngest and most metal-rich member of Fornax's unique GC system.

We present deep HST photometry for 26 globular clusters (GCs) residing in three Local Group dwarf galaxies: the isolated dwarf irregular (dIrr) galaxy, NGC6822, and the M31 dwarf elliptical (dE) satellites, NGC147 and NGC185. From their colour-magnitude diagrams (CMDs), we quantify their red giant branch (RGB) and horizontal branch (HB) morphologies, and employ new empirical relationships to derive measurements of metallicity, reddening and distance. Additionally, we measure sizes and V-band magnitudes from their integrated light profiles. We find that the clusters span a range of metallicities, from [Fe/H] ∼-0.7 to ≲-2; however, the three dwarfs have very similar mean GC metallicities (∼-1.7 dex) despite their very different evolutionary histories. In contrast, we find that almost all of NGC6822's GCs exhibit red HB morphologies, whereas those in the two dEs are predominantly blue. We highlight three outlying GCs in NGC6822 that have very low metallicities yet very red HBs; they are also the most extended clusters in our sample, with half-light radii of ∼ 13-17pc. If these clusters are young, their origin is difficult to explain given their very remote locations, and metallicities which are significantly lower than those of the old and intermediate-age stars in NGC6822. These three clusters are strikingly similar to several GCs linked to substructure in the outer halo of M31, suggesting either the recent accretion of an NGC6822-like dwarf by M31, or that both M31 and NGC6822 have recently accreted a similar low-mass system. Additionally, we find compelling evidence that another NGC6822 cluster, SC7, is tidally distorted.

We present Hubble Space Telescope Advanced Camera for Surveys observations of Gep I, a globular cluster (GC) candidate in the low-mass (M_⋆∼ 6.5 × 10^5 M_⊙) M31 dwarf spheroidal (dSph) satellite Andromeda XXV (And XXV). We confirm the nature of this object and provide the first detailed characterisation of its resolved stellar populations using a colour-magnitude diagram (CMD) that reaches 2 magnitudes below the horizontal branch. We compare Gep I's metallicity and distance with those of the surrounding And XXV stellar population, and find them to be strikingly similar, consistent with a physical association between the GC and the dSph. Gep I is very extended (R_h=24^+5_-4 pc) and faint (M_V = -4.5 ± 0.2 mag), similar to the star clusters residing in other low-mass dwarf galaxies. It is characterised by a very low metallicity ([Fe/H] = -2.4^+0.3_-0.4dex) and a red horizontal branch morphology, a combination also seen in suspected accreted GCs in the M31 halo and in the Local Group dwarf irregular galaxy NGC 6822. While Gep I is most likely a genuine star cluster, the current data do not exclude the tantalising possibility that it consists of And XXV stars temporarily captured by a dark subhalo orbiting within the dSph's potential well.

As ancient stellar systems, globular clusters (GCs) offer valuable insights into the dynamical histories of large galaxies. Previous studies of GC populations in the inner and outer regions of the Andromeda Galaxy (M31) have revealed intriguing subpopulations with distinct kinematic properties. Here, we build upon earlier studies by employing Bayesian modelling to investigate the kinematics of the combined inner and outer GC populations of M31. Given the heterogeneous nature of the data, we examine subpopulations defined by GCs' metallicity and by associations with substructure, in order to characterise possible relationships between the inner and outer GC populations. We find that lower-metallicity GCs and those linked to substructures exhibit a common, more rapid rotation, whose alignment is distinct from that of higher-metallicity and non-substructure GCs. Furthermore, the higher-metallicity GCs rotate in alignment with Andromeda's stellar disk. These pronounced kinematic differences reinforce the idea that different subgroups of GCs were accreted to M31 at distinct epochs, shedding light on the complex assembly history of the galaxy.

From >1000 orbits of HST imaging, we present deep homogeneous resolved star color–magnitude diagrams that reach the oldest main-sequence turnoff and uniformly measured star formation histories (SFHs) of 36 dwarf galaxies (−6 ≥ M V ≥ −17) associated with the M31 halo, and for 10 additional fields in M31, M33, and the Giant Stellar Stream. From our SFHs, we find: (i) The median stellar age and quenching epoch of M31 satellites correlate with galaxy luminosity and galactocentric distance. Satellite luminosity and present-day distance from M31 predict the satellite quenching epoch to within 1.8 Gyr at all epochs. This tight relationship highlights the fundamental connection between satellite halo mass, environmental history, and star formation duration. (ii) There is no difference between the median SFH of galaxies on and off the great plane of Andromeda satellites. (iii) ~50% of our M31 satellites show prominent ancient star formation (>12 Gyr ago) followed by delayed quenching (8–10 Gyr ago), which is not commonly observed among the MW satellites. (iv) A comparison with TNG50 and FIRE-2 simulated satellite dwarfs around M31-like hosts shows that some of these trends (dependence of SFH on satellite luminosity) are reproduced in the simulations while others (dependence of SFH on galactocentric distance, presence of the delayed-quenching population) are weaker or absent. We provide all photometric catalogs and SFHs as High-Level Science Products on MAST.

M31 hosts a rich population of outer halo (R-proj > 25 kpc) globular clusters (GCs), many of which show strong evidence for spatial and/or kinematical associations with large-scale tidal debris features. We present deep Hubble Space Telescope photometry of 48 halo GCs, including 18 with clear ties to stellar streams and 13 with potential associations. Using the colour-magnitude diagrams (CMDs), we quantify the horizontal branch (HB) morphologies and employ new empirical relationships, calibrated on Milky Way (MW) GCs, to consistently derive metallicities and line-of-sight extinctions. We find a remarkable correlation between HB morphology and the presence of substructure: GCs with very red HBs are almost exclusively associated with substructure, while 'non-substructure' GCs have extended blue HBs. This provides the first direct evidence that red HB halo clusters originate from satellite accretion, a notion introduced nearly 50 yr ago from MW studies which has remained unconfirmed until now. In addition to a more metal-rich tail, the substructure GC sample also contains a few clusters with very low metallicities and red HBs, unlike any objects known in the MW. We suggest these are recently accreted young clusters, supporting the growing evidence that M31 has experienced a more prolonged accretion history than the MW.

ABSTRACT We present detailed abundance results based on Ultraviolet and Visual Echelle Spectrograph high-dispersion spectra for seven very and extremely metal-poor stars in the Large Magellanic Cloud (LMC). We confirm that all seven stars, two of which have [Fe/H] ≤ −3.0, are the most metal-poor stars discovered so far in the Magellanic Clouds. The element abundance ratios are generally consistent with Milky Way halo stars of similar [Fe/H] values. We find that two of the more metal-rich stars in our sample are enhanced in r-process elements. This result contrasts with the literature, where all nine metal-poor LMC stars with higher [Fe/H] values than our sample were found to be rich in r-process elements. The absence of r-process enrichment in stars with lower [Fe/H] values is consistent with a minimum delay time-scale of ∼100 Myr for the neutron star binary merger process to generate substantial r-process enhancements in the LMC. We find that the occurrence rate of r-process enhancement (r-I or r-II) in our sample of very and extremely metal-poor stars is statistically indistinguishable from that found in the Milky Way’s halo, although including stars from the literature sample hints at a larger r-II frequency the LMC. Overall, our results shed light on the earliest epochs of star formation in the LMC that may be applicable to other galaxies of LMC-like mass.

We present the lifetime star formation histories (SFHs) for six ultrafaint dwarf (UFD; M V > − 7.0, 4.9 < log 10 ( M * ( z = 0 ) / M ⊙ ) < 5.5 ) satellite galaxies of M31 based on deep color–magnitude diagrams constructed from Hubble Space Telescope imaging. These are the first SFHs obtained from the oldest main-sequence turnoff of UFDs outside the halo of the Milky Way (MW). We find that five UFDs formed at least 50% of their stellar mass by z = 5 (12.6 Gyr ago), similar to known UFDs around the MW, but that 10%–40% of their stellar mass formed at later times. We uncover one remarkable UFD, And xiii , which formed only 10% of its stellar mass by z = 5, and 75% in a rapid burst at z ∼ 2–3, a result that is robust to choices of underlying stellar model and is consistent with its predominantly red horizontal branch. This “young” UFD is the first of its kind and indicates that not all UFDs are necessarily quenched by reionization, which is consistent with predictions from several cosmological simulations of faint dwarf galaxies. SFHs of the combined MW and M31 samples suggest reionization did not homogeneously quench UFDs. We find that the least-massive MW UFDs ( M * ( z = 5) ≲ 5 × 10 4 M ⊙ ) are likely quenched by reionization, whereas more-massive M31 UFDs ( M * ( z = 5) ≳ 10 5 M ⊙ ) may only have their star formation suppressed by reionization and quench at a later time. We discuss these findings in the context of the evolution and quenching of UFDs.

ABSTRACT Young star clusters enable us to study the effects of stellar rotation on an ensemble of stars of the same age and across a wide range in stellar mass and are therefore ideal targets for understanding the consequences of rotation on stellar evolution. We combine MUSE spectroscopy with HST photometry to measure the projected rotational velocities (Vsin i) of 2184 stars along the split main sequence and on the main sequence turn-off (MSTO) of the 100 Myr-old massive ($10^5\, {\rm M_{\odot }}$) star cluster NGC 1850 in the Large Magellanic Cloud. At fixed magnitude, we observe a clear correlation between Vsin i and colour, in the sense that fast rotators appear redder. The average Vsin i values for stars on the blue and red branches of the split main sequence are $\sim \! 100\, {\rm km\, s^{-1}}$ and $\sim \! 200\, {\rm km\, s^{-1}}$, respectively. The values correspond to about $25-30{{\ \rm per\ cent}}$ and $50-60{{\ \rm per\ cent}}$ of the critical rotation velocity and imply that rotation rates comparable to those observed in field stars of similar masses can explain the split main sequence. Our spectroscopic sample contains a rich population of ∼200 fast rotating Be stars. The presence of shell features suggests that 23 per cent of them are observed through their decretion discs, corresponding to a disc opening angle of 15 degrees. These shell stars can significantly alter the shape of the MSTO, hence care should be taken when interpreting this photometric feature. Overall, our findings impact our understanding of the evolution of young massive clusters and provide new observational constraints for testing stellar evolutionary models.

ABSTRACT The advent of large-scale stellar spectroscopic surveys naturally leads to the implementation of machine learning techniques to isolate, for example, small sub-samples of potentially interesting stars from the full data set. A recent example is the application of the t-SNE statistical method to ∼600 000 stellar spectra from the GALAH survey in order to identify a sample of candidate extremely metal-poor (EMP; [Fe/H] ≤ –3) stars. We report the outcome of low-resolution spectroscopic follow-up of 83 GALAH EMP candidates that lack any previous metallicity estimates. Overall, the statistical selection is found to be efficient (∼one-third of the candidates have [Fe/H] ≤ –2.75) with low contamination (<10 per cent have [Fe/H] > –2), and with a metallicity distribution function that is consistent with previous work. Five stars are found to have [Fe/H] ≤ –3.0, one of which is a main sequence turnoff star. Two other stars are revealed as likely carbon-enhanced metal-poor (CEMP) stars of type CEMP-s, and a known carbon star is re-identified. The results indicate that the statistical selection approach employed was successful, and therefore it can be applied to forthcoming even larger stellar spectroscopic surveys with the expectation of similar positive outcomes.

We present results of a search for extremely metal-poor (EMP) stars in the Large Magellanic Cloud, which can provide crucial information about the properties of the first stars as well as on the formation conditions prevalent during the earliest stages of star formation in dwarf galaxies. Our search utilised SkyMapper photometry, together with parallax and proper motion cuts (from Gaia), colour-magnitude cuts (by selecting the red giant branch region) and finally a metallicity-sensitive cut. Low-resolution spectra of a sample of photometric candidates were taken using the ANU 2.3m telescope/WiFeS spectrograph, from which 7 stars with [Fe/H] $\leq$ -2.75 were identified, two of which have [Fe/H] $\leq$ -3. Radial velocities, derived from the CaII triplet lines, closely match the outer rotation curve of the LMC for the majority of the candidates in our sample. Therefore, our targets are robustly members of the LMC based on their 6D phase-space information (coordinates, spectrophotometric distance, proper motions and radial velocities), and they constitute the most metal-poor stars so far discovered in this galaxy.

Based upon a metallicity selection, we identify a significant sub-population of the inner halo globular clusters in the Andromeda Galaxy which we name the Dulais Structure. It is distinguished as a co-rotating group of 10-20 globular clusters which appear to be kinematically distinct from, and on average more metal-poor than, the majority of the inner halo population. Intriguingly, the orbital axis of this Dulais Structure is closely aligned with that of the younger accretion event recently identified using a sub-population of globular clusters in the outer halo of Andromeda, and this is strongly suggestive of a causal relationship between the two. If this connection is confirmed, a natural explanation for the kinematics of the globular clusters in the Dulais Structure is that they trace the accretion of a substantial progenitor ( 10^11 Msun) into the halo of Andromeda during the last few billion years, that may have occurred as part of a larger group infall.

We use data from the Magellanic Edges Survey (MagES) in combination with Gaia EDR3 to study the extreme southern outskirts of the Small Magellanic Cloud (SMC), focussing on a field at the eastern end of a long arm-like structure which wraps around the southern periphery of the Large Magellanic Cloud (LMC). Unlike the remainder of this structure, which is thought to be comprised of perturbed LMC disk material, the aggregate properties of the field indicate a clear connection with the SMC. We find evidence for two stellar populations in the field: one having properties consistent with the outskirts of the main SMC body, and the other significantly perturbed. The perturbed population is on average 0.2 dex more metal-rich, and is located 7 kpc in front of the dominant population with a total space velocity relative to the SMC centre of 230 km/s broadly in the direction of the LMC. We speculate on possible origins for this perturbed population, the most plausible of which is that it comprises debris from the inner SMC that has been recently tidally stripped by interactions with the LMC.

ABSTRACT We present the results of a spectroscopic survey of the outskirts of four globular – 1261, NGC 4590, NGC 1904, and NGC 1851 – covering targets within 1° from the cluster centres, with 2dF/AAOmega on the Anglo-Australian Telescope (AAT) and FLAMES on the very large telescope (VLT). We extracted chemo-dynamical information for individual stars, from which we estimated the velocity dispersion profile and the rotation of each cluster. The observations are compared to direct N-body simulations and appropriate limepy/spes models for each cluster to interpret the results. In NGC 1851, the detected internal rotation agrees with existing literature, and NGC 1261 shows some rotation signal beyond the truncation radius, likely coming from the escaped stars. We find that the dispersion profiles for both the observations and the simulations for NGC 1261, NGC 1851, and NGC 1904 do not decrease as the limepy/spes models predict beyond the truncation radius, where the N-body simulations show that escaped stars dominate; the dispersion profile of NGC 4590 follows the predictions of the limepy/spes models, though the data do not effectively extend beyond the truncation radius. The increasing/flat dispersion profiles in the outskirts of NGC 1261, NGC 1851, and NGC 1904, are reproduced by the simulations. Hence, the increasing/flat dispersion profiles of the clusters in question can be explained by the tidal interaction with the galaxy without introducing dark matter.

It has recently been shown that the Large Magellanic Cloud (LMC) has a substantial effect on the Milky Way's stellar halo and stellar streams. Here, we explore how deformations of the Milky Way and LMC's dark matter haloes affect stellar streams, and whether these effects are observable. In particular, we focus on the Orphan-Chenab (OC) stream which passes particularly close to the LMC and spans a large portion of the Milky Way's halo. We represent the Milky Way-LMC system using basis function expansions that capture their evolution in an N-body simulation. We present the properties of this system, such as the evolution of the densities and force fields of each galaxy. The OC stream is evolved in this time-dependent, deforming potential, and we investigate the effects of the various moments of the Milky Way and the LMC. We find that the simulated OC stream is strongly influenced by the deformations of both the Milky Way and the LMC and that this effect is much larger than current observational errors. In particular, the Milky Way dipole has the biggest impact on the stream, followed by the evolution of the LMC's monopole, and the LMC's quadrupole. Detecting these effects would confirm a key prediction of collisionless, cold dark matter, and would be a powerful test of alternative dark matter and alternative gravity models.

ABSTRACT We present detailed C, O, Na, Mg, Si, Ca, Ti, V, Fe, Zr, Ba, and Eu abundance measurements for 20 red giant branch (RGB) stars in the Large Magellanic Cloud (LMC) star cluster NGC 1846 ([Fe/H] = −0.59). This cluster is 1.95 Gyr old and lies just below the supposed lower age limit (2 Gyr) for the presence of multiple populations in massive star clusters. Our measurements are based on high- and low-resolution VLT/FLAMES spectra combined with photometric data from Hubble Space Telescope (HST). Corrections for non-local thermodynamic equilibrium effects are also included for O, Na, Mg, Si, Ca, Fe, and Ba. Our results show that there is no evidence for multiple populations in this cluster based on the lack of any intrinsic star-to-star spread in the abundances of Na and O: We place 95 per cent confidence limits on the intrinsic dispersion for these elements of ≤0.07 and ≤0.09 dex, respectively. However, we do detect a significant spread in the carbon abundances, indicating varying evolutionary mixing occurring on the RGB that increases with luminosity. Overall, the general abundance patterns for NGC 1846 are similar to those seen in previous studies of intermediate-age LMC star clusters and field stars.

ABSTRACT We explore the structural and kinematic properties of the outskirts of the Large Magellanic Cloud (LMC) using data from the Magellanic Edges Survey (MagES) and Gaia EDR3. Even at large galactocentric radii (8° < R < 11°), we find the north-eastern LMC disc is relatively unperturbed: its kinematics are consistent with a disc of inclination ∼36.5° and line-of-nodes position angle ∼145° east of north. In contrast, fields at similar radii in the southern and western disc are significantly perturbed from equilibrium, with non-zero radial and vertical velocities, and distances significantly in front of the disc plane implied by our north-eastern fields. We compare our observations to simple dynamical models of the Magellanic or Milky Way system which describe the LMC as a collection of tracer particles within a rigid potential, and the Small Magellanic Cloud (SMC) as a rigid Hernquist potential. A possible SMC crossing of the LMC disc plane ∼400 Myr ago, in combination with the LMC’s infall to the Milky Way potential, can qualitatively explain many of the perturbations in the outer disc. Additionally, we find the claw-like and arm-like structures south of the LMC have similar metallicities to the outer LMC disc ([Fe/H] ∼ −1), and are likely comprised of perturbed LMC disc material. The claw-like substructure is particularly disturbed, with out-of-plane velocities >60 km s−1 and apparent counter-rotation relative to the LMC’s disc motion. More detailed N-body models are necessary to elucidate the origin of these southern features, potentially requiring repeated interactions with the SMC prior to ∼1 Gyr ago.

We report the discovery of J0624-6948, a low-surface brightness radio ring, lying between the Galactic Plane and the Large Magellanic Cloud (LMC). It was first detected at 888 MHz with the Australian Square Kilometre Array Pathfinder (ASKAP), and with a diameter of 196 arcsec. This source has phenomenological similarities to Odd Radio Circles (ORCs). Significant differences to the known ORCs - a flatter radio spectral index, the lack of a prominent central galaxy as a possible host, and larger apparent size - suggest that J0624-6948 may be a different type of object. We argue that the most plausible explanation for J0624-6948 is an intergalactic supernova remnant due to a star that resided in the LMC outskirts that had undergone a single-degenerate type Ia supernova, and we are seeing its remnant expand into a rarefied, intergalactic environment. We also examine if a massive star or a white dwarf binary ejected from either galaxy could be the supernova progenitor. Finally, we consider several other hypotheses for the nature of the object, including the jets of an active galactic nucleus (AGN) or the remnant of a nearby stellar super-flare.

The Fermi bubbles are giant, γ-ray-emitting lobes emanating from the nucleus of the Milky Way discovered in ~1–100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi bubbles that has been interpreted as a signature of collimated outflows from the Galaxy’s supermassive black hole. Here we show via a spatial template analysis that much of the γ-ray emission associated with the brightest region of substructure—the so-called cocoon—is probably due to the Sagittarius dwarf spheroidal galaxy (dSph). This large Milky Way satellite is viewed through the Fermi bubbles from the position of the Solar System. As a tidally and ram-pressure stripped remnant, the Sagittarius dSph has no ongoing star formation, but we nevertheless demonstrate that the dwarf’s millisecond pulsar population can plausibly supply the γ-ray signal that our analysis associates with its stellar template. The measured spectrum is naturally explained by inverse Compton scattering of cosmic microwave background photons by high-energy electron–positron pairs injected by millisecond pulsars belonging to the Sagittarius dSph, combined with these objects’ magnetospheric emission. This finding plausibly suggests that millisecond pulsars produce significant γ-ray emission among old stellar populations, potentially confounding indirect dark-matter searches in regions such as the Galactic Centre, the Andromeda galaxy and other massive Milky Way dSphs. A bright patch in the Fermi bubbles, previously attributed to a jet launched by the Galaxy’s central black hole, is actually due to gamma-ray emission by millisecond pulsars in a background, satellite galaxy of the Milky Way.