# SPHERE adaptive optics performance for faint targets

M I Jones, J. Milli, I. Blanchard, Z. Wahhaj, R J de Rosa, C. Romero, N.

Ihanec

► To cite this version:

M I Jones, J. Milli, I. Blanchard, Z. Wahhaj, R J de Rosa, et al.. SPHERE adaptive optics performance for faint targets. Astronomy & Astrophysics - A&A, 2022, 667, pp.A114. 10.1051/0004-6361/202243820. [hal-04975562](https://hal.science/hal-04975562v1)

HAL Id: hal-04975562

https://hal.science/hal-04975562v1

Submitted on 4 Mar 2025

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>
# SPHERE adaptive optics performance for faint targets

M. I. Jones<fnref n="1" />, J. Milli<fnref n="2" />, I. Blanchard<fnref n="1" />, Z. Wahhaj<fnref n="1" />, R. J. De Rosa<fnref n="1" />, C. Romero<fnref n="1" />, and N. Ihanec<fnref n="1" />,<fnref n="3" />

Received 20 April 2022 / Accepted 27 August 2022

## ABSTRACT

*Context.* High-contrast imaging is a powerful technique for detecting and characterizing planetary companions at orbital separations $\ge 100$ mas from their parent stars.

*Aims.* We aim to study the limiting magnitude of the VLT/SPHERE adaptive optics (AO) system and the corresponding instrument performance for faint targets ($G \ge 11.0$ mag).

*Methods.* We computed the coronagraphic H-band raw contrast and the full width at half maximum (FWHM) of the non-coronagraphic point spread function (PSF), for a total of 111 different stars observed between 2016 and 2022 with IRDIS. For this, we processed a large number of individual frames that were obtained under different atmospheric conditions. We then compared the resulting raw contrast and the PSF shape as a function of the visible wave front sensor (WFS) instant flux, which scales with the G-band stellar magnitude. We repeated this analysis for the top 10% (TCAT10) and top 30% (TCAT30) best turbulence conditions on Cerro Paranal.

*Results.* We found a strong decrease in the coronagraphic contrast for stars fainter than $G \sim 12.5$ mag, even under the best atmospheric conditions. In this regime, the AO correction is dominated by the read-out noise of the WFS detector. In particular we found roughly a factor of 10 decrease in the raw contrast ratio between stars with $G \sim 12.5$ and $G \sim 14.0$ mag. Similarly, we observed a sharp increase in the FWHM of the non-coronagraphic PSF beyond $G \sim 12.5$ mag, and a corresponding decrease in the strehl ratio from $\sim 0.5$ to $\sim 0.2$ for the faintest stars. The decrease in the contrast ratio and PSF sharpness is slightly more pronounced for TCAT30 than for TCAT10.

**Key words.** instrumentation: adaptive optics – atmospheric effects

## 1. Introduction

High-contrast direct imaging (DI) is a powerful technique for detecting and characterizing young and self-luminous giant planets, reflected light planets, highly polarized disks, Solar System objects, among other interesting objects. Thanks to its extreme adaptive optics (AO) system, called SAXO (Fusco et al. 2006), the Spectro-Polarimetric High-contrast Exoplanet REsearch instrument (SPHERE; Beuzit et al. 2019) is one of the most powerful instruments in the world for direct imaging science, and after more than 6 yr of operations, it has led to the detection and characterization of a variety of (forming) planetary systems and protoplanet disks. Notable examples of these include $\beta$ Pictoris b (Lagrange et al. 2018), the PDS 70 planet-forming system (Müller et al. 2018; Keppler et al. 2018), the multiplanet systems around HD 8799 (Zurlo et al. 2016) and YSES 1 (Bohn et al. 2020), and the TW Hydrae disk (van Boekel et al. 2017), among others.

While different long-term SPHERE DI programs have mainly studied relatively bright stars (e.g., SHARDDS: Milli et al. 2017a; SHINE: Desidera et al. 2021; YSES: Bohn et al. 2021; BEAST: Janson et al. 2021), there is still a large population of faint stars located in young (age $\le 10–20$ Myr) stellar associations, which are ideal laboratories for studying in-situ planet formation inside protoplanetary disks. However, these stars are embedded in gas- and dust-rich disks, which absorb a large fraction of the stellar visible light, making these objects very red (see a discussion in Boccaletti et al. 2020). This effect is even stronger for M-type stars, which are intrinsically fainter in the

optical when compared to the near infrared (NIR). The main disadvantage of targeting these kinds of stars with SPHERE is that the instant flux recorded in the visible wave front sensor (WFS) is strongly reduced, and thus the AO correction is highly degraded. As increasing number of such faint stars (in visible light) are routinely being observed with SPHERE, we are able to better characterize the SAXO limiting magnitude. Motivated by this idea, we studied the SPHERE performance for targets fainter than $G = 11.0$ mag, where $G$ corresponds to the Gaia Early Data Release 3 (Gaia EDR3) magnitude (Riello et al. 2021).

This paper is organized as follows: in Sect. 2 we present the data used for the analysis and we describe our methods in detail. In Sect. 3, we discuss about the instant flux received by the WFS. In Sect. 4, we present the results of the raw contrast (from coronagraphic images) and the full width at half maximum (FWHM) of the non-coronagraphic point spread function (PSF), as a function of the WFS flux and the stellar G-band magnitude. Similarly, we computed the strehl ratio (SR) for a handful of selected stars and we compare the resulting values with the FWHM previously measured. In Sect. 5, we compare the contrast from the raw images to that obtained in post-processing data. Finally, the summary, discussion and conclusions are presented in Sect. 6.

## 2. Data and method

We analyzed a total of 111 different stars observed with SPHERE between 2016 and 2022, all of them fainter than $G = 11.0$ mag, which are listed in Table A.1. Figure 1 shows the Gaia EDR3 color-magnitude diagram (CMD) with the positions of all of

<footnote n="1">¹ European Southern Observatory, Alonso de Córdova 3107, Vitacura, Casilla 19001, Santiago, Chile
e-mail: mjones@eso.org</footnote>

<footnote n="2">² Université Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France</footnote>

<footnote n="3">³ Astronomical Observatory, University of Warsaw, Al. Ujazdowskie 4, 00-478 Warsaw, Poland</footnote>
<table><caption>**Table 1.** Header keywords used to filter the data for the flux and coronagraphic observations.</caption><thead><tr><th>Keyword</th><th>Flux sequence</th><th>Coronagraphic sequence</th></tr></thead><tbody><tr><td>DPR TYPE</td><td>OBJECT, FLUX</td><td>OBJECT</td></tr><tr><td>DET SEQ1 DIT</td><td>≥0.87 s</td><td>≥2 s</td></tr><tr><td>DPR TECH</td><td>IMAGE/POLARIMETRY</td><td>IMAGE/POLARIMETRY</td></tr><tr><td>INS1 FILT NAME</td><td>B_H</td><td>B_H</td></tr><tr><td>INS COMB IFLT</td><td>DB_H23 / DP_0_BB_H</td><td>DB_H23/DP_0_BB_H</td></tr><tr><td>AOS VISWFS MODE</td><td>Gain_1000_FREQUENCY_300 Hz</td><td>Gain_1000_FREQUENCY_300 Hz</td></tr><tr><td>INS4 FILT3 NAME</td><td>OPEN</td><td>OPEN</td></tr><tr><td>INS4 OPTI22 ID</td><td>SMALL/MEDIUM/LARGE</td><td>SMALL/MEDIUM/LARGE</td></tr></tbody></table>

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Gaia EDR3 color-magnitude diagram of all 111 stars included in the analysis. The typical uncertainties are smaller than the symbol size. The position of the reddest star in the sample (AX Cha) is also labeled. The corresponding histogram of the G-band magnitude distribution is overplotted on the right side.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Same as for Fig. 1, but for the J and H bands.</figcaption></figure>

the targets. Similarly, Fig. 2 shows a Two Micron All Sky Survey (2MASS; Cutri et al. 2003) NIR CMD for our stars. The median G-band and H-band magnitudes of the sample are 12.4 mag and 8.9 mag, respectively. We included in the analysis only H-band ($\lambda_c$ ~ 1.63 µm) observations taken with the visible (WFS) in the high-gain (M = 1000) low-frequency mode (300 Hz), with the WFS spectral filter in the OPEN position. We also included

both imaging and polarimetric observations. In the case of the coronagraphic observations, we included individual frames with detector integration time (DIT) longer than 2 s. We note that the vast majority of the faint star observations were made using these setups. The list of header keywords used to filter the data are summarized in Table 1. We note that all of the raw data (non-coronagraphic FLUX and coronagraphic OBJECT data cubes) and their corresponding AO files (FITS binary tables containing atmospheric and SAXO telemetry data; see Milli et al. 2017b) were downloaded from the ESO archive. For the analysis, we developed an IDL-based code, which computes the contrast for each individual raw frame, following the steps listed below.

First, we fitted a 2D Gaussian function to each individual non-coronagraphic FLUX frames. From the fit we estimated the peak flux ($f_{nc}$) and the corresponding FWHM in the x- and y-pixels direction. The final $f_{nc}$ value was computed from the mean of all individual flux frames. We then used the mean $f_{nc}$ to predict the peak flux ($f_c$) in each individual coronagraphic OBJECT frame, by scaling the DIT, neutral density filters and other optical elements with different transmissions (e.g. when a sequence is obtained with and without the apodizer). We then estimated the noise level ($\sigma_{300}$) at 24 pixels from the coronagraph center (corresponding to ~300 mas), by computing the standard deviation in a 5-pixels-wide annulus. The resulting $5\sigma$ raw contrast at 300 mas (RC$_{300}$) was obtained from:

$$RC_{300} = \frac{5 \times \sigma_{300}}{f_c} \qquad (1)$$

We note that for the broad-band observations (B_H filter) and the differential polarimetry (DP_0_BB_H filter), we used only the left side of the Infra-Red Dual-beam Imaging and Spectroscopy (IRDIS; Dohlen et al. 2008) detector, while for the dual-imaging mode (DB_H23 filter) we independently analyzed the left ($\lambda_c$ ~ 1.59 µm) and the right ($\lambda_c$ ~ 1.67 µm) sides of the detector, and we then computed the mean value from the two channels. Finally, to study the effect of the atmospheric conditions on the contrast and the PSF, we retrieved the Differential Image Motion Monitor (DIMM) seeing and coherence time ($\tau_0$) at the beginning of each exposure, directly from the image header.

### 3. Instant flux on the visible WFS

The SPHERE instrument is equipped with a 40 × 40 Shack-Hartmann WFS, (with a total of 1240 effective subapertures imaging the SPHERE entrance pupil). The light reaching the WFS is recorded by a 240 × 240 pixel electron multiplying charge-coupled device (EMCCD) detector, sensitive to visible
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** H-band 5σ raw contrast at 300 mas as a function of the visible WFS flux. The black and red lines correspond to atmospheric turbulence categories TCAT10 and TCAT30, respectively.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Mean FWHM as a function of the visible WFS flux. The black and red lines correspond to the atmospheric turbulence categories TCAT10 and TCAT30, respectively.</figcaption></figure>

light (~0.45–0.95 µm). The main advantage of this detector is that it can operate at very high frequencies (up to 1380 Hz), low read-out noise (<0.2 e-/pixel/frame), and high multiplication gains (up to M = 1000; Sauvage et al. 2016b). However, as mentioned before, for our analysis we included only data taken with the EMCCD operating at 300 Hz and M = 1000, which is the mode automatically selected by the real-time computing system (SPARTA; Fedrigo et al. 2006) for faint targets (R > 10.5 mag).

The instant flux in the WFS EMCCD detector ($f_{WFS}$) is recorded by SPARTA, and it is saved in the AO files. These files (among other telemetry data) contain the total flux per pupil (i.e., integrated in all 1240 subapertures). This flux (tagged FLUX_AVG in the binary tables) is thus expressed in units of [ADU/pupil/frame]. To convert FLUX_AVG into [e-1/subaperture/frame], we used:

$$f_{WFS} [e^{-1}/\text{subaperture/frame}] = \frac{g \times \text{FLUX_AVG}}{n_s \times M}, \quad (2)$$

where g = 16.5 [e-/ADU], n_s = 1240 [subapertures/pupil], and M = 1000 (unit-less).

# 4. Results

## 4.1. Raw contrast and FWHM as a function of the WFS flux

After computing the raw contrast, we cleaned up the sample by removing those frames that lead to very deviant values. Typical examples of these included cases with the stars out of the coronagraph, AO loop open, contamination from close binary companions, extreme cases of low-wind effect (Sauvage et al. 2016a; Milli et al. 2018), among others. We then compared the raw contrast as a function of $f_{WFS}$. Figure 3 shows the RC300 versus the $f_{WFS}$ binned in 1 [e-/subaperture/frame] flux bins. We note that only frames with header values listed in Table 1 were included (as explained in Sect. 2). The black dots represent the atmospheric turbulence category TCAT10, corresponding to a seeing <0.6 arcsec and a coherence time ($\tau_0$) >5.2 ms, and the red dots represent TCAT30 (corresponding to a seeing <0.8 arcsec and a $\tau_0$ > 4.1 ms). The error bars correspond to a 1σ equal-tailed confidence interval for each bin. In total, we included ~12 000 individual OBJECT frames in Fig. 3.

As can be seen, there is a sharp decrease in the achievable contrast for $f_{WFS}$ below ~3 [e-/subaperture/frame]. In this regime, the AO correction is degraded due to the lack of photons collected by the WFS. Above this value the RC300 is nearly constant. We also note that the RC300 achieved under TCAT10 conditions is slightly better than that under TCAT30. Similarly, Fig. 4 shows the mean FWHM (from the x- and y-direction) measured on non-coronagraphic FLUX frames as a function of $f_{WFS}$. Again, this quantity was computed on the left side of the IRDIS detectors in the dual-band mode (see Table 1), and includes a total of ~9000 individual frames. It can be seen that for $f_{WFS}$ level above ~3 [e-/subaperture/frame] the PSF is nearly flat, with a FWHM of ~53 mas, for both TCAT10 and TCAT30. On the other hand, for lower $f_{WFS}$ values, the AO correction degrades, and thus the FWHM increases with a decreasing flux level, reaching a mean value larger than ~100 mas, when reaching the subelectron per frame regime, for both TCAT10 and TCAT30.

## 4.2. Correlation between stellar magnitude and WFS flux

To understand the SPHERE performance not only as a function of $f_{WFS}$, but also as a function of the target magnitude, we investigated the correlation between $f_{WFS}$ and the stellar G-band magnitude. We chose this particular band because the transmission curve of the WFS detector is very close to the G-band one, and also since nearly all of the targets observed by SPHERE show a G-band entry in the EDR3 catalog. Figure 5 shows the normalized transmission of the Gaia EDR3 bands<fnref n="1" />, the Johnson R broad-band filter (Bessell 1990), and the WFS response, including all of the upstream optical elements and with the WFS spectral filter in the OPEN position (which is the case for faint stars). It can be seen that the G-band transmission function is the closest one to the WFS response function, although the latter one is slightly more sensitive to longer wavelengths. For this reason, to predict the expected $f_{WFS}$ value during an observing sequence of a faint target, the G-band should be used, instead of the R-band (as currently implemented in SPARTA). In fact, even though the value of these two magnitudes are generally very close, there are cases where the difference is as large as one magnitude (corresponding to a factor of ~2.5 in flux). For this analysis, we first corrected the flux in the WFS for the target airmass, using an

<footnote n="1">¹ https://www.cosmos.esa.int/web/gaia/edr3-passbands</footnote>
<figure><img src="image_6.png" /><figcaption>Fig. 5. Normalized transmission curves for the Gaia EDR3 G (solid black line), $G_B$ (dotted blue line), $G_R$ (dotted red line) and Johnson R (solid green line) broad-band filters. The total transmission curve (including upstream optics) of the WFS is overplotted (solid black dots).</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 6. Stellar G-band magnitude versus the logarithm of the WFS flux. The solid line corresponds to the best linear fit between these two quantities. The position of AX cha is also shown (red dot). The corresponding $f_{WFS}$ values expressed in linear scale and units of [$e^{-1}$/subaperture/frame] are labeled in the upper x-axis.</figcaption></figure>

extinction coefficient of 0.12 (Patat et al. 2011). We then compared this value with the stellar G-band magnitude. We note that to minimize the stellar color effect, we only included stars with $1 < G_B - G_R < 3$ (see Fig. 1). The results are shown in Fig. 6. The error bars in the WFS flux correspond to the standard deviation from all individual measurements. The solid line corresponds to the WFS magnitude scale ($m_{WFS}$), given by:

$$m_{WFS} = -2.5 \log(f_{WFS}) + 25.7. \quad (3)$$

We note that we removed outliers beyond $2.5 \sigma$ from the fit (18 stars in total), leading to a scatter of 0.18 mag. This value is mainly explained by observations performed under different observing conditions (e.g., sky transparency, seeing) and, more importantly, due to the fact that we included observations using different spatial filters (Fusco et al. 2016), as listed in Table 1, meaning that a different amount of the light was blocked for different observations. Ideally, this analysis should be repeated

<figure><img src="image_8.png" /><figcaption>Fig. 7. H-band $5\sigma$ raw contrast at 300 mas as a function of the stellar G-band magnitude. The black and red lines correspond to the atmospheric turbulence categories TCAT10 and TCAT30, respectively.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Fig. 8. Mean FWHM as a function of the stellar G-band magnitude. The solid black and red dots (lines) correspond to the atmospheric turbulence categories TCAT10 and TCAT30, respectively.</figcaption></figure>

with data taken under photometric conditions and using the WFS spatial filter in the LARGE position, as opposed to the majority of the observations analyzed here that were done using either the SMALL or the MEDIUM spatial filter. In addition, variable stars or intrinsically red stars, such as AX Cha, might introduce extra scatter to this correlation.

Using this correlation, we can directly compare the stellar G magnitude and the expected raw contrast, under different atmospheric conditions. We therefore replicated Figs. 3 and 4, but this time we replaced $f_{WFS}$ by the stellar G magnitude. The results are shown in Figs. 7 and 8. As can be seen, there is a decrease of roughly one order of magnitude in the contrast ratio beyond G ~ 12.5. For stars brighter than this, the contrast is nearly constant, meaning that it is limited by the atmospheric conditions and by the ability of the AO system to correct for the atmospheric turbulence. In this regime (G ~ 10.5–12 mag), the medium AO frequency mode (600 Hz) should be used, while currently the 300 Hz frequency mode is set automatically by the real-time computer. This is confirmed by recent SPHERE technical observations of stars with G between 11 and 12, which yielded a systematic improvement in SR and contrast, even in atmospheric conditions worse than TCAT30. Similarly, there is an increase
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Normalized diffraction-limited theoretical PSF of the VLT pupil at 1.63 µm, as seen by IRDIS/SPHERE.</figcaption></figure>

in the FWHM of the PSF from ~53 mas to ≥100 mas between G ~ 12.5 mag and ~14 mag.

### 4.3. Strehl ratio

In addition to the FWHM values, we investigated the dependence of the strehl ratio (SR) on $f_{\text{WFS}}$ and the stellar magnitude. For this, we first computed a theoretical diffraction-limited PSF for the VLT as seen by SPHERE, (i.e., including the effect of the Lyot stop and the apodizer in the optical path) at 1.63 µm. The PSF was computed over an area of 29 × 29 pixels (including light up to the third airy ring) with the IRDIS plate scale of 12.27 mas pix⁻¹, and was then normalized by the total integrated flux in the whole 29 × 29 pixels area. A 3D view of the resulting normalized PSF is shown in Fig. 9. We then fitted a 2D Moffat profile (Moffat 1969) of the form:

$$I(x, y) = I_0/(U + 1)^\beta, \quad (4)$$

where U corresponds to:

$$U = [(x - x_0)/A_x]^2 + [(y - y_0)/A_y]^2. \quad (5)$$

The resulting coefficients are listed in Table 2. We then repeated the Moffat fit, but this time for a total of 332 observed FLUX frames. We note that, in this case, we also included brighter stars (leading to a higher SR) that were not part of the sample presented in Sect. 2, only to cover the full range of FWHM and SR values. The resulting SR was simply estimated by taking the ratio between the peak of the Moffat fit (corresponding to $I_0$ in Eq. (4)) of the observed FLUX frames and that from the theoretical PSF. An example of an observed PSF leading to SR = 0.78 and its corresponding Moffat fit is shown in Fig. 10.

In addition, we measured the FWHM to the FLUX frames (as explained in Sect. 4.1), and we compared them with the measured SR. The results are presented in Fig. 11. We finally obtained a correlation between the SR and FWHM by fitting a function of the form:

$$\text{SR} = a + b/(\text{FWHM} + c)^\alpha \quad (6)$$

<table><caption>**Table 2.** Resulting coefficients for the Moffat fit of the theoretical PSF at 1.63 µm, following Eqs. (4) and (5).</caption><thead><tr><th>Coefficient</th><th>Value</th></tr></thead><tbody><tr><td>$I_0$</td><td>0.058</td></tr><tr><td>$\beta$</td><td>26710.9</td></tr><tr><td>$A_x$</td><td>347.5</td></tr><tr><td>$A_y$</td><td>345.7</td></tr><tr><td>$x_0$</td><td>14.0</td></tr><tr><td>$y_0$</td><td>14.0</td></tr></tbody></table>

<table><caption>**Table 3.** Resulting coefficients for the SR-FWHM correlation presented in Eq. (6).</caption><thead><tr><th>Coefficient</th><th>Value</th></tr></thead><tbody><tr><td>a</td><td>-0.099</td></tr><tr><td>b</td><td>5.296</td></tr><tr><td>c</td><td>-33.236</td></tr><tr><td>α</td><td>0.716</td></tr></tbody></table>

The resulting coefficients are listed in Table 3. The RMS of the post-fit residuals is 0.023. The best-fit is also over plotted in Fig. 11. Using this correlation, we then converted the FWHM values into SR values. We note that this conversion is more accurate for faint stars, which leads to relatively low SR values. On the other hand, for brighter targets that lead to a SR ≥ 0.5, this conversion is less suited, since in this regime the FWHM varies very slowly, asymptotically approaching the full diffraction limit value of ~ 43 mas. We therefore replicated Figs. 4 and 8, but this time we replaced the FWHM by the SR. The results are presented in Figs. 12 and 13. As expected, we observed a similar result to that seen in the raw contrast and FWHM, that is, a plateau for stars brighter than G ~ 12.5 mag, and a sharp decrease in the SR beyond this value. By comparing our findings with the expected AO performance for SPHERE and commissioning data (see Fig. 13 in Fusco et al. 2016 and Fig. 7 in Beuzit et al. 2019), we found that, although the overall trend in the SR with the stellar magnitude is also predicted, we see two main differences. Firstly, we found a nearly constant SR between G ~ 11 and G ~ 12.5 mag, and then a sharp decrease, at the 35% per magnitude level, beyond this value. In contrast, they found a nearly constant decrease in SR between R ~ 11–16 mag, at the ~20% per magnitude level. Secondly, we measured an overall lower SR in the full range of stellar flux. However, this is partly explained by the different methods used to compute this quantity (see Roberts et al. 2004). In addition, we note that they used the R-band instead of the G-band, which might introduce some small systematic differences, although it cannot solely explain this trend. Moreover, other parameters such as the WFS gain, the frequency or spatial filter used, as well as the atmospheric conditions during the commissioning (which are not discussed in their paper), might be partially responsible for these differences.

## 5. Final contrast

To investigate the achievable final contrast as a function of the stellar magnitude, we compared the raw contrast with the contrast measured on post-processing data. We included in this analysis 19 different stars, leading to RC₃₀₀ values covering roughly
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Observed PSF of a star leading to a Strehl ratio of 0.78 (**left**) and its corresponding Moffat fit (**right**).</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Strehl ratio versus the FWHM measured from FLUX frame. The solid line corresponds to the fit following Eq. (6).</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Strehl ratio as a function of the visible WFS flux, using the SR-FWHM correlation. The black and red lines correspond to the atmospheric turbulence categories TCAT10 and TCAT30, respectively.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Strehl ratio as a function of the stellar G-band magnitude, using the SR-FWHM correlation. The solid black and red dots (lines) correspond to the atmospheric turbulence categories TCAT10 and TCAT30, respectively.</figcaption></figure>

two orders of magnitude (between ~2 × 10⁻² and 2 × 10⁻⁴). For this we retrieved the post-processing contrast curves (corrected for the throughput of the algorithm) produced by the SPHERE Data Center (Delorme et al. 2017), for the classical Angular Differential Imaging (cADI; Marois et al. 2006) and the Template Local Optimized Combination of Images (TLOCI; Marois et al. 2010; Galicher et al. 2018) techniques. We note that we only included relatively long observing sequences, leading to a total field of view (FoV) rotation >15 deg, to avoid strong self-subtraction, and thus avoid a low throughput of the reduction algorithm. These results are presented in Fig. 14. The black dots and blue stars correspond to cADI and TLOCI reduction techniques. The dashed black and blue lines correspond to the median conversion factors, which were 2.7 and 12.4 for the cADI and TLOCI techniques, respectively. We could then use this result to estimate the final contrast on post-processing data by scaling RC₃₀₀ using these two values (to convert Figs. 3 and 7 into final contrast). Finally, we repeated this conversion, but this
<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Comparison between the contrast measured at 300 mas from individual raw frames and the post-reduction full observing sequence. The black dots and blue stars correspond to cADI and TLOCI reduction techniques, respectively.</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Fig. 15.** Same as for Fig. 14, but for the final contrast at 500 mas.</figcaption></figure>

time from $RC_{300}$ to final contrast at 500 mas. The results are presented in Fig. 15. Similarly, we obtained median conversion factors of 11.0 and 43.5 for the cADI and TLOCI techniques, respectively.

# 6. Summary and discussion

In this paper, we have analyzed a sample of 111 relatively faint stars ($G \ge 11.0$ mag) observed with IRDIS/SPHERE between 2016 and 2022, with the aim of investigating how the AO correction (and thus the final performance in terms of achievable contrast) degrades when approaching the limiting magnitude regime. For this, we measured the H-band $RC_{300}$ from coronagraphic frames and we compared this quantity to the instant flux received by the visible WFS. We found that the $RC_{300}$ is at the ~$10^{-3}$ level for $f_{\text{WFS}} \ge 3$ [e$^{-1}$/subaperture/frame]. For $f_{\text{WFS}}$ below this value, the AO correction strongly degrades, decreasing the $RC_{300}$ by a factor of roughly ten when reaching the 1 [e$^{-1}$/subaperture/frame] level. In terms of stellar magnitude, this turning point is observed at $G \sim 12.5$ mag. At $G \sim 14$ mag, the measured $RC_{300}$ is at the $10^{-2}$ level. We observed this behavior in the data taken under both TCAT10 and TCAT30 atmospheric conditions, with the former performing slightly better.

Similarly, we computed the FWHM and SR from the FLUX frames. We observed a very similar trend above and below the ~3 [e$^{-1}$/subaperture/frame] limit. In particular, we measured an increase of a factor of roughly two and three in the FWHM from this limit and the ~1 [e$^{-1}$/subaperture/frame] $f_{\text{WFS}}$ level, for TCAT10 and TCAT30, respectively. Similarly, the SR drops from above 0.5 to less than 0.2, for both TCAT10 and TCAT30.

Finally, we computed the final contrast at 300 mas and 500 mas, in a selected sample of stars, using the post-reduction contrast curves. We obtained median conversion factors of roughly three and 12, at 300 mas, for the cADI and TLOCI techniques, respectively. Similarly, we found scaling factors of ~11 and 44, to 500 mas, for cADI and TLOCI, respectively. Using these conversion factors, we can then estimate the final contrast achievable at 300 mas and 500 mas, by scaling the results presented in this work.

As an example, we can consider a star with the exact same luminosity and colors of the well-known planet host star PDS70 (Keppler et al. 2018). We can also consider that such a star is located at ~213 pc from the Earth, so that it has: $G = 13.0$ mag and $H = 11.2$ mag. Following Fig. 7, we can expect a $RC_{300}$ of ~$4 \times 10^{-3}$ under good atmospheric conditions (TCAT30 or better). If we perform a relatively long (FoV rotation $\ge$15 deg) ADI sequence, we can then estimate a $5\sigma$ final contrast at 500 mas using the TLOCI technique of ~$10^{-4}$. If we compare these numbers with the AMES-COND models (Allard et al. 2001), this limit corresponds to a 10-Myr planet with a ~4.8 $M_J$ mass, located at a physical separation of ~107 AU from the host star.

Improving the contrast on faint targets and reaching fainter or redder stars represent one of the primary scientific drivers for the instrument upgrade called SPHERE+ (Boccaletti et al. 2022). Indeed, ALMA surveys such as DSHARP (Andrews et al. 2018) revealed dozens of star-forming disks with structures such as gaps, rings, or spirals, suggesting planet formation is on-going. PDS70 represents one of the brightest examples of this population, the bulk of which is currently out of reach for SPHERE (see Fig. 12 of Boccaletti et al. 2020). Coupling the current AO correction stage to a more sensitive second-stage AO in cascade represents the main technical improvement of SPHERE+ (Cerpa-Urra et al. 2022), and the characterization of the current performance on faint targets presented here is crucial in order to make informed technical choices in the design of the second stage AO for SPHERE+.

**Acknowledgements.** This work has made use of the SPHERE Data Centre, jointly operated by OSUG/IPAG (Grenoble), PYTHEAS/LAM/CESAM (Marseille), OCA/Lagrange (Nice), Observatoire de Paris/LESIA (Paris), and Observatoire de Lyon. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. 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

Allard, F., Hauschildt, P. H., Alexander, D. R., Tamanai, A., & Schweitzer, A. 2001, ApJ, 556, 357

Andrews, S. M., Huang, J., Pérez, L. M., et al. 2018, ApJ, 869, L41

Bessell, M. S. 1990, PASP, 102, 1181

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

Boccaletti, A., Chauvin, G., Mouillet, D., et al. 2020, ESO, submitted [arXiv:2003.05714]

Boccaletti, A., Chauvin, G., Wildi, F., et al. 2022, ArXiv e-prints [arXiv:2209.02092]

Bohn, A. J., Kenworthy, M. A., Ginski, C., et al. 2020, ApJ, 898, L16
Bohn, A. J., Ginski, C., Kenworthy, M. A., et al. 2021, A&A, 648, A73

Cerpa-Urra, N., Kasper, M., Kulcsár, C., Raynaud, H.-F., & Taïssir Heritier, C. 2022, J. Astron. Telescopes Instrum. Syst., 8, 019001

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, 2MASS All Sky Catalog of point sources.

Delorme, P., Meunier, N., Albert, D., et al. 2017, in SF2A-2017: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, eds. C. Reylé, P. Di Matteo, F. Herpin, E. Lagadec, A. Lançon, Z. Meliani, & F. Royer, Di

Desidera, S., Chauvin, G., Bonavita, M., et al. 2021, A&A, 651, A70

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

Fedrigo, E., Donaldson, R., Soenke, C., et al. 2006, SPIE Conf. Ser., 6272, 627210

Fusco, T., Petit, C., Rousset, G., et al. 2006, SPIE Conf. Ser., 6272, 62720K

Fusco, T., Sauvage, J. F., Mouillet, D., et al. 2016, SPIE Conf. Ser., 9909, 99090U

Galicher, R., Boccaletti, A., Mesa, D., et al. 2018, A&A, 615, A92

Janson, M., Squicciarini, V., Delorme, P., et al. 2021, A&A, 646, A164

Keppler, M., Benisty, M., Müller, A., et al. 2018, A&A, 617, A44

Lagrange, A. M., Keppler, M., Meunier, N., et al. 2018, A&A, 612, A108

Marois, C., Lafrenière, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006, ApJ, 641, 556

Marois, C., Macintosh, B., & Véran, J.-P. 2010, SPIE Conf. Ser., 7736, 77361J

Milli, J., Kasper, M., Bourget, P., et al. 2018, SPIE Conf. Ser., 10703, 107032A

Milli, J., Hibon, P., Christiaens, V., et al. 2017a, A&A, 597, A2

Milli, J., Mouillet, D., Fusco, T., et al. 2017b, Adaptive Optics for Extremely Large Telescopes 5, Conference Proceeding, Tenerife, Canary Islands, Spain https://doi.org/10.26698/A04ELT5.0034

Moffat, A. F. J. 1969, A&A, 3, 455

Müller, A., Keppler, M., Henning, T., et al. 2018, A&A, 617, A2

Patat, F., Moehler, S., O'Brien, K., et al. 2011, A&A, 527, A91

Riello, M., De Angeli, F., Evans, D. W., et al. 2021, A&A, 649, A3

Roberts, Lewis C. J., Perrin, M. D., Marchis, F., et al. 2004, SPIE Conf. Ser., 5490, 504

Sauvage, J.-F., Fusco, T., Lamb, M., et al. 2016a, SPIE Conf. Ser., 9909, 990916

Sauvage, J.-F., Fusco, T., Petit, C., et al. 2016b, J. Astron. Telescopes Instrum. Syst., 2, 025003

van Boekel, R., Henning, T., Menu, J., et al. 2017, ApJ, 837, 132

Zurlo, A., Vigan, A., Galicher, R., et al. 2016, A&A, 587, A57
Appendix A: Targets' apparent magnitudes

<table><caption>Table A.1. Gaia EDR3 and 2MASS JHK stellar magnitudes, listed in ascending order.</caption><thead><tr><th>G</th><th>G<sub>B</sub></th><th>G<sub>R</sub></th><th>J</th><th>H</th><th>K</th><th>Star name</th></tr></thead><tbody><tr><td>10.99</td><td>11.81</td><td>10.00</td><td>8.4</td><td>7.7</td><td>7.2</td><td>2MASS-J16274028-2422040</td></tr><tr><td>11.08</td><td>11.51</td><td>10.43</td><td>9.7</td><td>9.2</td><td>9.1</td><td>2MASS-J13015435-4249422</td></tr><tr><td>11.14</td><td>11.88</td><td>10.06</td><td>8.7</td><td>8.1</td><td>7.9</td><td>TYC-8083-455-1</td></tr><tr><td>11.16</td><td>11.88</td><td>10.32</td><td>9.4</td><td>8.7</td><td>8.5</td><td>2MASS-J10560422-6152054</td></tr><tr><td>11.20</td><td>11.83</td><td>10.42</td><td>9.3</td><td>8.5</td><td>8.0</td><td>2MASS-J11122441-7637064</td></tr><tr><td>11.24</td><td>11.86</td><td>10.47</td><td>9.5</td><td>8.9</td><td>8.8</td><td>2MASS-J12383556-5916438</td></tr><tr><td>11.25</td><td>11.68</td><td>10.39</td><td>9.4</td><td>8.8</td><td>8.7</td><td>2MASS-J13121764-5508258</td></tr><tr><td>11.34</td><td>12.67</td><td>10.19</td><td>...</td><td>...</td><td>...</td><td>HIP20990</td></tr><tr><td>11.34</td><td>12.06</td><td>10.50</td><td>8.6</td><td>7.9</td><td>7.7</td><td>2MASS-J00172353-6645124</td></tr><tr><td>11.45</td><td>12.14</td><td>10.62</td><td>8.6</td><td>8.0</td><td>7.8</td><td>2MASS-J12072738-3247002</td></tr><tr><td>11.45</td><td>12.88</td><td>10.27</td><td>9.4</td><td>8.6</td><td>8.2</td><td>2MASS-J11175186-6402056</td></tr><tr><td>11.46</td><td>12.34</td><td>10.51</td><td>9.1</td><td>8.3</td><td>8.0</td><td>DN-Tau</td></tr><tr><td>11.53</td><td>12.90</td><td>10.54</td><td>8.7</td><td>7.7</td><td>7.1</td><td>HQ-Tau</td></tr><tr><td>11.54</td><td>12.35</td><td>10.62</td><td>9.1</td><td>8.3</td><td>7.8</td><td>CHX-18N</td></tr><tr><td>11.56</td><td>11.60</td><td>11.46</td><td>11.3</td><td>11.3</td><td>11.3</td><td>2MASS-J21012926-0622148</td></tr><tr><td>11.59</td><td>12.27</td><td>10.77</td><td>9.6</td><td>9.0</td><td>8.9</td><td>2MASS-J16015918-3612555</td></tr><tr><td>11.61</td><td>12.39</td><td>10.71</td><td>9.6</td><td>8.8</td><td>8.5</td><td>CD-40-8434</td></tr><tr><td>11.62</td><td>12.36</td><td>10.77</td><td>9.7</td><td>9.0</td><td>8.8</td><td>2MASS-J12510556-5253121</td></tr><tr><td>11.63</td><td>12.20</td><td>10.90</td><td>9.9</td><td>9.4</td><td>9.3</td><td>2MASS-J12404664-5211046</td></tr><tr><td>11.65</td><td>12.87</td><td>10.55</td><td>9.1</td><td>8.4</td><td>8.2</td><td>2MASS-J22025453-6440441</td></tr><tr><td>11.65</td><td>12.35</td><td>10.73</td><td>8.8</td><td>7.8</td><td>6.9</td><td>2MASS-J04470620+1658428</td></tr><tr><td>11.67</td><td>12.55</td><td>10.73</td><td>9.3</td><td>8.6</td><td>8.4</td><td>2MASS-J04331003+2433433</td></tr><tr><td>11.68</td><td>12.65</td><td>10.70</td><td>9.5</td><td>8.7</td><td>8.6</td><td>2MASS-J04440099-6624036</td></tr><tr><td>11.72</td><td>12.41</td><td>10.90</td><td>9.5</td><td>8.7</td><td>8.2</td><td>2MASS-J18425797-3532427</td></tr><tr><td>11.73</td><td>12.41</td><td>10.79</td><td>9.3</td><td>8.6</td><td>8.3</td><td>2MASS-J04551098+3021595</td></tr><tr><td>11.73</td><td>12.47</td><td>10.71</td><td>9.3</td><td>8.4</td><td>7.8</td><td>SZ-cha</td></tr><tr><td>11.75</td><td>12.71</td><td>10.75</td><td>9.2</td><td>8.4</td><td>8.0</td><td>2MASS-J15392776-3446171</td></tr><tr><td>11.77</td><td>12.02</td><td>11.35</td><td>10.9</td><td>10.6</td><td>10.5</td><td>2MASS-J00190738-3712352</td></tr><tr><td>11.77</td><td>12.90</td><td>10.71</td><td>9.1</td><td>8.3</td><td>8.0</td><td>2MASS-J16303563-2434188</td></tr><tr><td>11.80</td><td>12.59</td><td>10.89</td><td>9.7</td><td>8.9</td><td>8.7</td><td>2MASS-J11040909-7627193</td></tr><tr><td>11.82</td><td>13.61</td><td>10.54</td><td>8.7</td><td>8.1</td><td>7.8</td><td>TX-PsA</td></tr><tr><td>11.84</td><td>12.54</td><td>11.01</td><td>10.0</td><td>9.3</td><td>9.2</td><td>RX-J1555.6-3709</td></tr><tr><td>11.86</td><td>12.59</td><td>11.01</td><td>10.0</td><td>9.3</td><td>9.2</td><td>GSC-06214-00210</td></tr><tr><td>11.88</td><td>12.86</td><td>10.77</td><td>8.7</td><td>7.7</td><td>7.0</td><td>DG-Tau</td></tr><tr><td>11.90</td><td>13.26</td><td>10.73</td><td>9.1</td><td>8.5</td><td>8.3</td><td>2MASS-J02014693+0117059</td></tr><tr><td>11.91</td><td>13.07</td><td>10.81</td><td>9.2</td><td>8.3</td><td>7.8</td><td>DE-Tau</td></tr><tr><td>11.91</td><td>12.94</td><td>10.88</td><td>9.5</td><td>8.8</td><td>8.6</td><td>2MASS-J08422372-7904030</td></tr><tr><td>11.92</td><td>12.39</td><td>10.85</td><td>9.3</td><td>8.5</td><td>8.1</td><td>2MASS-J02581122+2030037</td></tr><tr><td>11.97</td><td>12.57</td><td>11.21</td><td>10.2</td><td>9.7</td><td>9.5</td><td>GSC-08584-01898</td></tr><tr><td>11.99</td><td>13.03</td><td>10.89</td><td>9.1</td><td>8.1</td><td>7.5</td><td>GK-Tau</td></tr><tr><td>12.02</td><td>12.39</td><td>11.47</td><td>10.8</td><td>10.4</td><td>10.3</td><td>2MASS-J06281844-7621467</td></tr><tr><td>12.04</td><td>12.92</td><td>11.11</td><td>9.8</td><td>9.1</td><td>8.9</td><td>2MASS-J16093030-2104589</td></tr><tr><td>12.05</td><td>12.81</td><td>11.10</td><td>9.5</td><td>8.6</td><td>8.0</td><td>DS-Tau</td></tr><tr><td>12.07</td><td>12.93</td><td>11.12</td><td>9.6</td><td>8.8</td><td>8.4</td><td>2MASS-J16035767-2031055</td></tr><tr><td>12.10</td><td>13.02</td><td>11.12</td><td>9.7</td><td>9.0</td><td>8.9</td><td>2MASS-J16071370-3839238</td></tr><tr><td>12.14</td><td>13.01</td><td>11.11</td><td>9.6</td><td>8.7</td><td>8.0</td><td>DL-Tau</td></tr><tr><td>12.18</td><td>13.23</td><td>11.11</td><td>9.5</td><td>8.4</td><td>7.8</td><td>CI-Tau</td></tr><tr><td>12.20</td><td>12.50</td><td>11.68</td><td>10.6</td><td>9.7</td><td>8.8</td><td>2MASS-J05113654-0222484</td></tr><tr><td>12.22</td><td>14.11</td><td>10.89</td><td>8.1</td><td>6.9</td><td>6.2</td><td>Haro-1-6</td></tr><tr><td>12.23</td><td>13.09</td><td>11.21</td><td>9.5</td><td>8.7</td><td>8.0</td><td>V*-V1279-Sco</td></tr><tr><td>12.23</td><td>13.51</td><td>11.06</td><td>8.7</td><td>7.2</td><td>6.1</td><td>WW-Cha</td></tr><tr><td>12.25</td><td>13.06</td><td>11.19</td><td>9.4</td><td>8.3</td><td>7.9</td><td>2MASS-J16153456-2242421</td></tr><tr><td>12.26</td><td>13.28</td><td>11.24</td><td>9.5</td><td>8.6</td><td>8.1</td><td>2MASS-J16221852-2321480</td></tr><tr><td>12.36</td><td>13.21</td><td>11.48</td><td>10.1</td><td>9.3</td><td>8.7</td><td>2MASS-J05383368-0244141</td></tr><tr><td>12.36</td><td>13.26</td><td>11.39</td><td>9.9</td><td>9.1</td><td>8.8</td><td>2MASS-J13042410-7650012</td></tr><tr><td>12.37</td><td>13.28</td><td>11.39</td><td>9.9</td><td>9.1</td><td>8.6</td><td>2MASS-J15591647-4157102</td></tr></tbody></table>

Continue on the next page
<table><caption>**Table A.1.** Gaia EDR3 and 2MASS JHK stellar magnitudes, listed in ascending order.</caption><thead><tr><th>G</th><th>G<sub>B</sub></th><th>G<sub>R</sub></th><th>J</th><th>H</th><th>K</th><th>Star name</th></tr></thead><tbody><tr><td>12.37</td><td>13.56</td><td>11.26</td><td>...</td><td>...</td><td>...</td><td>UZ-Tau</td></tr><tr><td>12.37</td><td>13.28</td><td>11.39</td><td>9.9</td><td>9.1</td><td>8.6</td><td>WRAY-15-1400</td></tr><tr><td>12.38</td><td>13.34</td><td>11.28</td><td>9.5</td><td>8.2</td><td>7.3</td><td>DO-Tau</td></tr><tr><td>12.38</td><td>13.33</td><td>11.40</td><td>9.8</td><td>9.1</td><td>8.5</td><td>2MASS-J10180498-5816263</td></tr><tr><td>12.42</td><td>13.43</td><td>11.41</td><td>9.8</td><td>8.9</td><td>8.3</td><td>IP-Tau</td></tr><tr><td>12.42</td><td>12.54</td><td>12.19</td><td>12.0</td><td>11.9</td><td>11.8</td><td>TYC-9337-2511-1</td></tr><tr><td>12.46</td><td>13.69</td><td>11.34</td><td>9.5</td><td>8.5</td><td>8.0</td><td>DQ-Tau</td></tr><tr><td>12.46</td><td>13.47</td><td>11.42</td><td>9.8</td><td>8.8</td><td>8.2</td><td>DH-Tau</td></tr><tr><td>12.49</td><td>13.41</td><td>11.51</td><td>10.1</td><td>9.3</td><td>9.0</td><td>2MASS-J16100501-2132318</td></tr><tr><td>12.52</td><td>13.50</td><td>11.06</td><td>9.1</td><td>8.1</td><td>7.5</td><td>Haro-1-4</td></tr><tr><td>12.52</td><td>14.20</td><td>11.27</td><td>8.7</td><td>7.5</td><td>6.9</td><td>2MASS-J16283266-2422449</td></tr><tr><td>12.53</td><td>14.00</td><td>11.27</td><td>10.7</td><td>9.6</td><td>9.0</td><td>V409-Tau</td></tr><tr><td>12.54</td><td>13.54</td><td>11.53</td><td>9.9</td><td>9.0</td><td>8.7</td><td>V1094-Sco</td></tr><tr><td>12.55</td><td>13.53</td><td>11.56</td><td>10.0</td><td>9.2</td><td>8.5</td><td>PDS-65</td></tr><tr><td>12.57</td><td>13.58</td><td>11.55</td><td>10.2</td><td>9.5</td><td>9.3</td><td>2MASS-J13342461-6517473</td></tr><tr><td>12.60</td><td>13.88</td><td>11.47</td><td>9.9</td><td>9.1</td><td>8.8</td><td>CX-Tau</td></tr><tr><td>12.61</td><td>13.54</td><td>11.65</td><td>10.4</td><td>9.7</td><td>9.6</td><td>2MASS-J15545286-3827566</td></tr><tr><td>12.61</td><td>13.87</td><td>11.51</td><td>9.3</td><td>8.4</td><td>7.9</td><td>GI-Tau</td></tr><tr><td>12.63</td><td>13.58</td><td>11.61</td><td>10.1</td><td>9.4</td><td>9.2</td><td>2MASS-J16064385-1908056</td></tr><tr><td>12.64</td><td>13.71</td><td>11.60</td><td>9.9</td><td>9.1</td><td>8.6</td><td>V836-Tau</td></tr><tr><td>12.66</td><td>14.33</td><td>11.40</td><td>8.7</td><td>7.5</td><td>6.7</td><td>2MASS-J16271027-2419127</td></tr><tr><td>12.72</td><td>13.40</td><td>11.73</td><td>9.9</td><td>9.1</td><td>8.5</td><td>2MASS-J05065551-0321132</td></tr><tr><td>12.72</td><td>14.06</td><td>11.57</td><td>9.9</td><td>9.2</td><td>8.9</td><td>FP-Tau</td></tr><tr><td>12.73</td><td>13.80</td><td>11.70</td><td>10.1</td><td>9.5</td><td>9.3</td><td>2MASS-J18443116-3723347</td></tr><tr><td>12.74</td><td>12.84</td><td>12.55</td><td>12.4</td><td>12.3</td><td>12.3</td><td>WD-1202-232</td></tr><tr><td>12.76</td><td>13.67</td><td>11.80</td><td>10.1</td><td>9.2</td><td>8.6</td><td>TW-Cha</td></tr><tr><td>12.77</td><td>14.15</td><td>11.56</td><td>9.4</td><td>8.4</td><td>7.8</td><td>2MASS-J16262367-2443138</td></tr><tr><td>12.77</td><td>13.69</td><td>11.55</td><td>9.6</td><td>8.2</td><td>7.1</td><td>CW-Tau</td></tr><tr><td>12.79</td><td>14.02</td><td>11.68</td><td>9.6</td><td>8.6</td><td>7.6</td><td>PDS-198</td></tr><tr><td>12.80</td><td>14.26</td><td>11.61</td><td>9.4</td><td>8.3</td><td>7.9</td><td>2MASS-J16244839-2359160</td></tr><tr><td>12.85</td><td>13.70</td><td>11.96</td><td>10.6</td><td>9.8</td><td>9.4</td><td>2MASS-J05401286-0222020</td></tr><tr><td>12.89</td><td>14.02</td><td>11.83</td><td>10.2</td><td>9.5</td><td>9.2</td><td>2MASS-J16090075-1908526</td></tr><tr><td>12.92</td><td>14.12</td><td>11.82</td><td>10.2</td><td>9.5</td><td>9.2</td><td>2MASS-J16111534-1757214</td></tr><tr><td>12.97</td><td>13.09</td><td>12.75</td><td>12.7</td><td>12.6</td><td>12.6</td><td>LAWD-25</td></tr><tr><td>12.98</td><td>14.64</td><td>11.73</td><td>...</td><td>...</td><td>...</td><td>V710-Tau</td></tr><tr><td>13.03</td><td>14.32</td><td>11.78</td><td>9.5</td><td>8.5</td><td>7.6</td><td>HP-Tau</td></tr><tr><td>13.04</td><td>14.02</td><td>12.02</td><td>10.3</td><td>9.6</td><td>9.2</td><td>Ass-Cha-T-1-32</td></tr><tr><td>13.06</td><td>13.08</td><td>13.05</td><td>13.1</td><td>13.1</td><td>12.7</td><td>2MASS-J23284760+0514540</td></tr><tr><td>13.08</td><td>13.19</td><td>12.88</td><td>...</td><td>...</td><td>...</td><td>CRTS-J232812.6-395523</td></tr><tr><td>13.10</td><td>14.07</td><td>12.04</td><td>10.4</td><td>9.8</td><td>9.5</td><td>DM-Tau</td></tr><tr><td>13.10</td><td>14.34</td><td>11.96</td><td>10.1</td><td>9.2</td><td>8.6</td><td>2MASS-J15464473-3430354</td></tr><tr><td>13.13</td><td>13.84</td><td>12.32</td><td>11.3</td><td>10.7</td><td>10.6</td><td>2MASS-J06000075-3100276</td></tr><tr><td>13.19</td><td>14.28</td><td>12.06</td><td>10.0</td><td>8.9</td><td>8.4</td><td>2MASS-J16154416-1921171</td></tr><tr><td>13.19</td><td>14.49</td><td>12.02</td><td>9.4</td><td>8.4</td><td>7.8</td><td>IQ-Tau</td></tr><tr><td>13.33</td><td>14.65</td><td>12.18</td><td>10.2</td><td>9.3</td><td>8.9</td><td>2MASS-J13045571-7739495</td></tr><tr><td>13.34</td><td>14.51</td><td>12.17</td><td>10.5</td><td>9.8</td><td>9.4</td><td>ET-Cha</td></tr><tr><td>13.38</td><td>14.10</td><td>12.49</td><td>10.7</td><td>9.5</td><td>8.4</td><td>HN-Tau</td></tr><tr><td>13.45</td><td>14.26</td><td>12.53</td><td>11.3</td><td>10.4</td><td>10.0</td><td>V499-Ori</td></tr><tr><td>13.47</td><td>17.06</td><td>11.39</td><td>7.3</td><td>6.3</td><td>5.8</td><td>2MASS-J12450626-7720133</td></tr><tr><td>13.50</td><td>14.76</td><td>12.35</td><td>10.6</td><td>9.8</td><td>9.3</td><td>SZ-72</td></tr><tr><td>13.52</td><td>14.36</td><td>12.59</td><td>11.1</td><td>10.3</td><td>9.8</td><td>Haro-5-34</td></tr><tr><td>13.54</td><td>15.13</td><td>12.13</td><td>10.0</td><td>8.9</td><td>8.4</td><td>IRAS-12535-7623</td></tr><tr><td>13.66</td><td>14.43</td><td>12.79</td><td>11.5</td><td>10.6</td><td>10.0</td><td>2MASS-J05394017-0220480</td></tr><tr><td>13.66</td><td>14.49</td><td>12.75</td><td>11.7</td><td>10.8</td><td>10.2</td><td>V543-ORI</td></tr><tr><td>13.95</td><td>15.45</td><td>12.67</td><td>10.2</td><td>8.7</td><td>7.6</td><td>2MASS-J16394544-2402039</td></tr></tbody></table>