# Probing the innermost region of the AU Microscopii debris disc

A. Gallenne, C. Desgrange, J. Milli, J. Sanchez-Bermudez, G. Chauvin, S.

Kraus, J. H. Girard, A. Boccaletti

► To cite this version:

A. Gallenne, C. Desgrange, J. Milli, J. Sanchez-Bermudez, G. Chauvin, et al. Probing the innermost region of the AU Microscopii debris disc. Astronomy and Astrophysics - A&A, 2022, 665, 10.1051/0004-6361/202244226. insu-03860310

HAL Id: insu-03860310

https://insu.hal.science/insu-03860310

Submitted on 19 Nov 2022

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
# Probing the innermost region of the AU Microscopii debris disc*

A. Gallenne<sup>1,2</sup>●, C. Desgrange<sup>3,4</sup>●, J. Milli<sup>3</sup>●, J. Sanchez-Bermudez<sup>4,5</sup>●, G. Chauvin<sup>2,3</sup>●, S. Kraus<sup>6</sup>,
J. H. Girard<sup>7</sup>●, and A. Boccaletti<sup>8</sup>●

1. Universidad de Concepción, Departamento de Astronomía, Casilla 160-C, Concepción, Chile
e-mail: agallenne@astro-udec.cl

2. Unidad Mixta Internacional Franco-Chilena de Astronomía (CNRS UMI 3386), Departamento de Astronomía, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile

3. Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France

4. Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany

5. Instituto de Astronomía, Universidad Nacional Autónoma de México, Apdo. Postal 70264, Ciudad de México 04510, Mexico

6. School of Physics and Astronomy, University of Exeter, Stocker Road, Exeter EX4 4QL, UK

7. Space Telescope Science Institute, Baltimore, MD 21218, USA

8. LESIA, Observatoire de Paris, PSL Research Univ., CNRS, Univ. Paris Diderot, Sorbonne Paris Cité(c), UPMC Paris 6, Sorbonne Univ., 5 place Jules Janssen, 92195 Meudon, France

Received 9 June 2022 / Accepted 29 June 2022

## ABSTRACT

**Context.** AU Mic is a young and nearby M-dwarf star harbouring a circumstellar debris disc and one recently discovered planet on an eight-day orbit. Large-scale structures within the disc were also discovered and are moving outwards at high velocity.

**Aims.** We aim to study this system with the highest spatial resolution in order to probe the innermost regions and to search for additional low-mass companions or set detection limits.

**Methods.** The star was observed with two different high-angular resolution techniques probing complementary spatial scales. We obtained new $K_s$-band sparse aperture masking observations with VLT/SPHERE, which we combined with data from VLT/NACO, VLTI/PIONIER and VLTI/GRAVITY.

**Results.** We did not detect additional close companions within the separation range 0.02–7 au from the parent star. We determined magnitude upper limits for companions of $H \sim 9.8$ mag within 0.02–0.5 au, $K_s \sim 11.2$ mag within 0.4–2.4 au, and $L' \sim 10.7$ mag within 0.7–7 au. Using theoretical isochrones, we converted these magnitudes into upper limits on the mass of ~17 $M_{\text{Jup}}$, ~12 $M_{\text{Jup}}$, and ~9 $M_{\text{Jup}}$, respectively. The PIONIER observations also allowed us to determine the angular diameter of AU Mic, $\theta_{\text{LD}} = 0.825 \pm 0.033_{\text{stat}} \pm 0.038_{\text{sys}}$ mas, which converts to a linear radius $R = 0.862 \pm 0.052 R_{\odot}$ when combined with the Gaia parallax.

**Conclusions.** We did not detect the newly discovered planets orbiting AU Mic ($M < 0.2 M_{\text{Jup}}$), but we derived upper limit masses for the innermost region of AU Mic. We do not have any detection with a significance beyond $3\sigma$, the most significant signal with PIONIER being $2.9\sigma$ and that with SPHERE being $1.6\sigma$. We applied the pyMESS2 code to estimate the detection probability of companions by combining radial velocities, multi-band SPHERE imaging, and our interferometric detection maps. We show that 99% of the companions down to ~0.5 $M_{\text{Jup}}$ can be detected within 0.02 au or 1 $M_{\text{Jup}}$ down to 0.2 au. The low-mass planets orbiting at $\le 0.11$ au (≤11 mas) from the star will not be directly detectable with the current adaptive optics (AO) and interferometric instruments because of its close orbit and very high contrast (~$10^{-10}$ K). It will also be below the angular resolution and contrast limit of the next Extremely Large Telescope Infrared (ELT IR) imaging instruments.

**Key words.** instrumentation: high angular resolution – planets and satellites: formation – stars: individual: AU Mic

## 1. Introduction

AU Microscopii (AU Mic, HD 197481, $V = 8.63$, $K_s = 4.53$, $L' = 4.32$) is a young nearby M-type star, and a member of the 18.5 Myr old $\beta$ Pictoris Moving Group (Miret-Roig et al. 2020) located at $d = 9.714 \pm 0.002$ pc (Gaia Collaboration 2021, 2016; Lindegren et al. 2021). It hosts one of the best-studied edge-on debris discs and has been observed at X-ray to millimetre wavelengths (see e.g. Kalas et al. 2004; Quillen et al. 2007; Graham et al. 2007; Schneider & Schmitt 2010; Wilner et al. 2012; MacGregor et al. 2013; Matthews et al. 2015; Boccaletti et al. 2015; Daley et al. 2019; Esposito et al. 2020) with both space- and ground-based telescopes. Early

modelling suggested the existence of a ring of colliding parent bodies located at ~40 au, generating micron-sized dust particles later put on eccentric orbits through stellar wind and radiation forces, or migrating inwards through Poynting Robertson drag (Augereau & Beust 2006; Strubbe & Chiang 2006). Submillimetre ALMA observations were later used to infer that this birth ring lies between ~8 and 40 au (MacGregor et al. 2013; Matthews et al. 2015), although the inner radius is poorly constrained. An additional halo component made of particles down to submicron sizes is seen both in scattered light and thermal imaging and extends up to ~210 au (Kalas et al. 2004; Krist et al. 2005; Liu 2004; Schneider et al. 2014) with a complex morphology (Fitzgerald et al. 2007). Disc features such as gaps or asymmetries created by planetesimal or infant planets make debris discs good laboratories with which to study planet formation and evolution. AU Mic is a particularly interesting

<footnote>* Based on observations made with ESO telescopes at Paranal and La Silla observatory under program IDs 0101.C-0218(A), 087.C-0450(A), 087.C-0450(B) 087.C-0750(A), 088.C-0358(A) and 60.A-9165(A).</footnote>
system as several morphological features are present in the disc. The disc is asymmetrical with the northwest side brighter than the southwest side, and present some substructures possibly created by larger unseen orbiting bodies (Liu 2004). Brightness enhancements have been detected and attributed to dust-density variations in the range 25–40 au. A prominent bump at a projected separation of ~13 au was also identified on the southeast side above the disc plane. The disc also exhibits time variability with large-scale arc-like structures located at projected separations of ~5–10 au on the southeast side. Those arches are moving away from the star with projected speeds that are inconsistent with a Keplerian velocity, and some even reach speeds higher than the escape velocity (Boccaletti et al. 2015, 2018). One of the scenarios to explain those features would be an unseen point-source parent body releasing dust sequentially (Sezestre et al. 2017) or dust avalanches (Chiang & Fung 2017). Recent 1.3 mm ALMA observations (Daley et al. 2019) set upper limits on the largest planetesimals stirring the disc, and ruled out a perturbing body larger than a Neptune-size object within the outer belt.

Several high-contrast and high-spatial resolution observations have been reported, either using the *Hubble* Space Telescope (HST) or extreme adaptive optics (AO) systems, but none of these have yet been able to directly detect a planet-size body. Krist et al. (2005) performed the first coronagraphic observations of AU Mic but only suggested an unseen companion that would perturb the disc to explain the various detected features. Subsequent coronagraphic observations revealed more and finer details about the disc dynamic and dust distribution (see e.g. Wang et al. 2015; Boccaletti et al. 2015, 2018; Lomax et al. 2018), but still no detection of this hypothetical companion. Gauchet et al. (2016) observed the innermost region of the AU Mic disc in the *L'* band using the sparse aperture masking (SAM) technique, which provides a spatial resolution below the diffraction limit (~λ/2D). These latter authors did not find any companion brighter than Δ*L'* = 5.7 mag (5σ limit) within 30–500 mas (0.3–5 au). This would correspond to a companion mass of larger than 10 *M*<sub>jup</sub>. From angular differential imaging with a vortex coronagraph, Launhardt et al. (2020) reported a similar 5σ contrast limit in the *L'* band between 100 mas and 250 mas (6.2 mag). Detection limits were also provided by Lannier et al. (2017) who combined radial velocity (RV) measurements and direct imaging. These authors constrained the presence of planets with masses ≥2 *M*<sub>jup</sub> from a fraction of an astronomical unit (au) of to several tens of au.

Recently, a super-Neptune planet transiting AU Mic was detected by Plavchan et al. (2020) using NASA's Transiting Exoplanet Survey Satellite (TESS). This planet orbits the star in 8.46 days with a semi-major axis of 0.066 au (6.74 mas). Such a small separation partly explains why it was not detected before with high-contrast imaging, being hidden behind the coronagraph. From RV analysis, these latter authors inferred an upper limit for the planet of 0.18 *M*<sub>jup</sub> (= 3.4 *M*<sub>nep</sub>). Klein et al. (2021) confirmed the planet detection at a 3.9σ level with the SPIROU instrument, and inferred a mass of 0.054 *M*<sub>jup</sub>. Plavchan et al. (2020) suspected a possible second transiting planet, which was later confirmed by Martioli et al. (2021) with new TESS observations. This second planet of ~0.08 *M*<sub>jup</sub> orbits its star with a period of 18.9 d and a semi-major axis of 0.11 au (~11.3 mas).

In this paper, we aim at setting detection limits for additional orbiting companions within the range 0.04–7 au (4–730 mas) using high-angular-resolution observations, including infrared (IR) long-baseline interferometry and SAM. Section 2 presents new SAM observations acquired in 2018 with the SPHERE/IRDIS instrument mounted at the Very

Large Telescope (VLT), together with literature data from the VLTI/GRAVITY and VLT/NACO instruments. We present the details of our search for a companion and estimate our detection limits in Sect. 3 using the CANDID algorithm. We then discuss our results and provide conclusions in Sect. 4.

## 2. Observations and data reduction

### 2.1. SPHERE/SAM observations

We observed AU Mic with the extreme AO instrument SPHERE (Spectro-Polarimetric High-contrast Exoplanet REsearch, Beuzit et al. 2019) and the IRDIS/SAM mode (Cheetham et al. 2016). This is an interferometric method (Tuthill et al. 2000) which uses a mask with non-redundant holes placed in the pupil plane. Each pair of subapertures (or baseline) will form an interference fringe pattern. Observables that are closure phases (sum of the phases between three baselines forming a triangle) and squared visibilities (fringe contrast) can then be extracted, allowing us to reach an angular resolution down to $\lambda/2D$ ($D$ is the diameter of the primary mirror).

The observations were carried out on 2018 August 18 with the IRDIS (Dohlen et al. 2008) in the *Ks*-band filter ($\lambda$ = 2.10 µm, $\Delta\lambda$ = 0.30 µm) and a seven-hole aperture mask. This provides 21 baselines, that is 21 fringes with a unique spatial frequency and direction. Data were acquired in a data cube mode with tens of frames of 1 s exposure time. Similar sequences are executed on an empty sky region for background subtraction and on a calibrator star. The full log is summarised in Table 1. AU Mic observations were interleaved with the calibrator star HD 196919 (K4III, $R$ = 6.94 mag, $K$ = 4.31 mag, $\theta_{LD}$ = 0.717 ± 0.023 mas) in order to monitor the interferometric transfer function and calibrate the AU Mic visibilities. HD 196919 was chosen with the *SearchCal* software<fnref n="1" /> to be as close as possible to AU Mic in the sky and with a similar brightness.

Our data reduction steps were as follows. First, each image was corrected for detector cosmetics using a standard procedure, including sky subtraction, flat-fielding, and bad-pixel correction. Figure 1 shows an example of the interferogram and power spectrum of AU Mic. Subwindows of 200 × 200 px (2.45" × 2.45") of the two images (left and right channels of the IRDIS detector) in each frame of the cubes were then extracted, stored in new cubes of 450 interferograms each, and recentred at a subpixel level. To extract complex visibilities, we used our own Python algorithm based on the formalism detailed in Tuthill et al. (2000). Briefly, for each cube complex visibilities were computed frame by frame from the Fourier transform of the interferograms. Squared visibilities ($V^2$) were then calculated from the squared modulus, while the bispectrum was obtained by multiplying three complex visibilities corresponding to a triangle of holes in the mask. Closure phases (*CPs*) were then obtained from the argument of the bispectrum. Calibrated $V^2$ were obtained by dividing the AU Mic visibilities by those of the calibrator, while calibrated *CPs* were derived by subtracting those of the calibrator. The transfer functions are displayed in Figs. A.1 and A.2. We have a total of 294 $V^2$ and 490 *CPs* calibrated measurements, which are displayed in Fig. B.1.

Our SAM observations with SPHERE allow us to probe the spatial scale ~0.6–5.3 au (60–550 mas) from the star, with the lower limit set by the longest baselines and the upper

<footnote n="1">¹ https://www.jmmc.fr/english/tools/proposal-preparation/search-cal/</footnote>
<table><caption>**Table 1.** Log of our SAM observations with SPHERE/IRDIS and long-baseline interferometry with GRAVITY, together with some information about the observation of the literature data.</caption><thead><tr><th>Star</th><th>Date</th><th>MJD (day)</th><th>DIT (s)</th><th>NDIT</th><th>NEXP</th><th>Sp. Chan.</th><th>⟨Seeing⟩ (′′)</th><th>⟨τ<sub>0</sub>⟩ (ms)</th><th>Baselines</th><th>#</th></tr></thead><tbody><tr><td colspan="11" style="text-align:center; font-weight:bold;">SPHERE</td></tr><tr><td>AU Mic</td><td>2018-08-12</td><td>58342.22527</td><td>1</td><td>90</td><td>5</td><td>1</td><td>0.51</td><td>6.1</td><td>—</td><td>—</td></tr><tr><td>Sky</td><td>2018-08-12</td><td>58342.23169</td><td>1</td><td>10</td><td>1</td><td>1</td><td>0.47</td><td>6.0</td><td>—</td><td>—</td></tr><tr><td>AU Mic</td><td>2018-08-12</td><td>58342.23624</td><td>1</td><td>90</td><td>5</td><td>1</td><td>0.53</td><td>5.5</td><td>—</td><td>—</td></tr><tr><td>HD 196919</td><td>2018-08-12</td><td>58342.24860</td><td>1</td><td>40</td><td>6</td><td>1</td><td>0.50</td><td>6.0</td><td>—</td><td>—</td></tr><tr><td>Sky</td><td>2018-08-12</td><td>58342.25213</td><td>1</td><td>10</td><td>1</td><td>1</td><td>0.52</td><td>5.8</td><td>—</td><td>—</td></tr><tr><td>AU Mic</td><td>2018-08-12</td><td>58342.26256</td><td>1</td><td>30</td><td>15</td><td>1</td><td>0.48</td><td>5.8</td><td>—</td><td>—</td></tr><tr><td>Sky</td><td>2018-08-12</td><td>58342.26847</td><td>1</td><td>10</td><td>1</td><td>1</td><td>0.51</td><td>5.1</td><td>—</td><td>—</td></tr><tr><td>AU Mic</td><td>2018-08-12</td><td>58342.27387</td><td>1</td><td>30</td><td>15</td><td>1</td><td>0.51</td><td>5.6</td><td>—</td><td>—</td></tr><tr><td>HD 196919</td><td>2018-08-12</td><td>58342.28662</td><td>1</td><td>40</td><td>6</td><td>1</td><td>0.54</td><td>6.4</td><td>—</td><td>—</td></tr><tr><td>Sky</td><td>2018-08-12</td><td>58342.29015</td><td>1</td><td>10</td><td>1</td><td>1</td><td>0.51</td><td>6.4</td><td>—</td><td>—</td></tr><tr><td>AU Mic</td><td>2018-08-12</td><td>58342.29770</td><td>1</td><td>30</td><td>12</td><td>1</td><td>0.47</td><td>6.5</td><td>—</td><td>—</td></tr><tr><td>Sky</td><td>2018-08-12</td><td>58342.30257</td><td>1</td><td>10</td><td>1</td><td>1</td><td>0.41</td><td>7.2</td><td>—</td><td>—</td></tr><tr><td colspan="11" style="text-align:center; font-weight:bold;">NACO</td></tr><tr><td>AU Mic</td><td>2011-08-31</td><td>55804.15287</td><td>0.04</td><td>600</td><td>104</td><td>1</td><td>0.82</td><td>3.1</td><td>—</td><td>—</td></tr><tr><td colspan="11" style="text-align:center; font-weight:bold;">PIONIER</td></tr><tr><td>AU Mic</td><td>2011-08-09</td><td>55782.30522</td><td>0.21</td><td>100</td><td>1</td><td>3</td><td>0.64</td><td>6.9</td><td>A1-G1-I1-K0</td><td>1</td></tr><tr><td>AU Mic</td><td>2011-08-09</td><td>55782.31373</td><td>0.21</td><td>200</td><td>5</td><td>3</td><td>0.78</td><td>5.7</td><td>A1-G1-I1-K0</td><td>1</td></tr><tr><td>AU Mic</td><td>2011-08-09</td><td>55782.32517</td><td>0.21</td><td>100</td><td>1</td><td>3</td><td>0.66</td><td>6.8</td><td>A1-G1-I1-K0</td><td>1</td></tr><tr><td>AU Mic</td><td>2011-08-09</td><td>55782.34457</td><td>0.21</td><td>100</td><td>1</td><td>3</td><td>0.63</td><td>7.1</td><td>A1-G1-I1-K0</td><td>1</td></tr><tr><td>AU Mic</td><td>2011-11-03</td><td>55868.01837</td><td>0.39</td><td>100</td><td>5</td><td>7</td><td>0.94</td><td>2.3</td><td>D0-G1-H0-I1</td><td>2</td></tr><tr><td>AU Mic</td><td>2011-11-03</td><td>55868.03855</td><td>0.39</td><td>100</td><td>5</td><td>7</td><td>0.64</td><td>3.4</td><td>D0-G1-H0-I1</td><td>2</td></tr><tr><td>AU Mic</td><td>2011-11-03</td><td>55868.05373</td><td>0.39</td><td>100</td><td>5</td><td>7</td><td>0.97</td><td>2.3</td><td>D0-G1-H0-I1</td><td>2</td></tr><tr><td>AU Mic</td><td>2011-11-03</td><td>55868.06742</td><td>0.39</td><td>100</td><td>5</td><td>7</td><td>1.01</td><td>2.2</td><td>D0-G1-H0-I1</td><td>2</td></tr><tr><td>AU Mic</td><td>2011-11-03</td><td>55868.07997</td><td>0.39</td><td>100</td><td>5</td><td>7</td><td>1.17</td><td>1.9</td><td>D0-G1-H0-I1</td><td>2</td></tr><tr><td>AU Mic</td><td>2011-11-03</td><td>55868.09390</td><td>0.39</td><td>100</td><td>5</td><td>7</td><td>0.81</td><td>2.8</td><td>D0-G1-H0-I1</td><td>2</td></tr><tr><td>AU Mic</td><td>2011-11-03</td><td>55868.10677</td><td>0.39</td><td>100</td><td>5</td><td>7</td><td>0.99</td><td>2.2</td><td>D0-G1-H0-I1</td><td>2</td></tr><tr><td>AU Mic</td><td>2011-11-03</td><td>55868.12090</td><td>0.39</td><td>100</td><td>5</td><td>7</td><td>0.99</td><td>2.3</td><td>D0-G1-H0-I1</td><td>2</td></tr><tr><td>AU Mic</td><td>2013-08-09</td><td>56513.10329</td><td>0.39</td><td>100</td><td>5</td><td>3</td><td>0.92</td><td>2.5</td><td>A1-B2-C1-D0</td><td>3</td></tr><tr><td>AU Mic</td><td>2013-08-09</td><td>56513.11264</td><td>0.39</td><td>100</td><td>5</td><td>3</td><td>0.76</td><td>3.0</td><td>A1-B2-C1-D0</td><td>3</td></tr><tr><td>AU Mic</td><td>2013-08-09</td><td>56513.12271</td><td>0.39</td><td>100</td><td>5</td><td>3</td><td>0.75</td><td>3.0</td><td>A1-B2-C1-D0</td><td>3</td></tr><tr><td>AU Mic</td><td>2013-08-11</td><td>56515.07933</td><td>0.39</td><td>100</td><td>5</td><td>3</td><td>1.24</td><td>3.5</td><td>A1-B2-C1-D0</td><td>4</td></tr><tr><td>AU Mic</td><td>2013-08-11</td><td>56515.08979</td><td>0.39</td><td>100</td><td>5</td><td>3</td><td>1.42</td><td>4.0</td><td>A1-B2-C1-D0</td><td>4</td></tr><tr><td>AU Mic</td><td>2013-08-11</td><td>56515.09935</td><td>0.39</td><td>100</td><td>5</td><td>3</td><td>1.12</td><td>4.4</td><td>A1-B2-C1-D0</td><td>4</td></tr><tr><td>AU Mic</td><td>2014-10-12</td><td>56942.09731</td><td>0.38</td><td>100</td><td>5</td><td>3</td><td>1.42</td><td>1.7</td><td>A1-B2-C1-D0</td><td>5</td></tr><tr><td>AU Mic</td><td>2014-10-12</td><td>56942.10788</td><td>0.38</td><td>100</td><td>5</td><td>3</td><td>1.32</td><td>1.8</td><td>A1-B2-C1-D0</td><td>5</td></tr><tr><td>AU Mic</td><td>2014-10-12</td><td>56942.11896</td><td>0.38</td><td>100</td><td>5</td><td>3</td><td>1.61</td><td>1.5</td><td>A1-B2-C1-D0</td><td>5</td></tr><tr><td>AU Mic</td><td>2014-10-12</td><td>56942.13064</td><td>0.38</td><td>100</td><td>5</td><td>3</td><td>1.48</td><td>1.6</td><td>A1-B2-C1-D0</td><td>5</td></tr><tr><td colspan="11" style="text-align:center; font-weight:bold;">GRAVITY</td></tr><tr><td>AU Mic</td><td>2016-06-16</td><td>57555.43405</td><td>30</td><td>10</td><td>5</td><td>1742</td><td>1.18</td><td>3.7</td><td>A0-G1-J2-K0</td><td>—</td></tr><tr><td>AU Mic</td><td>2016-06-16</td><td>57555.44252</td><td>30</td><td>10</td><td>5</td><td>1742</td><td>1.12</td><td>3.0</td><td>A0-G1-J2-K0</td><td>—</td></tr><tr><td>AU Mic</td><td>2016-06-16</td><td>57555.45238</td><td>30</td><td>10</td><td>5</td><td>1742</td><td>1.35</td><td>3.0</td><td>A0-G1-J2-K0</td><td>—</td></tr></tbody></table>

**Notes.** Sp. Chan. is the number of spectral channel, ⟨Seeing⟩ is the average seeing as provided by the Differential Motion Image Monitor (DIMM) at 0.5 µm, and ⟨τ₀⟩ the coherence time provided by the MASS-DIMM and Baselines is the telescope configuration. The calibrators used are HD 197540 (G8III, V = 6.50 mag, H = 4.32 mag, θ₉₉ = 0.645 ± 0.046 mas), HD 201865 (K2III, V = 7.14 mag, H = 4.04 mag, θ₉₉ = 0.793 ± 0.075 mas) and HD 202775 (K4/5III, V = 7.49 mag, H = 4.05 mag, θ₉₉ = 0.844 ± 0.072 mas) for #1, HD 202293 (K5III, V = 7.92 mag, H = 3.74 mag, θ₉₉ = 0.994 ± 0.088 mas), HD 196387 (K4III, V = 7.28 mag, H = 3.52 mag, θ₉₉ = 1.110 ± 0.106 mas) and HD 202775 for #2, HD 197540, HD 194432 (K1III, V = 7.16 mag, H = 4.87 mag, θ₉₉ = 0.555 ± 0.013 mas) and HD 200750 (K1III, V = 6.81 mag, H = 4.53 mag, θ₉₉ = 0.627 ± 0.047 mas) for #3, HD 197540, HD 194432 and HD 200750 for #4, and HD 196387 for #5.

limit set by the bandwidth smearing. To improve the detection limits on a wider angular scale, we collected additional high-angular-resolution data from the VLTI/GRAVITY (Eisenhauer et al. 2011), VLTI/PIONIER (Le Bouquin et al. 2011) and VLT / NACO (Rousset et al. 2003; Lenzen et al. 2003) instruments.

## 2.2. NACO/SAM observations

We retrieved the reduced data published by Gauchet et al. (2016). These are also SAM mode observations and were made on 2011 August 04 with NACO in the L'-band filter (λ = 3.80 µm, Δλ = 0.62 µm) with a seven-hole mask. The observation log
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Interferogram of AU Mic trimmed within 2.45" × 2.45" (left image) and its power spectrum (squared visibilities, right image).</figcaption></figure>

from Gauchet et al. (2016) is reported in Table 1. Data reduction (flat-field, bad-pixels, and sky subtraction) were performed in the same way as for the SPHERE/SAM observations (see more details in Gauchet et al. 2016). The direct fringe-fitting method was applied instead of Fourier transform as we did previously; however both methods provide similar results. To calibrate the data, Gauchet et al. (2016) used HD 197339 (K2III, R = 6.47 mag, K = 4.47 mag, $\theta_{LD}$ = 0.646 ± 0.021 mas). These data allow to probe the spatial scale ~0.7–7 au (70–730 mas). Although these authors already published detection limits with this data set, we included it so that we can re-analyse the data with our CANDID tool, and perform a uniform analysis with the other datasets. The observations contain 546 $V^2$ and 910 CPs measurements, and are displayed in Fig. B.1.

## 2.3. VLTI/PIONIER observations

PIONIER calibrated data were retrieved from the Optical interferometry Database (OiDB<fnref n="2" />) developed by the Jean-Marie Marotti Center (JMMC<fnref n="3" />). PIONIER (Le Bouquin et al. 2011) is an interferometric instrument at the VLTI, combining the light coming from the four auxiliary telescopes (AT) in the H band, either in a broad band mode or with a low spectral resolution, where the light is dispersed across a few spectral channels.

The observations were performed on several dates from 2011 to 2014 for different observing programs. Data before 2014 were collected when PIONIER was a visitor instrument at the VLTI, and therefore the calibrated data are only available in the OiDB (not in the ESO archive). We listed these observations in Table 1. To monitor the instrumental and atmospheric contributions, the standard observational procedure was used which consists of interleaving the science target with reference stars. These stars are listed in Table 1 and taken from SearchCal.

The data were reduced with the `pndrs` package described in Le Bouquin et al. (2011, version v2.3 for 2011 data, v2.591 for 2013 and v2.71 for 2014). The main procedure is to compute squared visibilities and triple products for each baseline and spectral channel, and to correct for photon and readout noises. The combination of four telescopes provides six visibility and four closure phase measurements per spectral channel. Although the PIONIER observations were performed mostly with small telescope configurations, the baselines range from 10 m to 128 m, allowing us to probe the spatial scale 0.02–2 au (2–200 mas; where the interferometric FoV for our data set defined from the optical path difference separation of the fringe packet of two components). Data are displayed in Fig. B.2.

## 2.4. VLTI/GRAVITY observations

AU Mic was observed with GRAVITY and the Auxiliary Telescopes (ATs) on 2016 June 16 during the ESO science verification run of the instrument. GRAVITY is a near-IR interferometric beam combiner operating in the K-band ($\lambda$ = 1.990–2.450 µm). The observations were carried out using the highest spectral resolution (R ~ 4000) and the telescope baselines A0-G1-J2-K0 (projected baselines range from 48 m to 110 m). Three datasets of 300 s each (10 × 30 s) interleaved by sky exposures were recorded. To monitor the interferometric transfer function, the calibrator HD 155276 (K0.5IIICN0.5, K = 3.93 mag, $\theta_{LD}$ = 0.799 ± 0.078 mas) was observed (taken from SearchCal). However, HD 155276 was observed about 6 h before AU Mic, and therefore the calibration of the squared visibilities is likely unreliable. However, the closure phases can be used as they are independent of telescope-specific phase shifts induced by the atmosphere or optics (see e.g. Monnier 2007). Observables (visibilities and closure phases) were extracted using the provided GRAVITY pipeline (v1.1.2, Lapeyrere et al. 2014). The projected baselines allow us to probe the inner region of 0.03–2.4 au (3.5–250 mas) from the central star (FoV limited by the fiber FoV and not the bandwidth smearing thanks to the high spectral resolution). Data are displayed in Fig. B.2.

# 3. Companion search and detection limits

## 3.1. The CANDID algorithm

To search for a possible companion and set detection limits, we used the interferometric tool CANDID<fnref n="4" /> (Gallenne et al. 2015). The main function allows a systematic search for companions (point-sources in this case) performing an $N \times N$ grid of fits, the minimum required grid resolution of which is estimated a posteriori in order to find the global minimum in $\chi^2$. The tool delivers the binary parameters, namely the flux ratio $f$, and the relative astrometric separation ($\Delta\alpha, \Delta\delta$). The uniform disc (UD) angular diameters of both components can also be fitted, but in our case we do not expect a spatially resolved companion. Furthermore, the angular diameter $\theta_{UD}$ of AU Mic is only resolved by GRAVITY and PIONIER. The significance of the detection is also given, taking into account the reduced $\chi^2$ and the number of degrees of freedom.

The second main function of CANDID incorporates a robust method to set a $n\sigma$ detection limit on the flux ratio for undetected components, which is based on an analytical injection of a fake companion at each point in the grid. We refer the reader to Gallenne et al. (2015) for more details about CANDID and its formalism. Two-dimensional detection-limit maps are expressed in contrast ratio in the corresponding bandpass and spatial range of the instrument. To have a quantitative estimate of the sensitivity limit with respect to the separation, we estimated a radial profile using the 90% completeness level (i.e. 90% of all possible positions) from the cumulated histogram in rings for all azimuths.

## 3.2. Angular diameter of AU Mic

To measure the diameter, we only used the #1 PIONIER observations because this is the only dataset with long enough baselines to allow a reliable measurement. In addition, PIONIER operates in the H band and therefore provides a better spatial resolution

<footnote n="2">² http://www.jmmc.fr/english/tools/data-bases/oidb/</footnote>

<footnote n="3">³ https://www.jmmc.fr</footnote>

<footnote n="4">⁴ Available at https://github.com/amerand/CANDID and https://github.com/agallenne/GUIcandid for a GUI version.</footnote>
than GRAVITY. There is no bright companion orbiting AU Mic reported, and therefore we do not expect the angular diameter measurement to be biased by the presence of a low-contrast component.

We first estimated the uniform disc diameter $\theta_{UD}$ by fitting the calibrated squared visibilities with the following formula:

$$|V(u, v)|^2 = \left| \frac{2 J_1(x)}{x} \right|^2 \qquad (1)$$

where $J_1(x)$ is the first-order Bessel function, $x = \pi\theta_{UD}B/\lambda$, $\lambda$ the wavelength, $B = \sqrt{u^2 + v^2}$ the baseline length and $(u, v)$ are the spatial frequencies. We measured $\theta_{UD} = 0.797 \pm 0.032_{\text{stat}} \pm 0.038_{\text{sys}}$ mas. The statistical error was determined using the boot-strapping technique (with replacement) on all baselines. The given diameter corresponds to the median of the distribution, while the maximum value between the 16th and 84th percentile was chosen as the statistical uncertainty. The systematic error is the quadratic sum of two systematic uncertainties. The first comes from the calibrator uncertainties; we took the quadratic average of the calibrator uncertainties used. The second is a 0.35% systematic error due to the uncertainty of the wavelength calibration of PIONIER (Gallenne et al. 2012). We used this value in the following CANDID analysis. Our measurement is in good agreement with the $0.78 \pm 0.023$ mas reported by White et al. (2015, we assume they report the uniform disc diameter as there is no reference to any limb-darkening calculation).

We additionally tested a UD model including an IR emission from the hot dust in the inner disc. Such emission is usually difficult to detect because the flux ratio $f_d = f_{\text{disc}}/f_\star$ is of the order of 1% which needs high-precision visibilities and small baselines. We fitted a simplified visibility model $(1-2f_d) * |V(u, v)|^2$ which assumes the disc is fully resolved resulting in a visibility deficit compared to the values expected from the star alone. We found $\theta_{UD} = 0.679 \pm 0.12$ mas and $f_d = 1.09 \pm 0.67\%$. The angular diameter is smaller but within $1\sigma$ with the single star value. In addition, $f_d$ is barely constrained with a $\chi^2$ difference of 0.2 only, corresponding to a significance level of $1.8\sigma$. We can therefore assume that the single star model provides an unbiased angular diameter.

We determine the limb-darkened angular diameters, $\theta_{LD}$. Assuming a circular symmetry, we fitted the calibrated squared visibilities following the formalism of Mérand et al. (2015), which consists in extracting the radial intensity profile $I(r)$ of the spherical SATLAS models (Neilson & Lester 2013), which was converted to a visibility profile using a Hankel transform:

$$V_{\lambda}(x) = \frac{\int_{0}^{1} I_{\lambda}(r) J_{0}(rx) r dr}{\int_{0}^{1} I_{\lambda}(r) r dr}, \qquad (2)$$

where $\lambda$ is the wavelength, $x = \pi\theta_{LD}B/\lambda$, $B$ is the interferometric baseline projected onto the sky, $J_0$ the Bessel function of the first kind, and $r = \sqrt{1 - \mu^2}$, with $\mu = \cos(\theta)$, $\theta$ being the angle between the line of sight and a surface element of the star.

The SATLAS grid models span effective temperatures from 3000 to 8000 K in steps of 100 K, effective gravities from -1 to 3 in steps of 0.25, and masses from 0.5 to 20 $M_\odot$. We chose the model with the closest parameters, that is with $T_{\text{eff}} = 3700$ K, $\log g = 4.5$, and $M = 0.5 M_\odot$ (Plavchan et al. 2020, with $\log g \sim 4.4$ dex estimated from the mass, temperature and luminosity). We estimated $\theta_{LD} = 0.825 \pm 0.033_{\text{stat}} \pm 0.038_{\text{sys}}$ mas and errors were estimated in the same way as above. Changing the models with $T_{\text{eff}} \pm 200$ K, $\log g \pm 0.5$, and $M \pm 2.5 M_\odot$ changes the

<figure><img src="image_3.png" /><figcaption>Fig. 2. Best-fit solution of the limb-darkened disc model (red dots) to the PIONIER squared visibilities of AU Mic (blue dots).</figcaption></figure>

diameters by at most 0.1%, which we also quadratically added. Combining our estimate with the Gaia parallax of 102.943 $\pm 0.023$ mas provides a linear radius of $R = 0.862 \pm 0.052 R_\odot$. The best-fit solution for the angular diameter is shown in Fig. 2.

We also checked the angular diameter using the squared visibilities from GRAVITY. We found $\theta_{LD} = 1.078 \pm 0.255_{\text{stat}} \pm 0.078_{\text{sys}}$ mas (systematic error from the uncertainty of the calibrator). This value is higher than the one determined with PIONIER, although within $1\sigma$. The squared visibilities are not very accurate or precise, with a standard deviation of $\sim 10\%$ (see Fig. B.2). This is likely linked to a poor calibration, which was somewhat expected as the calibrator was observed 6 h before. However, as stated previously, the closure phases are independent of atmosphere seeing, and therefore we use only this observable for the following analyses.

### 3.3. Search for orbiting companion

Every instrument has a different minimum and maximum spatial range $r$ within which we can search for a companion. The loss of coherence caused by spectral smearing mainly limits the maximum separation for a reliable detection ($\le R\lambda/B$), while the smallest spatial scale is set by the angular resolution ($\sim \lambda/2B$). The ranges are mentioned in Sect. 2; we determined $r \sim 60-550$ mas for SPHERE, $r \sim 70-730$ mas for NACO, $r \sim 2-200$ mas for PIONIER and $r \sim 3.5-250$ mas for GRAVITY. In the following, we searched for companions using CPs only. As explained by Gallenne et al. (2015), the CPs are more sensitive to faint off-axis companions (although depending on its location and the $(u, v)$ coverage), and are also less affected by instrumental and atmospheric perturbations. The angular diameter is fixed in the fitting process to the uniform disc value listed above.

In the following, the contrasts have been converted to masses using the AMES-COND models (Baraffe et al. 2003).

**SPHERE data.** We did not detect any companion at more than $1.9\sigma$. Figure 3 shows the 2D detection maps given by
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** χ² map of the local minima for the SPHERE observations (**left panel**), detection level map (**middle panel**), and the (u, v) plane coverage of the observations (**right panel**). The yellow lines represent the convergence from the starting points to the final fitted position. The maps were re-interpolated in a regular grid for clarity. The yellow star denotes AU Mic, while the orange arrow in the middle indicates the disc orientation.</figcaption></figure>

**CANDID.** The highest detection is at ρ ~ 72 mas and PA ~ -17° with a flux ratio *f*<sub>K</sub> = 0.041%, which would correspond to a companion mass of ~7.3 *M*<sub>j</sub>. However, even if the position seems slightly above the disc plane (PA = 128.41°, MacGregor et al. 2013), this detection is not significant. The (u, v) coverage, that is, the sampling of the Fourier plane, is displayed in Fig. 3.

**NACO data.** A clear location at ρ ~ 139 mas and PA ~ -156° is detected in the χ² map but the detection level is only 1.6σ (see maps in Appendix C). The flux ratio is *f*<sub>L</sub> = 0.07% (~4.2 *M*<sub>j</sub>) but its location perpendicular to the disc plane makes this detection unlikely, although similar in brightness with the previous highest detection of SPHERE.

**PIONIER data.** We analysed all dates independently. There is also no detection with a statistical significance larger than 2.9σ (see maps in Appendix C). The highest peak has a flux ratio of *f*<sub>H</sub> = 0.49% (~14.1 *M*<sub>j</sub>) and is located at ρ ~ 44 mas and PA ~ -109°. We can see additional peaks but these are spuriously produced by an incomplete (u, v) coverage, and their significance is less than 2.9σ.

**GRAVITY data.** Only *CPs* from the science combiner have been used as it offers a much greater number of spectral channels than the fringe tracker (five channels only). No additional sources are detected at more than 2.6σ. The highest peak is at ρ ~ 14.8 mas, PA ~ -128°, and with a flux ratio of *f*<sub>K</sub> = 0.42% (~13.5 *M*<sub>j</sub>). This companion is unlikely as it is too bright, and the flux ratio is not consistent with the SPHERE estimate.

We therefore did not significantly detect (i.e. > 5σ) any component orbiting AU Mic. The highest probable locations given by the PIONIER and GRAVITY data are unlikely as they would correspond to a brown dwarf. Although such objects would not have been detected from AO imaging because of the limited angular resolution and size of the coronagraph, it would have been detected from the extensive observing RV observations (from Kepler's law and assuming our measured projected separation as a lower limit for the angular semi-major axis, the minimum RV amplitude would be ~0.8 km s⁻¹).

### 3.4. 5σ detection limits

CANDID has also implemented a robust method to derive the dynamic range we can reach with a given set of data. It consists in injecting a fake companion into the data at each astrometric position with different flux ratios. As we inject a companion, we therefore know that the binary model should be the true model.

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Contrast limit at 5σ for our SPHERE observations.</figcaption></figure>

We then compare the χ² with that of a single star model (uniform disc model) to obtain the probability of the binary model being the true model. We set the significance level on the flux ratios at 5σ, meaning that lower flux ratios are not significantly detected. Doing this for all points in the grid, we then have a 5σ detection limit map for the flux ratio. We refer the reader to Gallenne et al. (2015) for more information about the method. To obtain a quantitative estimate of the sensitivity limit with respect to the separation, we estimated a radial profile using the 90% completeness level (i.e. 90% of all possible positions) from the cumulated histogram in rings for all azimuths. The radial profile for the detection limit of our SPHERE observations is displayed in Fig. 4.

We list three different values in Table 2: the average limit for the separation ranges *r* < *r*<sub>mid</sub>, *r*<sub>mid</sub> *r* < *r* < *r*<sub>max</sub>, and the full range *r* < *r*<sub>max</sub>. All of the final 5σ contrast limits, Δ*m*<sub>5σ</sub> expressed in magnitude, are conservative as they correspond to the mean plus the standard deviation for the given radius range. Upper limit masses, that is, the lowest masses we can safely exclude, were estimated from the contrast limits calculated for the full separation range.

## 4. Discussion and conclusion

With our new SPHERE observations, we do not detect additional objects brighter than *K*<sub>s</sub> ~ 11.2 mag within 0.4–2.4 au, which
<table><caption>**Table 2.** Average $5\sigma$ contrast limits ($\Delta m_1$).</caption><thead><tr><th></th><th>Date</th><th colspan="3">Δm (mag)</th><th>Upper limit mass<br>M<sub>jup</sub></th></tr></thead><tbody><tr><td>SPHERE (K<sub>s</sub>)</td><td>2018-08</td><td>0.4 &lt; r &lt; 1 au<br>6.67</td><td>1 &lt; r &lt; 2.4 au<br>6.72</td><td>0.4 &lt; r &lt; 2.4 au<br>6.68</td><td>12</td></tr><tr><td>NACO (L')</td><td>2011-08</td><td>0.7 &lt; r &lt; 3.9 au<br>6.51</td><td>3.9 &lt; r &lt; 7 au<br>6.32</td><td>0.7 &lt; r &lt; 7 au<br>6.33</td><td>9</td></tr><tr><td>PIONIER (H)</td><td>2011-11</td><td>0.02 &lt; r &lt; 0.25 au<br>4.87</td><td>0.25 &lt; r &lt; 0.5 au<br>5.10</td><td>0.02 &lt; r &lt; 0.5 au<br>5.02</td><td>17</td></tr><tr><td>GRAVITY (K)</td><td>2016-06</td><td>0.03 &lt; r &lt; 1.4 au<br>5.25</td><td>1.2 &lt; r &lt; 2.4 au<br>5.52</td><td>0.03 &lt; r &lt; 2.4 au<br>5.39</td><td>16</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Detection probability maps of AU Mic assuming coplanarity with the disc (left panel) or with no hypothesis on the inclination (right panel), combining RV, interferometric, and direct-imaging observations. Blue colour-coded areas denote the detectability probability. The red area shows where the newly detected companions should be located (Martioli et al. 2021).</figcaption></figure>

would correspond to a mass limit for a planet of $12.3 \pm 0.5 M_{\text{jup}}$. The mass limit was estimated using a cubic spline interpolation of the AMES-COND isochrones encompassing 18.5 Ma (i.e. 9, 10, 20, and 30 Ma). We then interpolated the mass for the corresponding magnitude limit and age of $18.5 \pm 2$ Ma. The best mass sensitivity is achieved with NACO because the $L'$ band is favourable at young ages, providing an upper limit on the mass of $\sim 8.3 \pm 0.8 M_{\text{jup}}$ within 0.7–7.3 au. This is similar to the upper mass limit of $10.1 M_{\text{jup}}$ determined by Gauchet et al. (2016) with the same data set, but updated with our more robust CANDID tool. The SAM observations enable us to reach higher contrast thanks to a better AO correction and a more stable point spread function. In addition, the light beam in SPHERE is prone to a higher throughput than the whole VLTI optical path. Long-baseline interferometric observations enable us to probe smaller spatial scales but the contrast limit reached is not as good. Within 0.02–0.5 au, PIONIER provides a detectable mass limit of $\sim 17.1 \pm 1.0 M_{\text{jup}}$.

We added additional constraints on the mass limits by combining our calculated $5\sigma$ detection maps with the existing RV data from Lannier et al. (2017) and SPHERE contrast maps estimated from the data of Boccaletti et al. (2018, determined using a PCA analysis of the IFS and IRDIS images taken in the J, H, H23, and K12 bands). We used the Multi-epoch

multi-purpose Exoplanet Simulation System algorithm (MESS2, Lannier et al. 2017) which generates populations of synthetic planets through a Monte Carlo simulation for a given range of mass and separation. MESS2 grids are then compared to our combined data sets to explore the detection probability. Our contrast limit maps were converted into masses using an age of 18.5 Ma (Miret-Roig et al. 2020) and the COND-2003 evolutionary models (Baraffe et al. 2003). We considered two inclination distributions for the synthetic planets: a uniform distribution between 0° and 90°, and an inclination of $90 \pm 1°$, that is, assuming the synthetic planets orbit in the same plane as the resolved disc. We also took an ad hoc eccentricity distribution with $e < 0.6$. The probability maps are displayed in Fig. 5. As expected, assuming a uniform prior on the inclination significantly reduces the sensitivity to companions. This is particularly true within 1 au where the sensitivity is driven by the RV technique, which is more sensitive to edge-on systems. In the case of co-planar orbits, 99% of additional companions with $\ge 4 M_{\text{jup}}$ can be excluded within 0.2–1 au, and $\ge 2 M_{\text{jup}}$ within 1–8 au. These detection limits are not better than the ones derived by Lannier et al. (2017) combining RVs and direct imaging at separation $r > 0.2$ au. This is because direct imaging has a better sensitivity than interferometry at these angular separations. However, for the innermost region, at $\le 0.2$ au, which is the
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Same as Fig. 5 but only interferometric and direct-imaging observations.</figcaption></figure>

spatial scale probed by interferometry, we find that 99% of the companions down to ~0.5 $M_{\text{jup}}$ can be detected. This is actually similar to using RVs alone because of the lack of high-dynamic range of the current interferometric instruments and our limited datasets. In this specific case (i.e. for co-planar orbits of the edge-on disc), interferometry does not bring additional constraints when precise RVs are available. In a more general case where no precise RVs are available, interferometry is complementary to direct imaging by providing better contrast in the innermost region, as we can see in Fig. 6. We also notice that the coplanarity hypothesis mostly affects the RVs, and has a limited effect when only interferometry and direct imaging are used.

The newly detected planets have a mass of < 0.08 $M_{\text{jup}}$ for a maximum separation of 11.3 mas (0.11 au), well below our contrast limits (red area in Fig. 5). This would correspond to a lower limit K-band magnitude of ~33 mag, corresponding to a contrast of ~4 × 10⁻¹² ($\Delta K$ ~ 28.5 mag, extrapolated from the AMES-COND models). Such dynamic range is beyond the capabilities of the current imaging and interferometric instruments, and also beyond future facilities in the near future. Only space-based coronagraph like Habex or Louvoir (see e.g. Mennesson et al. 2016; Bolcar et al. 2017) could reach such contrasts but they would not have the required angular resolution. Direct imaging of these planets with AO is not possible due to its proximity with the host star. Only long-baseline interferometry can probe the innermost regions, but the current reachable contrast is about $\Delta H$ ~ 6.5 mag (Gallenne et al. 2015; Roettenbacher et al. 2015). The dual-field mode of GRAVITY reached a contrast of $\Delta K$ = 10.7 mag by directly measuring a spectrum of an exoplanet (GRAVITY Collaboration 2019), but this is a special case as the location of the companion must be known a priory to be able to position the fibre on it. A 'blind search' was tested on AU Mic with the Unit Telescopes by dithering the fibre along the southeastern part of the debris disc (GRAVITY Collaboration 2019). No planet was detected but they achieved a 5σ dynamic range of 11 mag at 1.2 au and 13.5 mag at 2.4 au. The astrometric mode of GRAVITY, which can provide a differential astrometric precision of a few tens of µas, could also be an alternative used to obtain the astrometric orbit of the more massive planet; however the astrometric wobble produced by this planet on the central star is too small (~1.7 µas, taking the upper limit of 0.08 $M_{\text{jup}}$ at 0.11 au).

The upcoming ELT with a SAM mode will likely also not detect these planets. Assuming a maximum baseline of 38 m (for the possibly two most distant mask holes) and observations performed in the K band, planets detectable down to ~7 mas would need an optimal AO correction and strongly improved sensitivity limits (although near the SAM resolution limit). Unfortunately, even the next generation of extreme-AO instruments for the ELT will not provide a contrast of better than 10⁻⁷ in such close proximity to the star.

**Acknowledgements.** A.G. acknowledges the support of the French Agence Nationale de la Recherche (ANR), under grant ANR-15-CE31-0012-01 (project UnlockCepheids). A.G. acknowledges financial support from "Programme National de Physique Stellaire" (PNPS) of CNRS/INSU, France. A.G. also acknowledges financial support from the ANID-ALMA fund No. ASTRO20-0059. The research leading to these results has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme under grant agreement No 695099 (project CepBin), grant agreement no. 951549 (project UniverScale) and from the National Science Center, Poland grants MAESTRO UMO-2017/26/A/ST9/00446 and BEETHOVEN UMO-2018/31/G/ST9/03050. A.G. acknowledges support from the IdP II 2015 0002 64 and DIR/WK/2018/09 grants of the Polish Ministry of Science and Higher Education, from the BASAL Centro de Astrofísica y Tecnologías Afines (CATA) AFB-170002, and from the Millennium Institute of Astrophysics (MAS) of the Iniciativa Científica Milenio del Ministerio de Economía, Fomento y Turismo de Chile, project IC120009. S.K. acknowledges support from an ERC Consolidator Grant (Grant Agreement ID 101003096). The GRAVITY data have been recorded as part of the Science Verification program and we thank the SV team, which is composed of ESO employees and GRAVITY consortium members, for carrying out the observations (https://www.eso.org/sci/activities/vltsv/gravitysv.html). This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

## References

Augereau, J. C., & Beust, H. 2006, A&A, 455, 987

Baraffe, I., Chabrier, G., Barman, T. S., Allard, F., & Hauschildt, P. H. 2003, A&A, 402, 701

Beuzit, J. L., Vigan, A., Mouillet, D., et al. 2019, A&A, 631, A155

Boccaletti, A., Thalmann, C., Lagrange, A.-M., et al. 2015, Nature, 526, 230

Boccaletti, A., Sezestre, E., Lagrange, A.-M., et al. 2018, A&A, 614, A52

Bolcar, M. R., Aloezos, S., Bly, V. T., et al. 2017, SPIE Conf. Ser., 10398, 1039809
Cheetham, A. C., Girard, J., Lacour, S., et al. 2016, SPIE Conf. Ser., 9907, 99072T

Chiang, E., & Fung, J. 2017, ApJ, 848, 4

Daley, C., Hughes, A. M., Carter, E. S., et al. 2019, ApJ, 875, 87

Dohlen, K., Langlois, M., Saisse, M., et al. 2008, SPIE Conf. Ser., 7014, 70143L

Eisenhauer, F., Perrin, G., Brandner, W., et al. 2011, The Messenger, 143, 16

Esposito, T. M., Kalas, P., Fitzgerald, M. P., et al. 2020, AJ, 160, 24

Fitzgerald, M. P., Kalas, P. G., Duchêne, G., Pinte, C., & Graham, J. R. 2007, ApJ, 670, 536

Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1

Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1

Gallenne, A., Kervella, P., Mérand, A., et al. 2012, A&A, 541, A87

Gallenne, A., Mérand, A., Kervella, P., et al. 2015, A&A, 579, A68

Gauchet, L., Lacour, S., Lagrange, A. M., et al. 2016, A&A, 595, A31

Graham, J. R., Kalas, P. G., & Matthews, B. C. 2007, ApJ, 654, 595

GRAVITY Collaboration (Abuter, R., et al.) 2019, The Messenger, 178, 47

GRAVITY Collaboration (Lacour, S., et al.) 2019, A&A, 623, A11

Kalas, P., Liu, M. C., & Matthews, B. C. 2004, Science, 303, 1990

Klein, B., Donati, J.-F., Moutou, C., et al. 2021, MNRAS, 502, 188

Krist, J. E., Ardila, D. R., Golimowski, D. A., et al. 2005, AJ, 129, 1008

Lannier, J., Lagrange, A. M., Bonavita, M., et al. 2017, A&A, 603, A54

Lapeyrere, V., Kervella, P., Lacour, S., et al. 2014, SPIE Conf. Ser., 9146, 91462D

Launhardt, R., Henning, T., Quirrenbach, A., et al. 2020, A&A, 635, A162

Le Bouquin, J.-B., Berger, J.-P., Lazareff, B., et al. 2011, A&A, 535, A67

Lenzen, R., Hartung, M., Brandner, W., et al. 2003, SPIE Conf. Ser., 4841, 944

Lindegren, L., Klioner, S. A., Hernández, J., et al. 2021, A&A, 649, A2

Liu, M. C. 2004, Science, 305, 1442

Lomax, J. R., Wisniewski, J. P., Roberge, A., et al. 2018, AJ, 155, 62

MacGregor, M. A., Wilner, D. J., Rosenfeld, K. A., et al. 2013, ApJ, 762, L21

Martioli, E., Hébrard, G., Correia, A. C. M., Laskar, J., & Lecavelier des Etangs, A. 2021, A&A, 649, A177

Matthews, B. C., Kennedy, G., Sibthorpe, B., et al. 2015, ApJ, 811, 100

Mennesson, B., Gaudí, S., Seager, S., et al. 2016, SPIE Conf. Ser., 9904, 99040L

Mérand, A., Kervella, P., Breitfelder, J., et al. 2015, A&A, 584, A80

Miret-Roig, N., Galli, P. A. B., Brandner, W., et al. 2020, A&A, 642, A179

Monnier, J. D. 2007, New Astron. Rev., 51, 604

Neilson, H. R., & Lester, J. B. 2013, A&A, 554, A98

Plavchan, P., Barclay, T., Gagné, J., et al. 2020, Nature, 582, 497

Quillen, A. C., Morbidelli, A., & Moore, A. 2007, MNRAS, 380, 1642

Roettenbacher, R. M., Monnier, J. D., Fekel, F. C., et al. 2015, ApJ, 809, 159

Rousset, G., Lacombe, F., Puget, P., et al. 2003, SPIE Conf. Ser., 4839, 140

Schneider, P. C., & Schmitt, J. H. M. M. 2010, A&A, 516, A8

Schneider, G., Grady, C. A., Hines, D. C., et al. 2014, AJ, 148, 59

Sezestre, É., Augereau, J. C., Boccaletti, A., & Thébault, P. 2017, A&A, 607, A65

Strubbe, L. E., & Chiang, E. I. 2006, ApJ, 648, 652

Tuthill, P. G., Monnier, J. D., Danchi, W. C., Wishnow, E. H., & Haniff, C. A. 2000, PASP, 112, 555

Wang, J. J., Graham, J. R., Pueyo, L., et al. 2015, ApJ, 811, L19

White, R. J., Schaefer, G., Ten Brummelaar, T., et al. 2015, BAAS, 225, 348.12

Wilner, D. J., Andrews, S. M., MacGregor, M. A., & Hughes, A. M. 2012, ApJ, 749, L27
Appendix A: Transfer functions of our
SPHERE/SAM observations

<figure><img src="image_8.png" /><figcaption>Fig. A.1. Transfer function calibration for the squared visibilities. Calibrator and science target observations are represented with black and blue dots, respectively. The solid line denotes the interpolation of the transfer function, while the grey area represents the standard deviation with the calibrators.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Fig. A.2. Same as Fig. A.1 but for the closure phases.</figcaption></figure>
Appendix B: Squared visibilities and closure
phases for all observations

<figure><img src="image_10.png" /><figcaption>Fig. B.1. Squared visibilities and closure phases for the SPHERE (left) and NACO (right) observations (blue). The red dots correspond to the best model found with CANDID.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. B.2.** Squared visibilities and closure phases for the PIONIER and GRAVITY observations (blue). The red dots correspond to the best model found with CANDID.</figcaption></figure>
Appendix C: Detection maps for the NACO,
PIONIER and GRAVITY observations

<figure><img src="image_12.png" /><figcaption>**Fig. C.1.** χ² map of the local minima for the NACO observations (left panel), detection level map (middle panel), and the (u, v) plane coverage of the observations (right panel). The yellow lines represent the convergence from the starting points to the final fitted position. The maps were re-interpolated in a regular grid for clarity. The yellow star denotes AU Mic, while the orange arrow in the middle indicates the disc orientation.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>Fig. C.2. Same as Fig. C.1 but for the PIONIER observations of 2011 and 2013.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>Fig. C.3. Same as Fig. C.1 but for the 2014 PIONIER (top) and GRAVITY (bottom) observations.</figcaption></figure>