Multi-Wavelength High-Resolution Spectroscopy for Exoplanet Detection: Motivation, Instrumentation and First Results

Exoplanet research has shown an incessant growth since the first claim of a hot giant planet around a solar-like star in the mid-1990s. Today, the new facilities are working to spot the first habitable rocky planets around low-mass stars as a forerunner for the detection of the long-awaited Sun-Earth analog system. All~the achievements in this field would not have been possible without the constant development of the technology and of new methods to detect more and more challenging planets. After the consolidation of a top-level instrumentation for high-resolution spectroscopy in the visible wavelength range, a huge effort is now dedicated to reaching the same precision and accuracy in the near-infrared. Actually, observations in this range present several advantages in the search for exoplanets around M dwarfs, known to be the most favorable targets to detect possible habitable planets. They are also characterized by intense stellar activity, which hampers planet detection, but~its impact on the radial velocity modulation is mitigated in the infrared. Simultaneous observations in the visible and near-infrared ranges appear to be an even more powerful technique since they provide combined and complementary information, also useful for many other exoplanetary science~cases.


Introduction
The radial velocity (or Doppler) method allowed detecting the first extrasolar planet around a solar-like star, 51 Peg b [1], opening the era of the quest for planets. The radial velocity (RV) technique is an "indirect" method to detect exoplanets since it reveals the reflex motion of the host star due to a hidden companion when both revolve around the common center of mass. The periodical variation of the stellar RV, obtained through the measurement of the Doppler shift of its spectral lines, is dependent on the characteristic of its small-mass companion. The modulation can be represented by a Keplerian function with a period equal to the orbital period of the planet, P, and a semi-amplitude, K, defined as follows: where M is the stellar mass, e the orbital eccentricity of the planet and m p sin i is the lower limit of the planetary mass since the real mass depends on the orbital inclination, i, unknown with the RV information alone. The actual value of the planet mass can be obtained if the planet transits in front of its host star and a proper modeling of the observed transit light curve is performed. Just like the shape of the RV variation allows measuring a set of orbital parameters with Equation (1), the same principle is valid for the transit light modulation, which can be modeled by specific functions presented in [2]. The shape of the light curve is dependent, among other things, on the inclination of the planetary orbital plane through the so-called "impact parameter", b (i.e., the projected separation between the center of the stellar disk and the center of the planet disk), according to the relation b = a cos i R , where a is the semi-major axis of the orbit and R is the stellar radius (see, e.g., [3]). Once the inclination is transits survey performed by the NASA-Kepler satellite [6,7], the RV was the most productive method for exoplanets' detection (e.g., [8]). The technology developed for the optical instrumentation showed an increasing improvement through the years: the fiber-fed ELODIE (installed at the 1.93-m telescope of Observatoire de Haute-Provence -OHP, France) spectrograph, which was used to find 51 Peg b in 1995, already ensured a good RV accuracy at that time (∼13 m s −1 ), by using the simultaneous thorium reference technique [9]. With this method, a scientific fiber feeds the instrument with the light of the target, while a second fiber is illuminated by a thorium-argon (ThAr) lamp, aiming to obtain a simultaneous wavelength calibration able to track and correct any instrumental drifts occurring during the observing night. Other spectrographs, such as HIRES(High-Resolution Echelle Spectrometer, mounted at the 10-m Keck I telescope, USA), adopted the Iodine cell technique [10], in which a gas absorption cell is inserted in the optical path of the instrument, obtaining a superimposition of the I 2 lines on the scientific spectra, producing thus a stable wavelength reference. This technique allows one to reach an RV precision up to 3 m s −1 (e.g., [11]). A significant step forward in the accuracy of the RV measurement in the optical band has been achieved with HARPS (High Accuracy Radial velocity Planet Searcher) [12], a high-resolution fiber-fed echelle spectrograph mounted at the ESO 3.6-m Telescope in La Silla (Chile). This spectrograph operates in a vacuum enclosure in controlled temperature and pressure environment, located in an isolated room to avoid vibrations connected to the telescope building. In order to obtain even more precise RV measurements, new calibration systems were studied and routinely adopted for HARPS and HARPS-N (its twin in the Northern Hemisphere; see Section 3.2), such as the Fabry-Perot (FP) interferometers [13]. As for the ThAr method, a calibration source (e.g., a lamp or laser) overimposes well-spaced emission lines on the detector, close to the stellar spectrum. As an example, Figure 1 displays a 2D spectrum of HARPS showing the traces of the FP compared to the less rich and unevenly-spaced lines of the ThAr lamp [13], demonstrating the high potential of this type of calibration. A similar, but more sophisticated and expensive technique, the laser frequency comb has been developed and tested, as well [14,15]. This method still needs a better optimization before its full exploitation, but thanks to the excellent calibration performances, it should contribute to providing an RV precision significantly lower than 1 m s −1 . Thanks to the state-of-the-art calibration techniques and to an optimized RV extraction [12], HARPS has been able to reach an RV accuracy of 1 m s −1 or less (depending on the SNR of the spectra and the stellar properties) and gave a significant contribution to exoplanet research in the last decade [16][17][18]. HARPS and HARPS-N have represented the state-of-the-art of the RV technique for a long time: now, at the beginning of the scientific operations of ESPRESSO (Echelle SPectrograph for Rocky Exoplanet and Stable Spectroscopic Observations, [19], https://obswww.unige.ch/Instruments/ espresso/) at the Very Large Telescope (VLT, in Chile), the incredible regime of 10 cm s −1 is waiting to be covered. The interested reader is referred to [20] for a detailed overview of the RV surveys of the last decades, as well as a resume of the status of the RV method.
With the previous discussion in mind, it should be emphasized that the planet distribution obtained by the RV surveys [8,21,22] suffers from a sample selection bias. Indeed, for years, the optical surveys have been mostly focused on "quiet" stars, i.e., stars with low or negligible stellar activity. Several processes can be ascribed to the stellar activity, such as the pulsations (timescale of a few tens of minutes for solar-type stars), granulation (hours) and magnetic activity (up to a few tens of days). All of them are responsible for the intrinsic RV variation of a star, together with all the other sources of stellar noise with no periodic patterns, called "jitter". In particular, magnetic activity induces a variable number of evolving photospheric inhomogeneities (e.g., cool spots and bright plage), the visibility of which is modulated by the stellar rotation. These phenomena produce subtle time-dependent changes in the line profiles of the star, resulting in an RV modulation that can be erroneously interpreted as a Keplerian motion [18,[23][24][25]. For this reason, stellar activity represents a real problem in the search for exoplanets, and several studies have been proposed to identify useful diagnostics helping in discriminating the real nature of an RV variation. First of all, a line deformation due to activity provokes an asymmetry of the line profile that can be observed through the line bisector, as shown, e.g., in [26]. In order to quantify the extent of the asymmetry, it is useful to evaluate the difference between average values of velocities in regions at the top and at the bottom parts of the line profile (see Figure 5 in [27]), which is defined as the bisector velocity span (BVS [28]). When photospheric inhomogeneities are responsible for the RV variation, the associated line profile variation produces a correlation between the RV and bisector span [27]. Besides the bisector, other activity indicators have been proposed, aiming to find false-positives in the RV time series: some of them take into account the asymmetry of the line profiles [29,30] or measure the emission of the H and K lines of the Ca II [31,32], centered at 3968.47 and 3933.66 Å, respectively. Unfortunately, the lack of correlation does not imply the presence of a planet [33], and the correlation is not always linear or easily detectable [34]. The treatment of the stellar activity can be very challenging, especially when one looks for low-mass planets having a small RV semi-amplitude comparable to or lower than the typical value of the RV variation due to the stellar activity. Only in recent years have specific tools been presented and heavily used to disentangle the planetary signal from the activity one; for instance, the Gaussian processes [35]. This is the reason why several stellar classes have been generally not considered in the past for the planet searching surveys, with a consequent lack of information about the planet frequency around all those targets showing high levels of stellar activity (e.g., young and intermediate-age stars, M dwarfs, etc.). Today, the treatment of stellar activity requires a substantial number of data, an optimized observing strategy with intense monitoring and a good degree of precision of the measurements. In principle, these conditions can be reached in the visible band (e.g., [36,37]), as discussed above, while a less favorable situation takes place in the case of the NIR facilities.

The Role of the NIR
A drastic turnaround occurred when the planet occurrence rate was assessed for low-mass stars in the Kepler sample, according to which the formation of small planets is highly favorable with respect to FGK stars [38,39], with an estimate of 2.5 ± 0.2 planets per M dwarfs with radii 1-4 R ⊕ and periods shorter than 200 days [40]; see Figure 2, left panel. Since the RV semi-amplitude K (Equation (1)) scales inversely with the stellar mass and the planet orbital period, planets with similar properties induce a larger RV amplitude on an M dwarf with respect to more massive stars.  [40]. Right panel: Planet occurrence rate vs. insolation. The green box indicates the typical range of insolation of the habitable zone. The green dashed lines indicate the limits for the maximum greenhouse and moist greenhouse insolation, while the magenta dotted-dashed lines represent the less conservative recent Venus and early Mars limits. Credit [40] ©AAS. Reproduced with permission.
The lower luminosity of M dwarfs also implies that their habitable zone (HZ), i.e., the circumstellar region where water is allowed to exist in liquid state on a planet surface, lies at shorter orbital distances [41]. As a result, even planets with relatively short orbital periods could be located inside the HZ of an M dwarf [42]. Their occurrence rate has been estimated as 0.16 +0. 17 −0.07 for Earth-sized planets (0.5-1.4 R ⊕ ) and 0.12 +0.1 −0.05 for Super Earths (1.4-4 R ⊕ ) in the case of the moist greenhouse inner and maximum greenhouse outer limits; 0.24 +0. 18 −0.08 for Earth-sized planets and 0.21 +0.11 −0.06 for Super Earths for the recent Venus and early Mars limits, as defined in [41]; see the right panel of Figure 2. Thanks to their characteristics and their ubiquity in the Universe [43], M dwarfs appeared as the ideal shortcut to find the first habitable Earth-like planet. As a downside, M dwarfs are extremely faint in the VIS wavelength and show increasing levels of stellar activity from earlier to later spectral types. This is where the NIR wavelength band comes to the rescue [44]: the M dwarfs' spectral energy distribution peaks between 1 and 2.5 µm, and the observation in the NIR range provides important advantages.
The main benefit when observing in the NIR band is that the stellar RV variation due to activity is expected to be lower than in the VIS [44,45], because of the lower intensity contrast between the stellar surface and the spots [46,47]. In the framework of the search for exoplanets, this represents an outstanding advantage: while the activity signal is wavelength dependent, the RV modulation due to a planetary companion must show the same amplitude regardless of the band. As a matter of fact, this allows discriminating between signals of different origin.
The search for habitable rocky planets around M dwarfs is not the only science case that benefits from observation in the NIR. The extension of the high-precision radial velocity measurements to this spectral range allows one to search for planetary companions to all those objects that were usually not considered due to the intrinsic difficulties, like stars with a moderate or high level of stellar activity, intermediate-age stars in open clusters or stars in very young associations (a few Myrs). In the latter case, exoplanets could be observed when their formation was ongoing, and the determination of their orbital parameters allows one to put constraints on the planet migration mechanisms [48][49][50]. According to the recent results [51,52] the frequency of planets around young stars is higher with respect to the older ones. A systematic investigation of these objects and the verification of this finding could help to understand the evolution of the systems. Comparison between the planetary mass distribution also allows inferring the photo-evaporation of the planet's atmosphere [53,54]. With the support of the NIR band, it is also possible to search for the presence of exoplanets around red giant stars. The typical pulsation periods of red giants range between two regimes: a short-time periodicity driven by pressure modes (with a timescale of the order of a few hours up to days, e.g., [55]) and a long-time periodicity (hundreds of days) probably due to oscillatory convective modes [56] (see [57] for a description of stellar oscillations in red giants). Obviously, this could hamper the detection of planets with similar orbital periods [58], but the oscillations being chromatic phenomena [59], it is possible to distinguish them with combined VIS-NIR observations.
Finally, alongside the RV search for exoplanets, high-resolution NIR spectroscopy also allows the atmospheric characterization of close-in giant planets, representing nowadays the training ground for the future detection of biosignatures of exo-Earths [60]. Actually, the main molecular constituents of a planet atmosphere show their absorption bands in the near-and mid-infrared (water, oxygen, methane, carbon monoxide, etc.; see, e.g., [61]) and can be detected through the transmission and emission spectroscopy method. In the former, the light of the host star is filtered by the atmosphere of a giant gas planet during the transit phase, impressing its extremely weak features on the stellar spectrum [62][63][64]. In the latter, the thermal emission from the planet, coming from internal heat and/or the reflected starlight, can be extrapolated from the stellar spectra by following part of the orbital phase, even for non-transiting exoplanets [65][66][67]. This review is focused on the planet search with the RV method, so the atmospheric characterization will be marginally discussed; anyway, the interested reader can refer to [68,69] for a detailed introduction to this argument.
It is now clear that the science achievable with the extension of the high-resolution spectroscopy to the NIR domain is really noteworthy. Unfortunately, all the previous efforts to develop techniques and build high-precision instruments working efficiently in the optical range could not be easily applied in another band. In the NIR range, the thermal noise represents an annoying source of disturbance in the spectra, especially toward the redder parts. For this reason, the instruments must be cooled down, which implies different technology and more costs. However, the main challenge to obtain high-precision RV measurements in the NIR is to define a reliable procedure for an accurate wavelength calibration. The simultaneous ThAr technique, so efficient in the VIS, is not suitable because in the NIR range, thorium shows a small number of faint spectral lines, which are generally contaminated by the scattered light of the argon emission lines, sensitive to changes in the environmental conditions (see [70] and the references therein). The calibration with a laser comb or FP is studied for NIR, as well. Both the hardware and a corresponding optimized technique were not available up to a few years ago, but after the promising results of dedicated studies [71,72], more efforts were performed and will be/are adopted in some of the current and future NIR facilities [73,74] (see also Sections 3.2 and 3.3).
Many studies have been carried out about the optimization of the NIR gas absorption cell technique [70,[75][76][77], similar to the I 2 cell in the visible. This approach has been demonstrated to be quite effective (RV precision, from ∼20 m s −1 [72] up to 3-5 m s −1 [76]), even if the emission lines of a selected gas are usually available for specific wavelength bands only (e.g., J or H) and do not cover the whole spectral range typically offered by a spectrograph. In the last few years, dedicated studies have been performed aiming to define suitable gas mixtures to be inserted in the cells, able to span larger wavelength ranges [78]. Finally, another approach to obtain a reliable wavelength reference for RV measurement is to exploit the rich forest of telluric lines that dominate the NIR spectra [79] (see, e.g., Figure 3). These lines are due to the light absorption of molecular species of the Earth's atmosphere (mainly OH, H 2 O, CO 2 ) that impress their features on the spectra just like an absorption gas cell. This technique produced satisfying results, with RV precision up to 5-10 m/s [80]. The drawback of this method is that the telluric lines are dependent on atmospheric conditions such as temperature, humidity and wind [81]. On the other hand, the presence of the dense telluric features reduces the number of spectral orders useful for the RV measurement. Therefore, several methods have been proposed to correct for the effect of those lines from the stellar spectra, allowing one to work with a cleaner and optimized spectrum (e.g., [82][83][84]).

Spectrographs for the NIR Spectral Domain
Starting in the 1990s, the first examples of high-resolution echelle spectrographs in the NIR band became available, allowing the pioneering studies of, e.g., young stars, M and L dwarfs and red giants (Section 2.3). However, the RV uncertainties were of the order of hundreds of m s −1 , preventing the detection of possible low-mass companions. Among those instruments, the Cryogenic Echelle Spectrograph (CSHELL) [86] at the NASA 3-m Infrared Telescope Facility (IRTF, Maunakea Observatory, USA) covered a wavelength range between 0.95 and 5.4 µm, allowing low-, intermediate and high-resolution modes (up to 40,000). The same spectral coverage is also offered by the Near InfraRed SPECtrograph (NIRSPEC) [87] at the Keck II 10-m telescope, providing a high-resolution of 25,000 (a low-resolution mode of R ∼2200 is also available). The Phoenix spectrograph [88] was previously mounted at the telescopes of the Kitt Peak Observatory before it was adapted for the 8 m Gemini South Telescope (Cerro Pachon, Chile). Its resolving power ranges between 50,000 and 80,000, and it covers the NIR range between 1 and 5 µm. For years, the state-of-the-art of NIR spectrographs were represented by CRIRES (CRyogenic high-resolution InfraRed Echelle Spectrograph) [89] at the ESO 8-m at VLT, operational since April 2007. CRIRES provided R = 100,000 in the wavelength range 0.95-5.38 µm. However, since it was a single-order spectrograph, the single exposure was limited to a narrow spectral range (∼1/70 of the central wavelength), resulting in a demanding observing request to obtain a sufficient wavelength coverage for the measure of reliable RVs. Nevertheless, several authors obtained remarkable results in the field of exoplanets thanks to very good RV precision by using different techniques [76,80]. To fulfill the increasing request of RV search for planets observing programs, a few years ago, a significant refurbishment has been proposed for CRIRES in order to improve the observing efficiency, giving rise to the CRIRES+ project [90]. The introduction of a full cross-dispersing element, as well as larger detectors will ensure an extension of the wavelength coverage for a single spectrum (∼1/7 of the central wavelength), by maintaining both the previous spectral coverage and the high resolution. Moreover, the new equipment foresees the insertion of an absorption gas cell filled with a suitable mixture of gas molecular species (acetylene, ammonia and an isotopologue of the methane) able to cover both the H and K bands. Since the region of maximum spectral emission of M dwarfs is located in the K band, the CRIRES+ cell is particularly optimized in that range ( Figure 4). According to the simulations, the use of the cell would allow reaching an RV measurement precision of ∼3 m s −1 [78].

First Multi-Wavelength Investigations
All the spectrographs previously mentioned operate in the VIS "or" in the NIR bands. Since the second half of the 2000s, some authors proposed the exploitation of the combined VIS-NIR information to obtain a more comprehensive description of the stellar activity. Actually, the science cases exploiting the NIR high-resolution spectroscopy illustrated in Section 2.2 obtain an even greater advance when studied with both VIS and NIR spectral bands. Those studies show the actual potential of this method, even if the first NIR spectrographs were not optimized for high-precision RV measurement.
The first investigation [91] that compared VIS and NIR RV datasets involved the observation of the brown dwarf LP 944-20. VIS data, obtained with UVES (UV-Visual Echelle Spectrograph) at VLT showed the significant RV dispersion of ∼3.5 km s −1 and a periodic pattern every 3.7 h, suggesting the presence of a giant planet companion very close to its host ( Figure 5, black circles). Anyway, both the RV amplitude and the period were unambiguously ruled-out by NIR data from NIRSPEC: their flat behavior indicates that the observed variation is not due to a reflex motion, but it is rather related to the fast rotation period of the brown dwarf. The presented NIR RV dispersion is ∼360 m s −1 , in full agreement with the prediction that activity-related processes are wavelength dependent and less prominent in the NIR with respect to the VIS band. On the other hand, the observations presented by [91] were not simultaneous, so they could not attest to the actual situation of the stellar activity at the time of the measurements, but they are surely conclusive about the fact that the signal observed in the VIS is not consistent with a giant planet orbiting the target. As Figure 5 shows, the uncertainties on the single RV are quite large (≥500 m s −1 for both the spectrographs), since the two instruments are not optimized for high-precision RV measurements. Anyway, due to the large amplitude of that specific signal, these uncertainties are appropriate to perform the proposed analysis. A similar approach was applied to a small sample of three young (∼2 Myrs) T Tauri stars [92] by using the VIS Coudé echelle spectrograph at the 2.7-m telescope at the McDonald Observatory and the CSHELL NIR spectrograph at the NASA IRTF 3-m telescope. One star of the sample shows a significant correlation between RV and bisector span, suggesting that the VIS RV modulation is due to the typical prominent activity of a young star, while the remaining targets show no correlation between RV and bisector. As in [91], the NIR data excluded planetary companions around the three targets, since the comparison between multi-wavelength RVs shows the same behavior as in Figure 5. This work demonstrates that the lack of correlation between RV and bisector span in this kind of objects is not sufficient to claim the presence of a planet, validating the simulation performed by [34] that demonstrates the dependence of such a correlation from the stellar rotation velocity and the resolution of the spectrograph used. Similar studies on young stars followed, including dedicated projects to search planets in their early stage of formation and evolution [93] or showing that in the particular case of very large and long-lived cool spots, the RV amplitude observed in the K-band can be significantly larger than in the VIS [94,95].
Among those studies, the one that became the most emblematic example of the power of the multi-wavelength observations method was the retreat of the hot Jupiter claimed around the classical T Tauri star TW Hya (about 10 Myrs old). In this case, the giant planet was claimed with data from the visible spectrograph FEROS (Fibre-fed Optical Echelle Spectrograph, at the 2.2-m Max-Planck-Gesellschaft ESO telescope in La Silla, Chile), which suggested the presence of such a companion with a period of 3.56 days and a mass of about 10 M Jup [96]. Even in this case, the suspect RV variation was found to be related to the stellar activity by NIR observations with CRIRES [97].
As previously mentioned, combined VIS-NIR observations are also suitable to investigate the presence of exoplanets around intermediate-mass stars, in order to understand how stellar evolution affects planetary systems. The optical RV survey on K giants presented by [98] allowed selecting low-mass planetary/brown dwarfs candidates, confirmed after NIR observations with CRIRES.

The Goal of the Simultaneity
The studies presented in the previous section showed how the combination between VIS and NIR data provides useful information on the origin of RV signals. All of them were necessarily performed with data obtained from observations collected at different observatories, where the VIS and NIR instruments are available, possibly with a time separation of months or even years. This means that a single research team should submit at least two different proposals to obtain telescope time with VIS and NIR spectrographs. Most important is the consequent loss of observing efficiency, since the VIS-NIR dataset of the same target must be collected twice. from the scientific point of view, the impossibility to obtain simultaneous VIS-NIR data prevents a correct interpretation of the behavior of stellar activity, since the two datasets monitor two different periods of the stellar activity cycle. The huge potential of the combined VIS-NIR observations has been caught in the last few years by several national and international consortia that started to design new facilities allowing simultaneous VIS-NIR high-resolution spectroscopy aiming to answer the new key questions of exoplanet research. CARMENES, GIARPS and NIRPS + HARPS, described in the next sections, are the first examples of the technological efforts made to fulfill the new scientific requirements.
The first result obtained from partially-simultaneous VIS-NIR observations regards the detection of a substellar companion with minimum mass m p sin i = 10.78 ± 0.12 M J and an orbital period P = 101.54 ± 0.05 days around the K giant star TYC 4282-605-1 [99]. This target was at first monitored with the VIS spectrograph HARPS-N at TNG, showing an RV modulation that could have been caused by a planet or stellar oscillations. A second dedicated monitoring has been arranged combining observations with HARPS-N and the NIR high-resolution spectrograph of TNG, GIANO, in its first configuration (see Section 3.2), obtaining thus quasi-simultaneous observations. Since the two spectrographs were mounted at the two opposite Nasmyth foci of the telescope, a certain amount of time was necessary to change the telescope configuration, so the time difference between the VIS and the NIR spectra in the same night was about a few hours, which is a very short time compared to the period of the candidate. The two datasets showed that the VIS and NIR RV semi-amplitudes were compatible (Figure 6), providing a "quick and easy" detection of the planet, avoiding performing two separate observing runs with a VIS and an NIR instrument. This work is the very first example of cooperation between two spectrographs that would have been part of the GIARPS instrument, introduced in Section 3.2.

New Breakthrough Facilities
A comprehensive overview of the current and future facilities providing high-precision RV in the visible and near-infrared bands (or both) is presented in [100] and partially resumed in Appendix A of this manuscript. In the following, I will be focused on those instruments that allow simultaneous multi-band observations. Two of them are already operational, i.e., CARMENES (Section 3.1) and GIARPS (Section 3.2), while others are ongoing, like HARPS + NIRPS (Section 3.3). All of them represent the groundwork of future facilities exploiting the large collecting area of new generation ELTs (Extremely Large Telescopes; Section 3.4).

CARMENES at the 3.5-m Telescope at Calar Alto
The primacy of the high-resolution spectrograph CARMENES (Calar Alto high-Resolution search for M dwarfs with exo-Earths with Near-infrared and optical Echelle Spectrographs) [101] is to be the first instrument to be specifically designed to provide simultaneous VIS-NIR observations. Moreover, the main science case of this spectrograph, i.e., the detection of low-mass rocky planets in the habitable zone of M-dwarfs, is openly declared in its name. The final concept of CARMENES had been established a few years after the first proposal in 2005. The first study considered only an NIR spectrograph, called NAHUAL (Near-infrAred High-resolUtion spectrogrAph for pLanet hunting) [102], presented as one of the instruments for the 10.4-m GTC (Gran Telescopio de Canarias) in La Palma (Canary Islands). After several modifications and the merging with a similar spectrograph and science driver, but in the VIS range, the CARMENES Project, equally led by a Spanish -German collaboration (the list of all the Institutions belonging to the CARMENES Consortium is available here: http://carmenes.caha.es/ext/consortium/index.html), was proposed and selected in 2009 for the 3.5-m telescope of the CAHA Observatory (Centro Astronómico Hispano-Alemán de Calar Alto, Almeria -Spain). A more detailed overview of the CARMENES history is provided in [103]. The commissioning of CARMENES occurred at the end of 2015, lasting for six weeks [104].
The instrument is composed of two similar highly-stabilized, high-resolution echelle spectrographs, with an optical and optomechanical design presented in [105]. The visible channel covers the wavelength range from 0.55-0.95 µm with a spectral resolution R ∼94,000, while the NIR channel covers the 0.95-1.7 µm interval with R ∼80,000. The thermal stability is guaranteed by a dedicated cooling system ensuring constant temperature inside vacuum tanks in which the two spectrographs are located, within few hundredths of a degree in 24 h [106]. While the VIS instrument operates at room temperature, the NIR one requires a cooler environment of about 140 K.
The light coming from the target is split with a dichroic, which feeds the VIS and NIR channels through fibers. Besides the scientific fiber, a second one is used to feed the spectrographs with a simultaneous calibration unit, by using an emission-line lamp or an FP (or the sky background in the case of faint objects).
Thanks to the monitored and stabilized environment and the simultaneous calibration technique, CARMENES is able to reach an RV precision of 1 m s −1 , necessary to reach the objectives of its primary science, as verified in the commissioning run [104]. Figure 7 shows the full simultaneous VIS and NIR spectrum of an M-dwarf obtained during the CARMENES commissioning. This plot demonstrates the huge information content of such a kind of instrument. The data reduction is assigned to the pipeline CARACAL (CARMENES Reduction And Calibration) included in the CARMENES data flow [107], while the RVs from CARMENES spectra are obtained with the dedicated code, SERVAL (SpEctrum Radial Velocity Analyser) [108], which creates a template spectrum for each star, starting from a number of observed spectra, and performs an iterative least-squares fitting between the template and the single spectrum. After the commissioning, in 2016, the Guaranteed Time Observations (GTO) program, which has been assigned 750 nights over five years, has started its operations to monitor a sample of more than 300 M dwarfs [109] (see Section 4.1 for a summary of the first results). Since the main goal of CARMENES is the detection of low-mass rocky planets in the HZ of M dwarfs, a wide recognition of literature data and specific analyses have been necessary for a careful selection of the sample. The coordinated effort for the characterization of the input catalog included: (i) the analysis of low resolution spectra to derive spectral types, surface gravity, metallicity and chromospheric activity for all the stars [110]; (ii) the census of binary companions to the stars through high-resolution images [111]; (iii) the measurement of stellar rotation velocity and stellar activity indices [112]. The complete database, named CARMENCITA (the CARMENES Cool dwarf Information and daTa Archive), is composed by ∼2200 M dwarfs [110], of which ∼300 were discarded due to stellar multiplicity. Finally, the GTO sample includes 324 stars [109] from the CARMENCITA catalog.   [20]. At the same time, a second high-resolution echelle spectrograph, GIANO [114], has been mounted at TNG, as well. It works in the NIR range, from 0.95-2.45 µm, at a spectral resolution of 50,000. Despite it being designed to be mounted at the Nasmyth-B focus receiving a direct light feed from the telescope, it was finally installed in the other focus chamber (Nasmyth-A) due to the scheduling constraints of TNG. In that configuration, the cryogenic vessel containing GIANO was detached from the telescope, which implied the use of ZBLAN fibers to feed the spectrograph, maintaining the wide wavelength range of the instrument, from the Y-K band. As a drawback, the adopted fiber introduced a modal noise [115] in the acquired spectra, which significantly degraded the overall instrumental efficiency [116].
A few years later, the coordinators of the Progetto Premiale WOW, "A Way to Other Worlds" (a funding scheme of the Italian Ministry of Education, University and Research promoting the cooperation and the connection of the Italian extrasolar planets community through the realization of several projects.), proposed a full exploitation of the potential of HARPS-N and GIANO, aiming to obtain a facility providing simultaneous VIS-NIR high-resolution spectroscopy. To reach this goal, it was necessary to join the two spectrographs by moving GIANO and allowing it to share the same focus of HARPS-N. A feasibility study to allow the direct feeding of GIANO from the telescope has been performed and included the definition of a new pre-slit system [73] and the building of a platform attached to the telescope fork to host the cryogenic dewar. Thanks to the new setup and the consequent removal of the fibers, the instrument, now called "GIANO-B", shows the expected efficiency [117]. After this operation, GIANO-B and HARPS-N are able to work together in the GIARPS (GIAno-b and haRPS-n, see Figure 8) [118] configuration thanks to the insertion of a dichroic in the light path that splits the entrance beam into the VIS and NIR components. An entrance slider allows choosing the preferred observing mode: HARPS-N only, GIANO-B only and GIARPS.
A uranium-neon lamp is used to obtain the wavelength calibration of the GIANO-B spectra, but the optimization of an FP interferometer is ongoing to provide simultaneous calibration and better accuracy. A preliminary study to build a gas absorption cell for GIANO-B has also been performed [119], with similar characteristics with respect to the gas cell for CRIRES+. The data reduction is performed in real time by the pipeline GOFIO [120], which processes calibration images and scientific spectra during the observing night. Due to the versatility of an instrument such as GIARPS, several science cases have been identified [118], but the main scientific driver is the search and characterization of exoplanets. Together with the TNG and the Observatories of Arcetri and Padova (the full list of institutes belonging to the GIARPS Consortium are listed on the web page: http://www.tng.iac.es/news/2017/04/04/ giarps/), the Italian community GAPS (Global Architecture of Planetary Systems; see, e.g., [121]) had a significant role in the building of GIARPS, as well as in its scientific validation. In order to exploit the capabilities of the instrument, an observing large program has been proposed by GAPS, aiming to study the frequency of hot young planets around stars having a young (age < 0.1 Gyrs) and intermediate age (<0.7 Gyrs) and to characterize hot Jupiters'/Neptunes' atmospheres. First observations are ongoing in parallel with the instrument characterization. In this framework, the instrumental stability of the red arm of GIARPS (i.e., GIANO-B) has been tested, since the one of HARPS-N is already well known and robust. To this aim, systematic observations of an RV standard have been collected over one semester in GIARPS mode. HD 3765 is a bright quiet star (K2V, V mag = 7.36, H mag = 5.27) with a low intrinsic RV variation of ∼2.4 m s −1 [32]. The NIR RVs are obtained as the weighted average RV for each exposure and its corresponding uncertainty, producing an internal error of 23 m s −1 (see also [122] for a detailed description of the RV extraction from GIANO spectra and Section 3.1 in [123] for an update), while the VIS counterpart shows internal errors of 0.3 m s −1 . This gap must be mainly ascribed to the different instrumental stability and different wavelength calibration methods, which are extremely optimized for HARPS-N. Figure 9 shows the resulting RVs, both in the VIS (blue open diamonds) and NIR (red dots) ranges. While the VIS RVs present an rms scatter of 2.3 m s −1 , in agreement with previous observations, the NIR RVs present an rms scatter of 17 m s −1 over a semester, representing the preliminary long-term stability of GIANO-B. It must be taken into account that the external contributions to HARPS-N RV noise, e.g., guiding errors, photon noise, wavelength calibration, are generally lower or much lower than 1 m s −1 , while these effects have more impact on the GIANO-B stability. In this case, some temperature drift can occur, and the wavelength calibration, obtained from the U-Ne lamp, cannot provide the same RV precision as for HARPS-N. The forthcoming use of an FP should help to limit this source of noise.  [124] is going to be the NIR counterpart of the HARPS spectrograph at the ESO 3.6-m telescope. The leadership of the building consortium is shared by the Geneva University (to exploit the expertise gained in the building of HARPS) and Montréal University (the full list of the NIRPS consortium can be found at http://www.eso.org/public/teles-instr/lasilla/ 36/nirps/). Similarly to CARMENES, NIRPS was initially proposed for a different telescope, as an answer to the 3.58-m ESO NTT (New Technology Telescope in La Silla, Chile) call for new instruments in 2014. At that time, another project was selected (the NIR medium-resolution spectrograph SOXS (Son Of X-Shooter [125]), but considering the high scientific impact of an instrument such as NIRPS, in particular, if coupled with HARPS, the ESO Board asked the proposers to adapt the instrumental design from the Ritchey-Chretien configuration of the NTT to the Cassegrain of 3.6 m. The wavelength coverage of NIRPS ranges between 0.98 and 1.80 µm (from the Y-H band), renouncing thus to the K band to favor the simplicity of the design, even if a possible extension to that range is not excluded in the future.

NIRPS (Near Infra Red Planet Searcher)
Just like HARPS, NIRPS will be installed inside a cryostat ( Figure 10) with controlled temperature and pressure, ensuring an accurate monitoring of the instrumental drift. NIRPS will benefit from a dedicated adaptive optic (AO) system [126] that allows reducing the light beam size and, as a direct consequence, the size of the instrument due to a more simple optical design. The nominal configuration foresees the concentration of the starlight into a fiber corresponding to 0.4 arcsec on the sky, producing a spectral resolution R = 100,000. In the case of faint targets or bad viewing condition, it will be possible to use a fiber with a larger field of view, 0.9, for a more efficient photon collection, with R = 75,000. This is the first practical high-performance attempt to combine AO with a fiber-fed high-resolution spectrograph. The cooperation between the dedicated AO system (that will scramble the modes at the fiber entrance) and the combination of three scrambling techniques (octagonal fibers, fiber stretcher and a double scrambler) allows mitigating the typical modal noise [124] that affected, e.g., GIANO in its previous configuration (Section 3.2). Users could choose between a laser frequency comb and a stabilized FP calibration system to obtain very high-precision RV. Indeed, the expected RV accuracy of NIRPS will not be very different from HARPS, reaching 1 m s −1 in less than 30 min for an M3 star with an H magnitude of 9 [74]. The overall wavelength coverage of the NIRPS + HARPS configuration will be 0.4-1.8 µm, except for a small window between 0.7 and 0.95 µm (i.e., the gap between the two instruments) that is used for the wave-front sensing of the AO system. NIRPS's first light is expected by 2019. The three science cases identified for the five years of GTO of NIRPS (combined with HARPS) are: (i) the characterization of transiting Earth-mass planets around M dwarfs; (ii) the selection of candidates for high-contrast imaging performed by ELT class telescopes; (iii) atmospheric characterization of exoplanets. The NIRPS + HARPS combination is then going to play an outstanding role in the context of the full characterization of potentially habitable planets in the near future. As soon as the first M dwarf planet candidates from the NASA TESS satellite (Transiting Exoplanet Survey Satellite, see https://tess.gsfc.nasa.gov/science.html [127]) are released, an RV follow-up with NIRPS + HARPS would provide their masses. Since the planet radius is measured with the transit observation, it is then possible to place those values in the mass-radius diagram (see, e.g., Figure 4 in [128]), which is the first step to reveal the planet composition, even if a very small uncertainty is required to avoid degeneracy among different models. Those planets could subsequently be subject for atmospheric characterization with JWST (James Webb Space Telescope, https://jwst.stsci.edu/science/sciencecorner/white-papers), and their host stars could be investigated to search for further companions on wider orbits through imaging observations or through the astrometry method by using the huge Gaia database (http://sci.esa.int/gaia/).
Finally, investigations of molecular features of hot Jupiters and hot Neptunes in the combined wavelength domains will be performed. As in the case of the other mentioned facilities, these studies will allow us to create a sort of chemical species database, useful for future investigations, as well as to validate laboratory experiments [129]. Anyway, the most important aim of such a type of detection is to explore and/or refine the techniques to recover atmospheric planetary signals waiting to investigate the atmosphere of potentially habitable rocky planets by using high-resolution VIS-NIR spectrographs installed at telescopes with larger collecting area.

Future Instrumentation
The exploitation of the huge collecting area of the future Extremely Large Telescopes (ELTs) represents an unprecedented opportunity and feeds great expectations for the search for exoplanets and characterization science cases. The quest for an Earth-analog planet around a solar-like star is the primary objective of the new-generation instruments, and the synergy among the next big international projects is going to make the difference. Among those instruments, the Giant Magellan Telescope (GMT, built by a consortium including American, Australian and Korean institutions, see https://www.gmto.org/) [130] is a 25.4-mdiameter segmented mirror telescope that will be based at Las Campanas Observatory (Chile), while the European ELT (E-ELT, an ESO facility, see https://www.eso.org/sci/facilities/eelt/) is a 39-m diameter segmented mirror telescope, which will be located at Cerro Armazones (Chile). Both of them are expected to be commissioned in 2024 and are planned to mount VIS and NIR high-resolution spectrographs.

•
G-CLEF is a fiber-fed, cross-dispersed echelle spectrograph built for general purposes, but designed to provide also high precision RVs (up to 10 cm s −1 [131]). The wavelength range covered by G-CLEF spans between 0.35 and 1 µm to satisfy the requirements of the main science cases. The bluer part is included to perform studies on very metal-poor stars since the iron content is evaluated through very faint features between 3500 and 3900 Å, undetectable from the existing facilities in terms of high resolution and telescope aperture. The wavelength coverage is pushed up to 1 µm to exploit the spectral information of the M dwarfs and obtain precise RVs for these stars. Among the other exoplanet science cases, it is noteworthy to mention the search for biomarkers (in particular, the features of molecular oxygen between 7600 and 7700 Å in the atmospheres of the exoplanets) and the follow-up of the planet candidates detected by the TESS satellite.
The spectral resolution can be chosen according to the science program. G-CLEF offers both high-resolution modes (R ∼105,000), available with or without an image scrambler that enhances the RV precision and medium resolution modes (35,000 and 19,000). As for GMTNIRS, particular attention is paid to providing the maximum performance in terms of throughput, to exploit the benefit of the huge telescope collecting area. To optimize the incoming stellar flux, the light beam will be split through a dichroic in two channels, red and blue (see Figure 11, upper panel), each one equipped with specific detectors and coatings. • GMTNIRS will cover the wavelength range between 1.15 and 5.3 µm making use of five spectrograph units (see Figure 11, lower panel), each one dedicated to a specific atmospheric window: J, H, K, L and M bands. The light that feeds the instrument will be subdivided into the five channels thanks to a series of dichroic elements. The spectral resolution is expected to be ∼60,000 in the JHK bands and ∼85,000 in the L and M bands. In order to test the innovative technical elements planned to equip GMTNIRS, the University of Texas and the Korea Astronomy and Space Science Institute (KASI) joined to build a sort of forerunner instrument, the cross-dispersed NIR spectrograph IGRINS (Immersion Grating Infrared Spectrometer, see http://www.as.utexas.edu/mcdonald/facilities/2.7m/igrins.html) [132]. IGRINS, which has operated since 2014, has been installed alternatively at the Harlan J. Smith 2.7-m telescope (McDonald Observatory, USA), at the 4.3-m Discovery Channel Telescope (Lowell Observatory, USA) and at the 8.1-m Gemini South Telescope. It covers the H and K windows, from 1.45-2.5 µm in a single acquisition with a resolving power of R = 45,000. The design of this instrument was particularly focused to optimize the throughput rather than to reach extreme precision in RV, and the same concept will be adopted for GMTNIRS, as well. The instrumental design originally proposed for GMTNIRS [133] is currently evolving [134,135], aiming to match the requirements of the main scientific driver emerging through the years, i.e., the observation of exoplanetary atmospheres, as for HIRES (see the discussion in Section 3.4.2). Being part of a large and international project as the GMT, the GMTNIRS spectrograph will not be focused on exoplanet science alone; it will provide a significant contribution to the study of young stellar objects, debris disks, protoplanetary systems, stellar evolution, interstellar medium and the star formation history of the Galaxy.
These two instruments alone show the peculiarity of having a larger spectral coverage with respect to "traditional" VIS or NIR spectrographs. For this reason, a sort of multi-wavelength comparison can be performed even with data from a single instrument; anyway, their simultaneous use will be extremely interesting in the framework of exoplanet study.

HIRES at E-ELT
The main scientific driver of HIRES, the High Resolution Spectrograph [136] for the E-ELT, is the atmospheric characterization of Neptune-like down to Earth-like planets that will be revealed by PLATO (PLAnetary Transits and Oscillations of stars, an ESA M-class space mission dedicated to the transit search of thousands of planets, including the terrestrial ones, and asteroseismic study of the host stars, with an expected launch date in 2026) [137]. Just like CARMENES, GIARPS and HARPS + NIRPS, the well-known potential of the combined VIS-NIR observations additionally coupled to a 40-m class telescope represents a great opportunity for exoplanets' characterization. The required accuracy for the RV precision is 10 cm s −1 . HIRES is the result of the combination between the two AO-assisted high-resolution spectrographs proposed for the Phase A study of E-ELT: CODEX (COsmic Dynamics and EXo-earth experiment) [138], operating in the VIS band (R = 120,000), and SIMPLE [139], operating in the NIR (R = 130,000). The new project thus foresees a unique instrument working in the wavelength range between 0.37 and 2.5 µm with a resolving power of 100,000.
Actually, as Figure 12 shows, HIRES should be composed by different modules, allowing the subdivision of the costs according to the resources. According to the preliminary design, there will be four fiber-fed cross-dispersed spectrograph modules, each one dedicated to a specific wavelength coverage: U-B, V-R-I, Y-J-H and K bands (the latter two would be cryogenic). HIRES will surely be the most demanding effort in the building of high-resolution spectrographs, and as Table 2 in [136] shows, it will be accomplished only with the collection of all the expertise gained in the production of the high-performance spectrographs realized in the last few years. HIRES will be a second-generation instrument, so it will not be available at the E-ELT first light. The HIRES Consortium includes institutes from several countries (see https://www.arcetri.astro.it/~hires/consortium.html), and INAF is the leading technical institute. The key science of HIRES will be the full investigation of planetary atmospheres (chemical composition, stratification and weather) of Neptune-like to Earth-like planets, including those in the HZ of their host stars [140]. The detection of biosignatures is also expected. Since a typical RV survey is very demanding from the point of view of time request, the use of E-ELT is not justified, since it will be available to a huge community having heterogeneous lists of relevant science cases. The high pressure on E-ELT instrumentation should relax the time request of other facilities. Therefore, the RV monitoring for the search of exoplanets or characterization should proceed, e.g., with ESPRESSO or other future available spectrographs at 10-m class telescopes (see Appendix A). Anyway, the extreme RV precision ensured by HIRES could be exploited in the measurement of the Rossiter-McLaughlin (RM) effect, an RV anomaly in the RV curve detectable during the transit of a planet. This effect, which represents an independent confirmation of the actual presence of the planet, allows deriving the obliquity of the planetary orbit with respect to the stellar spin axis, which is an important observable to understand the planet migration history and/or constrain the star-planet tidal interactions [141].

First Results
As discussed in Section 2.3, combined VIS-NIR observations have already been proposed by several authors, who exploited the multi-wavelength information to investigate different issues on different types of targets, even if with non-simultaneous observations. In the following paragraphs, I will present the first results of CARMENES and GIARPS, the above-mentioned simultaneous VIS-NIR facilities.

The CARMENES M Dwarf Survey
The science operations of the CARMENES GTO have been running since January 2016, and several works have been presented, both on the spectroscopic characterization of the sample [109,142,143] and on the analysis of the first RV time series. In [144], the first RV release of the VIS channel is provided. The mean RV precision of 287 targets out of 324 (in that work, only the stars with at least five RV epochs were considered) is 1.7 m s −1 , in agreement with the expectations. For targets with a spectral type later than M6, the mean RV scatter is typically larger (∼10 m s −1 with respect to the mean value of the whole sample, 3-4 m s −1 ) as the v sin i increases. An interesting result is the null-correlation between the RVs and the activity indicators for 50% of the considered sample. This could be explained by complex patterns of the dark spots on the stellar surface or related to the magnetic fields associated with those inhomogeneities or even a possible larger temperature contrast between the photosphere and the spots.
The first planetary companion detected by CARMENES is orbiting around the M0 star HD 147379 [145], located at 10.7 pc. This discovery has been presented by analyzing the VIS RV only, since the typical RV uncertainty reached in the NIR channel (8.6 m s −1 ) is not suitable to recover this kind of signal, because of the paucity of spectral lines of early M dwarfs in that wavelength range. Therefore, the VIS RV semi-amplitude is ∼5 m s −1 , providing a minimum mass of ∼25 M ⊕ (1.5-times the mass of Neptune). The orbital period is 88 days, which placed this planet inside the temperate zone around the star.
A test of the actual performance and capability of the instrument has been presented by [146], aiming to recover a number of known rocky planets around a small sample of seven M dwarfs (single-planet systems: GJ 15A, GJ 176, GJ 436, GJ 536, GJ 1148; multi-planet systems: GJ 581, GJ 876) with RV from the CARMENES VIS channel. In this work, the orbital parameters of the known systems are updated, and a new Neptune-mass companion with a period larger than 500 days in an eccentric orbit has been detected around GJ 1148. Anyway, the more intriguing result is the non-detection of the super-Earth (m sin i = 5.35 M ⊕ ) orbiting GJ 15A with a period of 11.44 days, previously claimed with HIRES data [147]. According to the data obtained during the first 16 months of monitoring with CARMENES (∼100 RVs used), no evidence of such a signal is shown (Figure 13, left panel). This fact is corroborated by the analysis of the combined HIRES and CARMENES datasets (spanning more than 7300 days): the authors conclude that this periodicity is caused by the stellar activity since the power of the generalized Lomb-Scargle periodogram (GLS) [148] of the corresponding peak decreases through the years. However, the presence of GJ 15Ab has been subsequently confirmed in a recent investigation of HARPS-N + HIRES data [149], also containing the detection of a further long period planet (P ∼7600 days, suggested, but not confirmed by [146]). In [149], the simultaneous modeling of the activity and planetary signals is performed through Gaussian processes, which retrieve the debated Super Earth (Figure 13, right panel), even if with a smaller value of the RV semi-amplitude with respect the previous finding [147], and therefore of the planet mass. According to the authors, this could be explained with a more suitable sampling of the periodicity, a higher RV accuracy of the analyzed dataset (new release of the HIRES data [11] and HARPS-N) and a better treatment of the activity contribution. Similar arguments are proposed to explain the non-detection from CARMENES data. Actually, the fitted RV semi-amplitude K b = 1.68 +0. 17 −0.18 m s −1 appears to be beyond the detectability limit of CARMENES, since its typical RV uncertainty is about 1.7 m s −1 with respect to the 0.62 m s −1 for HARPS-N. A further investigation is certainly required to define which analysis is better at interpreting the observed data. The computation of the RVs in different spectral regions is proposed [149] as an additional test to evaluate the behavior of the RVs as a function of the color, following the approach proposed by [150]. This could be done for the HARPS-N dataset since it should not be affected by a severe increase of the RV uncertainties when only a number of spectral orders are used to evaluate the RV measurement. In the case of CARMENES data, this could return a non-fully-reliable test. On the other hand, it would be very interesting to see the contribution of the NIR RVs of CARMENES.
CARMENES also performed an RV follow-up of the target K2-18b, a transiting Super Earth orbiting an M2.5 dwarf star, revealed by K2 photometry [151] (K2 is the update of the Kepler mission after its severe failure occurred in 2013) and then confirmed with the NASA Spitzer observations [152] (http://www.spitzer.caltech.edu/). Transit analysis returned a radius of about 2 R ⊕ . By using the VIS RVs of CARMENES (K b = 3.5 m s −1 ), a mass of 9 M ⊕ has been measured [153]. This implies a Super Earth having a density consistent with a solid core and a volatile-rich envelope in the temperate zone of an early M dwarf, receiving thus a similar stellar irradiation as the Earth. As a consequence, this target is a perfect candidate for the atmospheric characterization of a rocky planet in the HZ of its host star with JWST. Finally, CARMENES data tend to exclude the presence of a shorter period planet claimed with HARPS-N data [154], finding a signal dependent on time and wavelength. As the GLS periodograms in Figure 14 show, while the periodicity of the 33-day planet (0.3 1/d in the frequency domain) is clearly visible in both the first and the second halves of the data, the one at nine days (0.11 1/d) vanishes in the second part of the RV monitoring. This behavior is typical when the periodicity is due to the stellar activity and not to a Keplerian signal.  is the retreat of a controversial exoplanet around the young K5 star, BD+20 1790 (considered as a member of the AB Dor moving group, ∼150 Myrs) by [123]. This target shows an intense and peculiar stellar activity which the responsible for the significant RV modulation in the VIS range reported by [155,156]. This modulation shows an RV semi-amplitude K 1 km s −1 , which was interpreted as a signature of a massive hot Jupiter with a period of 7.8 days. In 2010, some doubts about the real presence of this companion were raised [157], because the modulation was not compatible with other spectroscopic data obtained in the VIS range in a different period, thus attributing the RV variations to photospheric processes. The new NIR RV dataset was obtained by using telluric lines from standard stars as a reference, following the approach by [122]. The comparison between the VIS archive data of BD+20 1790 and the GIANO/GIANO-B and IGRINS spectrographs ( Figure 15) shows that the average RV amplitude in the NIR is significantly lower than the VIS data (about 1/4). Moreover, the comparison between the previous optical data and the three RV measurements obtained with HARPS-N in the GIARPS configuration (black asterisks in the lower panel of Figure 15) shows that even the VIS modulation has changed with time, suggesting that it is caused by stellar activity. To have a more comprehensive indication of the photospheric activity, during the observing run of GIANO, a quasi-simultaneous photometric monitoring of BD+20 1790 had been performed with REMIR (REM InfraRed) and ROSS2 (REM Optical Slitless Spectrograph 2), two imagers in the NIR and VIS range, respectively, mounted at the REM telescope (Rapid Eye Mount, a 60-cm robotic telescope at the ESO-La Silla Observatory). ROSS2 and REMIR can operate simultaneously thanks to a dichroic. Surprisingly, the VIS light curve from ROSS2 shows a smaller amplitude with respect to the NIR one obtained with REMIR (see Figure 5 in [123]), which is not an intuitive result, according to the arguments presented in this review. Moreover, the two modulations are in anti-phase with each other. The authors propose a photometric model that includes a mixture of cool and hot spots in the same active region able to explain the peculiarity of this star. However, a similar behavior has been reported by [94,95] in the RVs of young stars (Section 2.3).
This work also demonstrates that the support of simultaneous photometry can be particularly useful both to interpret the RV curves and to add further information on the characteristics of the stellar activity.

Conclusions and Perspectives
The new scientific questions in the framework of the exoplanet search and characterization can be summarized as follow: (i) the search for Earth-mass rocky planets in the HZ of M dwarfs; (ii) the identification of the origin of planetary system diversity through the detection of planets around young stars; (iii) the characterization of the hot gas giant planets atmospheres as a laboratory for the future characterization of rocky habitable planets with ELTs or space-based telescopes. The simultaneous multi-band observations in the VIS and NIR ranges are expected to be the forthcoming parameter space for these new issues. The full monitoring of stellar activity represents the key point in the detection of exoplanets around targets with enhanced and peculiar activity such as the ones considered in (i) and (ii). On the other hand, the opportunity to explore molecular features of the planetary atmospheres, present over a wide interval of wavelengths, is now possible thanks to instruments covering the range between the bluer part of the VIS up to the edge of the NIR (iii). Instruments like CARMENES, GIARPS and the forthcoming NIRPS + HARPS configuration are expected to yield a significant contribution in these and in many other fields of astrophysics. They also represent the starting point for future generation spectrographs that are going to equip new astronomical facilities.
If we observe a planetary (period-mass) diagram (or, alternatively, a (semi-major axis-mass) diagram, e.g., Figure 1 in [158]) comparing the properties of the known exoplanets with the planets of our Solar System, it is clear that the current technology is still not sufficiently adequate to find planets like our Venus, Mars, Neptune, etc., or the Earth itself. The need to reach the 1-m s −1 (and less) RV accuracy is driven by the goal to find an Earth-like planet around a Sun-like star. As mentioned in this review, there are several impact factors on the achievable RV precision with high-resolution spectroscopy (both VIS and NIR): changes in the environmental conditions require a series of technical devices ensuring a high-level of instrumental stability; the wavelength calibration must rely on a very stable and reliable wavelength reference to reach the required precision; RV extraction methods must be more and more sophisticated. Other important technical elements can be the availability of state-of-the-art detectors and devices that allow a uniform illumination of the spectrograph slit. Last, but not least, from the observational point of view, the capability to conduct a proper data sampling of the RV signal helps to avoid the aliasing phenomenon that hampers the frequency (period) analysis. We are quite confident to say that all of these issues have been successfully addressed, at least in the visible range, at least to obtain the 1-m s −1 regime. As previously stated, the technology transfer from VIS to NIR is far from immediate, but many efforts have been already made to allow precise RV also in the NIR domain, starting from the lessons learned in the last few decades.
At this moment, we are approaching the new milestone of the Doppler method, since we are waiting for the precise NIR RV measurement from CARMENES (1 m s −1 expected), on the one hand, and the outstanding RV precision of 10 cm s −1 from ESPRESSO, on the other. On the scientific side, the challenge is to be able to manage all the other detectable effects, first of all, the stellar activity. Even in this case, valuable examples of data treatment are presented in literature, and it will be very interesting to see how they will work in the new ultra-precise RV domain. Reaching the cm s −1 regime is a necessary step of the journey, but it is also crucial to know how to exploit this opportunity. In the near future, these new facilities will necessarily be more and more specialized for RV monitoring for the search for exoplanets (some of them are already employed in such a way); just think about the follow-up of the forthcoming space satellite dedicated to transit detection. Actually, the cooperation among different detection methods (and thus, instruments and projects) will probably be the best way to characterize a star-planet system. In this context, an important role will be played by all the tools that allow obtaining a more precise estimate of the stellar parameters like mass and radius (e.g., stellar modeling and asteroseismology), since they are crucial to constrain the global planet properties.
Funding: This research received no external funding. Acknowledgments: S.B. is grateful to A. F. Lanza, I. Carleo and R. Claudi for the useful comments on this manuscript and the support with the specific topics, to C. Dressing for her kind advice and to M. Endl for the update on HRS-2. S.B. also wishes to thank the two anonymous referees that provided valid and useful suggestions that allowed this paper to be more complete and accurate. Finally, S.B. is grateful to the authors and the journal editorial boards for the permission to reproduce their plots and figures in this review.

Conflicts of Interest:
As Project Manager of the GIARPS instrument installed at the INAF TNG (Telescopio Nazionale Galileo), the author was deeply involved in its technical and scientific development, in the commissioning phase and science exploitation of the data. Nevertheless, in this manuscript, she tried to be as impartial as possible, just reporting the main characteristics and interesting results of both GIARPS and CARMENES instruments now available, as well as the forthcoming projects. The only aim of this paper is to show the outstanding potential of the multi-wavelength observations in the framework of exoplanet detection and characterization. No founding sponsors had a role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript; and in the decision to publish the results. Table A1. List of high-resolution spectrographs in the VIS or NIR range recently integrated and commissioned, or close to their first light.

Spectrograph Telescope and Site
Expected Link Availability