# A temperate super-Jupiter imaged with JWST in the mid-infrared

E. C. Matthews, A. L. Carter, P. Pathak, C. V. Morley, M. W. Phillips, S.

Krishanth P. M., F. Feng, M. J. Bonse, L. A. Boogaard, J. A. Burt, et al.

► To cite this version:

E. C. Matthews, A. L. Carter, P. Pathak, C. V. Morley, M. W. Phillips, et al. A temperate super-Jupiter imaged with JWST in the mid-infrared. Nature, 2024, 633, pp.789-792. 10.1038/s41586-024-07837-8. insu-04853413

HAL Id: insu-04853413

https://insu.hal.science/insu-04853413v1

Submitted on 24 Dec 2024

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A temperate super-Jupiter imaged with JWST in the mid-infrared

https://doi.org/10.1038/s41586-024-07837-8

Received: 5 April 2024

Accepted: 16 July 2024

Published online: 24 July 2024

Open access

<figure><img src="image_2.png" /></figure>

Check for updates

E. C. Matthews<fnref n="1" />,²,³, A. L. Carter², P. Pathak³, C. V. Morley<fnref n="4" />, M. W. Phillips⁵, S. Krishanth P. M.⁶, F. Feng⁷, M. J. Bonse⁸,<fnref n="9" />, L. A. Boogaard<fnref n="1" />, J. A. Burt¹⁰, I. J. M. Crossfield¹¹, E. S. Douglas⁶, Th. Henning<fnref n="1" />, J. Hom⁶, C.-L. Ko⁶, M. Kasper¹², A.-M. Lagrange¹³, D. Petit dit de la Roche¹⁴ & F. Philipot¹³

Of the approximately 25 directly imaged planets to date, all are younger than 500 Myr, and all but six are younger than 100 Myr (ref. 1). Eps Ind A (HD209100, HIP108870) is a K5V star of roughly solar age (recently derived as 3.7–5.7 Gyr (ref. 2) and 3.5$^{+0.8}_{-1.3}$ Gyr (ref. 3)). A long-term radial-velocity trend<fnref n="4" />,⁵ and an astrometric acceleration⁶,⁷ led to claims of a giant planet²,⁸,<fnref n="9" /> orbiting the nearby star (3.6384 ± 0.0013 pc; ref. 10). Here we report JWST coronagraphic images which reveal a giant exoplanet that is consistent with these radial and astrometric measurements but inconsistent with the previously claimed planet properties. The new planet has a temperature of approximately 275 K and is remarkably bright at 10.65 and 15.50 µm. Non-detections between 3.5 and 5.0 µm indicate an unknown opacity source in the atmosphere, possibly suggesting a high-metallicity, high carbon-to-oxygen ratio planet. The best-fitting temperature of the planet is consistent with theoretical thermal evolution models, which were previously untested at this temperature range. The data indicate that this is probably the only giant planet in the system, and therefore we refer to it as b, despite it having significantly different orbital properties than the previously claimed planet b.

We observed Eps Ind A with the mid-infrared instrument (MIRI)¹¹ coronagraph onboard the James Webb Space Telescope (JWST) on 3 July 2023 using two narrowband filters (10.65 and 15.50 µm). Figure 1 shows a bright point source detected in the north-east quadrant at a separation of 4.11". This is the opposite quadrant than expected based on previous orbital solutions²,⁸,<fnref n="9" />. The source is unresolved, has apparent magnitudes 13.16 and 11.20 mag at 10.65 and 15.50 µm, and is consistent with a cold, Jupiter-sized object at the host star distance. We ruled out chance-aligned background objects and recovered the companion at a low signal-to-noise ratio in archival observations, confirming that it is physically associated with the host star. The source has significantly different orbital properties than the literature solutions for Eps Ind Ab, but we conclude that it is probably the only massive planet in the system and is responsible for the previously observed radial velocity and astrometric acceleration of the host star.

Owing to the unexpected location of the source, we first assessed whether it could be a chance-aligned background object. This is statistically unfavoured. Source-counting studies with the MIRI broadband filters¹² indicate an expected background density of approximately 45 sources per square degree at least as bright as the point source, corresponding to just a 0.027% likelihood of finding such a source within 5" of the target. Although occurrence rates for giant, long-period planets are low¹³,¹⁴, the probability of a chance-aligned background object is much lower. A background contaminant was further constrained

based on archival data. The extremely high proper motion of Eps Ind A (4,708.15 ± 0.13 mas yr⁻<fnref n="1" />; ref. 10) means that the background location is unobstructed by the star in sufficiently old data. We studied previous observations of Eps Ind A, and no stationary background contaminant consistent with the point source was detected. A particularly strong constraint comes from archival observations by the infrared array камера (IRAC) onboard Spitzer¹⁵, which rules out companions brighter than 16.0 mag at 8 µm. The point source is 13.16 mag at 10.65 µm. Further, the source is unresolved, and most galaxies of this brightness would be spatially resolved at MIRI wavelengths. It is unlikely that a stationary background object would have evaded detection in all archival observations while also reproducing the candidate planet properties. A Solar System source (for example, an asteroid) can be excluded as it should have moved between the F1065C and F1550C observations. A transient source (for example, a burst) cannot be excluded based on archival data alone but is strongly disfavoured statistically.

We reanalysed archival data collected with the VISIR/NEAR instrument¹⁶,¹⁷ at the Very Large Telescope (VLT). These coronagraphic observations of Eps Ind A were collected over three nights during September 2019 using a broadband filter spanning approximately 10.0–12.5 µm. At approximately 15 au, a planet would have an orbital period of at least several decades, corresponding to little orbital motion in the approximately 3.8 yr between the VISIR and JWST observations, whereas a background source would have shifted by 18" over this time baseline.

<footnote n="1">¹Max Planck Institute for Astronomy, Heidelberg, Germany. ²Space Telescope Science Institute, Baltimore, MD, USA. ³Department of SPASE, Indian Institute of Technology Kanpur, Kanpur, India.</footnote>

<footnote n="4">⁴Department of Astronomy, University of Texas at Austin, Austin, TX, USA. ⁵Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh, UK. ⁶Steward Observatory and Department of Astronomy, Tucson, AZ, USA. ⁷Shanghai Jiao Tong University, Shanghai, People's Republic of China. ⁸ETH Zurich, Institute for Particle Physics & Astrophysics, Zurich, Switzerland.</footnote>

<footnote n="9">⁹Max Planck Institute for Intelligent Systems, Tübingen, Germany. ¹⁰Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA. ¹¹Department of Physics and Astronomy, University of Kansas, Lawrence, KS, USA. ¹²European Southern Observatory, Munich, Germany. ¹³LESIA, Observatoire de Paris, Université PSL, CNRS, Meudon, France. ¹⁴Department of Astronomy, University of Geneva, Versoix, Switzerland. ⁵e-mail: matthews@mpia.de</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** | Point source detected in JWST/MIRI coronagraphic images of Eps Ind A. The target was observed at 10.65 and 15.50 µm. Starlight was removed with RDI. Only the central portion of the field of view is shown. The stellar position is marked with an orange star. A bright point source was detected in the upper left corner of these images at a projected separation 4.11" (15.0 au at the distance of Eps Ind A). Top scale bar, 2.75" = 10.0 au; bottom scale bar, 2" = 7.28 au.</figcaption></figure>

Figure 2 shows the VISIR/NEAR observations. A faint point source at a signal-to-noise ratio of approximately 3 was observed. It is consistent with the expected flux and location of the JWST point source. We concluded that this is an archival redetection of the same object, which confirms its common proper motion and conclusively demonstrates that it is a planet (hereafter Eps Ind Ab).

The measured mid-infrared photometry (10.65 and 15.50 µm) of Eps Ind Ab is consistent with atmosphere models for an approximately 300 K planet, as demonstrated in Fig. 3. However, these models typically predict significant emission between 3.5 and 5.0 µm. Archival observations by the NaCo instrument at the Very Large Telescope (VLT)¹⁷,¹⁸ failed to detect Eps Ind Ab. Atmospheric models predict a 'window' for flux at these wavelengths that is between absorption features from molecular species including CH₄ and H₂O (ref. 19). This is observed for cold brown dwarfs²⁰,²¹ and expected for cold, giant planets but is in tension with NaCo upper limits (highlighted in Fig. 3). The NaCo non-detection suggests that there is another opacity source in the 3.5–5.0 µm region, such as molecular absorption or significant cloud or haze high in the atmosphere. Although some model grids explore cloudy scenarios²²,²³, these models require patchy clouds and the 3.5–5.0 µm flux is dominated by the cloud-free regions of the atmosphere. This, in turn, leads to a pronounced emission peak, similar to that seen in the ATMO models (Fig. 3). Sonora Elf Owl models²⁴ with temperature 275 K, high metallicity ([M/H] = 1.0), high carbon-to-oxygen ratio (×2.5 solar) and strong disequilibrium chemistry²⁵ (high log(K<sub>zz</sub>), where K<sub>zz</sub> is the eddy diffusion coefficient), are compatible with all in-hand observational constraints. This model includes significant CH₄, CO₂ and CO absorption, which suppresses the 3.5–5.0 µm flux, and a modest ammonia absorption feature that reproduces the 10.65 µm flux. A high metallicity is somewhat

surprising for a super-Jupiter planet; the C/O ratio is consistent with predictions for a planet formed by core accretion beyond the CO₂ ice line²⁶. Further photometric and spectroscopic characterization of the planet is crucial to confirm these preliminary indications of the atmospheric C/O ratio and to confirm whether chemistry or haze particles are the true cause of the suppressed 3.5–5.0 µm flux and allow a comparison to formation models.

To constrain the dynamical mass of Eps Ind Ab, we fitted a single-planet orbit to all published radial-velocity data (493 points from four facilities over 29 yr) for the two imaging epochs and to the Hipparcos–Gaia astrometry. This gave a mass of 6.31±0.60 Jupi masses ($M_J$) and a semimajor axis of $28.4_{-7.2}^{+10.0}$ au. The fitted orbit is eccentric ($e = 0.40_{-0.18}^{+0.15}$) and observed at apastron, and this is able to explain all in-hand data. Although previous studies claimed the presence of a $3.0 \pm 0.1 M_J$ planet with semimajor axis $8.8_{-0.1}^{+0.2}$ au (refs. 2,8,9), based on the radial-velocity data and astrometric information, the imaged planet imparts a significant acceleration on the host star, and so it is clear that there is no inner exoplanet with the previously published parameters. Although the data do not exclude the possibility that there is a second planet in the system, Eps Ind Ab is the dominant accelerator of the star. Future works should explore why the previous radial-velocity analyses led to underestimated mass and semimajor axis constraints for this system, which may be due to overfitting or if there are, in fact, several planets in the system.

Evolutionary models predict a temperature of approximately 280 K and a luminosity log(L/L<sub>☉</sub>) of approximately –7.2 for a 3.5 Gyr, 6.3 $M_J$ planet²⁷. Figure 4 shows evolutionary models²⁷ for two planets: a 6.3 $M_J$ object (matching the dynamical mass) and an 8.6 $M_J$ planet (matching the observed photometry), considering both equilibrium and disequilibrium chemistry. Based only on the age and MIRI photometry for Eps Ind Ab, the planet is overluminous relative to the models. However, the extent of the flux suppression at 3.5–5.0 µm remains unclear, and this impacts the derived temperature of the planet. The best-fitting model in Fig. 3 incorporates this suppressed flux. It gives a temperature of 275 K, which is very close to the temperature predictions from evolutionary models. Further photometry between 3.5 and 5.0 µm is crucial for confirming the extent of flux suppression, better constraining the temperature of Eps Ind Ab and, thereby, allowing a more thorough comparison to evolutionary models.

Eps Ind Ab is the coldest exoplanet to have been directly imaged. It has a temperature close to that of the coldest field brown dwarf (WISE 0855)²¹. The system is also co-moving with a widely separated brown dwarf binary²⁸,²⁹, making it a particularly valuable laboratory for comparative studies of substellar objects with a shared age and formation location. The exoplanet detection highlights the power of using indirect evidence to target direct detection efforts. Even though the detected planet does not match the previously claimed properties for an exoplanet orbiting this star, the long-term radial-velocity information was crucial in motivating the value of imaging this target. The photometry of Eps Ind Ab suggests it has a high metallicity

<figure><img src="image_4.png" /><figcaption>**Fig. 2** | Point source detected at a consistent location in VISIR/NEAR images of Eps Ind A, a, b, Co-added images (main panels) and PCA results for a small patch around the expected target location (insets). Panel b shows the images in a after convolution with a top-hat function.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 3** | Eps Ind Ab is consistent with theoretical atmosphere models with **suppressed 3.5–5.0 µm flux.** Measured photometry (squares and circles with uncertainties of $1\sigma$) and $5\sigma$ upper limits (triangles), compared to out-of-the-box model spectra<sup>24,27</sup> (coloured lines; circles indicate the integrated flux for each observed filter). Normalized filter profiles are overlaid as grey lines. Planet radii are scaled to match the measured 15.50 µm photometry, and labels</figcaption></figure>

indicate the temperature $T$, surface gravity $g$, log eddy diffusion coefficient (log(zz)), log-metallicity ([M/H]), and carbon-to-oxygen ratio (C/O). The ATMO models include non-equilibrium (NEQ) chemistry. A high metallicity and carbon-to-oxygen ratio are required to suppress the 3.5–5.0 µm flux to below the observed upper limit.

and carbon-to-oxygen ratio, and future works should aim to measure the planetary flux between 3.5 and 5.0 µm to determine the extent of flux suppression. Future works should also aim to expand the sample of cold exoplanets and determine whether other cold planets show the same low 3.5–5.0 µm flux, indicating high metallicity and high carbon-to-oxygen ratio. Two candidate companions to white dwarf stars<sup>30</sup> and several upcoming JWST observations of nearby accelerating stars<sup>31,32</sup> may mean that a small sample of cold, solar-age exoplanets can be assembled. The impact of atmospheric absorption is also critical when deriving mass detection limits, which are used to calculate occurrence rates, and when designing future observations

of cold exoplanets. The bright flux and wide separation of Eps Ind Ab mean that the planet is ideally suited to spectroscopic characterization efforts, which may more accurately constrain its metallicity and carbon-to-oxygen ratio.

## Online content

Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41586-024-07837-8.

<figure><img src="image_6.png" /><figcaption>**Fig. 4** | Comparison of Eps Ind Ab to evolutionary models. Cooling curves for a 6.3 $M_J$ planet and an 8.6 $M_J$ planet for models with equilibrium and with non-equilibrium chemistry<sup>27</sup>. The mid-infrared companion photometry is consistent with a more massive planet, although these models do not incorporate the high metallicity or carbon-to-oxygen ratio of the planet nor the correspondingly low 3.5–5.0 µm flux.</figcaption></figure>

1. NASA Exoplanet Archive. IPAC https://exoplanetarchive.ipac.caltech.edu/ (accessed 22 June 2024).

2. Feng, F. et al. Detection of the nearest Jupiter analogue in radial velocity and astrometry data. *Mon. Not. R. Astron. Soc.* **490**, 5002–5016 (2019).

3. Chen, M. et al. Precise dynamical masses of Eps Indi Ba and Bb: evidence of slowed cooling at the L/T transition. *Astron. J.* **163**, 288 (2022).

4. Endl, M. et al. The planet search program at the ESO coudé echelle spectrometer. III. The complete long camera survey results. *Astron. Astrophys.* **392**, 671–690 (2002).

5. Zechmeister, M. et al. The planet search programme at the ESO CES and HARPS. IV. The search for Jupiter analogues around solar-like stars. *Astron. Astrophys.* **552**, A78 (2013).

6. Brandt, T. D. The Hipparcos–Gaia catalog of accelerations. *Astrophys. J. Suppl. Ser.* **239**, 31 (2018).

7. Kervella, P., Arenou, F., Mignard, F. & Thévenin, F. Stellar and substellar companions of nearby stars from Gaia DR2. Binarity from proper motion anomaly. *Astron. Astrophys.* **623**, A72 (2019).

8. Feng, F., Butler, R. P., Vogt, S. S., Holden, B. & Rui, Y. Revised orbits of the two nearest Jupiters. *Mon. Not. R. Astron. Soc.* **525**, 607–619 (2023).

9. Philipot, F., Lagrange, A. M., Rubini, P., Kiefer, F. & Chomez, A. Updated characterization of long-period single companion by combining radial velocity, relative astrometry, and absolute astrometry. *Astron. Astrophys.* **670**, A65 (2023).
10. Gaia Collaboration et al. Gaia Data Release 3. Summary of the content and survey properties. *Astron. Astrophys.* **674**, A1 (2023).

11. Rieke, G. H. et al. The mid-infrared instrument for the James Webb Space Telescope. I. Introduction. *Publ. Astron. Soc. Pac.* **127**, 584 (2015).

12. Wu, C. K. W. et al. Source counts at 7.7–21µm in CEERS field with JWST. *Mon. Not. R. Astron. Soc.* **523**, 5187–5197 (2023).

13. Fulton, B. J. et al. California Legacy Survey. II. Occurrence of giant planets beyond the ice line. *Astrophys. J. Suppl. Ser.* **255**, 14 (2021).

14. Vigan, A. et al. The SPHERE infrared survey for exoplanets (SHINE). III. The demographics of young giant exoplanets below 300au with SPHERE. *Astron. Astrophys.* **651**, A72 (2021).

15. Marengo, M. et al. A Spitzer IRAC search for substellar companions of the debris disk star Eps Eridani. *Astrophys. J.* **647**, 1437–1451 (2006).

16. Pathak, P. et al. High-contrast imaging at ten microns: a search for exoplanets around Eps Indi A, Eps Eri, Tau Ceti, Sirius A, and Sirius B. *Astron. Astrophys.* **652**, A121 (2021).

17. Viswanath, G. et al. Constraints on the nearby exoplanet Eps Indi Ab from deep near- and mid-infrared imaging limits. *Astron. Astrophys.* **651**, A89 (2021).

18. Janson, M. et al. Imaging search for the unseen companion to Eps Ind A – improving the detection limits with 4 µm observations. *Mon. Not. R. Astron. Soc.* **399**, 377–384 (2009).

19. Burrows, A. et al. A nongray theory of extrasolar giant planets and brown dwarfs. *Astrophys. J.* **491**, 856–875 (1997).

20. Cushing, M. C. et al. The discovery of Y dwarfs using data from the Wide-field Infrared Survey Explorer (WISE). *Astrophys. J.* **743**, 50 (2011).

21. Luhman, K. L. et al. JWST/NIRSpec observations of the coldest known brown dwarf. *Astron. J.* **167**, 5 (2024).

22. Morley, C. V. et al. Water clouds in Y dwarfs and exoplanets. *Astrophys. J.* **787**, 78 (2014).

23. Lacy, B. & Burrows, A. Self-consistent models of Y dwarf atmospheres with water clouds and disequilibrium chemistry. *Astrophys. J.* **950**, 8 (2023).

24. Mukherjee, S. et al. The Sonora substellar atmosphere models. IV. Elf Owl: atmospheric mixing and chemical disequilibrium with varying metallicity and C/O ratios. *Astrophys. J.* **963**, 73 (2024).

25. Fegley, J. & Lodders, K. Atmospheric chemistry of the brown dwarf Gliese 229B: thermochemical equilibrium predictions. *Astrophys. J. Lett.* **472**, L37 (1996).

26. Öberg, K. I., Murray-Clay, R. & Bergin, E. A. The effects of snowlines on C/O in planetary atmospheres. *Astrophys. J. Lett.* **743**, L16 (2011).

27. Phillips, M. W. et al. A new set of atmosphere and evolution models for cool T-Y brown dwarfs and giant exoplanets. *Astron. Astrophys.* **637**, A38 (2020).

28. Scholz, R. D., McCaughrean, M. J., Lodieu, N. & Kuhlbrodt, B. ë Indi B: a new benchmark T dwarf. *Astron. Astrophys.* **398**, L29–L33 (2003).

29. McCaughrean, M. J. et al. Eps Indi Ba,Bb: the nearest binary brown dwarf. *Astron. Astrophys.* **413**, 1029–1036 (2004).

30. Mullally, S. E. et al. JWST directly images giant planet candidates around two metal-polluted white dwarf stars. *Astrophys. J. Lett.* **962**, L32 (2024).

31. Bardalez Gagliuffi, D. C. et al. Solving a solar neighborhood crime scene by imaging 14 Her c. *JWST Proposal Cycle 2*, 3337 (2023).

32. Matthews, E. C. et al. Super-Jupiters in our backyard: MIRI coronagraphic imaging of a massive planet/brown dwarf orbiting an M dwarf at 12pc. *JWST Proposal Cycle 3*, 5229 (2024).

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024
# Methods

## JWST observations and data reduction

We observed Eps Ind A with the JWST/MIRI coronagraphic imager¹¹,³³,³⁴ on 3 July 2023 (JWST General Observer programme 2243; ref. 35) in a single sequence. Our observations consisted of science, reference and background images, as detailed in Extended Data Table 1. We collected images with two narrowband coronagraphic filters, F1065C and F1550C, which were integrated for 3,772 and 3,922 s, respectively. Each coronagraphic filter has a dedicated four-quadrant phase-mask coronagraph³⁶ to suppress the stellar point spread function (PSF). Dedicated background observations of an empty patch of sky were made at the start of the sequence to allow subtraction of the 'glow stick' stray-light feature, first identified in MIRI commissioning³⁴. For each filter and each of the science and reference targets, we collected background observations using the same integration as a single science/reference dither, and we repeated this process over two dither positions. We used the same empty field for all background observations.

We collected PSF reference images of the star DI Tuc (HD211055) for reference differential imaging (RDI)³⁷. We collected five reference images for each filter using the five-point small-grid dither technique³⁸,³⁹ to account for the imperfect pointing of the JWST (commissioning data suggest a pointing stability of approximately 5–10 mas; ref. 40). The integration times per dither position were 2,631 and 2,724 s for the F1065C and F1550C observations, respectively. These integrations were chosen to match the peak flux of the Eps Ind A observations at each dither position. In advance of the JWST observations, we also screened the reference star for faint companions to ensure it was suitable for differential imaging. We observed DI Tuc with the Spectro-Polarimetric High-Contrast Exoplanet Research instrument (SPHERE) at VLT on 13 April 2023 (programme 110.25BR, P.I.E.C.M.) and collected 140 images, each with detector integration time (DIT) = 0.837 s (total exposure time 117 s), in the broadband H filter (λ = 1.6255 µm and Δλ = 0.291 µm). Images were dark- and flat-fielded and bad pixels corrected, and we then aligned all images based on the peak flux position and co-added across the stack of observations. No companions were detected anywhere in the SPHERE field of view (approximately 12" × 12"). Two companions are in the Gaia catalogue and the JWST field of view but outside the SPHERE field of view (separations 6.5" and 8.2"). These are sufficiently widely separated that they did not impact the PSF subtraction, though the brighter of these is clearly visible in the reduced data (Extended Data Fig. 1).

We used the spaceKLIP pipeline⁴¹,⁴² for the JWST data analysis. This pipeline provides various high-contrast imaging-specific functionalities for JWST data and, in particular, makes use of the jwst pipeline⁴³ and the pyKLIP Python package⁴⁴ for several key reduction steps. Our reduction follows the steps outlined in detail in ref. 42 and is briefly summarized here. We started from the stage 0 files ("uncal.fits") generated by the jwst pipeline. We performed ramp-fitting, calibrated the images from detector units to physical units, flat-fielded and subtracted the background images. We then subtracted the stellar PSF using RDI. We explored various parameters. The reductions shown in this work were made using principal component analysis (PCA) with a single subsection and a single annulus. Eight PCA modes were removed.

This process provided excellent starlight suppression for the F1550C images. The contrast performance is provided in Extended Data Fig. 2. The stellar PSF is less well matched between the science and reference images for F1065C, probably due to the larger time baseline between these observations (the order of observations is F1065C reference, F1550C reference, F1550C science and F1065C science, as listed in Extended Data Table 1). This reduced the planet sensitivity at close separations, and some residual starlight is clearly visible in these images. The full field of view for the 10.65 µm images is shown in Extended Data Fig. 1, and the central portion of both images is included in Fig. 1. Three astrophysical point sources were detected at F1065C:

1. A bright source to the north-east of Eps Ind A was not previously known. This source was confirmed to be the exoplanet Eps Ind Ab in the current work.

2. A faint source to the west of Eps Ind A is a background star. This source corresponds to Gaia DR3 6412595290591430784 and was also detected in archival Spitzer/IRAC observations.

3. A negative source to the south-east of Eps Ind A is a widely separated companion to DI Tuc. It was subtracted in our RDI reduction. This source corresponds to Gaia DR3 6411654761473726464.

We determined the companion properties using the forward-model PSF fitting routine provided by pyKLIP and implemented in spaceKLIP. We generated a PSF template at the location of the companion using webbpsf_ext. We then fitted the location and intensity of this PSF to the data to measure the photometry and astrometry of the observations; these values are listed in Extended Data Table 2.

We also reduced the data using non-negative matrix factorization⁴⁵,⁴⁶ to verify our results. Non-negative matrix factorization involves decomposing an image matrix into two non-negative matrices, one has features and the other weights for reconstructing the original matrix. The features matrix extracted from the RDI PSF library was used to model the stellar PSF in the target image frames. PSF subtraction was then performed to image the exoplanet. We performed the non-negative matrix factorization reduction using all frames from the DI Tuc PSF library to construct components. For F1065C, this involved using 605 frames to construct 605 components. For F1550C, the number of frames and components was 65. We tried alternative approaches, including using the entire PSF library with fewer components and using a smaller PSF library with components selected using their closeness to the target frames (using the Euclidean norm), but using the entire library with the maximum number of components consistently produced the best images.

## Background constraints derived from archival data products

Our analysis of archival data provided tight constraints on a possible background contaminant. If the point source were an extragalactic contaminant, it would have no proper motion. It would be found 4.1" to the north-east of the JWST Eps Ind A position, namely at right ascension 22 h 03 min 33.17 s and declination -56° 48' 06.0". For a background object, these absolute right ascension and declination coordinates would not change as Eps Ind A moves across the sky at a rate of 4.7" yr⁻¹ and would be well resolved from the star in sufficiently old archival data (allowing for a sufficient movement of the foreground star, relative to the background position). Even if a background contaminant were a distant object within the Milky Way, it would have a proper motion of at most a few milliarcseconds per year and at most a few tens of milliarcseconds over the approximately 2 decades spanned by the data products presented here. We searched the Gaia¹⁰ and 2MASS⁴⁷ catalogues and reanalysed data from the IRAC and the multiband imaging photometer (MIPS) onboard Spitzer, but we did not identify any background sources consistent with the point source.

Spitzer observations provide excellent background constraints and bracket the wavelengths of our JWST observations. Several Spitzer observations of Eps Ind A were carried out during 2004. Eps Ind A had moved approximately 90" in the 19 yr between the Spitzer and JWST observations. Eps Ind A was observed with Spitzer/IRAC on 1 May 2004 (programme ID 90; P.I. M. Werner), with data collected in all four IRAC channels (nominally 3.6, 4.5, 5.8 and 8.0 µm). The star was used primarily as a PSF reference star for a programme studying the Eps Eri debris disk¹⁵,⁴⁸,⁴⁹. We reanalysed these data, starting from the science-ready mosaic files from the Spitzer Heritage Archive for each of the four wavelength channels. Even though the background position was well separated from the star at this epoch, there was still a significant signal in the PSF wings (Extended Data Fig. 3). We used the Eps Eri observations from the same programme as a PSF reference for
RDI and for subtracting the stellar PSF wings, thereby allowing the best
possible contrast at the background location. We used a grid-based
approach to optimize the spatial (x, y) location and flux scaling of the
Eps Eri images and, thereby, to provide the best match to the PSF wing
structure of the Eps Ind images. We then subtracted this best-matching
image from the Eps Ind data to produce the RDI-subtracted images.
The science-ready images and RDI-subtracted images are shown in
Extended Data Fig. 3. The position of Eps Ind A during the Spitzer and
JWST epochs is highlighted. We calculated the flux limits based on
the pixel variance in the RDI-subtracted image, with known sources
masked, and validated our limits by performing aperture photometry
on identified sources in the image.

Eps Ind A was also observed with Spitzer/MIPS on 13 October 2004 (programme ID 41; PI G. Rieke)⁵⁰. Data were collected in all three MIPS channels (nominally 24, 70 and 160 µm). We analysed the science-ready mosaic files from the Spitzer Heritage Archive. Eps Ind A is unsaturated in the MIPS images, and the JWST position is well resolved from the star (Extended Data Fig. 4) and within the field of view for the 24 and 70 µm channels. We calculated the flux upper limit based on the standard deviation of pixel values in a large box around the position of Eps Ind A, excluding pixels within 1 arcmin of the star. No sources were observed within 10" of the JWST Eps Ind A position in the Spitzer IRAC 8.0 µm images or in any of the MIPS images. One source was identified approximately 9" west of Eps Ind A in the 3.6, 4.5 and 5.8 µm images. This source was also identified in the JWST/MIRI images and is labelled 2 in Extended Data Fig. 1. The location and magnitude measurements correspond to that of the star Gaia DR3 6412595290591430784. This is a background object and is not associated with Eps Ind A.

A background spectrum for the point source is shown in Extended Data Fig. 5. This indicates the flux and upper limits that would apply to the object if it were a stationary background object or a background object that had moved by less than a few hundred milliarcseconds per year. That is, these upper limits do not apply to the planet, which has the same proper motion as Eps Ind A and was not present at this location in the background observations. The 8.0 µm images provide a particularly stringent constraint on the background scenario. We ruled out background objects to 16.07 mag at this wavelength, whereas the point source is 13.16 mag at 10.65 µm. A background object would have to be extremely red to be compatible with both measurements.

Reanalysis of archival VISIR/NEAR data

We revisited archival VISIR/NEAR⁵¹,⁵² observations of Eps Ind A. We started from previously published reduced datacubes¹⁶, for which frames were calibrated and aligned. We removed low-quality frames, and then binned and spatially high-pass filtered the remaining frames. The NEAR sensitivity was optimized for the small separation (0.5"–3.0") around the stars. Ground-based thermal infrared observations require fast chopping to quickly subtract the varying sky background and excess low-frequency detector noise. With NEAR, chopping was done using the deformable secondary mirror of the VLT, which provides a maximum chop throw of about 4.5". The reduction with angular differential imaging (ADI) relies on subtracting images taken at a different field orientation, which leaves residuals of real sources in the azimuthal direction. As the separation of the imaged planet from the host star was close to the chop throw, these residuals masked the planet if ADI were used for data reduction. To avoid this effect, we simply co-added data without performing ADI, as the expected separation was well within the background-limited regime rather than the contrast-limited regime. We also performed a locally optimized combination of images (LOCI)⁵³ PCA ADI analysis with three principal components in a small patch around the expected companion location. This helped to improve the planet signal and suppress the thermal background. In each case, we also convolved the final images with a top hat to enhance any true signals in the data. Both images are shown in Fig. 2. To confirm the source, we also looked at the companion position relative to the off-axis PSFs.

We centred, aligned and derotated the frames relative to this position. We found that the planet signal was also present in the final averaged convolved frames relative to the off-axis PSFs, negating the possibility of the source being a speckle or a ghost from internal reflections. We assumed an error of half the full-width at half-maximum on the companion position, due to the high-spatial resolution and sampling of the NEAR. The planet was at a separation of 4.82 ± 0.16" and at a position angle of 40–45° anticlockwise of north. To estimate the flux, we used planet injection and recovery tests. As the planet signal was limited by the thermal background, we put a conservative range of (4–8) × 10⁻⁵ on contrast and a range of 0.18–0.35 mJy on flux. These values are also included in Extended Data Table 2.

Orbital constraints

Our goal in targeting Eps Ind A with JWST coronagraphic imaging was to detect a known, massive planet for which dynamical mass information was accessible. However, the imaged planet is inconsistent with the previously claimed planet in the system. The planet has a significantly larger mass and semimajor axis than the claimed Eps Ind Ab (refs. 2,8,9) and a strikingly different on-sky position than predicted. It is clear that the imaged planet imparts a significant, long-term, radial-velocity acceleration on the host star, although it was initially unclear whether this is the only planet in the system.

We attempted to fit the orbit of the imaged planet by assuming that it is, indeed, the only planet in the system and that it is responsible for all the dynamical measurements (radial-velocity data and astrometry) of the host star. We used radial-velocity data from the long camera (LC) and very long camera (VLC) of the coudé echelle spectrometer at the European Southern Observatory (ESO)⁵, from the ultraviolet and visual echelle spectrograph (UVES)² and from the high accuracy radial-velocity planet searcher (HARPS), using the data reduction presented in ref. 54. Eps Ind A has been targeted intensively with HARPS, and we binned any data points collected within a single night with HARPS. We treated the HARPS data as three separate instruments, to account for any possible baseline shifts during an instrument upgrade in 2015 and a telescope shutdown in 2020. This gave us 493 data points, from six independent instruments, covering from 3 November 1992 to 29 December 2021. We fitted orbits using orvaras⁵⁵, considering the 493 radial-velocity points, the two direct imaging epochs and the Hipparcos–Gaia astrometry of the host star, and using a host star mass prior of 0.76 ± 0.04 M☉. This process converged to an orbit with mass 6.31±0.60 M☉, semimajor axis 28.4±7.2 au, eccentricity 0.40±0.18 and inclination 103.7 ± 2.3° (Extended Data Fig. 6). Without the NEAR astrometry, we found an orbit with mass 6.17±0.65–0.57 M☉, semimajor axis 31±13 au, eccentricity 0.47±0.15 and inclination 104.0 ± 2.4°. These orbits are highly consistent with each other, further validating that the NEAR epoch is a positive redetection of the companion. These orbits are also consistent with all in-hand dynamical data. Curiously, several previous works had derived properties of the claimed planet Eps Ind Ab and found consistent results²,⁸,⁹, but these results are inconsistent with the planet observed in this work. This may be due to overfitting of the in-hand data, as fitting accurate orbits with insufficient orbital phase coverage is notoriously hard⁵⁶, or it may hint at another component in the system that biased the previous one-planet fits. The consistency between the fits with and without the NEAR astronomy points towards a one-planet system. We leave a detailed exploration of the discrepancy between the imaged planet and the previous radial-velocity solutions to a future work, but note that the dynamical mass derived in this work is valid only if Eps Ind Ab is the only massive planet present.

Stellar photometry

We used BT-NextGen models⁵⁷,⁵⁸ to predict the photometry of the host star Eps Ind A in the JWST/MIRI and VLT/NaCo filters. We fitted spectra to the in-hand photometry from TYCHO, Gaia, 2MASS, WISE, Spitzer/ MIPS and the photoconductor array camera and spectrometer
(PACS) onboard Herschel<sup>10,47,50,59,60</sup>. We allowed the effective temperature, log(g), metallicity and radius to vary and held the stellar distance fixed at 3.639 pc. We sampled the posterior with a Monte Carlo Markov chain (using emcee; ref. 61) to derive the best-fitting spectrum and uncertainty. The best-fitting values are effective temperature $T_{\text{eff}} = 4,760 \pm 15$ K, log(g) = $5.25^{+0.18}_{-0.34}$, [Fe/H] = 0.22 ± 0.12 and R = 0.679 ± 0.004 R<sub>☉</sub>, although we caution that these values are driven by photometry only, and spectroscopic fits of the stellar parameters should be used over these values when interpreting the metallicity and log(g) of the host star. We then derived model photometry and uncertainties by selecting samples from the Markov chain Monte Carlo (MCMC) posterior and integrating over the filter profiles. The best-fitting spectrum as well as the measured (green and pink) and model (blue) photometry of Eps Ind A is shown in Extended Data Fig. 7. We used this model to derive all companion photometry and flux limits in this work, including reconverting archival contrast curves to flux upper limits with this model.

The WISE observations of Eps Ind A appear to be unreliable. Eps Ind A is significantly saturated in the WISE W1, W2 and W3 observations (the WISE source catalogue<sup>60</sup> lists saturation fractions of w1sat = 0.253, w2sat = 0.239 and w3sat = 0.144). For this reason, we did not include the W1, W2 and W3 magnitudes directly in our fitting process. Our derived stellar magnitudes show some discrepancy with those used to derive previous L' limits<sup>17</sup>, which used the WISE W1 magnitude of 2.9 as a proxy for the L' magnitude of Eps Ind A. Our stellar fit indicates a somewhat brighter magnitude for the star of 2.119 ± 0.008 and 2.115 ± 0.009 for the W1 and L' filters respectively. The only published L-band magnitude for Eps Ind is L = 2.12 (ref. 62), consistent with our derived W1 and L' magnitudes.

## Data availability

All observational data in this work are available through public data archives. In particular, the JWST data were collected through the General Observer programme 2243 (PI E.C.M.) and are available through the Mikulski Archive for Space Telescopes (https://mast.stsci.edu). Spitzer data were collected through programmes 90 (PI M. Werner) and 41 (PI G. Rieke) and are available through the Spitzer Heritage Archive (https://irsa.ipac.caltech.edu/applications/Spitzer/SHA/). VLT data from VISIR/NEAR and SPHERE were collected through programmes 60.A-9107/103.201H (PI M. Meyer) and 110.25BR (PI E.C.M.), respectively, and are available through the ESO science archive (http://archive.eso.org/cms.html). Radial-velocity data and atmospheric models were taken from previously published literature, as indicated in the text.

## Code availability

Data reduction was carried out using the publicly available spaceKLIP package (https://github.com/kammerje/spaceKLIP), which also uses functions from the jwst pipeline (https://jwst-pipeline.readthedocs.io), the pyKLIP Python package (https://pyklip.readthedocs.io) and the webbpsf_ext package (https://webbpsf.readthedocs.io; https://github.com/JarronL/webbpsf_ext). Orbit fitting was carried out using the publicly available orvara package (https://github.com/t-brandt/orvara). We used various functions from astropy, pandas, matplotlib and seaborn in the analysis and to create figures.

**Acknowledgements** We thank M. Janson, J. Faherty and T. Brandt for discussions that aided the preparation of this paper. This work is based on observations made with the NASA/ESA/CSA JWST. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with JWST programme 2243. This work makes use of observations collected at the ESO (Programmes 110.25BR and 60.A-9107 (103.201H)). The NEAR project was made possible through contributions from the Breakthrough Foundation and Breakthrough Watch programme, as well as through contributions from the ESO and the NEAR consortium. This work makes use of observations made with the Spitzer Space Telescope, which was operated by the Jet Propulsion Laboratory, California Institute of Technology, under a contract with NASA. This research has made use of the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with NASA under the Exoplanet Exploration Program. Support for programme 2243 was provided by NASA through a grant from the Space Telescope Science Institute. Portions of this research used high-performance computing resources supported by the University of Arizona Technology and Research Initiative Fund (TRIF), University Information Technology Services (UITS), and Research, Innovation, and Impact (RII), and maintained by the University of Arizona's Research Technologies department. Part of this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with NASA. The contributions of D.P. were made within the National Centre of Competence in Research PlanetS framework supported by the Swiss National Science Foundation (Grant Nos. 51NF40_182901 and 51NF40_205606). This project has received funding from the European Research Council under the European Union's Horizon 2020 research and innovation programme (COBREX; Grant Agreement No. 885593). A.M.L. and F.P. received funding from Paris Sciences et Lettres University/l'Université Côte d'Azur (PSL/OCAV).

**Author contributions** E.C.M. designed the science programme, performed the analysis of JWST/MIRI, VLT/SPHERE, Spitzer/IRAC and Spitzer/MIPS data, carried out the statistical analysis of the background and companion hypotheses, contributed to orbit fitting, carried out the comparison to atmospheric and evolutionary models, reanalysed the stellar spectrum, led the interpretation of these results and wrote the manuscript. A.C. designed the MIRI observing sequence and supported the JWST/MIRI data reduction. P.P. and M.K.

<footnote n="33">33. Boccaletti, A. et al. The mid-infrared instrument for the James Webb Space Telescope. V. Predicted performance of the MIRI coronographs. *Publ. Astron. Soc. Pac.* **127**, 633 (2015).</footnote>

<footnote n="34">34. Boccaletti, A. et al. JWST/MIRI coronographic performances as measured on-sky. *Astron. Astrophys.* **667**, A165 (2022).</footnote>

<footnote n="35">35. Matthews, E. C. et al. A direct detection of the closest Jupiter analog with JWST/MIRI. *JWST Proposal Cycle 1*, 2243 (2021).</footnote>

<footnote n="36">36. Rouan, D., Riaud, P., Boccaletti, A., Clénet, Y. & Labeyrie, A. The four-quadrant phase-mask coronograph. I. Principle. *Publ. Astron. Soc. Pac.* **112**, 1479–1486 (2000).</footnote>

<footnote n="37">37. Smith, B. A. & Terrile, R. J. A circumstellar disk around β Pictoris. *Science* **226**, 1421–1424 (1984).</footnote>

<footnote n="38">38. Soummer, R., Hines, D. & Perrin, M. *Simulations of Target Acquisition with MIRI Four-quadrant Phase Mask Coronagraph (I)*, Report JWST-STSc-003063, SM-12 (STSI, 2013).</footnote>

<footnote n="39">39. Lajoie, C.-P. et al. Small-grid dithers for the JWST coronographs. *Proc. SPIE* **9904**, 99045K (2016).</footnote>

<footnote n="40">40. Rigby, J. et al. The science performance of JWST as characterized in commissioning. *Publ. Astron. Soc. Pac.* **135**, 048001 (2023).</footnote>

<footnote n="41">41. Kammerer, J. et al. Performance of near-infrared high-contrast imaging methods with JWST from commissioning. *Proc. SPIE* **12180**, 121803N (2022).</footnote>

<footnote n="42">42. Carter, A. L. et al. The JWST early release science program for direct observations of exoplanetary systems. I. High-contrast imaging of the exoplanet HIP 65426 b from 2 to 16 µm. *Astrophys. J. Lett.* **951**, L20 (2023).</footnote>

<footnote n="43">43. Bushouse, H. et al. JWST calibration pipeline. Zenodo https://doi.org/10.5281/zenodo.6984365 (2024).</footnote>

<footnote n="44">44. Wang, J. J. et al. pyKLIP: PSF subtraction for exoplanets and disks. Astrophysics Source Code Library, record ascl:1506.001 (ASCL, 2015).</footnote>

<footnote n="45">45. Ren, B., Pueyo, L., Zhu, G. B., Debes, J. & Duchêne, G. Non-negative matrix factorization: robust extraction of extended structures. *Astrophys. J.* **852**, 104 (2018).</footnote>

<footnote n="46">46. Krishanth P. M. S. et al. NMF-based GPU accelerated coronagraphy pipeline. *Proc. SPIE* **12680**, 1268021 (2023).</footnote>

<footnote n="47">47. Skrutskie, M. F. et al. The Two Micron All Sky Survey (2MASS). *Astron. J.* **131**, 1163–1183 (2006).</footnote>

<footnote n="48">48. Marengo, M. et al. Spitzer/Infrared array camera limits to planetary companions of Fomalhaut and Epsilon Eridani. *Astrophys. J.* **700**, 1647–1657 (2009).</footnote>

<footnote n="49">49. Backman, D. et al. Epsilon Eridani's planetary debris disk: structure and dynamics based on Spitzer and Caltech submillimeter observatory observations. *Astrophys. J.* **690**, 1522–1538 (2009).</footnote>

<footnote n="50">50. Gáspár, A., Rieke, G. H. & Balog, Z. The collisional evolution of debris disks. *Astrophys. J.* **768**, 25 (2013).</footnote>

<footnote n="51">51. Lagage, P. O. et al. Successful commissioning of VISIR: the mid-infrared VLT instrument. *Messenger* **117**, 12–16 (2004).</footnote>

<footnote n="52">52. Kasper, M. et al. NEAR: low-mass planets in a Cen with VISIR. *Messenger* **169**, 16–20 (2017).</footnote>

<footnote n="53">53. Lafrenière, D, Marois, C., Doyon, R., Nadeau, D. & Artigau, É. A new algorithm for point-spread function subtraction in high-contrast imaging: a demonstration with angular differential imaging. *Astrophys. J.* **660**, 770 (2007).</footnote>

<footnote n="54">54. Trifonov, T. et al. Public HARPS radial velocity database corrected for systematic errors. *Astron. Astrophys.* **636**, A74 (2020).</footnote>

<footnote n="55">55. Brandt, T. D. et al. orvara: an efficient code to fit orbits using radial velocity, absolute, and/or relative astrometry. *Astron. J.* **162**, 186 (2021).</footnote>

<footnote n="56">56. Lagrange, A. M. et al. Radial distribution of giant exoplanets at Solar System scales. *Astron. Astrophys.* **677**, A71 (2023).</footnote>

<footnote n="57">57. Allard, F., Homeier, D. & Freytag, B. Model atmospheres from very low mass stars to brown dwarfs. In *Proc. 16th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun* (eds Johns-Krull, C., Browning, M. K. & West, A. A.) **448**, 91 (2011).</footnote>

<footnote n="58">58. Allard, F., Homeier, D. & Freytag, B. Models of very-low-mass stars, brown dwarfs and exoplanets. *Philos. Trans. R. Soc. Lond. Ser. A* **370**, 2765–2777 (2012).</footnote>

<footnote n="59">59. Høg, E. et al. The Tycho-2 catalogue of the 2.5 million brightest stars. *Astron. Astrophys.* **355**, L27–L30 (2000).</footnote>

<footnote n="60">60. Cutri, R. M. et al. VizieR Online Data Catalog: WISE Preliminary Data Release (Cutri+ 2011) (ADS, 2012).</footnote>

<footnote n="61">61. Foreman-Mackey, D., Hogg, D. W., Lang, D. & Goodman, J. emcee: the MCMC hammer. *Publ. Astron. Soc. Pac.* **125**, 306 (2013).</footnote>

<footnote n="62">62. Morel, M. & Magnenat, P. UBVRIJKLMNH photoelectric photometric catalogue. (Magnetic tape). *Astron. Astrophys. Suppl. Ser.* **34**, 477 (1978).</footnote>
carried out the VISIR/NEAR data reduction. C.V.M. and M.W.P. provided a theoretical
interpretation of the photometric constraints and upper limits. S.K.P.M., E.D., J.H. and C.L.K.
provided an alternative JWST/MIRI data reduction. F.F., J.B., A.M.L. and F.P. contributed to the
orbit fitting and dynamical mass estimates for the target. M.J.B. reanalysed the VLT/NaCo
data. L.A.B. interpreted possible Galactic contaminants. I.J.M.C. and Th.H. contributed to
the interpretation of the results. D.P. contributed to the ground-based vetting of the
reference target. The work is based on a JWST proposal written by E.C.M. with contributions
from A.C., F.F., M.W.P., C.V.M., I.J.M.C., E.D. and J.B. All authors discussed the results and
commented on the manuscript.

**Funding** Open access funding provided by Max Planck Society.

**Competing interests** The authors declare no competing interests.

**Additional information**

**Correspondence and requests for materials** should be addressed to E. C. Matthews.

**Peer review information** *Nature* thanks Ignas Snellen and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

**Reprints and permissions information** is available at http://www.nature.com/reprints.
<figure><img src="image_8.png" /><figcaption>Extended Data Fig. 1 | Full field-of-view JWST/MIRI coronagraphic images of Eps Ind A in the 10.65 µm filter. Two visual companions to Eps Ind A, and one visual companion to DI Tuc, are identified: (1) is Eps Ind Ab, (2) is the</figcaption></figure>

background star Gaia DR3 6412595290591430784 (also identified in the Spitzer 8 µm observations of the background field) and (3) is the background star Gaia DR3 6411654761473726464 in the field of the DI Tuc observations.
<figure><img src="image_9.png" /><figcaption>Extended Data Fig. 2 | Contrast sensitivity of the MIRI observations. Curves are 5σ sensitivities calculated with spaceKLIP for each coronagraphic filter.</figcaption></figure>
<figure><img src="image_10.png" /></figure>

**Extended Data Fig. 3** | Spitzer IRAC images of Eps Ind A. The three columns show the science-ready Eps Ind A image from the Spitzer archive (left), an RDI processed image with Eps Eri subtracted from Eps Ind A (centre), and the RDI processed image zoomed at the JWST epoch position of Eps Ind A (right). For the RDI processed image, purple (negative) sources are in the field of the reference star, and orange (positive) sources are in the field of Eps Ind A. Blue

and green circles have radii 10", and are centred on the Spitzer and JWST epoch positions of Eps Ind A respectively. If Eps Ind Ab were a background object, it would be ~4" from the JWST epoch position of Eps Ind A. This is well within the blue circle. One source is identified within this circle at 3.6 µm, 4.5 µm and marginally at 5.8 µm: this is a background star also identified in the JWST/MIRI images (see text for details).
<figure><img src="image_11.png" /><figcaption>**Extended Data Fig. 4** | Spitzer MIPS 24µm and 70µm images of Eps Ind A.
Blue and green circles have radii 20", and are centred on the Spitzer and JWST epoch positions of Eps Ind A respectively; no sources are detected at the</figcaption></figure>

<figure><img src="image_12.png" /></figure>

background position (the JWST epoch position of Eps Ind A). The 70 µm data
cover only a narrow strip around the Eps Ind A location, and white regions at the
edge of the box indicate regions outside the telescope field of view.
<figure><img src="image_13.png" /></figure>

**Extended Data Fig. 5** | Spectrum (measured flux and upper limits) of the point source if it were a background object. These upper limits apply only if the point source were a stationary background at RA = 22:03:33.17 and Dec = -56:48:06.0. Constraints for the planetary case (which require high-contrast facilities to resolve a source that moves with same proper

motion as the host star) are given in Fig. 3. Green circles indicate the measured JWST photometry, while pink triangles indicate upper limits from Gaia, 2MASS, Spitzer/IRAC and Spitzer/MIPS. Most stationary background contaminants would not reproduce the 10.65 µm and 15.50 µm flux while evading detection at other wavelengths.
<figure><img src="image_14.png" /><figcaption>Extended Data Fig. 6 | Orbit fits including RV, Hipparcos-Gaia astrometry, and both imaging points (left) or the JWST imaging only (right), assuming Eps Ind Ab is the only massive planet in the system. Rows represent (1) the on-sky orbit, (2) the RV of the target (3) the proper motion in right ascension of</figcaption></figure>

the host star and (4) the proper motion in declination of the host star. Orbits are very similar in both cases, and explain all in-hand data. RVs are from LC (green), VLC (pink), HARPS03 (cyan), UVES (yellow), HARPS15 (blue) and HARPS20 (pink).
<figure><img src="image_15.png" /></figure>

**Extended Data Fig. 7** | **Best-fit spectrum for Eps Ind A.** The spectrum is shown in black and integrated to R-100. Measured photometry is shown in green if included in the fit, and pink otherwise, while the model photometry is indicated

in blue. The WISE W1, W2, and W3 channels are saturated, and in particular the W1 and W3 magnitudes⁶⁰ are several σ fainter than the best-fit model.
Extended Data Table 1 | JWST/MIRI observations for Eps Ind A

<table><thead><tr><th>Target</th><th>Telescope Coordinates</th><th>Filter</th><th>Readout</th><th>Groups/Int</th><th>Ints/Exp</th><th>Dithers</th><th>Total Time [s]</th></tr></thead><tbody><tr><td>background</td><td>22:17:54.000 -57:28:50.00</td><td>F1550C</td><td>FASTR1</td><td>811</td><td>14</td><td>2</td><td>5448.885</td></tr><tr><td>background</td><td>22:17:54.000 -57:28:50.00</td><td>F1550C</td><td>FASTR1</td><td>1168</td><td>14</td><td>2</td><td>7844.726</td></tr><tr><td>background</td><td>22:17:54.000 -57:28:50.00</td><td>F1065C</td><td>FASTR1</td><td>89</td><td>122</td><td>2</td><td>5262.893</td></tr><tr><td>background</td><td>22:17:54.000 -57:28:50.00</td><td>F1065C</td><td>FASTR1</td><td>128</td><td>122</td><td>2</td><td>7543.688</td></tr><tr><td>DI Tuc</td><td>22:16:08.195 -57:34:05.06</td><td>F1065C</td><td>FASTR1</td><td>89</td><td>122</td><td>5</td><td>13157.234</td></tr><tr><td>DI Tuc</td><td>22:16:08.195 -57:34:05.06</td><td>F1550C</td><td>FASTR1</td><td>811</td><td>14</td><td>5</td><td>13622.213</td></tr><tr><td>Eps Ind A</td><td>22:03:33.001 -56:48:09.13</td><td>F1550C</td><td>FASTR1</td><td>1168</td><td>14</td><td>1</td><td>3922.363</td></tr><tr><td>Eps Ind A</td><td>22:03:33.001 -56:48:09.13</td><td>F1065C</td><td>FASTR1</td><td>128</td><td>122</td><td>1</td><td>3771.844</td></tr></tbody></table>

All observations were carried out sequentially on 3 July 2023, and have one exposure per dither position.
Extended Data Table 2 | Astrometry, photometry and upper limits for Eps Ind Ab

<table><thead><tr><th>Instrument</th><th>Filter</th><th>Date</th><th>ρ<br>[arcsec]</th><th>θ<br>[deg]</th><th>Flux<br>[W/m²/μm]</th><th>5σ Limit<br>[W/m²/μm]</th></tr></thead><tbody><tr><td>JWST/MIRI</td><td>F1065C</td><td>2023-07-03</td><td>4.095±0.010<sup>(a)</sup></td><td>38.17±0.44<sup>(a)</sup></td><td>(4.98 ± 0.15) × 10<sup>-18</sup></td><td>–</td></tr><tr><td>JWST/MIRI</td><td>F1550C</td><td>2023-07-03</td><td>4.114±0.010</td><td>37.39±0.43</td><td>(6.55 ± 0.20) × 10<sup>-18</sup></td><td>–</td></tr><tr><td rowspan="3">VLT/VISIR</td><td rowspan="3">NEAR</td><td>2019-09-14</td><td rowspan="3">4.82 ± 0.164</td><td rowspan="3">42.5 ± 2.5</td><td rowspan="3">(6.6 ± 2.2) × 10<sup>-18</sup></td><td rowspan="3">–</td></tr><tr><td>2019-09-15</td></tr><tr><td>2019-09-17</td></tr><tr><td rowspan="4">VLT/NaCo</td><td rowspan="4">L'</td><td>2018-10-12</td><td rowspan="4">–</td><td rowspan="4">–</td><td rowspan="4">–</td><td rowspan="4">3.78×10<sup>-18</sup></td></tr><tr><td>2018-10-26</td></tr><tr><td>2018-11-03</td></tr><tr><td>2018-11-04</td></tr><tr><td rowspan="2">VLT/NaCo</td><td rowspan="2">NB4.05</td><td>2008-10-31</td><td rowspan="2">–</td><td rowspan="2">–</td><td rowspan="2">–</td><td rowspan="2">4.16×10<sup>-17</sup></td></tr><tr><td>2008-09-02</td></tr></tbody></table>

JWST and VLT/VISIR constraints are from this work; VLT/NaCo constraints are from ref. 17 (L', their "special LOCI" reduction with our stellar calibration from Section 2.5) and ref. 18 (NB4.05, we quote a median sensitivity of the observations between 3" and 8"). Values for multi-epoch visits (VLT/VISIR, VLT/NaCo) correspond to a combined analysis of all visits. a: Note that for F1065C the star is not well centred behind the coronagraph, and the astrometry should be considered unreliable.