---
title: PRIMAger Hyperspectral Imaging (PHI)
url: https://www.emergentmind.com/topics/primager-hyperspectral-imaging-phi
type: topic
---

# PRIMAger Hyperspectral Imaging (PHI)

Searching arXiv for the cited PRIMAger/PHI papers to ground the article in current literature.
PRIMAger Hyperspectral Imaging (PHI) is the hyperspectral imaging capability of PRIMAger, the infrared camera of the PRobe far-Infrared Mission for Astrophysics (PRIMA). In the Phase A instrument and General Observer descriptions, PRIMAger will offer hyperspectral imaging in medium resolution bands (\(R \sim 8\), using a linear variable filter) from 25 to 80 \(\mu\)m, with two contiguous LVF-equipped modules, PHI1 and PHI2, while a distinct mode provides broad band photometric and polarimetric imaging at longer wavelengths [2509.01727] [2509.17646]. Earlier confusion-noise studies describe PRIMAger as providing unique hyperspectral imaging simultaneously covering 25–235 \(\mu\)m and use that capability to test deblending far below the classical confusion limit [2404.06935]. This suggests evolving instrument specifications between mission-concept simulations and later Phase A descriptions.

## 1. Mission context and scope

PRIMA is an infrared observatory for the next decade, currently in Phase A, with a 1.8m telescope actively cooled to 4.5K. On board, an infrared camera, PRIMAger, equipped with ultra-sensitive kinetic inductance detector (KID) arrays, will provide observers with coverage of mid-infrared to far-infrared wavelengths from 24 to 264 microns. PRIMAger is being developed by an international collaboration bringing together French institutes (Laboratoire d'Astrophysique de Marseille and CEA) through the center National d'Etudes Spatiales (CNES, France), the Netherlands Institute for Space Research (SRON, Netherlands), and the Cardiff University (UK) in Europe, as well as the Jet Propulsion Laboratory (JPL) and Goddard Space Flight Center (GSFC) in the USA [2509.01727].

A recurring point of confusion in the literature is the spectral extent assigned to PHI or to PRIMAger’s hyperspectral capability. The distinction is explicit in the published documents.

| Source | Hyperspectral description |
|---|---|
| [2404.06935] | unique hyperspectral imaging simultaneously covering 25-235 \(\mu\)m |
| [2509.01727] | Hyperspectral mode will cover the 24-84 microns wavelength range with a spectral resolution \(R=8\) |
| [2509.17646] | hyperspectral imaging in medium resolution bands (\(R \sim 8\), using a linear variable filter) from 25 to 80 \(\mu\)m |

The later Phase A descriptions also partition PRIMAger into two imaging modes: the Hyperspectral mode will cover the 24-84 microns wavelength range with a spectral resolution \(R=8\), while the Polarimetric mode will provide polarimetric imaging in 4 broad bands, from 80 to 264 microns [2509.01727]. A plausible implication is that the 2024–2025 confusion analyses model a broader concept-stage hyperspectral configuration, whereas the instrument paper and GO program paper describe the Phase A implementation.

## 2. Optical architecture and detector system

PHI shares PRIMA’s 1.8 m, 4.5 K telescope focus via a 50 mm collimated input beam. A two-axis, cryogenic beam-steering mirror (BSM) at the telescope exit pupil, with \(\pm 10'\) throw, \(\lesssim 0.33''\) RMS accuracy, and \(\leq 3\) mW dissipation at 4.5 K, directs the beam into PRIMAger without moving the spacecraft. The BSM allows rapid (\(\leq 30'/\)s) 2D scanning patterns, including raster, Lissajous, and boustrophedon, essential both to build spectra and to mitigate \(1/f\) detector noise [2509.01727].

The cold optical train uses a single, bare-Al common mirror to refocus the BSM beam to an intermediate focal plane inside the 1 K instrument box. Two field mirrors define cold stops for PHI and PPI, each followed by dichroic-free optics. The PHI channel is relayed by two further aspherized fold mirrors to achieve an effective focal ratio of \(f/21\) at the detector. All optics are aluminum 6061-T6, machined in a homothetic design so that cooldown-induced contraction preserves alignment to \(\lesssim 5\ \mu\)m [2509.01727].

Instead of wheels or dichroics, PHI uses wafer-scale Linear Variable Filters (LVFs) bonded directly above each KID array. Each LVF’s resonant metal-mesh pattern produces a smoothly varying band-centre along one detector axis, and no moving filter mechanism is required. PHI employs lens-absorber-coupled, hybrid Microwave Kinetic Inductance Detectors, with arrays fabricated on SiN membranes and mated to monolithic Si lenslets via precision glue pillars that maintain a \(4\ \mu\)m micro-lens gap [2509.01727].

| Parameter | PHI1 | PHI2 |
|---|---|---|
| Wavelength range | 24–45 \(\mu\)m | 45–84 \(\mu\)m |
| Array size | \(61 \times 24\) pixels | \(34 \times 14\) pixels |
| Field of view | \(\approx 3.9' \times 1.3'\) | \(\approx 3.9' \times 1.4'\) |
| Pixel scale | \(3.8''\) | \(6.8''\) |

The instantaneous spatial sampling is set by \(f\lambda\) optics and the pixel pitch. PHI1 has \(3.8''\) pixels, approximately \(0.8\ f\lambda\) at 24 \(\mu\)m, and PHI2 has \(6.8''\) pixels, approximately \(0.8\ f\lambda\) at 45 \(\mu\)m. Because this is coarser than Nyquist (\(0.5\ f\lambda\)), the BSM or spacecraft scan fills in the PSF [2509.01727].

## 3. Spectral formation and observing strategy

PHI1 covers \(24\ \mu\mathrm{m} \leq \lambda \leq 45\ \mu\mathrm{m}\), and PHI2 covers \(45\ \mu\mathrm{m} \leq \lambda \leq 84\ \mu\mathrm{m}\). Each pixel’s band-centre \(\lambda_c\) varies along the long axis of the LVF, and the FWHM of the passband is \(\Delta \lambda \approx \lambda_c / R\), with \(R \equiv \lambda_c / \Delta \lambda = 8\). In practice, \(R \geq 8\) everywhere in PHI1 and PHI2 [2509.01727].

PHI has no discrete filters. Spectra are formed by scanning a point source across rows of pixels in the spectral direction. As the BSM moves the source through the LVF gradient, successive pixels measure contiguous \(\Delta \lambda\) slices. A full PHI spectrum from 24 to 84 \(\mu\)m is built by scanning first across PHI1’s long axis and then PHI2’s, or vice versa, in a single observation, with typical scan length matching the \(3.9'\) array dimension on the sky [2509.01727].

The observing modes are tied to this spectral-formation strategy. Large maps use spacecraft slews in boustrophedon patterns, with scan legs up to a few degrees and perpendicular spacing approximately equal to the pixel pitch, combined with slow BSM modulation to fill in \(f\lambda\) sampling and mitigate \(1/f\). Small maps use purely BSM-driven 2D scans over \(\lesssim 30'\) fields at speeds up to \(30'/\)s. PHI always requires scan modulation to reconstruct spectra; there is no staring mode [2509.01727].

In the \(\pi\)-IR survey concept, each scan imparts a spectral-slit mapping: at each sky coordinate, subsequent visits with different LVF columns sample different \(\lambda_0\). After 4–6 orthogonal scan passes the full 25–80 \(\mu\)m range is sampled at \(R=8\), and the software assembles 3D cubes \((\mathrm{RA}, \mathrm{Dec}, \lambda)\) on regular spatial and spectral grids [2509.17646].

## 4. Calibration and data products

The detector modules are housed at 125 mK and are thermally linked to a continuous Adiabatic Demagnetization Refrigerator. Eight independent frequency-multiplexed RF chains read out 2.6–4.9 GHz, sampling at 5 Gsps. A Xilinx Kintex KU060 FPGA per chain synthesizes a comb of tones via a polyphase filter-bank IFFT, sends them to the detectors via low-loss coax, and digitizes the returned tones through a digital spectrometer (WOLA PFB). The on-board tone-tracking firmware delivers per-pixel \(I/Q\) timestreams at approximately 100 Hz [2509.01727].

Calibration is based on both internal and celestial references. A set of on-board, cold (\(\sim 4\) K) blackbody sources can be inserted via a small flip-in mirror for detector gain and linearity checks; these calibrators are used at the start and end of each observing block. Regular observations of well-characterized FIR standard stars, including asteroids and planetary nebulae, provide absolute flux calibration to \(\lesssim 5\%\). Flat-fielding is accomplished by nodding a bright, spatially uniform field, exemplified by the Moon’s limb, across the LVF gradient to map pixel gains versus wavelength. Detector responsivities are stable to \(<1\%\) over 24 hr [2509.01727].

The raw data format consists of time-ordered \(I/Q\) streams for each KID at approximately 100 Hz, tagged with BSM angle, telescope attitude, and temperature monitors, and stored in PRIMA Level 0 packets. The Level 1 to Level 2 pipeline comprises deglitching and cosmic-ray removal; channelization, assigning each timestamp to a sky wavelength based on LVF pixel \(\lambda_c\) and BSM position; spectral extraction; flat-field and dark subtraction using internal calibrator data; and cube reconstruction, gridding \(I(\alpha,\delta,\lambda)\) into a 3D FITS cube with user-selectable \(\Delta \lambda\) bins, default \(R=8\), and spatial pixels matching the 3.8–6.8 arcsec physical pitch [2509.01727].

The FITS data cube is specified as a multi-extension file with a PRIMARY header containing observation metadata, an extension “WAVELENGTH”, an extension “MAP” containing calibrated surface brightness in MJy/sr, an extension “ERROR” containing \(1\sigma\) uncertainties, and an extension “DQ” containing a data-quality mask per voxel [2509.01727].

## 5. Confusion noise and Bayesian deblending

The confusion-limited use case of PRIMAger has been developed most explicitly in the simulation papers by Donnellan and collaborators. In the 2024 study, images representing a deep, \(1500\ \mathrm{hr}\ \mathrm{deg}^{-2}\) PRIMAger survey were synthesised with realistic instrumental and confusion noise [2404.06935]. In the 2025 follow-up, a 1 \(\mathrm{deg}^2\) deep-field survey was simulated by populating the sky with sources drawn from the SIDES extragalactic model of Bethermin et al. 2017, convolving each source list at each band with the appropriate PSF, adding Gaussian instrumental noise, and adding confusion noise as spatially correlated background fluctuations using HELP cirrus+clustering templates. Source positions were taken from a Euclid-like near-IR catalogue with \(\rho \approx 150\ \mathrm{arcmin}^{-2}\) and negligible positional error \(<0.2''\) [2512.13682].

The forward model at wavelength \(\lambda\) is
\[
d_\lambda = P_\lambda F_\lambda + n_\lambda,
\]
where \(d_\lambda\) is the vector of observed pixel values, \(F_\lambda\) is the vector of source fluxes, \(P_\lambda\) is the design matrix whose columns are the normalized PSF at each source position, and \(n_\lambda \sim N(0,N_\lambda)\) is Gaussian noise with covariance \(N_\lambda\). Donnellan et al. present a new Bayesian modelling approach, XID+stepwise, that exploits PRIMAger's hyperspectral imaging to derive self-consistent, informative flux priors by sequentially propagating constraints from short to long wavelengths [2512.13682].

Fluxes are inferred in ascending-\(\lambda\) order. At the first band no flux prior is used; once the posterior mean and covariance at band \(\lambda_1\) are obtained, they are pinned and used as priors for all subsequent \(\lambda > \lambda_1\). In practice a narrow Gaussian prior is used instead of a strict delta function to preserve error propagation. At each band the joint posterior over \(F_\lambda\) is sampled, and the resulting marginal flux estimates and uncertainties inform the next band [2512.13682].

This stepwise construction is not the only route considered. The 2024 study shows that catalogues of galaxies with a high purity (\(>95\) per cent) can be constructed at a source density of \(42\mathrm{k}\ \mathrm{deg}^{-2}\) using PRIMAger data alone, and the 2025 study further shows that positional priors from blind source detection followed by deblending via XID+ enables PRIMAger to achieve sensitivity beyond the confusion limits using PRIMAger data alone [2404.06935] [2512.13682].

## 6. Sensitivity and flux-recovery performance

In the instrument performance budget, PHI achieves a 5\(\sigma\) point-source sensitivity of 1.2–4.1 mJy per 10 per cent \(\Delta \lambda\) channel in a raster mapping mode over 1 \(\mathrm{deg}^2\) in 10 h. The noise model includes dark NEP \(\lesssim 3 \times 10^{-20}\ \mathrm{W\ Hz^{-1/2}}\), photon NEP of order \(1\)–\(2 \times 10^{-19}\ \mathrm{W\ Hz^{-1/2}}\) per pixel, single-mode throughput per pixel \(A\Omega=\lambda^2\), coupling efficiency \(\langle \eta \rangle \approx 0.7\), and end-to-end transmission excluding LVF \(\geq 45\%\). Detector \(1/f\) knee is at \(f_\mathrm{knee} \approx 0.1\) Hz and the KID time constant is \(\tau_\mathrm{det} \approx 20\) ms [2509.01727].

For the \(\pi\)-sr survey parameters, with scan speed \(250''\,\mathrm{s}^{-1}\), total on-sky time \(\simeq 2059\) h, and SNR \(=5\), the 5\(\sigma\) point-source sensitivities are 2.54 mJy at 34.3 \(\mu\)m and 3.44 mJy at 64.5 \(\mu\)m; PHI2 channels at 45–84 \(\mu\)m would lie in the 2–5 mJy range accordingly [2509.17646].

The confusion-noise studies quantify how far Bayesian deblending can push below the classical confusion floor. The 2024 deep-survey simulation shows fluxes measured with an accuracy better than 20 per cent to flux levels of 0.16, 0.80, 9.7, and 15 mJy at 47.4, 79.7, 172, and 235 \(\mu\)m, respectively. These levels are a factor of \(\sim 2\) and \(\sim 3\) fainter than the classical confusion limits for 72–96 \(\mu\)m and 126–235 \(\mu\)m, respectively. At \(1.5 \leq z \leq 2\), 8–10 of the 10 channels covering 47–235 \(\mu\)m are detected and accurately measured for sources with \(2 \leq \log(\mathrm{SFR}) \leq 2.5\), a 0.5 dex improvement on what might be expected from the classical confusion limit [2404.06935].

The 2025 XID+stepwise analysis strengthens this result. With Euclid-like prior source positions, fluxes are recovered to within 20 per cent to 0.2–0.7 mJy across 45–84 \(\mu\)m, corresponding to factors of 1.3–3.4 fainter than the confusion limit. In the most confusion-dominated channels, accurate fluxes are measured to 0.9, 2.5, 7.6, and 14.8 mJy at 92, 126, 183, and 235 \(\mu\)m, respectively, which are factors of 3–5 better than the confusion limit. Using a deeper Euclid-based prior catalogue and weak ancillary flux priors at 25 \(\mu\)m yields further improvements, reaching up to a factor \(\sim 7\) fainter than the confusion limit at 96 \(\mu\)m [2512.13682].

## 7. Survey programs and scientific applications

The scientific rationale for PHI is broad. PRIMAger’s observational capabilities have been tailored to answer fundamental astrophysical questions such as black hole and star-formation co-evolution in galaxies, the evolution of small dust grains over a wide range of redshifts, and the effects of interstellar magnetic fields in various environments, while also opening discovery space with versatile photometric and polarimetric capabilities [2509.01727].

Within the General Observer framework, the \(\pi\)-IR survey covers roughly \(\pi\) sr, or \(\simeq 10{,}300\ \mathrm{deg}^2\), and would exploit PRIMAger’s hyperspectral and polarimetric modes to collect data on about \(8 \times 10^6\) galaxies to \(z \sim 4\). The \(R=8\) spectral resolution of the PRIMAger Hyperspectral Imaging filters will enable users to study the emission of polycyclic aromatic hydrocarbon (PAH), and a large sample of galaxies will be observed with the polarimetric bands of PRIMAger [2509.17646].

The instrument paper states that PHI’s \(R=8\) sampling recovers broad PAH features, for example the 34 \(\mu\)m band at \(z \approx 1.2\), and continuum breaks, delivering redshifts to \(\Delta z \approx 0.1\) for dusty star-forming galaxies in confusion-limited deep fields. A 1 \(\mathrm{deg}^2\) deep survey with \(t=1500\) h achieves \(5\sigma \approx 1.2 \times 10^{-17}\ \mathrm{W\ m^{-2}}\) line sensitivity, integrated over \(\Delta \lambda\), sufficient to map the [O IV] 25.9 \(\mu\)m line in \(L^\star\) galaxies at \(z \approx 2\) [2509.01727].

Additional use cases extend from Galactic structure to compact dusty systems. PHI can spectrally map warm dust and PAH emission in giant molecular clouds such as Orion A, resolving substructures at 0.1 pc scales at 400 pc with \(\Delta \lambda \approx 3\ \mu\)m bins. Simulations show the ability to disentangle multiple temperature components and resolve 34 \(\mu\)m PAH emission from continuum in a \(4' \times 1'\) strip in \(\lesssim 1\) hr. Other key targets listed for PHI include protoplanetary disks, via solid-state features such as 33.6 \(\mu\)m forsterite, and AGN torus studies, tracing the warm dust continuum break near 30 \(\mu\)m at \(z \lesssim 0.4\) in moderate-luminosity Seyferts [2509.01727].

Across these science cases, one theme is consistent: hyperspectral imaging is used not only to provide \(R \approx 8\) mid-infrared spectral sampling, but also to mitigate source confusion by adding the spectral dimension. The confusion-noise studies indicate that, when employing XID+, confusion noise will not limit the key science from PRIMA extragalactic imaging surveys [2512.13682].

Source: https://www.emergentmind.com/topics/primager-hyperspectral-imaging-phi