Bursty star formation in dwarf galaxies can slowly transform a steep dark matter cusp into a constant density core. We explore the possibility that globular clusters (GCs) retain a dynamical memory of this transformation. To test this, we use the nbody6df code to simulate the dynamical evolution of GCs, including stellar evolution, orbiting in static and time-varying potentials for a Hubble time. We find that GCs orbiting within a cored dark matter halo, or within a halo that has undergone a cusp-core transformation, grow to a size that is substantially larger ($R_{\rm eff} > 10$ pc) than those in a static cusped dark matter halo. They also produce much less tidal debris. We find that the cleanest signal of an historic cusp-core transformation is the presence of large GCs with tidal debris. However, the effect is small and will be challenging to observe in real galaxies. Finally, we qualitatively compare our simulated GCs with the observed GC populations in the Fornax, NGC 6822, IKN and Sagittarius dwarf galaxies. We find that the GCs in these dwarf galaxies are systematically larger ($\langle R_{\rm eff}\rangle \simeq 7.8$ pc), and have substantially more scatter in their sizes, than in-situ metal rich GCs in the Milky Way and young massive star clusters forming in M83 ($\langle R_{\rm eff} \rangle \simeq 2.5$ pc). We show that the size, scatter and survival of GCs in dwarf galaxies are all consistent with them having evolved in a constant density core, or a potential that has undergone a cusp-core transformation, but not in a dark matter cusp.

We present a model for the concurrent formation of globular clusters (GCs) and supermassive stars (SMSs, $>10^3\,{\rm M}_\odot$) to address the origin of the HeCNONaMgAl abundance anomalies in GCs. GCs form in converging gas flows and accumulate low-angular momentum gas, which accretes onto protostars. This leads to an adiabatic contraction of the cluster and an increase of the stellar collision rate. A SMS can form via runaway collisions if the cluster reaches sufficiently high density before two-body relaxation halts the contraction. This condition is met if the number of stars $\gtrsim10^6$ and the gas accretion rate $\gtrsim10^5\,{\rm M}_\odot$/Myr, reminiscent of GC formation in high gas-density environments, such as -- but not restricted to -- the early Universe. The strong SMS wind mixes with the inflowing pristine gas, such that the protostars accrete diluted hot-hydrogen burning yields of the SMS. Because of continuous rejuvenation, the amount of processed material liberated by the SMS can be an order of magnitude higher than its maximum mass. This `conveyor-belt' production of hot-hydrogen burning products provides a solution to the mass budget problem that plagues other scenarios. Additionally, the liberated material is mildly enriched in helium and relatively rich in other hot-hydrogen burning products, in agreement with abundances of GCs today. Finally, we find a super-linear scaling between the amount of processed material and cluster mass, providing an explanation for the observed increase of the fraction of processed material with GC mass. We discuss open questions of this new GC enrichment scenario and propose observational tests.

We present a new technique to probe the central dark matter (DM) density profile of galaxies that harnesses both the survival and observed properties of star clusters. As a first application, we apply our method to the `ultra-faint' dwarf Eridanus II (Eri II) that has a lone star cluster ~45 pc from its centre. Using a grid of collisional $N$-body simulations, incorporating the effects of stellar evolution, external tides and dynamical friction, we show that a DM core for Eri II naturally reproduces the size and the projected position of its star cluster. By contrast, a dense cusped galaxy requires the cluster to lie implausibly far from the centre of Eri II (>1 kpc), with a high inclination orbit that must be observed at a particular orbital phase. Our results, therefore, favour a dark matter core. This implies that either a cold DM cusp was `heated up' at the centre of Eri II by bursty star formation, or we are seeing an evidence for physics beyond cold DM.

The central parsec of the Milky Way hosts two puzzlingly young stellar populations, a tight isotropic distribution of B stars around SgrA* (the S-stars) and a disk of OB stars extending to ~0.5pc. Using a modified version of Sverre Aarseth's direct summation code NBODY6 we explore the scenario in which a young star cluster migrates to the Galactic Centre within the lifetime of the OB disk population via dynamical friction. We find that star clusters massive and dense enough to reach the central parsec form a very massive star via physical collisions on a mass segregation timescale. We follow the evolution of the merger product using the most up to date, yet conservative, mass loss recipes for very massive stars. Over a large range of initial conditions, we find that the very massive star expels most of its mass via a strong stellar wind, eventually collapsing to form a black hole of mass 20 - 400 M_Sun, incapable of bringing massive stars to the Galactic Centre. No massive intermediate mass black hole can form in this scenario. The presence of a star cluster in the central ~10 pc within the last 15 Myr would also leave a ~2 pc ring of massive stars, which is not currently observed. Thus, we conclude that the star cluster migration model is highly unlikely to be the origin of either young population, and in-situ formation models or binary disruptions are favoured.

We present a dynamical friction model based on Chandrasekhar's formula that reproduces the fast inspiral and stalling experienced by satellites orbiting galaxies with a large constant density core. We show that the fast inspiral phase does not owe to resonance. Rather, it owes to the background velocity distribution function for the constant density cores being dissimilar from the usually-assumed Maxwellian distribution. Using the correct background velocity distribution function and the semi-analytic model from Petts et al. (2015), we are able to correctly reproduce the infall rate in both cored and cusped potentials. However, in the case of large cores, our model is no longer able to correctly capture core-stalling. We show that this stalling owes to the tidal radius of the satellite approaching the size of the core. By switching off dynamical friction when rt(r) = r (where rt is the tidal radius at the satellite's position) we arrive at a model which reproduces the N-body results remarkably well. Since the tidal radius can be very large for constant density background distributions, our model recovers the result that stalling can occur for Ms/Menc << 1, where Ms and Menc are the mass of the satellite and the enclosed galaxy mass, respectively. Finally, we include the contribution to dynamical friction that comes from stars moving faster than the satellite. This next-to-leading order effect becomes the dominant driver of inspiral near the core region, prior to stalling.

We present a new semi-analytic model for dynamical friction based on Chandrasekhar's formalism. The key novelty is the introduction of physically motivated, radially varying, maximum and minimum impact parameters. With these, our model gives an excellent match to full N-body simulations for isotropic background density distributions, both cuspy and shallow, without any fine-tuning of the model parameters. In particular, we are able to reproduce the dramatic core-stalling effect that occurs in shallow/constant density cores, for the first time. This gives us new physical insight into the core-stalling phenomenon. We show that core stalling occurs in the limit in which the product of the Coulomb logarithm and the local fraction of stars with velocity lower than the infalling body tends to zero. For cuspy backgrounds, this occurs when the infalling mass approaches the enclosed background mass. For cored backgrounds, it occurs at larger distances from the centre, due to a combination of a rapidly increasing minimum impact parameter and a lack of slow moving stars in the core. This demonstrates that the physics of core-stalling is likely the same for both massive infalling objects and low-mass objects moving in shallow density backgrounds. We implement our prescription for dynamical friction in the direct summation code NBODY6 as an analytic correction for stars that remain within the Roche volume of the infalling object. This approach is computationally efficient, since only stars in the inspiralling system need to be evolved with direct summation. Our method can be applied to study a variety of astrophysical systems, including young star clusters orbiting near the Galactic Centre; globular clusters moving within the Galaxy; and dwarf galaxies orbiting within dark matter halos.