# Planet transit and stellar granulation detection with interferometry

A. Chiavassa, R. Ligi, Z. Magic, R. Collet, M. Asplund, D. Mourard

► To cite this version:

A. Chiavassa, R. Ligi, Z. Magic, R. Collet, M. Asplund, et al. Planet transit and stellar granulation detection with interferometry. Astronomy and Astrophysics - A&A, 2014, 567, pp.A115. 10.1051/0004-6361/201323207. [hal-03473879](https://hal.science/hal-03473879)

HAL Id: hal-03473879

https://hal.science/hal-03473879

Submitted on 17 Dec 2021

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# Planet transit and stellar granulation detection with interferometry

## Using the three-dimensional stellar atmosphere STAGGER-grid simulations

A. Chiavassa<fnref n="1" />, R. Ligi<fnref n="1" />, Z. Magic<fnref n="2" />, R. Collet<fnref n="3" />, M. Asplund<fnref n="3" />, and D. Mourard<fnref n="1" />

Received 6 December 2013 / Accepted 29 April 2014

## ABSTRACT

**Context.** Stellar activity and, in particular, convection-related surface structures, potentially cause bias in planet detection and characterisation. In the latter, interferometry can help disentangle the signal of the transiting planet.

**Aims.** We used realistic three-dimensional (3D) radiative hydrodynamical (RHD) simulations from the STAGGER-grid and synthetic images computed with the radiative transfer code OPTIM3D to provide interferometric observables to extract the signature of stellar granulation and transiting planets.

**Methods.** We computed intensity maps from RHD simulations and produced synthetic stellar disk images as a nearby observer would see, thereby accounting for the centre-to-limb variations. We did this for twelve interferometric instruments covering wavelengths ranging from optical to infrared. We chose an arbitrary date and arbitrary star with coordinates, and this ensures observability throughout the night. This optimisation of observability allows for a broad coverage of spatial frequencies. The stellar surface asymmetries in the brightness distribution mostly affect closure phases, because of either convection-related structures or a faint companion. We then computed closure phases for all images and compared the system star with a transiting planet and the star alone. We considered the impact of magnetic spots with the construction of a hypothetical starspot image and compared the resulting closure phases with the system star that has a transiting planet.

**Results.** We analysed the impact of convection at different wavelengths. All the simulation depart from the axisymmetric case (closure phases not equal to 0 or $\pm\pi$) at all wavelengths. The levels of asymmetry and inhomogeneity of stellar disk images reach high values with stronger effects from the 3rd visibility lobe on. We present two possible targets (Beta Com and Procyon) either in the visible or in the infrared and find that departures up to 16° can be detected on the 3rd lobe and higher. In particular, MIRC is the most appropriate instrument because it combines good UV coverage and long baselines. Moreover, we explored the impact of convection on interferometric planet signature for three prototypes of planets with sizes corresponding to one hot Jupiter, one hot Neptune, and a terrestrial planet. The signature of the transiting planet in the closure phase is mixed with the signal due to the convection-related surface structures, but it is possible to disentangle it at particular wavelengths (either in the infrared or in the optical) by comparing the closure phases of the star at difference phases of the planetary transit. It must be noted that starspots caused by the magnetic field may pollute the granulation and the transiting planet signals. However, it is possible to differentiate the transiting planet signal because the time scale of a planet crossing the stellar disk is much smaller than the typical rotational modulation of a star.

**Conclusions.** Detection and characterisation of planets must be based on a comprehensive knowledge of the host star, and this includes the detailed study of the stellar surface convection with interferometric techniques. In this context, RHD simulations are crucial for this aim. We emphasise that interferometric observations should be pushed at high spatial frequencies by accumulating observations on closure phases at short and long baselines.

**Key words.** stars: atmospheres – hydrodynamics – radiative transfer – techniques: interferometric – planetary systems

## 1. Introduction

Two very successful ways to find exoplanets orbiting stars are the transit and radial velocity methods. The transit happens when a planet passes between the exoplanet and its host star. The planet then blocks some of the starlight during the transit and creates a periodic dip in the brightness of the star. Observations taken during both the primary and secondary transits can be used to deduce the composition of the planet's atmosphere.

As the star moves in the small orbit resulting from the pull of the exoplanet, it will move towards the planet and then away as it completes an orbit. Regular periodic changes in the star's

radial velocity (i.e., the velocity of the star along the line of sight of an observer on Earth) depend on the planet's mass and the inclination of its orbit to our line of sight. Measurements of the Doppler-shifted spectra give a minimum value for the mass of the planet.

A potential complication to planet detection, however, may be caused by stellar surface inhomogeneities (due to the presence of stellar granulation, magnetic spots, dust, etc.) of the host star. In this article we investigate the particular problem of stellar granulation. It was first observed on the Sun by Herschel (1801), and today modern telescopes provide direct observations (e.g., Carlsson et al. 2004). However, the best observational evidence

<footnote n="1">¹ Laboratoire Lagrange, UMR 7293, CNRS, Observatoire de la Côte d'Azur, Université de Nice Sophia-Antipolis, 06108 Nice, France
e-mail: andrea.chiavassa@oca.eu</footnote>

<footnote n="2">² Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1, 85741 Garching, Germany</footnote>

<footnote n="3">³ Research School of Astronomy & Astrophysics, Australian National University, Cotter Road, Weston ACT 2611, Australia</footnote>
comes from unresolved spectral line in terms of widths, shapes, and strengths that, when combined with numerical models of convection, allow quite robust results to be extracted from the simulations (Nordlund et al. 2009; Asplund et al. 2000). For this purpose, large efforts have been made in recent decades to use theoretical modelling of stellar atmospheres to solve multidimensional radiative hydrodynamic equations in which convection emerges naturally. These simulations take surface inhomogeneities into account (e.g., granulation pattern) and velocity fields. The widths of spectral lines are heavily influenced by the amplitude of the convective velocity field, which overshoots into the stable layers of the photosphere where the lines are formed. This results in characteristic asymmetries of spectral lines, as well as in net blueshifts (e.g., Dravins 1987). The observation and interpretation of unresolved stellar granulation is not limited to the Sun (Nordlund et al. 2009) because numerical simulations cover a substantial portion of the Hertzsprung-Russell diagram (Magic et al. 2013; Trampedach et al. 2013; Ludwig et al. 2009), including the evolutionary phases from the main sequence over the turnoff up to the red-giant branch for low-mass stars.

Since the discovery of 51 Peg (Mayor & Queloz 1995), various studies have looked at starspots. For instance, Saar et al. (1998) proposed the first quantitative impact of starspots on radial-velocity measurements. The authors studied the impact of these surface structures on the bisector (i.e., measurement of the spectral line asymmetries) global slope and found that convection leads to bisector variations up to a few tens m s⁻<fnref n="1" />. Saar & Donahue (1997) pointed out they can lead to even larger radial-velocity variations for G2V-type stars. It should be expected that, in the case of F dwarfs or K giants, the velocity fields would be even larger. Paulson et al. (2004) measured star-to-star variations of 50 m s⁻<fnref n="1" /> due to stellar activity in a sample of Hyades dwarfs. Desort et al. (2007) discuss the possibility that, in F-K type stars, radial-velocity variations may be due to either spots or planets. Sanchis-Ojeda & Winn (2011) show that the transit data of a super-Neptune planet exhibit numerous anomalies that they interpret as passages over dark spots.

The role of long-baseline interferometric observations in planet hunting is a complement to the radial velocity and adaptive optics surveys. Thanks to the higher angular resolution, interferometry is the ideal tool for exploring separations in the range 1 to 50 mas (Le Bouquin & Absil 2012). This is achieved by observing the closure phase measurements directly associated with the asymmetries in the brightness distribution and, as a consequence, off-axis detection of a companion. Long-baseline interferometry bridges the gap between the use of direct imaging, which finds wide companions, and the use of RV measurements, which detect close companions (Le Bouquin & Absil 2012). Several attempts and discussions regarding prospective ideas towards this end have already been carried out, in particular for hot Jupiter planets with the MIRC instrument at CHARA telescope (e.g., Zhao et al. 2008, 2011; van Belle 2008) or the AMBER, MIDI, PIONIER instruments at VLTI (e.g., Matter et al. 2010; Absil et al. 2011; Chiavassa et al. 2012; Lachaume & Berger 2013). However, extracting the planetary signal from the interferometric observables is a difficult task that requires very accurate precision levels, possible only with proportionate increase in the data signal-to-noise ratio.

In this work, we present interferometric predictions obtained from 3D surface convection simulations run for stars spanning different effective temperatures, surface gravities, and metallicities. Furthermore, we present results from a study of the impact of granulation on the detection of transiting planet for three

prototypes of planets of different sizes corresponding to a hot Jupiter, a hot Neptune, and a terrestrial planet.

## 2. Stellar model atmospheres

Magic et al. (2013) describe the large STAGGER-grid of realistic 3D radiative hydrodynamical (RHD) simulations of stellar convection for cool stars using STAGGER-code (originally developed by Nordlund & Galsgaard 1995¹, and continuously improved over the years by its user community), which is a state-of-the-art (magneto)hydrodynamic code that solves the time-dependent equations for conservation of mass, momentum, and energy. The code uses periodic boundary conditions horizontally and open boundaries vertically. At the bottom of the simulation, the inflows have constant entropy and pressure. The outflows are not constrained and are free to pass through the boundary. The code is based on a sixth-order explicit finite-difference scheme and a fifth-order interpolation. The considered large number over wavelength points is merged into 12 opacity bins (Nordlund 1982; Skartlien 2000). The equation-of-state accounts for ionisation, recombination, and dissociation (Mihalas et al. 1988). The opacities include continuous absorption and scattering coefficients as listed in Hayek et al. (2010) and the line opacities as described in Gustafsson et al. (2008), in turn based on the VALD-2 database (Stempels et al. 2001) of atomic lines and the SCAN-base (Jørgensen 1997) of molecular lines.

For the solar abundances, the authors employed the latest chemical composition by Asplund et al. (2009), which is based on a solar simulation performed with the same code and atomic physics as in Magic et al. (2013).

## 3. Three-dimensional radiative transfer

We used pure-LTE radiative transfer OPTIM3D (Chiavassa et al. 2009) to compute synthetic images from the snapshots of the RHD simulations of the STAGGER-grid (see Fig. 1 of Magic et al. 2013). The code takes the Doppler shifts occurring due to convective motions into account. The radiative transfer equation is solved monochromatically using pretabulated extinction coefficients as a function of temperature, density, and wavelength.

The lookup tables were computed for the same chemical compositions (Asplund et al. 2009) as the RHD simulations using the same extensive atomic and molecular continuum and line opacity data as the latest generation of MARCS models (Gustafsson et al. 2008). We assume zero microturbulence and model the non-thermal Doppler broadening of spectral lines using only the self-consistent velocity fields issued from the 3D simulations. The temperature and density ranges spanned by the tables are optimized for the values encountered in the RHD simulations. The detailed methods used in the code are explained in Chiavassa et al. (2009). OPTIM3D has already been employed in synergy with STAGGER-code within several works (Chiavassa et al. 2010, 2011, 2012) concerning either the extraction of interferometric observables or synthetic spectra.

## 4. Interferometric observable construction

The aim of the present work is to present a survey of the convective pattern ranging from the optical to the far infrared and to

<footnote n="1">¹ http://www.astro.ku.dk/~kg/Papers/MHD_code.ps.gz</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Synthetic stellar disk images of the RHD simulations of Table 1 (columns). The images correspond to a representative wavelength for each interferometric instrument of Table 2 from the optical (top row) to the far infrared (bottom row). The averaged intensity (×10⁵ erg cm⁻² s⁻¹ Å⁻¹) is reported in the lower left corner of each image.</figcaption></figure>

evaluate its effect on detection of planet transit. We chose representative simulations in the STAGGER-grid partially covering the *Kepler* planets to study the effect of the metallicity across the HR-diagram to cover typical *Kepler* planets, and included different metallicities for the solar model (Fig. 3). Our statistical approach aims to present results that can be extrapolated to other stars in the Hertzsprung-Russel. More detailed analysis with respect to particular stellar parameters can be conducted using specific simulations of the STAGGER grid.

We used OPTIM3D to compute intensity maps from the snapshots of the RHD simulations of Table 1 for different inclinations with respect to the vertical, $\mu \equiv \cos(\theta) = [1.000, 0.989, 0.978, 0.946, 0.913, 0.861, 0.809, 0.739, 0.669, 0.584, 0.500, 0.404, 0.309, 0.206, 0.104]$ (these angles have already been used in the previous works of Chiavassa et al. 2012, 2010), and for a representative series of ten snapshots adequately spaced apart so as

to capture several convective turnovers for each simulation. The wavelength range is between 4000 and 52 000 Å with a spectral resolution $\lambda/\Delta\lambda = 20 000$.

### 4.1. From a small portion of the stellar surface to spherical tile images

The computational domain of the RHD simulations is limited to a small representative volume located in the stellar photosphere including the top of the stellar convective envelope, the horizontal directions chosen so as to be large enough to cover an area corresponding to about ten granular cells. The intensity maps computed with OPTIM3D are limited to a small portion of the stellar surface (see, e.g., Fig. 1 of Chiavassa et al. 2010), thus to overcome this limitation, we applied the same method as explained in Chiavassa et al. (2010) to tile a spherical surface
<table><caption>**Table 1.** 3D simulations from STAGGER-grid used in this work.</caption><thead><tr><th>⟨T<sub>eff</sub>⟩<sup>a</sup><br>[K]</th><th>[Fe/H]</th><th>log g<br>[cgs]</th><th>x, y, z-dimensions<br>[Mm]</th><th>x, y, z-resolution<br>[grid points]</th><th>M<sub>*</sub><br>[M<sub>☉</sub>]</th><th>R<sub>*</sub><br>[R<sub>☉</sub>]</th><th>Number of tiles<br>over the diameter</th></tr></thead><tbody><tr><td>5768.51 (Sun)</td><td>0.0</td><td>4.4</td><td>3.33 × 3.33 × 2.16</td><td>240 × 240 × 240</td><td>1.0</td><td>1.0</td><td>286</td></tr><tr><td>5764.13</td><td>-1.0</td><td>4.4</td><td>3.12 × 3.12 × 1.63</td><td>240 × 240 × 240</td><td>1.0</td><td>1.0</td><td>305</td></tr><tr><td>5781.04</td><td>-2.0</td><td>4.4</td><td>2.75 × 2.75 × 1.67</td><td>240 × 240 × 240</td><td>1.0</td><td>1.0</td><td>347</td></tr><tr><td>5780.06</td><td>-3.0</td><td>4.4</td><td>3.00 × 3.00 × 1.61</td><td>240 × 240 × 240</td><td>1.0</td><td>1.0</td><td>318</td></tr><tr><td>4569.23</td><td>0.0</td><td>2.0</td><td>1000 × 1000 × 1288</td><td>240 × 240 × 240</td><td>1.3<sup>b</sup></td><td>18.9</td><td>17</td></tr><tr><td>5001.35</td><td>0.0</td><td>3.5</td><td>27.08 × 27.08 × 24.49</td><td>240 × 240 × 240</td><td>1.15<sup>c</sup></td><td>3.1</td><td>121</td></tr><tr><td>5993.42</td><td>0.0</td><td>4.0</td><td>10.83 × 10.83 × 5.66</td><td>240 × 240 × 240</td><td>1.0<sup>c</sup></td><td>1.6</td><td>266</td></tr><tr><td>5998.93</td><td>0.0</td><td>4.5</td><td>2.92 × 2.92 × 1.76</td><td>240 × 240 × 240</td><td>1.15<sup>c</sup></td><td>0.99</td><td>312</td></tr></tbody></table>

**Notes.** (a) Horizontally and temporal average of the emergent effective temperatures from Magic et al. (2013). (b) Averaged value from Fig. 4 of Mosser et al. (2012). (c) Averaged value from Fig. 2 of Silva Aguirre et al. (2011).

accounting for limb-darkened effects. The computed value of the $\theta$-angle used to generate each map depended on the position (longitude and latitude) of the tile on the sphere and was linearly interpolated among the inclination angles.

In addition to this, the statistical tile-to-tile fluctuations (i.e., number of granules, shape, and size) is considered by selecting random snapshots within each simulation's time series. As a consequence, the simulation assumption of periodic boundary conditions resulted in a tiled spherical surface displaying an artificial periodic granulation pattern. However, Chiavassa et al. (2010) proved that the signal artificially introduced into the interferometric observables is weaker than the signal caused by the inhomogeneities of the stellar surface.

We estimated a stellar radius based on an applied mass taken from the literature (Table 1, 6th column), then we computed the number ($N_{\text{tile}}$) of tiles needed to cover half a circumference from side to side on the sphere $N_{\text{tile}} = \frac{\pi R_{\odot}}{x, y\text{-dimension}}$, where $R_{\odot}$ (transformed in Mm) and the $x, y$-dimension come from Table 1.

The final result is an orthographic projection of the tiled spheres (Fig. 1). It must be noted that our method of constructing realisations of stellar disk images inevitably introduces some discontinuities between neighbouring tiles by randomly selecting temporal snapshots and by cutting intensity maps at high latitudes and longitudes. The figure shows that the centre-to-limb variations are more pronounced in the optical instruments than in the infrared ones. This effect has already been found in Chiavassa et al. (2012, 2010) for some K-giant and sub-giant stars and is explained by different sensitivities of the source (Planck) function at optical and at infrared wavelengths.

## 4.2. Choice of interferometric instruments

Actual interferometers ensure the wavelength coverage from optical to far infrared with a series of instruments mounted at different sites. Table 2 displays the instruments we chose, where they are mounted, and the number of telescopes recombined, as well as the wavelengths probed. We used the online Astronomical Software for Preparing Observations (ASPRO2) of the JMMC<fnref n="2" /> to extract an OIFITS file with the telescope real positions in the UV-plane, telescope configurations, and observing wavelength. Afterwards, we performed a top-hat average over the whole set of disk images to obtain one synthetic image for each observing wavelength. Even if a wavelength dependence exists on the interferometric observables, for simplicity, we assume a representative wavelength for each instrument in the rest of the work.

In Sect. 5, we introduce the closure phase observable and study its potentiality for the detection of surface related convective structures. We do not aim to interpret/observe a particular star and thus, for each instrument, we chose an arbitrary date and arbitrary star with coordinates that ensure observability for the whole night. This choice is taken to accommodate a broad coverage of the spatial frequencies up to the fifth to sixth lobe when possible (Fig. 4). Owing to the sparse selections of baselines (i.e., different apparent size of the targets), using this approach it is not possible to directly compare the instruments; however, in this section, we aim to present a closure phase survey of the convective pattern from the optical to the far infrared.

Finally, in Sect. 6, we investigate a more concrete scenario with the choice of two real targets either in the visible and in the infrared. Thanks to our considering fixed targets for visible and infrared instruments, we can directly compare the results among the different instruments and propose the best instrument and/or interferometric facility to detect the stellar granulation.

Owing to either convection-related structures or a faint companion, the stellar surface asymmetries in the brightness distribution affect closure phases. The sum of three phases around a closed triangle of baselines is the closure phase: this procedure removes the atmospheric contribution, leaving the phase information of the object visibility unaltered (Monnier 2007, 2003). Closure phases have the main advantage of being uncorrupted by telescope-specific phase errors, including pointing errors, atmospheric piston, and longitudinal dispersion due to air and water vapour (Le Bouquin & Absil 2012). Closure phase errors, when known, are reported in Table 2. Owing to the sparse structure of the point spread function associated with the diluted aperture of an interferometer (see Fig. 4), however, the position and the morphology of these surface inhomogeneities depend on their relative orientation and on the interferometric baselines.

Figure 1 shows irregular stellar surfaces with convection-related surface structures, whose sizes depend on the stellar parameters of the simulations. There are pronounced centre-to-limb variations in the optical (VEGA to NPOI instruments), while these are less noticeable in the infrared. This is mainly due to the differences in Planck functions in the optical range and in the infrared region.

Starting from these synthetic images, we followed the method described in Chiavassa et al. (2009) to calculate the discrete complex Fourier transform $\mathcal{F}$ for each image, with particular interest only in the closure phases. From the OIFITS files for each instrument, we know the set of all baseline vectors (i.e., UV-plane) of all the telescopes (see Sect. 4.2), and we matched their frequencies in the UV plane with the

<footnote n="2">² http://www.jmmc.fr/aspro_page.htm</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Enlargement of the synthetic stellar disk images of Fig. 1 for the VEGA instrument (Table 2).</figcaption></figure>

corresponding points in the Fourier transform of the synthetic images. The phase for each telescope is $\tan\varphi = \mathfrak{J}(\mathcal{F})/\mathfrak{R}(\mathcal{F})$, where $\mathfrak{J}(\mathcal{F})$ and $\mathfrak{R}(\mathcal{F})$ are the imaginary and real parts of the complex number $\mathcal{F}$, respectively. Finally, the closure phase is the sum of all phase differences between closed triangles of telescope baselines: e.g., for three telescopes: $CP(1-2-3) = \Phi_{1-2} + \Phi_{2-3} + \Phi_{3-1}$, where $CP(1-2-3)$ is the closure phase, and $\Phi_{1-2}$ the arctan of the Fourier phases $\tan\varphi$ for the telescopes $1-2$.

**5. Closure phase as an indicator of the stellar inhomogeneity**

In our survey, we used the setup of different instruments (Table 2) characterised by a number of telescopes ($N$) varying from 3 to 6. Monnier (2003) show that the possible closed triangles of baselines (i.e., one closed triangle gives one closure phase) is $\binom{N}{3} = \frac{(N)(N-1)(N-2)}{3(2)}$, but the independent Fourier phases are given by $\binom{N}{2} = \frac{(N)(N-1)}{(2)}$, and thus not all the closure phases are independent but only $\binom{N-1}{2} = \frac{(N-1)(N-2)}{2}$.

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** RHD simulations from STAGGER-grid used in this work with the aimed effective temperature (see also Table 1) over-plotted to the Kepler planets (confirmed and candidates) in fall 2013 from http://exoplanets.org (Wright et al. 2011).</figcaption></figure>

The number of independent closure phases is always less than the number of phases one would like to determine, but the percent of phase information retained by the closure phases improves as the number of telescopes in the array increases (Monnier 2003).

To sum up:

- with three telescopes one obtains one closed triangle of baselines (i.e., 1 closure phase), three Fourier phases, and one independent closure phase;

- with four telescopes one obtains four closed triangles of baselines (4 closure phases), six Fourier phases, and three independent closure phases;

- with six telescopes one obtains 20 closed triangles of baselines (20 closure phases), 15 Fourier phases, and 10 independent closure phases.

Chiavassa et al. (2012, 2010) demonstrate that, in the case of Procyon and K-giant stars, the synthetic visibility curves produced by the RHD simulations are systematically different from spherical symmetric modelling, with an impact on the radius, effective temperature, and departures from symmetry. This was noticeable at higher spatial frequencies and mostly affecting the signal of the closure phases. The authors interpreted this as the signature linked to the convection-related surface structures. Starting from these remarks, we decided to concentrate our survey study only on the closure phases.

Figure 5 displays closure phases deviating from the axisymmetric case that particularly occur in optical filters, where the dispersion is larger (e.g., VEGA, NPOI, and PAVO instruments). Depending on the instruments and spatial frequency spanned, the departures from symmetry may be large or not. However, it is apparent that the convection-related surface structures have a signature on the closure phases.

The characterisation of the granulation signature is analysed in Fig. 6, where the departure from axisymmetric case (i.e.,
<figure><img src="image_5.png" /><figcaption>Fig. 4. Typical UV coverage in metres for the different instruments of Table 2 over-plotted to the Fourier transform of the RHD simulation of the Sun (Table 1). Red dots correspond to the telescope baseline positions during the arbitrary observation we prepared. The observability is assured for approximatively a whole night with a large enough number of Earth rotation aperture synthesis baseline points. The arbitrary apparent size of the star is reported in the top right corner of each panel.</figcaption></figure>

closure phase different from zero or $\pm\pi$) of all the RHD simulations and interferometric instruments is displayed. We proceed as follows:

- for all the instruments, we determined the difference, $\psi_i$, between the synthetic closure phases and the axisymmetric values zero or $\pm\pi$ from all spatial frequency $i$;

- based on the UV coverage of Fig. 4 and the synthetic visibilities, we identified the frequency limits of the lobes (Fig. 5);

- we averaged $\bar{\psi}_i$ over the frequencies $i$ falling inside the lobe's limits;

- in case of multiple closed triangles (i.e., multiple values of closure phases such as for instruments working with four or six telescopes), we selected the largest $\bar{\psi}_i$.

It is remarkable that all the simulations show departure from the axisymmetric case at all the wavelengths. At least for the chosen configurations, it is difficult to determine clear differences among the stellar parameters and, in particular, for the different metallicities of the solar simulations. In addition to this, it must be noted that the averaged $\bar{\psi}_i$ may smooth out the differences even if these observables are useful for pointing out the signature of the convection-related surface structures.

The levels of asymmetry and inhomogeneity of stellar disk images reach very high values of several tens of degrees with stronger effects from third visibility lobe on. In this work we assumed precise values for the arbitrary observations (see Fig. 4). To estimate the baseline needed for other stellar sizes, the following equation can be used to retrieve second zero of the visibility curve (i.e., the third lobe):

$$B[m] = 2.23 \cdot \frac{\lambda}{\theta[\text{rad}]} = 2.23 \cdot \frac{\lambda \cdot 206\,265 \cdot 1000}{\theta[\text{mas}]} \quad (1)$$

where $\theta$ is the apparent diameter, $\lambda$ the wavelength in meters, and $B$ the baseline in meters. For example, a star observed at wavelength 0.7 µm (e.g., VEGA instrument) with $\theta = [0.5, 2, 5, 10]$ mas the third lobe will be probed with a baselines of about [645, 165, 65, 33] m. The same example in the H band 1.6 µm (e.g., PIONIER instrument) gives baselines of about [1472, 368, 148, 74] m, respectively, and for the L band at 30 µm with baselines of [27600, 2900, 2760, 1380] m, respectively.

Finally, to estimate the distance of the observed star the following equation can be used:

$$d \text{ [pc]} = \frac{R[R_{\odot}] \cdot 9.305}{\theta[\text{mas}]} \quad (2)$$

where $d$ is the distance of the observed star in parsec, $R[R_{\odot}]$ the radius of the star in solar radii, and $\theta$ the apparent diameter. The optical and the near infrared wavelengths are more affordable in terms of baseline lengths because they are less than 400 m for stellar sizes over 2 mas, it is however more complicated for the mid-infrared wavelengths where the baseline lengths become kilometric. The signature on the closure phases can be evaluated by accumulating observations at short and long baselines. This can be ensured by the fact that the instrumental errors on closure phases are much smaller than the expected closure phase departures (see Table 2). It is, however, important to note that probing high frequencies, the signal-to-noise ratio of the measurements would be very low due to low fringe visibilities, greatly deteriorating the closure phase precision and affecting the instrument capability.
<table><caption>**Table 2.** Interferometric instruments/configurations used in this work.</caption><thead><tr><th>Name</th><th>Location</th><th>Max<br>baseline [m]</th><th>Wavelength<br>range [Å]</th><th>Closure<br>phase error [°]</th><th>Configuration<br>chosen</th><th>Number of<br>telescopes</th><th>Active<br>since/from</th><th>Reference</th></tr></thead><tbody><tr><td>VEGA</td><td>CHARA<sup>a</sup></td><td>~331</td><td>~6400–8800</td><td>–</td><td>W2-E2-S1-E1</td><td>4</td><td>2009</td><td>1</td></tr><tr><td>PAVO</td><td>CHARA</td><td>~331</td><td>~6500–8000</td><td>–</td><td>S1-W1-W2</td><td>3</td><td>2011</td><td>2</td></tr><tr><td rowspan="2">VISION</td><td rowspan="2">NPOI<sup>b</sup></td><td rowspan="2">~432</td><td rowspan="2">~5500–8500</td><td rowspan="2">–</td><td>AC0-AW0-AN0-E06</td><td>4</td><td>2014?</td><td>3</td></tr><tr><td>AC0-AE0-AW0-AN0-E06-W07</td><td>6</td><td></td><td>3</td></tr><tr><td>CLIMB</td><td>CHARA</td><td>~331</td><td>12862/21349<br>broad-band</td><td>–</td><td>S1-W1-W2</td><td>3</td><td>2005</td><td>4</td></tr><tr><td>MIRC</td><td>CHARA</td><td>~331</td><td>~14 000–18 000</td><td>0.1–0.2<sup>c</sup></td><td>S1-S2-W1-W2-E1-E2</td><td>6</td><td>2007</td><td>5</td></tr><tr><td>AMBER</td><td>VLTI<sup>d</sup></td><td>130</td><td>~21 200–25 000</td><td>0.20–0.37<sup>e</sup></td><td>A1-G1-J3</td><td>3</td><td>2004</td><td>6</td></tr><tr><td>PIONIER</td><td>VLTI</td><td>130</td><td>16810/20510<br>broad-band</td><td>0.25–1<sup>f</sup></td><td>A1-G1-K0-J3</td><td>4</td><td>2010</td><td>7</td></tr><tr><td>GRAVITY</td><td>VLTI</td><td>130</td><td>~20 000–24 000</td><td>1<sup>g</sup></td><td>A1-G1-K0-J3</td><td>4</td><td>2015?</td><td>8</td></tr><tr><td>MATISSE</td><td>VLTI</td><td>130</td><td>~28 600–52 000</td><td>&lt;1.16<sup>h</sup></td><td>A1-G1-K0-J3</td><td>4</td><td>2016?</td><td>9</td></tr></tbody></table>

**Notes.** Some of the instruments may cover other wavelength ranges or be used with different configurations, but what we chose is a good representation for the purpose of the present work. All the information about the instrument/configuration/wavelengths have been retrieved on ASPRO2, except for VISION (Garcia, priv. comm.). (a) ten Brummelaar et al. (2005); (b) Armstrong et al. (2014); (c) Zhao et al. (2011, 2010, 2008); (d) Haguenauer et al. (2008); (e) Absil et al. (2010) for medium resolution; (f) Le Bouquin et al. (2011); Absil et al. (2011); (g) Final Design Review, priv. comm.; (h) Lopez (2012).

**References.** (1) Mourard et al. (2009); (2) Ireland et al. (2008); (3) Ghasempour et al. (2012); (4) ten Brummelaar et al. (2005); (5) Monnier et al. (2004); (6) Petrov et al. (2007); (7) Le Bouquin et al. (2011); (8) Eisenhauer et al. (2008); (9) Lopez et al. (2008).

# 6. Applications

## 6.1. Study of the granulation on two real targets

In the previous section we showed that the detection of closure phase departures from symmetry needs stars resolved up to the third and fourth lobes, as well as in some cases, the fifth or sixth lobes. This reduces the number of targets that can be observed with the actual interferometric baselines but foresees the need for an extension of the next-generation interferometric infrastructures.

In this section, we performed the closure phase analysis for two real targets: Beta Com and Procyon (Table 3). Beta Com has been chosen for its smaller angular diameter so as to illustrate observations in the visible, while the large diameter of Procyon ensures a good UV coverage in the infrared. As before, we prepared adapted OIFITS files for each instrument using real observability coverage.

Figure 7 displays smaller departures in closure phases with respect to the configurations taken in Fig. 6. In general all the instruments (except for MATISSE and NPOI, which do not probe frequencies greater than the first lobe) show closure phases departures (ψ) of a few degrees with the highest values of ~16° (to be compared with instrumental errors of Table 2). Since in this instance we only considered one single target for all the visible instruments and another for the infrared instruments, we can directly compare the results among the different instruments. Both PAVO and VEGA give values lower than 0.5° with a closure phase signal starting from the lower lobe for VEGA. More interesting is the infrared region with VLTI's instruments showing departures already from the second lobe: AMBER and GRAVITY with values lower than 0.8°, PIONIER in the H band with values of 4.3° (2nd lobe) 6.4° (3rd lobe), and PIONIER in the K band with 2.9° (2nd lobe). CHARA's instruments do not show departures on the second lobe, but they probe higher frequencies up to the sixth lobe with MIRC (13.8°, 15.3°, 16.4°, 13.1° for the 3rd, 4th, 5th, 6th lobes), CLIMB J band (6.5°, 0.2°, 5.1° for the 3rd, 4th, 5th lobes), and CLIMB K band (12.3°, 6.3° for the 3rd, 4th lobes).

The actual instruments and telescopes, with the errors on closure phases reported in Table 2, allow, in principle, the detection

<table><caption>**Table 3.** Reference targets and associated RHD simulations of Table 1.</caption><thead><tr><th>Name</th><th>Spectral type</th><th>Angular diameter [mas]</th><th>RHD [⟨T<sub>eff</sub>⟩]</th></tr></thead><tbody><tr><td>Beta Com (HD 114710)</td><td>G0V</td><td>1.1<sup>(a)</sup></td><td>Sun</td></tr><tr><td>Procyon (HD 62421)</td><td>F5IV 5.4</td><td>5.4<sup>(b)</sup></td><td>5993.42</td></tr></tbody></table>

**Notes.** (a) Richichi et al. (2005); (b) Chiavassa et al. (2012).

of the granulation. The closure phase signal is already more pronounced in the infrared for the second lobe and may be detected with very good weather and instrumental conditions, but it is certainly easier to detect from the third lobe on. The long baseline set of CHARA telescopes is more advantageous as higher frequencies are probed. Moreover, MIRC instrument with six telescope recombination is the most appropriate instrument because it combines good UV coverage and long baselines. Our analysis is based on the assumption of very good conditions, but in the context of non-zero exposure times and the presence of atmospheric turbulence, the accuracy on closure phase measurements is also affected by piston noise. The statistical uncertainty on the closure phases depends on the atmospheric conditions with the piston-noise contribution decreasing with the square root of the integration time (Le Bouquin & Absil 2012). Optimal observing strategy could, however, be defined to reach the needed accuracy as already demonstrated by some instruments of Table 2.

## 6.2. Transiting planet

Chiavassa et al. (2012) explored the impact of the convection on interferometric planetary signature around a RHD simulation of a sub-giant star. The authors estimated the impact of the granulation noise on a hot Jupiter detection using closure phases and found that there is a non-negligible and detectable contamination to the signal of the hot Jupiter due to the granulation from spatial frequencies longward of the third lobe. In this work, we extended this analysis to all the simulations of Table 1 using the following procedure:

- we chose three prototypes of planets representing different sizes and compositions (Table 4). However, the purpose of
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Scatter plot of closure phases of 20 000 random baseline triangles (black dots) as a function of the maximum linear extension corresponding to the configuration chosen for each instrument of Table 2 and for the RHD simulation of the Sun (Table 1). The coloured symbols over-plotted display the closure phases for the configuration chosen (see the UV-planes of Fig. 4), and the vertical dashed red lines give the approximate positions of the different lobes.</figcaption></figure>

this work is not to reproduce the exact conditions of the
planet-star system detected but to employ a statistical ap-
proach to the interferometric signature for different stellar
parameters hosting planets with different sizes;

- we simulated the transit of those planets for stars with stellar parameters of RHD simulations of Table 1. A representative example is reported in Fig. 8;

- we computed the closure phases for three planet-star sys-
tem images corresponding to three particular planet tran-
sit phases. The different selections Instrument+Wavelength
(configurations reported in Sect. 5) of the synthetic images
is a representative choice among the numerous possibilities;

- we determined the difference between the planet-star system
  and the star alone (Fig. 9).

For modelling the flux of the irradiated planet, we use the same
prescriptions as Chiavassa et al. (2012). Since our main interest
is related to the impact of the planet size on the interferomet-
ric observables, and the flux of the planet is much smaller than
the stellar flux, we used the same model for the planetary flux
(Barman et al. 2001), in particular, the spectra of hot extrasolar
planet around a generic cool star.

Our interest is related to the closure-phase signature due
to the planet with respect to the stellar granulation. Figure 8
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Differences, $\psi$, between the closure phases and zero or $\pm\pi$ (i.e., departure from axisymmetric case) for all the simulations of Table 1 (horizontal axis) and all the instruments of Table 2 (vertical axis). $\psi$ has been averaged over the spatial frequencies corresponding to the different lobes spanned by the instruments configuration (see text), and for each simulation, 6 lobes are displayed: black for the 1st lobe, violet for the 2nd, light blue for the 3rd, green for the 4th, yellow for the 5th, and red for the 6th. Only the lobes spanned in the UV planes of Fig. 4 are plotted. The symbols correspond to different values, in degrees.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Same as in Fig. 6 but for two real targets: Beta Com and Procyon (see Table 3 for the associated RHD simulation). MATISSE and NPOI do not probe frequencies higher than the first lobe and have not been reported. The sizes of the circles correspond to different values in degrees. The values of $\psi$ are lower than in Fig. 6. The crossed circle means that there are no detected departures.</figcaption></figure>

displays the geometrical configurations of the planet-star system for the representative example of the Sun. As already reported in Chiavassa et al. (2012), the ratio between the stellar intensity and the planet's integrated intensity is stronger in the infrared with respect to the optical. The stellar intensity, $I_{\text{star}}$, at its centre ($\mu = 1$) for the synthetic images of the simulations of Fig. 8,

<table><caption>**Table 4.** Prototypes of planets chosen to represent the planet transit phases of Fig. 8.</caption><thead><tr><th>Name</th><th>Jupiter mass</th><th>Jupiter radius</th><th>Semi-axis [AU]</th><th>Real hosting star T<sub>eff</sub> [K]/log g</th></tr></thead><tbody><tr><td>Kepler-11 f<sup>a</sup></td><td>0.006</td><td>0.222</td><td>0.2504</td><td>5663/4.37</td></tr><tr><td>HD 149026 b<sup>b</sup></td><td>0.360</td><td>0.654</td><td>0.0431</td><td>6160/4.28</td></tr><tr><td>CoRoT-14 b<sup>c</sup></td><td>7.570</td><td>1.090</td><td>0.0269</td><td>6040/4.45</td></tr></tbody></table>

Notes. (a) Lopez & Fortney (2013); (b) Sato et al. (2005); (c) Tingley et al. (2011).

and the planet integrated intensity, $I_{\text{planet}}$, at the wavelength corresponding to the instruments of Table 2 are reported in Table 5.

We considered three particular planet transition phases (Fig. 8) corresponding to the ingress and egress of the transition, as well as to the planet at the centre of the stellar disk. The resulting absolute differences (in degrees) between the closure phases of the planet-star system with the ones of the star alone are in Fig. 9. The figure shows that, for all the instruments, the absolute difference scales with the size of the planet considered: the smaller planet returns smaller differences. There is, however, an exception for the *H* band-MIRC and MATISSE where there is not a clear distinction between the different planets because the baselines probe very high spatial frequencies (Fig. 4) and thus finer details.

Moreover, the closure phase differences are greater in the optical wavelengths where the stellar surface is not flat but rather “corrugated”, owing to the larger fluctuations and the higher contrast of granulation than in the infrared. Finally, while for some instruments it is not possible to disentangle the transition phase of the planet (because of the configurations chosen and/or the
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Synthetic stellar disk images in PIONIER H band, together with three planet transiting phases (black colour) for the Sun of Table 1. The prototypes of planet are from Table 4 and are Kepler-11 f prototype planet (**left column**), HD 149026 b prototype (**central column**), and CoRoT-14 b prototype (**right column**).</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Absolute closure-phase differences (in degrees) between the star with a transiting planet (Fig. 8) and the star alone (Fig. 1) for all the instruments of Table 2. The black colour corresponds to the smallest prototype planet Kepler 11-f of Table 4, the red to the intermediate planet HD 149026 b, and the blue to the largest planet CoRoT 14-b. The star symbols connected with solid lines correspond to the planet phase entering the stellar disk (see Fig. 8), the circle symbols connected with dotted line to the planet at the centre of the stellar disk, and the triangles connected with dashed line to the planet exiting the stellar disk.</figcaption></figure>

spatial frequencies spanned), for others (VEGA, PAVO, CLIMB, and AMBER), it is clear that different transit positions have different effects on the closure phases (also shown in Ligi & Mourard, in prep.).

The signature of the transiting planet on the closure phase is mixed with the signal due to the convection-related surface structures. The time scale of granulation depends on the stellar parameters and varies from minutes or tens of minutes for solar type stars and sub-giants to hours for more evolved red giant stars. If the transit is longer that the granulation time scale (which is the case for most of main sequence stars), it is possible to disentangle its signal from convection by observing at particular wavelengths (either in the infrared or in the optical) and

measuring the closure phases for the star at difference phases of the planetary transit.

For this purpose, it is very important to have comprehensive knowledge of the host star to detect and characterise the orbiting planet, and RHD simulations are very important to reach this aim.

## 6.3. Closure phases impact: granulation versus limb darkened law

We show in this section that the planet detection with closure phases is strongly influenced by the intrinsic stellar granulation presented in Sect. 5. For this purpose, we computed images without stellar granulation and used the limb-darkened law
<table><caption>**Table 5.** Stellar intensity, *I*<sub>star</sub>, at its centre (μ = 1) for the synthetic images of the simulation of Fig. 8, and the planet integrated intensity, *I*<sub>planet</sub>, at the wavelength corresponding to the instruments of Table 2.</caption><thead><tr><th>Instrument</th><th>VEGA</th><th>PAVO</th><th>NPOI<br>4 Tel</th><th>NPOI<br>6 Tel</th><th>CLIMB<br>J band</th><th>MIRC<br>6 Tel</th></tr></thead><tbody><tr><td>Representative λ used</td><td>7312 Å</td><td>6400 Å</td><td>5669 Å</td><td>5817 Å</td><td>12 862 Å</td><td>15 940 Å</td></tr><tr><td><i>I</i><sub>star</sub>/<i>I</i><sub>planet</sub></td><td>1702</td><td>5524</td><td>13 228</td><td>82 619</td><td>148</td><td>29</td></tr></tbody><thead><tr><th>Instrument</th><th>AMBER</th><th>PIIONIER<br>H band</th><th>GRAVITY</th><th>PIIONIER<br>K band</th><th>CLIMB<br>K band</th><th>MATISSE<br>LM band</th></tr></thead><tbody><tr><td>Representative λ used</td><td>23 015 Å</td><td>16 810 Å</td><td>22 000 Å</td><td>20 510 Å</td><td>21 350 Å</td><td>28 675 Å</td></tr><tr><td><i>I</i><sub>star</sub>/<i>I</i><sub>planet</sub></td><td>84</td><td>15</td><td>16</td><td>22</td><td>15</td><td>13</td></tr></tbody></table>

and coefficients of Claret (2000). We proceeded using appropriate limb-darkened coefficients for the wavelength range of the interferometric instruments of Table 2 and for the same stellar parameters of RHD simulations as in Table 1. Then, we simulated the planet transitions and computed the resulting closure phases using the same approach and prototypes of planets as in the previous section.

Figure 10 (top) displays a typical example for the limb-darkened image of the Sun. The bottom panel of the figure shows the comparison between the closure phases of the synthetic images from the RHD simulation of the Sun and the ones from the limb-darkened image with the transiting planet. The closure phases of a limb disk without the presence of inhomogeneities on its surface is zero or ±π, while the transiting planet causes very small departures from spherical symmetry. It is evident that the departure from zero or ±π due to the convection-related surface structures are much larger than what it is expected by the transiting planet on axisymmetric images. This result is similar for all instruments and stellar parameters employed in this work.

It is essential to use reliable RHD hydrodynamical simulation for preparing and interpreting observations aimed to detect and characterise planets.

## 6.4. Magnetic starspots impact on closure phases

An immediate problem for detecting transiting planets is signal contamination from starspots caused by the magnetic field of the star. Starspots are created by local magnetic field on the stellar surface and they appear as cool (and therefore dark) regions as compared to the surrounding surface. This is due to the inhibition of the convective motions by a strong enough magnetic field that blocks or redirects the energy flow from the stellar interior (Strassmeier 2009). We used the intensity map from the RHD simulation of the Sun (Table 1) in the MIRC instrument wavelengths (Table 2) to construct a hypothetical starspot image. We chose to put four spots at different longitudes and distances from the centre (Table 6). The difference between the photosphere and the spot temperatures is up to 2000 K for F and early G stars and down to 200 K for late M stars (Berdyugina 2005). In our case we assumed temperatures for the spots of <2000 K. We used spot size values between ~0.1% and ~10% of the stellar radius based on the large compilation of detected stellar spots with Doppler imaging (Strassmeier 2009).

Figure 11 (top) displays the resulting stellar disk image. The apparent size of the spots should be compared to the apparent sizes of the transiting planets of Fig. 8. The closure phase signal for the RHD simulation considering only the granulation and the one with starspots show non-negligible differences (bottom left panel of the figure), even though it seems difficult to differentiate from the granulation signal due to its chaotic behaviour.

Moreover, it is also visible in the bottom right-hand panel that the starspot signal on closure phases can nearly be the same as the transiting planet signal. Consequently, the planet signal may be contaminated.

Starspots caused by the magnetic field may pollute the granulation and the transiting planet signals, at least for the starspots configuration we considered. However, it should be possible to differentiate the transiting planet signal since the time scale of a planet crossing stellar disk is much smaller than the typical rotational modulation of the star. A more detailed analysis will be reported in a forthcoming paper.

## 7. Conclusions

We presented an application of the STAGGER-grid of realistic, state-of-the-art, time-dependent, radiative-hydrodynamic stellar atmosphere. We used the simulations to provide synthetic images from the optical to the infrared and extract interferometric observables aimed to study stellar convection, as well as its impact on planet detection and characterisation. RHD simulations are essential for a proper quantitative analysis of interferometric observations and crucial for extracting the signal.

We analysed the impact of convection at different wavelengths using the closure phases. Closure phase is the interferometric observable with intrinsic and unaltered information about the stellar surface asymmetries in the brightness distribution, either owing to convection-related structures or to a faint companion. We made our predictions as real as possible using actual interferometric instruments and configurations. All the simulations show departure from the axisymmetric case (closure phases not equal to 0 or ±π) for all the wavelengths, but at least for the chosen configurations, it is difficult to determine clear differences among the stellar parameters and, in particular, for the different metallicities of the solar simulations. The levels of asymmetry and inhomogeneity of stellar disk images reach very high values of several tens of degrees with stronger effects from the third visibility lobe on. We explored the possibility of detecting the granulation pattern on two real targets (Beta Com and Procyon). We found that the detection on the second lobe is possible either in the visible or in the near infrared with closure phase departures of less than 1°; detections on the third, fourth, fifth, and sixth lobes (with departures up to 16°) are possible using CHARA's instruments, and, in particular, MIRC is the most appropriate instrument because it combines good UV coverage and long baselines. In general, interferometers probing optical and near infrared wavelengths are more adapted to reaching higher spatial frequencies because the third visibility lobe can be probed with baseline lengths less than 400 m for stellar sizes larger than 2 mas. It is more complicated for the mid-infrared wavelengths where the baselines become kilometric. We emphasise that stars
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Top: Limb-darkened image (made using the law and coefficients of Claret 2000) for AMBER instrument, together with three planet transiting phases (black) for a star with the Sun stellar parameters. Bottom: scatter plot of closure phases computed for the Sun with transiting planets using limb-darkening unidimensional models without granulation (coloured symbols) versus closure phases of the corresponding RHD simulation of Table 1 (black symbols). The dashed line indicates the zero degree.</figcaption></figure>

<table><caption>**Table 6.** Parameters for the starspots.</caption><thead><tr><th></th><th>Spot 1</th><th>Spot 2</th><th>Spot 3</th><th>Spot 4</th></tr></thead><tbody><tr><td>Size [% of stellar radius]</td><td>10.0</td><td>1.1</td><td>1.5</td><td>2.2</td></tr><tr><td>Temperature [K]</td><td>3800</td><td>3900</td><td>4100</td><td>3900</td></tr><tr><td>Longitude [°]</td><td>-45.</td><td>10.</td><td>30.</td><td>-10.</td></tr></tbody></table>

should be observed at high spatial frequencies by accumulating observations on closure phases at short and long baselines.

We explored the impact of convection on interferometric planet signature for three prototypes of planets with sizes corresponding to one hot Jupiter, one hot Neptune, and a terrestrial one. Considering three particular planet transition phases, we compared the closure phases of the star with the transiting planet and the star alone. The signature of the transiting planet on the closure phase is mixed with the signal due to the convection-related surface structure, but it is possible to distinguish it at particular wavelengths (either in the infrared or in the optical). It can be achieved by measuring the closure phases

<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Top panel: synthetic stellar disk image of the Sun (Table 1) for the MIRC instrument with four darker starspots (see text) and parameters reported in Table 6. Central panel: scatter plot of closure phases computed for the Sun with starspots (black stars) and for the Sun (red triangles). Bottom panel: same as in bottom left panel but for the Sun with transiting planets.</figcaption></figure>
for the star at different phases of the transit. Starspots caused by the magnetic field of the star may masquerade as planets for interferometric observations. We showed that the starspot signal on closure phases can be nearly the same as the transiting planet signal (at least in the example configuration we considered). However, it should be possible to differentiate between them because the time scale of a planet crossing the stellar disk is much smaller than the typical rotational modulation of the star. It is, however, important to note that when probing high spatial frequencies, the signal-to-noise ratio of the measurements would be very low owing to low fringe visibilities, greatly deteriorating the closure phase precision and affecting the instrument capability. Moreover, this would influence the capability and sensitivity of detecting the signatures of granulation and disentangling the planetary signal.

The detection and characterisation of planets must be based on comprehensive knowledge of the host star, and this includes the detailed study of the stellar surface convection. In this context, RHD simulations are crucial for reaching this aim.

***Acknowledgements.*** R.C. is the recipient of an Australian Research Council Discovery Early Career Researcher Award (project number DE120102940). This research has made use of the Exoplanet Orbit Database and the Exoplanet Data Explorer at exoplanets.org.

## References

* Absil, O., Le Bouquin, J.-B., Lebreton, J., et al. 2010, A&A, 520, L2

* Absil, O., Le Bouquin, J.-B., Berger, J.-P., et al. 2011, A&A, 535, A68

* Armstrong, J. T., Restaino, S. R., Clark, J. H., et al. 2014, in AAS Meeting Abstracts, 223, 202.02

* Asplund, M., Ludwig, H., Nordlund, Å., & Stein, R. F. 2000, A&A, 359, 669

* Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481

* Barman, T. S., Hauschildt, P. H., & Allard, F. 2001, ApJ, 556, 885

* Berdyugina, S. V. 2005, Liv. Rev. Sol. Phys., 2, 8

* Carlsson, M., Stein, R. F., Nordlund, Å., & Scharmer, G. B. 2004, ApJ, 610, L137

* Chiavassa, A., Plez, B., Josselin, E., & Freytag, B. 2009, A&A, 506, 1351

* Chiavassa, A., Collet, R., Casagrande, L., & Asplund, M. 2010, A&A, 524, A93

* Chiavassa, A., Bigot, L., Thévenin, F., et al. 2011, J. Phys. Conf. Ser., 328, 012012

* Chiavassa, A., Bigot, L., Kervella, P., et al. 2012, A&A, 540, A5

* Claret, A. 2000, A&A, 363, 1081

* Desort, M., Lagrange, A.-M., Galland, F., Udry, S., & Mayor, M. 2007, A&A, 473, 983

* Dravins, D. 1987, A&A, 172, 211

* Eisenhauer, F., Perrin, G., Brandner, W., et al. 2008, in SPIE Conf. Ser., 7013, 69

* Ghasempour, A., Muterspaugh, M., Hutter, D., et al. 2012, in AAS Meeting Abstracts, 219, 446.13

* Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951

* Haguenauer, P., Abuter, R., Alonso, J., et al. 2008, in SPIE Conf. Ser., 7013, 11

* Hayek, W., Asplund, M., Carlsson, M., et al. 2010, A&A, 517, A49

* Herschel, W. 1801, Roy. Soc. Lond. Philosoph. Trans. Ser. I, 91, 265

* Ireland, M. J., Mérand, A., ten Brummelaar, T. A., et al. 2008, in SPIE Conf. Ser., 7013, 63

* Jørgensen, U. G. 1997, ed. E. F. van Dishoeck, IAU Symp., 178, 441

* 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

* Lopez, B. 2012, Publications de l'Observatoire Astronomique de Beograd, 91, 129

* Lopez, E. D., & Fortney, J. J. 2013, ApJ, submitted [arXiv:1311.0329]

* Lopez, B., Antonelli, P., Wolf, S., et al. 2008, in SPIE Conf. Ser., 7013, 70

* Ludwig, H., Caffau, E., Steffen, M., et al. 2009, Mem. Soc. Astron. It., 80, 711

* Magic, Z., Collet, R., Asplund, M., et al. 2013, A&A, 557, A26

* Matter, A., Vannier, M., Morel, S., et al. 2010, A&A, 515, A69

* Mayor, M., & Queloz, D. 1995, Nature, 378, 355

* Mihalas, D., Dappen, W., & Hummer, D. G. 1988, ApJ, 331, 815

* Monnier, J. D. 2003, Rep. Prog. Phys., 66, 789

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

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

* Mosser, B., Elsworth, Y., Hekker, S., et al. 2012, A&A, 537, A30

* Mourard, D., Clausse, J. M., Marcotto, A., et al. 2009, A&A, 508, 1073

* Nordlund, A. 1982, A&A, 107, 1

* Nordlund, Å., Stein, R. F., & Asplund, M. 2009, Liv. Rev. Sol. Phys., 6, 2

* Paulson, D. B., Cochran, W. D., & Hatzes, A. P. 2004, AJ, 127, 3579

* Petrov, R. G., Malbet, F., Weigelt, G., et al. 2007, A&A, 464, 1

* Richichi, A., Percheron, I., & Khristoforova, M. 2005, A&A, 431, 773

* Saar, S. H., & Donahue, R. A. 1997, ApJ, 485, 319

* Saar, S. H., Butler, R. P., & Marcy, G. W. 1998, ApJ, 498, L153

* Sanchis-Ojeda, R., & Winn, J. N. 2011, ApJ, 743, 61

* Sato, B., Fischer, D. A., Henry, G. W., et al. 2005, ApJ, 633, 465

* Silva Aguirre, V., Chaplin, W. J., Ballot, J., et al. 2011, ApJ, 740, L2

* Skartlien, R. 2000, ApJ, 536, 465

* Stempels, H. C., Piskunov, N., & Barklem, P. S. 2001, in 11th Cambridge Workshop on Cool Stars, Stellar Systems and the Sun, eds. R. J. Garcia Lopez, R. Rebolo, & M. R. Zapaterio Osorio, ASP Conf. Ser., 223, 878

* Strassmeier, K. G. 2009, A&ARv, 17, 251

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

* Tingley, B., Endl, M., Gazzano, J.-C., et al. 2011, A&A, 528, A97

* Trampedach, R., Asplund, M., Collet, R., Nordlund, Å., & Stein, R. F. 2013, ApJ, 769, 18

* van Belle, G. T. 2008, PASP, 120, 617

* Wright, J. T., Fakhouri, O., Marcy, G. W., et al. 2011, PASP, 123, 412

* Zhao, M., Monnier, J. D., ten Brummelaar, T., Pedretti, E., & Thureau, N. D. 2008, in SPIE Conf. Ser., 7013

* Zhao, M., Monnier, J. D., Che, X., et al. 2010, in SPIE Conf. Ser., 7734

* Zhao, M., Monnier, J. D., Che, X., et al. 2011, PASP, 123, 964