Source author record

Zeljko Ivezic

Zeljko Ivezic 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

65works
8topics
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

65 published item(s)

preprint2026arXiv

Report on LEO satellite impacts on ground-based optical astronomy for the Rubin Observatory LSST

In August 2025 a workshop was convened to bring together experts to better understand steps that can be taken to mitigate the impact of satellite constellations on astronomical observations. At the time, just over 12,000 operational satellites were in low-Earth orbit (LEO). Although reflected sunlight and emissions all across the electromagnetic spectrum from artificial satellites impact scientific observations and the sky, the workshop focused on reflected sunlight in the wavelength range 330 nm to 1100 nm. This aligns with the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) planned imaging observations over the coming decade. Among other conclusions, we affirm previous recommendations that tracked satellite apparent magnitudes should be no brighter than 7th AB mag. The workshop participants discussed over 30 publications, reports, and presentations, and arrived at the Findings and Recommendations presented here. During the workshop, in addition to affirming many existing recommendations and best practices, the group discovered new issues and devised possible mitigations. These were nearly equally divided between advice to satellite builders and operators and to the observational astronomy community. While the workshop prioritized considerations for LSST, our hope is that many of the Findings and Recommendations will also apply to other observatories and constellations, and that all satellite companies will continue to engage in dialog with sky observers across the globe.

preprint2020arXiv

An Updated Version of the Unified Radio Catalog: A Multi-Wavelength Radio and Optical Catalog of Quasars and Radio Galaxies

We present an updated version of the Kimball & Ivezic (2008) radio catalog using updated radio and optical data, and new low-frequency radio data. The catalog, containing millions of radio sources, was created by consolidating large-area radio and optical surveys GB6 (6cm), FIRST (20cm), NVSS (20cm), WENSS (92cm), VLSSr (4m), and SDSS DR9 (optical). The region where all surveys overlap covers 3269 square degrees in the North Galactic Cap, and contains >160,000 20-cm sources, with about 12,000 detected in all five radio surveys and over one-third detected optically. Combining parameters from the sky surveys allows easy and efficient classification by radio and optical morphology and radio spectral index. The catalog is available at http://www.aoc.nrao.edu/~akimball/radiocat.shtml .

preprint2016arXiv

Radio Loud and Radio Quiet Quasars

We discuss 6 GHz JVLA observations covering a volume-limited sample of 178 low redshift ($0.2 < z < 0.3$) optically selected QSOs. Our 176 radio detections fall into two clear categories: (1) About $20$\% are radio-loud QSOs (RLQs) having spectral luminosities $L_6 \gtrsim 10^{\,23.2} \mathrm{~W~Hz}^{-1}$ primarily generated in the active galactic nucleus (AGN) responsible for the excess optical luminosity that defines a \emph{bona fide} QSO. (2) The radio-quiet QSOs (RQQs) have $10^{\,21} \lesssim L_6 \lesssim 10^{\,23.2} \mathrm{~W~Hz}^{-1}$ and radio sizes $\lesssim 10 \mathrm{~kpc}$, and we suggest that the bulk of their radio emission is powered by star formation in their host galaxies. "Radio silent" QSOs ($L_6 \lesssim 10^{\,21} \mathrm{~W~Hz}^{-1}$) are rare, so most RQQ host galaxies form stars faster than the Milky Way; they are not "red and dead" ellipticals. Earlier radio observations did not have the luminosity sensitivity $L_6 \lesssim 10^{\,21} \mathrm{~W~Hz}^{-1}$ needed to distinguish between such RLQs and RQQs. Strong, generally double-sided, radio emission spanning $\gg 10 \mathrm{~kpc}$ was found associated with 13 of the 18 RLQ cores having peak flux densities $S_\mathrm{p} > 5 \mathrm{~mJy~beam}^{-1}$ ($log(L) \gtrsim 24$). The radio luminosity function of optically selected QSOs and the extended radio emission associated with RLQs are both inconsistent with simple "unified" models that invoke relativistic beaming from randomly oriented QSOs to explain the difference between RLQs and RQQs. Some intrinsic property of the AGNs or their host galaxies must also determine whether or not a QSO appears radio loud.

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

Finding, characterizing and classifying variable sources in multi-epoch sky surveys: QSOs and RR Lyrae in PS1 3$π$ data

In area and depth, the Pan-STARRS1 (PS1) 3$π$ survey is unique among many-epoch, multi-band surveys and has enormous potential for all-sky identification of variable sources. PS1 has observed the sky typically seven times in each of its five bands ($grizy$) over 3.5 years, but unlike SDSS not simultaneously across the bands. Here we develop a new approach for quantifying statistical properties of non-simultaneous, sparse, multi-color lightcurves through light-curve structure functions, effectively turning PS1 into a $\sim 35$-epoch survey. We use this approach to estimate variability amplitudes and timescales $(ω_r, τ)$ for all point-sources brighter than $r_{\mathrm{P1}}=21.5$ mag in the survey. With PS1 data on SDSS Stripe 82 as ``ground truth", we use a Random Forest Classifier to identify QSOs and RR Lyrae based on their variability and their mean PS1 and WISE colors. We find that, aside from the Galactic plane, QSO and RR Lyrae samples of purity $\sim$75\% and completeness $\sim$92\% can be selected. On this basis we have identified a sample of $\sim 1,000,000$ QSO candidates, as well as an unprecedentedly large and deep sample of $\sim$150,000 RR Lyrae candidates with distances from $\sim$10 kpc to $\sim$120 kpc. Within the Draco dwarf spheroidal, we demonstrate a distance precision of 6\% for RR Lyrae candidates. We provide a catalog of all likely variable point sources and likely QSOs in PS1, a total of $25.8\times 10^6$ sources.

preprint2015arXiv

Periodograms for Multiband Astronomical Time Series

This paper introduces the multiband periodogram, a general extension of the well-known Lomb-Scargle approach for detecting periodic signals in time-domain data. In addition to advantages of the Lomb-Scargle method such as treatment of non-uniform sampling and heteroscedastic errors, the multiband periodogram significantly improves period finding for randomly sampled multiband light curves (e.g., Pan-STARRS, DES and LSST). The light curves in each band are modeled as arbitrary truncated Fourier series, with the period and phase shared across all bands. The key aspect is the use of Tikhonov regularization which drives most of the variability into the so-called base model common to all bands, while fits for individual bands describe residuals relative to the base model and typically require lower-order Fourier series. This decrease in the effective model complexity is the main reason for improved performance. We use simulated light curves and randomly subsampled SDSS Stripe 82 data to demonstrate the superiority of this method compared to other methods from the literature, and find that this method will be able to efficiently determine the correct period in the majority of LSST's bright RR Lyrae stars with as little as six months of LSST data. A Python implementation of this method, along with code to fully reproduce the results reported here, is available on GitHub.

preprint2015arXiv

Synergy between the Large Synoptic Survey Telescope and the Square Kilometre Array

We provide an overview of the science benefits of combining information from the Square Kilometre Array (SKA) and the Large Synoptic Survey Telescope (LSST). We first summarise the capabilities and timeline of the LSST and overview its science goals. We then discuss the science questions in common between the two projects, and how they can be best addressed by combining the data from both telescopes. We describe how weak gravitational lensing and galaxy clustering studies with LSST and SKA can provide improved constraints on the causes of the cosmological acceleration. We summarise the benefits to galaxy evolution studies of combining deep optical multi-band imaging with radio observations. Finally, we discuss the excellent match between one of the most unique features of the LSST, its temporal cadence in the optical waveband, and the time resolution of the SKA.

preprint2014arXiv

Introduction to astroML: Machine Learning for Astrophysics

Astronomy and astrophysics are witnessing dramatic increases in data volume as detectors, telescopes and computers become ever more powerful. During the last decade, sky surveys across the electromagnetic spectrum have collected hundreds of terabytes of astronomical data for hundreds of millions of sources. Over the next decade, the data volume will enter the petabyte domain, and provide accurate measurements for billions of sources. Astronomy and physics students are not traditionally trained to handle such voluminous and complex data sets. In this paper we describe astroML; an initiative, based on Python and scikit-learn, to develop a compendium of machine learning tools designed to address the statistical needs of the next generation of students and astronomical surveys. We introduce astroML and present a number of example applications that are enabled by this package.

preprint2014arXiv

The meaning of WISE colours - I. The Galaxy and its satellites

Through matches with the Sloan Digital Sky Survey (SDSS) catalogue we identify the location of various families of astronomical objects in WISE colour space. We identify reliable indicators that separate Galactic/local from extragalactic sources and concentrate here on the objects in our Galaxy and its closest satellites. We develop colour and magnitude criteria that are based only on WISE data to select asymptotic giant branch (AGB) stars with circumstellar dust shells, and separate them into O-rich and C-rich classes. With these criteria we produce an all-sky map for the count ratio of the two populations. The map reveals differences between the Galactic disc, the Magellanic Clouds and the Sgr Dwarf Spheroidal galaxy, as well as a radial gradient in the Large Magellanic Cloud (LMC) disc. We find that the C:O number ratio for dusty AGB stars increases with distance from the LMC centre about twice as fast as measured for near-IR selected samples of early AGB stars. Detailed radiative transfer models show that WISE colours are well explained by the emission of centrally heated dusty shells where the dust has standard properties of interstellar medium (ISM) grains. The segregation of different classes of objects in WISE colour space arises from differences in properties of the dust shells: those around young stellar objects have uniform density distributions while in evolved stars they have steep radial profiles.

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.

preprint2014arXiv

The SDSS-2MASS-WISE Ten Dimensional Stellar Color Locus

We present the fiducial main sequence stellar locus traced by 10 photometric colors observed by SDSS, 2MASS, and WISE. Median colors are determined using 1,052,793 stars with r-band extinction less than 0.125. We use this locus to measure the dust extinction curve relative to the r-band, which is consistent with previous measurements in the SDSS and 2MASS bands. The WISE band extinction coefficients are larger than predicted by standard extinction models. Using 13 lines of sight, we find variations in the extinction curve in H, Ks, and WISE bandpasses. Relative extinction decreases towards Galactic anti-center, in agreement with prior studies. Relative extinction increases with Galactic latitude, in contrast to previous observations. This indicates a universal mid-IR extinction law does not exist due to variations in dust grain size and chemistry with Galactocentric position. A preliminary search for outliers due to warm circumstellar dust is also presented, using stars with high signal-to-noise in the W3-band. We find 199 such outliers, identified by excess emission in Ks-W3. Inspection of SDSS images for these outliers reveals a large number of contaminants due to nearby galaxies. Six sources appear to be genuine dust candidates, yielding a fraction of systems with infrared excess of 0.12$\pm$0.05%.

preprint2013arXiv

AGN and Starburst Radio Emission from Optically Selected QSOs

We used the 1.4 GHz NVSS to study radio sources in two color-selected QSO samples: a volume-limited sample of 1313 QSOs defined by M_i < -23 in the redshift range 0.2 < z < 0.45 and a magnitude-limited sample of 2471 QSOs with m_r < 18.5 and 1.8 < z < 2.5. About 10% were detected above the 2.4 mJy NVSS catalog limit and are powered primarily by AGNs. The space density of the low-redshift QSOs evolves as rho proportional to (1+z)^6. In both redshift ranges the flux-density distributions and luminosity functions of QSOs stronger than 2.4 mJy are power laws, with no features to suggest more than one kind of radio source. Extrapolating the power laws to lower luminosities predicts the remaining QSOs should be extremely radio quiet, but they are not. Most were detected statistically on the NVSS images with median peak flux densities S_p(mJy/beam) ~ 0.3 and 0.05 in the low- and high-redshift samples, corresponding to 1.4 GHz spectral luminosities log[L(W/Hz)] ~ 22.7$ and 24.1, respectively. We suggest that the faint radio sources are powered by star formation at rates ~20 M_sun per year in the moderate luminosity (median M_i ~ -23.4) low-redshift QSOs and ~500 M_sun per year in the very luminous (M_i} ~ -27.5) high-redshift QSOs. Such luminous starbursts [ log(L / L_sun) ~ 11.2 and 12.6, respectively] are consistent with "quasar mode" accretion in which cold gas flows fuel both AGN and starburst.

preprint2013arXiv

AGN torus properties with WISE

The Wide-field Infrared Survey Explorer (WISE) has scanned the entire sky with unprecedented sensitivity in four infrared bands, at 3.4, 4.6, 12, and 22 micron. The WISE Point Source Catalog contains more than 560 million objects, among them hundreds of thousands of galaxies with Active Nuclei (AGN). While type 1 AGN, owing to their bright and unobscured nature, are easy to detect and constitute a rather complete and unbiased sample, their type 2 counterparts, postulated by AGN unification, are not as straightforward to identify. Matching the WISE catalog with known QSOs in the Sloan Digital Sky Survey we confirm previous identification of the type 1 locus in the WISE color space. Using a very large database of the popular CLUMPY torus models, we find the colors of the putative type 2 counterparts, and also, for the first time, predict their number vs. flux relation that can be expected to be observed in any given WISE color range. This will allow us to put statistically very significant constraints on the torus parameters. Our results are a successful test of the AGN unification scheme.

preprint2013arXiv

All Weather Calibration of Wide Field Optical and NIR Surveys

The science goals for ground-based large-area surveys, such as the Dark Energy Survey, Pan-STARRS, and the Large Synoptic Survey Telescope, require calibration of broadband photometry that is stable in time and uniform over the sky to precisions of a per cent or better. This performance will need to be achieved with data taken over the course of many years, and often in less than ideal conditions. This paper describes a strategy to achieve precise internal calibration of imaging survey data taken in less than photometric conditions, and reports results of an observational study of the techniques needed to implement this strategy. We find that images of celestial fields used in this case study with stellar densities of order one per arcmin-squared and taken through cloudless skies can be calibrated with relative precision of 0.5 per cent (reproducibility). We report measurements of spatial structure functions of cloud absorption observed over a range of atmospheric conditions, and find it possible to achieve photometric measurements that are reproducible to 1 per cent in images that were taken through cloud layers that transmit as little as 25 per cent of the incident optical flux (1.5 magnitudes of extinction). We find, however, that photometric precision below 1 per cent is impeded by the thinnest detectable cloud layers. We comment on implications of these results for the observing strategies of future surveys.

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.

preprint2013arXiv

Optical selection of quasars: SDSS and LSST

Over the last decade, quasar sample sizes have increased from several thousand to several hundred thousand, thanks mostly to SDSS imaging and spectroscopic surveys. LSST, the next-generation optical imaging survey, will provide hundreds of detections per object for a sample of more than ten million quasars with redshifts of up to about seven. We briefly review optical quasar selection techniques, with emphasis on methods based on colors, variability properties and astrometric behavior.

preprint2013arXiv

Optical variability of quasars: a damped random walk

A damped random walk is a stochastic process, defined by an exponential covariance matrix that behaves as a random walk for short time scales and asymptotically achieves a finite variability amplitude at long time scales. Over the last few years, it has been demonstrated, mostly but not exclusively using SDSS data, that a damped random walk model provides a satisfactory statistical description of observed quasar variability in the optical wavelength range, for rest-frame timescales from 5 days to 2000 days. The best-fit characteristic timescale and asymptotic variability amplitude scale with the luminosity, black hole mass, and rest wavelength, and appear independent of redshift. In addition to providing insights into the physics of quasar variability, the best-fit model parameters can be used to efficiently separate quasars from stars in imaging surveys with adequate long-term multi-epoch data, such as expected from LSST.

preprint2012arXiv

A Description of Quasar Variability Measured Using Repeated SDSS and POSS Imaging

We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along Stripe 82 (S82), whereas the remaining ~25,000 have 2-3 observations due to scan overlaps. The observed time lags span the range from a day to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 4 years, and at rest-frame wavelengths from 1000A to 6000A. We publicly release a user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Survey, and we use it to analyze the ensemble properties of quasar variability. Based on a damped random walk (DRW) model defined by a characteristic time scale and an asymptotic variability amplitude that scale with the luminosity, black hole mass, and rest wavelength for individual quasars calibrated in S82, we can fully explain the ensemble variability statistics of the non-S82 quasars such as the exponential distribution of large magnitude changes. All available data are consistent with the DRW model as a viable description of the optical continuum variability of quasars on time scales of ~5-2000 days in the rest frame. We use these models to predict the incidence of quasar contamination in transient surveys such as those from PTF and LSST.

preprint2012arXiv

Characterizing the Optical Variability of Bright Blazars: Variability-Based Selection of Fermi AGN

We investigate the use of optical photometric variability to select and identify blazars in large-scale time-domain surveys, in part to aid in the identification of blazar counterparts to the ~30% of gamma-ray sources in the Fermi 2FGL catalog still lacking reliable associations. Using data from the optical LINEAR asteroid survey, we characterize the optical variability of blazars by fitting a damped random walk model to individual light curves with two main model parameters, the characteristic timescales of variability (tau), and driving amplitudes on short timescales (sigma). Imposing cuts on minimum tau and sigma allows for blazar selection with high efficiency E and completeness C. To test the efficacy of this approach, we apply this method to optically variable LINEAR objects that fall within the several-arcminute error ellipses of gamma-ray sources in the Fermi 2FGL catalog. Despite the extreme stellar contamination at the shallow depth of the LINEAR survey, we are able to recover previously-associated optical counterparts to Fermi AGN with E > 88% and C = 88% in Fermi 95% confidence error ellipses having semimajor axis r < 8'. We find that the suggested radio counterpart to Fermi source 2FGL J1649.6+5238 has optical variability consistent with other gamma-ray blazars, and is likely to be the gamma-ray source. Our results suggest that the variability of the non-thermal jet emission in blazars is stochastic in nature, with unique variability properties due to the effects of relativistic beaming. After correcting for beaming, we estimate that the characteristic timescale of blazar variability is ~3 years in the rest-frame of the jet, in contrast with the ~320 day disk flux timescale observed in quasars. The variability-based selection method presented will be useful for blazar identification in time-domain optical surveys, and is also a probe of jet physics.

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

Disclosing the Radio Loudness Distribution Dichotomy in Quasars: An Unbiased Monte Carlo Approach Applied to the SDSS-FIRST Quasar Sample

We investigate the dichotomy in the radio loudness distribution of quasars by modelling their radio emission and various selection effects using a Monte Carlo approach. The existence of two physically distinct quasar populations, the radio-loud and radio-quiet quasars, is controversial and over the last decade a bimodal distribution of radio loudness of quasars has been both affirmed and disputed. We model the quasar radio luminosity distribution with simple unimodal and bimodal distribution functions. The resulting simulated samples are compared to a fiducial sample of 8,300 quasars drawn from the SDSS DR7 Quasar Catalog and combined with radio observations from the FIRST survey. Our results indicate that the SDSS-FIRST sample is best described by a radio loudness distribution which consists of two components, with 12+/-1 % of sources in the radio-loud component. On the other hand, the evidence for a local minimum in the loudness distribution (bimodality) is not strong and we find that previous claims for its existence were probably affected by the incompleteness of the FIRST survey close to its faint limit. We also investigate the redshift and luminosity dependence of the radio loudness distribution and find tentative evidence that at high redshift radio-loud quasars were rarer, on average "louder", and exhibited a smaller range in radio loudness. In agreement with other recent work, we conclude that the SDSS-FIRST sample strongly suggests that the radio loudness distribution of quasars is not a universal function, and that more complex models than presented here are needed to fully explain available observations.

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

Monitoring Quasar Colour Variability in Stripe 82

Broad Absorption Line (BAL) trough variability is predominantly due to cloud motion transverse to our line of sight. The rate at which the variability occurs indicates the velocity of the cloud, which can provide constraints on the cloud's distance from the central source. This requires detailed spectroscopy during a variability event. Such spectra have proven elusive, suggesting either the timescale of variability is too short to be caught, or too long to notice until a sufficient amount of time has passed. Photometric monitoring of BAL quasar colours may potentially be used as an early warning system to trigger time resolved spectroscopic monitoring of BAL variability. Towards this end, we are analyzing both BAL and non-BAL colour variability using time series photometry from Stripe 82 in the Sloan Digital Sky Survey.

preprint2012arXiv

The Sloan Digital Sky Survey quasar catalog: ninth data release

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $Λ$CDM cosmology with H0 = 70 km/s/Mpc, $Ω_{\rm M}$ = 0.3, and $Ω_Λ$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

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

Dusty Tori of Luminous Type 1 Quasars at z \sim 2

We present Spitzer infrared spectra and ultra-violet to mid-infrared spectral energy distributions (SEDs) of 25 luminous type 1 quasars at z \sim 2. In general, the spectra show a bump peaking around 3 μm, and the 10 μm silicate emission feature. The 3 μm emission is identified with hot dust emission at its sublimation temperature. We explore two approaches to modeling the SED: (i) using the Clumpy model SED from Nenkova et al. (2008a), and (ii) the Clumpy model SED, and an additional blackbody component to represent the 3 μm emission. In the first case, a parameter search of \sim 1.25 million Clumpy models shows: (i) if we ignore the UV-to-near-IR SED, models fit the 2-8 μm region well, but not the 10 μm feature; (ii) if we include the UV-to-near-IR SED in the fit, models do not fit the 2-8 μm region. The observed 10 μm features are broader and shallower than those in the best-fit models in the first approach. In the second case, the shape of the 10 μm feature is better reproduced by the Clumpy models. The additional blackbody contribution in the 2-8 μm range allows Clumpy models dominated by cooler temperatures (T < 800K) to better fit the 8-12μm SED. A centrally concentrated distribution of a small number of torus clouds is required in the first case, while in the second case the clouds are more spread out radially. The temperature of the blackbody component is ~ 1200 K as expected for graphite grains.

preprint2011arXiv

Efficient use of simultaneous multi-band observations for variable star analysis

The luminosity changes of most types of variable stars are correlated in the different wavelengths, and these correlations may be exploited for several purposes: for variability detection, for distinction of microvariability from noise, for period search or for classification. Principal component analysis is a simple and well-developed statistical tool to analyze correlated data. We will discuss its use on variable objects of Stripe 82 of the Sloan Digital Sky Survey, with the aim of identifying new RR Lyrae and SX Phoenicis-type candidates. The application is not straightforward because of different noise levels in the different bands, the presence of outliers that can be confused with real extreme observations, under- or overestimated errors and the dependence of errors on the magnitudes. These particularities require robust methods to be applied together with the principal component analysis. The results show that PCA is a valuable aid in variability analysis with multi-band data.

preprint2011arXiv

Formation and Evolution of the Disk System of the Milky Way: [alpha/Fe] Ratios and Kinematics of the SEGUE G-Dwarf Sample

We employ measurements of the [alpha/Fe] ratio derived from low-resolution (R~2000) spectra of 17,277 G-type dwarfs from the SEGUE survey to separate them into likely thin- and thick-disk subsamples. Both subsamples exhibit strong gradients of orbital rotational velocity with metallicity, of opposite signs, -20 to -30 km/s/dex for the thin-disk and +40 to +50 km/s/dex for the thick-disk population. The rotational velocity is uncorrelated with Galactocentric distance for the thin-disk subsample, and exhibits a small trend for the thick-disk subsample. The rotational velocity decreases with distance from the plane for both disk components, with similar slopes (-9.0 {\pm} 1.0 km/s/kpc). Thick-disk stars exhibit a strong trend of orbital eccentricity with metallicity (about -0.2/dex), while the eccentricity does not change with metallicity for the thin-disk subsample. The eccentricity is almost independent of Galactocentric radius for the thin-disk population, while a marginal gradient of the eccentricity with radius exists for the thick-disk population. Both subsamples possess similar positive gradients of eccentricity with distance from the Galactic plane. The shapes of the eccentricity distributions for the thin- and thick-disk populations are independent of distance from the plane, and include no significant numbers of stars with eccentricity above 0.6. Among several contemporary models of disk evolution we consider, radial migration appears to have played an important role in the evolution of the thin-disk population, but possibly less so for the thick disk, relative to the gas-rich merger or disk heating scenarios. We emphasize that more physically realistic models and simulations need to be constructed in order to carry out the detailed quantitative comparisons that our new data enable.

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 Case for the Dual Halo of the Milky Way

Carollo et al. have recently resolved the stellar population of the Milky Way halo into at least two distinct components, an inner halo and an outer halo. This result has been criticized by Schoenrich et al., who claim that the retrograde signature associated with the outer halo is due to the adoption of faulty distances. We refute this claim, and demonstrate that the Schoenrich et al. photometric distances are themselves flawed because they adopted an incorrect main-sequence absolute magnitude relationship from the work of Ivezić et al. When compared to the recommended relation from Ivezić et al., which is tied to a Milky Way globular cluster distance scale and accounts for age and metallicity effects, the relation adopted by Schoenrich et al. yields up to 18% shorter distances for stars near the main-sequence turnoff (TO). Use of the correct relationship yields agreement between the distances assigned by Carollo et al. and Ivezić et al. for low-metallicity dwarfs to within 6-10%. Schoenrich et al. also point out that intermediate-gravity stars (3.5 <= log g <= 4.0) with colors redder than the TO region are likely misclassified, with which we concur. We implement a new procedure to reassign luminosity classifications for the TO stars that require it. New derivations of the rotational behavior demonstrate that the retrograde signature and high velocity dispersion of the outer-halo population remains. We summarize additional lines of evidence for a dual halo, including a test of the retrograde signature based on proper motions alone, and conclude that the preponderance of evidence strongly rejects the single-halo interpretation.

preprint2011arXiv

The SDSS Coadd: 275 deg^2 of Deep SDSS Imaging on Stripe 82

We present details of the construction and characterization of the coaddition of the Sloan Digital Sky Survey Stripe 82 \ugriz\ imaging data. This survey consists of 275 deg$^2$ of repeated scanning by the SDSS camera of $2.5\arcdeg$ of $δ$ over $-50\arcdeg \le α\le 60\arcdeg$ centered on the Celestial Equator. Each piece of sky has $\sim 20$ runs contributing and thus reaches $\sim2$ magnitudes fainter than the SDSS single pass data, i.e. to $r\sim 23.5$ for galaxies. We discuss the image processing of the coaddition, the modeling of the PSF, the calibration, and the production of standard SDSS catalogs. The data have $r$-band median seeing of 1.1\arcsec, and are calibrated to $\le 1%$. Star color-color, number counts, and psf size vs modelled size plots show the modelling of the PSF is good enough for precision 5-band photometry. Structure in the psf-model vs magnitude plot show minor psf mis-modelling that leads to a region where stars are being mis-classified as galaxies, and this is verified using VVDS spectroscopy. As this is a wide area deep survey there are a variety of uses for the data, including galactic structure, photometric redshift computation, cluster finding and cross wavelength measurements, weak lensing cluster mass calibrations, and cosmic shear measurements.

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).

preprint2011arXiv

The Two-Component Radio Luminosity Function of QSOs: Star Formation and AGN

Despite decades of study, it remains unclear whether there are distinct radio-loud and radio-quiet populations of quasi-stellar objects (QSOs). Early studies were limited by inhomogeneous QSO samples, inadequate sensitivity to probe the radio-quiet population, and degeneracy between redshift and luminosity for flux-density-limited samples. Our new 6 GHz EVLA observations allow us for the first time to obtain nearly complete (97%) radio detections in a volume-limited color-selected sample of 179 QSOs more luminous than M_i = -23 from the Sloan Digital Sky Survey (SDSS) Data Release Seven in the narrow redshift range 0.2 < z < 0.3. The dramatic improvement in radio continuum sensitivity made possible with the new EVLA allows us, in 35 minutes of integration, to detect sources as faint as 20 microJy, or log[L_6 (W/Hz)] ~ 21.5 at z = 0.25, well below the radio luminosity, log[L_6 (W/Hz)] ~ 22.5, that separates star-forming galaxies from radio-loud active galactic nuclei (AGNs) driven by accretion onto a super-massive black hole. We calculate the radio luminosity function (RLF) for these QSOs using three constraints: (a) EVLA 6 GHz observations for log[L_6 (W/Hz)] < 23.5, (b) NRAO-VLA Sky Survey (NVSS) observations for log[L_6 (W/Hz)] > 23.5, and (c) the total number of SDSS QSOs in our volume-limited sample. We show that the RLF can be explained as a superposition of two populations, dominated by AGNs at the bright end and star formation in the QSO host galaxies at the faint end.

preprint2010arXiv

A closer view of the radio-FIR correlation: disentangling the contributions of star formation and AGN activity

We extend the Unified Radio Catalog, a catalog of sources detected by various (NVSS, FIRST, WENSS, GB6) radio surveys, and SDSS, to IR wavelengths by matching it to the IRAS Point and Faint Source catalogs. By fitting each NVSS-selected galaxy's NUV-NIR spectral energy distribution (SED) with stellar population synthesis models we add to the catalog star formation rates, stellar masses, and attenuations.We further add information about optical emission line properties for NVSS-selected galaxies with available SDSS spectroscopy. Using an NVSS 20cm (F_{1.4GHz} ge 2.5mJy) selected sample, matched to the SDSS spectroscopic ("main" galaxy and quasar) catalogs and IRAS data (0.04<z le 0.2) we perform an in depth analysis of the radio-FIR correlation for various types of galaxies, separated into i) quasars, ii) star forming, iii) composite, iv) Seyfert, v) LINER and vi) absorption line galaxies using the standard optical spectroscopic diagnostic tools. We utilize SED-based star formation rates to independently quantify the source of radio and FIR emission in our galaxies. Our results show that Seyfert galaxies have FIR/radio ratios lower than, but still within the scatter of, the canonical value due to an additional (likely AGN) contribution to their radio continuum emission. Furthermore, IR-detected absorption and LINER galaxies are on average strongly dominated by AGN activity in both their FIR and radio emission; however their average FIR/radio ratio is consistent with that expected for star forming galaxies. In summary, we find that most AGN-containing galaxies in our NVSS-IRAS-SDSS sample have FIR/radio flux ratios indistinguishable from those of the star-forming galaxies that define the radio-FIR correlation. Thus, attempts to separate AGNs from star-forming galaxies by their FIR/radio flux ratios alone can separate only a small fraction of the AGNs, such as the radio-loud quasars.

preprint2010arXiv

Detection of Periodic Variability in Simulated QSO Light Curves

Periodic light curve behavior predicted for some binary black hole systems might be detected in large samples, such as the multi-million quasar sample expected from the Large Synoptic Survey Telescope (LSST). We investigate the false-alarm probability for the discovery of a periodic signal in light curves simulated using damped random walk (DRW) model. This model provides a good description of observed light curves, and does not include periodic behavior. We used the Lomb-Scargle periodogram to search for a periodic signal in a million simulated light curves that properly sample the DRW parameter space, and the LSST cadence space. We find that even a very conservative threshold for the false-alarm probability still yields thousands of "good" binary black hole candidates. We conclude that the future claims for binary black holes based on Lomb-Scargle analysis of LSST light curves will have to be interpreted with caution.

preprint2010arXiv

Dusty winds II. Observational Implications

We compare observations of AGB stars and predictions of the Elitzur & Ivezic (2001) steady-state radiatively driven dusty wind model. The model results are described by a set of similarity functions of a single independent variable, and imply general scaling relations among the system parameters. We find that the model properly reproduces various correlations among the observed quantities and demonstrate that dust drift through the gas has a major impact on the structure of most winds. From data for nearby oxygen-rich and carbon-rich mass-losing stars we find that (1) the dispersion in grain properties within each group is rather small; (2) both the dust cross-section per gas particle and the dust-to-gas mass ratio are similar for the two samples even though the stellar atmospheres and grain properties are very different; (3) the dust abundance in both outflows is significantly below the Galactic average, indicating that most of the Galactic dust is not stardust - contrary to popular belief, but in support of Draine (2009). Our model results can be easily applied to recent massive data sets, such as the Spitzer SAGE survey of the Large Magellanic Cloud, and incorporated in galaxy evolution models.

preprint2010arXiv

Halo Velocity Groups in the Pisces Overdensity

We report spectroscopic observations with the Gemini South Telescope of 5 faint V~20 RR Lyrae stars associated with the Pisces overdensity. At a heliocentric and galactocentric distance of ~80 kpc, this is the most distant substructure in the Galactic halo known to date. We combined our observations with literature data and confirmed that the substructure is composed of two different kinematic groups. The main group contains 8 stars and has <V_{gsr}> = 50 km/s, while the second group contains four stars at a velocity of <V_{gsr}> = -52 km/s, where V_{gsr} is the radial velocity in the galactocentric standard of rest. The metallicity distribution of RR Lyrae stars in the Pisces overdensity is centered on [Fe/H]=-1.5 dex and has a width of 0.3 dex. The new data allowed us to establish that both groups are spatially extended making it very unlikely that they are bound systems, and are more likely to be debris of a tidally disrupted galaxy or galaxies. Due to small sky coverage, it is still unclear whether these groups have the same or different progenitors.

preprint2010arXiv

Photometric Response Functions of the SDSS Imager

The monochromatic illumination system is constructed to carry out in situ measurements of the response function of the mosaicked CCD imager used in the Sloan Digital Sky Survey (SDSS). The system is outlined and the results of the measurements, mostly during the first 6 years of the SDSS, are described. We present the reference response functions for the five colour passbands derived from these measurements, and discuss column to column variations and variations in time, and also their effects on photometry. We also discuss the effect arising from various, slightly different response functions of the associated detector systems that were used to give SDSS photometry. We show that the calibration procedures of SDSS remove these variations reasonably well with the resulting final errors from variant response functions being unlikely to be larger than 0.01 mag for g, r, i, and z bands over the entire duration of the survey. The considerable aging effect is uncovered in the u band, the response function showing a 30% decrease in the throughput in the short wavelength side during the survey years, which potentially causes a systematic error in photometry. The aging effect is consistent with variation of the instrumental sensitivity in u-band, which is calibrated out. The expected colour variation is consistent with measured colour variation in the catalog of repeated photometry. The colour variation is delta (u-g) ~ 0.01 for most stars, and at most delta (u-g) ~ 0.02 mag for those with extreme colours. We verified in the final catalogue that no systematic variations in excess of 0.01 mag are detected in the photometry which can be ascribed to aging and/or seasonal effects except for the secular u-g colour variation for stars with extreme colours.

preprint2010arXiv

Principal Component Analysis of SDSS Stellar Spectra

We apply Principal Component Analysis (PCA) to ~100,000 stellar spectra obtained by the Sloan Digital Sky Survey (SDSS). In order to avoid strong non-linear variation of spectra with effective temperature, the sample is binned into 0.02 mag wide intervals of the g-r color (-0.20<g-r<0.90, roughly corresponding to MK spectral types A3 to K3), and PCA is applied independently for each bin. In each color bin, the first four eigenspectra are sufficient to describe the observed spectra within the measurement noise. We discuss correlations of eigencoefficients with metallicity and gravity estimated by the Sloan Extension for Galactic Understanding and Exploration (SEGUE) Stellar Parameters Pipeline. The resulting high signal-to-noise mean spectra and the other three eigenspectra are made publicly available. These data can be used to generate high quality spectra for an arbitrary combination of effective temperature, metallicity, and gravity within the parameter space probed by the SDSS. The SDSS stellar spectroscopic database and the PCA results presented here offer a convenient method to classify new spectra, to search for unusual spectra, to train various spectral classification methods, and to synthesize accurate colors in arbitrary optical bandpasses.

preprint2010arXiv

Structure and Kinematics of the Stellar Halos and Thick Disks of the Milky Way Based on Calibration Stars from SDSS DR7

The structure and kinematics of the recognized stellar components of the Milky Way are explored, based on well-determined atmospheric parameters and kinematic quantities for 32360 calibration stars from the Sloan Digital Sky Survey (SDSS) and its first extension, (SDSS-II), which included the sub-survey SEGUE: Sloan Extension for Galactic Understanding and Exploration. Full space motions for a sub-sample of 16920 stars, exploring a local volume within 4 kpc of the Sun, are used to derive velocity ellipsoids for the inner- and outer-halo components of the Galaxy, as well as for the canonical thick-disk and proposed metal-weak thick-disk populations. We first examine the question of whether the data require the presence of at least a two-component halo in order to account for the rotational behavior of likely halo stars in the local volume, and whether more than two components are needed. We also address the question of whether the proposed metal-weak thick disk is kinematically and chemically distinct from the canonical thick disk. In addition, we consider the fractions of each component required to understand the nature of the observed kinematic behavior of the stellar populations of the Galaxy as a function of distance from the plane. Scale lengths and scale heights for the thick-disk and metal-weak thick-disk components are determined. Spatial density profiles for the inner- and outer-halo populations are inferred from a Jeans Theorem analysis. The full set of calibration stars (including those outside the local volume) is used to test for the expected changes in the observed stellar metallicity distribution function with distance above the Galactic plane in-situ, due to the changing contributions from the underlying stellar populations. [abridged]

preprint2010arXiv

The Genesis of the Milky Way's Thick Disk via Stellar Migration

The separation of the Milky Way disk into a thin and thick component is supported by differences in the spatial, kinematic and metallicity distributions of their stars. These differences have led to the view that the thick disk formed early via a cataclysmic event and constitutes fossil evidence of the hierarchical growth of the Milky Way. We show here, using N-body simulations, how a double-exponential vertical structure, with stellar populations displaying similar dichotomies can arise purely through internal evolution. In this picture, stars migrate radially, while retaining nearly circular orbits, as described by Sellwood & Binney (2002). As stars move outwards they populate a thickened component. Such stars found at the present time in the solar neighborhood formed early in the disk's history at smaller radii where stars are more metal-poor and alpha-enhanced, leading to the properties observed for thick disk stars. Classifying stars as members of the thin or thick disk by either velocity or metallicity leads to an apparent separation in the other property as observed. This scenario is supported by the SDSS observation that stars in the transition region do not show any correlation between rotational velocity and metallicity. The good qualitative agreement between our simulation and observations in the Milky Way hint that the thick disk may be a ubiquitous galaxy feature generated by stellar migration. Nonetheless, we cannot exclude that some fraction of the thick disk is a fossil of a past more violent history, nor can this scenario explain thick disks in all galaxies.

preprint2010arXiv

The Luminosity and Mass Functions of Low-Mass Stars in the Galactic Disk: II. The Field

We report on new measurements of the luminosity function (LF) and mass function (MF) of field low-mass dwarfs derived from Sloan Digital Sky Survey (SDSS) Data Release 6 (DR6) photometry. The analysis incorporates ~15 million low-mass stars (0.1 Msun < M < 0.8 Msun), spread over 8,400 square degrees. Stellar distances are estimated using new photometric parallax relations, constructed from ugriz photometry of nearby low-mass stars with trigonometric parallaxes. We use a technique that simultaneously measures Galactic structure and the stellar LF from 7 < M_r < 16. We compare the LF to previous studies and convert to a MF using the mass-luminosity relations of Delfosse et al., 2000. The system MF, measured over -1.0 < log M/Msun < -0.1, is well-described by a log-normal distribution with Mo = 0.25 Msun. We stress that our results should not be extrapolated to other mass regimes. Our work generally agrees with prior low-mass stellar MFs and places strong constraints on future star-formation studies of the Milky Way.

preprint2010arXiv

The Sloan Digital Sky Survey Quasar Catalog V. Seventh Data Release

We present the fifth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog, which is based upon the SDSS Seventh Data Release. The catalog, which contains 105,783 spectroscopically confirmed quasars, represents the conclusion of the SDSS-I and SDSS-II quasar survey. The catalog consists of the SDSS objects that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s or have interesting/complex absorption features, are fainter than i > 15.0 and have highly reliable redshifts. The catalog covers an area of 9380 deg^2. The quasar redshifts range from 0.065 to 5.46, with a median value of 1.49; the catalog includes 1248 quasars at redshifts greater than four, of which 56 are at redshifts greater than five. The catalog contains 9210 quasars with i < 18; slightly over half of the entries have i< 19. For each object the catalog presents positions accurate to better than 0.1" rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800-9200 Ang. at a spectral resolution R = 2000 the spectra can be retrieved from the SDSS public database using the information provided in the catalog. Over 96% of the objects in the catalog were discovered by the SDSS. We also include a supplemental list of an additional 207 quasars with SDSS spectra whose archive photometric information is incomplete.

preprint2009arXiv

Baryon Acoustic Oscillations in the Sloan Digital Sky Survey Data Release 7 Galaxy Sample

The spectroscopic Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample represents the final set of galaxies observed using the original SDSS target selection criteria. We analyse the clustering of galaxies within this sample, including both the Luminous Red Galaxy (LRG) and Main samples, and also include the 2-degree Field Galaxy Redshift Survey (2dFGRS) data. Baryon Acoustic Oscillations are observed in power spectra measured for different slices in redshift; this allows us to constrain the distance--redshift relation at multiple epochs. We achieve a distance measure at redshift z=0.275, of r_s(z_d)/D_V(0.275)=0.1390+/-0.0037 (2.7% accuracy), where r_s(z_d) is the comoving sound horizon at the baryon drag epoch, D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3), D_A(z) is the angular diameter distance and H(z) is the Hubble parameter. We find an almost independent constraint on the ratio of distances D_V(0.35)/D_V(0.2)=1.736+/-0.065, which is consistent at the 1.1sigma level with the best fit Lambda-CDM model obtained when combining our z=0.275 distance constraint with the WMAP 5-year data. The offset is similar to that found in previous analyses of the SDSS DR5 sample, but the discrepancy is now of lower significance, a change caused by a revised error analysis and a change in the methodology adopted, as well as the addition of more data. Using WMAP5 constraints on Omega_bh^2 and Omega_ch^2, and combining our BAO distance measurements with those from the Union Supernova sample, places a tight constraint on Omega_m=0.286+/-0.018 and H_0 = 68.2+/-2.2km/s/Mpc that is robust to allowing curvature and non-Lambda dark energy. This result is independent of the behaviour of dark energy at redshifts greater than those probed by the BAO and supernova measurements. (abridged)

preprint2009arXiv

Cosmological Constraints from the Clustering of the Sloan Digital Sky Survey DR7 Luminous Red Galaxies

We present the power spectrum of the reconstructed halo density field derived from a sample of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Seventh Data Release (DR7). The halo power spectrum has a direct connection to the underlying dark matter power for k <= 0.2 h/Mpc, well into the quasi-linear regime. This enables us to use a factor of ~8 more modes in the cosmological analysis than an analysis with kmax = 0.1 h/Mpc, as was adopted in the SDSS team analysis of the DR4 LRG sample (Tegmark et al. 2006). The observed halo power spectrum for 0.02 < k < 0.2 h/Mpc is well-fit by our model: chi^2 = 39.6 for 40 degrees of freedom for the best fit LCDM model. We find Ω_m h^2 * (n_s/0.96)^0.13 = 0.141^{+0.009}_{-0.012} for a power law primordial power spectrum with spectral index n_s and Ω_b h^2 = 0.02265 fixed, consistent with CMB measurements. The halo power spectrum also constrains the ratio of the comoving sound horizon at the baryon-drag epoch to an effective distance to z=0.35: r_s/D_V(0.35) = 0.1097^{+0.0039}_{-0.0042}. Combining the halo power spectrum measurement with the WMAP 5 year results, for the flat LCDM model we find Ω_m = 0.289 +/- 0.019 and H_0 = 69.4 +/- 1.6 km/s/Mpc. Allowing for massive neutrinos in LCDM, we find \sum m_ν < 0.62 eV at the 95% confidence level. If we instead consider the effective number of relativistic species Neff as a free parameter, we find Neff = 4.8^{+1.8}_{-1.7}. Combining also with the Kowalski et al. (2008) supernova sample, we find Ω_{tot} = 1.011 +/- 0.009 and w = -0.99 +/- 0.11 for an open cosmology with constant dark energy equation of state w.

preprint2009arXiv

HI Selected Galaxies in the Sloan Digital Sky Survey I: Optical Data

We present the optical data for 195 HI-selected galaxies that fall within both the Sloan Digital Sky Survey (SDSS) and the Parkes Equatorial Survey (ES). The photometric quantities have been independently recomputed for our sample using a new photometric pipeline optimized for large galaxies, thus correcting for SDSS's limited reliability for automatic photometry of angularly large or low surface brightness (LSB) galaxies. We outline the magnitude of the uncertainty in the SDSS catalog-level photometry and derive a quantitative method for correcting the over-sky subtraction in the SDSS photometric pipeline. The main thrust of this paper is to present the ES/SDSS sample and discuss the methods behind the improved photometry, which will be used in future scientific analysis. We present the overall optical properties of the sample and briefly compare to a volume-limited, optically-selected sample. Compared to the optically-selected SDSS sample (in the similar volume), HI-selected galaxies are bluer and more luminous (fewer dwarf ellipticals and more star formation). However, compared to typical SDSS galaxy studies, which have their own selection effects, our sample is bluer, fainter and less massive.

preprint2009arXiv

HI Selected Galaxies in the Sloan Digital Sky Survey II: The Colors of Gas-Rich Galaxies

We utilize color information for an HI-selected sample of 195 galaxies to explore the star formation histories and physical conditions that produce the observed colors. We show that the HI selection creates a significant offset towards bluer colors that can be explained by enhanced recent bursts of star formation. There is also no obvious color bimodality, because the HI selection restricts the sample to bluer, actively star forming systems, diminishing the importance of the red sequence. Rising star formation rates are still required to explain the colors of galaxies bluer than g-r < 0.3. We also demonstrate that the colors of the bluest galaxies in our sample are dominated by emission lines and that stellar population synthesis models alone (without emission lines) are not adequate for reproducing many of the galaxy colors. These emission lines produce large changes in the r-i colors but leave the g-r color largely unchanged. In addition, we find an increase in the dispersion of galaxy colors at low masses that may be the result of a change in the star formation process in low-mass galaxies.

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.

preprint2009arXiv

Photometric constraints on white dwarfs and the identification of extreme objects

It is possible to reliably identify white dwarfs (WDs) without recourse to spectra, instead using photometric and astrometric measurements to distinguish them from Main Sequence stars and quasars. WDs' colours can also be used to infer their intrinsic properties (effective temperature, surface gravity, etc.), but the results obtained must be interpreted with care. The difficulties stem from the existence of a solid angle degeneracy, as revealed by a full exploration of the likelihood, although this can be masked if a simple best-fit approach is used. Conversely, this degeneracy can be broken if a Bayesian approach is adopted, as it is then possible to utilise the prior information on the surface gravities of WDs implied by spectroscopic fitting. The benefits of such an approach are particularly strong when applied to outliers, such as the candidate halo and ultra-cool WDs identified by Vidrih et al. (2007). A reanalysis of these samples confirms their results for the latter sample but suggests that that most of the halo candidates are thick disk WDs in the tails of the photometric noise distribution.

preprint2009arXiv

Probing Spectroscopic Variability of Galaxies & Narrow-Line Active Galactic Nuclei in the Sloan Digital Sky Survey

Under the unified model for active galactic nuclei (AGNs), narrow-line (Type 2) AGNs are, in fact, broad-line (Type 1) AGNs but each with a heavily obscured accretion disk. We would therefore expect the optical continuum emission from Type 2 AGN to be composed mainly of stellar light and non-variable on the time-scales of months to years. In this work we probe the spectroscopic variability of galaxies and narrow-line AGNs using the multi-epoch data in the Sloan Digital Sky Survey (SDSS) Data Release 6. The sample contains 18,435 sources for which there exist pairs of spectroscopic observations (with a maximum separation in time of ~700 days) covering a wavelength range of 3900-8900 angstrom. To obtain a reliable repeatability measurement between each spectral pair, we consider a number of techniques for spectrophotometric calibration resulting in an improved spectrophotometric calibration of a factor of two. From these data we find no obvious continuum and emission-line variability in the narrow-line AGNs on average -- the spectroscopic variability of the continuum is 0.07+/-0.26 mag in the g band and, for the emission-line ratios log10([NII]/Halpha) and log10([OIII]/Hbeta), the variability is 0.02+/-0.03 dex and 0.06+/-0.08 dex, respectively. From the continuum variability measurement we set an upper limit on the ratio between the flux of varying spectral component, presumably related to AGN activities, and that of host galaxy to be ~30%. We provide the corresponding upper limits for other spectral classes, including those from the BPT diagram, eClass galaxy classification, stars and quasars.

preprint2008arXiv

AGN Dusty Tori: II. Observational Implications of Clumpiness

From extensive radiative transfer calculations we find that clumpy torus models with \No \about 5--15 dusty clouds along radial equatorial rays successfully explain AGN infrared observations. The dust has standard Galactic composition, with individual cloud optical depth \tV \about 30--100 at visual. The models naturally explain the observed behavior of the 10\mic silicate feature, in particular the lack of deep absorption features in AGN of any type. The weak 10\mic emission feature tentatively detected in type 2 QSO can be reproduced if in these sources \No drops to \about 2 or \tV exceeds \about 100. The clouds angular distribution must have a soft-edge, e.g., Gaussian profile, the radial distribution should decrease as $1/r$ or $1/r^2$. Compact tori can explain all observations, in agreement with the recent interferometric evidence that the ratio of the torus outer to inner radius is perhaps as small as \about 5--10. Clumpy torus models can produce nearly isotropic IR emission together with highly anisotropic obscuration, as required by observations. In contrast with strict variants of unification schemes where the viewing-angle uniquely determines the classification of an AGN into type 1 or 2, clumpiness implies that it is only a probabilistic effect; a source can display type 1 properties even from directions close to the equatorial plane. The fraction of obscured sources depends not only on the torus angular thickness but also on the cloud number \No. The observed decrease of this fraction at increasing luminosity can be explained with a decrease of either torus angular thickness or cloud number, but only the latter option explains also the possible emergence of a 10\mic emission feature in QSO2.

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.

preprint2005arXiv

The UV, Optical, and IR Properties of SDSS Sources Detected by GALEX

We discuss the UV, optical, and IR properties of the SDSS sources detected by GALEX as part of its All-sky Imaging Survey Early Release Observations. Virtually all of the GALEX sources in the overlap region are detected by SDSS. GALEX sources represent ~2.5% of all SDSS sources within these fields and about half are optically unresolved. Most unresolved GALEX/SDSS sources are bright blue turn-off thick disk stars and are typically detected only in the GALEX near-UV band. The remaining unresolved sources include low-redshift quasars, white dwarfs, and white dwarf/M dwarf pairs, and these dominate the optically unresolved sources detected in both GALEX bands. Almost all the resolved SDSS sources detected by GALEX are fainter than the SDSS 'main' spectroscopic limit. These sources have colors consistent with those of blue (spiral) galaxies (u-r<2.2), and most are detected in both GALEX bands. Measurements of their UV colors allow much more accurate and robust estimates of star-formation history than are possible using only SDSS data. Indeed, galaxies with the most recent (<20 Myr) star formation can be robustly selected from the GALEX data by requiring that they be brighter in the far-UV than in the near-UV band. However, older starburst galaxies have UV colors similar to AGN, and thus cannot be selected unambiguously on the basis of GALEX fluxes alone. With the aid of 2MASS data, we construct and discuss median 10 band UV-optical-IR spectral energy distributions for turn-off stars, hot white dwarfs, low-redshift quasars, and spiral and elliptical galaxies. We point out the high degree of correlation between the UV color and the contribution of the UV flux to the UV-optical-IR flux of galaxies detected by GALEX.

preprint2004arXiv

Andromeda IX: A New Dwarf Spheroidal Satellite of M31

We report the discovery of a new dwarf spheroidal satellite of M31, Andromeda IX, based on resolved stellar photometry from the Sloan Digital Sky Survey (SDSS). Using both SDSS and public archival data we have estimated its distance and other physical properties, and compared these to the properties of a previously known dwarf spheroidal companion, Andromeda V, also observed by SDSS. Andromeda IX is the lowest surface brightness galaxy found to date (mu_{V} ~ 26.8 mag arcsec^-2), and at the distance we estimate from the position of the tip of Andromeda IX's red giant branch, (m - M)_{0} ~ 24.5 (805 kpc), Andromeda IX would also be the faintest galaxy known (M_{V} ~ -8.3).

preprint2004arXiv

Spectroscopic Properties of Cool Stars in the SDSS: An Analysis of Magnetic Activity and a Search for Subdwarfs

We present a spectroscopic analysis of nearly 8000 late-type dwarfs in the Sloan Digital Sky Survey. Using the Halpha emission line as an activity indicator, we investigate the fraction of active stars as a function of spectral type and find a peak near type M8, confirming previous results. In contrast to past findings, we find that not all M7-M8 stars are active. We show that this may be a selection effect of the distance distributions of previous samples, as the active stars appear to be concentrated near the Galactic Plane. We also examine the activity strength (ratio of the luminosity emitted in Halpha to the bolometric luminosity) for each star, and find that the mean activity strength is constant over the range M0-M5 and declines at later types. The decline begins at a slightly earlier spectral type than previously found. We explore the effect that activity has on the broadband photometric colors and find no significant differences between active and inactive stars. We also carry out a search for subdwarfs using spectroscopic metallicity indicators, and find 60 subdwarf candidates. Several of these candidates are near the extreme subdwarf boundary. The spectroscopic subdwarf candidates are redder by \~0.2 magnitudes in g-r compared to disk dwarfs at the same r-i color.

preprint2003arXiv

Hdelta-Selected Galaxies in the Sloan Digital Sky Survey I: The Catalog

[Abridged] We present here a new and homogeneous sample of 3340 galaxies selected from the Sloan Digital Sky Survey (SDSS) based solely on the observed strength of their Hdelta absorption line. These galaxies are commonly known as ``post-starburst'' or ``E+A'' galaxies, and the study of these galaxies has been severely hampered by the lack of a large, statistical sample of such galaxies. In this paper, we rectify this problem by selecting a sample of galaxies which possess an absorption Hdelta equivalent width of EW(Hdelta_max) - Delta EW(Hdelta_max) > 4A from 106682 galaxies in the SDSS. We have performed extensive tests on our catalog including comparing different methodologies of measuring the Hdelta absorption and studying the effects of stellar absorption, dust extinction, emission-filling and measurement error. The measured abundance of our Hdelta-selected (HDS) galaxies is 2.6 +/- 0.1% of all galaxies within a volume-limited sample of 0.05<z<0.1 and M(r*)<-20.5, which is consistent with previous studies of such galaxies in the literature. We find that only 25 of our HDS galaxies in this volume-limited sample (3.5+/-0.7%) show no evidence for OII and Halpha emission, thus indicating that true E+A (or k+a) galaxies are extremely rare objects at low redshift, i.e., only 0.09+/-0.02% of all galaxies in this volume-limited sample are true E+A galaxies. In contrast, 89+/-5% of our HDS galaxies in the volume-limited sample have significant detections of the OII and Halpha emission lines. We find 27 galaxies in our volume-limited HDS sample that possess no detectable OII emission, but do however possess detectable Halpha emission. These galaxies may be dusty star-forming galaxies. We provide the community with this new catalog of Hdelta-selected galaxies to aid in the understanding of these galaxies.

preprint2002arXiv

Astrometric Calibration of the Sloan Digital Sky Survey

The astrometric calibration of the Sloan Digital Sky Survey is described. For point sources brighter than r ~ 20 the astrometric accuracy is 45 milliarcseconds (mas) rms per coordinate when reduced against the USNO CCD Astrograph Catalog, and 75 mas rms when reduced against Tycho-2, with an additional 20 - 30 mas systematic error in both cases. The rms errors are dominated by anomalous refraction and random errors in the primary reference catalogs. The relative astrometric accuracy between the r filter and each of the other filters (u g i z) is 25 - 35 mas rms. At the survey limit (r ~ 22), the astrometric accuracy is limited by photon statistics to approximately 100 mas rms for typical seeing. Anomalous refraction is shown to contain components correlated over two or more degrees on the sky.

preprint2001arXiv

The Ghost of Sagittarius and Lumps in the Halo of the Milky Way

We identify new structures in the halo of the Milky Way Galaxy from positions, colors and magnitudes of five million stars detected in the Sloan Digital Sky Survey. Most of these stars are within 1.26 degrees of the celestial equator. We present color-magnitude diagrams (CMDs) for stars in two previously discovered, tidally disrupted structures. The CMDs and turnoff colors are consistent with those of the Sagittarius dwarf galaxy, as had been predicted. In one direction, we are even able to detect a clump of red stars, similar to that of the Sagittarius dwarf, from stars spread across 110 square degrees of sky. Focusing on stars with the colors of F turnoff objects, we identify at least five additional overdensities of stars. Four of these may be pieces of the same halo structure, which would cover a region of the sky at least 40 degrees in diameter, at a distance of 11 kpc from the Sun (18 kpc from the center of the Galaxy). The turnoff is significantly bluer than that of thick disk stars, and closer to the Galactic plane than a power-law spheroid. We suggest two models to explain this new structure. One possibility is that this new structure could be a new dwarf satellite of the Milky Way, hidden in the Galactic plane, and in the process of being tidally disrupted. The other possibility is that it could be part of a disk-like distribution of stars which is metal-poor, with a scale height of approximately 2 kpc and a scale length of approximately 10 kpc. The fifth overdensity, which is 20 kpc away, is some distance from the Sagittarius dwarf streamer orbit and is not associated with any known structure in the Galactic plane. It is likely that there are many smaller streams of stars in the Galactic halo.

preprint1996arXiv

Dust Emission from IRC+10216

Infrared emission from the dust shell around IRC+10216 is analysed in detail, employing a self-consistent model for radiatively driven winds around late-type stars that couples the equations of motion and radiative transfer in the dust. The resulting model provides agreement with the wealth of available data, including the spectral energy distribution in the range 0.5--1000 \mic, and visibility and array observations. Previous conclusions about two dust shells, derived from modelling the data with a few single-temperature components of different radii, are not supported by our results. The extended, continuous temperature and density distributions derived from our model obviate the need for such discrete shells. The IR properties vary with the stellar phase, reflecting changes in both the dust condensation radius $r_1$ and the overall optical depth $τ$ -- as the luminosity increases from minimum to maximum, $r_1$ increases while $τ$ decreases. We find that the angular size of the dust condensation zone varies from 0.3 arcsec at minimum light to 0.5 arcsec at maximum. The shortage of flux at short wavelengths encountered in previous studies is resolved by employing a grain size distribution that includes grains larger than \about\ 0.1 \mic, required also for the visibility fits. This distribution is in agreement with the one recently proposed by Jura in a study that probed the outer regions of the envelope. Since our constraints on the size distribution mostly reflect the envelope's inner regions, the agreement of these independent studies is evidence against significant changes in grain sizes through effects like sputtering or grain growth after the initial formation at the dust condensation zone.

preprint1996arXiv

Infrared Imaging of Late-Type Stars

Infrared imaging properties of dusty winds around late-type stars are investigated in detail, employing a self-consistent model that couples the equations of motion and radiative transfer. Because of general scaling properties, the angular profiles of surface brightness are self-similar. In any given star, the profile shape is determined essentially by overall optical depth at each wavelength and it is self-similarly scaled by the size of the dust condensation zone. We find that the mid-IR is the best wavelength range to measure directly the angular size of this zone, and from {\it IRAS} data we identify the 15 best candidates for such future observations. We also show that the visibility function at short wavelengths ($\la$ 2 \mic) directly determines the scattering optical depth, and produce theoretical visibility curves for various characteristic wavelengths and the entire parameter range relevant to late-type stars. The infrared emission should display time variability because of cyclical changes in overall optical depth, reflecting luminosity-induced movement of the dust condensation point. Calculations of the wavelength dependence of photometric amplitudes and time variability of envelope sizes are in agreement with observations; envelopes are bigger and bluer at maximum light.

preprint1996arXiv

Self-similarity and scaling behavior of IR emission from radiatively heated dust: I. Theory

Dust infrared emission possesses scaling properties. Overall luminosity is never an input parameter of the radiative transfer problem, spectral shape is the only relevant property of the heating radiation when the inner boundary of the dusty region is controlled by dust sublimation. Similarly, the absolute scales of densities and distances are irrelevant; the geometry enters only through angles, relative thicknesses and aspect ratios, and the actual magnitudes of densities and distances enter only through one independent parameter, the overall optical depth. Dust properties enter only through dimensionless, normalized distributions that describe the spatial variation of density and the wavelength dependence of scattering and absorption efficiencies. Scaling enables a systematic approach to modeling and classification of IR spectra. We develop a new, fully scale-free method for solving radiative transfer, present exact numerical results, and derive approximate analytical solutions for spherical geometry, covering the entire range of parameter space relevant to observations. Scaling implies tight correlations among the SEDs of various members of the same class of sources such as young stellar objects, late-type stars, etc. In particular, all members of the same class occupy common, well defined regions in color-color diagrams. The observational data corroborate the existence of these correlations.

preprint1995arXiv

INFRARED EMISSION AND DYNAMICS OF OUTFLOWS IN LATE-TYPE STARS

The dynamical structure and infrared emission of winds around late-type stars are studied in a self-consistent model that couples the equations of motion and radiative transfer. Both the dynamics and IR spectrum of the solution are fully characterized by the flux averaged optical depth. Five types of dust grains are considered: astronomical silicate, SiC, crystalline olivine, graphite and amorphous carbon as well as mixtures. Analysis of infrared signatures provides constraints on the grain chemical composition and indications for the simultaneous existence of silicate and carbon grains. The region of parameter space where radiation pressure can support a given mass-loss rate is identified and a new method to derive mass loss rates from IR data is presented. IRAS LRS classes are associated with optical depth for various grain materials and the regions of color-color diagrams expected to be populated by late-type stars are identified. We show that cirrus emission can severely affect point source measurements to the extent that their listed IRAS long wavelength fluxes are unreliable. After accounting for cirrus contamination, essentially all IRAS point sources located in the relevant regions of the color-color diagrams can be explained as late-type stars.