We present the Vortex Image Processing (VIP) library, a python package dedicated to astronomical highcontrast 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 postprocessing 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.

Context. Transition disks offer the extraordinary opportunity to look for newly born planets and investigate the early stages of planet formation. Aims. In this context we observed the Herbig A5 star MWC 758 with the L'-band vector vortex coronagraph installed in the nearinfrared camera and spectrograph NIRC2 at the Keck II telescope, with the aim of unveiling the nature of the spiral structure by constraining the presence of planetary companions in the system. Methods. Our high-contrast imaging observations show a bright (∆L = 7.0 ± 0.3 mag) point-like emission, south of MWC 758 at a deprojected separation of ∼20 au (r = 0 . 111 ± 0 . 004) from the central star. We also recover the two spiral arms (south-east and north-west), already imaged by previous studies in polarized light, and discover a third one to the south-west of the star. No additional companions were detected in the system down to 5 Jupiter masses beyond 0 . 6 from the star. Results. We propose that the bright L'-band emission could be caused by the presence of an embedded and accreting protoplanet, although the possibility of it being an asymmetric disk feature cannot be excluded. The spiral structure is probably not related to the protoplanet candidate, unless on an inclined and eccentric orbit, and it could be due to one (or more) yet undetected planetary companions at the edge of or outside the spiral pattern. Future observations and additional simulations will be needed to shed light on the true nature of the point-like source and its link with the spiral arms.

Context. The direct detection of exoplanets with high-contrast imaging requires advanced data processing methods to disentangle potential planetary signals from bright quasi-static speckles. Among them, angular differential imaging (ADI) permits potential planetary signals with a known rotation rate to be separated from instrumental speckles that are either statics or slowly variable. The method presented in this paper, called ANDROMEDA for ANgular Differential OptiMal Exoplanet Detection Algorithm is based on a maximum likelihood approach to ADI and is used to estimate the position and the flux of any point source present in the field of view. Aims. In order to optimize and experimentally validate this previously proposed method, we applied ANDROMEDA to real VLT/NaCo data. In addition to its pure detection capability, we investigated the possibility of defining simple and efficient criteria for automatic point source extraction able to support the processing of large surveys. Methods. To assess the performance of the method, we applied ANDROMEDA on VLT/NaCo data of TYC-8979-1683-1 which is surrounded by numerous bright stars and on which we added synthetic planets of known position and flux in the field. In order to accommodate the real data properties, it was necessary to develop additional pre-processing and post-processing steps to the initially proposed algorithm. We then investigated its skill in the challenging case of a well-known target, β Pictoris, whose companion is close to the detection limit and we compared our results to those obtained by another method based on principal component analysis (PCA). Results. Application on VLT/NaCo data demonstrates the ability of ANDROMEDA to automatically detect and characterize point sources present in the image field. We end up with a robust method bringing consistent results with a sensitivity similar to the recently published algorithms, with only two parameters to be fine tuned. Moreover, the companion flux estimates are not biased by the algorithm parameters and do not require a posteriori corrections. Conclusions. ANDROMEDA is an attractive alternative to current standard image processing methods that can be readily applied to on-sky data.

An optical vortex coronagraph has been implemented within the NIRC2 camera on the Keck II telescope and used to carry out on-sky tests and observations. The development of this new L'-band observational mode is described, and an initial demonstration of the new capability is presented: a resolved image of the low-mass companion to HIP79124, which had previously been detected by means of interferometry. With HIP79124 B at a projected separation of 186.5 mas, both the small inner working angle of the vortex coronagraph and the related imaging improvements were crucial in imaging this close companion directly. Due to higher Strehl ratios and more relaxed contrasts in L' band versus H band, this new coronagraphic capability will enable high-contrast small-angle observations of nearby young exoplanets and disks on a par with those of shorterwavelength extreme adaptive optics coronagraphs.

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 frames 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.