$epsilon$~Eridani is a young planetary system hosting a complex multi-belt debris disk and a confirmed Jupiter-like planet orbiting at 3.48 AU from its host star. Its age and architecture are thus reminiscent of the early Solar System. The most recen
t study of Mawet et al. 2019, which combined radial velocity (RV) data and Ms-band direct imaging upper limits, started to constrain the planets orbital parameters and mass, but are still affected by large error bars and degeneracies. Here we make use of the most recent data compilation from three different techniques to further refine $epsilon$~Eridani~bs properties: RVs, absolute astrometry measurements from the Hipparcos~and Gaia~missions, and new Keck/NIRC2 Ms-band vortex coronagraph images. We combine this data in a Bayesian framework. We find a new mass, $M_b$ = $0.66_{-0.09}^{+0.12}$~M$_{Jup}$, and inclination, $i$ = $77.95_{-21.06}^{circ+28.50}$, with at least a factor 2 improvement over previous uncertainties. We also report updated constraints on the longitude of the ascending node, the argument of the periastron, and the time of periastron passage. With these updated parameters, we can better predict the position of the planet at any past and future epoch, which can greatly help define the strategy and planning of future observations and with subsequent data analysis. In particular, these results can assist the search for a direct detection with JWST and the Nancy Grace Roman Space Telescopes coronagraph instrument (CGI).
Using the Keck Planet Imager and Characterizer (KPIC), we obtained high-resolution (R$sim$35,000) $K$-band spectra of the four planets orbiting HR 8799. We clearly detected water{} and CO in the atmospheres of HR 8799 c, d, and e, and tentatively det
ected a combination of CO and water{} in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral resolution to date when considering both their angular separations and flux ratios. We developed a forward modeling framework that allows us to jointly fit the spectra of the planets and the diffracted starlight simultaneously in a likelihood-based approach and obtained posterior probabilities on their effective temperatures, surface gravities, radial velocities, and spins. We measured $vsin(i)$ values of $10.1^{+2.8}_{-2.7}$~km/s for HR 8799 d and $15.0^{+2.3}_{-2.6}$~km/s for HR 8799 e, and placed an upper limit of $< 14$~km/s of HR 8799 c. Under two different assumptions of their obliquities, we found tentative evidence that rotation velocity is anti-correlated with companion mass, which could indicate that magnetic braking with a circumplanetary disk at early times is less efficient at spinning down lower mass planets.
The Keck Planet Imager and Characterizer (KPIC) is an upgrade to the Keck II adaptive optics system and instrument suite with the goal of improving direct imaging and high-resolution spectroscopic characterization capabilities for giant exoplanets. K
PIC Phase I includes a fiber injection unit (FIU) downstream of a new pyramid wavefront sensor, coupling planet light to a single mode fiber fed into NIRSPEC, Kecks high-resolution infrared spectrograph. This enables high-dispersion spectroscopy (HDS) of directly imaged exoplanets at smaller separation and higher contrast, improving our spectral characterization capabilities for these objects. Here, we report performance results from the KPIC Phase I FIU commissioning, including analysis of throughput, stability, and sensitivity of the instrument.
We present $L$-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wavefront sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a mod
el of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semi-major axis of $20^{+3}_{-4}$~au and PDS 70 c to have a semi-major axis of $34^{+12}_{-6}$~au (95% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2-3~$R_{Jup}$. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 $M_{textrm{Jup}}$ and a mean mass accretion rate between $3 times 10^{-7}$ and $8 times 10^{-7}~M_{textrm{Jup}}/textrm{yr}$. For PDS 70 c, we computed a mass between 1 and 3 $M_{textrm{Jup}}$ and mean mass accretion rate between $1 times 10^{-7}$ and $5 times~10^{-7} M_{textrm{Jup}}/textrm{yr}$. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets SEDs.
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 Baades Window, a region of the galactic plane with relatively low extinction from i
nterstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($sim 3$ stars per square arcsecond with $H<22$), necessary because of the small field-of-view of the Gemini Planet Imager. Using nine epochs spanning over five years, we have fit a standard five-parameter astrometric model to the astrometry of seven background stars within close proximity to HD 165054 (angular separation $< 2$ arcsec). We achieved a proper motion precision of $sim 0.3$ mas/yr, and constrained the parallax of each star to be $lesssim 1$ mas. Our measured proper motions and parallax limits are consistent with the background stars being a part of the galactic bulge. Using these measurements we find no evidence of any systematic trend of either the plate scale or the north angle offset of GPI between 2014 and 2019. We compared our model describing the motions of the seven background stars to observations of the same field in 2014 and 2018 obtained with Keck/NIRC2, an instrument with an excellent astrometric calibration. We find that predicted position of the background sources is consistent with that measured by NIRC2, within the uncertainties of the calibration of the two instruments. In the future, we will use this field as a standard astrometric calibrator for the upgrade of GPI and potentially for other high-contrast imagers.
We present 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 imagi
ng, 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 planets Hill sphere between 2017 April 11 and 2018 February 16 (+/- 18 days). Following the recent radial velocity detection of a second planet in the system, beta Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass for the imaged planet beta Pic b, of 8.0 +/- 2.6 M_Jup, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4.
We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at
an intermediate inclination ($i=136_{-11}^{+10}$,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semi-major axis of $11.1_{-1.3}^{+4.2}$,au and a period of $28.1_{-4.9}^{+17.2}$,yr, assuming a mass of 1.75,M$_{odot}$ for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigated the potential of using absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51~Eri derived from a comparison of the {it Hipparcos} and {it Gaia} catalogues was found to be inconsistent at the 2--3$sigma$ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs or the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using {it Gaia} astrometry alone for a dynamical mass measurement of the planet by simulating {it Gaia} measurements of the motion of the photocenter of the system over the course of the extended eight-year mission. We find that such a measurement is only possible ($>98$% probability) given the most optimistic predictions for the {it Gaia} scan astrometric uncertainties for bright stars, and a high mass for the planet ($gtrsim3.6$,M$_{rm Jup}$).
The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the systems dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milli
arcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assuming the planets are near 1:2:4:8 period commensurabilities, or are coplanar, does not worsen the fit. We added the prior that the planets must have been stable for the age of the system (40 Myr) by running orbit configurations from our posteriors through $N$-body simulations and varying the masses of the planets. We found that only assuming the planets are both coplanar and near 1:2:4:8 period commensurabilities produces dynamically stable orbits in large quantities. Our posterior of stable coplanar orbits tightly constrains the planets orbits, and we discuss implications for the outermost planet b shaping the debris disk. A four-planet resonance lock is not necessary for stability up to now. However, planet pairs d and e, and c and d, are each likely locked in two-body resonances for stability if their component masses are above $6~M_{rm{Jup}}$ and $7~M_{rm{Jup}}$, respectively. Combining the dynamical and luminosity constraints on the masses using hot-start evolutionary models and a system age of $42 pm 5$~Myr, we found the mass of planet b to be $5.8 pm 0.5~M_{rm{Jup}}$, and the masses of planets c, d, and e to be $7.2_{-0.7}^{+0.6}~M_{rm{Jup}}$ each.
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-Loeve Image Processing (KLIP) algorithm with Angul
ar and Spectral Differential Imaging (ADI and SDI). The KLIP-induced distortion of the astrophysical signal is included in the matched filter template by computing a forward model of the PSF at every position in the image. To optimize the performance of the algorithm, we conduct extensive planet injection and recovery tests and tune the exoplanet spectra template and KLIP reduction aggressiveness to maximize the Signal-to-Noise Ratio (SNR) of the recovered planets. We show that only two spectral templates are necessary to recover any young Jovian exoplanets with minimal SNR loss. We also developed a complete pipeline for the automated detection of point source candidates, the calculation of Receiver Operating Characteristics (ROC), false positives based contrast curves, and completeness contours. We process in a uniform manner more than 330 datasets from the Gemini Planet Imager Exoplanet Survey (GPIES) and assess GPI typical sensitivity as a function of the star and the hypothetical companion spectral type. This work allows for the first time a comparison of different detection algorithms at a survey scale accounting for both planet completeness and false positive rate. We show that the new forward model matched filter allows the detection of $50%$ fainter objects than a conventional cross-correlation technique with a Gaussian PSF template for the same false positive rate.
A principal scientific goal of the Gemini Planet Imager (GPI) is obtaining milliarcsecond astrometry to constrain exoplanet orbits. However, astrometry of directly imaged exoplanets is subject to biases, systematic errors, and speckle noise. Here we
describe an analytical procedure to forward model the signal of an exoplanet that accounts for both the observing strategy (angular and spectral differential imaging) and the data reduction method (Karhunen-Lo`eve Image Projection algorithm). We use this forward model to measure the position of an exoplanet in a Bayesian framework employing Gaussian processes and Markov chain Monte Carlo (MCMC) to account for correlated noise. In the case of GPI data on $beta$ Pic b, this technique, which we call Bayesian KLIP-FM Astrometry (BKA), outperforms previous techniques and yields 1$sigma$-errors at or below the one milliarcsecond level. We validate BKA by fitting a Keplerian orbit to twelve GPI observations along with previous astrometry from other instruments. The statistical properties of the residuals confirm that BKA is accurate and correctly estimates astrometric errors. Our constraints on the orbit of $beta$ Pic b firmly rule out the possibility of a transit of the planet at 10-$sigma$ significance. However, we confirm that the Hill sphere of $beta$ Pic b will transit, giving us a rare chance to probe the circumplanetary environment of a young, evolving exoplanet. We provide an ephemeris for photometric monitoring of the Hill sphere transit event, which will begin at the start of April in 2017 and finish at the end of January in 2018.
