Application of a damped Locally Optimized Combination of Images method to the spectral characterization of faint companions using an Integral Field Spectrograph

High-contrast imaging instruments are now being equipped with integral field spectrographs (IFS) to facilitate the detection and characterization of faint substellar companions. Algorithms currently envisioned to handle IFS data, such as the Locally Optimized Combination of Images (LOCI) algorithm, rely upon aggressive point-spread-function (PSF) subtraction, which is ideal for initially identifying companions but results in significantly biased photometry and spectroscopy due to unwanted mixing with residual starlight. This spectro-photometric issue is further complicated by the fact that algorithmic color response is a function of the companion’s spectrum, making it difficult to calibrate the effects of the reduction without using iterations involving a series of injected synthetic companions. In this paper, we introduce a new PSF calibration method, which we call “damped LOCI”, that seeks to alleviate these concerns. By modifying the cost function that determines the weighting coefficients used to construct PSF reference images, and also forcing those coefficients to be positive, it is possible to extract companion spectra with a precision that is set by calibration of the instrument response and transmission of the atmosphere, and not by post-processing. We demonstrate the utility of this approach using on-sky data obtained with the Project 1640 IFS at Palomar. Damped-LOCI does not require any iterations on the underlying spectral type of the companion, nor does it rely upon priors involving the chromatic and statistical properties of speckles. It is a general technique that can readily be applied to other current and planned instruments that employ IFS’s.

💡 Research Summary

The paper addresses a fundamental problem in high‑contrast imaging with integral‑field spectrographs (IFS): while aggressive point‑spread‑function (PSF) subtraction algorithms such as LOCI (Locally Optimized Combination of Images) are excellent for detecting faint companions, they severely bias the photometry and spectroscopy of those companions. The bias arises because the linear combination of reference images, optimized to minimize residual speckle noise, also partially fits and removes the companion’s signal. This leads to a “flux depletion” that can be gray (wavelength‑independent) or, more problematically for IFS data, chromatic, causing spectral cross‑talk between adjacent wavelength channels. Existing mitigation strategies rely on iterative injections of synthetic companions and on priors about the companion’s spectrum, which are computationally expensive and can still leave significant systematic errors.

The authors propose “damped LOCI,” a modification of the LOCI cost function that simultaneously minimizes residual speckle noise and maximizes the retained companion flux. Two key changes are introduced: (1) a non‑negativity constraint on the LOCI coefficients (cₖ ≥ 0), preventing the algorithm from subtracting the companion’s flux; and (2) an additional penalty term that explicitly rewards flux preservation in a user‑defined subtraction zone (S‑zone) where a companion is expected. The optimization problem becomes

 min ‖W·(T − R·c)‖² + λ‖S·c‖² subject to c ≥ 0,

where T is the target image, R the set of reference images, W a weighting map for the optimization zone (O‑zone), S a mask for the subtraction zone, and λ a tunable regularization parameter. In matrix form this yields a regularized least‑squares problem that can be solved with non‑negative least squares (NNLS) algorithms. The penalty term effectively conditions the problem, reducing the ill‑posedness that arises in IFS data where many wavelength slices are highly correlated.

The method is tested on on‑sky data from the Project 1640 IFS at Palomar, which provides 23 wavelength channels across the J and H bands (R ≈ 40) behind a coronagraphic adaptive‑optics system. The authors compare three pipelines: (a) standard LOCI, (b) LOCI with binary masking of the companion region, and (c) damped LOCI. The results show that:

• For companions whose flux is comparable to the residual speckle level, damped LOCI recovers the spectral energy distribution (SED) to within the telluric calibration uncertainty (≈5 %). Standard LOCI underestimates the flux by 30–60 % across the band.
• For high‑contrast companions (contrast ≤ 10⁻³), damped LOCI eliminates the systematic gray loss seen in standard LOCI and dramatically reduces spectral cross‑talk, delivering accurate relative colors and absolute fluxes.
• The algorithm improves the overall contrast floor by 1–2 mag relative to standard LOCI, achieving detection limits near 10⁻⁵ at 1ʺ in 20 min exposures.
• The conditioning of the reference‑matrix M = RᵀWR is markedly better with the non‑negative constraint and λ‑penalty, making the solution robust even when the number of independent reference frames is smaller than the number of PSF cores in the optimization zone.

Importantly, damped LOCI does not require prior knowledge of the companion’s spectrum, nor does it need iterative synthetic‑planet injections to calibrate the algorithmic throughput. This simplifies the data‑reduction pipeline and makes real‑time processing feasible for future IFS‑based high‑contrast instruments such as GPI, SPHERE, and SCExAO. The authors argue that the technique is generic: any IFS that relies on chromatic diversity for speckle suppression can adopt damped LOCI to obtain unbiased spectra while retaining the aggressive speckle‑noise suppression needed for detection.

In conclusion, the paper presents a mathematically elegant yet practically implementable solution to the long‑standing trade‑off between detection sensitivity and spectro‑photometric fidelity in IFS high‑contrast imaging. By augmenting the LOCI cost function with positivity constraints and a flux‑preservation penalty, the authors achieve simultaneous speckle suppression and accurate spectral extraction, validated on real data. This work is likely to become a new standard for processing IFS data in exoplanet direct‑imaging surveys, enabling more reliable atmospheric characterizations of faint companions without the need for cumbersome post‑processing calibrations.