Issue |
A&A
Volume 506, Number 1, October IV 2009
The CoRoT space mission: early results
|
|
---|---|---|
Page(s) | 255 - 262 | |
Section | Stellar atmospheres | |
DOI | https://doi.org/10.1051/0004-6361/200811487 | |
Published online | 05 March 2009 |
The CoRoT space mission: early results
Photospheric activity and rotation of the planet-hosting star CoRoT-4a![[*]](/icons/foot_motif.png)
A. F. Lanza1 - S. Aigrain2 - S. Messina1 - G. Leto1 - I. Pagano1 - M. Auvergne3 - A. Baglin3 - P. Barge4 - A. S. Bonomo1,4,5 - A. Collier Cameron6 - G. Cutispoto1 - M. Deleuil4 - J. R. De Medeiros7 - B. Foing8 - C. Moutou4
1 - INAF - Osservatorio Astrofisico di Catania, via S. Sofia, 78, 95123 Catania, Italy
2 - School of Physics, University of Exeter, Stocker Road, Exeter, EX4 4QL, UK
3 - LESIA, CNRS UMR 8109, Observatoire de Paris, 5 place J. Janssen, 92195 Meudon, France
4 - Laboratoire d'Astrophysique de Marseille (UMR 6110),
Technopole de Château-Gombert,
38 rue Frédéric Joliot-Curie,
13388 Marseille Cedex 13, France
5 - Dipartimento di Fisica e Astronomia, Università degli Studi di Catania, via S. Sofia, 78, 95123 Catania, Italy
6 - School of Physics and Astronomy, University of St. Andrews,
North Haugh, St Andrews, Fife KY16 9SS, Scotland
7 - Departamento de Física, Universidade Federal do Rio Grande do Norte, 59072-970 Natal, RN, Brazil
8 - ESA/ESTEC/SRE-S, Postbus 299, 2200 AG Noordwijk, The Netherlands
Received 9 December 2008 / Accepted 28 January 2009
Abstract
Aims. The space experiment CoRoT has recently detected a transiting hot Jupiter in orbit around a moderately active F-type main-sequence star (CoRoT-4a). This planetary system is of particular interest because it has an orbital period of 9.202 days, the second longest one among the transiting planets known to date. We study the surface rotation and the activity of the host star during an uninterrupted sequence of optical observations of 58 days.
Methods. Our approach is based on a maximum entropy spot modelling technique extensively tested by modelling the variation in the total solar irradiance. It has been successfully applied to model the light curve of another active star with a transiting planet observed by CoRoT, i.e., CoRoT-2a. It assumes that stellar active regions consist of cool spots and bright faculae, analogous to sunspots and solar photospheric faculae, whose visibility is modulated by stellar rotation.
Results. The modelling of the light curve of CoRoT-4a reveals three main active longitudes with lifetimes between
and
days that rotate quasi-synchronously with the orbital motion of the planet. The different rotation rates of the active longitudes are interpreted in terms of surface differential rotation, and a lower limit of
is derived for its relative amplitude. The enhancement of activity observed close to the subplanetary longitude suggests a magnetic star-planet interaction, although the short duration of the time series prevents us from drawing definite conclusions.
Conclusions. The present work confirms the quasi-synchronicity between stellar rotation and planetary orbital motion in the CoRoT-4 system and provides a lower limit for the surface differential rotation of the star. This information can be important in trying to understand the formation and evolution of this highly interesting planetary system. Moreover, there is an indication of a possible star-planet magnetic interaction that needs to be confirmed by future studies.
Key words: stars: magnetic fields - stars: late-type - stars: activity - stars: rotation - planetary systems - stars: individual: CoRoT-4a
1 Introduction
CoRoT has recently discovered CoRoT-4b, a Jupiter-sized planet transiting the disc of an F-type main-sequence star with an orbital period of 9.202 days (Aigrain et al. 2008). This is the second longest period among the transiting planetary systems known to date, putting CoRoT-4b in a region of the mass-period parameter space that was previously empty (cf. Moutou et al. 2008). Moreover, the out-of-transit light curve shows a modulation with an amplitude of a few 0.001 mag that can be attributed to photospheric brightness inhomogeneities carried into and out of view by the rotation of the star. Given its late spectral type, those inhomogeneities can be considered analogous to the cool spots and bright faculae observed in the Sun, which owe their existence to photospheric magnetic fields.
The active regions of CoRoT-4a are stable enough to allow an estimate of its rotation period through an autocorrelation analysis, yielding a period of
days. This result indicates that the stellar rotation and the orbital motion of the planet are quasi-synchronized, which cannot be explained on the basis of the tidal theory given the large separation between the star and the planet (
)
and its low mass (
,
cf., Aigrain et al. 2008). Other synchronous or quasi-synchronous star-planet systems have recently been found, consisting of quite massive planets orbiting F-type stars (Catala et al. 2007; McCullough et al. 2009). Among those systems,
Bootis
has been studied in some detail revealing a surface differential rotation with a relative amplitude
,
and hints of an activity cycle of a few years (Donati et al. 2008; Catala et al. 2007).
However, those systems host more massive and closer planets that may explain the synchronization by tidal effects, assuming that only the outer stellar convection zone, having a mass of only a few 0.001 of the stellar mass, is synchronized (cf. Donati et al. 2008). On the other hand, the estimated time scale for the synchronization of the outer convective envelope of CoRoT-4a is
Gyr, assuming a mass of 0.001 of the stellar mass for the envelope.
In the present work, we present a modelling of the out-of-transit light curve of CoRoT-4a following the approach already applied to another star hosting a transiting hot Jupiter, i.e., CoRoT-2a. It is based on a solar analogy; i.e., it assumes that the stellar active regions responsible for the flux modulation consist of cool spots and bright faculae (Lanza et al. 2009).
Late-type stars accompanied by hot Jupiters sometimes show active regions rotating with the orbital period of the planet instead of the stellar rotation period. The phenomenon has been detected in some seasons in
Andromedae and HD 179949 that have shown chromospheric hot spots leading the planet by
and
,
respectively (Shkolnik et al. 2008). Moreover, the synchronous system
Boo has shown some evidence of
an active region rotating with the orbital period of its planet, but leading it by
during the period 2001-2005 (Walker et al. 2008). These phenomena suggest a magnetic interaction between the star and its hot Jupiter, as discussed by Lanza (2008). In the case of CoRoT-4a, the short duration of the observations prevents us from coming to definite conclusions about a possible star-planet interaction, but the longitude distributions of the active regions may be compared with the orbital position of the planet to see whether there is some analogy, say, with
Boo.
A detailed study of stellar activity by means of Doppler imaging techniques is not feasible owing to the low projected rotation velocity of CoRoT-4a (
km s-1) and its faintness (
). Therefore, a spot modelling approach is the most suitable for investigating the magnetic activity and rotation of this very interesting object.
2 Observations
CoRoT-4a was observed during the CoRoT Initial Run of 58 days, starting 6 February 2007. The time sampling was initially 512 s, then reduced to 32 s during the last three weeks after the detection of the transits by the so-called CoRoT alarm mode. CoRoT performed aperture photometry with a fixed mask. Since the star was brighter than R = 14.5, the flux was split along detector column boundaries into broad-band red, green, and blue channels.
The observations and data processing are described by Aigrain et al. (2008), to whom we refer the reader for details. The reduction pipeline applied corrections for the background and the pointing jitter of the satellite, which was particularly relevant during ingress and engress from the Earth shadow. Measurements during the crossing of the South Atlantic Anomaly of the Earth's magnetic field, which amounted to about 15-20 percent of a satellite orbit, were discarded. More information on the instrument, its operation, and performance can be found in Auvergne et al. (2009).
The light curve extracted from the CoRoT Mission Archive was further corrected for the effect of a hot pixel in the blue channel, as explained by Aigrain et al. (2008). To increase the signal-to-noise ratio and reduce residual systematic effects possibly present in individual channels, we summed up the flux in the red, green, and blue channels to obtain a light curve in a spectral range extending from 300 to 1100 nm. The oversampled section of this white-band light curve was rebinned with a regular 512 s sampling, to obtain an evenly sampled time series. It was further cleaned by applying a moving median 5-
clipping filter that allowed us to identify and discard residual outliers, as described by Aigrain et al. (2008). The standard deviation of the observations, averaged over 512 s bins, is
in relative flux units.
The contamination of the photometric aperture by stars other than CoRoT-4 was only
percent of the median of the measured flux, and it was subtracted to avoid dilution of the stellar light variations. Transits were removed by means of the ephemeris and parameters of Aigrain et al. (2008).
The flux variations related to stellar activity have time scales of a day or longer. Therefore, we rebinned the out-of-transit light curve by computing normal points along each orbital period of the satellite (6184 s). This has the advantage of removing tiny systematic variations that may still be associated with the orbital motion of the satellite (cf. Auvergne et al. 2009; Alonso et al. 2008). Since the orbital period is not a multiple of the individual exposures of 32 or 512 s, we performed a linear interpolation of the flux variation along the last 40 s of each orbital period to compute the normal points. In such a way we obtained a light curve consisting of 684 normal points ranging from HJD 2 454 135.0918 to HJD 2 454 192.7744, i.e., with a duration of 57.6826 days.
The light curve shows a long-term decreasing trend (cf. Fig. 2 in Aigrain et al. 2008) that may be instrumental in origin (Auvergne et al. 2009). Since it would produce only a monotonous increase of the spotted area in our model along the days of the present observations, we decided to subtract it before modelling the light curve by fitting a straight line to the data. Finally, the de-trended light curve was normalized to its maximum flux value observed at
HJD
=2 454 192.2730, which we assumed to represent the unspotted flux level of the star, whose true value is unknown.
3 Spot modelling
The reconstruction of the surface brightness distribution from the rotational modulation of the stellar flux is an ill-posed problem, because the variation of the flux vs. rotational phase contains information only on the distribution of the brightness inhomogeneities vs. longitude. The integration over the stellar disc cancels any latitudinal information, particularly when the inclination of the rotation axis along the line of sight is
,
as in the present case (see Sect. 4 and Lanza et al. 2009). Therefore, we need to include a priori information in the light curve inversion process to obtain a unique and stable map. This is done by computing a maximum entropy (hereinafter ME) map, which has been proven to best represent the case of the Sun (cf. Lanza et al. 2007). As a matter of fact, a sequence of seven maps was obtained for CoRoT-4a, covering its successive rotations during the 58 days of the observations. They will allow us to study the evolution of the active longitudes over the surface of the star.
Readers not interested in the details of our spot modelling approach may skip the remainder of this section and go on to Sect. 4 (or directly to Sect. 5) .
In our model the star is subdivided into several surface elements, in our case 200 squares of side
,
with each element containing unperturbed photosphere, cool spots, and facular areas. The fraction of an element covered by cool spots is indicated by the filling factor f, the fractional area of the faculae is Qf, and the fractional area of the unperturbed photosphere is 1-(Q+1)f. We fit the light curve by changing the value of f over the surface of the star, while Q is held uniform and constant. Even fixing the rotation period, the inclination, and the spot and facular contrasts (see Lanza et al. 2007, for details), the model has 200 free parameters and suffers from non-uniqueness and instability. To find a unique and stable spot map, we apply maximum entropy regularization, as described in Lanza et al. (2007), by minimizing a functional
,
which is a linear combination of the
and the entropy functional S; i.e.,
![]() |
(1) |
where





















Table 1: Parameters adopted for the spot modelling of CoRoT-4a.
In the case of the Sun, by assuming a fixed distribution of the filling factor, it is possible to obtain a good fit of the irradiance changes only for a limited time interval
,
not exceeding 14 days which is the lifetime of the largest sunspot groups dominating the irradiance variation. In the case of other active stars, the value of
must be determined from the observations themselves, looking for the maximum data extension that allows for a good fit with the applied model (see Sect. 4 for CoRoT-4a).
The optimal values of the spot and facular contrasts and of the facular-to-spotted area ratio Q in stellar active regions are unknown a priori. In our model the facular contrast
and the parameter Q enter as the product
,
so we can fix
and vary Q, estimating its best value
by minimizing the
of the model, as shown in Sect. 4. Since the number of free parameters of the ME model is large, for this specific application we make use of the model of Lanza et al. (2003), which fits the light curve by assuming only three active regions to model the rotational modulation of the flux plus a uniformly distributed background to account for the variations in the mean light level. This procedure is the same as was adopted to fix the value of Q in the case of CoRoT-2a
(cf. Lanza et al. 2009).
As for CoRoT-2a, we assume an inclination of the rotation axis of CoRoT-4a of
,
which implies that our model cannot provide information on the latitude distribution of the active regions. Such a limitation cannot be overcome by assuming an inclination
in the model, owing to the relatively low information content of the present light curve. Specifically, when we assume
,
we find that the ME regularization virtually puts all the spots at the sub-observer latitude (i.e.,
)
to minimize their area and maximize the entropy. Therefore, we are limited to mapping only the distribution of the active regions vs. longitude, which can be done with a resolution of at least
(cf. Lanza et al. 2009,2007). Our ignorance of the true value of the facular contribution to the light variations may produce systematic errors in the active region longitudes derived by our model, as discussed by Lanza et al. (2007) in the case of the Sun.
4 Model parameters
The fundamental stellar parameters are taken from Aigrain et al. (2008) and Moutou et al. (2008)
and are listed in
Table 1. A quadratic limb-darkening law is adopted for the stellar photosphere, viz.,
,
where
is the specific intensity at the limb position
,
with
the angle between the normal to a surface element and the line of sight (cf. Lanza et al. 2007,2003). The limb-darkening parameters have been derived from Kurucz (2000) model atmospheres for
K,
(cm s-2) and solar abundances, by adopting the CoRoT white-band transmission profile in Auvergne et al. (2009).
![]() |
Figure 1:
The ratio of the |
Open with DEXTER |
The polar flattening of the star due to the centrifugal potential is computed in the Roche approximation with a rotation period of 9.202 days (see below). The relative difference between the equatorial and the polar radii is
which may induce a flux variation about two orders of magnitude smaller for a spot coverage of
percent as a consequence of the gravity darkening of the equatorial regions of the star.
The inclination of the stellar rotation axis is difficult to constrain through the observation of the Rossiter-McLaughlin effect because of the small
of the star and its intrinsic line profile variations due to stellar activity (Moutou et al. 2008). Nevertheless, we assume that the stellar rotation axis is normal to the orbital plane of the planet, i.e., at an inclination of
from the line of sight (cf. Aigrain et al. 2008).
The rotation period adopted for our spot modelling is equal to the orbital period of the planet. This allows us to check the synchronization of stellar rotation to the planetary orbit through the study of the longitude drift of the active regions versus time (cf. Sect. 5). Another advantage of this choice is that the subplanetary longitude is fixed in our reference frame, which allows us to investigate the possible association between active longitudes and the orbital position of the planet. Nevertheless, considering that the autocorrelation analysis of the light modulation gives a rotation period of
days (Aigrain et al. 2008), it is worth also investigating spot models with a rotation period of 8.87 days.
![]() |
Figure 2: Upper panel: the ME-regularized composite best fit to the out-of-transit light curve of CoRoT-4a obtained for Q=4.5. The flux is in relative units, i.e., measured with respect to the maximum observed flux in the de-trended normal point light curve. Lower panel: the residuals from the composite best fit versus time. In both panels, the dashed vertical lines indicate the epochs of mid transits. |
Open with DEXTER |
The maximum time interval that our model can accurately fit with a fixed distribution of active regions has been determined by means of a series of tests and has been found to be
days. Therefore, we subdivided the time series into seven equal segments of duration
days, each of which was modelled with a fixed active region pattern.
To compute the spot contrast, we adopted the same mean temperature difference as derived for sunspot groups from their bolometric contrast, i.e. 560 K (Chapman et al. 1994). In other words, we assumed a spot effective temperature of 5630 K, yielding a constrast
in the CoRoT passband (cf. Lanza et al. 2007).
A different spot contrast changes the absolute spot coverages, but does not significantly affect their longitudes and their time evolution, as discussed in detail by Lanza et al. (2009). The facular contrast was assumed to be solar-like with
(Lanza et al. 2004).
The best value of the area ratio Q between the faculae and the spots in each active region has been estimated by means of the model by Lanza et al. (2003, cf. Sect. 3#. In Fig. 1, we plot the ratio
of the
total
of the composite best fit to the entire time series to
its minimum value
,
versus Q, and indicate the 99 percent confidence level as derived from the F-statistics (e.g., Lampton et al. 1976). The best value of Q turns out to be Q=4.5, with an acceptable range extending from
to
.
For comparison, the best value in the case of the Sun is
,
indicating a lower relative contribution of the faculae to the total light variation in CoRoT-4a.
The amplitude of the rotational modulation of our star is
mag, i.e., 2-3 times that of the Sun at the maximum of the eleven-year cycle (see Sect. 5). This suggests that CoRoT-4a is more active than the Sun, which may account for the reduced facular contribution to its light variations, as suggested by Radick et al. (1998) and Lockwood et al. (2007).
5 Results
The composite best fit obtained with the ME regularization is shown in Fig. 2, together with the residuals. The systematic negative values of the residuals are a consequence of the ME regularization that tends to reduce the spotted area as much as possible, thus systematically driving the fit above the observations (cf. Sect. 3).
The residuals of the unregularized composite best fit have a mean of
and a standard deviation
in relative flux units, i.e., assuming the maximum observed flux as the unit of measure.
In the absence of regularization, the standard deviation of the residuals is more than expected from the photon shot noise, i.e.,
(Aigrain et al. 2008). This may arise from intrinsic stellar microvariability on short time scales (up to a few days) or may be due to some residual instrumental effects.
On average, we have
normal points per fitted subset. The Lagrangian multipliers for the regularized ME models have been fixed to obtain a mean of the residuals of the composite fit
,
corresponding to
.
The optimal value
was established by considering the information content of the maps obtained for different values of
and minimizing it (i.e., maximizing the entropy) as much as possible, while still retaining an average acceptable fit. Lower values of
give an improvement of the fit over small sections of the light curve, but introduce several small spots that vary from one subset to the next. Therefore, they are actually not required by the data. On the other hand, by increasing the entropy of the maps as much as possible, we obtain smoother maps with a regular variation in the spot pattern from one subset to the next, which indicates that we are properly modelling the overall active region distribution on the star.
The two spikes of the best fit at HJD
and
stem from the small variations of the spot configuration when we move from one data subset to the next. The reader may wonder about the possibility of improving the regularized best fit by considering shorter time intervals
.
However, the improvement is negligible because the misfits at the matching points between one subset and the next are a consequence of the significant amount of regularization required by this light curve due to its relatively low information content. The effect is actually much greater in the present case than for CoRoT-2a, whose light curve had a significantly greater information content than the present one (cf. Lanza et al. 2009). Therefore, decreasing
in the present case does not alleviate the problem, but produces the opposite effect, owing to more matching points along the time series.
![]() |
Figure 3:
The distributions of the spotted area vs. longitude at the labelled times (
|
Open with DEXTER |
The distributions of the spotted area vs. longitude are plotted in Fig. 3 for the seven mean epochs of our individual subsets adopting a rotation period of 9.20205 days. The longitude zero corresponds to the point intercepted on the stellar photosphere by the line of sight to the centre of the star at HJD 2 454 135.0918, i.e., the sub-observer point at the initial epoch. The longitude increases in the same direction as the stellar rotation and the orbital motion of the planet, which has a fixed subplanetary point on the star in the adopted reference frame. This choice allows a direct comparison of the active region longitudes with the subplanetary longitude, although one cannot correlate the active regions to the dips in the light curves in a straightforward way.
Three main active regions can be identified in Fig. 3 and their migration has been traced with different straight lines. Specifically, their rotation rates with respect to the adopted reference frame are found by a linear best fit, i.e., assuming a constant migration rate. They are listed in Table 2 for three values of the facular-to-spotted area ratio Q, because the longitudes of the active regions derived from our model depend somewhat on that parameter (cf. Lanza et al. 2009,2007).
The relative amplitude of the surface differential rotation, estimated from the difference between the greatest and the lowest migration rates, is
,
,
and
for
Q=0, 4.5, and 7.0, respectively. For the sake of completeness, we have also computed ME spot models adopting Q=4.5 and a rotation period of 8.870 days, which comes from the autocorrelation analysis of the light curve (Aigrain et al. 2008). Again, three main active longitudes are seen in the longitude distribution of the spotted area. Their migration rates are listed in Table 2, and the derived relative amplitude for the differential rotation is
,
in very good agreement with the result obtained for a rotation period of 9.20205 days. It is interesting to note that a rotation period of 8.870 days makes the migration rate of AR2, i.e., the active longitude with the greatest spotted area, virtually zero. This agrees with the fact that the autocorrelation of the light curve is maximized for that rotation period.
Unfortunately, no information is available on the spot latitudes in our models, therefore, our
are only lower limits. These values suggest that CoRoT-4a has a surface differential rotation comparable to that of the Sun, for which we can estimate a relative amplitude of
0.04-0.05 considering active regions confined within the sunspot belt, i.e., within
from the equator.
The lifetimes of the active regions can be estimated from our spot models and range from 10 to 50 days, with longer lifetimes characterizing those having greater filling factors, as in the Sun.
Some active regions are located close to the subplanetary longitude most of the time, although the highest concentrations of spotted (and facular) areas are around a longitude following the planet by about
.
Conversely, the hemisphere centred around longitude
has the minimum active region coverage (cf. Fig. 3).
The association between the subplanetary longitude and an enhanced spot activity is particularly remarkable in the case of the models computed without facular contribution, i.e., for Q=0 (see Fig. 4). In six out of seven epochs, the ME models indicate a relative maximum of the spotted area within
from the subplanetary longitude. The time average of the spotted area in each longitude bin is plotted in Fig. 5 vs. longitude, together with the corresponding standard deviation. When Q=0 (Fig. 5, lower panel), a well-defined relative maximum appears close to the subplanetary longitude, with a height of
standard deviations, which is unlikely to be due to a chance fluctuation. Moreover, the subplanetary peak has an FWHM significantly lower than of the active longitude peak around
whose width reflects the migration of its active regions in time. A greater migration rate of the active regions is also responsible for smearing the active longitude around
and the subplanetary longitude, when we consider the average of the models computed with Q=4.5 (Fig. 5, upper panel).
The total spotted area varies only slightly with time (see Fig. 6) and may indicate a level of activity about 2-3 times that of the Sun at the maximum of the eleven-year cycle. The spotted area is determined here after removing the long-term decrease of the flux observed in the original light curve by Aigrain et al. (2008). Therefore, it does not show the long-term linear increase that one would expect in case the long-term trend were not removed. Moreover, the absolute value of the area depends on the adopted spot contrast
and the value of Q. For instance, a lower spot temperature would imply a stronger contrast and thus a smaller area, but the relative variations of the area are largely independent of
and Q (cf. Lanza et al. 2009).
Table 2:
Relative migration rates
of the active longitudes for different facular-to-spotted area ratio Q and rotation period
.
![]() |
Figure 4: The same as Fig. 3, but for Q=0.0. |
Open with DEXTER |
![]() |
Figure 5:
Upper panel: the spotted area per
|
Open with DEXTER |
![]() |
Figure 6: The total spotted area as derived from the regularized ME models vs. time for Q=4.5. |
Open with DEXTER |
6 Discussion
Our results confirm that the rotation of CoRoT-4a is quasi-synchronized with the orbital motion of its planet. Minimum spot migration rates are obtained for a rotation period of 8.87 days, which corresponds to the rotation period of the active longitude with the highest spot concentration. However, even when adopting that period, the other active longitudes continue to display a significant migration. Our interpretation of the different migration rates of the active longitudes is that the corresponding active regions are located at different latitudes on a differentially rotating star. Nevertheless, the different drift rates of the active longitudes may be produced by a phenomenon different from surface differential rotation. For instance, in the framework of a solar analogy, the rotation rate of sunspot groups during the early stages of their evolution (< 2-3 days from their appearance) is about 2 percent greater than that of recurrent groups and decreases steadily during their later evolution (e.g., Zappalà & Zuccarello 1991). However, we believe that it is unlikely that a similar phenomenon may explain the results obtained for CoRoT-4a because our active longitude drift lasts for several dozens of days. In conclusion, surface differential rotation in combination with active regions located at different latitudes seems to be the most likely explanation for the results presented in Sect. 5. Although the observed active longitudes rotate faster than the planet, the existence of a latitude exactly synchronized to the planetary orbit cannot be excluded because we have no information on the latitude range covered by the active regions.
The quasi-synchronous rotation status of CoRoT-4a might have been attained when the star was still magnetically locked to the inner parts of a circumstellar disc,
although already close to the ZAMS. If the planet had already migrated to the radius where the disc was truncated by the stellar magnetic field
(Romanova & Lovelace 2006), the system started its evolution in a rotational and orbital configuration close to the present status. The subsequent tidal interaction between the star and the planet could not significantly change the angular momentum of the planet and the star, as noticed by Aigrain et al. (2008). The magnetic braking of the star after the dissipation of the disc is expected to be negligible given its mid-F spectral type and an initial rotation period of
days
(Wolff & Simon 1997).
The surface differential rotation found through our spot modelling is not remarkably dependent on the adopted value of the facular-to-spotted area ratio Q and is at least 6-7 times that derived for CoRoT-2a that is rotating in 4.52 days and has spectral type G7V (Lanza et al. 2009).
For
Bootis, with
K and a rotation period of
days, Catala et al. (2007) found a relative differential rotation
between the equator and the pole by means of line-profile modelling. This compares well with our value, given that our determination is only a lower limit. The same method was applied by Reiners (2006) to a sample of F-type stars finding projected pole-equator relative amplitudes that are comparable to that of CoRoT-4a when one considers objects with a similar effective temperature and
km s-1. The small
of our star prevents the application of line profile analysis methods to measure its surface differential rotation (Reiners 2006), therefore,
a spot modelling approach is the only possibility of estimating its surface shear.
Theoretical predictions for the amplitude of the surface differential rotation in main-sequence F-type stars have recently been published by Küker & Rüdiger (2005). Their model is in general agreement with the present findings although the lack of a precise theory for the interaction between stellar rotation and turbulent convection makes a detailed quantitative comparison premature.
The possibility that stellar activity is influenced by the hot Jupiter in CoRoT-4a is suggested by
the present models, in particular by those assuming that the photometric modulation is only produced by cool spots (i.e., with Q=0). The plots in Fig. 5 suggest that the formation of spots may be triggered in an active region located at the subplanetary longitude; i.e., the emergence of magnetic flux may be promoted there in some way by the close planet.
It is important to note that tidal effects cannot account for this phenomenon because there is only one region of enhanced activity on the star most of the time, and not two separated by
,
as expected in that case.
The possibility that a hot Jupiter may in some way affect the dynamo action in the subsurface layers of the convection zone of its host star has been discussed by Lanza (2008), and the present results give some support to his speculations.
However, a word of caution is certainly appropriate in the present case, because the time interval covered by our data is rather limited. Therefore, the supposedly enhanced activity related to the planet might still be the result of statistical fluctuations in the appearance of small active regions over the star. Moreover, the relationship between enhanced activity and subplanetary longitude appears to be less clear when we consider models with Q > 0, also because the migration rates of the main active longitudes are greater.
7 Conclusions
Our analysis shows that the rotation of CoRoT-4a is quasi-synchronized with the orbital motion of the planet. Specifically, the active longitude that has the largest spotted area shows an angular velocity that is only 3-5 percent higher than the orbital angular velocity. Assuming that it traces the angular velocity of stellar rotation at some latitude, this translates into an upper limit for the difference between stellar and orbital angular velocities at that latitude.
Moreover, we find evidence of a significant differential rotation of the star, with a lower limit
,
when the optimal facular-to-spotted area ratio Q=4.5 is adopted.
It cannot be excluded that the planet may affect the longitude distribution of the active regions on the photosphere of the star in some way. Specifically, spot activity seems to be increasing close to the subplanetary longitude, particularly for the models computed without faculae. Unfortunately, the short duration of the present time series (only 58 days, i.e., 6.3 stellar rotations) does not allow us to derive definite conclusions on such a star-planet magnetic interaction.
A long-term monitoring of the star from the ground may help to clarify this issue, e.g., through measurements of the rotational modulation of the chromospheric Ca II H & K flux, which should be a good indicator of a possible star-planet interaction, as in the case of
Boo or
And (Walker et al. 2008; Shkolnik et al. 2008).
Acknowledgements
The authors are grateful to the Referee, Prof. Gordon Walker, for a careful reading of the manuscript and several valuable comments. This work has been partially supported by the Italian Space Agency (ASI) under contract ASI/INAF I/015/07/0, work package 3170. Active star research and exoplanetary studies at INAF-Osservatorio Astrofisico di Catania and Dipartimento di Fisica e Astronomia dell'Università degli Studi di Catania are funded by MIUR (Ministero dell'Istruzione, dell'Università e della Ricerca) and by Regione Siciliana, whose financial support is gratefully acknowledged. This research has made use of the ADS-CDS databases, operated at the CDS, Strasbourg, France.
References
- Aigrain, S., Collier Cameron, A., Ollivier, M., et al. 2008, A&A, 488, L43 [CrossRef] [EDP Sciences] (In the text)
- Alonso, R., Auvergne, M., Baglin, A., et al. 2008, A&A, 482, L21 [NASA ADS] [CrossRef] [EDP Sciences]
- Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411 [CrossRef] [EDP Sciences] (In the text)
- Catala, C., Donati, J.-F., Shkolnik, E., Bohlender, D., & Alecian, E. 2007, MNRAS, 374, L42 [NASA ADS]
- Chapman, G. A., Cookson, A. M., & Dobias, J. J. 1994, ApJ, 432, 403 [NASA ADS] [CrossRef] (In the text)
- Donati, J.-F., Moutou, C., Farès, R., et al. 2008, MNRAS, 385, 1179 [NASA ADS] [CrossRef]
- Küker, M., & Rüdiger, G. 2005, A&A, 433, 1023 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Kurucz, R. L. 2000, http://kurucz.harvard.edu/ (In the text)
- Lampton, M., Margon, B., & Bowyer, S. 1976, ApJ, 208, 177 [NASA ADS] [CrossRef] (In the text)
- Lanza, A. F. 2008, A&A, 487, 1163 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Lanza, A. F., Catalano, S., Cutispoto, G., Pagano, I., & Rodonò, M. 1998, A&A, 332, 541 [NASA ADS] (In the text)
- Lanza, A. F., Rodonò, M., Pagano, I., Barge P., & Llebaria, A. 2003, A&A, 403, 1135 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Lanza, A. F., Rodonò, M., & Pagano, I. 2004, A&A, 425, 707 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Lanza, A. F., Bonomo, A. S., & Rodonò, M. 2007, A&A, 464, 741 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Lanza, A. F., Pagano, I., Leto, G., et al. 2009, A&A, 493, 193 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Lockwood, G. W., Skiff, B. A., Henry, G. W., et al. 2007, ApJS, 171, 260 [NASA ADS] [CrossRef] (In the text)
- McCullough, P. R., Burke, C. J., Valenti, J. A., et al. 2009, ApJ, in press, [arXiv:0805.2921]
- Moutou, C., Bruntt, H., Guillot, T., et al. 2008, A&A, 488, L47 [CrossRef] [EDP Sciences] (In the text)
- Radick, R. R., Lockwood, G. W., Skiff, B. A., & Baliunas, S. L. 1998, ApJS, 118, 239 [NASA ADS] [CrossRef] (In the text)
- Reiners, A. 2006, A&A, 446, 267 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Romanova, M. M., & Lovelace, R. V. E. 2006, ApJ, 645, L73 [NASA ADS] [CrossRef] (In the text)
- Shkolnik, E., Bohlender, D. A., Walker, G. A. H., & Collier Cameron, A. 2008, ApJ, 676, 628 [NASA ADS] [CrossRef] (In the text)
- Walker, G. A. H., Croll, B., Matthews, J. M., et al. 2008, A&A, 482, 691 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Wolff, S. C., & Simon, T. 1997, PASP, 109, 759 [NASA ADS] [CrossRef] (In the text)
- Zappalà, R. A., & Zuccarello, F. 1991, A&A, 242, 480 [NASA ADS] (In the text)
Footnotes
- ... CoRoT-4a
- Based on observations obtained with CoRoT, a space project operated by the French Space Agency, CNES, with partecipation of the Science Programme of ESA, ESTEC/RSSD, Austria, Belgium, Brazil, Germany, and Spain.
- ... Archive
- http://idoc-corot.ias.u-psud.fr/
All Tables
Table 1: Parameters adopted for the spot modelling of CoRoT-4a.
Table 2:
Relative migration rates
of the active longitudes for different facular-to-spotted area ratio Q and rotation period
.
All Figures
![]() |
Figure 1:
The ratio of the |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Upper panel: the ME-regularized composite best fit to the out-of-transit light curve of CoRoT-4a obtained for Q=4.5. The flux is in relative units, i.e., measured with respect to the maximum observed flux in the de-trended normal point light curve. Lower panel: the residuals from the composite best fit versus time. In both panels, the dashed vertical lines indicate the epochs of mid transits. |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
The distributions of the spotted area vs. longitude at the labelled times (
|
Open with DEXTER | |
In the text |
![]() |
Figure 4: The same as Fig. 3, but for Q=0.0. |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Upper panel: the spotted area per
|
Open with DEXTER | |
In the text |
![]() |
Figure 6: The total spotted area as derived from the regularized ME models vs. time for Q=4.5. |
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.