# Accurate p-mode measurements of the G0V metal-rich CoRoT target HD 52265

J. Ballot, Laurent Gizon, R. Samadi, G. Vauclair, O. Benomar, H. Bruntt, B.

Mosser, T. Stahn, G.A. Verner, T.L. Campante, et al.

► To cite this version:

J. Ballot, Laurent Gizon, R. Samadi, G. Vauclair, O. Benomar, et al. Accurate p-mode measurements of the G0V metal-rich CoRoT target HD 52265. Astronomy & Astrophysics - A&A, 2011, 530, pp.97. 10.1051/0004-6361/201116547. [hal-00632626](https://hal.science/hal-00632626v1)

HAL Id: hal-00632626

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

Submitted on 24 Sep 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 CC BY 4.0 - Attribution - International License
Accurate p-mode measurements of the G0V metal-rich CoRoT* target HD 52265

J. Ballot<sup>1,2</sup>, L. Gizon<sup>3</sup>, R. Samadi<sup>4</sup>, G. Vauclair<sup>1,2</sup>, O. Benomar<sup>5,6</sup>, H. Bruntt<sup>4</sup>, B. Mosser<sup>4</sup>, T. Stahn<sup>3</sup>, G. A. Verner<sup>7,8</sup>, T. L. Campante<sup>9,10</sup>, R. A. García<sup>11</sup>, S. Mathur<sup>12</sup>, D. Salabert<sup>13,14,15</sup>, P. Gaulme<sup>5</sup>, C. Régulo<sup>13,14</sup>, I. W. Roxburgh<sup>7</sup>, T. Appourchaux<sup>5</sup>, F. Baudin<sup>5</sup>, C. Catala<sup>4</sup>, W. J. Chaplin<sup>8</sup>, S. Deheuvels<sup>4</sup>, E. Michel<sup>4</sup>, M. Bazot<sup>10</sup>, O. Creevey<sup>13,14</sup>, N. Dolez<sup>1,2</sup>, Y. Elsworth<sup>8</sup>, K. H. Sato<sup>11</sup>, S. Vauclair<sup>1,2</sup>, M. Auvergne<sup>4</sup>, and A. Baglin<sup>4</sup>

¹ Institut de Recherche en Astrophysique et Planétologie, CNRS, 14 avenue Edouard Belin, 31400 Toulouse, France
e-mail: jballot@astobs-mip.fr

² Université de Toulouse, UPS-OMP, IRAP, Toulouse, France

³ Max-Planck-Institut für Sonnensystemforschung, 37191 Katlenburg-Lindau, Germany

⁴ LESIA, UMR 8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon, France

⁵ Institut d'Astrophysique Spatiale, CNRS, Université Paris XI, 91405 Orsay, France

⁶ Sydney Institute for Astronomy, School of Physics, University of Sydney, NSW 2006, Australia

⁷ Astronomy Unit, Queen Mary, University of London Mile End Road, London E1 4NS, UK

⁸ School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, UK

⁹ Danish AstroSeismology Centre, Department of Physics and Astronomy, University of Aarhus, 8000 Aarhus C, Denmark

¹⁰ Centro de Astrofísica, DFA-Faculdade de Ciências, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal

¹¹ Laboratoire AIM, CEA/DSM, CNRS, Université Paris Diderot, IRFU/SAp, Centre de Saclay, 91191 Gif-sur-Yvette Cedex, France

¹² High Altitude Observatory, NCAR, PO Box 3000, Boulder, CO 80307, USA

¹³ Universidad de La Laguna, Dpto de Astrofísica, 38206 La Laguna, Tenerife, Spain

¹⁴ Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife, Spain

¹⁵ Université de Nice Sophia-Antipolis, CNRS, Observatoire de la Côte d'Azur, BP 4229, 06304 Nice Cedex 4, France

Received 20 January 2011 / Accepted 22 March 2011

ABSTRACT

**Context.** The star HD 52265 is a G0V metal-rich exoplanet-host star observed in the seismology field of the CoRoT space telescope from November 2008 to March 2009. The satellite collected 117 days of high-precision photometric data on this star, showing that it presents solar-like oscillations. HD 52265 was also observed in spectroscopy with the Narval spectrograph at the same epoch.

**Aims.** We characterise HD 52265 using both spectroscopic and seismic data.

**Methods.** The fundamental stellar parameters of HD 52265 were derived with the semi-automatic software VWA, and the projected rotational velocity was estimated by fitting synthetic profiles to isolated lines in the observed spectrum. The parameters of the observed p modes were determined with a maximum-likelihood estimation. We performed a global fit of the oscillation spectrum, over about ten radial orders, for degrees $l = 0$ to 2. We also derived the properties of the granulation, and analysed a signature of the rotation induced by the photospheric magnetic activity.

**Results.** Precise determinations of fundamental parameters have been obtained: $T_{\text{eff}} = 6100 \pm 60$ K, $\log g = 4.35 \pm 0.09$, [M/H] = 0.19 ± 0.05, as well as $v \sin i = 3.6_{-1.0}^{+0.3}$ km s⁻¹. We have measured a mean rotation period $P_{\text{rot}} = 12.3 \pm 0.15$ days, and find a signature of differential rotation. The frequencies of 31 modes are reported in the range 1500–2550 µHz. The large separation exhibits a clear modulation around the mean value $\Delta v = 98.3 \pm 0.1$ µHz. Mode widths vary with frequency along an S-shape with a clear local maximum around 1800 µHz. We deduce lifetimes ranging between 0.5 and 3 days for these modes. Finally, we find a maximal bolometric amplitude of about 3.96 ± 0.24 ppm for radial modes.

**Key words.** stars: individual: HD 52265 – asteroseismology – stars: fundamental parameters – stars: rotation – stars: solar-type

1. Introduction

Late-type stars oscillate when they have sufficiently deep convective envelopes. These oscillations are p modes, i.e., acoustic waves trapped in the stellar interior. Despite being normally stable, these modes are stochastically excited by near-surface turbulent motions (e.g. Goldreich & Keeley 1977). Their amplitudes remain small and produce tiny luminosity fluctuations. Observing them is now possible thanks to high-performance

dedicated instrumentation such as space-borne photometry. The current space-based missions CoRoT (COnvection, ROtation, and planetary Transits, Baglin et al. 2006) and Kepler (Koch et al. 2010) provide long uninterrupted time series of high-precision photometric data over months or years, allowing for the detection and the characterisation of p modes in late-type stars (e.g. Michel et al. 2008; Gilliland et al. 2010). These high-quality data give an exciting opportunity to conduct seismological analyses of solar-like stars. CoRoT has been successful in providing observations that allow detailed seismic analysis of late-type main-sequence (e.g. Appourchaux et al. 2008; Barban et al. 2009; García et al. 2009; Mosser et al. 2009b; Benomar et al. 2009b; Mathur et al. 2010a; Deheuvels et al. 2010) and

<footnote>* The CoRoT space mission, launched on December 27th 2006, has been developed and is operated by CNES, with the contribution of Austria, Belgium, Brazil, ESA (RSSD and Science Programme), Germany and Spain.</footnote>
giant (e.g. De Ridder et al. 2009; Hekker et al. 2009; Mosser et al. 2010) stars showing solar-like oscillations. Such accurate asteroseismic observations provide unprecedented constraints on the stellar structure of these classes of stars (e.g. Deheuvels & Michel 2010; Miglio et al. 2010). The Kepler mission is also currently providing interesting results on these two categories of stars (e.g. Chaplin et al. 2010; Huber et al. 2010; Kallinger et al. 2010; Mathur et al. 2011).

In this paper, we study HD 52265, a G0, rather metal-rich, main-sequence star hosting a planet that was independently discovered in 2000 by Butler et al. (2000) and Naef et al. (2001). HD 52265 was observed by CoRoT during almost four months from 13 November 2008 to 3 March 2009. CoRoT observations were complemented by ground-based observations achieved with the Narval spectrograph at the Pic du Midi observatory in December 2008 and January 2009, i.e. during the CoRoT observations.

Up to now, only a very few exoplanet-host stars (EHS) have been studied seismically. Ground-based observations of μ Ara, hosting two known planets, allowed Bouchy et al. (2005) to extract about 30 modes based on an eight-night observation and constrain its fundamental parameters (Bazot et al. 2005; Soriano & Vauclair 2010). Recently, similar constraints have been placed on another EHS, i Hor (Vauclair et al. 2008). CoRoT has also observed HD 46375, an unevolved K0 star hosting a Saturn-like planet, for which a mean large separation has been measured (Gaulme et al. 2010). Kepler also observed the EHS HAT-P-7, for which frequencies of 33 p modes were measured within 1.4 μHz accuracy, as well as the EHS HAT-P-11 and TrES-2, for which large separations have been estimated (Christensen-Dalsgaard et al. 2010). The opportunities to get seismic constraints on EHS are then rare. The high-quality observations of HD 52265 allow us to determine its seismic properties with an accuracy never obtained on EHS.

In Sect. 2, we first report the fundamental parameters of the star obtained with our spectroscopic study, including a detailed composition analysis. After describing the CoRoT photometric data in Sect. 3, we measure the rotation period and discuss the activity of the star (Sect. 4). We also measure granulation properties (Sect. 5) and extract the main p-mode characteristics – frequencies, amplitudes, and lifetimes – in Sect. 6, before concluding (Sect. 7).

## 2. Fundamental stellar parameters

### 2.1. Distance, luminosity, and chromospheric activity

The star HD 52265 (or HIP 33719, or HR 2622) was initially classified as G0III-IV in the Bright Star Catalogue (Hoffleit & Jaschek 1982) before finally being identified as a G0V dwarf thanks to Hipparcos parallax measurements (see Butler et al. 2000). It has a magnitude $V = 6.301$ and a parallax of $\pi = 34.54 \pm 0.40$ mas in the revised Hipparcos catalogue by van Leeuwen (2007) and this corresponds to a distance $d = 28.95 \pm 0.34$ pc.

The catalogue of van Leeuwen (2007) also gives the magnitude $HP = 6.4132 \pm 0.0006$. Using the bolometric correction $BC_{HP} = -0.17 \pm 0.01$ after Cayrel et al. (1997), one finds an absolute bolometric magnitude $M_{bol} = 3.93 \pm 0.06$; taking $M_{bol,\odot} = 4.75$, this gives luminosity of $L/L_{\odot} = 2.09 \pm 0.24$. For comparison, we also derived the luminosity by using the V magnitude and the bolometric correction $BC_{V} = -0.03 \pm 0.01$ (Flower 1996) and obtained a fully consistent value for L.

<table><caption>**Table 1.** Summary of general parameters of HD 52265.</caption><thead><tr><th></th><th>Ref.</th><th></th></tr></thead><tbody><tr><td>π</td><td>34.54 ± 0.40 mas</td><td>1</td></tr><tr><td>L/L<sub>☉</sub></td><td>2.09 ± 0.24</td><td>1, 2</td></tr><tr><td>log R'<sub>HK</sub></td><td>-5.02</td><td>3</td></tr><tr><td>T<sub>eff</sub></td><td>6100 ± 60 K</td><td>4</td></tr><tr><td>log g</td><td>4.35 ± 0.09</td><td>4</td></tr><tr><td>[M/H]</td><td>0.19 ± 0.05</td><td>4</td></tr><tr><td>ν sin i</td><td>3.6<sup>+0.3</sup><sub>-1.0</sub> km s<sup>-1</sup></td><td>4</td></tr></tbody></table>

**Notes.** See text for details.

**References.** (1) van Leeuwen (2007); (2) Cayrel et al. (1997); (3) Wright et al. (2004); (4) this work.

The magnetic activity of a star can be revealed by chromospheric emission lines, especially the very commonly used Ca II H and K lines, from which is derived the standard activity index $R'_{HK}$. For HD 52265, Wright et al. (2004) have measured log $R'_{HK} = -5.02$, which is consistent with the value of log $R'_{HK} = -4.99$ obtained by Butler et al. (2000). These values, comparable to the solar one, indicate HD 52265 can be classified as a magnetically quiet star. These parameters are summarised in Table 1.

### 2.2. Narval spectroscopic observations

Complementary ground-based observations were obtained with the Narval<fnref n="1" /> spectropolarimeter installed on the Bernard Lyot Telescope at the Pic du Midi Observatory (France). Spectra with a resolution of 65 000 were registered simultaneously in classical spectroscopy (Stokes I) and circularly polarised light (Stokes V) over the spectral domain 370–1000 nm in a single exposure. Nine high signal-to-noise ratio spectra were obtained in December 2008 and January 2009, i.e., during the CoRoT observation for this star: 1 spectrum recorded on 2008/12/20, 1 on 2008/12/21, 2 on 2009/1/10, and 5 on 2009/1/11. No Zeeman signature has been recovered in the polarised spectra. This gives us an upper limit of ~1 G for the average magnetic field by combining the five spectra registered on 2009/1/11 and ~2 G for the other days.

We analysed the observed spectrum of HD 52265 using the semi-automatic software package VWA (Bruntt et al. 2004; Bruntt 2009). We adopted atmospheric models interpolated in the MARCS grid (Gustafsson et al. 2008) and atomic parameters from VALD (Kupka et al. 1999). The abundances for more than 500 lines were calculated iteratively by fitting synthetic profiles to the observed spectrum using SYNTH (Valenti & Piskunov 1996). For each line, the abundances were calculated differentially with respect to the same line in a solar spectrum from Hinkle et al. (2000). The atmospheric parameters ($T_{eff}$, log g, [Fe/H], and $v_{micro}$) were determined by requiring that Fe lines give the same abundance independent of equivalent width (EW), excitation potential or ionisation stage. Only weak lines were used ($EW < 90$ mÅ) for this part of the analysis, but stronger lines ($EW < 140$ mÅ) were used for the calculation of the final mean abundances. The uncertainties on the atmospheric parameters and the abundances were determined by perturbing the best-fitting atmospheric parameters as described by Bruntt et al. (2008).

The results are $T_{eff} = 6100 \pm 31 \pm 50$ K, log g = 4.35 ± 0.04 ± 0.08, $v_{micro} = 1.18 \pm 0.04 \pm 0.05$ km s⁻<fnref n="1" />. We give two

<footnote n="1">¹ http://www.ast obs-mip.fr/projets/narval/v1/</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Abundance pattern of HD 52265 for 23 elements. Circles and box symbols are used for the mean abundance from neutral and singly ionised lines, respectively. The yellow horizontal bar marks the mean metallicity with $1-\sigma$ uncertainty range. The chemical symbols in blue correspond to α elements.</figcaption></figure>

<table><caption>**Table 2.** Abundances in HD 52265 relative to the Sun (log $\epsilon/\epsilon_\odot$) and number of lines (N) used in spectral analysis.</caption><thead><tr><th>Elem.</th><th>log ε/ε<sub>☉</sub><br>[dex]</th><th>N</th><th>Elem.</th><th>log ε/ε<sub>☉</sub><br>[dex]</th><th>N</th></tr></thead><tbody><tr><td>Li I</td><td>1.91</td><td>1</td><td>V I</td><td>0.18</td><td>12</td></tr><tr><td>C I</td><td>0.14</td><td>6</td><td>Cr I</td><td>0.20</td><td>21</td></tr><tr><td>O I</td><td>0.11</td><td>3</td><td>Cr II</td><td>0.22</td><td>4</td></tr><tr><td>Na I</td><td>0.22</td><td>4</td><td>Mn I</td><td>0.15</td><td>7</td></tr><tr><td>Mg I</td><td>0.14</td><td>2</td><td>Fe I</td><td>0.19</td><td>268</td></tr><tr><td>Mg II</td><td>0.17</td><td>2</td><td>Fe II</td><td>0.21</td><td>23</td></tr><tr><td>Al I</td><td>0.17</td><td>3</td><td>Co I</td><td>0.17</td><td>9</td></tr><tr><td>Si I</td><td>0.19</td><td>34</td><td>Ni I</td><td>0.19</td><td>72</td></tr><tr><td>Si II</td><td>0.18</td><td>2</td><td>Cu I</td><td>0.16</td><td>2</td></tr><tr><td>S I</td><td>0.20</td><td>2</td><td>Zn I</td><td>0.17</td><td>2</td></tr><tr><td>Ca I</td><td>0.19</td><td>12</td><td>Y II</td><td>0.20</td><td>4</td></tr><tr><td>Ca II</td><td>0.32</td><td>1</td><td>Zr I</td><td>0.18</td><td>1</td></tr><tr><td>Sc II</td><td>0.26</td><td>6</td><td>Ce II</td><td>0.13</td><td>1</td></tr><tr><td>Ti I</td><td>0.18</td><td>27</td><td>Nd II</td><td>0.01</td><td>1</td></tr><tr><td>Ti II</td><td>0.27</td><td>10</td><td></td><td></td><td></td></tr></tbody></table>

uncertainties: the first one is the intrinsic error and the second is our estimated systematic error. The systematic errors were estimated from a large sample of F5-K1 type stars (Bruntt et al. 2010). These two errors are independent and should be added quadratically.

We have also fitted profiles to 68 isolated lines to find the best combination of $v \sin i$ and macroturbulence. Lines with high $\chi^2$ of the best fit have been discarded, and 49 of the best lines remain. We found $v \sin i = 3.6^{+0.3}_{-1.0}$ km s$^{-1}$.

We calculated the mean metallicity from the six metals with at least 10 lines (Si, Ti, V, Cr, Fe, and Ni): [M/H] = +0.19 ± 0.05, where the uncertainty includes the errors on the atmospheric parameters. In Table 2 we list the abundances of the 23 elements that are also shown in Fig. 1. The estimated internal error is 0.04 dex for all elements.

HD 52265 is relatively metal rich (≈50% more heavy elements than the Sun), but there is no evidence of any strong enhancement of α elements. As seen in Fig. 1, almost all elements are enhanced by a similar amount. We can then use a global scaling of the metallicity of +0.19. Nevertheless, it is seen that the star has a high lithium abundance (1.91 dex). For the region around the Li I line at 6707.8 Å we used the line list from Ghezzi et al. (2009) but did not include the CN molecular lines. Considering the abundance log $\epsilon(\text{Li})_\odot = 1.05 \pm 0.1$ for the Sun (Asplund et al. 2009), we got log $\epsilon(\text{Li}) = 2.96 \pm 0.15$ for HD 52265. This measure is consistent with results of Israeli et al. (2004). According to their sample, such a high abundance of lithium is not unusual for a star with $T_{\text{eff}} = 6100$ K.

The fundamental parameters are summarised in Table 1. These were compared to previous work, especially to Valenti & Fischer (2005). We recover consistent results for $T_{\text{eff}}$, [M/H],

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** (Top) Light curve of HD 52265 observed by CoRoT, starting on 2008/11/13. Points taken across the SAA have been removed. The green curve is resampled within 8-h bins. (Bottom) The fluctuations around the average obtained after detrending with a 3rd-order polynomial.</figcaption></figure>

and log g, within the error bars. The major change concerns $v \sin i$, which is significantly reduced from 4.7 ± 0.5 km s$^{-1}$ to $3.6^{+0.3}_{-1.0}$ km s$^{-1}$, probably owing to the inclusion of macroturbulence in the present work.

## 3. CoRoT photometric observations

### 3.1. CoRoT data

HD 52265 was observed by CoRoT in the seismology field (see Auvergne et al. 2009) during 117 days starting on 2008 November 13 to 2009 March 3, during the second long run in the galactic anti-centre direction (LRa02). We used time series that are regularly spaced in the heliocentric frame with a 32 s sampling rate, i.e. the so-called Helreg level 2 (N2) data (Samadi et al. 2007a). The overall duty cycle of the time series is $r_{\text{dc}} = 94.6\%$. However, the extra noise on data taken across the south-Atlantic anomaly (SAA) creates strong harmonics of the satellite orbital period and of the day (for more details, see Auvergne et al. 2009). When all of the data registered across the SAA are removed, the duty cycle becomes $r_{\text{dc}} = 90.2\%$.

The light curve of HD 52265 from N2 data is plotted in Fig. 2-top. The light curve shows a slowly decreasing trend due to instrumental effects, and especially due to the ageing of CCDs. This trend is removed with a third-order polynomial fit. Residual fluctuations are plotted in the bottom of Fig. 2.

The power spectral density (PSD) is computed with a standard fast Fourier transform (FFT) algorithm, and we normalised
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Power spectrum of the observation window with a zoom around the zero frequency in the insert.</figcaption></figure>

it as the so-called one-sided power spectral density (Press et al. 1992). We interpolate the gaps produced by SAA with parabola fitting the points around the gap, in a similar way to the usual interpolation performed in N2 data for missing points.

To verify that the interpolation method has no influence on the results presented in this paper, we performed similar analyses on spectra where:

1. the PSD is computed as the FFT of the N2 data, and main harmonics and aliases induced by the SAA noise are removed in the spectrum;

2. all of the gaps and missing points are interpolated with the in-painting method described in Mathur et al. (2010a) and Sato et al. (2010).

We have also verified that the order of the polynomial used to detrend the data has no impact on the results, by using first to ninth order polynomials.

## 3.2. Influence of the observation window

At the precision reached by the current analysis, windowing effects must be taken into account when normalising the PSD. The first consequence of windowing is a lack of power in the PSD by a factor of $r_{dc}$. This factor should be taken into account a posteriori to avoid a general underestimation of the power by around 10%. Of course, the interpolation we performed injected power into the spectrum, but only at very low frequency. Thus, the PSD is properly calibrated by dividing it by $r_{dc}$.

Moreover, owing to the repetitive pattern of the gaps, the power is spread in side lobes. To estimate the resulting leakage, we computed the power spectrum of the observation window (Fig. 3). On a logarithmic scale we clearly see a forest of peaks corresponding to the orbital frequency (161.7 µHz), its harmonics, as well as daily aliases. These peaks have lower amplitudes than the central peak (smaller than 1%), and there is no peak close to the central one. As a result, all of the p modes, which are narrow structures, should be modelled well by Lorentzian profiles (see Sect. 6). The mode aliases then have very small amplitudes and do not perturb the other modes significantly, but a non-negligible part of the mode power is spread in the background, and the fitted mode heights will be underestimated by a factor of $r_{dc}$.

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** PSD of HD 52265 at low frequency. Dashed line shows an oversampling of the PSD by a factor of 10.</figcaption></figure>

However, we neglect the leakage effects on very broad structures, such as the background ones (see Sect. 5), because their widths are at least as broad as the range where the significant peaks of the window spectrum are found.

# 4. Activity, spots, and rotation period

The light curve (Fig. 2) exhibits clear variations on characteristic time scales of a few days and with amplitudes around 100 ppm, which we attribute to stellar activity. As shown in Sect. 2.1, this star is magnetically quiet, much like the Sun. The amplitude of the modulation is then compatible with those produced by spots in the photosphere. The modulation appears quasi-periodic with clear repetitive patterns.

This modulation produces a significant peak in the power spectrum at 1.05 µHz and a second very close one, at 0.91 µHz (Fig. 4). The signature is then not a single peak, but a slightly broadened structure. We interpret the broadening as the signature of differential rotation, since spots may appear at different latitudes that rotate with different periods. Peaks around 2 and 3 µHz are simply the harmonics of the rotation frequency.

A wavelet analysis (e.g. Torrence & Compo 1998) of the curve was performed to verify the persistence of this characteristic frequency during the complete length of observation. Interestingly, this method allowed us to monitor the evolution with time of the power spectral distribution, as well as to resolve the uncertainty between the rotation period and the first harmonic that could be observed in the power spectrum (see Mathur et al. 2010b). We used the Morlet wavelet (a moving Gaussian envelope convolved with a sinusoid function with a varying frequency) to produce the wavelet power spectrum shown Fig. 5. It represents the correlation between the wavelet with a given frequency along time. The frequency around 1 µHz previously identified as the rotation rate is very persistent during the whole time series, while the peaks at 2 and 3 µHz are not. Thus, the 1 µHz signature does not come from strong localised artefacts, but is the consequence of a continuous modulation. This confirms that the rotation period is around 12 days.

Finally, we performed spot-modelling of the light curve following the method described by Mosser et al. (2009a). The best-fit model is presented in Fig. 6. The fit is obtained with an inclination of $i = 30 \pm 10^{\circ}$. When we assumed a solid rotation,
<figure><img src="image_6.png" /><figcaption>Fig. 5. Wavelet power spectrum for HD 52265 as a function of the frequency of the Morlet wavelet (y-axis) and time (x-axis). The colour scale is logarithmic, white corresponds to the highest values, and red hatching indicates the cone of influence delimiting the reliable periodicity.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 6. Spot modelling of HD 52265. Relative luminosity fluctuations (d$\phi$/$\phi$) of HD 52265 are plotted as a function of time. The dots in the upper panel represent the data binned every CoRoT orbit (6184 s) and the solid grey curve is the best-fit model. Dots in the bottom panel represent residuals. Vertical grey lines indicate the equatorial rotation period.</figcaption></figure>

we found a rotation period $P_{rot} = 12.3 \pm 0.15$ days. However, considering differential rotation, with a law

$$P_{rot}(\theta) = P_{eq}/(1 - K \sin^2 \theta) \quad (1)$$

for the rotation as a function of the latitude $\theta$, provides a better fit. We derive an equatorial rotation period $P_{eq} = 11.7_{-0.2}^{+0.6}$ days and a differential rotation rate $K = 0.25_{-0.2}^{+0.05}$. The spot life time is short, around $\tau_{spot} = 8.0 \pm 1.5$ days, i.e. about 2/3 of the rotation period. As for similar stars, the modelling is not able to reproduce the sharpest features of the light curve. Fitting them would require introducing too many spots, which lowers the precision of the best-fit parameters. Moreover, some of these features, such as the jumps around $t = 32, 54, \text{ or } 64$ days, certainly do not have a stellar origin, but are probably introduced by instrumental glitches or cosmic events.

<figure><img src="image_8.png" /><figcaption>Fig. 7. Power spectral density of HD 52265 at full resolution (grey curve) and rebinned with a factor of 100 (black curve). The green line shows the spectrum smoothed by a box car with a width equal to the mean large separation $\Delta\nu$. The fitted background is the solid red line, and its three components ($W, B_g, B_a$, see text) are plotted as red dotted lines. P-mode power excess, fitted as a Gaussian profile, is also represented with a red dashed line.</figcaption></figure>

## 5. Fitting the stellar background

The PSD of HD 52265 is shown in Fig. 7. The p modes are clearly visible around 2000 µHz. The dominant peak at 323 µHz is the persisting first harmonic of the satellite orbital period. Before analysing the p-mode pattern, we determine the properties of the background, which gives useful information about the stellar granulation that constrains, for example, 3D models of superficial convection (e.g. Ludwig et al. 2009).

The background $B(\nu)$ of the spectrum is modelled as the sum of three components:

1. a white noise W modelling the photon shot noise, dominating the spectrum at high frequency;

2. a profile of Harvey (1985) to model the granulation spectrum:

$$B_g(\nu) = \frac{4\tau_g\sigma_g^2}{1 + (2\pi\tau_g\nu)^{\alpha_g}} \quad (2)$$

where $\tau_g$ and $\sigma_g$ are characteristic time scale and amplitude of the granulation and $\alpha_g$ an exponent characterising the temporal coherence of the phenomenon;

3. a third component taking the slow drifts or modulations due to the stellar activity, the instrument, etc., into account which we modelled either by an extra Harvey profile $B_a(\nu) = 4\tau_a\sigma_a^2/[1 + (2\pi\tau_a\nu)^{\alpha_a}]$ or a power law $B_a(\nu) = P_a\nu^{-\alpha_a}$.

We fit the PSD between 1 µHz and the Nyquist frequency, excluding the region where p modes are visible ([1200–3200] µHz). Fittings were performed with a maximum-likelihood estimation (MLE), or alternatively with Markov chains Monte Carlo (MCMC), by assuming that the noise follows a $\chi_2$ distribution. Results for $B_g$ and W do not depend significantly on the choice made for $B_a$. Table 3 lists the fitted parameters with statistical formal errors, obtained by inverting the Hessian matrix. These errors do not include any systematics and assume the model is correct.

Complementary fits including the p-mode region were performed by modelling their contribution as a Gaussian function.
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Power spectral density of HD 52265 in the p-mode frequency range at full resolution (grey curve) and smoothed by a 11-bin wide boxcar (black curve). The red line corresponds to the fitted spectrum.</figcaption></figure>

<table><caption>**Table 3.** White noise (W) and parameters of the granulation background (Eq. (2)) fitted on HD 52265 power spectrum.</caption><thead><tr><th><i>W</i><br>[ppm<sup>2</sup> µHz<sup>-1</sup>]</th><th>σ<sub>g</sub><br>[ppm]</th><th>τ<sub>g</sub><br>[s]</th><th>α<sub>g</sub></th></tr></thead><tbody><tr><td>0.542 ± 0.002</td><td>49.8 ± 0.4</td><td>242 ± 2</td><td>2.65 ± 0.15</td></tr></tbody></table>

The results for $B_g$ and W have not been modified. The p-mode Gaussian profile $B_p$ is characterised by a height $H_p^{(max)} = 0.426 \pm 0.016$ ppm$^2$ µHz$^{-1}$, a width $\Delta_p = 390 \pm 20$ µHz, and a central frequency $\nu_p^{(max)} = 2090 \pm 20$ µHz. We then redid fits over a frequency range starting at 5, 10, or 100 µHz, instead of 1 µHz, without any impact on the results. We also verified that removing the two bins still contaminated by the strong orbit harmonic does not influence the results.

## 6. P-mode analysis

Around 2000 µHz, the PSD presents a conspicuous comb structure that is typical of solar-like oscillations. A zoom on the PSD in the region is plotted in Fig. 8. The small separation between modes of degrees $l = 0$ and $2$ is also easily seen, making the mode identification simple. The autocorrelation of the spectrum between 1700 and 2400 µHz provides first estimates of the large and small separations: $\bar{\Delta v} \approx 98.5$ µHz and $\bar{\delta}_{02} \approx 8$ µHz. These values are recomputed from fitted frequencies in Sect. 6.2. We build an échelle diagram using $\bar{\Delta v}$ as the folding frequency. This is plotted in Fig. 9. This échelle diagram makes the mode identification even more obvious. On the right-hand side, the clear single ridge corresponds to $l = 1$ modes, whereas, on the left-hand side, the two ridges are identified as the $l = 2$ and $l = 0$ modes. We completely exclude the possibility that the two ridges on the left are $l = 1$ modes split by rotation for the following reasons. First, the two ridges show clear asymmetry in their power; second, the rotation needed to generate such high splittings is totally incompatible both with the low-frequency signature (Sect. 4) and the spectroscopic observations (Sect. 2.2). We are able to identify modes presenting a significant signal-to-noise ratio for about ten consecutive orders. Their characteristics are extracted with classical methods, described in the following section.

### 6.1. Fitting the spectrum

To estimate the mode parameters, we applied usual spectrum fitting techniques, as in previous similar works (e.g. Appourchaux et al. 2008; Barban et al. 2009; Deheuvels et al. 2010). To summarise the procedure, we assume that the observed PSD is

<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Échelle diagram of HD 52265, plotted with a folding frequency $\Delta v = 98.5$ µHz. The colour map shows the observed power spectrum normalised by the background (S/B), smoothed over 11 bins. Green symbols with error bars indicate the fitted modes listed in Table 4: plus signs, diamonds, and squares indicate $l = 0, 1$, and $2$ modes, respectively. White crosses indicate orbital harmonics.</figcaption></figure>

distributed around a mean profile spectrum $S(v)$ and follows a 2-degree-of-freedom $\chi^2$ statistics (e.g. Duvall & Harvey 1986). $S(v)$ is modelled as the sum of the background $B(v)$, described in Sect. 5, and p-mode profiles $P_{l,n}(v)$ for each considered degree $l$ and radial order $n$. $P_{l,n}(v)$ are multiplets of Lorentzian profiles and read as

$$P_{l,n}(v) = \sum_{m=-l}^{l} \frac{a_{l,m}(i)H_{l,n}}{1 + [2(v - v_{l,n} + mv_s)/\Gamma_{l,n}]^2}, \quad (3)$$

where $m$ is the azimuthal order, $H_{l,n}$, $v_{l,n}$, and $\Gamma_{l,n}$ are the height, the frequency, and the width of the mode, $v_s$ is the rotational splitting, common to all modes, and $a_{l,m}(i)$, the height ratios of multiplet components, are geometrical terms depending only on the inclination angle $i$ (e.g., see Gizon & Solanki 2003). Owing to cancellation effects of averaging the stellar flux over the whole disc, only modes with $l \le 2$ are considered in the present analysis. It is worth noticing that we analyse the full light curve disregarding any possible variations in the p-mode parameters due to any magnetic activity effects as the ones recently uncovered in HD 49933 (García et al. 2010).

The different parameters have been independently fitted by ten different groups, using either MLE (e.g. Appourchaux et al. 1998) or Bayesian priors and performing maximum a posteriori estimation (MAP, e.g. Gaulme et al. 2009) or MCMC (e.g.
Benomar et al. 2009a; Handberg & Campante 2011). Generally, global fits of all of the p modes were performed. Nevertheless, one group has used a local approach by fitting sequences of successive *l* = 2, 0, and 1 modes, according to the original CoRoT recipes (Appourchaux et al. 2006). Between 8 and 15 orders have been included in the fit, depending on the group.

To reduce the dimension of the parameter space, other constraints were used by taking advantage of the smooth variation of heights and widths. As a consequence, only one height and one width, $H_{0,n}$ and $\Gamma_{0,n}$, are fitted for each order. We linked the height (resp., width) of *l* = 2 mode to the height (width) of the nearby *l* = 0 mode through a relation

$$H_{2,n} = r_2^2 H_{0,n+1} \quad (\Gamma_{2,n} = \Gamma_{0,n+1}). \qquad (4)$$

We denote $r_2$ the visibility of the *l* = 2 mode relative to *l* = 0.

Concerning *l* = 1 modes, we can link them either to the previous or to the following *l* = 0 modes, leading to two fitting situations:

$$\text{fit A: } H_{1,n} = r_1^2 H_{0,n+1} \quad (\Gamma_{1,n} = \Gamma_{0,n+1}); \qquad (5)$$

$$\text{fit B: } H_{1,n} = r_1^2 H_{0,n} \quad (\Gamma_{1,n} = \Gamma_{0,n}), \qquad (6)$$

where $r_1$ the visibility of the *l* = 1 mode relative to *l* = 0. Even if the variations in heights and widths are limited, they can still be significant over half the large separation. Thus, results can be different for both cases.

The visibility factors $r_l$ depend mainly on the stellar limb-darkening profile (e.g. Gizon & Solanki 2003). It is then possible to fix them to theoretical values, deduced from stellar atmosphere models, or leave them as a free parameters. Six groups have fixed them, whereas four others left them free.

In the next two sections, we present the results of one of the fits used as reference. The reference fitting is based on MLE/MAP, and the fitted frequency range is [1430, 2610] µHz. The background component $B_a$ is fixed, whereas Bayesian priors for $B_g$ and $W$ are derived from the background fit (Sect. 5). There are no Bayesian priors for the mode parameters. The marginal probabilities derived from MCMC derived by other groups are compatible with the error bars presented hereafter.

## 6.2. Mode frequencies

For the two fitting configurations (A or B), the fitted frequencies are in very good agreement. Table 4 provides modes frequencies obtained with fit A. The fitted spectrum is plotted over the observation in Fig. 8. The determined frequencies are also plotted over the échelle diagram (Fig. 9). The frequencies listed in the table were found by at least eight groups out of ten within the error bars. The three modes labelled with an exponent *a* in the table correspond to modes fitted by only four groups. Nevertheless, the four groups have independently obtained consistent frequencies for these modes, so we report them, but they could be less reliable. An extra *l* = 1, around 1455 µHz, seems to be present in the échelle diagram and has been fitted by a few groups. Nevertheless, it is very close to one of the orbital harmonics (see Fig. 9) and has been rejected.

The radial order *n* provided in Table 4 is only relative and could be shifted by ±1. It is obtained by fitting the relation

$$v_{l,n} = \overline{\Delta v}(n + l/2 + 1/4 + \alpha) - l(l + 1)\overline{\delta}_{02}/6 \qquad (7)$$

derived from the asymptotic development for p modes (e.g. Tassoul 1980).

<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Large separations $\Delta v_{l,n} = v_{l,n} - v_{l,n-1}$ for *l* = 0, 1, and 2 modes plotted as a function of the frequency $(v_{l,n} + v_{l,n-1})/2$ for HD 52265 (symbols with error bars). A solid line shows the variations in $\Delta v$ obtained with the envelope autocorrelation function. The grey area indicates 1-$\sigma$ error bars. The inset shows the cosine filter with a full width at half maximum equal to 2 times the mean large separation used for computing the EACF.</figcaption></figure>

<table><caption>**Table 4.** Fitted mode frequencies $v_{l,n}$ for HD 52265.</caption><thead><tr><th>n</th><th>l</th><th>v<sub>l,n</sub><br>[µHz]</th><th>l</th><th>v<sub>l,n</sub><br>[µHz]</th><th>l</th><th>v<sub>l,n</sub><br>[µHz]</th></tr></thead><tbody><tr><td>14</td><td>0</td><td>1509.17 ± 0.06<sup>a</sup></td><td>1</td><td>1554.33 ± 0.42<sup>a</sup></td><td>...</td><td>...</td></tr><tr><td>15</td><td>0</td><td>1606.54 ± 0.33</td><td>1</td><td>1652.80 ± 0.34</td><td>2</td><td>1696.88 ± 0.53<sup>a</sup></td></tr><tr><td>16</td><td>0</td><td>1704.88 ± 0.33</td><td>1</td><td>1749.85 ± 0.48</td><td>2</td><td>1793.02 ± 1.13</td></tr><tr><td>17</td><td>0</td><td>1801.15 ± 0.56</td><td>1</td><td>1845.74 ± 0.51</td><td>2</td><td>1890.81 ± 0.77</td></tr><tr><td>18</td><td>0</td><td>1898.22 ± 0.38</td><td>1</td><td>1943.94 ± 0.25</td><td>2</td><td>1988.38 ± 0.40</td></tr><tr><td>19</td><td>0</td><td>1996.32 ± 0.20</td><td>1</td><td>2041.81 ± 0.23</td><td>2</td><td>2086.48 ± 0.45</td></tr><tr><td>20</td><td>0</td><td>2094.92 ± 0.25</td><td>1</td><td>2141.32 ± 0.20</td><td>2</td><td>2186.20 ± 0.31</td></tr><tr><td>21</td><td>0</td><td>2193.75 ± 0.23</td><td>1</td><td>2240.31 ± 0.26</td><td>2</td><td>2284.02 ± 0.33</td></tr><tr><td>22</td><td>0</td><td>2292.86 ± 0.24</td><td>1</td><td>2338.11 ± 0.38</td><td>2</td><td>2382.56 ± 0.77</td></tr><tr><td>23</td><td>0</td><td>2389.78 ± 0.76</td><td>1</td><td>2437.27 ± 0.44</td><td>2</td><td>2479.76 ± 1.14</td></tr><tr><td>24</td><td>0</td><td>2489.85 ± 0.71</td><td>1</td><td>2536.07 ± 1.19</td><td>...</td><td>...</td></tr></tbody></table>

**Notes.** (a) Less reliable frequencies. See text for details.

By using these frequency determinations, we can compute and plot large separations

$$\Delta v_{l,n} = v_{l,n} - v_{l,n-1} \qquad (8)$$

in Fig. 10. There are clear variations around the mean value $\Delta v = 98.3 \pm 0.1$ µHz, recomputed as the average of $\Delta v_{l,n}$. The variation of the large separation can also be measured from the autocorrelation of the time series computed as the power spectrum of the power spectrum windowed with a narrow filter, as proposed by Roxburgh & Vorontsov (2006) and Roxburgh (2009a). Mosser & Appourchaux (2009) have defined the envelope autocorrelation function (EACF) and explain its use as an automated pipeline. We applied this pipeline to the spectrum of HD 52265 to find the frequency variation $\Delta v(v)$ of the large separation. We used a narrow cosine filter, with a full width at half maximum equal to $2\Delta v$. The result is plotted in Fig. 10. The agreement between $\Delta v(v)$ and the fitted values $\Delta v_{l,n}$ are quite good. We recover the main variations, particularly the large oscillation, visible as a bump around 2100 µHz. Such oscillations in the frequencies and seismic variables can be created by sharp features in the stellar structure (e.g. Vorontsov 1988; Gough 1990). The observed oscillations are probably the signature of
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Small separations $\delta_{02,n}$ plotted as a function of the frequencies $v_{0,n}$ for HD 52265. The dashed line corresponds to the mean value.</figcaption></figure>

the helium's second ionisation zone located below the surface of the star. Such a signature, predicted by models and observed on the Sun, should help us to put constraints on the helium abundance of the star (e.g. Basu et al. 2004; Piau et al. 2005).

We then computed small separations

$$ \delta_{02,n} = v_{0,n} - v_{2,n-1}. \qquad (9) $$

Results are shown in Fig. 11. The error bars take the correlations between the frequency determination of the $l = 0$ and $2$ modes into account. These correlations are small. It is worth noticing that these small separations do not decrease with frequency as for the Sun, but remain close to their mean value, $\delta_{02} = 8.1 \pm 0.2$ µHz, re-evaluated from fitted frequencies.

The small separations between $l = 0$ and $1$ modes, defined by Roxburgh (1993) as

$$ \delta_{01,n} = v_{0,n} - \frac{1}{2}(v_{1,n} + v_{1,n-1}) \quad \text{and} \qquad (10) $$

$$ \delta_{10,n} = \frac{1}{2}(v_{0,n} + v_{0,n+1}) - v_{1,n}, \qquad (11) $$

are also interesting seismic diagnosis tools (e.g. Roxburgh & Vorontsov 2003; Roxburgh 2005). Results for HD 52265 are shown in Fig. 12. As with the $\delta_{02}$ separations they do not display the decrease in frequency seen in the solar values, and have an average of $3.28 \pm 0.09$ µHz. These are additional diagnostics of the interior and may even be showing an oscillation due to the radius of an outer convective zone and the signature of a convective core (e.g. Roxburgh 2009b; Silva Aguirre et al. 2011).

The results for the parameters $v_s$ and $i$ obtained by the different groups are spread but consistent. They appear to be sensitive to the prior for the background and the frequency range of the fit. A detailed study dedicated to determining $v_s$ and $i$ will be presented in a separate paper (Gizon et al., in prep.). The results for the projected splitting $v_s^* = v_s \sin i$, which is easier to determine than $v_s$ and $i$ separately (Ballot et al. 2006), are also consistent among the different groups. For the analysis presented in this work we find $v_s^* = 0.45 \pm 0.10$ µHz for HD 52265. Moreover, Ballot et al. (2008) show that the frequency estimates are not correlated with estimates for other parameters, in particular $v_s$ and $i$, so, the frequency estimates are robust.

<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Small separations $\delta_{01,n}$ and $\delta_{10,n}$ plotted as a function of the frequencies $v_{0,n}$ and $v_{1,n}$ for HD 52265. The dashed line corresponds to the mean value.</figcaption></figure>

<table><caption>**Table 5.** Fitted widths ($\Gamma$) and amplitudes (A) of modes ($l = 0, n$) for HD 52265, obtained with fits A and B.</caption><thead><tr><th>Fit</th><th>n</th><th>Γ<br>[μHz]</th><th>+err</th><th>-err</th><th>A<br>[ppm]</th><th>+err</th><th>-err</th></tr></thead><tbody><tr><td>A</td><td>14</td><td>0.11<sup>a</sup></td><td>+0.22</td><td>-0.08</td><td>1.34<sup>a</sup></td><td>+0.34</td><td>-0.27</td></tr><tr><td>B</td><td>14</td><td>0.10<sup>a</sup></td><td>+0.03</td><td>-0.02</td><td>1.24<sup>a</sup></td><td>+0.30</td><td>-0.24</td></tr><tr><td>A</td><td>15</td><td>1.23</td><td>+1.45</td><td>-0.67</td><td>1.63</td><td>+0.29</td><td>-0.25</td></tr><tr><td>B</td><td>15</td><td>1.56</td><td>+1.25</td><td>-0.69</td><td>1.81</td><td>+0.27</td><td>-0.23</td></tr><tr><td>A</td><td>16</td><td>1.45</td><td>+0.68</td><td>-0.46</td><td>1.97</td><td>+0.22</td><td>-0.20</td></tr><tr><td>B</td><td>16</td><td>1.46</td><td>+0.61</td><td>-0.43</td><td>2.00</td><td>+0.23</td><td>-0.20</td></tr><tr><td>A</td><td>17</td><td>2.96</td><td>+1.31</td><td>-0.91</td><td>2.38</td><td>+0.25</td><td>-0.23</td></tr><tr><td>B</td><td>17</td><td>3.87</td><td>+0.92</td><td>-0.74</td><td>2.88</td><td>+0.22</td><td>-0.21</td></tr><tr><td>A</td><td>18</td><td>3.14</td><td>+0.69</td><td>-0.57</td><td>3.11</td><td>+0.22</td><td>-0.20</td></tr><tr><td>B</td><td>18</td><td>2.08</td><td>+0.55</td><td>-0.43</td><td>2.99</td><td>+0.22</td><td>-0.21</td></tr><tr><td>A</td><td>19</td><td>1.53</td><td>+0.40</td><td>-0.32</td><td>3.07</td><td>+0.22</td><td>-0.21</td></tr><tr><td>B</td><td>19</td><td>1.33</td><td>+0.30</td><td>-0.25</td><td>3.26</td><td>+0.23</td><td>-0.22</td></tr><tr><td>A</td><td>20</td><td>1.91</td><td>+0.38</td><td>-0.32</td><td>3.72</td><td>+0.24</td><td>-0.22</td></tr><tr><td>B</td><td>20</td><td>2.03</td><td>+0.42</td><td>-0.35</td><td>3.86</td><td>+0.24</td><td>-0.23</td></tr><tr><td>A</td><td>21</td><td>1.61</td><td>+0.36</td><td>-0.30</td><td>3.69</td><td>+0.24</td><td>-0.22</td></tr><tr><td>B</td><td>21</td><td>1.86</td><td>+0.38</td><td>-0.31</td><td>3.68</td><td>+0.24</td><td>-0.22</td></tr><tr><td>A</td><td>22</td><td>1.92</td><td>+0.38</td><td>-0.31</td><td>3.70</td><td>+0.23</td><td>-0.22</td></tr><tr><td>B</td><td>22</td><td>2.29</td><td>+0.42</td><td>-0.36</td><td>3.74</td><td>+0.23</td><td>-0.22</td></tr><tr><td>A</td><td>23</td><td>4.02</td><td>+0.76</td><td>-0.64</td><td>3.58</td><td>+0.22</td><td>-0.20</td></tr><tr><td>B</td><td>23</td><td>5.08</td><td>+1.15</td><td>-0.94</td><td>3.51</td><td>+0.22</td><td>-0.21</td></tr><tr><td>A</td><td>24</td><td>4.39</td><td>+1.21</td><td>-0.95</td><td>2.84</td><td>+0.21</td><td>-0.20</td></tr><tr><td>B</td><td>24</td><td>5.88</td><td>+2.56</td><td>-1.78</td><td>2.30</td><td>+0.24</td><td>-0.22</td></tr></tbody></table>

**Notes.** (a) Less reliable determinations. See text for details.

## 6.3. Lifetimes and amplitudes of modes

To consider the window effects that affect p modes, we divide all the fitted mode heights by the factor $r_{dc}$, as mentioned in Sect. 3.2. As expected, the values of the fitted heights and widths $H_{0,n}$ and $\Gamma_{0,n}$ change slightly depending on the fitting configuration (A or B) since the values of fitted height and width average the contribution of the modes ($l = 0, n$) and ($l = 1, n-1$) in case A, or of the modes ($l = 0, n$) and ($l = 1, n$) in case B.

The mode width is a direct measurement of the mode lifetime $\tau_{0,n} = (\pi\Gamma_{0,n})^{-1}$. The widths of fitted modes are listed in Table 5 and plotted in Fig. 13. The two lowest values are quite low and close to the spectral resolution. Moreover, the values of fitted widths obtained for these modes by the different groups are rather widespread, compared to the other widths, which are
<figure><img src="image_14.png" /><figcaption>Fig. 13. Full width at half maximum of modes as a function of frequency for HD 52265. Crosses correspond to values obtained by fitting the same widths for modes ($l = 1, n - 1$) and ($l = 0, n$), i.e. fit A, whereas plus signs correspond to values obtained by fitting the same widths for modes ($l = 1, n$) and ($l = 0, n$), i.e. fit B. Widths are plotted as a function of the mean frequency of the two modes. The dotted line indicates the spectral resolution. The dashed line shows predictions by Chaplin et al. (2009). The y-axis on the right-hand side shows the corresponding lifetimes.</figcaption></figure>

consistent. As a consequence, these two widths are not considered as reliable and must be rejected. If we exclude these two points, the mode widths correspond to lifetimes ranging from ~0.5 days at the highest frequency to ~3 days at the lowest. The variation of widths is not monotonic and shows an S shape. The widths progressively increase until a local maximum of ~4 µHz ($\tau$ ~ 1 days) around 1850 µHz corresponding to the mode ($l = 1, n = 17$) that actually appears to be significantly wider than its neighbours. Then, the widths decrease until a plateau of ~2 µHz ($\tau$ ~ 2 days) covering ~4 orders, and finally increase again.

HD 52265 mode lifetimes are shorter than the solar ones, but significantly longer than lifetimes previously observed in F stars, such as HD 49933 (e.g. Appourchaux et al. 2008). Mode lifetimes clearly decrease as the effective temperature increases. By using theoretical models and former observations, Chaplin et al. (2009) find that the mean lifetimes of the most excited modes scale as $T_{eff}^{-4}$, which leads to $\tau \approx 2.6$ days for $T_{eff} = 6100$ K. This value is indicated in Fig. 13. We notice a qualitative agreement with the observations, but the lifetimes in the plateau are all shorter than the value deduced from this scaling law. Recently, Baudin et al. (2011) have used a sample of more accurate seismic observations – including these observations of HD 52265 – to show that the dependence of lifetimes on the effective temperature is even stronger and that it scales with $T_{eff}^{-14}$.

From $H_{0,n}$ and $\Gamma_{0,n}$, it is possible to recover the rms amplitude $A_{0,n}$ of a mode. Amplitudes are always better determined than heights and widths themselves. This is true regardless of the values of $i$ and $v_s$ (Ballot et al. 2008). The amplitudes follow the relation

$$A_{0,n} = \sqrt{\frac{\pi}{2}\Gamma_{0,n}H_{0,n}} \qquad (12)$$

The mode amplitudes are listed in Table 5 and plotted in Fig. 14. The errors take the correlations between the determinations of $H_{0,n}$ and $\Gamma_{0,n}$ into account. The amplitudes increase almost regularly until reaching a maximum (3.86 ppm) around the frequency

<figure><img src="image_15.png" /><figcaption>Fig. 14. Amplitudes of radial modes as a function of frequency for HD 52265. Crosses correspond to values obtained by fitting the same intrinsic amplitudes for modes ($l = 1, n - 1$) and ($l = 0, n$), i.e. fit A, whereas plus signs correspond to values obtained by fitting same intrinsic amplitudes for modes ($l = 1, n$) and ($l = 0, n$), i.e. fit B. Amplitudes are plotted as a function of the mean frequency of the two modes. The dashed line shows the amplitudes of radial modes deduced from the fit as a Gaussian profile of the p-mode power excess (see Sect. 5). Grey area indicates the associated $1-\sigma$ error bars. The green curve shows the amplitudes of radial modes deduced from the spectrum smoothed by a $\Delta\nu$-wide boxcar, after subtracting the background.</figcaption></figure>

2100 µHz. Then, the mode amplitudes remain close to this maximum before sharply dropping above 2450 µHz.

It is also possible to derive the amplitude of radial modes from the smoothed power spectrum after subtracting the background (see Kjeldsen et al. 2008), or, similarly, by using the fitted Gaussian profile $B_p(v)$ of the p-mode power excess derived in Sect. 5. By using Michel et al. (2009), we recover the amplitude of radial modes through a relation

$$A = \frac{R_{l=0}}{R_{\text{osc}}} \sqrt{B_p \Delta\nu}, \qquad (13)$$

where $R_{l=0}$ and $R_{osc}$ are the CoRoT response functions for $l = 0$ modes and for the sum of all modes of degree $l = 0$ to 4 (see definitions in Michel et al. 2009, notice that $R_{l=0}$ is identical to the granulation response $R_g$ defined in that article). We used the approximated formulations of the responses they derived, and find $R_{l=0} = 3.89$ and $R_{osc} = 6.86$, by considering $T_{eff} \approx 6100$ K for HD 52265.

When the spectrum is barely resolved, this procedure is a common way to extract the mode amplitudes. In the present situation, we are thus able to compare these estimates to the individually fitted mode amplitudes. In Fig. 14, we have also plotted the radial mode amplitudes deduced from the smoothed spectrum and from the profile $B_p$. Although the amplitudes of radial modes do not follow a perfect Gaussian profile, the Gaussian profile is a reasonable fit to the smoothed spectrum. The Gaussian profile reaches the maximum value $A^{(\text{max})} = \frac{R_{l=0}}{R_{\text{osc}}} \sqrt{H_p^{(\text{max})} \Delta\nu} \approx 3.7$ ppm. This is consistent with the maximum fitted amplitude ($A^{(\text{max})}_{l=0} = 3.86 \pm 0.24$ ppm). We get very good agreement at high frequency, but there is a clear departure at low frequency between the amplitudes obtained from the smoothed spectrum and from the fitted modes. The fitted background is probably underestimated around ~1200–1800 µHz. This could come from a
<table><caption>**Table 6.** Fitted amplitude mode ratios $r_l = A_l/A_0$ compared to theoretical values.</caption><thead><tr><th></th><th>Fit A</th><th>Fit B</th><th>Theo.</th></tr></thead><tbody><tr><td>$r_1$</td><td>1.22 ± 0.06</td><td>1.21 ± 0.06</td><td>1.22</td></tr><tr><td>$r_2$</td><td>0.76 ± 0.05</td><td>0.74 ± 0.05</td><td>0.72</td></tr><tr><td>$1 + r_1^2 + r_2^2$</td><td>3.07 ± 0.19</td><td>3.02 ± 0.19</td><td>3.01</td></tr></tbody></table>

component missing in our model, like a contribution of faculae or from an excess of power due to leakage effects. Thus, when the granulation is fitted in Sect. 5, the extra power increases the apparent contribution from p modes; however, when we fit p modes, in a restricted range (Sect. 6.1), the extra power excess is included in the background.

The amplitudes are measured in the CoRoT spectral band and must be converted into bolometric amplitudes. According to Michel et al. (2009), the bolometric correction factor is $c_{bol} = 4/R_{l=0}$ and is $c_{bol} = 1.03$ for HD 52265. The maximum bolometric amplitude of radial mode is then

$$A_{l=0,\text{bol}}^{(\text{max})} = 3.96 \pm 0.24 \text{ ppm.} \qquad (14)$$

By combining the adiabatic relation proposed by Kjeldsen & Bedding (1995) to relate mode amplitudes in intensity to mode amplitudes in velocity and the scaling law proposed by Samadi et al. (2007b), we get the relation

$$A^{(\text{max})} = A_{\odot}^{(\text{max})} \left( \frac{L/L_{\odot}}{M/M_{\odot}} \right)^{0.7} \sqrt{\frac{T_{\text{eff}}}{T_{\text{eff},\odot}}}. \qquad (15)$$

Moreover, if we assume that the frequency of maximum p-mode amplitude $\nu^{(\text{max})}$ scales with the acoustic cut-off frequency (e.g., see Bedding & Kjeldsen 2003), we obtain the relation (Deheuvels et al. 2010):

$$A^{(\text{max})} = A_{\odot}^{(\text{max})} \left( \frac{T_{\text{eff}}}{T_{\text{eff},\odot}} \right)^{1.95} \left( \frac{\nu^{(\text{max})}}{\nu_{\odot}^{(\text{max})}} \right)^{-0.7}. \qquad (16)$$

By using $A_{l=0,\text{bol},\odot}^{(\text{max})} = 2.53 \pm 0.11 \text{ ppm}$ (Michel et al. 2009), $T_{\text{eff},\odot} = 5777$ K, and $\nu_{\odot}^{(\text{max})} = 3050$ µHz for the Sun and by considering for HD 52265 our estimations $T_{\text{eff}} = 6100 \pm 60$ K, and $\nu^{(\text{max})} = 2090 \pm 20$ µHz, we obtain $A_{l=0,\text{bol}}^{(\text{max})} = 3.7 \pm 0.2 \text{ ppm}$. This predicted value is slightly lower than the value measured for HD 52265, but still consistent within the error bars. In contrast to previous observations of F stars (e.g. Michel et al. 2008) which have smaller amplitudes than expected, the observations for this G0V star are close to predictions.

The mode visibility $r_1$ and $r_2$ are free parameters in the reference fitting. The values obtained are given in Table 6 for fits A and B. Both results are consistent. These values are compatible with theoretical values computed with the CoRoT limb-darkening laws by Sing (2010). As a further test, we considered the sum of visibilities of modes over an interval $\Delta\nu$ that is equal to $1 + r_1^2 + r_2^2$. The fitted values are listed in Table 6. According to Michel et al. (2009), this is comparable to the quantity $(R_{osc}/R_{l=0})^2 = 3.10$. The second includes modes up to $l = 4$, whereas the first includes modes only up to $l = 2$. Nevertheless, the agreement between fits A and B and the theoretical expectation is met within the error bars.

## 7. Conclusion

Using photometric observations from the CoRoT space telescope spanning 117 days and a duty cycle of 90%, we measured a rotation period $P_{\text{rot}} = 12.3 \pm 0.15$ days for the G0 main sequence star HD 52265, thanks to a modulation in the light curve induced by photospheric activity. We have clearly detected solar-like oscillations in HD 52265, and characterised 31 p modes in the range 1500–2550 µHz with degrees $l = 0, 1$, and 2. We observed lifetimes for these modes ranging between 0.5 and 3 days. HD 52265 mode lifetimes are shorter than the solar ones, but significantly longer than the lifetimes previously observed in F stars, confirming that mode lifetimes decrease as the effective temperature increases. Moreover, we observed for HD 52265 that the variation in lifetimes with frequency is not monotonic and shows a clear S shape.

The fitted maximum bolometric amplitude for radial modes is $3.96 \pm 0.24 \text{ ppm}$, which is marginally higher than the theoretical models of Samadi et al. (2007b) but still compatible within the error bars. In the past, several analyses have shown smaller amplitudes than predicted; however, they were F stars, whereas HD 52265 is a G0 star, hence more like the Sun. Nevertheless, this star is over metallic. Samadi et al. (2010b) have shown that the surface metallicity has a strong impact on the efficiency of the mode driven by turbulent convection: the lower the surface metal abundance, the weaker the driving. These authors have precisely quantified this effect for the CoRoT target HD 49933, which is a rather metal-poor star compared to the Sun since for this star [Fe/H] = -0.37. Samadi et al. (2010a) found that ignoring the surface metal abundance of this target results in a significant underestimation of the observed mode amplitudes. The theoretical scaling law by Samadi et al. (2007b) (see Eq. (15)) was obtained on the basis of a series of 3D hydrodynamical models with a solar metal abundance. Given the result of Samadi et al. (2010a,b), we would expect for HD 52265 higher theoretical mode amplitudes than predicted by the theoretical scaling law of Samadi et al. (2007b). However, the amount by which the theoretical amplitudes are expected to increase remains to be precisely quantified and compared to the uncertainties associated with the present seismic data.

For HD 52265, we found that the mean large and small separations are $\bar{\Delta}\nu = 98.3 \pm 0.1$ µHz and $\bar{\delta}_{02} = 8.1 \pm 0.2$ µHz and that $\delta_{02}$ does not significantly decrease with frequency. These quantities are typical of a ~1.2-solar-mass star that is still on the main sequence (see preparatory models by Soriano et al. 2007). Moreover, HD 52265 does not show any mixed modes, as is the case for the more evolved G0 star HD 49385 (Deheuvels et al. 2010).

The variation of $\bar{\Delta}\nu$ with frequency shows an oscillation that we interpret as a possible signature of the second helium-ionisation region. Thanks to accurate eigenfrequency measurements (the error is about 0.2 µHz for the modes with the highest amplitudes) and to fine measurements of fundamental parameters obtained with the Narval spectrograph, HD 52265 is a very promising object for stellar modelling. These observations should especially help in determining whether HD 52265 has a convective core or not, which would put useful constraints on the mixing processes in such stars. These seismic inferences will also help improve our knowledge of the planet hosted by HD 52265.

Acknowledgements. J.B. acknowledges the support of the Agence National de la Recherche through the SIROCO project. I.W.R., G.A.V., W.J.C., and Y.E. acknowledge support from the UK Science and Technology Facilities Council
(STFC). This work has been partially supported by the CNES/GOLF grant at the Service d'Astrophysique (CEA/Saclay) and the grant PNAyA2007-62650 from the Spanish National Research Plan. NCAR is supported by the National Science Foundation. This work benefited from the support of the International Space Science Institute (ISSI), by funding the AstroFLAG international team. It was also partly supported by the European Helio- and Asteroseismology Network (HELAS), a major international collaboration funded by the European Commission's Sixth Framework Programme, and by the French PNPS programme. This research made use of the SIMBAD database, operated at the CDS, Strasbourg, France.

## References

Appourchaux, T., Gizon, L., & Rabello-Soares, M. 1998, A&AS, 132, 107

Appourchaux, T., Berthomieu, G., Michel, E., et al. 2006, in The CoRoT Mission: Pre-Launch Status, ed. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, ESA Publications Division, Noordwijk, ESA SP, 1306, 377

Appourchaux, T., Michel, E., Auvergne, M., et al. 2008, A&A, 488, 705

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

Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411

Baglin, A., Auvergne, M., Barge, P., et al. 2006, in The CoRoT Mission: Pre-Launch Status, ed. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, ESA Publications Division, Noordwijk, ESA SP, 1306, 33

Ballot, J., García, R. A., & Lambert, P. 2006, MNRAS, 369, 1281

Ballot, J., Appourchaux, T., Toutain, T., & Guittet, M. 2008, A&A, 486, 867

Barban, C., Deheuvels, S., Baudin, F., et al. 2009, A&A, 506, 51

Basu, S., Mazumdar, A., Antia, H. M., & Demarque, P. 2004, MNRAS, 350, 277

Baudin, F., Barban, C., Belkacem, K., et al. 2011, A&A, 529, A84

Bazot, M., Vauclair, S., Bouchy, F., & Santos, N. C. 2005, A&A, 440, 615

Bedding, T. R., & Kjeldsen, H. 2003, PASA, 20, 203

Benomar, O., Appourchaux, T., & Baudin, F. 2009a, A&A, 506, 15

Benomar, O., Baudin, F., Campante, T. L., et al. 2009b, A&A, 507, L13

Bouchy, F., Bazot, M., Santos, N. C., Vauclair, S., & Sosnowska, D. 2005, A&A, 440, 609

Bruntt, H. 2009, A&A, 506, 235

Bruntt, H., Bikmaev, I. F., Catala, C., et al. 2004, A&A, 425, 683

Bruntt, H., De Cat, P., & Aerts, C. 2008, A&A, 478, 487

Bruntt, H., Bedding, T. R., Quirion, P., et al. 2010, MNRAS, 405, 1907

Butler, R. P., Vogt, S. S., Marcy, G. W., et al. 2000, ApJ, 545, 504

Cayrel, R., Castelli, F., Katz, D., et al. 1997, in Proc. of the ESA Symposium Hipparcos - Venise '97, ed. M. A. C. Perryman, P. L. Bernacca, & B. Battrick, ESA Publications Division, Noordwijk, ESA SP, 402, 433

Chaplin, W. J., Houdek, G., Karoff, C., Elsworth, Y., & New, R. 2009, A&A, 500, L21

Chaplin, W. J., Appourchaux, T., Elsworth, Y., et al. 2010, ApJ, 713, L169

Christensen-Dalsgaard, J., Kjeldsen, H., Brown, T. M., et al. 2010, ApJ, 713, L164

De Ridder, J., Barban, C., Baudin, F., et al. 2009, Nature, 459, 398

Deheuvels, S., & Michel, E. 2010, Ap&SS, 328, 259

Deheuvels, S., Bruntt, H., Michel, E., et al. 2010, A&A, 515, A87

Duvall, Jr., T. L., & Harvey, J. W. 1986, in NATO ASIC Proc., Seismology of the Sun and the Distant Stars, ed. D. O. Gough (Dordrecht: D. Reidel Publishing Co.), 169, 105

Flower, P. J. 1996, ApJ, 469, 355

García, R. A., Régulo, C., Samadi, R., et al. 2009, A&A, 506, 41

García, R. A., Mathur, S., Salabert, D., et al. 2010, Science, 329, 1032

Gaulme, P., Appourchaux, T., & Boumier, P. 2009, A&A, 506, 7

Gaulme, P., Deheuvels, S., Weiss, W. W., et al. 2010, A&A, 524, A47

Ghezzi, L., Cunha, K., Smith, V. V., et al. 2009, ApJ, 698, 451

Gilliland, R. L., Brown, T. M., Christensen-Dalsgaard, J., et al. 2010, PASP, 122, 131

Gizon, L., & Solanki, S. K. 2003, ApJ, 589, 1009

Goldreich, P., & Keeley, D. A. 1977, ApJ, 212, 243

Gough, D. O. 1990, in Progress of Seismology of the Sun and Stars, ed. Y. Osaki, & H. Shibahashi (Berlin, Heidelberg, New-York: Springer-Verlag), Lect. Notes Phys., 367, 283

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

Handberg, R., & Campante, T. L. 2011, A&A, 527, A56

Harvey, J. 1985, in Future Missions in Solar, Heliospheric & Space Plasma Physics, ed. E. Rolfe, & B. Battrick, ESA Publications Division, Noordwijk, ESA SP, 235, 199

Hekker, S., Kallinger, T., Baudin, F., et al. 2009, A&A, 506, 465

Hinkle, K., Wallace, L., Valenti, J., & Harmer, D. 2000, Visible and Near Infrared Atlas of the Arcturus Spectrum 3727-9300 Å (San Francisco: ASP)

Hoffleit, D., & Jaschek, C. 1982, The Bright Star Catalogue, 4th ed. (New Haven: Yale University Observatory)

Huber, D., Bedding, T. R., Stello, D., et al. 2010, ApJ, 723, 1607

Israelian, G., Santos, N. C., Mayor, M., & Rebolo, R. 2004, A&A, 414, 601

Kallinger, T., Mosser, B., Hekker, S., et al. 2010, A&A, 522, A1

Kjeldsen, H., & Bedding, T. R. 1995, A&A, 293, 87

Kjeldsen, H., Bedding, T. R., Arentoft, T., et al. 2008, ApJ, 682, 1370

Koch, D. G., Borucki, W. J., Basri, G., et al. 2010, ApJ, 713, L79

Kupka, F., Piskunov, N., Ryabchikova, T. A., Stempels, H. C., & Weiss, W. W. 1999, A&AS, 138, 119

Ludwig, H., Samadi, R., Steffen, M., et al. 2009, A&A, 506, 167

Mathur, S., García, R. A., Catala, C., et al. 2010a, A&A, 518, A53

Mathur, S., García, R. A., Régulo, C., et al. 2010b, A&A, 511, A46

Mathur, S., Handberg, R., Campante, T. L., et al. 2011, ApJ, accepted [arXiv:1103.4085]

Michel, E., Baglin, A., Auvergne, M., et al. 2008, Science, 322, 558

Michel, E., Samadi, R., Baudin, F., et al. 2009, A&A, 495, 979

Miglio, A., Montalbán, J., Carrier, F., et al. 2010, A&A, 520, L6

Mosser, B., & Appourchaux, T. 2009, A&A, 508, 877

Mosser, B., Baudin, F., Lanza, A. F., et al. 2009a, A&A, 506, 245

Mosser, B., Michel, E., Appourchaux, T., et al. 2009b, A&A, 506, 33

Mosser, B., Belkacem, K., Goupil, M., et al. 2010, A&A, 517, A22

Naef, D., Mayor, M., Pepe, F., et al. 2001, A&A, 375, 205

Piau, L., Ballot, J., & Turck-Chièze, S. 2005, A&A, 430, 571

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in FORTRAN, The art of scientific computing, 2nd ed. (Cambridge: Cambridge University Press)

Roxburgh, I. W. 1993, in PRISMA, Report on the Phase-A Study, ed. T. Appourchaux et al., ESA SCI 93(3) (Paris: ESA), 31

Roxburgh, I. W. 2005, A&A, 434, 665

Roxburgh, I. W. 2009a, A&A, 506, 435

Roxburgh, I. W. 2009b, A&A, 493, 185

Roxburgh, I. W., & Vorontsov, S. V. 2003, A&A, 411, 215

Roxburgh, I. W., & Vorontsov, S. V. 2006, MNRAS, 369, 1491

Samadi, R., Fialho, F., Costa, J. E. S., et al. 2007a [arXiv:astro-ph/0703354]

Samadi, R., Georgobiani, D., Trampedach, R., et al. 2007b, A&A, 463, 297

Samadi, R., Ludwig, H., Belkacem, K., et al. 2010a, A&A, 509, A16

Samadi, R., Ludwig, H., Belkacem, K., Goupil, M. J., & Dupret, M. 2010b, A&A, 509, A15

Sato, K. H., García, R. A., Pires, S., et al. 2010, Astron. Nachr., 931, P06

Silva Aguirre, V., Ballot, J., Serenelli, A., & Weiss, A. 2011, A&A, 529, A63

Sing, D. K. 2010, A&A, 510, A21

Soriano, M., & Vauclair, S. 2010, A&A, 513, A49

Soriano, M., Vauclair, S., Vauclair, G., & Laymand, M. 2007, A&A, 471, 885

Tassoul, M. 1980, ApJS, 43, 469

Torrence, C., & Compo, G. P. 1998, Bull. Am. Meteor. Soc., 79, 61

Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141

Valenti, J. A., & Piskunov, N. 1996, A&AS, 118, 595

van Leeuwen, F. 2007, Hipparcos, the New Reduction of the Raw Data (Dordrecht: Springer), Astrophys. Space Sci. Lib., 350

Vauclair, S., Laymand, M., Bouchy, F., et al. 2008, A&A, 482, L5

Vorontsov, S. V. 1988, in Advances in Helio- and Asteroseismology, ed. J. Christensen-Dalsgaard, & S. Frandsen (D. Reidel Publishing Co., Dordrecht), IAU Symp., 123, 151

Wright, J. T., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2004, ApJS, 152, 261