Search

Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gravity, cosmology, and galaxy-formation code SWIFT. The software package exploits hybrid shared- and distributed-memory task-based parallelism, asynchronous communications, and domain-decomposition algorithms based on balancing the workload, rather than the data, to efficiently exploit modern high-performance computing cluster architectures. Gravity is solved for using a fast-multipole-method, optionally coupled to a particle mesh solver in Fourier space to handle periodic volumes. For gas evolution, multiple modern flavours of Smoothed Particle Hydrodynamics are implemented. SWIFT also evolves neutrinos using a state-of-the-art particle-based method. Two complementary networks of sub-grid models for galaxy formation as well as extensions to simulate planetary physics are also released as part of the code. An extensive set of output options, including snapshots, light-cones, power spectra, and a coupling to structure finders are also included. We describe the overall code architecture, summarise the consistency and accuracy tests that were performed, and demonstrate the excellent weak-scaling performance of the code using a representative cosmological hydrodynamical problem with $\approx$$300$ billion particles. The code is released to the community alongside extensive documentation for both users and developers, a large selection of example test problems, and a suite of tools to aid in the analysis of large simulations run with SWIFT., Comment: 43 pages, 20 figures, accepted for publication in MNRAS. Code, documentation, and examples available at www.swiftsim.com

Cosmological simulations require the use of a multiple time-stepping scheme. Without such a scheme, cosmological simulations would be impossible due to their high level of dynamic range; over eleven orders of magnitude in density. Such a large dynamic range leads to a range of over four orders of magnitude in time-step, which presents a significant load-balancing challenge. In this work, the extreme adaptivity that cosmological simulations present is tackled in three main ways through the use of the code SWIFT. First, an adaptive mesh is used to ensure that only the relevant particles are interacted in a given time-step. Second, task-based parallelism is used to ensure efficient load-balancing within a single node, using pthreads and SIMD vectorisation. Finally, a domain decomposition strategy is presented, using the graph domain decomposition library METIS, that bisects the work that must be performed by the simulation between nodes using MPI. These three strategies are shown to give SWIFT near-perfect weak-scaling characteristics, only losing 25% performance when scaling from 1 to 4096 cores on a representative problem, whilst being more than 30x faster than the de-facto standard Gadget-2 code.

In particle-based simulations, neighbour finding (i.e finding pairs of particles to interact within a given range) is the most time consuming part of the computation. One of the best such algorithms, which can be used for both Molecular Dynamics (MD) and Smoothed Particle Hydrodynamics (SPH) simulations, is the pseudo-Verlet list algorithm. This algorithm, however, does not vectorise trivially, and hence makes it difficult to exploit SIMD-parallel architectures. In this paper, we present several novel modifications as well as a vectorisation strategy for the algorithm which lead to overall speed-ups over the scalar version of the algorithm of 2.24x for the AVX instruction set (SIMD width of 8), 2.43x for AVX2, and 4.07x for AVX-512 (SIMD width of 16)., Comment: 10 pages, 3 figures. Proceedings of the ParCo 2017 conference, Bologna, Italy, September 12-15th, 2017

We present a new 3D map of interstellar dust reddening, covering three quarters of the sky (declinations greater than -30 degrees) out to a distance of several kiloparsecs. The map is based on high-quality stellar photometry of 800 million stars from Pan-STARRS 1 and 2MASS. We divide the sky into sightlines containing a few hundred stars each, and then infer stellar distances and types, along with the line-of-sight dust distribution. Our new map incorporates a more accurate average extinction law and an additional 1.5 years of Pan-STARRS 1 data, tracing dust to greater extinctions and at higher angular resolutions than our previous map. Out of the plane of the Galaxy, our map agrees well with 2D reddening maps derived from far-infrared dust emission. After accounting for a 15% difference in scale, we find a mean scatter of 10% between our map and the Planck far-infrared emission-based dust map, out to a depth of 0.8 mag in E(r-z), with the level of agreement varying over the sky. Our map can be downloaded at http://argonaut.skymaps.info, or by its DOI: 10.7910/DVN/LCYHJG., Comment: Submitted to MNRAS. 17 pages, 16 figures

We present a comprehensive and precise description of the Sagittarius (Sgr) stellar stream's 3D geometry as traced by its old stellar population. This analysis draws on the sample of ${\sim}44,000$ RR Lyrae (RRab) stars from the Pan-STARRS1 (PS1) 3$\pi$ survey (Hernitschek et al. 2016,Sesar et al. 2017b), which is ${\sim}80\%$ complete and ${\sim}90\%$ pure within 80~kpc, and extends to ${\gtrsim} 120$~kpc with a distance precision of ${\sim} 3\%$. A projection of RR Lyrae stars within $|\tilde{B}|_{\odot}<9^\circ$ of the Sgr stream's orbital plane reveals the morphology of both the leading and the trailing arms at very high contrast, across much of the sky. In particular, the map traces the stream near-contiguously through the distant apocenters. We fit a simple model for the mean distance and line-of-sight depth of the Sgr stream as a function of the orbital plane angle $\tilde{\Lambda}_{\odot}$, along with a power-law background-model for the field stars. This modeling results in estimates of the mean stream distance precise to ${\sim}1\%$ and it resolves the stream's line-of-sight depth. These improved geometric constraints can serve as new constraints for dynamical stream models.

RR Lyrae stars may be the best practical tracers of Galactic halo (sub-)structure and kinematics. The PanSTARRS1 (PS1) $3\pi$ survey offers multi-band, multi-epoch, precise photometry across much of the sky, but a robust identification of RR Lyrae stars in this data set poses a challenge, given PS1's sparse, asynchronous multi-band light curves ($\lesssim 12$ epochs in each of five bands, taken over a 4.5-year period). We present a novel template fitting technique that uses well-defined and physically motivated multi-band light curves of RR Lyrae stars, and demonstrate that we get accurate period estimates, precise to 2~sec in $>80\%$ of cases. We augment these light curve fits with other {\em features} from photometric time-series and provide them to progressively more detailed machine-learned classification models. From these models we are able to select the widest ($3/4$ of the sky) and deepest (reaching 120 kpc) sample of RR Lyrae stars to date. The PS1 sample of $\sim 45,000$ RRab stars is pure (90\%), and complete (80\% at 80 kpc) at high galactic latitudes. It also provides distances precise to 3\%, measured with newly derived period-luminosity relations for optical/near-infrared PS1 bands. With the addition of proper motions from {\em Gaia} and radial velocity measurements from multi-object spectroscopic surveys, we expect the PS1 sample of RR Lyrae stars to become the premier source for studying the structure, kinematics, and the gravitational potential of the Galactic halo. The techniques presented in this study should translate well to other sparse, multi-band data sets, such as those produced by the Dark Energy Survey and the upcoming Large Synoptic Survey Telescope Galactic plane sub-survey., Comment: 19 pages, 13 figures, 5 tables, accepted to AJ. The the PS1 catalog of RR Lyrae stars will become publicly available on Nov 1 2017. For collaborations on projects and earlier access to the catalog, please contact the first author

We present a panoramic map of the entire Milky Way halo north of dec~-30 degrees (~30,000 deg^2), constructed by applying the matched-filter technique to the Pan-STARRS1 3Pi Survey dataset. Using single-epoch photometry reaching to g~22, we are sensitive to stellar substructures with heliocentric distances between 3.5 and ~35 kpc. We recover almost all previously-reported streams in this volume and demonstrate that several of these are significantly more extended than earlier datasets have indicated. In addition, we also report five new candidate stellar streams. One of these features appears significantly broader and more luminous than the others and is likely the remnant of a dwarf galaxy. The other four streams are consistent with a globular cluster origin, and three of these are rather short in projection (<10 degrees), suggesting that streams like Ophiuchus may not be that rare. Finally, a significant number of more marginal substructures are also revealed by our analysis; many of these features can also be discerned in matched-filter maps produced by other authors from SDSS data, and hence they are very likely to be genuine. However, the extant 3Pi data is currently too shallow to determine their properties or produce convincing CMDs. The global view of the Milky Way provided by Pan-STARRS1 provides further evidence for the important role of both globular cluster disruption and dwarf galaxy accretion in building the Milky Way's stellar halo., Comment: 12 pages, 6 figures. MNRAS, in press. The maps in FITS format for the 26 distance slices are made available to the community at http://dx.doi.org/10.5281/zenodo.60518, while full sky colour maps in various projections are provided at http://www.roe.ac.uk/~ejb/streams.html

We present a new open-source cosmological code, called SWIFT, designed to solve the equations of hydrodynamics using a particle-based approach (Smooth Particle Hydrodynamics) on hybrid shared/distributed-memory architectures. SWIFT was designed from the bottom up to provide excellent strong scaling on both commodity clusters (Tier-2 systems) and Top100-supercomputers (Tier-0 systems), without relying on architecture-specific features or specialized accelerator hardware. This performance is due to three main computational approaches: (1) Task-based parallelism for shared-memory parallelism, which provides fine-grained load balancing and thus strong scaling on large numbers of cores. (2) Graph-based domain decomposition, which uses the task graph to decompose the simulation domain such that the work, as opposed to just the data, as is the case with most partitioning schemes, is equally distributed across all nodes. (3) Fully dynamic and asynchronous communication, in which communication is modelled as just another task in the task-based scheme, sending data whenever it is ready and deferring on tasks that rely on data from other nodes until it arrives. In order to use these approaches, the code had to be re-written from scratch, and the algorithms therein adapted to the task-based paradigm. As a result, we can show upwards of 60% parallel efficiency for moderate-sized problems when increasing the number of cores 512-fold, on both x86-based and Power8-based architectures., Comment: 9 pages, 7 figures. Code, scripts and examples available at http://icc.dur.ac.uk/swift/

Using the Pan-STARRS1 survey, we derive limiting magnitude, spatial completeness and density maps that we use to probe the three dimensional structure and estimate the stellar mass of the so-called Monoceros Ring. The Monoceros Ring is an enormous and complex stellar sub-structure in the outer Milky Way disk. It is most visible across the large Galactic Anticenter region, 120 < l < 240 degrees, -30 < b < +40 degrees. We estimate its stellar mass density profile along every line of sight in 2 X 2 degree pixels over the entire 30,000 square degree Pan-STARRS1 survey using the previously developed MATCH software. By parsing this distribution into a radially smooth component and the Monoceros Ring, we obtain its mass and distance from the Sun along each relevant line of sight. The Monoceros Ring is significantly closer to us in the South (6 kpc) than in the North (9 kpc). We also create 2D cross sections parallel to the Galactic plane that show 135 degrees of the Monoceros Ring in the South and 170 degrees of the Monoceros Ring in the North. We show that the Northern and Southern structures are also roughly concentric circles, suggesting that they may be a wave rippling from a common origin. Excluding the Galactic plane, we observe an excess stellar mass of 4 million solar masses across 120 < l < 240 degrees. If we interpolate across the Galactic plane, we estimate that this region contains 8 million solar masses. If we assume (somewhat boldly) that the Monoceros Ring is a set of two Galactocentric rings, its total stellar mass is 60 million solar masses. Finally, if we assume that it is a set of two circles centered at a point 4 kpc from the Galactic center in the anti-central direction, as our data suggests, we estimate its stellar mass to be 40 million solar masses., Comment: 24 pages, 23 figures, accepted by ApJ

We present a recalibration of the Sloan Digital Sky Survey (SDSS) photometry with new flat fields and zero points derived from Pan-STARRS1 (PS1). Using PSF photometry of 60 million stars with $16 < r < 20$, we derive a model of amplifier gain and flat-field corrections with per-run RMS residuals of 3 millimagnitudes (mmag) in $griz$ bands and 15 mmag in $u$ band. The new photometric zero points are adjusted to leave the median in the Galactic North unchanged for compatibility with previous SDSS work. We also identify transient non-photometric periods in SDSS ("contrails") based on photometric deviations co-temporal in SDSS bands. The recalibrated stellar PSF photometry of SDSS and PS1 has an RMS difference of {9,7,7,8} mmag in $griz$, respectively, when averaged over $15'$ regions., Comment: 13 pages, 7 figures, ApJ, in press. "Hypercalibration" refers to using repeat measurements of many stars from multiple surveys to constrain calibration parameters

We present the discovery of three new Milky Way satellites from our search for compact stellar overdensities in the photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS 1, or PS1) 3pi survey. The first satellite, Laevens 3, is located at a heliocentric distance of d=67+/-3 kpc. With a total magnitude of Mv=-4.4+/-0.3 and a half-light radius rh=7+/-2 pc, its properties resemble those of outer halo globular clusters. The second system, Draco II/Laevens 4 (Dra II), is a closer and fainter satellite (d~20 kpc, Mv =-2.9+/-0.8), whose uncertain size (rh = 19 +8/-6 pc) renders its classification difficult without kinematic information; it could either be a faint and extended globular cluster or a faint and compact dwarf galaxy. The third satellite, Sagittarius II/Laevens 5 (Sgr II), has an ambiguous nature as it is either the most compact dwarf galaxy or the most extended globular cluster in its luminosity range (rh = 37 +9/-8 pc and Mv=-5.2+/-0.4). At a heliocentric distance of 67+/-5 kpc, this satellite lies intriguingly close to the expected location of the trailing arm of the Sagittarius stellar stream behind the Sagittarius dwarf spheroidal galaxy (Sgr dSph). If confirmed through spectroscopic follow up, this connection would locate this part of the trailing arm of the Sagittarius stellar stream that has so far gone undetected. It would further suggest that Sgr II was brought into the Milky Way halo as a satellite of the Sgr dSph., Comment: Published in ApJ

We present a three-dimensional map of interstellar dust reddening, covering three-quarters of the sky out to a distance of several kiloparsecs, based on Pan-STARRS 1 and 2MASS photometry. The map reveals a wealth of detailed structure, from filaments to large cloud complexes. The map has a hybrid angular resolution, with most of the map at an angular resolution of 3.4' to 13.7', and a maximum distance resolution of ~25%. The three-dimensional distribution of dust is determined in a fully probabilistic framework, yielding the uncertainty in the reddening distribution along each line of sight, as well as stellar distances, reddenings and classifications for 800 million stars detected by Pan-STARRS 1. We demonstrate the consistency of our reddening estimates with those of two-dimensional emission-based maps of dust reddening. In particular, we find agreement with the Planck 353 GHz optical depth-based reddening map to within 0.05 mag in E(B-V) to a depth of 0.5 mag, and explore systematics at reddenings less than E(B-V) ~ 0.08 mag. We validate our per-star reddening estimates by comparison with reddening estimates for stars with both SDSS photometry and SEGUE spectral classifications, finding per-star agreement to within 0.1 mag out to a stellar E(B-V) of 1 mag. We compare our map to two existing three-dimensional dust maps, by Marshall et al. (2006) and Lallement et al. (2013), demonstrating our finer angular resolution, and better distance resolution compared to the former within ~3 kpc. The map can be queried or downloaded at http://argonaut.skymaps.info. We expect the three-dimensional reddening map presented here to find a wide range of uses, among them correcting for reddening and extinction for objects embedded in the plane of the Galaxy, studies of Galactic structure, calibration of future emission-based dust maps and determining distances to objects of known reddening., Comment: Accepted for publication in ApJ

We present the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest spectroscopic survey to explicitly target variable objects, avoiding pre-selection on the basis of colors or detailed modeling of specific variability characteristics. Kernel Density Estimate (KDE) analysis of our target population performed on SDSS Stripe 82 data suggests our target sample will be 95% pure (meaning 95% of objects we select have genuine luminosity variability of a few magnitudes or more). Our final spectroscopic sample will contain roughly 135,000 quasars and 85,000 stellar variables, approximately 4,000 of which will be RR Lyrae stars which may be used as outer Milky Way probes. The variability-selected quasar population has a smoother redshift distribution than a color-selected sample, and variability measurements similar to those we develop here may be used to make more uniform quasar samples in large surveys. The stellar variable targets are distributed fairly uniformly across color space, indicating that TDSS will obtain spectra for a wide variety of stellar variables including pulsating variables, stars with significant chromospheric activity, cataclysmic variables and eclipsing binaries. TDSS will serve as a pathfinder mission to identify and characterize the multitude of variable objects that will be detected photometrically in even larger variability surveys such as LSST., Comment: 22 pages, 18 figures, accepted by ApJ

The extensible N-Dimensional Data Format (NDF) was designed and developed in the late 1980s to provide a data model suitable for use in a variety of astronomy data processing applications supported by the UK Starlink Project. Starlink applications were used extensively, primarily in the UK astronomical community, and form the basis of a number of advanced data reduction pipelines today. This paper provides an overview of the historical drivers for the development of NDF and the lessons learned from using a defined hierarchical data model for many years in data reduction software, data pipelines and in data acquisition systems., Comment: 19 pages, 7 figures, submitted to the Astronomy & Computing special issue on astronomy data formats

We present the discovery of 61 wide (>5 arcsecond) separation, low-mass (stellar and substellar) companions to stars in the solar neighborhood identified from Pan-STARRS\,1 (PS1) data and the spectral classification of 27 previously known companions. Our companions represent a selective subsample of promising candidates and span a range in spectral type of K7-L9 with the addition of one DA white dwarf. These were identified primarily from a dedicated common proper motion search around nearby stars, along with a few as serendipitous discoveries from our Pan-STARRS1 brown dwarf search. Our discoveries include 24 new L dwarf companions and one known L dwarf not previously identified as a companion. The primary stars around which we searched for companions come from a list of bright stars with well-measured parallaxes and large proper motions from the Hipparcos catalog (8583 stars, mostly A-K~dwarfs) and fainter stars from other proper motion catalogues (79170 stars, mostly M~dwarfs). We examine the likelihood that our companions are chance alignments between unrelated stars and conclude that this is unlikely for the majority of the objects that we have followed-up spectroscopically. We also examine the entire population of ultracool (>M7) dwarf companions and conclude that while some are loosely bound, most are unlikely to be disrupted over the course of $\sim$10 Gyr. Our search increases the number of ultracool M dwarf companions wider than 300 AU by 88% and increases the number of L dwarf companions in the same separation range by 96%. Finally, we resolve our new L dwarf companion to HIP 6407 into a tight (0.13 arcsecond, 7.4 AU) L1+T3 binary, making the system a hierarchical triple. Our search for these key benchmarks against which brown dwarf and exoplanet atmosphere models are tested has yielded the largest number of discoveries to date., Comment: 74 pages, 17 figures, 13 tables, accepted to ApJ, updated with corrected version of Table 13

SPLAT-VO is a powerful graphical tool for displaying, comparing, modifying and analyzing astronomical spectra, as well as searching and retrieving spectra from services around the world using Virtual Observatory (VO) protocols and services. The development of SPLAT-VO started in 1999, as part of the Starlink StarJava initiative, sometime before that of the VO, so initial support for the VO was necessarily added once VO standards and services became available. Further developments were supported by the Joint Astronomy Centre, Hawaii until 2009. Since end of 2011 development of SPLAT-VO has been continued by the German Astrophysical Virtual Observatory, and the Astronomical Institute of the Academy of Sciences of the Czech Republic. From this time several new features have been added, including support for the latest VO protocols, along with new visualization and spectra storing capabilities. This paper presents the history of SPLAT-VO, it's capabilities, recent additions and future plans, as well as a discussion on the motivations and lessons learned up to now., Comment: 15 pages, 6 figures, accepted for publication in Astronomy & Computing

We use a WISE-2MASS-Pan-STARRS1 galaxy catalog to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. We obtain photometric redshifts using our multicolor data set to create a tomographic map of the galaxy distribution. The radial density profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic Microwave Background results, we test for underdensities within two angular radii, $5^\circ$, and $15^\circ$. Our data, combined with an earlier measurement by Granett et al 2010, are consistent with a large $R_{\rm void}=(192 \pm 15)h^{-1} Mpc $ $(2\sigma)$ supervoid with $\delta \simeq -0.13 \pm 0.03$ centered at $z=0.22\pm0.01$. Such a supervoid, constituting a $\sim3.5 \sigma$ fluctuation in the $\Lambda CDM$ model, is a plausible cause for the Cold Spot., Comment: 4 pages, 2 figures, Proceedings of IAU 306 Symposium: Statistical Challenges in 21st Century Cosmology

We present a panoptic view of the stellar structure in the Galactic disk's outer reaches commonly known as the Monoceros Ring, based on data from Pan-STARRS1. These observations clearly show the large extent of the stellar overdensities on both sides of the Galactic disk, extending between b = -25 and b = +35 degrees and covering over 130 degrees in Galactic longitude. The structure exhibits a complex morphology with both stream-like features and a sharp edge to the structure in both the north and the south. We compare this map to mock observations of two published simulations aimed at explaining such structures in the outer stellar disk, one postulating an origin as a tidal stream and the other demonstrating a scenario where the disk is strongly distorted by the accretion of a satellite. These morphological comparisons of simulations can link formation scenarios to observed structures, such as demonstrating that the distorted-disk model can produce thin density features resembling tidal streams. Although neither model produces perfect agreement with the observations--the tidal stream predicts material at larger distances which is not detected while in the distorted disk model the midplane is warped to an excessive degree--future tuning of the models to accommodate these latest data may yield better agreement., Comment: 10 pages, 7 figures. Accepted to ApJ. Key figures and data can be found at http://dept.astro.lsa.umich.edu/~ctslater/MRi/

Standard inflationary hot big bang cosmology predicts small fluctuations in the Cosmic Microwave Background (CMB) with isotropic Gaussian statistics. All measurements support the standard theory, except for a few anomalies discovered in the Wilkinson Microwave Anisotropy Probe maps and confirmed recently by the Planck satellite. The Cold Spot is one of the most significant of such anomalies, and the leading explanation of it posits a large void that imprints this extremely cold area via the linear Integrated Sachs-Wolfe (ISW) effect due to the decay of gravitational potentials over cosmic time, or via the Rees-Sciama (RS) effect due to late-time non-linear evolution. Despite several observational campaigns targeting the Cold Spot region, to date no suitably large void was found at higher redshifts $z > 0.3$. Here we report the detection of an $R =(192 \pm 15) h^{-1}Mpc$ size supervoid of depth $\delta = -0.13 \pm 0.03$, and centred at redshift $z = 0.22$. This supervoid, possibly the largest ever found, is large enough to significantly affect the CMB via the non-linear RS effect, as shown in our Lemaitre-Tolman-Bondi framework. This discovery presents the first plausible explanation for any of the physical CMB anomalies, and raises the possibility that local large-scale structure could be responsible for other anomalies as well., Comment: 8 pages, 4 figures, to appear in the proceedings of the Moriond Cosmology Conference 2014

We present the fiducial sequences of a sample of Galactic star clusters in the five bands of the Pan-STARRS1 (PS1) photometric system (g_P1, r_P1, i_P1, z_P1, y_P1). These empirical sequences -- which include the red giant and sub-giant branches, the main sequence, and the horizontal branch -- were defined from deep colour-magnitude diagrams reaching below the oldest main-sequence turn-offs of 13 globular and 3 old open clusters covering a wide range of metallicities (-2.4 < [Fe/H] < +0.4). We find excellent agreement for the nine clusters in common with previous studies in similar photometric systems when transformed to the PS1 system. Because the photometric and spectroscopic properties of these stellar populations are accurately known, the fiducials provide a solid basis for the interpretation of observations in the PS1 system, as well as valuable constraints to improve the empirical colour--$T_{eff}$ relations., Comment: 6 pages, 5 figures (+15 figures in Appendix). Re-submitted to MNRAS after addressing the referee's comments. Some figures degraded due to size limitations

We report the discovery of a thin stellar stream found in Pan-STARRS1 photometry near the Galactic bulge in the constellation of Ophiuchus. It appears as a coherent structure in the colour-selected stellar density maps produced to search for tidal debris around nearby globular clusters. The stream is exceptionally short and narrow; it is about 2.5{\deg} long and 6' wide in projection. The colour-magnitude diagram of this object, which harbours a blue horizontal-branch, is consistent with an old and relatively metal-poor population ([Fe/H]~-1.3) located 9.5 +/- 0.9 kpc away at (l,b) ~ (5{\deg},+32{\deg}), and 5.0 +/- 1.0 kpc from the Galactic centre. These properties argue for a globular cluster as progenitor. The finding of such a prominent, nearby stream suggests that many streams could await discovery in the more densely populated regions of our Galaxy., Comment: 5 pages, 3 figures. MNRAS, in press

We use the WISE-2MASS infrared galaxy catalog matched with Pan-STARRS1 (PS1) galaxies to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. Our imaging catalog has median redshift $z\simeq 0.14$, and we obtain photometric redshifts from PS1 optical colours to create a tomographic map of the galaxy distribution. The radial profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic Microwave Background results, we test for underdensities within two angular radii, $5^\circ$, and $15^\circ$. The counts in photometric redshift bins show significantly low densities at high detection significance, $\gtrsim 5 \sigma$ and $\gtrsim 6 \sigma$, respectively, for the two fiducial radii. The line-of-sight position of the deepest region of the void is $z\simeq 0.15-0.25$. Our data, combined with an earlier measurement by Granett et al. 2010, are consistent with a large $R_{\rm void}=(220 \pm 50) h^{-1}Mpc $ supervoid with $\delta_{m} \simeq -0.14 \pm 0.04$ centered at $z=0.22\pm0.03$. Such a supervoid, constituting at least a $\simeq 3.3\sigma$ fluctuation in a Gaussian distribution of the $\Lambda CDM$ model, is a plausible cause for the Cold Spot., Comment: 8 pages, 6 figures, accepted for publication by MNRAS