Free Access
Issue
A&A
Volume 613, May 2018
Article Number A50
Number of page(s) 12
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/201732350
Published online 30 May 2018

© ESO 2018

1 Introduction

The detection of giant planets around young stars can address key questions in the astrophysics of planetary formation and migration. In fact, looking at young stars provides an opportunity to observe the architecture of planets in their infancy. A variety of physical processes, such as planet-disc interaction, the Kozai mechanism (Kozai 1962), planet-planet scattering (Baruteau et al. 2014), or in-situ formation (Batygin et al. 2016), responsible for generating hot Jupiters (HJs), are expected to produce observable effects, including differences in orbital parameters (eccentricity and/or obliquity), migration timescales, an age-dependent frequency of such systems, and differences in the atmospheric composition of hot Jupiters across stellar ages. The investigation of the formation and migration histories may be carried out for example via intensive radial velocity (RV) monitoring of a sample of young stars, finding new planets and testing possible differences in the frequency of hot Jupiters with age. However, the high level of stellar activity characterizing young stars induces RV variations able to mimic planetary signals, demanding a specific processing in order to investigate its actual contribution. Recently, noteworthy results have been obtained through the modelling of stellar activity as correlated noise in a set of data with the Gaussian processes regression (see e.g. Haywood et al. 2014; Damasso et al. 2018). The first pioneering detection of a young HJ with the RV technique was recently announced around the weak-line T Tauri star V830 Tau using spectropolarimetry (Donati et al. 2016; Yu et al. 2017; Donati et al. 2017), after years of attempts marred by difficulties in dealing with the identification of planetary signals (both spectroscopic and photometric) in the presence of very high levels of activity of their young hosts. Simultaneous multi-band high resolution spectroscopy (HRS), for example in the visible (VIS) and near-infrared (NIR, 700–2500 nm) ranges, is a powerful tool for disentangling wavelength-dependent activity-induced RV signals (the impact of activity on RVs is expected to be typically three times lower in the NIR than in the VIS; see Prato et al. 2008; Mahmud et al. 2011; Crockett et al. 2012)from those Keplerian that are achromatic. An emblematic example of this principle is represented by the case of the HJ around TW Hya claimed by Setiawan et al. (2008), which was ruled out by NIR observations (Figueira et al. 2010a).

The building of new instrumentation allowing high precision HRS in the NIR range has been mainly pushed by the RV search of potentially habitable rocky planets around M-dwarf stars (see e.g. Reiners et al. 2010; Bonfils et al. 2013 and references therein). Some of those instruments, such as CARMENES (520–1710 nm, with spectral resolutions R = 80 000–100 000, Quirrenbach et al. 2014; Reiners et al. 2018) or HARPS+NIRPS (0.95–1.8 nm, R ~ 100 000, Conod et al. 2016), also guarantee contemporary VIS-NIR observations. In this framework, the GIARPS (GIANO–B + HARPS–N) project (Claudi et al. 2017) has been conceived precisely to exploit the full potential of the simultaneous VIS and NIR HRS technique, allowing us to have the two high resolution spectrographs of the Telescopio Nazionale Galileo (TNG, La Palma) working simultaneously: HARPS–N (High Accuracy Radial velocity Planet Searcher for the Northern hemisphere) in the visible and GIANO–B in the NIR (see details in Sect. 2). The RV precisions achievable by the VIS and NIR channels of GIARPS (see more details in Sect. 2) allow us to identify giant planets signals as well as to discriminate the high activity levels of the host stars, enabling an overall monitoring of stellar activity and simultaneous detection of atomic and molecular species in planetary atmospheres.

In this paper, we present the RV measurements of a debated HJ around BD+20 1790 (V429 Gem, K5Ve, V = 9.9, Jeffries 1995), a very active star (logRHK=3.7$\log R'_{HK}=-3.7$, Hernán-Obispo et al. 2015, hereafter HO15) and probable member of the AB Dor moving group (Torres et al. 2008) with an adopted age of 14919+51$149^{+51}_{-19}$ Myr (Bell et al. 2015). After dedicated spectroscopic observations, Hernán-Obispo et al. (2010; hereafter HO10) interpreted the RV variation of this target as being due to the presence of a massive HJ with a period of 7.8 days. Figueira et al. (2010b) questioned the planet, providing CORALIE RVs showing a clear correlation with the bisector span (BIS), thus attributing the RV variations to photospheric processes. Later, HO15 reported a new RV and activity study exploiting a larger spectroscopic and photometric dataset supporting again the presence of a HJ around BD+20 1790. Finally, Gagné et al. (2016) published a small number of CSHELL NIR RVs with no conclusive results.

This paper is organized as follows: we describe the instruments used in Sect. 2, then we present the observations and the characteristics of our dataset in Sect. 3. In Sect. 4 we show the analysis of both the spectroscopic and photometric time series, suggesting a theoretical model to explain the observed stellar behaviour in the VIS and NIR bands in Sect. 5. Finally we draw our conclusions in Sect. 6.

2 Instruments

In this section we summarize the characteristics of the instrumentation used to collect the data of BD+20 1790, with particular emphasis on the new TNG facility GIARPS. In 2012 the high resolution spectrograph HARPS–N (Cosentino et al. 2014) was installed at the Nasmyth–B focal station of the TNG. HARPS–N works in the visible range (0.39−0.68 μm) with a resolution of R = 115 000. In the same year, the NIR high resolution spectrograph GIANO (Oliva et al. 2006), working in the wavelength range from 0.97 to 2.45 μm at a resolution of R 50 000, was installed and commissioned at the Nasmyth – A focal station of the TNG in 2014. Although designed for direct light feed from the telescope, in its first phase GIANO was fed by fibers. Through the “WOW” (a Way to Others Worlds) Progetto Premiale funding scheme, the Italian National Institute for Astrophysics (INAF) proposed to move GIANO to the Nasmyth–B (re-naming the instrument GIANO–B) and to carry out its simultaneous use with HARPS–N aiming to achieve high resolution spectroscopy in a wide wavelength range (0.383−2.45 μm with a small gap between ~0.7 and 0.950 μm) obtained in a single exposure. The focus change, made in 2016, also allowed us to restore the original direct feeding from the telescope (Tozzi et al. 2016) and the coupling with HARPS–N. This was the beginning of the observing mode called GIARPS (GIANO–B & HARPS–N, Claudi et al. 2017). The two spectrographs are still able to work separately, so it is possible to have three different observing modes: a) HARPS–N only; b) GIANO–B only; c) GIARPS, by splitting the light with a dichroic. GIARPS uses both the instruments for high precision RV measurements, exploiting the simultaneous reference technique with HARPS–N in the visible (providing a RV precision σRV of ~ 1.0 m s−1) and a cross-correlation function (CCF) method based on telluric lines with GIANO–B (σRV ~10 m s−1, Carleo et al. 2016). The introduction of an absorption gas cell (Seeman et al. 2014), foreseen at the end of 2018, will allow us to reach a better precision (~3 m s−1) in the NIR, since the spectral features of a gas cell are more reliable in comparison with the instability of the telluric spectrum. Because of its characteristics, GIARPS can be considered the first and unique worldwide instrument providing not only high resolution but also high precision radial velocity measurements in such a large wavelength range between B and K bands.

The Immersion Grating Infrared Spectrometer (IGRINS; Mace et al. 2016) is the cross-dispersed NIR spectrograph, mounted alternatively at the Harlan J. Smith 2.7 m telescope (McDonald Observatory, TX, USA) and at the 4.3 m Discovery Channel Telescope (Lowell Observatory, AZ, USA), with a resolving power of R = 45 000 (Yuk et al. 2010; Park et al. 2014). It covers the H and K windows, from 1.45 to 2.5 μm in a single acquisition. REM (Rapid Eye Mount; Chincarini et al. 2003) is a 60 cm robotic telescope located at the La Silla station of the European Southern Observatory (ESO, Chile). The telescope hosts two instruments: REMIR, an infrared imaging camera, andROS2, a visible imager. The two cameras can also observe simultaneously the same field of view thanks to a dichroic placed in front of the telescope focus. The ROS2 camera is equipped with a back-illuminated CCD (2048 × 2048 pixels, 13.5 micron pixels size, 0.58 arcsec∕pixel plate scale), which has a corrected 9.8 × 9.8 arcmin field of view, and observes simultaneously through the four Sloan∕SDSS g,r,i,z filters. The REMIR camera is equipped with a Hawaii I CCD (512 × 512 pixels, 1.2 arcsec∕pixel plate scale), which has a corrected 10 × 10 arcmin field of view, and it is equipped with zJHK filters.

Table 1

Summary of the spectroscopic data presented in this work.

3 Observations and data reduction

In this section we describe the data acquired from each instrument. The spectroscopic datasets are presented in Sects. 3.13.3, and a complete summary of them, including the uncertainties of the RVs measurements (obtained as explained before), the typical signal-to-noise ratio (S/N), the RV rms, and the peak-to-valley amplitude, is reported in Table 1. The photometric observations are described in Sect. 3.4.

3.1 GIANO/GIANO–B

We collected spectra of BD+20 1790 with TNG instrumentation in three different observing campaigns (see Table 1). The first dataset consists of 18 observations with GIANO (March 28–April 3, 2016), supported by quasi-simultaneous photometry with the REM telescope (see below). The second dataset (seven spectra) was acquired during the commissioning of GIANO–B in November 2016. The third one contains five spectra of GIANO–B obtained during the GIARPS commissioning in March 2017. Two of them are acquired in GIARPS mode in order to test the simultaneity of GIANO–B and HARPS–N observations. One additional HARPS–N spectrum was later collected without the NIR counterpart because GIANO–B was temporarily unavailable.

GIANO data were reduced with the IRAF package ECHELLE and the dedicated scripts collected in the GIANO_TOOLS1 package, while GIANO–B spectra are processed with the dedicated pipeline GOFIO (Rainer et al., in prep.). We obtained a set of 30 NIR RVs over one year (listed in Table A.1), with the method described in Carleo et al. (2016), in which the telluric linesare used as wavelength reference and the CCF method is used to determine the stellar RV. For this purpose, we constructed two suitable digital masks that include about 2000 stellar lines and a similar number of telluric lines. After the correction of the spectra to the barycentre of the solar system, the procedure performs the cross correlation of individual orders of the normalized spectra with the appropriate masks (both stellar and telluric), with the derivation of individual CCF. Through a weighted sum of the CCFs of the individual orders, we obtain the final stellar and telluric CCFs. These are fitted with Gaussian profiles to derive the stellar and telluric RV, respectively. The latter are finally subtracted from the former, providing the relative stellar RVs. The uncertainties are then evaluated taking into account the photon statistics. In the present analysis we consider a slightly different approach, which takes into account the weight of the single orders, some of them being affected by telluric lines and thus contributing in different ways to the determination of the RVs. Starting from the RV of individual orders, we calculate the weighted average RV for each exposure and its corresponding error. As a final step, we derive the bisector velocity span (BIS, as in Carleo et al. 2016) of the CCF and we calculate the uncertainties on this quantity by considering the fractions of the CCF used for the derivation of the BIS, resulting in 10σRV$\sqrt{10}{\sigma_{\textrm{RV}}}$, where σRV is the RV error.

3.2 HARPS–N

HARPS–N RVs are extracted with the usual data reduction software (DRS, Pepe et al. 2002) by cross-correlating the observed spectrum with a numerical mask that depicts the spectral features of a K5 star. To work with uniform RV values, we processed our HARPS–N spectra and the ones collected by HO15 (except for one spectrum with very low S/N at JD 2456681) with the current HARPS–N DRS through the YABI workflow (Hunter et al. 2012) installed at IA22 at the INAF Observatoryof Trieste. Since this facility allows us to customize the re-processing, we enlarge the width of the CCF to take into account the quite large vsini of this star (10.03±0.47 km s−1, López-Santiago 2005) and the consequent line broadening. The resulting RVs are listed in Table A.2, together with the CaII activity indicator, log RHK, obtained with the dedicated tool of the HARPS–N DRS (the method is provided in Lovis et al. 2011) also available on YABI, and the BIS (see e.g. Queloz et al. 2001) estimated as in Lanza et al. (2018), starting from the computed CCF of the DRS. The usually adopted uncertainty for the BIS is twice the value of the RV uncertainty, on the basis that the BIS slope is calculated using the top and lower half of a single line measurement. However, this occurs in the ideal case in which the bins of the CCF used for the estimation of the bisector are not correlated, therefore, in our case the multiplying factor is increased to 2.5 to avoid an underestimation of the uncertainties (see the justification in Lanza et al. 2018).

3.3 IGRINS

The 29 spectra of BD+20 1790 collected with IGRINS over one year, from April 2016 to March 2017, are reduced with the IGRINS Pipeline Package3. We acquired three different datasets, since the run with IGRINS installed at the Discovery Channel Telescope (DCT) occurred in between the two runs at the McDonald (McD) Observatory (see Table 1). The RV measurements (Table A.3) are obtained with the same procedure used for GIANO.

3.4 REM

During the first run with GIANO we obtained quasi-simultaneous photometry with REM from May 30 to April 4, 2016. We used tasks within IRAF for bias correction and flat fielding, and the technique of aperture photometry to extract magnitude time series for BD+20 1790 and for other stars detected in the frames, which were selected as candidate comparison stars. In particular we identified two stars that were found to be non-variable and were used as comparison (C) and check (CK) stars (C: 2MASS 07233899+2025102, J = 11.23 mag, H = 11.43 mag, K = 11.17 mag; CK: 2MASS 07234597+2025328, J = 12.15 mag, H = 11.97 mag, K = 11.91 mag). We measured a standard deviation σCCK = 0.014 mag in their differential light curve. The visible and infrared magnitudes of BD+20 1790 were computed differentially with respect to the comparison star. After averaging the six consecutive differential magnitudes obtained on each night with ROS2, we obtained a time seriesof six average griz-band differential magnitudes and four JHK-band differential magnitudes for the subsequent analysis. The averagestandard deviation associated with the nightly averaged magnitudes, which we consider as our photometric precision, is σ < 0.015 mag for all filters, and σ ≃ 0.04 mag forthe H and K magnitudes.

4 Data analysis

4.1 Spectroscopic data

First, we reproduced the orbital fit of the VIS RVs as presented in HO15 (with SARG, FOCES, and HERMES data only) with a Keplerian function4 as in Desidera et al. (2011). Our model is displayed in the upper panel of Fig. 1 as a dashed line together with the VIS data used for the fit, represented by grey dots. The obtained RV semi-amplitude is 926.8 ± 34.6 ms−1 and the period is 7.7827 ± 10−4 days, in good agreement with the result proposed by HO15. As in HO15, we were not able to obtain a solution including their HARPS–N data so we only over-imposed them to the fitting function. We observe a match that by itself casts some doubt on a Keplerian interpretation of the RV variation.

Our phase-folded NIR RVs are shown as well in Fig. 1 (lower panel): GIANO and GIANO–B RVs (red dots) show an rms scatter of 129.6 ms−1, while IGRINS RVs (light blue dots) show an rms scatter of 109.2 ms−1. The rms scatter of the whole NIR dataset is 119 ms−1. Finally, the black asterisks represent new HARPS–N data (two of them have the NIR simultaneous counterpart). Figure 1 shows that the amplitude of NIR RVs (calculated as the difference between maximum and minimum RVs) is 437.3 ms−1, four times lower with respect to the optical one reported in HO15 and interpreted as a signature of a hot Jupiter. Therefore, according to our data we can exclude any companion with those characteristics, ascribing the observed variation to phenomena of stellar origin.

Actually, the Spearman rank correlation between RVs and BIS, in this case for the whole HARPS–N dataset (HO15 and the three epochs presented here) is − 0.92, with a very high statistical significance (p-value = 2 × 10−8, evaluated through the IDL Astronomy Library routine SAFE_CORRELATE; see Table 2 for a summary of the measured correlations between RVs and activity indices of the whole dataset), showing an unambiguous linear trend (Fig. 2). After subtracting this correlation from the HARPS–N RVs time series, the resulting residuals show an rms of 61.1 ms−1 and the Generalized Lomb-Scargle (GLS) periodogram (Zechmeister & Kürster 2009) does not show any significant periodicity (i.e. with an amplitude larger than four times the standard deviation of the power spectrum, corresponding to a false alarm probability (FAP) larger than 0.01). As further evidence that the original HARPS–N RVs are modulated by the stellar rotation of 2.801 ± 0.001 d (the photometric period by HO10), in Fig. 3 we fit them with a Keplerian function by using as a first guess that period, obtaining a very good agreement with our model (P =2.80150 ± 2 × 10−5 d, K = 0.449 ± 7.370 km s−1, e = 0.065 ± 0.018).

We also investigated the other activity indices available from HARPS–N, for example logRHK$\log R&#x0027;_{HK}$ as described in Sect. 3, and the CCF asymmetry indices: the Δ V (a measure of the RV shift produced by the exclusive contribution of the asymmetry of the CCF; see Nardetto et al. 2006 and Figueira et al. 2013) and the Vasy(mod) (the modified version of the index Vasy by Figueira et al. 2013 for which the dependence from the RV shift is removed), presented in Lanza et al. (2018).5 Finally, we also checked the correlation with the Hα index derived as in Sissa et al. (2016). As for the BIS, the other asymmetry indices show significant correlation with RVs, while we find weak correlations with logRHK$\log R&#x0027;_{HK}$ and Hα indicators (Table 2).

We then focused our attention on the NIR data. We computed the GLS for both GIANO and IGRINS RVs even if the sampling is not suitable for a proper periodogram analysis (too few and too sparse points) or for a proper resolution of the known photometric period. The periodograms are quite noisy and do not exhibit any significant periodicity. Only for GIANO did we investigate the FAP (estimated by generating 10 000 artificial RV curves obtained from the real one, keeping the epochs of observations fixed but making random permutations of the RV values) of a peak at7.74 d responsible for the apparent signal of the GIANO data when phase-folded with the orbital period proposed by HO15 (red dots in Fig. 1, lower panel). This test returned a FAP of 24.4%, so this periodicity is probably produced by random noise. We also noticed a small amount of power corresponding with the photometric period of the star, but it appears to be related to the GIANO data sampling according to the analysis of the window function. No prominent periodicity is found either when examining the GLS of the whole NIR dataset, which is obtained by applying a quite negligible RV offset between the two instruments (5 m s−1).

In order to investigate a possible correlation between NIR RVs and activity indicators, we first measured the Spearman coefficient for GIANO and GIANO–B data, in particular between the BIS and the RVs, the BIS and HeI index at 1.083 μm (extracted as in Robertson et al. 2016), and between the HeI index and the RVs, but we found no strong correlation (Table 2), mainly due to the uncertainties in the measurements. A similar analysis was performed for the IGRINS data. Since the HeI line is out of the spectral coverage of IGRINS, we investigated the impact of the activity on RVs through the Brackett γ (Brγ) emission line at 2.16 μm (not available in the GIANO spectra because of the discontinuity among the orders), extracted as in Robertson et al. (2016). As in the case of GIANO, no significant correlation was found between RVs and these indicators (Table 2). Apparently our NIR radial velocities for BD+20 1790 are not highly sensitive to activity.

thumbnail Fig. 1

Orbital fit at 7.78 days found by HO15 compared to phase-folded visible and NIR RVs. Top panel: orbital fit (black dashed line)obtained with the visible data (FOCES, SARG, and HERMES RVs from HO15, grey dots) and HARPS-N 2015 RVs (green dots). Bottom panel: orbital fit (black dashed line), GIANO/GIANO-B (red dots), IGRINS (light blue dots), and HARPS-N 2017 (black asterisks, two acquired in GIARPS mode) RVs.

Table 2

Summary of the values of Spearman (ρ) correlation coefficients and corresponding p-values between RVs and activity indicators.

thumbnail Fig. 2

Correlation between RVs and BIS for HARPS–N data (HO15 and the three new epochs in this work).

thumbnail Fig. 3

Phase-folded HARPS–N RVs (2015, reprocessed from HO15 dataset) at stellar rotational period.

4.2 Photometric data

We used theGLS and the CLEAN (Roberts et al. 1987) periodogram analyses to search for significant periodicities in the BD+20 1790 photometric time series related to its rotation period. As an example, the GLS and CLEAN periodograms are plotted in Fig. 4 for the case of the g filter. The solid black line represents the normalized power versus period, whereas the dotted red line is the spectral window function. The horizontal red dashed line represents the power level corresponding to a FAP = 0.01 (confidence level of 99%, obtained with 1000 moke light curves as in Sect. 4.1).

Our periodogram analysis of the photometric variation confirmed the already known stellar rotational period P = 2.76± 0.04 days (corresponding to the most powerful and significant peak in the periodogram, indicated with a red mark), slightly lower than the literature value possibly because of the short baseline of the observations (six consecutive nights) or the effect of differential rotation. The uncertainty on the period is calculated following the prescription of Lamm et al. (2004). We note that the secondary power peak in the GLS periodogram is absent in the CLEAN periodogram, which has the capability of effectively removing beat frequencies arising from the data sampling.

In Fig. 5 we plot the differential light curves of BD+20 1790 phased with the rotation period P = 2.76 d, using different colours for different filters. The solid lines represent sinusoidal fits to the phased magnitude computed using the rotation period. The peak-to-peak amplitudes of the sinusoidal fits measured in different filters are reported as labels in the figure. We observe that the REM light curves show different amplitudes in the differentbands: unexpectedly, the amplitudes at longer wavelengths are greater than those at shorter wavelengthsand the NIR modulation is almost in anti-phase with respect to the optical modulation in the gri passbands (the maximum of the light curves are located at ~0.8 for the gri filters, at ~ 0.1 for z and at ~ 0.2 for JHK). Moreover, the amplitudes of the optical light curves are a factor of approximately seven smaller than in other seasons (cf. the light curves in Fig. 3 in HO15). This could be due to a transient phase of peculiar activity during the REM observing season (that covered a very short period) with respect to what was observed in past campaigns. This particular variability is also supported by the comparison between our HARPS–N RVs and the ones in HO15, as noticed also by Figueira et al. (2010b) with CORALIE data.

thumbnail Fig. 4

Top panel: generalized Lomb-Scargle periodogram of the g-band photometric time series of BD+20 1790. The solid black line isthe normalized power versus period, the dotted red line is the spectral window, and the horizontal dashed line indicates the power level corresponding to a FAP = 0.01. Bottom panel: CLEAN periodogram. The red mark is the most powerful and significant peak in the periodogram.

5 Discussion

The present investigation adds another important piece of evidence to the multi-wavelength characterization of the RV variations in late-type stars. Our results show that the amplitude of the NIR RV modulation in BD+20 1790 is generally a factor of 2−3 smaller than in the optical band, in agreement with the results obtained by Crockett et al. (2012) for very active young stars. This indicates that brightness inhomogeneities, whose contrast is generally smaller in the NIR than in the optical passband, are mainly responsible for the RV variations via line profile distortions. Other effects, such as quenching of convective blueshifts (e.g. Lanza et al. 2011) or line profile distortions produced by the Zeeman effect (Reiners et al. 2013), are probably less important in these very active and rapidly rotating (v sin i ≥ 10 km s−1) stars. Nevertheless, an intriguing result is the small amplitude of the RV variations in the NIR as measured by GIANO contemporaneously to the REM observations that show a remarkable NIR rotational modulation in the J, H, and K passbands. To interpret these results, we first consider a simple model for the wide-band photometric variations, including the effects of both dark and bright spots. Several models for the simultaneous photometric and RV variations of late-type stars have been proposed (e.g. Boisse et al. 2012; Dumusque et al. 2014; Herrero et al. 2016). They include the effects of solar-like faculae whose contrast increases towards the limb. Here, we consider a hot spot in the photosphere that has a constant contrast at different limb positions, similarly to the behaviour generally assumed for a cool spot.

From the stellar vsini, radius, and rotation period (from HO10), we estimated an inclination of the stellar spin axis to the line of sight of ~ 50°. With a simulation we then reproduced the sinusoidal shapes, amplitude ratio, and phase difference of the optical and the NIR light curves by assuming a circumpolar active region, always in view, consisting of two co-spatial components (as observed in e.g. V410 Tau, Rice et al. 2011): a cool feature (hereinafter “cool spot”) covering a fraction fs of its total area and a hot feature (hereinafter “hot spot”) covering the remaining fraction 1 − fs.

Therefore, the average brightness of the active region at wavelength λ can be written as Ia=fsB(Tc,λ)+(1fs)B(Th,λ),\begin{equation*} I_{\textrm{a}} = f_{\textrm{s}} B(T_{\textrm{c}}, \lambda)&#x002B; (1-f_{\textrm{s}}) B(T_{\textrm{h}}, \lambda), \end{equation*}(1)

where, for simplicity, we assume that the brightness of each component is given by a Planck function B(T, λ), with Tc being the temperature of the cool spot and Th that of the hot spot. Those temperatures verify the inequality: Tc < Tphot < Th, where Tphot is the temperature of the unperturbed photosphere. The contrast of the active region is Cs = 1 − (IaIphot), where Iphot is the brightness of the unperturbed photosphere. For an active region dominated by the cool spot, Cs > 0 because Ia < Iphot, while for an active region dominated by the hot spot, Cs < 0.

For BD+20 1790, we assume Tphot = 4410 K (HO15), while for the cool and the hot spots we assume temperatures Tc = Tphot − 1000 K and Th = Tphot + 1000 K, respectively. Those temperature differences are typical of young and active stars such as Weak T-Tauri stars (cf. Rice et al. 2011; Koen 2016). The contrast is plotted versus the fraction of the active region covered by the cool spot in Fig. 6 for the optical wavelength λopt = 636 nm and the NIR wavelength λNIR = 1705 nm. Those values correspond approximately to the mean wavelengths of the gri bands and of the JHK bands, respectively.

We see that for a restricted range of fs, that is 0.65 < fs < 0.67, the contrast in the optical is negative and small, while that in the NIR is positive and remarkably larger. This leads to a rotational modulation of the optical flux remarkably smaller than, and in anti-phase with, the rotational modulation in the NIR as illustrated by the synthetic light curves in Fig. 7 computed with the model in Sect. 3.2 of Lanza (2016). Specifically, these light curves were computed for an inclination of the spin axis to the line of sight i = 50°, considering an active region of an area of 20% of that of the star’s disc, centred at a latitude of 60°. The quadratic limb-darkening coefficients at the two wavelengths were taken from Claret et al. (2012), while Cs (λopt) = −0.02 and Cs (λNIR) = 0.15 corresponding to fs ≃ 0.665 in our simple irradiance model (cf. Fig. 6). For simplicity, we assumed that the spot contrasts did not depend on the position on the disc.

The present model is simply illustrative. The amplitudes of the synthesized light curves are remarkably smaller than those observed in BD+20 1790, which may require a larger filling factor of the active region and/or larger temperature contrasts. For example, in the case of LkCa 4, Gully-Santiago et al. (2017) found a filling factor as large as 86%, which would imply an amplitude larger by a factor of approximately four in the case of the present model. We finally notice that the variable characteristics of the activity of this star with time might explain why past studies obtained quite different values for the rms scatter of the RVs, and then different interpretations of the nature of this object.

As previously mentioned, the small amplitude of the RV modulation in the NIR as measured by GIANO contemporaneously with REM photometry is intriguing. Assuming that the depth of the spectral lines relative to their adjacent continuum is constant and considering a spot with a contrast Cs at latitude ϕ with a filling factor fs, we expect an RV modulation approximately of Csfsvsinicosϕ (cf. Saar et al. 1997; Desort et al. 2007) that is ≈ 0.75−1.0 km s−1 for a spot at ϕ = 60°, which is remarkably higher than what has been observed. A cool spot at a higher latitude would reduce the amplitude of the NIR wide-band flux modulations; a quenching of the convective shifts or the Zeeman effect also do not appear to be viable explanations because they increase the effect of a cool spot on the RV at NIR wavelengths (Reiners et al. 2013). Nevertheless, the variation of the relative depths of the spectral lines in the NIR cannot be neglected and it is the dominant effect in the cool spot responsible for the large photometric modulation in the infrared. The relative line depths are strong functions of the continuum opacity and of the degree of element ionization both remarkably varying in a cool spot area with respect to the unperturbed photosphere. In general, these effects produce a remarkable increase of the relative depths of the spectral lines in the cool spot. This compensates for the decrease of the continuum intensity in the spot, reducing the distortions of the spectral line profiles and yielding an RV variation in the NIR smaller than expected from the wide-band photometric variation in the case of constant relative line depths. In any case, a quantitative analysis is not warranted by our data, since a larger number of observations would be required. This scenario suggests a need to better investigate this kind of target, since they might go through specific activity phases during which the VIS and NIR RV amplitudes are similar, possibly resulting in false positives. Looking at the curve phase shifts might give crucial information in these cases (see e.g. the recent result by Hatzes et al. 2018 for the K-giant γ Draconis).

thumbnail Fig. 5

Differential lightcurves of BD+20 1790 phased with the rotation period P = 2.76 d in different photometric bands. The solid line is a sinusoidal fit to the data with the same period. Labels show the peak-to-peak amplitudes of the lightcurves.

thumbnail Fig. 6

Contrast of an active region consisting of co-spatial cool and hot spots versus the fraction of its area covered by the cool spot. Two wavelengths are considered: λopt = 636 nm (solid line) and λNIR = 1705 nm (dashed line).

thumbnail Fig. 7

Light curves showing the rotational modulation of the flux at two wavelengths λopt = 636 nm (solid line) and λNIR = 1705 nm (dashed line). The amplitudes correspond to a spot having an area of 0.2 of the stellar disc.

6 Conclusions

In this paper we present the analysis of RV measurements of the star BD+20 1790, in order to resolve the debate on the presence of a hot Jupiter with an orbital period of 7.8 days, claimed by HO10 and questioned by Figueira et al. (2010b). Since all the previous RV measurements were provided by visible spectrographs, we observed this target with the GIANO (in a dedicated run and during the GIANO–B/GIARPS Commissioning) and IGRINS NIR spectrographs to disentangle the origin of the RV variation. The NIR data show a peak-to-valley amplitude of 437.3 ms−1, significantly lower than the VIS one (1853.6 ms−1), demonstrating that the detected RV modulation is due to activity. Therefore, we dismiss the presence of the claimed hot Jupiter around BD+20 1790, which is a very peculiar target from the point of view of stellar activity, as confirmed by our photometric monitoring with REM almost contemporary to the first spectroscopic run with GIANO. Surprisingly, REM light curves show larger amplitudes at longer wavelengths and the NIR curve is almost in anti-phase with respect to the visible one. This has been explained with a photometric model that includes a mixture of cool and hot spots in the same active region.

Multi-band spectroscopy is known to be crucial in the search for exoplanets around young and active stars: here we present the first contribution of a new facility, GIARPS, which adds further value to this method thanks to the simultaneous observations in VIS and NIR bands. Our result clears the current census of hot Jupiters from the only previously known case orbiting a star between 20 to 200 Myr old (Rizzuto et al. 2017; David et al. 2018). The paucity of hot Jupiters in this age range might be explained by the time dependence of planet migration mechanisms (Mann et al. 20186) or simply by small number statistics. Further investigations are needed to achieve a firm conclusion on this key issue for the evolution of planetary systems.

Acknowledgements

The authors are grateful to the anonymous referee for the careful review and the useful comments that helped to improve the quality of the paper. The authors acknowledge support from INAF through the “WOW Premiale” funding scheme of the Italian Ministry of Education, University, and Research. Based on observations made with the Italian Telescopio Nazionale Galileo (TNG) operated on the island of La Palma by the Fundación Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the Spanish Observatorio del Roque de los Muchachos (ORM) of the IAC. This work used the Immersion Grating Infrared Spectrometer (IGRINS) that was developed under a collaboration between the University of Texas at Austin and the Korea Astronomy and Space Science Institute (KASI) with the financial support of the US National Science Foundation under grant AST-1229522, of the University of Texas at Austin, and of the Korean GMT Project of KASI. These results made use of the Discovery Channel Telescope (DCT) at Lowell Observatory. Lowell is a private, non-profit institution dedicated to astrophysical research and public appreciation of astronomy, and it operates the DCT in partnership with Boston Univ., the Univ. of Maryland, the Univ. of Toledo, Northern Arizona Univ., and Yale Univ. This paper includes data taken with the 2.7-m H. J. Smith Telescope at The McDonald Observatoryof The University of Texas at Austin. Based on observations made with the REM Telescope, INAF Chile.

Appendix A Additional tables

Table A.1

Time series of BD+20 1790 from GIANO and GIANO–B data.

Table A.2

Time series of BD+20 1790 from HARPS–N data, from HO15 and GIARPS commissioning, uniformly reduced with the new HARPS–N DRS version.

Table A.3

Time series of BD+20 1790 from IGRINS data.

Table A.4

Time series of BD+20 1790 from REM data with different filters.

References

  1. Baruteau, C., Crida, A., Paardekooper, S.-J., et al. 2014, Protostars and Planets VI, 667 [Google Scholar]
  2. Batygin, K., Bodenheimer, P. H., & Laughlin, G. P. 2016, ApJ, 829, 114 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bell, C. P. M., Mamajek, E. E., & Naylor, T. 2015, MNRAS, 454, 593 [NASA ADS] [CrossRef] [Google Scholar]
  4. Boisse, I., Bonfils, X., & Santos, N. C. 2012 A&A, 545, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Bonfils, X., Delfosse, X., Udry, S., et al. 2013, A&A, 549, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Carleo, I., Sanna, N., Gratton, R., et al. 2016, Exp. Astron., 41, 351 [NASA ADS] [CrossRef] [Google Scholar]
  7. Chapman, G. A., & McGuire, T. E. 1977, ApJ, 217, 657 [NASA ADS] [CrossRef] [Google Scholar]
  8. Chincarini, G., Zerbi, F., Antonelli, A., et al. 2003, The Messenger, 113, 40 [NASA ADS] [Google Scholar]
  9. Claret, A., Hauschildt, P. H., & Witte, S. 2012, A&A, 546, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Claudi, R., Benatti, S., Carleo, I., et al. 2017, EPJP, 132, 364 [NASA ADS] [Google Scholar]
  11. Conod, U., Blind, N., Wildi, F., et al. 2016, Proc. SPIE, 9909, 990941 [CrossRef] [Google Scholar]
  12. Cosentino, R., Lovis, C., Pepe, F., et al. 2014, Proc. SPIE, 9147, 91478 [Google Scholar]
  13. Crockett, C. J., Mahmud, N. I., Prato, L., et al. 2012, ApJ, 761, 164 [NASA ADS] [CrossRef] [Google Scholar]
  14. Damasso, M., Bonomo, A. S., Astudillo-Defru, N., et al. 2018, A&A, in press, DOI:10.1051/0004-6361/201732459 [Google Scholar]
  15. David, T. J., Mamajek, E. E, Vanderburg, A., et al. 2018, MNRAS 444 1157 [Google Scholar]
  16. Desidera, S., Carolo, E., Gratton, R., et al. 2011, A&A, 533, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Desort, M., Lagrange, A.-M., Galland, F., et al. 2007, A&A, 473, 983 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Donati, J. F., Moutou, C., Malo, L., et al. 2016, Nature, 534, 662 [NASA ADS] [CrossRef] [Google Scholar]
  19. Donati, J. F., Yu, L., Moutou, C., et al. 2017, MNRAS, 465, 3343 [NASA ADS] [CrossRef] [Google Scholar]
  20. Dumusque, X., Boisse, I., & Santos, N. C. 2014, ApJ, 796, 132 [NASA ADS] [CrossRef] [Google Scholar]
  21. Figueira, P., Pepe, F., Santos, N. C., et al. 2010a, EAS Pub. Ser., 42, 125 [CrossRef] [Google Scholar]
  22. Figueira, P., Marmier, M., Bonfils, X., et al. 2010b, A&A, 513, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Figueira, P., Santos, N. C., Pepe, F., et al. 2013 A&A, 557, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Gagné, J., Plavchan, P., Gao, P., et al. 2016, ApJ, 822, 402 [Google Scholar]
  25. González-Álvarez, E., Affer, L., Micela, G., et al. 2017, A&A, 606, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Gully-Santiago, M. A., Herczeg, G. J., Czekala, I., et al. 2017, ApJ, 836, 200 [NASA ADS] [CrossRef] [Google Scholar]
  27. Hatzes, A. P., Endl, M., Cochran, W. D., et al. 2018, AJ, 155, 120 [NASA ADS] [CrossRef] [Google Scholar]
  28. Haywood, R. D., Collier Cameron, A., Queloz, D., et al. 2014, MNRAS, 443, 2517 [NASA ADS] [CrossRef] [Google Scholar]
  29. Hernán-Obispo, M., Gálvez-Ortiz, M. C., et al. 2010, A&A, 512, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Hernán-Obispo, M., Tuomi, M., Gálvez-Ortizi et al. 2015, A&A, 576, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Herrero, E., Ribas, I., Jordi, C., et al. 2016, A&A, 586, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Hunter, A., Macgregor, A. B., Szabo, T. O., et al. 2012, Source Code Biol. Med., 7, 1 [CrossRef] [Google Scholar]
  33. Jeffries, R. D. 1995, MNRAS, 273, 559 [NASA ADS] [Google Scholar]
  34. Koen, C. 2016, MNRAS, 463, 4383 [NASA ADS] [CrossRef] [Google Scholar]
  35. Kozai, Y. 1962, ApJ, 67, 592 [Google Scholar]
  36. Lamm, M. H., Bailer-Jones, C. A. L., Mundt, R., et al. 2004, A&A, 417, 557 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Lanza, A. F. 2016, Lect. Notes Phys., (Berlin: Springer Verlag), 914, 43 [NASA ADS] [CrossRef] [Google Scholar]
  38. Lanza, A. F., Boisse, I., Bouchy, F., et al. 2011, A&A, 533, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Lanza A. F., Malavolta L., Benatti S., et al. 2018, A&A in press, [arXiv:1804.07039] [Google Scholar]
  40. López-Santiago, J. 2005, Ph.D. Thesis, Universidad Complutense de Madrid [Google Scholar]
  41. Lovis, C., Dumusque, X., Santos, N. C., et al. 2011, ArXiv eprints [arXiv:1107.5325] [Google Scholar]
  42. Mace, G., Kim, H., Jaffe, D. T., et al. 2016, Proc. SPIE, 9908, 99080C [Google Scholar]
  43. Mahmud, N. I., Crockett, C. J., Johns-Krull, C. M., et al. 2011, ApJ, 736, 123 [NASA ADS] [CrossRef] [Google Scholar]
  44. Najita, J., Carr, J. S., & Tokunaga, A. T. 1996, ApJ, 456, 292 [NASA ADS] [CrossRef] [Google Scholar]
  45. Nardetto, N., Mourard, D., Kervella, P., et al. 2006, A&A, 453, 309 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Oliva, E., Origlia, L., Baffa, C., et al. 2006, Proc. SPIE, 6269, 626919 [CrossRef] [Google Scholar]
  47. Park, C., Jaffe, D. T., Yuk, I., et al. 2014, Proc. SPIE, 9147, 91471 [Google Scholar]
  48. Penn, M. J. 2014, Liv. Rev. Sol. Phys., 11, 2 [Google Scholar]
  49. Pepe, F., Mayor, M., Galland, F., et al. 2002, A&A, 388, 632 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Prato, L., Huerta, M., Johns-Krull, C. M., et al. 2008, ApJ, 687, L103 [NASA ADS] [CrossRef] [Google Scholar]
  51. Queloz, D., Henry, G. W., Sivan, J. P., et al. 2001, A&A, 379, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Quirrenbach, A., Amado, P. J., et al. 2014, Proc. SPIE, 9147, 91471F [Google Scholar]
  53. Reiners, A., Bean, J. L., Huber, K. F., et al. 2010, ApJ, 710, 432 [NASA ADS] [CrossRef] [Google Scholar]
  54. Reiners, A., Shulyak, D., Anglada-Escudé, G., et al. 2013, A&A, 552, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Reiners, A., Zechmeister, M. Caballero, J. A., et al. 2018, A&A, 612, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Rice, J. B., Strassmeier, K. G., & Kopf, M., 2011, ApJ, 728, 69 [NASA ADS] [CrossRef] [Google Scholar]
  57. Rizzuto, A. C., Mann, A. W., Vanderburg, A., et al. 2017, AJ, 154, 224 [NASA ADS] [CrossRef] [Google Scholar]
  58. Roberts, D. H., Lehar, J., & Dreher, J. W. 1987, AJ, 93, 968 [NASA ADS] [CrossRef] [Google Scholar]
  59. Robertson, P., Bender, C., Mahadevan, S., et al. 2016, ApJ, 832, 112 [NASA ADS] [CrossRef] [Google Scholar]
  60. Saar, S. H., & Donahue, R. A. 1997, ApJ, 485, 319 [NASA ADS] [CrossRef] [Google Scholar]
  61. Sánchez Cuberes, M., Vázquez, M., Bonet, J. A., et al. 2002, ApJ, 570, 886 [NASA ADS] [CrossRef] [Google Scholar]
  62. Scargle, J. D. 1982, ApJ, 263, 835 [NASA ADS] [CrossRef] [Google Scholar]
  63. Seemann U., Anglada – Escudé, G., Baade, D., et al. 2014, Proc. SPIE Conf. Ser., 9147, 91475 [Google Scholar]
  64. Setiawan, J., Henning, T., Launhardt, R., et al. 2008, Nature, 451, 38 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  65. Sissa, E., Gratton, R., Desidera, A. F., et al. 2016, A&A, 596, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Torres, C. A. O., Quast, G. R., Melo, C. H. F., & Sterik, M. F. 2008, in Handbook of Star Forming Regions, ed. Bo Reipurth, Vol. II: The Southern Sky ASP Monograph Publications, 5, 757 [NASA ADS] [Google Scholar]
  67. Tozzi, A., Oliva, E., Iuzzolino, M., et al. 2016, Proc. SPIE, 9908, 238 [Google Scholar]
  68. Yu, L., Donati, J.-F., Hébrard, E. M., et al. 2017 MNRAS, 467, 1342 [NASA ADS] [Google Scholar]
  69. Yuk, I., Jaffe, D. T.; Barnes, S., et al. 2010, Proc. SPIE, 7735, 77351 [CrossRef] [Google Scholar]
  70. Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

4

We adopted the IDL least-squares MPFIT package available at http://purl.com/net/mpfit

5

See a description of the Vasy(mod) in the poster “Line asymmetry indicators to detect stellar activity effects in radial velocity measurements" by Lanza et al., available at: https://sites.google.com/a/yale.edu/eprv-posters/home.

All Tables

Table 1

Summary of the spectroscopic data presented in this work.

Table 2

Summary of the values of Spearman (ρ) correlation coefficients and corresponding p-values between RVs and activity indicators.

Table A.1

Time series of BD+20 1790 from GIANO and GIANO–B data.

Table A.2

Time series of BD+20 1790 from HARPS–N data, from HO15 and GIARPS commissioning, uniformly reduced with the new HARPS–N DRS version.

Table A.3

Time series of BD+20 1790 from IGRINS data.

Table A.4

Time series of BD+20 1790 from REM data with different filters.

All Figures

thumbnail Fig. 1

Orbital fit at 7.78 days found by HO15 compared to phase-folded visible and NIR RVs. Top panel: orbital fit (black dashed line)obtained with the visible data (FOCES, SARG, and HERMES RVs from HO15, grey dots) and HARPS-N 2015 RVs (green dots). Bottom panel: orbital fit (black dashed line), GIANO/GIANO-B (red dots), IGRINS (light blue dots), and HARPS-N 2017 (black asterisks, two acquired in GIARPS mode) RVs.

In the text
thumbnail Fig. 2

Correlation between RVs and BIS for HARPS–N data (HO15 and the three new epochs in this work).

In the text
thumbnail Fig. 3

Phase-folded HARPS–N RVs (2015, reprocessed from HO15 dataset) at stellar rotational period.

In the text
thumbnail Fig. 4

Top panel: generalized Lomb-Scargle periodogram of the g-band photometric time series of BD+20 1790. The solid black line isthe normalized power versus period, the dotted red line is the spectral window, and the horizontal dashed line indicates the power level corresponding to a FAP = 0.01. Bottom panel: CLEAN periodogram. The red mark is the most powerful and significant peak in the periodogram.

In the text
thumbnail Fig. 5

Differential lightcurves of BD+20 1790 phased with the rotation period P = 2.76 d in different photometric bands. The solid line is a sinusoidal fit to the data with the same period. Labels show the peak-to-peak amplitudes of the lightcurves.

In the text
thumbnail Fig. 6

Contrast of an active region consisting of co-spatial cool and hot spots versus the fraction of its area covered by the cool spot. Two wavelengths are considered: λopt = 636 nm (solid line) and λNIR = 1705 nm (dashed line).

In the text
thumbnail Fig. 7

Light curves showing the rotational modulation of the flux at two wavelengths λopt = 636 nm (solid line) and λNIR = 1705 nm (dashed line). The amplitudes correspond to a spot having an area of 0.2 of the stellar disc.

In the text

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.