Source author record

Željko Ivezić

Željko Ivezić 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

26works
7topics
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

26 published item(s)

preprint2022arXiv

MUSSES2020J: The Earliest Discovery of a Fast Blue Ultraluminous Transient at Redshift 1.063

In this Letter, we report the discovery of an ultraluminous fast-evolving transient in rest-frame UV wavelengths, MUSSES2020J, soon after its occurrence by using the Hyper Suprime-Cam (HSC) mounted on the 8.2 m Subaru telescope. The rise time of about 5 days with an extremely high UV peak luminosity shares similarities to a handful of fast blue optical transients whose peak luminosities are comparable with the most luminous supernovae while their timescales are significantly shorter (hereafter "fast blue ultraluminous transient," FBUT). In addition, MUSSES2020J is located near the center of a normal low-mass galaxy at a redshift of 1.063, suggesting a possible connection between the energy source of MUSSES2020J and the central part of the host galaxy. Possible physical mechanisms powering this extreme transient such as a wind-driven tidal disruption event and an interaction between supernova and circumstellar material are qualitatively discussed based on the first multiband early-phase light curve of FBUTs, although whether the scenarios can quantitatively explain the early photometric behavior of MUSSES2020J requires systematical theoretical investigations. Thanks to the ultrahigh luminosity in UV and blue optical wavelengths of these extreme transients, a promising number of FBUTs from the local to the high-z universe can be discovered through deep wide-field optical surveys in the near future.

preprint2022arXiv

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

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

preprint2021arXiv

Optimization of the Observing Cadence for the Rubin Observatory Legacy Survey of Space and Time: a pioneering process of community-focused experimental design

Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the Solar System, exploring the transient optical sky, and mapping the Milky Way. The survey's massive data throughput will be transformational for many other astrophysics domains and Rubin's data access policy sets the stage for a huge potential users' community. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue.

preprint2021arXiv

The LSST DESC DC2 Simulated Sky Survey

We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep (WFD) area of approximately 300 deg^2 as well as a deep drilling field (DDF) of approximately 1 deg^2. We simulate 5 years of the planned 10-year survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the dataset to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic testbed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time-domain cosmology.

preprint2020arXiv

Improving Damped Random Walk parameters for SDSS Stripe 82 Quasars with Pan-STARRS1

We use the Panoramic Survey Telescope and Rapid Response System 1 Survey (Pan-STARRS1, PS1) data to extend the Sloan Digital Sky Survey (SDSS) Stripe 82 quasar light curves. Combining PS1 and SDSS light curves provides a 15 yr baseline for 9248 quasars - 5 yr longer than prior studies that used only SDSS. We fit the light curves with the damped random walk (DRW) model model - a statistical description of their variability. We correlate the resulting DRW model parameters: asymptotic variability amplitude SF$_{\infty}$, and characteristic timescale $τ$, with quasar physical properties - black hole mass, bolometric luminosity, and redshift. Using simulated light curves, we find that a longer baseline allows us to better constrain the DRW parameters. After adding PS1 data, the variability amplitude is a stronger function of the black hole mass, and has a weaker dependence on quasar luminosity. In addition, the characteristic timescale $τ$ dependence on quasar luminosity is marginally weaker. We also make predictions for the fidelity of DRW model parameter retrieval when light curves will be further extended with Zwicky Transient Facility (ZTF) and the Rubin Observatory Legacy Survey of Space and Time (LSST) data. Finally, we show how updated DRW parameters offer an independent method of discovering changing-look quasar candidates (CLQSOs). The candidates are outliers in terms of differences in magnitude and scatter between SDSS and PS1 segments. We identify 40 objects (35 newly reported) with tenfold increase in variability timescale between SDSS and SDSS--PS1 data, which is due to a large change in brightness (over 0.5 mag) - characteristic for CLQSOs.

preprint2020arXiv

Photometric Redshifts with the LSST II: The Impact of Near-Infrared and Near-Ultraviolet Photometry

Accurate photometric redshift (photo-$z$) estimates are essential to the cosmological science goals of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). In this work we use simulated photometry for mock galaxy catalogs to explore how LSST photo-$z$ estimates can be improved by the addition of near-infrared (NIR) and/or ultraviolet (UV) photometry from the Euclid, WFIRST, and/or CASTOR space telescopes. Generally, we find that deeper optical photometry can reduce the standard deviation of the photo-$z$ estimates more than adding NIR or UV filters, but that additional filters are the only way to significantly lower the fraction of galaxies with catastrophically under- or over-estimated photo-$z$. For Euclid, we find that the addition of ${JH}$ $5σ$ photometric detections can reduce the standard deviation for galaxies with $z>1$ ($z>0.3$) by ${\sim}20\%$ (${\sim}10\%$), and the fraction of outliers by ${\sim}40\%$ (${\sim}25\%$). For WFIRST, we show how the addition of deep ${YJHK}$ photometry could reduce the standard deviation by ${\gtrsim}50\%$ at $z>1.5$ and drastically reduce the fraction of outliers to just ${\sim}2\%$ overall. For CASTOR, we find that the addition of its ${UV}$ and $u$-band photometry could reduce the standard deviation by ${\sim}30\%$ and the fraction of outliers by ${\sim}50\%$ for galaxies with $z<0.5$. We also evaluate the photo-$z$ results within sky areas that overlap with both the NIR and UV surveys, and when spectroscopic training sets built from the surveys' small-area deep fields are used.

preprint2020arXiv

Predicting the Accuracy of Asteroid Size Estimation with Data from the Rubin Observatory Legacy Survey of Space and Time

Recent work has shown that the correlation between SDSS colors and optical albedo can be used to estimate asteroid sizes from optical data alone. We revisit a correlation between SDSS colors and optical albedo for asteroids, with the albedo derived using WISE-based size estimates. Moeyens, Myhrvold \& Ivezić (2020) showed that this correlation can be used to estimate asteroid sizes with optical data alone, with a precision of about 17\% relative to WISE-based size estimates. We present here several more sophisticated data-driven models for the variation of optical albedo with colors and estimate the contribution of SDSS photometric errors to the albedo and size estimate uncertainties. We use the results of our analysis to predict that LSST data will enable asteroid size precision of about 15\% relative to WISE-based size estimates. Compared to the accuracy of WISE-based size estimates of 15-20\%, the implied accuracy of optical size estimates, in the range 21-25\%, is thus only a factor of 1.3 to 1.4 worse. This size estimation accuracy is significantly better than commonly assumed for optical data and is due to accurate and homogeneous multi-band photometry delivered by modern digital sky surveys.

preprint2020arXiv

Proper motion measurements for stars up to $100$ kpc with Subaru HSC and SDSS Stripe 82

We present proper motion measurements for more than $0.55$ million main-sequence stars, by comparing astrometric positions of matched stars between the multi-band imaging datasets from the Hyper Suprime-Cam (HSC) Survey and the SDSS Stripe 82. In doing this we use $3$ million galaxies to recalibrate the astrometry and set up a common reference frame between the two catalogues. The exquisite depth and the nearly $12$ years of time baseline between HSC and SDSS enable high-precision measurements of statistical proper motions for stars down to $i\simeq 24$. A validation of our method is demonstrated by the agreement with the $Gaia$ proper motions, to the precision better than $0.1$ mas yr$^{-1}$. To retain the precision, we make a correction of the subtle effects due to the differential chromatic refraction in the SDSS images based on the comparison with the $Gaia$ proper motions against colour of stars, which is validated using the SDSS spectroscopic quasars. Combining with the photometric distance estimates for individual stars based on the precise HSC photometry, we show a significant detection of the net proper motions for stars in each bin of distance out to $100$ kpc. The two-component tangential velocities after subtracting the apparent motions due to our own motion display rich phase-space structures including a clear signature of the Sagittarius stream in the halo region of distance range $[10,\ 35]$ kpc. We also measure the tangential velocity dispersion in the distance range $5-20$ kpc and find that the data are consistent with a constant isotropic dispersion of $80\pm 10\ {\rm km/s}$. More distant stars appear to have random motions with respect to the Galactic centre on average.

preprint2019arXiv

Morphological Star-Galaxy Separation

We discuss the statistical foundations of morphological star-galaxy separation. We show that many of the star-galaxy separation metrics in common use today (e.g. by SDSS or SExtractor) are closely related both to each other, and to the model odds ratio derived in a Bayesian framework by Sebok (1979). While the scaling of these algorithms with the noise properties of the sources varies, these differences do not strongly differentiate their performance. We construct a model of the performance of a star-galaxy separator in a realistic survey to understand the impact of observational signal-to-noise ratio (or equivalently, 5-sigma limiting depth) and seeing on classification performance. The model quantitatively demonstrates that, assuming realistic densities and angular sizes of stars and galaxies, 10% worse seeing can be compensated for by approximately 0.4 magnitudes deeper data to achieve the same star-galaxy classification performance. We discuss how to probabilistically combine multiple measurements, either of the same type (e.g., subsequent exposures), or differing types (e.g., multiple bandpasses), or differing methodologies (e.g., morphological and color-based classification). These methods are increasingly important for observations at faint magnitudes, where the rapidly rising number density of small galaxies makes star-galaxy classification a challenging problem. However, because of the significant role that the signal-to-noise ratio plays in resolving small galaxies, surveys with large-aperture telescopes, such as LSST, will continue to see improving star-galaxy separation as they push to these fainter magnitudes.

preprint2015arXiv

Quasar Classification Using Color and Variability

We conduct a pilot investigation to determine the optimal combination of color and variability information to identify quasars in current and future multi-epoch optical surveys. We use a Bayesian quasar selection algorithm (Richards et al. 2004) to identify 35,820 type 1 quasar candidates in a 239 square degree field of the Sloan Digital Sky Survey (SDSS) Stripe 82, using a combination of optical photometry and variability. Color analysis is performed on 5-band single- and multi-epoch SDSS optical photometry to a depth of r ~22.4. From these data, variability parameters are calculated by fitting the structure function of each object in each band with a power law model using 10 to >100 observations over timescales from ~1 day to ~8 years. Selection was based on a training sample of 13,221 spectroscopically-confirmed type-1 quasars, largely from the SDSS. Using variability alone, colors alone, and combining variability and colors we achieve 91%, 93%, and 97% quasar completeness and 98%, 98%, and 97% efficiency respectively, with particular improvement in the selection of quasars at 2.7<z<3.5 where quasars and stars have similar optical colors. The 22,867 quasar candidates that are not spectroscopically confirmed reach a depth of i ~22.0; 21,876 (95.7%) are dimmer than coadded i-band magnitude of 19.9, the cut off for spectroscopic follow-up for SDSS on Stripe 82. Brighter than 19.9, we find 5.7% more quasar candidates without confirming spectra in sky regions otherwise considered complete. The resulting quasar sample has sufficient purity (and statistically correctable incompleteness) to produce a luminosity function comparable to those determined by spectroscopic investigations. We discuss improvements that can be made to the process in preparation for performing similar photometric selection and science on data from post-SDSS sky surveys.

preprint2015arXiv

Revealing the Nature of Extreme Coronal-line Emitter SDSS J095209.56+214313.3

Extreme coronal-line emitter (ECLE) SDSSJ095209.56+214313.3, known by its strong, fading, high ionization lines, has been a long standing candidate for a tidal disruption event, however a supernova origin has not yet been ruled out. Here we add several new pieces of information to the puzzle of the nature of the transient that powered its variable coronal lines: 1) an optical light curve from the Lincoln Near Earth Asteroid Research (LINEAR) survey that serendipitously catches the optical flare, and 2) late-time observations of the host galaxy with the Swift Ultraviolet and Optical Telescope (UVOT) and X-ray telescope (XRT) and the ground-based Mercator telescope. The well-sampled, $\sim10$-year long, unfiltered LINEAR light curve constrains the onset of the flare to a precision of $\pm5$ days and enables us to place a lower limit on the peak optical magnitude. Difference imaging allows us to estimate the location of the flare in proximity of the host galaxy core. Comparison of the \textsl{GALEX} data (early 2006) with the recently acquired Swift UVOT (June 2015) and Mercator observations (April 2015) demonstrate a decrease in the UV flux over a $\sim 10$ year period, confirming that the flare was UV-bright. The long-lived UV-bright emission, detected 1.8 rest-frame years after the start of the flare, strongly disfavors a SN origin. These new data allow us to conclude that the flare was indeed powered by the tidal disruption of a star by a supermassive black hole and that TDEs are in fact capable of powering the enigmatic class of ECLEs.

preprint2015arXiv

The Gaia-LSST Synergy

We discuss the synergy of Gaia and the Large Synoptic Survey Telescope (LSST) in the context of Milky Way studies. LSST can be thought of as Gaia's deep complement because the two surveys will deliver trigonometric parallax, proper-motion, and photometric measurements with similar uncertainties at Gaia's faint end at $r=20$, and LSST will extend these measurements to a limit about five magnitudes fainter. We also point out that users of Gaia data will have developed data analysis skills required to benefit from LSST data, and provide detailed information about how international participants can join LSST.

preprint2013arXiv

Exploring the Variable Sky with LINEAR. II. Halo Structure and Substructure Traced by RR Lyrae Stars to 30 kpc

We present a sample of ~5,000 RR Lyrae stars selected from the recalibrated LINEAR dataset and detected at heliocentric distances between 5 kpc and 30 kpc over ~8,000 deg^2 of sky. The coordinates and light curve properties, such as period and Oosterhoff type, are made publicly available. We find evidence for the Oosterhoff dichotomy among field RR Lyrae stars, with the ratio of the type II and I subsamples of about 1:4. The number density distribution of halo RRab stars as a function of galactocentric distance can be described as an oblate ellipsoid with the axis ratio q=0.63 and with either a single or a double power law with a power-law index in the range -2 to -3. Using a group-finding algorithm EnLink, we detected seven candidate halo groups, only one of which is statistically spurious. Three of these groups are near globular clusters (M53/NGC 5053, M3, M13), and one is near a known halo substructure (Virgo Stellar Stream); the remaining three groups do not seem to be near any known halo substructures or globular clusters, and seem to have a higher ratio of Oosterhoff type II to Oosterhoff type I RRab stars than what is found in the halo. The extended morphology and the position (outside the tidal radius) of some of the groups near globular clusters is suggestive of tidal streams possibly originating from globular clusters. Spectroscopic followup of detected halo groups is encouraged.

preprint2013arXiv

Exploring the Variable Sky with LINEAR. III. Classification of Periodic Light Curves

We describe the construction of a highly reliable sample of approximately 7,000 optically faint periodic variable stars with light curves obtained by the asteroid survey LINEAR across 10,000 sq.deg of northern sky. Majority of these variables have not been cataloged yet. The sample flux limit is several magnitudes fainter than for most other wide-angle surveys; the photometric errors range from ~0.03 mag at $r=15$ to ~0.20 mag at r=18. Light curves include on average 250 data points, collected over about a decade. Using SDSS-based photometric recalibration of the LINEAR data for about 25 million objects, we selected ~200,000 most probable candidate variables and visually confirmed and classified approximately 7,000 periodic variables using phased light curves. The reliability and uniformity of visual classification across eight human classifiers was calibrated and tested using a SDSS Stripe 82 region variable star catalog, and verified using an unsupervised machine learning approach. The resulting sample of periodic LINEAR variables is dominated by 3,900 RR Lyrae stars and 2,700 eclipsing binary stars of all subtypes, and includes small fractions of relatively rare populations such as asymptotic giant branch stars and SX Phoenicis stars. We discuss the distribution of these mostly uncataloged variables in various diagrams constructed with optical-to-infrared SDSS, 2MASS and WISE photometry, and with LINEAR light curve features. An interesting side result is a robust and precise quantitative description of a strong correlation between the light-curve period and color/spectral type for close and contact eclipsing binary stars. These large samples of robustly classified variable stars will enable detailed statistical studies of the Galactic structure and physics of binary and other stars, and we make them publicly available.

preprint2013arXiv

New Organizations to Support Astroinformatics and Astrostatistics

In the past two years, the environment within which astronomers conduct their data analysis and management has rapidly changed. Working Groups associated with international societies and Big Data projects have emerged to support and stimulate the new fields of astroinformatics and astrostatistics. Sponsoring societies include the Intenational Statistical Institute, International Astronomical Union, American Astronomical Society, and Large Synoptic Survey Telescope project. They enthusiastically support cross-disciplinary activities where the advanced capabilities of computer science, statistics and related fields of applied mathematics are applied to advance research on planets, stars, galaxies and the Universe. The ADASS community is encouraged to join these organizations and to explore and engage in their public communication Web site, the Astrostatistics and Astroinformatics Portal (http://asaip.psu.edu).

preprint2013arXiv

Variability-based AGN selection using image subtraction in the SDSS and LSST era

With upcoming all sky surveys such as LSST poised to generate a deep digital movie of the optical sky, variability-based AGN selection will enable the construction of highly-complete catalogs with minimum contamination. In this study, we generate $g$-band difference images and construct light curves for QSO/AGN candidates listed in SDSS Stripe 82 public catalogs compiled from different methods, including spectroscopy, optical colors, variability, and X-ray detection. Image differencing excels at identifying variable sources embedded in complex or blended emission regions such as Type II AGNs and other low-luminosity AGNs that may be omitted from traditional photometric or spectroscopic catalogs. To separate QSOs/AGNs from other sources using our difference image light curves, we explore several light curve statistics and parameterize optical variability by the characteristic damping timescale ($τ$) and variability amplitude. By virtue of distinguishable variability parameters of AGNs, we are able to select them with high completeness of 93.4% and efficiency (i.e., purity) of 71.3%. Based on optical variability, we also select highly variable blazar candidates, whose infrared colors are consistent with known blazars. One third of them are also radio detected. With the X-ray selected AGN candidates, we probe the optical variability of X-ray detected optically-extended sources using their difference image light curves for the first time. A combination of optical variability and X-ray detection enables us to select various types of host-dominated AGNs. Contrary to the AGN unification model prediction, two Type II AGN candidates (out of 6) show detectable variability on long-term timescales like typical Type I AGNs. This study will provide a baseline for future optical variability studies of extended sources.

preprint2012arXiv

The Stellar Metallicity Distribution Function of the Galactic Halo from SDSS Photometry

We explore the stellar metallicity distribution function of the Galactic halo based on SDSS ugriz photometry. A set of stellar isochrones is calibrated using observations of several star clusters and validated by comparisons with medium-resolution spectroscopic values over a wide range of metal abundance. We estimate distances and metallicities for individual main-sequence stars in the multiply scanned SDSS Stripe 82, at heliocentric distances in the range 5 - 8 kpc and |b| > 35 deg, and find that the in situ photometric metallicity distribution has a shape that matches that of the kinematically-selected local halo stars from Ryan & Norris. We also examine independent kinematic information from proper-motion measurements for high Galactic latitude stars in our sample. We find that stars with retrograde rotation in the rest frame of the Galaxy are generally more metal poor than those exhibiting prograde rotation, which is consistent with earlier arguments by Carollo et al. that the halo system comprises at least two spatially overlapping components with differing metallicity, kinematics, and spatial distributions. The observed photometric metallicity distribution and that of Ryan & Norris can be described by a simple chemical evolution model by Hartwick (or by a single Gaussian distribution); however, the suggestive metallicity-kinematic correlation contradicts the basic assumption in this model that the Milky Way halo consists primarily of a single stellar population. When the observed metallicity distribution is deconvolved using two Gaussian components with peaks at [Fe/H] ~ -1.7 and -2.3, the metal-poor component accounts for ~20% - 35% of the entire halo population in this distance range.

preprint2011arXiv

Correlations of Quasar Optical Spectra with Radio Morphology

(Abridged) With the largest homogeneous quasar sample with optical spectra and robust radio morphology classifications assembled to date, we investigate quasar radio and optical properties with unprecedented statistical power. The data consist of 4714 radio quasars from FIRST with S_{20}>2mJy and SDSS spectra. Radio morphology classes include core-only (core), core-lobe (lobe), core-jet (jet), lobe-core-lobe (triple), and double-lobe. We examine the optical colors of radio-morphology subsamples and find that radio quasars with core emission unresolved by FIRST (on 5" scale) have a redder color distribution than radio-quiet quasars (S_{20}<1mJy); other classes of radio quasars have optical color distributions similar to the radio-quiet quasars. This analysis also suggests that optical colors of z<2.7 SDSS quasars are not strongly (<0.1mag) biased blue. We show that the radio core-to-lobe flux density ratio (R) and the radio-to-optical (i-band) ratio of the quasar core (RI) are correlated, suggesting that both parameters are indicative of line-of-sight orientation. We investigate spectral line equivalent widths as a function of R and RI, including the [OIII] narrow line doublet and the CIV λ1549 and MgII λ2799 broad lines. We find that the rest equivalent widths (EWs) of the broad lines correlate positively with RI at the 4-8σ level. But we find no strong dependence of EW on R, in contrast to previous results. One interpretation is that EWs increase as the line-of-sight angle to the radio-jet axis decreases. These results are in stark contrast to commonly accepted orientation-based theories, which suggest that continuum emission should increase as the angle to the radio-jet axis decreases, resulting in smaller EW of emission lines (assumed isotropic). Finally, we find that the Baldwin effect in our sample does not depend on quasar radio morphology.

preprint2011arXiv

Exploring the Variable Sky with LINEAR. I. Photometric Recalibration with SDSS

We describe photometric recalibration of data obtained by the asteroid survey LINEAR. Although LINEAR was designed for astrometric discovery of moving objects, the dataset described here contains over 5 billion photometric measurements for about 25 million objects, mostly stars. We use SDSS data from the overlapping ~10,000 deg^2 of sky to recalibrate LINEAR photometry, and achieve errors of 0.03 mag for sources not limited by photon statistics, with errors of 0.2 mag at r~18. With its 200 observations per object on average, LINEAR data provide time domain information for the brightest 4 magnitudes of SDSS survey. At the same time, LINEAR extends the deepest similar wide-area variability survey, the Northern Sky Variability Survey, by 3 mag. We briefly discuss the properties of about 7,000 visually confirmed periodic variables, dominated by roughly equal fractions of RR Lyrae stars and eclipsing binary stars, and analyze their distribution in optical and infra-red color-color diagrams. The LINEAR dataset is publicly available from the SkyDOT website (http://skydot.lanl.gov).

preprint2011arXiv

Milky Way Tomography IV: Dissecting Dust

We use SDSS photometry of 73 million stars to simultaneously obtain best-fit main-sequence stellar energy distribution (SED) and amount of dust extinction along the line of sight towards each star. Using a subsample of 23 million stars with 2MASS photometry, whose addition enables more robust results, we show that SDSS photometry alone is sufficient to break degeneracies between intrinsic stellar color and dust amount when the shape of extinction curve is fixed. When using both SDSS and 2MASS photometry, the ratio of the total to selective absorption, $R_V$, can be determined with an uncertainty of about 0.1 for most stars in high-extinction regions. These fits enable detailed studies of the dust properties and its spatial distribution, and of the stellar spatial distribution at low Galactic latitudes. Our results are in good agreement with the extinction normalization given by the Schlegel et al. (1998, SFD) dust maps at high northern Galactic latitudes, but indicate that the SFD extinction map appears to be consistently overestimated by about 20% in the southern sky, in agreement with Schlafly et al. (2010). The constraints on the shape of the dust extinction curve across the SDSS and 2MASS bandpasses support the models by Fitzpatrick (1999) and Cardelli et al. (1989). For the latter, we find an $R_V=3.0\pm0.1$(random) $\pm0.1$(systematic) over most of the high-latitude sky. At low Galactic latitudes (|b|<5), we demonstrate that the SFD map cannot be reliably used to correct for extinction as most stars are embedded in dust, rather than behind it. We introduce a method for efficient selection of candidate red giant stars in the disk, dubbed "dusty parallax relation", which utilizes a correlation between distance and the extinction along the line of sight. We make these best-fit parameters, as well as all the input SDSS and 2MASS data, publicly available in a user-friendly format.

preprint2010arXiv

The Blue Tip of the Stellar Locus: Measuring Reddening with the SDSS

We present measurements of reddening due to dust using the colors of stars in the Sloan Digital Sky Survey (SDSS). We measure the color of main sequence turn-off stars by finding the "blue tip" of the stellar locus: the prominent blue edge in the distribution of stellar colors. The method is sensitive to color changes of order 18, 12, 7, and 8 mmag of reddening in the colors u-g, g-r, r-i, and i-z, respectively, in regions measuring 90' by 14'. We present maps of the blue tip colors in each of these bands over the entire SDSS footprint, including the new dusty southern Galactic cap data provided by the SDSS-III. The results disfavor the best fit O'Donnell (1994) and Cardelli et al. (1989) reddening laws, but are well described by a Fitzpatrick (1999) reddening law with R_V = 3.1. The SFD dust map is found to trace the dust well, but overestimates reddening by factors of 1.4, 1.0, 1.2, and 1.4 in u-g, g-r, r-i, and i-z, largely due to the adopted reddening law. In select dusty regions of the sky, we find evidence for problems in the SFD temperature correction. A dust map normalization difference of 15% between the Galactic north and south sky may be due to these dust temperature errors.

preprint2009arXiv

Spectroscopic Confirmation of the Pisces Overdensity

We present spectroscopic confirmation of the "Pisces Overdensity", also known as "Structure J", a photometric overdensity of RR Lyrae stars discovered by the Sloan Digital Sky Survey (SDSS) at an estimated photometric distance of ~85kpc. We measure radial velocities for 8 RR Lyrae stars within Pisces. We find that 5 of the 8 stars have heliocentric radial velocities within a narrow range of -87 km/s < v < -67 km/s, suggesting that the photometric overdensity is mainly due to a physically associated system, probably a dwarf galaxy or a disrupted galaxy. Two of the remaining 3 stars differ from one another by only 9 km/s, but it would be premature to identify them as a second system.

preprint2005arXiv

Active Galactic Nuclei in the Sloan Digital Sky Survey: I. Sample Selection

We have compiled a large sample of low-redshift active galactic nuclei (AGN) identified via their emission line characteristics from the spectroscopic data of the Sloan Digital Sky Survey. Since emission lines are often contaminated by stellar absorption lines, we developed an objective and efficient method of subtracting the stellar continuum from every galaxy spectrum before making emission line measurements. The distribution of the measured H$α$ Full Width at Half Maxima values of emission line galaxies is strongly bimodal, with two populations separated at about 1,200km s$^{-1}$. This feature provides a natural separation between narrow-line and broad-line AGN. The narrow-line AGN are identified using standard emission line ratio diagnostic diagrams. 1,317 broad-line and 3,074 narrow-line AGN are identified from about 100,000 galaxy spectra selected over 1151 square degrees. This sample is used in a companion paper to determine the emission-line luminosity function of AGN.

preprint2005arXiv

Active Galactic Nuclei in the Sloan Digital Sky Survey: II. Emission-Line Luminosity Function

The emission line luminosity function of active galactic nuclei (AGN) is measured from about 3000 AGN included in the main galaxy sample of the Sloan Digital Sky Survey within a redshift range of $0<z<0.15$. The $\Ha$ and [OIII]$λ5007$ luminosity functions for Seyferts cover luminosity range of $10^{5-9}$$L_\odot$ in H$α$ and the shapes are well fit by broken power laws, without a turnover at fainter nuclear luminosities. Assuming a universal conversion from emission line strength to continuum luminosity, the inferred B band magnitude luminosity function is comparable both to the AGN luminosity function of previous studies and to the low redshift quasar luminosity function derived from the 2dF redshift survey. The inferred AGN number density is approximately 1/5 of all galaxies and about $6\times 10^{-3}$ of the total light of galaxies in the $r$-band comes from the nuclear activity. The numbers of Seyfert 1s and Seyfert 2s are comparable at low luminosity, while at high luminosity, Seyfert 1s outnumber Seyfert 2s by a factor of 2-4. In making the luminosity function measurements, we assumed that the nuclear luminosity is independent of the host galaxy luminosity, an assumption we test {\it a posteriori}, and show to be consistent with the data. Given the relationship between black hole mass and host galaxy bulge luminosity, the lack of correlation between nuclear and host luminosity suggests that the main variable that determines the AGN luminosity is the Eddington ratio, not the black hole mass. This appears to be different from luminous quasars, which are most likely to be shining near the Eddington limit.

preprint2005arXiv

The Colours of Elliptical galaxies

We have compiled a sample of 2728 nearby ($z<0.08$) elliptical galaxies with photometry in the $g,r,i,z$ bands from the Sloan Digital Sky Survey (SDSS) and $J,H,K$ photometry from the Two Micron All Sky Survey (2MASS). Stellar masses, stellar velocity dispersions and structural parameters such as sizes and surface mass densities are also available for these objects. In order to correct the aperture mismatch between SDSS and 2MASS, we correct the SDSS magnitudes to the isophotal circular radius where the 2MASS magnitudes are measured. We compare the correlations between optical, optical-infrared and infrared colours and galaxy luminosity, stellar mass, velocity dispersion and surface mass density. We find that all galaxy colours correlate more strongly with stellar mass and velocity dispersion than with any other structural parameter. The dispersion about these two relations is also smaller. We also study the correlations between a variety of stellar absorption line indices and the same set of galaxy parameters and we reach very similar conclusions. Finally, we analyze correlations between absorption line indices and colour. Our results suggest that the optical colours of elliptical galaxies are sensitive to a combination of age, metallicity and $α$-enhancement, while the optical-infrared colours are sensitive to metallicity and to $α$-enhancement, but are somewhat less sensitive to age.