• We present a new framework (MADE) that produces distance and age estimates by applying a Bayesian isochrone pipeline to a combination of photometric, astrometric and spectroscopic data. For giant stars, the framework can supplement these observational constraints with posterior predictive distributions for mass from a new Bayesian spectroscopic mass estimator. The new mass estimator is a Bayesian artificial neural network (ANN) that learns the relationship between a specified set of inputs and outputs based on a training set. Posterior predictive distributions for the outputs given new inputs are computed, taking into account input uncertainties, and uncertainties in the parameters of the ANN. MADE trains the ANN on stars with spectroscopic and asteroseismology data to enable posterior predictive distributions for present masses of giant stars to be evaluated given spectroscopic data. We apply MADE to $\sim10\,000$ red giants in the overlap between APO Galactic Evolution Experiment (APOGEE) and the Tycho-Gaia astrometric solution (TGAS). The ANN is trained on a subsample of these stars with new asteroseismology determinations of mass, and is able to predict the masses to a similar degree of uncertainty as the measurement uncertainty. In particular, it is able to reduce the uncertainty on those with the highest measurement uncertainty. Using these masses in the Bayesian isochrone pipeline along with photometric and astrometric data, we are able to obtain distance estimates with uncertainties of order $\sim 10\%$ and age estimates with uncertainties of order $\sim 20\%$. Our resulting catalogue clearly demonstrates the expected thick and thin disc components in the [M/H]-[$\alpha$/M] plane when examined by age.
  • In this paper we propose a novel method that provides contrastive explanations justifying the classification of an input by a black box classifier such as a deep neural network. Given an input we find what should be minimally and sufficiently present (viz. important object pixels in an image) to justify its classification and analogously what should be minimally and necessarily \emph{absent} (viz. certain background pixels). We argue that such explanations are natural for humans and are used commonly in domains such as health care and criminology. What is minimally but critically \emph{absent} is an important part of an explanation, which to the best of our knowledge, has not been touched upon by current explanation methods that attempt to explain predictions of neural networks. We validate our approach on three real datasets obtained from diverse domains; namely, a handwritten digits dataset MNIST, a large procurement fraud dataset and an fMRI brain imaging dataset. In all three cases, we witness the power of our approach in generating precise explanations that are also easy for human experts to understand and evaluate.
  • The Theia Collaboration: Celine Boehm, Alberto Krone-Martins, Antonio Amorim, Guillem Anglada-Escude, Alexis Brandeker, Frederic Courbin, Torsten Ensslin, Antonio Falcao, Katherine Freese, Berry Holl, Lucas Labadie, Alain Leger, Fabien Malbet, Gary Mamon, Barbara McArthur, Alcione Mora, Michael Shao, Alessandro Sozzetti, Douglas Spolyar, Eva Villaver, Conrado Albertus, Stefano Bertone, Herve Bouy, Michael Boylan-Kolchin, Anthony Brown, Warren Brown, Vitor Cardoso, Laurent Chemin, Riccardo Claudi, Alexandre C. M. Correia, Mariateresa Crosta, Antoine Crouzier, Francis-Yan Cyr-Racine, Mario Damasso, Antonio da Silva, Melvyn Davies, Payel Das, Pratika Dayal, Miguel de Val-Borro, Antonaldo Diaferio, Adrienne Erickcek, Malcolm Fairbairn, Morgane Fortin, Malcolm Fridlund, Paulo Garcia, Oleg Gnedin, Ariel Goobar, Paulo Gordo, Renaud Goullioud, Nigel Hambly, Nathan Hara, David Hobbs, Erik Hog, Andrew Holland, Rodrigo Ibata, Carme Jordi, Sergei Klioner, Sergei Kopeikin, Thomas Lacroix, Jacques Laskar, Christophe Le Poncin-Lafitte, Xavier Luri, Subhabrata Majumdar, Valeri Makarov, Richard Massey, Bertrand Mennesson, Daniel Michalik, Andre Moitinho de Almeida, Ana Mourao, Leonidas Moustakas, Neil Murray, Matthew Muterspaugh, Micaela Oertel, Luisa Ostorero, Angeles Perez-Garcia, Imants Platais, Jordi Portell i de Mora, Andreas Quirrenbach, Lisa Randall, Justin Read, Eniko Regos, Barnes Rory, Krzysztof Rybicki, Pat Scott, Jean Schneider, Jakub Scholtz, Arnaud Siebert, Ismael Tereno, John Tomsick, Wesley Traub, Monica Valluri, Matt Walker, Nicholas Walton, Laura Watkins, Glenn White, Dafydd Wyn Evans, Lukasz Wyrzykowski, Rosemary Wyse
    July 2, 2017 astro-ph.IM
    In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.
  • Using 17 chemical elements as a proxy for stellar DNA, we present a full phylogenetic study of stars in the solar neighbourhood. This entails applying a clustering technique that is widely used in molecular biology to construct an evolutionary tree from which three branches emerge. These are interpreted as stellar populations which separate in age and kinematics and can be thus attributed to the thin disk, the thick disk, and an intermediate population of probable distinct origin. We further find six lone stars of intermediate age that could not be assigned to any population with enough statistical significance. Combining the ages of the stars with their position on the tree, we are able to quantify the mean rate of chemical enrichment of each of the populations, and thus show in a purely empirical way that the star formation rate in the thick disk is much higher than in the thin disk. We are also able to estimate the relative contribution of dynamical processes such as radial migration and disk heating to the distribution of chemical elements in the solar neighbourhood. Our method offers an alternative approach to chemical tagging methods with the advantage of visualising the behaviour of chemical elements in evolutionary trees. This offers a new way to search for `common ancestors' that can reveal the origin of solar neighbourhood stars.
  • The distribution of Milky Way halo blue horizontal-branch (BHB) stars is examined using action-based extended distribution functions (EDFs) that describe the locations of stars in phase space, metallicity, and age. The parameters of the EDFs are fitted using stars observed in the Sloan Extension for Galactic Understanding and Exploration-II (SEGUE-II) survey that trace the phase-space kinematics and chemistry out to ~70 kpc. A maximum a posteriori probability (MAP) estimate method and a Markov Chain Monte Carlo method are applied, taking into account the selection function in positions, distance, and metallicity for the survey. The best-fit EDF declines with actions less steeply at actions characteristic of the inner halo than at the larger actions characteristic of the outer halo, and older ages are found at smaller actions than at larger actions. In real space, the radial density profile steepens smoothly from -2 at ~2 kpc to -4 in the outer halo, with an axis ratio ~0.7 throughout. There is no indication for rotation in the BHBs, although this is highly uncertain. A moderate level of radial anisotropy is detected, with $\beta_s$ varying from isotropic to between ~0.1 and ~0.3 in the outer halo depending on latitude. The BHB data are consistent with an age gradient of -0.03 Gyr kpc$^{-1}$, with some uncertainty in the distribution of the larger ages. These results are consistent with a scenario in which older, larger systems contribute to the inner halo, whilst the outer halo is primarily comprised of younger, smaller systems.
  • We fit an Extended Distribution Function (EDF) to K giants in the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey. These stars are detected to radii ~80 kpc and span a wide range in [Fe/H]. Our EDF, which depends on [Fe/H] in addition to actions, encodes the entanglement of metallicity with dynamics within the Galaxy's stellar halo. Our maximum-likelihood fit of the EDF to the data allows us to model the survey's selection function. The density profile of the K giants steepens with radius from a slope ~-2 to ~-4 at large radii. The halo's axis ratio increases with radius from 0.7 to almost unity. The metal-rich stars are more tightly confined in action space than the metal-poor stars and form a more flattened structure. A weak metallicity gradient ~-0.001 dex/kpc, a small gradient in the dispersion in [Fe/H] of ~0.001 dex/kpc, and a higher degree of radial anistropy in metal-richer stars result. Lognormal components with peaks at ~-1.5 and ~-2.3 are required to capture the overall metallicity distribution, suggestive of the existence of two populations of K giants. The spherical anisotropy parameter varies between 0.3 in the inner halo to isotropic in the outer halo. If the Sagittarius stream is included, a very similar model is found but with a stronger degree of radial anisotropy throughout.
  • We create dynamical models of the massive elliptical galaxy, NGC 4649, using the N-body made-to-measure code, NMAGIC, and kinematic constraints from long-slit and planetary nebula (PN) data. We explore a range of potentials based on previous determinations from X-ray observations and a dynamical model fitting globular cluster (GC) velocities and a stellar density profile. The X-ray mass distributions are similar in the central region but have varying outer slopes, while the GC mass profile is higher in the central region and on the upper end of the range further out. Our models cannot differentiate between the potentials in the central region, and therefore if non-thermal pressures or multi-phase components are present in the hot gas, they must be smaller than previously inferred. In the halo, we find that the PN velocities are sensitive tracers of the mass, preferring a less massive halo than that derived from the GC mass profile, but similar to one of the mass distributions derived from X-rays. Our results show that the GCs may form a dynamically distinct system, and that the properties of the hot gas derived from X-rays in the outer halo have considerable uncertainties that need to be better understood. Estimating the mass in stars using photometric information and a stellar population mass-to-light ratio, we infer a dark matter mass fraction in NGC 4649 of ~0.39 at 1Re (10.5 kpc) and ~0.78 at 4Re. We find that the stellar orbits are isotropic to mildly radial in the central ~6 kpc depending on the potential assumed. Further out, the orbital structure becomes slightly more radial along R and more isotropic along z, regardless of the potential assumed. In the equatorial plane, azimuthal velocity dispersions dominate over meridional velocity dispersions, implying that meridional velocity anisotropy is the mechanism for flattening the stellar system.
  • Using a counter-dispersed slitless spectroscopy technique, we detect and measure the line-of-sight velocities of 187 planetary nebulae (PNe) around one of the nearest cD galaxies, NGC 1399, with FORS1 on the VLT. We describe the method for identifying and classifying the emission-line sources and the procedure for computing their J2000 coordinates and velocities. The number of PN detections and the errors in the velocity measurements (37 km/s indicate that this technique is comparable to other methods, such as that described by Teodorescu et al. (2005). We present the spatial distribution of the PNe and a basic analysis of their velocities. The PN two-dimensional velocity field shows marginal rotation consistent with other studies. We also find a low-velocity substructure in the halo and a flatter velocity-dispersion profile compared to previous observations that extends to ~400 arcsec. The detection of a low-velocity subcomponent underscores the importance of discrete velocity tracers for the detection of un-mixed components. The new velocity-dispersion profile is in good agreement with revised velocity dispersions for the red globular clusters in NGC 1399, using the data of Schuberth et al. (2009). The outer parts of this profile are consistent with one of the dynamical models of Kronawitter et al. (2000), which corresponds to a circular velocity of ~340 km/s and a rescaled B-band mass-to-light ratio of ~20 at 7' radius. These measurements trace the kinematics of the outer halo and disentangle the heterogenous populations in the Fornax Cluster core. The new data set the stage for a revised dynamical model of the outer halo of NGC 1399.
  • We use a new non-parametric Bayesian approach to obtain the most probable mass distributions and circular velocity curves along with their confidence ranges, given deprojected density and temperature profiles of the hot gas surrounding X-ray bright elliptical galaxies. For a sample of six X-ray bright ellipticals, we find that all circular velocity curves are rising in the outer parts due to a combination of a rising temperature profile and a logarithmic pressure gradient that increases in magnitude. Comparing the circular velocity curves we obtain from X-rays to those obtained from dynamical models, we find that the former are often lower in the central ~10 kpc. This is probably due to a combination of: i) Non-thermal contributions of up to ~35% in the pressure (with stronger effects in NGC 4486), ii) multiple-temperature components in the hot gas, iii) incomplete kinematic spatial coverage in the dynamical models, and iv) mass profiles that are insufficiently general in the dynamical modelling. Complementing the total mass information from the X-rays with photometry and stellar population models to infer the dark matter content, we find evidence for massive dark matter haloes with dark matter mass fractions of ~35-80% at 2Re, rising to a maximum of 80-90% at the outermost radii. We also find that the six galaxies follow a Tully-Fisher relation with slope ~4 and that their circular velocities at 1Re correlate strongly with the velocity dispersion of the local environment. As a result, the galaxy luminosity at 1Re also correlates with the velocity dispersion of the environment. These relations suggest a close link between the properties of central X-ray bright elliptical galaxies and their environments (abridged).
  • We present high resolution FLAMES/VLT spectroscopy of intracluster planetary nebula (ICPN) candidates, targeting three new fields in the Virgo cluster core with surface brightness down to mu_B = 28.5. Based on the projected phase space information we separate the old and 12 newly-confirmed PNs into galaxy and intracluster components. The M87 PNs are confined to the extended stellar envelope of M87, within a projected radius of ~ 160 kpc, while the ICPNs are scattered across the whole surveyed region between M87 and M86. The velocity dispersions determined from the M87 PNs at projected radii of 60 kpc and 144 kpc show that the galaxy's velocity dispersion profile decreases in the outer halo, down to 78 +/- 25 km/s. A Jeans model for the M87 halo stars in the gravitational potential traced by the X-ray emission fits the observed velocity dispersion profile only if the stellar orbits are strongly radially anisotropic (beta ~= 0.4 at r ~= 10 kpc increasing to 0.8 at the outer edge), and if additionally the stellar halo is truncated at ~= 150 kpc average elliptical radius. From the spatial and velocity distribution of the ICPNs we infer that M87 and M86 are falling towards each other and that we may be observing them just before the first close pass. The inferred luminosity-specific PN numbers for the M87 halo and the ICL are in the range of values observed for old (> 10 Gyr) stellar populations (abridged).