Point Spread Function Subtraction Research Articles

MINDS: JWST/NIRCam imaging of the protoplanetary disk PDS 70

Context. Two protoplanets have recently been discovered within the PDS 70 protoplanetary disk. JWST/NIRCam offers a unique opportunity to characterize them and their birth environment at wavelengths that are difficult to access from the ground. Aims. We image the circumstellar environment of PDS 70 at 1.87 μm and 4.83 μm, assess the presence of Pa-α emission due to accretion onto the protoplanets, and probe any IR excess indicative of heated circumplanetary material. Methods. We obtained noncoronagraphic JWST/NIRCam images of PDS 70 within the MIRI mid-INfrared Disk Survey (MINDS) program. We leveraged the Vortex Image Processing (VIP) package for data reduction, and we developed dedicated routines for optimal stellar point spread function subtraction, unbiased imaging of the disk, and protoplanet flux measurement in this type of dataset. A radiative transfer model of the disk was used to separate the contributions from the disk and the protoplanets. Results. We redetect both protoplanets and identify extended emission after subtracting a disk model, including a large-scale spiral-like feature. We interpret its signal in the direct vicinity of planet c as tracing the accretion stream that feeds its circumplanetary disk, while the outer part of the feature may rather reflect asymmetric illumination of the outer disk. We also report a bright signal that is consistent with a previously proposed protoplanet candidate enshrouded in dust near the 1:2:4 mean-motion resonance with planets b and c. The 1.87 μm flux of planet b is consistent with atmospheric model predictions, but the flux of planet c is not. We discuss potential origins for this discrepancy, including significant Pa-α line emission. The 4.83 μm fluxes of planets b and c suggest enshrouding dust or heated CO emission from their circumplanetary environment. Conclusions. The use of image-processing methods that are optimized for extended disk signals on high-sensitivity and high-stability from JWST can uniquely identify signatures of planet–disk interactions and enable accurate photometry of protoplanets at wavelengths that are difficult to probe from the ground. Our results indicate that more protoplanets can be identified and characterized in other JWST datasets.

JWST MIRI MRS Observations of T Cha: Discovery of a Spatially Resolved Disk Wind

Understanding when and how circumstellar disks disperse is crucial to constrain planet formation and migration. Thermal winds powered by high-energy stellar photons have long been theorized to drive disk dispersal. However, evidence for these winds is currently based only on small (∼3–6 km s−1) blueshifts in [Ne ii] 12.81 μm lines, which does not exclude MHD winds. We report JWST MIRI MRS spectro-imaging of T Cha, a disk with a large dust gap (∼30 au in radius) and blueshifted [Ne ii] emission. We detect four forbidden noble gas lines, [Ar ii], [Ar iii], [Ne ii], and [Ne iii], of which [Ar iii] is the first detection in any protoplanetary disk. We use line flux ratios to constrain the energy of the ionizing photons and find that argon is ionized by extreme ultraviolet, whereas neon is most likely ionized by X-rays. After performing continuum and point-spread function subtraction on the integral field unit cube, we discover a spatial extension in the [Ne ii] emission off the disk continuum emission. This is the first spatially resolved [Ne ii] disk wind emission. The mostly ionic spectrum of T Cha, in combination with the extended [Ne ii] emission, points to an evolved stage for any inner MHD wind and is consistent with the existence of an outer thermal wind ionized and driven by high-energy stellar photons. This work acts as a pathfinder for future observations aiming at investigating disk dispersal using JWST.

Probing Disk Ice Content and Polycyclic Aromatic Hydrocarbon Emission through Multiband MagAO+Clio Images of HD 141569

We present resolved images of the inner disk component around HD 141569A using the Magellan adaptive optics system with the Clio2 1–5 μm camera, offering a glimpse of a complex system thought to be in a short evolutionary phase between protoplanetary and debris disk stages. We use a reference star along with the Karhunen–Loéve image projection (KLIP) algorithm for point-spread function subtraction to detect the disk inward to about 0.″24 (∼25 au assuming a distance of 111 pc) at high signal-to-noise ratios at L′ (3.8 μm), Ls (3.3 μm), and narrowband Ice (3.1 μm). We identify an arc or spiral arm structure at the southeast extremity, consistent with previous studies. We implement forward modeling with a simple disk model within the framework of a Markov Chain Monte Carlo sampler to better constrain the geometrical attributes and photometry using our KLIP-reduced disk images. We then leverage these modeling results to facilitate a comparison of the measured brightness in each passband to find a reduction in scattered light from the disk in the Ice filter, implying significant absorption due to water ice in the dust. Additionally, our best-fit disk models exhibit peak brightness in the southwestern, back-scattering region of the disk, which we suggest to be possible evidence of 3.3 μm polycyclic aromatic hydrocarbon emission. However, we point out the need for additional observations with bluer filters and more complex modeling to confirm these hypotheses.

An imaged 15 MJup companion within a hierarchical quadruple system

Context. Since 2019, the direct imaging B-star Exoplanet Abundance STudy (BEAST) at SPHERE@VLT has been scanning the surroundings of young B-type stars in order to ascertain the ultimate frontiers of giant planet formation. Recently, the 17−4+3 Myr HIP 81208 was found to host a close-in (∼50 au) brown dwarf and a wider (∼230 au) late M star around the central 2.6 M⊙ primary. Aims. Alongside the continuation of the survey, we are undertaking a complete reanalysis of archival data aimed at improving detection performances so as to uncover additional low-mass companions. Methods. We present here a new reduction of the observations of HIP 81208 using the patch covariance algorithm (PACO), a recent and powerful algorithm dedicated to processing high-contrast imaging datasets, as well as more classical algorithms and a dedicated point spread function subtraction approach. The combination of different techniques allowed for a reliable extraction of astrometric and photometric parameters. Results. A previously undetected source was recovered at a short separation from the C component of the system. Proper motion analysis provided robust evidence for the gravitational bond of the object to HIP 81208 C. Orbiting C at a distance of ∼20 au, this 15 MJup brown dwarf becomes the fourth object of the hierarchical HIP 81208 system. Conclusions. Among the several BEAST stars which are being found to host substellar companions, HIP 81208 stands out as a particularly striking system. As the first stellar binary system with substellar companions around each component ever found by direct imaging, it yields exquisite opportunities for thorough formation and dynamical follow-up studies.

Speckle Space–Time Covariance in High-contrast Imaging

We introduce a new framework for point-spread function subtraction based on the spatiotemporal variation of speckle noise in high-contrast imaging data where the sampling timescale is faster than the speckle evolution timescale. One way that space–time covariance arises in the pupil is as atmospheric layers translate across the telescope aperture and create small, time-varying perturbations in the phase of the incoming wavefront. The propagation of this field to the focal plane preserves some of that space–time covariance. To utilize this covariance, our new approach uses a Karhunen–Loève transform on an image sequence, as opposed to a set of single reference images as in previous applications of Karhunen–Loève Image Processing (KLIP) for high-contrast imaging. With the recent development of photon-counting detectors, such as microwave kinetic inductance detectors, this technique now has the potential to improve contrast when used as a post-processing step. Preliminary testing on simulated data shows this technique can improve contrast by at least 10%–20% from the original image, with significant potential for further improvement. For certain choices of parameters, this algorithm may provide larger contrast gains than spatial-only KLIP.

The Giant Accreting Protoplanet Survey (GAPlanetS): Optimization Techniques for Robust Detections of Protoplanets

High-contrast imaging has afforded astronomers the opportunity to study light directly emitted by adolescent (tens of megayears) and “proto” (<10 Myr) planets still undergoing formation. Direct detection of these planets is enabled by empirical point-spread function (PSF) modeling and removal algorithms. The computational intensity of such algorithms, as well as their multiplicity of tunable input parameters, has led to the prevalence of ad hoc optimization approaches to high-contrast imaging results. In this work, we present a new, systematic approach to optimization vetted using data of the high-contrast stellar companion HD 142527 B from the Magellan Adaptive Optics Giant Accreting Protoplanet Survey (GAPlanetS). More specifically, we present a grid search technique designed to explore three influential parameters of the PSF subtraction algorithm pyKLIP: annuli, movement, and KL modes. We consider multiple metrics for postprocessed image quality in order to optimally recover at Hα (656 nm) synthetic planets injected into contemporaneous continuum (643 nm) images. These metrics include peak (single-pixel) signal-to-noise ratio (S/N), average (multipixel average) S/N, 5σ contrast, and false-positive fraction. We apply continuum-optimized KLIP reduction parameters to six Hα direct detections of the low-mass stellar companion HD 142527 B and recover the companion at a range of separations. Relative to a single-informed, nonoptimized set of KLIP parameters applied to all data sets uniformly, our multimetric grid search optimization led to improvements in companion S/N of up to 1.2σ, with an average improvement of 0.6σ. Since many direct imaging detections lie close to the canonical 5σ threshold, even such modest improvements may result in higher yields in future imaging surveys.

ISPY: NACO Imaging Survey for Planets around Young stars

Context. Planet formation is a frequent process, but little observational constraints exist about the mechanisms involved, especially for giant planets at large separation. The NaCo-ISPY large program is a 120 night L′-band direct imaging survey aimed at investigating the giant planet population on wide orbits (a &gt; 10 au) around stars hosting disks. Aims. Here we present the statistical analysis of a subsample of 45 young stars surrounded by protoplanetary disks (PPDs). This is the largest imaging survey uniquely focused on PPDs to date. Our goal is to search for young forming companions embedded in the disk material and to constrain their occurrence rate in relation to the formation mechanism. Methods. We used principal component analysis based point spread function subtraction techniques to reveal young companions forming in the disks. We calculated detection limits for our datasets and adopted a black-body model to derive temperature upper limits of potential forming planets. We then used Monte Carlo simulations to constrain the population of forming gas giant companions and compare our results to different types of formation scenarios. Results. Our data revealed a new binary system (HD 38120) and a recently identified triple system with a brown dwarf companion orbiting a binary system (HD 101412), in addition to 12 known companions. Furthermore, we detected signals from 17 disks, two of which (HD 72106 and T CrA) were imaged for the first time. We reached median detection limits of L′ = 15.4 mag at 2″.0, which were used to investigate the temperature of potentially embedded forming companions. We can constrain the occurrence of forming planets with semi-major axis a in [20–500] au and Teff in [600–3000] K to be 21.2-13.6+24.3%, 14.8-9.6+17.5%, and 10.8-7.0+12.6% for Rp = 2, 3, 5 RJ, which is in line with the statistical results obtained for more evolved systems from other direct imaging surveys. These values are obtained under the assumption that extinction from circumstellar and circumplanetary material does not affect the companion signal, but we show the potential impact these factors might have on the detectability of forming objects. Conclusions. The NaCo-ISPY data confirm that massive bright planets accreting at high rates are rare. More powerful instruments with better sensitivity in the near- to mid-infrared are likely required to unveil the wealth of forming planets sculpting the observed disk substructures.

NICMOS Kernel-phase Interferometry. I. Catalogue of Brown Dwarfs Observed in F110W and F170M

Filling out the dearth of detections between direct imaging and radial velocity surveys will test theories of planet formation and (sub)stellar binarity across the full range of semimajor axes, connecting formation of close- to wide-separation gas giants and substellar companions. Direct detection of close-in companions is notoriously difficult: coronagraphs and point-spread function subtraction techniques fail near the λ/D diffraction limit. We present a new faint companion detection pipeline called Argus that analyzes kernel phases, an interferometric observable analogous to closure phases from nonredundant aperture masking but utilizing the full unobscured telescope aperture. We demonstrate the pipeline, and the power of interferometry, by performing a companion search on the entire HST/NICMOS F110W and F170M image archive of 114 nearby brown dwarfs (observed in seven different programs). Our pipeline is able to detect companions down to flux ratios of ∼102 at half the classical diffraction limit. We discover no new companions, but recover and refine astrometry of 19 previous imaging companions (two with multiple epochs) and confirm two previous kernel-phase detections. We discuss the limitations of this technique with respect to nondetections of previously confirmed or proposed companions. We present contrast curves to enable population studies to leverage nondetections and to demonstrate the strength of this technique at separations inaccessible to classical imaging techniques. The binary fraction of our sample (%) is consistent with previous binary surveys, even with sensitivity to much tighter separation companions.

Hubble Space Telescope Reveals Spectacular Light Echoes Associated with the Stripped-envelope Supernova 2016adj in the Iconic Dust Lane of Centaurus A

We present a multiband sequence of Hubble Space Telescope images documenting the emergence and evolution of multiple light echoes (LEs) linked to the stripped-envelope supernova (SN) 2016adj located in the central dust lane of Centaurus A. Following point-spread function subtraction, we identify the earliest LE emission associated with an SN at only +34 days past the epoch of the B-band maximum. Additional HST images extending through +578 days cover the evolution of LE1 taking the form of a ring, while images taken on +1991 days reveal not only LE1 but also segments of a new inner LE ring (LE2) as well as two additional outer LE rings (LE3 and LE4). Adopting the single-scattering formalism, the angular radii of the LEs suggest they originate from discrete dust sheets in the foreground of the SN. This information, combined with measurements of color and optical depth of the scattering surfaces, informs a scenario with multiple sheets of clumpy dust characterized by a varying degree of holes. In this case, the larger the LE’s angular radii, the farther in the foreground of the SN its dust sheet is located. However, an exception to this is LE2, which is formed by a dust sheet located in closer proximity to the SN than the dust sheets producing LE1, LE3, and LE4. The delayed appearance of LE2 can be attributed to its dust sheet having a significant hole along the line of sight between the SN and Earth.

StraKLIP: A Novel Pipeline for Detection and Characterization of Close-in Faint Companions through the Karhunen–Loéve Image Processing Algorithm

We present a new pipeline developed to detect and characterize faint astronomical companions at small angular separation from the host star using sets of wide-field imaging observations not specifically designed for high-contrast imaging analysis. The core of the pipeline relies on Karhunen–Loéve truncated transformation of the reference point-spread function (PSF) library to perform PSF subtraction and identify candidates. Tests of reliability of detections and characterization of companions are made through simulation of binaries and generation of receiver operating characteristic curves for false-positive/true-positive analysis. The algorithm has been successfully tested on large HST/ACS and WFC3 data sets acquired for two HST Treasury Programs on the Orion Nebula Cluster. Based on these extensive numerical experiments we find that, despite being based on methods designed for observations of a single star at a time, our pipeline performs very well on mosaic space-based data. In fact, we are able to detect brown-dwarf-mass companions almost down to the planetary-mass limit. The pipeline is able to reliably detect signals at separations as close as ≳0.″1 with a completeness of ≳10%, or ∼0.″2 with a completeness of ∼30%. This approach can potentially be applied to a wide variety of space-based imaging surveys, from data in the existing HST archive to near-future JWST mosaics and future wide-field Roman images.

Companion mass limits for 17 binary systems obtained with binary differential imaging and MagAO/Clio

ABSTRACT Improving direct detection capability close to the star through improved star subtraction and post-processing techniques is vital for discovering new low-mass companions and characterizing known ones at longer wavelengths. We present results of 17 binary star systems observed with the Magellan adaptive optics system (MagAO) and the Clio infrared camera on the Magellan Clay Telescope using binary differential imaging (BDI). BDI is an application of reference differential imaging (RDI) and angular differential imaging (ADI) applied to wide binary star systems (2 arcsec &amp;lt;Δρ &amp;lt; 10 arcsec) within the isoplanatic patch in the infrared. Each star serves as the point spread function (PSF) reference for the other, and we performed PSF estimation and subtraction using principal component analysis. We report contrast and mass limits for the 35 stars in our initial survey using BDI with MagAO/Clio in L′ and 3.95 µm bands. Our achieved contrasts varied between systems, and spanned a range of contrasts from 3.0 to 7.5 magnitudes and a range of separations from 0.2 to 2 arcsec. Stars in our survey span a range of masses, and our achieved contrasts correspond to late-type M-dwarf masses down to ∼10 MJup. We also report detection of a candidate companion signal at 0.2 arcsec (18 au) around HIP 67506 A (SpT G5V, mass ∼1.2 M⊙), which we estimate to be $\sim 60-90 \, \rm{M_{Jup}}$. We found that the effectiveness of BDI is highest for approximately equal brightness binaries in high-Strehl conditions.

Detection of Near-infrared Water Ice at the Surface of the (Pre)Transitional Disk of AB Aur: Informing Icy Grain Abundance, Composition, and Size

Abstract We present near-infrared Large Binocular Telescope LMIRCam imagery of the disk around the Herbig Ae/Be star AB Aurigae. A comparison of the surface brightness at K s (2.16 μm), H2O narrowband (3.08 μm), and L ′ (3.7 μm) allows us to probe the presence of icy grains in this (pre)transitional disk environment. By applying reference differential imaging point-spread function subtraction, we detect the disk at high signal-to-noise ratios in all three bands. We find strong morphological differences between the bands, including asymmetries consistent with the observed spiral arms within 100 au in L ′ . An apparent deficit of scattered light at 3.08 μm relative to the bracketing wavelengths (K s and L ′ ) is evocative of ice absorption at the disk surface layer. However, the Δ(K s − H2O) color is consistent with grains with little to no ice (0%–5% by mass). The Δ ( H 2 O − L ′ ) color, conversely, suggests grains with a much higher ice mass fraction (∼0.68), and the two colors cannot be reconciled under a single grain population model. Additionally, we find that the extremely red Δ ( K s − L ′ ) disk color cannot be reproduced under conventional scattered light modeling with any combination of grain parameters or reasonable local extinction values. We hypothesize that the scattering surfaces at the three wavelengths are not colocated, and that the optical depth effects in each wavelength result from probing the grain population at different disk surface depths. The morphological similarity between K s and H2O suggests that their scattering surfaces are near one another, lending credence to the Δ(K s − H2O) disk color constraint of &lt;5% ice mass fraction for the outermost scattering disk layer.

Circumnuclear Dust in AP Librae and the Source of Its VHE Emission

Abstract The broad high-energy spectral component in blazars is usually attributed to various inverse Compton scattering processes in the relativistic jet, but has not been clearly identified in most cases due to degeneracies in physical models. AP Librae, a low-synchrotron-peaking BL Lac object (LBL) detected in 2015 by H.E.S.S. at very high energies (VHE; &gt;0.5 TeV), has an extremely broad high-energy spectrum, covering ∼9 decades in energy. Standard synchrotron self-Compton models generally fail to reproduce the VHE emission, which has led to the suggestion that it might arise not from the blazar core, but on kiloparsec scales from inverse Compton (IC) scattering of cosmic microwave background (CMB) photons by a still-relativistic jet (IC/CMB). IC/CMB models for the TeV emission of AP Librae in prior works have implied a high level of infrared emission from the kiloparsec-scale jet. With newly obtained Hubble Space Telescope (HST) imaging, we obtain a deep upper limit on the kiloparsec-scale jet emission at 1.6 μm, well below the expected level. High-resolution Atacama Large Millimeter/submillimeter Array imaging in bands 3–9 reveals a residual dust-disk signature after core subtraction, with a clearly thermal spectrum, and an extent (∼500 pc) that matches with a nonjet residual emission seen after point-spread function subtraction in our 1.6 μm HST imaging. We find that the unusually broad GeV and VHE emission in AP Librae can be reproduced through the combined IC scattering of photons from the CMB and the dust disk, respectively, by electrons in both the blazar core and subkiloparsec jet.

Auto-RSM: An automated parameter-selection algorithm for the RSM map exoplanet detection algorithm

Context. Most of the high-contrast imaging (HCI) data-processing techniques used over the last 15 years have relied on the angular differential imaging (ADI) observing strategy, along with subtraction of a reference point spread function (PSF) to generate exoplanet detection maps. Recently, a new algorithm called regime switching model (RSM) map has been proposed to take advantage of these numerous PSF-subtraction techniques; RSM uses several of these techniques to generate a single probability map. Selection of the optimal parameters for these PSF-subtraction techniques as well as for the RSM map is not straightforward, is time consuming, and can be biased by assumptions made as to the underlying data set. Aims. We propose a novel optimisation procedure that can be applied to each of the PSF-subtraction techniques alone, or to the entire RSM framework. Methods. The optimisation procedure consists of three main steps: (i) definition of the optimal set of parameters for the PSF-subtraction techniques using the contrast as performance metric, (ii) optimisation of the RSM algorithm, and (iii) selection of the optimal set of PSF-subtraction techniques and ADI sequences used to generate the final RSM probability map. Results. The optimisation procedure is applied to the data sets of the exoplanet imaging data challenge, which provides tools to compare the performance of HCI data-processing techniques. The data sets consist of ADI sequences obtained with three state-of-the-art HCI instruments: SPHERE, NIRC2, and LMIRCam. The results of our analysis demonstrate the interest of the proposed optimisation procedure, with better performance metrics compared to the earlier version of RSM, as well as to other HCI data-processing techniques.

Characterizing the Discovery of a New Trans-Neptunian Object Binary in a Trailed Point-spread Function Search

Abstract The Latitude Density Search utilized the Hyper Suprime-Cam on the Subaru Telescope to discover 60 moving objects in the outer solar system, 54 of which have semimajor axes beyond 30 au. The images were acquired in exceptional seeing (0.″4) and reached a detection limit of m r ≃ 25.2. The two night arcs were used to calculate orbits that are poorly constrained; however, the distance and inclination are the parameters best constrained by short arcs, and a reasonable determination can be made of which objects are cold classical trans-Neptunian objects (TNOs) and which are dynamically excited. We identify 10 objects as likely cold classical objects. We searched all of the detections for binary sources using a trailed point-spread function (PSF) subtraction method and identified one binary object with a separation of 0.″34 and a secondary with 17% the brightness of the primary (2.0 magnitudes fainter). This is the brightest TNO in the sample, the previously known object (471165) 2010 HE79, which has a dynamically excited orbit. Because of the excellent seeing, this search was sensitive to secondaries with 0.″34 separation and a brightness of ≥50% the primary brightness for seven objects, including one cold classical. This gives an intrinsic binary fraction of 17 − 10 + 19 % (one of six) for the dynamically excited objects within our detection limits. The trailed PSF subtraction method used in the Latitude Density Search to identify binaries, fit the two components, and determine the sensitivity limits is a useful tool that could be more broadly applied to identify binary TNOs and track known binary TNO orbits.

SCExAO/MEC and CHARIS Discovery of a Low-mass, 6 au Separation Companion to HIP 109427 Using Stochastic Speckle Discrimination and High-contrast Spectroscopy*

Abstract We report the direct imaging discovery of a low-mass companion to the nearby accelerating A star, HIP 109427, with the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument coupled with the Microwave Kinetic Inductance Detector Exoplanet Camera (MEC) and CHARIS integral field spectrograph. CHARIS data reduced with reference star point spread function (PSF) subtraction yield 1.1–2.4 μm spectra. MEC reveals the companion in Y and J band at a comparable signal-to-noise ratio using stochastic speckle discrimination, with no PSF subtraction techniques. Combined with complementary follow-up L p photometry from Keck/NIRC2, the SCExAO data favors a spectral type, effective temperature, and luminosity of M4–M5.5, 3000–3200 K, and log 10 ( L / L ⊙ ) = − 2.28 − 0.04 + 0.04 , respectively. Relative astrometry of HIP 109427 B from SCExAO/CHARIS and Keck/NIRC2, and complementary Gaia–Hipparcos absolute astrometry of the primary favor a semimajor axis of 6.55+3.0 −0.48 au, an eccentricity of 0.54 − 0.15 + 0.28 , an inclination of 66.7 − 14 + 8.5 degrees, and a dynamical mass of 0.280 − 0.059 + 0.18 M ⊙. This work shows the potential for extreme AO systems to utilize speckle statistics in addition to widely used postprocessing methods to directly image faint companions to nearby stars near the telescope diffraction limit.

Absorption-selected galaxies trace the low-mass, late-type, star-forming population at z ∼ 2–3

ABSTRACT We report on the stellar content, half-light radii and star formation rates of a sample of 10 known high-redshift (z ≳ 2) galaxies selected on strong neutral hydrogen (H i) absorption ($\log ({\rm N_{H\, \rm {I}}\: /\: cm}^{-2})\,\gt\, 19$) towards background quasars. We use observations from the Hubble Space Telescope Wide Field Camera 3 in three broad-band filters to study the spectral energy distribution (SED) of the galaxies. Using careful quasar point spread function subtraction, we study their galactic environments, and perform the first systematic morphological characterization of such absorption-selected galaxies at high redshifts. Our analysis reveals complex, irregular hosts with multiple star-forming clumps. At a spatial sampling of 0.067 arcsec per pixel (corresponding to 0.55 kpc at the median redshift of our sample), 40 per cent of our sample requires multiple Sérsic components for an accurate modelling of the observed light distributions. Placed on the mass–size relation and the ‘main sequence’ of star-forming galaxies, we find that absorption-selected galaxies at high redshift extend known relations determined from deep luminosity-selected surveys to an order of magnitude lower stellar mass, with objects primarily composed of star-forming, late-type galaxies. We measure half-light radii in the range r1/2 ∼ 0.4 to 2.6 kpc based on the reddest band (F160W) to trace the oldest stellar populations, and stellar masses in the range log (M⋆/M⊙) ∼ 8 to 10 derived from fits to the broad-band SED. Spectroscopic and SED-based star formation rates are broadly consistent, and lie in the range $\log (\mathrm{SFR}/{\rm M}_{\odot }\, {\rm yr}^{-1}) \sim 0.0$ to 1.7.

Discovery of a directly imaged planet to the young solar analog YSES 2

Context. To understand the origin and formation pathway of wide-orbit gas giant planets, it is necessary to expand the limited sample of these objects. The mass of exoplanets derived with spectrophotometry, however, varies strongly as a function of the age of the system and the mass of the primary star. Aims. By selecting stars with similar ages and masses, the Young Suns Exoplanet Survey (YSES) aims to detect and characterize planetary-mass companions to solar-type host stars in the Scorpius-Centaurus association. Methods. Our survey is carried out with VLT/SPHERE with short exposure sequences on the order of 5 min per star per filter. The subtraction of the stellar point spread function (PSF) is based on reference star differential imaging using the other targets (with similar colors and magnitudes) in the survey in combination with principal component analysis. Two astrometric epochs that are separated by more than one year are used to confirm co-moving companions by proper motion analysis. Results. We report the discovery of YSES 2b, a co-moving, planetary-mass companion to the K1 star YSES 2 (TYC 8984-2245-1, 2MASS J11275535-6626046). The primary has a Gaia EDR3 distance of 110 pc, and we derive a revised mass of 1.1 M⊙ and an age of approximately 14 Myr. We detect the companion in two observing epochs southwest of the star at a position angle of 205° and with a separation of ~1.′′05, which translates to a minimum physical separation of 115 au at the distance of the system. Photometric measurements in the H and Ks bands are indicative of a late L spectral type, similar to the innermost planets around HR 8799. We derive a photometric planet mass of 6.3−0.9+1.6 MJup using AMES-COND and AMES-dusty evolutionary models; this mass corresponds to a mass ratio of q = (0.5 ± 0.1)% with the primary. This is the lowest mass ratio of a direct imaging planet around a solar-type star to date. We discuss potential formation mechanisms and find that the current position of the planet is compatible with formation by disk gravitational instability, but its mass is lower than expected from numerical simulations. Formation via core accretion must have occurred closer to the star, yet we do not find evidence that supports the required outward migration, such as via scattering off another undiscovered companion in the system. We can exclude additional companions with masses greater than 13 MJup in the full field of view of the detector (0.′′15&lt;ρ&lt;5.′′50), at 0.′′5 we can rule out further objects that are more massive than 6 MJup, and for projected separations ρ &gt;2′′ we are sensitive to planets with masses as low as 2 MJup. Conclusions. YSES 2b is an ideal target for follow-up observations to further the understanding of the physical and chemical formation mechanisms of wide-orbit Jovian planets. The YSES strategy of short snapshot observations (≤5 min) and PSF subtraction based on a large reference library proves to be extremely efficient and should be considered for future direct imaging surveys.

New mid-infrared imaging constraints on companions and protoplanetary disks around six young stars

Context. Mid-infrared (mid-IR) imaging traces the sub-micron and micron-sized dust grains in protoplanetary disks and it offers constraints on the geometrical properties of the disks and potential companions, particularly if those companions have circumplanetary disks. Aims. We use the VISIR instrument and its upgrade NEAR on the VLT to take new mid-IR images of five (pre-)transition disks and one circumstellar disk with proposed planets and obtain the deepest resolved mid-IR observations to date in order to put new constraints on the sizes of the emitting regions of the disks and the presence of possible companions. Methods. We derotated and stacked the data to find the disk properties. Where available, we compare the data to PRODIMO (Protoplanetary Disk Model) radiation thermo-chemical models to achieve a deeper understanding of the underlying physical processes within the disks. We applied the circularised point spread function subtraction method to find upper limits on the fluxes of possible companions and model companions with circumplanetary disks. Results. We resolved three of the six disks and calculated position angles, inclinations, and (upper limits to) sizes of emission regions in the disks, improving upper limits on two of the unresolved disks. In all cases the majority of the mid-IR emission comes from small inner disks or the hot inner rims of outer disks. We refined the existing PRODIMO HD 100546 model spectral energy distribution (SED) fit in the mid-IR by increasing the PAH abundance relative to the ISM, adopting coronene as the representative PAH, and increasing the outer cavity radius to 22.3 AU. We produced flux estimates for putative planetary-mass companions and circumplanetary disks, ruling out the presence of planetary-mass companions with L &gt; 0.0028 L⊙ for a &gt; 180 AU in the HD 100546 system. Upper limits of 0.5–30 mJy are obtained at 8–12 μm for potential companions in the different disks. We rule out companions with L &gt; 10−2 L⊙ for a &gt; 60 AU in TW Hydra, a &gt; 110 AU in HD 169142, a &gt; 150 AU in HD 163296, and a &gt; 160 AU in HD 36112. Conclusions. The mid-IR emission comes from the central regions and traces the inner areas of the disks, including inner disks and inner rims of outer disks. Planets with mid-IR luminosities corresponding to a runaway accretion phase can be excluded from the HD 100546, HD 169142, TW Hydra, and HD 36112 systems at separations &gt;1′′. We calculated an upper limit to the occurrence rate of wide-orbit massive planets with circumplanetary disks of 6.2% (68% confidence). Future observations with METIS on the ELT will be able to achieve a factor of 10 better sensitivity with a factor of five better spatial resolution. MIRI on JWST will be able to achieve 250 times better sensitivity. Both will possibly detect the known companions to all six targets.

Improving the RSM map exoplanet detection algorithm

Context. High-contrast imaging is one of the most challenging techniques for exoplanet detection. It relies on sophisticated data processing to reach high contrasts at small angular separations. Most data processing techniques of this type are based on the angular differential imaging observing strategy to perform the subtraction of a reference point spread function (PSF). In addition, such techniques generally make use of signal-to-noise (S/N) maps to infer the existence of planetary signals via thresholding. Aims. An alternative method for generating the final detection map was recently proposed with the regime-switching model (RSM) map, which uses a regime-switching framework to generate a probability map based on cubes of residuals generated by different PSF subtraction techniques. In this paper, we present several improvements to the original RSM map, focusing on novel PSF subtraction techniques and their optimal combinations, as well as a new procedure for estimating the probabilities involved. Methods. We started by implementing two forward-model versions of the RSM map algorithm based on the LOCI and KLIP PSF subtraction techniques. We then addressed the question of optimally selecting the PSF subtraction techniques to optimise the overall performance of the RSM map. A new forward-backward approach was also implemented to take into account both past and future observations to compute the RSM map probabilities, leading to improved precision in terms of astrometry and lowering the background speckle noise. Results. We tested the ability of these various improvements to increase the performance of the RSM map based on data sets obtained with three different instruments: VLT/NACO, VLT/SPHERE, and LBT/LMIRCam via a computation of receiver operating characteristic curves. These results demonstrate the benefits of these proposed improvements. Finally, we present a new framework to generate contrast curves based on probability maps. The contrast curves highlight the higher performance of the RSM map compared to a standard S/N map at small angular separations.