We present the Vortex Image Processing (VIP) library, a python package dedicated to astronomical high-contrast imaging. Our package relies on the extensive python stack of scientific libraries and ...aims to provide a flexible framework for high-contrast data and image processing. In this paper, we describe the capabilities of VIP related to processing image sequences acquired using the angular differential imaging (ADI) observing technique. VIP implements functionalities for building high-contrast data processing pipelines, encompassing pre- and post-processing algorithms, potential source position and flux estimation, and sensitivity curve generation. Among the reference point-spread function subtraction techniques for ADI post-processing, VIP includes several flavors of principal component analysis (PCA) based algorithms, such as annular PCA and incremental PCA algorithms capable of processing big datacubes (of several gigabytes) on a computer with limited memory. Also, we present a novel ADI algorithm based on non-negative matrix factorization, which comes from the same family of low-rank matrix approximations as PCA and provides fairly similar results. We showcase the ADI capabilities of the VIP library using a deep sequence on HR 8799 taken with the LBTI/LMIRCam and its recently commissioned L-band vortex coronagraph. Using VIP, we investigated the presence of additional companions around HR 8799 and did not find any significant additional point source beyond the four known planets. VIP is available at http://github.com/vortex-exoplanet/VIP and is accompanied with Jupyter notebook tutorials illustrating the main functionalities of the library.
We present the first observational evidence for a circumplanetary disk (CPD) around the protoplanet PDS 70 b, based on a new spectrum in the K-band acquired with Very Large Telescope/SINFONI. We ...tested three hypotheses to explain the spectrum: atmospheric emission from the planet with either (1) a single value of extinction, (2) a variable extinction, and (3) a combined atmospheric and CPD model. Goodness-of-fit indicators favor the third option, suggesting that circumplanetary material contributes excess thermal emission-most prominent at λ 2.3 m. Inferred accretion rates (∼10−7.8-10−7.3 MJ yr−1) are compatible with observational constraints based on the H and Brγ lines. For the planet, we derive an effective temperature of 1500-1600 K, surface gravity , radius ∼1.6RJ, mass ∼10MJ, and possible thick clouds. Models with variable extinction lead to slightly worse fits. However, the amplitude (ΔAV 3 mag) and timescale of variation ( years) required for the extinction would also suggest circumplanetary material.
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.
ABSTRACT
Planet–disc interactions build up local pressure maxima that may halt the radial drift of protoplanetary dust, and pile it up in rings and crescents. ALMA observations of the HD 135344B disc ...revealed two rings in the thermal continuum stemming from ∼mm-sized dust. At higher frequencies the inner ring is brighter relative to the outer ring, which is also shaped as a crescent rather than a full ring. In near-IR scattered light images, the disc is modulated by a two-armed grand-design spiral originating inside the ALMA inner ring. Such structures may be induced by a massive companion evacuating the central cavity, and by a giant planet in the gap separating both rings, that channels the accretion of small dust and gas through its filamentary wakes while stopping the larger dust from crossing the gap. Here we present ALMA observations in the J = (2 − 1) CO isotopologue lines and in the adjacent continuum, with up to 12 km baselines. Angular resolutions of ∼0${_{.}^{\prime\prime}}$03 reveal the tentative detection of a filament connecting both rings, and which coincides with a local discontinuity in the pitch angle of the IR spiral, proposed previously as the location of the protoplanet driving this spiral. Line diagnostics suggests that turbulence, or superposed velocity components, is particularly strong in the spirals. The 12CO(2-1) 3D rotation curve points at stellocentric accretion at radii within the inner dust ring, with a radial velocity of up to ${\sim}5{{\ \rm per\ cent}}\pm 0.5{{\ \rm per\ cent}}$ Keplerian, which corresponds to an excessively large accretion rate of ${\sim}2\times 10^{-6}\, M_\odot \,$yr−1 if all of the CO layer follows the 12CO(2-1) kinematics. This suggests that only the surface layers of the disc are undergoing accretion, and that the line broadening is due to superposed laminar flows.
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
Over the past decade, hundreds of nights have been spent on the world’s largest telescopes to search for and directly detect new exoplanets using high-contrast imaging (HCI). Thereby, two ...scientific goals are of central interest: first, to study the characteristics of the underlying planet population and distinguish between different planet formation and evolution theories. Second, to find and characterize planets in our immediate solar neighborhood. Both goals heavily rely on the metric used to quantify planet detections and nondetections. Current standards often rely on several explicit or implicit assumptions about noise. For example, it is often assumed that the residual noise after data postprocessing is Gaussian. While being an inseparable part of the metric, these assumptions are rarely verified. This is problematic as any violation of these assumptions can lead to systematic biases. This makes it hard, if not impossible, to compare results across data sets or instruments with different noise characteristics. We revisit the fundamental question of how to quantify detection limits in HCI. We focus our analysis on the error budget resulting from violated assumptions. To this end, we propose a new metric based on bootstrapping that generalizes current standards to non-Gaussian noise. We apply our method to archival HCI data from the NACO instrument at the Very Large Telescope and derive detection limits for different types of noise. Our analysis shows that current standards tend to give detection limits that are about one magnitude too optimistic in the speckle-dominated regime. That is, HCI surveys may have excluded planets that can still exist.
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.
Abstract Recent Atacama Large Millimeter/submillimeter Array (ALMA) observations of protoplanetary disks in the millimeter continuum have shown a variety of radial gaps, cavities, and spiral ...features. These substructures may be signposts for ongoing planet formation, and therefore these systems are promising targets for direct imaging planet searches in the near-infrared. To this end, we present results from a deep imaging survey in the L ′ band (3.8 μ m) with the Keck/NIRC2 vortex coronagraph to search for young planets in 43 disks with resolved features in the millimeter continuum or evidence for gaps/central cavities from their spectral energy distributions. Although we do not detect any new point sources, using the vortex coronagraph allows for high sensitivity to faint sources at small angular separations (down to ∼0.″1), allowing us to place strong upper limits on the masses of potential gas giant planets. We compare our mass sensitivities to the masses of planets derived using ALMA observations, and while we are sensitive to ∼1 M Jup planets in the gaps in some of our systems, we are generally not sensitive to planets of the masses expected from the ALMA observations. In addition to placing upper limits on the masses of gas giant planets that could be interacting with the dust in the disks to form the observed millimeter substructures, we are also able to map the micron-sized dust as seen in scattered light for 8 of these systems. Our large sample of systems also allows us to investigate limits on planetary accretion rates and disk viscosities.
We present L′-band Keck/NIRC2 imaging and H-band Subaru/AO188+HiCIAO polarimetric observations of the CQ Tau disk with a new spiral arm. Apart from the spiral feature, our observations could not ...detect any companion candidates. We traced the spiral feature from the r2-scaled High-Contrast Coronographic Imager for Adaptive Optics (HiCIAO) polarimetric intensity image and the fitted result is used for forward modeling to reproduce the ADI-reduced NIRC2 image. We estimated the original surface brightness after throughput correction in the L′ band to be ∼126 mJy arcsec−2 at most. We suggest that the grain temperature of the spiral may be heated up to ∼200 K in order to explain both of the H- and L′-band results. The H-band emission at the location of the spiral originates from the scattering from the disk surface while both scattering and thermal emission may contribute to the L′-band emission. If the central star is only the light source of scattered light, the spiral emission at the L′ band should be thermal emission. If an inner disk also acts as the light source, the scattered light and the thermal emission may equally contribute to the L′-band spiral structure.
ABSTRACT Debris dust in the habitable zones of stars-otherwise known as exozodiacal dust-comes from extrasolar asteroids and comets and is thus an expected part of a planetary system. Background flux ...from the solar system's zodiacal dust and the exozodiacal dust in the target system is likely to be the largest source of astrophysical noise in direct observations of terrestrial planets in the habitable zones of nearby stars. Furthermore, dust structures like clumps, thought to be produced by dynamical interactions with exoplanets, are a possible source of confusion. In this article, we qualitatively assess the primary impact of exozodiacal dust on high-contrast direct imaging at optical wavelengths, such as would be performed with a coronagraph. Then we present the sensitivity of previous, current, and near-term facilities to thermal emission from debris dust at all distances from nearby solar-type stars, as well as our current knowledge of dust levels from recent surveys. Finally, we address the other method of detecting debris dust, through high-contrast imaging in scattered light. This method is currently far less sensitive than thermal emission observations, but provides high spatial resolution for studying dust structures. This article represents the first report of NASA's Exoplanet Exploration Program Analysis Group (ExoPAG).