The Asteroid Terrestrial impact Last Alert System (ATLAS) system consists of two 0.5m Schmidt telescopes with cameras covering 29 square degrees at plate scale of 1.86 arcsec per pixel. Working in tandem, the telescopes routinely survey the whole sky visible from Hawaii (above $\delta > -50^{\circ}$) every two nights, exposing four times per night, typically reaching $o < 19$ magnitude per exposure when the moon is illuminated and $c < 19.5$ per exposure in dark skies. Construction is underway of two further units to be sited in Chile and South Africa which will result in an all-sky daily cadence from 2021. Initially designed for detecting potentially hazardous near earth objects, the ATLAS data enable a range of astrophysical time domain science. To extract transients from the data stream requires a computing system to process the data, assimilate detections in time and space and associate them with known astrophysical sources. Here we describe the hardware and software infrastructure to produce a stream of clean, real, astrophysical transients in real time. This involves machine learning and boosted decision tree algorithms to identify extragalactic and Galactic transients. Typically we detect 10-15 supernova candidates per night which we immediately announce publicly. The ATLAS discoveries not only enable rapid follow-up of interesting sources but will provide complete statistical samples within the local volume of 100 Mpc. A simple comparison of the detected supernova rate within 100 Mpc, with no corrections for completeness, is already significantly higher (factor 1.5 to 2) than the current accepted rates.

Evolved stars near the tip of the red giant branch (TRGB) show solar-like oscillations with periods spanning hours to months and amplitudes ranging from $\sim$1 mmag to $\sim$100 mmag. The systematic detection of the resulting photometric variations with ground-based telescopes would enable the application of asteroseismology to a much larger and more distant sample of stars than is currently accessible with space-based telescopes such as \textit{Kepler} or the ongoing Transiting Exoplanet Survey Satellite (\textit{TESS}) mission. We present an asteroseismic analysis of 493 M giants using data from two ground-based surveys: the Asteroid Terrestrial-impact Last Alert System (ATLAS) and the All-Sky Automated Survey for Supernovae (ASAS-SN). By comparing the extracted frequencies with constraints from \textit{Kepler}, the Sloan Digital Sky Survey Apache Point Observatory Galaxy Evolution Experiment (APOGEE), and Gaia we demonstrate that ground-based transient surveys allow accurate distance measurements to oscillating M giants with a precision of $\sim$15$\%$. Using stellar population synthesis models we predict that ATLAS and ASAS-SN can provide asteroseismic distances to $\sim$2$\times$10$^{6}$ galactic M giants out to typical distances of $20-50 \; \rm{kpc}$, vastly improving the reach of Gaia and providing critical constraints for Galactic archaeology and galactic dynamics.

We present here the c-o colors for identified Flora, Vesta, Nysa-Polana, Themis, and Koronis family members within the historic data set (2015-2018) of the Asteroid Terrestrial-impact Last Alert System (ATLAS). The Themis and Koronis families are known to be relatively pure C- and S-type Bus-DeMeo taxonomic families, respectively, and the extracted color data from the ATLAS broadband c- and o-filters of these two families is used to demonstrate that the ATLAS c-o color is a sufficient parameter to distinguish between the C- and S-type taxonomies. The Vesta and Nysa-Polana families are known to display a mixture of taxonomies possibly due to Vesta's differentiated parent body origin and Nysa-Polana actually consisting of two nested families with differing taxonomies. Our data show that the Flora family also displays a large degree of taxonomic mixing and the data reveal a substantial H-magnitude dependence on color. We propose and exclude several interpretations for the observed taxonomic mix. Additionally, we extract rotation periods of all of the targets reported here and find good agreement with targets that have previously reported periods.

Digital tracking enables telescopes to detect asteroids several times fainter than conventional techniques. We describe our optimized methodology to acquire, process, and interpret digital tracking observations, and we apply it to probe the apparent magnitude distribution of main belt asteroids fainter than any previously detected from the ground. All-night integrations with the Dark Energy Camera (DECam) yield 95\% completeness at $R$ magnitude 25.0, and useful sensitivity to $R=25.6$ mag when we use an analytical detection model to correct flux overestimation bias. In a single DECam field observed over two nights, we detect a total of 3234 distinct asteroids, of which 3123 are confirmed on both nights. At opposition from the Sun, we find a sky density of $697 \pm 15$ asteroids per square degree brighter than $R = 25.0$ mag, and $1031 \pm 23$ brighter than $R = 25.6$ mag. We agree with published results for the sky density and apparent magnitude distribution of asteroids brighter than $R=23$ mag. For a power law defined by $dN/dR \propto 10^{\alpha R}$, we find marginally acceptable fits with a constant slope $\alpha = 0.28 \pm 0.02$ from $R=20$ to 25.6 mag. Better fits are obtained for a broken power law with $\alpha=0.218 \pm 0.026$ for $R=20$ to 23.5 mag, steepening to $\alpha=0.340 \pm 0.025$ for $R = 23.5$ to 25.6 mag. The constant or steepening power law indicates asteroids fainter than $R = 23.5$ mag are abundant, contrary to some previous claims but consistent with theory.

We present the discovery of PS18kh, a tidal disruption event (TDE) discovered at the center of SDSS J075654.53+341543.6 ($d\simeq322$ Mpc) by the Pan-STARRS Survey for Transients. Our dataset includes pre-discovery survey data from Pan-STARRS, the All-Sky Automated Survey for Supernovae (ASAS-SN), and the Asteroid Terrestrial-impact Last Alert System (ATLAS) as well as high-cadence, multi-wavelength follow-up data from ground-based telescopes and Swift, spanning from 56 days before peak light until 75 days after. The optical/UV emission from PS18kh is well-fit as a blackbody with temperatures ranging from $T\simeq12000$ K to $T\simeq25000$ K and it peaked at a luminosity of $L\simeq8.8\times10^{43}$ ergs s$^{-1}$. PS18kh radiated $E=(3.45\pm0.22)\times10^{50}$ ergs over the period of observation, with $(1.42\pm0.20)\times10^{50}$ ergs being released during the rise to peak. Spectra of PS18kh show a changing, boxy/double-peaked H$\alpha$ emission feature, which becomes more prominent over time. We use models of non-axisymmetric accretion disks to describe the profile of the H$\alpha$ line and its evolution. We find that at early times the high accretion rate leads the disk to emit a wind which modifies the shape of the line profile and makes it bell-shaped. At late times, the wind becomes optically thin, allowing the non-axisymmetric perturbations to show up in the line profile. The line-emitting portion of the disk extends from $r_{\rm in}\sim60r_{\rm g}$ to an outer radius of $r_{\rm out}\sim1400r_{\rm g}$ and the perturbations can be represented either as an eccentricity in the outer rings of the disk or as a spiral arm in the inner disk.

On 2019 January 5 a streamer associated with the 4--10 km main-belt asteroid (6478)~Gault was detected by the ATLAS sky survey, a rare discovery of activity around a main-belt asteroid. Archival data from ATLAS and Pan-STARRS1 show the trail in early December 2018, but not between 2010 and January 2018. The feature has significantly changed over one month, perfectly matching predictions of pure dust dynamical evolution and changes in observing geometry for a short release of dust around 2018 October 28. Follow-up observations with HST show a second narrow trail corresponding to a brief release of dust on 2018 December 30. Both releases occurred with negligible velocity. We find the dust grains to be fairly large, with power-law size distributions in the $10^{-5} - 10^{-3}$~m range and power-law indices of $\sim -1.5$. Three runs of ground-based data find a signature of $\sim 2\,\rm h$ rotation, close to the rotational limit, suggesting that the activity is the result of landslides or reconfigurations after YORP spin-up.

We present observations of supernova (SN) 2017ens, discovered by the ATLAS survey and identified as a hot blue object through the GREAT program. The redshift z=0.1086 implies a peak brightness of M_g=-21.1 mag, placing the object within the regime of superluminous supernovae. We observe a dramatic spectral evolution, from initially being blue and featureless, to later developing features similar to those of the broadlined Type Ic SN 1998bw, and finally showing ~2000 km s^-1 wide H-alpha and H-beta emission. Relatively narrow Balmer emission (reminiscent of a SN IIn) is present at all times. We also detect coronal lines, indicative of a dense circumstellar medium. We constrain the progenitor wind velocity to ~50-60 km s^-1 based on P-Cygni profiles, which is far slower than those present in Wolf-Rayet stars. This may suggest that the progenitor passed through a luminous blue variable phase, or that the wind is instead from a binary companion red supergiant star. At late times we see the ~2000 km s^-1 wide H-alpha emission persisting at high luminosity (~3x10^40 erg s^-1) for at least 100 day, perhaps indicative of additional mass loss at high velocities that could have been ejected by a pulsational pair instability.

Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically-confirmed type Ia supernova (SN Ia) observed in the $Kepler$ field. The $Kepler$ data revealed an excess emission in its early light curve, allowing to place interesting constraints on its progenitor system (Dimitriadis et al. 2018, Shappee et al. 2018b). Here, we present extensive optical, ultraviolet, and near-infrared photometry, as well as dense sampling of optical spectra, for this object. SN 2018oh is relatively normal in its photometric evolution, with a rise time of 18.3$\pm$0.3 days and $\Delta$m$_{15}(B)=0.96\pm$0.03 mag, but it seems to have bluer $B - V$ colors. We construct the "uvoir" bolometric light curve having peak luminosity as 1.49$\times$10$^{43}$erg s$^{-1}$, from which we derive a nickel mass as 0.55$\pm$0.04M$_{\odot}$ by fitting radiation diffusion models powered by centrally located $^{56}$Ni. Note that the moment when nickel-powered luminosity starts to emerge is +3.85 days after the first light in the Kepler data, suggesting other origins of the early-time emission, e.g., mixing of $^{56}$Ni to outer layers of the ejecta or interaction between the ejecta and nearby circumstellar material or a non-degenerate companion star. The spectral evolution of SN 2018oh is similar to that of a normal SN Ia, but is characterized by prominent and persistent carbon absorption features. The C II features can be detected from the early phases to about 3 weeks after the maximum light, representing the latest detection of carbon ever recorded in a SN Ia. This indicates that a considerable amount of unburned carbon exists in the ejecta of SN 2018oh and may mix into deeper layers.

We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unusual two-component shape, where the flux rises with a steep linear gradient for the first few days, followed by a quadratic rise as seen for typical SNe Ia. This "flux excess" relative to canonical SN Ia behavior is confirmed in our $i$-band light curve, and furthermore, SN 2018oh is especially blue during the early epochs. The flux excess peaks 2.14$\pm0.04$ days after explosion, has a FWHM of 3.12$\pm0.04$ days, a blackbody temperature of $T=17,500^{+11,500}_{-9,000}$ K, a peak luminosity of $4.3\pm0.2\times10^{37}\,{\rm erg\,s^{-1}}$, and a total integrated energy of $1.27\pm0.01\times10^{43}\,{\rm erg}$. We compare SN 2018oh to several models that may provide additional heating at early times, including collision with a companion and a shallow concentration of radioactive nickel. While all of these models generally reproduce the early K2 light curve shape, we slightly favor a companion interaction, at a distance of $\sim$$2\times10^{12}\,{\rm cm}$ based on our early color measurements, although the exact distance depends on the uncertain viewing angle. Additional confirmation of a companion interaction in future modeling and observations of SN 2018oh would provide strong support for a single-degenerate progenitor system.

On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System (ATLAS). The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN~Ia light curve, and it unambiguously shows a ~4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a non-degenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing, published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe~Ia light curves. Finally, we use Swift X-ray non-detections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant density CSM these non-detections constrain rho<4.5 * 10^5 cm^-3 at a radius of 4 *10^15 cm from the progenitor star. Assuming a wind-like environment, we place mass-loss limits of Mdot< 8 * 10^-6 M_sun yr^-1 for v_w=100 km s^-1, ruling out some symbiotic progenitor systems.

We present the discovery of ASASSN-18ey (MAXI J1820+070), a new black hole low-mass X-ray binary discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN). A week after ASAS-SN discovered ASASSN-18ey as an optical transient, it was detected as an X-ray transient by MAXI/GCS. Here, we analyze ASAS-SN and Asteroid Terrestrial-impact Last Alert System (ATLAS) pre-outburst optical light curves, finding evidence of intrinsic variability for several years prior to the outburst. While there was no long-term rise leading to outburst, as has been seen in several other systems, the start of the outburst in the optical preceded that in the X-rays by $7.20\pm0.97~\rm days$. We analyze the spectroscopic evolution of ASASSN-18ey from pre-maximum to $> 100~\rm days$ post-maximum. The spectra of ASASSN-18ey exhibit broad, asymmetric, double-peaked H$\alpha$ emission. The Bowen blend ($\lambda \approx 4650$\AA) in the post-maximum spectra shows highly variable double-peaked profiles, likely arising from irradiation of the companion by the accretion disk, typical of low-mass X-ray binaries. The optical and X-ray luminosities of ASASSN-18ey are consistent with black hole low-mass X-ray binaries, both in outburst and quiescence.

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.

The Asteroid Terrestrial-impact Last Alert System (ATLAS) observes most of the sky every night in search of dangerous asteroids. Its data are also used to search for photometric variability, where sensitivity to variability is limited by photometric accuracy. Since each exposure spans 7.6 deg corner to corner, variations in atmospheric transparency in excess of 0.01 mag are common, and 0.01 mag photometry cannot be achieved by using a constant flat field calibration image. We therefore have assembled an all-sky reference catalog of approximately one billion stars to m~19 from a variety of sources to calibrate each exposure's astrometry and photometry. Gaia DR2 is the source of astrometry for this ATLAS Refcat2. The sources of g, r, i, z photometry include Pan-STARRS DR1, the ATLAS Pathfinder photometry project, ATLAS re-flattened APASS data, SkyMapper DR1, APASS DR9, the Tycho-2 catalog, and the Yale Bright Star Catalog. We have attempted to make this catalog at least 99% complete to m<19, including the brightest stars in the sky. We believe that the systematic errors are no larger than 5 millimag RMS, although errors are as large as 20 millimag in small patches near the galactic plane.

Super-luminous supernovae (SLSNe) are rare events defined as being significantly more luminous than normal terminal stellar explosions. The source of the extra powering needed to achieve such luminosities is still unclear. Discoveries in the local Universe (i.e. $z<0.1$) are scarce, but afford dense multi-wavelength observations. Additional low-redshift objects are therefore extremely valuable. We present early-time observations of the type I SLSN ASASSN-18km/SN~2018bsz. These data are used to characterise the event and compare to literature SLSNe and spectral models. Host galaxy properties are also analysed. Optical and near-IR photometry and spectroscopy were analysed. Early-time ATLAS photometry was used to constrain the rising light curve. We identified a number of spectral features in optical-wavelength spectra and tracked their time evolution. Finally, we used archival host galaxy photometry together with HII region spectra to constrain the host environment. ASASSN-18km/SN~2018bsz is found to be a type I SLSN in a galaxy at a redshift of 0.0267 (111 Mpc), making it the lowest-redshift event discovered to date. Strong CII lines are identified in the spectra. Spectral models produced by exploding a Wolf-Rayet progenitor and injecting a magnetar power source are shown to be qualitatively similar to ASASSN-18km/SN~2018bsz, contrary to most SLSNe-I that display weak/non-existent CII lines. ASASSN-18km/SN~2018bsz displays a long, slowly rising, red 'plateau' of $>$26 days, before a steeper, faster rise to maximum. The host has an absolute magnitude of --19.8 mag ($r$), a mass of M$_{*}$ = 1.5$^{+0.08}_{-0.33}$ $\times$10$^{9}$ M$_{\odot}$ , and a star formation rate of = 0.50$^{+2.22}_{-0.19}$ M$_{\odot}$ yr$^{-1}$. A nearby HII region has an oxygen abundance (O3N2) of 8.31$\pm$0.01 dex.

We present the ATLAS discovery and initial analysis of the first 18 days of the unusual transient event, ATLAS18qqn/AT2018cow. It is characterized by a high peak luminosity ($\sim$1.7 $\times$ 10$^{44}$ erg s$^{-1}$), rapidly evolving light curves ($>$5 mag rise to peak in $\sim$3.5 days), and hot blackbody spectra, peaking at $\sim$27000 K that are relatively featureless and unchanging over the first two weeks. The bolometric light curve cannot be powered by radioactive decay under realistic assumptions. The detection of high-energy emission may suggest a central engine as the powering source. Using a magnetar model, we estimated an ejected mass of $0.1-0.4$ \msol, which lies between that of low-energy core-collapse events and the kilonova, AT2017gfo. The spectra cooled rapidly from 27000 to 15000 K in just over 2 weeks but remained smooth and featureless. Broad and shallow emission lines appear after about 20 days, and we tentatively identify them as He I although they would be redshifted from their rest wavelengths. We rule out that there are any features in the spectra due to intermediate mass elements up to and including the Fe-group. The presence of r-process elements cannot be ruled out. If these lines are due to He, then we suggest a low-mass star with residual He as a potential progenitor. Alternatively, models of magnetars formed in neutron-star mergers give plausible matches to the data.

Technology has advanced to the point that it is possible to image the entire sky every night and process the data in real time. The sky is hardly static: many interesting phenomena occur, including variable stationary objects such as stars or QSOs, transient stationary objects such as supernovae or M dwarf flares, and moving objects such as asteroids and the stars themselves. Funded by NASA, we have designed and built a sky survey system for the purpose of finding dangerous near-Earth asteroids (NEAs). This system, the "Asteroid Terrestrial-impact Last Alert System" (ATLAS), has been optimized to produce the best survey capability per unit cost, and therefore is an efficient and competitive system for finding potentially hazardous asteroids (PHAs) but also for tracking variables and finding transients. While carrying out its NASA mission, ATLAS now discovers more bright ($m < 19$) supernovae candidates than any ground based survey, frequently detecting very young explosions due to its 2 day cadence. ATLAS discovered the afterglow of a gamma-ray burst independent of the high energy trigger and has released a variable star catalogue of 5$\times10^{6}$ sources. This, the first of a series of articles describing ATLAS, is devoted to the design and performance of the ATLAS system. Subsequent articles will describe in more detail the software, the survey strategy, ATLAS-derived NEA population statistics, transient detections, and the first data release of variable stars and transient lightcurves.

SN 2017dio shows both spectral characteristics of a type-Ic supernova (SN) and signs of a hydrogen-rich circumstellar medium (CSM). Prominent, narrow emission lines of H and He are superposed on the continuum. Subsequent evolution revealed that the SN ejecta are interacting with the CSM. The initial SN Ic identification was confirmed by removing the CSM interaction component from the spectrum and comparing with known SNe Ic, and reversely, adding a CSM interaction component to the spectra of known SNe Ic and comparing them to SN 2017dio. Excellent agreement was obtained with both procedures, reinforcing the SN Ic classification. The light curve constrains the pre-interaction SN Ic peak absolute magnitude to be around $M_g = -17.6$ mag. No evidence of significant extinction is found, ruling out a brighter luminosity required by a SN Ia classification. These pieces of evidence support the view that SN 2017dio is a SN Ic, and therefore the first firm case of a SN Ic with signatures of hydrogen-rich CSM in the early spectrum. The CSM is unlikely to have been shaped by steady-state stellar winds. The mass loss of the progenitor star must have been intense, $\dot{M} \sim 0.02$ $(\epsilon_{H\alpha}/0.01)^{-1}$ $(v_\textrm{wind}/500$ km s$^{-1}$) $(v_\textrm{shock}/10 000$ km s$^{-1})^{-3}$ $M_\odot$~yr$^{-1}$, peaking at a few decades before the SN. Such a high mass loss rate might have been experienced by the progenitor through eruptions or binary stripping.

We report the discovery and multi-wavelength data analysis of the peculiar optical transient, ATLAS17aeu. This transient was identified in the skymap of the LIGO gravitational wave event GW170104 by our ATLAS and Pan-STARRS coverage. ATLAS17aeu was discovered 23.1hrs after GW170104 and rapidly faded over the next 3 nights, with a spectrum revealing a blue featureless continuum. The transient was also detected as a fading x-ray source by Swift and in the radio at 6 and 15 GHz. A gamma ray burst GRB170105A was detected by 3 satellites 19.04hrs after GW170104 and 4.10hrs before our first optical detection. We analyse the multi-wavelength fluxes in the context of the known GRB population and discuss the observed sky rates of GRBs and their afterglows. We find it statistically likely that ATLAS17aeu is an afterglow associated with GRB170105A, with a chance coincidence ruled out at the 99\% confidence or 2.6$\sigma$. A long, soft GRB within a redshift range of $1 \lesssim z \lesssim 2.9$ would be consistent with all the observed multi-wavelength data. The Poisson probability of a chance occurrence of GW170104 and ATLAS17aeu is $p=0.04$. This is the probability of a chance coincidence in 2D sky location and in time. These observations indicate that ATLAS17aeu is plausibly a normal GRB afterglow at significantly higher redshift than the distance constraint for GW170104 and therefore a chance coincidence. However if a redshift of the faint host were to place it within the GW170104 distance range, then physical association with GW170104 should be considered.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

We present results from a two-night R~4000 0.9-2.5 micron spectroscopic monitoring campaign of Luhman 16AB (L7.5 + T0.5). We assess the variability amplitude as a function of pressure level in the atmosphere of Luhman 16B: the more variable of the two components. The amplitude decreases monotonically with decreasing pressure, indicating that the source of variability - most likely patchy clouds - lies in the lower atmosphere. An unexpected result is that the strength of the K I absorption is higher in the faint state of Luhman 16B and lower in the bright state. We conclude that either the abundance of K I increases when the clouds roll in, potentially because of additional K I in the cloud itself, or that the temperature-pressure profile changes. We reproduce the change in K I absorption strengths with combinations of spectral templates to represent the bright and the faint variability states. These are dominated by a warmer L8 or L9 component, with a smaller contribution from a cooler T1 or T2 component. The success of this approach argues that the mechanism responsible for brown dwarf variability is also behind the diverse spectral morphology across the L-to-T transition. We further suggest that the L9-T1 part of the sequence represents a narrow but random ordering of effective temperatures and cloud fractions, obscured by the monotonic progression in methane absorption strength.

Supernova (SN) 2016bdu is an unusual transient resembling SN 2009ip. SN 2009ip-like events are characterized by a long-lasting phase of erratic variability which ends with two luminous outbursts a few weeks apart. The second outburst is significantly more luminous (about 3 mag) than the first. In the case of SN 2016bdu, the first outburst (Event A) reached an absolute magnitude M(r) ~ -15.3 mag, while the second one (Event B) occurred over one month later and reached M(r) ~ -18 mag. By inspecting archival data, a faint source at the position of SN 2016bdu is detectable several times in the past few years. We interpret these detections as signatures of a phase of erratic variability, similar to that experienced by SN 2009ip between 2008 and mid-2012, and resembling the currently observed variability of the luminous blue variable SN 2000ch in NGC 3432. Spectroscopic monitoring of SN 2016bdu during the second peak initially shows features typical of a SN IIn. One month after the Event B maximum, the spectra develop broad Balmer lines with P Cygni profiles and broad metal features. At these late phases, the spectra resemble those of a typical Type II SN. All members of this SN 2009ip-like group are remarkably similar to the Type IIn SN 2005gl. For this object, the claim of a terminal SN explosion is supported by the disappearance of the progenitor star. The similarity with SN 2005gl suggests that all members of this family may finally explode as genuine SNe, although the unequivocal detection of nucleosynthesised elements in their nebular spectra is still missing.

Recent photometric studies have revealed that surface spots that produce flux variations are present on virtually all L and T dwarfs. Their likely magnetic or dusty nature has been a much-debated problem, the resolution to which has been hindered by paucity of diagnostic multi-wavelength observations. To test for a correlation between magnetic activity and photometric variability, we searched for H$\alpha$ emission among eight L3$-$T2 ultra-cool dwarfs with extensive previous photometric monitoring, some of which are known to be variable at 3.6 $\mu$m or 4.5 $\mu$m. We detected H$\alpha$ only in the non-variable T2 dwarf 2MASS J12545393$-$0122474. The remaining seven objects do not show H$\alpha$ emission, even though six of them are known to vary photometrically. Combining our results with those for 86 other L and T dwarfs from the literature show that the detection rate of H$\alpha$ emission is very high (94$\%$) for spectral types between L0 and L3.5 and much smaller (20$\%$) for spectral types $\ge$L4, while the detection rate of photometric variability is approximately constant (30$\%-$55$\%$) from L0 to T8 dwarfs. We conclude that chromospheric activity, as evidenced by H$\alpha$ emission, and large-amplitude photometric variability are not correlated. Consequently, dust clouds are the dominant driver of the observed variability of ultra-cool dwarfs at spectral types at least as early as L0.

SN 2016gkg is a nearby Type IIb supernova discovered shortly after explosion. Like several other Type IIb events with early-time data, SN 2016gkg displays a double-peaked light curve, with the first peak associated with the cooling of a low-mass extended progenitor envelope. We present unprecedented intranight-cadence multi-band photometric coverage of the first light-curve peak of SN 2016gkg obtained from the Las Cumbres Observatory Global Telescope network, the Asteroid Terrestrial-impact Last Alert System, the Swift satellite and various amateur-operated telescopes. Fitting these data to analytical shock-cooling models gives a progenitor radius of ~25-140 solar radii with ~2-30 x 10^-2 solar masses of material in the extended envelope (depending on the model and the assumed host-galaxy extinction). Our radius estimates are broadly consistent with values derived independently (in other works) from HST imaging of the progenitor star. However, the shock-cooling model radii are on the lower end of the values indicated by pre-explosion imaging. Hydrodynamical simulations could refine the progenitor parameters deduced from the shock-cooling emission and test the analytical models.

In an earlier study we reported nearly 100 previously unknown dusty debris disks around Hipparcos main sequence stars within 75 pc by selecting stars with excesses in individual WISE colors. Here, we further scrutinize the Hipparcos 75 pc sample to (1) gain sensitivity to previously undetected, fainter mid-IR excesses and (2) to remove spurious excesses contaminated by previously unidentified blended sources. We improve upon our previous method by adopting a more accurate measure of the confidence threshold for excess detection, and by adding an optimally-weighted color average that incorporates all shorter-wavelength WISE photometry, rather than using only individual WISE colors. The latter is equivalent to spectral energy distribution fitting, but only over WISE band passes. In addition, we leverage the higher resolution WISE images available through the unWISE.me image service to identify contaminated WISE excesses based on photocenter offsets among the W3- and W4-band images. Altogether, we identify 19 previously unreported candidate debris disks. Combined with the results from our earlier study, we have found a total of 107 new debris disks around 75 pc Hipparcos main sequence stars using precisely calibrated WISE photometry. This expands the 75 pc debris disk sample by 22% around Hipparcos main-sequence stars and by 20% overall (including non-main sequence and non-Hipparcos stars).

We present a search for an electromagnetic counterpart of the gravitational wave source GW151226. Using the Pan-STARRS1 telescope we mapped out 290 square degrees in the optical i_ps filter starting 11.5hr after the LIGO information release and lasting for a further 28 days. The first observations started 49.5hr after the time of the GW151226 detection. We typically reached sensitivity limits of i_ps = 20.3-20.8 and covered 26.5% of the LIGO probability skymap. We supplemented this with ATLAS survey data, reaching 31% of the probability region to shallower depths of m~19. We found 49 extragalactic transients (that are not obviously AGN), including a faint transient in a galaxy at 7Mpc (a luminous blue variable outburst) plus a rapidly decaying M-dwarf flare. Spectral classification of 20 other transient events showed them all to be supernovae. We found an unusual transient, PS15dpn, with an explosion date temporally coincident with GW151226 which evolved into a type Ibn supernova. The redshift of the transient is secure at z=0.1747 +/- 0.0001 and we find it unlikely to be linked, since the luminosity distance has a negligible probability of being consistent with that of GW151226. In the 290 square degrees surveyed we therefore do not find a likely counterpart. However we show that our survey strategy would be sensitive to NS-NS mergers producing kilonovae at D < 100 Mpc which is promising for future LIGO/Virgo searches.

We searched for an optical counterpart to the first gravitational wave source discovered by LIGO (GW150914), using a combination of the Pan-STARRS1 wide-field telescope and the PESSTO spectroscopic follow-up programme. As the final LIGO sky maps changed during analysis, the total probability of the source being spatially coincident with our fields was finally only 4.2 per cent. Therefore we discuss our results primarily as a demonstration of the survey capability of Pan-STARRS and spectroscopic capability of PESSTO. We mapped out 442 square degrees of the northern sky region of the initial map. We discovered 56 astrophysical transients over a period of 41 days from the discovery of the source. Of these, 19 were spectroscopically classified and a further 13 have host galaxy redshifts. All transients appear to be fairly normal supernovae and AGN variability and none is obviously linked with GW150914. We illustrate the sensitivity of our survey by defining parameterised lightcurves with timescales of 4, 20 and 40 days and use the sensitivity of the Pan-STARRS1 images to set limits on the luminosities of possible sources. The Pan-STARRS1 images reach limiting magnitudes of i = 19.2, 20.0 and 20.8 respectively for the three timescales. For long timescale parameterised lightcurves (with FWHM=~40d) we set upper limits of M_i <= -17.2 -0.9/+1.4 if the distance to GW150914 is D = 400 +/- 200Mpc. The number of type Ia SN we find in the survey is similar to that expected from the cosmic SN rate, indicating a reasonably complete efficiency in recovering supernova like transients out to D = 400 +/- 200 Mpc.

Using data from the WISE All-Sky Survey, we have found >100 new infrared excess sources around main-sequence Hipparcos stars within 75pc. Our empirical calibration of WISE photospheric colors and removal of non-trivial false-positive sources are responsible for the high confidence (>99.5%) of detections, while our corrections to saturated W1 (3.4um) and W2 (4.6um) photometry have for the first time allowed us to search for new infrared excess sources around bright field stars in WISE. The careful calibration and filtering of the WISE data have allowed us to probe excess fluxes down to roughly 8% of the photospheric emission at 22um around saturated stars in WISE. We expect that the increased sensitivity of our survey will not only aid in understanding the evolution of debris disks, but will also benefit future studies using WISE.

We present a method for calculating precise distances to asteroids using only two nights of data from a single location --- far too little for an orbit --- by exploiting the angular reflex motion of the asteroids due to Earth's axial rotation. We refer to this as the rotational reflex velocity method. While the concept is simple and well-known, it has not been previously exploited for surveys of main-belt asteroids. We offer a mathematical development, estimates of the errors of the approximation, and a demonstration using a sample of 197 asteroids observed for two nights with a small, 0.9-meter telescope. This demonstration used digital tracking to enhance detection sensitivity for faint asteroids, but our distance determination works with any detection method. Forty-eight asteroids in our sample had known orbits prior to our observations, and for these we demonstrate a mean fractional error of only 1.6% between the distances we calculate and those given in ephemerides from the Minor Planet Center. In contrast to our two-night results, distance determination by fitting approximate orbits requires observations spanning 7--10 nights. Once an asteroid's distance is known, its absolute magnitude and size (given a statistically-estimated albedo) may immediately be calculated. Our method will therefore greatly enhance the efficiency with which 4-meter and larger telescopes can probe the size distribution of small (e.g. 100 meter) main belt asteroids. This distribution remains poorly known, yet encodes information about the collisional evolution of the asteroid belt --- and hence the history of the Solar System.

We describe digital tracking, a method for asteroid searches that greatly increases the sensitivity of a telescope to faint unknown asteroids. It has been previously used to detect faint Kuiper Belt objects using the Hubble Space Telescope and large ground-based instruments, and to find a small, fast-moving asteroid during a close approach to Earth. We complement this earlier work by developing digital tracking methodology for detecting asteroids using large-format CCD imagers. We demonstrate that the technique enables the ground-based detection of large numbers of new faint asteroids. Our methodology resolves or circumvents all major obstacles to the large-scale application of digital tracking for finding main belt and near-Earth asteroids. We find that for both asteroid populations, digital tracking can deliver a factor of ten improvement over conventional searches. Digital tracking has long been standard practice for deep Kuiper Belt surveys, but even there our methodology enables deeper integrations than have yet been attempted. Our search for main belt asteroids using a one-degree imager on the 0.9m WIYN telescope on Kitt Peak validates our methodology, delivers sensitivity to asteroids in a regime previously probed only with 4-meter and larger instruments, and leads to the detection of 156 previously unknown asteroids and 59 known objects in a single field. Digital tracking has the potential to revolutionize searches for faint moving objects ranging from the Kuiper Belt through main belt and near-Earth asteroids, and perhaps even anthropogenic space debris in low Earth orbit.

We have monitored twelve T dwarfs with the Kitt Peak 2.1m telescope using an f814w filter (0.7-0.95 microns) to place in context the remarkable 10-20% variability exhibited by the nearby T dwarf Luhman 16B in this wavelength regime. The motivation was the poorly known red optical behavior of T dwarfs, which have been monitored almost exclusively at infrared wavelengths, where variability amplitudes greater than 10% have been found to be very rare. We detect highly significant variability in two T dwarfs. The T2.5 dwarf 2MASS 13243559+6358284 shows consistent ~17% variability on two consecutive nights. The T2 dwarf 2MASS J16291840+0335371 exhibits ~10% variability that may evolve from night to night, similarly to Luhman 16B. Both objects were previously known to be variable in the infrared, but with considerably lower amplitudes. We also find evidence for variability in the T6 dwarf J162414.37+002915.6, but since it has lower significance, we conservatively refrain from claiming this object as a variable. We explore and rule out various telluric effects, demonstrating that the variations we detect are astrophysically real. We suggest that high-amplitude photometric variability for T dwarfs is likely more common in the red optical than at longer wavelengths. The two new members of the growing class of high-amplitude variable T dwarfs offer excellent prospects for further study of cloud structures and their evolution.

We present time-resolved near-infrared spectroscopy of two L5 dwarfs, 2MASS J18212815+1414010 and 2MASS J15074759-1627386, observed with the Wide Field Camera 3 instrument on the Hubble Space Telescope (HST). We study the wavelength dependence of rotation-modulated flux variations between 1.1 $\mu$m and 1.7 $\mu$m. We find that the water absorption bands of the two L5 dwarfs at 1.15 $\mu$m and 1.4 $\mu$m vary at similar amplitudes as the adjacent continuum. This differs from the results of previous HST observations of L/T transition dwarfs, in which the water absorption at 1.4 $\mu$m displays variations of about half of the amplitude at other wavelengths. We find that the relative amplitude of flux variability out of the water band with respect to that in the water band shows a increasing trend from the L5 dwarfs toward the early T dwarfs. We utilize the models of Saumon & Marley (2008) and find that the observed variability of the L5 dwarfs can be explained by the presence of spatially varying high-altitude haze layers above the condensate clouds. Therefore, our observations show that the heterogeneity of haze layers - the driver of the variability - must be located at very low pressures, where even the water opacity is negligible. In the near future, the rotational spectral mapping technique could be utilized for other atomic and molecular species to probe different pressure levels in the atmospheres of brown dwarfs and exoplanets and uncover both horizontal and vertical cloud structures.

We present results from the "Weather on Other Worlds" Spitzer Exploration Science program to investigate photometric variability in L and T dwarfs, usually attributed to patchy clouds. We surveyed 44 L3-T8 dwarfs, spanning a range of $J-K_s$ colors and surface gravities. We find that 14/23 (61%; 95% confidence interval: 41%-78%) of our single L3-L9.5 dwarfs are variable with peak-to-peak amplitudes between 0.2% and 1.5%, and 5/16 (31%; 95% confidence interval: 14%-56%) of our single T0-T8 dwarfs are variable with amplitudes between 0.8% and 4.6%. After correcting for sensitivity, we find that 80% (95% confidence interval: 53%-100%) of L dwarfs vary by >0.2%, and 36% (95% confidence interval: 19%-52%) of T dwarfs vary by >0.4%. Given viewing geometry considerations, we conclude that photospheric heterogeneities causing >0.2% 3-5-micron flux variations are present on virtually all L dwarfs, and probably on most T dwarfs. A third of L dwarf variables show irregular light curves, indicating that L dwarfs may have multiple spots that evolve over a single rotation. Also, approximately a third of the periodicities are on time scales >10 h, suggesting that slowly-rotating brown dwarfs may be common. We observe an increase in the maximum amplitudes over the entire spectral type range, revealing a potential for greater temperature contrasts in T dwarfs than in L dwarfs. We find a tentative association (92% confidence) between low surface gravity and high-amplitude variability among L3-L5.5 dwarfs. Although we can not confirm whether lower gravity is also correlated with a higher incidence of variables, the result is promising for the characterization of directly imaged young extrasolar planets through variability.

We present a sensitive search for WISE W3 (12um) and W4 (22um) excesses from warm optically thin dust around Hipparcos main sequence stars within 75pc from the Sun. We use contemporaneously measured photometry from WISE, remove sources of contamination, and derive and apply corrections to saturated fluxes to attain optimal sensitivity to >10um excesses. We use data from the WISE All-Sky Survey Catalog rather than the AllWISE release, because we find that its saturated photometry is better behaved, allowing us to detect small excesses even around saturated stars in WISE. Our new discoveries increase by 45% the number of stars with warm dusty excesses and expand the number of known debris disks (with excess at any wavelength) within 75pc by 29%. We identify 220 Hipparcos debris disk-host stars, 108 of which are new detections at any wavelength. We present the first measurement of a 12um and/or 22um excess for 10 stars with previously known cold (50-100 K) disks. We also find five new stars with small but significant W3 excesses, adding to the small population of known exozodi, and we detect evidence for a W2 excess around HIP96562 (F2V), indicative of tenuous hot (780 K) dust. As a result of our WISE study, the number of debris disks with known 10-30um excesses within 75pc (379) has now surpassed the number of disks with known >30um excesses (289, with 171 in common), even if the latter have been found to have a higher occurrence rate in unbiased samples.

Photometric monitoring from warm Spitzer reveals that the L3 dwarf DENIS-P J1058.7-1548 varies sinusoidally in brightness with a period of 4.25 +0.26/-0.16 hours and an amplitude of 0.388 +/- 0.043% (peak-to-valley) in the 3.6 micron band, confirming the reality of a 4.31 +/- 0.31 hour periodicity detected in J-band photometry from the SOAR telescope. The J-band variations are a factor of 2.17 +/- 0.35 larger in amplitude than those at 3.6 microns, while 4.5 micron Spitzer observations yield a 4.5 micron/3.6 micron amplitude ratio of only 0.23 +/- 0.15, consistent with zero 4.5 micron variability. This wide range in amplitudes indicates rotationally modulated variability due to magnetic phenomena and/or inhomogeneous cloud cover. Weak H-alpha emission indicates some magnetic activity, but it is difficult to explain the observed amplitudes by magnetic phenomena unless they are combined with cloud inhomogeneities (which might have a magnetic cause). However, inhomogenous cloudcover alone can explain all our observations, and our data align with theory in requiring that the regions with the thickest clouds also have the lowest effective temperature. Combined with published v sin(i) results, our rotation period yields a 95% confidence lower limit of 0.111 solar radii on the radius of DENIS-P J1058.7-1548, suggesting upper limits of 320 Myr and 0.055 solar masses on the age and mass. These limits should be regarded cautiously because of 3-sigma inconsistencies with other data; however, a lower limit of 45 degrees on the inclination is more secure. DENIS-P J1058.7-1548 is only the first of nearly two dozen low-amplitude variables discovered and analyzed by the Weather on Other Worlds project.

We present the observational results of an L' and M band Adaptive Optics (AO) imaging survey of 54 nearby, sunlike stars for extrasolar planets, carried out using the Clio camera on the MMT. We have concentrated more strongly than all other planet imaging surveys to date on very nearby F, G, and K stars, prioritizing stellar proximity higher than youth. Ours is also the first survey to include extensive observations in the M band, which supplement the primary L' observations. Models predict much better planet/star flux ratios at the L' and M bands than at more commonly used shorter wavelengths (i.e. the H band). We have carried out extensive blind simulations with fake planets inserted into the raw data to verify our sensitivity, and to establish a definitive relationship between source significance in $\sigma$ and survey completeness. We find 97% confident-detection completeness for 10$\sigma$ sources, but only 46% for 7$\sigma$ sources -- raising concerns about the standard procedure of assuming high completeness at 5$\sigma$, and demonstrating that blind sensitivity tests to establish the significance-completeness relation are an important analysis step for all planet-imaging surveys. We discovered a previously unknown, approximately 0.15 solar-mass stellar companion to the F9 star GJ 3876, at a projected separation of about 80 AU. Twelve additional candidate faint companions are detected around other stars. Of these, eleven are confirmed to be background stars, and one is a previously known brown dwarf. We obtained sensitivity to planetary-mass objects around almost all of our target stars, with sensitivity to objects below 3 Jupiter masses in the best cases. Constraints on planet populations based on this null result are presented in our Modeling Results paper.

We have carried out an L' and M band Adaptive Optics (AO) extrasolar planet imaging survey of 54 nearby, sunlike stars using the Clio camera at the MMT. Our survey concentrates more strongly than all others to date on very nearby F, G, and K stars, in that we have prioritized proximity higher than youth. Our survey is also the first to include extensive observations in the M band, which supplemented the primary L' observations. These longer wavelength bands are most useful for very nearby systems in which low temperature planets with red IR colors (i.e. H - L', H - M) could be detected. The survey detected no planets, but set interesting limits on planets and brown dwarfs in the star systems we investigated. We have interpreted our null result by means of extensive Monte Carlo simulations, and constrained the distributions of extrasolar planets in mass $M$ and semimajor axis $a$. If planets are distributed according to a power law with $dN \propto M^{\alpha} a^{\beta} dM da$, normalized to be consistent with radial velocity statistics, we find that a distribution with $\alpha = -1.1$ and $\beta = -0.46$, truncated at 110 AU, is ruled out at the 90% confidence level. These particular values of $\alpha$ and $\beta$ are significant because they represent the most planet-rich case consistent with current statistics from radial velocity observations. With 90% confidence no more than 8.1% of stars like those in our survey have systems with three widely spaced, massive planets like the A-star HR 8799. Our observations show that giant planets in long-period orbits around sun-like stars are rare, confirming the results of shorter-wavelength surveys, and increasing the robustness of the conclusion.

We present direct imaging observations at wavelengths of 3.3, 3.8 (L',band), and 4.8 (M band) microns, for the planetary system surrounding HR 8799. All three planets are detected at L'. The c and d component are detected at 3.3 microns, and upper limits are derived from the M band observations. These observations provide useful constraints on warm giant planet atmospheres. We discuss the current age constraints on the HR 8799 system, and show that several potential co-eval objects can be excluded from being co-moving with the star. Comparison of the photometry is made to models for giant planet atmospheres. Models which include non-equilibrium chemistry provide a reasonable match to the colors of c and d. From the observed colors in the thermal infrared we estimate T_eff < 960 K for b, and T_eff=1300 and 1170 K for c and d, respectively. This provides an independent check on the effective temperatures and thus masses of the objects from the Marois 2008 results.

A candidate < 3 Jupiter mass, extrasolar planet was recently imaged by Kalas et al. (2008) using HST/ACS at 12.7" (96 AU) separation from the nearby (d = 7.7 pc) young (~200 Myr) A2V star Fomalhaut. Here we report results from M-band (4.8 micron) imaging of Fomalhaut on 5 Dec 2006 using the Clio IR imager on the 6.5-m MMT with the adaptive secondary mirror. Our images are sensitive to giant planets at orbital radii comparable to the outer solar system (~5-40 AU). Comparing our 5-sigma M-band photometric limits to theoretical evolutionary tracks for substellar objects, our results rule out the existence of planets with masses greater than 2 Jupiter masses, from ~13-40 AU and objects greater than 13 Jupiter masses from ~8-40 AU.

We have obtained deep Adaptive Optics (AO) images of Vega and epsilon Eri to search for planetary-mass companions. We observed at the MMT in the L' (3.8 micron) and M (4.8 micron) bands using Clio, a recently commissioned imager optimized for these wavelengths. Observing at these long wavelengths represents a departure from the H band (1.65 microns) more commonly used for AO imaging searches for extrasolar planets. The long wavelengths offer better predicted planet/star flux ratios and cleaner (higher Strehl) AO images, at the cost of lower diffraction limited resolution and higher sky background. We have not detected any planets or planet candidates around Vega or epsilon Eri. We report the sensitivities obtained around both stars, which correspond to upper limits on any planetary companions which may exist. The sensitivities of our L' and M band observations are comparable to those of the best H-regime observations of these stars. For epsilon Eri our M band observations deliver considerably better sensitivity to close-in planets than any previously published results, and we show that the M band is by far the best wavelength choice for attempts at ground-based AO imaging of the known planet epsilon Eri b. The Clio camera itself with MMTAO may be capable of detecting epsilon Eri b at its 2010 apastron, given a multi-night observing campaign. Clio appears to be the only currently existing AO imager that has a realistic possibility of detecting epsilon Eri b.

We present the first astronomical observations obtained with an Apodizing Phase Plate (APP). The plate is designed to suppress the stellar diffraction pattern by 5 magnitudes from 2-9 lambda/D over a 180 degree region. Stellar images were obtained in the M' band (4.85 microns) at the MMTO 6.5m telescope, with adaptive wavefront correction made with a deformable secondary mirror designed for low thermal background observations. The measured PSF shows a halo intensity of 0.1% of the stellar peak at 2 lambda/D (0.36 arcsec), tapering off as r^{-5/3} out to radius 9 lambda/D. Such a profile is consistent with residual errors predicted for servo lag in the AO system. We project a 5 sigma contrast limit, set by residual atmospheric fluctuations, of 10.2 magnitudes at 0.36 arcsec separation for a one hour exposure. This can be realised if static and quasi-static aberrations are removed by differential imaging, and is close to the sensitivity level set by thermal background photon noise for target stars with M'>3. The advantage of using the phase plate is the removal of speckle noise caused by the residuals in the diffraction pattern that remain after PSF subtraction. The APP gives higher sensitivity over the range 2-5 lambda/D compared to direct imaging techniques.

We are carrying out a survey to search for giant extrasolar planets around nearby, moderate-age stars in the mid-infrared L' and M bands (3.8 and 4.8 microns, respectively), using the Clio camera with the adaptive optics system on the MMT telescope. To date we have observed 7 stars, of a total 50 planned, including GJ 450 (distance about 8.55pc, age about 1 billion years, no real companions detected), which we use as our example here. We report the methods we use to obtain extremely high contrast imaging in L', and the performance we have obtained. We find that the rotation of a celestial object over time with respect to a telescope tracking it with an altazimuth mount can be a powerful tool for subtracting telescope-related stellar halo artifacts and detecting planets near bright stars. We have carried out a thorough Monte Carlo simulation demonstrating our ability to detect planets as small as 6 Jupiter masses around GJ 450. The division of a science data set into two independent parts, with companions required to be detected on both in order to be recognized as real, played a crucial role in detecting companions in this simulation. We mention also our discovery of a previously unknown faint stellar companion to another of our survey targets, HD 133002. Followup is needed to confirm this as a physical companion, and to determine its physical properties.

Clio is an adaptive-optics camera mounted on the 6.5 meter MMT optimized for diffraction-limited L' and M-band imaging over a ~15'' field. The instrument was designed from the ground up with a large well-depth, fast readout thermal infrared (~3-5 micron) 320 by 256 pixel InSb detector, cooled optics, and associated focal plane and pupil masks (with the option for a coronograph) to minimize the thermal background and maximize throughput. When coupled with the MMT's adaptive secondary AO (two warm reflections) system's low thermal background, this instrument is in a unique position to image nearby warm planets, which are the brightest in the L' and M-band atmospheric windows. We present the current status of this recently commissioned instrument that performed exceptionally during first light. Our instrument sensitivities are impressive and are sky background limited: for an hour of integration, we obtain an L'-band 5 sigma detection limit of of 17.0 magnitudes (Strehl ~80%) and an M-band limit of 14.5 (Strehl ~90%). Our M-band sensitivity is lower due to the increase in thermal sky background. These sensitivities translate to finding relatively young planets five times Jupiter mass at 10 pc within a few AU of a star. Presently, a large Clio survey of nearby stellar systems is underway including a search for planets around solar-type stars, M dwarfs, and white dwarfs. Even with a null result, we can place strong constraints on planet distribution models.

Vega may have a massive companion in a wide orbit, as evidenced by structure in its cold dust debris. We have tested this hypothesis by direct imaging with adaptive optics in the M band. The observations were made with a newly commissioned thermal infrared camera, Clio, on the 6.5 MMT AO system with low-background deformable secondary. The observations constrain a planet to be less than 7 M$_J$ at the approximate position angle expected from the dust structure and at a radius $>$ 20AU (2.5 arcsec) . This result is more stringent than similar previous near-infrared observations of Vega, that achieve limits of 20 and 10 M$_J$ at separations of 7 arcsec. The higher sensitivity is due both to the more favorable contrast of gas giant planets at M band and to the higher Strehl and more stable point spread function at longer wavelengths. Future L' or M band observations could provide a powerful approach for wide separation planet detection, especially for cooler, and thus older or less massive planets. The natural best targets are nearby stars where planets in the range of 5-15 M$_J$ and as old as several Gyr are expected to be detectable with this technique.