Researcher profile

Geraint F. Lewis

Geraint F. Lewis contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
37works
0followers
10topics
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

37 published item(s)

preprint2026arXiv

Multiple Populations and a CH Star Found in the 300S Globular Cluster Stellar Stream

Milky Way globular clusters (GCs) display chemical enrichment in a phenomenon called multiple stellar populations (MSPs). While the enrichment mechanism is not fully understood, there is a correlation between a cluster's mass and the fraction of enriched stars found therein. However, present-day GC masses are often smaller than their masses at the time of formation due to dynamical mass loss. In this work, we explore the relationship between mass and MSPs using the stellar stream 300S. We present the chemical abundances of eight red giant branch member stars in 300S with high-resolution spectroscopy from Magellan/MIKE. We identify one enriched star characteristic of MSPs and no detectable metallicity dispersion, confirming that the progenitor of 300S was a globular cluster. The fraction of enriched stars (12.5\%) observed in our 300S stars is less than the 50\% of stars found enriched in Milky Way GCs of comparable present-day mass ($\sim10^{4.5}$\msun). We calculate the mass of 300S's progenitor and compare it to the initial masses of intact GCs, finding that 300S aligns well with the trend between the system mass at formation and enrichment. 300S's progenitor may straddle the critical mass threshold for the formation of MSPs and can therefore serve as a benchmark for the stellar enrichment process. Additionally, we identify a CH star, with high abundances of \textit{s}-process elements, probably accreted from a binary companion. The rarity of such binaries in intact GCs may imply stellar streams permit the survival of binaries that would otherwise be disrupted.

preprint2026arXiv

Revisiting The Cosmological Time Dilation of Distant Quasars: Influence of Source Properties and Evolution

After decades of searching, cosmological time dilation was recently identified in the timescale of variability seen in distant quasars. Here, we expand on the previous analysis to disentangle this cosmological signal from the influence of the properties of the source population, specifically the quasar bolometric luminosity and the rest-frame emission wavelength at which the variability was observed. Furthermore, we consider the potential influence of the evolution of the quasar population over cosmic time. We find that a significant intrinsic scatter of 0.288 +- 0.021 dex in the variability timescales, which was not considered in the previous analysis, is favoured by the data. This slightly increases the uncertainty in the results. However, the expected cosmological dependence of the variability timescales is confirmed to be robust to changes in the underlying assumptions. We find that the variability timescales increase smoothly with both wavelength and bolometric luminosity, and that black hole mass has no effect on the variability timescale once rest wavelength and bolometric luminosity are accounted for. Moreover, if the standard cosmological model is correct, governed by relativistic expansion, we also find very little cosmological evolution in the intrinsic variability timescales of distant quasars.

preprint2026arXiv

Rotation of the Globular Cluster Population of the Dark Matter Deficient Galaxy NGC 1052-DF4: Implication for the Total Mass

We explore the globular cluster population of NGC 1052-DF4, a dark matter deficient galaxy, using Bayesian inference to search for the presence of rotation. The existence of such a rotating component is relevant to the estimation of the mass of the galaxy, and therefore the question of whether NGC 1052-DF4 is truly deficient of dark matter, similar to NGC 1052-DF2 another galaxy in the same group. The rotational characteristics of seven globular clusters in NGC 1052-DF4 were investigated, finding that a non-rotating kinematic model has a higher Bayesian evidence than a rotating model, by a factor of approximately 2.5. In addition, we find that under the assumption of rotation, its amplitude must be small. This distinct lack of rotation strengthens the case that, based on its intrinsic velocity dispersion, NGC 1052-DF4 is a truly dark matter deficient galaxy.

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

LensingETC: a tool to optimize multi-filter imaging campaigns of galaxy-scale strong lensing systems

Imaging data is the principal observable required to use galaxy-scale strong lensing in a multitude of applications in extragalactic astrophysics and cosmology. In this paper, we develop Lensing Exposure Time Calculator (LensingETC) to optimize the efficiency of telescope time usage when planning multi-filter imaging campaigns for galaxy-scale strong lenses. This tool simulates realistic data tailored to specified instrument characteristics and then automatically models them to assess the power of the data in constraining lens model parameters. We demonstrate a use case of this tool by optimizing a two-filter observing strategy (in IR and UVIS) within the limited exposure time per system allowed by a Hubble Space Telescope (HST) Snapshot program. We find that higher resolution is more advantageous to gain constraining power on the lensing observables, when there is a trade-off between signal-to-noise ratio and resolution; e.g., between the UVIS and IR filters of the HST. We also find that, whereas a point spread function (PSF) with sub-Nyquist sampling allows the sample mean for a model parameter to be robustly recovered for both galaxy-galaxy and point-source lensing systems, a sub-Nyquist sampled PSF introduces a larger scatter than a Nyquist sampled one in the deviation from the ground truth for point-source lens systems.

preprint2022arXiv

Metallicity distribution of the progenitor of the Giant Stellar Stream in the Andromeda Galaxy

The Giant Stellar Stream (GSS) is a prominent tidal feature in the halo of the Andromeda Galaxy (M31), representing the ongoing destruction of a satellite galaxy. In this paper, we investigate the formation of the GSS through detailed numerical simulations of the tidal disruption of a progenitor system. Assuming that the stream was created in a single merger event between M31 and a dwarf spheroidal galaxy with stellar mass of $10^{9}M_{\odot}$, we successfully reproduce the dynamical properties of the GSS. As the metallicity distribution along the stream has been well determined from the observations (PAndAS and AMIGA data sets), we use Monte Carlo simulations to reconstruct the original metallicity distribution of the dwarf progenitor. We find that a progenitor dwarf galaxy with a negative radial metallicity gradient, $Δ$FeH = -$0.3 \pm 0.2$, reproduces the observed GSS properties at a time between 2.4 and 2.9 Gyrs into the merger. We also show that the observed double peak metallicity distribution along the stream is a transitory structure caused by unique merger circumstances where two groups of streaming stars are moving in opposite directions, intersecting to produce the peaks.

preprint2022arXiv

Observational window effects on multi-object Reverberation Mapping

Contemporary reverberation mapping campaigns are employing wide-area photometric data and high-multiplex spectroscopy to efficiently monitor hundreds of active galactic nuclei (AGN). However, the interaction of the window function(s) imposed by the observation cadence with the reverberation lag and AGN variability time scales (intrinsic to each source over a range of luminosities) impact our ability to recover these fundamental physical properties. Time dilation effects due to the sample source redshift distribution introduces added complexity. We present comprehensive analysis of the implications of observational cadence, seasonal gaps and campaign baseline duration (i.e., the survey window function) for reverberation lag recovery. We find the presence of a significant seasonal gap dominates the efficacy of any given campaign strategy for lag recovery across the parameter space, particularly for those sources with observed-frame lags above 100 days. Using the OzDES survey as a baseline, we consider the implications of this analysis for the 4MOST/TiDES campaign providing concurrent follow-up of the LSST deep-drilling fields, as well as upcoming programs. We conclude that the success of such surveys will be critically limited by the seasonal visibility of some potential field choices, but show significant improvement from extending the baseline. Optimising the sample selection to fit the window function will improve survey efficacy.

preprint2022arXiv

Ringing the universe with cosmic emptiness: void properties through a combined analysis of stacked weak gravitational and Doppler lensing

An essential aspect of cosmic voids is that these underdense regions provide complementary information about the properties of our Universe. Unlike dense regions, voids are avoided by matter and are less contaminated by baryonic processes. The first step to understanding the properties of cosmic voids is to correctly infer their mass profiles. In the literature, various techniques have been implemented. In this paper, we review them and implement a new technique that is based on Doppler lensing. We use a relativistic $N$-body code, \textsc{Gevolution}, to generate cosmological mass perturbations and implement a three-dimensional ray-tracing technique, which follows the evolution of a ray-bundles. We focus on the various properties of cosmic voids (e.g. void size function, 2-point correlation function, and the density profile of voids), and compare the results with their universal trends. We show that when weak-lensing is combined with the Doppler lensing we obtain even tighter constraints than weak-lensing alone. We also obtain better agreement between density profiles within central parts of voids inferred from lensing and density profiles inferred from halo tracers. The implication of the result relevant to the ongoing and prospective low-redshift spectroscopic surveys is briefly discussed.

preprint2022arXiv

The GALAH Survey: A New Sample of Extremely Metal-Poor Stars Using A Machine Learning Classification Algorithm

Extremely Metal-Poor (EMP) stars provide a valuable probe of early chemical enrichment in the Milky Way. Here we leverage a large sample of $\sim600,000$ high-resolution stellar spectra from the GALAH survey plus a machine learning algorithm to find 54 candidates with estimated [Fe/H]~$\leq$~-3.0, 6 of which have [Fe/H]~$\leq$~-3.5. Our sample includes $\sim 20 \%$ main sequence EMP candidates, unusually high for \emp surveys. We find the magnitude-limited metallicity distribution function of our sample is consistent with previous work that used more complex selection criteria. The method we present has significant potential for application to the next generation of massive stellar spectroscopic surveys, which will expand the available spectroscopic data well into the millions of stars.

preprint2022arXiv

The GALAH Survey: Chemical tagging and chrono-chemodynamics of accreted halo stars with GALAH+ DR3 and $Gaia$ eDR3

Since the advent of $Gaia$ astrometry, it is possible to identify massive accreted systems within the Galaxy through their unique dynamical signatures. One such system, $Gaia$-Sausage-Enceladus (GSE), appears to be an early &#34;building block&#34; given its virial mass $> 10^{10}\,\mathrm{M_\odot}$ at infall ($z\sim1-3$). In order to separate the progenitor population from the background stars, we investigate its chemical properties with up to 30 element abundances from the GALAH+ Survey Data Release 3 (DR3). To inform our choice of elements for purely chemically selecting accreted stars, we analyse 4164 stars with low-$α$ abundances and halo kinematics. These are most different to the Milky Way stars for abundances of Mg, Si, Na, Al, Mn, Fe, Ni, and Cu. Based on the significance of abundance differences and detection rates, we apply Gaussian mixture models to various element abundance combinations. We find the most populated and least contaminated component, which we confirm to represent GSE, contains 1049 stars selected via [Na/Fe] vs. [Mg/Mn] in GALAH+ DR3. We provide tables of our selections and report the chrono-chemodynamical properties (age, chemistry, and dynamics). Through a previously reported clean dynamical selection of GSE stars, including $30 < \sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} < 55$, we can characterise an unprecedented 24 abundances of this structure with GALAH+ DR3. Our chemical selection allows us to prevent circular reasoning and characterise the dynamical properties of the GSE, for example mean $\sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} = 26_{-14}^{+9}$. We find only $(29\pm1)\%$ of the GSE stars within the clean dynamical selection region. Our methodology will improve future studies of accreted structures and their importance for the formation of the Milky Way.

preprint2022arXiv

The Hierarchical Structure of Galactic Haloes: Generalised N-Dimensional Clustering with CluSTAR-ND

We present CluSTAR-ND, a fast hierarchical galaxy/(sub)halo finder that produces {\bf Clu}stering {\bf S}tructure via {\bf T}ransformative {\bf A}ggregation and {\bf R}ejection in {\bf N}-{\bf D}imensions. It is designed to improve upon Halo-OPTICS -- an algorithm that automatically detects and extracts significant astrophysical clusters from the 3D spatial positions of simulation particles -- by decreasing run-times, possessing the capability for metric adaptivity, and being readily applicable to data with any number of features. We directly compare these algorithms and find that not only does CluSTAR-ND produce a similarly robust clustering structure, it does so in a run-time that is at least $3$ orders of magnitude faster. In optimising CluSTAR-ND&#39;s clustering performance, we have also carefully calibrated $4$ of the $7$ CluSTAR-ND parameters which -- unless specified by the user -- will be automatically and optimally chosen based on the input data. We conclude that CluSTAR-ND is a robust astrophysical clustering algorithm that can be leveraged to find stellar satellite groups on large synthetic or observational data sets.

preprint2022arXiv

The K2 Galactic Archaeology Program Data Release 3: Age-abundance patterns in C1-C8 and C10-C18

We present the third and final data release of the K2 Galactic Archaeology Program (K2 GAP) for Campaigns C1-C8 and C10-C18. We provide asteroseismic radius and mass coefficients, $κ_R$ and $κ_M$, for $\sim 19,000$ red giant stars, which translate directly to radius and mass given a temperature. As such, K2 GAP DR3 represents the largest asteroseismic sample in the literature to date. K2 GAP DR3 stellar parameters are calibrated to be on an absolute parallactic scale based on Gaia DR2, with red giant branch and red clump evolutionary state classifications provided via a machine-learning approach. Combining these stellar parameters with GALAH DR3 spectroscopy, we determine asteroseismic ages with precisions of $\sim 20-30\%$ and compare age-abundance relations to Galactic chemical evolution models among both low- and high-$α$ populations for $α$, light, iron-peak, and neutron-capture elements. We confirm recent indications in the literature of both increased Ba production at late Galactic times, as well as significant contribution to r-process enrichment from prompt sources associated with, e.g., core-collapse supernovae. With an eye toward other Galactic archaeology applications, we characterize K2 GAP DR3 uncertainties and completeness using injection tests, suggesting K2 GAP DR3 is largely unbiased in mass/age and with uncertainties of $2.9\%\,(\rm{stat.})\,\pm0.1\%\,(\rm{syst.})$ & $6.7\%\,(\rm{stat.})\,\pm0.3\%\,(\rm{syst.})$ in $κ_R$ & $κ_M$ for red giant branch stars and $4.7\%\,(\rm{stat.})\,\pm0.3\%\,(\rm{syst.})$ & $11\%\,(\rm{stat.})\,\pm0.9\%\,(\rm{syst.})$ for red clump stars. We also identify percent-level asteroseismic systematics, which are likely related to the time baseline of the underlying data, and which therefore should be considered in TESS asteroseismic analysis.

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 GALAH survey: tracing the Galactic disk with Open Clusters

Open clusters are unique tracers of the history of our own Galaxy&#39;s disk. According to our membership analysis based on \textit{Gaia} astrometry, out of the 226 potential clusters falling in the footprint of GALAH or APOGEE, we find that 205 have secure members that were observed by at least one of the survey. Furthermore, members of 134 clusters have high-quality spectroscopic data that we use to determine their chemical composition. We leverage this information to study the chemical distribution throughout the Galactic disk of 21 elements, from C to Eu. The radial metallicity gradient obtained from our analysis is $-$0.076$\pm$0.009 dex kpc$^{-1}$, which is in agreement with previous works based on smaller samples. Furthermore, the gradient in the [Fe/H] - guiding radius (r$_{\rm guid}$) plane is $-$0.073$\pm$0.008 dex kpc$^{-1}$. We show consistently that open clusters trace the distribution of chemical elements throughout the Galactic disk differently than field stars. In particular, at given radius, open clusters show an age-metallicity relation that has less scatter than field stars. As such scatter is often interpreted as an effect of radial migration, we suggest that these differences are due to the physical selection effect imposed by our Galaxy: clusters that would have migrated significantly also had higher chances to get destroyed. Finally, our results reveal trends in the [X/Fe]$-$r$_{\rm guid}$$-$age space, which are important to understand production rates of different elements as a function of space and time.

preprint2020arXiv

A unified framework for 21cm tomography sample generation and parameter inference with Progressively Growing GANs

Creating a database of 21cm brightness temperature signals from the Epoch of Reionisation (EoR) for an array of reionisation histories is a complex and computationally expensive task, given the range of astrophysical processes involved and the possibly high-dimensional parameter space that is to be probed. We utilise a specific type of neural network, a Progressively Growing Generative Adversarial Network (PGGAN), to produce realistic tomography images of the 21cm brightness temperature during the EoR, covering a continuous three-dimensional parameter space that models varying X-ray emissivity, Lyman band emissivity, and ratio between hard and soft X-rays. The GPU-trained network generates new samples at a resolution of $\sim 3&#39;$ in a second (on a laptop CPU), and the resulting global 21cm signal, power spectrum, and pixel distribution function agree well with those of the training data, taken from the 21SSD catalogue \citep{Semelin2017}. Finally, we showcase how a trained PGGAN can be leveraged for the converse task of inferring parameters from 21cm tomography samples via Approximate Bayesian Computation.

preprint2020arXiv

Abundances in the Milky Way across five nucleosynthetic channels from 4 million LAMOST stars

Large stellar surveys are revealing the chemodynamical structure of the Galaxy across a vast spatial extent. However, the many millions of low-resolution spectra observed to date are yet to be fully exploited. We employ The Cannon, a data-driven approach to estimating abundances, to obtain detailed abundances from low-resolution (R = 1800) LAMOST spectra, using the GALAH survey as our reference. We deliver five (for dwarfs) or six (for giants) estimated abundances representing five different nucleosynthetic channels, for 3.9 million stars, to a precision of 0.05 - 0.23 dex. Using wide binary pairs, we demonstrate that our abundance estimates provide chemical discriminating power beyond metallicity alone. We show the coverage of our catalogue with radial, azimuthal and dynamical abundance maps, and examine the neutron capture abundances across the disk and halo, which indicate different origins for the in-situ and accreted halo populations. LAMOST has near-complete Gaia coverage and provides an unprecedented perspective on chemistry across the Milky Way.

preprint2020arXiv

Gravitational Microlensing Time Delays at High Optical Depth: Image Parities and the Temporal Properties of Fast Radio Bursts

Due to differing gravitational potentials and path lengths, gravitational lensing induces time delays between multiple images of a source which, for solar mass objects, is of order $\sim10^{-5}$ seconds. If an astrophysically compact source, such as a Fast Radio Burst (FRB), is observed through a region with a high optical depth of such microlensing masses, this gravitational lensing time delay can be imprinted on short timescale transient signals. In this paper, we consider the impact of the parity of the macroimage on the resultant microlensing time delays. It is found that this parity is directly imprinted on the microlensing signal, with macroimages formed at minima of the time arrival surface beginning with the most highly magnified microimages and then progressing to the fainter microimages. At macroimages at the maxima of the time arrival surface, this situation is reversed, with fainter images observed first and finishing with the brightest microimages. For macroimages at saddle-points, the signal again begins with fainter images, followed by brighter images before again fading through the fainter microimages. The growing populations of cosmologically distant bursty transient sources will undoubtedly result in the discovery of strong lensed, multiply imaged FRBs, which will be susceptible to microlensing by compact masses. With the temporal resolution being offered my modern and future facilities, the detection of microlensing induced time delays will reveal the parities of the gravitational lens macroimages, providing additional constraints on macrolensing mass models and improving the efficacy of these transient sources as a cosmological probes.

preprint2020arXiv

Mapping the tilt of the Milky Way bulge velocity ellipsoids with ARGOS and $Gaia$ DR2

Until the recent advent of $Gaia$ Data Release 2 (DR2) and deep multi-object spectroscopy, it has been difficult to obtain 6-D phase space information for large numbers of stars beyond 4 kpc, in particular towards the Galactic centre, where dust and crowding effects are significant. In this study we combine line-of-sight velocities from the Abundances and Radial velocity Galactic Origins Survey (ARGOS) spectroscopic survey with proper motions from $Gaia$ DR2, to obtain a sample of $\sim$ 7,000 red clump stars with 3-D velocities. We perform a large scale stellar kinematics study of the Milky Way (MW) bulge to characterize the bulge velocity ellipsoids. We measure the tilt $l_{v}$ of the major-axis of the velocity ellipsoid in the radial-longitudinal velocity plane in 20 fields across the bulge. The tilt or vertex deviation, is characteristic of non-axisymmetric systems and a significant tilt is a robust indicator of non-axisymmetry or bar presence. We compare the observations to the predicted kinematics of an N-body boxy-bulge model formed from dynamical instabilities. In the model, the $l_{v}$ values are strongly correlated with the angle ($α$) between the bulge major-axis and the Sun-Galactic centre line-of-sight. We use a maximum likelihood method to obtain an independent measurement of $α$, from bulge stellar kinematics alone. The most likely value of $α$ given our model is $α= (29 \pm 3)^{\circ}$. In the Baade&#39;s window, the metal-rich stars display a larger vertex deviation ($l_{v} = -40^{\circ}$) than the metal-poor stars ($l_{v} = 10^{\circ}$) but we do not detect significant $l_{v}-$metallicity trends in the other fields.

preprint2020arXiv

The GALAH Survey: A new constraint on cosmological lithium and Galactic lithium evolution from warm dwarf stars

Lithium depletion and enrichment in the cosmos is not yet well understood. To help tighten constraints on stellar and Galactic evolution models, we present the largest high-resolution analysis of Li abundances A(Li) to date, with results for over 100 000 GALAH field stars spanning effective temperatures $5900\,\mathrm{K} \lesssim \rm{T_{eff}} \lesssim7000\,\mathrm{K}$ and metallicities $-3 \lesssim \rm[Fe/H] \lesssim +0.5$. We separated these stars into two groups, on the warm and cool side of the so-called Li-dip, a localised region of the Kiel diagram wherein lithium is severely depleted. We discovered that stars in these two groups show similar trends in the A(Li)-[Fe/H] plane, but with a roughly constant offset in A(Li) of 0.4 dex, the warm group having higher Li abundances. At $\rm[Fe/H]\gtrsim-0.5$, a significant increasing in Li abundance with increasing metallicity is evident in both groups, signalling the onset of significant Galactic production. At lower metallicity, stars in the cool group sit on the Spite plateau, showing a reduced lithium of around 0.4 dex relative to the primordial value predicted from Big Bang nucleosynthesis (BBN). However, stars in the warm group between [Fe/H] = -1.0 and -0.5, form an elevated plateau that is largely consistent with the BBN prediction. This may indicate that these stars in fact preserve the primordial Li produced in the early Universe.

preprint2020arXiv

The GCE in a New Light: Disentangling the $γ$-ray Sky with Bayesian Graph Convolutional Neural Networks

A fundamental question regarding the Galactic Center Excess (GCE) is whether the underlying structure is point-like or smooth. This debate, often framed in terms of a millisecond pulsar or annihilating dark matter (DM) origin for the emission, awaits a conclusive resolution. In this work we weigh in on the problem using Bayesian graph convolutional neural networks. In simulated data, our neural network (NN) is able to reconstruct the flux of inner Galaxy emission components to on average $\sim$0.5%, comparable to the non-Poissonian template fit (NPTF). When applied to the actual $\textit{Fermi}$-LAT data, we find that the NN estimates for the flux fractions from the background templates are consistent with the NPTF; however, the GCE is almost entirely attributed to smooth emission. While suggestive, we do not claim a definitive resolution for the GCE, as the NN tends to underestimate the flux of point-sources peaked near the 1$σ$ detection threshold. Yet the technique displays robustness to a number of systematics, including reconstructing injected DM, diffuse mismodeling, and unmodeled north-south asymmetries. So while the NN is hinting at a smooth origin for the GCE at present, with further refinements we argue that Bayesian Deep Learning is well placed to resolve this DM mystery.

preprint2020arXiv

The Hierarchical Structure of Galactic Haloes: Classification and characterisation with Halo-OPTICS

We build upon Ordering Points To Identify Clustering Structure (OPTICS), a hierarchical clustering algorithm well-known to be a robust data-miner, in order to produce Halo-OPTICS, an algorithm designed for the automatic detection and extraction of all meaningful clusters between any two arbitrary sizes. We then apply Halo-OPTICS to the 3D spatial positions of halo particles within four separate synthetic Milky Way type galaxies, classifying the stellar and dark matter structural hierarchies. Through visualisation of the Halo-OPTICS output, we compare its structure identification to the state-of-the-art galaxy/(sub)halo finder VELOCIraptor, finding excellent agreement even though Halo-OPTICS does not consider kinematic information in this current implementation. We conclude that Halo-OPTICS is a robust hierarchical halo finder, although its determination of lower spatial-density features such as the tails of streams could be improved with the inclusion of extra localised information such as particle kinematics and stellar metallicity into its distance metric.

preprint2020arXiv

The One-Way Speed of Light and the Milne Universe

In Einstein&#39;s Special Theory of Relativity, all observers measure the speed of light, $c$, to be the same. However, this refers to the round trip speed, where a clock at the origin times the outward and return trip of light reflecting off a distant mirror. Measuring the one-way speed of light is fraught with issues of clock synchronisation, and, as long as the average speed of light remains $c$, the speeds on the outward and return legs could be different. One objection to this anisotropic speed of light is that views of the distant universe would be different in different directions, especially with regards to the ages of observed objects and the smoothness of the Cosmic Microwave Background. In this paper, we explore this in the Milne universe, the limiting case of a Friedmann-Robertson-Walker universe containing no matter, radiation or dark energy. Given that this universe is empty, it can be mapped onto flat Minkowski space-time, and so can be explored in terms of the one-way speed of light. The conclusion is that the presence of an anisotropic speed of light leads to anisotropic time dilation effects, and hence observers in the Milne universe would be presented with an isotropic view of the distant cosmos.

preprint2020arXiv

The Pristine Inner Galaxy Survey (PIGS) II: Uncovering the most metal-poor populations in the inner Milky Way

Metal-poor stars are important tools for tracing the early history of the Milky Way, and for learning about the first generations of stars. Simulations suggest that the oldest metal-poor stars are to be found in the inner Galaxy. Typical bulge surveys, however, lack low metallicity ([Fe/H] < -1.0) stars because the inner Galaxy is predominantly metal-rich. The aim of the Pristine Inner Galaxy Survey (PIGS) is to study the metal-poor and very metal-poor (VMP, [Fe/H] < -2.0) stars in this region. In PIGS, metal-poor targets for spectroscopic follow-up are selected from metallicity-sensitive CaHK photometry from the CFHT. This work presents the ~250 deg^2 photometric survey as well as intermediate-resolution spectroscopic follow-up observations for ~8000 stars using AAOmega on the AAT. The spectra are analysed using two independent tools: ULySS with an empirical spectral library, and FERRE with a library of synthetic spectra. The comparison between the two methods enables a robust determination of the stellar parameters and their uncertainties. We present a sample of 1300 VMP stars -- the largest sample of VMP stars in the inner Galaxy to date. Additionally, our spectroscopic dataset includes ~1700 horizontal branch stars, which are useful metal-poor standard candles. We furthermore show that PIGS photometry selects VMP stars with unprecedented efficiency: 86%/80% (lower/higher extinction) of the best candidates satisfy [Fe/H] < -2.0, as do 80%/63% of a larger, less strictly selected sample. We discuss future applications of this unique dataset that will further our understanding of the chemical and dynamical evolution of the innermost regions of our Galaxy.

preprint2020arXiv

The SAMI Galaxy Survey: Gas velocity dispersions in low-$z$ star-forming galaxies and the drivers of turbulence

We infer the intrinsic ionised gas kinematics for 383 star-forming galaxies across a range of integrated star-formation rates (SFR $\in [10^{-3}, 10^2]$ M$_\odot$ yr$^{-1}$) at $z \lesssim 0.1$ using a consistent 3D forward-modelling technique. The total sample is a combination of galaxies from the SAMI Galaxy Survey and DYNAMO survey. For typical low-$z$ galaxies taken from the SAMI Galaxy Survey, we find the vertical velocity dispersion ($σ_{v, z}$) to be positively correlated with measures of star-formation rate, stellar mass, HI gas mass, and rotational velocity. The greatest correlation is with star-formation rate surface density ($Σ_\text{SFR}$). Using the total sample, we find $σ_{v, z}$ increases slowly as a function of integrated star-formation rate in the range SFR $\in$ [$10^{-3}$, 1] M$_\odot$ yr$^{-1}$ from $17\pm3$ km s$^{-1}$ to $24\pm5$ km s$^{-1}$ followed by a steeper increase up to $σ_{v, z}$ $\sim 80$ km s$^{-1}$ for SFR $\gtrsim 1$ M$_\odot$ yr$^{-1}$. This is consistent with recent theoretical models that suggest a $σ_{v, z}$ floor driven by star-formation feedback processes with an upturn in $σ_{v, z}$ at higher SFR driven by gravitational transport of gas through the disc.

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 data-driven model of nucleosynthesis with chemical tagging in a lower-dimensional latent space

Chemical tagging seeks to identify unique star formation sites from present-day stellar abundances. Previous techniques have treated each abundance dimension as being statistically independent, despite theoretical expectations that many elements can be produced by more than one nucleosynthetic process. In this work we introduce a data-driven model of nucleosynthesis where a set of latent factors (e.g., nucleosynthetic yields) contribute to all stars with different scores, and clustering (e.g., chemical tagging) is modelled by a mixture of multivariate Gaussians in a lower-dimensional latent space. We use an exact method to simultaneously estimate the factor scores for each star, the partial assignment of each star to each cluster, and the latent factors common to all stars, even in the presence of missing data entries. We use an information-theoretic Bayesian principle to estimate the number of latent factors and clusters. Using the second Galah data release we find that six latent factors are preferred to explain N = 2,566 stars with 17 chemical abundances. We identify the rapid- and slow-neutron capture processes, as well as latent factors consistent with Fe-peak and α-element production, and another where K and Zn dominate. When we consider N ~ 160,000 stars with missing abundances we find another 7 factors, as well as 16 components in latent space. Despite these components showing separation in chemistry that is explained through different yield contributions, none show significant structure in their positions or motions. We argue that more data, and joint priors on cluster membership that are constrained by dynamical models, are necessary to realise chemical tagging at a galactic-scale. We release software that allows for model parameters to be optimised in seconds given a fixed number of latent factors, components, and $10^7$ abundance measurements.

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.

preprint2019arXiv

Comparing the Quenching Times of Faint M31 and Milky Way Satellite Galaxies

We present the star formation histories (SFHs) of 20 faint M31 satellites ($-12 \lesssim M_V \lesssim -6$) that were measured by modeling sub-horizontal branch (HB) depth color-magnitude diagrams constructed from Hubble Space Telescope (HST) imaging. Reinforcing previous results, we find that virtually all galaxies quenched between 3 and 9 Gyr ago, independent of luminosity, with a notable concentration $3-6$ Gyr ago. This is in contrast to the Milky Way (MW) satellites, which are generally either faint with ancient quenching times or luminous with recent ($<3$ Gyr) quenching times. We suggest that systematic differences in the quenching times of M31 and MW satellites may be a reflection of the varying accretion histories of M31 and the MW. This result implies that the formation histories of low-mass satellites may not be broadly representative of low-mass galaxies in general. Among the M31 satellite population we identify two distinct groups based on their SFHs: one with exponentially declining SFHs ($τ\sim 2$ Gyr) and one with rising SFHs with abrupt quenching. We speculate how these two groups could be related to scenarios for a recent major merger involving M31. The Cycle 27 HST Treasury survey of M31 satellites will provide well-constrained ancient SFHs to go along with the quenching times we measure here. The discovery and characterization of M31 satellites with $M_V \gtrsim -6$ would help quantify the relative contributions of reionization and environment to quenching of the lowest-mass satellites.

preprint2019arXiv

Cosmological Signatures of Dark Sector Physics: The Evolution of Haloes and Spin Alignment

The standard cosmological paradigm currently lacks a detailed account of physics in the dark sector, the dark matter and energy that dominate cosmic evolution. In this paper, we consider the distinguishing factors between three alternative models - warm dark matter, quintessence and coupled dark matter-energy - and \lcdm{} through numerical simulations of cosmological structure formation. Key halo statistics - halo spin/velocity alignment between large-scale structure and neighboring haloes, halo formation time and migration - were compared across cosmologies within the redshift range 0$\leq$$z$$\leq$2.98. We found the alignment of halo motion and spin to large-scale structures and neighbouring haloes to be similar in all cosmologies for a range of redshifts. The search was extended to low density regions, avoiding non-linear disturbances of halo spins, yet very similar alignment trends were found between cosmologies which are difficult to characterize and use as a probe of cosmology. We found haloes in quintessence cosmologies form earlier than their \lcdm{} counterparts. Relating this to the fact that such haloes originate in high density regions, such findings could hold clues to distinguishing factors for the quintessence cosmology from the standard model. Although in general, halo statistics are not an accurate probe of the dark sector physics.

preprint2019arXiv

Discovery of a 21 Myr old stellar population in the Orion complex

The Orion complex is arguably the most studied star-forming region in the Galaxy. While stars are still being born in the Orion nebula, the oldest part was believed to be no more than 13 Myr old. In order to study the full hierarchy of star formation across the Orion complex, we perform a clustering analysis of the Ori OB1a region using new stellar surveys and derive robust ages for each identified stellar aggregate. We use Gaia DR2 parameters supplemented with radial velocities from the GALAH and APOGEE surveys to perform clustering of the Ori OB1a association. Five overdensities are resolved in a six-dimensional parameter space (positions, distance, proper motions, and radial velocity). Most correspond to previously known structures (ASCC 16, 25 Orionis, ASCC 20, ASCC 21). We use Gaia DR2, Pan-STARRS1 and 2MASS photometry to fit isochrones to the colour-magnitude diagrams of the identified clusters. The ages of the clusters can thus be measured with ~10% precision. While four of the clusters have ages between 11 and 13 Myr, the ASCC 20 cluster stands out at an age of 21 $\pm$ 3 Myr. This is significantly greater than the age of any previously known component of the Orion complex. To some degree, all clusters overlap in at least one of the six phase-space dimensions. We argue that the formation history of the Orion complex, and its relation to the Gould belt, must be reconsidered. A significant challenge in reconstructing the history of the Ori OB1a association is to understand the impact of the newly discovered 21 Myr old population on the younger parts of the complex, including their formation.

preprint2019arXiv

Microlensing and Photon Bunching: The impact of decoherence

Gravitational microlensing within the Galaxy offers the prospect of probing the details of distant stellar sources, as well as revealing the distribution of compact (and potentially non-luminous) masses along the line-of-sight. Recently, it has been suggested that additional constraints on the lensing properties can be determined through the measurement of the time delay between images through the correlation of the bunching of photon arrival times; an application of the Hanbury-Brown Twiss effect. In this paper, we revisit this analysis, examining the impact of decoherence of the radiation from a spatially extended source along the multiple paths to an observer. The result is that, for physically reasonable situations, such decoherence completely erases any correlation that could otherwise be used to measure the gravitational lensing time delay. Indeed, the divergent light paths traverse extremely long effective baselines at the lens plane, corresponding to extremes of angular resolving power well beyond those attainable with any terrestrial technologies; the drawback being that few conceivable celestial objects would be sufficiently compact with high enough surface brightness to yield usable signals.

preprint2019arXiv

On the origin of the asymmetric dwarf galaxy distribution around Andromeda

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_{esc}^{2} - V_{init}^{2} \sim 200^2$ -- $300^2\ \rm{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.

preprint2019arXiv

The Canada-France Imaging Survey: Reconstructing the Milky Way Star Formation History from its White Dwarf Population

As the remnants of stars with initial masses $\lesssim$ 8 M$_{\odot}$, white dwarfs contain valuable information on the formation histories of stellar populations. In this paper, we use deep, high-quality, u-band photometry from the Canada France Imaging Survey (CFIS), griz photometry from Pan-STARRS 1 (PS1), as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over $\sim$4500 deg$^2$ using a reduced proper motion diagram. We develop a new white dwarf population synthesis code that returns mock observations of the Galactic field white dwarf population for a given star formation history, while simultaneously taking into account the geometry of the Milky Way, survey parameters, and selection effects. We use this model to derive the star formation histories of the thin disk, thick disk, and stellar halo. Our results show that the Milky Way disk began forming stars (11.3 $\pm$ 0.5) Gyr ago, with a peak rate of (8.8 $\pm$ 1.4) M$_{\odot}$yr$^{-1}$ at (9.8 $\pm$ 0.4) Gyr, before a slow decline to a constant rate until the present day --- consistent with recent results suggesting a merging event with a satellite galaxy. Studying the residuals between the data and best-fit model shows evidence for a slight increase in star formation over the past 3 Gyr. We fit the local fraction of helium-atmosphere white dwarfs to be (21 $\pm$ 3) %. Incorporating this methodology with data from future wide-field surveys such as LSST, Euclid, CASTOR, and WFIRST should provide an unprecedented view into the formation of the Milky Way at its earliest epoch through its white dwarfs.

preprint2019arXiv

The GALAH Survey: Temporal Chemical Enrichment of the Galactic Disk

We present isochrone ages and initial bulk metallicities ($\rm [Fe/H]_{bulk}$, by accounting for diffusion) of 163,722 stars from the GALAH Data Release 2, mainly composed of main sequence turn-off stars and subgiants ($\rm 7000 K>T_{eff}>4000 K$ and $\rm log g>3$ dex). The local age-metallicity relationship (AMR) is nearly flat but with significant scatter at all ages; the scatter is even higher when considering the observed surface abundances. After correcting for selection effects, the AMR appear to have intrinsic structures indicative of two star formation events, which we speculate are connected to the thin and thick disks in the solar neighborhood. We also present abundance ratio trends for 16 elements as a function of age, across different $\rm [Fe/H]_{bulk}$ bins. In general, we find the trends in terms of [X/Fe] vs age from our far larger sample to be compatible with studies based on small ($\sim$ 100 stars) samples of solar twins but we now extend it to both sub- and super-solar metallicities. The $α$-elements show differing behaviour: the hydrostatic $α$-elements O and Mg show a steady decline with time for all metallicities while the explosive $α$-elements Si, Ca and Ti are nearly constant during the thin disk epoch (ages $\lessapprox $ 12 Gyr). The s-process elements Y and Ba show increasing [X/Fe] with time while the r-process element Eu have the opposite trend, thus favouring a primary production from sources with a short time-delay such as core-collapse supernovae over long-delay events such as neutron star mergers.

preprint2019arXiv

The Globular Cluster Population of NGC 1052-DF2: Evidence for Rotation

Based upon the kinematics of ten globular clusters, it has recently been claimed that the ultra-diffuse galaxy, NCD 1052-DF2, lacks a significant quantity of dark matter. Dynamical analyses have generally assumed that this galaxy is pressure supported, with the relatively small velocity dispersion of the globular cluster population indicating the deficit of dark matter. However, the presence of a significant rotation of the globular cluster population could substantially modify this conclusion. Here we present the discovery of such a signature of rotation in the kinematics of NGC 1052-DF2&#39;s globular clusters, with a velocity amplitude of $\sim12.44^{+4.40}_{-5.16}$ km/s, which, through Bayesian model comparison, represents a marginally better fit to the available kinematic data; note that this rotation is distinct from, and approximately perpendicular to, the recently identified rotation of the stellar component of NGC 1052-DF2. Assuming this truly represents an underlying rotation, it is shown that the determined mass depends upon the inclination of the rotational component and, with a moderate inclination, the resultant mass to light ratio can exceed $M/L\sim10$.