M. C. Smith

THE MILKY WAY'S CIRCULAR VELOCITY CURVE TO 60 kpc AND AN ESTIMATE OF THE DARK MATTER HALO MASS FROM THE KINEMATICS OF ∼2400 SDSS BLUE HORIZONTAL-BRANCH STARS

We derive new constraints on the mass of the Milky Ways dark matter halo, based on 2401 rigorously selected blue horizontal-branch halo stars from SDSS DR6. This sample enables construction of the full line-of-sight velocity distribution at different galactocentric radii. To interpret these distributions, we compare them to matched mock observations drawn from two different cosmological galaxy formation simulations designed to resemble the Milky Way. This procedure results in an estimate of the Milky Ways circular velocity curve to ~60 kpc, which is found to be slightly falling from the adopted value of 220 km s?1 at the Suns location, and implies -->M( Vcir(r) , derived in statistically independent bins, is found to be consistent with the expectations from an NFW dark matter halo with the established stellar mass components at its center. If we assume that an NFW halo profile of characteristic concentration holds, we can use the observations to estimate the virial mass of the Milky Ways dark matter halo, -->Mvir = 1.0+ 0.3?0.2 ? 1012 M?, which is lower than many previous estimates. We have checked that the particulars of the cosmological simulations are unlikely to introduce systematics larger than the statistical uncertainties. This estimate implies that nearly 40% of the baryons within the virial radius of the Milky Ways dark matter halo reside in the stellar components of our Galaxy. A value for -->Mvir of only ~ -->1 ? 1012 M? also (re)opens the question of whether all of the Milky Ways satellite galaxies are on bound orbits.

Substructure revealed by RR Lyraes in SDSS Stripe 82

We present an analysis of the substructure revealed by RR Lyraes in Sloan Digital Sky Survey Stripe 82, which covers 2°5 in declination on the celestial equator over the right ascension range α = 20 h .7 to 3 h .3. We use the new public archive of light-motion curves in Stripe 82, published by Bramich et al. in 2008, to identify a set of high-quality RR Lyrae candidates. Period estimates are determined to high accuracy using a string-length method. A subset of 178 RR Lyraes with spectrally derived metallicities are employed to derive metallicity-period-amplitude relations, which are then used, together with archive magnitude data and light-curve Fourier decomposition, to estimate metallicities and hence distances for the entire sample. The RR Lyraes lie 5-115 kpc from the Galactic Centre, with distance estimates accurate to ∼8 per cent. The RR Lyraes are further divided into subsets of 316 RRab types and 91 RRc types based on their period, colour and metallicity. We fit a smooth density law to the distribution as a simple representation of the data. For Galactocentric radii 5-25 kpc the number density of RR Lyraes falls as r ―2,4 , but beyond 25 kpc, the number density falls much more steeply, as r ―4.5 . However, we stress that in practice the density distribution is not smooth, but dominated by clumps and substructure. Samples of 55 and 237 RR Lyraes associated with the Sagittarius Stream and the Hercules-Aquila Cloud, respectively, are identified. Hence, ∼70 per cent of the RR Lyraes in Stripe 82 belong to known substructure, and the sharp break in the density law reflects the fact that the dominant substructure in Stripe 82 - the Hercules-Aquila Cloud and the Sagittarius Stream -lie within 40 kpc. In fact, almost 60 per cent of all the RR Lyraes in Stripe 82 are associated with the Hercules-Aquila Cloud alone, which emphasizes the clouds pre-eminence. Additionally, evidence of a new and distant substructure - the Pisces Overdensity ― is found, consisting of 28 faint RRLyraes centred on Galactic coordinates (l ≈ 80°, b ≈ -55°), with distances of ∼80 kpc. The total stellar mass in the Pisces Overdensity is ∼10 4 M ⊙ and its metallicity is [Fe/H]∼ ―1.5.

Light and motion in SDSS Stripe 82 : the catalogues

We present a new public archive of light-motion curves in Sloan Digital Sky Survey (SDSS) Stripe 82, covering 99° in right ascension from a = 20 h 7 to 3 h .3 and spanning 2.°52 in declination from δ = -1°.26 to 1°.26, for a total sky area of ∼249 deg 2 . Stripe 82 has been repeatedly monitored in the u, g, r, i and z bands over a seven-year baseline. Objects are cross-matched between runs, taking into account the effects of any proper motion. The resulting catalogue contains almost 4 million light-motion curves of stellar objects and galaxies. The photometry are recalibrated to correct for varying photometric zero-points, achieving ∼20 and 30 mmag rms accuracy down to 18 mag in the g, r, i and z bands for point sources and extended sources, respectively. The astrometry are recalibrated to correct for inherent systematic errors in the SDSS astrometric solutions, achieving ∼32 and 35 mas rms accuracy down to 18 mag for point sources and extended sources, respectively. For each light-motion curve, 229 photometric and astrometric quantities are derived and stored in a higher level catalogue. On the photometric side, these include mean exponential and point spread function (PSF) magnitudes along with uncertainties, rms scatter, Χ 2 per degree of freedom, various magnitude distribution percentiles, object type (stellar or galaxy), and eclipse, Stetson and Vidrih variability indices. On the astrometric side, these quantities include mean positions, proper motions as well as their uncertainties and Χ 2 per degree of freedom. The light-motion curve catalogue presented here is complete down to r ∼ 21.5 and is at present the deepest large-area photometric and astrometric variability catalogue available.

Kinematics of SDSS subdwarfs: Structure and substructure of the Milky Way halo

We construct a new sample of ∼ 1700 solar neighbourhood halo subdwarfs from the Sloan Digital Sky Survey (SDSS), selected using a reduced proper-motion diagram. Radial velocities come from the SDSS spectra and proper motions from the light-motion curve catalogue of Bramich et al. Using a photometric parallax relation to estimate distances gives us the full phase-space coordinates. Typical velocity errors are in the range 30-50 km s ―1 . This halo sample is one of the largest constructed to date and the disc contamination is at a level of ≲1 per cent. This enables us to calculate the halo velocity dispersion to excellent accuracy. We find that the velocity dispersion tensor is aligned in spherical polar coordinates and that (σ r ,σ φ , σ θ ) = (143 ± 2, 82 ± 2, 77 ± 2) km s ―1 . The stellar halo exhibits no net rotation, although the distribution of υ φ shows tentative evidence for asymmetry. The kinematics are consistent with a mildly flattened stellar density falling with distance like r ―3.75 . Using the full phase-space coordinates, we look for signs of kinematic substructure in the stellar halo. We find evidence for four discrete overdensities localized in angular momentum and suggest that they may be possible accretion remnants. The most prominent is the solar neighbourhood stream previously identified by Helmi et al., but the remaining three are new. One of these overdensities is potentially associated with a group of four globular clusters (NGC 5466, NGC 6934, M2 and M13) and raises the possibility that these could have been accreted as part of a much larger progenitor.

The OGLE view of microlensing towards the Magellanic Clouds – I. A trickle of events in the OGLE-II LMC data

We present the results from the Optical Gravitational Lensing Experiment II (OGLE-II) survey (1996–2000) towards the Large Magellanic Cloud (LMC), which has the aim of detecting the microlensing phenomena caused by dark matter compact objects in the Galactic halo [massive compact halo objects (MACHOs)]. n n n nWe use high-resolution Hubble Space Telescope images of the OGLE fields and derive the correction for the number of monitored stars in each field. This also yields blending distributions which we use in ‘catalogue-level’ Monte Carlo simulations of the microlensing events in order to calculate the detection efficiency of the events. n n n nWe detect two candidates for microlensing events in the All Stars Sample, which translates into an optical depth of 0.43 ± 0.33 × 10−7. If both events were due to MACHO, the fraction of mass of compact dark matter objects in the Galactic halo would be 8 ± 6 per cent. This optical depth, however, along with the characteristics of the events seems to be consistent with the self-lensing scenario, i.e. self-lensing alone is sufficient to explain the observed microlensing signal. Our results indicate the non-detection of MACHOs lensing towards the LMC with an upper limit on their abundance in the Galactic halo of 19 per cent for M= 0.4xa0M⊙ and 10 per cent for masses between 0.01 and 0.2xa0M⊙.

New ultracool and halo white dwarf candidates in SDSS Stripe 82

Abstract : A 2.5 100 region along the celestial equator (Stripe 82) has been imaged repeatedly from 1998 to 2005 by the Sloan Digital Sky Survey (SDSS). A new catalogue of ~4 million light-motion curves, together with over 200 derived statistical quantities, for objects in Stripe 82 brighter than r~21.5 has been constructed by combining these data by Bramich, et al. This catalogue is at present the deepest catalogue of its kind. Extracting ~130,000 objects with highest signal-to-noise ratio proper motions, we build a reduced proper motion diagram to illustrate the scientific promise of the catalogue. In this diagram, disc and halo sub-dwarfs are well-separated from the cool white dwarf sequence. Our sample of 1049 cool white dwarf candidates includes at least eight and possibly 21 new ultracool H-rich white dwarfs (Teff < 4000 K) and one new ultracool He-rich white dwarf candidate identified from their SDSS optical and UKIDSS infrared photometry. At least 10 new halo white dwarfs are also identified from their kinematics.

The OGLE View of Microlensing towards the Magellanic Clouds. II. OGLE-II SMC data

The primary goal of this paper is to provide the evidence that can either prove or falsify the hypothesis that dark matter in the Galactic halo can clump into stellar-mass compact objects. If such objects existed, they would act as lenses to external sources in the Magellanic Clouds, giving rise to an observable effect of microlensing. We present the results of our search for such events, based on the data from the second phase of the OGLE survey (1996-2000) towards the SMC. The data set we used is comprised of 2.1 million monitored sources distributed over an area of 2.4 square degrees. We found only one microlensing event candidate, however its poor quality light curve limited our discussion on the exact distance to the lensing object. nGiven a single event, taking the blending (crowding of stars) into account for the detection efficiency simulations, and deriving the HST-corrected number of monitored stars, the microlensing optical depth is tau=(1.55+-1.55)10e-7. This result is consistent with the expected SMC self-lensing signal, with no need of introducing dark matter microlenses. Rejecting the unconvincing event leads to the upper limit on the fraction of dark matter in the form of MACHOs to f<20 per cent for deflectors masses around 0.4 Msun and f<11 per cent for masses between 0.003 and 0.2 Msun (95 per cent confidence limit). Our result indicates that the Milky Ways dark matter is unlikely to be clumpy and form compact objects in the sub-solar-mass range.

The OGLE view of microlensing towards the Magellanic Clouds – II. OGLE-II Small Magellanic Cloud data

The primary goal of this paper is to provide the evidence that can either prove or falsify the hypothesis that dark matter in the Galactic halo can clump into stellar-mass compact objects. If such objects existed, they would act as lenses to external sources in the Magellanic Clouds, giving rise to an observable effect of microlensing. We present the results of our search for such events, based on the data from the second phase of the OGLE survey (1996-2000) towards the SMC. The data set we used is comprised of 2.1 million monitored sources distributed over an area of 2.4 square degrees. We found only one microlensing event candidate, however its poor quality light curve limited our discussion on the exact distance to the lensing object. nGiven a single event, taking the blending (crowding of stars) into account for the detection efficiency simulations, and deriving the HST-corrected number of monitored stars, the microlensing optical depth is tau=(1.55+-1.55)10e-7. This result is consistent with the expected SMC self-lensing signal, with no need of introducing dark matter microlenses. Rejecting the unconvincing event leads to the upper limit on the fraction of dark matter in the form of MACHOs to f<20 per cent for deflectors masses around 0.4 Msun and f<11 per cent for masses between 0.003 and 0.2 Msun (95 per cent confidence limit). Our result indicates that the Milky Ways dark matter is unlikely to be clumpy and form compact objects in the sub-solar-mass range.

A strip search for new very wide halo binaries

We report on a search for new wide halo binary stars in Sloan Digital Sky Survey (SDSS) Stripe 82. A list of new halo wide binary candidates which satisfy common proper motion and photometric constraints is provided. The projected separations of the sample lie between 0.007 and 0.25 pc. Although the sample is not large enough to improve constraints on dark matter in the halo, we find the wide binary angular separation function is broadly consistent with past work. We discuss the significance of the new sample for a number of astrophysical applications, including as a testbed for ideas about wide binary formation. For the subset of candidates which have radial velocity information, we make use of integrals of motion to investigate one such scheme in which the origin of Galactic wide binaries is associated with the accretion/disruption of stellar systems in the Galaxy. Additional spectroscopic observations of these candidate binaries will strengthen their usefulness in many of these respects. Based on our search experience in Stripe 82 we estimate that the upcoming Pan-STARRS survey will increase the sample size of wide halo binaries by over an order of magnitude.

Supernova feedback in numerical simulations of galaxy formation: separating physics from numerics

While feedback from massive stars exploding as supernovae (SNe) is thought to be one of the key ingredients regulating galaxy formation, theoretically it is still unclear how the available energy couples to the interstellar medium and how galactic scale outflows are launched. We present a novel implementation of six sub-grid SN feedback schemes in the moving-mesh code Arepo, including injections of thermal and/or kinetic energy, two parametrizations of delayed cooling feedback and a `mechanical feedback scheme that injects the correct amount of momentum depending on the relevant scale of the SN remnant resolved. All schemes make use of individually time-resolved SN events. Adopting isolated disk galaxy setups at different resolutions, with the highest resolution runs reasonably resolving the Sedov-Taylor phase of the SN, we aim to find a physically motivated scheme with as few tunable parameters as possible. As expected, simple injections of energy overcool at all but the highest resolution. Our delayed cooling schemes result in overstrong feedback, destroying the disk. The mechanical feedback scheme is efficient at suppressing star formation, agrees well with the Kennicutt-Schmidt relation and leads to converged star formation rates and galaxy morphologies with increasing resolution without fine tuning any parameters. However, we find it difficult to produce outflows with high enough mass loading factors at all but the highest resolution, indicating either that we have oversimplified the evolution of unresolved SN remnants, require other stellar feedback processes to be included, require a better star formation prescription or most likely some combination of these issues.