# Robust high-contrast companion detection from interferometric observations

Alexandre Gallenne, Antoine Mérand, Pierre Kervella, J. Monnier, G.

Schaefer, F. Baron, J. Breitfelder, J. Le Bouquin, R. Roettenbacher, W. Gieren, et al.

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

## ► To cite this version:

Alexandre Gallenne, Antoine Mérand, Pierre Kervella, J. Monnier, G. Schaefer, et al. Robust high-contrast companion detection from interferometric observations. Astronomy & Astrophysics - A&A, 2015, 579, pp.A68. 10.1051/0004-6361/201525917. hal-02505418

HAL Id: hal-02505418

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

Submitted on 9 Jan 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Robust high-contrast companion detection from interferometric observations

## The CANDID algorithm and an application to six binary Cepheids*

A. Gallenne<fnref n="1" />, A. Mérand<fnref n="2" />, P. Kervella<fnref n="3" />,<fnref n="4" />, J. D. Monnier<fnref n="5" />, G. H. Schaefer<fnref n="6" />, F. Baron<fnref n="6" />,<fnref n="7" />, J. Breitfelder<fnref n="3" />,
J. B. Le Bouquin<fnref n="8" />, R. M. Roettenbacher<fnref n="5" />, W. Gieren<fnref n="1" />,<fnref n="10" />, G. Pietrzyński<fnref n="1" />,<fnref n="9" />, H. McAlister<fnref n="6" />, T. ten Brummelaar<fnref n="6" />,
J. Sturmann<fnref n="6" />, L. Sturmann<fnref n="6" />, N. Turner<fnref n="6" />, S. Ridgway<fnref n="11" />, and S. Kraus<fnref n="12" />

Received 17 February 2015 / Accepted 4 May 2015

## ABSTRACT

**Context.** Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By combining several telescopes together and spectrally dispersing the light, it is possible to detect faint components around bright stars in a few hours of observations.

**Aims.** We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, and estimate the dynamic range from interferometric observations.

**Methods.** We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit using all interferometric observables, i.e., the squared visibilities, closure phases and bispectrum amplitudes. The search procedure is made on a $N \times N$ grid of fit, whose minimum needed resolution is estimated a posteriori. It includes a tool to estimate the detection level of the companion in the number of sigmas. The code CANDID also incorporates a robust method to set a $3\sigma$ detection limit on the flux ratio, which is based on an analytical injection of a fake companion at each point in the grid. Our injection method also allows us to analytically remove a detected component to 1) search for a second companion; and 2) set an unbiased detection limit.

**Results.** We used CANDID to search for the companions around the binary Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul. First, we showed that our previous discoveries of the components orbiting V1334 Cyg and AX Cir were detected at $>25\sigma$ and $>13\sigma$, respectively. The astrometric positions and flux ratios provided by CANDID for these two stars are in good agreement with our previously published values. The companion around AW Per is detected at more than $15\sigma$ with a flux ratio of $f = 1.22 \pm 0.30\%$, and it is located at $\rho = 32.16 \pm 0.29$ mas and PA = $67.1 \pm 0.3^\circ$. We made a possible detection of the companion orbiting RT Aur with $f = 0.22 \pm 0.11\%$, and at $\rho = 2.10 \pm 0.23$ mas and PA = $-136 \pm 6^\circ$. It was detected at $3.8\sigma$ using the closure phases only, and so more observations are needed to confirm the detection. No companions were detected around SU Cas and T Vul. We also set the detection limit for possible undetected companions around these stars. We found that there is no companion with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V orbiting the Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively. This work also demonstrates the capabilities of the MIRC and PIONIER instruments, which can reach a dynamic range of 1:200, depending on the angular distance of the companion and the $(u, v)$ plane coverage. In the future, we plan to work on improving the sensitivity limits for realistic data through better handling of the correlations.

**Key words.** stars: variables: Cepheids – binaries: close – techniques: high angular resolution – techniques: interferometric

## 1. Introduction

Long-baseline interferometry (LBI) enables us to spatially resolve components in close orbits, providing astrometric

positions at <50 milli-arcsec (mas) scale with micro-arcsecond accuracy. When combined with spectroscopic radial velocities, we can obtain model independent estimates of the stellar masses and orbital parallaxes, which are fundamental parameters that help us study stellar properties and evolution. However, so far LBI is limited to bright stars ($H < 7$ mag) with typical magnitude differences $\Delta H < 6$ mag. Adaptive optics imaging with

<footnote n="1">¹ Universidad de Concepción, Departamento de Astronomía, Casilla 160-C, Concepción, Chile
e-mail: agallenne@astro-udec.cl</footnote>

<footnote n="2">² European Southern Observatory, Alonso de Córdova 3107, Casilla 19001, Santiago 19, Chile</footnote>

<footnote n="3">³ LESIA, Observatoire de Paris, CNRS UMR 8109, UPMC, Université Paris Diderot, 5 place Jules Janssen, 92195 Meudon, France</footnote>

<footnote n="4">⁴ Unidad Mixta Internacional Franco-Chilena de Astronomía, CNRS/INSU, France (UMI 3386) and Departamento de Astronomía, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile</footnote>

<footnote n="5">⁵ Astronomy Department, University of Michigan, 941 Dennison Bldg, Ann Arbor, MI 48109-1090, USA</footnote>

<footnote n="6">⁶ The CHARA Array of Georgia State University, Mount Wilson, CA 91023, USA</footnote>

<footnote n="7">⁷ Department of Physics & Astronomy, Georgia State University, 25 Park Place NE, Atlanta, GA 30303-2911, USA</footnote>

<footnote n="8">⁸ Université Grenoble Alpes, IPAG, 38000 Grenoble, France and CNRS, IPAG, 38000 Grenoble, France</footnote>

<footnote n="9">⁹ Warsaw University Observatory, Al. Ujazdowskie 4, 00-478 Warsaw, Poland</footnote>

<footnote n="10">¹⁰ Millenium Institute of Astrophysics, Santiago, Chile</footnote>

<footnote n="11">¹¹ National Optical Astronomy Observatories, 950 North Cherry Avenue, Tucson, AZ 85719, USA</footnote>

<footnote n="12">¹² School of Physics, University of Exeter, Stocker Road, Exeter EX4 4QL, UK</footnote>

<footnote>* The current version of the code is available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/579/A68</footnote>
single-dish telescopes reaches better contrasts, down to $\Delta H \sim 12$ mag (Zurlo et al. 2014), but the angular resolution is limited to $0.2''$ at these detection levels. Long-baseline interferometry is therefore a complementary technique in terms of spatial scale by probing the innermost regions.

LBI can detect components down to ~1 mas in the infrared, but the main limitation is sensitivity to high-contrast companions. We roughly know the performances reachable by current interferometric combiners. So far, the faintest companion detected with LBI has a flux ratio of 0.75% and orbits a Cepheid star of magnitude $H = 3.85$ mag (Gallenne et al. 2014). Objects, such as brown dwarfs and hot giant planets, are still inaccessible because of a lack of sensitivity and accuracy of the instruments. Absil et al. (2011) demonstrated a possible dynamic range of 1:500 with the VLTI/PIONIER instrument (Le Bouquin et al. 2011), but this range has not yet been reached. To achieve this detection level, several hours of multitelescope observations are required to obtain as many simultaneous interferometric measurements as possible.

Deriving the sensitivity limit from imaging can be determined directly from the noise level, however, this is not the case for interferometric observations. There are some papers in the literature that discuss detection limits and methods to search for companions from interferometric data (Absil et al. 2011; Le Bouquin & Absil 2012), but these studies have some shortcomings: the searching method is not formalized, the sigma detection is not robust, and they do not take the bandwidth smearing into account. Therefore, a robust implementation to search for components does not exist thus far. This kind of method is particularly critical to detect faint companions, as they can be at the sensitivity limit of the instrument or even not be a statistically significant detection (i.e., $<3\sigma$). We therefore created CANDID (Companion Analysis and Non-Detection in Interferometric Data) to address these aspects. This is a suite of Python tools, which contains two main functions: 1) one to perform a systematic search for faint companions (Sect. 2); and 2) one to estimate the detection limit from long-baseline interferometric observations (Sect. 3). This tool is made available to the community<fnref n="1" />.

CANDID is made for detecting high-contrast, point-source companions orbiting a spatially resolved primary star, although it also works with an unresolved primary and contrast <50%. In this paper, we used CANDID to look for companions in binary Cepheids. We present the first main function to search for companions in Sect. 2, verify and clarify the detection level of our previously detected faint companions, and we report new detections for other Cepheids. In Sect. 3, we present the second main function of CANDID and explain our robust method to set detection limits from interferometric data. In Sect. 4, we then use this formalism to our set of Cepheid observations to derive the detection limits. We finally present our conclusions in Sect. 5.

## 2. Searching for companions

Our grid search is traditionally performed in a three-dimensional space. Specifically, we vary the position of the component ($\Delta\alpha, \Delta\delta$) and the companion/star flux ratio $f$, and then compute the $\chi^2$ for each of these positions and flux ratios (the third dimension being $f$). The weakness of this method is the resolution of the grid, i.e., if the grid is to too coarse, the detection can be missed.

In addition, for high-contrast binaries and/or low accuracy data, a $\chi^2$ map can show fake or nonsignificant detections. Sometimes, some authors show the estimated $\chi^2$ map with the most probable location of the component without any detection level mentioned, only the reduced $\chi^2$ variation is given. This parameter, however, is not optimal for checking the detection level because it depends on the number of degrees of freedom (d.o.f.). The large quantity of data required to detect a faint companion using LBI results in a large dof. For instance, we can have the most probable location with the lowest reduced chi-square, $\chi^2$, equal to 1.0, and the highest in the whole map equal to 1.1, i.e., with only a variation of 0.1. With this information, however, there is no way to know if the detection is statistically significant or not. In this section, we therefore present a systematic approach to search for components and set the detection level using the $\chi^2$ and the number of dof. We then applied it to the case of binary Cepheids. We verify our previous detections for the Cepheids V1334 Cyg and AX Cir first (Gallenne et al. 2013b, 2014), and then report new detections for RT Aur and AW Per.

### 2.1. Detection method

A more rigorous approach is to perform a grid of fit using a least-squares algorithm, with a starting grid spacing that is guaranteed to find the global minimum. The grid in question is the 2D grid of starting points for the companion position. For each starting position, a multiparameter fit is performed: the companion position and its flux ratio (possibly the stellar diameters) are adjusted. Each position of the grid leads to a local minimum. Ideally, if the starting grid is fine enough, multiple starting points lead to the same local minima: this guarantees that all the local minima are explored and that the global minimum is indeed the global minimum. Hence, the criteria to decide if the global minimum is global is a posteriori. In CANDID, we require (a posteriori) that on average, each unique minimum is reached from two starting points of the starting grid. We also provide statistics on the “traveling” distance of the fit (between the starting and end points), compared to the size of the grid. We require that the median travel distance should be less than $\sqrt{2}/2$ of the size of the square grid. The 10 and 90% percentiles are also provided to the user to assess the typical travel distance.

A systematic search using a grid of fit is an iterative process. First, a coarse starting grid is chosen and the fits are performed. CANDID estimates, based on the traveling distances of the fits, how to refine the starting grid. A second series of fits are run, using the finer grid, which might take much longer, but the global minimum can be trusted to be global (at least within the area searched).

We searched for components with a maximum distance to the main star of 50 mas. For a wider range, the loss of coherence caused by spectral smearing of the companion is the main limitation because it degrades the dynamic range. A complete discussion about the search region is presented by Absil et al. (2011) and Le Bouquin & Absil (2012). The main limitation is the spectral sampling compared to the relative position of the component, i.e., to avoid a significant smearing we need $R > \rho B_p / \lambda$, with $R$, $B_p$ and $\rho$ the spectral resolution, the projected baseline, and the separation, respectively. For PIONIER, with $R = 18$ and a mean projected baseline $B_p = 100$ m (used for our observations), we have $\rho < 60$ mas, and for MIRC with $R = 42$ and $B = 200$ m, we have $\rho < 70$ mas. Zhao et al. (2007) recommended a more stringent criteria $\rho < R\lambda/(5B)$ to assure uncorrupted data, leading to $\rho \le 15$ mas for PIONIER and MIRC.

<footnote n="1">¹ Available at https://github.com/amerand/CANDID</footnote>
However, a companion can still be detected at separations larger than 15 mas by taking bandwidth smearing effects into account. We chose 50 mas as our limiting range to avoid making the grid search too stringent. In addition, companions located at more than 50 mas are detected more efficiently using adaptive optics on a single-dish telescope (through imaging or sparse aperture masking).

Each point in the grid is fitted with the following binary model, representing a spatially resolved primary star with a point source component:

$$ \tilde{V}(u, v) = \frac{V_{\star}(u, v) + G(\zeta) f \tilde{V}_{c}(u, v)}{1 + f}, \quad (1) $$

with,

$$ V_{\star}(u, v) = \frac{2J_{1}(x)}{x}, \quad (2) $$

$$ \tilde{V}_{c}(u, v) = \exp[-2i\pi(u\Delta\alpha + v\Delta\delta)/\lambda], \quad (3) $$

$$ G(\zeta) = \left| \frac{\sin \zeta}{\zeta} \right| \quad \text{with} \quad \zeta = \frac{\pi(u\Delta\alpha + v\Delta\delta)}{R\lambda}, \quad (4) $$

where $J_1(x)$ is the first-order Bessel function, $x = \pi\theta_{UD}\sqrt{u^2 + v^2}/\lambda$, $(u, v)$ the spatial frequencies, $\theta_{UD}$ the uniform disk angular diameter of the primary star, $\lambda$ the wavelength, $f$ the flux ratio between the companion and the primary star, $(\Delta\alpha, \Delta\delta)$ the relative position of the component with respect to the primary, and $R = \lambda/\Delta\lambda$ the spectral resolution. The function $G(\zeta)$ is a corrective term to overcome the effect of bandwidth smearing (Lachaume & Berger 2013). The fitted parameters are $\theta_{UD}$, $f$, $\Delta\alpha$, and $\Delta\delta$. Although the bandwidth smearing can be fit in CANDID, we kept it fixed for this work because we noticed that it cannot be constrained by these data.

The interferometric observables, i.e., the squared visibility $V^2$, the closure phase CP, and the bispectrum amplitude $B_{amp}$, are then estimated from the squared modulus and the bispectrum in closed triangles as follows:

$$ V^2 = |\tilde{V}|^2 \quad \text{and} \quad \tilde{B} = \tilde{V}_{12}\tilde{V}_{23}\tilde{V}_{31}^*, \quad (5) $$

which provides the bispectrum amplitude and the closure phase from the definition $\tilde{B} = B_{amp}e^{-iCP}$:

$$ B_{amp} = |\tilde{B}| \quad \text{and} \quad CP = \arg(\tilde{B}). \quad (6) $$

The data are then fitted simultaneously using a Levenberg-Marquardt least-squares minimization algorithm with

$$ \chi^2 = \sum (\text{CP}_o - \text{CP}_m)^2 / \sigma_{\text{CP}}^2 + \sum (B_{\text{amp},o} - B_{\text{amp},m})^2 / \sigma_{B_{\text{amp}}}^2 + \sum (V_o^2 - V_m^2)^2 / \sigma_{V^2}^2, \quad (7) $$

where the indexes o and m denote the data and the model, respectively. CANDID can fit all observables or just one, depending on the data available and the user. We then divided by the number of degrees of freedom to obtain a map of the chi-square minima, and interpolated these minima on a regular grid.

The grid resolution is critical as it depends on the distance explored from the initial and final positions. This is why we implemented an estimate of the optimum grid resolution in our code. It is worth mentioning that a denser starting grid is not necessarily better as the fit would not be improved, and we would lose in computation time (CANDID performs a 30×30 grid with 1 mas steps in 50 s using six cores and a data set with 879 degrees of freedom). To make the grid search faster, CANDID was developed for parallel processing on multicore machines.

Once the map of the minima is computed, we can check the variation in the whole map, however this does not tell us about the statistical significance: specifically, whether the most probable location is detected at $1\sigma$ or more. Although it might not be important for low-contrast companions for which the variation is large enough, this is critical for components with a flux ratio <5% with a variation of the minima of a few percent. The number of dof is an important parameter in that context. Assuming that the data follow Gaussian statistics, we implemented in CANDID an estimate of the number of sigma for each point in the grid in order to obtain a $n\sigma$ detection map. The formalism we used is based on the probability $P$ (or confidence interval) with $\nu$ degrees of freedom, as already employed by Absil et al. (2011). The number of sigma demonstrates how our binary model is significant compared to a uniform disk model (i.e., a single star). We used the following formula for the probability:

$$ P(\Delta\alpha, \Delta\delta) = 1 - \text{CDF}_{\nu} \left( \frac{\nu\chi_{\text{UD}}^2}{\chi_{\text{r,bin}}^2(\Delta\alpha, \Delta\delta)} \right), \quad (8) $$

where $\chi_{\text{bin}}^2$ and $\chi_{\text{UD}}^2$ are the minimum chi-square for the binary and the uniform disk model (i.e., fitting Eqs. (1) and (2), respectively), and CDF denotes the $\chi^2$ cumulative probability distribution function with $\nu$ degrees of freedom. We then convert the probability into the number of sigmas, $n\sigma$ (e.g., 99.73% = $3\sigma$, 99.99% = $4\sigma$, ...). To avoid big float numbers, we limited the maximum value to 50$\sigma$.

This formalism therefore provides a $\chi^2$ map to find the most probable location of a companion, if any, and a $n\sigma$ map giving the detection level at each point in the grid.

## 2.2. Previously published detections

To validate our method, we computed the maps of the $\chi^2$ minima and estimated the detection level for our previously detected companions around the Cepheids V1334 Cyg and AX Cir (Gallenne et al. 2013b, 2014, observed with the instruments CHARA/MIRC and VLTI/PIONIER, respectively). The journal of these previous observations are reported in Table 1. The maps are presented in Figs. 1 and 2, for which all observables were fitted, except for PIONIER where only CP and $V^2$ were used because there is no good estimator of $B_{amp}$ so far. The central part has been hidden to improve the clarity of the intensity map level, which can be biased by the primary star. The companion orbiting V1334 Cyg is detected at more than $25\sigma$, and at more than $13\sigma$ for AX Cir. We summarized the detection levels for these two stars in Table 2, including fitting only the closure phase signal. The CP is more sensitive to faint off-axis companions and is also less affected by instrumental and atmospheric perturbations than the other observables (i.e., $V^2$ and $B_{amp}$). Fitting all of the observables can improve the detection level because we add more information, but it can also affect the results, depending on the magnitude of the biases. We notice from Table 2 that including $V^2$ and $B_{amp}$ degrades the detection level for V1334 Cyg, although it is still significant, and including these observables improves the detection for AX Cir. A possible explanation that addresses why the detection level decreases sometimes when we add the $V^2$ is that we add correlated noise, which is not consistent with our hypothesis of uncorrelated noise. When only the CP is used, the angular diameter of the primary is first determined by fitting only a uniform disk model to the square visibilities, and then kept fixed during the grid search.
<table><caption>**Table 1.** Journal of the observations.</caption><thead><tr><th>UT</th><th>MJD</th><th>Star</th><th>N<sub>bracket</sub></th><th>Configuration</th><th>N<sub>V2</sub></th><th>N<sub>CP</sub></th><th>Calibrators</th></tr></thead><tbody><tr><td>2012 Jul. 26</td><td>56 135.449</td><td>V1334 Cyg</td><td>2</td><td>S1-S2-E1-E2-W2</td><td>48</td><td>42</td><td>HD 200577, HD 214200</td></tr><tr><td>2012 Sep. 30</td><td>56 201.221</td><td>V1334 Cyg</td><td>3</td><td>S1-S2-E1-E2-W1-W2</td><td>62</td><td>68</td><td>HD 185395, HD 199956, HD 218470, HD 207978</td></tr><tr><td>2013 Jul. 14</td><td>56 487.983</td><td>AX Cir</td><td>6</td><td>D0-H0-G1-I1</td><td>300</td><td>200</td><td>HD 133869, HD 129462</td></tr><tr><td>2012 Jul. 26</td><td>56 134.354</td><td>T Vul</td><td>2</td><td>S1-S2-E1-E2-W1-W2</td><td>199</td><td>260</td><td>HD 192518, HD 205852</td></tr><tr><td>2012 Sept. 30</td><td>56 200.228</td><td>T Vul</td><td>2</td><td>S1-S2-E1-E2-W1-W2</td><td>175</td><td>210</td><td>HD 189849, HD 198692, HD 207978</td></tr><tr><td>2012 Sept. 30</td><td>56 200.434</td><td>SU Cas</td><td>4</td><td>S1-S2-E1-E2-W1</td><td>259</td><td>255</td><td>HD 12216, HD 19267, HD 34200</td></tr><tr><td>2012 Oct. 01</td><td>56 201.428</td><td>AW Per</td><td>1</td><td>S1-S2-E1-E2-W1-W2</td><td>210</td><td>280</td><td>HD 19845, HD 30825, HD 35940</td></tr><tr><td>2012 Oct. 01</td><td>56 201.507</td><td>RT Aur</td><td>1</td><td>S1-S2-E1-E2-W1-W2</td><td>105</td><td>140</td><td>HD 48682</td></tr></tbody></table>

**Notes.** $N_{bracket}$: number of data blocks. $N_{V2}$ and $N_{CP}$: number of squared visibilities and closure phase. Adopted calibrator diameters: HD 200577 = 0.758 ± 0.052 mas, HD 214200 = 0.790 ± 0.050 mas, HD 185395 = 0.726 ± 0.014 mas, HD 199956 = 0.603 ± 0.043 mas, HD 218470 = 0.477 ± 0.033 mas, HD 207978 = 0.571 ± 0.040 mas, HD 133869 = 1.043 ± 0.015 mas, HD 129462 = 0.857 ± 0.061 mas, HD 192518 = 0.418 ± 0.029 mas, HD 205852 = 0.461 ± 0.032 mas, HD 189849 = 0.510 ± 0.036 mas, HD 198692 = 0.660 ± 0.056 mas, HD 207978 = 0.571 ± 0.040 mas, HD 12216 = 0.467 ± 0.033 mas, HD 19267 = 0.586 ± 0.042 mas, HD 34200 = 0.652 ± 0.046 mas, HD 48682 = 0.616 ± 0.043 mas, HD 19845 = 0.788 ± 0.056 mas, HD 30825 = 0.564 ± 0.040 mas, HD 35940 = 0.615 ± 0.044 mas.

<figure><img src="image_3.png" /><figcaption>**Fig. 1.** $\chi^2$ map of the local minima (left) and detection level map (right) of V1334 Cyg for the observations on 2012-10-01. The yellow lines represent the convergence from the starting points to the final fitted position. The maps were reinterpolated in a regular grid for clarity. The axis limit was chosen according to the location of the companion.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Same as Fig. 1, except for AX Cir for the observations on 2013-07-14.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>Fig. 3. Same as Fig. 1, except for RT Aur using only the closure phase.</figcaption></figure>

<table><caption>**Table 2.** Summary of the detection level for our previous published data and new detections.</caption><thead><tr><th>Star</th><th>CP only<br>(nσ)</th><th>ALL <sup>a</sup><br>(nσ)</th><th>d.o.f.</th><th>Instrument</th></tr></thead><tbody><tr><td>V1334 Cyg</td><td>42.1</td><td>25.4</td><td>1583</td><td>MIRC</td></tr><tr><td>AX Cir</td><td>13.6</td><td>19.8</td><td>1499</td><td>PIONIER</td></tr><tr><td>RT Aur</td><td>3.8</td><td>2.3</td><td>879</td><td>MIRC</td></tr><tr><td>AW Per</td><td>19.6</td><td>15.0</td><td>1757</td><td>MIRC</td></tr></tbody></table>

**Notes.** (a) All means $CP + V^2 + B_{amp}$, and $B_{amp}$ is only used with MIRC data.

The resulting fitted parameters (i.e., the astrometric position, flux ratio, and the angular diameter) are in good agreement with the values determined in Gallenne et al. (2013b, 2014).

## 2.3. New detections of the companions around the Cepheids RT Aur and AW Per

Our interferometric program on Galactic binary Cepheids, which started two years ago, is promising to directly measure the dynamical masses. We have obtained several observing nights in both hemispheres with the multitelescope combiners CHARA/MIRC and VLTI/PIONIER to detect the close companions of a few Cepheids (Gallenne et al. 2014, 2013b,a). Here we report new detections for the Cepheids AW Per and RT Aur.

The observations were performed in 2012 using the Michigan InfraRed Combiner (MIRC) installed at the CHARA Array (ten Brummelaar et al. 2005), located on Mount Wilson, California. The CHARA Array consists of six 1 m aperture telescopes in a Y-shaped configuration (two telescopes on each branch), oriented to the east (E1, E2), west (W1,W2) and south (S1, S2), providing good coverage of the (u, v) plane. The baselines range from 34 m to 331 m, providing high angular resolution down to 0.5 mas at H. The MIRC instrument (Monnier et al. 2004, 2010) is an image-plane combiner, which enables us to combine the light coming from all six telescopes in H or K. MIRC also offers three spectral resolutions ($R$ = 42, 150 and 400), which provide 15 visibility and 20 closure phase measurements across a range of spectral channels.

We observed the short-period Cepheids RT Aur (HD 45412, $P_{puls}$ = 3.73 d) and AW Per (HD 30282, $P_{puls}$ = 6.46 d) with all six telescopes. We used the H-band filter with the lowest spectral resolution in which the light is split into eight spectral channels. Table 1 lists the journal of our observations. We followed a standard observing procedure in which we monitored the interferometric transfer function by observing a calibrator before and/or after the Cepheids. The calibrators, listed in Table 1, were selected using the SearchCal<fnref n="2" /> software (Bonneau et al. 2006, 2011) provided by the Jean-Marie Mariotti Center<fnref n="3" />.

We reduced the data using the standard MIRC pipeline (Monnier et al. 2007), which consists of computing the squared visibilities and triple products for each baseline and spectral channel, and correcting for photon and readout noise. Squared visibilities are estimated using Fourier transforms, while the triple products are evaluated from the amplitudes and phases between three baselines forming a closed triangle.

We used CANDID to search for a component within ±50 mas. For RT Aur, we might have detected a companion at $3.8\sigma$ using only the closure phases, while we only have a $2.3\sigma$ detection using all observables. The possible companion is detected at $\rho$ = 2.1 mas and PA = -136° with a flux ratio $f$ = 0.22%. The component orbiting AW Per is detected at >15$\sigma$ at $\rho$ = 32 mas and PA = 67° with a flux ratio $f$ = 1.22%. The grids of fit and detection level maps are shown in Fig 3 and 4, and the final fitted parameters are listed in Table 3. For the uncertainties, we used the conservative formalism of Boffin et al. (2014) for all the fitted parameters, i.e.,

$$ \sigma_X^2 = N_{sp}\sigma_{stat}^2 + 0.0001X^2, \qquad (9) $$

where $N_{sp}$ is the number of spectral channels and $X$ denotes the fitted parameters (i.e., $\Delta x, \Delta y, f$ and $\theta_{UD}$). The first term takes into account that the spectral channels are almost perfectly correlated, and the second term comes from the fact that the wavelength calibration is only precise at a 1% level. The parameter $\sigma_{stat}$ is the statistical error from the bootstrapping technique (bootstrap on the calibrated data with replacement) using 10 000 bootstrap samples (also included in CANDID). We then took from the distributions the median and the maximum

<footnote n="2">² Available at http://www.jmmc.fr/searchcal</footnote>

<footnote n="3">³ http://www.jmmc.fr</footnote>
<figure><img src="image_6.png" /><figcaption>Fig. 4. Same as Fig. 1, except for AW Per.</figcaption></figure>

<table><caption>**Table 3.** Final best-fit parameters.</caption><thead><tr><th></th><th>RT Aur<br>φ = 0.32</th><th>AW Per<br>φ = 0.52</th><th>SU Cas<br>φ = 0.77</th><th>T Vul<br>φ = 0.27 / φ = 0.12</th></tr></thead><tbody><tr><td>θ<sub>UD</sub> (mas)</td><td>0.699 ± 0.011</td><td>0.627 ± 0.018</td><td>0.609 ± 0.043</td><td>0.608 ± 0.013/0.635 ± 0.018</td></tr><tr><td>f (%)</td><td>0.22 ± 0.11</td><td>1.22 ± 0.30</td><td>—</td><td>—</td></tr><tr><td>Δα (mas)</td><td>-1.458 ± 0.238</td><td>29.624 ± 0.305</td><td>—</td><td>—</td></tr><tr><td>Δδ (mas)</td><td>-1.506 ± 0.224</td><td>12.523 ± 0.147</td><td>—</td><td>—</td></tr></tbody></table>

**Notes.** $θ_{UD}$: uniform disk angular diameter, respectively. f, Δx, Δy: flux ratio and position of the companion.

value between the 16% and 84% percentiles as the uncertainty (although the distributions were roughly symmetrical).

There are additional significant peaks in the RT Aur maps, with a detection level >3σ, however, these are spuriously produced by the (u, v) coverage and the presence of the companion. It is worth mentioning that RT Aur was only observed for one hour (one sequence), and we need more data to confirm the presence of the companion. A more complete discussion about the detected companions is presented in Sect. 4.

## 2.4. Undetected companions

We also observed with MIRC the Cepheids SU Cas (HD 17463, *P* = 1.95 day) and T Vul (HD 198726, *P* = 4.44 day). We used five and six telescope configurations, using the same instrument setup and calibration procedures as explained in the previous section. Table 1 lists the journal of these observations.

We did not detect any significant companions around these stars, either using all of the observables or only the CP, i.e., no more than 1.9σ and 2.9σ for SU Cas and T Vul, respectively. However, we were able to measure their angular diameters, which are listed in Table 3. A more detailed discussion of these binary systems is presented in Sect. 4.

## 3. Detection limit of high-contrast binaries

It is important to be able to check the dynamic range that can be reached with a given set of data and any interferometric combiner. Absil et al. (2011) already presented a method to

set detection limits for the VLTI/PIONIER instrument, but it has some shortcomings, that is why we propose a more robust formalism.

**Absil's method:** their method is based on comparing a uniform disk model with a binary model for each position (Δα, Δδ) in the grid. They then checked whether the probability of the binary model is consistent with the data using,

$$P(\Delta\alpha, \Delta\delta) = 1 - \text{CDF}_{\nu} \left( \frac{\nu\chi_{r,\text{bin}}^2 (\Delta\alpha, \Delta\delta)}{\chi_{\text{UD}}^2} \right). \quad (10)$$

We notice that this equation has a different ratio in the CDF than Eq. (8). This is because we assumed that the binary model is the true model, while Absil et al. (2011) assumed the uniform disk as the true model. In theory, both equations should lead to the same results, however, as we see in subsequent sections, their method is more sensitive to biased data and can sometimes lead to under- or overestimated detection limits.

**Our method:** we suggest an alternative method, which is based on the injection of a 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, and we can use Eq. (8) to obtain the probability of the binary model to be the true model. We introduced this method because we think it is more robust, as we demonstrate in the next section.

In this section, we first introduce approximate formulae for high-contrast companions. We then explain how we inject an additional component into the data and derive the detection limits.
<figure><img src="image_7.png" /><figcaption>Fig. 5. Same as Fig. 3, except we analytically removed the companion.</figcaption></figure>

### 3.1. High-contrast approximation

The complex visibility for a binary system composed of a resolved primary star and an unresolved component is given by Eq. (1). The squared visibility is given by,

$$V^2 = |\tilde{V} \cdot \tilde{V}^*| = \frac{1}{1 + f^2} |(V_* + G(\zeta)e^{i\varphi})(V_* + G(\zeta)e^{-i\varphi})| \\ = \left| \frac{V_*^2 + 2G(\zeta)fV_* \cos \varphi + G(\zeta)^2 f^2}{(1 + f)^2} \right|, \quad (11)$$

where we kept the $f^2$ term to avoid having an (small) offset (see next section). For a high-contrast companion, i.e., for $f \ll 1$, we can approximate the bispectrum at the first order in $f$ as,

$$\tilde{B} = \frac{(V_1 + G(\zeta_1)e^{i\varphi_1})(V_2 + G(\zeta_2)e^{i\varphi_2})(V_3 + G(\zeta_3)e^{-i\varphi_3})}{(1 + f)^3} \\ \sim \frac{\tilde{B}_*}{(1 + f)^3} \left[ 1 + f \left( \frac{G(\zeta_1)e^{-i\varphi_1}}{\tilde{V}_{1*}} + \frac{G(\zeta_2)e^{-i\varphi_2}}{\tilde{V}_{2*}} + \frac{G(\zeta_3)e^{i\varphi_3}}{\tilde{V}_{3*}^*} \right) \right] \\ \sim \frac{\tilde{B}_*}{(1 + f)^3} Z, \quad (12)$$

where $\tilde{B}_* = \tilde{V}_{1*} \tilde{V}_{2*} \tilde{V}_{3*}^*$.

The bispectrum amplitude and the closure phase are then estimated as,

$$B_{\text{amp}} = \frac{1}{(1 + f)^3} |\tilde{B}_{\star}| |Z| \quad (13)$$

$$CP = \arg(\tilde{B}_{\star}) - \arg(Z). \quad (14)$$

### 3.2. Adding/removing a component

If no companion is detected in the data, we can assume that the measured values only represent the primary star, i.e., a uniform disk, and we can substitute the index “★” in the previous equations by the index “obs” (note that the function $Z$ also depends on the primary star visibility). It is now simple, from Eqs. (11), (13) and (14), to inject a companion to

the observed data, which corresponds to the following equations for the following new observables:

$$V^2 = \frac{V_{\text{obs}}^2 + 2G(\zeta)fV_{\text{obs}} \cos \varphi + G(\zeta)^2 f^2}{(1 + f)^2}, \quad (15)$$

$$B_{\text{amp}} = \frac{1}{(1 + f)^3} B_{\text{amp,obs}} |Z|, \quad (16)$$

$$CP = CP_{\text{obs}} - \arg(Z). \quad (17)$$

The original “oifits” files provide all the necessary spatial and spectral information to reconstruct the individual phases and visibilities. We also take the coherence loss effect caused by the spectral smearing of the companion into account.

Inversely, if a component is detected, we can use this formula to analytically remove the companion (with a negative flux ratio), and then check for another possible fainter component or estimate the detection limit to rule out any other companions. This step is critical to obtain unbiased detection limits. As an example in Fig. 5, we analytically removed the detected component orbiting RT Aur, and notice that there is no significant detection of a third component and that the other $\chi^2$ minima are not additional companions.

To check our approximation, we created a model of a single star with a uniform disk of 1 mas (3h of PIONIER observations with a point every 30 min), for which we added a companion at a position $\Delta\alpha = \Delta\delta = 50$ mas with a flux ratio $f = 5\%$. We then compare this model with a true binary model using the same parameters (Eq. (1)). The difference between the approximation and the true model is shown in Fig. 6. We found for the amplitude of the bispectrum a standard deviation of the relative error <0.3%, and a difference <0.005° for the closure phase, which is lower than the achievable interferometric accuracy. Closer companions with the same or higher contrast (i.e., <5%) lead to smaller errors.

We kept the $f^2$ term in the squared visibility, otherwise it gives an offset in the relative error of about 0.3% for $f = 5\%$, and decreases with decreasing flux ratios. Although this value is negligible compared to the current possible data accuracy ($\Delta V^2/V^2 = \Delta B_{\text{amp}}/B_{\text{amp}} \sim 2\%$, $\Delta CP \sim 0.5^\circ$), the offset is larger for lower contrasts (~1.2% with $f = 10\%$). As no approximation is made for $V^2$, the formula is valid for any flux ratio. We therefore also checked the bispectrum approximation
<figure><img src="image_8.png" /><figcaption>**Fig. 6.** Difference between the true binary model and the approximation model. For the bispectrum, the abscissa is the maximum of the three spatial frequencies.</figcaption></figure>

<table><caption>**Table 4.** Validity domain of our approximation.</caption><thead><tr><th>f<br>(%)</th><th>ΔCP<br>(°)</th><th>ΔB<sub>amp</sub>/B<sub>amp</sub><br>(%)</th></tr></thead><tbody><tr><td>5</td><td>&lt;0.01</td><td>&lt;0.5</td></tr><tr><td>10</td><td>&lt;0.02</td><td>&lt;1.9</td></tr><tr><td>50</td><td>&lt;0.25</td><td>&lt;67</td></tr><tr><td>75</td><td>&lt;0.7</td><td>&gt;10</td></tr><tr><td>100</td><td>&lt;3.5</td><td>&gt;100</td></tr></tbody></table>

for brighter components, and found that we can use it up to a flux ratio of 50% only using CP and $V^2$. We listed in Table 4 the relative error for different flux ratios and various positions ($1 < \rho < 50$ mas and $0 < \text{PA} < 2\pi$). We see that CANDID can also be used for low-contrast companions.

### 3.3. Estimating the detection limit

We implemented in CANDID a tool to estimate the detection limit for a given set of interferometric data. The method is to compare the $\chi^2$ obtained for a model without a companion and the $\chi^2$ obtained for a model with an injected companion. The method does not work, a priori, with an already detected companion, first because the previous approximation might not be valid and also because the resulting detection limit would be biased by a systematic offset related to the flux ratio between the components. This means that any detected component has to be analytically removed first.

The method is based on an $N \times N$ grid with a range of ±50 mas, for which the minimum $N$ depends on the optimum resolution estimated from the $\chi^2$ map (see Sect. 2). At each point in the grid, we inject a companion with various flux ratios and we compute the $\chi^2$. As we know that the true model is the binary model (because we injected a companion), we used our previous equation, Eq. (8), to estimate the number of sigmas for each flux ratio. We then interpolated the flux ratio values at $3\sigma$, which we set as the significance level. This means that lower flux ratios are not detected significantly. Doing this for all points in the grid, we then have a $3\sigma$ detection limit map for the flux ratio. To have a quantitative estimate of the sensitivity limit with respect to the separation, we estimated a radial profile, $f_{3\sigma}(r)$, using the 90% completeness level (i.e., 90% of all possible positions) from

<figure><img src="image_9.png" /><figcaption>**Fig. 7.** Comparison between the CANDID detection limit tool and the formalism of Absil et al. (2011) for a uniform disk model with Gaussian noise.</figcaption></figure>

the cumulated histogram in rings for all azimuths. This tool also includes parallel processing to make the calculation faster.

In theory, for uncorrelated data with Gaussian statistics, our method should lead to the same results as Absil et al. (2011) in terms of the detection limit. However, real data are often biased by different sources (atmospheric turbulence, mechanical vibrations, ...). We performed two tests to compare both methods, one including uncorrelated Gaussian noise and another with a noise model, i.e., correlated non-Gaussian noise. We used all observables and three data sets for each test. The Aspro2 software<fnref n="4" /> were used to create these synthetic data sets, however, it does not have implemented the bandwidth smearing effect. We therefore did not take it into account in CANDID (this would not change the conclusion of our test).

**First test:** we created the first “ideal” synthetic data set (i.e., without noise) representing a uniform disk of 1 mas (3 h observation with three spectral channels with the PIONIER instrument). We then added uncorrelated Gaussian noise and estimated the detection limit using the Absil’s method and our formalism. As expected, we see in Fig. 7 that both methods provide the same results. We then created a second data set by adding a component at $\Delta\alpha = \Delta\delta = 5$ mas with $f = 2\%$ as a source of bias to the previous model. The companion creates a departure from the ideal measurements (like noise) but in such a way that the observations at different spatial frequencies have a correlated departure (hence, correlated noise). We then estimated the detection limits for this second data set, and although the trend between the two methods are similar, there is an offset of ~2% because of the presence of the faint companion. Finally, we analytically removed the companion to get our third data set. Comparing the three data sets, i.e., estimating the total variation of the detection limit for each radius, we noticed that both methods vary in a similar way. This was expected as all of the noise sources inserted have Gaussian statistics. However, it is worth mentioning that the scatter tends to be a bit larger for the Absil’s method for increasing flux ratios and wider separations of the companion (due to the correlated noise introduced by the component).

**Second test:** we created two new data sets similar to the first test, i.e., one for a uniform disk of 1 mas and another for a binary system with $\Delta\alpha = \Delta\delta = 5$ mas and $f = 2\%$ (with average atmospheric conditions, which is an option for the noise model). The only difference is that the noise is no longer Gaussian; it is represented by a more complicated model that includes the instrument

<footnote n="4">⁴ Available at http://www.jmmc.fr/aspro_page.htm</footnote>
<figure><img src="image_10.png" /><figcaption>**Fig. 8.** Comparison between the CANDID detection limit tool and the formalism of Absil et al. (2011) for a uniform disk model with a noise model, biased by the presence of a faint companion.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 9.** Total variation (minimum and maximum values) between the three synthetic data sets with a noise model for both methods: a single star, a single star + an unresolved component as source of bias, and the previous model with the component analytically removed. The variation has been normalized to unity.</figcaption></figure>

response, atmospheric turbulence, photon, and detector noise<fnref n="5" />. As shown in Fig. 8, the absil's method is more sensitive to the presence of the companion, which in this case gives a higher sensitivity limit. The third data set is created by removing the component using our approximation. We then compared the variation of the detection limit between all of the data sets, as described previously, i.e., the minimum and maximum value between the three synthetic data sets at each separation. Figure 9 shows that the scatter from the Absil's method is larger. We performed additional tests with lower flux ratios and at different positions, and the scatter from the Absil's method is larger than our formalism most of the time. We therefore conclude that our method is more robust to biased data when estimating detection limits.

In our tests, the limits derived from the injection method are lower than when we use Absil's method, but this is not always the case for real data. The injection method occasionally results in higher sensitivity limits, depending on the magnitude of the biases. We find that Absil's method may under- or overestimate the detection limits depending on the data set. However, in some cases when the data are not affected by significant biases, both methods give similar results.

<figure><img src="image_12.png" /><figcaption>**Fig. 10.** Flux ratio detection limit at $3\sigma$ for a second companion around V1334 Cyg (for 2012 Oct. 01).</figcaption></figure>

## 4. Detection limit of our binary Cepheids

In this section, we set the detection limit for the sample of Cepheids previously presented. For V1334 Cyg, AX Cir, RT Aur, and AW Per, where a companion is detected, we removed analytically the companion first. Although both methods are estimated by CANDID, all of the given detection limits are derived from our injection method. We present two detection limits for each star, one using all of the observables, and one using only the closure phases. We listed three different values, the average for $r < 50$ mas, $r < 25$ mas and $r > 25$ mas, which can be relevant when the limit increases with $r$. All of the final detection limits are listed in Table 5, where the values are conservative as they correspond to the mean plus the standard deviation for the given radius range.

From an evolutionary timescale point of view, most of the companions should be stars close to the main sequence. We therefore set upper limits for the spectral type of the companion assuming it is on the main sequence, based on their H-band luminosities.

V1334 Cyg: this binary system contains a visual and a spectroscopic component. While the visual companion is $>150$ mas, the close component was spatially resolved using interferometry by Gallenne et al. (2013b), and has a flux ratio of ~3.1% in H (more details on this companion are presented in Gallenne et al. 2013b). There is no evidence of a third component so far. After removing the close companion analytically, we estimated the dynamic range as explained in the previous section. Between the two epochs, we reached a maximum average sensitivity limit of $f_{3\sigma} = 1.54\%$. The limit is lower for the first epoch because the atmospheric conditions were better than the second epoch. The contrast upper limit at $3\sigma$ is shown in Fig. 10, with the average values for the two epochs in Table 5. This converts to magnitude difference of $\Delta m_H > -2.5 \log f_{3\sigma} = 4.5$ mag. Using the distance $d = 683$ pc from the K-band P–L relation for first overtone (FO) pulsators (Bono et al. 2002) and the average magnitudes of the Cepheid $H = 4.66$ mag (Cutri et al. 2003), we can exclude the presence of additional companions with a spectral type earlier than B7V stars.

AX Cir: the spectroscopic companion was first detected by Lloyd Evans (1982), and was spatially resolved for the first time with LBI by Gallenne et al. (2014, with a more detailed discussion about this companion). The average detection limits are listed in Table 5. We reached a dynamic range of $f_{3\sigma} = 0.36\%$, corresponding to $\Delta m_H > 6.1$ mag if another companion

<footnote n="5">⁵ Details are explained in http://www.jmmc.fr/doc/approved/JMMC-MEM-2800-0001.pdf</footnote>
<table><caption>**Table 5.** 3σ average detection limits of the flux ratio.</caption><thead><tr><th rowspan="2"></th><th colspan="2">All observables<br>r &lt; 50 mas</th><th colspan="2">Only CP<br>r &lt; 50 mas</th><th rowspan="2">Sp. type<br>upper limit</th><th rowspan="2">Instrument</th></tr><tr><th>r &lt; 25 mas</th><th>r &gt; 25 mas</th><th>r &lt; 25 mas</th><th>r &gt; 25 mas</th></tr></thead><tbody><tr><td>V1334 Cyg<br>(2012 Jul. 27)</td><td>1.71%<br>1.45%</td><td>1.72%<br>1.72%</td><td>1.54%<br>1.19%</td><td>1.57%<br>1.57%</td><td>B7V</td><td>MIRC</td></tr><tr><td>V1334 Cyg<br>(2012 Oct. 01)</td><td>2.59%<br>2.15%</td><td>2.60%<br>2.60%</td><td>2.44%<br>1.86%</td><td>2.42%<br>2.42%</td><td>-</td><td>-</td></tr><tr><td>AX Cir<br>(2013 Jul. 14)</td><td>0.36%<br>0.40%</td><td>0.34%<br>0.34%</td><td>0.73%<br>0.80%</td><td>0.69%<br>0.69%</td><td>A5V</td><td>PIONIER</td></tr><tr><td>RT Aur</td><td>0.64%<br>0.44%</td><td>0.64%<br>0.64%</td><td>0.47%<br>0.33%</td><td>0.47%<br>0.47%</td><td>F0V</td><td>MIRC</td></tr><tr><td>AW Per</td><td>0.62%<br>0.54%</td><td>0.62%<br>0.62%</td><td>0.72%<br>0.52%</td><td>0.73%<br>0.73%</td><td>B9V</td><td>MIRC</td></tr><tr><td>SU Cas</td><td>1.43%<br>1.06%</td><td>1.44%<br>1.44%</td><td>1.37%<br>0.99%</td><td>1.37%<br>1.37%</td><td>A0V</td><td>MIRC</td></tr><tr><td>T Vul<br>(2012 Jul. 26)</td><td>1.04%<br>0.84%</td><td>1.05%<br>1.05%</td><td>1.13%<br>0.81%</td><td>1.14%<br>1.14%</td><td>B9V</td><td>MIRC</td></tr><tr><td>T Vul<br>(2012 Sep. 30)</td><td>1.21%<br>0.89%</td><td>1.21%<br>1.21%</td><td>1.29%<br>0.92%</td><td>1.30%<br>1.30%</td><td>-</td><td>-</td></tr></tbody></table>

is present. With $d = 500$ pc (from the K-band P–L relation of Storm et al. 2011) and $H = 4.66$ mag for the Cepheid (Cutri et al. 2003), we can rule out any other component with a spectral type earlier than A5V star.

**RT Aur:** the binary nature of this short-period Cepheid is still uncertain, however we might have detected it for the first time. Some authors suggested the presence of an early-type companion (Janot-Pacheco 1976; Balona 1977), while others did not see evidence of an additional component (Harris 1981; Gieren 1985). Leonard & Turner (1986) summarized various studies from that time about the possible companion orbiting RT Aur. Based on spectra from the International Ultraviolet Explorer satellite (IUE), Evans (1992a) did not report any detection and showed that any main-sequence secondary has to be cooler than an A4 star. Radial velocity measurements do not show any variations from orbital motion. Recently, Turner et al. (2007) studied long-term photometric light curves and reported a sinusoidal trend consistent with a light travel time effect in the binary system. The companion we might have resolved has a flux ratio of $0.21 \pm 0.12\%$, i.e., $\Delta m_H = 6.7 \pm 0.6$ mag, leading to $m_H = 10.6 \pm 1.0$ mag (using for the Cepheid at our given phase $m_H = 3.94 \pm 0.01$ mag from Monson & Pierce 2011). Using the K-band P–L relation to get $d = 428$ pc (Storm et al. 2011), we estimate its spectral type to be later than an F1 star. This is compatible with the A4V star upper limit determined by Evans (1992a) from the International Ultraviolet Explorer (IUE) spectra. However, we need more observations to confirm the existence of the companion as we are at the sensitivity limit of this data set.

Removing this possible companion from our interferometric data, we estimated the detection limit. We reached a minimum dynamic range of 0.47%, corresponding to $\Delta m_H > 5.8$ mag. This is a 90% completeness azimuthal value, which is a conservative value, and the component at a given position with a slightly higher contrast might be detected (as in this case). This limit allows us to exclude the presence of an additional companion with a spectral type earlier than F0V.

**AW Per:** this Cepheid is a spectroscopic binary with an orbital period of ~40 yr. First discovered by Miller & Preston (1964), it took several years to derive the first orbit from radial velocity

data (Evans 1989). It is likely that this companion is itself a binary, as argued by Evans (1989), Evans et al. (2000), because the magnitude difference between the Cepheid and its companion is not consistent with equal masses and predictions from evolutionary tracks. Unfortunately, we do not have enough angular resolution and sensitivity to detect this third companion. The properties of the brightest companion were studied based on IUE spectra by Evans (1994, 1995), who found its spectral type to be B8.3V. This is in agreement with our detection with a flux ratio $f = 1.22 \pm 0.30\%$, i.e., a spectral type in the range B6-B9V (using $d = 853$ pc from a P–L relation and $K = 4.63$ mag for the Cepheid). We estimated for the companion $m_H = 9.6 \pm 0.3$ mag (using for the Cepheid at our given phase $m_H = 4.84 \pm 0.01$ mag from Monson & Pierce 2011). Massa & Evans (2008) also determined the angular separation of the component for another epoch, which allows us later, with more astrometric points from interferometry, to estimate all the orbital elements, including the inclination and the semi-major axis (see also Gallenne et al. 2013a).

We derived a maximum sensitivity limit at $3\sigma$ of 0.62%, i.e., $\Delta m_H > 5.5$ mag. We can therefore exclude any additional components with a spectral type earlier than B9V.

**SU Cas:** the binary nature of this Cepheid is still ambiguous. A component was first detected by Evans (1985) based on studying the CAII H and K lines, and then a spectral type of B9.5V was determined from IUE spectra (Evans 1991). Although the location of SU Cas in a two color diagram is consistent with the presence of a companion, the radial velocity measurements do not show convincing evidence. Szabados (1991) found four possible orbital periods from the observations available at that time. Later, (Gorynya et al. 1996) compiled more data and suggested SU Cas as a possible spectroscopic binary. They derived an orbital period of 408 days with an eccentricity of $e = 0.43$, but Groenewegen (2008) could not confirm the eccentricity with a larger data set and found a period of 407 days, fixing $e = 0$. Recently, Evans et al. (2013) reanalyzed special dates of Gorynya data, where the velocity difference was supposed to be the largest from the derived orbit, but they concluded that the orbital motion from radial velocity data do not show any convincing detection.
Using the distance *d* = 392 pc (from the K-band P–L relation of first overtone pulsators of Bono et al. 2002) and a Cepheid magnitude *H* = 4.27 mag (Cutri et al. 2003), a B9.5V companion should give a flux ratio ~1.8% in the *H* band, which should be detectable with MIRC. Taking the velocity amplitude of *K*₁ = 1 km s⁻¹ from Groenewegen (2008, while Gorynya et al. 1996 derived a velocity amplitude of 3 km s⁻¹), we estimated *a* sin *i* > 0.1 mas. Our nondetection could be explained if the companion was located at < 0.5 mas because it would not have been spatially resolved, but this type of close component would have an effect on the radial velocities (unless the orbit is face-on, but the probability of this kind of an orbit is low). The other possibility would be that this companion has a wider orbit. We therefore searched within a 100 mas radius range, but we did not find a significant detection. We then estimated the sensitivity limit for *r* < 100 mas, and found a maximum 3σ flux ratio of 1.65%. This means that if the B9.5V companion was within 100 mas, we would have detected it. Therefore, if this companion exists, it should have a wider orbit.

The average sensitivity limits for a given radius range are tabulated in Table 5. We reached a mean contrast of 1.37%, which means that $\Delta m_H > 4.6$ mag. This converts to an upper limit of the spectral type of an A0V star.

T Vul: as with SU Cus, the radial velocities of this Cepheid do not show any signature of an orbiting companion, while a hot A0.8V component was detected by Evans (1992b). In the literature, we found contradictory estimates of the orbital period from radial velocity measurements; for instance, Kovacs et al. (1990) found a long period modulation of 738 d, Szabados (1991) estimated *P* = 1745 d from a larger data set, while Bersier et al. (1994) did not find any orbital motion larger 0.55 km s⁻¹ using additional more accurate data. They showed that the long period of Szabados (1991) is not compatible and argued that the 738 d period might be an artifact of the time sampling because the observations were only made in autumn. Kiss & Vinkó (2000) reached the same conclusion with additional measurements showing no signature of orbital motion in the radial velocity curve.

From our interferometric observations, we did not detect any companion within 50 mas. The A0.8V component detected from IUE spectra would correspond to a flux ratio in *H* of ~0.7%. According to our estimated interferometric detection limit, listed in Table 5, this kind of a component is below our detection level. We reached an average sensitivity limit of 1.04%, and we can therefore exclude any other possible companion with a spectral type earlier than a B9V star.

## 5. Conclusion

We presented an overview of CANDID, a new tool to search for point-source companions and estimate the sensitivity level from interferometric observations using the squared visibilities, closure phases, and amplitude of the bispectrum, when available. CANDID allows us to:

- efficiently detect companions using a grid of fit and determine the detection level by giving the number of sigmas;

- set the detection limit for a companion in data where a companion has not been detected; and

- set the detection limit for a tertiary companion, in the case where a companion has been detected.

We used CANDID to investigate a sample of binary Cepheids. We first determined the detection level for our previous detections (Gallenne et al. 2013b, 2014) and showed that the

components were detected at >13σ for AX Cir and >25σ for V1334 Cyg. We also reported a new detection for AW Per, with a detection level >15σ; the companion is located at ρ = 32 mas and PA = 67°, with a flux ratio of *f* = 1.22%. The companion orbiting RT Aur might have been detected at 3.8σ, using only the closure phase signal, however, more observations are needed to confirm the presence of this component. Any additional companions were not detected significantly (i.e., with a detection level >3σ) around these stars. Likewise, no companions were detected around SU Cas and T Vul. From these interferometric data, we were able to set upper limits for the spectral types; we found no components with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V for V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively.

The fitting procedure of CANDID also allowed us to measure the uniform disk angular diameters of the new Cepheids observed. We found for RT Aur $\theta_{UD} = 0.699 \pm 0.011$ mas (at the pulsation phase $\phi = 0.32$), for AW Per $\theta_{UD} = 0.627 \pm 0.018$ mas (at $\phi = 0.52$), for SU Cas $\theta_{UD} = 0.609 \pm 0.043$ mas (at $\phi = 0.77$), and $\theta_{UD} = 0.608 \pm 0.013$ mas and $\theta_{UD} = 0.635 \pm 0.018$ mas for T Vul (at $\phi = 0.27$ and $\phi = 0.12$).

We demonstrated that the approximation we used to analytically inject a companion and estimate the detection limits is valid (i.e., error <0.5%) for contrasts *f* < 5% if we use all of the observables, and up to *f* < 50% using only the squared visibilities and the closure phases. This makes CANDID a useful tool for analyzing long-baseline interferometric observations of binary star systems.

Finally, this work demonstrates the capabilities of the MIRC and PIONIER instruments, which can reach a dynamic range of 1:200, depending on the angular distance of the companion and the (*u*, *v*) plane coverage.

*Acknowledgements.* The authors would like to thank the CHARA Array and Mount Wilson Observatory staff for their support. Research conducted at the CHARA Array is funded by the National Science Foundation through NSF grant AST-0908253, by Georgia State University, the W. M. Keck Foundation, the Packard Foundation, and the NASA Exoplanet Science Institute. The authors also thank all the people involved in the VLTI project. A.G. acknowledges support from FONDECYT grant 3130361. J.D.M. acknowledges funding from the NSF grants AST-0707927, AST-0807577, and AST 1108963. We also acknowledge the STSCI grant HST-GO13841.006-A and HST-GO-13454.07-A. W.G. and G.P. gratefully acknowledge financial support for this work from the BASAL Centro de Astrofísica y Tecnologías Afines (CATA) PFB-06/2007. Support from the Polish National Science Centre grant MAESTRO DEC-2012/06/A/ST9/00269 and the Polish Ministry of Science grant Ideas Plus (awarded to G. P.) is also acknowledge. W.G. also acknowledges financial support from the Millenium Institute of Astrophysics (MAS) of the Iniciativa Científica Milenio del Ministerio de Economía, Fomento y Turismo de Chile, project IC120009. A.G. also acknowledge computational resources from FONDECYT grant 1130521 and 1130721. We also acknowledge support from the ECOS/Conicyt grant C13U01. PIONIER is funded by the Université Joseph Fourier (UJF), the Institut de Planétologie et d'Astrophysique de Grenoble (IPAG), the Agence Nationale pour la Recherche (ANR-06-BLAN-0421, ANR-10-BLAN-0505, ANR-10-LABX56), and the Institut National des Science de l'Univers (INSU PNP and PNPS). The integrated optics beam combiner is the result of a collaboration between IPAG and CEA-LETI based on CNES R&T funding. This research received the support of PHASE, the high angular resolution partnership between ONERA, Observatoire de Paris, CNRS, and University Denis Diderot Paris 7. This work made use of the SIMBAD and VIZIER astrophysical database from CDS, Strasbourg, France and the bibliographic informations from the NASA Astrophysics Data System. This research has made use of the Jean-Marie Mariotti Center SearchCal and ASPRO services, co-developed by FIZEAU and LAOG/IPAG, and of CDS Astronomical Databases SIMBAD and VIZIER.

## References

Absil, O., Le Bouquin, J.-B., Berger, J.-P., et al. 2011, A&A, 535, A68
Balona, L. A. 1977, MNRAS, 178, 231
Bersier, D., Burki, G., Mayor, M., & Duquennoy, A. 1994, A&AS, 108, 25

Boffin, H. M. J., Hillen, M., Berger, J. P., et al. 2014, A&A, 564, A1

Bonneau, D., Clausse, J.-M., Delfosse, X., et al. 2006, A&A, 456, 789

Bonneau, D., Delfosse, X., Mourard, D., et al. 2011, A&A, 535, A53

Bono, G., Groenewegen, M. A. T., Marconi, M., & Caputo, F. 2002, ApJ, 574, L33

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246

Evans, N. 1985, in Cepheids: Theory and Observation, ed. B. F. Madore, IAU Colloq., 82, 79

Evans, N. R. 1989, AJ, 97, 1737

Evans, N. R. 1991, ApJ, 372, 597

Evans, N. R. 1992a, ApJ, 384, 220

Evans, N. R. 1992b, AJ, 104, 216

Evans, N. R. 1994, ApJ, 436, 273

Evans, N. R. 1995, ApJ, 445, 393

Evans, N. R., Vinko, J., & Wahlgren, G. M. 2000, AJ, 120, 407

Evans, N. R., Bond, H. E., Schaefer, G. H., et al. 2013, AJ, 146, 93

Gallenne, A., Kervella, P., Mérand, A., et al. 2013a, in EAS Pub. Ser., 64, eds. A. Pavlovski, K. Tkachenko, & G. Torres, 197

Gallenne, A., Monnier, J. D., Mérand, A., et al. 2013b, A&A, 552, A21

Gallenne, A., Mérand, A., Kervella, P., et al. 2014, A&A, 561, L3

Gieren, W. P. 1985, A&A, 148, 138

Gorynya, N. A., Rastorguev, A. S., & Samus, N. N. 1996, Astron. Lett., 22, 33

Groenewegen, M. A. T. 2008, A&A, 488, 25

Harris, H. C. 1981, AJ, 86, 707

Janot-Pacheco, E. 1976, A&AS, 25, 159

Kiss, L. L., & Vinkó, J. 2000, MNRAS, 314, 420

Kovacs, G., Kisvarsanyi, E. G., & Buchler, J. R. 1990, ApJ, 351, 606

Lachaume, R., & Berger, J.-P. 2013, MNRAS, 435, 2501

Le Bouquin, J.-B., & Absil, O. 2012, A&A, 541, A89

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

Leonard, P. J. T., & Turner, D. G. 1986, JRASC, 80, 240

Lloyd Evans, T. 1982, MNRAS, 199, 925

Massa, D., & Evans, N. R. 2008, MNRAS, 383, 139

Miller, J., & Preston, G. 1964, PASP, 76, 47

Monson, A. J., & Pierce, M. J. 2011, ApJS, 193, 12

Monnier, J. D., Berger, J.-P., Millan-Gabet, R., & ten Brummelaar, T. A. 2004, SPIE Conf. Ser., ed. W. A. Traub, 5491, 1370

Monnier, J. D., Zhao, M., Pedretti, E., et al. 2007, Science, 317, 342

Monnier, J. D., Anderson, M., Baron, F., et al. 2010, SPIE Conf. Ser., 7734 eds. W. C. Danchi, F. Delplancke, & J. K. Rajagopal, 77340G

Storm, J., Gieren, W., Fouqué, P., et al. 2011, A&A, 534, A94

Szabados, L. 1991, Commun. of the Konkoly Observatory Hungary, 96, 123

ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., et al. 2005, ApJ, 628, 453

Turner, D. G., Bryukhanov, I. S., Balyuk, I. I., et al. 2007, PASP, 119, 1247

Zhao, M., Monnier, J. D., Torres, G., et al. 2007, ApJ, 659, 626

Zurlo, A., Vigan, A., Mesa, D., et al. 2014, A&A, 572, A85