The Arizona Robotic Telescope Network (ARTN) project is a long term effort to develop a system of telescopes to carry out a flexible program of PI observing, survey projects, and time domain astrophysics including monitoring, rapid response, and transient/target-of-opportunity followup. Steward Observatory operates and shares in several 1-3m class telescopes with quality sites and instrumentation, largely operated in classical modes. Science programs suited to these telescopes are limited by scheduling flexibility and people-power of available observers. Our goal is to adapt these facilities for multiple co-existing queued programs, interrupt capability, remote/robotic operation, and delivery of reduced data. In the long term, planning for the LSST era, we envision an automated system coordinating across multiple telescopes and sites, where alerts can trigger followup, classification, and triggering of further observations if required, such as followup imaging that can trigger spectroscopy. We are updating telescope control systems and software to implement this system in stages, beginning with the Kuiper 61'' and Vatican Observatory 1.8-m telescopes. The Kuiper 61'' and its Mont4K camera can now be controlled and queue-scheduled by the RTS2 observatory control software, and operated from a remote room at Steward. We discuss science and technical requirements for ARTN, and some of the challenges in adapting heterogenous legacy facilities, scheduling, data pipelines, and maintaining capabilities for a diverse user base.

We conduct a statistical analysis of a combined sample of direct imaging data, totalling nearly 250 stars. The stars cover a wide range of ages and spectral types, and include five detections ($\kappa$ And b, two $\sim$60 M$_{\rm J}$ brown dwarf companions in the Pleiades, PZ Tel B, and CD$-$35 2722B). For some analyses we add a currently unpublished set of SEEDS observations, including the detections GJ 504b and GJ 758B. We conduct a uniform, Bayesian analysis of all stellar ages using both membership in a kinematic moving group and activity/rotation age indicators. We then present a new statistical method for computing the likelihood of a substellar distribution function. By performing most of the integrals analytically, we achieve an enormous speedup over brute-force Monte Carlo. We use this method to place upper limits on the maximum semimajor axis of the distribution function derived from radial-velocity planets, finding model-dependent values of $\sim$30--100 AU. Finally, we model the entire substellar sample, from massive brown dwarfs to a theoretically motivated cutoff at $\sim$5 M$_{\rm Jup}$, with a single power law distribution. We find that $p(M, a) \propto M^{-0.65\pm0.60} a^{-0.85\pm0.39}$ (1$\sigma$ errors) provides an adequate fit to our data, with 1.0--3.1\% (68\% confidence) of stars hosting 5--70 $M_{\rm Jup}$ companions between 10 and 100 AU. This suggests that many of the directly imaged exoplanets known, including most (if not all) of the low-mass companions in our sample, formed by fragmentation in a cloud or disk, and represent the low-mass tail of the brown dwarfs.

We present results from the first three years of observations of moving group targets in the SEEDS high-contrast imaging survey of exoplanets and disks using the Subaru telescope. We achieve typical contrasts of ~10^5 at 1" and ~10^6 beyond 2" around 63 proposed members of nearby kinematic moving groups. We review each of the kinematic associations to which our targets belong, concluding that five, \beta Pictoris (~20 Myr), AB Doradus (~100 Myr), Columba (~30 Myr), Tucana-Horogium (~30 Myr), TW Hydrae (~10 Myr), are sufficiently well-defined to constrain the ages of individual targets. Somewhat less than half of our targets are high-probability members of one of these moving groups. For all of our targets, we combine proposed moving group membership with other age indicators where available, including Ca II HK emission, X-ray activity, and rotation period, to produce a posterior probability distribution of age. SEEDS observations discovered a substellar companion to one of our targets, \kappa And, a late B star. We do not detect any other substellar companions, but do find seven new close binary systems, of which one still needs to be confirmed. A detailed analysis of the statistics of this sample, and of the companion mass constraints given our age probability distributions and exoplanet cooling models, will be presented in a forthcoming paper.

We describe Algorithms for Calibration, Optimized Registration, and Nulling the Star in Angular Differential Imaging (ACORNS-ADI), a new, parallelized software package to reduce high-contrast imaging data, and its application to data from the SEEDS survey. We implement several new algorithms, including a method to register saturated images, a trimmed mean for combining an image sequence that reduces noise by up to ~20%, and a robust and computationally fast method to compute the sensitivity of a high-contrast observation everywhere on the field-of-view without introducing artificial sources. We also include a description of image processing steps to remove electronic artifacts specific to Hawaii2-RG detectors like the one used for SEEDS, and a detailed analysis of the Locally Optimized Combination of Images (LOCI) algorithm commonly used to reduce high-contrast imaging data. ACORNS-ADI is written in python. It is efficient and open-source, and includes several optional features which may improve performance on data from other instruments. ACORNS-ADI requires minimal modification to reduce data from instruments other than HiCIAO. It is freely available for download at www.github.com/t-brandt/acorns-adi under a BSD license.

We present Subaru/HiCIAO H-band high-contrast images of the debris disk around HIP 79977, whose pres- ence was recently inferred from an infrared excess. Our images resolve the disk for the first time, allowing characterization of its shape, size, and dust grain properties. We use angular differential imaging (ADI) to reveal the disk geometry in unpolarized light out to a radius of ~2", as well as polarized differential imaging (PDI) to measure the degree of scattering polarization out to ~1.5". In order to strike a favorable balance between suppression of the stellar halo and conservation of disk flux, we explore the application of principal component analysis (PCA) to both ADI and reference star subtraction. This allows accurate forward modeling of the effects of data reduction on simulated disk images, and thus direct comparison with the imaged disk. The resulting best-fit values and well-fitting intervals for the model parameters are a surface brightness power-law slope of S_out = -3.2 [-3.6,-2.9], an inclination of i = 84{\deg} [81{\deg},86{\deg}], a high Henyey-Greenstein forward-scattering parameter of g = 0.45 [0.35, 0.60], and a non-significant disk-star offset of u = 3.0 [-1.5, 7.5] AU = 24 [-13, 61] mas along the line of nodes. Furthermore, the tangential linear polarization along the disk rises from ~10% at 0.5" to ~45% at 1.5". These measurements paint a consistent picture of a disk of dust grains produced by collisional cascades and blown out to larger radii by stellar radiation pressure.

We present high-contrast images of HR 4796 A taken with Subaru/HiCIAO in H-band, resolving the debris disk in scattered light. The application of specialized angular differential imaging methods (ADI) allows us to trace the inner edge of the disk with high precision, and reveals a pair of "streamers" extending radially outwards from the ansae. Using a simple disk model with a power-law surface brightness profile, we demonstrate that the observed streamers can be understood as part of the smoothly tapered outer boundary of the debris disk, which is most visible at the ansae. Our observations are consistent with the expected result of a narrow planetesimal ring being ground up in a collisional cascade, yielding dust with a wide range of grain sizes. Radiation forces leave large grains in the ring and push smaller grains onto elliptical, or even hyperbolic trajectories. We measure and characterize the disk's surface brightness profile, and confirm the previously suspected offset of the disk's center from the star's position along the ring's major axis. Furthermore, we present first evidence for an offset along the minor axis. Such offsets are commonly viewed as signposts for the presence of unseen planets within a disk's cavity. Our images also offer new constraints on the presence of companions down to the planetary mass regime (~9 Jupiter masses at 0.5", ~3 Jupiter masses at 1").

GJ 758 B is a cold (~600K) companion to a Sun-like star at 29 AU projected separation, which was recently detected with high-contrast imaging. Here we present photometry of the companion in seven photometric bands from Subaru/HiCIAO, Gemini/NIRI and Keck/NIRC2, providing a rich sampling of the spectral energy distribution in the 1-5 micron wavelength range. A clear detection at 1.58 micron combined with an upper limit at 1.69 micron shows methane absorption in the atmosphere of the companion. The mass of the companion remains uncertain, but an updated age estimate indicates that the most likely mass range is ~30-40 Mjup. In addition, we present an updated astrometric analysis that imposes tighter constraints on GJ 758 B's orbit and identifies the proposed second candidate companion, "GJ 758 C", as a background star.

We report high-resolution 1.6 $\micron$ polarized intensity ($PI$) images of the circumstellar disk around the Herbig Ae star AB Aur at a radial distance of 22 AU ($0."15$) up to 554 AU (3.$"$85), which have been obtained by the high-contrast instrument HiCIAO with the dual-beam polarimetry. We revealed complicated and asymmetrical structures in the inner part ($\lesssim$140 AU) of the disk, while confirming the previously reported outer ($r$ $\gtrsim$200 AU) spiral structure. We have imaged a double ring structure at $\sim$40 and $\sim$100 AU and a ring-like gap between the two. We found a significant discrepancy of inclination angles between two rings, which may indicate that the disk of AB Aur is warped. Furthermore, we found seven dips (the typical size is $\sim$45 AU or less) within two rings as well as three prominent $PI$ peaks at $\sim$40 AU. The observed structures, including a bumpy double ring, a ring-like gap, and a warped disk in the innermost regions, provide essential information for understanding the formation mechanism of recently detected wide-orbit ($r$ $>$20 AU) planets.

We present H- and Ks-band imaging data resolving the gap in the transitional disk around LkCa 15, revealing the surrounding nebulosity. We detect sharp elliptical contours delimiting the nebulosity on the inside as well as the outside, consistent with the shape, size, ellipticity, and orientation of starlight reflected from the far-side disk wall, whereas the near-side wall is shielded from view by the disk's optically thick bulk. We note that forward-scattering of starlight on the near-side disk surface could provide an alternate interpretation of the nebulosity. In either case, this discovery provides confirmation of the disk geometry that has been proposed to explain the spectral energy distributions (SED) of such systems, comprising an optically thick outer disk with an inner truncation radius of ~46 AU enclosing a largely evacuated gap. Our data show an offset of the nebulosity contours along the major axis, likely corresponding to a physical pericenter offset of the disk gap. This reinforces the leading theory that dynamical clearing by at least one orbiting body is the cause of the gap. Based on evolutionary models, our high-contrast imagery imposes an upper limit of 21 Jupiter masses on companions at separations outside of 0.1" and of 13 Jupiter masses outside of 0.2". Thus, we find that a planetary system around LkCa 15 is the most likely explanation for the disk architecture.

We present results of direct imaging observations for HAT-P-7 taken with the Subaru HiCIAO and the Calar Alto AstraLux. Since the close-in transiting planet HAT-P-7b was reported to have a highly tilted orbit, massive bodies such as giant planets, brown dwarfs, or a binary star are expected to exist in the outer region of this system. We show that there are indeed two candidates for distant faint stellar companions around HAT-P-7. We discuss possible roles played by such companions on the orbital evolution of HAT-P-7b. We conclude that as there is a third body in the system as reported by Winn et al. (2009, ApJL, 763, L99), the Kozai migration is less likely while planet-planet scattering is possible.