Researcher profile

Dougal Mackey

Dougal Mackey contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2026arXiv

Rotational Kinematics in the Globular Cluster System of M31: Insights from Bayesian Inference

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.

preprint2022arXiv

$S^5$: The Orbital and Chemical Properties of One Dozen Stellar Streams

We report the kinematic, orbital, and chemical properties of 12 stellar streams with no evident progenitors, using line-of-sight velocities and metallicities from the Southern Stellar Stream Spectroscopic Survey ($S^5$), proper motions from $Gaia$ EDR3, and distances derived from distance tracers or the literature. This data set provides the largest homogeneously analyzed set of streams with full 6D kinematics and metallicities. All streams have heliocentric distances between ${\sim}10-50$ kpc. The velocity and metallicity dispersions show that half of the stream progenitors were disrupted dwarf galaxies (DGs), while the other half originated from disrupted globular clusters (GCs), hereafter referred to as DG and GC streams. Based on the mean metallicities of the streams and the mass-metallicity relation, the luminosities of the progenitors of the DG streams range between Carina and Ursa Major I ($-9.5\lesssim M_V\lesssim-5.5$). Four of the six GC streams have mean metallicities of [Fe/H]$< -2$, more metal-poor than typical Milky Way (MW) GCs at similar distances. Interestingly, the 300S and Jet GC streams are the only streams on retrograde orbits in our dozen stream sample. Finally, we compare the orbital properties of the streams with known DGs and GCs in the MW, finding several possible associations. Some streams appear to have been accreted with the recently discovered Gaia-Enceladus-Sausage system, and others suggest that GCs were formed in and accreted together with the progenitors of DG streams whose stellar masses are similar to Draco to Carina ($\sim10^5-10^6M_\odot$).

preprint2022arXiv

A search for stellar structures around nine outer halo globular clusters in the Milky Way

We use deep imaging from the Dark Energy Camera to explore the peripheral regions of nine globular clusters in the outer halo of the Milky Way. Apart from Whiting 1 and NGC 7492, which are projected against the Sagittarius stream, we see no evidence for adjacent stellar populations to indicate any of these clusters is associated with coherent tidal debris from a destroyed host dwarf. We also find no evidence for tidal tails around any of the clusters in our sample; however, both NGC 1904 and 6981 appear to possess outer envelopes. Motivated by a slew of recent Gaia-based discoveries, we compile a sample of clusters with robust detections of extra-tidal structure, and search for correlations with orbital properties. While we observe that clusters with tidal tails are typically on moderately or very eccentric orbits that are highly inclined to the Galactic plane and often retrograde, these are neither necessary nor sufficient conditions for the formation of extra-tidal structure. That many objects with tidal tails appear to be accreted leads us to speculate that this lack of consistency may stem from the inhomogeneous dynamical history of the Milky Way globular cluster system. Finally, we note that clusters with prominent stellar envelopes detected in ground-based imaging (such as NGC 1851 and 7089) are now all known from Gaia to possess long tidal tails -- experimental confirmation that the presence of an extended envelope is indicative of tidal erosion.

preprint2022arXiv

Gamma-ray emission from the Sagittarius Dwarf Spheroidal galaxy due to millisecond pulsars

The Fermi Bubbles are giant, gamma-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&#39;s super-massive black hole. Here we show via a spatial template analysis that much of the gamma-ray emission associated to the brightest region of substructure -- the so-called cocoon -- is likely due to the Sagittarius dwarf spheroidal (Sgr dSph) galaxy. 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 Sgr dSph has no on-going star formation, but we nevertheless demonstrate that the dwarf&#39;s millisecond pulsar (MSP) population can plausibly supply the gamma-ray signal that our analysis associates to 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 MSPs belonging to the Sgr dSph, combined with these objects&#39; magnetospheric emission. This finding plausibly suggests that MSPs produce significant gamma-ray emission amongst old stellar populations, potentially confounding indirect dark matter searches in regions such as the Galactic Centre, the Andromeda galaxy, and other massive Milky Way dwarf spheroidals.

preprint2022arXiv

The PAndAS View of the Andromeda Satellite System. III. Dwarf galaxy detection limits

We determine the detection limits of the search for dwarf galaxies in the Pan-Andromeda Archaeological Survey (PAndAS) using the algorithm developed by the PAndAS team. The recovery fractions of artificial dwarf galaxies are, as expected, a strong function of physical size and luminosity and, to a lesser extent, distance. We show that these recovery fractions vary strongly with location in the surveyed area because of varying levels of contamination from both the Milky Way foreground stars and the stellar halo of Andromeda. We therefore provide recovery fractions that are a function of size, luminosity, and location within the survey on a scale of 1 square degree. Overall, the effective surface brightness for a 50-percent detection rate range between 28 and 30 mag per square arcsecond. This is in line with expectations for a search that relies on photometric data that are as deep as the PAndAS survey. The derived detection limits are an essential ingredient on the path to constraining the global properties of Andromeda&#39;s system of satellite dwarf galaxies and, more broadly, to provide constraints on dwarf galaxy formation and evolution in a cosmological context.

preprint2021arXiv

Measuring the Mass of the Large Magellanic Cloud with Stellar Streams Observed by ${S}^5$

Stellar streams are excellent probes of the underlying gravitational potential in which they evolve. In this work, we fit dynamical models to five streams in the Southern Galactic hemisphere, combining observations from the Southern Stellar Stream Spectroscopic Survey (${S}^5$), Gaia EDR3, and the Dark Energy Survey (DES), to measure the mass of the Large Magellanic Cloud (LMC). With an ensemble of streams, we find a mass of the LMC ranging from 14 to $19 \times 10^{10}\ \mathrm{M}_{\odot}$, probed over a range of closest approach times and distances. With the most constraining stream (Orphan-Chenab), we measure an LMC mass of $18.8^{+ 3.5}_{- 4.0} \times 10^{10}\ \mathrm{M}_{\odot}$, probed at a closest approach time of 310 Myr and a closest approach distance of 25.4 kpc. This mass is compatible with previous measurements, showing that a consistent picture is emerging of the LMC&#39;s influence on structures in the Milky Way. Using this sample of streams, we find that the LMC&#39;s effect depends on the relative orientation of the stream and LMC at their point of closest approach. To better understand this, we present a simple model based on the impulse approximation and we show that the LMC&#39;s effect depends both on the magnitude of the velocity kick imparted to the stream and the direction of this kick.

preprint2021arXiv

The dynamics of the globular cluster NGC3201 out to the Jacobi radius

As part of a chemo-dynamical survey of five nearby globular clusters with 2dF/AAOmega on the AAT, we have obtained kinematic information for the globular cluster NGC3201. Our new observations confirm the presence of a significant velocity gradient across the cluster which can almost entirely be explained by the high proper motion of the cluster. After subtracting the contribution of this perspective rotation, we found a remaining rotation signal with an amplitude of $\sim1\ km/s$ around a different axis to what we expect from the tidal tails and the potential escapers, suggesting that this rotation is internal and can be a remnant of its formation process. At the outer part, we found a rotational signal that is likely a result from potential escapers. The proper motion dispersion at large radii reported by Bianchini et al. has previously been attributed to dark matter. Here we show that the LOS dispersion between 0.5-1 Jacobi radius is lower, yet above the predictions from an N-body model of NGC3201 that we ran for this study. Based on the simulation, we find that potential escapers cannot fully explain the observed velocity dispersion. We also estimate the effect on the velocity dispersion of different amounts of stellar-mass black holes and unbound stars from the tidal tails with varying escape rates and find that these effects cannot explain the difference between the LOS dispersion and the N-body model. Given the recent discovery of tidal tail stars at large distances from the cluster, a dark matter halo is an unlikely explanation. We show that the effect of binary stars, which is not included in the N-body model, is important and can explain part of the difference in dispersion. We speculate that the remaining difference must be the result of effects not included in the N-body model, such as initial cluster rotation, velocity anisotropy and Galactic substructure.

preprint2020arXiv

Exploring the nature and synchronicity of early cluster formation in the Large Magellanic Cloud V: Multiple Populations in ancient Globular Clusters

We examine four ancient Large Magellanic Cloud (LMC) globular clusters (GCs) for evidence of multiple stellar populations using the Advanced Camera for Surveys and Wide Field Camera 3 on the Hubble Space Telescope Programme GO-14164. NGC 1466, NGC 1841, and NGC 2257 all show evidence for a redder, secondary population along the main-sequence. Reticulum does not show evidence for the presence of a redder population, but this GC has the least number of stars and Monte Carlo simulations indicate that the sample of main sequence stars is too small to robustly infer whether a redder population exists in this cluster. The second, redder, population of the other three clusters constitutes $\sim30-40\%$ of the total population along the main-sequence. This brings the total number of ancient LMC GCs with known split or broadened main-sequences to five. However, unlike for Hodge 11 and NGC 2210 (see arXiv:1904.01434), none of the clusters show evidence for multiple populations in the horizontal branch. We also do not find evidence of a second population along the Red Giant Branch (RGB).

preprint2020arXiv

How stellar rotation shapes the colour magnitude diagram of the massive intermediate-age star cluster NGC 1846

We present a detailed study of stellar rotation in the massive 1.5 Gyr old cluster NGC 1846 in the Large Magellanic Cloud. Similar to other clusters at this age, NGC 1846 shows an extended main sequence turn-off (eMSTO), and previous photometric studies have suggested it could be bimodal. In this study, we use MUSE integral-field spectroscopy to measure the projected rotational velocities (vsini) of around 1400 stars across the eMSTO and along the upper main sequence of NGC 1846. We measure vsini values up to ~250 km/s and find a clear relation between the vsini of a star and its location across the eMSTO. Closer inspection of the distribution of rotation rates reveals evidence for a bimodal distribution, with the fast rotators centred around vsini = 140 km/s and the slow rotators centred around vsini = 60 km/s. We further observe a lack of fast rotating stars along the photometric binary sequence of NGC 1846, confirming results from the field that suggest that tidal interactions in binary systems can spin down stars. However, we do not detect a significant difference in the binary fractions of the fast and slowly rotating sub-populations. Finally, we report on the serendipitous discovery of a planetary nebula associated with NGC 1846.

preprint2020arXiv

Milky Way Tomography with the SkyMapper Southern Survey. II. Photometric Re-calibration of SMSS DR2

We apply the spectroscopy-based stellar-color regression (SCR) method to perform an accurate photometric re-calibration of the second data release from the SkyMapper Southern Survey (SMSS DR2). From comparison with a sample of over 200,000 dwarf stars with stellar atmospheric parameters taken from GALAH+ DR3 and with accurate, homogeneous photometry from $Gaia$ DR2, zero-point offsets are detected in the original photometric catalog of SMSS DR2, in particular for the gravity- and metallicity-sensitive $uv$ bands. For $uv$ bands, the zero-point offsets are close to zero at very low extinction, and then steadily increase with $E (B - V)$, reaching as large as 0.174 and 0.134 mag respectively, at $E (B - V) \sim 0.5$ mag. These offsets largely arise from the adopted dust term in the transformations used by SMSS DR2 to construct photometric calibrators from the ATLAS reference catalog. For the $gr$ bands, the zero-point offsets exhibit negligible variations with SFD $E(B - V )$, due to their tiny coefficients on the dust term in the transformation. Our study also reveals small, but significant, spatial variations of the zero-point offsets in all $uvgr$ bands. External checks using Strömgren photometry, WD loci and the SDSS Stripe 82 standard-star catalog independently confirm the zero-points found by our revised SCR method.

preprint2020arXiv

SkyMapper Southern Survey: Second Data Release (DR2)

We present the second data release (DR2) of the SkyMapper Southern Survey, a hemispheric survey carried out with the SkyMapper Telescope at Siding Spring Observatory in Australia, using six optical filters: $u,v,g,r,i,z$. DR2 is the first release to go beyond the $\sim$18mag (10$σ$) limit of the Shallow Survey released in DR1, and includes portions of the sky at full survey depth that reach >21mag in $g$ and $r$ filters. The DR2 photometry has a precision as measured by internal reproducibility of 1% in $u$ and $v$, and 0.7% in $griz$. More than 21 000 deg$^2$ have data in some filters (at either Shallow or Main Survey depth) and over 7 000 deg$^2$ have deep Main Survey coverage in all six filters. Finally, about 18 000 deg$^2$ have Main Survey data in $i$ and $z$ filters, albeit not yet at full depth. The release contains over 120 000 images, as well as catalogues with over 500 million unique astrophysical objects and nearly 5 billion individual detections. It also contains cross-matches with a range of external catalogues such as Gaia DR2, Pan-STARRS1 DR1, GALEX GUVcat, 2MASS, and AllWISE, as well as spectroscopic surveys such as 2MRS, GALAH, 6dFGS, and 2dFLenS.

preprint2020arXiv

The tidal remnant of an unusually metal-poor globular cluster

Globular clusters are some of the oldest bound stellar structures observed in the Universe. They are ubiquitous in large galaxies and are believed to trace intense star formation events and the hierarchical build-up of structure. Observations of globular clusters in the Milky Way, and a wide variety of other galaxies, have found evidence for a `metallicity floor&#39;, whereby no globular clusters are found with chemical (`metal&#39;) abundances below approximately 0.3 to 0.4 per cent of that of the Sun. The existence of this metallicity floor may reflect a minimum mass and a maximum redshift for surviving globular clusters to form, both critical components for understanding the build-up of mass in the universe. Here we report measurements from the Southern Stellar Streams Spectroscopic Survey of the spatially thin, dynamically cold Phoenix stellar stream in the halo of the Milky Way. The properties of the Phoenix stream are consistent with it being the tidally disrupted remains of a globular cluster. However, its metal abundance ([Fe/H] = -2.7) is substantially below that of the empirical metallicity floor. The Phoenix stream thus represents the debris of the most metal-poor globular cluster discovered so far, and its progenitor is distinct from the present-day globular cluster population in the local Universe. Its existence implies that globular clusters below the metallicity floor have probably existed, but were destroyed during Galactic evolution.

preprint2019arXiv

A SkyMapper view of the Large Magellanic Cloud: The dynamics of stellar populations

We present the first SkyMapper stellar population analysis of the Large Magellanic Cloud (hereafter LMC),including the identification of 3578 candidate Carbon Stars through their extremely red $g-r$ colours. Coupled with Gaia astrometry, we analyse the distribution and kinematics of this Carbon Star population, finding the LMC to be centred at $(R.A., Dec.) = (80.90^{\circ}\pm{0.29}, -68.74^{\circ}\pm{0.12})$, with a bulk proper motion of $(μ_α,μ_δ) = (1.878\pm0.007,0.293\pm0.018) \mathrm{mas\ yr^{-1}}$ and a disk inclination of $i = 25.6^{\circ}\pm1.1$ at position angle $θ= 135.6^{\circ}\pm 3.3^{\circ}$. We complement this study with the identification and analysis of additional stellar populations, finding that the dynamical centre for Red Giant Branch (RGB) stars is similar to that seen for the Carbon Stars, whereas for young stars the dynamical centre is significantly offset from the older populations. This potentially indicates that the young stars were formed as a consequence of a strong tidal interaction, probably with the Small Magellanic Cloud (SMC). In terms of internal dynamics, the tangential velocity profile increases linearly within $\sim3\ \mathrm{kpc}$, after which it maintains an approximately constant value of $V_{rot} = 83.6\pm 1.7 \mathrm{km\ s^{-1}}$ until $\sim7 \mathrm{kpc}$. With an asymmetric drift correction, we estimate the mass within $7 \mathrm{kpc}$ to be $M_{\rm LMC}(<7\mathrm{kpc}) = (2.5\pm0.1)\times10^{10}\ \mathrm{M}_{\odot}$ and within the tidal radius ($\sim 30\ \mathrm{kpc}$) to be $M_{\rm LMC}(<30\mathrm{kpc}) = (1.06 \pm 0.32)\times10^{11}\ \mathrm{M}_{\odot}$, consistent with other recent measurements.