Source author record

Geraint F. Lewis

Geraint F. Lewis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

125works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

125 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 "building block" 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'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'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's influence on structures in the Milky Way. Using this sample of streams, we find that the LMC'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'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'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'$ 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'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'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', whereby no globular clusters are found with chemical (`metal') 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'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$.

preprint2016arXiv

A computational approach to the twin paradox in curved spacetime

Despite being a major component in the teaching of special relativity, the twin `paradox' is generally not examined in courses on general relativity. Due to the complexity of analytical solutions to the problem, the paradox is often neglected entirely, and students are left with an incomplete understanding of the relativistic behaviour of time. This article outlines a project, undertaken by undergraduate physics students at the University of Sydney, in which a novel computational method was derived in order to predict the time experienced by a twin following a number of paths between two given spacetime coordinates. By utilising this method, it is possible to make clear to students that following a geodesic in curved spacetime does not always result in the greatest experienced proper time.

preprint2016arXiv

Exploring the reality of density substructures in the Palomar 5 stellar stream

We present an analysis of the presence of substructures in the stellar stream of the Palomar 5 globular cluster, as derived from Sloan Digital Sky Survey data. Using a matched filter technique, we recover the positions and sizes of overdensities reported in previous studies. To explore the reality of these structures, we also create an artificial model of the stream, in which we construct a realistic background on top of which we add a perfectly smooth stream structure, taking into account the effects of photometric completeness and interstellar extinction. We find that the smooth artificial stream then shows similarly-pronounced substructures as the real structure. Interestingly, our best-fit N-body simulation does display real projected density variations linked to stellar epicyclic motions, but these become less significant when taking into account the SDSS star-count constraints. The substructures found when applying our matched filter technique to the N-body particles converted into observable stars are thus mostly unrelated to these epicyclic motions. This analysis suggests that the majority of the previously-detected substructures along the tidal tail of Palomar 5 are artefacts of observational inhomogeneities.

preprint2016arXiv

Major Substructure in the M31 Outer Halo: Distances and Metallicities along the Giant Stellar Stream

We present a renewed look at M31's Giant Stellar Stream along with the nearby structures Stream C and Stream D, exploiting a new algorithm capable of fitting to the red giant branch (RGB) of a structure in both colour and magnitude space. Using this algorithm, we are able to generate probability distributions in distance, metallicity and RGB width for a series of subfields spanning these structures. Specifically, we confirm a distance gradient of approximately 20 kpc per degree along a 6 degree extension of the Giant Stellar Stream, with the farthest subfields from M31 lying ~ 120 kpc more distant than the inner-most subfields. Further, we find a metallicity that steadily increases from -0.7^{+0.1}_{-0.1} dex to -0.2^{+0.2}_{-0.1} dex along the inner half of the stream before steadily dropping to a value of -1.0^{+0.2}_{-0.2} dex at the farthest reaches of our coverage. The RGB width is found to increase rapidly from 0.4^{+0.1}_{-0.1} dex to 1.1^{+0.2}_{-0.1} dex in the inner portion of the stream before plateauing and decreasing marginally in the outer subfields of the stream. In addition, we estimate Stream C to lie at a distance between 794 and 862 kpc and Stream D between 758 kpc and 868 kpc. We estimate the median metallicity of Stream C to lie in the range -0.7 to -1.6 dex and a metallicity of -1.1^{+0.3}_{-0.2} dex for Stream D. RGB widths for the two structures are estimated to lie in the range 0.4 to 1.2 dex and 0.3 to 0.7 dex respectively. In total, measurements are obtained for 19 subfields along the Giant Stellar Stream, 4 along Stream C, 5 along Stream D and 3 general M31 spheroid fields for comparison. We thus provide a higher resolution coverage of the structures in these parameters than has previously been available in the literature.

preprint2016arXiv

MARZ: Manual and Automatic Redshifting Software

The Australian Dark Energy Survey (OzDES) is a 100-night spectroscopic survey underway on the Anglo-Australian Telescope using the fibre-fed 2-degree-field (2dF) spectrograph. We have developed a new redshifting application Marz with greater usability, flexibility, and the capacity to analyse a wider range of object types than the Runz software package previously used for redshifting spectra from 2dF. Marz is an open-source, client-based, Javascript web-application which provides an intuitive interface and powerful automatic matching capabilities on spectra generated from the AAOmega spectrograph to produce high quality spectroscopic redshift measurements. The software can be run interactively or via the command line, and is easily adaptable to other instruments and pipelines if conforming to the current FITS file standard is not possible. Behind the scenes, a modified version of the Autoz cross-correlation algorithm is used to match input spectra against a variety of stellar and galaxy templates, and automatic matching performance for OzDES spectra has increased from 54% (Runz) to 91% (Marz). Spectra not matched correctly by the automatic algorithm can be easily redshifted manually by cycling automatic results, manual template comparison, or marking spectral features.

preprint2016arXiv

Matter in the beam: Weak lensing, substructures and the temperature of dark matter

Warm Dark Matter (WDM) models offer an attractive alternative to the current Cold Dark Matter (CDM) cosmological model. We present a novel method to differentiate between WDM and CDM cosmologies, namely using weak lensing; this provides a unique probe as it is sensitive to all the "matter in the beam", not just dark matter haloes and the galaxies that reside in them, but also the diffuse material between haloes. We compare the weak lensing maps of CDM clusters to those in a WDM model corresponding to a thermally produced $0.5$~keV dark matter particle. Our analysis clearly shows that the weak lensing magnification, convergence and shear distributions can be used to distinguish between CDM and WDM models. WDM models {\em increase} the probability of weak magnifications, with the differences being significant to $\gtrsim5σ$, while leaving no significant imprint on the shear distribution. WDM clusters analysed in this work are more homogeneous than CDM ones, and the fractional decrease in the amount of material in haloes is proportional to the average increase in the magnification. This difference arises from matter that would be bound in compact haloes in CDM being smoothly distributed over much larger volumes at lower densities in WDM. Moreover, the signature does not solely lie in the probability distribution function but in the full spatial distribution of the convergence field.

preprint2016arXiv

On the dynamical state of galaxy clusters: insights from cosmological simulations II

Using a suite of cosmology simulations of a sample of $> 120$ galaxy clusters with $\log(M_{DM, vir}) \le 14.5$. We compare clusters that form in purely dark matter run and their counterparts in hydro runs and investigate 4 independent parameters, that are normally used to classify dynamical state. We find that the virial ratio $η$ in hydro-dynamical runs is $\sim 10$ per cent lower than in the DM run, and there is no clear separation between the relaxed and unrelaxed clusters for any parameter. Further, using the velocity dispersion deviation parameter $ζ$, which is defined as the ratio between cluster velocity dispersion $σ$ and the theoretical prediction $σ_t = \sqrt{G M_{total}/R}$, we find that there is a linear correlation between the virial ratio $η$ and this $ζ$ parameter. We propose to use this $ζ$ parameter, which can be easily derived from observed galaxy clusters, as a substitute of the $η$ parameter to quantify the cluster dynamical state.

preprint2016arXiv

On The Relativity of Redshifts: Does Space Really "Expand"?

In classes on cosmology, students are often told that photons stretch as space expands, but just how physical is this picture? Does space really expand? In this article, we explore the notion of the redshift of light with Einstein's general theory of relativity, showing that the core underpinning principles reveal that redshifts are both simpler and more complex than you might naively think. This has significant implications for the observed redshifting of photons as they travel across the universe, often refereed to as the cosmological redshift, and for the idea of expanding space.

preprint2016arXiv

On the Stability of Satellite Planes I: Effects of Mass, Velocity, Halo Shape and Alignment

The recently discovered vast thin plane of dwarf satellites orbiting the Andromeda Galaxy (M31) adds to the mystery of the small scale distribution of the Local Group's galaxy population. Such well defined planar structures are apparently rare occurrences in cold dark matter cosmological simulations, and we lack a coherent explanation of their formation and existence. In this paper, we explore the long-term survivability of thin planes of dwarfs in galactic halos, focusing, in particular, on systems mimicking the observed Andromeda distribution. The key results show that, in general, planes of dwarf galaxies are fragile, sensitive to the shape of the dark matter halo and other perturbing effects. In fact, long lived planes of satellites only exist in polar orbits in spherical dark matter halos, presenting a challenge to the observed Andromeda plane which is significantly tilted with respect to the optical disk. Our conclusion is that, in standard cosmological models, planes of satellites are generally short lived, and hence we must be located at a relatively special time in the evolution of the Andromeda Plane, lucky enough to see its coherent pattern.

preprint2016arXiv

Primordial Nucleosynthesis in the Rh = ct cosmology: Pouring cold water on the Simmering Universe

Primordial nucleosynthesis is rightly hailed as one of the great successes of the standard cosmological model. Here we consider the initial forging of elements in the recently proposed Rh = ct universe, a cosmology that demands linear evolution of the scale factor. Such a universe cools extremely slowly compared to standard cosmologies, considerably depleting the available neutrons during nucleosynthesis; this has significant implications for the resultant primordial abundances of elements, predicting a minuscule quantity of helium which is profoundly at odds with observations. The production of helium can be enhanced in such a "simmering universe" by boosting the baryon to photon ratio, although more than an order of magnitude increase is required to bring the helium mass fraction into accordance with observations. However, in this scenario, the prolonged period of nucleosynthesis results of the efficient cooking of lighter into heavier elements, impacting the resultant abundances of all elements so that, other than hydrogen and helium, there are virtually no light elements present in the universe. Without the addition of substantial new physics in the early universe, it is difficult to see how the Rh = ct universe can be considered a viable cosmological model.

preprint2016arXiv

Spurious Small-Scale Structure & Discreteness-Driven Relaxation in Cosmological Simulations

There is strong evidence that cosmological N-body simulations dominated by Warm Dark Matter (WDM) contain spurious or unphysical haloes, most readily apparent as regularly spaced low-mass haloes strung along filaments. We show that spurious haloes are a feature of traditional N-body simulations of cosmological structure formation models, including WDM and Cold Dark Matter (CDM) models, in which gravitational collapse proceeds in an initially anisotropic fashion, and arises naturally as a consequence of discreteness-driven relaxation. We demonstrate this using controlled N-body simulations of plane-symmetric collapse and show that spurious haloes are seeded at shell crossing by localised velocity perturbations induced by the discrete nature of the density field, and that their characteristic separation should be approximately the mean inter-particle separation of the N-body simulation, which is fixed by the mass resolution within the volume. Using cosmological N-body simulations in which particles are split into two collisionless components of fixed mass ratio, we find that the spatial distribution of the two components show signatures of discreteness-driven relaxation in their spatial distribution on both large and small scales. Adopting a spline kernel gravitational softening that is of order the comoving mean inter-particle separation helps to suppress the effect of discreteness-driven relaxation, but cannot eliminate it completely. These results provide further motivation for recent developments of new algorithms, which include, for example, revisions of the traditional N-body approach by means of spatially adaptive anistropric gravitational softenings or explicit solutions for the evolution of dark matter in phase space.

preprint2016arXiv

The GALAH survey: The data reduction pipeline

We present the data reduction procedures being used by the GALAH survey, carried out with the HERMES fibre-fed, multi-object spectrograph on the 3.9~m Anglo-Australian Telescope. GALAH is a unique survey, targeting 1 million stars brighter than magnitude V=14 at a resolution of 28,000 with a goal to measure the abundances of 29 elements. Such a large number of high resolution spectra necessitates the development of a reduction pipeline optimized for speed, accuracy, and consistency. We outline the design and structure of the Iraf-based reduction pipeline that we developed, specifically for GALAH, to produce fully calibrated spectra aimed for subsequent stellar atmospheric parameter estimation. The pipeline takes advantage of existing Iraf routines and other readily available software so as to be simple to maintain, testable and reliable. A radial velocity and stellar atmospheric parameter estimator code is also presented, which is used for further data analysis and yields a useful verification of the reduction quality. We have used this estimator to quantify the data quality of GALAH for fibre cross-talk level ($\lesssim0.5$%) and scattered light ($\sim5$ counts in a typical 20 minutes exposure), resolution across the field, sky spectrum properties, wavelength solution reliability (better than $1$ $\mathrm{km\ s^{-1}}$ accuracy) and radial velocity precision.

preprint2016arXiv

The PAndAS view of the Andromeda satellite system - II. Detailed properties of 23 M31 dwarf spheroidal galaxies

We present a comprehensive analysis of the structural properties and luminosities of the 23 dwarf spheroidal galaxies that fall within the footprint of the Pan-Andromeda Archaeological Survey (PAndAS). These dwarf galaxies represent the large majority of Andromeda's known satellite dwarf galaxies and cover a wide range in luminosity ($-11.6<M_V<-5.8$ or $10^{4.2}< L <10^{6.5} L_\odot$) and surface brightness ($25.1<μ_0<29.3$ mag/arcsec$^2$). We confirm most previous measurements, but find And XIX to be significantly larger than before ($r_h=3065^{+1065}_{-935}$ pc, $M_V=-10.1^{+0.8}_{-0.4}$) and cannot derive parameters for And XXVII as it is likely not a bound stellar system. We also significantly revise downward the luminosities of And~XV and And~XVI, which are now $M_V\sim-7.5$ or $L\sim10^5 L_\odot$. Finally, we provide the first detailed analysis of Cas II/And XXX, a fairly faint system ($M_V=-8.0^{+0.4}_{-0.3}$) of typical size ($r_h=270\pm50$ pc), located in close proximity to the two bright elliptical dwarf galaxies NGC 147 & 185. Combined with the set of homogeneous distances published in an earlier contribution, our analysis dutifully tracks all relevant sources of uncertainty in the determination of the properties of the dwarf galaxies from the PAndAS photometric catalogue. We further publish the posterior probability distribution functions of all the parameters we fit for in the form of MCMC chains available online; these inputs should be used in any analysis that aims to remain truthful to the data and properly account for covariance between parameters.

preprint2016arXiv

Ultracompact minihalos as probes of inflationary cosmology

Cosmological inflation generates primordial density perturbations on all scales, including those far too small to contribute to the cosmic microwave background. At these scales, isolated ultracompact minihalos of dark matter can form well before standard structure formation, if the perturbations have sufficient amplitude. Minihalos affect pulsar timing data and are potentially bright sources of gamma rays. The resulting constraints significantly extend the observable window of inflation in the presence of cold dark matter, coupling two of the key problems in modern cosmology.

preprint2015arXiv

Eppur si muove: Positional and kinematic correlations of satellite pairs in the low Z universe

We have recently shown (Ibata et al. 2014) that pairs of satellite galaxies located diametrically opposite each other around their host possess predominantly anti-correlated velocities. This is consistent with a scenario in which $\sim 50$% of satellite galaxies belong to kinematically-coherent rotating planar structures, similar to those detected around the giant galaxies of the Local Group. Here we extend this analysis, examining the incidence of satellites of giant galaxies drawn from an SDSS photometric redshift catalog. We find that there is a $\sim 17$% overabundance ($> 3 σ$ significance) of candidate satellites at positions diametrically opposite a spectroscopically confirmed satellite. We show that cosmological simulations do not possess this property when the contamination is included, and that there are in fact, after subtracting contamination, 2 to 3 times more satellites diametrically opposed to a spectroscopically confirmed satellite than at $90°$ from it. We also examine the correlation between the satellite pair positions and the orientation of the host galaxy major axis. We find that those satellite pairs with anti-correlated velocities have a strong preference ($\sim 3:1$) to align with the major axis of the host whereas those with correlated velocities display the opposite behavior. This correlation of the satellite alignments appears to be stronger than the well-documented preference of satellites to be located close to the major axis of their host. We finally show that repeating a similar analysis to Ibata et al. (2014) with same-side satellites is generally hard to interpret, but is not inconsistent with our previous results when strong quality-cuts are applied on the sample. All these unexpected correlations strongly suggest that a substantial fraction of satellite galaxies are causally-linked in their formation and evolution.

preprint2015arXiv

Feeling the pull, a study of natural Galactic accelerometers. I: photometry of the delicate stellar stream of the Palomar 5 globular cluster

We present an analysis of wide-field photometric surveys of the Palomar 5 globular cluster and its stellar stream, based on g- and r-band measures together with narrow-band DDO51 photometry. In this first study, we use the deep (g,r) data to measure the incidence of gaps and peaks along the stream. Examining the star-counts profile of the stream plus contaminating populations, we find no evidence for significant under-densities, and find only a single significant over-density. This is at odds with earlier studies based on matched-filter maps derived from shallower SDSS data if the contaminating population possesses plausible spatial properties. The lack of substantial sub-structure along the stream may be used in future dynamical simulations to examine the incidence of dark matter sub-halos in the Galactic halo. We also present a measurement of the relative distances along the stream which we use to create the deepest wide-field map of this system to date.

preprint2015arXiv

Hidden from view: Coupled Dark Sector Physics and Small Scales

We study cluster mass dark matter haloes, their progenitors and surroundings in an coupled Dark Matter-Dark Energy model and compare it to quintessence and $Λ$CDM models with adiabatic zoom simulations. When comparing cosmologies with different expansions histories, growth functions & power spectra, care must be taken to identify unambiguous signatures of alternative cosmologies. Shared cosmological parameters, such as $σ_8$, need not be the same for optimal fits to observational data. We choose to set our parameters to $Λ$CDM $z=0$ values. We find that in coupled models, where DM decays into DE, haloes appear remarkably similar to $Λ$CDM haloes despite DM experiencing an additional frictional force. Density profiles are not systematically different and the subhalo populations have similar mass, spin, and spatial distributions, although (sub)haloes are less concentrated on average in coupled cosmologies. However, given the scatter in related observables ($V_{\rm max},R_{V_{\rm max}}$), this difference is unlikely to distinguish between coupled and uncoupled DM. Observations of satellites of MW and M31 indicate a significant subpopulation reside in a plane. Coupled models do produce planar arrangements of satellites of higher statistical significance than $Λ$CDM models, however, in all models these planes are dynamically unstable. In general, the nonlinear dynamics within and near large haloes masks the effects of a coupled dark sector. The sole environmental signature we find is that small haloes residing in the outskirts are more deficient in baryons than their $Λ$CDM counterparts. The lack of a pronounced signal for a coupled dark sector strongly suggests that such a phenomena would be effectively hidden from view.

preprint2015arXiv

How does our choice of observable influence our estimation of the centre of a galaxy cluster? Insights from cosmological simulations

Galaxy clusters are an established and powerful test-bed for theories of both galaxy evolution and cosmology. Accurate interpretation of cluster observations often requires robust identification of the location of the centre. Using a statistical sample of clusters drawn from a suite of cosmological simulations in which we have explored a range of galaxy formation models, we investigate how the location of this centre is affected by the choice of observable - stars, hot gas, or the full mass distribution as can be probed by the gravitational potential. We explore several measures of cluster centre: the minimum of the gravitational potential, which would expect to define the centre if the cluster is in dynamical equilibrium; the peak of the density; the centre of BCG; and the peak and centroid of X-ray luminosity. We find that the centre of BCG correlates more strongly with the minimum of the gravitational potential than the X-ray defined centres, while AGN feedback acts to significantly enhance the offset between the peak X-ray luminosity and minimum gravitational potential. These results highlight the importance of centre identification when interpreting clusters observations, in particular when comparing theoretical predictions and observational data.

preprint2015arXiv

NGC 147, NGC 185 and CassII: a genetic approach to orbital properties, star formation and tidal debris

NGC147, NGC185 and CassiopeiaII (CassII) have similar positions in the sky, distances and measured line of sight velocities. This proximity in phase space suggests that these three satellites of M31 form a subgroup within the Local Group. Nevertheless, the differences in their star formation history and interstellar medium, and the recent discovery of a stellar stream in NGC~147, combined with the lack of tidal features in the other two satellites, are all indications of complex and diverse interactions between M31 and these three satellites. We use a genetic algorithm to explore the different orbits that these satellites can have and select six sets of orbits that could best explain the observational features of the NGC147, NGC185 and CassII satellites. The parameters of these orbits are then used as a starting point for N-body simulations. We present models for which NGC147, NGC185 and CassII are a bound group for a total time of at least one Gyr but still undergo different interactions with M31 and as a result NGC147 has a clear stellar stream whereas the other two satellites have no significant tidal features. This result shows that it is possible to find solutions that reproduce the contrasting properties of the satellites and for which NGC147-NGC185-CassII have been gravitationally bound.

preprint2015arXiv

The Spatially-Resolved Star Formation History of the M31 Outer Disc

We present deep Hubble Space Telescope Advanced Camera for Surveys observations of the stellar populations in two fields lying at 20 and 23 kpc from the centre of M31 along the south-west semi-major axis. These data enable the construction of colour-magnitude diagrams reaching the oldest main-sequence turn-offs (~13 Gyr) which, when combined with another field at 25 kpc from our previous work, we use to derive the first precision constraints on the spatially-resolved star formation history of the M31 disc. The star formation rates exhibit temporal as well as field-to-field variations, but are generally always within a factor of two of their time average. There is no evidence of inside-out growth over the radial range probed. We find a median age of ~7.5 Gyr, indicating that roughly half of the stellar mass in the M31 outer disc was formed before z ~ 1. We also find that the age-metallicity relations (AMRs) are smoothly increasing from [Fe/H]~-0.4 to solar metallicity between 10 and 3 Gyr ago, contrary to the flat AMR of the Milky Way disc at a similar number of scale lengths. Our findings provide insight on the roles of stellar feedback and radial migration in the formation and evolution of large disc galaxies.

preprint2015arXiv

Trans-Dimensional Bayesian Inference for Gravitational Lens Substructures

We introduce a Bayesian solution to the problem of inferring the density profile of strong gravitational lenses when the lens galaxy may contain multiple dark or faint substructures. The source and lens models are based on a superposition of an unknown number of non-negative basis functions (or "blobs") whose form was chosen with speed as a primary criterion. The prior distribution for the blobs' properties is specified hierarchically, so the mass function of substructures is a natural output of the method. We use reversible jump Markov Chain Monte Carlo (MCMC) within Diffusive Nested Sampling (DNS) to sample the posterior distribution and evaluate the marginal likelihood of the model, including the summation over the unknown number of blobs in the source and the lens. We demonstrate the method on two simulated data sets: one with a single substructure, and one with ten. We also apply the method to the g-band image of the "Cosmic Horseshoe" system, and find evidence for more than zero substructures. However, these have large spatial extent and probably only point to misspecifications in the model (such as the shape of the smooth lens component or the point spread function), which are difficult to guard against in full generality.

preprint2014arXiv

A Genetic Approach to the History of the Magellanic Clouds

The history of the Magellanic Clouds is investigated using N-body hydrodynamic simulations where the initial conditions are set by a genetic algorithm. This technique allows us to identify possible orbits for the Magellanic Clouds around the Milky Way, by directly comparing the simulations with observational constraints. We explore the parameter space of the interaction between the Magellanic Clouds and the Milky Way, considering as free parameters the proper motions of the Magellanic Clouds, the virial mass and the concentration parameter ($c$) of the Galactic dark matter halo. The best orbital scenarios presented here are considered with two different sets of parameters for the Milky Way disc and bulge components. The total circular velocity at the Sun's position ($\rm{R_{\odot}}=8.5\,\rm {kpc}$) is directly calculated from the rotation curve of the corresponding Galactic mass model. Our analysis suggests that the Magellanic Clouds have orbited inside the virial radius of the Milky Way for at least $3\,\rm{Gyr}$, even for low mass haloes. However, this is possible only with high values for the concentration parameter ($c\ge20$). In both orbital models presented here, the mutual interaction between the Magellanic Clouds is able to reproduce the observed features of the Magellanic System.

preprint2014arXiv

A thousand shadows of Andromeda: rotating planes of satellites in the Millennium-II cosmological simulation

In a recent contribution, Bahl \& Baumgardt investigated the incidence of planar alignments of satellite galaxies in the Millennium-II simulation, and concluded that vast thin planes of dwarf galaxies, similar to that observed in the Andromeda galaxy (M31), occur frequently by chance in $Λ$-Cold Dark Matter cosmology. However, their analysis did not capture the essential fact that the observed alignment is simultaneously radially extended, yet thin, and kinematically unusual. With the caveat that the Millennium-II simulation may not have sufficient mass resolution to identify confidently simulacra of low-luminosity dwarf galaxies, we re-examine that simulation for planar structures, using the same method as employed by Ibata et al. (2013) on the real M31 satellites. We find that 0.04\% of host galaxies display satellite alignments that are at least as extreme as the observations, when we consider their extent, thickness and number of members rotating in the same sense. We further investigate the angular momentum properties of the co-planar satellites, and find that the median of the specific angular momentum derived from the line of sight velocities in the real M31 structure ($1.3\times10^4$ km/s kpc) is very high compared to systems drawn from the simulations. This analysis confirms that it is highly unlikely that the observed structure around the Andromeda galaxy is due to a chance occurrence. Interestingly, the few extreme systems that are similar to M31 arise from the accretion of a massive sub-halo with its own spatially-concentrated entourage of orphan satellites.

preprint2014arXiv

Another coincidence problem for $Λ$CDM?

Over the last nine years of cosmic microwave background observations, the Wilkinson Microwave Anisotropy Probe ($WMAP$) results were consistent with a $Λ$CDM cosmological model in which the age of the Universe is one Hubble time, and the time-averaged value of the deceleration parameter is consistent with zero. This curious observation has been put forward as a new coincidence problem for the $Λ$CDM concordance cosmology, which is in fact a `greater' coincidence than the near equality of the density parameters of matter and the cosmological constant. At the moment of writing these conference proceedings, the Planck Collaboration has released its first cosmological data, which revealed a small shift in the $Λ$CDM cosmological parameters when compared to $WMAP$. We show that under the assumption of a spatially flat $Λ$CDM cosmology, Planck's results remove this coincidence problem for $Λ$CDM at greater than 99\% confidence level.

preprint2014arXiv

Comparing the observable properties of dwarf galaxies on and off the Andromeda plane

The thin, extended planes of satellite galaxies detected around both the Milky Way and Andromeda are not a natural prediction of the LCDM paradigm. Galaxies in these distinct planes may have formed and evolved in a different way (e.g., tidally) to their off-plane neighbours. If this were the case, one would expect the on- and off-plane dwarf galaxies in Andromeda to have experienced different evolutionary histories, which should be reflected by the chemistries, dynamics, and star formation histories of the two populations. In this work, we present new, robust kinematic observations for 2 on-plane M31 dSphs (And XVI and XVII) and compile and compare all available observational metrics for the on- and off-plane dwarfs to search for a signal that would corroborate such a hypothesis. We find that, barring their spatial alignment, the on- and off-plane Andromeda dwarf galaxies are indistinguishable from one another, arguing against vastly different formative and evolutionary histories for these two populations.

preprint2014arXiv

Hydrodynamical simulations of coupled and uncoupled quintessence models I: Halo properties and the cosmic web

We present the results of a series of adiabatic hydrodynamical simulations of several quintessence models (both with a free and an interacting scalar field) in comparison to a standard \LCDM\ cosmology. For each we use $2\times1024^3$ particles in a $250$\hMpc\ periodic box assuming WMAP7 cosmology. In this work we focus on the properties of haloes in the cosmic web at $z=0$. The web is classified into \emph{voids}, \emph{sheets}, \emph{filaments} and \emph{knots} depending on the eigenvalues of the velocity shear tensor, which are an excellent proxy for the underlying overdensity distribution. We find that the properties of objects classified according to their surrounding environment shows a substantial dependence on the underlying cosmology; for example, while $V_{\rm max}$ shows average deviations of $\approx5$ per cent across the different models when considering the full halo sample, comparing objects classified according to their environment, the size of the deviation can be as large as $20$ per cent. We also find that halo spin parameters are positively correlated to the coupling, whereas halo concentrations show the opposite behaviour. Furthermore, when studying the concentration-mass relation in different environments, we find that in all cosmologies underdense regions have a larger normalization and a shallower slope. While this behaviour is found to characterize all the models, differences in the best-fit relations are enhanced in (coupled) dark energy, thus providing a clearer prediction for this class of models.

preprint2014arXiv

Hydrodynamical simulations of coupled and uncoupled quintessence models II: Galaxy clusters

We study the $z=0$ properties of clusters (and large groups) of galaxies within the context of interacting and non-interacting quintessence cosmological models, using a series of adiabatic SPH simulations. Initially, we examine the average properties of groups and clusters, quantifying their differences in LCDM, uncoupled Dark Energy (\ude) and coupled Dark Energy (\cde) cosmologies. In particular, we focus upon radial profiles of the gas density, temperature and pressure, and we also investigate how the standard hydrodynamic equilibrium hypothesis holds in quintessence cosmologies. While we are able to confirm previous results about the distribution of baryons, we also find that the main discrepancy (with differences up to $20\%$) can be seen in cluster pressure profiles. We then switch attention to individual structures, mapping each halo in quintessence cosmology to its \LCDM\ counterpart. We are able to identify a series of small correlations between the coupling in the dark sector and halo spin, triaxiality and virialization ratio. When looking at spin and virialization of dark matter haloes, we find a weak ($5\%$) but systematic deviation in fifth force scenarios from \LCDM.

preprint2014arXiv

On the Shoulders of Giants: Properties of the Stellar Halo and the Milky Way Mass Distribution

Halo stars orbit within the potential of the Milky Way and hence their kinematics can be used to understand the underlying mass distribution. However, the inferred mass distribution depends sensitively upon assumptions made on the density and the velocity anisotropy profiles of the tracers. Also, there is a degeneracy between the parameters of the halo and that of the disk or bulge. Here, we decompose the Galaxy into bulge, disk and dark matter halo and then model the kinematic data of the halo BHB and K-giants from the SEGUE. Additionally, we use the gas terminal velocity curve and the Sgr A$^*$ proper motion. With $R_\odot = 8.5$kpc, our study reveals that the density of the stellar halo has a break at $17.2^{+1.1}_{-1.0}$ kpc, and an exponential cut-off in the outer parts starting at $97.7^{+15.6}_{-15.8}$kpc. Also, we find the velocity anisotropy is radially biased with $β_s= 0.4\pm{0.2}$ in the outer halo. We measure halo virial mass $M_{\text{vir}} = 0.80^{+0.31}_{-0.16} \times 10^{12} M_{\odot}$, concentration $c=21.1^{+14.8}_{-8.3}$, disk mass of $0.95^{+0.24}_{-0.30}\times10^{11} M_{\odot}$, disk scale length of $4.9^{+0.4}_{-0.4}$ kpc and bulge mass of $0.91^{+0.31}_{-0.38} \times 10^{10} M_{\odot}$. The mass of halo is found to be small and this has important consequences. The giant stars reveal that the outermost halo stars have low velocity dispersion interestingly suggesting a truncation of the stellar halo density rather than a small overall mass of the Galaxy. Our estimates of local escape velocity $v_{\rm esc} = 550.9^{+32.4}_{-22.1}$ kms$^{-1}$ and dark matter density $ρ^{\rm DM}_{\odot} = 0.0088^{+0.0024}_{-0.0018} M_{\odot} {\rm pc^{-3}} $ ($0.35^{+0.08}_{-0.07}$ GeV cm$^{-3}$) are in good agreement with recent estimates. Some of the above estimates are depended on the adopted value of $R_\odot$ and of outer power-law index of the tracer number density.

preprint2014arXiv

Probing spatial homogeneity with LTB models: a detailed discussion

Do current observational data confirm the assumptions of the cosmological principle, or is there statistical evidence for deviations from spatial homogeneity on large scales? To address these questions, we developed a flexible framework based on spherically symmetric, but radially inhomogeneous Lemaitre-Tolman-Bondi (LTB) models with synchronous Big Bang. We expanded the (local) matter density profile in terms of flexible interpolation schemes and orthonormal polynomials. A Monte Carlo technique in combination with recent observational data was used to systematically vary the shape of these profiles. In the first part of this article, we reconsider giant LTB voids without dark energy to investigate whether extremely fine-tuned mass profiles can reconcile these models with current data. While the local Hubble rate and supernovae can easily be fitted without dark energy, however, model-independent constraints from the Planck 2013 data require an unrealistically low local Hubble rate, which is strongly inconsistent with the observed value; this result agrees well with previous studies. In the second part, we explain why it seems natural to extend our framework by a non-zero cosmological constant, which then allows us to perform general tests of the cosmological principle. Moreover, these extended models facilitate explorating whether fluctuations in the local matter density profile might potentially alleviate the tension between local and global measurements of the Hubble rate, as derived from Cepheid-calibrated type Ia supernovae and CMB experiments, respectively. We show that current data provide no evidence for deviations from spatial homogeneity on large scales. More accurate constraints are required to ultimately confirm the validity of the cosmological principle, however.

preprint2014arXiv

Resolving the mass--anisotropy degeneracy of the spherically symmetric Jeans equation I: theoretical foundation

A widely employed method for estimating the mass of stellar systems with apparent spherical symmetry is dynamical modelling using the spherically symmetric Jeans equation. Unfortunately this approach suffers from a degeneracy between the assumed mass density and the second order velocity moments. This degeneracy can lead to significantly different predictions for the mass content of the system under investigation, and thus poses a barrier for accurate estimates of the dark matter content of astrophysical systems. In a series of papers we describe an algorithm that removes this degeneracy and allows for unbiased mass estimates of systems of constant or variable mass-to-light ratio. The present contribution sets the theoretical foundation of the method that reconstructs a unique kinematic profile for some assumed free functional form of the mass density. The essence of our method lies in using flexible B-spline functions for the representation of the radial velocity dispersion in the spherically symmetric Jeans equation. We demonstrate our algorithm through an application to synthetic data for the case of an isotropic King model with fixed mass-to-light ratio, recovering excellent fits of theoretical functions to observables and a unique solution. The mass-anisotropy degeneracy is removed to the extent that, for an assumed functional form of the potential and mass density pair $(Φ,ρ)$, and a given set of line-of-sight velocity dispersion $σ_{los}^2$ observables, we recover a unique profile for $σ_{rr}^2$ and $σ_{tt}^2$. Our algorithm is simple, easy to apply and provides an efficient means to reconstruct the kinematic profile.

preprint2014arXiv

Resolving the mass--anisotropy degeneracy of the spherically symmetric Jeans equation II: optimum smoothing and model validation

The spherical Jeans equation is widely used to estimate the mass content of a stellar systems with apparent spherical symmetry. However, this method suffers from a degeneracy between the assumed mass density and the kinematic anisotropy profile, $β(r)$. In a previous work, we laid the theoretical foundations for an algorithm that combines smoothing B-splines with equations from dynamics to remove this degeneracy. Specifically, our method reconstructs a unique kinematic profile of $σ_{rr}^2$ and $σ_{tt}^2$ for an assumed free functional form of the potential and mass density $(Φ,ρ)$ and given a set of observed line-of-sight velocity dispersion measurements, $σ_{los}^2$. In Paper I (submitted to MNRAS: MN-14-0101-MJ) we demonstrated the efficiency of our algorithm with a very simple example and we commented on the need for optimum smoothing of the B-spline representation; this is in order to avoid unphysical variational behaviour when we have large uncertainty in our data. In the current contribution we present a process of finding the optimum smoothing for a given data set by using information of the behaviour from known ideal theoretical models. Markov Chain Monte Carlo methods are used to explore the degeneracy in the dynamical modelling process. We validate our model through applications to synthetic data for systems with constant or variable mass-to-light ratio $Υ$. In all cases we recover excellent fits of theoretical functions to observables and unique solutions. Our algorithm is a robust method for the removal of the mass-anisotropy degeneracy of the spherically symmetric Jeans equation for an assumed functional form of the mass density.

preprint2014arXiv

Stochastic Star Formation & Feedback: Mapping Low-Mass Galaxies to Dark Matter Haloes

Comparison of observed satellite galaxies of the Milky Way (hereafter MW) with dark matter subhaloes in cosmological $N$-body simulations of MW-mass haloes suggest that such subhaloes, if they exist, are occupied by satellites in a stochastic fashion. We examine how inefficient massive star formation and associated supernova feedback in high-redshift progenitors of present-day low-mass subhaloes might contribute to this stochasticity. Using a Monte Carlo approach to follow the assembly histories of present-day low-mass haloes with $10^7 \lesssim M \leq 10^{10}$ ${\rm M}_{\odot}$, we identify when cooling and star formation is likely to proceed, and observe that haloes with present-day masses $\lesssim 10^9 {\rm M}_{\odot}$ never grow sufficiently massive to support atomic hydrogen line cooling. Noting that the star formation timescale decreases sharply with stellar mass as $t_{\rm PMS} \propto m_{\ast}^{-2.5}$, we argue that, should the conditions for high mass star formation arise in low-mass haloes, the ensuing supernovae are likely to disrupt ongoing lower-mass star formation and unbind gas within the halo. This potentially star-forming gas is unlikely to be replenished in lower mass haloes because of, e.g. cosmological reionization, and so we expect galaxy formation to be stymied in a manner that depends on host halo assembly history and the efficiency and timing of star formation in proto-galaxies, which we illustrate using a Monte Carlo model. Based on these simple physical arguments, we assert that stochasticity of star formation and feedback is an essential but overlooked ingredient in modelling galaxy formation on the smallest scales.

preprint2014arXiv

The Nature and Origin of Substructure in the Outskirts of M31 -- II. Detailed Star Formation Histories

While wide-field surveys of M31 have revealed much substructure at large radii, understanding the nature and origin of this material is not straightforward from morphology alone. Using deep HST/ACS data, we have derived further constraints in the form of quantitative star formation histories (SFHs) for 14 inner halo fields which sample diverse substructures. In agreement with our previous analysis of colour-magnitude diagram morphologies, we find the resultant behaviours can be broadly separated into two categories. The SFHs of 'disc-like' fields indicate that most of their mass has formed since z~1, with one quarter of the mass formed in the last 5 Gyr. We find 'stream-like' fields to be on average 1.5 Gyr older, with <10 percent of their stellar mass formed within the last 5 Gyr. These fields are also characterised by an age--metallicity relation showing rapid chemical enrichment to solar metallicity by z=1, suggestive of an early-type progenitor. We confirm a significant burst of star formation 2 Gyr ago, discovered in our previous work, in all the fields studied here. The presence of these young stars in our most remote fields suggests that they have not formed in situ but have been kicked-out from through disc heating in the recent past.

preprint2014arXiv

The PAndAS Field of Streams: stellar structures in the Milky Way halo toward Andromeda and Triangulum

We reveal the highly structured nature of the Milky Way stellar halo within the footprint of the PAndAS photometric survey from blue main sequence and main sequence turn-off stars. We map no fewer than five stellar structures within a heliocentric range of ~5 to 30 kpc. Some of these are known (the Monoceros Ring, the Pisces/Triangulum globular cluster stream), but we also uncover three well-defined stellar structures that could be, at least partly, responsible for the so-called Triangulum/Andromeda and Triangulum/Andromeda 2 features. In particular, we trace a new faint stellar stream located at a heliocentric distance of ~17 kpc. With a surface brightness of Σ_V ~ 32-32.5 mag/arcsec^2, it follows an orbit that is almost parallel to the Galactic plane north of M31 and has so far eluded surveys of the Milky Way halo as these tend to steer away from regions dominated by the Galactic disk. Investigating our follow-up spectroscopic observations of PAndAS, we serendipitously uncover a radial velocity signature from stars that have colors and magnitudes compatible with the stream. From the velocity of eight likely member stars, we show that this stellar structure is dynamically cold, with an unresolved velocity dispersion that is lower than 7.1 km/s at the 90-percent confidence level. Along with the width of the stream (300-650 pc), its dynamics points to a dwarf-galaxy-accretion origin. The numerous stellar structures we can map in the Milky Way stellar halo between 5 and 30 kpc and their varying morphology is a testament to the complex nature of the stellar halo at these intermediate distances.

preprint2014arXiv

Velocity anti-correlation of diametrically opposed galaxy satellites in the low redshift universe

Recent work has shown that both the Milky Way and the Andromeda galaxies possess the unexpected property that their dwarf satellite galaxies are aligned in thin and kinematically coherent planar structures. It is now important to evaluate the incidence of such planar structures in the larger galactic population, since the Local Group may not be a sufficiently representative environment. Here we report that the measurement of the velocity of pairs of diametrically opposed galaxy satellites provides a means to determine statistically the prevalence of kinematically coherent planar alignments. In the local universe (redshift $z<0.05$), we find that such satellite pairs out to a galactocentric distance of $150$ kpc are preferentially anti-correlated in their velocities (99.994% confidence level), and that the distribution of galaxies in the larger scale environment (beyond 150 kpc and up to $\approx 2$ Mpc) is strongly elongated along the axis joining the inner satellite pair ($>7σ$ confidence). Our finding may indicate that co-rotating planes of satellites, similar to that seen around the Andromeda galaxy, are ubiquitous in nature, while their coherent motion also suggests that they are a significant repository of angular momentum on $\sim 100$ kpc scales.

preprint2014arXiv

Warm Dark Haloes Accretion Histories and their Gravitational Signatures

We study clusters in Warm Dark Matter (WDM) models of a thermally produced dark matter particle $0.5$ keV in mass. We show that, despite clusters in WDM cosmologies having similar density profiles as their Cold Dark Matter (CDM) counterparts, the internal properties, such as the amount of substructure, shows marked differences. This result is surprising as clusters are at mass scales that are {\em a thousand times greater} than that at which structure formation is suppressed. WDM clusters gain significantly more mass via smooth accretion and contain fewer substructures than their CDM brethren. The higher smooth mass accretion results in subhaloes which are physically more extended and less dense. These fine-scale differences can be probed by strong gravitational lensing. We find, unexpectedly, that WDM clusters have {\em higher} lensing efficiencies than those in CDM cosmologies, contrary to the naive expectation that WDM clusters should be less efficient due to the fewer substructures they contain. Despite being less dense, the larger WDM subhaloes are more likely to have larger lensing cross-sections than CDM ones. Additionally, WDM subhaloes typically reside at larger distances, which radially stretches the critical lines associated with strong gravitational lensing, resulting in excess in the number of clusters with large radial cross-sections at the $\sim2σ$ level. Though lensing profile for an individual cluster vary significantly with the line-of-sight, the radial arc distribution based on a sample of $\gtrsim100$ clusters may prove to be the crucial test for the presence of WDM.

preprint2013arXiv

A kinematic study of the Andromeda dwarf spheroidal system

We present a homogeneous kinematic analysis of red giant branch stars within 18 of the 28 Andromeda dwarf spheroidal (dSph) galaxies, obtained using the Keck I LRIS and Keck II DEIMOS spectrographs. Based on their g-i colors (taken with the CFHT MegaCam imager), physical positions on the sky, and radial velocities, we assign probabilities of dSph membership to each observed star. Using this information, the velocity dispersions, central masses and central densities of the dark matter halos are calculated for these objects, and compared with the properties of the Milky Way dSph population. We also measure the average metallicity ([Fe/H]) from the co-added spectra of member stars for each M31 dSph and find that they are consistent with the trend of decreasing [Fe/H] with luminosity observed in the Milky Way population. We find that three of our studied M31 dSphs appear as significant outliers in terms of their central velocity dispersion, And XIX, XXI and XXV, all of which have large half-light radii (>700 pc) and low velocity dispersions (sigma_v<5 km/s). In addition, And XXV has a mass-to-light ratio within its half-light radius of just [M/L]_{half}=10.3^{+7.0}_{-6.7}, making it consistent with a simple stellar system with no appreciable dark matter component within its 1 sigma uncertainties. We suggest that the structure of the dark matter halos of these outliers have been significantly altered by tides.

preprint2013arXiv

A Vast Thin Plane of Co-rotating Dwarf Galaxies Orbiting the Andromeda Galaxy

Dwarf satellite galaxies are thought to be the remnants of the population of primordial structures that coalesced to form giant galaxies like the Milky Way. An early analysis noted that dwarf galaxies may not be isotropically distributed around our Galaxy, as several are correlated with streams of HI emission, and possibly form co-planar groups. These suspicions are supported by recent analyses, and it has been claimed that the apparently planar distribution of satellites is not predicted within standard cosmology, and cannot simply represent a memory of past coherent accretion. However, other studies dispute this conclusion. Here we report the existence (99.998% significance) of a planar sub-group of satellites in the Andromeda galaxy, comprising approximately 50% of the population. The structure is vast: at least 400 kpc in diameter, but also extremely thin, with a perpendicular scatter <14.1 kpc (99% confidence). Radial velocity measurements reveal that the satellites in this structure have the same sense of rotation about their host. This finding shows conclusively that substantial numbers of dwarf satellite galaxies share the same dynamical orbital properties and direction of angular momentum, a new insight for our understanding of the origin of these most dark matter dominated of galaxies. Intriguingly, the plane we identify is approximately aligned with the pole of the Milky Way's disk and is co-planar with the Milky Way to Andromeda position vector. The existence of such extensive coherent kinematic structures within the halos of massive galaxies is a fact that must be explained within the framework of galaxy formation and cosmology.

preprint2013arXiv

Dynamical Modeling of NGC 6809: Selecting the best model using Bayesian Inference

The precise cosmological origin of globular clusters remains uncertain, a situation hampered by the struggle of observational approaches in conclusively identifying the presence, or not, of dark matter in these systems. In this paper, we address this question through an analysis of the particular case of NGC 6809. While previous studies have performed dynamical modeling of this globular cluster using a small number of available kinematic data, they did not perform appropriate statistical inference tests for the choice of best model description; such statistical inference for model selection is important since, in general, different models can result in significantly different inferred quantities. With the latest kinematic data, we use Bayesian inference tests for model selection and thus obtain the best fitting models, as well as mass and dynamic mass-to-light ratio estimates. For this, we introduce a new likelihood function that provides more constrained distributions for the defining parameters of dynamical models. Initially we consider models with a known distribution function, and then model the cluster using solutions of the spherically symmetric Jeans equation; this latter approach depends upon the mass density profile and anisotropy $β$ parameter. In order to find the best description for the cluster we compare these models by calculating their Bayesian evidence. We find smaller mass and dynamic mass-to-light ratio values than previous studies, with the best fitting Michie model for a constant mass-to-light ratio of $Υ= 0.90^{+0.14}_{-0.14}$ and $M_{\text{dyn}}=6.10^{+0.51}_{-0.88} \times 10^4 M_{\odot}$. We exclude the significant presence of dark matter throughout the cluster, showing that no physically motivated distribution of dark matter can be present away from the cluster core.

preprint2013arXiv

Inferring the Andromeda Galaxy's mass from its giant southern stream with Bayesian simulation sampling

M31 has a giant stream of stars extending far to the south and a great deal of other tidal debris in its halo, much of which is thought to be directly associated with the southern stream. We model this structure by means of Bayesian sampling of parameter space, where each sample uses an N-body simulation of a satellite disrupting in M31's potential. We combine constraints on stellar surface densities from the Isaac Newton Telescope survey of M31 with kinematic data and photometric distances. This combination of data tightly constrains the model, indicating a stellar mass at last pericentric passage of log(M_s / Msun) = 9.5+-0.1, comparable to the LMC. Any existing remnant of the satellite is expected to lie in the NE Shelf region beside M31's disk, at velocities more negative than M31's disk in this region. This rules out the prominent satellites M32 or NGC 205 as the progenitor, but an overdensity recently discovered in M31's NE disk sits at the edge of the progenitor locations found in the model. M31's virial mass is constrained in this model to be log(M200) = 12.3+-0.1, alleviating the previous tension between observational virial mass estimates and expectations from the general galactic population and the timing argument. The techniques used in this paper, which should be more generally applicable, are a powerful method of extracting physical inferences from observational data on tidal debris structures.

preprint2013arXiv

Matter Matters: Unphysical Properties of the Rh = ct Universe

It is generally agreed that there is matter in the universe and, in this paper, we show that the existence of matter is extremely problematic for the proposed Rh = ct universe. Considering a dark energy component with an equation of state of w=-1/3, it is shown that the presence of matter destroys the strict expansion properties that define the evolution of Rh = ct cosmologies, distorting the observational properties that are touted as its success. We further examine whether an evolving dark energy component can save this form of cosmological expansion in the presence of matter by resulting in an expansion consistent with a mean value of <w> = -1/3, finding that the presence of mass requires unphysical forms of the dark energy component in the early universe. We conclude that matter in the universe significantly limits the fundamental properties of the Rh = ct cosmology, and that novel, and unphysical, evolution of the matter component would be required to save it. Given this, Rh = ct cosmology is not simpler or more accurate description of the universe than prevailing cosmological models, and its presentation to date possesses significant flaws.

preprint2013arXiv

Phantom Energy and the Cosmic Horizon: Rh is still not a horizon!

There has been a recent spate of papers on the Cosmic Horizon, an apparently fundamental, although unrecognised, property of the universe. The misunderstanding of this horizon, it is claimed, demonstrates that our determination of the cosmological makeup of the universe is incorrect, although several papers have pointed out key flaws in these arguments. Here, we identify additional flaws in the most recent claims of the properties of the Cosmic Horizon in the presence of phantom energy, simply demonstrating that it does not act as a horizon, and that its limiting of our view of the universe is a trivial statement.

preprint2013arXiv

The large-scale structure of the halo of the Andromeda Galaxy Part I: global stellar density, morphology and metallicity properties

We present an analysis of the large-scale structure of the halo of the Andromeda galaxy, based on the Pan-Andromeda Archeological Survey (PAndAS), currently the most complete map of resolved stellar populations in any galactic halo. Despite copious substructure, the global halo populations follow closely power law profiles that become steeper with increasing metallicity. We divide the sample into stream-like populations and a smooth halo component. Fitting a three-dimensional halo model reveals that the most metal-poor populations ([Fe/H]<-1.7) are distributed approximately spherically (slightly prolate with ellipticity c/a=1.09+/-0.03), with only a relatively small fraction (42%) residing in discernible stream-like structures. The sphericity of the ancient smooth component strongly hints that the dark matter halo is also approximately spherical. More metal-rich populations contain higher fractions of stars in streams (86% for [Fe/H]>-0.6). The space density of the smooth metal-poor component has a global power-law slope of -3.08+/-0.07, and a non-parametric fit shows that the slope remains nearly constant from 30kpc to 300kpc. The total stellar mass in the halo at distances beyond 2 degrees is 1.1x10^10 Solar masses, while that of the smooth component is 3x10^9 Solar masses. Extrapolating into the inner galaxy, the total stellar mass of the smooth halo is plausibly 8x10^9 Solar masses. We detect a substantial metallicity gradient, which declines from [Fe/H]=-0.7 at R=30kpc to [Fe/H]=-1.5 at R=150kpc for the full sample, with the smooth halo being 0.2dex more metal poor than the full sample at each radius. While qualitatively in-line with expectations from cosmological simulations, these observations are of great importance as they provide a prototype template that such simulations must now be able to reproduce in quantitative detail.

preprint2013arXiv

The masses of Local Group dwarf spheroidal galaxies: The death of the universal mass profile

We investigate the claim that all dwarf spheroidal galaxies (dSphs) reside within halos that share a common, universal mass profile as has been derived for dSphs of the Galaxy. By folding in kinematic information for 25 Andromeda dSphs, more than doubling the previous sample size, we find that a singular mass profile can not be found to fit all the observations well. Further, the best-fit dark matter density profile measured for solely the Milky Way dSphs is marginally discrepant (at just beyond the 1 sigma level) with that of the Andromeda dSphs, where a profile with lower maximum circular velocity, and hence mass, is preferred. The agreement is significantly better when three extreme Andromeda outliers, And XIX, XXI and XXV, all of which have large half-light radii (>600pc) and low velocity dispersions (sigma_v < 5km/s) are omitted from the sample. We argue that the unusual properties of these outliers are likely caused by tidal interactions with the host galaxy.

preprint2013arXiv

The PAndAS view of the Andromeda satellite system - I. A Bayesian search for dwarf galaxies using spatial and color-magnitude information

We present a generic algorithm to search for dwarf galaxies in photometric catalogs and apply it to the Pan-Andromeda Archaeological Survey (PAndAS). The algorithm is developed in a Bayesian framework and, contrary to most dwarf-galaxy-search codes, makes use of both the spatial and color-magnitude information of sources in a probabilistic approach. Accounting for the significant contamination from the Milky Way foreground and from the structured stellar halo of the Andromeda galaxy, we recover all known dwarf galaxies in the PAndAS footprint with high significance, even for the least luminous ones. Some Andromeda globular clusters are also recovered and, in one case, discovered. We publish a list of the 143 most significant detections yielded by the algorithm. The combined properties of the 39 most significant isolated detections show hints that at least some of these trace genuine dwarf galaxies, too faint to be individually detected. Follow-up observations by the community are mandatory to establish which are real members of the Andromeda satellite system. The search technique presented here will be used in an upcoming contribution to determine the PAndAS completeness limits for dwarf galaxies. Although here tuned to the search of dwarf galaxies in the PAndAS data, the algorithm can easily be adapted to the search for any localised overdensity whose properties can be modeled reliably in the parameter space of any catalog.

preprint2013arXiv

The Three-Dimensional Structure of the M31 Satellite System; Strong Evidence for an Inhomogeneous Distribution of Satellites

We undertake an investigation into the spatial structure of the M31 satellite system utilizing the distance distributions presented in a previous publication. These distances make use of the unique combination of depth and spatial coverage of the Pan-Andromeda Archaeological Survey (PAndAS) to provide a large, homogeneous sample consisting of 27 of M31's satellites, as well as M31 itself. We find that the satellite distribution, when viewed as a whole, is no more planar than one would expect from a random distribution of equal size. A disk consisting of 15 of the satellites is however found to be highly significant, and strikingly thin, with a root-mean-square thickness of just $12.34^{+0.75}_{-0.43}$ kpc. This disk is oriented approximately edge on with respect to the Milky Way and almost perpendicular to the Milky Way disk. It is also roughly orthogonal to the disk like structure regularly reported for the Milky Way satellite system and in close alignment with M31's Giant Stellar Stream. A similar analysis of the asymmetry of the M31 satellite distribution finds that it is also significantly larger than one would expect from a random distribution. In particular, it is remarkable that 20 of the 27 satellites most likely lie on the Milky Way side of the galaxy, with the asymmetry being most pronounced within the satellite subset forming the aforementioned disk. This lopsidedness is all the more intriguing in light of the apparent orthogonality observed between the satellite disk structures of the Milky Way and M31.

preprint2012arXiv

A Bayesian Approach to Locating the Red Giant Branch Tip Magnitude (Part II); Distances to the Satellites of M31

In `A Bayesian Approach to Locating the Red Giant Branch Tip Magnitude (PART I),' a new technique was introduced for obtaining distances using the TRGB standard candle. Here we describe a useful complement to the technique with the potential to further reduce the uncertainty in our distance measurements by incorporating a matched-filter weighting scheme into the model likelihood calculations. In this scheme, stars are weighted according to their probability of being true object members. We then re-test our modified algorithm using random-realization artificial data to verify the validity of the generated posterior probability distributions (PPDs) and proceed to apply the algorithm to the satellite system of M31, culminating in a 3D view of the system. Further to the distributions thus obtained, we apply a satellite-specific prior on the satellite distances to weight the resulting distance posterior distributions, based on the halo density profile. Thus in a single publication, using a single method, a comprehensive coverage of the distances to the companion galaxies of M31 is presented, encompassing the dwarf spheroidals Andromedas I - III, V, IX-XXVII and XXX along with NGC147, NGC 185, M33 and M31 itself. Of these, the distances to Andromeda XXIV - XXVII and Andromeda XXX have never before been derived using the TRGB. Object distances are determined from high-resolution tip magnitude posterior distributions generated using the Markov Chain Monte Carlo (MCMC) technique and associated sampling of these distributions to take into account uncertainties in foreground extinction and the absolute magnitude of the TRGB as well as photometric errors. The distance PPDs obtained for each object both with, and without the aforementioned prior are made available to the reader in tabular form...

preprint2012arXiv

Does the Sagittarius Stream constrain the Milky Way halo to be triaxial?

Recent analyses of the stellar stream of the Sagittarius dwarf galaxy have claimed that the kinematics and three-dimensional location of the M-giant stars in this structure constrain the dark matter halo of our Galaxy to possess a triaxial shape that is extremely flattened, being essentially an oblate ellipsoid oriented perpendicular to the Galactic disk. Using a new stream-fitting algorithm, based on a Markov Chain Monte Carlo procedure, we investigate whether this claim remains valid if we allow the density profile of the Milky Way halo greater freedom. We find stream solutions that fit the leading and trailing arms of this structure even in a spherical halo, although this would need a rising Galactic rotation curve at large Galactocentric radius. However, the required rotation curve is not ruled out by current constraints. It appears therefore that for the Milky Way, halo triaxiality, despite its strong theoretical motivation, is not required to explain the Sagittarius stream. This degeneracy between triaxiality and the halo density profile suggests that in future endeavors to model this structure, it will be advantageous to relax the strict analytic density profiles that have been used to date.

preprint2012arXiv

How does the Hubble Sphere limit our view of the Universe?

It has recently been claimed that the Hubble Sphere represents a previously unknown limit to our view of the universe, with light we detect today coming from a proper distance less than this "Cosmic Horizon" at the present time. By considering the paths of light rays in several cosmologies, we show that this claim is not generally true. In particular, in cosmologies dominated by phantom energy (with an equation of state of ω< -1) the proper distance to the Hubble Sphere decreases, and light rays can cross it more than once in both directions; such behaviour further diminishes the claim that the Hubble Sphere is a fundamental, but unrecognised, horizon in the universe.

preprint2012arXiv

Kinematics of the stellar halo and the mass distribution of the Milky Way using BHB stars

Here we present a kinematic study of the Galactic halo out to a radius of $\sim$ 60 kpc, using 4664 blue horizontal branch (BHB) stars selected from the SDSS/SEGUE survey, to determine key dynamical properties. Using a maximum likelihood analysis, we determine the velocity dispersion profiles in spherical coordinates ($σ_{r}$, $σ_θ$, $σ_ϕ$) and the anisotropy profile ($β$). The radial velocity dispersion profile ($σ_{r}$) is measured out to a galactocentric radius of $r \sim 60$ kpc, but due to the lack of proper-motion information, $σ_θ$, $σ_ϕ$ and $β$ could only be derived directly out to $r \sim25$ kpc. From a starting value of $β\approx 0.5$ in the inner parts ($9<r/\kpc<12$), the profile falls sharply in the range $r \approx 13-18$ kpc, with a minimum value of $β=-1.2$ at $r=17$ kpc, rising sharply at larger radius. In the outer parts, in the range $25<r/\kpc<56$, we predict the profile to be roughly constant with a value of $β\approx 0.5$. The newly discovered kinematic anomalies are shown not to arise from halo substructures. We also studied the anisotropy profile of simulated stellar halos formed purely by accretion and found that they cannot reproduce the sharp dip seen in the data. From the Jeans equation, we compute the stellar rotation curve ($v_{\rm circ}$) of the Galaxy out to $r \sim 25$ kpc. The mass of the Galaxy within $r \lesssim 25$ kpc is determined to be $2.1 \times 10^{11}$ $M_{\sun}$, and with a 3-component fit to $v_{\rm circ}(r)$, we determine the virial mass of the Milky Way dark matter halo to be $M_{\rm vir} = 0.9 ^{+0.4}_{-0.3} \times 10^{12}$ $M_{\sun}$ ($R_{\rm vir} = 249^{+34}_{-31}$ kpc).

preprint2012arXiv

Mapping Growth and Gravity with Robust Redshift Space Distortions

Redshift space distortions caused by galaxy peculiar velocities provide a window onto the growth rate of large scale structure and a method for testing general relativity. We investigate through a comparison of N-body simulations to various extensions of perturbation theory beyond the linear regime, the robustness of cosmological parameter extraction, including the gravitational growth index, γ. We find that the Kaiser formula and some perturbation theory approaches bias the growth rate by 1-sigma or more relative to the fiducial at scales as large as k > 0.07 h/Mpc. This bias propagates to estimates of the gravitational growth index as well as Ω_m and the equation of state parameter and presents a significant challenge to modelling redshift space distortions. We also determine an accurate fitting function for a combination of line of sight damping and higher order angular dependence that allows robust modelling of the redshift space power spectrum to substantially higher k.

preprint2012arXiv

PAndAS in the mist: The stellar and gaseous mass within the halos of M31 and M33

Large scale surveys of the prominent members of the Local Group have provided compelling evidence for the hierarchical formation of massive galaxies, revealing a wealth of substructure that is thought to be the debris from ancient and on-going accretion events. In this paper, we compare two extant surveys of the M31-M33 subgroup of galaxies; the Pan-Andromeda Archaeological Survey (PAndAS) of the stellar structure, and a combination of observations of the HI gaseous content, detected at 21cm. Our key finding is a marked lack of spatial correlation between these two components on all scales, with only a few potential overlaps between stars and gas.The paucity of spatial correlation significantly restricts the analysis of kinematic correlations, although there does appear to the HI kinematically associated with the Giant Stellar Stream where it passes the disk of M31. These results demonstrate that that different processes must significantly influence the dynamical evolution of the stellar and HI components of substructures, such as ram pressure driving gas away from a purely gravitational path. Detailed modelling of the offset between the stellar and gaseous substructure will provide a determination of the properties of the gaseous halo of M31 and M33.

preprint2012arXiv

The Alcubierre Warp Drive: On the Matter of Matter

The Alcubierre warp drive allows a spaceship to travel at an arbitrarily large global velocity by deforming the spacetime in a bubble around the spaceship. Little is known about the interactions between massive particles and the Alcubierre warp drive, or the effects of an accelerating or decelerating warp bubble. We examine geodesics representative of the paths of null and massive particles with a range of initial velocities from -c to c interacting with an Alcubierre warp bubble travelling at a range of globally subluminal and superluminal velocities on both constant and variable velocity paths. The key results for null particles match what would be expected of massive test particles as they approach +/- c. The increase in energy for massive and null particles is calculated in terms of v_s, the global ship velocity, and v_p, the initial velocity of the particle with respect to the rest frame of the origin/destination of the ship. Particles with positive v_p obtain extremely high energy and velocity and become "time locked" for the duration of their time in the bubble, experiencing very little proper time between entering and eventually leaving the bubble. When interacting with an accelerating bubble, any particles within the bubble at the time receive a velocity boost that increases or decreases the magnitude of their velocity if the particle is moving towards the front or rear of the bubble respectively. If the bubble is decelerating, the opposite effect is observed. Thus Eulerian matter is unaffected by bubble accelerations/decelerations. The magnitude of the velocity boosts scales with the magnitude of the bubble acceleration/deceleration.

preprint2012arXiv

The cosmic history of the spin of dark matter haloes within the large scale structure

We use N-body simulations to investigate the evolution of the orientation and magnitude of dark matter halo angular momentum within the large scale structure since z=3. We look at the evolution of the alignment of halo spins with filaments and with each other, as well as the spin parameter, which is a measure of the magnitude of angular momentum. It was found that the angular momentum vectors of dark matter haloes at high redshift have a weak tendency to be orthogonal to filaments and high mass haloes have a stronger orthogonal alignment than low mass haloes. Since z=1, the spins of low mass haloes have become weakly aligned parallel to filaments, whereas high mass haloes kept their orthogonal alignment. This recent parallel alignment of low mass haloes casts doubt on tidal torque theory as the sole mechanism for the build up of angular momentum. We see evidence for bulk flows and the broadening of filaments over time in the alignments of halo spin and velocities. We find a significant alignment of the spin of neighboring dark matter haloes only at very small separations, $r<0.3$Mpc/h, which is driven by substructure. A correlation of the spin parameter with halo mass is confirmed at high redshift.

preprint2012arXiv

The Lives of High Redshift Mergers

We present a comparative study of recent works on merger-timescales with dynamical friction and find a strong contrast between idealized/isolated mergers (Boylan-Kolchin et al. 2008) and mergers from a cosmological volume (Jiang et al. 2008). Our study measures the duration of mergers in a cosmological N-body simulation of dark matter, with emphasis on higher redshifts (z < 10) and a lower mass range. In our analysis we consider and compare two merger definitions; tidal disruption and coalescence. We find that the merger-time formula proposed by Jiang et al. (2008) describes our results well and conclude that cosmologically motivated merger-time formulae provide a more versatile and statistically robust approximation for practical applications such as semi-analytic/hybrid models.

preprint2012arXiv

The X-ray Transient 2XMMi J003833.3+402133: A Candidate Magnetar at High Galactic Latitude

We present detailed analysis of the transient X-ray source 2XMMi J003833.3+402133 detected by XMM-Newton in January 2008 during a survey of M 31. The X-ray spectrum is well fitted by either a steep power law plus a blackbody model or a double blackbody model. Prior observations with XMM-Newton, Chandra, Swift and ROSAT spanning 1991 to 2007, as well as an additional Swift observation in 2011, all failed to detect this source. No counterpart was detected in deep optical imaging with the Canada France Hawaii Telescope down to a 3sigma lower limit of g = 26.5 mag. This source has previously been identified as a black hole X-ray binary in M 31. While this remains a possibility, the transient behaviour, X-ray spectrum, and lack of an optical counterpart are equally consistent with a magnetar interpretation for 2XMMi J003833.3+402133. The derived luminosity and blackbody emitting radius at the distance of M 31 argue against an extragalactic location, implying that if it is indeed a magnetar it is located within the Milky Way but 22deg out of the plane. The high Galactic latitude could be explained if 2XMMi J003833.3+402133 were an old magnetar, or if its progenitor was a runaway star that traveled away from the plane prior to going supernova.

preprint2012arXiv

Unearthing Foundations of a Cosmic Cathedral: Searching the Stars for M33's Halo

We use data from the Pan-Andromeda Archaeological Survey (PAndAS) to search for evidence of an extended halo component belonging to M33 (the Triangulum Galaxy). We identify a population of red giant branch (RGB) stars at large radii from M33's disk whose connection to the recently discovered extended "disk substructure" is ambiguous, and which may represent a "bona-fide" halo component. After first correcting for contamination from the Milky Way foreground population and misidentified background galaxies, we average the radial density of RGB candidate stars over circular annuli centered on the galaxy and away from the disk substructure. We find evidence of a low-luminosity, centrally concentrated component that is everywhere in our data fainter than mu_V ~ 33 mag arcsec^(-2). The scale length of this feature is not well constrained by our data, but it appears to be of order r_exp ~ 20 kpc; there is weak evidence to suggest it is not azimuthally symmetric. Inspection of the overall CMD for this region that specifically clips out the disk substructure reveals that this residual RGB population is consistent with an old population with a photometric metallicity of around [Fe/H] ~ -2 dex, but some residual contamination from the disk substructure appears to remain. We discuss the likelihood that our findings represent a bona-fide halo in M33, rather than extended emission from the disk substructure. We interpret our findings in terms of an upper limit to M33's halo that is a few percent of its total luminosity, although its actual luminosity is likely much less.

preprint2011arXiv

A New Collisional Ring Galaxy at z = 0.111: Auriga's Wheel

We report the serendipitous discovery of a collision ring galaxy, identified as 2MASX J06470249+4554022, which we have dubbed 'Auriga's Wheel', found in a SUPRIME-CAM frame as part of a larger Milky Way survey. This peculiar class of galaxies is the result of a near head-on collision between typically, a late type and an early type galaxy. Subsequent GMOS-N long-slit spectroscopy has confirmed both the relative proximity of the components of this interacting pair and shown it to have a redshift of 0.111. Analysis of the spectroscopy reveals that the late type galaxy is a LINER class Active Galactic Nuclei while the early type galaxy is also potentially an AGN candidate, this is very uncommon amongst known collision ring galaxies. Preliminary modeling of the ring finds an expansion velocity of ~200 kms^-1 consistent with our observations, making the collision about 50 Myr old. The ring currently has a radius of about 10 kpc and a bridge of stars and gas is also visible connecting the two galaxies.

preprint2011arXiv

AAOmega spectroscopy of 29 351 stars in fields centered on ten Galactic globular clusters

Galactic globular clusters have been pivotal in our understanding of many astrophysical phenomena. Here we publish the extracted stellar parameters from a recent large spectroscopic survey of ten globular clusters. A brief review of the project is also presented. Stellar parameters have been extracted from individual stellar spectra using both a modified version of the Radial Velocity Experiment (RAVE) pipeline and a pipeline based on the parameter estimation method of RAVE. We publish here all parameters extracted from both pipelines. We calibrate the metallicity and convert this to [Fe/H] for each star and, furthermore, we compare the velocities and velocity dispersions of the Galactic stars in each field to the Besançon Galaxy model. We find that the model does not correspond well with the data, indicating that the model is probably of little use for comparisons with pencil beam survey data such as this.

preprint2011arXiv

Gravitational Lensing with Three-Dimensional Ray Tracing

High redshift sources suffer from magnification or demagnification due to weak gravitational lensing by large scale structure. One consequence of this is that the distance-redshift relation, in wide use for cosmological tests, suffers lensing-induced scatter which can be quantified by the magnification probability distribution. Predicting this distribution generally requires a method for ray-tracing through cosmological N-body simulations. However, standard methods tend to apply the multiple thin-lens approximation. In an effort to quantify the accuracy of these methods, we develop an innovative code that performs ray-tracing without the use of this approximation. The efficiency and accuracy of this computationally challenging approach can be improved by careful choices of numerical parameters; therefore, the results are analysed for the behaviour of the ray-tracing code in the vicinity of Schwarzschild and Navarro-Frenk-White lenses. Preliminary comparisons are drawn with the multiple lens-plane ray-bundle method in the context of cosmological mass distributions for a source redshift of $z_{s}=0.5$.

preprint2011arXiv

PAndAS' progeny: extending the M31 dwarf galaxy cabal

We present the discovery of five new dwarf galaxies, Andromeda XXIII-XXVII, located in the outer halo of M31. These galaxies were found in the second year of data from the Pan-Andromeda Archaeological Survey (PAndAS) of the M31/M33 subgroup. This survey now provides an almost complete panoramic view of the M31 halo out to an average projected radius of ~150 kpc. The metal-poor stellar density map for this whole region serves, not only as an illustration of the discovery space for satellite galaxies, but also gives a birds-eye view of the ongoing assembly process of an L* disk galaxy. Four of the new dwarfs appear as well-defined spatial over-densities of stars lying on the expected locus of metal-poor red giant branch stars at the distance of M31. The fifth over-density, And XXVII, is embedded in an extensive stream of such stars and is possibly the remnant of a strong tidal disruption event. All five satellites have metallicities and luminosities typical of dwarf spheroidal galaxies and continue the trend whereby the brighter dwarf spheroidal satellites of M31 generally have much larger half-light radii than their Milky Way counterparts. With an extended sample of M31 satellite galaxies we also revisit the spatial distribution of this population and find that, within the current limits of the PAndAS survey, the surface density of satellites is essentially constant out to 150 kpc. This corresponds to a radial density distribution of satellites varying as 1/r, a result seemingly in conflict with the predictions of cosmological simulations. (Abridged)

preprint2011arXiv

Stellar Streams as Probes of Dark Halo Mass and Morphology: A Bayesian Reconstruction

Tidal streams provide a powerful tool by means of which the matter distribution of the dark matter halos of their host galaxies can be studied. However, the analysis is not straightforward because streams do not delineate orbits, and for most streams, especially those in external galaxies, kinematic information is absent. We present a method wherein streams are fit with simple corrections made to possible orbits of the progenitor, using a Bayesian technique known as Parallel Tempering to efficiently explore the parameter space. We show that it is possible to constrain the shape of the host halo potential or its density distribution using only the projection of tidal streams on the sky, if the host halo is considered to be axisymmetric. By adding kinematic data or the circular velocity curve of the host to the fitting data, we are able to recover other parameters of the matter distribution such as its mass and profile. We test our method on several simulated low mass stellar streams and also explore the cases for which additional data are required.

preprint2011arXiv

The M33 Globular Cluster System with PAndAS Data: The Last Outer Halo Cluster?

We use CFHT/MegaCam data to search for outer halo star clusters in M33 as part of the Pan-Andromeda Archaeological Survey (PAndAS). This work extends previous studies out to a projected radius of 50 kpc and covers over 40 square degrees. We find only one new unambiguous star cluster in addition to the five previously known in the M33 outer halo (10 kpc <= r <= 50 kpc). Although we identify 2440 cluster candidates of various degrees of confidence from our objective image search procedure, almost all of these are likely background contaminants, mostly faint unresolved galaxies. We measure the luminosity, color and structural parameters of the new cluster in addition to the five previously-known outer halo clusters. At a projected radius of 22 kpc, the new cluster is slightly smaller, fainter and redder than all but one of the other outer halo clusters, and has g' ~ 19.9, (g'-i') ~ 0.6, concentration parameter c ~ 1.0, a core radius r_c ~ 3.5 pc, and a half-light radius r_h ~ 5.5 pc. For M33 to have so few outer halo clusters compared to M31 suggests either tidal stripping of M33's outer halo clusters by M31, or a very different, much calmer accretion history of M33.

preprint2011arXiv

The Star Formation History and Dust Content in the Far Outer Disc of M31

We present a detailed analysis of two fields located 26 kpc (~5 scalelengths) from the centre of M31. One field samples the major axis populations--the Outer Disc field--while the other is offset by ~18' and samples the Warp in the stellar disc. The CMDs based on HST/ACS imaging reach old main-sequence turn-offs (~12.5 Gyr). We apply the CMD-fitting technique to the Warp field to reconstruct the star formation history (SFH). We find that after undergoing roughly constant SF until about 4.5 Gyr ago, there was a rapid decline in activity and then a ~1.5 Gyr lull, followed by a strong burst lasting 1.5 Gyr and responsible for 25% of the total stellar mass in this field. This burst appears to be accompanied by a decline in metallicity which could be a signature of the inflow of metal-poor gas. The onset of the burst (~3 Gyr ago) corresponds to the last close passage of M31 and M33 as predicted by detailed N-body modelling, and may have been triggered by this event. We reprocess the deep M33 outer disc field data of Barker et al. (2011) in order to compare consistently-derived SFHs. This reveals a similar duration burst that is exactly coeval with that seen in the M31 Warp field, lending further support to the interaction hypothesis. The complex SFHs and the smoothly-varying age-metallicity relations suggest that the stellar populations observed in the far outer discs of both galaxies have largely formed in situ rather than migrated from smaller galactocentric radii. The strong differential reddening affecting the CMD of the Outer Disc field prevents derivation of the SFH. Instead, we quantify this reddening and find that the fine-scale distribution of dust precisely follows that of the HI gas. This indicates that the outer HI disc of M31 contains a substantial amount of dust and therefore suggests significant metal enrichment in these parts, consistent with inferences from our CMD analysis.

preprint2010arXiv

AAOmega Observations of 47 Tucanae: Evidence for a Past Merger?

The globular cluster 47 Tucanae is well studied but it has many characteristics that are unexplained, including a significant rise in the velocity dispersion profile at large radii, indicating the exciting possibility of two distinct kinematic populations. In this Letter we employ a Bayesian approach to the analysis of the largest available spectral dataset of 47 Tucanae to determine whether this apparently two-component population is real. Assuming the two models were equally likely before taking the data into account, we find that the evidence favours the two-component population model by a factor of ~3x10^7. Several possible explanations for this result are explored, namely the evaporation of low-mass stars, a hierarchical merger, extant remnants of two initially segregated populations, and multiple star formation epochs. We find the most compelling explanation for the two-component velocity distribution is that 47 Tuc formed as two separate populations arising from the same proto-cluster cloud which merged <7.3 +/- 1.5 Gyr ago. This may also explain the extreme rotation, low mass-to-light ratio and mixed stellar populations of this cluster.

preprint2010arXiv

Halo globular clusters observed with AAOmega: dark matter content, metallicity and tidal heating

Globular clusters have proven to be essential to our understanding of many important astrophysical phenomena. Here we analyse spectroscopic observations of ten Halo globular clusters to determine their dark matter content, their tidal heating by the Galactic disc and halo, describe their metallicities and the likelihood that Newtonian dynamics explain their kinematics. We analyse a large number of members in all clusters, allowing us to address all these issues together, and we have included NGC 288 and M30 to overlap with previous studies. We find that any flattening of the velocity dispersion profiles in the outer regions of our clusters can be explained by tidal heating. We also find that all our GCs have M/L_V < 5, therefore, we infer the observed dynamics do not require dark matter, or a modification of gravity. We suggest that the lack of tidal heating signatures in distant clusters indicates the Halo is not triaxial. The isothermal rotations of each cluster are measured, with M4 and NGC 288 exhibiting rotation at a level of 0.9 +/- 0.1 km/s and 0.25 +/- 0.15 km/s, respectively. We also indirectly measure the tidal radius of NGC 6752, determining a more realistic figure for this cluster than current literature values. Lastly, an unresolved and intriguing puzzle is uncovered with regard to the cooling of the outer regions of all ten clusters.

preprint2010arXiv

Modelling of the Complex CASSOWARY/SLUGS Gravitational Lenses

We present the first high-resolution images of CSWA 31, a gravitational lens system observed as part of the SLUGS (Sloan Lenses Unravelled by Gemini Studies) program. These systems exhibit complex image structure with the potential to strongly constrain the mass distribution of the massive lens galaxies, as well as the complex morphology of the sources. In this paper, we describe the strategy used to reconstruct the unlensed source profile and the lens galaxy mass profiles. We introduce a prior distribution over multi-wavelength sources that is realistic as a representation of our knowledge about the surface brightness profiles of galaxies and groups of galaxies. To carry out the inference computationally, we use Diffusive Nested Sampling, an efficient variant of Nested Sampling that uses Markov Chain Monte Carlo (MCMC) to sample the complex posterior distributions and compute the normalising constant. We demonstrate the efficacy of this approach with the reconstruction of the group-group gravitational lens system CSWA 31, finding the source to be composed of five merging spiral galaxies magnified by a factor of 13.

preprint2010arXiv

No evidence for internal rotation in the remnant core of the Sagittarius dwarf

We have conducted a spectroscopic survey of the inner regions of the Sagittarius (Sgr) dwarf galaxy using the AAOmega spectrograph on the Anglo-Australian Telescope. We determine radial velocities for over 1800 Sgr star members in 6 fields that cover an area 18.84 deg^2, with a typical accuracy of ~2 km/s. Motivated by recent numerical models of the Sgr tidal stream that predict a substantial amount of rotation in the dwarf remnant core, we compare the kinematic data against N-body models that simulate the stream progenitor as (i) a pressure-supported, mass-follows-light system, and (ii) a late-type, rotating disc galaxy embedded in an extended dark matter halo. We find that the models with little, or no intrinsic rotation clearly yield a better match to the mean line-of-sight velocity in all surveyed fields, but fail to reproduce the shape of the line-of-sight velocity distribution. This result rules out models wherein the prominent bifurcation observed in the leading tail of the Sgr stream was caused by a transfer from intrinsic angular momentum from the progenitor satellite into the tidal stream. It also implies that the trajectory of the young tidal tails has not been affected by internal rotation in the progenitor system. Our finding indicates that new, more elaborate dynamical models, in which the dark and luminous components are treated independently, are necessary for simultaneously reproducing both the internal kinematics of the Sgr dwarf and the available data for the associated tidal stream.

preprint2010arXiv

The photometric properties of a vast stellar substructure in the outskirts of M33

We have surveyed $\sim40$sq.degrees surrounding M33 with CFHT MegaCam in the g and i filters, as part of the Pan-Andromeda Archaeological Survey. Our observations are deep enough to resolve the top 4mags of the red giant branch population in this galaxy. We have previously shown that the disk of M33 is surrounded by a large, irregular, low-surface brightness substructure. Here, we quantify the stellar populations and structure of this feature using the PAndAS data. We show that the stellar populations of this feature are consistent with an old population with $<[Fe/H]>\sim-1.6$dex and an interquartile range in metallicity of $\sim0.5$dex. We construct a surface brightness map of M33 that traces this feature to $μ_V\simeq33$mags\,arcsec$^{-2}$. At these low surface brightness levels, the structure extends to projected radii of $\sim40$kpc from the center of M33 in both the north-west and south-east quadrants of the galaxy. Overall, the structure has an "S-shaped" appearance that broadly aligns with the orientation of the HI disk warp. We calculate a lower limit to the integrated luminosity of the structure of $-12.7\pm0.5$mags, comparable to a bright dwarf galaxy such as Fornax or AndII and slightly less than $1\$ of the total luminosity of M33. Further, we show that there is tentative evidence for a distortion in the distribution of young stars near the edge of the HI disk that occurs at similar azimuth to the warp in HI. The data also hint at a low-level, extended stellar component at larger radius that may be a M33 halo component. We revisit studies of M33 and its stellar populations in light of these new results, and we discuss possible formation scenarios for the vast stellar structure. Our favored model is that of the tidal disruption of M33 in its orbit around M31.

preprint2010arXiv

Through the Looking Glass: Why the "Cosmic Horizon" is not a horizon

The present standard model of cosmology, $Λ$CDM, contains some intriguing coincidences. Not only are the dominant contributions to the energy density approximately of the same order at the present epoch, but we note that contrary to the emergence of cosmic acceleration as a recent phenomenon, the time averaged value of the deceleration parameter over the age of the universe is nearly zero. Curious features like these in $Λ$CDM give rise to a number of alternate cosmologies being proposed to remove them, including models with an equation of state w = -1/3. In this paper, we examine the validity of some of these alternate models and we also address some persistent misconceptions about the Hubble sphere and the event horizon that lead to erroneous conclusions about cosmology.

preprint2009arXiv

Fractal Bubble Cosmology: A concordant cosmological model?

The Fractal Bubble model has been proposed as a viable cosmology that does not require dark energy to account for cosmic acceleration, but rather attributes its observational signature to the formation of structure. In this paper it is demonstrated that, in contrast to previous findings, this model is not a good fit to cosmological supernovae data; there is significant tension in the best fit parameters obtained from different samples, whereas LCDM is able to fit all datasets consistently. Furthermore, the concordance between galaxy clustering scales and data from the cosmic microwave background is not achieved with the most recent supernova compilations. The validity of the FB formalism as a sound cosmological model is further challenged as it is shown that previous studies of this model achieve concordance by requiring a value for the present day Hubble constant that is derived from supernovae data containing an arbitrary distance normalisation.

preprint2009arXiv

Lyman alpha absorbers in motion: consequences of gravitational lensing for the cosmological redshift drift experiment

The evolution of the expansion rate of the Universe results in a drift in the redshift of distant sources over time. A measurement of this drift would provide us with a direct probe of expansion history. The Lyman alpha forest has been recognized as the best candidate for this experiment, but the signal would be weak and it will take next generation large telescopes coupled with ultra-stable high resolution spectrographs to reach the cm/s resolution required. One source of noise that has not yet been assessed is the transverse motion of Lyman alpha absorbers, which varies the gravitational potential in the line of sight and subsequently shifts the positions of background absorption lines. We examine the relationship between the pure cosmic signal and the observed redshift drift in the presence of moving Lyman alpha clouds, particularly the collapsed structures associated with Lyman limit systems (LLSs) and damped Lyman alpha systems (DLAs). Surprisingly, the peculiar velocities and peculiar accelerations both enter the expression, although the acceleration term stands alone as an absolute error, whilst the velocity term appears as a fractional noise component. An estimate of the magnitude of the noise reassures us that the motion of the Lyman alpha absorbers will not pose a threat to the detection of the signal.

preprint2009arXiv

PAndAS' cubs: discovery of two new dwarf galaxies in the surroundings of the Andromeda and Triangulum galaxies

We present the discovery of two new dwarf galaxies, Andromeda XXI and Andromeda XXII, located in the surroundings of the Andromeda and Triangulum galaxies (M31 and M33). These discoveries stem from the first year data of the Pan-Andromeda Archaeological Survey (PAndAS), a photometric survey of the M31/M33 group conducted with the Megaprime/MegaCam wide-field camera mounted on the Canada-France-Hawaii Telescope. Both satellites appear as spatial overdensities of stars which, when plotted in a color-magnitude diagram, follow metal-poor, [Fe/H]=-1.8, red giant branches at the distance of M31/M33. Andromeda XXI is a moderately bright dwarf galaxy (M_V=-9.9+/-0.6), albeit with low surface brightness, emphasizing again that many relatively luminous M31 satellites still remain to be discovered. It is also a large satellite, with a half-light radius close to 1 kpc, making it the fourth largest Local Group dwarf spheroidal galaxy after the recently discovered Andromeda XIX, Andromeda II and Sagittarius around the Milky Way, and supports the trend that M31 satellites are larger than their Milky Way counterparts. Andromeda XXII is much fainter (M_V=-6.5+/-0.8) and lies a lot closer in projection to M33 than it does to M31 (42 vs. 224 kpc), suggesting that it could be the first Triangulum satellite to be discovered. Although this is a very exciting possibility in the context of a past interaction of M33 with M31 and the fate of its satellite system, a confirmation will have to await a good distance estimate to confirm its physical proximity to M33. Along with the dwarf galaxies found in previous surveys of the M31 surroundings, these two new satellites bring the number of dwarf spheroidal galaxies in this region to 20.

preprint2009arXiv

Testing Newtonian Gravity with AAOmega: Mass-To-Light Profiles and Metallicity Calibrations From 47 Tuc and M55

Globular clusters are an important test bed for Newtonian gravity in the weak-acceleration regime, which is vital to our understanding of the nature of the gravitational interaction. Recent claims have been made that the velocity dispersion profiles of globular clusters flatten out at large radii, despite an apparent paucity of dark matter in such objects, indicating the need for a modification of gravitational theories. We continue our investigation of this claim, with the largest spectral samples ever obtained of 47 Tucanae and M55. Furthermore, this large sample allows for an accurate metallicity calibration based on the equivalent widths of the calcium triplet lines and K band magnitude of the Tip of the Red Giant Branch. Assuming an isothermal distribution, the rotations of each cluster are also measured with both clusters exhibiting clear rotation signatures. The global velocity dispersions of NGC 121 and Kron 3, two globular clusters in the Small Magellanic Cloud, are also calculated. By applying a simple dynamical model to the velocity dispersion profiles of 47 Tuc and M55, we calculate their mass-to-light profiles, total masses and central velocity dispersions. We find no statistically significant flattening of the velocity dispersion at large radii for M55, and a marked increase in the profile of 47 Tuc for radii greater than approximately half the tidal radius. We interpret this increase as an evaporation signature, indicating that 47 Tuc is undergoing, or has undergone, core-collapse, but find no requirement for dark matter or a modification of gravitational theories in either cluster.

preprint2009arXiv

Testing Newtonian Gravity with AAOmega: Mass-to-Light Profiles of Four Globular Clusters

Testing Newtonian gravity in the weak-acceleration regime is vital to our understanding of the nature of the gravitational interaction. It has recently been claimed that the velocity dispersion profiles of several globular clusters flatten out at large radii, reminiscent of galaxy rotation curves, even though globular clusters are thought to contain little or no dark matter. We investigate this claim, using AAOmega observations of four globular clusters, namely M22, M30, M53 and M68. M30, one such cluster that has had this claim made for its velocity dispersion, was included for comparison with previous studies. We find no statistically significant flattening of the velocity dispersion at large radii for any of our target clusters and therefore we infer the observed dynamics do not require that globular clusters are dark matter dominated, or a modification of gravity. Furthermore, by applying a simple dynamical model we determine the radial mass-to-light profiles for each cluster. The isothermal rotations of each cluster are also measured, with M22 exhibiting clear rotation, M68 possible rotation and M30 and M53 lacking any rotation, within the uncertainties.

preprint2009arXiv

The Adventures of the Rocketeer: Accelerated Motion Under the Influence of Expanding Space

It is well known that interstellar travel is bounded by the finite speed of light, but on very large scales any rocketeer would also need to consider the influence of cosmological expansion on their journey. This paper examines accelerated journeys within the framework of Friedmann- Lemaitre-Robertson-Walker universes, illustrating how the duration of a fixed acceleration sharply divides exploration over interstellar and intergalactic distances. Furthermore, we show how the universal expansion increases the difficulty of intergalactic navigation, with small uncertainties in cosmological parameters resulting in significantly large deviations. This paper also shows that, contrary to simplistic ideas, the motion of any rocketeer is indistinguishable from Newtonian gravity if the acceleration is kept small.

preprint2009arXiv

The remnants of galaxy formation from a panoramic survey of the region around M31

In hierarchical cosmological models, galaxies grow in mass through the continual accretion of smaller ones. The tidal disruption of these systems is expected to result in loosely bound stars surrounding the galaxy, at distances that reach $10 - 100$ times the radius of the central disk. The number, luminosity and morphology of the relics of this process provide significant clues to galaxy formation history, but obtaining a comprehensive survey of these components is difficult because of their intrinsic faintness and vast extent. Here we report a panoramic survey of the Andromeda galaxy (M31). We detect stars and coherent structures that are almost certainly remnants of dwarf galaxies destroyed by the tidal field of M31. An improved census of their surviving counterparts implies that three-quarters of M31's satellites brighter than $M_V < -6$ await discovery. The brightest companion, Triangulum (M33), is surrounded by a stellar structure that provides persuasive evidence for a recent encounter with M31. This panorama of galaxy structure directly confirms the basic tenets of the hierarchical galaxy formation model and reveals the shared history of M31 and M33 in the unceasing build-up of galaxies.

preprint2008arXiv

A Molecular Einstein Ring at z=4.12: Imaging the Dynamics of a Quasar Host Galaxy Through a Cosmic Lens

We present high-resolution (0.3") Very Large Array (VLA) imaging of the molecular gas in the host galaxy of the high redshift quasar PSS J2322+1944 (z=4.12). These observations confirm that the molecular gas (CO) in the host galaxy of this quasar is lensed into a full Einstein ring, and reveal the internal dynamics of the molecular gas in this system. The ring has a diameter of ~1.5", and thus is sampled over ~20 resolution elements by our observations. Through a model-based lens inversion, we recover the velocity gradient of the molecular reservoir in the quasar host galaxy of PSS J2322+1944. The Einstein ring lens configuration enables us to zoom in on the emission and to resolve scales down to ~1 kpc. From the model-reconstructed source, we find that the molecular gas is distributed on a scale of 5 kpc, and has a total mass of M(H2)=1.7 x 10^10 M_sun. A basic estimate of the dynamical mass gives M_dyn = 4.4 x 10^10 (sin i)^-2 M_sun, that is, only ~2.5 times the molecular gas mass, and ~30 times the black hole mass (assuming that the dynamical structure is highly inclined). The lens configuration also allows us to tie the optical emission to the molecular gas emission, which suggests that the active galactic nucleus (AGN) does reside within, but not close to the center of the molecular reservoir. Together with the (at least partially) disturbed structure of the CO, this suggests that the system is interacting. Such an interaction, possibly caused by a major `wet' merger, may be responsible for both feeding the quasar and fueling the massive starburst of 680 M_sun/yr in this system, in agreement with recently suggested scenarios of quasar activity and galaxy assembly in the early universe.

preprint2008arXiv

Can Early Dark Energy be Detected in Non-Linear Structure?

We present the first study of early dark energy cosmologies using N-body simulations to investigate the formation of non-linear structure. In contrast to expectations from semi-analytic approaches, we find that early dark energy does not imprint a unique signature on the statistics of non-linear structures. Investigating the non-linear power spectra and halo mass functions, we show that universal mass functions hold for early dark energy, making its presence difficult to distinguish from $Λ$CDM. Since early dark energy biases the baryon acoustic oscillation scale, the lack of discriminating power is problematic.

preprint2008arXiv

Cosmological Radar Ranging in an Expanding Universe

While modern cosmology, founded in the language of general relativity, is almost a century old, the meaning of the expansion of space is still being debated. In this paper, the question of radar ranging in an expanding universe is examined, focusing upon light travel times during the ranging; it has recently been claimed that this proves that space physically expands. We generalize the problem into considering the return journey of an accelerating rocketeer, showing that while this agrees with expectations of special relativity for an empty universe, distinct differences occur when the universe contains matter. We conclude that this does not require the expansion of space to be a physical phenomenon, rather that we cannot neglect the influence of matter, seen through the laws of general relativity, when considering motions on cosmic scales.

preprint2008arXiv

Halo Mass Functions in Early Dark Energy Cosmologies

We examine the linear density contrast at collapse time, $δ_c$ for large-scale structure in dynamical dark energy cosmologies, including models with early dark energy. Contrary to previous results, we find that as long as dark energy is homogeneous on small scales, $δ_c$ is insensitive to dark energy properties for parameter values fitting current data, including the case of early dark energy. This is significant since using the correct $δ_c$ is crucial for accurate Press-Schechter prediction of the halo mass function. Previous results have found an apparent failing of the extended Press-Schechter approach (Sheth-Tormen) for early dark energy. Our calculations demonstrate that with the correct $δ_c$ the accuracy of this approach is restored. We discuss the significance of this result for the halo mass function and examine what dark energy physics would be needed to cause significant change in $δ_c$, and the observational signatures this would leave.

preprint2005arXiv

A Keck/DEIMOS Kinematic Study of Andromeda IX: dark matter on the smallest galactic scales

We present the results of a kinematic survey of the dwarf spheroidal satellite of M31, And-IX, which appears to be the lowest surface brightness and also the faintest galaxy (M_V = -8.3) found to date. Using Keck/DEIMOS spectroscopic data, we have measured its velocity relative to M31, its velocity dispersion, and its metallicity. It exhibits a significant velocity dispersion sigma_v = 6.8+3.0-2.0km/s, which coupled with the low luminosity implies a very high mass to V-band light ratio, M/L ~ 93+120-50msun/lsun (M/L > 17msun/lsun at 99% confidence). Unless strong tidal forces have perturbed this system, this smallest of galaxies is a highly dark matter dominated system.

preprint2005arXiv

The Influence of Evolving Dark Energy on Cosmology

Observational evidence indicating that the expansion of the universe is accelerating has surprised cosmologists in recent years. Cosmological models have sought to explain this acceleration by incorporating `dark energy', of which the traditional cosmological constant is just one possible candidate. Several cosmological models involving an evolving equation of state of the dark energy have been proposed, as well as possible energy exchange to other components, such as dark matter. This paper summarises the forms of the most prominent models and discusses their implications for cosmology and astrophysics. Finally, this paper examines the current and future observational constraints on the nature of dark energy.

preprint2003arXiv

Nanolensing of gamma-ray bursts

All quasars vary in their optical flux on a time-scale of years, and it has been proposed that these variations are principally due to gravitational lensing by a cosmologically distributed population of planetary mass objects. This interpretation has implications for the observable properties of gamma-ray bursts (GRBs) -- as a source expands across the nano-arcsecond caustic network, variability is expected -- and data on GRBs can be used to test the proposed model of quasar variability. Here we employ an ultra-relativistic blast-wave model of the source, with no intrinsic variations, to study the effects of nanolensing on GRBs. Taken in isolation the light-curves of the caustic crossings are predictable, and we find that a subset of the predicted light-curves (the image-annihilating fold crossings) resemble the ``pulses'' which are commonly seen in long GRBs. Furthermore, for sources at high redshift the expected time between caustic crossings is of order seconds, comparable to the observed time between pulses. These points suggest that it might be possible to model some of the observed variations of GRBs in terms of nanolensing; however, our simulated light-curves exhibit a small depth of modulation compared to what is observed. This means that the GRB data do not significantly constrain the quasar nanolensing model; it also means that the simplest nanolensing model cannot explain the observed GRB ``pulses''. Viable nanolensing models for pulses probably require a large external beam shear. If a viable model can be constructed it would effect a considerable simplification in source modelling and, ironically, it would explain why no macro-lensed GRBs have been identified to date. (Abridged)

preprint2002arXiv

An investigation of gravitational lens determinations of H_o in quintessence cosmologies

There is growing evidence that the majority of the energy density of the universe is not baryonic or dark matter, rather it resides in an exotic component with negative pressure. The nature of this `quintessence' influences our view of the universe, modifying angular diameter and luminosity distances. Here, we examine the influence of a quintessence component upon gravitational lens time delays. As well as a static quintessence component, an evolving equation of state is also considered. It is found that the equation of state of the quintessence component and its evolution influence the value of the Hubble's constant derived from gravitational lenses. However, the differences between evolving and non-evolving cosmologies are relatively small. We undertake a suite of Monte Carlo simulations to examine the potential constraints that can be placed on the universal equation of state from the monitoring of gravitational lens system, and demonstrate that at least an order of magnitude more lenses than currently known will have to be discovered and analysed to accurately probe any quintessence component.

preprint2002arXiv

Cosmology on a Mesh

An adaptive multi grid approach to simulating the formation of structure from collisionless dark matter is described. MLAPM (Multi-Level Adaptive Particle Mesh) is one of the most efficient serial codes available on the cosmological 'market' today. As part of Swinburne University's role in the development of the Square Kilometer Array, we are implementing hydrodynamics, feedback, and radiative transfer within the MLAPM adaptive mesh, in order to simulate baryonic processes relevant to the interstellar and intergalactic media at high redshift. We will outline our progress to date in applying the existing MLAPM to a study of the decay of satellite galaxies within massive host potentials.

preprint2001arXiv

Was SN1997ff at z~1.7 magnified by gravitational lensing?

The quest for the cosmological parameters has come to fruition with the identification of a number of supernovae at a redshift of $z\sim1$. Analyses of the brightness of these standard candles reveal that the Universe is dominated by a large cosmological constant. The recent identification of the $z\sim1.7$ SN1997ff in the northern Hubble Deep Field has provided further evidence for this cosmology. Here we examine the case for gravitational lensing of SN1997ff due to the presence of galaxies lying along our line of sight. We find that, while the alignment of SN1997ff with foreground masses was not favorable for it to be multiply imaged and strongly magnified, two galaxies did lie close enough to result in significant magnification: $μ\sim1.4$ for the case where these elliptical galaxies have velocity dispersion $200 {\rm km/s}$. Given the small difference between supernova brightnesses in different cosmologies, detailed modeling of the gravitational lensing properties of the intervening matter is therefore required before the true cosmological significance of SN1997ff can be deduced.

preprint1999arXiv

KECK HIRES Spectroscopy of APM 08279+5255

With an optical R-band magnitude of 15.2, the recently discovered z=3.911 BAL quasar APM 08279+5255 is an exceptionally bright high redshift source. Its brightness has allowed us to acquire a high signal-to-noise ratio (~100), high resolution (~6 km/s) spectrum using the HIRES echelle spectrograph on the 10-m Keck I telescope. Given the quality of the data, these observations provide an unprecedented view of associated and intervening absorption systems. Here we announce the availability of this spectrum to the general astronomical community and present a brief analysis of some of its main features.

preprint1997arXiv

Two Close-Separation Quasar-Quasar Pairs in the Large Bright Quasar Survey

We present photometric and spectroscopic observations of two close quasar-quasar pairs found in the Large Bright Quasar Survey (LBQS) (Hewett et al. 1995, AJ, 109, 1498). The two components of the 2153-2056 pair ($z=1.845$, $Δθ=7\farcs8$, $B=17.9$ and 21.3) have the same redshifts within the relatively large uncertainty, $Δv_{A-B} =-1100\pm1500 km s^{-1}$, of the observations. The quasars are most likely spatially coincident although the possibility that the pair is the result of gravitational lensing cannot be ruled out. The two components of 1148+0055 ($z=1.879$, $B=18.5$ and $z=1.409$, $B=21.1$, $Δθ=3\farcs9$) have disjoint redshifts and the pair has attracted some attention in the context of gravitational lensing (e.g. Wampler 1997, ApJ, 476, L55) following the independent discovery of the pair by Surdej and collaborators. Four close, $Δθ\le10\farcs0$, quasar-quasar pairs have now been discovered in the LBQS and we discuss the probability of identifying pairs with disjoint redshifts and of locating spatially coincident pairs from the systematic investigation of the well-defined quasar catalogue.