Source author record

Paul Kalas

Paul Kalas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2022arXiv

A Multi-Wavelength Study of the Highly Asymmetrical Debris Disk Around HD 111520

We observed the nearly edge-on debris disk system HD 111520 at $J$, $H$, & $K1$ near infrared (NIR) bands using both the spectral and polarization modes of the Gemini Planet Imager (GPI). With these new observations, we have performed an empirical analysis in order to better understand the disk morphology and its highly asymmetrical nature. We find that the disk features a large brightness and radial asymmetry, most prominent at shorter wavelengths. We also find that the radial location of the peak polarized intensity differs on either side of the star by 11 AU, suggesting that the disk may be eccentric, although, such an eccentricity does not fully explain the large brightness and radial asymmetry observed. Observations of the disk halo with HST also show the disk to be warped at larger separations, with a bifurcation feature in the northwest, further suggesting that there may be a planet in this system creating an asymmetrical disk structure. Measuring the disk color shows that the brighter extension is bluer compared to the dimmer extension, suggesting that the two sides have different dust grain properties. This finding, along with the large brightness asymmetry, are consistent with the hypothesis that a giant impact occurred between two large bodies in the northern extension of the disk, although confirming this based on NIR observations alone is not feasible. Follow-up imaging with ALMA to resolve the asymmetry in the dust mass distribution is essential in order to confirm this scenario.

preprint2022arXiv

Search for Stellar Flybys in the Sco-Cen OB Association with the Gaia DR2

High-contrast imaging studies of debris disks have revealed a significant diversity in their morphologies, including large-scale asymmetries. Theories involving stellar flybys, an external source of gravitational disturbance, have offered a plausible explanation for the origin of these morphological variations. Our study is an experiment to gain empirical evidence that has been lacking from such theories. We explore this paradigm by using astrometric and radial velocity measurements from the Gaia DR2 and ground-based observations to trace the trajectories of 625 stars in the Sco-Cen OB Association from 5 Myr in the past to 2 Myr in the future. We identified 119 stars that had at least one past flyby event occurring within one Hill radius, and 23 of these experienced flybys within 0.5 Hill radius. We found no evidence of a significant correlation between the presence of flyby events and infrared excess detections, although the sample is not uniformly sensitive to infrared excess emission. Ten stars that had past flyby events host resolved circumstellar disks that appear relatively symmetric in the existing data except for the circumbinary disk surrounding HD 106906. We determined the trajectory and relative velocity of each of these flyby events, and compared these to the geometry of the spatially-resolved disks. Future work is needed to measure the kinematics of lower mass stars and to improve sensitivity to circumstellar disks for the entire sample.

preprint2020arXiv

Debris Disk Results from the Gemini Planet Imager Exoplanet Survey's Polarimetric Imaging Campaign

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 $μ$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 °$). 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.

preprint2020arXiv

HD 165054: an astrometric calibration field for high-contrast imagers in Baade's Window

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.

preprint2020arXiv

Multiband GPI Imaging of the HR 4796A Debris Disk

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.

preprint2020arXiv

Multiband Polarimetric Imaging of HR 4796A with the Gemini Planet Imager

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.

preprint2020arXiv

On the Chemical Abundance of HR 8799 and the Planet c

Comparing chemical abundances of a planet and the host star reveals the origin and formation path. Stellar abundance is measured with high-resolution spectroscopy. Planet abundance, on the other hand, is usually inferred from low-resolution data. For directly imaged exoplanets, the data are available from a slew of high-contrast imaging/spectroscopy instruments. Here, we study the chemical abundance of HR 8799 and its planet c. We measure stellar abundance using LBT/PEPSI (R=120,000) and archival HARPS data: stellar [C/H], [O/H], and C/O are 0.11$\pm$0.12, 0.12$\pm$0.14, and 0.54$^{+0.12}_{-0.09}$, all consistent with solar values. We conduct atmospheric retrieval using newly obtained Subaru/CHARIS data together with archival Gemini/GPI and Keck/OSIRIS data. We model the planet spectrum with petitRADTRANS and conduct retrieval using PyMultiNest. Retrieved planetary abundance can vary by $\sim$0.5 dex, from sub-stellar to stellar C and O abundances. The variation depends on whether strong priors are chosen to ensure a reasonable planet mass. Moreover, comparison with previous works also reveals inconsistency in abundance measurements. We discuss potential issues that can cause the inconsistency, e.g., systematics in individual data sets and different assumptions in the physics and chemistry in retrieval. We conclude that no robust retrieval can be obtained unless the issues are fully resolved.

preprint2020arXiv

Revised Astrometric Calibration of the Gemini Planet Imager

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.

preprint2020arXiv

The Breakthrough Listen Search for Intelligent Life: A 3.95-8.00 GHz Search for Radio Technosignatures in the Restricted Earth Transit Zone

We report on a search for artificial narrowband signals of 20 stars within the restricted Earth Transit Zone as a part of the ten-year Breakthrough Listen (BL) search for extraterrestrial intelligence. The restricted Earth Transit Zone is the region of the sky from which an observer would see the Earth transit the Sun with an impact parameter of less than 0.5. This region of the sky is geometrically unique, providing a potential way for an extraterrestrial intelligence to discover the Solar System. The targets were nearby (7-143 pc) and the search covered an electromagnetic frequency range of 3.95-8.00 GHz. We used the Robert C. Byrd Green Bank Telescope to perform these observations with the standard BL data recorder. We searched these data for artificial narrowband ($\sim$Hz) signals with Doppler drift rates of $\pm 20$ Hz s$^{-1}$. We found one set of potential candidate signals on the target HIP 109656 which was then found to be consistent with known properties of anthropogenic radio frequency interference. We find no evidence for radio technosignatures from extraterrestrial intelligence in our observations. The observing campaign achieved a minimum detectable flux which would have allowed detections of emissions that were $10^{-3}$ to $0.88$ times as powerful as the signaling capability of the Arecibo radar transmitter, for the nearest and furthest stars respectively. We conclude that at least $8\%$ of the systems in the restricted Earth Transit Zone within 150 pc do not possess the type of transmitters searched in this survey. To our knowledge, this is the first targeted search for extraterrestrial intelligence of the restricted Earth Transit Zone. All data used in this paper are publicly available via the Breakthrough Listen Public Data Archive (http://seti.berkeley.edu/bldr2).

preprint2020arXiv

The Gemini Planet Imager view of the HD 32297 debris disk

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.

preprint2019arXiv

An updated visual orbit of the directly-imaged exoplanet 51 Eridani b and prospects for a dynamical mass measurement with Gaia

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$σ$ 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}$).

preprint2019arXiv

First Resolved Scattered-Light Images of Four Debris Disks in Scorpius-Centaurus with the Gemini Planet Imager

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.

preprint2019arXiv

The Gemini Planet Imager Exoplanet Survey: Dynamical Mass of the Exoplanet beta Pictoris b from Combined Direct Imaging and Astrometry

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.

preprint2019arXiv

The Possible Astrometric Signature of a Planetary-mass Companion to the Nearby Young Star TW Piscis Austrini (Fomalhaut B): Constraints from Astrometry, Radial Velocities, and Direct Imaging

We present constraints on the presence of substellar companions to the nearby ($d\sim7.6$ pc), young ($440\pm40$ Myr) K4Ve star TW Piscis Austrini, the wide ($\sim$0.3 pc) companion to the A4V star Fomalhaut. We combined absolute astrometry from Hipparcos and Gaia with literature radial velocity measurements and dedicated high-contrast imaging observations obtained with Keck/NIRC2 to achieve sensitivity to brown dwarf and planetary-mass companions ($\gtrsim2$ M$_{\rm Jup}$) over many decades of orbital period ($\lesssim 10^3$ yr). The significant astrometric acceleration measured between the Hipparcos and Gaia catalogues, reported previously in the literature, cannot be explained by the orbital motion of TW PsA around the barycenter of the Fomalhaut triple system. Instead, we find that it is consistent with the reflex motion induced by an orbiting substellar companion. The combination of astrometry, radial velocities, and a deep $L^{\prime}$ imaging dataset leads to a constraint on the companion mass of $1.2_{-0.6}^{+0.7}$ M$_{\rm Jup}$. However, the period of the companion is poorly constrained, with a highly multi-modal period posterior distribution due to aliasing with the 24.25 year baseline between Hipparcos and Gaia. If confirmed through continued astrometric or spectroscopic monitoring, or via direct detection, the companion to TW PsA would represent a choice target for detailed atmospheric characterization with high-contrast instruments on the upcoming James Webb Space Telescope and Wide Field Infrared Survey Telescope.

preprint2016arXiv

Bringing "The Moth" to Light: A Planet-Sculpting Scenario for the HD 61005 Debris Disk

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.

preprint2016arXiv

Constraints on the architecture of the HD 95086 planetary system with the Gemini Planet Imager

We present astrometric monitoring of the young exoplanet HD 95086 b obtained with the Gemini Planet Imager between 2013 and 2016. A small but significant position angle change is detected at constant separation; the orbital motion is confirmed with literature measurements. Efficient Monte Carlo techniques place preliminary constraints on the orbital parameters of HD 95086 b. With 68% confidence, a semimajor axis of 61.7^{+20.7}_{-8.4} au and an inclination of 153.0^{+9.7}_{-13.5} deg are favored, with eccentricity less than 0.21. Under the assumption of a co-planar planet-disk system, the periastron of HD 95086 b is beyond 51 au with 68% confidence. Therefore HD 95086 b cannot carve the entire gap inferred from the measured infrared excess in the SED of HD 95086. We use our sensitivity to additional planets to discuss specific scenarios presented in the literature to explain the geometry of the debris belts. We suggest that either two planets on moderately eccentric orbits or three to four planets with inhomogeneous masses and orbital properties are possible. The sensitivity to additional planetary companions within the observations presented in this study can be used to help further constrain future dynamical simulations of the planet-disk system.

preprint2016arXiv

Discovery of a Substellar Companion to the Nearby Debris Disk Host HR 2562

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σ$. 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.

preprint2016arXiv

Imaging an 80 AU Radius Dust Ring Around the F5V Star HD 157587

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.

preprint2016arXiv

Spectroscopic characterization of HD 95086 b with the Gemini Planet Imager

We present new $H$ (1.5-1.8 $μ$m) photometric and $K_1$ (1.9-2.2 $μ$m) spectroscopic observations of the young exoplanet HD 95086 b obtained with the Gemini Planet Imager. The $H$-band magnitude has been significantly improved relative to previous measurements, whereas the low resolution $K_1$ ($λ/δλ\approx 66$) spectrum is featureless within the measurement uncertainties, and presents a monotonically increasing pseudo-continuum consistent with a cloudy atmosphere. By combining these new measurements with literature $L^{\prime}$ photometry, we compare the spectral energy distribution of the planet to other young planetary-mass companions, field brown dwarfs, and to the predictions of grids of model atmospheres. HD 95086 b is over a magnitude redder in $K_1-L^{\prime}$ color than 2MASS J12073346-3932539 b and HR 8799 c and d, despite having a similar $L^{\prime}$ magnitude. Considering only the near-infrared measurements, HD 95086 b is most analogous to the brown dwarfs 2MASS J2244316+204343 and 2MASS J21481633+4003594, both of which are thought to have dusty atmospheres. Morphologically, the spectral energy distribution of HD 95086 b is best fit by low temperature ($T_{\rm eff} =$ 800-1300 K), low surface gravity spectra from models which simulate high photospheric dust content. This range of effective temperatures is consistent with field L/T transition objects, but the spectral type of HD 95086 b is poorly constrained between early L and late T due to its unusual position the color-magnitude diagram, demonstrating the difficulty in spectral typing young, low surface gravity substellar objects. As one of the reddest such objects, HD 95086 b represents an important empirical benchmark against which our current understanding of the atmospheric properties of young extrasolar planets can be tested.

preprint2016arXiv

The Orbit and Transit Prospects for $β$ Pictoris b constrained with One Milliarcsecond Astrometry

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 (MCMC) 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 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 $β$ 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.

preprint2016arXiv

The PDS 66 Circumstellar Disk as seen in Polarized Light with the Gemini Planet Imager

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.

preprint2016arXiv

The Peculiar Debris Disk of HD 111520 as Resolved by the Gemini Planet Imager

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.

preprint2015arXiv

$β$ Pictoris' inner disk in polarized light and new orbital parameters for $β$ Pictoris b

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 ($ϕ=85.27°^{+0.26}_{-0.19}$) with a position angle $θ_{PA}=30.35°^{+0.29}_{-0.28}$ (slightly offset from the main outer disk, $θ_{PA}\approx29°$), 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 $β$ 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 $Ω=31.75°\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 $β$ Pic to $1.60\pm0.05 M_{\odot}$. 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 $4°$, 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 $β$ Pic $b$ or self-stirring by the disk's parent bodies.

preprint2015arXiv

Gemini Planet Imager Observations of the AU Microscopii Debris Disk: Asymmetries within One Arcsecond

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.

preprint2015arXiv

The AU Mic Debris Disk: far-infrared and submillimeter resolved imaging

We present far-infrared and submillimeter maps from the Herschel Space Observatory and the James Clerk Maxwell Telescope of the debris disk host star AU Microscopii. Disk emission is detected at 70, 160, 250, 350, 450, 500 and 850 micron. The disk is resolved at 70, 160 and 450 micron. In addition to the planetesimal belt, we detect thermal emission from AU Mic's halo for the first time. In contrast to the scattered light images, no asymmetries are evident in the disk. The fractional luminosity of the disk is $3.9 \times 10^{-4}$ and its mm-grain dust mass is 0.01 MEarth (+/- 20%). We create a simple spatial model that reconciles the disk SED as a blackbody of 53 +/- 2 K (a composite of 39 and 50 K components) and the presence of small (non-blackbody) grains which populate the extended halo. The best fit model is consistent with the "birth ring" model explored in earlier works, i.e., an edge-on dust belt extending from 8.8-40 AU, but with an additional halo component with an $r^{-1.5}$ surface density profile extending to the limits of sensitivity (140 AU). We confirm that AU Mic does not exert enough radiation force to blow out grains. For stellar mass loss rates of 10-100x solar, compact (zero porosity) grains can only be removed if they are very small, consistently with previous work, if the porosity is 0.9, then grains approaching 0.1 micron can be removed via corpuscular forces (i.e., the stellar wind).

preprint2014arXiv

Gemini Planet Imager Observational Calibrations V: Astrometry and Distortion

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$°$. 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.

preprint2014arXiv

Gemini Planet Imager Observational Calibrations VIII: Characterization and Role of Satellite Spots

The Gemini Planet Imager (GPI) combines extreme adaptive optics, an integral field spectrograph, and a high performance coronagraph to directly image extrasolar planets in the near-infrared. Because the coronagraph blocks most of the light from the star, it prevents the properties of the host star from being measured directly. Instead, satellite spots, which are created by diffraction from a square grid in the pupil plane, can be used to locate the star and extract its spectrum. We describe the techniques implemented into the GPI Data Reduction Pipeline to measure the properties of the satellite spots and discuss the precision of the reconstructed astrometry and spectrophotometry of the occulted star. We find the astrometric precision of the satellite spots in an $H$-band datacube to be $0.05$ pixels and is best when individual satellite spots have a signal to noise ratio (SNR) of $> 20$. In regards to satellite spot spectrophotometry, we find that the total flux from the satellite spots is stable to $\sim 7\%$ and scales linearly with central star brightness and that the shape of the satellite spot spectrum varies on the $2\%$ level.

preprint2014arXiv

Hubble Space Telescope High Resolution Imaging of Kepler Small and Cool Exoplanet Host Stars

High resolution imaging is an important tool for follow-up study of exoplanet candidates found via transit detection with the Kepler Mission. We discuss here HST imaging with the WFC3 of 23 stars that host particularly interesting Kepler planet candidates based on their small size and cool equilibrium temperature estimates. Results include detections, exclusion of background stars that could be a source of false positives for the transits, and detection of physically-associated companions in a number of cases providing dilution measures necessary for planet parameter refinement. For six KOIs, we find that there is ambiguity in which star hosts the transiting planet(s), with potentially strong implications for planetary characteristics. Our sample is evenly distributed in G, K, and M spectral types. Albeit with a small sample size, we find that physically-associated binaries are more common than expected at each spectral type, reaching a factor of 10 frequency excess at M. We document the program detection sensitivities, detections, and deliverables to the Kepler follow-up program archive.

preprint2014arXiv

Insights on the dynamical history of the Fomalhaut system - Investigating the Fom c hypothesis

The eccentric shape of the debris disk observed around Fomalhaut was first attributed to Fom b, a companion detected near the belt inner-edge, but new constraints on its orbit revealed that it is belt-crossing, highly eccentric $(e \sim 0.6-0.9)$, and can hardly account for the shape of the belt. The best scenario to explain this paradox is that there is another massive body in this system, Fom c, which drives the debris disk shape. The resulting planetary system is highly unstable, which hints at a dynamical scenario involving a recent scattering of Fom b on its current orbit, potentially with the putative Fom c. Our goal is to give insights on the probability for Fom b to have been set on its highly eccentric orbit by a close-encounter with the putative Fom c. We aim to study in particular the part played by mean-motion resonances with Fom c, which could have brought Fom b sufficiently close to Fom c for it to be scattered on its current orbit, but also delay this scattering event. Using N-body simulations, we found that the generation of orbits similar to that of Fom b, either in term of dimensions or orientation, is a robust process involving a scattering event and a further secular evolution of inner material with an eccentric massive body such as the putative Fom c. We found in particular that mean-motion resonances can delay scattering events, and thus the production of Fom b-like orbits, on timescales comparable to the age of the system, thus explaining the witnessing of an unstable configuration. We conclude that Fom b probably originated from an inner resonance with Fom c, which is at least Neptune-Saturn size, and was set on its current orbit by a scattering event with Fom c. Since Fom b could not have formed from material in resonance, our scenario also hints at former migration processes in this planetary system.

preprint2014arXiv

The First H-band Spectrum of the Massive Gas Giant Planet beta Pictoris b with the Gemini Planet Imager

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.

preprint2013arXiv

An independent determination of Fomalhaut b's orbit and the dynamical effects on the outer dust belt

The nearby star Fomalhaut harbours a cold, moderately eccentric dust belt with a sharp inner edge near 133 au. A low-mass, common proper motion companion (Fom b), was discovered near the inner edge and was identified as a planet candidate that could account for the belt morphology. However, the most recent orbit determination based on four epochs of astrometry over eight years reveals a highly eccentric orbit that appears to cross the belt in the sky plane projection. We perform here a full orbital determination based on the available astrometric data to independently validate the orbit estimates previously presented. Adopting our values for the orbital elements and their associated uncertainties, we then study the dynamical interaction between the planet and the dust ring, to check whether the proposed disk sculpting scenario by Fom b is plausible. We used a dedicated MCMC code to derive the statistical distributions of the orbital elements of Fom b. Then we used symplectic N-body integration to investigate the dynamics of the dust belt, as perturbed by a single planet. Different attempts were made assuming different masses for Fom b. We also performed a semi-analytical study to explain our results. Our results are in good agreement with others regarding the orbit of Fom b. We find that the orbit is highly eccentric, is close to apsidally aligned with the belt, and has a moderate mutual inclination relative to the belt plane of. If coplanar, this orbit crosses the disk. Our dynamical study then reveals that the observed planet could sculpt a transient belt configuration with a similar eccentricity to what is observed, but it would not be simultaneously apsidally aligned with the planet. This transient configuration only occurs a short time after the planet is placed on such an orbit (assuming an initially circular disk), a time that is inversely proportional to the planet's mass, and that is in any case much less than the 440 Myr age of the star. We constrain how long the observed dust belt could have survived with Fom b on its current orbit, as a function of its possible mass. This analysis leads us to conclude that Fom b is likely to have low mass, that it is unlikely to be responsible for the sculpting of the belt, and that it supports the hypothesis of a more massive, less eccentric planet companion Fom c.

preprint2013arXiv

Modeling Self-Subtraction in Angular Differential Imaging: Application to the HD 32297 Debris Disk

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.

preprint2013arXiv

STIS Coronagraphic Imaging of Fomalhaut: Main Belt Structure and the Orbit of Fomalhaut b

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.

preprint2012arXiv

Millimeter Emission Structure in the first ALMA Image of the AU Mic Debris Disk

We present 1.3 millimeter ALMA Cycle 0 observations of the edge-on debris disk around the nearby, ~10 Myr-old, M-type star AU Mic. These observations obtain 0.6 arcsec (6 AU) resolution and reveal two distinct emission components: (1) the previously known dust belt that extends to a radius of 40 AU, and (2) a newly recognized central peak that remains unresolved. The cold dust belt of mass about 1 lunar mass is resolved in the radial direction with a rising emission profile that peaks sharply at the location of the outer edge of the "birth ring" of planetesimals hypothesized to explain the midplane scattered light gradients. No significant asymmetries are discerned in the structure or position of this dust belt. The central peak identified in the ALMA image is ~6 times brighter than the stellar photosphere, which indicates an additional emission process in the inner regions of the system. Emission from a stellar corona or activity may contribute, but the observations show no signs of temporal variations characteristic of radio-wave flares. We suggest that this central component may be dominated by dust emission from an inner planetesimal belt of mass about 0.01 lunar mass, consistent with a lack of emission shortward of 25 microns and a location <3 AU from the star. Future millimeter observations can test this assertion, as an inner dust belt should be readily separated from the central star at higher angular resolution.

preprint2012arXiv

Resolved Debris Discs Around A Stars in the Herschel DEBRIS Survey

The majority of debris discs discovered so far have only been detected through infrared excess emission above stellar photospheres. While disc properties can be inferred from unresolved photometry alone under various assumptions for the physical properties of dust grains, there is a degeneracy between disc radius and dust temperature that depends on the grain size distribution and optical properties. By resolving the disc we can measure the actual location of the dust. The launch of Herschel, with an angular resolution superior to previous far-infrared telescopes, allows us to spatially resolve more discs and locate the dust directly. Here we present the nine resolved discs around A stars between 20 and 40 pc observed by the DEBRIS survey. We use these data to investigate the disc radii by fitting narrow ring models to images at 70, 100 and 160 μm and by fitting blackbodies to full spectral energy distributions. We do this with the aim of finding an improved way of estimating disc radii for unresolved systems. The ratio between the resolved and blackbody radii varies between 1 and 2.5. This ratio is inversely correlated with luminosity and any remaining discrepancies are most likely explained by differences to the minimum size of grain in the size distribution or differences in composition. We find that three of the systems are well fit by a narrow ring, two systems are borderline cases and the other four likely require wider or multiple rings to fully explain the observations, reflecting the diversity of planetary systems.

preprint2011arXiv

Radio Interferometric Planet Search II: Constraints on sub-Jupiter-Mass Companions to GJ 896A

We present results from the Radio Interferometric Planet (RIPL) search for compan- ions to the nearby star GJ 896A. We present 11 observations over 4.9 years. Fitting astrometric parameters to the data reveals a residual with peak-to-peak amplitude of ~ 3 mas in right ascension. This residual is well-fit by an acceleration term of 0.458 \pm 0.032 mas/y^2. The parallax is fit to an accuracy of 0.2 mas and the proper motion terms are fit to accuracies of 0.01 mas/y. After fitting astrometric and acceleration terms residuals are 0.26 mas in each coordinate, demonstrating that stellar jitter does not limit the ability to carry out radio astrometric planet detection and characterization. The acceleration term originates in part from the companion GJ 896B but the amplitude of the acceleration in declination is not accurately predicted by the orbital model. The acceleration sets a mass upper limit of 0.15 MJ at a semi-major axis of 2 AU for a planetary companion to GJ 896A. For semi-major axes between 0.3 and 2 AU upper limits are determined by the maximum angular separation; the upper limits scale from the minimum value in proportion to the inverse of the radius. Upper limits at larger radii are set by the acceleration and scale as the radius squared. An improved solution for the stellar binary system could improve the exoplanet mass sensitivity by an order of magnitude.

preprint2009arXiv

Orbital Constraints on the beta Pic Inner Planet Candidate with Keck Adaptive Optics

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.

preprint2009arXiv

Radio Interferometric Planet Search I: First Constraints on Planetary Companions for Nearby, Low-Mass Stars from Radio Astrometry

Radio astrometry of nearby, low-mass stars has the potential to be a powerful tool for the discovery and characterization of planetary companions. We present a Very Large Array survey of 172 active M dwarfs at distances of less than 10 pc. Twenty nine stars were detected with flux densities greater than 100 microJy. We observed 7 of these stars with the Very Long Baseline Array at milliarcsecond resolution in three separate epochs. With a detection threshold of 500 microJy in images of sensitivity 1 sigma ~ 100 microJy, we detected three stars three times (GJ 65B, GJ896A, GJ 4247), one star twice (GJ 285), and one star once (GJ 803). Two stars were undetected (GJ 412B and GJ 1224). For the four stars detected in multiple epochs, residuals from the optically-determined proper motions have an rms deviation of ~0.2 milliarcseconds, consistent with statistical noise limits. Combined with previous optical astrometry, these residuals provide acceleration upper limits that allow us to exclude planetary companions more massive than 3-6 M_Jup at a distance of ~1 AU with a 99% confidence level.

preprint2001arXiv

Discovery of Reflection Nebulosity Around Five Vega-like Stars

Coronagraphic optical observations of six Vega-like stars reveal reflection nebulosities, five of which were previously unknown. The nebulosities illuminated by HD 4881, HD 23362, HD 23680, HD 26676, and HD 49662 resemble that of the Pleiades, indicating an interstellar origin for dust grains. The reflection nebulosity around HD 123160 has a double-arm morphology, but no disk-like feature is seen as close as 2.5 arcsec from the star in K-band adaptive optics data. We demonstrate that uniform density dust clouds surrounding HD 23362, HD 23680 and HD 123160 can account for the observed 12-100 micron spectral energy distributions. For HD 4881, HD 26676, and HD 49662 an additional emission source, such as from a circumstellar disk or non-equilibrium grain heating, is required to fit the 12-25 micron data. These results indicate that in some cases, particularly for Vega-like stars located beyond the Local Bubble (>100 pc), the dust responsible for excess thermal emission may originate from the interstellar medium rather than from a planetary debris system.