• ### GPI spectra of HR 8799 c, d, and e from 1.5 to 2.4$\mu$m with KLIP Forward Modeling(1804.07774)

April 20, 2018 astro-ph.EP, astro-ph.IM
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.
• ### GPI Spectroscopy of the Mass, Age, and Metallicity Benchmark Brown Dwarf HD 4747 B(1801.06192)

Jan. 18, 2018 astro-ph.SR, astro-ph.EP
The physical properties of brown dwarf companions found to orbit nearby, solar-type stars can be benchmarked against independent measures of their mass, age, chemical composition, and other parameters, offering insights into the evolution of substellar objects. The TRENDS high-contrast imaging survey has recently discovered a (mass/age/metallicity) benchmark brown dwarf orbiting the nearby (d=18.69+/-0.19 pc), G8V/K0V star HD 4747. We have acquired follow-up spectroscopic measurements of HD 4747 B using the Gemini Planet Imager to study its spectral type, effective temperature, surface gravity, and cloud properties. Observations obtained in the H-band and K1-band recover the companion and reveal that it is near the L/T transition (T1+/-2). Fitting atmospheric models to the companion spectrum, we find strong evidence for the presence of clouds. However, spectral models cannot satisfactorily fit the complete data set: while the shape of the spectrum can be well-matched in individual filters, a joint fit across the full passband results in discrepancies that are a consequence of the inherent color of the brown dwarf. We also find a $2\sigma$ tension in the companion mass, age, and surface gravity when comparing to evolutionary models. These results highlight the importance of using benchmark objects to study "secondary effects" such as metallicity, non-equilibrium chemistry, cloud parameters, electron conduction, non-adiabatic cooling, and other subtleties affecting emergent spectra. As a new L/T transition benchmark, HD 4747 B warrants further investigation into the modeling of cloud physics using higher resolution spectroscopy across a broader range of wavelengths, polarimetric observations, and continued Doppler radial velocity and astrometric monitoring.
• ### Automated data processing architecture for the Gemini Planet Imager Exoplanet Survey(1801.01902)

Jan. 5, 2018 astro-ph.EP, astro-ph.IM
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.
• ### Type Ia Supernova Distances at z > 1.5 from the Hubble Space Telescope Multi-Cycle Treasury Programs: The Early Expansion Rate(1710.00844)

Oct. 2, 2017 astro-ph.CO
We present an analysis of 15 Type Ia supernovae (SNe Ia) at redshift z > 1 (9 at 1.5 < z < 2.3) recently discovered in the CANDELS and CLASH Multi-Cycle Treasury programs using WFC3 on the Hubble Space Telescope. We combine these SNe Ia with a new compilation of 1050 SNe Ia, jointly calibrated and corrected for simulated survey biases to produce accurate distance measurements. We present unbiased constraints on the expansion rate at six redshifts in the range 0.07 < z < 1.5 based only on this combined SN Ia sample. The added leverage of our new sample at z > 1.5 leads to a factor of ~3 improvement in the determination of the expansion rate at z = 1.5, reducing its uncertainty to ~20%, a measurement of H(z=1.5)/H0=2.67 (+0.83,-0.52). We then demonstrate that these six measurements alone provide a nearly identical characterization of dark energy as the full SN sample, making them an efficient compression of the SN Ia data. The new sample of SNe Ia at z > 1 usefully distinguishes between alternative cosmological models and unmodeled evolution of the SN Ia distance indicators, placing empirical limits on the latter. Finally, employing a realistic simulation of a potential WFIRST SN survey observing strategy, we forecast optimistic future constraints on the expansion rate from SNe Ia.
• ### Evidence that the Directly-Imaged Planet HD 131399 Ab is a Background Star(1705.06851)

Sept. 7, 2017 astro-ph.SR, astro-ph.EP
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.
• ### Improving and Assessing Planet Sensitivity of the GPI Exoplanet Survey with a Forward Model Matched Filter(1705.05477)

May 15, 2017 astro-ph.EP, astro-ph.IM
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}$.
• ### Complex Spiral Structure in the HD 100546 Transitional Disk as Revealed by GPI and MagAO(1704.06260)

April 19, 2017 astro-ph.EP
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.
• ### Bringing "The Moth" to Light: A Planet-Sculpting Scenario for the HD 61005 Debris Disk(1605.06161)

Sept. 16, 2016 astro-ph.EP
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.
• ### Imaging an 80 AU Radius Dust Ring Around the F5V Star HD 157587(1609.00382)

Sept. 6, 2016 astro-ph.EP
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.
• ### First images of debris disks around TWA 7, TWA 25, HD 35650, and HD 377(1512.02220)

Dec. 7, 2015 astro-ph.SR, astro-ph.EP
We present the first images of four debris disks observed in scattered light around the young (4--250 Myr old) M dwarfs TWA 7 and TWA 25, the K6 star HD 35650, and the G2 star HD 377. We obtained these images by reprocessing archival Hubble Space Telescope NICMOS coronagraph data with modern post-processing techniques as part of the Archival Legacy Investigation of Circumstellar Environments (ALICE) program. All four disks appear faint and compact compared with other debris disks resolved in scattered light. The disks around TWA 25, HD 35650, and HD 377 appear very inclined, while TWA 7's disk is viewed nearly face-on. The surface brightness of HD 35650's disk is strongly asymmetric. These new detections raise the number of disks resolved in scattered light around M and late-K stars from one (the AU Mic system) to four. This new sample of resolved disks enables comparative studies of heretofore scarce debris disks around low-mass stars relative to solar-type stars.
• ### Discovery of a Companion Candidate in the HD169142 Transition Disk and the Possibility of Multiple Planet Formation(1408.0813)

Aug. 4, 2014 astro-ph.EP
We present L' and J-band high-contrast observations of HD169142, obtained with the VLT/NACO AGPM vector vortex coronagraph and the Gemini Planet Imager, respectively. A source located at 0".156+/-0".032 north of the host star (PA=7.4+/-11.3 degrees) appears in the final reduced L' image. At the distance of the star (~145 pc), this angular separation corresponds to a physical separation of 22.7+/-4.7 AU, locating the source within the recently resolved inner cavity of the transition disk. The source has a brightness of L'=12.2+/-0.5 mag, whereas it is not detected in the J band (J>13.8 mag). If its L' brightness arose solely from the photosphere of a companion and given the J-L' color constraints, it would correspond to a 28-32 MJupiter object at the age of the star, according to the COND models. Ongoing accretion activity of the star suggests, however, that gas is left in the inner disk cavity from which the companion could also be accreting. In this case the object could be lower in mass and its luminosity enhanced by the accretion process and by a circumplanetary disk. A lower mass object is more consistent with the observed cavity width. Finally, the observations enable us to place an upper limit on the L'-band flux of a second companion candidate orbiting in the disk annular gap at ~50 AU, as suggested by millimeter observations. If the second companion is also confirmed, HD169142 might be forming a planetary system, with at least two companions opening gaps and possibly interacting with each other.
• ### The First H-band Spectrum of the Massive Gas Giant Planet beta Pictoris b with the Gemini Planet Imager(1407.4469)

July 16, 2014 astro-ph.EP
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.
• ### Gemini Planet Imager Observational Calibrations IX: Least-Squares Inversion Flux Extraction(1407.2309)

July 9, 2014 astro-ph.EP, astro-ph.IM
The Gemini Planet Imager (GPI) is an instrument designed to directly image planets and circumstellar disks from 0.9 to 2.5 microns (the $YJHK$ infrared bands) using high contrast adaptive optics with a lenslet-based integral field spectrograph. We develop an extraction algorithm based on a least-squares method to disentangle the spectra and systematic noise contributions simultaneously. We utilize two approaches to adjust for the effect of flexure of the GPI optics which move the position of light incident on the detector. The first method is to iterate the extraction to achieve minimum residual and the second is to cross-correlate the detector image with a model image in iterative extraction steps to determine an offset. Thus far, this process has made clear qualitative improvements to the cube extraction by reducing the Moir\'{e} pattern. There are also improvements to the automated routines for finding flexure offsets which are reliable to with $\sim0.5$ pixel accuracy compared to pixel accuracy prior. Further testing and optimization will follow before implementation into the GPI pipeline.
• ### Gemini Planet imager Observational Calibrations X: Non-Redundant Masking on GPI(1407.2310)

July 9, 2014 astro-ph.EP, astro-ph.IM
The Gemini Planet Imager (GPI) Extreme Adaptive Optics Coronograph contains an interferometric mode: a 10-hole non-redundant mask (NRM) in its pupil wheel. GPI operates at $Y, J, H$, and $K$ bands, using an integral field unit spectrograph (IFS) to obtain spectral data at every image pixel. NRM on GPI is capable of imaging with a half resolution element inner working angle at moderate contrast, probing the region behind the coronagraphic spot. The fine features of the NRM PSF can provide a reliable check on the plate scale, while also acting as an attenuator for spectral standard calibrators that would otherwise saturate the full pupil. NRM commissioning data provides details about wavefront error in the optics as well as operations of adaptive optics control without pointing control from the calibration system. We compare lab and on-sky results to evaluate systematic instrument properties and examine the stability data in consecutive exposures. We discuss early on-sky performance, comparing images from integration and tests with the first on-sky images, and demonstrate resolving a known binary. We discuss the status of NRM and implications for future science with this mode.
• ### Type Ia Supernova Rate Measurements to Redshift 2.5 from CANDELS : Searching for Prompt Explosions in the Early Universe(1401.7978)

April 25, 2014 astro-ph.CO
The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) was a multi-cycle treasury program on the Hubble Space Telescope (HST) that surveyed a total area of ~0.25 deg^2 with ~900 HST orbits spread across 5 fields over 3 years. Within these survey images we discovered 65 supernovae (SN) of all types, out to z~2.5. We classify ~24 of these as Type Ia SN (SN Ia) based on host-galaxy redshifts and SN photometry (supplemented by grism spectroscopy of 6 SN). Here we present a measurement of the volumetric SN Ia rate as a function of redshift, reaching for the first time beyond z=2 and putting new constraints on SN Ia progenitor models. Our highest redshift bin includes detections of SN that exploded when the universe was only ~3 Gyr old and near the peak of the cosmic star-formation history. This gives the CANDELS high-redshift sample unique leverage for evaluating the fraction of SN Ia that explode promptly after formation (<500 Myr). Combining the CANDELS rates with all available SN Ia rate measurements in the literature we find that this prompt SN Ia fraction is fP=0.53 +0.09 -0.10 (stat) +0.10 -0.26 (sys), consistent with a delay time distribution that follows a simple t^{-1} power law for all times t>40 Myr. However, a mild tension is apparent between ground-based low-z surveys and space-based high-z surveys. In both CANDELS and the sister HST program CLASH, we find a low rate of SN Ia at z>1. This could be a hint that prompt progenitors are in fact relatively rare, accounting for only ~20% of all SN Ia explosions -- though further analysis and larger samples will be needed to examine that suggestion.
• ### Five Debris Disks Newly Revealed in Scattered Light from the HST NICMOS Archive(1404.5614)

April 22, 2014 astro-ph.SR
We have spatially resolved five debris disks (HD 30447, HD 35841, HD 141943, HD 191089, and HD 202917) for the first time in near-infrared scattered light by reanalyzing archival Hubble Space Telescope (HST)/NICMOS coronagraphic images obtained between 1999 and 2006. One of these disks (HD 202917) was previously resolved at visible wavelengths using HST/Advanced Camera for Surveys. To obtain these new disk images, we performed advanced point-spread function subtraction based on the Karhunen-Loeve Image Projection (KLIP) algorithm on recently reprocessed NICMOS data with improved detector artifact removal (Legacy Archive PSF Library And Circumstellar Environments Legacy program). Three of the disks (HD 30447, HD 35841, and HD 141943) appear edge-on, while the other two (HD 191089 and HD 202917) appear inclined. The inclined disks have been sculpted into rings; in particular, the disk around HD 202917 exhibits strong asymmetries. All five host stars are young (8-40 Myr), nearby (40-100 pc) F and G stars, and one (HD 141943) is a close analog to the young sun during the epoch of terrestrial planet formation. Our discoveries increase the number of debris disks resolved in scattered light from 19 to 23 (a 21% increase). Given their youth, proximity, and brightness (V = 7.2 to 8.5), these targets are excellent candidates for follow-up investigations of planet formation at visible wavelengths using the HST/STIS coronagraph, at near-infrared wavelengths with the Gemini Planet Imager (GPI) and Very Large Telescope (VLT)/SPHERE, and at thermal infrared wavelengths with the James Webb Space Telescope NIRCam and MIRI coronagraphs.
• ### Neptune on tiptoes: dynamical histories that preserve the cold classical Kuiper belt(1112.1954)

Feb. 27, 2012 astro-ph.EP
The current dynamical structure of the Kuiper belt was shaped by the orbital evolution of the giant planets, especially Neptune, during the era following planet formation, when the giant planets may have undergone planet-planet scattering and/or planetesimal-driven migration. Numerical simulations of this process, while reproducing many properties of the belt, fail to generate the high inclinations and eccentricities observed for some objects while maintaining the observed dynamically "cold" population. We present the first of a three-part parameter study of how different dynamical histories of Neptune sculpt the planetesimal disk. Here we identify which dynamical histories allow an in situ planetesimal disk to remain dynamically cold, becoming today's cold Kuiper belt population. We find that if Neptune undergoes a period of elevated eccentricity and/or inclination, it secularly excites the eccentricities and inclinations of the planetesimal disk. We demonstrate that there are several well-defined regimes for this secular excitation, depending on the relative timescales of Neptune's migration, the damping of Neptune's orbital inclination and/or eccentricity, and the secular evolution of the planetesimals. We model this secular excitation analytically in each regime, allowing for a thorough exploration of parameter space. Neptune's eccentricity and inclination can remain high for a limited amount of time without disrupting the cold classical belt. In the regime of slow damping and slow migration, if Neptune is located (for example) at 20 AU, then its eccentricity must stay below 0.18 and its inclination below 6{\deg}.