Source author record

Mario Juric

Mario Juric 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

20works
11topics
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

20 published item(s)

preprint2022arXiv

Simulating the Legacy Survey of Space and Time stellar content with TRILEGAL

We describe a large simulation of the stars to be observed by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The simulation is based on the TRILEGAL code, which resorts to large databases of stellar evolutionary tracks, synthetic spectra, and pulsation models, added to simple prescriptions for the stellar density and star formation histories of the main structures of the Galaxy, to generate mock stellar samples through a population synthesis approach. The main bodies of the Magellanic Clouds are also included. A complete simulation is provided for single stars, down to the $r=27.5$ mag depth of the co-added wide-fast-deep survey images. A second simulation is provided for a fraction of the binaries, including the interacting ones, as derived with the BinaPSE module of TRILEGAL. We illustrate the main properties and numbers derived from these simulations, including: comparisons with real star counts; the expected numbers of Cepheids, long-period variables and eclipsing binaries; the crowding limits as a function of seeing and filter; the star-to-galaxy ratios, etc. Complete catalogs are accessible through the NOIRLab Astro Data Lab, while the stellar density maps are incorporated in the LSST metrics analysis framework (MAF).

preprint2021arXiv

Checkpoint, Restore, and Live Migration for Science Platforms

We demonstrate a fully functional implementation of (per-user) checkpoint, restore, and live migration capabilities for JupyterHub platforms. Checkpointing -- the ability to freeze and suspend to disk the running state (contents of memory, registers, open files, etc.) of a set of processes -- enables the system to snapshot a user's Jupyter session to permanent storage. The restore functionality brings a checkpointed session back to a running state, to continue where it left off at a later time and potentially on a different machine. Finally, live migration enables moving running Jupyter notebook servers between different machines, transparent to the analysis code and w/o disconnecting the user. Our implementation of these capabilities works at the system level, with few limitations, and typical checkpoint/restore times of O(10s) with a pathway to O(1s) live migrations. It opens a myriad of interesting use cases, especially for cloud-based deployments: from checkpointing idle sessions w/o interruption of the user's work (achieving cost reductions of 4x or more), execution on spot instances w. transparent migration on eviction (with additional cost reductions up to 3x), to automated migration of workloads to ideally suited instances (e.g. moving an analysis to a machine with more or less RAM or cores based on observed resource utilization). The capabilities we demonstrate can make science platforms fully elastic while retaining excellent user experience.

preprint2021arXiv

Galactic Mass Estimates using Dwarf Galaxies as Kinematic Tracers

New mass estimates and cumulative mass profiles with Bayesian credible regions (c.r.) for the Milky Way (MW) are found using the Galactic Mass Estimator (GME) code and dwarf galaxy (DG) kinematic data from multiple sources. GME takes a hierarchical Bayesian approach to simultaneously estimate the true positions and velocities of the DGs, their velocity anisotropy, and the model parameters for the Galaxy's total gravitational potential. In this study, we incorporate meaningful prior information from past studies and simulations. The prior distributions for the physical model are informed by the results of Eadie & Juric (2019), which used globular clusters instead of DGs, as well as by the subhalo distributions of the Ananke Gaia-like surveys from Feedback In Realistic Environments-2 (Fire-2) cosmological simulations (see Sanderson et al. 2020). Using DGs beyond 45 kpc, we report median and 95% c.r estimates for $r_{200}$ = 212.8 (191.12,238.44) kpc, and for the total enclosed mass $M_{200}$ = 1.19 (0.87,1.68)$\times10^{12}M_{\odot}$ (adopting $Δ_c=200$). Median mass estimates at specific radii are also reported (e.g., $M(<50\text{ kpc})=0.52\times10^{12}M_{\odot}$ and $M(100\text{ kpc})=0.78\times10^{12}M_{\odot}$). Estimates are comparable to other recent studies using GAIA DR2 and DGs, but notably different from the estimates of Eadie & Juric (2019). We perform a sensitivity analysis to investigate whether individual DGs and/or a more massive Large Magellanic Cloud (LMC) on the order of $10^{11}M_{\odot}$ may be affecting our mass estimates. We find possible supporting evidence for the idea that some DGs are affected by a massive LMC and are not in equilibrium with the MW.

preprint2020arXiv

Characterization of the Nucleus, Morphology and Activity of Interstellar Comet 2I/Borisov by Optical and Near-Infrared GROWTH, Apache Point, IRTF, ZTF and Keck Observations

We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium.

preprint2015arXiv

Asteroid Discovery and Characterization with the Large Synoptic Survey Telescope (LSST)

The Large Synoptic Survey Telescope (LSST) will be a ground-based, optical, all-sky, rapid cadence survey project with tremendous potential for discovering and characterizing asteroids. With LSST's large 6.5m diameter primary mirror, a wide 9.6 square degree field of view 3.2 Gigapixel camera, and rapid observational cadence, LSST will discover more than 5 million asteroids over its ten year survey lifetime. With a single visit limiting magnitude of 24.5 in r-band, LSST will be able to detect asteroids in the Main Belt down to sub-kilometer sizes. The current strawman for the LSST survey strategy is to obtain two visits (each visit being a pair of back-to-back 15s exposures) per field, separated by about 30 minutes, covering the entire visible sky every 3-4 days throughout the observing season, for ten years. The catalogs generated by LSST will increase the known number of small bodies in the Solar System by a factor of 10-100 times, among all populations. The median number of observations for Main Belt asteroids will be on the order of 200-300, with Near Earth Objects receiving a median of 90 observations. These observations will be spread among ugrizy bandpasses, providing photometric colors and allowing sparse lightcurve inversion to determine rotation periods, spin axes, and shape information. These catalogs will be created using automated detection software, the LSST Moving Object Processing System (MOPS), that will take advantage of the carefully characterized LSST optical system, cosmetically clean camera, and recent improvements in difference imaging. Tests with the prototype MOPS software indicate that linking detections (and thus discovery) will be possible at LSST depths with our working model for the survey strategy, but evaluation of MOPS and improvements in the survey strategy will continue. All data products and software created by LSST will be publicly available.

preprint2015arXiv

Hypercalibration: A Pan-STARRS1-based recalibration of the Sloan Digital Sky Survey

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.

preprint2014arXiv

Radio Astronomy in LSST Era

A community meeting on the topic of "Radio Astronomy in the LSST Era" was hosted by the National Radio Astronomy Observatory in Charlottesville, VA (2013 May 6--8). The focus of the workshop was on time domain radio astronomy and sky surveys. For the time domain, the extent to which radio and visible wavelength observations are required to understand several classes of transients was stressed, but there are also classes of radio transients for which no visible wavelength counterpart is yet known, providing an opportunity for discovery. From the LSST perspective, the LSST is expected to generate as many as 1 million alerts nightly, which will require even more selective specification and identification of the classes and characteristics of transients that can warrant follow up, at radio or any wavelength. The LSST will also conduct a deep survey of the sky, producing a catalog expected to contain over 38 billion objects in it. Deep radio wavelength sky surveys will also be conducted on a comparable time scale, and radio and visible wavelength observations are part of the multi-wavelength approach needed to classify and understand these objects. Radio wavelengths are valuable because they are unaffected by dust obscuration and, for galaxies, contain contributions both from star formation and from active galactic nuclei. The workshop touched on several other topics, on which there was consensus including the placement of other LSST "Deep Drilling Fields," inter-operability of software tools, and the challenge of filtering and exploiting the LSST data stream. There were also topics for which there was insufficient time for full discussion or for which no consensus was reached, which included the procedures for following up on LSST observations and the nature for future support of researchers desiring to use LSST data products.

preprint2014arXiv

The Milky Way Tomography with SDSS. V. Mapping the Dark Matter Halo

We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding 10 kpc and galactocentric distances exceeding 20 kpc. Our analysis utilizes Jeans equations to generate two-dimensional acceleration maps throughout the volume; this approach is thoroughly tested on a cosmologically derived N-body+SPH simulation of a MW-like galaxy. We show that the known accelerations (gradients of the gravitational potential) can be successfully recovered in such a realistic system. Leveraging the baryonic gravitational potential derived by Bovy & Rix (2013), we show that the gravitational potential implied by the SDSS observations cannot be explained, assuming Newtonian gravity, by visible matter alone: the gravitational force experienced by stars at galactocentric distances of 20 kpc is as much as three times stronger than what can be attributed to purely visible matter. We also show that the SDSS data provide a strong constraint on the shape of the dark matter halo potential. Within galactocentric distances of 20 kpc, the dark matter halo potential is well described as an oblate halo with axis ratio qDM=0.7+/-0.1; this corresponds to an axis ratio qDM=0.4+/-0.1 for the dark matter density distribution. Because of our precise two-dimensional measurements of the acceleration of the halo stars, we can reject several MOND models as an explanation of the observed behavior.

preprint2013arXiv

Galactic Stellar Populations in the Era of SDSS and Other Large Surveys

Studies of stellar populations, understood to mean collections of stars with common spatial, kinematic, chemical, and/or age distributions, have been reinvigorated during the last decade by the advent of large-area sky surveys such as SDSS, 2MASS, RAVE, and others. We review recent analyses of these data that, together with theoretical and modeling advances, are revolutionizing our understanding of the nature of the Milky Way, and galaxy formation and evolution in general. The formation of galaxies like the Milky Way was long thought to be a steady process leading to a smooth distribution of stars. However, the abundance of substructure in the multi-dimensional space of various observables, such as position, kinematics, and metallicity, is by now proven beyond doubt, and demonstrates the importance of mergers in the growth of galaxies. Unlike smooth models that involve simple components, the new data reviewed here clearly show many irregular structures, such as the Sagittarius dwarf tidal stream and the Virgo and Pisces overdensities in the halo, and the Monoceros stream closer to the Galactic plane. These recent developments have made it clear that the Milky Way is a complex and dynamical structure, one that is still being shaped by the merging of neighboring smaller galaxies. We also briefly discuss the next generation of wide-field sky surveys, such as SkyMapper, Pan-STARRS, Gaia and LSST, which will improve measurement precision manyfold, and comprise billions of individual stars. The ultimate goal, development of a coherent and detailed story of the assembly and evolutionary history of the Milky Way and other large spirals like it, now appears well within reach.

preprint2012arXiv

Constraints on the Shape of the Milky Way Dark Matter Halo from Jeans Equations Applied to SDSS Data

We search for evidence of dark matter in the Milky Way by utilizing the stellar number density distribution and kinematics measured by the Sloan Digital Sky Survey (SDSS) to heliocentric distances exceeding ~10 kpc. We employ the cylindrically symmetric form of Jeans equations and focus on the morphology of the resulting acceleration maps, rather than the normalization of the total mass as done in previous, mostly local, studies. Jeans equations are first applied to a mock catalog based on a cosmologically derived N-body + SPH simulation, and the known acceleration (gradient of gravitational potential) is successfully recovered. The same simulation is also used to quantify the impact of dark matter on the total acceleration. We use Galfast, a code designed to quantitatively reproduce SDSS measurements and selection effects, to generate a synthetic stellar catalog. We apply Jeans equations to this catalog and produce two-dimensional maps of stellar acceleration. These maps reveal that in a Newtonian framework, the implied gravitational potential cannot be explained by visible matter alone. The acceleration experienced by stars at galactocentric distances of ~20 kpc is three times larger than what can be explained by purely visible matter. The application of an analytic method for estimating the dark matter halo axis ratio to SDSS data implies an oblate halo with q_DM = 0.47 +/- 0.14 within the same distance range. These techniques can be used to map the dark matter halo to much larger distances from the Galactic center using upcoming deep optical surveys, such as LSST.

preprint2012arXiv

Ensemble Properties of Comets in the Sloan Digital Sky Survey

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u,g,r,i,z) photometry is used to determine the comets' colors, sizes, surface brightness profiles, and rates of dust production in terms of the Afρ formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(< H) \propto 10(0.49\pm0.05)H for H < 18, with evidence of a much shallower fit N(< H) \propto 10(0.19\pm0.03)H for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 \pm 0.05 in g - r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets.

preprint2012arXiv

Swarm-NG: a CUDA Library for Parallel n-body Integrations with focus on Simulations of Planetary Systems

We present Swarm-NG, a C++ library for the efficient direct integration of many n-body systems using highly-parallel Graphics Processing Unit (GPU), such as NVIDIA's Tesla T10 and M2070 GPUs. While previous studies have demonstrated the benefit of GPUs for n-body simulations with thousands to millions of bodies, Swarm-NG focuses on many few-body systems, e.g., thousands of systems with 3...15 bodies each, as is typical for the study of planetary systems. Swarm-NG parallelizes the simulation, including both the numerical integration of the equations of motion and the evaluation of forces using NVIDIA's "Compute Unified Device Architecture" (CUDA) on the GPU. Swarm-NG includes optimized implementations of 4th order time-symmetrized Hermite integration and mixed variable symplectic integration, as well as several sample codes for other algorithms to illustrate how non-CUDA-savvy users may themselves introduce customized integrators into the Swarm-NG framework. To optimize performance, we analyze the effect of GPU-specific parameters on performance under double precision. Applications of Swarm-NG include studying the late stages of planet formation, testing the stability of planetary systems and evaluating the goodness-of-fit between many planetary system models and observations of extrasolar planet host stars (e.g., radial velocity, astrometry, transit timing). While Swarm-NG focuses on the parallel integration of many planetary systems,the underlying integrators could be applied to a wide variety of problems that require repeatedly integrating a set of ordinary differential equations many times using different initial conditions and/or parameter values.

preprint2012arXiv

The Metallicity of the Monoceros Stream

We present low-resolution MMT Hectospec spectroscopy of 594 candidate Monoceros stream member stars. Based on strong color-magnitude diagram overdensities, we targeted three fields within the stream's footprint, with 178 deg < l < 203 deg and -25 deg < b < 25 deg. By comparing the measured iron abundances with those expected from smooth Galactic components alone, we measure, for the first time, the spectroscopic metallicity distribution function for Monoceros. We find the stream to be chemically distinct from both the thick disk and halo, with [Fe/H] = -1, and do not detect a trend in the stream's metallicity with Galactic longitude. Passing from b = +25 deg to b = -25 deg the median Monoceros metallicity trends upward by 0.1 dex, though uncertainties in modeling sample contamination by the disk and halo make this a marginal detection. In each field, we find Monoceros to have an intrinsic [Fe/H] dispersion of 0.10-0.22 dex. From the CaII K line, we measure [Ca/Fe] for a subsample of metal poor program stars with -1.1 < [Fe/H] < -0.5. In two of three fields, we find calcium deficiencies qualitatively similar to previously reported [Ti/Fe] underabundances in Monoceros and the Sagittarius tidal stream. Further, using 90 spectra of thick disk stars in the Monoceros pointings with b ~ +/-25 deg, we detect a 0.22 dex north/south metallicity asymmetry coincident with known stellar density asymmetry at R_GC ~ 12 kpc and |Z| ~ 1.7 kpc. Our median Monoceros [Fe/H] = -1.0 and its relatively low dispersion naturally fit the expectation for an appropriately luminous M_V ~ -13 dwarf galaxy progenitor.

preprint2012arXiv

Update on the Nature of Virgo Overdensity

We use the Eighth Data Release of Sloan Digital Sky Survey (SDSS DR8) catalog with its additional sky coverage of the southern Galactic hemisphere, to measure the extent and study the nature of the Virgo Overdensity (VOD; Juric et al. 2008). The data show that the VOD extends over no less than 2000 deg^2, with its true extent likely closer to 3000 deg^2. We test whether the VOD can be attributed to a tilt in the stellar halo ellipsoid with respect to the plane of the Galactic disk and find that the observed symmetry of the north-south Galactic hemisphere star counts excludes this possibility. We argue that the Virgo Overdensity, in spite of its wide area and cloud-like appearance, is still best explained by a minor merger. Its appearance and position is qualitatively similar to a near perigalacticon merger event and, assuming that the VOD and the Virgo Stellar Stream share the same progenitor, consistent with the VSS orbit determined by Casetti-Dinescu et al. (2009).

preprint2011arXiv

SDSS, LSST, and Gaia: Lessons and Synergies

The advent of deep, wide, accurate, digital photometric surveys exemplified by the Sloan Digital Sky Survey (SDSS) has had a profound impact on studies of the Milky Way. In the past decade, we have transitioned from a scarcity to an (over)abundance of precise, well calibrated, observations of stars over a large fraction of the Galaxy. The avalanche of data will continue throughout this decade, culminating with Gaia and LSST. This new reality will necessitate changes in methodology, habits, and expectations both on the side of the large survey projects as well as the astrophysics community at large. We argue, based on the experience with SDSS, that surveys should release data as early and often as possible incorporating incremental improvements in each subsequent release, as opposed to holding off for a single, big, final release. The scientific community will need to reciprocate by performing analyses and (re-analyses) appropriate to the current fidelity of the released data, understanding that these are continually evolving and improving products.

preprint2011arXiv

The Shape and Profile of the Milky Way Halo as Seen by the CFHT Legacy Survey

We use Canada-France-Hawaii Telescope Legacy Survey data for 170 deg^2, recalibrated and transformed to the Sloan Digital Sky Survey ugri photometric system, to study the distribution of near-turnoff main-sequence stars in the Galactic halo along four lines of sight to heliocentric distances of ~35 kpc. We find that the halo stellar number density profile becomes steeper at Galactocentric distances greater than R_{gal}~28 kpc, with the power law index changing from n_{inner}=-2.62+-0.04 to n_{outer}=-3.8+-0.1. In particular, we test a series of single power law models and find them to be strongly disfavored by the data. The parameters for the best-fit Einasto profile are n=2.2+-0.2 and R_e=22.2+-0.4 kpc. We measure the oblateness of the halo to be q=c/a=0.70+-0.01 and detect no evidence of it changing across the range of probed distances. The Sagittarius stream is detected in the l=173 and b=-62 direction as an overdensity of [Fe/H] ~ -1.5 dex stars at R_{gal}~32 kpc, providing a new constraint for the Sagittarius stream and dark matter halo models. We also detect the Monoceros stream as an overdensity of [Fe/H] > -1.5 dex stars in the l=232 and b=26 direction at R_{gal}<25 kpc. In the two sightlines where we do not detect significant substructure, the median metallicity is found to be independent of distance within systematic uncertainties ([Fe/H] ~ -1.5+-0.1 dex).

preprint2009arXiv

Light Curve Templates and Galactic Distribution of RR Lyrae Stars from Sloan Digital Sky Survey Stripe 82

We present an improved analysis of halo substructure traced by RR Lyrae stars in the SDSS stripe 82 region. With the addition of SDSS-II data, a revised selection method based on new ugriz light curve templates results in a sample of 483 RR Lyrae stars that is essentially free of contamination. The main result from our first study persists: the spatial distribution of halo stars at galactocentric distances 5--100 kpc is highly inhomogeneous. At least 20% of halo stars within 30 kpc from the Galactic center can be statistically associated with substructure. We present strong direct evidence, based on both RR Lyrae stars and main sequence stars, that the halo stellar number density profile significantly steepens beyond a Galactocentric distance of ~30 kpc, and a larger fraction of the stars are associated with substructure. By using a novel method that simultaneously combines data for RR Lyrae and main sequence stars, and using photometric metallicity estimates for main sequence stars derived from deep co-added u-band data, we measure the metallicity of the Sagittarius dSph tidal stream (trailing arm) towards R.A.2h-3h and Dec~0 deg to be 0.3 dex higher ([Fe/H]=-1.2) than that of surrounding halo field stars. Together with a similar result for another major halo substructure, the Monoceros stream, these results support theoretical predictions that an early forming, smooth inner halo, is metal poor compared to high surface brightness material that have been accreted onto a later-forming outer halo. The mean metallicity of stars in the outer halo that are not associated with detectable clumps may still be more metal-poor than the bulk of inner-halo stars, as has been argued from other data sets.

preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

preprint2008arXiv

The Milky Way Tomography with SDSS: I. Stellar Number Density Distribution

Abridged: We estimate the distances to ~48 million stars detected by the Sloan Digital Sky Survey and map their 3D number density distribution in 100 < D < 20 kpc range over 6,500 deg^2 of sky. The data show strong evidence for a Galaxy consisting of an oblate halo, a disk component, and a number of localized overdensities with exponential disk parameters (bias-corrected for an assumed 35% binary fraction) H_1 = 300 pc, L_1 = 2600 pc, H_2 = 900 pc, L_2 = 3600 pc, and local density normalization of 12%. We find the halo to be oblate, with best-fit axis ratio c/a = 0.64, r^{-2.8} profile, and the local halo-to-thin disk normalization of 0.5%. We estimate the errors of derived model parameters to be no larger than ~20% (disk scales) and ~10% (thick disk normalization). While generally consistent with the above model, the density distribution shows a number of statistically significant localized deviations. We detect two overdensities in the thick disk region at (R, Z) ~ (6.5, 1.5)kpc and (R, Z) ~ (9.5, 0.8) kpc, and a remarkable density enhancement in the halo covering >1000deg^2 of sky towards the constellation of Virgo, at distances of ~6-20 kpc. Compared to a region symmetric with respect to the l=0 line, the Virgo overdensity is responsible for a factor of 2 number density excess and may be a nearby tidal stream or a low-surface brightness dwarf galaxy merging with the Milky Way. After removal of the resolved overdensities, the remaining data are consistent with a smooth density distribution; we detect no evidence of further unresolved clumpy substructure at scales ranging from ~50pc in the disk, to ~1 - 2 kpc in the halo.

preprint2007arXiv

The Virgo stellar over-density: Mapping the infall of the Sagittarius tidal stream onto the Milky Way disk

The recently discovered Virgo stellar over-density, which expands over \~1000deg^2 perpendicularly to the Galactic disk plane (7< Z <15 kpc, R~7 kpc), is the largest clump of tidal debris ever detected in the outer halo and is likely related with the accretion of a nearby dwarf galaxy by the Milky Way. We carry out N-body simulations of the Sagittarius stream to show that this giant stellar over-density is a confirmation of theoretical model predictions for the leading tail of the Sagittarius stream to cross the Milky Way plane in the Solar neighborhood. Radial velocity measurements are needed to confirm this association and to further constrain the shape of the Milky Way dark matter halo through a new generation of theoretical models. If the identification of Virgo over-density and the Sagittarius leading arm is correct, we predict highly negative radial velocities for the stars of Virgo over-density. The detection of this new portion of the Sagittarius tidal stream would represent an excellent target for the on-going and future kinematic surveys and for dark matter direct detection experiments in the proximity of the Sun.