Systematic investigation of chemical abundances derived using IR spectra obtained with GIANO

Detailed chemical abundances of Galactic stars are needed in order to improve our knowledge of the formation and evolution of our galaxy, the Milky Way. We took advantage of the GIANO archive spectra to select a sample of Galactic disc stars in order to derive their chemical inventory and to compare the abundances we derived from these infrared spectra to the chemical pattern derived from optical spectra. We analysed high-quality spectra of 40 stars observed with GIANO. We derived the stellar parameters from the photometry and the Gaia data-release 2 (DR2) parallax; the chemical abundances were derived with the code MyGIsFOS. For a subsample of stars we compared the chemical pattern derived from the GIANO spectra with the abundances derived from optical spectra. We derived P abundances for all 40 stars, increasing the number of Galactic stars for which phosphorus abundance is known. We could derive abundances of 14 elements, 8 of which are also derived from optical spectra. The comparison of the abundances derived from infrared and optical spectra is very good. The chemical pattern of these stars is the one expected for Galactic disc stars and is in agreement with the results from the literature. GIANO is providing the astronomical community with an extremely useful instrument, able to produce spectra with high resolution and a wide wavelength range in the infrared.


Introduction
The knowledge of the chemical composition of the stars is a necessary ingredient to understand and to model the formation and the evolution of the Milky Way and of the Local-Group Galaxies. In order to derive the chemical pattern of a star, astronomers compare the strengths of the atomic and molecular lines in an observed stellar spectrum to theoretical predictions. Typically, spectra with medium/high resolving power (R > 15000), a good signal-to-noise ratio (S/N > 50), and a wide-enough wavelength range allow to derive precise abundances for several elements. Until a few decades ago, securing stellar spectra was the bottleneck in the chemical investigation of stars. Any single star had to be observed for a reasonable amount of time, depending on its luminosity, the size of the telescope, and the efficiency of the spectrograph, in order to obtain a spectrum with the quality necessary to derive the desired precision in the chemical abundances. For instance, with a four-meter-class telescope, a spectrum with a typical S/N of about 50 at 520 nm is obtained in one hour for a V = 12 mag star (e.g. HARPS Mayor et al., 2003) while with an eight-meter-class telescope, a spectrum of ⋆ GIANO programme A31TAC the same S/N is obtained in one hour for a star 300 times fainter in the V band (e.g. FLAMES Pasquini et al., 2002). The advent of the multi-object spectrographs (e.g. FLAMES Pasquini et al., 2002) allowed one to observe of the order of one hundred stars at the same time, considerably increasing the number of stars in the Galaxy for which the chemical composition is known. The Gaia-ESO Survey (Gilmore et al., 2012) is an example of the great efficiency that can be reached by using multi-objects spectrographs. In the near future multi-object spectrographs with a capacity to observe about ten times more stars at the same time are expected (e.g. WEAVE: Dalton et al. 2016 and4MOST: de Jong et al. 2014).
Before the commissioning of the near-infrared (NIR) high-resolution spectrograph GIANO  at Telescopio Nazionale Galileo (TNG), the detailed chemical investigations of stars were based mainly on optical spectra. Other IR spectrographs appeared before GIANO, such as for example CRIRES (Kaeufl et al., 2004). However, their wavelength coverage was small, useful only to add one element on the chemical pattern of a star, and did not allow us to perform a systematic investigation of the IR range. GIANO is a crossdispersed echelle spectrograph, and covers a large wavelength range (950-2450 nm) allowing the abundance of several elements to be derived. The multi-object spectrograph APOGEE (Allende Prieto et al., 2008), with a wavelength range of about 200 nm (1510-1700 nm), has been shown to be able to provide detailed chemical composition of Galactic stars. In the future we expect the Multi Object Optical and Near-infrared Spectrograph (MOONS, Cirasuolo et al., 2014), that will be able to observe several hundred objects in three arms with two resolution options. MOONS is a high-multiplicity, multi-object, fibre-fed spectrograph for the ESO VLT. It will be capable of both low (R∼ 5 000) and high (R∼ 20 000) resolution and will deploy 1000 fibres on the 20 ′ field of view of the VLT. Its high-resolution H arm will observe the wavelength range 1521-1641 nm, covered by the orders 47-50 of GIANO.
The present work has two main drivers. On the one hand, we are deeply involved in the preparation of MOONS and wish to be prepared for the millions of spectra it will provide. Therefore a detailed study of the IR range is in order. On the other hand, we are interested in the Galactic evolution of phosphorus, which is still a very little studied element, and of sulphur, both elements being available in the wavelength range of GIANO. Our main goals in this work are: (i) to investigate the IR range and see which elements can be derived and (ii) to derive the abundance of phosphorus. We also wish to provide detailed abundances for several elements for those stars in this sample for which only the abundances of a few elements are known. We took the advantage of the GIANO spectra present in the Italian centre for Astronomical Archives (Molinaro et al., 2018) and we selected a number of dwarf stars (the IR P i lines are contaminated in the spectrum of a giant star) of slightly sub-solar metallicity to super-solar metallicity. For a subsample of stars, SOPHIE (Bouchy & Sophie Team, 2006) spectra of good quality were available, allowing us to compare for some elements the abundances derived from the IR and the optical spectra.

Observations
Stars were observed with GIANO, the high-resolution (R≃50 000) IR (950-2450 nm) spectrometer of the TNG . The instrument was designed for direct feeding of light at a dedicated focus of the TNG. In 2012 the spectrometer was provisionally positioned on the rotating building and fed via a pair of IR (ZBLAN) fibres connected to another focal station . In 2016 it was eventually moved to the originally foreseen configuration where it can also be used for simultaneous observations with HARPS-N (Tozzi et al., 2016). The spectra presented here were collected in March 2015, when the spectrometer was still in its provisional configuration and fed via fibres. The raw data were downloaded from the public archive of the TNG. The stars were observed by nodding on fibre, that is, target and sky were taken in pairs and alternatively acquired on the two fibers (A,B) for an optimal subtraction of the detector noise and background. Integration time was 5 minutes per A,B position. Dark-subtracted frames of flat (halogen lamp) and wavelength (U-Ne lamp) calibration sources taken at the beginning and/or at the end of each night were used to reduce the data. For each nodding pair we created a two-dimensional (2D) frame by subtracting the raw frames (A-B) and dividing the difference by the normalised flat. Spectral extraction and wavelength calibration were performed using a physical model of the spectrometer that accurately matches instrumental effects such as variable slit tilt and orders curvature over the echellogram (Oliva et al., 2018). The telluric absorption features were corrected using the spectra of two telluric standards (O-type stars) taken at different airmasses during the same nights. The normalised spectra of the telluric standard taken at low and high airmass values were combined with different weights to match the depth of the telluric lines in the stellar spectra.
In Fig. 1 the H-band is shown in the case of the star HD 24040; this range corresponds exactly to the wavelength range that will be observed at high resolution by MOONS.
The SOPHIE spectra have been retrieved from the SOPHIE archive 1 . Several spectra for each star were present in the archive, but only the observations with a corresponding sky spectrum were retained. We subtracted the sky, shifted the spectra at laboratory wavelength and, in the case of several observations, added them together. For one star, HD 76909, we had a single observation.

Stellar parameters
Within the stars observed by GIANO available in the archive we selected those which had Gaia DR2 (Gaia Collaboration et al., 2018b,a) parallax as well as 2MASS (Cutri et al., 2003) and Tycho (ESA, 1997) photometry. To derive effective temperatures we used the colour-temperature calibrations of González Hernández & Bonifacio (2009) that are based on temperatures derived through the IR flux method for a sample of calibrators. To do so we used the colour transformations given by equations 1.3.20 and 1.3.26 in the Hipparcos and Tycho catalogues (ESA, 1997) 2 to obtain V an (B − V) in the Johnson system. All the stars are fairly near and little reddened, nevertheless we derived the reddening for each star from the reddening maps of Capitanio et al. (2017) and corrected the colours for reddening. We then applied the equations given in table 5 of González Hernández & Bonifacio (2009) and (V − K) that were averaged among the different calibrations. The colours involving the 2MASS photometry were retained only if the corresponding quality flag was A. This implies that not all stars have the same number of colour determinations. In Table A.1, T eff is provided with an uncertainty which is the scatter among the temperatures derived by the five calibrations. For one star, HD 190360, T eff is derived from one single calibration (that using the B − V colour), and we assumed 250 K as uncertainty in the temperature; for two stars, HD 28005 and HD 34575, four calibrations were used (those based on colours B − V, V − J, V − K and J − K), and for all the other stars the five calibrations were used. The effective temperature of this sample of stars spans from 5006 K to 6292 K, with T eff = (5727 ± 245) K.
For each calibration used to derive the effective temperature, we derived the corresponding surface gravity from the Stefan-Boltzmann relation using the distances from the parallaxes provided by Gaia DR2 and the stellar absolute luminosity, derived from the Gaia G band and a bolometric correction based on three-dimensional (3D) model atmospheres (Bonifacio et al., 2018). We assumed the mass of each star according to its tem- perature, by looking at the isochrones. The uncertainty in the gravity is the scatter among the values we derived, and this value is on average below 0.05 and for all stars below 0.1. The gravities we derived spanned from 3.53 to 4.56 [c.g.s], with log g = 4.16 ± 0.25 (see Table A.1).
The stars we are analysing have a metallicity from slightly sub-solar to super-solar, meaning that the lines are strong and several lines are blended. Also, to derive the micro-turbulence in the IR range is not trivial since the lines we can detect are stronger than in the optical range. We then fixed the value for the microturbulence at 1.0 km/s, which is a good value for this sample of dwarf stars. Also for the Sun a microturbulence at 1.0 km/s is the preferred value (Steffen et al., 2013).
The chemical analysis has been done with the pipeline MyGIsFOS (Sbordone et al., 2014), in the complete wavelength range provided by GIANO, except the ranges that are too heavily contaminated by telluric absorption. This choice was also motivated by two facts: (i) the GIANO archive is providing more and more stars that can be analysed; (ii) the future MOONS-GTO as well as the future MOONS observations will imply several thousand spectra each night. The numbers of stars to analyse is too large to reasonably allow a by-hand analysis. To analyse the present sample of stars, MyGIsFOS operated at fixed T eff and log g and derived, from selected atomic lines, the chemical abundances by comparing the observed spectrum to a grid of synthetic spectra. The theoretical synthesis was computed by SYNTHE (see Kurucz, 2005;Sbordone et al., 2004) from a grid of ATLAS 12 models (Kurucz, 2005) with a step of 0.25 dex in metallicity. The atomic parameters of the lines are from the most recent line-list provided by R. Kurucz on his site 3 ; in the H band we used the atomic data provided by Shetrone et al. (2015). The stars we analysed all had metallicities from slightly sub-solar to super-solar, and from Fe i lines we derived for the sample of stars: [Fe/H] = +0.10 ± 0.11; for all the stars −0.14 ≤ [Fe/H] ≤ +0.33. Similar results were derived using Fe ii lines, −0.07 ≤ [Fe/H] ≤ +0.40, with an average difference of [Fe/H] from Fe i and Fe ii lines smaller than 0.01 dex, the largest one being just smaller than 0.10 dex. The difference of [Fe/H] from Fe i and Fe ii lines is in any case smaller than the uncertainties (see Table A.1).

Comparison optical -IR
For nine stars (see Table 1), we had optical SOPHIE spectra for which one fibre was on the sky. We could then subtract the sky and compare the analysis of GIANO and SOPHIE for the elements in common.
Eight elements are in common in the two analyses and the agreement is generally very good.  Fig. 2. Comparison of the abundances derived from IR and optical spectra.
-The agreement in the Fe abundance derived from optical and IR spectra is good: For one star in particular (HD 76909), the difference on A(Fe ii) from optical and IR spectra is large (0.18 dex) but still within uncertainties; in fact it is the optical spectrum that shows a larger disagreement in the A(Fe) from Fe i and Fe ii lines and the line-to-line scatter from the Fe ii lines is the largest among the optical spectra. -There is a systematic difference in A(Na) from IR and optical spectra: A(Na) IR − A(Na) opt = −0.23 ± 0.14. From the IR and optical spectra we have [Na/Fe] = −0.11 ± 0.09 and [Na/Fe] = +0.11 ± 0.08, respectively. This difference could be due to NLTE effects; in fact the IR lines form close to LTE, but this is not true for the optical lines, which are sensitive to NLTE. -The agreement in A(Al) is really good; in the case of IR spectra, in seven cases out of nine, A(Al) is derived from one single line and in the other two stars from two lines: A(Al) IR − A(Al) opt = +0.07 ± 0.06. -For A(Si) from Si i lines, the agreement of IR and optical spectra is good, well within the uncertainties: A(Si) IR − A(Si) opt = −0.08 ± 0.04. -Also the agreement for A(S) is extremely good: A(S) IR − A(S) opt = +0.03 ± 0.04 (see Fig. 2). -For A(Ca) from Ca i lines the agreement is also perfect: Fig. 2). -The Ti abundance, derived from Ti i and Ti ii lines is very close from IR and optical spectra: -For Ni we also find a good agreement between IR and optical: A(Ni) IR − A(Ni) opt = +0.02 ± 0.03 (see Fig. 2).
A summary on the comparison of the abundances derived from IR and optical spectra is presented in Table 2.
Previously published investigations exist for all nine stars from Table 1. Three of the papers that contain parameters for some of these stars are based on spectra (Valenti & Fischer, 2005;Petigura & Marcy, 2011;Brewer et al., 2016) and two are based on photometry (Masana et al., 2006;Casagrande et al., 2011). All the spectra used in Valenti & Fischer (2005)  In the following we discuss abundances in terms of A(X) = log(X/H)+12.

HD 34575
For this star we derived T eff = 5582 K, in very good agreement with 5529 K derived by Masana et al. (2006) and 5546 K derived by Casagrande et al. (2011). The analysis by Brewer et al. (2016) is in very good agreement with our investigation; they derive very close stellar parameters (T eff = 5551 K, log g = 4.25) and also the abundances (C, Na, Al, Si, Ca, Ti, Fe) are in agreement with our findings, well within the uncertainties. Valenti & Fischer (2005) derived a close effective temperature, 5651 K, and a gravity 0.23 dex larger than our adopted value. The abundances they derive (Na, Si, Ti, Fe, Ni) are mostly slightly larger, but well compatible within the uncertainties. With the stellar parameters from Valenti & Fischer (2005), Petigura & Marcy (2011) derived a C abundance for this star that is compatible with our value, but larger, as in the abundances from Valenti & Fischer (2005).

HD 67346
The temperature we derive for this star of 5953 K is in close agreement with T eff = 5941 K from Masana et al. (2006)

HD 69056
The effective temperature we derive (T eff = 5637 K) for this star is in very good agreement with the values from Masana et al. (2006) and Casagrande et al. (2011) (T eff of 5598 and 5635 K, respectively). Petigura & Marcy (2011) used T eff = 5490 K and log g of 4.27 and derived C and Ni abundances in agreement with our results.

HD 69960
Our T eff = 5655 K is in extremely good agreement with T eff = 5625 K by Casagrande et al. (2011) and also in agreement with the cooler value of Masana et al. (2006)  We are, by and large, in agreement with Chen et al. (2008), in spite of their effective temperature being about 200 K cooler.

HD 73933
The effective temperature we derive (6243 K) is about 100 K hotter than the values from Casagrande et al. (2011) and Masana et al. (2006). Petigura & Marcy (2011), with an effective temperature of 6076 K and log g = 4.39, derived a C abundance in agreement with our value, within uncertainties 4.6. HD 76909 The T eff = 5655 K we adopt is in very good agreement with Masana et al. (2006). Valenti & Fischer (2005) adopted an effective temperature about 100 K higher, a gravity 0.2 dex higher, and derive a metallicity about 0.1 dex higher. Our results are in very good agreement with those of Brewer et al. (2016) both for the stellar parameters and the detailed abundances (C, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni). Only our Na abundance is definitely larger. Petigura & Marcy (2011) adopted parameters compatible with ours (T eff = 5740 K and log g = 4.39) to derive a C abundance in perfect agreement with our value.
We remind the reader that we used only one SOPHIE spectrum for the optical analysis, the only one that has the sky.

HD 90681
For this star we adopted T eff = 5950 K and log g = 4.26. With very close stellar parameters (T eff = 5995 K and log g = 4.49), Petigura & Marcy (2011) derived the same C abundance.

HD 97645
We adopted T eff = 6127 K and log g = 4.10. With consistent stellar parameters (T eff = 6171 K and log g = 4.32), Petigura & Marcy (2011) derived a C abundance (A(C)=8.51) very close to our value of A(C)=8.48.

Granulation effects
We investigated whether the effects of granulation on abundance determinations are similar at optical and IR wavelengths. In order to assess this, we synthesised a number of fictitious Fe i and Fe ii lines of varying strength with excitation potentials between 3 and 5 eV at 500 nm (optical), 850 nm (NIR), and 1 600 nm (IR) using both 1D and 3D model atmospheres of solar-metallicity ([Fe/H] = +0.0) for stars with physical parameters representative of the present sample: T eff in the range 5 000 − 6 500 K and log g = 4.0 (cm s −2 ) plus the Sun (T eff = 5771 K and log g = 4.44). The 3D hydrodynamic and corresponding stationary 1D hydrostatic models were adopted from the staggergrid (Magic et al., 2013). We used these models in combination with the scate code (Hayek et al., 2011) to compute spectral line profiles under the assumption of local thermodynamic equilibrium (LTE). For the calculations with 1D models, we adopted a micro-turbulent parameter of ξ = 1.0 km s −1 . We derived the 3D−1D abundance corrections for the fictitious lines with varying strength by comparing the 3D and 1D Fe abundances required to match a given equivalent width (EW) value. We then quantified the differential granulation effects between different spectral regions by computing the 3D-1D abundance corrections of lines from the same ion at two different wavelengths, for example, optical and IR (but with identical excitation potential and identical normalised equivalent width EW/λ) and analysing the difference. We find the differential granulation effects for unsaturated Fe ii lines at optical and IR wavelengths to be generally small, ranging from about 0.0 dex at T eff = 5500 K to about 0.06 dex at T eff = 6500 K. Differential granulation effects for Fe i lines are typically 0.01 to 0.02 dex larger, but following a similar trend with effective temperature. For comparison, differential 3D−1D abundance corrections between NIR and IR lines show a flatter trend with values in the range 0.00 − 0.05 dex. Results are summarised in Fig. 3. Calculations performed with models from the CIFIST grid (Ludwig et al., 2009) show similar behaviour for the differential granulation effects (see also Caffau & Ludwig, 2007).

Detailed chemical investigation of the full sample
We were able to derive the abundance of 14 elements (C, Na, Mg, Al, Si, P, S, K, Ca, Ti, Cr, Fe, Ni, Sr) from the GIANO spectra. These are reported in Tables A.1-A.5. Some lines of Mn i are present in the range of GIANO (1297, 1328, 1329, 1521 nm), but the abundance we derive is not trustworthy, probably due to departure from LTE or granulation effects. Fig. 3. Differential granulation effects, i.e. differential 3D−1D abundance corrections δ ∆ 3D−1D A(Fe), between optical and IR and between NIR and IR wavelengths derived from moderatelyhigh-excitation unsaturated Fe lines for dwarf stars as a function of effective temperature. Filled and open symbols refer to Fe i and Fe ii lines, respectively. Circles and triangles represent lines with a lower-level excitation potential of 3 and 5 eV, respectively. Blue (green) symbols indicate differential corrections between optical and IR (NIR and IR) wavelengths. Abundance corrections are computed for fictitious lines with the same normalised equivalent width, EW/λ, and corresponding to a value of w λ = 4 pm at λ = 850 nm.

Phosphorus
Phosphorus has a single stable isotope 31 P, which is probably produced via neutron capture and the most likely site where it is produced is in the oxygen and neon burning shells during late stages of massive stars (see e.g. Caffau et al., 2011).
With this analysis, we added 40 points in the [P/Fe] versus [Fe/H] diagram. The P i lines we looked at are two of the lines of Mult. 1 at 1053.24 and 1058.44 nm (log gf of 0.20 and 0.48, respectively), which have already been used by Caffau et al. (2011Caffau et al. ( , 2016, and a line of Mult. 4 at 1008.70 nm (with log gf = 0.130). The two lines of Mult. 1 are shown in Fig. 4 in the case of HD 24040. The P abundance we derive from the 1008.70 nm line is on average larger (0.06 and 0.15 dex) than the one derived from the other two lines. This difference is not particularly large, the difference between the two lines of Mult. 1 being 0.09 dex.
In Fig.5, we compare these results with previous analyses. In particular, (Caffau et al., 2011(Caffau et al., , 2016) claimed a behaviour of P similar to an α-element, such as S. Unfortunately, the metallicity of this sample of stars is too high to be able to confirm this finding.

Sulphur
Sulphur is an α element produced by type-II supernovae. The fact that sulphur is a relatively volatile element, preventing it from being locked into dust grains in the interstellar medium, makes it a perfect element to compare α abundances in stars and in blue compact galaxies and damped Lyα systems. Details on S and an exhaustive summary of the results in the literature can be found in Duffau et al. (2017).
To derive A(S), we investigated the three S i lines of Mult. 3 at 1045 nm as we did in Caffau et al. (2016). We also looked at the triplet at 1541 nm, one single line in the multiplet at 1547 nm and the multiplet at 2261 nm (see Fig. 6 in the case of HD 24040). This multiplet is in a wavelength range contaminated by telluric absorptions that have to be removed when it is analysed.
We can compare the S abundance derived here with the metal-rich sample analysed by Ecuvillon et al. (2004) and with our previous investigations on sulphur. In Fig.7, we added the sample analysed here to the metal-rich range of the plot we showed in Duffau et al. (2017); the agreement with the previous investigations is very good. The fact that [S/Fe] analysed here is on average slightly higher than the result by Ecuvillon et al. (2004) could be explained by the fact that among the S i lines that we adopt to derive the S abundance, the lines of Mult. 3 at 1045 nm are affected by NLTE, which would slightly decrease A(S), while Ecuvillon et al. (2004) analysed the lines of Mult. 8 formed close to the LTE condition.

The other elements
We were able to find several atomic lines of C i in the GIANO wavelength range. This is extremely important because nitrogen can be derived from CN molecules once the C abundance is known. [C/Fe] is consistent with Galactic stars of slightly subsolar to super-solar metallicity ([C/Fe] = −0.04 ± 0.08) with just one star, HD 190360, with [C/Fe] > +0.3.
We used three lines of Na i to derive the Na abundances; according to our calculations, they all form close to LTE and therefore we do not expect large NLTE corrections (see in Fig. 8 the line at 1075 nm for HD 24040). The NLTE atomic model of sodium was presented by Korotin & Mishenina (1999) and then updated by Dobrovolskas et al. (2014). The updated sodium model currently consists of twenty energy levels of Na i and the ground level of Na ii. In total, 46 radiative transitions were taken into account for the calculation of the population of all levels. We carried out test calculations for a grid of stellar atmosphere models in the ranges: 5000 < T eff < 6000 K and 3.0 < log g < 4.5 with solar metallicity. For the IR sodium lines, NLTE effects lead to a slight increase in EW by 1 to 3%, while the EW in the optical lines increases by 10-18%. For all but one star we were able to derive the Na abundance with [Na/Fe] = +0.05 ± 0.17 in the sample. The large scatter is due to six stars with [Na/Fe] around +0.1 dex and several stars with negative [Na/Fe], down to about −0.2 dex.
The Mg i lines we selected to derive the Mg abundance are formed close to the LTE condition. In the case of Mg, we used the model atom from Mishenina et al. (2004), then updated by Cerniauskas et al. (2017). It consisted of 84 levels of Mg i, 12 levels of Mg ii, and the ground state of Mg iii. In the computation of departure coefficients, the radiative transitions between the first 59 levels of Mg i and the ground level of Mg ii were taken into account. We carried out test calculations within the same grid of star parameters as for sodium. It turned out that due to the influence of the effects of the deviation from the LTE, the EW of the magnesium lines used in this analysis vary by ±1%; that is, these lines are formed almost in the LTE. The average [Mg/Fe] we derived ( [Mg/Fe] = +0.02 ± 0.07) is the one expected for disc stars around solar metallicity.
For Al, we investigated two lines at 1087.59 and 1676.79 nm, which are formed close to LTE. Our Al atomic model is described in detail in Andrievsky et al. (2008). This model atom consists of 76 levels of Al i and 13 levels of Al ii. As shown by test calculations within that same grid star parameters as for the  other NLTE calculations, Al i lines 1087 and 1676 nm respond very poorly to NLTE effects. The effect on EW is within ±4%, depending on the parameters of the atmosphere. For all stars but three, A(Al) is based only on the Al i line at 1087.59 nm; this line is weak but is in a range that is not overly contaminated by telluric absorption (see Fig. 9). For the three stars for which also the line at 1676.79 nm is kept, the agreement in A(Al) derived from the two lines is good. For the 40 stars, we find [Al/Fe] = +0.10 ± 0.06, with some stars with [Al/Fe] of about +0.2 dex. Only one star, HD 136618, has a negative [Al/Fe]; this remains compatible within uncertainties with a solar ratio.
Five lines of K i have been analysed, at 1177. 28, 1177.50, 1243.56, 1252.55, and 1516.72 nm. The model atom of K was taken from Andrievsky et al. (2010) and consisted of 20 levels of K i and the ground level of K ii. In addition, another 15 levels of K i and seven levels of K ii were used to ensure particle number conservation. The total number of bound-bound radiative transitions taken into account was 62 (see Andrievsky et al., 2010, for further details). Infrared lines of potassium showed a strong susceptibility to the influence of NLTE effects. Therefore, due to deviations from the LTE, the EW of the lines used increased from 10% to 32%. We were able to derive A(K) for all but two stars, and [K/Fe] = +0.16 ± 0.06 for the 38 stars for which K has been derived. All the lines are sensitive to NLTE, and therefore the above-solar [K/Fe] value can easily be attributed to our LTE assumption.
A(Ca) from Ca i and Ca ii lines is in good agreement ( A(Ca i) − A(Ca ii) = −0.04±0.08) with a slightly higher value of A(Ca) when derived from the Ca ii lines, which can be explain by the fact that the Ca i lines we selected to derive the Ca abundance are close to LTE while some of the Ca ii lines are sensitive to NLTE and the A(Ca) derived in LTE overestimates the abundance. However, the A(Ca) abundance we derived for the sample stars is slightly high: from Ca i lines, [Ca/Fe] = +0.12 ± 0.05. In particular for four stars, HD 97645, HD 99492, HD 147231 and HD 190360, [Ca/Fe]   All three Sr ii lines already investigated in Caffau et al. (2016)  This could be due to an underestimated effective temperature, HD 190360 being the only star for which the temperature is derived from a single calibration. Increasing the temperature by 180 K to match for example the T eff of Mishenina et al. (2004) of 5606 K, the star no longer shows overabundance in C, P, and S, but still shows an overabundance in Ca. Also Bensby et al. (2003), with an effective temperature of 5490 K and a surface gravity of 4.23 (values very similar to the ones we adopted) found an overabundance with respect to Fe of several elements. In Bensby et al. (2014) however, with a T eff = 5572 K, this overabundance is no longer present. It is interesting to note that a linear regression of the mean T eff over the B − V T eff for the 37 stars for which five temperature determinations are available predicts T eff = 5593 K very close to the values of Mishenina et al. (2004) and Bensby et al. (2014).

Helium
The triplet line of neutral He arising from the high-excitation (19.72 eV) metastable level 2s 3 S 1 is observed in absorption in the Sun and in solar-type stars. It is of too high an excitation to be formed in the photosphere and therefore must be formed in the chromosphere and very likely excited by radiation from the corona. In spatially resolved spectra of the Sun it can be observed in absorption in active regions (plages). Following the activity cycle of the Sun, the EW measured in the flux-integrated spectrum varies depending on the number of active regions. Andretta & Giampapa (1995) suggested to use this line together with He i 587.6 nm to estimate the fraction of the stellar surface covered by active regions. In a pioneering work, Zirin (1982) obtained IR spectra at the coudé focus of the Palomar five-metre telescope using image cathode tubes and photographic plates and films as detectors for 455 stars. The EW of the He i 1083.0 nm line is well correlated with the soft X-ray luminosity and in spite of the technical advances even today, the "observational data on the λ10830 triplet line remains sparse" (verbatim from Smith, 2016). This line has not been studied in great detail because it requires high-resolution observations. The large spectral coverage provided by GIANO ensures that the spectral range is always observed, and it lies in a region that is contaminated by a few telluric lines that sometimes affect the He line. We inspected all the spectra of our sample and were able to detect this line for 16 stars. When the line was detected we measured its EW (Table 3). The typical uncertainty is 1 pm. For five stars in our sample we found X-ray measurements in the All-Sky Survey Faint source Catalogue (Voges et al., 2000). From these we derived log( f x / f bol ), also provided in Table 3. The apparent bolometric magnitude was derived from the Gaia G magnitude and the bolometric corrections of Bonifacio et al. (2018). The bolometric flux was then derived as recommended in IAU resolution 2015 B2 4 equation (4). These values are plotted in Fig.10 against the EW of the He i 1083.0 nm line and a very good correlation can be seen between the two quantities. GIANO observations offer a good opportunity to monitor the activity cycles of solar-type stars from the observation of this line. This line may also be used to monitor the activity of stars with candidate planets found in radial-velocity searches. The monitoring of the He i 1083.0 nm line should allow us to derive the activity cycle of the star. This may help decipher whether radial velocity variations are due to planets or stellar activity. A well known example is HD 166435 (Queloz et al., 2001).

Conclusions
In this work, we analysed 40 stars for which we provided detailed chemical abundances of 14 elements. These stars belong to the Galactic disc and the abundances we derive are consistent with the chemical behaviour expected for Galactic disc stars; they are also in line with previous investigations (e.g. see Mikolaitis et al., 2014;Smiljanic et al., 2014Smiljanic et al., , 2016Kordopatis et al., 2015;Duffau et al., 2017). For some of the stars investigated here there are no detailed chemical analyses in the literature.
We were also able to add 40 values for the P abundance in Galactic stars, which is still a rather poorly studied element. We also introduced a P i line of Mult. 4, the abundance of which is consistent with, though slightly higher than, the P abundance derived from the P i lines of Mult. 1, used in previous investigations (Caffau et al., 2011(Caffau et al., , 2016.
For a subsample of nine stars, we compared the abundances derived from the IR GIANO spectra and SOPHIE optical spectra. For eight chemical abundances common to both surveys agreement was found to be very good. Nevertheless, the line-toline scatter in the IR spectra is larger than in the optical spectra, probably due to a better knowledge on atomic data in the latter range and this is probably related to the longer history of use of optical ranges.
Both the good agreement on the abundances derived from IR and optical spectra and the good agreement with the literature are indicators that the abundances we derive from GIANO spectra are reliable. For the star HD 98736, we analysed only < 5 d a angular distance ′′ between the X-ray source and the star four orders of GIANO (47-50) which correspond to the highresolution mode of the H-band for MOONS. In this case we can derive abundances for seven elements (C, Si, S, Ti, Cr, Fe and Ni). We analysed this range also for the star HD 24040, the spectral range shown in Fig. 1, and we derive the abundance for C, Si, S, K, Ti, Cr, Fe, Ni. Clearly, by using only four orders of the GIANO spectra, the uncertainties in the abundances are much larger because the number of lines available is much lower than when analysing the complete GIANO range; for the majority of the elements abundances are based on a single line. On top of that, because of the lower resolution that MOONS will have with respect to GIANO, some of these elements will be difficult to be derived by MOONS, in particular Ti, Cr, and Ni.
For several stars in the sample we were able to detect the He i line at 1083.0 nm and confirm that its EW is strongly correlated with the logarithmic ratio of X-ray flux to bolometric flux and, hence, stellar activity. Since He measurements in dwarf stars are still limited in number, GIANO provides a great opportunity in regards to this topic.
GIANO provides us with a great chance to investigate Galactic stars and to derive their detailed chemical abundances. From GIANO spectra we could derive the abundances of some elements (C, P, K, Sr) that could not be derived from the SOPHIE optical spectra. GIANO's high resolution and large spectral range are fundamental to deriving trustable abundances for several elements. GIANO is a single-object spectrograph meaning that stars can only be observed one at a time, but in the next years its archive will fill and the number of stars available will be larger and larger.