Search

OSIRIS is a near-infrared (1.0-2.4 μm) integral field spectrograph operating behind the adaptive optics system at Keck Observatory and 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 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, bad-pixel mask, and 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.

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey. 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, semimajor axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M ∗ >1.5 M o more likely to host planets with masses between 2 and 13M Jup and semimajor axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semimajor axis (a) for planet populations around high-mass stars (M ∗ >1.5 M o) of the form , finding α = -2.4 +0.8 and β = -2.0 +0.5, and an integrated occurrence rate of % between 5-13M Jup and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with % of stars hosting a brown dwarf companion between 13-80M Jup and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semimajor axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semimajor axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the radial velocity method, our results are consistent with a peak in occurrence of giant planets between ∼1 and 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.

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 μm), we detect the highly inclined (i = 85°) 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°, 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 μm 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 the H and 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 data set. 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 and K spectra, but do not find any conclusive differences between d and e, nor between 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, as well as how differences between the properties of these planets can be further explored.

The Gemini Planet Imager Exoplanet Survey (GPIES) is a multiyear 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 (DRPs) 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 DRPs. 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.

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.

NIRSPEC is a 1-5 um echelle spectrograph in use on the Keck II Telescope since 1999. The spectrograph is capable of both moderate (R=λ/Δλ∼2000) 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 μm slit-viewing camera, detector characterization of the science and Engineering A grade arrays, electronics systems, and updated software design. The optomechanical design of the slit-viewing camera and replacement detector head assembly have both been assembled and cold-tested in our lab. We also show results from the GigE interface to the SAM/ASIC boards to control the H2RG. The upgrade will continue NIRSPEC's legacy as a powerful near-infrared spectrograph behind one of the world's most scientifically productive telescopes.

We present evidence that the recently discovered, directly imaged planet HD 131399 Ab is a background star with nonzero 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 reinvestigate 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-1. 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 spectrophotometry spanning 1-5 μm of 51 Eridani b, a 2-10 planet discovered by the Gemini Planet Imager Exoplanet Survey. In this study, we present new K1 (1.90-2.19 μm) and K2 (2.10-2.40 μm) spectra taken with the Gemini Planet Imager as well as an updated L P (3.76 μm) and new M S (4.67 μm) photometry from the NIRC2 Narrow camera. The new data were combined with J (1.13-1.35 μm) and H (1.50-1.80 μm) spectra from the discovery epoch with the goal of better characterizing the planet properties. The 51 Eri b photometry is redder than field brown dwarfs as well as known young T-dwarfs with similar spectral type (between T4 and 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. The model fits suggest that 51 Eri b has an effective temperature ranging between 605 and 737 K, a solar metallicity, and a surface gravity of log(g) = 3.5-4.0 dex, and the atmosphere requires a patchy cloud atmosphere to model the spectral energy distribution (SED). From the model atmospheres, we infer a luminosity for the planet of -5.83 to -5.93 (logL/L⊙), leaving 51 Eri b in the unique position of being one of the only directly imaged planets consistent with having formed via a cold-start scenario. Comparisons of the planet SED against warm-start models indicate that the planet luminosity is best reproduced by a planet formed via core accretion with a core mass between 15 and 127 M⊕.

Recent Atacama Large Millimeter/submillimeter Array observations present mounting evidence for the presence of exocometary gas released within Kuiper Belt analogs 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 colocated with dust emission from the cometary belt in the 440 Myr old Fomalhaut system. Through spectrospatial filtering, we achieve a 5.4σ detection and determine that the ring's sky-projected rotation axis matches that of the star. The CO mass derived ((0.65-42) × 10-7 M⊕) 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+CO2 mass fraction of exocomets around Fomalhaut to be between 4.6% and 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 interstellar medium inheritance. In addition, we find tentative evidence that (49 ± 27)% of the detected flux originates from a region near the eccentric belt's pericenter. If confirmed, the latter may be explained through a recent impact event or CO pericenter 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éve 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 (S/N) of the recovered planets. We show that only two spectral templates are necessary to recover any young Jovian exoplanets with minimal S/N loss. We also developed a complete pipeline for the automated detection of point-source candidates, the calculation of receiver operating characteristics (ROC), contrast curves based on false positives, and completeness contours. We process in a uniform manner more than 330 data sets from the Gemini Planet Imager Exoplanet Survey 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 ALMA mosaic observations at 1.3 mm (223 GHz) of the Fomalhaut system with a sensitivity of 14 μ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 an 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 ± 0.9 au and width of 13.5 ± 1.8 au. We determine a best-fit eccentricity of 0.12 ± 0.01. Assuming a size distribution power-law index of q = 3.46 ± 0.09, we constrain the dust absorptivity power-law index β to be 0.9 < β < 1.5. The geometry of the disk is robustly constrained with inclination 65.°6 ± 0.°3, position angle 337.°9 ± 0.°3, and argument of periastron 22.°5 ± 4.°3. Our observations do not confirm any of the azimuthal features found in previous imaging studies of the disk with Hubble Space Telescope, SCUBA, and ALMA. However, we cannot rule out structures ≤10 au in size or that only affect smaller grains. The central star is clearly detected with a flux density of 0.75 ± 0.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.

We present H band spectroscopic and Hμ photometric observations of HD 100546 obtained with the Gemini Planet Imager and the Magellan Visible AO camera. We detect H band emission at the location of the protoplanet HD 100546 b, but show that the 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 years, inconsistent at the 2σ 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 this is consistent with the spectrum of the star and the spectra extracted at other locations in the disk. A non-detection at the 5σ level of HD 100546 b in differential Hμ imaging places an upper limit, assuming the protoplanet lies in a gap free of extinction, on the accretion luminosity of 1.7 ? 10-4 L o and for 1 R Jup. These limits are comparable to the accretion luminosity and accretion rate of T-Tauri stars or LkCa 15 b. Taken together, these lines of evidence suggest that the H band source at the location of HD 100546 b is not emitted by a planetary photosphere or an accreting circumplanetary disk but is a disk feature enhanced by the point-spread function 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 100546 b 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α. 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 data sets, 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 HD 100546 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 (GPI) as a part of the GPI Exoplanet Survey campaign. Images of HD 984 B were obtained in the J (1.12-1.3 μm) and H (1.50-1.80 μm) 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°with a 68% confidence interval between 114°and 125°. To address the 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 them 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 K. Photometry of the companion yields a luminosity of log(Lbol/L⊙)=2.88 ± 0.07 dex with 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 located at Gemini South, we measured the near-infrared (1.0-2.4 μm) spectrum of the planetary companion to the nearby, young star β 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 ±0.2 , an effective temperature of 1724 ±15 K, a radius of 1.46 ±0.01 , and a surface gravity of [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 -4.0 [dex] depending upon the model. These values agree well with other publications and with "hot-start" predictions from planetary evolution models. Further, we find that the spectrum of β Pic b best matches a low surface gravity L2 ±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 a system, as it harbors a confirmed exterior 11 MJup 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.

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 β Pictoris (β Pic) 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 MAa = 1.10 ± 0.10 M⊙ and MAb = 0.290 ± 0.018 M⊙. 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 β Pic moving group. We derive a model-dependent age for the β 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.

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 ∼80-100 au. The new GPI data complement recent Hubble Space Telescope/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 ∼70° and a 1.5 ± 0.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.

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ève 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 to account for correlated noise. In the case of GPI data on β Pic b, this technique, which we call Bayesian KLIP-FM Astrometry (BKA), outperforms previous techniques and yields 1σ errors at or below the one milliarcsecond level. We validate BKA by fitting a Keplerian orbit to 12 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 β Pic b firmly rule out the possibility of a transit of the planet at 10-σ significance. However, we confirm that the Hill sphere of β 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.

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 μm that further constrains its outer morphology (projected separations of 27-135 au). We also present complementary Gemini Planet Imager 1.6 μm 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 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 ± 0.3 au (0.″618 ± 0. ″004) from the star. With the high astrometric precision afforded by GPI, we have confirmed, to more than 5?, the common proper motion of HR 2562B with the star, with only a month-long time baseline between observations. Spectral data in the J-, H-, and K-bands show a morphological similarity to L/T transition objects. We assign a spectral type of L7 ± 3 to HR 2562Band derive a luminosity of log(Lbol/L⊙) = -4.62 ± 0.12, corresponding to a mass of 30 ± 15 MJup 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 that is 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.

Using the Gemini Planet Imager, 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° along the spine of emission. A slight inclination and asymmetric warp are covariant and alter the interpretation of the observed disk emission. We employ three point-spread function 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 example 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.

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 pCL99.73% ≤ 2.4%. We discuss our results in the context of T dwarf cloud models and photometric variability.

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 the H band, almost three times closer to the star than the previous Hubble Space Telescope (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 that is inferred to be nearer to us. We also detect a lateral asymmetry in the south possibly due to shadowing from material within the IWA. This likely corresponds to a temporally variable azimuthal asymmetry observed in HST/STIS coronagraphic imaging.

An important technique for discovering and characterizing planets beyond our solar system relies upon measurement of weak Doppler shifts in the spectra of host stars induced by the influence of orbiting planets. A recent advance has been the introduction of optical frequency combs as frequency references. Frequency combs produce a series of equally spaced reference frequencies and they offer extreme accuracy and spectral grasp that can potentially revolutionize exoplanet detection. Here we demonstrate a laser frequency comb using an alternate comb generation method based on electro-optical modulation, with the comb centre wavelength stabilized to a molecular or atomic reference. In contrast to mode-locked combs, the line spacing is readily resolvable using typical astronomical grating spectrographs. Built using commercial off-the-shelf components, the instrument is relatively simple and reliable. Proof of concept experiments operated at near-infrared wavelengths were carried out at the NASA Infrared Telescope Facility and the Keck-II telescope.

We present the first scattered-light image of the debris disk around HD 131835 in the H band using the Gemini Planet Imager. HD 131835 is a ∼15 Myr old A2IV star at a distance of ∼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, in scattered light the disk shows similar orientation but different morphology. The scattered-light disk extends from ∼75 to ∼210 AU in the disk plane with roughly flat surface density. Our Monte Carlo radiative transfer model can 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σ level.

We present the first scattered light detections of the HD 106906 debris disk using the Gemini/Gemini Planet Imager in the infrared and Hubble Space Telescope (HST)/Advanced Camera for Surveys 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 ∼50 AU, and an outer extent >500 AU. The HST data show that the outer regions are highly asymmetric, resembling the "needle" morphology seen for the HD 15115 debris disk. The planet candidate is oriented ∼21° 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 new Gemini Planet Imager observations of the young exoplanet 51 Eridani b that 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 × 10-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 AU, corresponding to a period of years (assuming a mass of 1.75 Mo for the central star), and an inclination of deg. The remaining orbital elements are only marginally constrained by the current measurements. These preliminary values suggest an orbit that 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 compositions and luminosities, which are influenced by their formation mechanisms. Using the Gemini Planet Imager, we discovered a planet orbiting the ~20-million-year-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 (normalized by the luminosity of the Sun) of 1.6 to 4.0 × 10(-6) and an effective temperature of 600 to 750 kelvin. 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 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 (SE) and northwest (NW) sides. The SE 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 NW side at similar separations. This part of the disk is also vertically offset by 69 ± 30 mas to the northeast at 1″ when compared to the established disk midplane and is consistent with prior Atacama Large Millimeter/submillimeter Array and Hubble Space Telescope/Space Telescope Imaging Spectrograph observations. We see hints that the SE 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 NW side, where the disk midplane is offset in the opposite direction ∼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 ∼4 MJup 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 β 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 β 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 () with a position angle (PA) (slightly offset from the main outer disk, ), that extends from an inner disk radius of to well outside GPI's field of view. In addition, we present an updated orbit for β Pic b based on new astrometric measurements taken in GPI's spectroscopic mode spanning 14 months. The planet has a semimajor axis of , with an eccentricity The PA of the ascending node is offset from both the outer main disk and the inner disk seen in the GPI image. The orbital fit constrains the stellar mass of β Pic to Dynamical sculpting by β 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 β Pic b; (2) the mutual inclination of the inner disk and β Pic b is when it is expected to be closer to zero; and (3) the aspect ratio of the disk () is larger than expected from interactions with β Pic b or self-stirring by the disk's parent bodies.

Context. Thermal emission by stochastically heated dust grains (SHGs) plays an important role in the radiative transfer (RT) problem for a dusty medium. It is therefore essential to verify that RT codes properly calculate the dust emission before studying the effects of spatial distribution and other model parameters on the simulated observables. Aims. We define an appropriate problem for benchmarking dust emissivity calculations in the context of RT simulations, specifically including the emission from SHGs. Our aim is to provide a self-contained guide for implementors of such functionality and to offer insight into the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. Methods. 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. It also includes a series of analytically defined radiation fields to which the dust population is to be exposed and instructions for the desired output. We processed this problem using six RT codes participating in this benchmark effort and compared the results to a reference solution computed with the publicly available dust emission code DustEM. Results. 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 3 μm ≤1000 μm. Conclusions. We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of SHGs. This work can serve as a reference for implementors of dust RT codes, and it will pave the way for a more extensive benchmark effort focusing on the RT aspects of the various codes.

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-12.8+16.2 pc. The extended disk has been detected to ∼1.″5 (200 AU) at 11.7 and 18.3 μm with T-ReCS on Gemini South. The disk is inclined at an angle of ∼75° with the position angle of ∼61°. The flux of the HD 131835 system is 49.3 ± 7.6 and 84 ± 45 mJy at 11.7 and 18.3 μ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 present the first results from the polarimetry mode of the Gemini Planet Imager (GPI), which uses anew 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 (PSF) 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 thewest 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 Fring.

Using the recently installed Gemini Planet Imager (GPI), we have obtained the first H-band spectrum of the planetary companion to the nearby young star β 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 μm). The spectrum has a resolving power of ∼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 1600-1700K and a surface gravity of log(g) = 3.5-4.5 (cgs units). These values agree well with "hot-start" predictions from planetary evolution models for a gas giant with mass between 10 and 12 MJup and age between 10 and 20 Myr.

During the first-light run of the Gemini Planet Imager 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 μm 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.

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 μm Keck adaptive optics NIRC2 scattered-light observations of the HD 32297 debris disk reduced using a variant of the "locally optimized combination of images" 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 toward the northwest relative to a straight fiducial midplane. © 2014. The American Astronomical Society. All rights reserved.

We present new optical coronagraphic data of Fomalhaut obtained with HST/STIS in 2010 and 2012. 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 Fomalhaut b. The two epochs of Space Telescope Imaging Spectrograph (STIS) photometry exclude optical variability greater than 35%. A Markov chain Monte Carlo analysis demonstrates that the orbit of Fomalhaut b is highly eccentric, with e = 0.8 ± 0.1, a = 177 ± 68 AU, and q = 32 ± 24 AU. Fomalhaut b is apsidally aligned with the belt and 90% of allowed orbits have mutual inclination ≤36°. Fomalhaut b's orbit is belt crossing in the sky plane projection, but only 12% of possible orbits have ascending or descending nodes within a 25 AU wide belt annulus. The high eccentricity invokes a dynamical history where Fomalhaut 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 Tisserand parameter with respect to a hypothetical Fomalhaut planet at 30 AU or 120 AU lies in the range 2-3, similar to highly eccentric dwarf planets in our solar system. We argue that Fomalhaut 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. 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. © 2013. The American Astronomical Society. All rights reserved.

Aim: To consider the currently available knowledge and understanding of the symptom of urgency. Materials & Methods: Each faculty member reviewed the literature base of a different aspect of urgency and along with their personal clinical experience provided a base of evidence for discussion. Results: This overview summarises relevant published literature and the current clinical experience of the authors. Discussion: Whilst the mechanisms producing the sensation of urgency are still not fully understood and we are working within a definition that may complicate measurement and treatment, our pressing need is to effectively manage our patients for whom the practical nature of urgency can be all too apparent. Conclusion: Health care professionals have an important role to play today in helping to alleviate the widespread problem of urgency and its consequences. [ABSTRACT FROM AUTHOR]

Background: Women with urge urinary incontinence are commonly treated with antimuscarinic medications, but many discontinue therapy.Objective: To determine whether combining antimuscarinic drug therapy with supervised behavioral training, compared with drug therapy alone, improves the ability of women with urge incontinence to achieve clinically important reductions in incontinence episodes and to sustain these improvements after discontinuing drug therapy.Design: 2-stage, multicenter, randomized clinical trial conducted from July 2004 to January 2006.Setting: 9 university-affiliated outpatient clinics.Patients: 307 women with urge-predominant incontinence.Intervention: 10 weeks of open-label, extended-release tolterodine alone (n = 153) or combined with behavioral training (n = 154), followed by discontinuation of therapy and follow-up at 8 months.Measurements: The primary outcome, measured at 8 months, was no receipt of drugs or other therapy for urge incontinence and a 70% or greater reduction in frequency of incontinence episodes. Secondary outcomes were reduction in incontinence, self-reported satisfaction and improvement, and scores on validated questionnaires measuring symptom distress and bother and health-related quality of life. Study staff who performed outcome evaluations, but not participants and interventionists, were blinded to group assignment.Results: 237 participants completed the trial. According to life-table estimates, the rate of successful discontinuation of therapy at 8 months was the same in the combination therapy and drug therapy alone groups (41% in both groups; difference, 0 percentage points [95% CI, -12 to 12 percentage points]). A higher proportion of participants who received combination therapy than drug therapy alone achieved a 70% or greater reduction in incontinence at 10 weeks (69% vs. 58%; difference, 11 percentage points [CI, -0.3 to 22.1 percentage points]). Combination therapy yielded better outcomes over time on the Urogenital Distress Inventory and the Overactive Bladder Questionnaire (both P <0.001) at both time points for patient satisfaction and perceived improvement but not health-related quality of life. Adverse events were uncommon (12 events in 6 participants [3 in each group]).Limitations: Behavioral therapy components (daily bladder diary and recommendations for fluid management) in the group receiving drug therapy alone may have attenuated between-group differences. Assigned treatment was completed by 68% of participants, whereas 8-month outcome status was assessed on 77%.Conclusion: The addition of behavioral training to drug therapy may reduce incontinence frequency during active treatment but does not improve the ability to discontinue drug therapy and maintain improvement in urinary incontinence. Combination therapy has a beneficial effect on patient satisfaction, perceived improvement, and reduction of other bladder symptoms. [ABSTRACT FROM AUTHOR]

BACKGROUND AND PURPOSE: The purpose of this study was to describe how clinical pelvic-floor muscle (PFM) strength (force-generating capacity) is related to patient characteristics, lower urinary tract symptoms, and fecal incontinence symptoms. SUBJECTS: Data were obtained from 643 women who were participating in a randomized surgical trial for treatment of stress urinary incontinence. METHODS: Patient demographic variables, baseline urinary and fecal incontinence symptom questionnaires, urodynamic data and urinary diary data, pad test results, and standardized assessment of pelvic organ support were compared with PFM strength as described by the Brink scoring system. Bivariate analysis of factors associated with the Brink scale score was done using analysis of variance and linear regression. Multivariate analysis included patient variables that were significant on bivariate analysis. RESULTS: The mean Brink scale score was 9 (SD=2) and did not vary widely in this large, but highly select, patient sample. We found a weak, but statistically strong, relationship between age and Brink score. Brink scores were not related to diary and pad test measures of incontinence severity. DISCUSSION AND CONCLUSION: Overall, PFM strength was good in this sample of women with stress incontinence. Scores tended to be similar, and it is possible that the Brink scale does not reflect real clinical differences in PFM strength. [ABSTRACT FROM AUTHOR]