Professor Justin Read

We present a new non-parametric Jeans code, GravSphere, that recovers the density ρ(r) and velocity anisotropy β(r) of spherical stellar systems, assuming only that they are in a steady-state. Using a large suite of mock data, we confirm that with only line-of-sight velocity data, GravSphere provides a good estimate of the density at the projected stellar half mass radius, ρ(R1/2), but is not able to measure ρ(r) or β(r), even with 10,000 tracer stars. We then test three popular methods for breaking this ρ − β degeneracy: using multiple populations with different R1/2; using higher order ‘Virial Shape Parameters’ (VSPs); and including proper motion data. We find that two populations provide an excellent recovery of ρ(r) in-between their respective R1/2. However, even with a total of ∼7, 000 tracers, we are not able to well-constrain β(r) for either population. By contrast, using 1000 tracers with higher order VSPs we are able to measure ρ(r) over the range 0.5 < r/R1/2 < 2 and broadly constrain β(r). Including proper motion data for all stars gives an even better performance, with ρ and β well-measured over the range 0.25 < r/R1/2 < 4. Finally, we test GravSphere on a triaxial mock galaxy that has axis ratios typical of a merger remnant, [1: 0.8: 0.6]. In this case, GravSphere can become slightly biased. However, we find that when this occurs the data are poorly fit, allowing us to detect when such departures from spherical symmetry become problematic.

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 core being dissimilar from the usually-assumed Maxwellian distribution. Using the correct background velocity distribution function and the semi-analytic model from Petts, Gualandris & Read (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.

The dwarf galaxy distribution surrounding M31 is significantly anisotropic in nature. Of the 30 dwarf galaxies in this distribution, 15 form a disc-like structure and 23 are contained within the hemisphere facing the Milky Way. Using a realistic local potential, we analyse the conditions required to produce and maintain these asymmetries. We find that some dwarf galaxies are required to have highly eccentric orbits in order to preserve the presence of the hemispherical asymmetry with an appropriately large radial dispersion. Under the assumption that the dwarf galaxies originate from a single association or accretion event, we find that the initial size and specific energy of that association must both be relatively large in order to produce the observed hemispherical asymmetry. However if the association was large in physical size, the very high-energy required would enable several dwarf galaxies to escape from the M31 and be captured by the Milky Way. Furthermore, we find that associations that result in this structure have total specific energies concentrated around E = V_esc2 - V_init2 ̃ 200^2 - 300^2 {km^2 s^{-2}}, implying that the initial velocity and initial position needed to produce the structure are strongly correlated. The overlap of initial conditions required to produce the radial dispersion, angular dispersion, and the planar structure is small and suggests that either they did not originate from a single accretion event, or that these asymmetric structures are short-lived.

We demonstrate how the least luminous galaxies in the universe, ultra-faint dwarf galaxies, are sensitive to their dynamical mass at the time of cosmic reionization. We select a low-mass (~ ´ 1.5 10 M☉ 9 ) dark matter halo from a cosmological volume, and perform zoom hydrodynamical simulations with multiple alternative histories using “genetically modified” initial conditions. Earlier-forming ultra-faints have higher stellar mass today, due to a longer period of star formation before their quenching by reionization. Our histories all converge to the same final dynamical mass, demonstrating the existence of extended scatter (1 dex) in stellar masses at fixed halo mass due to the diversity of possible histories. One of our variants builds less than 2% of its final dynamical mass before reionization, rapidly quenching in situ star formation. The bulk of its final stellar mass is later grown by dry mergers, depositing stars in the galaxy’s outskirts and hence expanding its effective radius. This mechanism constitutes a new formation scenario for highly diffuse (r1 2 ~ 820 pc, ~ - 32 mag arcsec 2 ), metal-poor ([Fe H 2.9 ] = - ), ultra-faint (V = -5.7) dwarf galaxies within the reach of next-generation low surface brightness surveys.

We present a novel positive potential-density pair expansion for modelling galaxies, based on the Miyamoto–Nagai disc. By using three sets of such discs, each one of them aligned along each symmetry axis, we are able to reconstruct a broad range of potentials that correspond to density profiles from exponential discs to 3D power-law models with varying triaxiality (henceforth simply ‘twisted’ models). We increase the efficiency of our expansion by allowing the scalelength parameter of each disc to be negative. We show that, for suitable priors on the scalelength and scaleheight parameters, these ‘MNn discs’ (Miyamoto–Nagai negative) have just one negative density minimum. This allows us to ensure global positivity by demanding that the total density at the global minimum is positive. We find that at better than 10 per cent accuracy in our density reconstruction, we can represent a radial and vertical exponential disc over 0.1–10 scalelengths/scaleheights with four MNn discs; a Navarro, Frenk and White (NFW) profile over 0.1–10 scalelengths with four MNn discs; and a twisted triaxial NFW profile with three MNn discs per symmetry axis. Our expansion is efficient, fully analytic, and well suited to reproducing the density distribution and gravitational potential of galaxies from discs to ellipsoids.

We introduce the 'Engineering Dwarfs at Galaxy Formation's Edge' (EDGE) project to study the cosmological formation and evolution of the smallest galaxies in the Universe. In this first paper, we explore the effects of resolution and sub-grid physics on a single low-mass halo (M_halo=109{ M}_☉), simulated to redshift z = 0 at a mass and spatial resolution of ̃ 20{ M}_☉ and ∼3 pc. We consider different star formation prescriptions, supernova feedback strengths, and on-the-fly radiative transfer (RT). We show that RT changes the mode of galactic self-regulation at this halo mass, suppressing star formation by causing the interstellar and circumgalactic gas to remain predominantly warm (∼104 K) even before cosmic reionization. By contrast, without RT, star formation regulation occurs only through starbursts and their associated vigorous galactic outflows. In spite of this difference, the entire simulation suite (with the exception of models without any feedback) matches observed dwarf galaxy sizes, velocity dispersions, V-band magnitudes, and dynamical mass-to-light-ratios. This is because such structural scaling relations are predominantly set by the host dark matter halo, with the remaining model-to-model variation being smaller than the observational scatter. We find that only the stellar mass-metallicity relation differentiates the galaxy formation models. Explosive feedback ejects more metals from the dwarf, leading to a lower metallicity at a fixed stellar mass. We conclude that the stellar mass-metallicity relation of the very smallest galaxies provides a unique constraint on galaxy formation physics.

The Magellanic Bridge (MB) is a gaseous stream that links the Large (LMC) and Small (SMC) Magellanic Clouds. Current simulations suggest that the MB forms from a recent interaction between the Clouds. In this scenario, the MB should also have an associated stellar bridge formed by stars tidally stripped from the SMC by the LMC. There are several observational evidences for these stripped stars, from the presence of intermediate age populations in the MB and carbon stars, to the recent observation of an over-density of RR Lyrae stars offset from the MB. However, spectroscopic confirmation of stripped stars in the MB remains lacking. In this paper, we use medium resolution spectra to derive the radial velocities and metallicities of stars in two fields along the MB. We show from both their chemistry and kinematics that the bulk of these stars must have been tidally stripped from the SMC. This is the first spectroscopic evidence for a dwarf galaxy being tidally stripped by a larger dwarf.

Gravitational potential fluctuations driven by bursty star formation can kinematically ‘heat up’ dark matter at the centres of dwarf galaxies. A key prediction of such models is that, at a fixed dark matter halo mass, dwarfs with a higher stellar mass will have a lower central dark matter density. We use stellar kinematics and HI gas rotation curves to infer the inner dark matter densities of eight dwarf spheroidal and eight dwarf irregular galaxies with a wide range of star formation histories. For all galaxies, we estimate the dark matter density at a common radius of 150 pc, ρDM(150pc)⁠. We find that our sample of dwarfs falls into two distinct classes. Those that stopped forming stars over 6 Gyr ago favour central densities ρDM(150pc)>108 M⊙ kpc−3, consistent with cold dark matter cusps, while those with more extended star formation favour ρDM(150pc)

We employ the earlier published proper motions of the newly discovered Antlia 2 dwarf galaxy derived from Gaia data to calculate its orbital distribution in the cosmologically recent past. Using these observationally motivated orbits, we calculate the effect of the Antlia 2 dwarf galaxy on the outer H I disk of the Milky Way, using both test particle and smoothed particle hydrodynamics simulations. We find that orbits with low pericenters, ∼10 kpc, produce disturbances that match the observed outer H I disk perturbations. We have independently recalculated the proper motion of the Antlia 2 dwarf from Gaia data and found a proper motion of (μ α cosδ, μ δ ) = (−0.068, 0.032) ± (0.023, −0.031) mas yr−1, which agrees with results from Torrealba et al. within the errors, but gives lower mean pericenters, e.g., ∼15 kpc for our fiducial model of the Milky Way. We also show that the Sagittarius dwarf galaxy interaction does not match the observed perturbations in the outer gas disk. Thus, Antlia 2 may be the driver of the observed large perturbations in the outer gas disk of the Galaxy. The current location of the Antlia 2 dwarf galaxy closely matches that predicted by an earlier dynamical analysis of the dwarf galaxy that drove ripples in the outer Galaxy, and, in particular, its orbit is nearly coplanar to the Galactic disk. If the Antlia 2 dwarf galaxy is responsible for the perturbations in the outer Galactic disk, it would have a specific range of proper motions that we predict here; this can be tested soon with Gaia DR-3 and Gaia DR-4 data.

We introduce a novel abundance matching technique that produces a more accurate estimate of the pre-infall halo mass, M200, for satellite galaxies. To achieve this, we abundance match with the mean star formation rate, averaged over the time when a galaxy was forming stars, ⟨SFR⟩, instead of the stellar mass, M∗. Using data from the Sloan Digital Sky Survey, the GAMA survey and the Bolshoi simulation, we obtain a statistical ⟨SFR⟩−M200 relation in ΛCDM. We then compare the pre-infall halo mass, Mabund200, derived from this relation with the pre-infall dynamical mass, Mdyn200, for 21 nearby dSph and dIrr galaxies, finding a good agreement between the two. As a first application, we use our new ⟨SFR⟩−M200 relation to empirically measure the cumulative mass function of a volume-complete sample of bright Milky Way satellites within 280 kpc of the Galactic centre. Comparing this with a suite of cosmological 'zoom' simulations of Milky Way-mass halos that account for subhalo depletion by the Milky Way disc, we find no missing satellites problem above M200∼09M⊙ in the Milky Way. We discuss how this empirical method can be applied to a larger sample of nearby spiral galaxies.

We fit the rotation curves of isolated dwarf galaxies to directly measure the stellar mass-halo mass relation (M∗ − M200) over the mass range 5 × 105

We derive the local dark matter density by applying the integrated Jeans equation method from Silverwood et al. to SDSS-SEGUE G-dwarf data processed and presented by Büdenbender et al. We use the MULTINEST Bayesian nested sampling software to fit a model for the baryon distribution, dark matter, and tracer stars, including a model for the `tilt term' that couples the vertical and radial motions, to the data. The α-young population from Büdenbender et al. yields the most reliable result of ρ_dm= 0.46^{+0.07}_{-0.08} {GeV cm}^{-3}= 0.012^{+0.002}_{-0.002} M_{☉} pc^{-3}. Our analyses yield inconsistent results for the α-young and α-old data, pointing to problems in the tilt term and its modelling, the data itself, the assumption of a flat rotation curve, or the effects of disequilibria.

We use a new mass modelling method, GRAVSPHERE, to measure the central dark matter density profile of the Draco dwarf spheroidal galaxy. Draco's star formation shut down long ago, making it a prime candidate for hosting a `pristine' dark matter cusp, unaffected by stellar feedback during galaxy formation. We first test GRAVSPHERE on a suite of tidally stripped mock `Draco'-like dwarfs. We show that we are able to correctly infer the dark matter density profile of both cusped and cored mocks within our 95 per cent confidence intervals. While we obtain only a weak inference on the logarithmic slope of these density profiles, we are able to obtain a robust inference of the amplitude of the inner dark matter density at 150 pc, ρ _DM(150 pc). We show that, combined with constraints on the density profile at larger radii, this is sufficient to distinguish a Λ Cold Dark Matter (ΛCDM) cusp - that has ρ _DM(150 pc) ≳ 1.8 × 10^8 M_☉ kpc^{-3} - from alternative dark matter models that have lower inner densities. We then apply GRAVSPHERE to the real Draco data. We find that Draco has an inner dark matter density of ρ _DM(150 pc) = 2.4_{-0.6}^{+0.5} × 10^8 M_☉ kpc^{-3}, consistent with a ΛCDM cusp. Using a velocity-independent SIDM model, calibrated on ΛSIDM cosmological simulations, we show that Draco's high central density gives an upper bound on the SIDM cross-section of σ/m < 0.57 cm2 g-1 at 99 per cent confidence. We conclude that the inner density of nearby dwarf galaxies like Draco provides a new and competitive probe of dark matter models.

We study the impact of stellar feedback in shaping the density and velocity structure of neutral hydrogen (H I) in disc galaxies. For our analysis, we carry out ∼4.6 pc resolution N-body+adaptive mesh refinement hydrodynamic simulations of isolated galaxies, set up to mimic a Milky Way and a Large and Small Magellanic Cloud. We quantify the density and velocity structure of the interstellar medium using power spectra and compare the simulated galaxies to observed H I in local spiral galaxies from THINGS (The H I Nearby Galaxy Survey). Our models with stellar feedback give an excellent match to the observed THINGS H I density power spectra. We find that kinetic energy power spectra in feedback-regulated galaxies, regardless of galaxy mass and size, show scalings in excellent agreement with supersonic turbulence (E(k) ∝ k-2) on scales below the thickness of the H I layer. We show that feedback influences the gas density field, and drives gas turbulence, up to large (kpc) scales. This is in stark contrast to density fields generated by large-scale gravity-only driven turbulence. We conclude that the neutral gas content of galaxies carries signatures of stellar feedback on all scales.

We apply the GC3 stream-finding method to RR Lyrae stars (RRLSs) in the Catalina survey. We find 2 RRLS stream candidates at >4σ confidence and another 12 at >3.5σ confidence over the Galactocentric distance range 4 < D/kpc < 26. Of these, only two are associated with known globular clusters (NGC 1261 and Arp2). The remainder are candidate `orphan' streams, consistent with the idea that globular cluster streams are most visible close to dissolution. Our detections are likely a lower bound on the total number of dissolving globulars in the inner galaxy, since many globulars have few RRLSs, while only the brightest streams are visible over the Galactic RRLS background, particularly given the current lack of kinematical information. We make all of our candidate streams publicly available and provide a new galstreamsPYTHON library for the footprints of all known streams and overdensities in the Milky Way.

According to star formation histories (SFHs), Local Group dwarf galaxies can be broadly classified in two types: those forming most of their stars before z = 2 (fast) and those with more extended SFHs (slow). The most precise SFHs are usually derived from deep but not very spatially extended photometric data; this might alter the ratio of old to young stars when age gradients are present. Here, we correct for this effect and derive the mass formed in stars by z = 2 for a sample of 16 Local Group dwarf galaxies. We explore early differences between fast and slow dwarfs, and evaluate the impact of internal feedback by supernovae (SNe) on the baryonic and dark matter (DM) component of the dwarfs. Fast dwarfs assembled more stellar mass at early times and have larger amounts of DM within the half-light radius than slow dwarfs. By imposing that slow dwarfs cannot have lost their gas by z = 2, we constrain the maximum coupling efficiency of SN feedback to the gas and to the DM to be ̃10 per cent. We find that internal feedback alone appears insufficient to quench the SFH of fast dwarfs by gas deprivation, in particular for the fainter systems. Nonetheless, SN feedback can core the DM halo density profiles relatively easily, producing cores of the sizes of the half-light radius in fast dwarfs by z = 2 with very low efficiencies. Amongst the `classical' Milky Way satellites, we predict that the smallest cores should be found in Draco and Ursa Minor, while Sculptor and Fornax should host the largest ones.

The shape and wide diversity of dwarf galaxy rotation curves is at apparent odds with dark matter halos in LCDM. We generate mock rotation curve data from dwarf galaxy simulations to show that this owes to bursty star formation driven by stellar feedback. There are three main effects. Firstly, stellar feedback transforms dark matter cusps into cores. Ignoring such transformations leads to a poor fit of the rotation curve shape and a large systematic bias on the halo concentration parameter c. Secondly, if close to a recent starburst, large HI bubbles push the rotation curve out of equilibrium. This makes the gas rotational velocity a poor probe of the underlying potential, leading to a systematic error on the halo virial mass M200 of up to half a dex. Thirdly, when galaxies are viewed near face-on (i

Coalescing massive black hole binaries, formed during galaxy mergers, are expected to be a primary source of low-frequency gravitational waves. Yet in isolated gas-free spherical stellar systems, the hardening of the binary stalls at parsec-scale separations owing to the inefficiency of relaxation-driven loss-cone refilling. Repopulation via collisionless orbit diffusion in triaxial systems is more efficient, but published simulation results are contradictory. While sustained hardening has been reported in simulations of galaxy mergers with N ∼ 106 stars and in early simulations of rotating models, in isolated non-rotating triaxial models the hardening rate continues to fall with increasing N, a signature of spurious two-body relaxation. We present a novel approach for studying loss-cone repopulation in galactic nuclei. Since loss-cone repopulation in triaxial systems owes to orbit diffusion, it is a purely collisionless phenomenon and can be studied with an approximated force calculation technique, provided the force errors are well behaved and sufficiently small. We achieve this using an accurate fast multipole method and define a proxy for the hardening rate that depends only on stellar angular momenta. We find that the loss cone is efficiently replenished even in very mildly triaxial models (with axis ratios 1:0.9:0.8). Such triaxiality is unavoidable following galactic mergers and can drive binaries into the gravitational wave regime. We conclude that there is no ‘final parsec problem’.

The coalescence of massive black hole binaries (BHBs) in galactic mergers is the primary source of gravitational waves (GWs) at low frequencies. Current estimates of GW detection rates for the Laser Interferometer Space Antenna and the Pulsar Timing Array vary by three orders of magnitude. To understand this variation, we simulate the merger of equal-mass, eccentric, galaxy pairs with central massive black holes and shallow inner density cusps. We model the formation and hardening of a central BHB using the Fast Multiple Method as a force solver, which features a O¹Nº scaling with the number N of particles and obtains results equivalent to direct-summation simulations. At N 5105, typical for contemporary studies, the eccentricity of the BHBs can vary significantly for different random realisations of the same initial condition, resulting in a substantial variation of the merger timescale. This scatter owes to the stochasticity of stellar encounters with the BHB and decreases with increasing N. We estimate that N 107 within the stellar half-light radius suffices to reduce the scatter in the merger timescale to 10%. Our results suggest that at least some of the uncertainty in low-frequency GW rates owes to insufficient numerical resolution.

Massive black hole (MBH) binaries, formed as a result of galaxy mergers, are expected to harden by dynamical friction and three-body stellar scatterings, until emission of gravitational waves (GWs) leads to their final coalescence. According to recent simulations, MBH binaries can efficiently harden via stellar encounters only when the host geometry is triaxial, even if only modestly, as angular momentum diffusion allows an efficient repopulation of the binary loss cone. In this paper, we carry out a suite of N-body simulations of equal-mass galaxy collisions, varying the initial orbits and density profiles for the merging galaxies and running simulations both with and without central MBHs. We find that the presence of an MBH binary in the remnant makes the system nearly oblate, aligned with the galaxy merger plane, within a radius enclosing 100 MBH masses. We never find binary hosts to be prolate on any scale. The decaying MBHs slightly enhance the tangential anisotropy in the centre of the remnant due to angular momentum injection and the slingshot ejection of stars on nearly radial orbits. This latter effect results in about 1% of the remnant stars being expelled from the galactic nucleus. Finally, we do not find any strong connection between the remnant morphology and the binary hardening rate, which depends only on the inner density slope of the remnant galaxy. Our results suggest that MBH binaries are able to coalesce within a few Gyr, even if the binary is found to partially erase the merger-induced triaxiality from the remnant.

The origin of the gas in between the Magellanic Clouds (MCs), known as the Magellanic Bridge, has always been the subject of controversy. To shed light into this, we present the results from the MAGellanic Inter-Cloud II (MAGIC II) project aimed at probing the stellar populations in 10 large fields located perpendicular to the main ridge-line of H i in the Inter-Cloud region. We secured these observations of the stellar populations in between the MCs using the WFI (Wide Field Imager) camera on the 2.2 m telescope in La Silla. Using colour–magnitude diagrams, we trace stellar populations across the Inter-Cloud region. In good agreement with MAGIC I, we find significant intermediate-age stars in the Inter-Cloud region as well as young stars of a similar age to the last pericentre passage in between the MCs (∼200 Myr ago). We show here that the young, intermediate-age and old stars have distinct spatial distributions. The young stars correlate well with the H i gas suggesting that they were either recently stripped from the Small Magellanic Cloud (SMC) or formed in situ. The bulk of intermediate-age stars are located mainly in the Bridge region where the H i column density is higher, but they are more spread out than the young stars. They have very similar properties to stars located ∼2 kpc from the SMC centre, suggesting that they were tidally stripped from this region. Finally, the old stars extend to some 8 kpc from the SMC supporting the idea that all galaxies have a large extended metal-poor stellar halo.

We present a new method for determining the local dark matter density using kinematic data for a population of tracer stars. The Jeans equation in the z-direction is integrated to yield an equation that gives the velocity dispersion as a function of the total mass density, tracer density, and the ‘tilt’ term that describes the coupling of vertical and radial motions. We then fit a dark matter mass profile to tracer density and velocity dispersion data to derive credible regions on the vertical dark matter density profile. Our method avoids numerical differentiation, leading to lower numerical noise, and is able to deal with the tilt term while remaining one dimensional. In this study we present the method and perform initial tests on idealised mock data. We also demonstrate the importance of dealing with the tilt term for tracers that sample > ∼ 1 kpc above the disc plane. If ignored, this results in a systematic underestimation of the dark matter density.

We simulate star formation in two molecular clouds extracted from a larger disc-galaxy simulation with a spatial resolution of ∼0.1 pc, one exiting a spiral arm dominated by compression, and another in an inter-arm region more strongly affected by galactic shear. Treating the stars as ‘sink particles’, we track their birth angular momentum, and the later evolution of their angular momentum due to gas accretion. We find that in both clouds, the sinks have spin vectors that are aligned with one another, and with the global angular momentum vector of the star cluster. This alignment is present at birth, but enhanced by later gas accretion. In the compressive cloud, the sink-spins remain aligned with the gas for at least a free-fall time. By contrast, in the shear cloud, the increased turbulent mixing causes the sinks to rapidly misalign with their birth cloud on approximately a gas free-fall time. In spite of this, both clouds show a strong alignment of sink-spins at the end of our simulations, independently of environment.

Due to their close proximity, the Large and Small Magellanic Clouds (LMC/SMC) provide natural laboratories for understanding how galaxies form and evolve. With the goal of determining the structure and dynamical state of the SMC, we present new spectroscopic data for ∼3000 SMC red giant branch stars observed using the AAOmega spectrograph at the Anglo-Australian Telescope. We complement our data with further spectroscopic measurements from previous studies that used the same instrumental configuration as well as proper motions from the Gaia Data Release 2 catalogue. Analysing the photometric and stellar kinematic data, we find that the SMC centre of mass presents a conspicuous offset from the velocity centre of its associated H i gas, suggesting that the SMC gas is likely to be far from dynamical equilibrium. Furthermore, we find evidence that the SMC is currently undergoing tidal disruption by the LMC within 2 kpc of the centre of the SMC, and possibly all the way into the very core. This is revealed by a net outward motion of stars from the SMC centre along the direction towards the LMC and an apparent tangential anisotropy at all radii. The latter is expected if the SMC is undergoing significant tidal stripping, as we demonstrate using a suite of N-body simulations of the SMC/LMC system disrupting around the Milky Way. Our results suggest that dynamical models for the SMC that assume a steady state will need to be revisited.

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 (Reff ˃ 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 (⟨Reff⟩ ≃ 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 (⟨Reff⟩ ≃ 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.

Dwarf irregular galaxies (dIrrs) are the smallest stellar systems with extended H I discs. The study of the kinematics of such discs is a powerful tool to estimate the total matter distribution at these very small scales. In this work, we study the H I kinematics of 17 galaxies extracted from the ‘Local Irregulars That Trace Luminosity Extremes, The H I Nearby Galaxy Survey’ (LITTLE THINGS). Our approach differs significantly from previous studies in that we directly fit 3D models (two spatial dimensions plus one spectral dimension) using the software 3DBAROLO, fully exploiting the information in the H I data cubes. For each galaxy, we derive the geometric parameters of the H I disc (inclination and position angle), the radial distribution of the surface density, the velocity-dispersion (σv) profile and the rotation curve. The circular velocity (Vc), which traces directly the galactic potential, is then obtained by correcting the rotation curve for the asymmetric drift. As an initial application, we show that these dIrrs lie on a baryonic Tully–Fisher relation in excellent agreement with that seen on larger scales. The final products of this work are high-quality, ready-to-use kinematic data (Vc and σv) that we make publicly available. These can be used to perform dynamical studies and improve our understanding of these low-mass galaxies.

We present new FLAMES+GIRAFFE spectroscopy of 36 member stars in the isolated Local Group dwarf spheroidal galaxy Tucana. We measure a systemic velocity for the system of vTuc=216.7+2.9−2.8vTuc=216.7−2.8+2.9 km s−1, and a velocity dispersion of σv,Tuc=14.4+2.8−2.3σv,Tuc=14.4−2.3+2.8 km s−1. We also detect a rotation gradient of dvrdχ=7.6+4.2−4.3dvrdχ=7.6−4.3+4.2 km s−1 kpc−1, which reduces the systemic velocity to vTuc=215.2+2.8−2.7vTuc=215.2−2.7+2.8 km s−1 and the velocity dispersion to σv,Tuc=13.3+2.7−2.3σv,Tuc=13.3−2.3+2.7 km s−1. We perform Jeans modelling of the density profile of Tucana, using the line-of-sight velocities of the member stars. We find that it favours a high central density consistent with ‘pristine’ subhaloes in Λ cold dark matter, and a massive dark matter halo (∼1010 M⊙) consistent with expectations from abundance matching. Tucana appears to be significantly more centrally dense than other isolated Local Group dwarfs, making it an ideal laboratory for testing dark matter models.

We make the first attempt to find dwarf galaxies in eight Fermi-LAT extended, unassociated, source fields using Gaia DR2. After probing previously unexplored heliocentric distances of d ˂ 20 kpc with an extreme-deconvolution (XD) technique, we find no sign of a dwarf galaxy in any of these fields despite Gaia's excellent astrometric accuracy. Our detection limits are estimated by applying the XD method to mock data, obtaining a conservative limit on the stellar mass of M* ˂ 104 M☉ for d ˂ 20 kpc. Such a low stellar mass implies either a low-mass subhalo or a massive stripped-down subhalo. We use an analytic model for stripped subhaloes to argue that, given the sizes and fluxes of the Fermi-LAT sources, we can reject the hypothesis that they owe to dark matter annihilation.

We use high resolution simulations of isolated dwarf galaxies to study the physics of dark matter cusp-core transformations at the edge of galaxy formation: M200 = 107 109M .We work at a resolution ( 4 pc minimum cell size; 250M per particle) at which the impact from individual supernovae explosions can be resolved, becoming insensitive to even large changes in our numerical `sub-grid' parameters. We nd that our dwarf galaxies give a remarkable match to the stellar light pro le; star formation history; metallicity distribution function; and star/gas kinematics of isolated dwarf irregular galaxies. Our key result is that dark matter cores of size comparable to the stellar half mass radius r1=2 always form if star formation proceeds for long enough. Cores fully form in less than 4 Gyrs for the M200 = 108M and 14 Gyrs for the 109M dwarf. We provide a convenient two parameter `coreNFW' tting function that captures this dark matter core growth as a function of star formation time and the projected stellar half mass radius. Our results have several implications: (i) we make a strong prediction that if CDM is correct, then `pristine' dark matter cusps will be found either in systems that have truncated star formation and/or at radii r > r1=2; (ii) complete core formation lowers the projected velocity dispersion at r1=2 by a factor 2, which is su cient to fully explain the `too big to fail problem'; and (iii) cored dwarfs will be much more susceptible to tides, leading to a dramatic scouring of the subhalo mass function inside galaxies and groups.

Dark matter-only simulations of galaxy formation predict many more subhalos around a Milky Way-like galaxy than the number of observed satellites. Proposed solutions require the satellites to inhabit dark matter halos with masses 10(9)-10(10 )Msun at the time they fell into the Milky Way. Here we use a modelling approach, independent of cosmological simulations, to obtain a pre-infall mass of 3.6(-2.3)(+3.8) × 10(8) Msun for one of the Milky Way's satellites: Carina. This determination of a low halo mass for Carina can be accommodated within the standard model only if galaxy formation becomes stochastic in halos below ∼10(10 )Msun. Otherwise Carina, the eighth most luminous Milky Way dwarf, would be expected to inhabit a significantly more massive halo. The implication of this is that a population of 'dark dwarfs' should orbit the Milky Way: halos devoid of stars and yet more massive than many of their visible counterparts.

High-mass galaxies, with halo masses M200 ≥ 1010M⊙ , reveal a remarkable near-linear relation between their globular cluster (GC) system mass and their host galaxy halo mass. Extending this relation to the mass range of dwarf galaxies has been problematic due to the difficulty in measuring independent halo masses. Here we derive new halo masses based on stellar and H i gas kinematics for a sample of nearby dwarf galaxies with GC systems. We find that the GC system mass–halo mass relation for galaxies populated by GCs holds from halo masses of M200 ∼ 1014 M ⊙ ⊙ down to below M200 ∼109 M ⊙ ⊙ , although there is a substantial increase in scatter towards low masses. In particular, three well-studied ultradiffuse galaxies, with dwarf-like stellar masses, reveal a wide range in their GC-to-halo mass ratios. We compare our GC system–halo mass relation to the recent model of El Badry et al., finding that their fiducial model does not reproduce our data in the low-mass regime. This may suggest that GC formation needs to be more efficient than assumed in their model, or it may be due to the onset of stochastic GC occupation in low-mass haloes. Finally, we briefly discuss the stellar mass–halo mass relation for our low-mass galaxies with GCs, and we suggest some nearby dwarf galaxies for which searches for GCs may be fruitful.

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 imply 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.

Andromeda (And) XXV has previously been reported as a dwarf spheroidal galaxy (dSph) with little-to-no dark matter. However, the uncertainties on this result were significant. In this study, we double the number of member stars and re-derive the kinematics and mass of And XXV. We find that And XXV has a systemic velocity of $\nu_\mathrm{r}=-107.7\pm1.0 \mathrm{~km s}^{-1}$ and a velocity dispersion of $\sigma_\nu=4.5\pm1.0\mathrm{~km s}^{-1}$. With this better constrained velocity dispersion, we derive a mass contained within the half-light radius of $M(r< r_\mathrm{h})=6.9^{+3.2}_{-2.8}\times10^6\mathrm{~M}_\odot$. This mass corresponds to a mass-to-light ratio of $\mathrm{[M/L]}_\mathrm{r_\mathrm{h}}=37^{+17}_{-15}\mathrm{~M}_\odot/\mathrm{L}_\odot$, demonstrating, for the first time, that And XXV has an unambiguous dark matter component. We also measure the metallicity of And XXV to be $\mathrm{[Fe/H]}=-1.9\pm0.1$$\mathrm{~}$dex, which is in agreement with previous results. Finally, we extend the analysis of And XXV to include mass modelling using GravSphere. We find that And XXV has a low central dark matter density, $\rho_\mathrm{DM}(150\mathrm{pc})= 2.7^{+1.8}_{-1.6}\times10^7\mathrm{~M}_\odot\mathrm{kpc}^{-3}$, making And XXV a clear outlier when compared to other Local Group (LG) dSphs of the similar stellar mass. In a companion paper, we will explore whether some combination of dark matter cusp-core transformations and/or tides can explain And XXV's low density.

We apply four different mass modelling methods to a suite of publicly available mock data for spherical stellar systems. We focus on the recovery of the density and velocity anisotropy as a function of radius, either using line-of-sight velocity data only or adding proper motion data. All methods perform well on isotropic and tangentially anisotropic mock data, recovering the density and velocity anisotropy within their 95 per cent confidence intervals over the radial range 0.25 < R/R1/2 < 4, where R1/2 is the half-light radius. However, radially anisotropic mocks are more challenging. For line-of-sight data alone, only methods that use information about the shape of the velocity distribution function are able to break the degeneracy between the density profile and the velocity anisotropy, β, to obtain an unbiased estimate of both. This shape information can be obtained through directly fitting a global phase-space distribution function, by using higher order ‘virial shape parameters’ or by assuming a Gaussian velocity distribution function locally, but projecting it self-consistently along the line of sight. Including proper motion data yields further improvements, and in this case, all methods give a good recovery of both the radial density and velocity anisotropy profiles.

ABSTRACT The Eridanus II (EriII) ‘ultra-faint’ dwarf has a large (15 pc) and low-mass (4.3 × 103 M⊙) star cluster (SC) offset from its centre by 23 ± 3 pc in projection. Its size and offset are naturally explained if EriII has a central dark matter core, but such a core may be challenging to explain in a ΛCDM cosmology. In this paper, we revisit the survival and evolution of EriII’s SC, focusing for the first time on its puzzlingly large ellipticity ($0.31^{+0.05}_{-0.06}$). We perform a suite of 960 direct N-body simulations of SCs, orbiting within a range of spherical background potentials fit to ultra-faint dwarf (UFD) galaxy simulations. We find only two scenarios that come close to explaining EriII’s SC. In the first scenario, EriII has a low-density dark matter core (of size ${\sim}70\, \text{pc}$ and density $\lesssim 2\times 10^8\, \text{M}_{\odot }\, \text{kpc}^{-3}$). In this model, the high ellipticity of EriII’s SC is set at birth, with the lack of tidal forces in the core allowing its ellipticity to remain frozen for long times. In the second scenario, EriII’s SC orbits in a partial core, with its high ellipticity owing to its imminent tidal destruction. However, this latter model struggles to reproduce the large size of EriII’s SC, and it predicts substantial tidal tails around EriII’s SC that should have already been seen in the data. This leads us to favour the cored model. We discuss potential caveats to these findings, and the implications of the cored model for galaxy formation and the nature of dark matter.

Feedback to the interstellar medium (ISM) from ionising radiation, stellar winds and supernovae is central to regulating star formation in galaxies. Due to their low mass (M * < 10 9 M), dwarf galaxies are particularly susceptible to such processes, making them ideal sites to study the detailed physics of feedback. In this perspective, we summarise the latest observational evidences for feedback from star forming regions and how this drives the formation of 'superbubbles' and galaxy-wide winds. We discuss the important role of external ionising radiation – 'reionisation' – for the smallest galaxies. And, we discuss the observational evidences that this feedback directly impacts galaxy properties such as their star formation histories, metal content, colours, sizes, morphologies and even their inner dark matter densities. We conclude with a look to the future, summarising the key questions that remain unanswered and listing some of the outstanding challenges for galaxy formation theories.

Massive elliptical galaxies are typically observed to have central cores in their projected radial light profiles. Such cores have long been thought to form through ‘binary scouring’ as supermassive black holes (SMBHs), brought in through mergers, form a hard binary and eject stars from the galactic centre. However, the most massive cores, like the $\sim 3{\, \mathrm{kpc}}$ core in A2261-BCG, remain challenging to explain in this way. In this paper, we run a suite of dry galaxy merger simulations to explore three different scenarios for central core formation in massive elliptical galaxies: ‘binary scouring’, ‘tidal deposition’, and ‘gravitational wave (GW) induced recoil’. Using the griffin code, we self-consistently model the stars, dark matter, and SMBHs in our merging galaxies, following the SMBH dynamics through to the formation of a hard binary. We find that we can only explain the large surface brightness core of A2261-BCG with a combination of a major merger that produces a small $\sim 1{\, \mathrm{kpc}}$ core through binary scouring, followed by the subsequent GW recoil of its SMBH that acts to grow the core size. Key predictions of this scenario are an offset SMBH surrounded by a compact cluster of bound stars and a non-divergent central density profile. We show that the bright ‘knots’ observed in the core region of A2261-BCG are best explained as stalled perturbers resulting from minor mergers, though the brightest may also represent ejected SMBHs surrounded by a stellar cloak of bound stars.

Andromeda XXI (And XXI) has been proposed as a dwarf spheroidal galaxy with a central dark matter density that is lower than expected in the standard $\Lambda$ cold dark matter ($\Lambda$CDM) cosmology. In this work, we present dynamical observations for 77 member stars in this system, more than doubling previous studies to determine whether this galaxy is truly a low-density outlier. We measure a systemic velocity of $v_r=-363.4\pm 1.0{\rm \, km\, s^{-1}}$ and a velocity dispersion of $\sigma _v=6.1^{+1.0}_{-0.9}{\rm \, km\, s^{-1}}$, consistent with previous work and within $1\sigma$ of predictions made using the modified Newtonian dynamics framework. We also measure the metallicity of our member stars from their spectra, finding a mean value of ${\rm [Fe/H]}=-1.7\pm 0.1$ dex. We model the dark matter density profile of And XXI using an improved version of gravsphere, finding a central density of $\rho _{\rm DM}({\rm 150 pc})=2.6_{-1.5}^{+2.4} \times 10^7 \, {\rm M_\odot \, kpc^{-3}}$ at 68 per cent confidence, and a density at two half-light radii of $\rho _{\rm DM}({\rm 1.75 kpc})=0.9_{-0.2}^{+0.3} \times 10^6 \, {\rm M_\odot \, kpc^{-3}}$ at 68 per cent confidence. These are both a factor of${\sim }3\!-\!5$ lower than the densities expected from abundance matching in $\Lambda$CDM. We show that this cannot be explained by ‘dark matter heating’ since And XXI had too little star formation to significantly lower its inner dark matter density, while dark matter heating only acts on the profile inside the half-light radius. However, And XXI’s low density can be accommodated within $\Lambda$CDM if it experienced extreme tidal stripping (losing ${\gt}95{{\ \rm per\ cent}}$ of its mass), or if it inhabits a low-concentration halo on a plunging orbit that experienced repeated tidal shocks.

ABSTRACT Motivated by the observed bottom-light initial mass function (IMF) in faint dwarfs, we study how a metallicity-dependent IMF affects the feedback budget and observables of an ultra-faint dwarf galaxy. We model the evolution of a low-mass ($\approx 8 \, \times \, 10^{8} \, \rm M_{\odot }$) dark matter halo with cosmological, zoomed hydrodynamical simulations capable of resolving individual supernovae explosions, which we complement with an empirically motivated subgrid prescription for systematic IMF variations. In this framework, at the low gas metallicities typical of faint dwarfs, the IMF of newborn stellar populations becomes top-heavy, increasing the efficiency of supernova and photoionization feedback in regulating star formation. This results in a 100-fold reduction of the final stellar mass of the dwarf compared to a canonical IMF, at fixed dynamical mass. The increase in the feedback budget is none the less met by increased metal production from more numerous massive stars, leading to nearly constant iron content at z = 0. A metallicity-dependent IMF therefore provides a mechanism to produce low-mass ($\rm M_{\star }\sim 10^3 \rm M_{\odot }$), yet enriched ($\rm [Fe/H]\approx -2$) field dwarf galaxies, thus opening a self-consistent avenue to populate the plateau in $\rm [Fe/H]$ at the faintest end of the mass–metallicity relation.

We show how the interplay between feedback and mass-growth histories introduces scatter in the relationship between stellar and neutral gas properties of field faint dwarf galaxies (M-*less than or similar to 10(6) M-circle dot). Across a suite of cosmological, high-resolution zoomed simulations, we find that dwarf galaxies of stellar masses 10(5)

The density of dark matter near the Sun, ρ DM, ⊙ , is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way’s dark matter halo. Estimates to date have typically assumed that the Milky Way’s stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. We assess the impact of these assumptions on determinations of ρ DM, ⊙  by applying a free-form, steady-state, Jeans method to two different  N -body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized  Gaia –Sausage–Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans–Poisson equations on our ability to correctly extract ρ DM, ⊙  from the simulated data. We find that common approximations employed in the literature – axisymmetry and a locally flat rotation curve – can lead to significant systematic errors of up to a factor ∼1.5 in the recovered surface mass density ∼2 kpc above the disc plane, implying a fractional error on ρ DM, ⊙  of the order of unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρ DM, ⊙ , consistent with the true value within our 95 per cent confidence intervals for realistic 20 per cent uncertainties on the baryonic surface density of the disc. Other terms – the axial tilt, 2nd Poisson and time-dependent terms – contribute less than 10 per cent to ρ DM, ⊙  (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis.

The low dark matter density in the Fornax dwarf galaxy is often interpreted as being due to the presence of a constant density 'core', but it could also be explained by the effects of Galactic tides. The latter interpretation has been disfavoured because it is apparently inconsistent with the orbital parameters and star formation history of Fornax. We revisit these arguments with the help of the APOSTLE cosmological hydrodynamics simulations. We show that simulated dwarfs with similar properties to Fornax are able to form stars after infall, so that star formation is not necessarily a good tracer of infall time. We also examine the constraints on the pericentre of Fornax and point out that small pericentres (

We introduce a new method to mitigate numerical diffusion in adaptive mesh refinement (AMR) simulations of cosmological galaxy formation, and study its impact on a simulated dwarf galaxy as part of the 'EDGE' project. The target galaxy has a maximum circular velocity of 21 km s(-1) but evolves in a region that is moving at up to 90 km s(-1) relative to the hydrodynamic grid. In the absence of any mitigation, diffusion softens the filaments feeding our galaxy. As a result, gas is unphysically held in the circumgalactic medium around the galaxy for 320 Myr, delaying the onset of star formation until cooling and collapse eventually triggers an initial starburst at z = 9. Using genetic modification, we produce 'velocity-zeroed' initial conditions in which the grid-relative streaming is strongly suppressed; by design, the change does not significantly modify the large-scale structure or dark matter accretion history. The resulting simulation recovers a more physical, gradual onset of star formation starting at z = 17. While the final stellar masses are nearly consistent (4.8 x 10(6) M-circle dot and 4.4 x 10(6) M-circle dot for unmodified and velocity-zeroed, respectively), the dynamical and morphological structure of the z = 0 dwarf galaxies are markedly different due to the contrasting histories. Our approach to diffusion suppression is suitable for any AMR zoom cosmological galaxy formation simulations, and is especially recommended for those of small galaxies at high redshift.

Dwarf spheroidal galaxies are excellent systems to probe the nature of fermionic dark matter due to their high observed dark matter phase-space density. In this work, we review, revise, and improve upon previous phase-space considerations to obtain lower bounds on the mass of fermionic dark matter particles. The refinement in the results compared to previous works is realized particularly due to a significantly improved Jeans analysis of the galaxies. We discuss two methods to obtain phase-space bounds on the dark matter mass, one model-independent bound based on Pauli's principle, and the other derived from an application of Liouville's theorem. As benchmark examples for the latter case, we derive constraints for thermally decoupled particles and (non-)resonantly produced sterile neutrinos. Using the Pauli principle, we report a model-independent lower bound of m >= 0.18 keV at 68 per cent CL and m >= 0.13 keV at 95 per cent CL. For relativistically decoupled thermal relics, this bound is strengthened to m >= 0.59 keV at 68 per cent CL and m >= 0.41 keV at 95 per cent CL, while for non-resonantly produced sterile neutrinos the constraint is m >= 2.80 keV at 68 per cent CL and m >= 1.74 keV at 95 per cent CL. Finally, the phase-space bounds on resonantly produced sterile neutrinos are compared with complementary limits from X-ray, Lyman alpha, and big bang nucleosynthesis observations.

In the standard Lambda cold dark matter paradigm, pure dark matter simulations predict dwarf galaxies should inhabit dark matter haloes with a centrally diverging density 'cusp'. This is in conflict with observations that typically favour a constant density 'core'. We investigate this 'cusp-core problem' in 'ultra-faint' dwarf galaxies simulated as part of the 'Engineering Dwarfs at Galaxy formation's Edge' project. We find, similarly to previous work, that gravitational potential fluctuations within the central region of the simulated dwarfs kinematically heat the dark matter particles, lowering the dwarfs' central dark matter density. However, these fluctuations are not exclusively caused by gas inflow/outflow, but also by impulsive heating from minor mergers. We use the genetic modification approach on one of our dwarf's initial conditions to show how a delayed assembly history leads to more late minor mergers and, correspondingly, more dark matter heating. This provides a mechanism by which even ultra-faint dwarfs (M-star < 10(5) M-circle dot), in which star formation was fully quenched at high redshift, can have their central dark matter density lowered over time. In contrast, we find that late major mergers can regenerate a central dark matter cusp, if the merging galaxy had sufficiently little star formation. The combination of these effects leads us to predict significant stochasticity in the central dark matter density slopes of the smallest dwarfs, driven by their unique star formation and mass assembly histories.

ABSTRACT Using the cosmological zoom simulation VINTERGATAN, we present a new scenario for the onset of star formation at the metal-poor end of the low-[α/Fe] sequence in a Milky Way-like galaxy. In this scenario, the galaxy is fuelled by two distinct gas flows. One is enriched by outflows from massive galaxies, but not the other. While the former feeds the inner galactic region, the latter fuels an outer gas disc, inclined with respect to the main galactic plane, and with a significantly poorer chemical content. The first passage of the last major merger galaxy triggers tidal compression in the outer disc, which increases the gas density and eventually leads to star formation, at a metallicity 0.75 dex lower than the inner galaxy. This forms the first stars of the low-[α/Fe] sequence. These in situ stars have halo-like kinematics, similar to what is observed in the Milky Way, due to the inclination of the outer disc that eventually aligns with the inner one via gravitational torques. We show that this tilting disc scenario is likely to be common in Milky Way-like galaxies. This process implies that the low-[α/Fe] sequence is populated in situ, simultaneously from two formation channels, in the inner and the outer galaxy, with distinct metallicities. This contrasts with purely sequential scenarios for the assembly of the Milky Way disc and could be tested observationally.

ABSTRACT Using the VINTERGATAN cosmological zoom simulation, we explore the contributions of the in situ and accreted material, and the effect of galaxy interactions and mergers in the assembly of a Milky Way-like galaxy. We find that the initial growth phase of galaxy evolution, dominated by repeated major mergers, provides the necessary physical conditions for the assembly of a thick, kinematically hot disc populated by high-[α/Fe] stars, formed both in situ and in accreted satellite galaxies. We find that the diversity of evolutionary tracks followed by the simulated galaxy and its progenitors leads to very little overlap of the in situ and accreted populations for any given chemical composition. At a given age, the spread in [α/Fe] abundance ratio results from the diversity of physical conditions in VINTERGATAN and its satellites, with an enhancement in [α/Fe] found in stars formed during starburst episodes. Later, the cessation of the merger activity promotes the in situ formation of stars in the low-[α/Fe] regime, in a radially extended, thin and overall kinematically colder disc, thus establishing chemically bimodal thin and thick discs, in line with observations. We draw links between notable features in the [Fe/H]-[α/Fe] plane with their physical causes, and propose a comprehensive formation scenario explaining self-consistently, in the cosmological context, the main observed properties of the Milky Way.

ABSTRACT We test a non-parametric higher order Jeans analysis method, GravSphere, on 32 simulated dwarf galaxies comparable to classical Local Group dwarfs like Fornax. The galaxies are selected from A Project Of Simulating The Local Environment (APOSTLE) suite of cosmological hydrodynamics simulations with cold dark matter (CDM) and self-interacting dark matter (SIDM) models, allowing us to investigate cusps and cores in density distributions. We find that, for CDM dwarfs, the recovered enclosed mass profiles have a bias of no more than 10 per cent, with a 50 per cent scatter in the inner regions and a 20 per cent scatter near the half-light radius, consistent with standard mass estimators. The density profiles are also recovered with a bias of no more than 10 per cent and a scatter of 30 per cent in the inner regions. For SIDM dwarfs, the mass and density profiles are recovered within our 95 per cent confidence intervals but are biased towards cuspy dark matter distributions. This is mainly due to a lack of sufficient constraints from the data. We explore the sources of scatter in the accuracy of the recovered profiles and suggest a χ2 statistic to separate successful models from biased ones. Finally, we show that the uncertainties on the mass profiles obtained with GravSphere are smaller than those for comparable Jeans methods and that they can be further improved if stronger priors, motivated by cosmological simulations, are placed on the velocity anisotropy. We conclude that GravSphere is a promising Jeans-based approach for modelling dark matter distributions in dwarf galaxies.

We present an updated analysis of the gamma-ray flux from the directions of classical dwarf spheroidal galaxies, deriving new constraints on WIMP dark matter (DM) annihilation using a decade of Fermi-LAT data. Among the major novelties, we infer the dwarfs' J-factors by including new observations without imposing any a priori parametric profile for the DM distribution. While statistically compatible with results obtained from more conventional parameterisations, this procedure reduces the theoretical bias imposed on the data. Furthermore, we retain the full data-driven shape of the J-factors' empirical probability distributions when setting limits on DM, without imposing log-normality as is typically done. In conjunction with the data-driven J-factors, we improve on a new method for estimating the probability distribution function of the astrophysical background at the dwarf position [1], fully profiling over background uncertainties. We show that, for most "classical" dwarfs, the background systematic uncertainty dominates over the uncertainty on their J-factors. Raw distributions of J- and D-factors (the latter being the analogous of J-factors for decaying DM) are available upon request.

The coalescence of massive black hole binaries (BHBs) in galactic mergers is the primary source of gravitational waves (GWs) at low frequencies. Current estimates of GW detection rates for the Laser Interferometer Space Antenna and the Pulsar Timing Array vary by three orders of magnitude. To understand this variation, we simulate the merger of equal-mass, eccentric, galaxy pairs with central massive black holes and shallow inner density cusps. We model the formation and hardening of a central BHB using the fast multiple method as a force solver, which features a O(N) scaling with the number N of particles and obtains results equivalent to direct-summation simulations. At N similar to 5 x 10(5), typical for contemporary studies, the eccentricity of the BHBs can vary significantly for different random realizations of the same initial condition, resulting in a substantial variation of the merger time-scale. This scatter owes to the stochasticity of stellar encounters with the BHB and decreases with increasing N. We estimate that N similar to 10(7) within the stellar half-light radius suffices to reduce the scatter in the merger time-scale to similar to 10 per cent. Our results suggest that at least some of the uncertainty in low-frequency GW rates owes to insufficient numerical resolution.

ABSTRACT Spectroscopic surveys of the Milky Way’s stars have revealed spatial, chemical, and kinematical structures that encode its history. In this work, we study their origins using a cosmological zoom simulation, VINTERGATAN, of a Milky Way-mass disc galaxy. We find that in connection to the last major merger at z ∼ 1.5, cosmological accretion leads to the rapid formation of an outer, metal-poor, low-[α/Fe] gas disc around the inner, metal-rich galaxy containing the old high-[α/Fe] stars. This event leads to a bimodality in [α/Fe] over a range of [Fe/H]. A detailed analysis of how the galaxy evolves since z ∼ 1 is presented. We demonstrate the way in which inside-out growth shapes the radial surface density and metallicity profile and how radial migration preferentially relocates stars from the inner disc to the outer disc. Secular disc heating is found to give rise to increasing velocity dispersions and scale heights with stellar age, which together with disc flaring explains several trends observed in the Milky Way, including shallower radial [Fe/H] profiles above the mid-plane. We show how the galaxy formation scenario imprints non-trivial mappings between structural associations (i.e. thick and thin discs), velocity dispersions, α-enhancements, and ages of stars; e.g. the most metal-poor stars in the low-[α/Fe] sequence are found to have a scale height comparable to old high-[α/Fe] stars. Finally, we illustrate how at low spatial resolution, comparable to the thickness of the galaxy, the proposed pathway to distinct sequences in [α/Fe]–[Fe/H] cannot be captured.

ABSTRACT We use the Milky Way’s nuclear star cluster (NSC) to test the existence of a dark matter ‘soliton core’, as predicted in ultra-light dark matter (ULDM) models. Since the soliton core size is proportional to $m_{\rm DM}^{-1}$, while the core density grows as $m_{\rm DM}^{2}$, the NSC (dominant stellar component within ∼3 pc) is sensitive to a specific window in the dark matter particle mass, mDM. We apply a spherical isotropic Jeans model to fit the NSC line-of-sight velocity dispersion data, assuming priors on the precisely measured Milky Way’s supermassive black hole (SMBH) mass and the well-measured NSC density profile. We find that the current observational data reject the existence of a soliton core for a single ULDM particle with mass in the range 10−20.4 eV ≲ mDM ≲ 10−18.5 eV, assuming that the soliton core structure is not affected by the Milky Way’s SMBH. We test our methodology on mock data, confirming that we are sensitive to the same range in ULDM mass as for the real data. Dynamical modelling of a larger region of the Galactic centre, including the nuclear stellar disc, promises tighter constraints over a broader range of mDM. We will consider this in future work.