Source author record

B. Sesar

B. Sesar 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

22works
5topics
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

22 published item(s)

preprint2020arXiv

A Larger Extent for the Ophiuchus Stream

We present new kinematic data for the Ophiuchus stellar stream. Spectra have been taken of member candidates at the MMT telescope using Hectospec, Hectochelle and Binospec, which provide more than 1800 new velocities. Combined with proper motion measurements of stars in the field by the Gaia - DR2 catalog, we have derived stream membership probabilities, resulting in the detection of more than 200 likely members. These data show the stream extends to more than three times the length shown in the discovery data. A spur to the main stream is also detected. The high resolution spectra allow us to resolve the stellar velocity dispersion, found to be $1.6 \pm 0.3 $ km/sec.

preprint2016arXiv

A Global Astrometric Solution for Pan-STARRS referenced to ICRF2

We describe development and application of a Global Astrometric Solution (GAS) to the problem of Pan-STARRS1 (PS1) astrometry. Current PS1 astrometry is based on differential astrometric measurements using 2MASS reference stars, thus PS1 astrometry inherits the errors of the 2MASS catalog. The GAS, based on a single, least squares adjustment to approximately 750k grid stars using over 3000 extragalactic objects as reference objects, avoids this catalog-to-catalog propagation of errors to a great extent. The GAS uses a relatively small number of Quasi-Stellar Objects (QSOs, or distant AGN) with very accurate (<1 mas) radio positions, referenced to the ICRF2. These QSOs provide a hard constraint in the global least squares adjustment. Solving such a system provides absolute astrometry for all the stars simultaneously. The concept is much cleaner than conventional astrometry but is not easy to perform for large catalogs. In this paper we describe our method and its application to Pan-STARRS1 data. We show that large-scale systematic errors are easily corrected but our solution residuals for position (~60 mas) are still larger than expected based on simulations (~10 mas). We provide a likely explanation for the reason the small-scale residual errors are not corrected in our solution as would be expected.

preprint2016arXiv

iPTF Search for an Optical Counterpart to Gravitational Wave Trigger GW150914

The intermediate Palomar Transient Factory (iPTF) autonomously responded to and promptly tiled the error region of the first gravitational wave event GW150914 to search for an optical counterpart. Only a small fraction of the total localized region was immediately visible in the Northern night sky, due both to sun-angle and elevation constraints. Here, we report on the transient candidates identified and rapid follow-up undertaken to determine the nature of each candidate. Even in the small area imaged of 126 sq deg, after extensive filtering, 8 candidates were deemed worthy of additional follow-up. Within two hours, all 8 were spectroscopically classified by the Keck II telescope. Curiously, even though such events are rare, one of our candidates was a superluminous supernova. We obtained radio data with the Jansky Very Large Array and X-ray follow-up with the Swift satellite for this transient. None of our candidates appear to be associated with the gravitational wave trigger, which is unsurprising given that GW150914 came from the merger of two stellar-mass black holes. This end-to-end discovery and follow-up campaign bodes well for future searches in this post-detection era of gravitational waves.

preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

preprint2015arXiv

Evidence of Fanning in the Ophiuchus Stream

The Ophiuchus stellar stream presents a dynamical puzzle: its old stellar populations ($\sim 12$ Gyr) cannot be reconciled with (1) its orbit in a simple model for the Milky Way potential and (2) its short angular extent, both of which imply that the observed stream formed within the last $<1$ Gyr. Recent theoretical work has shown that streams on chaotic orbits may abruptly fan out near their apparent ends; stars in these fans are dispersed in both position and velocity and may be difficult to associate with the stream. Here we present the first evidence of such stream-fanning in the Ophiuchus stream, traced by four blue horizontal branch (BHB) stars beyond the apparent end of the stream. These stars stand out from the background by their high velocities ($v_{\rm los} > 230$ km s$^{-1}$) against $\sim 40$ other stars: their velocities are comparable to those of the stream, but would be exceptional if they were unrelated halo stars. Their positions and velocities are, however, inconsistent with simple extrapolation of the observed cold, high-density portion of the stream. These observations suggest that stream-fanning may be a real, observable effect and, therefore, that Ophiuchus may be on a chaotic orbit. They also show that the Ophiuchus stream is more extended and hence dynamically older than previously thought, easing the stellar population vs. dynamical age tension.

preprint2015arXiv

The fastest unbound star in our Galaxy ejected by a thermonuclear supernova

Hypervelocity stars (HVS) travel with velocities so high, that they exceed the escape velocity of the Galaxy. Several acceleration mechanisms have been discussed. Only one HVS (US 708, HVS 2) is a compact helium star. Here we present a spectroscopic and kinematic analysis of US\,708. Travelling with a velocity of $\sim1200\,{\rm km\,s^{-1}}$, it is the fastest unbound star in our Galaxy. In reconstructing its trajectory, the Galactic center becomes very unlikely as an origin, which is hardly consistent with the most favored ejection mechanism for the other HVS. Furthermore, we discovered US\,708 to be a fast rotator. According to our binary evolution model it was spun-up by tidal interaction in a close binary and is likely to be the ejected donor remnant of a thermonuclear supernova.

preprint2015arXiv

The Nature and Orbit of the Ophiuchus Stream

The Ophiuchus stream is a recently discovered stellar tidal stream in the Milky Way. We present high-quality spectroscopic data for 14 stream member stars obtained using the Keck and MMT telescopes. We confirm the stream as a fast moving ($v_{los}\sim290$ km s$^{-1}$), kinematically cold group ($σ_{v_{los}}\lesssim1$ km s$^{-1}$) of $α$-enhanced and metal-poor stars (${\rm [α/Fe]\sim0.4}$ dex, ${\rm [Fe/H]\sim-2.0}$ dex). Using a probabilistic technique, we model the stream simultaneously in line-of-sight velocity, color-magnitude, coordinate, and proper motion space, and so determine its distribution in 6D phase-space. We find that that the stream extends in distance from 7.5 to 9 kpc from the Sun; it is 50 times longer than wide, merely appearing highly foreshortened in projection. The analysis of the stellar population contained in the stream suggests that it is $\sim12$ Gyr old, and that its initial stellar mass was $\sim2\times10^4$ $M_{\odot}$ (or at least $\gtrsim7\times10^3$ $M_{\odot}$). Assuming a fiducial Milky Way potential, we fit an orbit to the stream which matches the observed phase-space distribution, except for some tension in the proper motions: the stream has an orbital period of $\sim350$ Myr, and is on a fairly eccentric orbit ($e\sim0.66$) with a pericenter of $\sim3.5$ kpc and an apocenter of $\sim17$ kpc. The phase-space structure and stellar population of the stream show that its progenitor must have been a globular cluster that was disrupted only $\sim240$ Myr ago. We do not detect any significant overdensity of stars along the stream that would indicate the presence of a progenitor, and conclude that the stream is all that is left of the progenitor.

preprint2014arXiv

A multi-wavelength investigation of the radio-loud supernova PTF11qcj and its circumstellar environment

We present the discovery, classification, and extensive panchromatic (from radio to X-ray) follow-up observations of PTF11qcj, a supernova discovered by the Palomar Transient Factory. PTF11qcj is located at a distance of dL ~ 124 Mpc. Our observations with the Karl G. Jansky Very Large Array show that this event is radio-loud: PTF11qcj reached a radio peak luminosity comparable to that of the famous gamma-ray-burst-associated supernova 1998bw (L_{5GHz} ~ 10^{29} erg/s/Hz). PTF11qcj is also detected in X-rays with the Chandra observatory, and in the infrared band with Spitzer. Our multi-wavelength analysis probes the supernova interaction with circumstellar material. The radio observations suggest a progenitor mass-loss rate of ~10^{-4} Msun/yr x (v_w/1000 km/s), and a velocity of ~(0.3-0.5)c for the fastest moving ejecta (at ~10d after explosion). However, these estimates are derived assuming the simplest model of supernova ejecta interacting with a smooth circumstellar material characterized by radial power-law density profile, and do not account for possible inhomogeneities in the medium and asphericity of the explosion. The radio light curve shows deviations from such a simple model, as well as a re-brightening at late times. The X-ray flux from PTF11qcj is compatible with the high-frequency extrapolation of the radio synchrotron emission (within the large uncertainties). An IR light echo from pre-existing dust is in agreement with our infrared data. Our analysis of pre-explosion data from the Palomar Transient Factory suggests that a precursor eruption of absolute magnitude M_r ~ -13 mag may have occurred ~ 2.5 yr prior to the supernova explosion. Based on our panchromatic follow-up campaign, we conclude that PTF11qcj fits the expectations from the explosion of a Wolf-Rayet star. Precursor eruptions may be a feature characterizing the final pre-explosion evolution of such stars.

preprint2014arXiv

Connection between optical and gamma-ray variability in blazars

We use optical data from the Palomar Transient Factory (PTF) and the Catalina Real-Time Transient Survey (CRTS) to study the variability of gamma-ray detected and non-detected objects in a large population of active galactic nuclei (AGN) selected from the Candidate Gamma-Ray Blazar Survey and Fermi Gamma-Ray Space Telescope catalogs. Our samples include 714 sources with PTF data and 1244 sources with CRTS data. We calculate the intrinsic modulation index to quantify the optical variability amplitude in these samples. We find the gamma-ray detected objects to be more variable than the non-detected ones. The flat spectrum radio quasars (FSRQs) are more variable than the BL Lac objects in our sample, but the significance of the difference depends on the sample used. When dividing the objects based on their synchrotron peak frequency, we find the low synchrotron peaked (LSP) objects to be significantly more variable than the high synchrotron peaked (HSP) ones, explaining the difference between the FSRQs and BL Lacs. This could be due to the LSPs being observed near their electron energy peak, while in the HSPs the emission is caused by lower energy electrons, which cool more slowly. We also find a significant correlation between the optical and gamma-ray fluxes that is stronger in the HSP BL Lacs than in the FSRQs. The FSRQs in our sample are also more Compton dominated than the HSP BL Lacs. These findings are consistent with models where the gamma-ray emission of HSP objects is produced by the synchrotron self-Compton mechanism, while the LSP objects need an additional external Compton component that increases the scatter in the flux-flux correlation.

preprint2014arXiv

Stacking the Invisibles: A Guided Search for Low-Luminosity Milky Way Satellites

Almost every known low-luminosity Milky Way dwarf spheroidal (dSph) satellite galaxy contains at least one RR Lyrae star. Assuming that a fraction of distant (60 < d_{helio} < 100 kpc) Galactic halo RR Lyrae stars are members of yet to be discovered low-luminosity dSph galaxies, we perform a {\em guided} search for these low-luminosity dSph galaxies. In order to detect the presence of dSph galaxies, we combine stars selected from more than 123 sightlines centered on RR Lyrae stars identified by the Palomar Transient Factory. We find that this method is sensitive enough to detect the presence of Segue 1-like galaxies (M_V= -1.5^{+0.6}_{-0.8}, r_h=30 pc) even if only ~20 sightlines were occupied by such dSph galaxies. Yet, when our method is applied to the SDSS DR10 imaging catalog, no signal is detected. An application of our method to sightlines occupied by pairs of close (<200 pc) horizontal branch stars, also did not yield a detection. Thus, we place upper limits on the number of low-luminosity dSph galaxies with half-light radii from 30 pc to 120 pc, and in the probed volume of the halo. Stronger constraints on the luminosity function may be obtained by applying our method to sightlines centered on RR Lyrae stars selected from the Pan-STARRS1 survey, and eventually, from LSST. In the Appendix, we present spectroscopic observations of an RRab star in the Boötes 3 dSph and a light curve of an RRab star near the Boötes 2 dSph.

preprint2013arXiv

High-precision 2MASS JHKs light curves and other data for RR Lyrae star SDSS J015450+001501: strong constraints for non-linear pulsation models

We present and discuss an extensive data set for the non-Blazhko ab-type RR Lyrae star SDSSJ015450+001501, including optical SDSS ugriz light curves and spectroscopic data, LINEAR and CSS unfiltered optical light curves, and infrared 2MASS JHKs and WISE W1 and W2 light curves. Most notably, light curves obtained by 2MASS include close to 9000 photometric measures collected over 3.3 years and provide exceedingly precise view of near-IR variability. These data demonstrate that static atmosphere models are insufficient to explain multi-band photometric light curve behavior and present strong constraints for non-linear pulsation models for RR Lyrae stars. It is a challenge to modelers to produce theoretical light curves that can explain data presented here, which we make publicly available.

preprint2012arXiv

Search for high-amplitude Delta Scuti and RR Lyrae stars in Sloan Digital Sky Survey Stripe 82 using principal component analysis

We propose a robust principal component analysis (PCA) framework for the exploitation of multi-band photometric measurements in large surveys. Period search results are improved using the time series of the first principal component due to its optimized signal-to-noise ratio.The presence of correlated excess variations in the multivariate time series enables the detection of weaker variability. Furthermore, the direction of the largest variance differs for certain types of variable stars. This can be used as an efficient attribute for classification. The application of the method to a subsample of Sloan Digital Sky Survey Stripe 82 data yielded 132 high-amplitude Delta Scuti variables. We found also 129 new RR Lyrae variables, complementary to the catalogue of Sesar et al., 2010, extending the halo area mapped by Stripe 82 RR Lyrae stars towards the Galactic bulge. The sample comprises also 25 multiperiodic or Blazhko RR Lyrae stars.

preprint2012arXiv

The complex case of V445 Lyr observed with Kepler: Two Blazhko modulations, a non-radial mode, possible triple mode RR Lyrae pulsation, and more

Rapid and strong changes in the Blazhko modulation of RR Lyrae stars, as they have recently been detected in high precision satellite data, have become a crucial topic in finding an explanation of the long-standing mystery of the Blazhko effect. We present here an analysis of the most extreme case detected so far, the RRab star V445 Lyr (KIC 6186029) which was observed with the Kepler space mission. V445 Lyr shows very strong cycle-to-cycle changes in its Blazhko modulation, which are caused both by a secondary long-term modulation period as well as irregular variations. In addition to the complex Blazhko modulation, V445 Lyr also shows a rich spectrum of additional peaks in the frequency range between the fundamental pulsation and the first harmonic. Among those peaks, the second radial overtone could be identified, which, combined with a metallicity estimate of [Fe/H]=-2.0 dex from spectroscopy, allowed to constrain the mass (0.55-0.65 M_sun) and luminosity (40-50 L_sun) of V445 Lyr through theoretical Petersen diagrams. A non-radial mode as well as possibly the first overtone are also excited. Furthermore, V445 Lyr shows signs of the period doubling phenomenon and a long term period change. A detailed Fourier analysis along with a study of the O-C variation of V445 Lyr is presented, and the origin of the additional peaks and possible causes of the changes in the Blazhko modulation are discussed. The results are then put into context with those of the only other star with a variable Blazhko effect for which a long enough set of high precision continuous satellite data has been published so far, the CoRoT star 105288363.

preprint2012arXiv

The Palomar Transient Factory photometric catalog 1.0

We construct a photometrically calibrated catalog of non-variable sources from the Palomar Transient Factory (PTF) observations. The first version of this catalog presented here, the PTF photometric catalog 1.0, contains calibrated R_PTF-filter magnitudes for about 21 million sources brighter than magnitude 19, over an area of about 11233 deg^2. The magnitudes are provided in the PTF photometric system, and the color of a source is required in order to convert these magnitudes into other magnitude systems. We estimate that the magnitudes in this catalog have typical accuracy of about 0.02 mag with respect to magnitudes from the Sloan Digital Sky Survey. The median repeatability of our catalog's magnitudes for stars between 15 and 16 mag, is about 0.01 mag, and it is better than 0.03 mag for 95% of the sources in this magnitude range. The main goal of this catalog is to provide reference magnitudes for photometric calibration of visible light observations. Subsequent versions of this catalog, which will be published incrementally online, will be extended to a larger sky area and will also include g_PTF-filter magnitudes, as well as variability and proper motion information.

preprint2012arXiv

Two Distant Halo Velocity Groups Discovered by the Palomar Transient Factory

We report the discovery of two new halo velocity groups (Cancer groups A and B) traced by 8 distant RR Lyrae stars and observed by the Palomar Transient Factory (PTF) survey at R.A.~129 deg, Dec~20 deg (l~205 deg, b~32 deg). Located at 92 kpc from the Galactic center (86 kpc from the Sun), these are some of the most distant substructures in the Galactic halo known to date. Follow-up spectroscopic observations with the Palomar Observatory 5.1-m Hale telescope and W. M. Keck Observatory 10-m Keck I telescope indicate that the two groups are moving away from the Galaxy at v_{gsr} = 78.0+-5.6 km/s (Cancer group A) and v_{gsr} = 16.3+-7.1 km/s (Cancer group B). The groups have velocity dispersions of σ_{v_{gsr}}=12.4+-5.0 km/s and σ_{v_{gsr}}=14.9+-6.2 km/s, and are spatially extended (about several kpc) making it very unlikely that they are bound systems, and are more likely to be debris of tidally disrupted dwarf galaxies or globular clusters. Both groups are metal-poor (median metallicities of [Fe/H] = -1.6 dex and [Fe/H] =-2.1 dex), and have a somewhat uncertain (due to small sample size) metallicity dispersion of ~0.4 dex, suggesting dwarf galaxies as progenitors. Two additional RR Lyrae stars with velocities consistent with those of the Cancer groups have been observed ~25 deg east, suggesting possible extension of the groups in that direction.

preprint2012arXiv

X-ray emission from supernovae in dense circumstellar matter environments: A search for collisionless shocks

(Abridged). The optical light curve of some SNe may be powered by the outward diffusion of the energy deposited by the explosion shock in optically thick circumstellar matter (CSM). Recently, it was shown that the radiation-mediated and -dominated shock in an optically thick wind must transform into a collisionless shock and can produce hard X-rays. The X-rays are expected to peak at late times, relative to maximum visible light. Here we report on a search, using Swift and Chandra, for X-ray emission from 28 SNe that belong to classes whose progenitors are suspected to be embedded in dense CSM (IIn/Ibn/SLSN-I). Two SNe in our sample have X-ray properties that are roughly consistent with the expectation for X-rays from a collisionless shock in optically thick CSM. Therefore, we suggest that their optical light curves are powered by shock breakout in CSM. We show that two other events were too X-ray bright during the SN maximum optical light to be explained by the shock breakout model. We conclude that the light curves of some, but not all, type-IIn/Ibn SNe are powered by shock breakout in CSM. For the rest of the SNe in our sample, including all the SLSN-I events, our X-ray limits are not deep enough and were typically obtained at too early times to conclude about their nature. We argue that the optical light curves of SNe, for which the X-ray emission peaks at late times, are likely powered by the diffusion of shock energy from a dense CSM. We comment about the possibility to detect some of these events in radio.

preprint2011arXiv

Periodic Variability of Low-Mass Stars in SDSS Stripe 82

We present a catalog of periodic stellar variability in the "Stripe 82" region of the Sloan Digital Sky Survey (SDSS). After aggregating and recalibrating catalog-level data from the survey, we ran a period-finding algorithm (Supersmoother) on all point-source lightcurves. We used color selection to identify systems that are likely to contain low-mass stars, in particular M dwarfs and white dwarfs. In total, we found 207 candidates, the vast majority of which appear to be in eclipsing binary systems. The catalog described in this paper includes 42 candidate M dwarf / white dwarf pairs, 4 white-dwarf pairs, 59 systems whose colors indicate they are composed of 2 M dwarfs and whose lightcurve shapes suggest they are in detached eclipsing binaries, and 28 M dwarf systems whose lightcurve shapes suggest they are in contact binaries. We find no detached systems with periods longer than 3 days, thus the majority of our sources are likely to have experienced orbital spin-up and enhanced magnetic activity. Indeed, twenty-six of twenty-seven M dwarf systems that we have spectra for show signs of chromospheric magnetic activity, far higher than the 24% seen in field stars of the same spectral type. We also find binaries composed of stars that bracket the expected boundary between partially and fully convective interiors, which will allow the measurement of the stellar mass-radius relationship across this transition. The majority of our contact systems have short orbital periods, with small variance (0.02 days) in the sample near the observed cutoff of 0.22 days. The accumulation of these stars at short orbital period suggests that the process of angular momentum loss, leading to period evolution, becomes less efficient at short periods. (Abridged)

preprint2011arXiv

The Palomar Transient Factory photometric calibration

The Palomar Transient Factory (PTF) provides multiple epoch imaging for a large fraction of the sky. Here we describe the photometric calibration of the PTF data products that relates the PTF magnitudes to other mag systems. The calibration process utilizes SDSS r~16 mag point source objects as photometric standards. During photometric conditions, this allows us to solve for the extinction coefficients and color terms, and to estimate the camera illumination correction. This also enables the calibration of fields that are outside the SDSS footprint. We test the precision and repeatability of the PTF photometric calibration. Given that PTF is observing in a single filter each night, we define a PTF calibrated magnitude system for R-band and g-band. We show that, in this system, 59% (47%) of the photometrically calibrated PTF R-band (g-band) data achieve a photometric precision of 0.02-0.04 mag. Given the objects' color, the PTF magnitude system can be converted to other systems. A night-by-night comparison of the calibrated magnitudes of individual stars observed on multiple nights shows that they are consistent to a level of ~0.02 mag. Most of the data that were taken under non-photometric conditions can be calibrated relative to other epochs of the same sky footprint obtained during photometric conditions. We provide a guide describing the use of the PTF photometric calibration data products, as well as the transformations between the PTF magnitude system and the SDSS and Johnson-Cousins systems. (abridged).

preprint2010arXiv

Modeling the Time Variability of SDSS Stripe 82 Quasars as a Damped Random Walk

We model the time variability of ~9,000 spectroscopically confirmed quasars in SDSS Stripe 82 as a damped random walk. Using 2.7 million photometric measurements collected over 10 years, we confirm the results of Kelly et al. (2009) and Kozłowski et al. (2010) that this model can explain quasar light curves at an impressive fidelity level (0.01-0.02 mag). The damped random walk model provides a simple, fast [O(N) for N data points], and powerful statistical description of quasar light curves by a characteristic time scale (tau) and an asymptotic rms variability on long time scales (SF_inf). We searched for correlations between these two variability parameters and physical parameters such as luminosity and black hole mass, and rest-frame wavelength. We find that tau increases with increasing wavelength with a power law index of 0.17, remains nearly constant with redshift and luminosity, and increases with increasing black hole mass with power law index of 0.21+/-0.07. The amplitude of variability is anti-correlated with the Eddington ratio, which suggests a scenario where optical fluctuations are tied to variations in the accretion rate. The radio-loudest quasars have systematically larger variability amplitudes by about 30%, when corrected for the other observed trends, while the distribution of their characteristic time scale is indistinguishable from that of the full sample. We do not detect any statistically robust differences in the characteristic time scale and variability amplitude between the full sample and the small subsample of quasars detected by ROSAT. Our results provide a simple quantitative framework for generating mock quasar light curves, such as currently used in LSST image simulations. (abridged)

preprint2010arXiv

Quasar Selection Based on Photometric Variability

We develop a method for separating quasars from other variable point sources using SDSS Stripe 82 light curve data for ~10,000 variable objects. To statistically describe quasar variability, we use a damped random walk model parametrized by a damping time scale, tau, and an asymptotic amplitude (structure function), SF_inf. With the aid of an SDSS spectroscopically confirmed quasar sample, we demonstrate that variability selection in typical extragalactic fields with low stellar density can deliver complete samples with reasonable purity (or efficiency, E). Compared to a selection method based solely on the slope of the structure function, the inclusion of the tau information boosts E from 60% to 75% while maintaining a highly complete sample (98%) even in the absence of color information. For a completeness of C=90%, E is boosted from 80% to 85%. Conversely, C improves from 90% to 97% while maintaining E=80% when imposing a lower limit on tau. With the aid of color selection, the purity can be further boosted to 96%, with C= 93%. Hence, selection methods based on variability will play an important role in the selection of quasars with data provided by upcoming large sky surveys, such as Pan-STARRS and the Large Synoptic Survey Telescope (LSST). For a typical (simulated) LSST cadence over 10 years and a photometric accuracy of 0.03 mag (achieved at i~22), C is expected to be 88% for a simple sample selection criterion of tau>100 days. In summary, given an adequate survey cadence, photometric variability provides an even better method than color selection for separating quasars from stars.

preprint2010arXiv

The Milky Way Tomography with SDSS: III. Stellar Kinematics

We study Milky Way kinematics using a sample of 18.8 million main-sequence stars with r<20 and proper-motion measurements derived from SDSS and POSS astrometry, including ~170,000 stars with radial-velocity measurements from the SDSS spectroscopic survey. Distances to stars are determined using a photometric parallax relation, covering a distance range from ~100 pc to 10 kpc over a quarter of the sky at high Galactic latitudes (|b|>20 degrees). We find that in the region defined by 1 kpc <Z< 5 kpc and 3 kpc <R< 13 kpc, the rotational velocity for disk stars smoothly decreases, and all three components of the velocity dispersion increase, with distance from the Galactic plane. In contrast, the velocity ellipsoid for halo stars is aligned with a spherical coordinate system and appears to be spatially invariant within the probed volume. The velocity distribution of nearby ($Z<1$ kpc) K/M stars is complex, and cannot be described by a standard Schwarzschild ellipsoid. For stars in a distance-limited subsample of stars (<100 pc), we detect a multimodal velocity distribution consistent with that seen by HIPPARCOS. This strong non-Gaussianity significantly affects the measurements of the velocity ellipsoid tilt and vertex deviation when using the Schwarzschild approximation. We develop and test a simple descriptive model for the overall kinematic behavior that captures these features over most of the probed volume, and can be used to search for substructure in kinematic and metallicity space. We use this model to predict further improvements in kinematic mapping of the Galaxy expected from Gaia and LSST.