Search

Changing-Look Active Galactic Nuclei (CLAGN) are characterised by extreme variations in line emission over short timescales, mostly affecting broad H$\beta$ lines. While a few hundred CLAGN are known, a complete sample of turn-on CLAGN is still elusive. Here, we present a search for turn-on CLAGN in a complete sample of galaxies, using archival spectra and recent light curves. We obtained light curves from the Asteroid Terrestrial Impact Last Alert System (ATLAS) for 16,232 emission line galaxies, including both star-forming and active galaxies, at $z<0.1$ with spectra from the Six-degree Field Galaxy Survey (6dFGS). We first establish typical variability behaviour for different AGN types, as recorded between 2001 and 2009, and then select outliers from the bulk behaviour as CLAGN candidates. We obtain new spectra for the candidates and identify 12 new turn-on CLAGN (appearing broad H$\beta$ line) and 19 new turn-off CLAGN (disappearing broad H$\beta$ line). We may have missed AGN that changed and reverted their state over the 15$-$20 years since 6dFGS spectra were taken, and thus our CLAGN rates of 1.7$\%$ for turn-on and 9.6$\%$ for turn-off are lower limits. The turn-on rate is naturally much lower as the type 1.9/2 sample is dominated by obscured AGN due to orientation, which are not expected to change. However, the number of turn-on (27) and turn-off (24) CLAGN we find are similar, suggesting that our parent AGN sample is reasonably complete in our search volume at $z<0.1$., Comment: 19 pages, 12 figures. Accepted by MNRAS

We utilize serendipitous observations from the Asteroid Terrestrial-impact Last Alert System (ATLAS) and the Zwicky Transient Facility (ZTF) in addition to targeted follow-up observations from the Las Cumbres Observatory (LCO) and Liverpool Telescope to analyze the first observed instance of cometary activity by the newly-discovered Jupiter-family comet C/2023 RN3 (ATLAS), whose orbital dynamics place it close to residing on a Centaur-like orbit. Across our 7-month baseline, we observe an epoch of cometary activity commencing in August 2023 with an increase in brightness of >5.4 mag. The lightcurve of 2023 RN3 indicates the presence of continuous cometary activity across our observations, suggesting the onset of a new period of sustained activity. We find no evidence of any outbursts on top of the observed brightening, nor do we find any significant color evolution across our observations. 2023 RN3 is visibly extended in LCO and Liverpool Telescope observations, indicating the presence of a spatially-extended coma. Numerical integration of 2023 RN3's orbit reveals the comet to have recently undergone a slight increase in semimajor axis due to a planetary encounter with Jupiter, however whether this orbital change could trigger 2023 RN3's cometary activity is unclear. Our estimate for the maximum dust production metric of Afrho ~400 cm is consistent with previous measurements for the Jupiter-family comet and Centaur populations., Comment: 20 pages, 9 figures

Centaurs are small Solar System objects on chaotic orbits in the giant planet region, forming an evolutionary continuum with the Kuiper belt objects and Jupiter-family comets. Some Centaurs are known to exhibit cometary activity, though unlike comets this activity tends not to correlate with heliocentric distance and the mechanism behind it is currently poorly understood. We utilize serendipitous observations from the Asteroid Terrestrial-impact Last Alert System (ATLAS), Zwicky Transient Facility (ZTF), Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), Dark Energy Survey (DES), and Gaia in addition to targeted follow-up observations from the Las Cumbres Observatory, TRAnsiting Planets and PlanetesImals Small Telescope South (TRAPPIST-South), and Gemini North telescope to analyze an unexpected brightening exhibited by the known active Centaur (2060) Chiron in 2021. This is highly indicative of a cometary outburst. As of 2023 February, Chiron has still not returned to its pre-brightening magnitude. We find Chiron's rotational lightcurve, phase curve effects, and possible high-albedo surface features to be unlikely causes of this observed brightening. We consider the most likely cause to be an epoch of either new or increased cometary activity, though we cannot rule out a possible contribution from Chiron's reported ring system, such as a collision of as-yet unseen satellites shepherding the rings. We find no evidence for coma in our Gemini or TRAPPIST-South observations, though this does not preclude the possibility that Chiron is exhibiting a coma that is too faint for observation or constrained to the immediate vicinity of the nucleus., Comment: 39 pages, 14 figures, 14 tables. Has been accepted for publication in PSJ

We present the discovery and early observations of the nearby Type II supernova (SN) 2024ggi in NGC 3621 at 6.64 +/- 0.3 Mpc. The SN was caught 5.8 (+1.9 -2.9) hours after its explosion by the ATLAS survey. Early-phase, high-cadence, and multi-band photometric follow-up was performed by the Kinder (Kilonova Finder) project, collecting over 1000 photometric data points within a week. The combined o- and r-band light curves show a rapid rise of 3.3 magnitudes in 13.7 hours, much faster than SN 2023ixf (another recent, nearby, and well-observed SN II). Between 13.8 and 18.8 hours after explosion SN 2024ggi became bluer, with u-g colour dropping from 0.53 to 0.15 mag. The rapid blueward evolution indicates a wind shock breakout (SBO) scenario. No hour-long brightening expected for the SBO from a bare stellar surface was detected during our observations. The classification spectrum, taken 17 hours after the SN explosion, shows flash features of high-ionization species such as Balmer lines, He I, C III, and N III. Detailed light curve modeling reveals critical insights into the properties of the circumstellar material (CSM). Our favoured model has an explosion energy of 2 x 10^51 erg, a mass-loss rate of 10^-3 solar_mass/yr (with an assumed 10 km/s wind), and a confined CSM radius of 6 x 10^14 cm. The corresponding CSM mass is 0.4 solar_mass. Comparisons with SN 2023ixf highlight that SN 2024ggi has a smaller CSM density, resulting in a faster rise and fainter UV flux. The extensive dataset and the involvement of citizen astronomers underscore that a collaborative network is essential for SBO searches, leading to more precise and comprehensive SN characterizations., Comment: 11 pages, 5 figures in manuscript, 6 pages in appendix, submitted to ApJL

We report an analysis of a sample of 186 spectroscopically confirmed Type II supernova (SN) light curves (LCs) obtained from a combination of Zwicky Transient Facility (ZTF) and Asteroid Terrestrial-impact Last Alert System (ATLAS) observations. We implement a method to infer physical parameters from these LCs using hydrodynamic models that take into account the progenitor mass, the explosion energy, and the presence of circumstellar matter (CSM). The CSM is modelled via the mass loss rate, wind acceleration at the surface of the progenitor star with a $\beta$ velocity law, and the CSM radius. We also infer the time of explosion, attenuation (A$_V$), and the redshift for each SN. Our results favor low-mass progenitor stars (M$_{ZAMS}$\,$<$14\,$M_\odot$) with a dense CSM ($\dot{M}$ $>$ 10$^{-3}$ [M$_\odot$ yr$^{-1}$], a CSM radius of $\sim$ 10$^{15}$ cm, and $\beta$ $>$ 2). Additionally, we find that the redshift inferred from the supernova LCs is significantly more accurate than that inferred using the host galaxy photometric redshift, suggesting that this method could be used to infer more accurate host galaxy redshifts from large samples of SNe II in the LSST era. Lastly, we compare our results with similar works from the literature.

Measurements of Type Ia Supernovae (SNe Ia) in the near-infrared (NIR) have been used both as an alternate path to cosmology compared to optical measurements and as a method of constraining key systematics for the larger optical studies. With the DEHVILS sample, the largest published NIR sample with consistent NIR coverage of maximum light across three NIR bands ($Y$, $J$, and $H$), we check three key systematics: (i) the reduction in Hubble residual scatter as compared to the optical, (ii) the measurement of a "mass step" or lack thereof and its implications, and (iii) the ability to distinguish between various dust models by analyzing slopes and correlations between Hubble residuals in the NIR and optical. We produce SN Ia simulations of the DEHVILS sample and find that it is $\textit{harder}$ to differentiate between various dust models than previously understood. Additionally, we find that fitting with the current SALT3-NIR model does not yield accurate wavelength-dependent stretch-luminosity correlations, and we propose a limited solution for this problem. From the data, we see that (i) the standard deviation of Hubble residual values from NIR bands treated as standard candles are 0.007-0.042 mag smaller than those in the optical, (ii) the NIR mass step is not constrainable with the current sample size of 47 SNe Ia from DEHVILS, and (iii) Hubble residuals in the NIR and optical are correlated in the data. We test a few variations on the number and combinations of filters and data samples, and we observe that none of our findings or conclusions are significantly impacted by these modifications., Comment: 22 pages, 10 figures. Accepted by A&A

We introduce the Hawai`i Supernova Flows project and present summary statistics of the first 1218 astronomical transients observed, 669 of which are spectroscopically classified Type Ia Supernovae (SNe Ia). Our project is designed to obtain systematics-limited distances to SNe Ia while consuming minimal dedicated observational resources. This growing sample will provide increasing resolution into peculiar velocities as a function of position on the sky and redshift, allowing us to more accurately map the structure of dark matter. This can be used to derive cosmological parameters such as $\sigma_8$ and can be compared with large scale flow maps from other methods such as luminosity-line width or luminosity-velocity dispersion correlations in galaxies. Additionally, our photometry will provide a valuable test bed for analyses of SNe Ia incorporating near-infrared data. In this survey paper, we describe the methodology used to select targets, collect and reduce data, and calculate distances., Comment: 33 pages, 23 figures

We present the fourth data release (DR4) of the SkyMapper Southern Survey (SMSS), the last major step in our hemispheric survey with six optical filters: u, v, g, r, i, z. SMSS DR4 covers 26,000 sq.deg from over 400,000 images acquired by the 1.3m SkyMapper telescope between 2014-03 and 2021-09. The 6-band sky coverage extends from the South Celestial Pole to Dec = +16deg, with some images reaching Dec ~ +28deg. In contrast to previous DRs, we include all good-quality images from the facility taken during that time span, not only those explicitly taken for the public Survey. From the image dataset, we produce a catalogue of nearly 13 billion detections made from ~700 million unique astrophysical objects. The typical 10sigma depths for each field range between 18.5 and 20.5 mag, depending on the filter, but certain sky regions include longer exposures that reach as deep as 22 mag in some filters. As with previous SMSS catalogues, we have cross-matched with a host of other imaging and spectroscopic datasets to facilitate additional science outcomes. SMSS DR4 is now available to the worldwide astronomical community., Comment: 32 pages. Submitted to the Publications of the Astronomical Society of Australia

Jupiter-family Comets (JFCs) exhibit a wide range of activity levels and mass-loss over their orbits. We analyzed high-cadence observations of 42 active JFCs with the wide-field Asteroid Terrestrial-impact Last Alert System (ATLAS) survey in 2020-2021. We measured dust production rates of the JFCs using the Af\rho parameter and its variation as a function of heliocentric distance. There is a tendency for our JFC sample to exhibit a maximum Af\rho after perihelion, with 254P/McNaught and P/2020 WJ5 (Lemmon) having their maximum Af\rho over a year after perihelion. On average, the rate of change of activity post-perihelion was shallower than that pre-perihelion. We also estimated the mass maximum loss rate for 17 of the JFCs in our sample, finding 4P/Faye to be the most active. We present a subset of comets whose measured Af\rho have been interpolated and extrapolated to a common distance of 2 au pre-perihelion and post-perihelion. From these measurements we found no correlation of intrinsic activity with current perihelion distance. For three of the JFCs in our sample, 6P/d'Arrest, 156P/Russell-LINEAR and 254P/McNaught, there was no visible coma but a constant absolute magnitude which we attributed to a probable detection of the nucleus. We derived upper limits for the nuclear radii of \leq 2.1 +/- 0.3 km, \leq 2.0 +/- 0.2 km and \leq 4.0 +/- 0.8 km respectively. Finally, we found that 4P/Faye, 108P/Ciffreo, 132P/Helin-Roman-Alu 2, 141P/Machholz 2, and 398P/Boattini experienced outbursts between 2020 and 2022., Comment: Accepted for publication in The Planetary Science Journal (PSJ)

The Kuiper belt objects, the Centaurs, and the Jupiter-family comets form an evolutionary continuum of small outer Solar System objects, and their study allows us to gain insight into the history and evolution of the Solar System. Broadband photometry can be used to measure their phase curves, allowing a first-order probe into the surface properties of these objects, though limited telescope time makes measuring accurate phase curves difficult. We make use of serendipitous broadband photometry from the long-baseline, high-cadence ATLAS survey to measure the phase curves for a sample of 18 Kuiper belt objects, Centaurs, and Jupiter-family comets with unprecedentedly large datasets. We find phase curves with previously reported negative slopes become positive with increased data and are thus due to insufficient sampling of the phase curve profile, and not a real physical effect. We search for correlations between phase curve parameters, finding no strong correlations between any parameter pair, consistent with the findings of previous studies. We search for instances of cometary activity in our sample, finding a previously reported outburst by Echeclus and a new epoch of increased activity by Chiron. Applying the main belt asteroid HG1G2 phase curve model to three Jupiter-family comets in our sample with large phase angle spans, we find their slope parameters imply surfaces more consistent with those of carbonaceous main belt asteroids than silicaceous ones., Comment: 42 pages, 24 figures, accepted for publication to the Planetary Science Journal. Supplementary materials will be available in the online journal version

While the sample of optical Type Ia Supernova (SN Ia) light curves (LCs) usable for cosmological parameter measurements surpasses 2000, the sample of published, cosmologically viable near-infrared (NIR) SN Ia LCs, which have been shown to be good "standard candles," is still $\lesssim$ 200. Here, we present high-quality NIR LCs for 83 SNe Ia ranging from $0.002 < z < 0.09$ as a part of the Dark Energy, H$_0$, and peculiar Velocities using Infrared Light from Supernovae (DEHVILS) survey. Observations are taken using UKIRT's WFCAM, where the median depth of the images is 20.7, 20.1, and 19.3 mag (Vega) for $Y$, $J$, and $H$-bands, respectively. The median number of epochs per SN Ia is 18 for all three bands ($YJH$) combined and 6 for each band individually. We fit 47 SN Ia LCs that pass strict quality cuts using three LC models, SALT3, SNooPy, and BayeSN and find scatter on the Hubble diagram to be comparable to or better than scatter from optical-only fits in the literature. Fitting NIR-only LCs, we obtain standard deviations ranging from 0.128-0.135 mag. Additionally, we present a refined calibration method for transforming 2MASS magnitudes to WFCAM magnitudes using HST CALSPEC stars that results in a 0.03 mag shift in the WFCAM $Y$-band magnitudes., Comment: 24 pages, 9 figures. Accepted by MNRAS

We analyze high-cadence data from the Transiting Exoplanet Survey Satellite (TESS) of the ambiguous nuclear transient (ANT) ASASSN-18el. The optical changing-look phenomenon in ASASSN-18el has been argued to be due to either a drastic change in the accretion rate of the existing active galactic nucleus (AGN) or the result of a tidal disruption event (TDE). Throughout the TESS observations, short-timescale stochastic variability is seen, consistent with an AGN. We are able to fit the TESS light curve with a damped-random-walk (DRW) model and recover a rest-frame variability amplitude of $\hat{\sigma} = 0.93 \pm 0.02$ mJy and a rest-frame timescale of $\tau_{DRW} = 20^{+15}_{-6}$ days. We find that the estimated $\tau_{DRW}$ for ASASSN-18el is broadly consistent with an apparent relationship between the DRW timescale and central supermassive black hole mass. The large-amplitude stochastic variability of ASASSN-18el, particularly during late stages of the flare, suggests that the origin of this ANT is likely due to extreme AGN activity rather than a TDE., Comment: Updated to reflect the accepted version in MNRAS

We present the discovery that ATLAS18mlw was a tidal disruption event (TDE) in the galaxy WISEA J073544.83+663717.3, at a luminosity distance of 334 Mpc. Initially discovered by the Asteroid Terrestrial Impact Last Alert System (ATLAS) on 2018 March 17.3, the TDE nature of the transient was uncovered only recently with the re-reduction of a SuperNova Integral Field Spectrograph (SNIFS) spectrum. This spectrum, taken by the Spectral Classification of Astronomical Transients (SCAT) survey, shows a strong blue continuum and a broad H$\alpha$ emission line. Here we present roughly six years of optical survey photometry beginning before the TDE to constrain AGN activity, optical spectroscopy of the transient, and a detailed study of the host galaxy properties through analysis of archival photometry and a host spectrum. ATLAS18mlw was detected in ground-based light curves for roughly two months. From a blackbody fit to the transient spectrum and bolometric correction of the optical light curve, we conclude that ATLAS18mlw is best explained by a low-luminosity TDE with a peak luminosity of log(L [erg s$^{-1}$]) = $43.5 \pm 0.2$. The TDE classification is further supported by the quiescent Balmer strong nature of the host galaxy. We also calculated the TDE decline rate from the bolometric light curve and find $\Delta L_{40} = -0.7 \pm 0.2$ dex, making ATLAS18mlw a member of the growing class of ``faint and fast'' TDEs with low peak luminosities and fast decline rates., Comment: Updated to reflect the accepted version in MNRAS

We present a novel technique to study Type Ia supernovae by constraining surviving companions of historical extragalactic SN by combining archival photographic plates and Hubble Space Telescope imaging. We demonstrate this technique for Supernova 1972E, the nearest known SN Ia in over 125 years. Some models of Type Ia supernovae describe a white dwarf with a non-degenerate companion that donates enough mass to trigger thermonuclear detonation. Hydrodynamic simulations and stellar evolution models show that these donor stars survive the explosion, and show increased luminosity for at least a thousand years. Thus, late-time observations of the exact location of a supernova after its ejecta have faded can constrain the presence of a surviving donor star and progenitor models. We find the explosion site of SN 1972E by analyzing 17 digitized photographic plates taken with the European Southern Observatory 1m Schmidt and 1 plate taken with the Cerro Tololo Inter-American Observatory 1.5m telescope from 1972-1974. Using the \textit{Gaia} eDR3 catalog to determine Supernova 1972E's equatorial coordinates yields: $\alpha$ = 13$^h$ 39$^m$ 52.708$^s$ $\pm$ 0.004$^s$ and $\delta$ = $-$31\degree 40' 8\farcs97 $\pm$ 0\farcs04 (ICRS). In 2005, HST/ACS imaged NGC 5253, the host galaxy of SN 1972E, with the $F435W$, $F555W$, and $F814W$ filters covering the explosion site. The nearest source detected is offset by 3.0 times our positional precision, and is inconsistent with the colors expected of a surviving donor star. Thus, the 2005 HST observation rules out all Helium-star companion models, and the most luminous main-sequence companion model currently in the literature. The remaining main-sequence companion models could be tested with deeper HST imaging., Comment: 13 pages, 5 figures, and 2 tables in main text, plus 2 pages of appendices. Submitted to MNRAS

We estimate the total population of near-Earth objects (NEOs) in the Solar System, using an extensive, `Solar System to pixels' fake-asteroid simulation to debias detections of real NEOs by the ATLAS survey. Down to absolute magnitudes $H=25$ and 27.6 (diameters of $\sim 34$ and 10 meters, respectively, for 15% albedo), we find total populations of $(3.72 \pm 0.49) \times 10^5$ and $(1.59 \pm 0.45) \times 10^7$ NEOs, respectively. Most plausible sources of error tend toward underestimation, so the true populations are likely larger. We find the distribution of $H$ magnitudes steepens for NEOs fainter than $H \sim 22.5$, making small asteroids more common than extrapolation from brighter $H$ mags would predict. Our simulation indicates a strong bias against detecting small but dangerous asteroids that encounter Earth with high relative velocities -- i.e., asteroids in highly inclined and/or eccentric orbits. Worldwide NEO discovery statistics indicate this bias affects global NEO detection capability, to the point that an observational census of small asteroids in such orbits is probably not currently feasible. Prompt and aggressive followup of NEO candidates, combined with closer collaborations between segments of the global NEO community, can increase detection rates for these dangerous objects., Comment: 23 pages, 11 figures, accepted by the Planetary Science Journal

We present the discovery that ASASSN-14ko is a periodically flaring AGN at the center of the galaxy ESO 253-G003. At the time of its discovery by the All-Sky Automated Survey for Supernovae (ASAS-SN), it was classified as a supernova close to the nucleus. The subsequent six years of V- and g-band ASAS-SN observations reveal that ASASSN-14ko has nuclear flares occurring at regular intervals. The seventeen observed outbursts show evidence of a decreasing period over time, with a mean period of $P_0 = 114.2 \pm 0.4$ days and a period derivative of $\dot{P} = -0.0017\pm0.0003$. The most recent outburst in May 2020, which took place as predicted, exhibited spectroscopic changes during the rise and a had a UV bright, blackbody spectral energy distribution similar to tidal disruption events (TDEs). The X-ray flux decreased by a factor of 4 at the beginning of the outburst and then returned to its quiescent flux after ~8 days. TESS observed an outburst during Sectors 4-6, revealing a rise time of $5.60 \pm 0.05$ days in the optical and a decline that is best fit with an exponential model. We discuss several possible scenarios to explain ASASSN-14ko's periodic outbursts, but currently favor a repeated partial TDE. The next outbursts should peak in the optical on UT 2020-09-7.4$ \pm $1.1 and UT 2020-12-26.5$ \pm $1.4., Comment: 26 pages, 15 figures, 7 tables. Will be submitted to ApJ. The latest flare is currently ongoing, as we predicted

A large sample of white dwarfs is selected by both proper motion and colours from the Pan-STARRS 1 3{\pi} Steradian Survey Processing Version 2 to construct the White Dwarf Luminosity Functions of the discs and halo in the solar neighbourhood. Four-parameter astrometric solutions were recomputed from the epoch data. The generalised maximum volume method is then used to calculate the density of the populations. After removal of crowded areas near the Galactic plane and centre, the final sky area used by this work is 7.833 sr, which is 83% of the 3{\pi} sky and 62% of the whole sky. By dividing the sky using Voronoi tessellation, photometric and astrometric uncertainties are recomputed at each step of the integration to improve the accuracy of the maximum volume. Interstellar reddening is considered throughout the work. We find a disc-to-halo white dwarf ratio of about 100., Comment: 18 pages, 13 figures, 2 tables in text, 4 tables in Appendix A, supplementary online materials are available on MNRAS. Accepted for publication on 3 Oct 2018

We report the discovery of the ultra-luminous QSO SMSS~J215728.21-360215.1 with magnitude $z=16.9$ and W4$=7.42$ at redshift 4.75. Given absolute magnitudes of $M_{145,\rm AB}=-29.3$, $M_{300,\rm AB}=-30.12$ and $\log L_{\rm bol}/L_{\rm bol,\odot} = 14.84$, it is the QSO with the highest unlensed UV-optical luminosity currently known in the Universe. It was found by combining proper-motion data from Gaia DR2 with photometry from SkyMapper DR1 and the Wide-field Infrared Survey Explorer (WISE). In the Gaia database it is an isolated single source and thus unlikely to be strongly gravitationally lensed. It is also unlikely to be a beamed source as it is not discovered in the radio domain by either NVSS or SUMSS. It is classed as a weak-emission-line QSO and possesses broad absorption line features. A lightcurve from ATLAS spanning the time from October 2015 to December 2017 shows little sign of variability., Comment: 6 pages, 4 figures, accepted for publication in PASA

The Asteroid Terrestrial-impact Last Alert System (ATLAS) carries out its primary planetary defense mission by surveying about 13000 deg^2 at least four times per night. The resulting data set is useful for the discovery of variable stars to a magnitude limit fainter than r~18, with amplitudes down to 0.01 mag for bright objects. Here we present a Data Release One catalog of variable stars based on analyzing 142 million stars measured at least 100 times in the first two years of ATLAS operations. Using a Lomb-Scargle periodogram and other variability metrics, we identify 4.7 million candidate variables which we analyze in detail. Through Space Telescope Science Institute, we publicly release lightcurves for all of them, together with a vector of 169 classification features for each star. We do this at the level of unconfirmed candidate variables in order to provide the community with a large set of homogeneously analyzed photometry and avoid pre-judging which types of objects others may find most interesting. We use machine learning to classify the candidates into fifteen different broad categories based on lightcurve morphology. About 10% (430,000 stars) pass extensive tests designed to screen out spurious variability detections: we label these as `probable' variables. Of these, 230,000 receive specific classifications as eclipsing binaries, pulsating, Mira-type, or sinusoidal variables: these are the `classified' variables. New discoveries among the probable variables number more than 300,000, while 150,000 of the classified variables are new, including about 10,000 pulsating variables, 2,000 Mira stars, and 70,000 eclipsing binaries., Comment: Accepted by AJ; gives instructions for querying ATLAS variable star database; this new version has nicer lightcurve figures

We characterize the spatial density of the Pan-STARRS1 (PS1) sample of RR Lyrae stars, to study the properties of the old Galactic stellar halo as traced by RRab stars. This sample of 44,403 sources spans Galactocentric radii of $0.55 \; \mathrm{kpc} \leq R_{\mathrm{gc}} \leq 141 \; \mathrm{kpc}$ with a distance precision of 3\% and thus is able to trace the halo out to larger distances than most previous studies. After excising stars that are attributed to dense regions such as stellar streams, the Galactic disc and bulge as well as halo globular clusters, the sample contains ${\sim}11,000$ sources within $20 \; \mathrm{kpc} \leq R_{\mathrm{gc}} \leq 131 \; \mathrm{kpc}$. We then apply forward modeling using ellipsoidal stellar density models $\rho(l,b,R_{\mathrm{gc}})$ both with a constant and a radius-dependent halo flattening $q(R_{\mathrm{gc}})$. Assuming constant flattening $q$, the distribution of the sources is reasonably well fit from $20 \; \mathrm{kpc}$ to $131 \; \mathrm{kpc}$ by a single power law with $n=4.40^{+0.05}_{-0.04}$ and $q=0.918^{+0.016}_{-0.014}$. The distance distribution is fit comparably well by an Einasto profile with $n=9.53^{+0.27}_{-0.28}$, an effective radius $r_{\mathrm{eff}}=1.07 \pm 0.10 \; \mathrm{kpc}$ and a halo flattening of $q=0.923 \pm 0.007$. If we allow for a radius-dependent flattening $q(R_{\mathrm{gc}})$, we find evidence for a distinct flattening of $q{\sim}0.8$ of the inner halo at ${\sim} 25 \; \mathrm{kpc}$. Additionally, we find that the south Galactic hemisphere is more flattened than the north Galactic hemisphere. The results of our work are largely consistent with many earlier results, e.g. \cite{Watkins2009}, \cite{Iorio2017}. We find that the stellar halo, as traced in RR Lyrae stars, exhibits a substantial number of further significant over- and underdensities, even after all known overdensities have been masked.

We present a new 3D map of interstellar dust reddening, covering three quarters of the sky (declinations greater than -30 degrees) out to a distance of several kiloparsecs. The map is based on high-quality stellar photometry of 800 million stars from Pan-STARRS 1 and 2MASS. We divide the sky into sightlines containing a few hundred stars each, and then infer stellar distances and types, along with the line-of-sight dust distribution. Our new map incorporates a more accurate average extinction law and an additional 1.5 years of Pan-STARRS 1 data, tracing dust to greater extinctions and at higher angular resolutions than our previous map. Out of the plane of the Galaxy, our map agrees well with 2D reddening maps derived from far-infrared dust emission. After accounting for a 15% difference in scale, we find a mean scatter of 10% between our map and the Planck far-infrared emission-based dust map, out to a depth of 0.8 mag in E(r-z), with the level of agreement varying over the sky. Our map can be downloaded at http://argonaut.skymaps.info, or by its DOI: 10.7910/DVN/LCYHJG., Comment: Submitted to MNRAS. 17 pages, 16 figures

A detection of the predicted anticorrelation between 21cm and either Ly-alpha or H-alpha from the Epoch of Reionization (EOR) would be a powerful probe of the first galaxies. While 3D intensity maps isolate foregrounds in low k_\parallel modes, infrared surveys cannot yet match the field of view and redshift resolution of radio intensity mapping experiments. In contrast, 2D (i.e., broad band) infrared intensity maps can be measured with current experiments and are limited by foregrounds instead of photon or thermal noise. We show 2D experiments can measure most of the 3D fluctuation power at k<0.2 Mpc^-1 while preserving its correlation properties. However, we show foregrounds pose two challenges: (1) simple geometric effects produce percent-level correlations between radio and infrared fluxes, even if their luminosities are uncorrelated; and (2) radio and infrared foreground residuals contribute sample variance noise to the cross spectrum. The first challenge demands better foreground masking and subtraction, while the second demands large fields of view to average away uncorrelated radio and infrared power. Using radio observations from the Murchison Widefield Array and near-infrared observations from the Asteroid Terrestrial-impact Last Alert System, we set an upper limit on residual foregrounds of the 21cm--Ly-alpha cross power spectrum at z\sim7 of \Delta^2<181 kJy/sr * mK (95\%) at \ell\sim800. We predict levels of foreground correlation and sample variance noise in future experiments, showing that higher resolution surveys such as LOFAR, SKA-LOW, and the Dark Energy Survey can start to probe models of the 21cm--Ly\alpha EOR cross spectrum., Comment: 24 pages, 13 f igures

Large modern surveys require efficient review of data in order to find transient sources such as supernovae, and to distinguish such sources from artefacts and noise. Much effort has been put into the development of automatic algorithms, but surveys still rely on human review of targets. This paper presents an integrated system for the identification of supernovae in data from Pan-STARRS1, combining classifications from volunteers participating in a citizen science project with those from a convolutional neural network. The unique aspect of this work is the deployment, in combination, of both human and machine classifications for near real-time discovery in an astronomical project. We show that the combination of the two methods outperforms either one used individually. This result has important implications for the future development of transient searches, especially in the era of LSST and other large-throughput surveys., Comment: 10 pages, 9 figures, submitted to MNRAS

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., Comment: 19 pages, 13 figures, 5 tables, accepted to AJ. The the PS1 catalog of RR Lyrae stars will become publicly available on Nov 1 2017. For collaborations on projects and earlier access to the catalog, please contact the first author

Although the majority of Centaurs are thought to have originated in the scattered disk, with the high-inclination members coming from the Oort cloud, the origin of the high inclination component of trans-Neptunian objects (TNOs) remains uncertain. We report the discovery of a retrograde TNO, which we nickname "Niku", detected by the Pan-STARRS 1 Outer Solar System Survey. Our numerical integrations show that the orbital dynamics of Niku are very similar to that of 2008 KV$_{42}$ (Drac), with a half-life of $\sim 500$ Myr. Comparing similar high inclination TNOs and Centaurs ($q > 10$ AU, $a < 100$ AU and $i > 60^\circ$), we find that these objects exhibit a surprising clustering of ascending node, and occupy a common orbital plane. This orbital configuration has high statistical significance: 3.8-$\sigma$. An unknown mechanism is required to explain the observed clustering. This discovery may provide a pathway to investigate a possible reservoir of high-inclination objects., Comment: 18 pages, 4 figures, 1 table, accepted for publication in ApJ Letters

Substellar members of young ($\lesssim$150 Myr) moving groups are valuable benchmarks to empirically define brown dwarf evolution with age and to study the low-mass end of the initial mass function. We have combined Pan-STARRS1 (PS1) proper motions with optical$-$IR photometry from PS1, 2MASS and $\textit{WISE}$ to search for substellar members of the AB Dor Moving Group within $\approx$50 pc and with spectral types of late-M to early-L, corresponding to masses down to $\approx$30 M$_{Jup}$ at the age of the group ($\approx$125 Myr). Including both photometry and proper motions allows us to better select candidates by excluding field dwarfs whose colors are similar to young AB~Dor Moving Group members. Our near-IR spectroscopy has identified six ultracool dwarfs (M6$-$L4; $\approx$30$-$100 M$_{Jup}$) with intermediate surface gravities (INT-G) as candidate members of the AB Dor Moving Group. We find another two candidate members with spectra showing hints of youth but consistent with field gravities. We also find four field brown dwarfs unassociated with the AB Dor Moving Group, three of which have INT-G gravity classification. While signatures of youth are present in the spectra of our $\approx$125 Myr objects, neither their $J-K$ nor $W1-W2$ colors are significantly redder than field dwarfs with the same spectral types, unlike younger ultracool dwarfs. We also determined PS1 parallaxes for eight of our candidates and one previously identified AB Dor Moving Group candidate. Although radial velocities (and parallaxes, for some) are still needed to fully assess membership, these new objects provide valuable insight into the spectral characteristics and evolution of young brown dwarfs., Comment: ApJ, accepted

The ESSENCE survey discovered 213 Type Ia supernovae at redshifts 0.1 < z < 0.81 between 2002 and 2008. We present their R and I-band photometry, measured from images obtained using the MOSAIC II camera at the CTIO 4 m Blanco telescope, along with rapid-response spectroscopy for each object. We use our spectroscopic follow-up observations to determine an accurate, quantitative classification and a precise redshift. Through an extensive calibration program we have improved the precision of the CTIO Blanco natural photometric system. We use several empirical metrics to measure our internal photometric consistency and our absolute calibration of the survey. We assess the effect of various potential sources of systematic bias on our measured fluxes, and we estimate that the dominant term in the systematic error budget from the photometric calibration on our absolute fluxes is ~1%., Comment: (40 pages, 21 figures, accepted for publication in ApJS)

We present a recalibration of the Sloan Digital Sky Survey (SDSS) photometry with new flat fields and zero points derived from Pan-STARRS1 (PS1). Using PSF photometry of 60 million stars with $16 < r < 20$, we derive a model of amplifier gain and flat-field corrections with per-run RMS residuals of 3 millimagnitudes (mmag) in $griz$ bands and 15 mmag in $u$ band. The new photometric zero points are adjusted to leave the median in the Galactic North unchanged for compatibility with previous SDSS work. We also identify transient non-photometric periods in SDSS ("contrails") based on photometric deviations co-temporal in SDSS bands. The recalibrated stellar PSF photometry of SDSS and PS1 has an RMS difference of {9,7,7,8} mmag in $griz$, respectively, when averaged over $15'$ regions., Comment: 13 pages, 7 figures, ApJ, in press. "Hypercalibration" refers to using repeat measurements of many stars from multiple surveys to constrain calibration parameters

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., Comment: Published in ApJ

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\"otes 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., Comment: 6 pages, 4 figures, Accepted for publication by ApjL

Supermassive black hole binaries (SMBHBs) should be an inevitable consequence of the hierarchical growth of massive galaxies through mergers, and the strongest sirens of gravitational waves (GWs) in the cosmos. And yet, their direct detection has remained elusive due to the compact (sub-parsec) orbital separations of gravitationally bound SMBHBs. Here we exploit a theoretically predicted signature of a SMBHB in the time domain: periodic variability caused by a mass accretion rate that is modulated by the binary's orbital motion. We report our first significant periodically varying quasar detection from the systematic search in the Pan-STARRS1 (PS1) Medium Deep Survey. Our SMBHB candidate, PSO J334.2028+01.4075, is a luminous radio-loud quasar at $z=2.060$, with extended baseline photometry from the Catalina Real-Time Transient Survey, as well as archival spectroscopy from the FIRST Bright Quasar Survey. The observed period ($542 \pm 15$ days) and estimated black hole mass ($\log (M_{\rm BH}/M_\odot) = 9.97 \pm 0.50$), correspond to an orbital separation of $7^{+8}_{-4}$ Schwarzschild radii ($\sim 0.006^{+0.007}_{-0.003}$ pc), assuming the rest-frame period of the quasar variability traces the orbital period of the binary. This SMBHB candidate, discovered at the peak redshift for SMBH mergers, is in a physically stable configuration for a circumbinary accretion disk, and within the regime of GW-driven orbital decay. Our search with PS1 is a benchmark study for the exciting capabilities of LSST, which will have orders of magnitude larger survey power, and will potentially pinpoint the locations of thousands of SMBHBs in the variable night sky., Comment: Published in ApJ Letters; matches published version

We analyze 760475 observations of 333026 main-belt objects obtained by the Pan-STARRS1 (PS1) survey telescope between 2012 May 20 and 2013 November 9, a period during which PS1 discovered two main-belt comets, P/2012 T1 (PANSTARRS) and P/2013 R3 (Catalina-PANSTARRS). PS1 comet detection procedures currently consist of the comparison of the point spread functions (PSFs) of moving objects to those of reference stars, and the flagging of objects that show anomalously large radial PSF widths. Based on the number of missed discovery opportunities among comets discovered by other observers, we estimate an upper limit comet discovery efficiency rate of ~70% for PS1. Additional analyses that could improve comet discovery yields in future surveys include linear PSF analysis, modeling of trailed stellar PSFs for comparison to trailed moving object PSFs, searches for azimuthally localized activity, comparison of point-source-optimized photometry to extended-source-optimized photometry, searches for photometric excesses in objects with known absolute magnitudes, and crowd-sourcing. Analysis of PS1 survey statistics indicates an expected fraction of 59 MBCs per 10^6 outer main-belt asteroids, and a 95% confidence upper limit of 96 MBCs per 10^6 outer main-belt asteroids. We note that more sensitive future surveys could detect many more MBCs than estimated here. We find an excess of high eccentricities (0.1 < e < 0.3) among all known MBCs relative to the background asteroid population. Theoretical calculations show that, given these eccentricities, the sublimation rate for a typical MBC is orders of magnitude larger at perihelion than at aphelion, providing a plausible physical explanation for the observed behavior of MBCs peaking in observed activity strength near perihelion., Comment: 38 pages, 14 figures; accepted for publication in Icarus

We present Keck II/DEIMOS spectroscopy of the three distant dwarf galaxies of M31 Lacerta I, Cassiopeia III, and Perseus I, recently discovered within the Pan-STARRS1 3\pi imaging survey. The systemic velocities of the three systems (v_{r,helio} = -198.4 +/- 1.1 km/s, -371.6 +/- 0.7 km/s, and -326 +/- 3 km/s, respectively) confirm that they are satellites of M31. In the case of Lacerta I and Cassiopeia III, the high quality of the data obtained for 126 and 212 member stars, respectively, yields reliable constraints on their global velocity dispersions (\sigma_{vr} = 10.3 +/- 0.9 km/s and 8.4 +/- 0.6 km/s, respectively), leading to dynamical-mass estimates for both of ~4x10^7 Msun within their half-light radius. These translate to V-band mass-to-light ratios of 15^{+12}_{-9} and 8^{+9}_{-5} in solar units. We also use our spectroscopic data to determine the average metallicity of the 3 dwarf galaxies ([Fe/H] = -2.0 +/- 0.1, -1.7 +/- 0.1, and -2.0 +/- 0.2, respectively). All these properties are typical of dwarf galaxy satellites of Andromeda with their luminosity and size., Comment: 6 pages, 4 figures. ApJL, in press

We present the discovery of 61 wide (>5 arcsecond) separation, low-mass (stellar and substellar) companions to stars in the solar neighborhood identified from Pan-STARRS\,1 (PS1) data and the spectral classification of 27 previously known companions. Our companions represent a selective subsample of promising candidates and span a range in spectral type of K7-L9 with the addition of one DA white dwarf. These were identified primarily from a dedicated common proper motion search around nearby stars, along with a few as serendipitous discoveries from our Pan-STARRS1 brown dwarf search. Our discoveries include 24 new L dwarf companions and one known L dwarf not previously identified as a companion. The primary stars around which we searched for companions come from a list of bright stars with well-measured parallaxes and large proper motions from the Hipparcos catalog (8583 stars, mostly A-K~dwarfs) and fainter stars from other proper motion catalogues (79170 stars, mostly M~dwarfs). We examine the likelihood that our companions are chance alignments between unrelated stars and conclude that this is unlikely for the majority of the objects that we have followed-up spectroscopically. We also examine the entire population of ultracool (>M7) dwarf companions and conclude that while some are loosely bound, most are unlikely to be disrupted over the course of $\sim$10 Gyr. Our search increases the number of ultracool M dwarf companions wider than 300 AU by 88% and increases the number of L dwarf companions in the same separation range by 96%. Finally, we resolve our new L dwarf companion to HIP 6407 into a tight (0.13 arcsecond, 7.4 AU) L1+T3 binary, making the system a hierarchical triple. Our search for these key benchmarks against which brown dwarf and exoplanet atmosphere models are tested has yielded the largest number of discoveries to date., Comment: 74 pages, 17 figures, 13 tables, accepted to ApJ, updated with corrected version of Table 13

We present a panoptic view of the stellar structure in the Galactic disk's outer reaches commonly known as the Monoceros Ring, based on data from Pan-STARRS1. These observations clearly show the large extent of the stellar overdensities on both sides of the Galactic disk, extending between b = -25 and b = +35 degrees and covering over 130 degrees in Galactic longitude. The structure exhibits a complex morphology with both stream-like features and a sharp edge to the structure in both the north and the south. We compare this map to mock observations of two published simulations aimed at explaining such structures in the outer stellar disk, one postulating an origin as a tidal stream and the other demonstrating a scenario where the disk is strongly distorted by the accretion of a satellite. These morphological comparisons of simulations can link formation scenarios to observed structures, such as demonstrating that the distorted-disk model can produce thin density features resembling tidal streams. Although neither model produces perfect agreement with the observations--the tidal stream predicts material at larger distances which is not detected while in the distorted disk model the midplane is warped to an excessive degree--future tuning of the models to accommodate these latest data may yield better agreement., Comment: 10 pages, 7 figures. Accepted to ApJ. Key figures and data can be found at http://dept.astro.lsa.umich.edu/~ctslater/MRi/

We present the fiducial sequences of a sample of Galactic star clusters in the five bands of the Pan-STARRS1 (PS1) photometric system (g_P1, r_P1, i_P1, z_P1, y_P1). These empirical sequences -- which include the red giant and sub-giant branches, the main sequence, and the horizontal branch -- were defined from deep colour-magnitude diagrams reaching below the oldest main-sequence turn-offs of 13 globular and 3 old open clusters covering a wide range of metallicities (-2.4 < [Fe/H] < +0.4). We find excellent agreement for the nine clusters in common with previous studies in similar photometric systems when transformed to the PS1 system. Because the photometric and spectroscopic properties of these stellar populations are accurately known, the fiducials provide a solid basis for the interpretation of observations in the PS1 system, as well as valuable constraints to improve the empirical colour--$T_{eff}$ relations., Comment: 6 pages, 5 figures (+15 figures in Appendix). Re-submitted to MNRAS after addressing the referee's comments. Some figures degraded due to size limitations

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{\deg} 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{\deg},+32{\deg}), 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., Comment: 5 pages, 3 figures. MNRAS, in press

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., Comment: 6 pages, 4 figures, accepted

The NGVS-IR project (Next Generation Virgo Survey - Infrared) is a contiguous near-infrared imaging survey of the Virgo cluster of galaxies. It complements the optical wide-field survey of Virgo (NGVS). The current state of NGVS-IR consists of Ks-band imaging of 4 deg^2 centered on M87, and J and Ks-band imaging of 16 deg^2 covering the region between M49 and M87. In this paper, we present the observations of the central 4 deg^2 centered on Virgo's core region. The data were acquired with WIRCam on the Canada-France-Hawaii Telescope and the total integration time was 41 hours distributed in 34 contiguous tiles. A survey-specific strategy was designed to account for extended galaxies while still measuring accurate sky brightness within the survey area. The average 5\sigma limiting magnitude is Ks=24.4 AB mag and the 50% completeness limit is Ks=23.75 AB mag for point source detections, when using only images with better than 0.7" seeing (median seeing 0.54"). Star clusters are marginally resolved in these image stacks, and Virgo galaxies with \mu_Ks=24.4 AB mag arcsec^-2 are detected. Combining the Ks data with optical and ultraviolet data, we build the uiK color-color diagram which allows a very clean color-based selection of globular clusters in Virgo. This diagnostic plot will provide reliable globular cluster candidates for spectroscopic follow-up campaigns needed to continue the exploration of Virgo's photometric and kinematic sub-structures, and will help the design of future searches for globular clusters in extragalactic systems. Equipped with this powerful new tool, future NGVS-IR investigations based on the uiK diagram will address the mapping and analysis of extended structures and compact stellar systems in and around Virgo galaxies., Comment: 23 pages, 18 figures. Accepted for publication in ApJS

We present the discovery of a new dwarf galaxy, Perseus I/Andromeda XXXIII, found in the vicinity of Andromeda (M31) in stacked imaging data from the Pan-STARRS1 3{\pi} survey. Located 27.9{\deg} away from M31, Perseus I has a heliocentric distance of 785 +/- 65 kpc, compatible with it being a satellite of M31 at 374 +14/-10 kpc from its host. The properties of Perseus I are typical for a reasonably bright dwarf galaxy (M_V = -10.3 +/- 0.7), with an exponential half-light radius of r_h = 1.7 +/- 0.4 arcminutes or r_h = 400 +105/-85 pc at this distance, and a moderate ellipticity (\epsilon = 0.43 +0.15/-0.17). The late discovery of Perseus I is due to its fairly low surface brightness (\mu_0=25.7 +1.0/-0.9 mag/arcsec^2), and to the previous lack of deep, high quality photometric data in this region. If confirmed to be a companion of M31, the location of Perseus I, far east from its host, could place interesting constraints on the bulk motion of the satellite system of M31., Comment: 6 pages, 4 figures. Accepted for publication in ApJL

We report the discovery of two new dwarf galaxies, Lacerta I/Andromeda XXXI (Lac I/And XXXI) and Cassiopeia III/Andromeda XXXII (Cas III/And XXXII), in stacked Pan-STARRS1 r_P1- and i_P1-band imaging data. Both are luminous systems (M_V ~ -12) located at projected distances of 20.3{\deg} and 10.5{\deg} from M31. Lac I and Cas III are likely satellites of the Andromeda galaxy with heliocentric distances of 756^{+44}_{-28} kpc and 772^{+61}_{-56} kpc, respectively, and corresponding M31-centric distances of 275+/-7 kpc and 144^{+6}_{-4} kpc . The brightest of recent Local Group member discoveries, these two new dwarf galaxies owe their late discovery to their large sizes (r_h = 4.2^{+0.4}_{-0.5} arcmin or 912^{+124}_{-93} pc for Lac I; r_h = 6.5^{+1.2}_{-1.0} arcmin or 1456+/-267 pc for Cas III), and consequently low surface brightness (\mu_0 ~ 26.0 mag/arcsec^2), as well as to the lack of a systematic survey of regions at large radii from M31, close to the Galactic plane. This latter limitation is now alleviated by the 3{\pi} Pan-STARRS1 survey, which could lead to the discovery of other distant Andromeda satellite dwarf galaxies., Comment: 7 pages, 5 figures. Accepted for publication in ApJ

We describe the Pan-STARRS Moving Object Processing System (MOPS), a modern software package that produces automatic asteroid discoveries and identifications from catalogs of transient detections from next-generation astronomical survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from a synthetic but realistic population of asteroids whose measurements were simulated for a Pan-STARRS4-class telescope. Additionally, using a non-physical grid population, we demonstrate that MOPS can detect populations of currently unknown objects such as interstellar asteroids. MOPS has been adapted successfully to the prototype Pan-STARRS1 telescope despite differences in expected false detection rates, fill-factor loss and relatively sparse observing cadence compared to a hypothetical Pan-STARRS4 telescope and survey. MOPS remains >99.5% efficient at detecting objects on a single night but drops to 80% efficiency at producing orbits for objects detected on multiple nights. This loss is primarily due to configurable MOPS processing limits that are not yet tuned for the Pan-STARRS1 mission. The core MOPS software package is the product of more than 15 person-years of software development and incorporates countless additional years of effort in third-party software to perform lower-level functions such as spatial searching or orbit determination. We describe the high-level design of MOPS and essential subcomponents, the suitability of MOPS for other survey programs, and suggest a road map for future MOPS development., Comment: 57 Pages, 26 Figures, 13 Tables

We report the discovery of a wide (approximately 400 AU projected separation), common proper motion companion to the nearby M dwarf LHS 2803 (PSO J207.0300-13.7422). This object was discovered during our census of the local T dwarf population using Pan-STARRS1 and 2MASS data. Using IRTF/SpeX near-infrared spectroscopy, we classify the secondary to be spectral type T5.5. University of Hawai`i 2.2m/SNIFS optical spectroscopy indicates the primary has a spectral type of M4.5, with approximately solar metallicity and no measurable H_alpha emission. We use this lack of activity to set a lower age limit for the system of 3.5 Gyr. Using a comparison with chance alignments of brown dwarfs and nearby stars, we conclude that the two objects are unlikely to be a chance association. The primary's photometric distance of 21 pc and its proper motion implies thin disk kinematics. Based on these kinematics and its metallicity, we set an upper age limit for the system of 10 Gyr. Evolutionary model calculations suggest the secondary has a mass of 72^+4_-7 M_Jup, temperature of 1120+-80 K, and log g=5.4+-0.1 dex. Model atmosphere fitting to the near-IR spectrum gives similar physical parameters of 1100 K and log g=5.0., Comment: Accepted to ApJ (1) Max Planck Institute for Astronomy, (2) Institute for Astronomy, University of Hawai`i, Manoa, (3) Visiting Astronomer at the Infrared Telescope Facility, (4) Facebook, (5) Institute for Astronomy, University of Hawai`i, Hilo, (6) Princeton University Observatory

Photometric calibration is currently the dominant source of systematic uncertainty in exploiting type Ia supernovae to determine the nature of the dark energy. We review our ongoing program to address this calibration challenge by performing measurements of both the instrumental response function and the optical transmission function of the atmosphere. A key aspect of this approach is to complement standard star observations by using NIST-calibrated photodiodes as a metrology foundation for optical flux measurements. We present our first attempt to assess photometric consistency between synthetic photometry and observations, by comparing predictions based on a NIST-diode-based determination of the PanSTARRS-1 instrumental response and empirical atmospheric transmission measurements, with fluxes we obtained from observing spectrophotometric standards., Comment: 11 pages, submitted to Proceedings of the Calibration and Standardization of Large Surveys and Missions in Astronomy and Astrophysics Conference, held at Fermilab April 16-19 2012

In this first paper in a series we present 1298 low-redshift (z\leq0.2) optical spectra of 582 Type Ia supernovae (SNe Ia) observed from 1989 through 2008 as part of the Berkeley SN Ia Program (BSNIP). 584 spectra of 199 SNe Ia have well-calibrated light curves with measured distance moduli, and many of the spectra have been corrected for host-galaxy contamination. Most of the data were obtained using the Kast double spectrograph mounted on the Shane 3 m telescope at Lick Observatory and have a typical wavelength range of 3300-10,400 Ang., roughly twice as wide as spectra from most previously published datasets. We present our observing and reduction procedures, and we describe the resulting SN Database (SNDB), which will be an online, public, searchable database containing all of our fully reduced spectra and companion photometry. In addition, we discuss our spectral classification scheme (using the SuperNova IDentification code, SNID; Blondin & Tonry 2007), utilising our newly constructed set of SNID spectral templates. These templates allow us to accurately classify our entire dataset, and by doing so we are able to reclassify a handful of objects as bona fide SNe Ia and a few other objects as members of some of the peculiar SN Ia subtypes. In fact, our dataset includes spectra of nearly 90 spectroscopically peculiar SNe Ia. We also present spectroscopic host-galaxy redshifts of some SNe Ia where these values were previously unknown. [Abridged], Comment: 34 pages, 11 figures, 11 tables, revised version, re-submitted to MNRAS. Spectra will be released in January 2013. The SN Database homepage (http://hercules.berkeley.edu/database/index_public.html) contains the full tables, plots of all spectra, and our new SNID templates