# Spectrally resolved interferometric observations of α Cephei and physical modeling of fast rotating stars

O. Delaa, J. Zorec, A. Domiciano de Souza, D. Mourard, K. Perraut, Ph.

Stee, Y. Frémat, J. Monnier, S. Kraus, X. Che, et al.

► To cite this version:

O. Delaa, J. Zorec, A. Domiciano de Souza, D. Mourard, K. Perraut, et al. Spectrally resolved interferometric observations of α Cephei and physical modeling of fast rotating stars. Astronomy and Astrophysics - A&A, 2013, 555, pp.A100. 10.1051/0004-6361/201220689. hal-03473885

HAL Id: hal-03473885

https://hal.science/hal-03473885

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
# Spectrally resolved interferometric observations of α Cephei and physical modeling of fast rotating stars

O. Delaa<fnref n="1" />,<fnref n="2" />, J. Zorec<fnref n="2" />, A. Domiciano de Souza<fnref n="1" />, D. Mourard<fnref n="1" />, K. Perraut<fnref n="3" />, Ph. Stee<fnref n="1" />, Y. Frémat<fnref n="4" />, J. Monnier<fnref n="7" />,
S. Kraus<fnref n="7" />, X. Che<fnref n="7" />, Ph. Bério<fnref n="1" />, D. Bonneau<fnref n="1" />, J. M. Clausse<fnref n="1" />, M. Challouf<fnref n="1" />, R. Ligi<fnref n="1" />, A. Meilland<fnref n="1" />, N. Nardetto<fnref n="1" />,
A. Spang<fnref n="1" />, H. McAlister<fnref n="5" />,<fnref n="6" />, T. ten Brummelaar<fnref n="6" />, J. Sturmann<fnref n="6" />, L. Sturmann<fnref n="6" />, N. Turner<fnref n="6" />,
C. Farrington<fnref n="6" />, and P. J. Goldfinger<fnref n="6" />

Received 2 November 2012 / Accepted 25 February 2013

## ABSTRACT

**Context.** When a given observational quantity depends on several stellar physical parameters, it is generally very difficult to obtain observational constraints for each of them individually. Therefore, we studied under which conditions constraints for some individual parameters can be achieved for fast rotators, knowing that their geometry is modified by the rapid rotation which causes a non-uniform surface brightness distribution.

**Aims.** We aim to study the sensitivity of interferometric observables on the position angle of the rotation axis (PA) of a rapidly rotating star, and whether other physical parameters can influence the determination of PA, and also the influence of the surface differential rotation on the determination of the $\beta$ exponent in the gravity darkening law that enters the interpretation of interferometric observations, using $\alpha$ Cep as a test star.

**Methods.** We used differential phases obtained from observations carried out in the $H\alpha$ absorption line of $\alpha$ Cep with the VEGA/CHARA interferometer at high spectral resolution, $R = 30\,000$ to study the kinematics in the atmosphere of the star.

**Results.** We studied the influence of the gravity darkening effect (GDE) on the determination of the PA of the rotation axis of $\alpha$ Cep and determined its value, $PA = -157^{+17}_{-10}^{\circ}$. We conclude that the GDE has a weak influence on the dispersed phases. We showed that the surface differential rotation can have a rather strong influence on the determination of the gravity darkening exponent. A new method of determining the inclination angle of the stellar rotational axis is suggested. We conclude that differential phases obtained with spectro-interferometry carried out on the $H\alpha$ line can in principle lead to an estimate of the stellar inclination angle $i$. However, to determine both $i$ and the differential rotation parameter $\alpha$, lines free from the Stark effect and that have collision-dominated source functions are to be preferred.

**Key words.** techniques: high angular resolution – techniques: interferometric – stars: individual: $\alpha$ Cephei – stars: rotation – stars: individual: Aderamin

## 1. Introduction

The centrifugal force of rapid rotators induces a strong deformation of the stellar geometry that has been nicely shown observationally, thanks to long-baseline optical interferometrists that today allow observations at high spatial resolution to be performed (van Belle et al. 2001, 2006; Domiciano de Souza et al. 2003; McAlister et al. 2005; Kervella & Domiciano de Souza 2006; Monnier et al. 2007; Zhao et al. 2009). Concomitant to the stellar flattening, there is a non-uniform effective temperature distribution, called the gravity darkening effect (GDE; von Zeipel 1924). For stars with conservative rotation laws (rigid rotation is a particular case of a conservative law), whose envelopes are in hydrostatic and radiative equilibrium and where the diffusion approximation to the radiative transfer equation is used, the

following relation for the effective temperature distribution as a function of the co-latitude $\theta$ holds:

$$T_{\text{eff}}(\theta) = T_{\text{p}} \left[ \frac{g_{\text{eff}}(\theta)}{g_{\text{p}}} \right]^{\beta}, \quad (1)$$

where $T_{\text{eff}}(\theta)$ and $g_{\text{eff}}(\theta)$ are the co-latitude $\theta$ dependent effective temperature and effective gravity $T_{\text{p}}$ is the polar effective temperature, and $g_{\text{p}}$ is the polar gravity. The canonical value $\beta = 0.25$ (i.e., the standard gravity-darkening law von Zeipel 1924) is thus valid only in the frame of the above-mentioned physical conditions and approximations. The value $\beta = 0.08$ holds for stars with convective envelopes (Lucy 1967; Reiners 2003). Tuominen (1972) also found that $\beta = 0.15$ for non-barotropic differentially rotating stars with masses $M \approx 2.5 M_{\odot}$. In a recent work, Zorec et al. (2011, and references therein)

<footnote n="1">¹ Lab. Lagrange, CNRS UMR 6525, Univ. de Nice-Sophia Antipolis, Obs. de la Côte d'Azur, Avenue Copernic, 06130 Grasse, France e-mail: Omar.Delaa@obs-azur.fr</footnote>

<footnote n="2">² Institut d'Astrophysique de Paris, UMR 7095 du CNRS, Université Pierre & Marie Curie, 98bis bd. Arago, 75014 Paris, France</footnote>

<footnote n="3">³ UJF - Grenoble I/CNRS - INSU, Institut de Planétologie et d'Astrophysique de Grenoble (IPAG) UMR 5274, 38041 Grenoble, France</footnote>

<footnote n="4">⁴ Royal Observatory of Belgium, 3 Av. Circulaire, 1180 Brussels, Belgium</footnote>

<footnote n="5">⁵ Georgia State University, PO Box 3969, Atlanta, GA 30302-3969, USA</footnote>

<footnote n="6">⁶ CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA</footnote>

<footnote n="7">⁷ University of Michigan, Astronomy Department, 941 Dennison Bldg, Ann Arbor, MI 48109-1090, USA</footnote>
commented that $\beta < 0.25$ may reveal a non-conservative rotational law in the stellar surface. Nevertheless, when more detailed expressions for the radiative flux are used and not simply the diffusion approximation, $\beta \le 0.25$ (Claret 2012). Che et al. (2011) suggested $\beta = 0.19$ for rapid rotators, and Espinosa Lara & Rieutord (2012) showed that $\beta$ ranges from 0.25 for slow rotators and decreases for rapid rotators, compatible with observations of Che et al. (2011).

Rotation generates hydrodynamical instabilities that redistribute the angular momentum in the stellar interior. Some of them were taken into account to study their effect on the stellar structure and evolution with rotation (cf. Endal & Sofia 1979; Pinsonneault et al. 1991; Heger & Langer 2000; Meynet & Maeder 2000). These instabilities also induce a turbulent diffusion (Zahn 1992) that, in massive objects, distributes the CNO-cycled material from the core to the envelope and changes the stellar atmospheric chemical composition profile (Meynet & Maeder 2000). Owing to the mass-compensation effect related to the rapid rotation (Sackmann 1970), the total radiative flux produced in the stellar core is lower than that of a parent rotationless star. The combined effect of the mass-compensation effect and the refueling of the stellar core by the meridian circulation contributes to extend the lifespan of stars in the main sequence (cf. Kiziloglu & Civelek 1996; Talon et al. 1997; Meynet & Maeder 2000). Stellar winds and the related continuum mass-loss rate are somewhat modified by the rotation (Lamers & Cassinelli 1999; Maeder et al. 2007). Maeder et al. (2008) have shown that in massive stars, rapid rotation considerably enlarges the Fe and He ionization-related convective zones in their envelope. This convection can then induce a non-conservative, albeit non-shellular, law of differential rotation in the envelope, on which rely both the stellar external geometry and its surface brightness distribution (Zorec et al. 2011).

This account of observational phenomena related with rotation shows that all stellar observables, in particular those we can obtain with interferometry, can be complex functions of several fundamental parameters, some of intrinsic stellar nature (mass, age, chemical composition, amount of angular momentum stored in the star, rotation law, etc.) and others of more geometrical or extrinsic character (inclination angle, position in the sky, etc.). The stellar geometry and the wavelength-dependent surface brightness distribution are both measurable with interferometric techniques. We can then ask whether by using observables, some of the quoted independent parameters can be tested individually, or whether they should all be determined and analyzed at once.

At a first glance, we may expect that due to their extrinsic nature, the position angle of the rotation axis and the stellar inclination can be determined independently of the stellar parameters of intrinsic character, even though the determination of the inclination angle relies on the surface brightness asymmetries, which in turn are controlled by intrinsic stellar quantities.

The aim of the present paper is the following: thanks to a spectral resolution that can reach 30 000, the VEGA/CHARA instrument (Mourard et al. 2009, 2010) offers a great possibility to perform precise measurements of differential phases of fringes to study the prevailing kinematics in the stellar atmospheres. By taking advantage of the link that exists between the phase of fringes and the position of the photo-center in a partially resolved star, we can carry out spectro-interferometric studies and provide measurements of the position angle (PA) of stars. A similar approach was recently used to measure the PA of the rotation axis of the fast rotators Fomalhaut (Le Bouquin et al. 2009) and Achernar (Domiciano de Souza et al. 2012a). In addition, this

work differs from previous studies of the PA (van Belle et al. 2006; Zhao et al. 2009), because we use information on the stellar atmosphere from the observed spectral lines. The method is tested for $\alpha$ Cep using the $H\alpha$ absorption line.

In principle, the combination of the high spectral and spatial resolution should enable us to extract information from the surface radial iso-velocity curves, whose properties reflect the kinematics in the stellar atmosphere. The apparent aspect of the radial iso-velocity curves depends, however, on the inclination angle $i$ of the rotation axis, and on the degree of the surface differential rotation. It is then worthwhile to discuss whether we can estimate the inclination angle $i$ and the degree of differential rotation independently of the intrinsic stellar parameters.

The present paper is organized as follows: in Sect. 2 we summarize previous studies of $\alpha$ Cep. In Sect. 3, we present our interferometric observations and reduction technique. In Sect. 4, we introduce our method for deriving the position angle PA of the stellar rotation axis. In this section, we also introduce the model used to interpret the interferometric measurements and present the results obtained with our best model. In Sect. 5 we study the effects on the von Zeipel coefficient $\beta$ determination that depend on the law of surface differential rotation and analyze the feasibility of the inclination angle $i$ determination. The conclusions are summarized in Sect. 6.

## 2. $\alpha$ Cephei

$\alpha$ Cephei (Alderamin, HD 203 280, HR 8162, $V = 2.46$, $d = 14$ pc) is the brightest star in the constellation of Cepheus. Originally classified as A2n (Douglas 1926), the spectral type of $\alpha$ Cephei has been re-estimated in Johnson & Morgan (1953), who classified it as A7IV-V. Gray et al. (2003) classified the star as A8V. Like Altair, $\alpha$ Cep has a chromosphere (Walter et al. 1995; Simon & Landsman 1997; Simon et al. 2002). According to various authors, the $V\sin i$ of $\alpha$ Cep ranges from $180 \text{ km s}^{-1}$ to $245 \text{ km s}^{-1}$ (Bernacca & Perinotto 1973a,b,c; Uesugi & Fukuda 1970; Royer et al. 2007; Abt & Morrell 1995).

Up to now, few interferometric observations have been carried out on $\alpha$ Cep. van Belle et al. (2006) used the Georgia State University CHARA array to observe this star in the infrared and estimated its oblateness, the rotational velocity, and its gravity darkening. Their data showed a non-circular projected disk brightness distribution from which they deduced the equatorial and polar diameter of $\alpha$ Cep using two different models. The first one, purely geometrical, represents an elliptical star whose brightness distribution is uniform. From this model, they found a polar diameter of $1.355 \pm 0.080$ mas, an equatorial diameter of $1.625 \pm 0.056$ mas and a position angle of the rotation axis of $3^\circ \pm 14^\circ$ (or $-177^\circ$, since there can be an indetermination of $180^\circ$) with a reduced $\chi^2$ of 1.08. The second, and more physical model, takes into account the shape of the photosphere under conditions of hydrostatic equilibrium, uniform rotation, Roche's approximation for the total potential at the surface, the limb darkening that responds to the geometrical deformation of the star, and the surface gravity darkening. The parameters derived from their best model are summarized in Table 1. Their estimate of the $\beta$ exponent in the gravitational darkening law gives $\beta = 0.084^{+0.026}_{-0.049}$, which is consistent with $\beta = 0.08$ expected for a star whose envelope is entirely convective (Lucy 1967; Reiners 2003).

Taking advantage of a better $(u, v)$-plane coverage than in van Belle et al. (2006), Zhao et al. (2009) used the Michigan Infra-Red Combiner (MIRC) at the CHARA Array to produce a model for $\alpha$ Cep as fast rotator. These authors calculated models of brightness distribution in fast rotating stars according to
<table><caption>**Table 1.** Parameters from best models obtained by van Belle et al. (2006) and Zhao et al. (2009).</caption><thead><tr><th>Model parameters</th><th>Van Belle et al. (2006)</th><th>Zhao et al. (2009)</th></tr></thead><tbody><tr><td><i>i</i></td><td>88.2°<sup>+1.8°</sup><sub>-13.3°</sub></td><td>55.70° ± 6.23°</td></tr><tr><td>PA</td><td>-163°<sup>+3.2°</sup><sub>-4.3°</sub></td><td>-178.84° ± 4.28°</td></tr><tr><td><i>T</i><sub>p</sub></td><td>8440<sup>+430</sup><sub>-700</sub> K</td><td>8588 ± 300 K</td></tr><tr><td><i>R</i><sub>p</sub>/<i>R</i><sub>☉</sub></td><td>2.175 ± 0.046</td><td>2.162 ± 0.036</td></tr><tr><td><i>T</i><sub>eq</sub></td><td>7600 K</td><td>6574 ± 200 K</td></tr><tr><td><i>R</i><sub>eq</sub>/<i>R</i><sub>☉</sub></td><td>2.823 ± 0.097</td><td>2.740 ± 0.044</td></tr><tr><td>Ω/Ω<sub>c</sub></td><td>0.8287<sup>+0.0482</sup><sub>-0.0232</sub></td><td>0.941 ± 0.020</td></tr><tr><td>β</td><td>0.084<sup>+0.026</sup><sub>-0.049</sub></td><td>0.216 ± 0.021</td></tr><tr><td>Vsin<i>i</i></td><td>283 ± 19 km s<sup>-1</sup></td><td>225 km s<sup>-1</sup></td></tr><tr><td>χ<sup>2</sup><sub>red</sub></td><td>2.16</td><td>1.18</td></tr></tbody></table>

prescriptions given by Aufdenberg et al. (2006, and references therein), which they fitted to the interferometric observations. In these models the authors assumed solid-body rotation, adopted a Roche potential to characterize the stellar surface, and used the von Zeipel gravity darkening law (von Zeipel 1924) to determine the latitude-dependent effective temperature profile. Hence, for α Cep they deduced the inclination angle *i*, the position angle PA of the rotation axis, the polar *R*<sub>pol</sub> and equatorial *R*<sub>eq</sub> radii, the polar effective temperature *T*<sub>p</sub>, and equatorial effective temperature *T*<sub>eq</sub>, the angular velocity ratio *ω* = Ω/Ω<sub>c</sub>, and the β exponent. These values are summarized in Table 1.

The exponent $\beta = 0.216 \pm 0.021$ found by Zhao et al. (2009) suggests that the envelope of α Cep is probably not convective as found by van Belle et al. (2006). Because their β exponent is lower than 0.25, at least three possibilities exist: 1) the stellar envelope is mostly radiative but interspersed with convective regions, which in the equatorial region case perhaps induced by the GDE where the effective temperature is low enough; 2) the surface rotation reflects a deeper non-conservative rotation law in the envelope; or 3) the diffusion approximation used to derive the gravity darkening law is not sufficient, and more detailed expressions for the emerging radiative flux must be taken into account (Claret 2012).

# 3. Observations and data reduction

## 3.1. Observations

The VEGA instrument (Mourard et al. 2009) at the CHARA Array (ten Brummelaar et al. 2005), located on the Mount Wilson Observatory (LA, California, USA), operates in the visible domain and benefits from a spectrograph and a polarimeter. The spectrograph is designed to sample in the visible band from 0.45 to 0.85 µm and is equipped with two photon-counting detectors that observe in two different spectral bands (the detector that observes at the shortest wavelength is denoted “blue”, the other is called “red”). Interferometric observations of α Cep with the VEGA/CHARA instrument were carried out in July 2008 and October 2010 using the high spectral resolution mode (R = 30 000). All observations were made with two telescopes only. Observations were conducted with the red and blue detectors centered at 656 nm (Hα line) and 632 nm, respectively. In this spectral resolution mode, the Hα line covers almost the entire surface of the red detector, which allows us to make a precise measurement of the differential phase. Because we based our study on a differential analysis, we did not need to observe

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** (u, v)-plane coverage by our observations (u is in abscissas).</figcaption></figure>

<table><caption>**Table 2.** Log of observations with VEGA/CHARA.</caption><thead><tr><th>Date (UTC)</th><th>Telescopes</th><th>Projected baseline Length (m)</th><th>PA (in °)</th><th>r<sub>0</sub> (cm)</th></tr></thead><tbody><tr><td>28/07/2008</td><td>03h44 S1S2<sub>a</sub></td><td>28.4</td><td>+29.3°</td><td>8</td></tr><tr><td>28/07/2008</td><td>09h26 S1S2<sub>b</sub></td><td>29.3</td><td>-20.1°</td><td>13</td></tr><tr><td>28/07/2008</td><td>11h48 S1S2<sub>c</sub></td><td>26.4</td><td>-42.1°</td><td>8</td></tr><tr><td>13/10/2010</td><td>03h02 E1E2<sub>a</sub></td><td>63.6</td><td>-116.4°</td><td>8</td></tr><tr><td>13/10/2010</td><td>03h52 E1E2<sub>b</sub></td><td>65.8</td><td>-125.7°</td><td>8</td></tr><tr><td>13/10/2010</td><td>04h59 E1E2<sub>c</sub></td><td>65.7</td><td>-145.0°</td><td>7</td></tr></tbody></table>

Notes. $r_0$ is the Fried turbulence parameter; sub-indices *a*, *b* and *c* identify the observing dates.

calibrators. The log of the observations is given in Table 2 and the (u, v)-plane coverage during these observation campaigns is shown in Fig. 1.

## 3.2. Data reduction

The data were reduced with the software developed in the VEGA group and described by Mourard et al. (2009). This software offers two data reduction modes, the first one is the spectral density analysis and the second one the cross-spectral analysis scheme. In the spectral density mode, we only obtain the squared visibility averaged over a spectral band centered at a given wavelength $\lambda_1$. With the cross-spectral analysis mode, we obtain a visibility modulus $V(\lambda)$ and a differential phase $\phi_{\text{diff}}$ within a narrow channel of a few angstroms, which is step-by-step translated over a large reference spectral band. In each step, we measure the product $V(\lambda_1) \times V(\lambda_2)$ and $\phi_{\text{diff}}(\lambda_1)$. In this study, we are only interested by the quantity $\phi_{\text{diff}}(\lambda_1)$. The spectral calibration was performed with known spectral lines of the ThAr spectral lamp. A polynomial fit was made to determine the dispersion law. The residual of the fit is typically two times lower than the spectral resolution.

Figure 2 presents the differential phases obtained after the data processing. In this figure we also superimposed (red solid line) the best-fit model. As explained in Sect. 4.4, the differential phases in the wings of the Hα line (outside 6555 Å to 6569 Å) can be considered to be equal to zero. We used the measurements
<figure><img src="image_3.png" /><figcaption>Fig. 2. Differential phases obtained after data processing. This figure also shows the best-fit (red line) models calculated as a function of the position angle of the rotation axis (PA = -157°) obtained as explained in Sect. 4.4.</figcaption></figure>

in this “wings” region to estimate the variance of the phase and decided to attribute this variance to the individual phase error. We then applied a fit of the model in the S-shaped region, where a signal is expected, i.e. from 6555 Å to 6569 Å.

# 4. Determination of the position angle

## 4.1. Principles

Interferometric observables were modeled from a computed brightness distribution (see Sect. 4.2) using the Zernike-van Cittert’s theorem. This theorem relates the complex visibility $Y_{\delta\lambda}$ with a stellar surface brightness distribution emitted in a given spectral channel of bandwidth $\delta\lambda$:

$$Y_{\delta\lambda}(u, v) = \frac{\int_{-\infty}^{\infty} \int_{-\infty}^{\infty} I_{\delta\lambda}(\alpha, \beta) \exp[-2\pi i(\alpha u + \beta u)] d\alpha d\beta}{\int_{-\infty}^{\infty} \int_{-\infty}^{\infty} I_{\delta\lambda}(\alpha, \beta) d\alpha d\beta}, \quad (2)$$

where $I_{\delta\lambda}(\alpha, \beta)$ is the intensity of the emitted energy by the source toward the observer caught in the spectral channel of bandwidth $\delta\lambda$ and centered on the wavelength $\lambda$; $(\alpha, \beta)$ are the angular coordinates of an emitting point in the source projected onto the sky, while $(u, v)$ are the Fourier conjugate spatial frequencies related to the baseline coordinates and measured in the observing plane. The modulus of the complex visibility, $V_{\delta\lambda} = |Y_{\delta\lambda}(u, v)|$, measures the contrast of the interferometric fringes, while its argument determines the phase of the fringes, $\phi_{\delta\lambda} = \arg[Y_{\delta\lambda}(u, v)]$. For an unresolved source, $\phi_{\delta\lambda}$ is directly related to the position of the photo-center of the brightness distribution at a given wavelength $\lambda$. Unfortunately, the atmospheric turbulence disrupts the phase of fringes when interferometric observations are performed with only two telescopes. To overcome this difficulty, we used the differential phase $\phi_{\text{diff}}$, defined as the difference between the phase measured in two spectral channels centered on the wavelengths $\lambda_1$ and $\lambda_2$:

$$\phi_{\text{diff}} = \phi_{\lambda_2} - \phi_{\lambda_1}. \quad (3)$$

The quantity $\phi_{\text{diff}}$ enables us to quantify the displacement of photo-centers at different wavelengths $\lambda$ of an unresolved source free from perturbations from the atmospheric turbulence. The spatial resolution reached with the S1S2 and E1E2 baselines used in this study are 5.5 and 2.7 mas, respectively, at 0.656 µm.

Hence, $\alpha$ Cep remained unresolved (see Table 1). The cross-spectrum method described in Berio et al. (1999) and Tatulli et al. (2007) used in the VEGA/CHARA data processing aims at measuring the differential phases $\phi_{\text{diff}}(\lambda)$ and the differential visibility $V_{\text{diff}}(\lambda)$ between two different spectral channels centered on the wavelengths $\lambda_1$ and $\lambda_2$.

## 4.2. Models of stellar surface brightness distribution

To calculate $\phi_{\text{diff}}$ and $V_{\text{diff}}(\lambda)$ we used the model described in Domiciano de Souza et al. (2002) and Domiciano de Souza et al. (2012b), which can produce the intensity maps in the $H\alpha$ line needed in this work. This model is based on the following assumptions: a) stars have rigid rotation; b) Roche-approximation is used to describe the stellar shape modified by the rigid rotation; c) at each point of the observed stellar hemisphere, the atmosphere is approximated by a plane-parallel model for the local effective temperature $T_{\text{eff}}(\theta)$ and gravity $g(\theta)$.

A modified version of the BRUCE/KYLIE calculation code (Townsend 1997) was used to divide the apparent hemisphere into elementary surface elements (blocks), each characterized by a given set of local physical parameters. The plane-parallel stellar atmosphere models were computed using the Kurucz stellar atmosphere models with solar abundances (Kurucz 1979). The SYNSPEC code (Hubeny & Lanz 1995) with chemical compositions compatible with metallicity $Z = 0.02$ was employed to calculate the local emergent specific radiation intensities directed toward the observer $I_{\delta\lambda}(\alpha, \beta)$ which determine the monochromatic intensity maps of the apparent photosphere in the wavelengths of the instrumental transmitting channels. The models depend on five intrinsic stellar parameters: polar radius $R_p$, the polar temperature $T_p$, the polar gravity $g_p$, the stellar rotational flattening $D = R_c/R_p$, and the gravity darkening exponent $\beta$. In addition, there are two extrinsic parameters: the stellar inclination angle $i$ and the position angle PA of the stellar rotation axis.

## 4.3. Method

We adopted a Cartesian reference system in the sky centered on the star, so that the x-axis lies along the equator of the star and the y-axis contains the plane of the stellar rotation axis. In the upper panel of Fig. 3, we show the computed displacement of the photo-center of the $H\alpha$ line as observed along the equatorial axis of $\alpha$ Cep. The corresponding differential phases $\phi_{\text{diff}}(\lambda)$ are shown in the lower panel of Fig. 3. Phases are calculated for a sixty-meter baseline parallel to the direction of the stellar equatorial axis. When the orientation of the baseline projected onto the sky changes with respect to the x-axis, we obtain different plots of $\phi_{\text{diff}}(\lambda)$, as illustrated in Fig. 4, where the line-styles change according to the orientation of the baseline. $\phi_{\text{diff}}(\lambda)$ is well spread for a baseline perpendicular to the rotational axis of the star (Mourard 1990), whereas for a baseline parallel to the rotational axis of the star, $\phi_{\text{diff}}(\lambda)$ is close to zero. So, to determine the position angle PA of the stellar rotation axis of $\alpha$ Cep, we carried out spectro-interferometric observations at different orientations of the baseline and took measurements of the respective variation of phases $\phi_{\text{diff}}(\lambda)$. Thus, for a given observed distribution $\phi_{\text{diff}}(\lambda)$ as a function of $\lambda$, we can estimate the PA by searching the best fit of the observed differential phases with calculated ones for different position angles.
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Top: intensity maps calculated for selected wavelengths inside the Hα line (from left to right): 6559.5 Å, 6561.3 Å, 6562.8 Å, and 6565.2 Å. Ordinates and abscissas are given in units of the equatorial radius R<sub>eq</sub>. The inclination angle used is i = 55°. Bottom: photo-center displacements along the y-axis (**left side**); differential phase associated with the photo-center shift (**right side**). All quantities are computed for a sixty-meter baseline parallel to the equatorial axis (**x-axis**). The photo-center displacements are given in units of the angular stellar radius.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Differential phases $\phi_{\text{diff}}(\lambda)$ in the H$\alpha$ line for different baseline orientations. Top: stellar disk deformed by the fast rotation and projected onto the sky where we have represented the baseline orientations: 0° (cross), 45° (star), 90° (solid line) and 160° (diamond). Bottom: evolution of $\phi_{\text{diff}}(\lambda)$ according to the baseline orientations represented on top and computed for a sixty meter-baseline. In this example we assumed that the position angle of the stellar rotation axis is equal to zero.</figcaption></figure>

4.4. Model fitting

The position angle PA of the rotation axis is an extrinsic pa-
rameter whose determination can be assumed independent of

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Spectrum of alpha Cep obtained in the high spectral resolution mode of VEGA/CHARA (green) and i th the ELODIE spectrograph (red).</figcaption></figure>

the choice of the intrinsic stellar parameters. The results ob-
tained will justify this assumption a posteriori. So, to calculate
the brightness distribution, we adopted the stellar quantities R<sub>p</sub>,
T<sub>p</sub>, g<sub>p</sub>, D, and i obtained by Zhao et al. (2009, see Table 1) be-
cause of their excellent (u, v)-plane coverage. Then, we calcu-
lated monochromatic intensity maps for wavelengths inside the
Hα line as detailed in Sect. 4.2, and for a series of angle PA
and β parameters. The corresponding differential phases φ<sub>diff</sub>(λ)
were estimated by calculating the FT of the monochromatic in-
tensity maps.

It is important to note that in the high spectral resolution
mode, the red detector of the VEGA/CHARA instrument cov-
ers only the spectral interval λλ6542.49–6582.68 Å of the ob-
served Hα line of α Cep, as shown in Fig. 5. Because of the
lack of measurements in the continuum spectrum around the
line, we calibrated the differential phases by assuming that both
line wings in the uncovered spectral ranges are symmetric. To
test this hypothesis, we computed images of the apparent stellar
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Surface brightness distribution in $\alpha$ Cep at four different wavelengths: 6548.7 Å (top, left), 6557.7 Å (top, right), $\lambda6562.8$ Å (bottom, left), and $\lambda6568.5$ Å, (bottom, right). Ordinates and abscissas are given in units of the equatorial radius $R_{\text{eq}}$. The inclination angle used is $i = 55^{\circ}$.</figcaption></figure>

<table><caption>**Table 3.** Differential phases as a function of different wavelengths inside the $H_{\alpha}$ line for the simulated intensity maps of Fig. 6.</caption><thead><tr><td>$\lambda$ (Å)</td><td>$\phi_{\text{diff}}(\lambda)$</td></tr></thead><tbody><tr><td>6548.7</td><td>3°</td></tr><tr><td>6557.7</td><td>4°</td></tr><tr><td>6562.8</td><td>0°</td></tr><tr><td>6568.5</td><td>4°</td></tr></tbody></table>

hemisphere and the respective differential phases at four different wavelengths: $\lambda$ 6548.7, 6557.7, 6562.8, and 6568.5 Å. A perfectly symmetric distribution of the stellar brightness in the given wavelengths should lead the corresponding photocenters to lie in the center of the disk and the respective differential phases to be zero. The images of the brightness distribution over the stellar disk at the chosen wavelength are shown in Fig. 6. The values of the differential phase for each image are summarized Table 3. Because the expected uncertainty on the estimate of differential phases is roughly 3°, we conclude that in the wavelength-intervals $\lambda\lambda6542.49-6557.7$ Å and $\lambda\lambda6568.5-6582.68$ Å our hypothesis of wing-symmetry is a reasonable first approximation.

According to the results shown in Table 3, we can assume that the differential phase in the wings of the $H\alpha$ line is equal to zero, because we were dealing with the continuum spectrum proper. The differential phases in these spectral ranges were fitted with a polynomial and then used to divide the observed differential phases to obtain the sought calibration.

To determine PA, a grid of models parameterized according to the exponent $\beta$, and the position angle of the rotation axis of $\alpha$ Cep was produced. The $\beta$ exponents used range from 0.08 to 0.25 by steps of 0.01, while the angle PA ranged from -130° to -180° by steps of 1°. This interval of PA was chosen according to previous indications on the PA value inferred in van Belle et al. (2006) and Zhao et al. (2009). Each model provides thus the observable $\phi_{\text{diff}}(\lambda)$, which is compared to the data and the fit controlled through a $\chi^2$ test. In Fig. 2 are

<figure><img src="image_8.png" /><figcaption>**Fig. 7.** $\chi^2_{\text{red}}$ as a function of the PA of the rotation axis of $\alpha$ Cep for different von Zeipel exponents ranging from $\beta = 0.08$ to $\beta = 0.25$.</figcaption></figure>

shown the respective best fits obtained (red lines) to the observed differential phases. We computed the reduced $\chi^2$ residuals (hereafter $\chi_{\text{red}}$) and plotted as a function of the PA angle, as shown in Fig. 7, where each curve corresponds to a different gravity darkening exponent. We note that the $\chi_{\text{red}}$ curves are calculated for $\lambda6560$ Å and each of them is an average over the six better fits obtained from the respective observed differential phases.

From the $\chi_{\text{red}}$ curves in Fig. 7 we can see that if the PA determination corresponds to the minimum value [$\chi^2(\beta)_{\text{min}}$], roughly the same angle PA will be obtained whatever the $\beta$ exponent. Because each $\beta$ value produces a different effective temperature contrast over the observed stellar hemisphere, with the consequent contrast of the emitted radiative fluxes, which mimic the effect of different sets of intrinsic stellar parameters, we conclude that the observed differential phases $\phi_{\text{diff}}(\lambda)$ against the position angle reflect the external stellar geometry quite independently of any explicit intrinsic physical parameter.

To avoid inconveniences introduced by possible correlations among data points, a reliable approach to the uncertainty affecting the estimate of the position angle PA can be inferred by using the random re-sampling of data or “bootstrap simulation” (Efron & Tibshirani 1986). We thus randomly re-sampled $B = 700$ times the six data sets of differential phases shown in Fig. 8 in the wavelength interval used to fit the model-differential phases. This leads to $B$ independent determinations of PA values shown in Fig. 8a against the bootstrap replication number. The position angles of Fig. 8a were calculated for the von Zeipel coefficient $\beta = 0.19$ that was suggested in Che et al. (2011) to analyze interferometric data. The position angles PA were then gathered into a histogram of frequencies where the class-step (bin-width) $\Delta PA = 2^\circ$ was chosen following the bin-width optimization method in Shimazaki & Shinomoto (2007). The histogram is shown in Fig. 8b, where we also show a smoothed version of the histogram obtained using the Gaussian-kernel-estimator method (Bowman & Azzalini 1997) (dashed curve). We can readily see that the estimate of PA with the highest probability is $PA = -157^\circ$, while more than 99% of estimates are in the interval from $PA = -167^\circ$ to $PA = -140^\circ$ ($PA = -157^{+17^\circ}_{-10^\circ}$). The whole fitting procedure was also repeated for $\beta = 0.08$ and $\beta = 0.25$. While the confidence interval remains in all cases the same, for $\beta = 0.08$ the maximum is slightly shifted to more negative values of PA and slightly to more positive PA for $\beta = 0.25$ than with $\beta = 0.19$, but the whole displacement implies $\Delta PA \approx 2^\circ$. The low value of $\Delta\beta$ demonstrates the high independence of the PA-angle determination on the chosen value of $\beta$.
<figure><img src="image_9.png" /><figcaption>Fig. 8. a) Simulated position angle PA against the number of “bootstrap” replications when $\beta = 0.19$. The colored line indicates the position angle PA of maximum probability. b) Frequency distribution of the simulated position angles PA (histogram and smoothed distribution). In the figure are also indicated the position angle PA of maximum probability and the 99% confidence interval.</figcaption></figure>

The contrast of the brightness distribution over the stellar surface was determined simultaneously: a) the rotation law, characterized in Sect. 5.1 by the parameter $\alpha$; b) the value of the $\beta$ exponent; and c) the equatorial angular velocity $\Omega_{eq}$. It can be expected that some degeneracy of $\alpha$, $\Omega_{eq}$, and $\beta$ may appear when they are determined together from the observed stellar surface brightness distributions. Accordingly, there may also appear some indetermination of the inclination angle, since the $V\sin i$, a key quantity for estimating $i$, depends on $V_{eq}$, which in turn is determined by $\Omega_{eq}$ and the true stellar equatorial radius $R_{eq}$. In this section we have taken a wide range of possible values of $\beta$, which can be thought of as blurring uncertainties on the surface brightness distribution carried by $\alpha$ and $\Omega_{eq}$. We can conclude that at least the PA angle can be safely determined even though other intrinsic stellar quantities are uncertain.

## 5. Inclination angle and surface differential rotation

In this section we explore whether the inclination angle $i$ can be estimated with interferometric observations and what effects are induced by the surface differential rotation that may hamper the angle $i$-determination.

Compared to the rigid rotation where the angular velocity $\Omega(\theta)$ is constant over the entire stellar surface, the differential rotation introduces a supplementary $\theta$-dependent centrifugal acceleration that changes the stellar geometry, its surface gravity, and the brightness distribution (Zorec et al. 2011). All these changes should be reflected in the differential phases of the spectral lines studied.

It is important to note that the gravity darkening law given by Eq. (1) is valid only when the proportionality between the local bolometric radiative flux and the effective gravity is constant over the stellar surface. As already noted in Sect. 1, this is possible only when 1) the rotation law in the stellar envelope is conservative (the angular velocity depends only on the distance to

the rotation axis); 2) the envelope is in radiative equilibrium; and 3) the diffusion approximation to the radiative transfer holds.

The aim of this section is to provide an exploratory study of the effects shown in spectral lines by the GDE and the differential rotation when stars are seen at different inclination angles. As this is mainly an exploratory exercise, depending on the needs, we proceed with different degrees of approximation that are detailed according to the effect or the spectral line studied.

### 5.1. The $\beta$ exponent

The surface stellar brightness distributions are strongly dependent on the adopted darkening law. In the present discussion we write this law as follows:

$$F_{\text{bol}}(\theta, \alpha) = c(\eta, \alpha) g_{\text{eff}}^{4\beta}(\theta, \alpha), \quad (4)$$

where $F_{\text{bol}}(\theta, \alpha) \propto T_{\text{eff}}^4(\theta)$ is the local bolometric radiative flux and $T_{\text{eff}}(\theta)$ is the local effective temperature; $c(\eta, \alpha)$ is a function that should normally depend on both $\eta$ and $\alpha$, where $\alpha$ is the constant that parameterizes the degree of differential rotation, but in general possible also on $\theta$. Here, $\eta$ is the relation between the centrifugal acceleration at the equator ($F_c$) and the gravitational acceleration ($F_g$) at the equator, and $\alpha$ parameterizes the degree of surface differential rotation. We can assume that $c$ is still a constant dependent only on $\eta$, but the effect of $\alpha$ passed entirely on $\beta$, which accordingly becomes $\beta = \beta(\alpha, \beta, \theta)$. To put forward more clearly the effect of $\alpha$, we adopt $\beta$ to be constant and think of it as an average quantity over the co-latitude $\theta$. To assess the influence of $\alpha$ on the stellar brightness distribution, we changed the surface effective gravity adopting the Maunder relation for the surface angular velocity

$$\Omega(\theta) = \Omega_{\text{eq}} (1 + \alpha \cos^2 \theta), \quad (5)$$

where $\Omega_{\text{eq}}$ is the angular velocity at the equator. For $\alpha > 0$ the angular velocity is accelerated toward the pole ($\partial\Omega/\partial\theta < 0$; $\alpha$ can has any positive value), while for $\alpha < 0$ angular velocity accelerates toward the equator ($\partial\Omega/\partial\theta > 0$; it must be $\alpha > -1.0$ to keep all stellar latitudes rotating in the same sense). The fit of the solar surface differential angular velocity $\Omega(\theta)$ law (Snodgrass 1984) with Eq. (5) produces $\alpha < -0.29$ (Zorec et al. 2011). The expression for the effective gravity $g_{\text{eff}}$ then becomes

$$\frac{g_{\text{eff}}(\theta)}{g_p} = \left(\frac{R_p}{R_{\text{eq}}}\right)^2 \frac{1}{r^2(\theta)} \left\{1 + \left[1 - \eta_{\alpha}(\theta)r^3(\theta)\right]^2 - 1\right\} \sin^2 \theta^{1/2} \\ \eta_{\alpha}(\theta) = \eta \times (1 + \alpha \cos^2 \theta)^2$$

with

$$\eta = \frac{F_c}{F_g} = \frac{\Omega_{\text{eq}}^2 R_{\text{eq}}^3}{GM_*} = \left(\frac{\Omega_{\text{eq}}}{\Omega_c}\right)^2 \left(\frac{R_{\text{eq}}}{R_c}\right)^3 = \left(\frac{V_{\text{eq}}}{V_c}\right)^2 \left(\frac{R_{\text{eq}}}{R_c}\right). \quad (6)$$

In Eq. (6), $r(\theta)$ is the radius-vector normalized to the equatorial radius $R_{\text{eq}}$ that describes the geometry of the stellar surface. In Eq. (6), $g_p$ is the polar gravity; $R_p$ and $R_c$ are the polar and equatorial critical radii, respectively; $M_*$ is the stellar mass; $V_{\text{eq}}$ and $V_c$ are the equatorial linear velocity at $\eta$ and critical velocity, respectively. The radius vector $r(\theta)$, radii $R_p$ and $R_{\text{eq}}$ must be calculated consistently as a function of $\eta$ and $\alpha$ as indicated in Zorec et al. (2011). The effect in Eq. (6) of the centrifugal acceleration as a function of the co-latitude $\theta$ is produced by $\eta_{\alpha}(\theta)$ and $r(\theta)$.

To see the influence of $\alpha$ on the stellar brightness distribution, we computed several maps with different values of $\eta$ and $\alpha$.
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Brightness distribution on the stellar disk in the continuum at $\lambda = 4437$ Å, near the conspicuous He I 4471 Å line, calculated for a star rotating at $\eta = 0.9$ and for different values of $\alpha$ (from left to right): -1.0, 0.0 and 1.5, seen at an inclination angle of 90°.</figcaption></figure>

<table><caption>**Table 4.** Apparent $\beta_{\text{app}}$ exponents as a function of the $\alpha$ and $\eta$ parameters.</caption><thead><tr><th>η</th><th>0.01</th><th>0.1</th><th>0.3</th><th>0.5</th><th>0.7</th><th>0.9</th></tr><tr><th>Ω/Ω<sub>c</sub></th><td>0.172</td><td>0.521</td><td>0.796</td><td>0.918</td><td>0.976</td><td>0.998</td></tr><tr><th>α</th><td></td><td></td><td></td><td></td><td></td><td></td></tr></thead><tbody><tr><th>-1.0</th><td>0.250</td><td>0.253</td><td>0.261</td><td>0.267</td><td>0.271</td><td>0.274</td></tr><tr><th>-0.5</th><td>0.250</td><td>0.252</td><td>0.256</td><td>0.260</td><td>0.262</td><td>0.263</td></tr><tr><th>0.0</th><td>0.250</td><td>0.250</td><td>0.250</td><td>0.250</td><td>0.250</td><td>0.250</td></tr><tr><th>0.5</th><td>0.250</td><td>0.247</td><td>0.242</td><td>0.239</td><td>0.236</td><td>0.234</td></tr><tr><th>1.0</th><td>0.250</td><td>0.243</td><td>0.233</td><td>0.226</td><td>0.222</td><td>0.215</td></tr><tr><th>1.5</th><td>0.250</td><td>0.239</td><td>0.223</td><td>0.214</td><td>0.208</td><td>0.204</td></tr></tbody></table>

The maps obtained for $\lambda = 4437$ Å, $\eta = 0.9$, $\beta = 0.25$, and for an object characterized by $T_p = 19000$ K, $\log g_p = 4.0$ seen at $i = 90^\circ$, and several values of $\alpha$ are shown in Fig. 9. In Fig. 9, it is clearly apparent that in rapid rotators, $\alpha$ has a strong effect on the surface brightness distribution. Thus, even with the value of $\beta = 0.25$ (adopted for the calculations), the distribution of the surface brightness can be mis-interpreted in terms of apparent parameters, $\beta_{\text{app}} \neq \beta$, if the original von Zeipel relation is used to account for the GDE. The apparent $\beta_{\text{app}}$ can be estimated by dividing the stellar surface into elementary elements labeled with sub-indices $(i, j)$ that indicate the respective co-latitude and longitude angles $(\theta_i, \phi_j)$. Under the above adopted simplifications for the gravity darkening law, we let the parameter $\beta$ vary over the stellar surface and adopted an average estimate of it, which is given by $\beta_{\text{app}} = \langle\beta_{i,j}\rangle$

$$ \beta_{\text{app}} = 0.25 \times \left\{ \frac{\sum_{i=1}^{n} \sum_{j=1}^{m} \{ \ln g_{i,j}^{(0)} + \ln [F_{i,j}^{(\alpha)} / F_{i,j}^{(0)}] \}}{\sum_{i=1}^{n} \sum_{j=1}^{m} \ln g_{i,j}^{(\alpha)}} \right\}. \quad (7) $$

where $g_{i,j}^{(0)}, g_{i,j}^{(\alpha)}$ are the effective gravities at $(\theta_i, \phi_j)$ for $\alpha = 0$ and $\alpha \neq 0$, respectively; while $F_{i,j}^{(0)}, F_{i,j}^{(\alpha)}$ are the corresponding bolometric radiative fluxes. Table 4 summarizes the values of the apparent $\beta_{\text{app}}$ calculated for the stellar example presented in Fig. 9, except that calculations were made for several values of $\eta$.

From Table 4 we can see that $\beta_{\text{app}} \ge 0.25$ when $\alpha \le 0.0$ and $\beta_{\text{app}} \le 0.25$ if $\alpha \ge 0.0$. The higher the value of $\eta$, the larger the differences between $\beta_{\text{app}}$ and $\beta$. These results show that the interpretation of $\beta_{\text{app}}$ deduced from observations deserves a thorough discussion before anything can be said regarding the law characterizing the stellar surface rotation, whether radiation or convection prevails on the stellar surface. All these calculations rely on the assumption that the stellar atmosphere is in radiative equilibrium and that the radiative flux responds to the temperature gradient as established by the diffusion approximation. Claret (2012) has shown, however, that this approximation may not hold. Under more detailed representations of the radiative flux we can also have a $\beta_{\text{app}} \le 0.25$ effect that certainly relates in some way to the changes foreseen above as a function of $\alpha$ and $\beta$.

## 5.2. The $H\alpha$ line as a function of the inclination angle $i$ and differential rotation parameter $\alpha$

Because we have at our disposal interferometric observations covering almost the entire $H\alpha$ line of $\alpha$ Cep, in the present section, we explore whether we can use this line to estimate the inclination angle $i$ of the rotation axis. Knowing that the inclination angle is in principle a parameter of extrinsic character, it may be interesting to inquire on how much its determination may depend on the simultaneous knowledge of other intrinsic stellar quantities. A thorough discussion of this question is far beyond the scope of the present paper. We shall, however, briefly skim its dependence on the inclination angle $i$ and on the differential rotation parameter $\alpha$.

First, we calculated the $H\alpha$ line emitted by the $i$-apparent stellar hemisphere by considering the stellar geometry controlled by a given $\eta$ and several values of the differential rotation parameter $\alpha$. Using the stellar geometry derived according to the rules given in Zorec et al. (2011), we assumed that each surface element is characterized by a plane-parallel stellar atmosphere corresponding to $\log g_{\text{eff}}(\theta)$, and $T_{\text{eff}}(\theta)$ derived from Eq. (4). All these requirements are taken into account in the calculation code FASTROT modified for surface differential rotation (Frémat et al. 2005; Zorec et al. 2011). Adopting the stellar fundamental parameters of $\alpha$ Cep obtained by Zhao et al. (2009), we calculated two series of $H\alpha$ line profiles assuming rigid rotation, i.e., $\alpha = 0.0$. In the first series, shown in Fig. 10, we adopted $V_{\text{eq}} = 267$ km s$^{-1} \approx [225 \pm 20]/\sin[55.7^\circ \pm 6.2^\circ]$ which implies $\eta \approx 0.57$, and varied the inclination angle in the interval $35^\circ < i < 90^\circ$. In this figure we clearly see the very low sensitivity of the $H\alpha$ line profile to the $V\sin i$ parameter. In the second series, also shown in Fig. 10, we have adopted $V\sin i = 225$ km s$^{-1}$ and varied the force-ratio in the interval $0.4 < \eta < 0.8$, which implies a change in the inclination angle from $46.7^\circ$ to $77.1^\circ$.

We see again that only the core of the line profile may show some significant changes as a function of $\eta$. For a third series of $H\alpha$ line profiles we assumed $\eta = 0.57$ for the rotation at the equator, so that $V\sin i = 225$ km s$^{-1}$ suited for $i \approx 56^\circ$ and adopted several values of the differential parameter in the interval $= -1.0 < \eta < +1.5$. The line profiles obtained are shown in Fig. 11, where we see that once a $V\sin i$ is adopted, the core of the line changes somewhat with $\alpha$.

In a second step, we calculated brightness distributions over the apparent stellar hemisphere in the wavelengths of the $H\alpha$ line by considering that the stellar geometry is controlled by a given $\eta$ and differential rotation parameters $\alpha$, seen at various inclination angles. The respective differential phases as a function of the wavelength $\lambda$ inside the $H\alpha$ line were calculated as before, with the physical parameters $T_p$, $\log g_p$, $\Omega/\Omega_c$, $R_{\text{eq}}/R_p$, $R_p/R_\odot$, and $\beta$ of $\alpha$ Cep found by Zhao et al. (2009) (see Table 1). We also considered the position angles or arrays during the observations carried out on $\alpha$ Cep from 2008 to 2010 and the corresponding baselines of VEGA/CHARA given in Table 2. For each baseline configuration and inclination angles $i$ ranging from $47.7^\circ$ to $75.7^\circ$ by a $3^\circ$ step, we obtained the differential phases shown in Fig. 12. If the method were sensitive enough
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Rotationally broadened Hα line profiles using the fundamental parameters of α Cep and assuming rigid rotation. a) Line profiles adopting η = 0.57 (or Veq = 267 km s⁻¹) and varying the inclination angle. For comparison, we have included the observed Hα line profile of α Cep. b) Line profiles adopting Veqsin i = 225 km s⁻¹ and varying η (or the inclination angle i).</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Table 5.** Strongest variation of φ<sub>diff</sub>(λ) across the Hα line of α Cep as a function inclination angles i ranging from 47.7° to i = -75.7°.</figcaption></figure>

<table><thead><tr><th>Array</th><th>Strongest variation of<br>φ<sub>diff</sub>(λ)</th></tr></thead><tbody><tr><td>S1S2<sub>a</sub></td><td>0.84°</td></tr><tr><td>S1S2<sub>b</sub></td><td>0.58°</td></tr><tr><td>S1S2<sub>c</sub></td><td>0.11°</td></tr><tr><td>E1E2<sub>a</sub></td><td>2.55°</td></tr><tr><td>E1E2<sub>b</sub></td><td>1.43°</td></tr><tr><td>E1E2<sub>c</sub></td><td>2.20°</td></tr></tbody></table>

<figure><img src="image_13.png" /><figcaption>**Fig. 11.** Rotationally broadened Hα line profiles using the fundamental parameters of α Cep and surface differential rotation. The line profiles are plotted for η = 0.57, Vsin i = 225 km s⁻¹ and for several values of α of the Maunder relation for surface differential rotation.</figcaption></figure>

to the inclination i, the interval chosen for the tested inclination angle i would, in principle, enable us to distinguish between the estimates of i by van Belle et al. (2006) and Zhao et al. (2009). The respective largest amplitudes of deviations of φ<sub>diff</sub>(λ) obtained are given in Table 5, which show that for the configurations of the interferometric array used, our observations of the Hα line can hardly lead to conclusive results.

In spite of these results, we repeated the calculation of differential phases for a CHARA Array optimal configuration,

represented by a baseline of 63 meters and position angle PA = 157°. Using the same stellar physical parameters and inclination angles 0° ≤ i ≤ 90° as before, we obtained the differential phases shown in Fig. 13. In the left panel are shown differential phases when Vsin i = 225 km s⁻¹ is conserved and the inclination angle varied from in the interval 35° < i < 90°. We can see in this figure that the maximum amplitude of the phase variation is slightly larger than the uncertainty on the observational differential phase estimate (σ<sub>φ</sub> ≈ 3°). In the left panel, models keep V(=267 km s⁻¹) constant but vary the inclination angle from i = 0° to i = 90°. Since this is another diagnostic way to interpret the observed differential phases, we see that the method leads to amplitude changes of phases that are larger than Δφ = 5 × σ<sub>φ</sub>. We also note that changes are also seen in the position of maxima larger than 10 km s⁻¹, which provide a supplementary tool for interpreting the observations. We may then expect that observations carried on the Hα line might lead us to some estimate of the inclination angle i if they are performed at optimal conditions, which unfortunately were not fulfilled in the present case. Finally, we note that model calculations should be made with fundamental parameters corresponding to the parent, non-rotating object. This implies that the estimate of i, α, and of the remaining stellar fundamental parameters have to be iterated. These problems will be discussed in greater detail in a forthcoming paper.
<figure><img src="image_14.png" /><figcaption>**Fig. 12.** Estimated variations of the differential phases as function of the inclination (from 47.7° to 75.7°) of α Cep across the Hα line for our six observations. Calculations were performed with PA = −158°.</figcaption></figure>

### 5.3. Properties of spectral lines and the determination of *i* and *α*

To understand better the lack of results on *i* and *α* obtained for α Cep from studies based on the Hα line, we considered the rotational broadening and differential phases as generalized convolutions between the line flux profile and a rotational broadening-like (RB) function. To this end, we note that the width of the rotational broadening function proper is about Δλ<sub>RB</sub> = λ<sub>Hα</sub>Vsin<i>i</i>/c = 5.6 Å, while the Hα line spreads over some Δλ<sub>Hα</sub> ≈ 80 Å. This means that the several intervening rotational broadening functions used should also behave like quasi-delta functions that nearly preserve the shape of the original spectral line profile, except perhaps in the line core wavelengths. Because the broad width of the Hα line is due to the Stark effect, to obtain more detailed information on the rotational broadening, in particular that produced by a surface differential rotation, spectral lines as free as possible from Stark broadening should be employed. The Hα line wings are sensitive to log *g*, so a variation of the line profile in rapidly rotating stars induced by the θ-dependent log *g*<sub>eff</sub> is expected. Nevertheless, the rapidly rotating equatorial regions where both log *g*<sub>eff</sub> and the effective temperature are low do not contribute efficiently to the global line flux. The Doppler broadening of the line is thus left to the higher stellar latitudes where log *g*<sub>eff</sub> is comparatively larger and more uniform.

Finally, the nature of the Hα line source function also deserves some attention. In stars with T<sub>eff</sub> ≳ 8000 K, the Hα-source

function is dominated by photo-ionizations, which make the line quite non-local in character (Thomas 1965, 1983; Mihalas 1978). Consequently, the Hα line formation is fairly independent of details concerning the non-uniform distributions of T<sub>eff</sub>(θ) and log *g*<sub>eff</sub>(θ) induced by the GDE. Because we expect that the surface differential rotation is correlated with physical condition changes in the line formation regions, to increase the degree of information on the external rotation law of stars, we should use spectral lines that not only carry signatures of the macroscopic broadening velocity fields, but also have source functions sensitive to changes of local physical conditions (T<sub>e</sub>, N<sub>e</sub>). This means that these spectral lines should have source functions that are collision-dominated.

### 5.4. Information expected from spectral lines sensitive to effects induced by the rotation

According to spectral type and luminosity class, several spectral lines can be good candidates to carry out such studies. Since the identification of these lines and the subsequent calculation of expected effects from all of them is far beyond the scope of the present paper, we limit our discussion to a hypothetical line that behaves like the He Iλ 4471 line, which fulfills the above quoted requirements in stars with T<sub>eff</sub> ≳ 10 000. Though this line is not useful for α Cep, our aim is to briefly sketch some expected effects of quite general character.

Because the He Iλ 4471 line is only a hypothetical model line, we do not proceed with detailed radiation transfer
<figure><img src="image_15.png" /><figcaption>**Fig. 13.** Differential phases $\phi_{\text{diff}}(\lambda)$ over the H$\alpha$ line for different inclination angles $i$ of $\alpha$ Cep. Computations were done for a baseline of 63 m and baseline orientation -158°. Left panel: differential phases for $V\sin i = 225$ km s$^{-1}$ fixed, but varying inclination angle from $i = 35^\circ$ to $i = 90^\circ$ by steps $\Delta i = 5^\circ$ and $\eta$ accordingly. Right panel: differential phases for $V = 267$ km s$^{-1}$ fixed ($\eta = 0.57$) and inclination angles varying from $i = 10^\circ$ to $i = 90^\circ$ by steps $\Delta i = 5^\circ$. The parent star without rotation has $T_{\text{eff}} = 7500$ K and $\log g = 4.0$.</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Fig. 14.** He I 4471-like line profiles as a function of different values of the differential rotation parameter $\alpha$. The star is characterized by $T_p = 19\,000$ K, $\log g_p = 4.0$, $\eta = 0.9$, $\beta = 0.25$ and $i = 45^\circ$.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Fig. 15.** Iso-radial velocity zone for $\lambda = 4438.78$ Å. Top, form left to right: $\alpha = -1.0$ and $-0.5$. Bottom, from left to right: $\alpha = 0.0$ and $1.0$. Computations were done for an inclination angle of $i = 45^\circ$.</figcaption></figure>

calculations, but represent the emergent monochromatic specific intensity $I(\lambda, \mu)$ directed toward the observer as follows:

$$I(\theta, \phi, \eta, \alpha, \lambda, \mu) = H \left[ \theta, \lambda - \lambda_o \frac{V_{\text{proj}}(\theta, \phi, \alpha, i)}{c} \right] I_c(\theta, \eta, \alpha, \mu), \quad (8)$$

where $\theta$ reflects the dependence of $I(\lambda, \mu)$ with the local $T_{\text{eff}}(\theta)$ and $\log g_{\text{eff}}(\theta)$, whose variation with $\theta$ is determined by $\eta$ and $\alpha$; the dependence of $I(\lambda, \mu)$ with the longitude $\phi$ is due to the law of the surface differential rotation that produces iso-radial velocity curves that are a function of $\theta$ and $\phi$ over the apparent stellar hemisphere; $\eta$ controls the equatorial rotational velocity, and $\alpha$ parameterizes the surface differential rotation according to the Maunder law; $I_c(\theta, \eta, \alpha, \mu)$ is the specific radiation intensity directed toward the observer in the continuum at the central wavelength $\lambda_c$ of the studied line; $\mu(\theta, \phi, i) = \cos(\hat{n}, \hat{\beta}) = \hat{n} \cdot \hat{\beta}$, where $\hat{n} = \bar{g}_{\text{eff}}(\theta)/|g_{\text{eff}}(\theta)|$ and $\hat{\beta} = \cos i\hat{e}_y + \sin i\hat{e}_z$ are the

unit normal vector to the stellar surface and the unit vector directed toward the observer, respectively. The specific intensity $I_c$ is represented as

$$I_c(\theta, \eta, \alpha, \mu) = \frac{F_c(\theta)}{\pi (1 - \epsilon_{\lambda_c} / 3)} (1 - \epsilon_{\lambda_c} + \epsilon_{\lambda_c} \cos \mu), \quad (9)$$

where $F_c(\theta)$ is the emergent continuum radiation flux $\lambda_c$ at $\lambda_c$; $\epsilon_{\lambda_c}$ is the limb darkening coefficient $\epsilon_\lambda$ at $\lambda_c$ of the continuum spectrum that we interpolated in Claret (2012) for the solar chemical abundance and the local $[T_{\text{eff}}(\theta), \log g_{\text{eff}}(\theta)]$. In a more detailed treatment of spectral lines, the limb darkening coefficient should vary with the wavelength inside the spectral line (Collins & Truax 1995; Chauville et al. 2001).

Let $(x, y, z)$ be a Cartesian reference system centered on the apparent stellar disk, where $(x, y)$ are in the plane of the sky, $z$ is directed towards the observer, $y$ is in the plane containing $z$ and
<figure><img src="image_18.png" /><figcaption>Fig. 16. Evolution of the differential phase as a function of the wavelength through the He I line for different values of α. Differential phases were computed for different orientations and a 60-m baseline. The position angle of the star is assumed to be equal to zero, the inclination angle of the rotation axis is i = 45°, and the stellar angular diameter is assumed to be 1.3 mas.</figcaption></figure>

the stellar rotation axis. The radial velocity projected towards the observer of a surface point (r, θ, φ) is then

$$ 
\begin{aligned}
V_{\text{proj}}(\theta, \phi, \eta, \alpha, i) &= \left[ \overline{\Omega(\theta)} \wedge \overline{r(\theta, \eta, \alpha)} \right] \hat{\beta} \\
&= \Omega(\theta, \alpha) R_{\text{eq}}(\eta) x(\eta, \theta, \phi) \sin i \\
x(\eta, \theta, \phi) &= r(\theta, \eta, \alpha) \sin \theta \sin \phi,
\end{aligned}
\quad (10) $$

where Ω(θ, α) is given by Eq. (5); i is the inclination angle of the stellar rotation axis; r(θ, η, α) is the radius vector a point in the stellar surface.

The calculation of the flux line profile F(λ) was made as in all preceding calculations: the star was divided into surface elements to which were associated with the respective T<sub>eff</sub>(θ), log g<sub>eff</sub>(θ) and V<sub>proj</sub>(θ, φ, α, i) that enable one to obtain the local Doppler-shifted specific intensity directed toward the observer. The integration over the visible stellar hemisphere of the projected specific intensities produces F(λ). The central intensity and equivalent widths dependent on [T<sub>eff</sub>, log g<sub>eff</sub>] needed to determine the intrinsic line profile H(λ) lower case In Eq. (8) were taken from Frémat et al. (2005). Equation (9) was used to obtain the local specific radiation intensity, either monochromatic or bolometric, directed toward the observer, which determines the sought distribution of the surface brightness. For the test examples calculated here, the input parameters were the polar effective temperature T<sub>p</sub> = 19 000 K, the polar gravity g<sub>p</sub> = 4.0, the inclination angle i = 45°, gravity darkening exponent of the gravity darkening law β = 0.25, the ratio η = 0.9, and several values of the differential rotation parameter α. The test He II 4471 line-like profiles obtained are shown in Fig. 14, where the strong variation of the line profiles according to the value of α is clearly apparent.

The corresponding maps of brightness distribution containing the iso-radial velocity zone for λ = 4439 Å, calculated for the same fundamental parameters, are shown in Fig. 15. At the top of this figure, from left to right, the images are for α = -1.0 and -0.5, while in the bottom panels we have from left to right, α = 0.0 and 1.0. We can clearly see in this figure that the regions crossed by a given iso-radial velocity curve are quite different. The iso-radial velocity curves thus run across zones with different (T<sub>eff</sub>, log g) contributing to the final line

flux profile with local line profiles having different intensities and Doppler broadenings. In Fig. 16 we note that according to the parameter α, not only the amplitude variations of differential phase are stronger than the expected uncertainties affecting the observed differential phases, but the strongest phase deviations are radial-velocity shifted, which are α-dependent. In this figure we notice that deviations of maxima in the differential phases exceed 10 km s<sup>-1</sup> and the amplitude variation of this maximum widely exceeds 3°. These limits correspond to the VEGA/CHARA sensitivity. A spectral line in α Cep that would behave like the witness spectral line treated here could lead us to reliable inferences of the differential parameter α. Looking for appropriate configurations of the interferometric array might improve the estimates even more.

## 6. Conclusions

Differential phases in the Hα line of α Cep observed during six observations at the VEGA/CHARA array in September 2008 and October 2010 were presented. We successfully applied a method for determining the position angle of the rotationally deformed star that takes advantage of the differential phases obtained at different configurations of the interferometric array, and is independent of the knowledge of the stellar fundamental parameters. We have thus obtained PA = -157<sup>+17°</sup><sub>-10°</sub>.

We discussed the feasibility of determining differential phases using the gravity darkening exponent β in von Zeipel's relation, the inclination angle i of the star, and the differential rotation parameter α of Maunders's expression for the stellar surface rotation law. The discussion showed that the exponent β is a function of the rotation rate η = F<sub>c</sub>/F<sub>g</sub> (F<sub>c</sub> = centrifugal acceleration at the equator; F<sub>g</sub> = gravitational acceleration) and of α. Its determination is probably not straightforward and needs a careful discussion of interferometric imaging of stellar surfaces.

From models of the Hα emitted by a rapidly rotating star with the fundamental parameters α Cep, we showed that for very ideal configurations of interferometric arrays it would be possible to estimate the inclination angle i. Nevertheless, a knowledge
of the stellar fundamental parameters is needed. They can be in-
ferred by iteration together with the inclination angle *i*. We fi-
nally concluded that spectral lines that are not strongly affected
by Stark broadening and have collision-dominated source func-
tions should be observed to estimate *i* and α with differential
phases. They are sensitive to these quantities in three different
ways: height of the deviation wings, the value of the maximum
deviation, and the spectral position of the maximum deviation.

Acknowledgements. VEGA is a collaboration between CHARA and OCA/LAOG/CRAL/LESIA that has been supported by the French programs PNPS and ASHRA, by INSU and by the Région PACA. The project has obviously benefitted from the strong support of the OCA and CHARA technical teams. The CHARA Array is operated with support from the National Science Foundation through grant AST-0908253, the W. M. Keck Foundation, the NASA Exoplanet Science Institute, and from Georgia State University. This work has made use of the BeSS database, operated at GEPI, Observatoire de Meudon, France: http://basebe.obspm.fr, use of the Jean-Marie Mariotti Center SearchCal service<fnref n="1" /> co-developed by FIZEAU and LAOG, and of CDS Astronomical Databases SIMBAD and VIZIER<fnref n="2" />. We are grateful to an anonymous referee for her/his valuable suggestions that helped to improve the presentation of our results.

References

Abt, H. A., & Morrell, N. I. 1995, ApJS, 99, 135

Aufdenberg, J. P., Mérand, A., Coudé du Foresto, V., et al. 2006, ApJ, 645, 664

Berio, P., Mourard, D., Bonneau, D., et al. 1999, J. Opt. Soc. Am. A, 16, 872

Bernacca, P. L., & Perinotto, M. 1973a, in Main sequence single stars, A cata-
logue of stellar rotational velocities (Universita di Padova), 1

Bernacca, P. L., & Perinotto, M. 1973b, in Main sequence spectroscopic binaries, A catalogue of stellar rotational velocities (Universita di Padova), 2

Bernacca, P. L., & Perinotto, M. 1973c, in Evolved stars, A catalogue of stellar rotational velocities (Universita di Padova), 4

Bowman, A. W., & Azzalini, A. 1997, Applied smoothing techniques for data analysis: the kernel approach with S-plus illustrations, Oxford statistical science series No. 18 (Clarendon Press)

Chauville, J., Zorec, J., Ballereau, D., et al. 2001, A&A, 378, 861

Che, X., Monnier, J. D., Zhao, M., et al. 2011, ApJ, 732, 68

Claret, A. 2012, A&A, 538, A3

Collins, II, G. W., & Truax, R. J. 1995, ApJ, 439, 860

Domiciano de Souza, A., Vakili, F., Jankov, S., Janot-Pacheco, E., & Abe, L. 2002, A&A, 393, 345

Domiciano de Souza, A., Kervella, P., Jankov, S., et al. 2003, A&A, 407, L47

Domiciano de Souza, A., Hadjara, M., Vakili, F., et al. 2012a, A&A, 545, A130

Domiciano de Souza, A., Zorec, J., & Vakili, F. 2012b, SF2A Proc., in press

Douglas, A. V. 1926, ApJ, 64, 262

Efron, B., & Tibshirani, R. 1986, Statist. Sci., 1, 54

Endal, A. S., & Sofia, S. 1979, ApJ, 232, 531

Espinosa Lara, F., & Rieutord, M. 2012, A&A, 547, A32

Frémat, Y., Zorec, J., Hubert, A.-M., & Floquet, M. 2005, A&A, 440, 305

Gray, R. O., Corbally, C. J., Garrison, R. F., McFadden, M. T., & Robinson, P. E. 2003, AJ, 126, 2048

Heger, A., & Langer, N. 2000, ApJ, 544, 1016

Hubeny, I., & Lanz, T. 1995, ApJ, 439, 875

Johnson, H. L., & Morgan, W. W. 1953, ApJ, 117, 313

Kervella, P., & Domiciano de Souza, A. 2006, A&A, 453, 1059

Kiziloglu, N., & Civelek, R. 1996, A&A, 307, 88

Kurucz, R. L. 1979, ApJS, 40, 1

Lamers, H. J. G. L. M., & Cassinelli, J. P. 1999, Introduction to Stellar Winds (Cambridge University Press)

Le Bouquin, J.-B., Absil, O., Benisty, M., et al. 2009, A&A, 498, L41

Lucy, L. B. 1967, ZAp, 65, 89

Maeder, A., Meynet, G., & Ekström, S. 2007, in From Stars to Galaxies: Building the Pieces to Build Up the Universe, eds. A. Vallenari, R. Tantalo, L. Portinari, & A. Moretti, ASP Conf. Ser., 374, 13

Maeder, A., Georgy, C., & Meynet, G. 2008, A&A, 479, L37

McAlister, H. A., ten Brummelaar, T. A., Gies, D. R., et al. 2005, ApJ, 628, 439

Meynet, G., & Maeder, A. 2000, A&A, 361, 101

Mihalas, D. 1978, Stellar atmospheres, 2nd Ed. (San Francisco: W. H. Freeman and Co.)

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

Mourard, D. 1990, Ph.D. Thesis, Université de Nice

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

Mourard, D., Tallon, M., Bério, P., et al. 2010, in SPIE Conf. Ser., 7734

Pinsonneault, M. H., Deliyannis, C. P., & Demarque, P. 1991, ApJ, 367, 239

Reiners, A. 2003, A&A, 408, 707

Royer, F., Zorec, J., & Gómez, A. E. 2007, A&A, 463, 671

Sackmann, I. J. 1970, A&A, 8, 76

Shimazaki, H., & Shinomoto, S. 2007, Neural Computation, 19, 1503

Simon, T., & Landsman, W. B. 1997, ApJ, 483, 435

Simon, T., Ayres, T. R., Redfield, S., & Linsky, J. L. 2002, ApJ, 579, 800

Snodgrass, H. B. 1984, Sol. Phys., 94, 13

Talon, S., Zahn, J.-P., Maeder, A., & Meynet, G. 1997, A&A, 322, 209

Tatulli, E., Millour, F., Chelli, A., et al. 2007, A&A, 464, 29

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

Thomas, R. N. 1965, Some aspects of non-equilibrium thermodynamics in the presence of a radiation field (University of Colorado Press)

Thomas, R. N. 1983, Stellar atmospheric structural patterns, 471 (NASA-CNRS)

Townsend, R. H. D. 1997, MNRAS, 284, 839

Tuominen, I. V. 1972, Ann. Acad. Sci. Fennicae, 392

Uesugi, A., & Fukuda, I. 1970, Catalogue of rotational velocities of the stars (University of Kyoto)

van Belle, G. T., Ciardi, D. R., Thompson, R. R., Akeson, R. L., & Lada, E. A. 2001, ApJ, 559, 1155

van Belle, G. T., Ciardi, D. R., ten Brummelaar, T., et al. 2006, ApJ, 637, 494

von Zeipel, H. 1924, MNRAS, 84, 665

Walter, F. M., Matthews, L. D., & Linsky, J. L. 1995, ApJ, 447, 353

Zahn, J.-P. 1992, A&A, 265, 115

Zhao, M., Monnier, J. D., Pedretti, E., et al. 2009, ApJ, 701, 209

Zorec, J., Frémat, Y., Domiciano de Souza, A., et al. 2011, A&A, 526, A87

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

<footnote n="2">² Available at http://cdsweb.u-strasbg.fr/</footnote>