Source author record

Branimir Sesar

Branimir 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

50works
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

50 published item(s)

preprint2020arXiv

The extended Gaia-PS1-SDSS (GPS1+) proper motion catalog

The GPS1 catalog was released in 2017. It delivered precise proper motions for around 350 million sources across three-fourths of the sky down to a magnitude of $r\sim20$\,mag. In this study, we present GPS1+ the extension GPS1 catalog down to $r\sim22.5$\,mag, based on {\it Gaia} DR2, PS1, SDSS and 2MASS astrometry. The GPS1+ totally provides proper motions for $\sim$400 million sources with a characteristic systematic error of less than 0.1\masyr. This catalog is divided into two sub-samples, i.e., the primary and secondary parts. The primary $\sim$264 million sources have either or both of the {\it Gaia} and SDSS astrometry, with a typical precision of 2.0-5.0 \masyr. In this part, $\sim$160 million sources have {\it Gaia} proper motions, we provide another new proper motion for each of them by building a Bayesian model. Relative to {\it Gaia}'s values, the precision is improved by $\sim$0.1\,dex on average at the faint end; $\sim$50 million sources are the objects whose proper motions are missing in {\it Gaia} DR2, we provide their proper motion with a precision of $\sim$4.5\masyr; the remaining $\sim$54 million faint sources are beyond {\it Gaia} detecting capability, we provide their proper motions for the first time with a precision of 7.0 \masyr. However, the secondary $\sim$136 million sources only have PS1 astrometry, the average precision is worse than 15.0 \masyr. All the proper motions have been validated using QSOs and the existing {\it Gaia} proper motions. The catalog will be released on-line and available via the VO-TAP Service, or via the National Astronomical Data Center serviced by China-VO: https://nadc.china-vo.org/data/data/gps1p/f.

preprint2016arXiv

A Synoptic Map of Halo Substructures from the Pan-STARRS1 3π Survey

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

preprint2016arXiv

Galactoseismology: Discovery of a cluster of receding, variable halo stars

A dynamical characterization of dark matter dominated dwarf galaxies from their observed effects on galactic disks (i.e. Galactoseismology) has remained an elusive goal. Here, we present preliminary results from spectroscopic observations of three clustered Cepheid candidates identified from $K$-band light curves towards Norma. The average heliocentric radial velocity of these stars is $\sim$ 156 km/s, which is large and distinct from that of the Galaxy's stellar disk. These objects at $l \sim 333 ^\circ$ and $b \sim -1 ^\circ$ are therefore halo stars; using the $3.6~\micron$ period-luminosity relation of Type I Cepheids, they are at $\sim$ 73 kpc. Our ongoing $I$-band photometry indicates variability on the same time scale as the period determined from the $K_{s}$-band light curve. Distances determined from the $K$-band period-luminosity relation and the 3.6 $\micron$ period-luminosity relation are comparable. The observed radial velocity of these stars agrees roughly with predictions from dynamical models. If these stars are indeed members of the predicted dwarf galaxy that perturbed the outer HI disk of the Milky Way, this would mark the first application of Galactoseismology.

preprint2016arXiv

Maximizing Science in the Era of LSST: A Community-Based Study of Needed US Capabilities

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

preprint2016arXiv

TGAS Error Renormalization from the RR Lyrae Period-Luminosity Relation

The Gaia team has applied a renormalization to their internally-derived parallax errors $σ_{\rm int}(π)$ $$ σ_{tgas}(π) = \sqrt{[Aσ_{int}(π)]^2 + σ_0^2}; \ \ \ \ (A,σ_0) = (1.4,0.20\ \rm mas) $$ based on comparison to Hipparcos astrometry. We use a completely independent method based on the RR Lyrae $K$-band period-luminosity relation to derive a substantially different result, with smaller ultimate errors $$ (A,σ_0) = (1.1,0.12\ \rm mas) \ \ \ \ (this\ paper). $$ We argue that our estimate is likely to be more accurate and therefore that the reported TGAS parallax errors should be reduced according to the prescription: $$ σ_{true}(π) = \sqrt{(0.79σ_{tgas}(π))^2 - (0.10\ \rm mas)^2}. $$

preprint2016arXiv

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

preprint2016arXiv

Why are rapidly rotating M dwarfs in the Pleiades so (infra)red? New period measurements confirm rotation-dependent color offsets from the cluster sequence

Stellar rotation periods measured in open clusters have proved to be extremely useful for studying stars' angular momentum content and rotationally driven magnetic activity, which are both age- and mass-dependent processes. While period measurements have been obtained for hundreds of solar-mass members of the Pleiades, period measurements exist for only a few low-mass ($<$0.5 M$_{\odot}$) members of this key laboratory for stellar evolution theory. To fill this gap, we report rotation periods for 132 low-mass Pleiades members (including nearly 100 with M $\leq$ 0.45 M$_{\odot}$), measured from photometric monitoring of the cluster conducted by the Palomar Transient Factory in late 2011 and early 2012. These periods extend the portrait of stellar rotation at 125 Myr to the lowest-mass stars and re-establish the Pleiades as a key benchmark for models of the transport and evolution of stellar angular momentum. Combining our new rotation periods with precise BVIJHK photometry reported by Stauffer et al. and Kamai et al., we investigate known anomalies in the photometric properties of K and M Pleiades members. We confirm the correlation detected by Kamai et al. between a star's rotation period and position relative to the main sequence in the cluster's color-magnitude diagram. We find that rapid rotators have redder (V-K) colors than slower rotators at the same V, indicating that rapid and slow rotators have different binary frequencies and/or photospheric properties. We find no difference in the photometric amplitudes of rapid and slow rotators, indicating that asymmetries in the longitudinal distribution of starspots do not scale grossly with rotation rate.

preprint2015arXiv

A New Faint Milky Way Satellite Discovered in the Pan-STARRS1 3 pi Survey

We present the discovery of a faint Milky Way satellite, Laevens 2/Triangulum II, found in the Panoramic Survey Telescope And Rapid Response System (Pan-STARRS 1) 3 pi imaging data and confirmed with follow-up wide-field photometry from the Large Binocular Cameras. The stellar system, with an absolute magnitude of M_V=-1.8 +/-0.5, a heliocentric distance of 30 +2/-2 kpc, and a half-mass radius of 34 +9/-8 pc, shows remarkable similarity to faint, nearby, small satellites such as Willman 1, Segue 1, Segue 2, and Boötes II. The discovery of Laevens 2/Triangulum II further populates the region of parameter space for which the boundary between dwarf galaxies and globular clusters becomes tenuous. Follow-up spectroscopy will ultimately determine the nature of this new satellite, whose spatial location hints at a possible connection with the complex Triangulum-Andromeda stellar structures.

preprint2015arXiv

A re-interpretation of the Triangulum-Andromeda stellar clouds: a population of halo stars kicked out of the Galactic disk

The Triangulum-Andromeda stellar clouds (TriAnd1 and TriAnd2) are a pair of concentric ring- or shell-like over-densities at large $R$ ($\approx$ 30 kpc) and $Z$ ($\approx$ -10 kpc) in the Galactic halo that are thought to have been formed from the accretion and disruption of a satellite galaxy. This paper critically re-examines this formation scenario by comparing the number ratio of RR Lyrae to M giant stars associated with the TriAnd clouds with other structures in the Galaxy. The current data suggest a stellar population for these over-densities ($f_{\rm RR:MG} < 0.38$ at 95% confidence) quite unlike any of the known satellites of the Milky Way ($f_{\rm RR:MG} \approx 0.5$ for the very largest and $f_{\rm RR:MG} >>1$ for the smaller satellites) and more like the population of stars born in the much deeper potential well inhabited by the Galactic disk ($f_{\rm RR:MG} < 0.01$). N-body simulations of a Milky-Way-like galaxy perturbed by the impact of a dwarf galaxy demonstrate that, in the right circumstances, concentric rings propagating outwards from that Galactic disk can plausibly produce similar over-densities. These results provide dramatic support for the recent proposal by Xu et al. (2015) that, rather than stars accreted from other galaxies, the TriAnd clouds could represent stars kicked-out from our own disk. If so, these would be the first populations of disk stars to be found in the Galactic halo and a clear signature of the importance of this second formation mechanism for stellar halos more generally. Moreover, their existence at the very extremities of the disk places strong constraints on the nature of the interaction that formed them.

preprint2015arXiv

Asteroid Spin-Rate Study using the Intermediate Palomar Transient Factory

Two dedicated asteroid rotation-period surveys have been carried out using data taken on January 6-9 and February 20-23 of 2014 by the Intermediate Palomar Transient Factory (iPTF) in the $R$~band with $\sim 20$-min cadence. The total survey area covered 174~deg$^2$ in the ecliptic plane. Reliable rotation periods for 1,438 asteroids are obtained from a larger data set of 6,551 mostly main-belt asteroids, each with $\geq 10$~detections. Analysis of 1751, PTF based, reliable rotation periods clearly shows the "spin barrier" at $\sim 2$~hours for "rubble-pile" asteroids. We also found a new large-sized super-fast rotator, 2005 UW163 (Chang et al., 2014), and other five candidates as well. Our spin-rate distributions of asteroids with $3 < D < 15$~km shows number decrease when frequency greater than 5 rev/day, which is consistent to that of the Asteroid Light Curve Database (LCDB, Warner et al., 2009) and the result of (Masiero et al., 2009). We found the discrepancy in the spin-rate distribution between our result and (Pravec et al., 2008, update 2014-04-20) is mainly from asteroids with $Δm < 0.2$ mag that might be primarily due to different survey strategies. For asteroids with $D \leq 3$~km, we found a significant number drop at $f = 6$ rev/day. The YORP effect timescale for small-sized asteroid is shorter that makes more elongate objets spun up to reach their spin-rate limit and results in break-up. The K-S test suggests a possible difference in the spin-rate distributions of C- and S-type asteroids. We also find that C-type asteroids have a smaller spin-rate limit than the S-type, which agrees with the general sense that the C-type has lower bulk density than the S-type.

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

Properties and Evolution of the Redback Millisecond Pulsar Binary PSR J2129-0429

PSR J2129-0429 is a "redback" eclipsing millisecond pulsar binary with an unusually long 15.2 hour orbit. It was discovered by the Green Bank Telescope in a targeted search of unidentified Fermi gamma-ray sources. The pulsar companion is optically bright (mean $m_R = 16.6$ mag), allowing us to construct the longest baseline photometric dataset available for such a system. We present ten years of archival and new photometry of the companion from LINEAR, CRTS, PTF, the Palomar 60-inch, and LCOGT. Radial velocity spectroscopy using the Double-Beam Spectrograph on the Palomar 200-inch indicates that the pulsar is massive: $1.74\pm0.18 M_\odot$. The G-type pulsar companion has mass $0.44\pm0.04 M_\odot$, one of the heaviest known redback companions. It is currently 95\% Roche-lobe filling and only mildly irradiated by the pulsar. We identify a clear 13.1 mmag yr$^{-1}$ secular decline in the mean magnitude of the companion as well as smaller-scale variations in the optical lightcurve shape. This behavior may indicate that the companion is cooling. Binary evolution calculations indicate that PSR J2129-0429 has an orbital period almost exactly at the bifurcation period between systems that converge into tighter orbits as black widows and redbacks and those that diverge into wider pulsar--white dwarf binaries. Its eventual fate may depend on whether it undergoes future episodes of mass transfer and increased irradiation.

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

Sagittarius II, Draco II and Laevens 3: three new Milky Way satellites discovered in the Pan-STARRS 1 3pi Survey

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

preprint2015arXiv

The Outer Galactic Halo As Probed By RR Lyr Stars From the Palomar Transient Facility + Keck

We present initial results from our study of the outer halo of the Milky Way using a large sample of RR Lyr(ab) variables datamined from the archives of the Palomar Transient Facility. Of the 464 RR Lyr in our sample with distances exceeding 50 kpc, 62 have been observed spectroscopically at the Keck Observatory. Radial velocities and sigma(vr) are given as a function of distance between 50 and 110 kpc, and a very preliminary rather low total mass for the Milky Way out to 110 kpc of ~7 (+-1.5) x 10**11 solar masses is derived from our data.

preprint2015arXiv

The radial profile and flattening of the Milky Way's stellar halo to $\rm 80~$kpc from the SEGUE K-giant Survey

We characterise the radial density, metallicity and flattening profile of the Milky Way's stellar halo, based on the large sample of 1757 spectroscopically confirmed giant stars from SDSS/SEGUE-2 after excising stars that were algorithmically attributed to apparent halo substructure (including the Sagittarius stream). Compared to BHB stars or RR Lyrae, giants are more readily understood tracers of the overall halo star population, with less bias in age or metallicity. The well-characterized selection function of the sample enables forward modelling of those data, based on ellipsoidal stellar density models, $ν_* (R,z)$, with Einasto profiles and (broken) power laws for their radial dependence, combined with a model for the metallicity gradient and the flattening profile. Among models with constant flattening, these data are reasonably well fit by an Einasto profile of $n=3.1\pm 0.5$ with an effective radius $\rm r_{eff} = 15\pm2~$kpc and a flattening of $q=0.7\pm 0.02$; or comparably well by an equally flattened broken power-law, with radial slopes of $α_{in}=2.1\pm 0.3$ and $α_{out}=3.8\pm 0.1$, with a break-radius of $r_{break}=18\pm1$~kpc; this is largely consistent with earlier work. We find a modest, but significant metallicity gradient within the "outer" stellar halo, $\rm [Fe/H]$ decreasing outward. If we allow for a variable flattening $q = f(r_{GC} )$, we find the distribution of halo giants to be considerably more flattened at small radii, $q({\rm 10~kpc})\sim 0.57$, compared to $q(>30{\rm kpc})\sim 0.8$. Remarkably, the data are then very well fit by a single power-law of index $\rm \sim 4.2\pm0.1$ of the variable $r_q\equiv\sqrt{R^2+(z/q(r))^2}$. In this simple and better fitting model, there is a break in flattening at $\sim 20$~kpc, instead of a break in the radial density function.

preprint2014arXiv

313 new asteroid rotation periods from Palomar Transient Factory observations

A new asteroid rotation period survey have been carried out by using the Palomar Transient Factory (PTF). Twelve consecutive PTF fields, which covered an area of 87 deg$^2$ in the ecliptic plane, were observed in $R$ band with a cadence of $\sim$20 min during February 15--18, 2013. We detected 2500 known asteroids with a diameter range of 0.5 km $\leq D \leq$ 200 km. Of these, 313 objects had highly reliable rotation periods and exhibited the "spin barrier" at $\sim2$ hours. In contrast to the flat spin rate distribution of the asteroids with 3 km $\leq D \leq$ 15 km shown by Pravec et al. (2008), our results deviated somewhat from a Maxwellian distribution and showed a decrease at the spin rate greater than 5 rev/day. One super-fast-rotator candidate and two possible binary asteroids were also found in this work.

preprint2014arXiv

A new distant Milky Way globular cluster in the Pan-STARRS1 3π survey

We present a new satellite in the outer halo of the Galaxy, the first Milky Way satellite found in the stacked photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS1) Survey. From follow-up photometry obtained with WFI on the MPG/ESO 2.2m telescope, we argue that the object, located at a heliocentric distance of 145+/-17 kpc, is the most distant Milky Way globular cluster yet known. With a total magnitude of M_V=-4.3+/-0.2 and a half-light radius of 20+/-2 pc, it shares the properties of extended globular clusters found in the outer halo of our Galaxy and the Andromeda galaxy. The discovery of this distant cluster shows that the full spatial extent of the Milky Way globular cluster system has not yet been fully explored.

preprint2014arXiv

Identification of the Optical Counterpart of Fermi Black Widow Millisecond Pulsar PSR J1544+4937

We report the optical identification of the companion to the {\it Fermi} black widow millisecond pulsar PSR J1544+4937. We find a highly variable source on Keck LRIS images at the nominal pulsar position, with 2 magnitude variations over orbital period in the B, g, R, and I bands. The nearly achromatic light curves are difficult to explain with a simply irradiated hemisphere model, and suggest that the optical emission is dominated by a nearly isothermal hot patch on the surface of the companion facing the pulsar. We roughly constrain the distance to PSR J1544+4937 to be between 2 and 5 kpc. A more reliable distance measurement is needed in order to constrain the composition of the companion.

preprint2014arXiv

IPAC Image Processing and Data Archiving for the Palomar Transient Factory

The Palomar Transient Factory (PTF) is a multi-epochal robotic survey of the northern sky that acquires data for the scientific study of transient and variable astrophysical phenomena. The camera and telescope provide for wide-field imaging in optical bands. In the five years of operation since first light on December 13, 2008, images taken with Mould-R and SDSS-g' camera filters have been routinely acquired on a nightly basis (weather permitting), and two different H-alpha filters were installed in May 2011 (656 nm and 663 nm). The PTF image-processing and data-archival program at the Infrared Processing and Analysis Center (IPAC) is tailored to receive and reduce the data, and, from it, generate and preserve astrometrically and photometrically calibrated images, extracted source catalogs, and coadded reference images. Relational databases have been deployed to track these products in operations and the data archive. The fully automated system has benefited by lessons learned from past IPAC projects and comprises advantageous features that are potentially incorporable into other ground-based observatories. Both off-the-shelf and in-house software have been utilized for economy and rapid development. The PTF data archive is curated by the NASA/IPAC Infrared Science Archive (IRSA). A state-of-the-art custom web interface has been deployed for downloading the raw images, processed images, and source catalogs from IRSA. Access to PTF data products is currently limited to an initial public data release (M81, M44, M42, SDSS Stripe 82, and the Kepler Survey Field). It is the intent of the PTF collaboration to release the full PTF data archive when sufficient funding becomes available.

preprint2014arXiv

Long-term Photometric Behavior of Outbursting AM CVn Systems

The AM CVn systems are a class of He-rich, post-period minimum, semi-detached, ultra-compact binaries. Their long-term light curves have been poorly understood due to the few systems known and the long (hundreds of days) recurrence times between outbursts. We present combined photometric light curves from the LINEAR, CRTS, and PTF synoptic surveys to study the photometric variability of these systems over an almost 10 yr period. These light curves provide a much clearer picture of the outburst phenomena that these systems undergo. We characterize the photometric behavior of most known outbursting AM CVn systems and establish a relation between their outburst properties and the systems' orbital periods. We also explore why some systems have only shown a single outburst so far and expand the previously accepted phenomenological states of AM CVn systems. We conclude that the outbursts of these systems show evolution with respect to the orbital period, which can likely be attributed to the decreasing mass transfer rate with increasing period. Finally, we consider the number of AM CVn systems that should be present in modeled synoptic surveys.

preprint2014arXiv

PTF1 J191905.19+481506.2 - A Partially Eclipsing AM CVn System Discovered in the Palomar Transient Factory

We report on PTF1 J191905.19+481506.2, a newly discovered, partially eclipsing, outbursting AM CVn system found in the Palomar Transient Factory synoptic survey. This is only the second known eclipsing AM CVn system. We use high-speed photometric observations and phase-resolved spectroscopy to establish an orbital period of 22.4559(3) min. We also present a long-term light curve and report on the normal and super-outbursts regularly seen in this system, including a super-outburst recurrence time of 36.8(4) d. We use the presence of the eclipse to place upper and lower limits on the inclination of the system and discuss the number of known eclipsing AM CVn systems versus what would be expected.

preprint2014arXiv

Serendipitous Discovery of a Thin Stellar Stream near the Galactic Bulge in the Pan-STARRS1 3Pi Survey

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

preprint2013arXiv

Detection of a Nearby Halo Debris Stream in the WISE and 2MASS Surveys

Combining the Wide-Field Infrared Survey Explorer All-Sky Release with the 2MASS Point Source Catalog, we detect a nearby, moderately metal-poor stellar debris stream spanning 24 degrees across the southern sky. The stream, which we designate Alpheus, is at an estimated distance of ~1.9 kpc. Its position, orientation, width, estimated metallicity, and to some extent its distance, are in approximate agreement with what one might expect of the leading tidal tail of the southern globular cluster NGC 288.

preprint2013arXiv

Discovery of a Cosmological, Relativistic Outburst via its Rapidly Fading Optical Emission

We report the discovery by the Palomar Transient Factory (PTF) of the transient source PTF11agg, which is distinguished by three primary characteristics: (1) bright, rapidly fading optical transient emission; (2) a faint, blue quiescent optical counterpart; and (3) an associated year-long, scintillating radio transient. We argue that these observed properties are inconsistent with any known class of Galactic transients, and instead suggest a cosmological origin. The detection of incoherent radio emission at such distances implies a large emitting region, from which we infer the presence of relativistic ejecta. The observed properties are all consistent with the population of long-duration gamma-ray bursts (GRBs), marking the first time such an outburst has been discovered in the distant universe independent of a high-energy trigger. We searched for possible high-energy counterparts to PTF11agg, but found no evidence for associated prompt emission. We therefore consider three possible scenarios to account for a GRB-like afterglow without a high-energy counterpart: an "untriggered" GRB (lack of satellite coverage), an "orphan" afterglow (viewing-angle effects), and a "dirty fireball" (suppressed high-energy emission). The observed optical and radio light curves appear inconsistent with even the most basic predictions for off-axis afterglow models. The simplest explanation, then, is that PTF11agg is a normal, on-axis long-duration GRB for which the associated high-energy emission was simply missed. However, we have calculated the likelihood of such a serendipitous discovery by PTF and find that it is quite small (~ 2.6%). While not definitive, we nonetheless speculate that PTF11agg may represent a new, more common (> 4 times the on-axis GRB rate at 90% confidence) class of relativistic outbursts lacking associated high-energy emission.

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

Main-belt comets in the Palomar Transient Factory survey: I. The search for extendedness

Cometary activity in main-belt asteroids probes the ice content of these objects and provides clues to the history of volatiles in the inner solar system. We search the Palomar Transient Factory (PTF) survey to derive upper limits on the population size of active main-belt comets (MBCs). From data collected March 2009 through July 2012, we extracted 2 million observations of 220 thousand known main-belt objects (40% of the known population, down to 1-km diameter) and discovered 626 new objects in multi-night linked detections. We formally quantify the extendedness of a small-body observation, account for systematic variation in this metric (e.g., due to on-sky motion) and evaluate this method's robustness in identifying cometary activity using observations of 115 comets, including two known candidate MBCs and six newly-discovered non-main-belt comets (two of which were originally designated as asteroids by other surveys). We demonstrate a 66% detection efficiency with respect to the extendedness distribution of the 115 sampled comets, and a 100% detection efficiency with respect to extendedness levels greater than or equal to those we observed in the known candidate MBCs P/2010 R2 (La Sagra) and P/2006 VW139. Using a log-constant prior, we infer 95% confidence upper limits of 33 and 22 active MBCs (per million main-belt asteroids down to 1-km diameter), for detection efficiencies of 66% and 100%, respectively. In a follow-up to this morphological search, we will perform a photometric (disk-integrated brightening) search for MBCs.

preprint2013arXiv

Metal Abundances, Radial Velocities and other Physical Characteristics for the RR Lyrae Stars in the Kepler Field

Spectroscopic iron-to-hydrogen ratios, radial velocities, atmospheric parameters, and new photometric analyses are presented for 41 RR Lyrae stars (and one probable high-amplitude delta Scuti star) located in the field-of-view of the Kepler space telescope. Thirty-seven of the RR Lyrae stars are fundamental-mode pulsators (i.e., RRab stars) of which 16 exhibit the Blazhko effect. Four of the stars are multiperiodic RRc pulsators oscillating primarily in the first-overtone mode. Spectroscopic [Fe/H] values for the 34 stars for which we were able to derive estimates range from -2.54 +/- 0.13 (NR Lyr) to -0.05 +/- 0.13 dex (V784 Cyg), and for the 19 Kepler-field non-Blazhko stars studied by Nemec et al.(2011) the abundances agree will with their photometric [Fe/H] values. Four non-Blazhko RR Lyrae stars that they identified as metal-rich (KIC 6100702, V2470 Cyg, V782 Cyg and V784 Cyg) are confirmed as such, and four additional stars (V839 Cyg, KIC 5520878, KIC 8832417, KIC 3868420) are also shown here to be metal-rich. Five of the non-Blazhko RRab stars are found to be more metal-rich than [Fe/H] -0.9 dex while all of the 16 Blazhko stars are more metal-poor than this value. New period-phi31s-[Fe/H] relationships are derived based on ~970 days of quasi-continuous high-precison Q0-Q11 long- and short-cadence Kepler photometry. With the exception of some Blazhko stars, the spectroscopic and photometric [Fe/H] values are in good agreement. Several stars with unique photometric characteristics are identified, including a Blazhko variable with the smallest known amplitude and frequency modulations (V838 Cyg).

preprint2013arXiv

R Coronae Borealis Stars in M31 from the Palomar Transient Factory

We report the discovery of R Coronae Borealis (RCB) stars in the Andromeda galaxy (M31) using the Palomar Transient Factory (PTF). RCB stars are rare hydrogen-deficient, carbon-rich supergiant variables, most likely the merger products of two white dwarfs. These new RCBs, including two confirmed ones and two candidates, are the first to be found beyond the Milky Way and the Magellanic Clouds. All of M31 RCBs showed >1.5 mag irregular declines over timescales of weeks to months. Due to the limiting magnitude of our data (R~21-22 mag), these RCB stars have R~19.5 to 20.5 mag at maximum light, corresponding to M_R= -4 to -5, making them some of the most luminous RCBs known. Spectra of two objects show that they are warm RCBs, similar to the Milky Way RCBs RY Sgr and V854 Cen. We consider these results, derived from a pilot study of M31 variables, as an important proof-of-concept for the study of rare bright variables in nearby galaxies with the PTF or other synoptic surveys.

preprint2013arXiv

Statistical Searches for Microlensing Events in Large, Non-Uniformly Sampled Time-Domain Surveys: A Test Using Palomar Transient Factory Data

Many photometric time-domain surveys are driven by specific goals, such as searches for supernovae or transiting exoplanets, which set the cadence with which fields are re-imaged. In the case of the Palomar Transient Factory (PTF), several sub-surveys are conducted in parallel, leading to non-uniform sampling over its $\sim$$20,000 \mathrm{deg}^2$ footprint. While the median $7.26 \mathrm{deg}^2$ PTF field has been imaged $\sim$40 times in \textit{R}-band, $\sim$$2300 \mathrm{deg}^2$ have been observed $>$100 times. We use PTF data to study the trade-off between searching for microlensing events in a survey whose footprint is much larger than that of typical microlensing searches, but with far-from-optimal time sampling. To examine the probability that microlensing events can be recovered in these data, we test statistics used on uniformly sampled data to identify variables and transients. We find that the von Neumann ratio performs best for identifying simulated microlensing events in our data. We develop a selection method using this statistic and apply it to data from fields with $>$10 $R$-band observations, $1.1\times10^9$ light curves, uncovering three candidate microlensing events. We lack simultaneous, multi-color photometry to confirm these as microlensing events. However, their number is consistent with predictions for the event rate in the PTF footprint over the survey's three years of operations, as estimated from near-field microlensing models. This work can help constrain all-sky event rate predictions and tests microlensing signal recovery in large data sets, which will be useful to future time-domain surveys, such as that planned with the Large Synoptic Survey Telescope.

preprint2013arXiv

The PTF Orion Project: a Possible Planet Transiting a T-Tauri Star

We report observations of a possible young transiting planet orbiting a previously known weak-lined T-Tauri star in the 7-10 Myr old Orion-OB1a/25-Ori region. The candidate was found as part of the Palomar Transient Factory (PTF) Orion project. It has a photometric transit period of 0.448413 +- 0.000040 days, and appears in both 2009 and 2010 PTF data. Follow-up low-precision radial velocity (RV) observations and adaptive optics imaging suggest that the star is not an eclipsing binary, and that it is unlikely that a background source is blended with the target and mimicking the observed transit. RV observations with the Hobby-Eberly and Keck telescopes yield an RV that has the same period as the photometric event, but is offset in phase from the transit center by approximately -0.22 periods. The amplitude (half range) of the RV variations is 2.4 km/s and is comparable with the expected RV amplitude that stellar spots could induce. The RV curve is likely dominated by stellar spot modulation and provides an upper limit to the projected companion mass of M_p sin i_orb < 4.8 +- 1.2 M_Jup; when combined with the orbital inclination, i orb, of the candidate planet from modeling of the transit light curve, we find an upper limit on the mass of the planetary candidate of M_p < 5.5 +- 1.4 M_Jup. This limit implies that the planet is orbiting close to, if not inside, its Roche limiting orbital radius, so that it may be undergoing active mass loss and evaporation.

preprint2013arXiv

Tracing the Orphan Stream to 55 kpc with RR Lyrae Stars

We report positions, velocities and metallicities of 50 ab-type RR Lyrae (RRab) stars observed in the vicinity of the Orphan stellar stream. Using about 30 RRab stars classified as being likely members of the Orphan stream, we study the metallicity and the spatial extent of the stream. We find that RRab stars in the Orphan stream have a wide range of metallicities, from -1.5 dex to -2.7 dex. The average metallicity of the stream is -2.1 dex, identical to the value obtained by Newberg et al. (2010) using blue horizontal branch stars. We find that the most distant parts of the stream (40-50 kpc from the Sun) are about 0.3 dex more metal-poor than the closer parts (within ~30 kpc), suggesting a possible metallicity gradient along the stream's length. We have extended the previous studies and have mapped the stream up to 55 kpc from the Sun. Even after a careful search, we did not identify any more distant RRab stars that could plausibly be members of the Orphan stream. If confirmed with other tracers, this result would indicate a detection of the end of the leading arm of the stream. We have compared the distances of Orphan stream RRab stars with the best-fit orbits obtained by Newberg et al. (2010). We find that model 6 of Newberg et al. (2010) cannot explain the distances of the most remote Orphan stream RRab stars, and conclude that the best fit to distances of Orphan stream RRab stars and to the local circular velocity is provided by potentials where the total mass of the Galaxy within 60 kpc is M_{60}~2.7x10^{11} Msun, or about 60% of the mass found by previous studies. More extensive modelling that would consider non-spherical potentials and the possibility of misalignment between the stream and the orbit, is highly encouraged.

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

Addressing Decadal Survey Science through Community Access to Highly Multiplexed Spectroscopy with BigBOSS on the KPNO Mayall Telescope

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

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

Five New Outbursting AM CVn Systems Discovered by the Palomar Transient Factory

We present five new outbursting AM CVn systems and one candidate discovered as part of an ongoing search for such systems using the Palomar Transient Factory (PTF). This is the first large-area, systematic search for AM CVn systems using only large-amplitude photometric variability to select candidates. Three of the confirmed systems and the candidate system were discovered as part of the PTF transient search. Two systems were found as part of a search for outbursts through the PTF photometric database. We discuss the observed characteristics of each of these systems, including the orbital periods of two systems. We also consider the position of these systems, selected in a colour-independent survey, in colour-colour space and compare to systems selected solely by their colours. We find that the colours of our newly discovered systems do not differ significantly from those of previously known systems, but significant errors preclude a definitive answer.

preprint2012arXiv

Multi-wavelength characterization of stellar flares on low-mass stars using SDSS and 2MASS time domain surveys

We present the first rates of flares from M dwarf stars in both red optical and near infrared (NIR) filters. We have studied ~50,000 M dwarfs from the SDSS Stripe 82 area, and 1,321 M dwarfs from the 2MASS Calibration Scan Point Source Working Database that overlap SDSS imaging fields. We assign photometric spectral types from M0 to M6 using (r-i) and (i-z) colors for every star in our sample. Stripe 82 stars each have 50-100 epochs of data, while 2MASS Calibration stars have ~1900 epochs. From these data we estimate the observed rates and theoretical detection thresholds for flares in eight photometric bands as a function of spectral type. Optical flare rates are found to be in agreement with previous studies, while the frequency per hour of NIR flare detections is found to be more than two orders of magnitude lower. An excess of small amplitude flux increases in all bands exhibits a power-law distribution, which we interpret as the result of flares below our detection thresholds. In order to investigate the recovery efficiency for flares in each filter, we extend a two-component flare model into the NIR. Quiescent M0-M6 spectral templates were used with the model to predict the photometric response of flares from u to Ks. We determine that red optical filters are sensitive to flares with u-band amplitudes >2 mag, and NIR filters to flares with delta u>4.5 mag. Our model predicts that M0 stars have the best color-contrast for J-band detections, but M4-M6 stars should show the highest rate of NIR flares with amplitudes of delta J ~0.01 mag. Characterizing flare rates and photometric variations at longer wavelengths is important for predicting the signatures of M dwarf variability in next-generation surveys, and we discuss their impact on surveys such as LSST.

preprint2012arXiv

Template RR Lyrae H alpha, H beta, and H gamma Velocity Curves

We present template radial velocity curves of $ab$-type RR Lyrae stars constructed from high-precision measurements of ${\rm Hα}$, ${\rm Hβ}$, and ${\rm Hγ}$ lines. Amplitude correlations between the Balmer line velocity curves, Johnson $V$-band, and SDSS $g$- and $r$-band light curves are also derived. Compared to previous methods, these templates and derived correlations reduce the uncertainty in measured systemic (center-of-mass) velocities of RR Lyrae stars by up to 15 {\kms}, and will be of particular interest to wide-area spectroscopic surveys such as the Sloan Digital Sky Survey (SDSS) and LAMOST Experiment for Galactic Understanding and Exploration (LEGUE).

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

Mapping the Galactic Halo with SDSS, LINEAR and PTF RR Lyrae Stars

We present an analysis of Galactic halo structure and substructure traced by ~500 RR Lyrae stars in the SDSS stripe 82 region. The main result of this study is a 2D map of the Galactic halo that reaches distances of 100 kpc and traces previously known and new halo substructures, such as the Sagittarius and Pisces tidal streams. We also present strong direct evidence, based on both RR Lyrae and main sequence stars, that the halo stellar number density profile significantly steepens beyond 30 kpc from the Galactic center. Using a novel photometric metallicity method that simultaneously combines data for RR Lyrae and main sequence stars, we show that the median metallicity of the Sagittarius trailing stream in the SDSS stripe 82 region is [Fe/H]}=-1.2 +- 0.1 dex. In addition to these results, we will present ongoing work on a 3D map of the Galactic halo constructed using a sample of ~4000 RR Lyrae stars that covers 8500 deg^2 of northern sky and probes up to 30 kpc from the Sun.

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.

preprint2011arXiv

Mining Databases for M Dwarf Variability

Time-resolved databases with large spatial coverage are quickly becoming a standard tool for all types of astronomical studies. We report preliminary results from our search for stellar flares in the 2MASS calibration fields. A sample of 4343 M dwarfs, spatially matched between the SDSS and the 2MASS calibration fields, each with hundreds to thousands of epochs in near infrared bandpasses, is analyzed using a modified Welch-Stetson index to characterize the variability. A Monte Carlo model was used to assess the noise of the variability index. We find significnat residuals above the noise with power-law slopes of -3.37 and -4.05 for our JH and HKs distributions respectively. This is evidence for flares being observed from M dwarfs in infrared photometry.

preprint2011arXiv

PTF1 J071912.13+485834.0: An outbursting AM CVn system discovered by a synoptic survey

We present extensive photometric and spectroscopic observations of PTF1 J071912.13+485834.0, an outbursting AM CVn system discovered by the Palomar Transient Factory (PTF). AM CVn systems are stellar binaries with some of the smallest separations known and orbital periods ranging from 5 to 65 minutes. They are believed to be composed of a white dwarf accretor and a (semi)-degenerate He-rich donor and are considered to be the helium equivalents of Cataclysmic Variables. We have spectroscopically and photometrically identified an orbital period of 26.77 \pm 0.02 minutes for PTF1 J071912.13+485834.0 and found a super-outburst recurrence time of greater than 65 days along with the presence of "normal" outbursts - rarely seen in AM CVn systems but well known in super-outbursting Cataclysmic Variables. We present a long-term light curve over two super-cycles as well as high cadence photometry of both outburst and quiescent stages, both of which show clear variability. We also compare both the outburst and quiescent spectra of PTF1 J071912.13+485834.0 to other known AM CVn systems, and use the quiescent phase-resolved spectroscopy to determine the origin of the photometric variability. Finally, we draw parallels between the different subclasses of SU UMa-type Cataclysmic Variables and outbursting AM CVn systems. We conclude by predicting that the Palomar Transient Factory may more than double the number of outbursting AM CVn systems known, which would greatly increase our understanding of AM CVn systems.

preprint2011arXiv

SN2010jp (PTF10aaxi): A Jet-Driven Type II Supernova

We present photometry and spectroscopy of the peculiar TypeII supernova (SN) 2010jp, also named PTF10aaxi. The light curve exhibits a linear decline with a relatively low peak absolute magnitude of only -15.9, and a low radioactive decay luminosity at late times that suggests a nickel mass below 0.003 $M_{\odot}$. Spectra of SN2010jp display an unprecedented triple-peaked H$α$ line profile, showing: (1) a narrow (800 km/s) central component that suggests shock interaction with dense CSM; (2) high-velocity blue and red emission features centered at -12600 and +15400 km/s; and (3) broad wings extending from -22000 to +25000 km/s. These features persist during 100 days after explosion. We propose that this line profile indicates a bipolar jet-driven explosion, with the central component produced by normal SN ejecta and CSM interaction at mid latitudes, while the high-velocity bumps and broad line wings arise in a nonrelativistic bipolar jet. Two variations of the jet interpretation seem plausible: (1) A fast jet mixes 56Ni to high velocities in polar zones of the H-rich envelope, or (2) the reverse shock in the jet produces blue and red bumps in Balmer lines when a jet interacts with dense CSM. Jet-driven SNeII are predicted for collapsars resulting from a wide range of initial masses above 25 $M_{\odot}$ at sub-solar metallicity. This seems consistent with the SN host environment, which is either an extremely low-luminosity dwarf galaxy or very remote parts of an interacting pair of star-forming galaxies. It also seems consistent with the low 56Ni mass that may accompany black hole formation. We speculate that the jet survives to produce observable signatures because the star's H envelope was mostly stripped away by previous eruptive mass loss.

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

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.

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

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.

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.