• (Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.
• ### A Study of the Point Spread Function in SDSS Images(1805.02845)

May 8, 2018 astro-ph.IM
We use SDSS imaging data in $ugriz$ passbands to study the shape of the point spread function (PSF) profile and the variation of its width with wavelength and time. We find that the PSF profile is well described by theoretical predictions based on von K\'{a}rm\'{a}n's turbulence theory. The observed PSF radial profile can be parametrized by only two parameters, the profile's full width at half maximum (FWHM) and a normalization of the contribution of an empirically determined "instrumental" PSF. The profile shape is very similar to the "double gaussian plus power-law wing" decomposition used by SDSS image processing pipeline, but here it is successfully modeled with two free model parameters, rather than six as in SDSS pipeline. The FWHM variation with wavelength follows the $\lambda^{\alpha}$ power law, where $\alpha \approx-0.3$ and is correlated with the FWHM itself. The observed behavior is much better described by von K\'{a}rm\'{a}n's turbulence theory with the outer scale parameter in the range 5$-$100 m, than by the Kolmogorov's turbulence theory. We also measure the temporal and angular structure functions for FWHM and compare them to simulations and results from literature. The angular structure function saturates at scales beyond 0.5$-$1.0 degree. The power spectrum of the temporal behavior is found to be broadly consistent with a damped random walk model with characteristic timescale in the range $\sim5-30$ minutes, though data show a shallower high-frequency behavior. The latter is well fit by a single power law with index in the range $-1.5$ to $-1.0$. A hybrid model is likely needed to fully capture both the low-frequency and high-frequency behavior of the temporal variations of atmospheric seeing.
• ### APO Time Resolved Color Photometry of Highly-Elongated Interstellar Object 1I/'Oumuamua(1711.04927)

Jan. 29, 2018 astro-ph.EP
We report on $g$, $r$ and $i$ band observations of the Interstellar Object 'Oumuamua (1I) taken on 2017 October 29 from 04:28 to 08:40 UTC by the Apache Point Observatory (APO) 3.5m telescope's ARCTIC camera. We find that 1I's colors are $g-r=0.41\pm0.24$ and $r-i=0.23\pm0.25$, consistent with the visible spectra of Masiero (2017), Ye et al. (2017) and Fitzsimmons et al. (2017), and most comparable to the population of Solar System C/D asteroids, Trojans, or comets. We find no evidence of any cometary activity at a heliocentric distance of 1.46 au, approximately 1.5 months after 1I's closest approach distance to the Sun. Significant brightness variability was seen in the $r$ observations, with the object becoming notably brighter towards the end of the run. By combining our APO photometric time series data with the Discovery Channel Telescope (DCT) data of Knight et al. (2017), taken 20 h later on 2017 October 30, we construct an almost complete light curve with a most probable lightcurve period of $P \simeq 4~{\rm h}$. Our results imply a double peaked rotation period of 8.1 $\pm$ 0.02 h, with a peak-to-peak amplitude of 1.5 - 2.1 mags. Assuming that 1I's shape can be approximated by an ellipsoid, the amplitude constraint implies that 1I has an axial ratio of 3.5 to 10.3, which is strikingly elongated. Assuming that 1I is rotating above its critical break up limit, our results are compatible with 1I having having modest cohesive strength and may have obtained its elongated shape during a tidal disruption event before being ejected from its home system. Astrometry useful for constraining 1I's orbit was also obtained and published in Weaver et al. (2017).
• ### The Large Synoptic Survey Telescope as a Near-Earth Object Discovery Machine(1711.10621)

Nov. 29, 2017 astro-ph.EP, astro-ph.IM
Using the most recent prototypes, design, and as-built system information, we test and quantify the capability of the Large Synoptic Survey Telescope (LSST) to discover Potentially Hazardous Asteroids (PHAs) and Near-Earth Objects (NEOs). We empirically estimate an expected upper limit to the false detection rate in LSST image differencing, using measurements on DECam data and prototype LSST software and find it to be about $450$~deg$^{-2}$. We show that this rate is already tractable with current prototype of the LSST Moving Object Processing System (MOPS) by processing a 30-day simulation consistent with measured false detection rates. We proceed to evaluate the performance of the LSST baseline survey strategy for PHAs and NEOs using a high-fidelity simulated survey pointing history. We find that LSST alone, using its baseline survey strategy, will detect $66\%$ of the PHA and $61\%$ of the NEO population objects brighter than $H=22$, with the uncertainty in the estimate of $\pm5$ percentage points. By generating and examining variations on the baseline survey strategy, we show it is possible to further improve the discovery yields. In particular, we find that extending the LSST survey by two additional years and doubling the MOPS search window increases the completeness for PHAs to $86\%$ (including those discovered by contemporaneous surveys) without jeopardizing other LSST science goals ($77\%$ for NEOs). This equates to reducing the undiscovered population of PHAs by additional $26\%$ ($15\%$ for NEOs), relative to the baseline survey.
• ### A hybrid type Ia supernova with an early flash triggered by helium-shell detonation(1710.01824)

Oct. 4, 2017 astro-ph.HE
Type Ia supernovae (SNe Ia) arise from the thermonuclear explosion of carbon-oxygen white dwarfs. Though the uniformity of their light curves makes them powerful cosmological distance indicators, long-standing issues remain regarding their progenitors and explosion mechanisms. Recent detection of the early ultraviolet pulse of a peculiar subluminous SN Ia has been claimed as new evidence for the companion-ejecta interaction through the single-degenerate channel. Here, we report the discovery of a prominent but red optical flash at $\sim$ 0.5 days after the explosion of a SN Ia which shows hybrid features of different SN Ia sub-classes: a light curve typical of normal-brightness SNe Ia, but with strong titanium absorptions, commonly seen in the spectra of subluminous ones. We argue that the early flash of such a hybrid SN Ia is different from predictions of previously suggested scenarios such as the companion-ejecta interaction. Instead it can be naturally explained by a SN explosion triggered by a detonation of a thin helium shell either on a near-Chandrasekhar-mass white dwarf ($\gtrsim$ 1.3 M$_{\odot}$) with low-yield $^{56}$Ni or on a sub-Chandrasekhar-mass white dwarf ($\sim$ 1.0 M$_{\odot}$) merging with a less massive white dwarf. This finding provides compelling evidence that one branch of the previously proposed explosion models, the helium-ignition scenario, does exist in nature, and such a scenario may account for explosions of white dwarfs in a wider mass range in contrast to what was previously supposed.
• The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.
• ### Linear feature detection algorithm for astronomical surveys - II. Defocusing effects on meteor tracks(1707.07223)

July 22, 2017 astro-ph.EP
Given the current limited knowledge of meteor plasma micro-physics and its interaction with the surrounding atmosphere and ionosphere, meteors are a highly interesting observational target for high-resolution wide-field astronomical surveys. Such surveys are capable of resolving the physical size of meteor plasma heads, but they produce large volumes of images that need to be automatically inspected for possible existence of long linear features produced by meteors. Here we show how big aperture sky survey telescopes detect meteors as defocused tracks with a central brightness depression. We derive an analytic expression for a defocused point source meteor track and use it to calculate brightness profiles of meteors modeled as uniform brightness disks. We apply our modeling to meteor images as seen by the SDSS and LSST telescopes. The expression is validated by Monte Carlo ray-tracing simulations of photons traveling through the atmosphere and the LSST telescope optics. We show that estimates of the meteor distance and size can be extracted from the measured FWHM and the strength of the central dip in the observed brightness profile. However, this extraction becomes difficult when the defocused meteor track is distorted by the atmospheric seeing or contaminated by a long lasting glowing meteor trail. The FWHM of satellites tracks is distinctly narrower than meteor values, which enables removal of a possible confusion between satellites and meteors.
• ### Photometric Redshifts with the LSST: Evaluating Survey Observing Strategies(1706.09507)

June 28, 2017 astro-ph.CO, astro-ph.GA
In this paper we present and characterize a nearest-neighbors color-matching photometric redshift estimator that features a direct relationship between the precision and accuracy of the input magnitudes and the output photometric redshifts. This aspect makes our estimator an ideal tool for evaluating the impact of changes to LSST survey parameters that affect the measurement errors of the photometry, which is the main motivation of our work (i.e., it is not intended to provide the "best" photometric redshifts for LSST data). We show how the photometric redshifts will improve with time over the 10-year LSST survey and confirm that the nominal distribution of visits per filter provides the most accurate photo-$z$ results. We also demonstrate how deep LSST imaging of a spectroscopic galaxy sample can significantly improve photo-$z$ quality, especially in the survey's early years. The LSST survey strategy naturally produces observations over a range of airmass, which offers the opportunity of using an SED- and $z$-dependent atmospheric affect on the observed photometry as a color-independent redshift indicator. We show that measuring this airmass effect and including it as a prior has the potential to improve the photometric redshifts and can ameliorate extreme outliers, but also find that it will only be adequately measured for the brightest galaxies, which limits its overall impact on LSST photometric redshifts. Ultimately, we intend for this work to serve as a guide for the expectations and preparations of the LSST science community with regards to the minimum quality of photo-$z$ as the survey progresses.
• ### Machine-Learned Identification of RR Lyrae Stars from Sparse, Multi-band Data: the PS1 Sample(1611.08596)

March 7, 2017 astro-ph.GA
RR Lyrae stars may be the best practical tracers of Galactic halo (sub-)structure and kinematics. The PanSTARRS1 (PS1) $3\pi$ survey offers multi-band, multi-epoch, precise photometry across much of the sky, but a robust identification of RR Lyrae stars in this data set poses a challenge, given PS1's sparse, asynchronous multi-band light curves ($\lesssim 12$ epochs in each of five bands, taken over a 4.5-year period). We present a novel template fitting technique that uses well-defined and physically motivated multi-band light curves of RR Lyrae stars, and demonstrate that we get accurate period estimates, precise to 2~sec in $>80\%$ of cases. We augment these light curve fits with other {\em features} from photometric time-series and provide them to progressively more detailed machine-learned classification models. From these models we are able to select the widest ($3/4$ of the sky) and deepest (reaching 120 kpc) sample of RR Lyrae stars to date. The PS1 sample of $\sim 45,000$ RRab stars is pure (90\%), and complete (80\% at 80 kpc) at high galactic latitudes. It also provides distances precise to 3\%, measured with newly derived period-luminosity relations for optical/near-infrared PS1 bands. With the addition of proper motions from {\em Gaia} and radial velocity measurements from multi-object spectroscopic surveys, we expect the PS1 sample of RR Lyrae stars to become the premier source for studying the structure, kinematics, and the gravitational potential of the Galactic halo. The techniques presented in this study should translate well to other sparse, multi-band data sets, such as those produced by the Dark Energy Survey and the upcoming Large Synoptic Survey Telescope Galactic plane sub-survey.
• ### Everything we'd like to do with LSST data, but we don't know (yet) how(1612.04772)

Dec. 14, 2016 astro-ph.IM
The Large Synoptic Survey Telescope (LSST), the next-generation optical imaging survey sited at Cerro Pachon in Chile, will provide an unprecedented database of astronomical measurements. The LSST design, with an 8.4m (6.7m effective) primary mirror, a 9.6 sq. deg. field of view, and a 3.2 Gigapixel camera, will allow about 10,000 sq. deg. of sky to be covered twice per night, every three to four nights on average, with typical 5-sigma depth for point sources of $r$=24.5 (AB). With over 800 observations in $ugrizy$ bands over a 10-year period, these data will enable a deep stack reaching $r$=27.5 (about 5 magnitudes deeper than SDSS) and faint time-domain astronomy. The measured properties of newly discovered and known astrometric and photometric transients will be publicly reported within 60 sec after observation. The vast database of about 30 trillion observations of 40 billion objects will be mined for the unexpected and used for precision experiments in astrophysics. In addition to a brief introduction to LSST, we discuss a number of astro-statistical challenges that need to be overcome to extract maximum information and science results from LSST dataset.

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

Dec. 29, 2015 astro-ph.GA, astro-ph.HE
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.
• The Large Synoptic Survey Telescope (LSST) is a large-aperture, wide-field, ground-based survey system that will image the sky in six optical bands from 320 to 1050 nm, uniformly covering approximately $18,000$deg$^2$ of the sky over 800 times. The LSST is currently under construction on Cerro Pach\'on in Chile, and expected to enter operations in 2022. Once operational, the LSST will explore a wide range of astrophysical questions, from discovering "killer" asteroids to examining the nature of Dark Energy. The LSST will generate on average 15 TB of data per night, and will require a comprehensive Data Management system to reduce the raw data to scientifically useful catalogs and images with minimum human intervention. These reductions will result in a real-time alert stream, and eleven data releases over the 10-year duration of LSST operations. To enable this processing, the LSST project is developing a new, general-purpose, high-performance, scalable, well documented, open source data processing software stack for O/IR surveys. Prototypes of this stack are already capable of processing data from existing cameras (e.g., SDSS, DECam, MegaCam), and form the basis of the Hyper-Suprime Cam (HSC) Survey data reduction pipeline.
• ### Finding, characterizing and classifying variable sources in multi-epoch sky surveys: QSOs and RR Lyrae in PS1 3$\pi$ data(1511.05527)

Nov. 17, 2015 astro-ph.GA
In area and depth, the Pan-STARRS1 (PS1) 3$\pi$ 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 $(\omega_r, \tau)$ 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.
• ### Asteroid Discovery and Characterization with the Large Synoptic Survey Telescope (LSST)(1511.03199)

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

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

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

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

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

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

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

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

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

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

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