Multi-band high resolution spectroscopy rules out thehot Jupiter BD+20 1790bBased 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 Observatory of The University of Texas at Austin.{}^{,}Based on observations made with the REM Telescope, INAF Chile.

Multi-band high resolution spectroscopy rules out the hot Jupiter BD+20 1790b - First data from the GIARPS Commissioning

First data from the GIARPS Commissioning
Key Words.:
Instrumentation: spectrographs – exoplanets – Spectroscopy – Active stars – Radial velocities – Stars: individual: BD+20 1790


Context:Stellar activity is currently challenging the detection of young planets via the radial velocity (RV) technique.

Aims:We attempt to definitively discriminate the nature of the RV variations for the young active K5 star BD+20 1790, for which visible (VIS) RV measurements show divergent results on the existence of a substellar companion.

Methods:We compare VIS data with high precision RVs in the near infrared (NIR) range by using the GIANO – B and IGRINS spectrographs. In addition, we present for the first time simultaneous VIS-NIR observations obtained with GIARPS (GIANO – B and HARPS – N) at Telescopio Nazionale Galileo (TNG). Orbital RVs are achromatic, so the RV amplitude does not change at different wavelengths, while stellar activity induces wavelength-dependent RV variations, which are significantly reduced in the NIR range with respect to the VIS.

Results:The NIR radial velocity measurements from GIANO – B and IGRINS show an average amplitude of about one quarter with respect to previously published VIS data, as expected when the RV jitter is due to stellar activity. Coeval multi-band photometry surprisingly shows larger amplitudes in the NIR range, explainable with a mixture of cool and hot spots in the same active region.

Conclusions:In this work, the claimed massive planet around BD+20 1790 is ruled out by our data. We exploited the crucial role of multi-wavelength spectroscopy when observing young active stars: thanks to facilities like GIARPS that provide simultaneous observations, this method can reach its maximum potential.

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. 2017) or HARPS+NIRPS (0.95–1.8 nm, R100,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 emisphere) 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 (, 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 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 ( m) with a resolution of R . In the same year, the NIR high resolution spectrograph GIANO (Oliva et al., 2006), working in the wavelength range from to m at a resolution of R , 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 ( - m with a small gap between and 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 of m  s) and a Cross-Correlation Function (CCF) method based on telluric lines with GIANO – B ( m  s, 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) 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.

IGRINS (Immersion Grating Infrared Spectrometer) (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= (Yuk et al., 2010; Park et al., 2014). It covers the H and K windows, from 1.45 to 2.5 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, and ROS2, 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 (20482048 pixels, 13.5 micron pixels size, 0.58 arcsecpixel plate scale), which has a corrected 9.89.8 arcmin field of view, and observes simultaneously through the four SloanSDSS g,r,i,z filters. The REMIR camera is equipped with a Hawaii I CCD (512512 pixels, 1.2 arcsecpixel plate scale), which has a corrected 1010 arcmin field of view, and it is equipped with zJHK filters.

3 Observations and data reduction

In this section we describe the data acquired from each instrument. The spectroscopic datasets are presented in Sects. 3.1, 3.2, and 3.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 r.m.s., and the peak-to-valley amplitude, is reported in Table 1. The photometric observations are described in Sec. 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).

Instrument N Spect. range SNR RV r.m.s. Peak-to-valley
(kms) (kms) (kms)
GIANO/GIANO – B 18 0.95 - 2.45 72 0.036 0.130 0.384
HARPS – N 20 0.38 - 0.69 35 0.029 0.280 1.036
IGRINS 29 1.45 - 2.5 134 0.059 0.109 0.384
Table 1: Summary of the spectroscopic data presented in this work. For each dataset we list the instrument used for the observations, the number of spectra, the spectral range, the typical SNR, the RV nominal internal error (), the RV r.m.s. scatter, and the peak-to-valley value of the RVs.

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_TOOLS6 package, while GIANO – B spectra are processed with the dedicated pipeline GOFIO (Rainer et al, SPIE, in prep). We obtained a set of 30 NIR RVs over one year (listed in Table LABEL:tab:bd201790_giano), with the method described in Carleo et al. (2016), in which the telluric lines are 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 , where 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 signal-to-noise ratio at JD 2456681) with the current HARPS – N DRS through the YABI workflow (Hunter et al., 2012) installed at IA27 at the INAF Observatory of 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 of this star (10.030.47 km  s, López-Santiago 2005) and the consequent line broadening. The resulting RVs are listed in Table LABEL:tab:bd201790_harpsn, together with the CaII activity indicator, R’, 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. (submitted), 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. (accepted)).

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 Package.8 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 LABEL:tab:bd201790_igrins) 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 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 series of six average griz-band differential magnitudes and four JHK-band differential magnitudes for the subsequent analysis. The average standard 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 for the 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 function9 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 and the period is 7.7827 10 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 r.m.s. scatter of 129.6 ms, while IGRINS RVs (light blue dots) show an r.m.s. scatter of 109.2 ms. The r.m.s. scatter of the whole NIR dataset is 119 ms. 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, 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.

Figure 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 -values between RVs and activity indicators.

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 = , 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 r.m.s. of 61.1 ms 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.801502 d, K=0.4497.370 kms, e=0.0650.018).

We also investigated the other activity indices available from HARPS – N, for example 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 V (the modified version of the index V by Figueira et al. 2013 for which the dependence from the RV shift is removed), presented in Lanza et al. (submitted).10 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 and H indicators (Table 2).

Figure 2: Correlation between RVs and BIS for HARPS – N data (HO15 and the three new epochs in this work).
Figure 3: Phase-folded HARPS – N RVs (2015, reprocessed from HO15 dataset) at stellar rotational period.

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 at 7.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).

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 (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 (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.

4.2 Photometric data

We used the GLS 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 1,000 moke light curves as in Sect. 4.1).

Figure 4: Top panel: Generalized Lomb-Scargle periodogram of the g-band photometric time series of BD+20 1790. The solid black line is the 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.

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 different bands: unexpectedly, the amplitudes at longer wavelengths are greater than those at shorter wavelengths and the NIR modulation is almost in anti-phase with respect to the optical modulation in the passbands (the maximum of the light curves are located at for the filters, at for and at 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.

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 BD20 1790 is generally a factor of 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 ( km s) 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 , radius, and rotation period (from HO10), we estimated an inclination of the stellar spin axis to the line of sight of . 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 of its total area and a hot feature (hereinafter ”hot spot”) covering the remaining fraction .
Therefore, the average brightness of the active region at wavelength can be written as


where, for simplicity, we assume that the brightness of each component is given by a Planck function with being the temperature of the cool spot and that of the hot spot. Those temperatures verify the inequality: , where is the temperature of the unperturbed photosphere. The contrast of the active region is , where is the brightness of the unperturbed photosphere. For an active region dominated by the cool spot, because , while for an active region dominated by the hot spot, .

For BD+20 1790, we assume  K (HO15), while for the cool and the hot spots we assume temperatures  K and  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  nm and the NIR wavelength  nm. Those values correspond approximately to the mean wavelengths of the bands and of the bands, respectively.

We see that for a restricted range of , that is , 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 , considering an active region of an area of 20% of that of the star’s disc, centred at a latitude of . The quadratic limb-darkening coefficients at the two wavelengths were taken from Claret et al. (2012), while and corresponding to 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.

Figure 5: Differential lightcurves of BD+20 1790 phased with the rotation period P = 2.76d 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.

Figure 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:  nm (solid line) and  nm (dashed line).

Figure 7: Light curves showing the rotational modulation of the flux at two wavelengths  nm (solid line) and  nm (dashed line). The amplitudes correspond to a spot having an area of 0.2 of the stellar disc.

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 r.m.s. 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 at latitude with a filling factor , we expect an RV modulation approximately of (cf. Saar et al. 1997; Desort et al. 2007) that is km s for a spot at , 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).

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, significantly lower than the VIS one (1853.6 ms), 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. 2018 11) 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.

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.


  1. thanks: 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.
  2. thanks: 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.
  3. thanks: 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.
  4. thanks: This paper includes data taken with the 2.7-m H. J. Smith Telescope at The McDonald Observatory of The University of Texas at Austin.
  5. thanks: Based on observations made with the REM Telescope, INAF Chile.
  6. Available at the TNG webpage:
  9. We adopted the IDL least-squares MPFIT package available at
  10. See a description of the V in the poster “Line asymmetry indicators to detect stellar activity effects in radial velocity measurements” by Lanza et al., available at:


  1. Baruteau, C., Crida, A., Paardekooper, S.-J., et al. 2014, Protostars and Planets VI, 667
  2. Batygin, K., Bodenheimer, P. H., Laughlin, G. P. 2016, ApJ, 829, 114
  3. Bell, C. P. M., Mamajek, E. E., Naylor, T. 2015, MNRAS, 454, 593
  4. Boisse, I., Bonfils, X. & Santos, N. C. 2012 A&A, 545, A109
  5. Bonfils, X., Delfosse, X., Udry, S., et al. 2013, A&A, 549, 109
  6. Carleo, I., Sanna, N., Gratton, R., et al. 2016, Experimental Astronomy, 41, 351
  7. Chapman, G. A. & McGuire, T. E. 1977, ApJ, 217, 657
  8. Chincarini, G., Zerbi, F., Antonelli, A. et al. 2003, The Messenger, 113, 40
  9. Claret, A., Hauschildt, P. H., & Witte, S. 2012, A&A, 546, A14
  10. Claudi, R., Benatti, S., Carleo, I., et al. 2017, EPJP, 132, 364
  11. Conod, U., Blind, N., Wildi, F., et al. 2016, Proc. SPIE, 9909, 990941
  12. Cosentino, R., Lovis, C., Pepe, F.,et al. 2014, Proc. SPIE, 9147, 91478
  13. Crockett, C. J., Mahmud, N. I., Prato, L., et al. 2012, ApJ, 761, 164
  14. Damasso, M., Bonomo, A. S., Astudillo-Defru, N., et al. 2018, ArXiv e-prints [arXiv:1802.08320]
  15. David, T. J., Mamajek, E. E, Vanderburg, A. et al. 2018, ArXiv e-prints [arXiv:1801.07320]
  16. Desidera, S., Carolo, E., Gratton, R., et al. 2011, A&A, 533, 90
  17. Desort, M., Lagrange, A.-M., Galland, F. et al. 2007, A&A, 473, 983
  18. Donati, J. F., Moutou, C., Malo, L. et al. 2016, Nature, 534, 662
  19. Donati, J. F., Yu, L., Moutou, C. et al. 2017, MNRAS, 465, 3343
  20. Dumusque, X., Boisse, I. and Santos, N. C. 2014, ApJ, 796, 132
  21. Figueira, P., Pepe, F., Santos, N. C., et al. 2010a, EAS Publications Series, 42, 125
  22. Figueira, P., Marmier, M., Bonfils, X., et al. 2010b, A&A, 513, L8
  23. Figueira, P., Santos, N. C., Pepe, F., et al. 2013 A&A, 557, A93
  24. Gagné, J., Plavchan, P., Gao, P., et al. 2016, ApJ, 822, 40 2
  25. González-Álvarez, E., Affer, L., Micela, G., et al. 2017, A&A, 606, A51
  26. Gully-Santiago, M. A., Herczeg, G. J., Czekala, I., et al. 2017, ApJ, 836, 200
  27. Hatzes, A. P., Endl, M., Cochran, W. D., et al. 2018, AJ, 155, 120
  28. Haywood, R. D., Collier Cameron, A., Queloz, D., et al. 2014, MNRAS, 443, 2517
  29. Hernán-Obispo, M., Gálvez-Ortiz, M. C., et al. 2010, A&A, 512, A45
  30. Hernán-Obispo, M., Tuomi, M., et al. 2015, A&A, 576, A66
  31. Herrero, E., Ribas, I., Jordi, C. et al. 2016, A&A, 586, A131
  32. Hunter, A., Macgregor, A.B., Szabo, T.O., et al. 2012, Source Code Biol. Med., 2012, vol. 7, 1
  33. Jeffries, R. D. 1995, MNRAS, 273, 559
  34. Koen, C. 2016, MNRAS, 463, 4383
  35. Kozai, Y. 1962, ApJ, 67, 592
  36. Lamm, M. H., Bailer-Jones, C. A. L. et al. 2004, A&A, 417, 557
  37. Lanza, A. F., Boisse, I., Bouchy, F. et al. 2011, A&A, 533, A44
  38. Lanza, A. F. 2016, Lecture Notes in Physics, Berlin Springer Verlag, 914, 43
  39. López-Santiago, J. 2005, PhD Thesis, Universidad Complutense de Madrid
  40. Lovis, C., Dumusque, X., Santos, N. C., et al. 2011, arXiv:1107.5325
  41. Mace, G., Kim, H., Jaffe, D. T., et al. 2016, Proc. SPIE, 9908, 99080C
  42. Mahmud, N. I., Crockett, C. J., Johns-Krull, C. M. et al., 2011, ApJ, 736, 123
  43. Najita, J., Carr, J. S. and Tokunaga, A. T. 1996, ApJ, 456, 292
  44. Nardetto, N., Mourard, D., Kervella, P., et al. 2006, A&A, 453, 309
  45. Oliva, E., Origlia, L., Baffa, C., et al. 2006, Proc. SPIE, 6269, 626919
  46. Park, C., Jaffe, D. T., Yuk, I. et al. 2014, Proc. SPIE, 9147, 91471
  47. Penn, M. J. 2014, Living Reviews in Solar Physics, 11, 2
  48. Pepe, F., Mayor, M., Galland, F. et al. 2002, A&A, 388, 632
  49. Prato, L., Huerta, M., Johns-Krull, C. M. et al, 2008, ApJ, 687, L103
  50. Queloz, D., Henry, G. W., Sivan, J. P. et al, 2001, A&A, 379, 279
  51. Quirrenbach, A., Amado, P. J. et al. 2014, Proc. SPIE, 9147, 91471F
  52. Reiners, A., Bean, J. L., Huber, K. F., et al. 2010, ApJ, 710, 432
  53. Reiners, A., Shulyak, D., Anglada-Escudé, G. et al. 2013, A&A, 552, A103
  54. Reiners, A., Zechmeister, M. Caballero, J. A., et al. 2017, arXiv171106576R
  55. Rice, J. B., Strassmeier, K. G., & Kopf, M., 2011, ApJ, 728, 69
  56. Rizzuto, A. C., Mann, A. W., Vanderburg, A. et al., 2017, AJ, 154,224
  57. Roberts, D. H., Lehar, J. and Dreher, J. W. 1987, AJ, 93, 968
  58. Robertson, P., Bender, C., Mahadevan, S., et al. 2016, ApJ, 832, 112
  59. Saar, S. H. and Donahue, R. A. 1997, ApJ, 485,319
  60. Sánchez Cuberes, M., Vázquez, M., Bonet, J. A., et al. 2002, ApJ, 570, 886
  61. Scargle, J. D., 1982, ApJ, 263, 835
  62. Seemann U., Anglada – Escudé G., Baade D. et al. 2014, Proc. SPIEConf. Ser., 9147, 91475
  63. Setiawan, J., Henning, T., Launhardt, R. Müller, A., Weise, P., Kürster, M., 2008, Nature, 451, 38
  64. Sissa, E., Gratton, R., Desidera, et al. 2016, A&A, 596A, 76S
  65. Torres, C. A. O., Quast, G. R., Melo, C. H. F., et al. 2008, in Handbook of Star Forming Regions, Volume II: The Southern Sky ASP Monograph Publications, Vol. 5., Ed. Bo Reipurth, p.757
  66. Tozzi, A., Oliva, E., Iuzzolino, M., et al. 2016, Proc. SPIE, 9908, 238
  67. Yu, L., Donati, J.-F., Hébrard, E. M. et al. 2017 MNRAS, 467, 1342
  68. Yuk, I., Jaffe, D. T.; Barnes, S., et al. 2010, Proc. SPIE, 7735, 77351
  69. Zechmeister, M. and Kürster, M. 2009, A&A, 496, 577
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description