Issue 
A&A
Volume 570, October 2014



Article Number  A54  
Number of page(s)  16  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201424326  
Published online  15 October 2014 
Asteroseismic inference on the spinorbit misalignment and stellar parameters of HATP7^{⋆}
^{1} Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark
email: mikkelnl@phys.au.dk
^{2} Sydney Institute for Astronomy (SIfA), School of Physics, University of Sydney, NSW 2006, Australia
^{3} Research School of Astronomy & Astrophysics, Mount Stromlo Observatory, The Australian National University, ACT 2611, Australia
^{4} School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
^{5} Department of Geoscience, Aarhus University, HøeghGuldbergs Gade 2, 8000, Aarhus C, Denmark
^{6} Institut für Astrophysik, GeorgAugustUniversität Göttingen, FriedrichHundPlatz 1, 37077 Göttingen, Germany
^{7} MaxPlanckInstitut für Sonnensystemforschung, JustusvonLiebigWeg 3, 37077 Göttingen, Germany
^{8} Instituut voor Sterrenkunde, Katholieke Universiteit Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
Received: 2 June 2014
Accepted: 28 July 2014
Context. The measurement of obliquities – the angle between the orbital and stellar rotation – in starplanet systems is of great importance for understanding planet system formation and evolution. The bright and wellstudied HATP7 (Kepler2) system is intriguing because several RossiterMcLaughlin (RM) measurements found a high projected obliquity in this system, but it was not possible so far to determine whether the orbit is polar and/or retrograde.
Aims. The goal of this study is to measure the stellar inclination and hereby the full 3D obliquity of the HATP7 system instead of only the 2D projection as measured by the RM effect. In addition, we provide an updated set of stellar parameters for the star.
Methods. We used the full set of available observations from Kepler spanning Q0Q17 to produce the power spectrum of HATP7. We extracted oscillationmode frequencies via an Markov chain Monte Carlo peakbagging routine and used the results from this to estimate the stellar inclination angle. Combining this with the projected obliquity from RM and the inclination of the orbital plane allowed us to determine the stellar obliquity. Furthermore, we used asteroseismology to model the star from the extracted frequencies using two different approaches to the modelling, for which either the stellar evolution codes MESA or GARSTEC were adopted.
Results. Our updated asteroseismic modelling shows, i.a., the following stellar parameters for HATP7: M_{⋆} = 1.51^{+ 0.04}_{0.05} M_{⊙}, R_{⋆} = 2.00^{+ 0.01}_{0.02} R_{⊙}, and age = 2.07^{+ 0.28}_{0.23} Gyr. The modelling offers a high precision on the stellar parameters, the uncertainty on age, for instance, is of the order ∼ 11%. For the stellar inclination we estimate i_{⋆}< 36.5^{°}, which translates into an obliquity of 83^{°}<ψ< 111^{°}. The planet HATP7b is likely retrograde in its orbit, and the orbit is close to being polar. The new parameters for the star give an updated planetary density of ρ_{p} = 0.65 ± 0.03 g cm^{3}, which is lower than previous estimates.
Key words: asteroseismology / planetary systems / stars: oscillations / stars: individual: HATP7 / methods: data analysis / stars: rotation
Appendices are available in electronic form at http://www.aanda.org
© ESO, 2014
1. Introduction
Asteroseismology can provide detailed information about stellar parameters such as mass, radius, and age (see Chaplin & Miglio 2013, and references therein). Furthermore, an estimate for the stellar inclination can be obtained for solarlike oscillators (Gizon & Solanki 2003), which in turn is needed in asserting the obliquity of planethosting systems. The obliquity of planetary systems, ψ, which is the angle between the stellar spinaxis and the angular momentum vector of the planetary orbit, is an important parameter for a better understanding of how these systems form and evolve (see, e.g., Nagasawa et al. 2008; Winn et al. 2010; Triaud et al. 2010; Morton & Johnson 2011; Rogers et al. 2012). The obliquity is especially interesting for systems that appear to have retrograde orbits from measurements of the skyprojected obliquity, λ. The reason is that while the orbit might indeed be retrograde, it is not known whether it is closer to polar than equatorial. This distinction makes a great difference for theories dealing with planetary system formation and evolution because they must be able to account for such a configuration (see, e.g., Chatterjee et al. 2008; Ford & Rasio 2008; Matsumura et al. 2010; Lai et al. 2011; Albrecht et al. 2012).
Stellar obliquities are difficult to measure because stars are unresolved by modern telescopes. Therefore no spacial information can be obtained^{1}. During planetary transits, however, parts of the stellar surface are covered, which breaks the degeneracy. Such information on the obliquity can be obtained, for instance, from studies of the anomalous effect in the radial velocity (RV) curve that is known as the RossiterMcLaughlin (RM) effect (Rossiter 1924; McLaughlin 1924; see Fabrycky & Winn 2009 for an overview). Unfortunately, only λ can be obtained from RM measurements (see Fig. 1). Other means of obtaining λ are spotcrossing anomalies that are observed during planetary transits (see, e.g., SanchisOjeda et al. 2011, 2013; Désert et al. 2011), Doppler tomography (see, e.g., Gandolfi et al. 2012), or the effects of gravity darkening (see, e.g., Barnes et al. 2011; Ahlers et al. 2014), among others.
The true obliquity can only be unequivocally determined if the stellar angle of inclination, i_{⋆}^{2}, can be measured, and combined with λ and the inclination of the planetary orbital plane, i_{p}. The orbital inclination can be estimated fairly easily from analysis of the photometric light curve if the planet happens to transit its host star.
Fig. 1
Configuration of the HATP7 system to scale using values for HATP7b (full black circle) from Van Eylen et al. (2013). Top: the observer’s view of the system from Earth, with the angular momentum vectors of the planetary orbit, (normal to the orbital plane), and stellar spin, , given as red arrows. The projected angle, λ, is indicated by the shaded magenta region and is found as the angle between and when these are projected onto the plane of the sky (dashed lines). This angle is obtained from RM measurements, and in this panel we have used λ = 155 ± 37° following Albrecht et al. (2012; the uncertainty on λ is not included in the figure). The stellar inclination, i_{⋆}, which is the parameter measured from asteroseismology, is set to 15°, and is given by the direct angle between the line of sight (midpoint of star) and . The inclination of the planetary orbit, i_{p}, is set to 83° following Van Eylen et al. (2013). The true angle, ψ, is the direct angle between and . Middle: side view of the system, with the observer’s viewpoint from the right, indicated by “line of sight”. To properly show i_{⋆} (shaded green) and i_{p} (shaded red), and not their projected values, we have set λ = 180° such that both and lie on the same plane as the line of sight. For illustrative purposes we have in addition decreased i_{p} to 63°. When adopting this configuration ψ (shaded blue) is given by the sum of i_{⋆} and i_{p}. Bottom: top view of the system, with the observer’s viewpoint from the bottom, indicated by “line of sight”. 
A measure of the stellar inclination angle can be obtained from combining v sin i_{⋆} from spectroscopy with the stellar rotation period from modulations of the light curve from stellar spots (see Hirano et al. 2012, 2014, for recent uses of this method for planetary systems). These estimates can be quite uncertain because of the difficulty of calibrating the spectroscopic v sin i_{⋆}, separating the rotational signal from other broadening effects, and the need for an estimate of the stellar radius R_{⋆}. A more direct method for obtaining i_{⋆} is that of asteroseismology, where the stellar inclination can be estimated by analysing solarlike acoustic (pmode) oscillations (Gizon & Solanki 2003). Another great advantage of using asteroseismology is that a detailed stellar model can be obtained with welldetermined parameters that are needed in simulations of planetary systems dynamics.
The high photometric quality that enabled the Kepler mission (see Borucki et al. 2010; Koch et al. 2010) to detect the transits of extrasolar planets, and thus allows for a determination of i_{p}, also makes the data ideal for asteroseismic analysis (Gilliland et al. 2010). Analyses of the obliquity using the asteroseismic method have to date only been performed in the systems Kepler50 and 65 by Chaplin et al. (2013), Kepler56 by Huber et al. (2013a), Kepler410 by Van Eylen et al. (2014), and 16 Cygni by Davies et al. (2014). However, for these systems the obliquity could only be assessed in a statistical sense, since the projected angle, λ, was unavailable.
Our aim in this paper is to use asteroseismology to determine precise stellar parameters of HATP7 and to determine the obliquity of the HATP7^{3} system with a F6V (Faedi et al. 2013) type star and a closein ~ 1.78 M_{J} planet (HATP7b) in a ~ 2.2 day orbit (Pál et al. 2008). From the very onset it is clear, however, that this is a challenging task because HATP7 is a lateFtype star. This spectral type is notorious for having short lifetimes of the pmode oscillations and consequently very broad (in frequency) oscillation modes, which highly obscures the potentially small effects imposed by rotation.
A fortuitous feature of the system is that not only i_{⋆} and i_{p} can be estimated from Kepler data, but the RM effect has been studied independently by Winn et al. (2009), Narita et al. (2009), and Albrecht et al. (2012) using HIRES and/or HDS data^{4}. Values for λ and v sin i_{⋆} from these studies are given in Table 3. Despite disagreement on the actual values, all three works agree that the system is misaligned, that the planetary orbit might be retrograde, and that the very low value measured for v sin i_{⋆} suggests a low i_{⋆}, which for a transiting planet implies a nearpolar orbit of the planet.
An interesting aspect of the system with regard to the obliquity is that a third body (an M5.5V dwarf known as HATP7B) is found to be associated with the system (Narita et al. 2010, 2012; Bergfors et al. 2013; Faedi et al. 2013), and that a fourth associated body, likely more massive than Jupiter, is thought to be likely based on an unexplained RV excess (Winn et al. 2009; Narita et al. 2012). The system is thus a prime candidate for obliquity studies and theories concerned with planetary system formation and evolution.
The paper is structured as follows: in Sect. 2 we describe the data used in our analysis. Section 3 explains the model we used to obtain the stellar inclination angle. In Sect. 4 we present our analysis and results, including stellar modelling using two different codes (Sect. 4.2), and the results obtained for the stellar inclination and rotation (Sect. 4.3). Our results for the obliquity are the topic of Sect. 5. In Sect. 6 we compare our result on the stellar rotation rate with gyrochronology, and Sect. 7 is concerned with possible activity signatures. Finally, we discuss our results in Sect. 8 and conclude in Sect. 9.
Fig. 2
Left: power spectrum of HATP7 (black) with the fitted model overplotted (red; cf. Eq. (1)). Right: power spectrum of HATP7 (black) overplotted with the best fit to the background (red; cf. Eq. (2)). The lightgrey part up to 100 μHz was not included in the fit. The fit includes a granulation component (P_{G}; green) and a white/shot noise (P_{S}; black) level in addition to the Gaussian envelope from pmodes centred around ν_{max} ≈ 1115 μHz. The dashed red line shows the background fit without the Gaussian envelope. 
2. Data
We extracted shortcadence (SC; Δt = 58.8 s) simple aperture photometry (SAP) data from target pixel files (TPFs) using the procedure of Steven Bloemen (private communication), which starts from the original Kepler mask and adds to or removes pixels from the aperture based on the amount of signal in each pixel. The outcome of this is a new mask that most often is slightly larger than the original one.
The data span quarters from Q0 to Q17 (~ 1470 days), with a duty cycle of ~ 90.4%. This constitutes the full amount of data available from Kepler in its normal mode of operation. Data were downloaded from the Mikulski Archive for Space Telescopes (MAST) and corrected using the procedure described in Handberg & Lund (2014). Briefly, two medianfiltered versions of the time series are computed with different filter windows. A weighted combination of the two are then used together with a filtered version of the planetary phasecurve to correct the time series for both instrumental features and the planetary signal. For the asteroseismic analysis we used filter windows of 1 (τ_{long}) and 0.028 (τ_{short}) days, while windows of 15 (τ_{long}) and 5 (τ_{short}) days were adopted in preparing the time series used in Sect. 7.2 to search for a lowfrequency imprint of rotational modulation. We refer to Handberg & Lund (2014) for further details on the filter and for a view of the corrected time series for HATP7.
For the transit parameters needed in both the correction of the time series (orbital period) and for the estimation of the obliquity (inclination of orbital plane) we used the results of Van Eylen et al. (2013). The power spectrum was calculated using a weighted sinewave fitting method (see, e.g., Kjeldsen 1992; Frandsen et al. 1995), normalised according to the amplitudescaled version of Parseval’s theorem (see Kjeldsen & Frandsen 1992), in which a sine wave of peak amplitude, A, will have a corresponding peak in the power spectrum of A^{2}.
3. Fitting the power spectrum
The harmonic eigenmodes of acoustic solarlike oscillations are characterised by their degree, l, which gives the number of nodal lines on the stellar surface, and their radial order, n, giving the number of radial nodes. In addition, an eigenmode is characterised by its azimuthal order, m, of which there are 2l + 1. Only in the case of broken spherical symmetry, broken, for instance, by rotation, will the degeneracy between different mvalues be lifted. The removal of this degeneracy makes it possible to measure the stellar rotation rate and inclination angle.
3.1. Modelling the power spectrum
We modelled the power spectral density of the oscillations with a series of standard Lorentzian functions. The limit spectrum (noise free) to be fit to the power spectrum can be expressed as follows: (1)Here n_{a} and n_{b} represent the first and last mode orders included from the power spectrum, while ν_{nlm} is the mode frequency including the effect of rotation, B(ν_{j}) describes the contribution from the stellar noise background at frequency ν_{j}, S_{nl} is the overall height of the multiplet, that is, the maximum power spectral density, and Γ_{nl} is the mode line width. A geometrical modulation of the relative visibility between components of a split multiplet is given by ℰ_{lm}(i_{⋆}). The fitted parameters are denoted by Θ. This fit can be seen in the left panel of Fig. 2.
The background signal B(ν) is described by a series of power laws (Harvey 1985), each of which relate to a specific physical phenomenon. The power laws included in this work describe the signals from granulation and faculae. The specific functional from for the background signal is that suggested by Karoff (2008): (2)In this equation, σ_{i} gives the flux rms variation in time and τ_{i} the characteristic time scales of the different phenomena. The constant B_{0} is a measure of the photon shotnoise. The fit to the background can be seen in the right panel of Fig. 2.
3.2. Optimisation procedure
The fit of Eq. (1) to the power spectrum is optimised in a Bayesian manner using the Markov chain Monte Carlo (MCMC) sampler emcee (ForemanMackey et al. 2013).
Frequencies for all modes are free parameters in the fit. Heights and widths are only free parameters for radial (l = 0) modes. For a nonradial mode the width is found from a linear interpolation between the two nearest radial modes. The same is true for the heights where the linearly interpolated value for the nearest radial modes is scaled via the visibility parameter (assumed constant as a function of frequency). In this work we keep the relative visibilities as free parameters. With this setup we have the following set of free parameters in the fitting: .
We refer to Appendix A for details on the fitted model and the adopted optimisation.
4. Analysis and results
4.1. General results from peakbagging
To estimate mode frequencies for the modelling we fitted Eq. (1) to the frequency range 600 − 1650 μHz, which is the range where we can visually identify modes. All estimated frequencies are given in Table D.1. Here we report the median of the marginalised posterior (kernel) probability distributions (PPDs) for the respective modes, while the uncertainties were obtained from the 68% highest probability density (HPD) credible region. For the frequency uncertainties used in the stellar modelling (see Sect. 4.2) we adopted the mean value of the (potentially asymmetric) uncertainties from the HPD credible region. For results on mode line widths and visibilities we refer to Appendix B.
To estimate the inclination and splitting parameters we fitted Eq. (1) to a smaller range in the power spectrum that included only the frequencies in the range 780 − 1400 μHz. The selection of this interval was based on the estimates for the mode line widths from the large fit and to obtain modes with a relatively high signaltonoise ratio (S/N). For the splitting we used a flat prior from −8 to 8 μHz and then used the absolute value of the splitting in Eq. (1). For the inclination we used a flat prior from −90 to 180°, and then folded values onto the interval from 0 to 90°. The symmetry of the priors was chosen to avoid potential boundary effects in the MCMC sampling.
We note that a component of Eq. (1) that might cause problems in the fitting is the noisebackground. Initially, we did not fix the background in the fit of Eq. (1), but instead set Gaussian priors on the background parameters from the posteriors of a backgroundonly fit. However, given that the background is very poorly constrained in the relatively small part of the power spectrum occupied by the oscillation modes, we found that not even the Gaussian priors were able to constrain the background, and especially the granulation time scale, τ_{g}, wandered to lower values. These lower values for τ_{g} were found to correlate with a wide range of high values for the splitting at a particular value for the inclination. Because of this apparent degeneracy we chose to fix the background in the fits.
4.2. Modelling of HATP7
Selected properties of our bestfit models compared with the first asteroseismic modelling by ChristensenDalsgaard et al. (2010, CD10).
A detailed modelling of HATP7 was first made by ChristensenDalsgaard et al. (2010) based on asteroseismic measurements of the solarlike pmode oscillations in the star. This work was based on the SC Q0Q1 data from the Kepler satellite. A fit of 33 individual mode frequencies, with values obtained from peaks of a smoothed power spectrum, was made to models computed using the Aarhus stellar evolution code (ASTEC; ChristensenDalsgaard 2008a), with adiabatic pulsation frequencies calculated using the Aarhus adiabatic oscillation package (ADIPLS; ChristensenDalsgaard 2008b).
Here we used the updated set of 50 frequencies from our peakbagging to model HATP7 and two different codes and modelling schemes to asses the robustness of our results. We note that this approach does not safeguard against other potential sources of systematics that can arise, for instance, from input physics that are not covered by the two codes or from assumptions made on certain quantities in the modelling. The effect of such systematics will be studied elsewhere (Silva Aguirre et al., in prep.).
The results from the modelling can be found in Table 1. The échelle diagram (Grec et al. 1983) of HATP7 is given in Fig. 3 after correcting for the background, and overlaid are both peakbagged and modelled frequencies. In the construction of the échelle diagram we plot on the ordinate the medium frequency of the respective Δνlength segments. Here Δν denotes the socalled large separation, computed as the frequency difference between consecutive radial orders of a given degree. To obtain a better representation of the ridges for illustrative purposes a constant value was added to the frequencies before taking the modulo, thus allowing a shift of the ridges on the abscissa (see, e.g., Bedding 2011). The Δν from Huber et al. (2013b) was used.
Fig. 3
Échelle diagram for HATP7 using Δν = 59.22 μHz. The grey scale ranges from white at low power to black at high power. Circles give the extracted frequencies with corresponding uncertainties, triangles connected by red lines give model frequencies from GARSTEC, while squares connected by green lines give model frequencies from MESA. The degree of each ridge is indicated in the top part of the plot. The radial order of the l = 0 modes is indicated by the numbers at the right side of the plot. 
4.2.1. MESA model
As a first approach, HATP7 was modelled using Modules for experiments in stellar astrophysics (MESA; Paxton et al. 2011, 2013) and ADIPLS (ChristensenDalsgaard 2008b). For the modelling we used the value of the effective temperature, T_{eff} = 6350 ± 126 K, the heavyelement abundance, [Fe/H] = + 0.26 ± 0.15 (both from Huber et al. 2014), and the frequencies given in Table D.1, except for the lowest l = 0 mode, which consistently agreed only extremely poorly with the models and was subsequently excluded.
For the input physics we chose to neglect diffusion and settling following ChristensenDalsgaard et al. (2010). We used the 2005 update of the OPAL EOS (Rogers et al. 1996; Rogers & Nayfonov 2002) and the NACRE nuclear reaction rates (Angulo et al. 1999) with the updated ^{14}N(p,γ)^{15}O reaction rate by Formicola et al. (2004) and the updated ^{12}C(α,γ)^{16}O reaction rate by Kunz et al. (2002). Furthermore, we used the OPAL opacities (Iglesias & Rogers 1996) assuming the solar chemical composition given by Asplund et al. (2009), supplemented by the Ferguson et al. (2005) opacities at low temperatures. We used the mixinglength theory of convection as formulated by BöhmVitense (1958). Finally, we chose to use the “simple photosphere” option in MESA for the atmospheric boundary condition, which constitutes a grey atmosphere with the optical depth, τ_{s}, to the base of the atmosphere of 2/3 (see Paxton et al. 2011, their Eq. (3)).
In matching the model frequencies to the observed frequencies, we corrected the model frequencies for nearsurface effects (Kjeldsen et al. 2008) using the prescription by Brandão et al. (2011) with b = 4.90 and a reference frequency of ν_{0} = 1100 μHz. The bestfitting model was found via a χ^{2}minimisation (see, e.g., ChristensenDalsgaard et al. 2010; Brandão et al. 2011; Doǧan et al. 2013). The total χ^{2} was found as the weighted average of a normalised frequency and a spectroscopic component with weights given as (following Paxton et al. 2013) and the normalisation given by the number of values entering the χ^{2} component. The component χ_{spectro} included the match of model and spectroscopic values for T_{eff} and [Fe/H]. The uncertainty of a given model parameter was found as the likelihoodweighted standard deviation of total χ^{2} values for all the computed models. Since we have chosen to use the model parameters for the bestfitting model instead of the likelihoodweighted mean values, we added the difference between the mean and the bestfitting values in quadrature to the uncertainties. The parameters of the bestfitting model and their uncertainties are listed in Table 1.
4.2.2. GARSTEC model
Fig. 4
Ratios r_{02} and r_{010} as a function of frequency. The lines show the ratios obtained for the bestfit GARSTEC (black solid) and MESA (dashed blue) models. Note that the relatively small uncertainties on frequencies (horizontal error bars) renders them indiscernible on this scale. 
In the second approach we used two grids of stellar models computed with the Garching stellar evolution code (GARSTEC; Weiss & Schlattl 2008).
The input physics is similar to the description given for the MESA model, with the difference that we used the mixinglength theory of convection of Kippenhahn et al. (2013), and Grevesse & Sauval (1998) solar abundances. The updated ^{12}C(α,γ)^{16}O reaction rate was not used. In one grid we included the effects of coreovershooting using an exponential decay of the convective velocities with an efficiency of f_{ov} = 0.016 (Magic et al. 2010).
The grids cover a mass range between 0.7 and 1.8 M_{⊙} in steps of 0.01 M_{⊙} and initial compositions of −0.65 < [Fe/H] < + 0.50 in steps of 0.05 dex. These were determined using a galactic chemical evolution law of ΔY/ ΔZ = 1.4 (see, e.g., Balser 2006) anchored to the Big Bang nucleosynthesis value of Y_{p} = 0.248 (Steigman 2010). For hundreds of models along each evolutionary track, we computed theoretical oscillation frequencies using ADIPLS (ChristensenDalsgaard 2008b). These allowed us to construct dense grids of models that covered the spectroscopic and asteroseismic parameter space of our target.
To determine the stellar parameters, we used the Bayesian approach described in Silva Aguirre et al. (in prep.). Briefly, we assumed a flat prior in [Fe/H] and age including only a strict cut on the latter at 15 Gyr, and a standard Salpeter IMF (Salpeter 1955). We computed the likelihood of the observed values given a set of model parameters assuming Gaussian distributed errors. In our case, the observables included were the spectroscopic T_{eff} and [Fe/H], and the frequency ratios defined as (Roxburgh & Vorontsov 2003) Here the d_{01} and d_{10} are the smooth fivepoint small frequency separations given as We refer to Silva Aguirre et al. (2013) for further details. In Fig. 4 we show the ratios obtained from the peakbagging as a function of frequency along with ratios from the bestfit model. Here we also show the corresponding ratios obtained from the bestfit MESA model. Note that the ratios were not fitted in the MESA modelling. The construction of the ratios introduces correlations as a function of frequency that need to be taken into account when calculating the likelihood. This is done by calculating the χ^{2} entering the likelihood function as (7)where C is the covariance matrix and the vectors x (of length N) give the observed and modelled values. As an illustration of C we show in Fig. 5 the Hinton diagram for the r_{010} ratios, which provides a qualitative view of the correlation matrix for these ratios. To construct ratios we used the PPDs from the individual frequencies. From these distributions we obtained our central values and uncertainties and used the Pearson standard correlation coefficient^{5} to compute the correlation matrix.
We found that the grid that includes coreovershooting provided the best results. The final set of parameters from the modelling are given in Table 1. These were obtained from the median of the posterior probability distribution function and the 16 and 84 percent values. For the sake of comparison with the MESA model, the frequencies of the bestfit model are shown in Fig. 3 after applying the Kjeldsen et al. (2008) method to correct for nearsurface effects. Note that the surface correction is not needed for this modelling approach because frequency ratios, which are not strongly affected by the surface layers (Roxburgh & Vorontsov 2003; Silva Aguirre et al. 2011a), were used instead of the actual frequencies. Consequently, any possible apparent misfit between observed and surface corrected model frequencies in the échelle diagram cannot simply be interpreted as being caused by a poor model, but might just as well be due to a poor surface correction. In the current case, the fit, qualitatively speaking, reproduces the observations satisfactorily.
Fig. 5
Hinton diagram for the correlation matrix of r_{010} ratios. The first part of the subscript denotes whether the ratio is of type 01 or 10, while the second gives the radial order of the central frequency. White (black) squares indicate positive (negative) covariances between the ratios in question, and the size gives the relative size of the correlation (one along the diagonal). 
4.2.3. Model comparison
An improvement in our analysis compared with other works using Kepler data is that now the full available dataset can be utilised. Furthermore, the approach we adopted to extract mode frequencies, that is, using an MCMC peakbagging scheme, is expected to provide more reliable estimates and uncertainties for the frequencies than the approach of assigning frequencies from peaks in a smoothed power spectrum.
In Table 1 we compare our results with the original asteroseismic results from ChristensenDalsgaard et al. (2010). It is noteworthy that the GARSTEC model agrees quite well with these values, where only data from Q0Q1 were used. The Ftype character of the star limits the frequency precision by the large mode line widths (see Appendix B). This contributes to the relatively small changes in model results from the addition of significantly more data. Our MESA modelling results in a slightly more massive star than the stars obtained by GARSTEC and ChristensenDalsgaard et al. (2010), but still agree within uncertainties despite the different fitting techniques.
The asteroseismic modelling by Van Eylen et al. (2012) resulted in a less massive (M_{⋆} = 1.36 M_{⊙}) and less metalrich ([Fe/H] = 0.13) star than the current modelling efforts. This is likely because the authors explored a limited metallicity range and no spectroscopic constraints were included in the optimisation. Our stellar parameters agree with those derived by Pál et al. (2008) from a combined analysis of stellar isochrones, photometry, and spectroscopy. With our asteroseismic analysis we reduce the uncertainties on mass, radius, and age to 4.2%, 1.1%, and 11.2%, respectively, from the original estimates of 6.4%, 13.9%, and 45.5% given in Pál et al. (2008). In this comparison asymmetric uncertainties were added in quadrature.
Since the GARSTEC values are computed using frequency ratios that are insensitive to nearsurface effects, this set of stellar parameters was adopted for the remainder of the paper. Furthermore, this Bayesian approach offers a direct estimation of the parameter uncertainties from the posterior distributions from a dense grid of models.
Finally, we note that the values of mass, radius, and age obtained from AME (Lundkvist et al. 2014); M/M_{⊙} = 1.55 ± 0.05, R/R_{⊙} = 1.99 ± 0.02, age = 1.9 ± 0.5 Gyr, are consistent with the results from this work.
Fig. 6
Top: marginalised posterior probability distribution (PPD) for the stellar inclination i_{⋆}. Here we have folded the full distribution (−90 to 180°) onto the the range from 0 to 90°. Bottom right: PPD for the rotational splitting. Bottom left: correlation map between the inclination and the rotational splitting. The 68% credible regions (highest posterior density credible regions) are indicated by the darkgrey parts, while lightgrey indicates the additional parts covered in a 95% credible region. Included are also lines of constant v sin i_{⋆}, computed using the radius estimate from our analysis and with v sin i_{⋆} values from the literature (see Table 3), arranged in the legend in order of increasing v sin i_{⋆} value (top to bottom). The dashed black horizontal line indicates the value of ν_{s} = 5 μHz above which walkers were excluded from the PPDs. The grey dashdotted lines give lines of constant v sin i_{⋆} from 8 to 48 km s^{1} in steps of 8 km s^{1}. 
4.2.4. Updated planetary parameters
With our new estimates for the stellar parameters we can update the mass and radius for HATP7b. The mass is found from (see, e.g., Winn 2010) (8)where M_{p} is the planet mass, e is the eccentricity, P_{orb} is the orbital period, and K_{⋆} is the stellar reflex velocity. As input we use i_{p} and P_{orb} from Van Eylen et al. (2013), e is set to zero following Husnoo et al. (2012), and K_{⋆} = 212.2 ± 3.2 m s^{1} from combining estimates of K_{⋆} by Winn et al. (2009) and Narita et al. (2009). Solving for the planet mass then gives^{6}M_{p} = 1.80 ± 0.05 M_{J}. For the planetary radius we derive R_{p} = 1.51 ± 0.02 R_{J} when using R_{p}/R_{⋆} from Van Eylen et al. (2013). Considering the discussion in Van Eylen et al. (2013) on the size of systematic effects on transit depth measurements, we adopted a one per cent uncertainty on R_{p}/R_{⋆}. This gives a planetary density of ρ_{p} = 0.65 ± 0.03 g cm^{3}, which is just within the uncertainty of the estimate from Pál et al. (2008) of . We note that the stellar parameters estimated from asteroseismology offers a greatly improved precision on the planetary parameters.
4.3. Splitting and inclination
In Fig. 6 we present the results for the splitting and inclination of HATP7 as the PPDs for these parameters, along with the 2D correlation map. These results were obtained from the small fit described above. The 68% and 95% HPD credible regions are indicated in Fig. 6.
The distributions for inclination and splitting for HATP7 are unfortunately not simple and Gaussian, but rather cover an extended region of parameter space. Therefore we refrain from using measures of central tendency, such as the median or mode of the distributions, but instead report values from the 68% HPD credible regions. The overall behaviour seen in most of the correlation map corresponds well to what might be expected given the relatively large line widths: Γ = 4.8 − 8.4 μHz for the range fitted (see Fig. B.1). At low inclinations, the central m = 0 component dominates the relative heights of azimuthal components of a rotationally split multiplet, while at high inclinations the sectoral (m = ± l) components dominate. For l = 2 the tesseral (0 <  m  <l) components dominate at an inclination of ~ 50°. Thereby, the splitting is generally less constrained at low inclinations because the m ≠ 0 components are small and harder to fit. Furthermore, a line width exceeding the splitting effectively hides these small components up to a splitting of half the line width (i.e., up to ~ 4 μHz). At higher inclinations the splitting is more easily discernible because the visible azimuthal components are the ones farthest apart, thus the line width becomes less of a problem.
There is a high density at i_{⋆} ~ 10° and ν_{s}> 5 μHz. This should not be interpreted as a good solution to the fit of the power spectrum but is rather the result of walkers trying to escape the parameter space, with the consequence that the walkers pile up at the boundary of the prior on the splitting. A contributing factor could also be that the small separation between l = 0 and l = 2 modes is of the order ~ 5 μHz, whereby l = 2 azimuthal components would be within the l = 0 mode profile. We see no correlations between these high splitting values and any of the other parameters of the fit and so these parameters are unaffected by this feature. These walkers will, however, contribute to the shape of the splitting and inclination PPDs, and therefore we chose to exclude walkers with ν_{s}> 5 μHz in the marginalisation of the i_{⋆} and ν_{s} PPDs. This is the same effect as observed when the background was kept free in the fit, but now much decreased. It is, however, not clear to us what about the inclination of i_{⋆} ~ 10° facilitates this escape of walkers.
From the 68% HPD credible regions of the inclination and splitting PPDs we find i_{⋆}< 36.5° and ν_{s}< 0.87 μHz (see Table 2). The value for the inclination agrees overall with previous reports that were based on the low value for v sin i_{⋆}. In addition, agreement is found with the statistically derived estimate by Schlaufman (2010) of i_{⋆} = 6.7 − 10° based on a simple P_{rot} ∝ t^{1/2} model (Weber & Davis 1967; Skumanich 1972) for the evolution of the stellar rotation period as a function of mass and age. The splitting, on the other hand, is lower than expected for an F6type star (see below). However, it is clear from the correlation map that at low inclinations, values of the splitting up to ~ 4 μHz are allowed.
Values related to the stellar inclination and rotation.
Literature values for parameters related to the HATP7 system.
Fig. 7
Distribution for v sin i_{⋆} constructed from the PPDs obtained for ν_{s} and i_{⋆} and assuming (see Eq. (9)). The dark and lightgrey regions of the distribution correspond to the 68% and 95% credible regions. Literature values for v sin i_{⋆} obtained from spectral and RM analysis (see Table 3) are given by vertical lines where the shaded regions above the distribution give the corresponding uncertainties. 
With our estimate for the stellar radius (see Sect. 4.2) we may convert the estimated ranges for the splitting and inclination to a measure of v sin i_{⋆} as (Chaplin et al. 2013) (9)Using the PPDs obtained for ν_{s} and i_{⋆}, while assuming a distribution for the radius as , we find value of v sin i_{⋆}< 2.21 km s^{1} (68% HPD credible region). In Fig. 7 we have plotted the full distribution for v sin i_{⋆} and indicate the values obtained from spectral and RM analysis (see Table 3). We find that the value of Narita et al. (2009) agrees within uncertainties with our results. Winn et al. (2009) found the highest value for v sin i_{⋆}. Their higher v sin i_{⋆} value results from the lower orbital inclination value (see Table 3) they used to model the RM effect. For HATP7 a larger impact parameter (lower inclination) requires a larger v sin i_{⋆}. If the analysis were repeated with the inclination derived from the Kepler light curve instead of the groundbased data that were available to them, this disagreement in v sin i_{⋆} would vanish.
After measuring v sin i_{⋆} and i_{⋆} and assuming solidbody rotation, we can now calculate the true rotation speed of HATP7, v< 7.66 km s^{1} (68% HPD credible region). This upperlimit estimate agrees with the values from Nielsen et al. (2013) for the rotation periods of Kepler stars of an approximate spectral type between F4 and 6.
5. Obliquity
With our estimate for the stellar inclination, i_{⋆}, λ, and the planetary orbital inclination, i_{p}, we are now able to calculate the system obliquity, ψ, from (Winn et al. 2005) (10)In Fig. 8 we show the distribution for ψ for the obtained distribution for i_{⋆} (see Fig. 6), while adopting i_{p} from Van Eylen et al. (2013), and λ from Albrecht et al. (2012). From the 68% HPD credible region we obtain 83°<ψ< 111°, consistent with a polar orbit. The corresponding results using λ from Winn et al. (2009) and Narita et al. (2009) are given in Table 3.
We refer to Albrecht et al. (2012) for a discussion on the different values for λ and the possible reasons for their disagreement (see also Albrecht et al. 2011). Here we note that regardless of which λvalue we used we found a polar orbit for HATP7b.
When the obliquity is assessed without knowing the stellar inclination, a flat distribution in cos i_{⋆} is generally assumed for the stellar orientation, because this results in an isotropic distribution for the stellar inclination. From this distribution it is a priori much more likely to observe a random star in an equatoron configuration. In Fig. 8 we show the distribution in ψ from adopting this isotropic distribution in i_{⋆}. In this way, Winn et al. (2009) estimated ψ> 86.3° with 99.73% confidence, while Narita et al. (2009) found ψ> 90° with 99.70% confidence (both used i_{p} from Pál et al. 2008). From this approach a retrograde orbit is thus strongly suggested, but the orbit has a higher probability of being more equatorial than polar. With our asteroseismic estimate for the inclination we can substantiate these statistical results because we find ψ> 90° with 68% credibility using λ from Albrecht et al. (2012), and now a nearpolar orbit is the most likely configuration for the system.
Fig. 8
Distribution of the true angle ψ using the distribution for i_{⋆} from the peakbagging, and with the assumption of normal distributions for the planetary inclination i_{p} and projected angle λ. For i_{p} the value from Van Eylen et al. (2013) was adopted, while λ was taken from Albrecht et al. (2012). The dark and lightgrey regions of the distribution correspond to the 68% and 95% credible regions. The dashed curve gives the distribution assuming instead an isotropic distribution for i_{⋆}, i.e., flat in cos i_{⋆}. We indicate also which values of ψ correspond to a retrograde or prograde orbit of HATP7b. 
6. Comparison with gyrochronology
Our result for the limits on the stellar rotation rate (see Table 2) can be compared with empirically calibrated gyrochronology relations. Here we use the form described by Barnes (2007) given as (11)where t is the stellar age in Myr, while a, b, and n are empirically determined coefficients that vary depending on the calibration set used (see, e.g., Epstein & Pinsonneault 2014). To compare with this relation, we first need an estimate for the dereddened colour (B − V)_{0}. The procedure used to derive (B − V)_{0} is described in Appendix C.
Using the estimate (B − V)_{0} = 0.495 ± 0.022 together with the age determined in Sect. 4.2 of t = 2.07 ± 0.36 Gyr (see Table 1; asymmetric uncertainties were added in quadrature), we can estimate the rotation period from the relation in Eq. (11). In Fig. 9 we show two versions of this relation, those by Barnes (2007) (B07), and those by Meibom et al. (2009; M09; see also Barnes & Kim 2010; Meibom et al. 2011). From these we compute periods of 9.9 ± 1.8 (B07) and 5.0 ± 4.0 days (M09). The upper value from the B07 relation is similar to our lowerlimit period estimate of about 13 days (see Table 2), which means that the agreement is not very convincing. On the other hand, the splittings at low inclinations match this range of rotation periods well because the level of ~ 4 μHz corresponds to a rotation period of ~ 2.9 days (see Fig. 6).
We note that a P_{rot} ∝ t^{1/2} law might provide a poor description of the rotational evolution for HATP7, because its (B − V)_{0} places it in close proximity to the socalled Kraft break (Kraft 1967) where loss of angular momentum via a stellar wind is inhibited by the lack of a sufficiently deep convection zone. If HATP7 is on the low (B − V)_{0} side of this break, the rotation rate becomes a strong function of the initial conditions (see, e.g., van Saders & Pinsonneault 2013), and a gyrochronology scaling is not applicable.
The closein hotJupiter HATP7b has potentially had an impact on the rotation rate of its host star. While a detailed dynamical analysis of this system is beyond the scope of this paper, the synchronisation of rotation and alignment clearly has not been reached yet. The time scale for circularisation has been reached, however, with an upper limit on the eccentricity of e< 0.038 given by Husnoo et al. (2012). While it is difficult to assess the past interaction between the planet and the star, we note that with the configuration we estimate for the system, that is, close to a polar orbit of the planet, the rotational synchronisation time will likely be very long because the angular momentum vectors are close to being orthogonal, thereby decreasing the tidal interaction. For discussions on the interaction between hot Jupiters and their host stars we refer to Hut (1981), Miller et al. (2009), Cohen et al. (2010), Matsumura et al. (2010), Mardling (2011), and Valsecchi & Rasio (2014), for example.
Fig. 9
Gyrochronology relations from Meibom et al. (2009; M09) and Barnes (2007; B07). The shaded region around each relation (with the same colour) represents the standard error of the relationship between (B − V)_{0} and rotation period from propagating the uncertainties reported for the coefficients entering Eq. (11) together with the uncertainty in our age estimate. The vertical line and shaded region give the value of (B − V)_{0} = 0.495 ± 0.022 of HATP7, while darker shaded regions around the relations vertically bound the corresponding uncertainty in the period from the uncertainty in (B − V)_{0}, age, and the gyrochronology relations. The horizontal lines indicate the limits on the 68% (dashed) and 95% (dotted) credible regions on the stellar rotation period from pmode splittings. The upwardpointing arrows indicate that these rotation periods are lower limits from the respective HPD credible regions. 
7. Activity signatures
Activity of stars is linked to the interaction between rotation, convection, and magnetic fields. The interaction can cause magnetic features, such as dark spots and bright faculae, to appear in the photosphere of the stars (see, e.g., Berdyugina 2005, for a review), and plages in the chromosphere. If a star is rotating, such dark or bright regions will induce a temporal modulation of the integrated stellar flux.
The photospheric features can also cause departures from radiative equilibrium in the stellar chromosphere in plages, and induce emission in the cores of specific spectral lines. For example, emission in cores of Ca ii H&K lines is an often used indicator of magnetic activity (Knaack et al. 2001; Wright et al. 2004; Frasca et al. 2011; Fröhlich et al. 2012) because it is believed to reflect the amount of nonthermal chromospheric heating above faculae (see, e.g., Hall 2008, for a review). For closein hotJupiter systems there is furthermore the possibility for a direct magnetic interaction between the planet and the star, where magnetic reconnections, similar to those seen in flares, can cause a heating of the chromosphere (see, e.g., Shkolnik et al. 2005).
7.1. Chromospheric activity
To test for signatures of chromospheric activity we searched the Ca ii H&K lines in the HIRES spectra from Winn et al. (2009), but found no signs of emission or high levels of activity. Because we see pmode oscillations in HATP7, it must have an outer convection zone, and the (B − V)_{0} colour derived in Sect. 6 places it above the limit from Simon & Landsman (1991) for the onset of activity. The absence of an emission signal could point towards an intrinsically low surface activity caused by a low rotation rate.
Another effect that would influence the signal is if HATP7 has a low angle of inclination. On the Sun, faculae are primarily located in the active latitude bands between around 5 and 40° latitude. As faculae and plages are believed to have the same magnetic driver the Ca ii H&K emission from plages will thus decrease with decreasing inclination from the decreasing projected area of the active regions (Knaack et al. 2001; Freire Ferrero et al. 2004; Chaplin et al. 2007).
Fig. 10
Top: autocorrelation function (ACF) of the HATP7 time series (black). The dashed horizontal line gives the expectation value for random independent and identically distributed values, while the grey curves give the largelag 95% confidence levels. Bottom: lowfrequency end of the power spectrum in units of period. 
Fig. 11
Morlet wavelet power spectrum as a function of time for HATP7, where the time series was binned by 31 points. The colour changes from red at low power to blue at high power – the colours are on a logarithmic scale. The crosshatched regions at high and low times indicate the cone of influence where edge effects become important (see Torrence & Compo 1998). No clear signatures from rotation are apparent. 
7.2. Lowfrequency region
We searched for a signal imparted by a temporal flux modulation at low frequencies in the power spectrum of HATP7 (see, e.g., Campante et al. 2011; Nielsen et al. 2013). Examples of such activity signals in Ftype Kepler stars can be seen, for instance, in Mathur et al. (2014).
The signal from the stellar rotation can also be estimated for example by the autocorrelation function (ACF) of the time series, as shown for instance by McQuillan et al. (2013). In Fig. 10 we show the ACF of the corrected time series together with the lowfrequency end of the power spectrum (in units of period for convenience). The ACF shows no clear sign of modulation. To test the hypothesis that a signal in the time series has died out at lag k, we used the largelag standard error (Anderson 1976), indicated in Fig. 10. From this it is clear that the signal seen in the ACF is not significant. However, the lowamplitude hump at around ~ 9 days in the ACF does seem to align with an increase in power in the power spectrum.
We investigated the presence of a modulation also via the Morlet wavelet transform (Torrence & Compo 1998) of the time series (see also Mathur et al. 2013), this is shown in Fig. 11. No clear periodicity is seen, and the signal that causes the ~ 9 day hump in the ACF and the power spectrum seems to be very intermittent. Finally, we checked the magnetic proxy from Campante et al. (2014), but we did not find indications of a signal here either. The fact that a stronger or more localised signal is missing for HATP7 might again be linked to a low stellar inclination. Indeed, if magnetic features such as stellar spots primarily reside near the equator of the star, while we view it close to poleon, a strong modulation in the light curve would not be expected.
We also note a collection of peaks in the power spectrum around a frequency of 6.7 ± 0.4 μHz (corresponding to a period of ~ 41.6 h). The presence of this power excess is quite robust and is consistently seen in random segments of the total time series and in the wavelet spectrum. We therefore rule out that it originates from random noise. First, we checked that there are no known artefacts at this frequency (Christiansen et al. 2013). To determine whether the signal might originate from another star in the vicinity of HATP7 (K_{p} = 10.463) we located all Kepler targets with available data within a radius of 3′ (12 were found). Of these, the star KIC 10666727 (K_{p} = 13.166) was found to have strong signatures of spot modulations in its time series and broadened power excess peaks in its power spectrum around ~ 2.2 μHz, ~ 4.6 μHz, and most interestingly, in the region 6.7 − 7.2 μHz. The angular separation between this target and HATP7 is around 155″, which places the two stars close enough for direct PRF (pixel response function) contamination to occur between them (Coughlin et al. 2014). However, according to Coughlin et al. (2014), stars have to be bright to contaminate over such a large separation. We do not see any significant power excess in HATP7 at other frequencies where KIC 10666727 shows an even stronger excess than at 6.7 ± 0.4 μHz. Therefore we find it possible, but unlikely, that KIC 10666727 is contaminating the light curve of HATP7.
8. Discussion
With asteroseismic modelling we have provided a precise and detailed model for HATP7, with parameters that can be used in tests of theories on the formation and evolution of planetary systems. These, together with the obliquity, are especially important to constrain when attempting to explain a system such as HATP7, which likely has a closeto polar orbit of its hotJupiter planet.
Our estimate of the obliquity of the HATP7 system supports the hypothesis described by Winn et al. (2010), who stated that hot Jupiters are born with a wide range of obliquities. Planetplanet scatterings (Chatterjee et al. 2008) and the effect of Kozai cycles and tidal friction (Fabrycky & Tremaine 2007) could for example create large initial obliquities. For cool dwarfs with deep convection zones, tidal dissipation would operate efficiently and result in aligned systems. Hot stars fail to align because they lack such a deep convection zone for most of their life on the mainsequence, if not altogether. This hypothesis was put forth based on the observed trend of λ against T_{eff}, where a broad distribution is seen in λ for T_{eff}> 6250 K, while predominantly low values for λ are seen below this temperature. This result was corroborated by Albrecht et al. (2012) based on a larger sample of measurements. These authors also found that the systems with respect to λ could be sorted into relative tidal timescales, which in turn depend on the planet to star mass ratio, and a/R_{⋆} (a being the semimajor axis), among others. While HATP7 falls on the hot side of the dividing temperature between the two regimes and tidal interactions are expected to be weak, it is worth noting that the temperature is quite close to this dividing line. Moreover, a/R_{⋆} is small, which increases tidal interaction (Albrecht et al. 2012). The relatively slow rotation we found suggests that some magnetic braking has taken place and is still ongoing. The lack of alignment is likely linked to the combined effects of a low mass of the convection zone (making tides ineffective), the closeto orthogonal alignment of the angular momentum vectors, and the relatively young age of the system. Indeed, the star might have started out with a more rapid rotation, where magnetic braking from the developing convection zone has been more effective in slowing down the star than has the tidal interaction in realigning the system.
With regard to the obliquity of the system, the fact that a third, and possibly a fourth, body is found to be associated (Winn et al. 2009; Narita et al. 2012) might support the scenario of fewbody dynamical interaction early in the life of the system, which resulted in the high obliquity of HATP7b.
9. Conclusions
Using asteroseismology, we have estimated the stellar inclination and provided new and precise stellar parameters for HATP7 based on the extracted mode frequencies. Mode frequencies were extracted using a Bayesian MCMC approach to peakbag the frequency power spectrum from corrected Kepler data, for which we utilised the full SC dataset from Q0Q17 available from the Kepler satellite. From this, information was obtained about stellar, planetary, and system parameters (age, mass, radius, composition, luminosity, mass of convective core) that are important ingredients, for instance, in dynamical simulations used to test theories for the evolution of planetary systems.
We found the star of the HATP7 system to have a low inclination, i_{⋆}< 36.5° with 68% credibility, meaning that the star is seen close to poleon. Combining this with estimates for the planetary orbital inclination and the projected obliquity from RM measurements, the closein hotJupiter planet is in a high obliquity and likely retrograde orbit (the retrograde solutions account for ~ 68% of the PDF in Fig. 8). Our estimate for the stellar rotation matches empirical findings for stars of the same spectral type, and using the age of the system from our modelling of the oscillation frequencies in combination with an improved estimate for the colour of the star (see Appendix C) yields that estimates from gyrochronology are not in conflict with our results. While the lack of signatures from activity are by no means proof of a low inclination for the star, they are not incompatible with a low inclination.
To our knowledge, this analysis is the first wherein asteroseismology has been able to provide an estimate for i_{⋆} that, together with λ from RM measurements and i_{p} from analysis of the transit profile, has allowed for a near complete description of the system geometry.
For theories attempting to explain the formation and evolution of planetary systems, the HATP7 system is highly interesting because any such theory must be able to offer an explanation for the system geometry. The analysis presented in this paper has shed more light on the obliquity of the system, and not just the projected obliquity that is normally used in obliquity studies. The result on the obliquity agrees with assumptions based on the high value of λ and the low value for v sin i_{⋆}. Furthermore, the high value for ψ corroborates the theory where the degree of alignment is connected to the tidal evolution of the system (see, e.g., Winn et al. 2010; Albrecht et al. 2012; Valsecchi & Rasio 2014).
An aspect of the HATP7 system that might be investigated to constrain the obliquity even better is the apparent asymmetry of the transit light curve, seen for instance in the phase curve presented in Esteves et al. (2013) and Van Eylen et al. (2013). Such an asymmetry is also seen in the KOI13 system, for example, and it was found by Szabó et al. (2011) and Barnes et al. (2011) to agree well with the predictions by Barnes (2009) for a planet crossing over the gravitybrightened polar region of its rapidly rotating host star on a highobliquity orbit.
KeckI/HIRES: High Resolution Spectrograph (Vogt et al. 1994); Subaru/HDS: High Dispersion Spectrograph (Noguchi et al. 2002).
Using the geweke module of PyMC (Patil et al. 2010).
Acknowledgments
The authors wish to thank the entire Kepler team, without whom these results would not be possible. We are thankful to O. Benomar and his collaborators for their approach to the contemporaneousness of our respective studies. We would like to thank Eric Stempels and Heidi Korhonen for useful discussions in the early phases of the project. Thanks to Joshua N. Winn and John A. Johnson for making their Keck spectra available to us. Funding for the Stellar Astrophysics Centre (SAC) is provided by The Danish National Research Foundation. The research is supported by the ASTERISK project (ASTERoseismic Investigations with SONG and Kepler) funded by the European Research Council (Grant agreement No. 267864). M.L. would like to thank the asteroseismology group at SIfA for their hospitality during a research stay where some of this work was carried out. M.L. also wishes to thank Niels Bohr Fondet for financial support for the research stay at SIfA. T.L.C., G.R.D., W.J.C., and R.H. acknowledge the support of the UK Science and Technology Facilities Council (STFC). C.K. acknowledges the support of the Villum Foundation. M.B.N. acknowledges research funding by Deutsche Forschungsgemeinschaft (DFG) under grant SFB 963/1 “Astrophysical flow instabilities and turbulence” (Project A18). The research leading to these results has received funding from the European Community’s Seventh Framework Programme (FP7/20072013) under grant agreement No. 269194. This research has made use of the following web resources: the SIMBAD database (simbad.ustrasbg.fr), operated at CDS, Strasbourg, France; NASAs Astrophysics Data System Bibliographic Services (adswww.harvard.edu); arxiv.org, maintained and operated by the Cornell University Library.
References
 Ahlers, J. P., Seubert, S. A., & Barnes, J. W. 2014, ApJ, 786, 131 [NASA ADS] [CrossRef] [Google Scholar]
 Albrecht, S., Winn, J. N., Johnson, J. A., et al. 2011, ApJ, 738, 50 [NASA ADS] [CrossRef] [Google Scholar]
 Albrecht, S., Winn, J. N., Johnson, J. A., et al. 2012, ApJ, 757, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, E. R., Duvall, Jr., T. L., & Jefferies, S. M. 1990, ApJ, 364, 699 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, O. 1976, Time series analysis and forecasting: the BoxJenkins approach (Butterworth) [Google Scholar]
 Angulo, C., Arnould, M., Rayet, M., et al. 1999, Nucl. Phys. A, 656, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Appourchaux, T., Benomar, O., Gruberbauer, M., et al. 2012, A&A, 537, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Ballot, J., Barban, C., & van’t VeerMenneret, C. 2011, A&A, 531, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Balser, D. S. 2006, AJ, 132, 2326 [NASA ADS] [CrossRef] [Google Scholar]
 Barnes, J. W. 2009, ApJ, 705, 683 [Google Scholar]
 Barnes, J. W., Linscott, E., & Shporer, A. 2011, ApJS, 197, 10 [NASA ADS] [CrossRef] [Google Scholar]
 Barnes, S. A. 2007, ApJ, 669, 1167 [NASA ADS] [CrossRef] [Google Scholar]
 Barnes, S. A., & Kim, Y.C. 2010, ApJ, 721, 675 [NASA ADS] [CrossRef] [Google Scholar]
 Bedding, T. R. 2011 [arXiv:1107.1723] [Google Scholar]
 Benomar, O., Appourchaux, T., & Baudin, F. 2009, A&A, 506, 15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Berdyugina, S. V. 2005, Liv. Rev. Sol. Phys., 2 [Google Scholar]
 Bergfors, C., Brandner, W., Daemgen, S., et al. 2013, MNRAS, 428, 182 [NASA ADS] [CrossRef] [Google Scholar]
 BöhmVitense, E. 1958, ZAp, 46, 108 [NASA ADS] [Google Scholar]
 Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Brandão, I. M., Doğan, G., ChristensenDalsgaard, J., et al. 2011, A&A, 527, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Campante, T. L., Handberg, R., Mathur, S., et al. 2011, A&A, 534, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Campante, T. L., Chaplin, W. J., Lund, M. N., et al. 2014, ApJ, 783, 123 [NASA ADS] [CrossRef] [Google Scholar]
 Casagrande, L., & VandenBerg, D. A. 2014, MNRAS, submitted [arXiv:1407.6095] [Google Scholar]
 Casagrande, L., Ramírez, I., Meléndez, J., Bessell, M., & Asplund, M. 2010, A&A, 512, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Casagrande, L., Ramírez, I., Meléndez, J., & Asplund, M. 2012, ApJ, 761, 16 [NASA ADS] [CrossRef] [Google Scholar]
 Casagrande, L., Silva Aguirre, V., Stello, D., et al. 2014, ApJ, 787, 110 [NASA ADS] [CrossRef] [Google Scholar]
 Chaplin, W. J., & Miglio, A. 2013, ARA&A, 51, 353 [NASA ADS] [CrossRef] [Google Scholar]
 Chaplin, W. J., Elsworth, Y., Houdek, G., & New, R. 2007, MNRAS, 377, 17 [Google Scholar]
 Chaplin, W. J., SanchisOjeda, R., Campante, T. L., et al. 2013, ApJ, 766, 101 [NASA ADS] [CrossRef] [Google Scholar]
 Chatterjee, S., Ford, E. B., Matsumura, S., & Rasio, F. A. 2008, ApJ, 686, 580 [NASA ADS] [CrossRef] [Google Scholar]
 ChristensenDalsgaard, J. 2008a, Ap&SS, 316, 13 [NASA ADS] [CrossRef] [Google Scholar]
 ChristensenDalsgaard, J. 2008b, Ap&SS, 316, 113 [NASA ADS] [CrossRef] [Google Scholar]
 ChristensenDalsgaard, J., Kjeldsen, H., Brown, T. M., et al. 2010, ApJ, 713, L164 [NASA ADS] [CrossRef] [Google Scholar]
 Christiansen, J. L., Jenkins, J. M., Caldwell, D. A., et al. 2013, Kepler Data Characteristics Handbook, (KSCI  19040  004), http://archive.stsci.edu/kepler/documents.html [Google Scholar]
 Cohen, O., Drake, J. J., Kashyap, V. L., Sokolov, I. V., & Gombosi, T. I. 2010, ApJ, 723, L64 [NASA ADS] [CrossRef] [Google Scholar]
 Coughlin, J. L., Thompson, S. E., Bryson, S. T., et al. 2014, AJ, 147, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Davies, G. R., Chaplin, W. J., Farr, W. M., et al. 2014, MNRAS, submitted [Google Scholar]
 Désert, J.M., Charbonneau, D., Demory, B.O., et al. 2011, ApJS, 197, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Doǧan, G., Metcalfe, T. S., Deheuvels, S., et al. 2013, ApJ, 763, 49 [NASA ADS] [CrossRef] [Google Scholar]
 Dziembowski, W. 1977, Acta Astron., 27, 203 [NASA ADS] [Google Scholar]
 Epstein, C. R., & Pinsonneault, M. H. 2014, ApJ, 780, 159 [NASA ADS] [CrossRef] [Google Scholar]
 Esteves, L. J., De Mooij, E. J. W., & Jayawardhana, R. 2013, ApJ, 772, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Fabrycky, D., & Tremaine, S. 2007, ApJ, 669, 1298 [NASA ADS] [CrossRef] [Google Scholar]
 Fabrycky, D. C., & Winn, J. N. 2009, ApJ, 696, 1230 [NASA ADS] [CrossRef] [Google Scholar]
 Faedi, F., Staley, T., Gómez Maqueo Chew, Y., et al. 2013, MNRAS, 433, 2097 [NASA ADS] [CrossRef] [Google Scholar]
 Ferguson, J. W., Alexander, D. R., Allard, F., et al. 2005, ApJ, 623, 585 [NASA ADS] [CrossRef] [Google Scholar]
 Ford, E. B., & Rasio, F. A. 2008, ApJ, 686, 621 [NASA ADS] [CrossRef] [Google Scholar]
 ForemanMackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [CrossRef] [Google Scholar]
 Formicola, A., Imbriani, G., Costantini, H., et al. 2004, Phys. Lett. B, 591, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Frandsen, S., Jones, A., Kjeldsen, H., et al. 1995, A&A, 301, 123 [NASA ADS] [Google Scholar]
 Frasca, A., Fröhlich, H.E., Bonanno, A., et al. 2011, A&A, 532, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Freire Ferrero, R., Frasca, A., Marilli, E., & Catalano, S. 2004, A&A, 413, 657 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fröhlich, H.E., Frasca, A., Catanzaro, G., et al. 2012, A&A, 543, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gandolfi, D., Collier Cameron, A., Endl, M., et al. 2012, A&A, 543, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Geweke, J. 1992, in Bayesian statistics, 4 (Oxford Univ. Press, New York), 169 [Google Scholar]
 Gilliland, R. L., Brown, T. M., ChristensenDalsgaard, J., et al. 2010, PASP, 122, 131 [NASA ADS] [CrossRef] [Google Scholar]
 Gizon, L., & Solanki, S. K. 2003, ApJ, 589, 1009 [NASA ADS] [CrossRef] [Google Scholar]
 Goldreich, P., Murray, N., & Kumar, P. 1994, ApJ, 424, 466 [NASA ADS] [CrossRef] [Google Scholar]
 Gough, D. O. 1977a, ApJ, 214, 196 [NASA ADS] [CrossRef] [Google Scholar]
 Gough, D. O. 1977b, in Problems of Stellar Convection, eds. E. A. Spiegel, & J.P. Zahn (Berlin: Springer Verlag), Lect. Notes Phys., 71, 57 [Google Scholar]
 Grec, G., Fossat, E., & Pomerantz, M. A. 1983, Sol. Phys., 82, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Grevesse, N., & Sauval, A. J. 1998, Space Sci. Rev., 85, 161 [NASA ADS] [CrossRef] [Google Scholar]
 Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hall, J. C. 2008, Liv. Rev. Sol. Phys., 5, 2 [Google Scholar]
 Handberg, R., & Campante, T. L. 2011, A&A, 527, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Handberg, R., & Lund, M. N. 2014, MNRAS, accepted [arXiv:1409.1366] [Google Scholar]
 Harvey, J. 1985, in Future Missions in Solar, Heliospheric & Space Plasma Physics, eds. E. Rolfe, & B. Battrick, ESA SP, 235, 199 [Google Scholar]
 Henden, A. A., Welch, D. L., Terrell, D., & Levine, S. E. 2009, in AAS Meeting Abstracts, 214, 407.02 [Google Scholar]
 Hirano, T., SanchisOjeda, R., Takeda, Y., et al. 2012, ApJ, 756, 66 [NASA ADS] [CrossRef] [Google Scholar]
 Hirano, T., SanchisOjeda, R., Takeda, Y., et al. 2014, ApJ, 783, 9 [NASA ADS] [CrossRef] [Google Scholar]
 Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27 [NASA ADS] [Google Scholar]
 Hou, F., Goodman, J., Hogg, D. W., Weare, J., & Schwab, C. 2012, ApJ, 745, 198 [NASA ADS] [CrossRef] [Google Scholar]
 Houdek, G. 2006, in Proc. SOHO 18/GONG 2006/HELAS I, Beyond the spherical Sun, ESA SP, 624, [Google Scholar]
 Houdek, G. 2012, in Progress in Solar/Stellar Physics with Helio and Asteroseismology, eds. H. Shibahashi, M. Takata, & A. E. LynasGray, ASP Conf. Ser., 462, 7 [Google Scholar]
 Houdek, G., & Gough, D. O. 2002, MNRAS, 336, L65 [NASA ADS] [CrossRef] [Google Scholar]
 Houdek, G., Balmforth, N. J., ChristensenDalsgaard, J., & Gough, D. O. 1999, A&A, 351, 582 [NASA ADS] [Google Scholar]
 Huber, D., Carter, J. A., Barbieri, M., et al. 2013a, Science, 342, 331 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Huber, D., Chaplin, W. J., ChristensenDalsgaard, J., et al. 2013b, ApJ, 767, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Huber, D., Silva Aguirre, V., Matthews, J. M., et al. 2014, ApJS, 211, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Husnoo, N., Pont, F., Mazeh, T., et al. 2012, MNRAS, 422, 3151 [NASA ADS] [CrossRef] [Google Scholar]
 Hut, P. 1981, A&A, 99, 126 [NASA ADS] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
 Karoff, C. 2008, Ph.D. Thesis, University of Aarhus, Denmark [Google Scholar]
 Kippenhahn, R., Weigert, A., & Weiss, A. 2013, Stellar Structure and Evolution [Google Scholar]
 Kjeldsen, H. 1992, Ph.D. Thesis, University of Aarhus, Denmark [Google Scholar]
 Kjeldsen, H., & Frandsen, S. 1992, PASP, 104, 413 [NASA ADS] [CrossRef] [Google Scholar]
 Kjeldsen, H., Bedding, T. R., & ChristensenDalsgaard, J. 2008, ApJ, 683, L175 [NASA ADS] [CrossRef] [Google Scholar]
 Knaack, R., Fligge, M., Solanki, S. K., & Unruh, Y. C. 2001, A&A, 376, 1080 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Koch, D. G., Borucki, W. J., Basri, G., et al. 2010, ApJ, 713, L79 [NASA ADS] [CrossRef] [Google Scholar]
 Kraft, R. P. 1967, ApJ, 150, 551 [NASA ADS] [CrossRef] [Google Scholar]
 Kunz, R., Fey, M., Jaeger, M., et al. 2002, ApJ, 567, 643 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Lai, D., Foucart, F., & Lin, D. N. C. 2011, MNRAS, 412, 2790 [NASA ADS] [CrossRef] [Google Scholar]
 Ledoux, P. 1951, ApJ, 114, 373 [NASA ADS] [CrossRef] [Google Scholar]
 Lund, M. N., Kjeldsen, H., ChristensenDalsgaard, J., Handberg, R., & Silva Aguirre, V. 2014, ApJ, 782, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Lundkvist, M., Kjeldsen, H., & Silva Aguirre, V. 2014, A&A, 566, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Magic, Z., Serenelli, A., Weiss, A., & Chaboyer, B. 2010, ApJ, 718, 1378 [NASA ADS] [CrossRef] [Google Scholar]
 Mardling, R. A. 2011, in IAU Symp. 276, eds. A. Sozzetti, M. G. Lattanzi, & A. P. Boss, 238 [Google Scholar]
 Mathur, S., García, R. A., Morgenthaler, A., et al. 2013, A&A, 550, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathur, S., García, R. A., Ballot, J., et al. 2014, A&A, 562, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Matsumura, S., Peale, S. J., & Rasio, F. A. 2010, ApJ, 725, 1995 [NASA ADS] [CrossRef] [Google Scholar]
 McLaughlin, D. B. 1924, ApJ, 60, 22 [NASA ADS] [CrossRef] [Google Scholar]
 McQuillan, A., Aigrain, S., & Mazeh, T. 2013, MNRAS, 432, 1203 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Meibom, S., Mathieu, R. D., & Stassun, K. G. 2009, ApJ, 695, 679 [NASA ADS] [CrossRef] [Google Scholar]
 Meibom, S., Mathieu, R. D., Stassun, K. G., Liebesny, P., & Saar, S. H. 2011, ApJ, 733, 115 [NASA ADS] [CrossRef] [Google Scholar]
 Miller, N., Fortney, J. J., & Jackson, B. 2009, ApJ, 702, 1413 [NASA ADS] [CrossRef] [Google Scholar]
 Morris, B. M., Mandell, A. M., & Deming, D. 2013, ApJ, 764, L22 [NASA ADS] [CrossRef] [Google Scholar]
 Morton, T. D., & Johnson, J. A. 2011, ApJ, 729, 138 [NASA ADS] [CrossRef] [Google Scholar]
 Nagasawa, M., Ida, S., & Bessho, T. 2008, ApJ, 678, 498 [NASA ADS] [CrossRef] [Google Scholar]
 Narita, N., Sato, B., Hirano, T., & Tamura, M. 2009, PASJ, 61, L35 [NASA ADS] [Google Scholar]
 Narita, N., Kudo, T., Bergfors, C., et al. 2010, PASJ, 62, 779 [NASA ADS] [Google Scholar]
 Narita, N., Takahashi, Y. H., Kuzuhara, M., et al. 2012, PASJ, 64, L7 [NASA ADS] [Google Scholar]
 Nielsen, M. B., Gizon, L., Schunker, H., & Karoff, C. 2013, A&A, 557, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Noguchi, K., Aoki, W., Kawanomoto, S., et al. 2002, PASJ, 54, 855 [NASA ADS] [CrossRef] [Google Scholar]
 Pál, A., Bakos, G. Á., Torres, G., et al. 2008, ApJ, 680, 1450 [NASA ADS] [CrossRef] [Google Scholar]
 Patil, A., Huard, D., & Fonnesbeck, C. J. 2010, Journal of Statistical Software, 35, 1 [Google Scholar]
 Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 [Google Scholar]
 Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Ramírez, I., Michel, R., Sefako, R., et al. 2012, ApJ, 752, 5 [NASA ADS] [CrossRef] [Google Scholar]
 Rogers, F. J., & Nayfonov, A. 2002, ApJ, 576, 1064 [Google Scholar]
 Rogers, F. J., Swenson, F. J., & Iglesias, C. A. 1996, ApJ, 456, 902 [NASA ADS] [CrossRef] [Google Scholar]
 Rogers, T. M., Lin, D. N. C., & Lau, H. H. B. 2012, ApJ, 758, L6 [NASA ADS] [CrossRef] [Google Scholar]
 Rossiter, R. A. 1924, ApJ, 60, 15 [NASA ADS] [CrossRef] [Google Scholar]
 Roxburgh, I. W., & Vorontsov, S. V. 2003, A&A, 411, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Salpeter, E. E. 1955, ApJ, 121, 161 [Google Scholar]
 SanchisOjeda, R., Winn, J. N., Holman, M. J., et al. 2011, ApJ, 733, 127 [NASA ADS] [CrossRef] [Google Scholar]
 SanchisOjeda, R., Winn, J. N., Marcy, G. W., et al. 2013, ApJ, 775, 54 [NASA ADS] [CrossRef] [Google Scholar]
 Schlaufman, K. C. 2010, ApJ, 719, 602 [NASA ADS] [CrossRef] [Google Scholar]
 Shkolnik, E., Walker, G. A. H., Bohlender, D. A., Gu, P.G., & Kürster, M. 2005, ApJ, 622, 1075 [NASA ADS] [CrossRef] [Google Scholar]
 Silva Aguirre, V., Ballot, J., Serenelli, A. M., & Weiss, A. 2011a, A&A, 529, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Silva Aguirre, V., Chaplin, W. J., Ballot, J., et al. 2011b, ApJ, 740, L2 [Google Scholar]
 Silva Aguirre, V., Casagrande, L., Basu, S., et al. 2012, ApJ, 757, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Silva Aguirre, V., Basu, S., Brandão, I. M., et al. 2013, ApJ, 769, 141 [NASA ADS] [CrossRef] [Google Scholar]
 Simon, T., & Landsman, W. 1991, ApJ, 380, 200 [Google Scholar]
 Skumanich, A. 1972, ApJ, 171, 565 [NASA ADS] [CrossRef] [Google Scholar]
 Spiegelhalter, D. J., Best, N. G., Carlin, B. P., & van der Linde, A. 2002, J. R. Stat. Soc. Ser. B Stat. Methodol., 64, 583 [Google Scholar]
 Steigman, G. 2010 [arXiv:1008.4765] [Google Scholar]
 Szabó, G. M., Szabó, R., Benkő, J. M., et al. 2011, ApJ, 736, L4 [NASA ADS] [CrossRef] [Google Scholar]
 Torrence, C., & Compo, G. P. 1998, Bull. Am. Meteorol. Soc., 79, 61 [Google Scholar]
 Torres, G., Fischer, D. A., Sozzetti, A., et al. 2012, ApJ, 757, 161 [NASA ADS] [CrossRef] [Google Scholar]
 Toutain, T., & Appourchaux, T. 1994, A&A, 289, 649 [NASA ADS] [Google Scholar]
 Triaud, A. H. M. J., Collier Cameron, A., Queloz, D., et al. 2010, A&A, 524, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Valsecchi, F., & Rasio, F. A. 2014, ApJ, 786, 102 [NASA ADS] [CrossRef] [Google Scholar]
 Van Eylen, V., Kjeldsen, H., ChristensenDalsgaard, J., & Aerts, C. 2012, Astron. Nachr., 333, 1088 [NASA ADS] [CrossRef] [Google Scholar]
 Van Eylen, V., Lindholm Nielsen, M., Hinrup, B., Tingley, B., & Kjeldsen, H. 2013, ApJ, 774, L19 [NASA ADS] [CrossRef] [Google Scholar]
 Van Eylen, V., Lund, M. N., Silva Aguirre, V., et al. 2014, ApJ, 782, 14 [NASA ADS] [CrossRef] [Google Scholar]
 van Saders, J. L., & Pinsonneault, M. H. 2013, ApJ, 776, 67 [NASA ADS] [CrossRef] [Google Scholar]
 Vogt, S. S., Allen, S. L., Bigelow, B. C., et al. 1994, in Instrumentation in Astronomy VIII, eds. D. L. Crawford & E. R. Craine, SPIE Conf. Ser., 2198, 362 [Google Scholar]
 Weber, E. J., & Davis, Jr., L. 1967, ApJ, 148, 217 [NASA ADS] [CrossRef] [Google Scholar]
 Weiss, A., & Schlattl, H. 2008, Ap&SS, 316, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Winn, J. N. 2010 [arXiv:1001.2010] [Google Scholar]
 Winn, J. N., Noyes, R. W., Holman, M. J., et al. 2005, ApJ, 631, 1215 [NASA ADS] [CrossRef] [Google Scholar]
 Winn, J. N., Johnson, J. A., Albrecht, S., et al. 2009, ApJ, 703, L99 [NASA ADS] [CrossRef] [Google Scholar]
 Winn, J. N., Fabrycky, D., Albrecht, S., & Johnson, J. A. 2010, ApJ, 718, L145 [NASA ADS] [CrossRef] [Google Scholar]
 Wright, J. T., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2004, ApJS, 152, 261 [NASA ADS] [CrossRef] [Google Scholar]
Online material
Appendix A: Powerspectrum modelling and optimisation
Appendix A.1: Modelling the power spectrum
To describe the observed power spectral density of a mode peak in the frequency power spectrum, we use a standard Lorentzian function (see, e.g., Anderson et al. 1990; Gizon & Solanki 2003) given by (A.1)The use of a Lorentzian function for the mode line profile comes from the nature of solarlike pmodes; the modes are stochastically driven by turbulent convection in the outer envelope after which they are intrinsically damped (see, e.g., Goldreich et al. 1994). In this equation, H_{nlm} is the mode height, ν_{nlm} is the resonance frequency of the mode, while Γ_{nl} is a measure of the damping rate of the mode and gives the FWHM of L_{nlm}(ν).
For slow stellar rotation the star is generally assumed to rotate as a rigid body and the modes will to first order be split as (Ledoux 1951) (A.2)Here m is the azimuthal order of the mode, Ω is the angular rotation rate of the star, and C_{nl} is a dimensionless constant that describes the effect of the Coriolis force (the Ledoux constant). For highorder lowdegree solarlike oscillations, like the ones we wish to analyse, this quantity is of the order C_{nl}< 10^{2}, and is therefore neglected. In this way we see that the splitting due to rotation between adjacent components of a multiplet will to a good approximation be given by ν_{s} = Ω/2π.
In assuming equipartition of power between the components of a multiplet (i.e. no assumed preference in the excitation for prograde over retrograde propagating modes), it is possible to calculate the geometrical modulation of the relative visibility between the 2l + 1 multiplet components as a function of i_{⋆} as (see, e.g., Dziembowski 1977; Gizon & Solanki 2003) (A.3)where are the associated Legendre functions.
With this, the limit spectrum (noisefree) to be fit to the power spectrum is as expressed in Eq. (1) (see Fig. 2). By comparing Eqs. (A.1) and (1), we see that the height is given by ℰ_{lm}(i)S_{nl}. By assuming equipartition of energy between different radial orders, this can be written as (A.4)The factor is a measure of the relative visibility in power (primarily set by partial cancellation) between nonradial and radial (l = 0) modes, while α_{l = 0}(ν) represents a (mainly) frequencydependent height modulation for the radial modes, generally represented by a Gaussian centred on the frequency of maximum oscillation power, ν_{max}.
Appendix A.2: Optimisation procedure
The fitting of Eq. (1) to the power spectrum is made in a Bayesian manner by mapping the posterior probability: (A.5)Here p(Θ  I) is the prior probability assigned to the parameters Θ from any prior information I, and p(D  Θ,I) is the likelihood of the observed data D given the parameters Θ. The posterior is approximated using the affine invariant MCMC sampler emcee (ForemanMackey et al. 2013; see also Hou et al. 2012). Using the emcee routine, the posterior distribution is mapped, after which parameter estimates are evaluated as the median of the respective marginalised distributions (see Sect. 4). In the sampling we enable the parallel tempering scheme of emcee and use five temperatures with tempering parameters set as β_{i} = 1.2^{1 − i} (Benomar et al. 2009; Handberg & Campante 2011). The affine invariant character of the emcee sampler ensures that it works efficiently in spite of linear parameter correlations, which are a problem for many MCMC algorithms. In our optimisation we make both a fit to a large (full fit) and a small (small fit) frequency range (see Sect. 4 for details). We employ 1500 (full fit) and 2000 (small fit) walkers, all initiated from a sampling of the prior distributions. Each walker is stopped after 10 000 steps, after which we thin the chains by a factor of 10 (full fit) or 5 (small fit). We cut away a burnin part of each chain based on the Geweke^{7} statistics (Geweke 1992), and check for good mixing using the autocorrelation time of the chains and by performing a visual inspection of the traces of walkers in parameter space. We refer to Handberg & Campante (2011) and references therein for further details on the MCMC nomenclature and ForemanMackey et al. (2013) for the specifics of the emcee sampler. To ensure better numerical stability we map the logarithm of the posterior with the description of the loglikelihood function from Anderson et al. (1990) and Toutain & Appourchaux (1994). With regard to priors, we use tophat priors for location parameters (e.g., ν_{nl}) and scale invariant modified Jeffryes’ priors for scale parameters (e.g., S_{n 0}). To decrease the computation time the limitspectrum (Eq. (1)) was only fit to the frequency range that includes the identified oscillation modes (see Sect. 4.1 for further details). To better constrain the stellar noisebackground in the relatively narrow range occupied by the oscillation modes, Eq. (2) was first fit to the power spectrum in the frequency range from 100 − 8496 μHz (the upper limit is the approximate Nyquist frequency of SC data) and included either one or two characteristic timescales corresponding to the contributions from granulation only or granulation and faculae (this lowerlimit frequency ensures that the activity component can be omitted). We also added a Gaussian function to Eq. (2) to account for the power excess from solarlike oscillations seen in HATP7. Using the deviance information criterion (DIC; Spiegelhalter et al. 2002), we found that only one component is needed to describe of the background. The medians of the posteriors from this fit were then used to fix the background in the fit of Eq. (1).
Appendix B: Line widths and visibilities
Appendix B.1: Line widths
The line widths can have a strong impact on the estimated splitting because a small splitting could be equally well fit by a slightly larger line width. This is especially a problem when the splitting is smaller than the line width. We show in Fig. B.1 our fitted line widths for the radial (l = 0) modes and their associated uncertainties. As a sanity check, we may first compare the mode line width at the frequency of maximum power, ν_{max}, with the estimate from Eq. (2) of Appourchaux et al. (2012) and using combined values from their Table 2. Using the spectroscopic temperature of T_{eff} = 6350 ± 126 K, yields Γ ≈ 3.5 ± 1.1 μHz. From the three l = 0 modes closest to ν_{max} we find (central) values of Γ between 5.5 and 6.4 μHz (see Fig. B.1), which is higher than expected from the Appourchaux et al. (2012) formulae.
Fig. B.1
Measured line widths for radialorder modes are given by the open symbols as a function of frequency (left axis), and plotted with associated errors. Theoretical linear damping rates, multiplied by two, are shown by the solid curve (right axis). The dotted curve shows smoothed damping rates, multiplied by a factor of 2.2. Inputs to the theoretical model were taken from the bestfit GARSTEC model, e.g. the radius at the base of the surface convection zone (R_{bcz}/R_{⋆} = 0.862). 
Additionally, we estimated linear damping rates, η, which we assumed to be approximately equal to half of the observed line widths, i.e. Γ ≃ ηπ^{1}, if η is in units of angular frequency. The outcome is shown in Fig. B.1. The computations included a full nonadiabatic treatment of the pulsations and convection dynamics. Both the convective heat and momentum (turbulent pressure) fluxes were treated consistently in the equilibrium and pulsation computations using the nonlocal generalisation of the timedependent convection model by Gough (1977a,b). The computations were carried out as described by Houdek et al. (1999) and Houdek & Gough (2002). For the nonlocal convection parameters we adopted the values a^{2} = 900 and b^{2} = 2000, and the mixinglength paramter was calibrated to obtain the same depth of the surface convection zone as in the best GARSTEC model. For the anisotropy parameter, Φ (see Houdek & Gough 2002), the value 2.50 was adopted.
After applying a median smoothing filter to ηπ^{1}, with a width in frequency corresponding to five radial modes, the result of 2.8 μHz at ν_{max} of the oscillation heights lies within the error bars of the observational scaling relation by Appourchaux et al. (2012). To fit the observations in Fig. B.1 we multiplied the mediansmoothed estimates by a factor of 2.2 (dotted curve in Fig. B.1), which agrees with previous comparisons between line width observations and model estimates for hotter solarlike stars (see, e.g., Houdek 2006, 2012).
Appendix B.2: Visibilities
For the visibilities (see Appendix A), we estimate and from the small fit. These agree reasonably well with the theoretical values of and estimated from the tables of Ballot et al. (2011). We do note, however, that this agreement is no guarantee for correct values because some stars deviate from the simple theoretical estimates (see, e.g., Lund et al. 2014). If we calculate the visibilities using the method described in Ballot et al. (2011), but adopt a quadratic limbdarkening (LD) law and measured LDparameters from fits to the planetary transit by Van Eylen et al. (2013) and Morris et al. (2013), we obtain values that are slightly lower than those from theoretical LD parameters: and . These values agree within the errors with the fit values, which is encouraging given the very simplified assumptions adopted in the Ballot et al. (2011) calculation, where for instance all nonadiabatic effects are neglected.
Appendix C: Determining (B – V)_{0}
We determined the photometric stellar parameters for HATP7 by combining asteroseismic results with the Infrared flux method (IRFM) (see Silva Aguirre et al. 2011b, 2012). We adopted our seismic log g and the spectroscopic metallicity from Huber et al. (2013b, 2014) and used the IRFM implementation described in Casagrande et al. (2014), where different threedimensional reddening maps are used to constrain extinction. At a distance of 320 pc (approximate distance to HATP7 determined by Pál et al. 2008), reddening varies between 0.02 <E^{(}B − V^{)}< 0.03.
Unfortunately, optical measurements of HATP7 are quite uncertain, and depending on whether the Tycho2 (Høg et al. 2000) or APASS (Henden et al. 2009) photometry is used, the resulting T_{eff} will vary anywhere between 6350 and 6650 K. At the magnitude of our star, Tycho2 photometry becomes increasingly uncertain (Høg et al. 2000) (although its T_{eff} would be in overall good agreement with the spectroscopic estimate of T_{eff} = 6350 ± 126 K), while the APASS (B − V)_{0} index (i.e., after correcting it for reddening) is almost as red as the solar one (Ramírez et al. 2012), thus suggesting a T_{eff} close to solar. The higher T_{eff} is, however, confirmed by the (reddening corrected) J − K_{s} index of HATP7 (indeed bluer than the solar one; Casagrande et al. 2012). From these considerations we thus discard the APASS (B − V)_{0} as faulty, and adopt a photometric T_{eff} = 6500 ± 150 K, where the generous errors account for the discussed uncertainties. With this T_{eff} and [Fe/H] = 0.26, we can invert the colourT_{eff}[Fe/H] relation of Casagrande et al. (2010), which returns an intrinsic (i.e., unreddened) colour of (B − V)_{0} = 0.455 ± 0.040 mag. For spectroscopic T_{eff} and uncertainty, its (B − V)_{0} = 0.495 ± 0.022.
Synthetic photometry offers an alternative way of assessing the (B − V)_{0} colour for HATP7 (all synthetic quantities are obtained by interpolating at the “known” physical parameters of the star, and are thus unaffected by reddening). We use the large grid of MARCS (Gustafsson et al. 2008) synthetic colours and interpolation routines provided by Casagrande & VandenBerg (2014) to infer the (B − V)_{0} index of HATP7, for the spectroscopic parameters and asteroseismic log g.
To estimate the uncertainty in the synthetic (B − V)_{0}, we also compute all possible T_{eff} and [Fe/H] combinations allowed by the spectroscopic uncertainties (while the seismic log g is so precisely known that changing it makes no difference). With this procedure, we obtain (B − V)_{0} = 0.498 ± 0.020, in excellent agreement with the estimate from the empirical colourT_{eff}[Fe/H] relation when using the spectroscopic T_{eff}.
For the gyrochronology calculation we used the value (B − V)_{0} = 0.495 ± 0.022 from the IRFM when using the spectroscopic T_{eff}.
Appendix D: Peakbagging results
Frequencies extracted from the MCMC peakbagging.
All Tables
Selected properties of our bestfit models compared with the first asteroseismic modelling by ChristensenDalsgaard et al. (2010, CD10).
All Figures
Fig. 1
Configuration of the HATP7 system to scale using values for HATP7b (full black circle) from Van Eylen et al. (2013). Top: the observer’s view of the system from Earth, with the angular momentum vectors of the planetary orbit, (normal to the orbital plane), and stellar spin, , given as red arrows. The projected angle, λ, is indicated by the shaded magenta region and is found as the angle between and when these are projected onto the plane of the sky (dashed lines). This angle is obtained from RM measurements, and in this panel we have used λ = 155 ± 37° following Albrecht et al. (2012; the uncertainty on λ is not included in the figure). The stellar inclination, i_{⋆}, which is the parameter measured from asteroseismology, is set to 15°, and is given by the direct angle between the line of sight (midpoint of star) and . The inclination of the planetary orbit, i_{p}, is set to 83° following Van Eylen et al. (2013). The true angle, ψ, is the direct angle between and . Middle: side view of the system, with the observer’s viewpoint from the right, indicated by “line of sight”. To properly show i_{⋆} (shaded green) and i_{p} (shaded red), and not their projected values, we have set λ = 180° such that both and lie on the same plane as the line of sight. For illustrative purposes we have in addition decreased i_{p} to 63°. When adopting this configuration ψ (shaded blue) is given by the sum of i_{⋆} and i_{p}. Bottom: top view of the system, with the observer’s viewpoint from the bottom, indicated by “line of sight”. 

In the text 
Fig. 2
Left: power spectrum of HATP7 (black) with the fitted model overplotted (red; cf. Eq. (1)). Right: power spectrum of HATP7 (black) overplotted with the best fit to the background (red; cf. Eq. (2)). The lightgrey part up to 100 μHz was not included in the fit. The fit includes a granulation component (P_{G}; green) and a white/shot noise (P_{S}; black) level in addition to the Gaussian envelope from pmodes centred around ν_{max} ≈ 1115 μHz. The dashed red line shows the background fit without the Gaussian envelope. 

In the text 
Fig. 3
Échelle diagram for HATP7 using Δν = 59.22 μHz. The grey scale ranges from white at low power to black at high power. Circles give the extracted frequencies with corresponding uncertainties, triangles connected by red lines give model frequencies from GARSTEC, while squares connected by green lines give model frequencies from MESA. The degree of each ridge is indicated in the top part of the plot. The radial order of the l = 0 modes is indicated by the numbers at the right side of the plot. 

In the text 
Fig. 4
Ratios r_{02} and r_{010} as a function of frequency. The lines show the ratios obtained for the bestfit GARSTEC (black solid) and MESA (dashed blue) models. Note that the relatively small uncertainties on frequencies (horizontal error bars) renders them indiscernible on this scale. 

In the text 
Fig. 5
Hinton diagram for the correlation matrix of r_{010} ratios. The first part of the subscript denotes whether the ratio is of type 01 or 10, while the second gives the radial order of the central frequency. White (black) squares indicate positive (negative) covariances between the ratios in question, and the size gives the relative size of the correlation (one along the diagonal). 

In the text 
Fig. 6
Top: marginalised posterior probability distribution (PPD) for the stellar inclination i_{⋆}. Here we have folded the full distribution (−90 to 180°) onto the the range from 0 to 90°. Bottom right: PPD for the rotational splitting. Bottom left: correlation map between the inclination and the rotational splitting. The 68% credible regions (highest posterior density credible regions) are indicated by the darkgrey parts, while lightgrey indicates the additional parts covered in a 95% credible region. Included are also lines of constant v sin i_{⋆}, computed using the radius estimate from our analysis and with v sin i_{⋆} values from the literature (see Table 3), arranged in the legend in order of increasing v sin i_{⋆} value (top to bottom). The dashed black horizontal line indicates the value of ν_{s} = 5 μHz above which walkers were excluded from the PPDs. The grey dashdotted lines give lines of constant v sin i_{⋆} from 8 to 48 km s^{1} in steps of 8 km s^{1}. 

In the text 
Fig. 7
Distribution for v sin i_{⋆} constructed from the PPDs obtained for ν_{s} and i_{⋆} and assuming (see Eq. (9)). The dark and lightgrey regions of the distribution correspond to the 68% and 95% credible regions. Literature values for v sin i_{⋆} obtained from spectral and RM analysis (see Table 3) are given by vertical lines where the shaded regions above the distribution give the corresponding uncertainties. 

In the text 
Fig. 8
Distribution of the true angle ψ using the distribution for i_{⋆} from the peakbagging, and with the assumption of normal distributions for the planetary inclination i_{p} and projected angle λ. For i_{p} the value from Van Eylen et al. (2013) was adopted, while λ was taken from Albrecht et al. (2012). The dark and lightgrey regions of the distribution correspond to the 68% and 95% credible regions. The dashed curve gives the distribution assuming instead an isotropic distribution for i_{⋆}, i.e., flat in cos i_{⋆}. We indicate also which values of ψ correspond to a retrograde or prograde orbit of HATP7b. 

In the text 
Fig. 9
Gyrochronology relations from Meibom et al. (2009; M09) and Barnes (2007; B07). The shaded region around each relation (with the same colour) represents the standard error of the relationship between (B − V)_{0} and rotation period from propagating the uncertainties reported for the coefficients entering Eq. (11) together with the uncertainty in our age estimate. The vertical line and shaded region give the value of (B − V)_{0} = 0.495 ± 0.022 of HATP7, while darker shaded regions around the relations vertically bound the corresponding uncertainty in the period from the uncertainty in (B − V)_{0}, age, and the gyrochronology relations. The horizontal lines indicate the limits on the 68% (dashed) and 95% (dotted) credible regions on the stellar rotation period from pmode splittings. The upwardpointing arrows indicate that these rotation periods are lower limits from the respective HPD credible regions. 

In the text 
Fig. 10
Top: autocorrelation function (ACF) of the HATP7 time series (black). The dashed horizontal line gives the expectation value for random independent and identically distributed values, while the grey curves give the largelag 95% confidence levels. Bottom: lowfrequency end of the power spectrum in units of period. 

In the text 
Fig. 11
Morlet wavelet power spectrum as a function of time for HATP7, where the time series was binned by 31 points. The colour changes from red at low power to blue at high power – the colours are on a logarithmic scale. The crosshatched regions at high and low times indicate the cone of influence where edge effects become important (see Torrence & Compo 1998). No clear signatures from rotation are apparent. 

In the text 
Fig. B.1
Measured line widths for radialorder modes are given by the open symbols as a function of frequency (left axis), and plotted with associated errors. Theoretical linear damping rates, multiplied by two, are shown by the solid curve (right axis). The dotted curve shows smoothed damping rates, multiplied by a factor of 2.2. Inputs to the theoretical model were taken from the bestfit GARSTEC model, e.g. the radius at the base of the surface convection zone (R_{bcz}/R_{⋆} = 0.862). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.