Researcher profile

Paul J. McMillan

Paul J. McMillan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

14 published item(s)

preprint2022arXiv

The velocity distribution of white dwarfs in Gaia EDR3

Using a penalised maximum likelihood we estimate, for the first time, the velocity distribution of white dwarfs in the Solar neighbourhood. Our sample consists of 129 675 white dwarfs within 500 pc in Gaia Early Data Release 3 The white dwarf velocity distributions reveal a similar structure to the rest of the Solar neighbourhood stars, reflecting that white dwarfs are subjected to the same dynamical processes. In the velocity distribution for three magnitude-binned subsamples we however find a novel structure at $(U, V) = (7, -19)$ km s$^{-1}$ in fainter samples, potentially related to the Coma Berenices stream. We also see a double-peaked feature in $U$-$W$ at $U \approx -30$ km s$^{-1}$ and in $V$-$W$ at $V \approx -20$ km s$^{-1}$ for fainter samples. We determine the velocity distribution and velocity moments as a function of absolute magnitude for two samples based on the bifurcation identified in Gaia Data Release 2 in the colour-magnitude diagram. The brighter, redder sequence has a larger velocity dispersion than the fainter, bluer sequence across all magnitudes. It is hard to reconcile this kinematic difference with a bifurcation caused purely by atmospheric composition, while it fits neatly with a significant age difference between the two sequences. Our results provide novel insights into the kinematic properties of white dwarfs and demonstrate the power of analytical techniques that work for the large fraction of stars that do not have measured radial velocities in the current era of large-scale astrometric surveys.

preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- I: Survey Description, Spectra and Radial Velocities

The Radial Velocity Experiment (RAVE) is a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A). The 6th and final data release (DR6 or FDR) is based on 518387 observations of 451783 unique stars. RAVE observations were taken between 12 April 2003 and 4 April 2013. Here we present the genesis, setup and data reduction of RAVE as well as wavelength-calibrated and flux-normalized spectra and error spectra for all observations in RAVE DR6. Furthermore, we present derived spectral classification and radial velocities for the RAVE targets, complemented by cross matches with Gaia DR2 and other relevant catalogs. A comparison between internal error estimates, variances derived from stars with more than one observing epoch and a comparison with radial velocities of Gaia DR2 reveals consistently that 68% of the objects have a velocity accuracy better than 1.4 km/s, while 95% of the objects have radial velocities better than 4.0 km/s. Stellar atmospheric parameters, abundances and distances are presented in subsequent publication. The data can be accessed via the RAVE Web (http://rave-survey.org) or the Vizier database.

preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- II: Stellar Atmospheric Parameters, Chemical Abundances and Distances

We present part 2 of the 6th and final Data Release (DR6 or FDR) of the Radial Velocity Experiment (RAVE), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A) and span the complete time frame from the start of RAVE observations on 12 April 2003 to their completion on 4 April 2013. In the second of two publications, we present the data products derived from 518387 observations of 451783 unique stars using a suite of advanced reduction pipelines focussing on stellar atmospheric parameters, in particular purely spectroscopically derived stellar atmospheric parameters (Teff, log(g), and the overall metallicity), enhanced stellar atmospheric parameters inferred via a Bayesian pipeline using Gaia DR2 astrometric priors, and asteroseismically calibrated stellar atmospheric parameters for giant stars based on asteroseismic observations for 699 K2 stars. In addition, we provide abundances of the elements Fe, Al, and Ni, as well as an overall [alpha/Fe] ratio obtained using a new pipeline based on the GAUGUIN optimization method that is able to deal with variable signal-to-noise ratios. The RAVE DR6 catalogs are cross matched with relevant astrometric and photometric catalogs, and are complemented by orbital parameters and effective temperatures based on the infrared flux method. The data can be accessed via the RAVE Web site (http://rave-survey.org) or the Vizier database.

preprint2016arXiv

Chemical separation of disc components using RAVE

We present evidence from the RAdial Velocity Experiment (RAVE) survey of chemically separated, kinematically distinct disc components in the solar neighbourhood. We apply probabilistic chemical selection criteria to separate our sample into $α$-low (`thin disc&#39;) and $α$-high (`thick disc&#39;) sequences. Using newly derived distances, which will be utilized in the upcoming RAVE DR5, we explore the kinematic trends as a function of metallicity for each of the disc components. For our $α$-low disc, we find a negative trend in the mean rotational velocity ($V_{\mathrmϕ}$) as a function of iron abundance ([Fe/H]). We measure a positive gradient $\partial V_{\mathrmϕ}$/$\partial$[Fe/H] for the $α$-high disc, consistent with results from high-resolution surveys. We also find differences between the $α$-low and $α$-high discs in all three components of velocity dispersion. We discuss the implications of an $α$-low, metal-rich population originating from the inner Galaxy, where the orbits of these stars have been significantly altered by radial mixing mechanisms in order to bring them into the solar neighbourhood. The probabilistic separation we propose can be extended to other data sets for which the accuracy in [$α$/Fe] is not sufficient to disentangle the chemical disc components a priori. For such datasets which will also have significant overlap with Gaia DR1, we can therefore make full use of the improved parallax and proper motion data as it becomes available to investigate kinematic trends in these chemical disc components.

preprint2016arXiv

The mass distribution and gravitational potential of the Milky Way

We present mass models of the Milky Way created to fit observational constraints and to be consistent with expectations from theoretical modelling. The method used to create these models is that demonstrated in McMillan (2011), and we improve on those models by adding gas discs to the potential, considering the effects of allowing the inner slope of the halo density profile to vary, and including new observations of maser sources in the Milky Way amongst the new constraints. We provide a best fitting model, as well as estimates of the properties of the Milky Way. Under the assumptions in our main model, we find that the Sun is $R_0 = (8.20\pm0.09)\,\mathrm{kpc}$ from the Galactic Centre, with the circular speed at the Sun being $v_0 = (232.8\pm3.0)\,\mathrm{km}\,\mathrm{s}^{-1}$; that the Galaxy has a total stellar mass of $(54.3\pm5.7)\times10^9\,{\rm M}_\odot$, a total virial mass of $(1.30 \pm 0.30)\times10^{12}\,{\rm M}_\odot$ and a local dark-matter density of $0.38\pm0.04\,\mathrm{GeV\,cm}^{-3}$, where the quoted uncertainties are statistical. These values are sensitive to our choice of priors and constraints. We investigate systematic uncertainties, which in some cases may be larger. For example, if we weaken our prior on $R_0$, we find it to be $(7.97\pm0.15)\,\mathrm{kpc}$ and that $v_0=(226.8\pm4.2)\,\mathrm{km}\,\mathrm{s}^{-1}$. We find that most of these properties, including the local dark-matter density, are remarkably insensitive to the assumed power-law density slope at the centre of the dark-matter halo. We find that it is unlikely that the local standard of rest differs significantly from that found under assumptions of axisymmetry. We have made code to compute the force from our potential, and to integrate orbits within it, publicly available.

preprint2015arXiv

Making action-angle disc models for Gaia

A brief review of recent work. I describe dynamical modelling of the Milky Way using action-angle coordinates. I explain what action-angle coordinates are, and what progress has been made in the past few years to ensuring they can be used in reasonably realistic Galactic potentials. I then describe recent modelling efforts, and progress they have made in constraining the potential of the Milky Way and the local dark matter density.

preprint2013arXiv

Analysing surveys of our Galaxy -- II. Determining the potential

We consider the problem of determining the Galaxy&#39;s gravitational potential from a star catalogue. We show that orbit-based approaches to this problem suffer from unacceptable numerical noise deriving from the use of only a finite number of orbits. An alternative approach, which requires an ability to determine the model&#39;s phase-space density at predetermined positions and velocities, has a level of numerical noise that lies well below the intrinsic uncertainty associated with the finite size of the catalogue analysed. A catalogue of 10000 stars brighter than V=17 and distributed over the sky at b>30 degrees enables us to determine the scaleheight of the disc that contributes to the potential with an uncertainty below 20pc if the catalogue gives proper motions, line-of-sight velocities and parallaxes with errors typical of the Gaia Catalogue, rising to 36pc if only proper motions are available. The uncertainty in the disc&#39;s scalelength is significantly smaller than 0.25kpc.

preprint2013arXiv

Extending the Hyades

We explore the implications of models of the Hyades moving group in which it has a resonant origin, for regions of the Galaxy beyond the Solar neighbourhood. We show that while models associated with different resonances can produce nearly identical substructure in the local velocity distribution, the velocity distribution away from the Solar neighbourhood has different properties for different models. In particular there is a variation between different models of where in Galactocentric radius the observed Hyades signal in velocity space is strongest, at a given Galactic azimuth. We note, however, that the uncertainties in currently available data, primarily due to uncertain distances to stars, hide these signatures rather effectively, meaning we are not yet able to determine which resonance is the cause of the Hyades.

preprint2011arXiv

Analysing surveys of our Galaxy I: basic astrometric data

We consider what is the best way to extract science from large surveys of the Milky Way galaxy. The diversity of data gathered in these surveys, together with our position within the Galaxy, imply that science must be extracted by fitting dynamical models to the data in the space of the observables. Models based on orbital tori promise to be superior for this task than traditional types of models, such as N-body models and Schwarzschild models. A formalism that allows such models to be fitted to data is developed and tested on pseudodata of varying richness.

preprint2011arXiv

Mass models of the Milky Way

We present a simple method for fitting parametrized mass models of the Milky Way to observational constraints. We take a Bayesian approach which allows us to take into account input from photometric and kinematic data, and expectations from theoretical modelling. This provides us with a best-fitting model, which is a suitable starting point for dynamical modelling. We also determine a probability density function on the properties of the model, which demonstrates that the mass distribution of the Galaxy remains very uncertain. For our choices of parametrization and constraints, we find disc scale lengths of 3.00 \pm 0.22 kpc and 3.29 \pm 0.56 kpc for the thin and thick discs respectively; a Solar radius of 8.29 \pm 0.16 kpc and a circular speed at the Sun of 239 \pm 5 km/s; a total stellar mass of 6.43 \pm 0.63 * 10^10 M_sun; a virial mass of 1.26 \pm 0.24 * 10^12 M_sun and a local dark matter density of 0.40 \pm 0.04 GeV/cm^3. We find some correlations between the best-fitting parameters of our models (for example, between the disk scale lengths and the Solar radius), which we discuss. The chosen disc scale-heights are shown to have little effect on the key properties of the model.

preprint2011arXiv

The Solar neighbourhood in angle coordinates: the Hyades moving group

I investigate the suggestion that the Hyades moving group in the Solar neighbourhood is the result of a recent inner Lindblad resonance. I use dynamical &#34;torus&#34; models of the Galaxy to understand the expected distribution of solar neighbourhood stars in angle coordinates for phase-mixed models and models which include a resonant component. I show that attempts to find the signatures of resonances in angle coordinates are strongly influenced by selection effects, including rather subtle effects associated with the relationship between action and angle for stars at a given point. These effects mean that one can not use simple tests to determine whether substructures seen in the Solar neighbourhood are associated with any given resonance.

preprint2009arXiv

The dangers of deprojection of proper motions

We re-examine the method of deprojection of proper motions, which has been used for finding the velocity ellipsoid of stars in the nearby Galaxy. This method is only legitimate if the lines of sight to the individual stars are uncorrelated with the stars&#39; velocities. Very simple models are used to show that spurious results similar to ones recently reported are obtained when velocity dispersion decreases with galactocentric radius in the expected way. A scheme that compensates for this bias is proposed.

preprint2009arXiv

The uncertainty in Galactic parameters

We reanalyse the measurements of parallax, proper motion, and line-of-sight velocity for 18 masers in high mass star-forming regions presented by Reid et al. (2009). We use a likelihood analysis to investigate the distance of the Sun from the Galactic centre, R_0, the rotational speed of the local standard of rest, v_0, and the peculiar velocity of the Sun, vsol, for various models of the rotation curve, and models which allow for a typical peculiar motion of the high mass star-forming regions. We find that these data are best fit by models with non-standard values for vsol or a net peculiar motion of the high mass star-forming regions. We argue that a correction to vsol is much more likely, and that these data support the conclusion of Binney (2009) that V_sol should be revised upwards from 5.2 km/s to 11 km/s. We find that the values of R_0 and v_0 that we determine are heavily dependent on the model we use for the rotation curve, with model-dependent estimates of R_0 ranging from 6.7 \pm 0.5kpc to 8.9 \pm 0.9kpc, and those of v_0 ranging from 200 \pm 20 km/s to 279 \pm 33 km/s. We argue that these data cannot be thought of as implying any particular values of R_0 or v_0. However, we find that v_0/R_0 is better constrained, lying in the range 29.9-31.6 km/s/kpc for all models but one.