Optimising reference library selection for reference-star differential imaging of discs with SPHERE/IRDIS (2509.03325v1)

Abstract: The direct detection of circumstellar discs through high-contrast imaging provides key insights into the history and dynamics of planetary systems. Pole-on discs, especially faint debris discs, are difficult to detect and require careful consideration during post-processing to remove stellar residuals from the data while preserving the disc signal. Reference-star differential imaging (RDI) serves as one of the primary post-processing methods for disc observations. We aim to develop a method of reference frame selection that is optimised for the reduction of pole-on discs. Method: We performed principal component analysis based RDI on seven known disc targets and 20 disc-free targets with varying observational conditions, using reference libraries built from frames preselected to best match different observational, atmospheric, and stellar parameters of the science frames. The contrast of the disc-free reductions was measured, and forward modelling was used to estimate the signal loss from over-subtraction using synthetic pole-on discs with two different widths and four different radii. The signal-to-noise ratio (S/N) of the real disc targets was measured. Results: Diverse reference libraries built using subsets of frames that closely matched different parameters achieved the best disc S/N and smallest deviation from the best contrast of each target, outperforming libraries built using a single criterion as a selection metric. Libraries built using frame-to-frame Pearson correlation coefficient alone as a selection criterion achieved the best mean contrast overall. Both selection metrics performed consistently well for all disc radii and observational conditions. We also found that reference libraries built using frames observed close in time to the science frame performed well for discs at small separations, giving the best contrast for ~30% of the targets at a radius of 20px.