We have obtained Gemini Planet Imager (GPI) J-, H-, K1-, and K2-Spec observations of the iconic debris ring around the young, main-sequence star HR 4796A. We applied several point-spread function (PSF) subtraction techniques to the observations (Mask-and-Interpolate, RDI-NMF, RDI-KLIP, and ADI-KLIP) to measure the geometric parameters and the scattering phase function for the disk. To understand the systematic errors associated with PSF subtraction, we also forward-modeled the observations using a Markov Chain Monte Carlo framework and a simple model for the disk. We found that measurements of the disk geometric parameters were robust, with all of our analyses yielding consistent results; however, measurements of the scattering phase function were challenging to reconstruct from PSF-subtracted images, despite extensive testing. As a result, we estimated the scattering phase function using disk modeling. We searched for a dependence of the scattering phase function with respect to the GPI filters but found none. We compared the H-band scattering phase function with that measured by Hubble Space Telescope STIS at visual wavelengths and discovered a blue color at small scattering angles and a red color at large scattering angles, consistent with predictions and laboratory measurements of large grains. Finally, we successfully modeled the SPHERE H2 HR 4796A scattered phase function using a distribution of hollow spheres composed of silicates, carbon, and metallic iron.

We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager high-contrast integral field spectrograph in $H$-band (1.6 $\mu$m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and three protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks but HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 \deg$). Based on post-survey statistics, we improved upon our pre-survey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing $H$-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.

HR4796A hosts a well-studied debris disk with a long history due to its high fractional luminosity and favorable inclination lending itself well to both unresolved and resolved observations. We present new J- and K1-band images of the resolved debris disk HR4796A taken in the polarimetric mode of the Gemini Planet Imager (GPI). The polarized intensity features a strongly forward scattered brightness distribution and is undetected at the far side of the disk. The total intensity is detected at all scattering angles and also exhibits a strong forward scattering peak. We use a forward modelled geometric disk in order to extract geometric parameters, polarized fraction and total intensity scattering phase functions for these data as well as H-band data previously taken by GPI. We find the polarized phase function becomes increasingly more forward scattering as wavelength increases. We fit Mie and distribution of hollow spheres grain (DHS) models to the extracted functions. We find that while it is possible to describe generate a satisfactory model for the total intensity using a DHS model, but not with a Mie model. We find that no single grain population of DHS or Mie grains of arbitrary composition can simultaneously reproduce the polarized fraction and total intensity scattering phase functions, indicating the need for more sophisticated grain models.

Due to the efforts by numerous ground-based surveys and NASA's Kepler and TESS, there will be hundreds, if not thousands, of transiting exoplanets ideal for atmospheric characterization via spectroscopy with large platforms such as JWST and ARIEL. However their next predicted mid-transit time could become so increasingly uncertain over time that significant overhead would be required to ensure the detection of the entire transit. As a result, follow-up observations to characterize these exoplanetary atmospheres would require less-efficient use of an observatory's time---which is an issue for large platforms where minimizing observing overheads is a necessity. Here we demonstrate the power of citizen scientists operating smaller observatories ($\le$1-m) to keep ephemerides "fresh", defined here as when the 1$\sigma$ uncertainty in the mid-transit time is less than half the transit duration. We advocate for the creation of a community-wide effort to perform ephemeris maintenance on transiting exoplanets by citizen scientists. Such observations can be conducted with even a 6-inch telescope, which has the potential to save up to $\sim$10,000~days for a 1000-planet survey. Based on a preliminary analysis of 14 transits from a single 6-inch MicroObservatory telescope, we empirically estimate the ability of small telescopes to benefit the community. Observations with a small-telescope network operated by citizen scientists are capable of resolving stellar blends to within 5''/pixel, can follow-up long period transits in short-baseline TESS fields, monitor epoch-to-epoch stellar variability at a precision 0.67\%$\pm$0.12\% for a 11.3 V-mag star, and search for new planets or constrain the masses of known planets with transit timing variations greater than two minutes.

We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of $\approx$100au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

We present a revision to the astrometric calibration of the Gemini Planet Imager (GPI), an instrument designed to achieve the high contrast at small angular separations necessary to image substellar and planetary-mass companions around nearby, young stars. We identified several issues with the GPI Data Reduction Pipeline (DRP) that significantly affected the determination of angle of north in reduced GPI images. As well as introducing a small error in position angle measurements for targets observed at small zenith distances, this error led to a significant error in the previous astrometric calibration that has affected all subsequent astrometric measurements. We present a detailed description of these issues, and how they were corrected. We reduced GPI observations of calibration binaries taken periodically since the instrument was commissioned in 2014 using an updated version of the DRP. These measurements were compared to observations obtained with the NIRC2 instrument on Keck II, an instrument with an excellent astrometric calibration, allowing us to derive an updated plate scale and north offset angle for GPI. This revised astrometric calibration should be used to calibrate all measurements obtained with GPI for the purposes of precision astrometry.

We present a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager. HD 165054 is a bright star within Baade's Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($\sim 3$ stars per square arcsecond with $H<22$), necessary because of the small field-of-view of the Gemini Planet Imager. Using nine epochs spanning over five years, we have fit a standard five-parameter astrometric model to the astrometry of seven background stars within close proximity to HD 165054 (angular separation $< 2$ arcsec). We achieved a proper motion precision of $\sim 0.3$ mas/yr, and constrained the parallax of each star to be $\lesssim 1$ mas. Our measured proper motions and parallax limits are consistent with the background stars being a part of the galactic bulge. Using these measurements we find no evidence of any systematic trend of either the plate scale or the north angle offset of GPI between 2014 and 2019. We compared our model describing the motions of the seven background stars to observations of the same field in 2014 and 2018 obtained with Keck/NIRC2, an instrument with an excellent astrometric calibration. We find that predicted position of the background sources is consistent with that measured by NIRC2, within the uncertainties of the calibration of the two instruments. In the future, we will use this field as a standard astrometric calibrator for the upgrade of GPI and potentially for other high-contrast imagers.

We present the first polarimetric detection of the inner disk component around the pre-main sequence B9.5 star HD 141569A. Gemini Planet Imager H-band (1.65 micron) polarimetric differential imaging reveals the highest signal-to-noise ratio detection of this ring yet attained and traces structure inwards to 0.25" (28 AU at a distance of 111 pc). The radial polarized intensity image shows the east side of the disk, peaking in intensity at 0.40" (44 AU) and extending out to 0.9" (100 AU). There is a spiral arm-like enhancement to the south, reminiscent of the known spiral structures on the outer rings of the disk. The location of the spiral arm is coincident with 12CO J=3-2 emission detected by ALMA, and hints at a dynamically active inner circumstellar region. Our observations also show a portion of the middle dusty ring at ~220 AU known from previous observations of this system. We fit the polarized H-band emission with a continuum radiative transfer Mie model. Our best-fit model favors an optically thin disk with a minimum dust grain size close to the blow-out size for this system: evidence of on-going dust production in the inner reaches of the disk. The thermal emission from this model accounts for virtually all of the far-infrared and millimeter flux from the entire HD 141569A disk, in agreement with the lack of ALMA continuum and CO emission beyond ~100 AU. A remaining 8-30 micron thermal excess a factor of ~2 above our model argues for a yet-unresolved warm innermost 5-15 AU component of the disk.

We present new observations of the planet beta Pictoris b from 2018 with GPI, the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit fits. We find a mass consistent with predictions of hot-start evolutionary models and previous works following similar methods, though with larger uncertainties: 12.8 [+5.3, -3.2] M_Jup. Our eccentricity determination of 0.12 [+0.04, -0.03] disfavors circular orbits. We consider orbit fits to several different imaging datasets, and find generally similar posteriors on the mass for each combination of imaging data. Our analysis underscores the importance of performing joint fits to the absolute and relative astrometry simultaneously, given the strong covariance between orbital elements. Time of conjunction is well constrained within 2.8 days of 2017 September 13, with the star behind the planet's Hill sphere between 2017 April 11 and 2018 February 16 (+/- 18 days). Following the recent radial velocity detection of a second planet in the system, beta Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass for the imaged planet beta Pic b, of 8.0 +/- 2.6 M_Jup, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4.

We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet-disk system. Uniquely, HD 98363 has a wide co-moving stellar companion Wray 15-788 with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one non-detection to one of the brightest disks in the region. With the new results, a total of 15 A/F-stars in Sco-Cen have resolved scattered light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F-stars at this important stage in the development of planetary systems.

We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semi-major axis of $11.1_{-1.3}^{+4.2}$\,au and a period of $28.1_{-4.9}^{+17.2}$\,yr, assuming a mass of 1.75\,M$_{\odot}$ for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigated the potential of using absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51~Eri derived from a comparison of the {\it Hipparcos} and {\it Gaia} catalogues was found to be inconsistent at the 2--3$\sigma$ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs or the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using {\it Gaia} astrometry alone for a dynamical mass measurement of the planet by simulating {\it Gaia} measurements of the motion of the photocenter of the system over the course of the extended eight-year mission. We find that such a measurement is only possible ($>98$\% probability) given the most optimistic predictions for the {\it Gaia} scan astrometric uncertainties for bright stars, and a high mass for the planet ($\gtrsim3.6$\,M$_{\rm Jup}$).

The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to be resolved with ground-based high contrast imaging instruments; long-period equal mass ratio stellar companions that are also consistent with the measured acceleration are excluded with previous imaging observations. The small but significant amplitude of the acceleration made HR 1645 a promising candidate for targeted searches for brown dwarf and planetary-mass companions around nearby, young stars. In this paper we explore the origin of the astrometric acceleration by modelling the signal induced by a wide-orbit M8 companion discovered with the Gemini Planet Imager, as well as the effects of an inner short-period spectroscopic companion discovered a century ago but not since followed-up. We present the first constraints on the orbit of the inner companion, and demonstrate that it is a plausible cause of the astrometric acceleration. This result demonstrates the importance of vetting of targets with measured astrometric acceleration for short-period stellar companions prior to conducting targeted direct imaging surveys for wide-orbit substellar companions.

An explanation for the origin of asymmetry along the preferential axis of the PSF of an AO system is developed. When phase errors from high altitude turbulence scintillate due to Fresnel propagation, wavefront amplitude errors may be spatially offset from residual phase errors. These correlated errors appear as asymmetry in the image plane under the Fraunhofer condition. In an analytic model with an open-loop AO system, the strength of the asymmetry is calculated for a single mode of phase aberration, which generalizes to two dimensions under a Fourier decomposition of the complex illumination. Other parameters included are the spatial offset of the AO correction, which is the wind velocity in the frozen flow regime multiplied by the effective AO time delay, and propagation distance or altitude of the turbulent layer. In this model, the asymmetry is strongest when the wind is slow and nearest to the coronagraphic mask when the turbulent layer is far away, such as when the telescope is pointing low towards the horizon. A great emphasis is made about the fact that the brighter asymmetric lobe of the PSF points in the opposite direction as the wind, which is consistent analytically with the clarification that the image plane electric field distribution is actually the inverse Fourier transform of the aperture plane. Validation of this understanding is made with observations taken from the Gemini Planet Imager, as well as being reproducible in end-to-end AO simulations.

The Keck Planet Imager and Characterizer (KPIC) is an upgrade to the Keck II adaptive optics system enabling high contrast imaging and high-resolution spectroscopic characterization of giant exoplanets in the mid-infrared (2-5 microns). The KPIC instrument will be developed in phases. Phase I entails the installation of an infrared pyramid wavefront sensor (PyWFS) based on a fast, low-noise SAPHIRA IR-APD array. The ultra-sensitive infrared PyWFS will enable high contrast studies of infant exoplanets around cool, red, and/or obscured targets in star forming regions. In addition, the light downstream of the PyWFS will be coupled into an array of single-mode fibers with the aid of an active fiber injection unit (FIU). In turn, these fibers route light to Keck's high-resolution infrared spectrograph NIRSPEC, so that high dispersion coronagraphy (HDC) can be implemented for the first time. HDC optimally pairs high contrast imaging and high-resolution spectroscopy allowing detailed characterization of exoplanet atmospheres, including molecular composition, spin measurements, and Doppler imaging. Here we provide an overview of the instrument, its science scope, and report on recent results from on-sky commissioning of Phase I. The instrument design and techniques developed will be key for more advanced instrument concepts needed for the extremely large telescopes of the future.

The fields of Astronomy and Astrophysics are technology limited, where the advent and application of new technologies to astronomy usher in a flood of discoveries altering our understanding of the Universe (e.g., recent cases include LIGO and the GRAVITY instrument at the VLTI). Currently, the field of astronomical spectroscopy is rapidly approaching an impasse: the size and cost of instruments, especially multi-object and integral field spectrographs for extremely large telescopes (ELTs), are pushing the limits of what is feasible, requiring optical components at the very edge of achievable size and performance. For these reasons, astronomers are increasingly looking for innovative solutions like photonic technologies that promote instrument miniaturization and simplification, while providing superior performance. Astronomers have long been aware of the potential of photonic technologies. The goal of this white paper is to draw attention to key photonic technologies and developments over the past two decades and demonstrate there is new momentum in this arena. We outline where the most critical efforts should be focused over the coming decade in order to move towards realizing a fully photonic instrument. A relatively small investment in this technology will advance astronomical photonics to a level where it can reliably be used to solve challenging instrument design limitations. For the benefit of both ground and space borne instruments alike, an endorsement from the National Academy of Sciences decadal survey will ensure that such solutions are set on a path to their full scientific exploitation, which may one day address a broad range of science cases outlined in the KSPs.

HISPEC (High-resolution Infrared Spectrograph for Exoplanet Characterization) is a proposed diffraction-limited spectrograph for the W.M. Keck Observatory, and a pathfinder for the MODHIS facility project (Multi-Object Diffraction-limited High-resolution Infrared Spectrograph) on the Thirty Meter Telescope. HISPEC/MODHIS builds on diffraction-limited spectrograph designs which rely on adaptively corrected single-mode fiber feeds. Seeing-limited high-resolution spectrographs, by virtue of the conservation of beam etendue, grow in volume following a D^3 power law (D is the telescope diameter), and are subject to daunting challenges associated with their large size. Diffraction-limited spectrographs fed by single mode fibers are decoupled from the telescope input, and are orders of magnitude more compact and have intrinsically stable line spread functions. Their efficiency is directly proportional to the performance of the adaptive optics (AO) system. AO technologies have matured rapidly over the past two decades and are baselined for future extremely large telescopes. HISPEC/MODHIS will take R>100,000 spectra of a few objects in a 10" field-of-view sampled at the diffraction limit (~10-50 mas), simultaneously from 0.95 to 2.4 microns (y-K). The scientific scope ranges from exoplanet infrared precision radial velocities, spectroscopy of transiting, close-in, and directly imaged exoplanets (atmospheric composition and dynamics, RM effect, spin measurements, Doppler imaging), brown dwarf characterization, stellar physics/chemistry, proto-planetary disk kinematics/composition, Solar system, extragalactic science, and cosmology. HISPEC/MODHIS features a compact, cost-effective design optimized to fully exploit the existing Keck-AO and future TMT-NFIRAOS infrastructures and boost the scientific reach of Keck Observatory and TMT soon after first light.

We have obtained Hubble Space Telescope STIS and NICMOS, and Gemini/GPI scattered light images of the HD 191089 debris disk. We identify two spatial components: a ring resembling Kuiper Belt in radial extent (FWHM: ${\sim}$25 au, centered at ${\sim}$46 au), and a halo extending to ${\sim}$640 au. We find that the halo is significantly bluer than the ring, consistent with the scenario that the ring serves as the "birth ring" for the smaller dust in the halo. We measure the scattering phase functions in the 30{\deg}-150{\deg} scattering angle range and find the halo dust is both more forward- and backward-scattering than the ring dust. We measure a surface density power law index of -0.68${\pm}$0.04 for the halo, which indicates the slow-down of the radial outward motion of the dust. Using radiative transfer modeling, we attempt to simultaneously reproduce the (visible) total and (near-infrared) polarized intensity images of the birth ring. Our modeling leads to mutually inconsistent results, indicating that more complex models, such as the inclusion of more realistic aggregate particles, are needed.

Over the past three decades instruments on the ground and in space have discovered thousands of planets outside the solar system. These observations have given rise to an astonishingly detailed picture of the demographics of short-period planets, but are incomplete at longer periods where both the sensitivity of transit surveys and radial velocity signals plummet. Even more glaring is that the spectra of planets discovered with these indirect methods are either inaccessible (radial velocity detections) or only available for a small subclass of transiting planets with thick, clear atmospheres. Direct detection can be used to discover and characterize the atmospheres of planets at intermediate and wide separations, including non-transiting exoplanets. Today, a small number of exoplanets have been directly imaged, but they represent only a rare class of young, self-luminous super-Jovian-mass objects orbiting tens to hundreds of AU from their host stars. Atmospheric characterization of planets in the <5 AU regime, where radial velocity (RV) surveys have revealed an abundance of other worlds, is technically feasible with 30-m class apertures in combination with an advanced AO system, coronagraph, and suite of spectrometers and imagers. There is a vast range of unexplored science accessible through astrometry, photometry, and spectroscopy of rocky planets, ice giants, and gas giants. In this whitepaper we will focus on one of the most ambitious science goals --- detecting for the first time habitable-zone rocky (<1.6 R_Earth) exoplanets in reflected light around nearby M-dwarfs

The Gemini Planet Imager (GPI) contains a 10-hole non-redundant mask (NRM), enabling interferometric resolution in complement to its coronagraphic capabilities. The NRM operates both in spectroscopic (integral field spectrograph, henceforth IFS) and polarimetric configurations. NRM observations were taken between 2013 and 2016 to characterize its performance. Most observations were taken in spectroscopic mode with the goal of obtaining precise astrometry and spectroscopy of faint companions to bright stars. We find a clear correlation between residual wavefront error measured by the AO system and the contrast sensitivity by comparing phase errors in observations of the same source, taken on different dates. We find a typical 5-$\sigma$ contrast sensitivity of $2-3~\times~10^{-3}$ at $\sim\lambda/D$. We explore the accuracy of spectral extraction of secondary components of binary systems by recovering the signal from a simulated source injected into several datasets. We outline data reduction procedures unique to GPI's IFS and describe a newly public data pipeline used for the presented analyses. We demonstrate recovery of astrometry and spectroscopy of two known companions to HR 2690 and HD 142527. NRM+polarimetry observations achieve differential visibility precision of $\sigma\sim0.4\%$ in the best case. We discuss its limitations on Gemini-S/GPI for resolving inner regions of protoplanetary disks and prospects for future upgrades. We summarize lessons learned in observing with NRM in spectroscopic and polarimetric modes.

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^\alpha a^\beta$, finding $\alpha$ = -2.4 $\pm$ 0.8 and $\beta$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.

Planned astronomical observatories of the 2020s will be capable of obtaining reflected light photometry and spectroscopy of cool extrasolar giant planets. Here we explain that such data are valuable both for understanding the origin and evolution of giant planets as a whole and for preparing for the interpretation of similar datasets from potentially habitable extrasolar terrestrial planets in the decades to follow.

We have now accumulated a wealth of observations of the planet-formation environment and of mature planetary systems. These data allow us to test and refine theories of gas-giant planet formation by placing constraints on the conditions and timescale of this process. Yet a number of fundamental questions remain unanswered about how protoplanets accumulate material, their photospheric properties and compositions, and how they interact with protoplanetary disks. While we have begun to detect protoplanet candidates during the last several years, we are presently only sensitive to the widest separation, highest mass / accretion rate cases. Current observing facilities lack the angular resolution and inner working angle to probe the few-AU orbital separations where giant planet formation is thought to be most efficient. They also lack the contrast to detect accretion rates that would form lower mass gas giants and ice giants. Instruments and telescopes coming online over the next decade will provide high contrast in the inner giant-planet-forming regions around young stars, allowing us to build a protoplanet census and to characterize planet formation in detail for the first time.

OSIRIS is a near-infrared (1.0--2.4 $\mu$m) integral field spectrograph operating behind the adaptive optics system at Keck Observatory, and is one of the first lenslet-based integral field spectrographs. Since its commissioning in 2005, it has been a productive instrument, producing nearly half the laser guide star adaptive optics (LGS AO) papers on Keck. The complexity of its raw data format necessitated a custom data reduction pipeline (DRP) delivered with the instrument in order to iteratively assign flux in overlapping spectra to the proper spatial and spectral locations in a data cube. Other than bug fixes and updates required for hardware upgrades, the bulk of the DRP has not been updated since initial instrument commissioning. We report on the first major comprehensive characterization of the DRP using on-sky and calibration data. We also detail improvements to the DRP including characterization of the flux assignment algorithm; exploration of spatial rippling in the reduced data cubes; and improvements to several calibration files, including the rectification matrix, the bad pixel mask, and the wavelength solution. We present lessons learned from over a decade of OSIRIS data reduction that are relevant to the next generation of integral field spectrograph hardware and data reduction software design.

The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assuming the planets are near 1:2:4:8 period commensurabilities, or are coplanar, does not worsen the fit. We added the prior that the planets must have been stable for the age of the system (40 Myr) by running orbit configurations from our posteriors through $N$-body simulations and varying the masses of the planets. We found that only assuming the planets are both coplanar and near 1:2:4:8 period commensurabilities produces dynamically stable orbits in large quantities. Our posterior of stable coplanar orbits tightly constrains the planets' orbits, and we discuss implications for the outermost planet b shaping the debris disk. A four-planet resonance lock is not necessary for stability up to now. However, planet pairs d and e, and c and d, are each likely locked in two-body resonances for stability if their component masses are above $6~M_{\rm{Jup}}$ and $7~M_{\rm{Jup}}$, respectively. Combining the dynamical and luminosity constraints on the masses using hot-start evolutionary models and a system age of $42 \pm 5$~Myr, we found the mass of planet b to be $5.8 \pm 0.5~M_{\rm{Jup}}$, and the masses of planets c, d, and e to be $7.2_{-0.7}^{+0.6}~M_{\rm{Jup}}$ each.

Direct detection, also known as direct imaging, is a method for discovering and characterizing the atmospheres of planets at intermediate and wide separations. It is the only means of obtaining spectra of non-transiting exoplanets. Characterizing the atmospheres of planets in the <5 AU regime, where RV surveys have revealed an abundance of other worlds, requires a 30-m-class aperture in combination with an advanced adaptive optics system, coronagraph, and suite of spectrometers and imagers - this concept underlies planned instruments for both TMT (the Planetary Systems Imager, or PSI) and the GMT (GMagAO-X). These instruments could provide astrometry, photometry, and spectroscopy of an unprecedented sample of rocky planets, ice giants, and gas giants. For the first time habitable zone exoplanets will become accessible to direct imaging, and these instruments have the potential to detect and characterize the innermost regions of nearby M-dwarf planetary systems in reflected light. High-resolution spectroscopy will not only illuminate the physics and chemistry of exo-atmospheres, but may also probe rocky, temperate worlds for signs of life in the form of atmospheric biomarkers (combinations of water, oxygen and other molecular species). By completing the census of non-transiting worlds at a range of separations from their host stars, these instruments will provide the final pieces to the puzzle of planetary demographics. This whitepaper explores the science goals of direct imaging on 30-m telescopes and the technology development needed to achieve them.

The Star-Planet Activity Research CubeSat (SPARCS) is a NASA-funded astrophysics mission, devoted to the study of the ultraviolet (UV) time-domain behavior in low-mass stars. Given their abundance and size, low-mass stars are important targets in the search for habitable-zone, exoplanets. However, not enough is known about the stars flare and quiescent emission, which powers photochemical reactions on the atmospheres of possible planets. Over its initial 1-year mission, SPARCS will stare at ~10 stars in order to measure short- (minutes) and long- (months) term variability simultaneously in the near-UV (NUV - lam = 280 nm) and far-UV (FUV - lam = 162 nm). The SPARCS payload consists of a 9-cm reflector telescope paired with two high-sensitivity 2D-doped CCDs. The detectors are kept passively cooled at 238K, in order to reduce dark-current contribution. The filters have been selected to provide strong rejection of longer wavelengths, where most of the starlight is emitted. The payload will be integrated within a 6U CubeSat to be placed on a Sun-synchronous terminator orbit, allowing for long observing stares for all targets. Launch is expected to occur not earlier than October 2021.

NIRSPEC is a 1-5 micron echelle spectrograph in use on the Keck II Telescope since 1999. The spectrograph is capable of both moderate (R~2,000) and high (R~25,000) resolution observations and has been a workhorse instrument across many astronomical fields, from planetary science to extragalactic observations. In the latter half of 2018, we will upgrade NIRSPEC to improve the sensitivity and stability of the instrument and increase its lifetime. The major components of the upgrade include replacing the spectrometer and slit-viewing camera detectors with Teledyne H2RG arrays and replacing all transputer-based electronics. We present detailed design, testing, and analysis of the upgraded instrument, including the finalized optomechanical design of the new 1-5 micron slit-viewing camera, final alignment and assembly of the science array, electronics systems, and updated software design.

We summarize the red channel (2-5 micron) of the Planetary Systems Imager (PSI), a proposed second-generation instrument for the TMT. Cold exoplanets emit the majority of their light in the thermal infrared, which means these exoplanets can be detected at a more modest contrast than at other wavelengths. PSI-Red will be able to detect and characterize a wide variety of exoplanets, including radial-velocity planets on wide orbits, accreting protoplanets in nearby star-forming regions, and reflected-light planets around the nearest stars. PSI-Red will feature an imager, a low-resolution lenslet integral field spectrograph, a medium-resolution lenslet+slicer integral field spectrograph, and a fiber-fed high-resolution spectrograph.

The success of ground-based, high contrast imaging for the detection of exoplanets in part depends on the ability to differentiate between quasi-static speckles caused by aberrations not corrected by adaptive optics (AO) systems, known as non-common path aberrations (NCPAs), and the planet intensity signal. Frazin (ApJ, 2013) introduced a post-processing algorithm demonstrating that simultaneous millisecond exposures in the science camera and wavefront sensor (WFS) can be used with a statistical inference procedure to determine both the series expanded NCPA coefficients and the planetary signal. We demonstrate, via simulation, that using this algorithm in a closed-loop AO system, real-time estimation and correction of the quasi-static NCPA is possible without separate deformable mirror (DM) probes. Thus the use of this technique allows for the removal of the quasi-static speckles that can be mistaken for planetary signals without the need for new optical hardware, improving the efficiency of ground-based exoplanet detection. In our simulations, we explore the behavior of the Frazin Algorithm (FA) and the dependence of its convergence to an accurate estimate on factors such as Strehl ratio, NCPA strength, and number of algorithm search basis functions. We then apply this knowledge to simulate running the algorithm in real-time in a nearly ideal setting. We then discuss adaptations that can be made to the algorithm to improve its real-time performance, and show their efficacy in simulation. A final simulation tests the technique's resilience against imperfect knowledge of the AO residual phase, motivating an analysis of the feasibility of using this technique in a real closed-loop Extreme AO system such as SCExAO or MagAO-X, in terms of computational complexity and the accuracy of the estimated quasi-static NCPA correction.

We present new high resolution imaging of a light-scattering dust ring and halo around the young star HD 35841. Using spectroscopic and polarimetric data from the Gemini Planet Imager in H-band (1.6 microns), we detect the highly inclined (i=85 deg) ring of debris down to a projected separation of ~12 au (~0.12") for the first time. Optical imaging from HST/STIS shows a smooth dust halo extending outward from the ring to >140 au (>1.4"). We measure the ring's scattering phase function and polarization fraction over scattering angles of 22-125 deg, showing a preference for forward scattering and a polarization fraction that peaks at ~30% near the ansae. Modeling of the scattered-light disk indicates that the ring spans radii of ~60-220 au, has a vertical thickness similar to that of other resolved dust rings, and contains grains as small as 1.5 microns in diameter. These models also suggest the grains have a low porosity, are more likely to consist of carbon than astrosilicates, and contain significant water ice. The halo has a surface brightness profile consistent with that expected from grains pushed by radiation pressure from the main ring onto highly eccentric but still bound orbits. We also briefly investigate arrangements of a possible inner disk component implied by our spectral energy distribution models, and speculate about the limitations of Mie theory for doing detailed analyses of debris disk dust populations.

We explore KLIP forward modeling spectral extraction on Gemini Planet Imager coronagraphic data of HR 8799, using PyKLIP and show algorithm stability with varying KLIP parameters. We report new and re-reduced spectrophotometry of HR 8799 c, d, and e in H & K bands. We discuss a strategy for choosing optimal KLIP PSF subtraction parameters by injecting simulated sources and recovering them over a range of parameters. The K1/K2 spectra for HR 8799 c and d are similar to previously published results from the same dataset. We also present a K band spectrum of HR 8799 e for the first time and show that our H-band spectra agree well with previously published spectra from the VLT/SPHERE instrument. We show that HR 8799 c and d show significant differences in their H & K spectra, but do not find any conclusive differences between d and e or c and e, likely due to large error bars in the recovered spectrum of e. Compared to M, L, and T-type field brown dwarfs, all three planets are most consistent with mid and late L spectral types. All objects are consistent with low gravity but a lack of standard spectra for low gravity limit the ability to fit the best spectral type. We discuss how dedicated modeling efforts can better fit HR 8799 planets' near-IR flux and discuss how differences between the properties of these planets can be further explored.

Detection of weak radial velocity shifts of host stars induced by orbiting planets is an important technique for discovering and characterizing planets beyond our solar system. Optical frequency combs enable calibration of stellar radial velocity shifts at levels required for detection of Earth analogs. A new chip-based device, the Kerr soliton microcomb, has properties ideal for ubiquitous application outside the lab and even in future space-borne instruments. Moreover, microcomb spectra are ideally suited for astronomical spectrograph calibration and eliminate filtering steps required by conventional mode-locked-laser frequency combs. Here, for the calibration of astronomical spectrographs, we demonstrate an atomic/molecular line-referenced, near-infrared soliton microcomb. Efforts to search for the known exoplanet HD 187123b were conducted at the Keck-II telescope as a first in-the-field demonstration of microcombs.

The Gemini Planet Imager Exoplanet Survey (GPIES) is a multi-year direct imaging survey of 600 stars to discover and characterize young Jovian exoplanets and their environments. We have developed an automated data architecture to process and index all data related to the survey uniformly. An automated and flexible data processing framework, which we term the Data Cruncher, combines multiple data reduction pipelines together to process all spectroscopic, polarimetric, and calibration data taken with GPIES. With no human intervention, fully reduced and calibrated data products are available less than an hour after the data are taken to expedite follow-up on potential objects of interest. The Data Cruncher can run on a supercomputer to reprocess all GPIES data in a single day as improvements are made to our data reduction pipelines. A backend MySQL database indexes all files, which are synced to the cloud, and a front-end web server allows for easy browsing of all files associated with GPIES. To help observers, quicklook displays show reduced data as they are processed in real-time, and chatbots on Slack post observing information as well as reduced data products. Together, the GPIES automated data processing architecture reduces our workload, provides real-time data reduction, optimizes our observing strategy, and maintains a homogeneously reduced dataset to study planet occurrence and instrument performance.

We present evidence that the recently discovered, directly-imaged planet HD 131399 Ab is a background star with non-zero proper motion. From new JHK1L' photometry and spectroscopy obtained with the Gemini Planet Imager, VLT/SPHERE, and Keck/NIRC2, and a reanalysis of the discovery data obtained with VLT/SPHERE, we derive colors, spectra, and astrometry for HD 131399 Ab. The broader wavelength coverage and higher data quality allow us to re-investigate its status. Its near-infrared spectral energy distribution excludes spectral types later than L0 and is consistent with a K or M dwarf, which are the most likely candidates for a background object in this direction at the apparent magnitude observed. If it were a physically associated object, the projected velocity of HD 131399 Ab would exceed escape velocity given the mass and distance to HD 131399 A. We show that HD 131399 Ab is also not following the expected track for a stationary background star at infinite distance. Solving for the proper motion and parallax required to explain the relative motion of HD 131399 Ab, we find a proper motion of 12.3 mas/yr. When compared to predicted background objects drawn from a galactic model, we find this proper motion to be high, but consistent with the top 4% fastest-moving background stars. From our analysis we conclude that HD 131399 Ab is a background K or M dwarf.

We present ALMA mosaic observations at 1.3 mm (223 GHz) of the Fomalhaut system with a sensitivity of 14 $\mu$Jy/beam. These observations provide the first millimeter map of the continuum dust emission from the complete outer debris disk with uniform sensitivity, enabling the first conclusive detection of apocenter glow. We adopt a MCMC modeling approach that accounts for the eccentric orbital parameters of a collection of particles within the disk. The outer belt is radially confined with an inner edge of $136.3\pm0.9$ AU and width of $13.5\pm1.8$ AU. We determine a best-fit eccentricity of $0.12\pm0.01$. Assuming a size distribution power law index of $q=3.46\pm 0.09$, we constrain the dust absorptivity power law index $\beta$ to be $0.9<\beta<1.5$. The geometry of the disk is robustly constrained with inclination $65.\!\!^\circ6\pm0.\!\!^\circ3$, position angle $337.\!\!^\circ9\pm0.\!\!^\circ3$, and argument of periastron $22.\!\!^\circ5\pm4.\!\!^\circ3$. Our observations do not confirm any of the azimuthal features found in previous imaging studies of the disk with HST, SCUBA, and ALMA. However, we cannot rule out structures $\leq10$ AU in size or which only affect smaller grains. The central star is clearly detected with a flux density of $0.75\pm0.02$ mJy, significantly lower than predicted by current photospheric models. We discuss the implications of these observations for the directly imaged Fomalhaut b and the inner dust belt detected at infrared wavelengths.

Recent ALMA observations present mounting evidence for the presence of exocometary gas released within Kuiper belt analogues around nearby main sequence stars. This represents a unique opportunity to study their ice reservoir at the younger ages when volatile delivery to planets is most likely to occur. We here present the detection of CO J=2-1 emission co-located with dust emission from the cometary belt in the 440 Myr-old Fomalhaut system. Through spectro-spatial filtering, we achieve a 5.4$\sigma$ detection and determine that the ring's sky-projected rotation axis matches that of the star. The CO mass derived ($0.65-42 \times10^{-7}$ M$_{\oplus}$) is the lowest of any circumstellar disk detected to date, and must be of exocometary origin. Using a steady state model, we estimate the CO+CO$_2$ mass fraction of exocomets around Fomalhaut to be between 4.6-76%, consistent with Solar System comets and the two other belts known to host exocometary gas. This is the first indication of a similarity in cometary compositions across planetary systems that may be linked to their formation scenario and is consistent with direct ISM inheritance. In addition, we find tentative evidence that $(49\pm 27)$% of the detected flux originates from a region near the eccentric belt's pericentre. If confirmed, the latter may be explained through a recent impact event or CO pericentre glow due to exocometary release within a steady state collisional cascade. In the latter scenario, we show how the azimuthal dependence of the CO release rate leads to asymmetries in gas observations of eccentric exocometary belts.

We present a new matched filter algorithm for direct detection of point sources in the immediate vicinity of bright stars. The stellar Point Spread Function (PSF) is first subtracted using a Karhunen-Lo\'eve Image Processing (KLIP) algorithm with Angular and Spectral Differential Imaging (ADI and SDI). The KLIP-induced distortion of the astrophysical signal is included in the matched filter template by computing a forward model of the PSF at every position in the image. To optimize the performance of the algorithm, we conduct extensive planet injection and recovery tests and tune the exoplanet spectra template and KLIP reduction aggressiveness to maximize the Signal-to-Noise Ratio (SNR) of the recovered planets. We show that only two spectral templates are necessary to recover any young Jovian exoplanets with minimal SNR loss. We also developed a complete pipeline for the automated detection of point source candidates, the calculation of Receiver Operating Characteristics (ROC), false positives based contrast curves, and completeness contours. We process in a uniform manner more than 330 datasets from the Gemini Planet Imager Exoplanet Survey (GPIES) and assess GPI typical sensitivity as a function of the star and the hypothetical companion spectral type. This work allows for the first time a comparison of different detection algorithms at a survey scale accounting for both planet completeness and false positive rate. We show that the new forward model matched filter allows the detection of $50\%$ fainter objects than a conventional cross-correlation technique with a Gaussian PSF template for the same false positive rate.

We present spectro-photometry spanning 1-5 $\mu$m of 51 Eridani b, a 2-10 M$_\text{Jup}$ planet discovered by the Gemini Planet Imager Exoplanet Survey. In this study, we present new $K1$ (1.90-2.19 $\mu$m) and $K2$ (2.10-2.40 $\mu$m) spectra taken with the Gemini Planet Imager as well as an updated $L_P$ (3.76 $\mu$m) and new $M_S$ (4.67 $\mu$m) photometry from the NIRC2 Narrow camera. The new data were combined with $J$ (1.13-1.35 $\mu$m) and $H$ (1.50-1.80 $\mu$m) spectra from the discovery epoch with the goal of better characterizing the planet properties. 51 Eri b photometry is redder than field brown dwarfs as well as known young T-dwarfs with similar spectral type (between T4-T8) and we propose that 51 Eri b might be in the process of undergoing the transition from L-type to T-type. We used two complementary atmosphere model grids including either deep iron/silicate clouds or sulfide/salt clouds in the photosphere, spanning a range of cloud properties, including fully cloudy, cloud free and patchy/intermediate opacity clouds. Model fits suggest that 51 Eri b has an effective temperature ranging between 605-737 K, a solar metallicity, a surface gravity of $\log$(g) = 3.5-4.0 dex, and the atmosphere requires a patchy cloud atmosphere to model the SED. From the model atmospheres, we infer a luminosity for the planet of -5.83 to -5.93 ($\log L/L_{\odot}$), leaving 51 Eri b in the unique position as being one of the only directly imaged planet consistent with having formed via cold-start scenario. Comparisons of the planet SED against warm-start models indicates that the planet luminosity is best reproduced by a planet formed via core accretion with a core mass between 15 and 127 M$_{\oplus}$.

We present H band spectroscopic and Halpha photometric observations of HD 100546 obtained with GPI and MagAO. We detect H band emission at the location of the protoplanet HD 100546b, but show that choice of data processing parameters strongly affects the morphology of this source. It appears point-like in some aggressive reductions, but rejoins an extended disk structure in the majority of the others. Furthermore, we demonstrate that this emission appears stationary on a timescale of 4.6 yrs, inconsistent at the 2sigma level with a Keplerian clockwise orbit at 59 au in the disk plane. The H band spectrum of the emission is inconsistent with any type of low effective temperature object or accreting protoplanetary disk. It strongly suggests a scattered light origin, as it is consistent with the spectrum of the star and the spectra extracted at other locations in the disk. A non detection at the 5sigma level of HD 100546b in differential Halpha imaging places an upper limit, assuming the protoplanet lies in a gap free of extinction, on the accretion luminosity and accretion rate of 1.7E-4 Lsun and MMdot<6.4E-7Mjup^2/yr for 1Rjup. These limits are comparable to the accretion luminosity and rate of TTauri-stars or LkCa 15b. Taken together, these lines of evidence suggest that the H band source at the location of HD 100546b is not emitted by a planetary photosphere or an accreting circumplanetary disk but is a disk feature enhanced by the PSF subtraction process. This non-detection is consistent with the non-detection in the K band reported in an earlier study but does not exclude the possibility that HD 100546b is deeply embedded.

We present optical and near-infrared high contrast images of the transitional disk HD 100546 taken with the Magellan Adaptive Optics system (MagAO) and the Gemini Planet Imager (GPI). GPI data include both polarized intensity and total intensity imagery, and MagAO data are taken in Simultaneous Differential Imaging mode at H{\alpha}. The new GPI H -band total intensity data represent a significant enhancement in sensitivity and field rotation compared to previous data sets and enable a detailed exploration of substructure in the disk. The data are processed with a variety of differential imaging techniques (polarized, angular, reference, and simultaneous differential imaging) in an attempt to identify the disk structures that are most consistent across wavelengths, processing techniques, and algorithmic parameters. The inner disk cavity at 15 au is clearly resolved in multiple datasets, as are a variety of spiral features. While the cavity and spiral structures are identified at levels significantly distinct from the neighboring regions of the disk under several algorithms and with a range of algorithmic parameters, emission at the location of HD 100546 c varies from point-like under aggressive algorithmic parameters to a smooth continuous structure with conservative parameters, and is consistent with disk emission. Features identified in the HD100546 disk bear qualitative similarity to computational models of a moderately inclined two-armed spiral disk, where projection effects and wrapping of the spiral arms around the star result in a number of truncated spiral features in forward-modeled images.

We present new observations of the low-mass companion to HD 984 taken with the Gemini Planet Imager as a part of the Gemini Planet Imager Exoplanet Survey campaign. Images of HD 984 B were obtained in the J (1.12--1.3 micron) and H (1.50--1.80 micron) bands. Combined with archival epochs from 2012 and 2014, we fit the first orbit to the companion to find an 18 AU (70 year) orbit with a 68% confidence interval between 14 and 28 AU, an eccentricity of 0.18 with a 68% confidence interval between 0.05 and 0.47, and an inclination of 119 degrees with a 68% confidence interval between 114 degrees and 125 degrees. To address considerable spectral covariance in both spectra, we present a method of splitting the spectra into low and high frequencies to analyze the spectral structure at different spatial frequencies with the proper spectral noise correlation. Using the split spectra, we compare to known spectral types using field brown dwarf and low-mass star spectra and find a best fit match of a field gravity M6.5+/-1.5 spectral type with a corresponding temperature of 2730+120 K. Photometry of the companion yields a luminosity of log(L_bol/L_sun) = -2.88+/-0.07 dex, using DUSTY models. Mass estimates, again from DUSTY models, find an age-dependent mass of 34+/-1 to 95+/-4 M_Jup. These results are consistent with previous measurements of the object.

Using the Gemini Planet Imager (GPI) located at Gemini South, we measured the near-infrared (1.0-2.4 micron) spectrum of the planetary companion to the nearby, young star $\beta$ Pictoris. We compare the spectrum obtained with currently published model grids and with known substellar objects and present the best matching models as well as the best matching observed objects. Comparing the empirical measurement of the bolometric luminosity to evolutionary models, we find a mass of $12.9\pm0.2$ $\mathcal{M}_\mathrm{Jup}$, an effective temperature of $1724\pm15$ K, a radius of $1.46\pm0.01$ $\mathcal{R}_\mathrm{Jup}$, and a surface gravity of $\log g = 4.18\pm0.01$ [dex] (cgs). The stated uncertainties are statistical errors only, and do not incorporate any uncertainty on the evolutionary models. Using atmospheric models, we find an effective temperature of $1700-1800$ K and a surface gravity of $\log g = 3.5$-$4.0$ [dex] depending upon model. These values agree well with other publications and with "hot-start" predictions from planetary evolution models. Further, we find that the spectrum of $\beta$ Pic b best matches a low-surface gravity L2$\pm$1 brown dwarf. Finally comparing the spectrum to field brown dwarfs we find the the spectrum best matches 2MASS J04062677-381210 and 2MASS J03552337+1133437.

Models of debris disk morphology are often focused on the effects of a planet orbiting interior to or within the disk. Nonetheless, an exterior planetary-mass perturber can also excite eccentricities in a debris disk, via Laplace-Lagrange secular perturbations in the coplanar case or Kozai-Lidov perturbations for mutually inclined companions and disks. HD 106906 is an ideal example of such a system, as it harbors a confirmed exterior 11 M_Jup companion at a projected separation of 650 au outside a resolved, asymmetric disk. We use collisional and dynamical simulations to investigate the interactions between the disk and the companion, and to use the disk's observed morphology to place constraints on the companion's orbit. We conclude that the disk's observed morphology is consistent with perturbations from the observed exterior companion. Generalizing this result, we suggest that exterior perturbers, as well as interior planets, should be considered when investigating the cause of observed asymmetries in a debris disk.

The power of robotic telescopes to transform science education has been voiced by multiple sources, since the 1980s. Since then, much technical progress has been made in robotic telescope provision to end users via a variety of different approaches. The educational transformation hoped for by the provision of this technology has, so far, yet to be achieved on a scale matching the technical advancements. In this paper, the history, definition, role and rationale of optical robotic telescopes with a focus on their use in education is provided. The current telescope access providers and educational projects and their broad uses in traditional schooling, undergraduate and outreach are then outlined. From this background, the current challenges to the field, which are numerous, are then presented. This review is concluded with a series of recommendations for current and future projects that are apparent and have emerged from the literature.

We present new spatially resolved astrometry and photometry from the Gemini Planet Imager of the inner binary of the young multiple star system V343 Normae, which is a member of the beta Pictoris moving group. V343 Normae comprises a K0 and mid-M star in a ~4.5 year orbit (AaAb) and a wide 10" M5 companion (B). By combining these data with archival astrometry and radial velocities we fit the orbit and measure individual masses for both components of M_Aa = 1.10 +/- 0.10 M_sun and M_Ab = 0.290 +/- 0.018 M_sun. Comparing to theoretical isochrones, we find good agreement for the measured masses and JHK band magnitudes of the two components consistent with the age of the beta Pic moving group. We derive a model-dependent age for the beta Pic moving group of 26 +/- 3 Myr by combining our results for V343 Normae with literature measurements for GJ 3305, which is another group member with resolved binary components and dynamical masses.

The Gemini Planet Imager (GPI) has been designed for the direct detection and characterization of exoplanets and circumstellar disks. GPI is equipped with a dual channel polarimetry mode designed to take advantage of the inherently polarized light scattered off circumstellar material to further suppress the residual seeing halo left uncorrected by the adaptive optics. We explore how recent advances in data reduction techniques reduce systematics and improve the achievable contrast in polarimetry mode. In particular, we consider different flux extraction techniques when constructing datacubes from raw data, division by a polarized flat-field and a method for subtracting instrumental polarization. Using observations of unpolarized standard stars we find that GPI's instrumental polarization is consistent with being wavelength independent within our errors. In addition, we provide polarimetry contrast curves that demonstrate typical performance throughout the GPIES campaign.

The HD 61005 debris disk ("The Moth") stands out from the growing collection of spatially resolved circumstellar disks by virtue of its unusual swept-back morphology, brightness asymmetries, and dust ring offset. Despite several suggestions for the physical mechanisms creating these features, no definitive answer has been found. In this work, we demonstrate the plausibility of a scenario in which the disk material is shaped dynamically by an eccentric, inclined planet. We present new Keck NIRC2 scattered-light angular differential imaging of the disk at 1.2-2.3 microns that further constrains its outer morphology (projected separations of 27-135 AU). We also present complementary Gemini Planet Imager 1.6 micron total intensity and polarized light detections that probe down to projected separations less than 10 AU. To test our planet-sculpting hypothesis, we employed secular perturbation theory to construct parent body and dust distributions that informed scattered-light models. We found that this method produced models with morphological and photometric features similar to those seen in the data, supporting the premise of a planet-perturbed disk. Briefly, our results indicate a disk parent body population with a semimajor axis of 40-52 AU and an interior planet with an eccentricity of at least 0.2. Many permutations of planet mass and semimajor axis are allowed, ranging from an Earth mass at 35 AU to a Jupiter mass at 5 AU.

We present $H$-band near-infrared polarimetric imaging observations of the F5V star HD~157587 obtained with the Gemini Planet Imager (GPI) that reveal the debris disk as a bright ring structure at a separation of $\sim$80$-$100~AU. The new GPI data complement recent HST/STIS observations that show the disk extending out to over 500~AU. The GPI image displays a strong asymmetry along the projected minor axis as well as a fainter asymmetry along the projected major axis. We associate the minor and major axis asymmetries with polarized forward scattering and a possible stellocentric offset, respectively. To constrain the disk geometry we fit two separate disk models to the polarized image, each using a different scattering phase function. Both models favor a disk inclination of $\sim 70\degr$ and a $1.5\pm0.6$ AU stellar offset in the plane of the sky along the projected major axis of the disk. We find that the stellar offset in the disk plane, perpendicular to the projected major axis is degenerate with the form of the scattering phase function and remains poorly constrained. The disk is not recovered in total intensity due in part to strong adaptive optics residuals, but we recover three point sources. Considering the system's proximity to the galactic plane and the point sources' positions relative to the disk, we consider it likely that they are background objects and unrelated to the disk's offset from the star.

We present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3$\pm$0.3 au (0.618$\pm$0.004") from the star. With the high astrometric precision afforded by GPI, we have confirmed common proper motion of HR 2562B with the star with only a month time baseline between observations to more than $5\sigma$. Spectral data in $J$, $H$, and $K$ bands show morphological similarity to L/T transition objects. We assign a spectral type of L7$\pm$3 to HR 2562B, and derive a luminosity of $\log$(L$_{\rm bol}$/L$_{\odot}$)=-4.62$\pm$0.12, corresponding to a mass of 30$\pm$15 M$_{\rm Jup}$ from evolutionary models at an estimated age of the system of 300-900 Myr. Although the uncertainty in the age of the host star is significant, the spectra and photometry exhibit several indications of youth for HR 2562B. The source has a position angle consistent with an orbit in the same plane as the debris disk recently resolved with Herschel. Additionally, it appears to be interior to the debris disk. Though the extent of the inner hole is currently too uncertain to place limits on the mass of HR 2562B, future observations of the disk with higher spatial resolution may be able to provide mass constraints. This is the first brown dwarf-mass object found to reside in the inner hole of a debris disk, offering the opportunity to search for evidence of formation above the deuterium burning limit in a circumstellar disk.

A principal scientific goal of the Gemini Planet Imager (GPI) is obtaining milliarcsecond astrometry to constrain exoplanet orbits. However, astrometry of directly imaged exoplanets is subject to biases, systematic errors, and speckle noise. Here we describe an analytical procedure to forward model the signal of an exoplanet that accounts for both the observing strategy (angular and spectral differential imaging) and the data reduction method (Karhunen-Lo\`eve Image Projection algorithm). We use this forward model to measure the position of an exoplanet in a Bayesian framework employing Gaussian processes and Markov chain Monte Carlo (MCMC) to account for correlated noise. In the case of GPI data on $\beta$ Pic b, this technique, which we call Bayesian KLIP-FM Astrometry (BKA), outperforms previous techniques and yields 1$\sigma$-errors at or below the one milliarcsecond level. We validate BKA by fitting a Keplerian orbit to twelve GPI observations along with previous astrometry from other instruments. The statistical properties of the residuals confirm that BKA is accurate and correctly estimates astrometric errors. Our constraints on the orbit of $\beta$ Pic b firmly rule out the possibility of a transit of the planet at 10-$\sigma$ significance. However, we confirm that the Hill sphere of $\beta$ Pic b will transit, giving us a rare chance to probe the circumplanetary environment of a young, evolving exoplanet. We provide an ephemeris for photometric monitoring of the Hill sphere transit event, which will begin at the start of April in 2017 and finish at the end of January in 2018.

Using the Gemini Planet Imager (GPI), we have resolved the circumstellar debris disk around HD 111520 at a projected range of ~30-100 AU in both total and polarized $H$-band intensity. The disk is seen edge-on at a position angle of ~165$^{\circ}$ along the spine of emission. A slight inclination or asymmetric warping are covariant and alters the interpretation of the observed disk emission. We employ 3 point spread function (PSF) subtraction methods to reduce the stellar glare and instrumental artifacts to confirm that there is a roughly 2:1 brightness asymmetry between the NW and SE extension. This specific feature makes HD 111520 the most extreme examples of asymmetric debris disks observed in scattered light among similar highly inclined systems, such as HD 15115 and HD 106906. We further identify a tentative localized brightness enhancement and scale height enhancement associated with the disk at ~40 AU away from the star on the SE extension. We also find that the fractional polarization rises from 10 to 40% from 0.5" to 0.8" from the star. The combination of large brightness asymmetry and symmetric polarization fraction leads us to believe that an azimuthal dust density variation is causing the observed asymmetry.

We present H and K band imaging polarimetry for the PDS 66 circumstellar disk obtained during the commissioning of the Gemini Planet Imager (GPI). Polarization images reveal a clear detection of the disk in to the 0.12'' inner working angle (IWA) in H band, almost 3 times as close to the star as the previous HST observations with NICMOS and STIS (0.35'' effective IWA). The centro-symmetric polarization vectors confirm that the bright inner disk detection is due to circumstellar scattered light. A more diffuse disk extends to a bright outer ring centered at 80 AU. We discuss several physical mechanisms capable of producing the observed ring + gap structure. GPI data confirm enhanced scattering on the East side of the disk which is inferred to be nearer to us. We also detect a lateral asymmetry in the South possibly due to shadowing from material within the inner working angle. This likely corresponds to a temporally variable azimuthal asymmetry observed in HST/STIS coronagraphic imaging.

Detecting polarized light from self-luminous exoplanets has the potential to provide key information about rotation, surface gravity, cloud grain size, and cloud coverage. While field brown dwarfs with detected polarized emission are common, no exoplanet or substellar companion has yet been detected in polarized light. With the advent of high contrast imaging spectro-polarimeters such as GPI and SPHERE, such a detection may now be possible with careful treatment of instrumental polarization. In this paper, we present 28 minutes of $H$-band GPI polarimetric observations of the benchmark T5.5 companion HD 19467 B. We detect no polarization signal from the target, and place an upper limit on the degree of linear polarization of $p_{\text{CL}99.73\%} \leq 2.4\%$. We discuss our results in the context of T dwarf cloud models and photometric variability.

We present the first scattered-light image of the debris disk around HD 131835 in $H$ band using the Gemini Planet Imager. HD 131835 is a $\sim$15 Myr old A2IV star at a distance of $\sim$120 pc in the Sco-Cen OB association. We detect the disk only in polarized light and place an upper limit on the peak total intensity. No point sources resembling exoplanets were identified. Compared to its mid-infrared thermal emission, the disk in scattered light shows similar orientation but different morphology. The scattered-light disk extends from $\sim$75 to $\sim$210 AU in the disk plane with roughly flat surface density. Our Monte Carlo radiative transfer model can well describe the observations with a model disk composed of a mixture of silicates and amorphous carbon. In addition to the obvious brightness asymmetry due to stronger forward scattering, we discover a weak brightness asymmetry along the major axis with the northeast side being 1.3 times brighter than the southwest side at a 3-{\sigma} level.

We present new Gemini Planet Imager observations of the young exoplanet 51 Eridani b which provide further evidence that the companion is physically associated with 51 Eridani. Combining this new astrometric measurement with those reported in the literature, we significantly reduce the posterior probability that 51 Eridani b is an unbound foreground or background T-dwarf in a chance alignment with 51 Eridani to $2\times10^{-7}$, an order of magnitude lower than previously reported. If 51 Eridani b is indeed a bound object, then we have detected orbital motion of the planet between the discovery epoch and the latest epoch. By implementing a computationally efficient Monte Carlo technique, preliminary constraints are placed on the orbital parameters of the system. The current set of astrometric measurements suggest an orbital semimajor axis of $14^{+7}_{-3}$ AU, corresponding to a period of $41^{+35}_{-12}$ years (assuming a mass of $1.75$ $M_{\odot}$ for the central star), and an inclination of $138^{+15}_{-13}$ deg. The remaining orbital elements are only marginally constrained by the current measurements. These preliminary values suggest an orbit which does not share the same inclination as the orbit of the distant M-dwarf binary, GJ 3305, which is a wide physically bound companion to 51 Eridani.

Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter.

We present the first scattered light detections of the HD 106906 debris disk using Gemini/GPI in the infrared and HST/ACS in the optical. HD 106906 is a 13 Myr old F5V star in the Sco-Cen association, with a previously detected planet-mass candidate HD 106906b projected 650 AU from the host star. Our observations reveal a near edge-on debris disk that has a central cleared region with radius $\sim$50 AU, and an outer extent $>$500 AU. The HST data show the outer regions are highly asymmetric, resembling the ''needle'' morphology seen for the HD 15115 debris disk. The planet candidate is oriented $\sim$21$\deg$ away from the position angle of the primary's debris disk, strongly suggesting non-coplanarity with the system. We hypothesize that HD 106906b could be dynamically involved in the perturbation of the primary's disk, and investigate whether or not there is evidence for a circumplanetary dust disk or cloud that is either primordial or captured from the primary. We show that both the existing optical properties and near-infrared colors of HD 106906b are weakly consistent with this possibility, motivating future work to test for the observational signatures of dust surrounding the planet.

We present Gemini Planet Imager (GPI) observations of AU Microscopii, a young M dwarf with an edge-on, dusty debris disk. Integral field spectroscopy and broadband imaging polarimetry were obtained during the commissioning of GPI. In our broadband imaging polarimetry observations, we detect the disk only in total intensity and find asymmetries in the morphology of the disk between the southeast and northwest sides. The southeast side of the disk exhibits a bump at 1$''$ (10 AU projected separation) that is three times more vertically extended and three times fainter in peak surface brightness than the northwest side at similar separations. This part of the disk is also vertically offset by 69$\pm$30 mas to the northeast at 1$''$ when compared to the established disk mid-plane and consistent with prior ALMA and Hubble Space Telescope/STIS observations. We see hints that the southeast bump might be a result of detecting a horizontal sliver feature above the main disk that could be the disk backside. Alternatively when including the morphology of the northwest side, where the disk mid-plane is offset in the opposite direction $\sim$50 mas between 0$.''$4 and 1$.''$2, the asymmetries suggest a warp-like feature. Using our integral field spectroscopy data to search for planets, we are 50% complete for $\sim$4 $M_\mathrm{Jup}$ planets at 4 AU. We detect a source, resolved only along the disk plane, that could either be a candidate planetary mass companion or a compact clump in the disk.

We present $H$-band observations of $\beta$ Pic with the Gemini Planet Imager's (GPI's) polarimetry mode that reveal the debris disk between ~0.3" (~6 AU) and ~1.7" (~33 AU), while simultaneously detecting $\beta$ Pic $b$. The polarized disk image was fit with a dust density model combined with a Henyey-Greenstein scattering phase function. The best fit model indicates a disk inclined to the line of sight ($\phi=85.27{\deg}^{+0.26}_{-0.19}$) with a position angle $\theta_{PA}=30.35{\deg}^{+0.29}_{-0.28}$ (slightly offset from the main outer disk, $\theta_{PA}\approx29{\deg}$), that extends from an inner disk radius of $23.6^{+0.9}_{-0.6}$ AU to well outside GPI's field of view. In addition, we present an updated orbit for $\beta$ Pic $b$ based on new astrometric measurements taken in GPI's spectroscopic mode spanning 14 months. The planet has a semi-major axis of $a=9.2^{+1.5}_{-0.4}$AU, with an eccentricity $e\leq 0.26$. The position angle of the ascending node is $\Omega=31.75{\deg}\pm0.15$, offset from both the outer main disk and the inner disk seen in the GPI image. The orbital fit constrains the stellar mass of $\beta$ Pic to $1.60\pm0.05 M_{\odot}$. Dynamical sculpting by $\beta$ Pic $b$ cannot easily account for the following three aspects of the inferred disk properties: 1) the modeled inner radius of the disk is farther out than expected if caused by $\beta$ Pic b; 2) the mutual inclination of the inner disk and $\beta$ Pic $b$ is $4{\deg}$, when it is expected to be closer to zero; and 3) the aspect ratio of the disk ($h_0 = 0.137^{+0.005}_{-0.006}$) is larger than expected from interactions with $\beta$ Pic $b$ or self-stirring by the disk's parent bodies.

We define an appropriate problem for benchmarking dust emissivity calculations in the context of radiative transfer (RT) simulations, specifically including the emission from stochastically heated dust grains. Our aim is to provide a self-contained guide for implementors of such functionality, and to offer insights in the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. The benchmark problem definition includes the optical and calorimetric material properties, and the grain size distributions, for a typical astronomical dust mixture with silicate, graphite and PAH components; a series of analytically defined radiation fields to which the dust population is to be exposed; and instructions for the desired output. We process this problem using six RT codes participating in this benchmark effort, and compare the results to a reference solution computed with the publicly available dust emission code DustEM. The participating codes implement different heuristics to keep the calculation time at an acceptable level. We study the effects of these mechanisms on the calculated solutions, and report on the level of (dis)agreement between the participating codes. For all but the most extreme input fields, we find agreement within 10% across the important wavelength range from 3 to 1000 micron. We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of stochastically heated dust grains.

Optical UBVRI photometric measurements using the Faulkes Telescope North were taken in early 2011 and combined with 2MASS JHK$_s$ and WISE infrared photometry as well as UCAC4 proper motion data in order to estimate the main parameters of the galactic open cluster NGC 2215 of which large uncertainty exists in the current literature. Fitting a King model we estimate a core radius of 1.12$'\pm$0.04$'$ (0.24$\pm$0.01pc) and a limiting radius of $4.3'\pm$0.5$'$ (0.94$\pm$0.11pc) for the cluster. The results of isochrone fits indicates an age of $log(t)=8.85\pm0.10$ with a distance of $d=790\pm90$pc, a metallicity of $[Fe/H]=-0.40\pm0.10$ dex and a reddening of $E(B-V)=0.26\pm0.04$. A proportion of the work in this study was undertaken by Australian and Canadian upper secondary school students involved in the Space to Grow astronomy education project, and is the first scientific publication to have utilized our star cluster photometry curriculum materials.

We report the discovery of the resolved disk around HD 131835 and present the analysis and modeling of its thermal emission. HD 131835 is a ~15 Myr A2 star in the Scorpius-Centaurus OB association at a distance of 122.7 +16.2 -12.8 parsec. The extended disk has been detected to ~1.5" (200 AU) at 11.7 {\mu}m and 18.3 {\mu}m with T-ReCS on Gemini South. The disk is inclined at an angle of ~75{\deg} with the position angle of ~61{\deg}. The flux of HD 131835 system is 49.3+-7.6 mJy and 84+-45 mJy at 11.7 {\mu}m and 18.3 {\mu}m respectively. A model with three grain populations gives a satisfactory fit to both the spectral energy distribution and the images simultaneously. This best-fit model is composed of a hot continuous power-law disk and two rings. We characterized the grain temperature profile and found that the grains in all three populations are emitting at temperatures higher than blackbodies. In particular, the grains in the continuous disk are unusually warm; even when considering small graphite particles as the composition.

We describe a successful effort to produce a laser comb around 1.55 $\mu$m in the astronomical H band using a method based on a line-referenced, electro-optical-modulation frequency comb. We discuss the experimental setup, laboratory results, and proof of concept demonstrations at the NASA Infrared Telescope Facility (IRTF) and the Keck-II telescope. The laser comb has a demonstrated stability of $<$ 200 kHz, corresponding to a Doppler precision of ~0.3 m/s. This technology, when coupled with a high spectral resolution spectrograph, offers the promise of $<$1 m/s radial velocity precision suitable for the detection of Earth-sized planets in the habitable zones of cool M-type stars.

During the first-light run of the Gemini Planet Imager (GPI) we obtained K-band spectra of exoplanets HR 8799 c and d. Analysis of the spectra indicates that planet d may be warmer than planet c. Comparisons to recent patchy cloud models and previously obtained observations over multiple wavelengths confirm that thick clouds combined with horizontal variation in the cloud cover generally reproduce the planets' spectral energy distributions. When combined with the 3 to 4 um photometric data points, the observations provide strong constraints on the atmospheric methane content for both planets. The data also provide further evidence that future modeling efforts must include cloud opacity, possibly including cloud holes, disequilibrium chemistry, and super-solar metallicity.

Since the early 1990s with the arrival of a variety of new technologies, the capacity for authentic astronomical research at the high school level has skyrocketed. This potential, however, has not realized the bright-eyed hopes and dreams of the early pioneers who expected to revolutionise science education through the use of telescopes and other astronomical instrumentation in the classroom. In this paper, a general history and analysis of these attempts is presented. We define what we classify as an Astronomy Research in the Classroom (ARiC) project and note the major dimensions on which these projects differ before describing the 22 major student research projects active since the early 1990s. This is followed by a discussion of the major issues identified that affected the success of these projects and provide suggestions for similar attempts in the future.

Using the recently installed Gemini Planet Imager (GPI), we have taken the first H-band spectrum of the planetary companion to the nearby young star beta Pictoris. GPI is designed to image and provide low-resolution spectra of Jupiter sized, self-luminous planetary companions around young nearby stars. These observations were taken covering the H-band (1.65 microns). The spectrum has a resolving power of $\sim$ 45 and demonstrates the distinctive triangular shape of a cool substellar object with low surface gravity. Using atmospheric models, we find an effective temperature of $1650 \pm 50$ K and a surface gravity of $\log(g) = 4.0 \pm 0.25$ (cgs units). These values agree well with predictions from planetary evolution models for a gas giant with mass between 10 and 12 $M_{\rm Jup}$ and age between 10 and 20 Myrs.

We present the first results from the polarimetry mode of the Gemini Planet Imager (GPI), which uses a new integral field polarimetry architecture to provide high contrast linear polarimetry with minimal systematic biases between the orthogonal polarizations. We describe the design, data reduction methods, and performance of polarimetry with GPI. Point spread function subtraction via differential polarimetry suppresses unpolarized starlight by a factor of over 100, and provides sensitivity to circumstellar dust reaching the photon noise limit for these observations. In the case of the circumstellar disk around HR 4796A, GPI's advanced adaptive optics system reveals the disk clearly even prior to PSF subtraction. In polarized light, the disk is seen all the way in to its semi-minor axis for the first time. The disk exhibits surprisingly strong asymmetry in polarized intensity, with the west side >9 times brighter than the east side despite the fact that the east side is slightly brighter in total intensity. Based on a synthesis of the total and polarized intensities, we now believe that the west side is closer to us, contrary to most prior interpretations. Forward scattering by relatively large silicate dust particles leads to the strong polarized intensity on the west side, and the ring must be slightly optically thick in order to explain the lower brightness in total intensity there. These findings suggest that the ring is geometrically narrow and dynamically cold, perhaps shepherded by larger bodies in the same manner as Saturn's F ring.

We present the results of both laboratory and on sky astrometric characterization of the Gemini Planet Imager (GPI). This characterization includes measurement of the pixel scale of the integral field spectrograph (IFS), the position of the detector with respect to north, and optical distortion. Two of these three quantities (pixel scale and distortion) were measured in the laboratory using two transparent grids of spots, one with a square pattern and the other with a random pattern. The pixel scale in the laboratory was also estimate using small movements of the artificial star unit (ASU) in the GPI adaptive optics system. On sky, the pixel scale and the north angle are determined using a number of known binary or multiple systems and Solar System objects, a subsample of which had concurrent measurements at Keck Observatory. Our current estimate of the GPI pixel scale is 14.14 $\pm$ 0.01 millarcseconds/pixel, and the north angle is -1.00 $\pm$ 0.03$\deg$. Distortion is shown to be small, with an average positional residual of 0.26 pixels over the field of view, and is corrected using a 5th order polynomial. We also present results from Monte Carlo simulations of the GPI Exoplanet Survey (GPIES) assuming GPI achieves ~1 milliarcsecond relative astrometric precision. We find that with this precision, we will be able to constrain the eccentricities of all detected planets, and possibly determine the underlying eccentricity distribution of widely separated Jovians.

We present on-sky polarimetric observations with the Gemini Planet Imager (GPI) obtained at straight Cassegrain focus on the Gemini South 8-m telescope. Observations of polarimetric calibrator stars, ranging from nearly unpolarized to strongly polarized, enable determination of the combined telescope and instrumental polarization. We find the conversion of Stokes $I$ to linear and circular instrumental polarization in the instrument frame to be $I \rightarrow (Q_{\rm IP}, U_{\rm IP}, P_{\rm IP}, V_{\rm IP}) = (-0.037 \pm 0.010\%, +0.4338 \pm 0.0075\%, 0.4354 \pm 0.0075\%, -6.64 \pm 0.56\%)$. Such precise measurement of instrumental polarization enables $\sim 0.1\%$ absolute accuracy in measurements of linear polarization, which together with GPI's high contrast will allow GPI to explore scattered light from circumstellar disk in unprecedented detail, conduct observations of a range of other astronomical bodies, and potentially even study polarized thermal emission from young exoplanets. Observations of unpolarized standard stars also let us quantify how well GPI's differential polarimetry mode can suppress the stellar PSF halo. We show that GPI polarimetry achieves cancellation of unpolarized starlight by factors of 100-200, reaching the photon noise limit for sensitivity to circumstellar scattered light for all but the smallest separations at which the calibration for instrumental polarization currently sets the limit.

The Gemini Planet Imager (GPI) has as its science instrument an infrared integral field spectrograph/polarimeter (IFS). Integral field spectrographs are scientificially powerful but require sophisticated data reduction systems. For GPI to achieve its scientific goals of exoplanet and disk characterization, IFS data must be reconstructed into high quality astrometrically and photometrically accurate datacubes in both spectral and polarization modes, via flexible software that is usable by the broad Gemini community. The data reduction pipeline developed by the GPI instrument team to meet these needs is now publicly available following GPI's commissioning. This paper, the first of a series, provides a broad overview of GPI data reduction, summarizes key steps, and presents the overall software framework and implementation. Subsequent papers describe in more detail the algorithms necessary for calibrating GPI data. The GPI data reduction pipeline is open source, available from planetimager.org, and will continue to be enhanced throughout the life of the instrument. It implements an extensive suite of task primitives that can be assembled into reduction recipes to produce calibrated datasets ready for scientific analysis. Angular, spectral, and polarimetric differential imaging are supported. Graphical tools automate the production and editing of recipes, an integrated calibration database manages reference files, and an interactive data viewer customized for high contrast imaging allows for exploration and manipulation of data.

The Gemini Planet Imager (GPI) is a complex optical system designed to directly detect the self-emission of young planets within two arcseconds of their host stars. After suppressing the starlight with an advanced AO system and apodized coronagraph, the dominant residual contamination in the focal plane are speckles from the atmosphere and optical surfaces. Since speckles are diffractive in nature their positions in the field are strongly wavelength dependent, while an actual companion planet will remain at fixed separation. By comparing multiple images at different wavelengths taken simultaneously, we can freeze the speckle pattern and extract the planet light adding an order of magnitude of contrast. To achieve a bandpass of 20%, sufficient to perform speckle suppression, and to observe the entire two arcsecond field of view at diffraction limited sampling, we designed and built an integral field spectrograph with extremely low wavefront error and almost no chromatic aberration. The spectrograph is fully cryogenic and operates in the wavelength range 1 to 2.4 microns with five selectable filters. A prism is used to produce a spectral resolution of 45 in the primary detection band and maintain high throughput. Based on the OSIRIS spectrograph at Keck, we selected to use a lenslet-based spectrograph to achieve an rms wavefront error of approximately 25 nm. Over 36,000 spectra are taken simultaneously and reassembled into image cubes that have roughly 192x192 spatial elements and contain between 11 and 20 spectral channels. The primary dispersion prism can be replaced with a Wollaston prism for dual polarization measurements. The spectrograph also has a pupil-viewing mode for alignment and calibration.

The Gemini Planet Imager (GPI) is a dedicated facility for directly imaging and spectroscopically characterizing extrasolar planets. It combines a very high-order adaptive optics system, a diffraction-suppressing coronagraph, and an integral field spectrograph with low spectral resolution but high spatial resolution. Every aspect of GPI has been tuned for maximum sensitivity to faint planets near bright stars. During first light observations, we achieved an estimated H band Strehl ratio of 0.89 and a 5-sigma contrast of $10^6$ at 0.75 arcseconds and $10^5$ at 0.35 arcseconds. Observations of Beta Pictoris clearly detect the planet, Beta Pictoris b, in a single 60-second exposure with minimal post-processing. Beta Pictoris b is observed at a separation of $434 \pm 6$ milli-arcseconds and position angle $211.8 \pm 0.5$ deg. Fitting the Keplerian orbit of Beta Pic b using the new position together with previous astrometry gives a factor of three improvement in most parameters over previous solutions. The planet orbits at a semi-major axis of $9.0^{+0.8}_{-0.4}$ AU near the 3:2 resonance with the previously-known 6 AU asteroidal belt and is aligned with the inner warped disk. The observations give a 4% posterior probability of a transit of the planet in late 2017.

We have conducted a detailed multi-wavelength study of the peculiar nebula Abell 48 and its central star. We classify the nucleus as a helium-rich, hydrogen-deficient star of type [WN4-5]. The evidence for either a massive WN or a low-mass [WN] interpretation is critically examined, and we firmly conclude that Abell 48 is a planetary nebula (PN) around an evolved low-mass star, rather than a Population I ejecta nebula. Importantly, the surrounding nebula has a morphology typical of PNe, and is not enriched in nitrogen, and thus not the `peeled atmosphere' of a massive star. We estimate a distance of 1.6 kpc and a reddening, E(B-V) = 1.90 mag, the latter value clearly showing the nebula lies on the near side of the Galactic bar, and cannot be a massive WN star. The ionized mass (~0.3 M_Sun) and electron density (700 cm^-3) are typical of middle-aged PNe. The observed stellar spectrum was compared to a grid of models from the Potsdam Wolf-Rayet (PoWR) grid. The best fit temperature is 71 kK, and the atmospheric composition is dominated by helium with an upper limit on the hydrogen abundance of 10 per cent. Our results are in very good agreement with the recent study of Todt et al., who determined a hydrogen fraction of 10 per cent and an unusually large nitrogen fraction of ~5 per cent. This fraction is higher than any other low-mass H-deficient star, and is not readily explained by current post-AGB models. We give a discussion of the implications of this discovery for the late-stage evolution of intermediate-mass stars. There is now tentative evidence for two distinct helium-dominated post-AGB lineages, separate to the helium and carbon dominated surface compositions produced by a late thermal pulse. Further theoretical work is needed to explain these recent discoveries.

We present a new technique for forward-modeling self-subtraction of spatially extended emission in observations processed with angular differential imaging (ADI) algorithms. High-contrast direct imaging of circumstellar disks is limited by quasi-static speckle noise and ADI is commonly used to suppress those speckles. However, the application of ADI can result in self-subtraction of the disk signal due to the disk's finite spatial extent. This signal attenuation varies with radial separation and biases measurements of the disk's surface brightness, thereby compromising inferences regarding the physical processes responsible for the dust distribution. To compensate for this attenuation, we forward-model the disk structure and compute the form of the self-subtraction function at each separation. As a proof of concept, we apply our method to 1.6 and 2.2 micron Keck AO NIRC2 scattered-light observations of the HD 32297 debris disk reduced using a variant of the "locally optimized combination of images" (LOCI) algorithm. We are able to recover disk surface brightness that was otherwise lost to self-subtraction and produce simplified models of the brightness distribution as it appears with and without self-subtraction. From the latter models, we extract radial profiles for the disk's brightness, width, midplane position, and color that are unbiased by self-subtraction. Our analysis of these measurements indicates a break in the brightness profile power law at r~110 AU and a disk width that increases with separation from the star. We also verify disk curvature that displaces the midplane by up to 30 AU towards the northwest relative to a straight fiducial midplane.

We present new optical coronagraphic data of the bright star Fomalhaut obtained with the HST in 2010/2012 using STIS. Fomalhaut b is recovered at both epochs to high significance. The observations include the discoveries of tenuous nebulosity beyond the main dust belt detected to at least 209 AU projected radius and a ~50 AU wide azimuthal gap in the belt northward of Fom b. The morphology of Fomalhaut b appears elliptical in the STIS detections. We show that residual noise in the processed data can plausibly result in point sources appearing extended. A MCMC analysis demonstrates that the orbit of Fom b is highly eccentric, with e=0.8+/-0.1, a=177+/-68 AU, and q = 32+/-24 AU. Fom b is apsidally aligned with the belt and 90% of allowed orbits have mutual inclination 36 deg or less. Fomalhaut b's orbit is belt-crossing in projection, but only 12% of possible orbits have nodes within a 25 AU wide belt annulus (133-158 AU). The high e invokes a dynamical history where Fom b may have experienced a significant dynamical interaction with a hypothetical planet Fomalhaut c, and the current orbital configuration may be relatively short-lived. The new value for the periastron distance diminishes the Hill radius of Fom b and any weakly bound satellite system surrounding a planet would be sheared and dynamically heated at periapse. We argue that Fom b's minimum mass is that of a dwarf planet in order for a circumplanetary satellite system to remain bound to a sufficient radius from the planet to be consistent with the dust scattered light hypothesis. Fom b may be optically bright because the recent passage through periapse and/or the ascending node has increased the erosion rates of planetary satellites. In the coplanar case, Fomalhaut b will collide with the main belt around 2032, and the subsequent emergent phenomena may help determine its physical nature.

We present a calibration of the metallicity and physical parameters (temperature, luminosity, gravity, mass, radius) for RR Lyrae stars using the ugriz SDSS photometric system. Our work is based on calculations of synthetic horizontal branches (HBs), fully taking into account evolutionary effects for a wide range in metallicities and HB morphologies. We provide analytical fits that are able to provide all quantities mentioned with very high (internal) precision, based solely on mean SDSS magnitudes and colors.

We present resolved Herschel images of a circumbinary debris disk in the 99 Herculis system. The primary is a late F-type star. The binary orbit is well characterised and we conclude that the disk is misaligned with the binary plane. Two different models can explain the observed structure. The first model is a ring of polar orbits that move in a plane perpendicular to the binary pericenter direction. We favour this interpretation because it includes the effect of secular perturbations and the disk can survive for Gyr timescales. The second model is a misaligned ring. Because there is an ambiguity in the orientation of the ring, which could be reflected in the sky plane, this ring either has near-polar orbits similar to the first model, or has a 30 degree misalignment. The misaligned ring, interpreted as the result of a recent collision, is shown to be implausible from constraints on the collisional and dynamical evolution. Because disk+star systems with separations similar to 99 Herculis should form coplanar, possible formation scenarios involve either a close stellar encounter or binary exchange in the presence of circumstellar and/or circumbinary disks. Discovery and characterisation of systems like 99 Herculis will help understand processes that result in planetary system misalignment around both single and multiple stars.

We present new imaging data and archival multiwavelength observations of the little studied emission nebula K 1-6 and its central star. Narrow-band images in H-alpha (+ [NII]) and [OIII] taken with the Faulkes Telescope North reveal a stratified, asymmetric, elliptical nebula surrounding a central star which has the colours of a late G- or early K-type subgiant or giant. GALEX ultraviolet images reveal a very hot subdwarf or white dwarf coincident in position with this star. The cooler, optically dominant star is strongly variable with a period of 21.312 +/- 0.008 days, and is possibly a high amplitude member of the RS CVn class, although an FK Com classification is also possible. Archival ROSAT data provide good evidence that the cool star has an active corona. We conclude that K 1-6 is most likely an old bona fide planetary nebula at a distance of ~1.0 kpc, interacting with the interstellar medium, and containing a binary or ternary central star. The observations and data analyses reported in this paper were conducted in conjunction with Year 11 high school students as part of an Australian Research Council Linkage Grant science education project, denoted Space To Grow, conducted jointly by professional astronomers, educational researchers, teachers, and high-school students.

We have observed Eta Carinae over 34 nights between 4th January 2009 and 27th March 2009 covering the estimated timeframe for a predicted spectroscopic event related to a suspected binary system concealed within the homunculus nebula. A photometric minimum feature was confirmed to be periodic and comparison to a previous event indicated that the period to within our error at 2022.6 +/-1.0 d. Using the E-region standard star system, the apparent V magnitudes determined for the local comparison stars were HD303308 8.14+/-0.02, HD 93205 7.77 +/-0.03 and HD93162 8.22 +/-0.05. The latter star was found to be dimmer than previously reported.

We present Hubble Space Telescope optical coronagraphic polarization imaging observations of the dusty debris disk HD 61005. The scattered light intensity image and polarization structure reveal a highly inclined disk with a clear asymmetric, swept back component, suggestive of significant interaction with the ambient interstellar medium. The combination of our new data with the published 1.1 micron discovery image shows that the grains are blue scattering with no strong color gradient as a function of radius, implying predominantly sub-micron sized grains. We investigate possible explanations that could account for the observed swept back, asymmetric morphology. Previous work has suggested that HD 61005 may be interacting with a cold, unusually dense interstellar cloud. However, limits on the intervening interstellar gas column density from an optical spectrum of HD 61005 in the Na I D lines render this possibility unlikely. Instead, HD 61005 may be embedded in a more typical warm, low-density cloud that introduces secular perturbations to dust grain orbits. This mechanism can significantly distort the ensemble disk structure within a typical cloud crossing time. For a counterintuitive relative flow direction--parallel to the disk midplane--we find that the structures generated by these distortions can very roughly approximate the HD 61005 morphology. Future observational studies constraining the direction of the relative interstellar medium flow will thus provide an important constraint for future modeling. Independent of the interpretation for HD 61005, we expect that interstellar gas drag likely plays a role in producing asymmetries observed in other debris disk systems, such as HD 15115 and Delta-Velorum.

A point source observed 8 AU in projection from beta Pictoris in L' (3.8 micron) imaging in 2003 has been recently presented as a planet candidate. Here we show the results of L'-band adaptive optics imaging obtained at Keck Observatory in 2008. We do not detect beta Pic b beyond a limiting radius of 0.29 arcsec, or 5.5 AU in projection, from the star. If beta Pic b is an orbiting planet, then it has moved >=0.12 arcsec (2.4 AU in projection) closer to the star in the five years separating the two epochs of observation. We examine the range of orbital parameters consistent with the observations, including likely bounds from the locations of previously inferred planetesimal belts. We find a family of low-eccentricity orbits with semimajor axes ~8-9 AU that are completely allowed, as well as a broad region of orbits with e<~0.2, a>~10 AU that are allowed if the apparent motion of the planet was towards the star in 2003. We compare this allowed space with predictions of the planetary orbital elements from the literature. Additionally, we show how similar observations in the next several years can further constrain the space of allowed orbits. Non-detections of the source through 2013 will exclude the interpretation of the candidate as a planet orbiting between the 6.4 and 16 AU planetesimal belts.

Fomalhaut is a bright star 7.7 parsecs (25 light years) from Earth that harbors a belt of cold dust with a structure consistent with gravitational sculpting by an orbiting planet. Here, we present optical observations of an exoplanet candidate, Fomalhaut b. In the plane of the belt, Fomalhaut b lies approximately 119 astronomical units (AU) from the star and 18 AU from the belt, matching predictions. We detect counterclockwise orbital motion using Hubble Space Telescope observations separated by 1.73 years. Dynamical models of the interaction between the planet and the belt indicate that the planet's mass is at most three times that of Jupiter for the belt to avoid gravitational disruption. The flux detected at 800 nm is also consistent with that of a planet with mass no greater than a few times that of Jupiter. The brightness at 600 nm and the lack of detection at longer wavelengths suggest that the detected flux may include starlight reflected off a circumplanetary disk, with dimension comparable to the orbits of the Galilean satellites. We also observed variability of unknown origin at 600 nm.

We present the first detection and mapping of the HD 32297 debris disk at 1.3 mm with the Combined Array for Research in Millimeter-wave Astronomy (CARMA). With a sub-arcsecond beam, this detection represents the highest angular resolution (sub)mm debris disk observation made to date. Our model fits to the spectral energy distribution from the CARMA flux and new Spitzer MIPS photometry support the earlier suggestion that at least two, possibly three, distinct grain populations are traced by the current data. The observed millimeter map shows an asymmetry between the northeast and southwest disk lobes, suggesting large grains may be trapped in resonance with an unseen exoplanet. Alternatively, the observed morphology could result from the recent breakup of a massive planetesimal. A similar-scale asymmetry is also observed in scattered light but not in the mid-infrared. This contrast between asymmetry at short and long wavelengths and symmetry at intermediate wavelengths is in qualitative agreement with predictions of resonant debris disk models. With resolved observations in several bands spanning over three decades in wavelength, HD 32297 provides a unique testbed for theories of grain and planetary dynamics, and could potentially provide strong multi-wavelength evidence for an exoplanetary system.

We have obtained a Gemini South T-ReCS Qa-band (18.3 micron) image and a Spitzer MIPS SED-mode observation of HD181327, an F5/F6V member of the ~12 Myr old beta Pictoris moving group. We resolve the disk in thermal-emission for the first time and find that the northern arm of the disk is 1.4x brighter than the southern arm. In addition, we detect a broad peak in the combined Spitzer IRS and MIPS spectra at 60 - 75 micron that may be produced by emission from crystalline water ice. We model the IRS and MIPS data using a size distribution of amorphous olivine and water ice grains (dn/da proportional to a^{-2.25} with a_{min} consistent with the minimum blow out size and a_{max} = 20 micron) located at a distance of 86.3 AU from the central star, as observed in previously published scattered-light images. Since the photo-desorption lifetime for the icy particles is ~1400 yr, significantly less than the estimated ~12 Myr age of the system, we hypothesize that we have detected debris that may be steadily replenished by collisions among icy Kuiper belt object-like parent bodies in a newly forming planetary system.

Using the Advanced Camera for Surveys aboard the Hubble Space Telescope, we have discovered dust-scattered light from the debris disk surrounding the F2V star HD 15745. The circumstellar disk is detected between 2.0" and 7.5" radius, corresponding to 128 - 480 AU radius. The circumstellar disk morphology is asymmetric about the star, resembling a fan, and consistent with forward scattering grains in an optically thin disk with an inclination of ~67 degrees to our line of sight. The spectral energy distribution and scattered light morphology can be approximated with a model disk composed of silicate grains between 60 and 450 AU radius, with a total dust mass of 10E-7 M_sun (0.03 M_earth) representing a narrow grain size distribution (1 - 10 micron). Galactic space motions are similar to the Castor Moving Group with an age of ~10E+8 yr, although future work is required to determine the age of HD 15745 using other indicators.

We report the detection of a ring of warm dust in the edge-on disk surrounding HD 32297 with the Gemini-N/MICHELLE mid-infrared imager. Our N'-band image shows elongated structure consistent with the orientation of the scattered-light disk. The Fnu(11.2 um) = 49.9+/-2.1 mJy flux is significantly above the 28.2+/-0.6 mJy photosphere. Subtraction of the stellar point spread function reveals a bilobed structure with peaks 0.5"-0.6" from the star. An analysis of the stellar component of the SED suggests a spectral type later than A0, in contrast to commonly cited literature values. We fit three-dimensional, single-size grain models of an optically thin dust ring to our image and the SED using a Markov chain Monte Carlo algorithm in a Bayesian framework. The best-fit effective grain sizes are submicron, suggesting the same dust population is responsible for the bulk of the scattered light. The inner boundary of the warm dust is located 0.5"-0.7" (~65 AU) from the star, which is approximately cospatial with the outer boundary of the scattered-light asymmetry inward of 0.5". The addition of a separate component of larger, cooler grains that provide a portion of the 60 um flux improves both the fidelity of the model fit and consistency with the slopes of the scattered-light brightness profiles. Previous indirect estimates of the stellar age (~30 Myr) indicate the dust is composed of debris. The peak vertical optical depths in our models (~0.3-1 x 1e-2) imply that grain-grain collisions likely play a significant role in dust dynamics and evolution. Submicron grains can survive radiation pressure blow-out if they are icy and porous. Similarly, the inferred warm temperatures (130-200 K) suggest that ice sublimation may play a role in truncating the inner disk.

(abridged) Debris disks around main sequence stars are produced by the erosion and evaporation of unseen parent bodies. AU Microscopii (GJ 803) is a compelling object to study in the context of disk evolution across different spectral types, as it is an M dwarf whose near edge-on disk may be directly compared to that of its A5V sibling beta Pic. We resolve the disk from 8-60 AU in the near-IR JHK' bands at high resolution with the Keck II telescope and adaptive optics, and develop a novel data reduction technique for the removal of the stellar point spread function. The point source detection sensitivity in the disk midplane is more than a magnitude less sensitive than regions away from the disk for some radii. We measure a blue color across the near-IR bands, and confirm the presence of substructure in the inner disk. Some of the structural features exhibit wavelength-dependent positions. The disk architecture and characteristics of grain composition are inferred through modeling. We approach the modeling of the dust distribution in a manner that complements previous work. Using a Monte Carlo radiative transfer code, we compare a relatively simple model of the distribution of porous grains to a broad data set, simultaneously fitting to midplane surface brightness profiles and the spectral energy distribution. Our model confirms that the large-scale architecture of the disk is consistent with detailed models of steady-state grain dynamics. Here, a belt of parent bodies from 35-40 AU is responsible for producing dust that is then swept outward by the stellar wind and radiation pressures. We infer the presence of very small grains in the outer region, down to sizes of ~0.05 micron. These sizes are consistent with stellar mass-loss rates Mdot_* << 10^2 Mdot_sun.

We report the first scattered light detection of a dusty debris disk surrounding the F2V star HD 15115 using the Hubble Space Telescope in the optical, and Keck adaptive optics in the near-infrared. The most remarkable property of the HD 15115 disk relative to other debris disks is its extreme length asymmetry. The east side of the disk is detected to ~315 AU radius, whereas the west side of the disk has radius >550 AU. We find a blue optical to near-infrared scattered light color relative to the star that indicates grain scattering properties similar to the AU Mic debris disk. The existence of a large debris disk surrounding HD 15115 adds further evidence for membership in the Beta Pic moving group, which was previously argued based on kinematics alone. Here we hypothesize that the extreme disk asymmetry is due to dynamical perturbations from HIP 12545, an M star 0.5 degrees (0.38 pc) east of HD 15115 that shares a common proper motion vector, heliocentric distance, galactic space velocity, and age.

We present the first scattered light images of debris disks around a K star (HD 53143) and an F star (HD 139664) using the coronagraphic mode of the Advanced Camera for Surveys (ACS) aboard the Hubble Space Telescope (HST). With ages 0.3 - 1 Gyr, these are among the oldest optically detected debris disks. HD 53143, viewed ~45 degrees from edge-on, does not show radial variation in disk structure and has width >55 AU. HD 139664 is seen close to edge-on and has belt-like morphology with a dust peak 83 AU from the star and a distinct outer boundary at 109 AU. We discuss evidence for significant diversity in the radial architecture of debris disks that appears unconnected to stellar spectral type or age. HD 139664 and possibly the solar system belong in a category of narrow belts 20-30 AU wide. HD 53143 represents a class of wide disk architecture with characteristic width >50 AU.

(abridged) Imaging observations are generally affected by a fluctuating background of speckles, a particular problem when detecting faint stellar companions at small angular separations. Knowing the distribution of the speckle intensities at a given location in the image plane is important for understanding the noise limits of companion detection. The speckle noise limit in a long-exposure image is characterized by the intensity variance and the speckle lifetime. In this paper we address the former quantity through the distribution function of speckle intensity. Previous theoretical work has predicted a form for this distribution function at a single location in the image plane. We developed a fast readout mode to take short exposures of stellar images corrected by adaptive optics at the ground-based UCO/Lick Observatory, with integration times of 5 ms and a time between successive frames of 14.5 ms ($\lambda=2.2$ $\mu$m). These observations temporally oversample and spatially Nyquist sample the observed speckle patterns. We show, for various locations in the image plane, the observed distribution of speckle intensities is consistent with the predicted form. Additionally, we demonstrate a method by which $I_c$ and $I_s$ can be mapped over the image plane. As the quantity $I_c$ is proportional to the PSF of the telescope free of random atmospheric aberrations, this method can be used for PSF calibration and reconstruction.