The scope of the project "A PSF-based Approach to TESS High Quality data Of Stellar clusters" (PATHOS) is the extraction and analysis of high-precision light curves of stars in stellar clusters and young associations for the identification of candidate exoplanets and variable stars. The cutting-edge tools used in this project allow us to measure the real flux of stars in dense fields, minimising the effects due to contamination by neighbour sources. We extracted about 200 000 light curves of stars in 645 open clusters located in the southern ecliptic hemisphere and observed by TESS during the first year of its mission. We searched for transiting signals and we found 33 objects of interest, 11 of them are strong candidate exoplanets. Because of the limited S/N, we did not find any Earth or super-Earth. We identified two Neptune-size planets orbiting stars with $R_{\star}<1.5\,R_{\odot}$, implying a frequency $f_{\star}=1.34 \pm 0.95\,\%$, consistent with the frequency around field stars. The 7 Jupiter candidates around stars with $R_{\star}<\,1.5R_{\odot}$ imply a frequency $f_{\star}=0.19\pm 0.07\,\%$, smaller than in the field. A more complete estimate of the survey completeness and false positive rate is needed to confirm these results. Light curves used in this work will be made available to the astronomical community on the Mikulski Archive for Space Telescope under the project PATHOS.

We present the first detection of atomic emission lines from the atmosphere of an exoplanet. We detect neutral iron lines from the day-side of KELT-9b (Teq $\sim$ 4, 000 K). We combined thousands of spectrally resolved lines observed during one night with the HARPS-N spectrograph (R $\sim$ 115, 000), mounted at the Telescopio Nazionale Galileo. We introduce a novel statistical approach to extract the planetary parameters from the binary mask cross-correlation analysis. We also adapt the concept of contribution function to the context of high spectral resolution observations, to identify the location in the planetary atmosphere where the detected emission originates. The average planetary line profile intersected by a stellar G2 binary mask was found in emission with a contrast of 84 $\pm$ 14 ppm relative to the planetary plus stellar continuum (40 $\pm$ 5$\%$ relative to the planetary continuum only). This result unambiguously indicates the presence of an atmospheric thermal inversion. Finally, assuming a modelled temperature profile previously published (Lothringer et al. 2018), we show that an iron abundance consistent with a few times the stellar value explains the data well. In this scenario, the iron emission originates at the $10^{-3}$-$10^{-5}$ bar level.

The TESS mission will survey ~85 % of the sky, giving us the opportunity of extracting high-precision light curves of millions of stars, including stellar cluster members. In this work, we present our project "A PSF-based Approach to TESS High quality data Of Stellar clusters" (PATHOS), aimed at searching and characterise candidate exoplanets and variable stars in stellar clusters using our innovative method for the extraction of high-precision light curves of stars located in crowded environments. Our technique of light-curve extraction involves the use of empirical Point Spread Functions (PSFs), an input catalogue and neighbour-subtraction. The PSF-based approach allows us to minimise the dilution effects in crowded environments and to extract high-precision photometry for stars in the faint regime (G>13). For this pilot project, we extracted, corrected, and analysed the light curves of 16641 stars located in a dense region centred on the globular cluster 47 Tuc. We were able to reach the TESS magnitude T~16.5 with a photometric precision of ~1 % on the 6.5-hour timescale; in the bright regime we were able to detect transits with depth of ~34 parts per million. We searched for variables and candidate transiting exoplanets. Our pipeline detected one planetary candidate orbiting a main sequence star in the Galactic field. We analysed the period-luminosity distribution for red-giant stars of 47 Tuc and the eclipsing binaries in the field. Light curves are uploaded on the Mikulski Archive for Space Telescopes under the project PATHOS.

The origin of the observed diversity of planetary system architectures is one of the main topic of the exoplanetary research. The detection of a statistically significant sample of planets around young stars allows us to study the early stages of planet formation and evolution, but only a handful of them is known so far. In this regard, a considerable contribution is expected from the NASA TESS satellite, which is now performing a survey of $\sim 85 \%$ of the sky to search for short-period transiting planets In its first month of operations, TESS found a planet candidate with an orbital period of 8.14 days around a member of the Tuc-Hor young association ($\sim$ 40 Myr), the G6V main component of the binary system DS\,Tuc. If confirmed, it would be the first transiting planet around a young star suitable for radial velocity and/or atmospheric characterization. We aim to validate the planetary nature of this companion and to measure its orbital and physical parameters. We obtain accurate planet parameters by coupling an independent reprocessing of the TESS light curve with improved stellar parameters and the dilution caused by the binary companion; we analyse high precision archival radial velocities to impose an upper limit of about 0.1 M$_{\rm Jup}$ on the planet mass; we finally rule out the presence of external companions beyond 40 au with adaptive optics images. We confirm the presence of a young, giant (${\rm R} = 0.50$ R$_{\rm Jup}$) planet having a not negligible possibility to be inflated (theoretical mass $\lesssim 20$ M$_{\oplus}$) around DS\,Tuc~A. We discuss the feasibility of mass determination, Rossiter-McLaughlin analysis and atmosphere characterization, allowed by the brightness of the star.

We investigated the discrepancy between planetary mass determination using the transit timing variations (TTVs) and radial velocities (RVs), by analysing the multi-planet system Kepler-9. Despite being the first system characterised with TTVs, there are several discrepant solutions in the literature, with those reporting lower planetary densities being apparently in disagreement with high-precision RV observations. To resolve this, we gathered HARPS-N RVs at epochs that maximised the difference between the predicted RV curves from discrepant solutions in the literature. We also re-analysed the full Kepler data-set and performed a dynamical fit, within a Bayesian framework, using the newly derived central and duration times of the transits. We compared these results with the RV data and found that our solution better describes the RV observations, despite the masses of the planets being nearly half that presented in the discovery paper. We therefore confirm that the TTV method can provide mass determinations that agree with those determined using high-precision RVs. The low densities of the planets place them in the scarcely populated region of the super-Neptunes / inflated sub-Saturns in the mass-radius diagram.

We report a detailed characterization of the Kepler-19 system. This star was previously known to host a transiting planet with a period of 9.29 days, a radius of 2.2 R$_\oplus$ and an upper limit on the mass of 20 M$_\oplus$. The presence of a second, non-transiting planet was inferred from the transit time variations (TTVs) of Kepler-19b, over 8 quarters of Kepler photometry, although neither mass nor period could be determined. By combining new TTVs measurements from all the Kepler quarters and 91 high-precision radial velocities obtained with the HARPS-N spectrograph, we measured through dynamical simulations a mass of $8.4 \pm 1.6$ M$_\oplus$ for Kepler-19b. From the same data, assuming system coplanarity, we determined an orbital period of 28.7 days and a mass of $13.1 \pm 2.7$ M$_\oplus$ for Kepler-19c and discovered a Neptune-like planet with a mass of $20.3 \pm 3.4$ M$_\oplus$ on a 63 days orbit. By comparing dynamical simulations with non-interacting Keplerian orbits, we concluded that neglecting interactions between planets may lead to systematic errors that could hamper the precision in the orbital parameters when the dataset spans several years. With a density of $4.32 \pm 0.87$ g cm$^{-3}$ ($0.78 \pm 0.16$ $\rho_\oplus$) Kepler-19b belongs to the group of planets with a rocky core and a significant fraction of volatiles, in opposition to low-density planets characterized by transit-time variations only and the increasing number of rocky planets with Earth-like density. Kepler-19 joins the small number of systems that reconcile transit timing variation and radial velocity measurements.

We present FIES@NOT, HARPS-N@TNG, and HARPS@ESO-3.6m radial velocity follow-up observations of K2-19, a compact planetary system hosting three planets, of which the two larger ones, namely K2-19b and K2-19c, are close to the 3:2 mean motion resonance. An analysis considering only the radial velocity measurements detects K2-19b, the largest and most massive planet in the system, with a mass of $54.8\pm7.5$~M${_\oplus}$ and provides a marginal detection of K2-19c, with a mass of M$_\mathrm{c}$=$5.9^{+7.6}_{-4.3}$ M$_\oplus$. We also used the TRADES code to simultaneously model both our RV measurements and the existing transit-timing measurements. We derived a mass of $54.4\pm8.9$~M${_\oplus}$ for K2-19b and of $7.5^{+3.0}_{-1.4}$~M${_\oplus}$ for K2-19c. A prior K2-19b mass estimated by Barros et al. 2015, based principally on a photodynamical analysis of K2-19's light-curve, is consistent with both analysis, our combined TTV and RV analysis, and with our analysis based purely on RV measurements. Differences remain mainly in the errors of the more lightweight planet, driven likely by the limited precision of the RV measurements and possibly some yet unrecognized systematics.

In the framework of the GAPS project, we are conducting an observational programme aimed at the determination of the orbital obliquity of known transiting exoplanets. The targets are selected to probe the obliquity against a wide range of stellar and planetary physical parameters. We exploit high-precision radial velocity (RV) measurements, delivered by the HARPS-N spectrograph at the 3.6m Telescopio Nazionale Galileo, to measure the Rossiter-McLaughlin (RM) effect in RV time-series bracketing planet transits, and to refine the orbital parameters determinations with out-of-transit RV data. We also analyse new transit light curves obtained with several 1-2m class telescopes to better constrain the physical fundamental parameters of the planets and parent stars. We report here on new transit spectroscopic observations for three very massive close-in giant planets: WASP43b, HATP20b and Qatar2b orbiting dwarf K-type stars with effective temperature well below 5000K. We find lambda = 3.5pm6.8 deg for WASP43b and lambda = -8.0pm6.9 deg for HATP20b, while for Qatar2, our faintest target, the RM effect is only marginally detected, though our best-fit value lambda = 15pm20 deg is in agreement with a previous determination. In combination with stellar rotational periods derived photometrically, we estimate the true spin-orbit angle, finding that WASP43b is aligned while the orbit of HATP20b presents a small but significant obliquity (Psi=36 _{-12}^{+10} deg). By analyzing the CaII H&K chromospheric emission lines for HATP20 and WASP43, we find evidence for an enhanced level of stellar activity which is possibly induced by star-planet interactions.

In the third paper of this series we continue the exploitation of Kepler/K2 data in dense stellar fields using our PSF-based method. This work is focused on a ~720-arcmin^2 region centred on the Solar-metallicity and Solar-age open cluster M 67. We extracted light curves for all detectable sources in the Kepler channels 13 and 14, adopting our technique based on the usage of a high-angular-resolution input catalogue and target-neighbour subtraction. We detrended light curves for systematic errors, and searched for variables and exoplanets using several tools. We found 451 variables, of which 299 are new detection. Three planetary candidates were detected by our pipeline in this field. Raw and detrended light curves, catalogues, and K2 stacked images used in this work will be released to the community.

In this work we keep pushing K2 data to a high photometric precision, close to that of the Kepler main mission, using a PSF-based, neighbour-subtraction technique, which also overcome the dilution effects in crowded environments. We analyse the open cluster M 44 (NGC 2632), observed during the K2 Campaign 5, and extract light curves of stars imaged on module 14, where most of the cluster lies. We present two candidate exoplanets hosted by cluster members and five by field stars. As a by-product of our investigation, we find 1680 eclipsing binaries and variable stars, 1071 of which are new discoveries. Among them, we report the presence of a heartbeat binary star. Together with this work, we release to the community a catalogue with the variable stars and the candidate exoplanets found, as well as all our raw and detrended light curves.

We confirm the planetary nature of Kepler-539b (aka Kepler object of interest K00372.01), a giant transiting exoplanet orbiting a solar-analogue G2 V star. The mass of Kepler-539b was accurately derived thanks to a series of precise radial velocity measurements obtained with the CAFE spectrograph mounted on the CAHA 2.2m telescope. A simultaneous fit of the radial-velocity data and Kepler photometry revealed that Kepler-539b is a dense Jupiter-like planet with a mass of Mp = 0.97 Mjup and a radius of Rp = 0.747 Rjup, making a complete circular revolution around its parent star in 125.6 days. The semi-major axis of the orbit is roughly 0.5 au, implying that the planet is at roughly 0.45 au from the habitable zone. By analysing the mid-transit times of the 12 transit events of Kepler-539b recorded by the Kepler spacecraft, we found a clear modulated transit time variation (TTV), which is attributable to the presence of a planet c in a wider orbit. The few timings available do not allow us to precisely estimate the properties of Kepler-539c and our analysis suggests that it has a mass between 1.2 and 3.6 Mjup, revolving on a very eccentric orbit (0.4<e<0.6) with a period larger than 1000 days. The high eccentricity of planet c is the probable cause of the TTV modulation of planet b. The analysis of the CAFE spectra revealed a relatively high photospheric lithium content, A(Li)=2.48 dex, which, together with both a gyrochronological and isochronal analysis, suggests that the parent star is relatively young.

Open cluster (OC) stars share the same age and metallicity, and, in general, their age and mass can be estimated with higher precision than for field stars. For this reason, OCs are considered an important laboratory to study the relation between the physical properties of the planets and those of their host stars, and the evolution of planetary systems. We started an observational campaign within the GAPS collaboration to search for and characterize planets in OCs We monitored the Praesepe member Pr0211 to improve the eccentricity of the Hot-Jupiter (HJ) already known to orbit this star and search for additional planets. An eccentric orbit for the HJ would support a planet-planet scattering process after its formation. From 2012 to 2015, we collected 70 radial velocity (RV) measurements with HARPS-N and 36 with TRES of Pr0211. Simultaneous photometric observations were carried out with the robotic STELLA telescope in order to characterize the stellar activity. We discovered a long-term trend in the RV residuals that we show to be due to the presence of a second, massive, outer planet. Orbital parameters for the two planets are derived by simultaneously fitting RVs and photometric light curves, with the activity signal modelled as a series of sinusoids at the rotational period of the star and its harmonics. We confirm that Pr0211b has a nearly circular orbit ($e = 0.02 \pm 0.01$), with an improvement of a factor two with respect to the previous determination of its eccentricity, and estimate that Pr0211c has a mass $M_p\sin i = 7.9 \pm 0.2 M_J$, a period $P>$3500 days and a very eccentric orbit ($e>$0.60). Such peculiar systems may be typical of open clusters if the planet-planet scattering phase leading to the formation of HJs is caused by stellar encounters rather than unstable primordial orbits. Pr0211 is the first multi-planet system discovered around an OC star. (abridged)

In this paper we continue the release of high-level data products from the multiyear photometric survey collected at the 67/92 cm Schmidt Telescope in Asiago. The primary goal of the survey is to discover and to characterise variable objects and exoplanetary transits in four fields containing five nearby open clusters spanning a broad range of ages. This second paper releases a photometric catalogue, in five photometric bands, of the Solar-age, Solar-metallicity open cluster M 67 (NGC 2682). Proper motions are derived comparing the positions observed in 2013 at the Asiago's Schmidt Telescope with those extracted from WFI@2.2m MPG/ESO images in 2000. We also analyse the variable sources within M 67. We detected 68 variables, 43 of which are new detection. Variable periods and proper-motion memberships of a large majority of sources in our catalogue are improved with respect to previous releases. The entire catalogue will be available in electronic format. Besides the general interest on an improved catalogue, this work will be particularly useful because of: (1) the imminent release of Kepler/K2 Campaign-5 data of this cluster, for which our catalogue will provide an excellent, high spatial resolution input list, and (2) characterisation of the M 67 stars which are targets of intense HARPS and HARPS-N radial-velocity surveys for planet search.

XO-2 is the first confirmed wide stellar binary system where the almost twin components XO-2N and XO-2S have planets. This stimulated a detailed characterization study of the stellar and planetary components based on new observations. We collected high-resolution spectra with the HARPS-N spectrograph and multi-band light curves. Spectral analysis led to an accurate determination of the stellar atmospheric parameters and characterization of the stellar activity. We collected 14 transit light curves of XO-2Nb used to improve the transit parameters. Photometry provided accurate magnitude differences between the stars and a measure of their rotation periods. The iron abundance of XO-2N was found to be +0.054 dex greater, within more than 3-sigma, than that of XO-2S. We confirm a long-term variation in the radial velocities of XO-2N, and we detected a turn-over with respect to previous measurements. We suggest the presence of a second massive companion in an outer orbit or the stellar activity cycle as possible causes of the observed acceleration. The latter explanation seems more plausible with the present dataset. We obtained an accurate value of the projected spin-orbit angle for the XO-2N system (lambda=7+/-11 degrees), and estimated the real 3-D spin-orbit angle (psi=27 +12/-27 degrees). We measured the XO-2 rotation periods, and found a value of P=41.6 days in the case of XO-2N, in excellent agreement with the predictions. The period of XO-2S appears shorter, with an ambiguity between 26 and 34.5 days that we cannot solve with the present dataset alone. XO-2N appears to be more active than the companion, and this could be due to the fact that we sampled different phases of their activity cycle, or to an interaction between XO-2N and its hot Jupiter that we could not confirm.

We present a multi-year survey aimed at collecting (1) high-precision (~5 milli-mag), (2) fast-cadence (~3 min), and (3) relatively long duration (~10 days) multi-band photometric series. The goal of the survey is to discover and characterize efficiently variable objects and exoplanetary transits in four fields containing five nearby open clusters spanning a broad range of ages. More in detail, our project will (1) constitute a preparatory survey for HARPS-N@TNG, which will be used for spectroscopic follow-up of any target of interest that this survey discovers or characterizes, (2) measure rotational periods and estimate the activity level of targets we are already monitoring with HARPS and HARPS-N for exoplanet transit search, and (3) long term characterization of selected targets of interest in open clusters within the planned K2 fields. In this first paper we give an overview of the project, and report on the variability of objects within the first of our selected fields, which contains two open clusters: M 35 and NGC 2158. We detect 519 variable objects, 273 of which are new discoveries, while the periods of most of the previously known variables are considerably improved.

Aims. With the purpose of determining the orbital parameters of exoplanetary systems from observational data, we have developed a software, named TRADES (TRAnsits and Dynamics of Exoplanetary Systems), to simultaneously fit observed radial velocities and transit times data. Methods. We implemented a dynamical simulator for N-body systems, which also fits the available data during the orbital integration and determines the best combination of the orbital parameters using grid search, $\chi^2$ minimization, genetic algorithms, particle swarm optimization, and bootstrap analysis. Results. To validate TRADES, we tested the code on a synthetic three-body system and on two real systems discovered by the Kepler mission: Kepler-9 and Kepler-11. These systems are good benchmarks to test multiple exoplanet systems showing transit time variations (TTVs) due to the gravitational interaction among planets. We have found that orbital parameters of Kepler-11 planets agree well with the values proposed in the discovery paper and with a a recent work from the same authors. We analyzed the first three quarters of Kepler-9 system and found parameters in partial agreement with discovery paper. Analyzing transit times (T0s) covering 12 quarters of Kepler data, that we have found a new best-fit solution. This solution outputs masses that are about 55% of the values proposed in the discovery paper; this leads to a reduced semi-amplitude of the radial velocities of about 12.80 m/s.

We present four new light curves of transiting exoplanets WASP-1b and HAT-P-20b, observed within the TASTE (The Asiago Search for Transit timing variations of Exoplanets) project. We re-analyzed light curves from the literature in a homogeneous way, calculating a refined ephemeris and orbital-physical parameters for both objects. WASP-1b does not show any significant Transit Timing Variation signal at the 120 s-level. As for HAT-P-20b, we detected a deviation from our re-estimated linear ephemeris that could be ascribed to the presence of a perturber or, more probably, to a previously unnoticed high level of stellar activity. The rotational period of HAT-P-20 A we obtained from archival data (P_rot ~ 14.5 days), combined with its optical variability and strong emission of CaII H&K lines, is consistent with a young stellar age (< 1 Gyr) and support the hypothesis that stellar activity may be responsible of the measured deviations of the transit times.

A long-term multi-purpose observational programme has started with HARPS-N@TNG aimed to characterise the global architectural properties of exoplanetary systems. In this first paper we fully characterise the transiting system Qatar-1. We exploit HARPS-N high-precision radial velocity measurements obtained during a transit to measure the Rossiter-McLaughlin effect in the Qatar-1 system, and out-of-transit measurements to redetermine the spectroscopic orbit. New photometric transit light-curves are analysed and a spectroscopic characterisation of the host star atmospheric parameters is performed based on various methods (line equivalent width ratios, spectral synthesis, spectral energy distribution). We achieved a significant improvement in the accuracy of the orbital parameters and derived the spin-orbit alignment of the system; this information, combined with the spectroscopic determination of the host star properties, allows us to derive the fundamental physical parameters for star and planet (masses and radii). The orbital solution for the Qatar-1 system is consistent with a circular orbit and the system presents a sky-projected obliquity of lambda = -8.4+-7.1 deg. The planet, with a mass of 1.33+-0.05 M_J, is found to be significantly more massive than previously reported. The host star is confirmed to be metal-rich ([Fe/H]= 0.20+-0.10) and slowly rotating (vsinI = 1.7+-0.3 km/s), though moderately active, as indicated by strong chromospheric emission in the Ca II H&K line cores (logR'_HK about -4.60). The system is well aligned and fits well within the general lambda vs Teff trend. We definitely rule out any significant orbital eccentricity. The evolutionary status of the system is inferred based on gyrochronology, and the present orbital configuration and timescale for orbital decay are discussed in terms of star-planet tidal interactions.

The TASTE project is searching for low-mass planets with the Transit Timing Variation (TTV) technique, by gathering high-precision, short-cadence light curves for a selected sample of transiting exoplanets. It has been claimed that the "hot Jupiter" WASP-3b could be perturbed by a second planet. Presenting eleven new light curves (secured at the IAC80 and UDEM telescopes) and re-analyzing thirty-eight archival light curves in a homogeneous way, we show that new data do not confirm the previously claimed TTV signal. However, we bring evidence that measurements are not consistent with a constant orbital period, though no significant periodicity can be detected. Additional dynamical modeling and follow-up observations are planned to constrain the properties of the perturber or to put upper limits to it. We provide a refined ephemeris for WASP-3b and improved orbital/physical parameters. A contact eclipsing binary, serendipitously discovered among field stars, is reported here for the first time.

TASTE (The Asiago Search for Transit timing variations of Exoplanets) project is collecting high-precision, short-cadence light curves for a selected sample of transiting exoplanets. It has been claimed that the hot jupiter HAT-P-13b suddenly deviated from a linear ephemeris by $\sim 20$ min, implying that there is a perturber in the system. Using five new transits, we discuss the plausibility of this transit time variation (TTV), and show that a periodic signal should not be excluded. More follow-up observations are required to constrain the mass and the orbit of the hypothetical perturber.