We present L'-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wave front 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 model 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 semimajor axis of au and PDS 70 c to have a semimajor axis of 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 RJup. 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 MJup and a mean mass accretion rate between 3 × 10−7 and 8 × 10−7 MJup/yr. For PDS 70 c, we computed a mass between 1 and 3 MJup and mean mass accretion rate between 1 × 10−7 and 5 × 10−7 MJup/yr. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets' SEDs.
To characterize circumstellar systems in high-contrast imaging, the fundamental step is to construct a best point-spread function (PSF) template for the noncircumstellar signals (i.e., starlight and ...speckles) and separate it from the observation. With existing PSF construction methods, the circumstellar signals (e.g., planets, circumstellar disks) are unavoidably altered by overfitting and/or self-subtraction, making forward modeling a necessity to recover these signals. We present a forward modeling-free solution to these problems with data imputation using sequential nonnegative matrix factorization (DI-sNMF), which first converts this signal separation problem to a "missing data" problem in statistics by flagging the regions that host circumstellar signals as missing data, then attributes PSF signals to these regions. We mathematically prove it to have negligible alteration to circumstellar signals when the imputation region is relatively small, which thus enables precise measurement of these circumstellar objects. We apply it to simulated point-source and circumstellar disk observations to demonstrate its proper recovery of them. We apply it to Gemini Planet Imager K1-band observations of the debris disk surrounding HR 4796A, finding a tentative trend that the dust is more forward scattering as the wavelength increases. We expect DI-sNMF to be applicable to other general scenarios where the separation of signals is needed.
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.
Reference star differential imaging (RDI) is a powerful strategy for high-contrast imaging. Using example observations taken with the vortex coronagraph mode of Keck/NIRC2 in L′ band, we demonstrate ...that RDI provides improved sensitivity to point sources at small angular separations compared to angular differential imaging (ADI). Applying RDI to images of the low-mass stellar companions HIP 79124 C (192 mas separation, ΔL′ = 4.01) and HIP 78233 B (141 mas separation, ΔL′ = 4.78), the latter a first imaging detection, increases the significance of their detections by up to a factor of 5 with respect to ADI. We compare methods for reference frame selection and find that pre-selection of frames improves detection significance of point sources by up to a factor of 3. In addition, we use observations of the circumstellar disks around MWC 758 and 2MASS J16042165−2130284 to show that RDI allows for accurate mapping of scattered light distributions without self-subtraction artifacts.
The NIRC2 vortex coronagraph is an instrument on Keck II designed to directly image exoplanets and circumstellar disks at mid-infrared bands L′ (3.4-4.1 m) and Ms (4.55-4.8 m). We analyze imaging ...data and corresponding adaptive optics telemetry, observing conditions, and other metadata over a three-year time period to characterize the performance of the instrument and predict the detection limits of future observations. We systematically process images from 359 observations of 304 unique stars to subtract residual starlight (i.e., the coronagraphic point-spread function) of the target star using two methods: angular differential imaging (ADI) and reference star differential imaging (RDI). We find that for the typical parallactic angle (PA) rotation of our data set (∼10°), RDI provides gains over ADI for angular separations smaller than 0 25. Furthermore, we find a power-law relation between the angular separation from the host star and the minimum PA rotation required for ADI to outperform RDI, with a power-law index of −1.18 0.08. Finally, we use random forest models to estimate ADI and RDI post-processed detection limits a priori. These models, which we provide publicly on a website, explain 70%-80% of the variance in ADI detection limits and 30%-50% of the variance in RDI detection limits. Averaged over a range of angular separations, our models predict both ADI and RDI contrast to within a factor of 2. These results illuminate important factors in high-contrast imaging observations with the NIRC2 vortex coronagraph, help improve observing strategies, and inform future upgrades to the hardware.
The Hubble Space Telescope NICMOS instrument was used from 1997 to 2008 to perform coronagraphic observations of about 400 targets. Most of them were part of surveys looking for substellar companions ...or resolved circumstellar disks to young nearby stars, making the NICMOS coronagraphic archive a valuable database for exoplanets and disks studies. As part of the Archival Legacy Investigations of Circumstellar Environments program, we have consistently reprocessed a large fraction of the NICMOS coronagrahic archive using advanced starlight subtraction methods. We present here the high-level science products of these re-analyzed data, which we delivered back to the community through the Mikulski Archive for Space Telescopes: doi:10.17909/T9W89V. We also present the second version of the HCI-FITS format (for High-Contrast Imaging FITS format), which we developed as a standard format for data exchange of imaging reduced science products. These re-analyzed products are openly available for population statistics studies, characterization of specific targets, or detected point-source identification.
ABSTRACT HD 141569 A is a pre-main sequence B9.5 Ve star surrounded by a prominent and complex circumstellar disk, likely still in a transition stage from protoplanetary to debris disk phase. Here, ...we present a new image of the third inner disk component of HD 141569 A made in the L′ band (3.8 m) during the commissioning of the vector vortex coronagraph that has recently been installed in the near-infrared imager and spectrograph NIRC2 behind the W.M. Keck Observatory Keck II adaptive optics system. We used reference point-spread function subtraction, which reveals the innermost disk component from the inner working distance of 23 au and up to 70 au. The spatial scale of our detection roughly corresponds to the optical and near-infrared scattered light, thermal Q, N, and 8.6 m PAH emission reported earlier. We also see an outward progression in dust location from the L′ band to the H band (Very Large Telescope/SPHERE image) to the visible (Hubble Space Telescope (HST)/STIS image), which is likely indicative of dust blowout. The warm disk component is nested deep inside the two outer belts imaged by HST-NICMOS in 1999 (at 406 and 245 au, respectively). We fit our new L′-band image and spectral energy distribution of HD 141569 A with the radiative transfer code MCFOST. Our best-fit models favor pure olivine grains and are consistent with the composition of the outer belts. While our image shows a putative very faint point-like clump or source embedded in the inner disk, we did not detect any true companion within the gap between the inner disk and the first outer ring, at a sensitivity of a few Jupiter masses.
ABSTRACT We present the first images of four debris disks observed in scattered light around the young (4-250 Myr old) M dwarfs TWA 7 and TWA 25, the K6 star HD 35650, and the G2 star HD 377. We ...obtained these images by reprocessing archival Hubble Space Telescope NICMOS coronagraph data with modern post-processing techniques as part of the Archival Legacy Investigation of Circumstellar Environments program. All four disks appear faint and compact compared with other debris disks resolved in scattered light. The disks around TWA 25, HD 35650, and HD 377 appear very inclined, while TWA 7's disk is viewed nearly face-on. The surface brightness of HD 35650's disk is strongly asymmetric. These new detections raise the number of disks resolved in scattered light around M and late-K stars from one (the AU Mic system) to four. This new sample of resolved disks enables comparative studies of heretofore scarce debris disks around low-mass stars relative to solar-type stars.
Radio and near-infrared observations have observed dozens of protoplanetary disks that host spiral arm features. Numerical simulations have shown that companions may excite spiral density waves in ...protoplanetary disks via companion–disk interaction. However, the lack of direct observational evidence for spiral-driving companions poses challenges to current theories of companion–disk interaction. Here we report multi-epoch observations of the binary system
HD 100453
with the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) facility at the Very Large Telescope. By recovering the spiral features via robustly removing starlight contamination, we measure spiral motion across 4 yr to perform dynamical motion analyses. The spiral pattern motion is consistent with the orbital motion of the eccentric companion. With this first observational evidence of a companion driving a spiral arm among protoplanetary disks, we directly and dynamically confirm the long-standing theory on the origin of spiral features in protoplanetary disks. With the pattern motion of companion-driven spirals being independent of companion mass, here we establish a feasible way of searching for hidden spiral-arm-driving planets that are beyond the detection of existing ground-based high-contrast imagers.
Large-scale spiral arms have been revealed in scattered light images of a few protoplanetary disks. Theoretical models suggest that such arms may be driven by and corotate with giant planets, which ...has called for remarkable observational efforts to look for them. By examining the rotation of the spiral arms for the MWC 758 system over a 10 year timescale, we are able to provide dynamical constraints on the locations of their perturbers. We present reprocessed Hubble Space Telescope (HST)/NICMOS F110W observations of the target in 2005, and the new Keck/NIRC2 L′-band observations in 2017. MWC 758's two well-known spiral arms are revealed in the NICMOS archive at the earliest observational epoch. With additional Very Large Telescope (VLT)/SPHERE data, our joint analysis leads to a pattern speed of at 3 for the two major spiral arms. If the two arms are induced by a perturber on a near-circular orbit, its best-fit orbit is at 89 au (0 59), with a 3 lower limit of 30 au (0 20). This finding is consistent with the simulation prediction of the location of an arm-driving planet for the two major arms in the system.