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 constructed a Hubble Space Telescope (HST) astro-photometric catalog of the central region of the Galactic globular cluster NGC 1261. This catalog, complemented with Gaia DR2 data sampling the external regions, has been used to estimate the structural parameters of the system (i.e., core, half-mass, tidal radii and concentration) from its resolved star density profile. We computed high-precision proper motions thanks to multi-epoch HST data and derived the cluster velocity dispersion profile in the plane of the sky for the innermost region, finding that the system is isotropic. The combination with line-of-sight information collected from spectroscopy in the external regions provided us with the cluster velocity dispersion profile along the entire radial extension. We also measured the absolute proper motion of NGC 1261 using a few background galaxies as a reference. The radial distribution of the Blue Straggler Star population shows that the cluster is in a low/intermediate phase of dynamical evolution.

A number of scenarios for the formation of multiple populations in globular clusters (GCs) predict that second generation (2G) stars form in a compact and dense subsystem embedded in a more extended first-generation (1G) system. If these scenarios are accurate, a consequence of the denser 2G formation environment is that 2G binaries should be more significantly affected by stellar interactions and disrupted at a larger rate than 1G binaries. The fractions and properties of binary stars can thus provide a dynamical fingerprint of the formation epoch of multiple-population GCs and their subsequent dynamical evolution. We investigate the connection between binaries and multiple populations in five GCs, NGC 288, NGC 6121 (M 4), NGC 6352, NGC 6362, and NGC 6838 (M 71). To do this, we introduce a new method based on the comparison of Hubble Space Telescope observations of binaries in the F275W, F336W, F438W, F606W and F814W filters with a large number of simulated binaries. In the inner regions probed by our data we do not find large differences between the local 1G and the 2G binary incidences in four of the studied clusters, the only exception being M 4 where the 1G binary incidence is about three times larger than the 2G incidence. The results found are in general agreement with the results of simulations predicting significant differences in the global 1G and 2G incidences and in the local values in the clusters' outer regions but similar incidences in the inner regions. The significant difference found in M 4 is consistent with simulations with a larger fraction of wider binaries. Our analysis also provides the first evidence of mixed (1G-2G) binaries, a population predicted by numerical simulations to form in a cluster's inner regions as a result of stellar encounters during which one component of a binary is replaced by a star of a different population.

In the present work we analyzed seven globular clusters selected from their location in the Galactic bulge and with metallicity values in the range $-1.30\lesssim\rm{[Fe/H]}\lesssim-0.50$. The aim of this work is first to derive cluster ages assuming single stellar populations, and secondly, to identify the stars from first (1G) and second generations (2G) from the main sequence, subgiant and red giant branches, and to derive their age differences. Based on a combination of UV and optical filters used in this project, we apply the Gaussian mixture models to distinguish the multiple stellar populations. Applying statistical isochrone fitting, we derive self-consistent ages, distances, metallicities, and reddening values for the sample clusters. An average of $12.3\pm0.4$ Gyr was obtained both using Dartmouth and BaSTI (accounting atomic diffusion effects) isochrones, without a clear distinction between the moderately metal-poor and the more metal-rich bulge clusters, except for NGC 6717 and the inner halo NGC 6362 with $\sim 13.5$ Gyr. We derived a weighted mean age difference between the multiple populations hosted by each globular cluster of $41\pm170$ Myr adopting canonical He abundances; whereas for higher He in 2G stars, this difference reduces to $17\pm170$ Myr, but with individual uncertainties of $500$ Myr.

The Wide-Field InfraRed Space Telescope (WFIRST) will be capable of delivering precise astrometry for faint sources over the enormous field of view of its main camera, the Wide-Field Imager (WFI). This unprecedented combination will be transformative for the many scientific questions that require precise positions, distances, and velocities of stars. We describe the expectations for the astrometric precision of the WFIRST WFI in different scenarios, illustrate how a broad range of science cases will see significant advances with such data, and identify aspects of WFIRST's design where small adjustments could greatly improve its power as an astrometric instrument.

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.

We present a detailed analysis of HARPS-N radial velocity observations of K2-100, a young and active star in the Praesepe cluster, which hosts a transiting planet with a period of 1.7 days. We model the activity-induced radial velocity variations of the host star with a multi-dimensional Gaussian Process framework and detect a planetary signal of $10.6 \pm 3.0 {\rm m\,s^{-1}}$, which matches the transit ephemeris, and translates to a planet mass of $21.8 \pm 6.2 M_\oplus$. We perform a suite of validation tests to confirm that our detected signal is genuine. This is the first mass measurement for a transiting planet in a young open cluster. The relatively low density of the planet, $2.04^{+0.66}_{-0.61} {\rm g\,cm^{-3}}$, implies that K2-100b retains a significant volatile envelope. We estimate that the planet is losing its atmosphere at a rate of $10^{11}-10^{12}\,{\rm g\,s^{-1}}$ due to the high level of radiation it receives from its host star.

Globular clusters are collisional systems, meaning that the stars inside them interact on timescales much shorter than the age of the Universe. These frequent interactions transfer energy between stars and set up observable trends that tell the story of a cluster's evolution. This contribution focuses on what we can learn by studying velocity anisotropy and energy equipartition in globular clusters with Hubble Space Telescope proper motions.

We report on the white dwarf cooling sequence of the old globular cluster NGC 6752, which is chemically complex and hosts a blue horizontal branch. This is one of the last globular cluster white dwarf (WD) cooling sequences accessible to imaging by the Hubble Space Telescope. Our photometry and completeness tests show that we have reached the peak of the luminosity function of the WD cooling sequence, at a magnitude m_F606W=29.4+/-0.1, which is consistent with a formal age of ~14Gyr. This age is also consistent with the age from fits to the main-sequence turnoff (13-14Gyr), reinforcing our conclusion that we observe the expected accumulation of white dwarfs along the cooling sequence.

We have constructed the Spectral Energy Distributions (SEDs) of a sample of Blue Straggler Stars (BSSs) in the core of the globular cluster 47 Tucanae, taking advantage of the large set of high resolution images, ranging from the ultraviolet to the near infrared, obtained with the ACS/HRC camera of the Hubble Space Telescope. Our final BSS sample consists of 22 objects, spanning the whole color and magnitude extension of the BSS sequence in 47 Tucanae. We fitted the BSS broadband SEDs with models to derive temperature, surface gravity, radius, luminosity and mass. We show that BSSs indeed define a mass sequence, where the mass increases for increasing luminosity. Interestingly, the BSS masses estimates from the SED fitting turn out to be comparable to those derived from the projection of the stellar position in the color-magnitude diagram onto standard star evolutionary tracks. We compare our results with previous, direct mass estimates of a few BSSs in 47 Tucanae. We also find a couple of supermassive BSS candidates, i.e., BSSs with a mass larger than twice the turn-off mass, the formation of which must have involved more than two progenitors.

We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observed only hundreds of millions of years after the Big Bang, and all theoretically viable paths to making supermassive BHs require a stage of "intermediate" mass. However, no BHs have yet been reliably detected in the 100-10}^5 M_sun mass range. Uncovering these intermediate-mass BHs of 10^3-10^5 M_sun is within reach in the coming decade. In this white paper we highlight the crucial role that 30-m class telescopes will play in dynamically detecting intermediate-mass black holes, should they exist.

As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion.

The outer regions of globular clusters can enable us to answer many fundamental questions concerning issues ranging from the formation and evolution of clusters and their multiple stellar populations to the study of stars near and beyond the hydrogen-burning limit and to the dynamics of the Milky Way. The outskirts of globular clusters are still uncharted territories observationally. A very efficient way to explore them is through high-precision proper motions of low-mass stars over a large field of view. The Wide Field InfraRed Survey Telescope (WFIRST) combines all these characteristics in a single telescope, making it the best observational tool to uncover the wealth of information contained in the clusters' outermost regions.

We present the analysis of the radial distributions and kinematic properties of the multiple stellar populations (mPOPs) hosted in the globular cluster (GC) NGC 6352 as part of the Hubble Space Telescope (HST) "UV Legacy Survey of Galactic Globular Clusters" program. NGC 6352 is one of the few GCs for which the mPOP tagging in appropriate color-magnitude diagrams is clear in all evolutionary sequences. We computed high-precision stellar proper motions for the stars from the cluster's core out to 75 arcsec ($\sim$1.5 core radii, or $\sim$0.6 half-light radii). We find that, in the region explored, first- and second-generation stars share the same radial distribution and kinematic properties. Velocity dispersions, anisotropy radial profiles, differential rotation, and level of energy equipartition, all suggest that NGC 6352 is probably in an advanced evolutionary stage, and any possible difference in the structural and kinematic properties of its mPOPs have been erased by dynamical processes in the core of the cluster. We also provide an estimate of the mass of blue stragglers and of main-sequence binaries through kinematics alone. In general, in order to build a complete dynamical picture of this and other GCs, it will be essential to extend the analyses presented in this paper to the GCs' outer regions where some memories of the initial differences in the mPOP properties, and those imprinted by dynamical processes, might still be present.

As part of a large Hubble Space Telescope investigation aiming at reaching the faintest stars in the Galactic globular cluster NGC 6752, an ACS/WFC field was the subject of deep optical observations reaching magnitudes as faint as V~30. In this field we report the discovery of Bedin I a dwarf spheroidal galaxy too faint and too close to the core of NGC 6752 for detection in earlier surveys. As it is of broad interest to complete the census of galaxies in the local Universe, in this Letter we provide the position of this new object along with preliminary assessments of its main param eters. Assuming the same reddening as for NGC 6752, we estimate a distance modulus of (m-M)_0 = 29.70+/-0.13 from the observed red giant branch, i.e., 8.7 (+0.5 -0.7) Mpc, and size of ~840x340pc, about 1/5 the size of the LMC. A comparison of the observed colour-magnitude diagram with synthetic counterparts that account for the galaxy distance modulus, reddening and photometric errors, suggests the presence of an old (~13Gyr) and metal poor ([Fe/H]~-1.3) population. This object is most likely a relatively isolated satellite dwarf spheroidal galaxy of the nearby great spiral NGC 6744, or potentially the most distant isolated dwarf spheroidal known with a secure distance.

Historically, multiple populations in Globular Clusters (GCs) have been mostly studied from ultraviolet and optical filters down to stars that are more massive than ~0.6 solar masses. Here we exploit deep near-infrared (NIR) photometry from the Hubble Space Telescope to investigate multiple populations among M-dwarfs in the GC NGC6752. We discovered that the three main populations (A, B and C), previously observed in the brightest part of the color-magnitude diagram, define three distinct sequences that run from the main-sequence (MS) knee towards the bottom of the MS (~0.15 solar masses). These results, together with similar findings on NGC2808, M4, and omega Centauri, demonstrate that multiple sequences of M-dwarfs are common features of the color-magnitude diagrams of GCs. The three sequences of low-mass stars in NGC6752 are consistent with stellar populations with different oxygen abundances. The range of [O/Fe] needed to reproduce the NIR CMD of NGC6752 is similar to the oxygen spread inferred from high-resolution spectroscopy of red-giant branch (RGB) stars. The relative numbers of stars in the three populations of M-dwarfs are similar to those derived among RGB and MS stars more massive than ~0.6 solar masses. As a consequence, the evidence that the properties of multiple populations do not depend on stellar mass is a constraint for the formation scenarios.

HP$\,$1 is an $\alpha$-enhanced and moderately metal-poor bulge globular cluster with a blue horizontal branch. These combined characteristics make it a probable relic of the early star formation in the innermost Galactic regions. Here we present a detailed analysis of a deep near-infrared (NIR) photometry of HP$\,$1 obtained with the NIR GSAOI+GeMS camera at the Gemini-South telescope. $J$ and $K_{\rm S}$ images were collected with an exquisite spatial resolution (FWHM $\sim 0.1$ arcsec), reaching stars at two magnitudes below the MSTO. We combine our GSAOI data with archival F606W-filter $HST$ ACS/WFC images to compute relative proper motions and select bona fide cluster members. Results from statistical isochrone fits in the NIR and optical-NIR colour-magnitude diagrams indicate an age of $12.8^{+0.9}_{-0.8}$ Gyr, confirming that HP$\,$1 is one of the oldest clusters in the Milky Way. The same fits also provide apparent distance moduli in the $K_{\rm S}$ and $V$ filters in very good agreement with the ones from 11 RR Lyrae stars. By subtracting the extinction in each filter, we recover a heliocentric distance of $6.59^{+0.17}_{-0.15}$ kpc. Furthermore, we refine the orbit of HP$\,$1 using this accurate distance and update and accurate radial velocities (from high resolution spectroscopy) and absolute proper motions (from Gaia DR2), reaching mean perigalactic and apogalactic distances of $\sim$0.12 and $\sim$3 kpc respectively.

In this paper we present the astro-photometric catalogues of 56 globular clusters and one open cluster. Astrometry and photometry are mainly based on images collected within the "HST Legacy Survey of Galactic Globular Clusters: Shedding UV Light on Their Populations and Formation" (GO-13297, PI:~Piotto), and the "ACS Survey of Galactic Globular Clusters" (GO-10775, PI:~Sarajedini). For each source in the catalogues for which we have reliable proper motion we also publish a membership probability for separation of field and cluster stars. These new catalogues, which we make public in Mikulski Archive for Space Telescopes, replace previous catalogues by Paper VIII of this series.

Line-of-sight kinematic studies indicate that many Galactic globular clusters have a significant degree of internal rotation. However, three-dimensional kinematics from a combination of proper motions and line-of-sight velocities are needed to unveil the role of angular momentum in the formation and evolution of these old stellar systems. Here we present the first quantitative study of internal rotation on the plane-of-the-sky for a large sample of globular clusters using proper motions from Gaia DR2. We detect signatures of rotation in the tangential component of proper motions for 11 out of 51 clusters at a $>$3-sigma confidence level, confirming the detection reported in Gaia collaboration et al. (2018) for 8 clusters, and additionally identify 11 GCs with a 2-sigma rotation detection. For the clusters with a detected global rotation, we construct the two-dimensional rotation maps and proper motion rotation curves, and we assess the relevance of rotation with respect to random motions ($V/\sigma\sim0.08-0.51$). We find evidence of a correlation between the degree of internal rotation and relaxation time, highlighting the importance of long-term dynamical evolution in shaping the clusters current properties. This is a strong indication that angular momentum must have played a fundamental role in the earliest phases of cluster formation. Finally, exploiting the spatial information of the rotation maps and a comparison with line-of-sight data, we provide an estimate of the inclination of the rotation axis for a subset of 8 clusters. Our work demonstrates the potential of Gaia data for internal kinematic studies of globular clusters and provides the first step to reconstruct their intrinsic three-dimensional structure.

We present an improved data-reduction technique to obtain high-precision proper motions (PMs) of globular clusters using Hubble Space Telescope data. The new reduction is superior to the one presented in the first paper of this series for the faintest sources in very crowded fields. We choose the globular cluster NGC 362 as a benchmark to test our new procedures. We measure PMs of 117 450 sources in the field, showing that we are able to obtain a PM precision better than 10 $\mu$as yr$^{-1}$ for bright stars. We make use of this new PM catalog of NGC 362 to study the cluster's internal kinematics. We investigate the velocity-dispersion profiles of the multiple stellar populations hosted by NGC 362 and find new pieces of information on the kinematics of first- and second-generation stars. We analyze the level of energy equipartition of the cluster and find direct evidence for its post-core-collapsed state from kinematic arguments alone. We refine the dynamical mass of the blue stragglers and study possible kinematic differences between blue stragglers formed by collisions and mass transfer. We also measure no significant cluster rotation in the plane of the sky. Finally, we measure the absolute PM of NGC 362 and of the background stars belonging to the Small Magellanic Cloud, finding a good agreement with previous estimates in the literature. We make the PM catalog publicly available.

The parameter A+, defined as the area enclosed between the cumulative radial distribution of blue straggler stars (BSSs) and that of a reference population, is a powerful indicator of the level of BSS central segregation. As part of the Hubble Space Telescope UV Legacy Survey of Galactic globular clusters (GCs), here we present the BSS population and the determination of A+ in 27 GCs observed out to about one half-mass radius. In combination with 21 additional clusters discussed in a previous paper this provides us with a global sample of 48 systems (corresponding to \sim 32\% of the Milky Way GC population), for which we find a strong correlation between A+ and the ratio of cluster age to the current central relaxation time. Tight relations have been found also with the core radius and the central luminosity density, which are expected to change with the long-term cluster dynamical evolution. An interesting relation is emerging between A+ and the ratio of the BSS velocity dispersion relative to that of main sequence turn-off stars, which measures the degree of energy equipartition experienced by BSSs in the cluster. These results provide further confirmation that BSSs are invaluable probes of GC internal dynamics and A+ is a powerful dynamical clock.

In the context of the Hubble Space Telescope UV Survey of Galactic Globular Clusters (GCs), we derived high-precision, multi-band photometry to investigate the multiple stellar populations in the massive and metal-poor GC M 15. By creating for red-giant branch (RGB) stars of the cluster a 'chromosome map', which is a pseudo two-colour diagram made with appropriate combination of F275W, F336W, F438W, and F814W magnitudes, we revealed colour spreads around two of the three already known stellar populations. These spreads cannot be produced by photometric errors alone and could hide the existence of (two) additional populations. This discovery increases the complexity of the multiple-population phenomenon in M 15. Our analysis shows that M 15 exhibits a faint sub-giant branch (SGB), which is also detected in colour-magnitude diagrams (CMDs) made with optical magnitudes only. This poorly-populated SGB includes about 5% of the total number of SGB stars and evolves into a red RGB in the m_F336W vs. m_F336W-m_F814W CMD, suggesting that M 15 belongs to the class of Type II GCs. We measured the relative number of stars in each population at various radial distances from the cluster centre, showing that all of these populations share the same radial distribution within statistic uncertainties. These new findings are discussed in the context of the formation and evolution scenarios of the multiple populations.

We present an analysis of a slightly eccentric ($e=0.05$), partially eclipsing long-period ($P = 69.73$ d) main sequence binary system (WOCS 12009, Sanders 1247) in the benchmark old open cluster M67. Using Kepler K2 and ground-based photometry along with a large set of new and reanalyzed spectra, we derived highly precise masses ($1.111\pm0.015$ and $0.748\pm0.005 M_\odot$) and radii ($1.071\pm0.008\pm0.003$ and $0.713\pm0.019\pm0.026 R_\odot$, with statistical and systematic error estimates) for the stars. The radius of the secondary star is in agreement with theory. The primary, however, is approximately $15\%$ smaller than reasonable isochrones for the cluster predict. Our best explanation is that the primary star was produced from the merger of two stars, as this can also account for the non-detection of photospheric lithium and its higher temperature relative to other cluster main sequence stars at the same $V$ magnitude. To understand the dynamical characteristics (low measured rotational line broadening of the primary star and the low eccentricity of the current binary orbit), we believe that the most probable (but not the only) explanation is the tidal evolution of a close binary within a primordial triple system (possibly after a period of Kozai-Lidov oscillations), leading to merger approximately 1Gyr ago. This star appears to be a future blue straggler that is being revealed as the cluster ages and the most massive main sequence stars die out.

Bulge globular clusters (GCs) with metallicities [Fe/H] <~ -1.0 and blue horizontal branches are candidates to harbor the oldest populations in the Galaxy. Based on the analysis of HST proper-motion-cleaned color-magnitude diagrams in filters F435W and F625W, we determine physical parameters for the old bulge globular clusters NGC 6522 and NGC 6626 (M28), both with well-defined blue horizontal branches. We compare these results with similar data for the inner halo cluster NGC 6362. These clusters have similar metallicities (-1.3 <= [Fe/H] <= -1.0) obtained from high resolution spectroscopy. We derive ages, distance moduli, and reddening values by means of statistical comparisons between observed and synthetic fiducial lines employing likelihood statistics and the Markov chain Monte Carlo method. The synthetic fiducial lines were generated using alpha-enhanced BaSTI and Dartmouth stellar evolutionary models, adopting both canonical (Y ~ 0.25) and enhanced (Y~ 0.30-0.33) helium abundances. RR Lyrae stars were employed to determine the HB magnitude level, providing an independent indicator to constrain the apparent distance modulus and the helium enhancement. The shape of the observed fiducial line could be compatible with some helium enhancement for NGC 6522 and NGC 6626, but the average magnitudes of RR Lyrae stars tend to rule out this hypothesis. Assuming canonical helium abundances, BaSTI and Dartmouth models indicate that all three clusters are coeval, with ages between ~ 12.5 and 13.0 Gyr. The present study also reveals that NGC 6522 has at least two stellar populations, since its CMD shows a significantly wide subgiant branch compatible with 14% +/- 2% and 86% +/- 5% for first and second generations, respectively.

In this second installment of the series, we look at the internal kinematics of the multiple stellar populations of the globular cluster $\omega$ Centauri in one of the parallel Hubble Space Telescope (HST) fields, located at about 3.5 half-light radii from the center of the cluster. Thanks to the over 15-year-long baseline and the exquisite astrometric precision of the HST cameras, well-measured stars in our proper-motion catalog have errors as low as $\sim 10\ \mu$as yr$^{-1}$, and the catalog itself extends to near the hydrogen-burning limit of the cluster. We show that second-generation (2G) stars are significantly more radially anisotropic than first-generation (1G) stars. The latter are instead consistent with an isotropic velocity distribution. In addition, 1G have excess systemic rotation in the plane of the sky with respect to 2G stars. We show that the six populations below the main-sequence (MS) knee identified in our first paper are associated to the five main population groups recently isolated on the upper MS in the core of cluster. Furthermore, we find both 1G and 2G stars in the field to be far from being in energy equipartition, with $\eta_{\rm 1G}=-0.007\pm0.026$ for the former, and $\eta_{\rm 2G}=0.074\pm0.029$ for the latter, where $\eta$ is defined so that the velocity dispersion $\sigma_\mu$ scales with stellar mass as $\sigma_\mu \propto m^{-\eta}$. The kinematical differences reported here can help constrain the formation mechanisms for the multiple stellar populations in $\omega$ Centauri and other globular clusters. We make our astro-photometric catalog publicly available.

In this paper we report a new estimate of the absolute proper motion (PM) of the globular cluster NGC 5139 ($\omega$ Cen) as part of the HST large program GO-14118+14662. We analyzed a field 17 arcmin South-West of the center of $\omega$ Cen and computed PMs with an epoch span of $\sim$15.1 years. We employed 45 background galaxies to link our relative PMs to an absolute reference-frame system. The absolute PM of the cluster in our field is: $(\mu_\alpha \cos\delta , \mu_\delta) = (-3.341 \pm 0.028 , -6.557 \pm 0.043)$ mas yr$^{-1}$. Upon correction for the effects of viewing perspective and the known cluster rotation, this implies that for the cluster center of mass $(\mu_\alpha \cos\delta , \mu_\delta) = (-3.238 \pm 0.028, -6.716 \pm 0.043)$ mas yr$^{-1}$. This measurement is direct and independent, has the highest random and systematic accuracy to date, and will provide an external verification for the upcoming Gaia Data Release 2. It also differs from most reported PMs for $\omega$ Cen in the literature by more than 5$\sigma$, but consistency checks compared to other recent catalogs yield excellent agreement. We computed the corresponding Galactocentric velocity, calculated the implied orbit of $\omega$ Cen in two different Galactic potentials, and compared these orbits to the orbits implied by one of the PM measurements available in the literature. We find a larger (by about 500 pc) perigalactic distance for $\omega$ Cen with our new PM measurement, suggesting a larger survival expectancy for the cluster in the Galaxy.

Located at ~2pc, the L7.5+T0.5 dwarfs system WISE J104915.57-531906.1 (Luhman16AB) is the third closest system known to Earth, making it a key benchmark for detailed investigation of brown dwarf atmospheric properties, thermal evolution, multiplicity, and planet-hosting frequency. In the first study of this series -- based on a multi-cycle Hubble Space Telescope (HST) program -- we provide an overview of the project and present improved estimates of positions, proper motions, annual parallax, mass ratio, and the current best assessment of the orbital parameters of the A-B pair. Our HST observations encompass the apparent periastron of the binary at 220.5+/-0.2 mas at epoch 2016.402. Although our data seem to be inconsistent with recent ground-based astrometric measurements, we also exclude the presence of third bodies down to Neptune masses and periods longer than a year.

We present the results of a seven-year-long radial velocity survey of a sample of 88 main-sequence and evolved stars to reveal signatures of Jupiter-mass planets in the solar-age and solar-metallicity open cluster M67. We aim at studying the frequency of giant planets in this cluster with respect to the field stars. In addition, our sample is also ideal to perform a long-term study to compare the chemical composition of stars with and without giant planets in detail. We analyzed precise radial velocity (RV) measurements obtained with five different instruments. We conducted Monte Carlo simulations to estimate the occurrence rate of giant planets in our radial velocity survey. All the planets previously announced in this RV campaign with their properties are summarized here: 3 hot Jupiters around the main-sequence stars YBP1194, YBP1514, and YBP401, and 1 giant planet around the evolved star S364. Two additional planet candidates around the stars YBP778 and S978 are also analyzed in the present work. We discuss stars that exhibit large RV variability or trends individually. For 2 additional stars, long-term trends are compatible with new binary candidates or substellar objects, which increases the total number of binary candidates detected in our campaign to 14. Based on the Doppler-detected planets discovered in this survey, we find an occurrence of giant planets of ~18.0%(+12.0/-8.0%) in the selected period-mass range. This frequency is slightly higher but consistent within the errors with the estimate for the field stars, which leads to the general conclusion that open cluster and field statistics agree. However, we find that the rate of hot Jupiters in the cluster (~5.7%(+5.5/-3.0%)) is substantially higher than in the field.

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 present a simultaneous, multi-wavelength campaign targeting the nearby (7.2 pc) L8/L9 (optical/near-infrared) dwarf WISEP J060738.65+242953.4 in the mid-infrared, radio, and optical. Spitzer Space Telescope observations show no variability at the 0.2% level over 10 hours each in the 3.6 and 4.5 micron bands. Kepler K2 monitoring over 36 days in Campaign 0 rules out stable periodic signals in the optical with amplitudes great than 1.5% and periods between 1.5 hours and 2 days. Non-simultaneous Gemini optical spectroscopy detects lithium, constraining this L dwarf to be less than ~2 Gyr old, but no Balmer emission is observed. The low measured projected rotation velocity (v sin i < 6 km/s) and lack of variability are very unusual compared to other brown dwarfs, and we argue that this substellar object is likely viewed pole-on. We detect quiescent (non-bursting) radio emission with the VLA. Amongst radio detected L and T dwarfs, it has the lowest observed L_nu and the lowest v sin i. We discuss the implications of a pole-on detection for various proposed radio emission scenarios.

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)

Kepler and K2 data analysis reported in the literature is mostly based on aperture photometry. Because of Kepler's large, undersampled pixels and the presence of nearby sources, aperture photometry is not always the ideal way to obtain high-precision photometry and, because of this, the data set has not been fully exploited so far. We present a new method that builds on our experience with undersampled HST images. The method involves a point-spread function (PSF) neighbour-subtraction and was specifically developed to exploit the huge potential offered by the K2 "super-stamps" covering the core of dense star clusters. Our test-bed targets were the NGC 2158 and M 35 regions observed during the K2 Campaign 0. We present our PSF modeling and demonstrate that, by using a high-angular-resolution input star list from the Asiago Schmidt telescope as the basis for PSF neighbour subtraction, we are able to reach magnitudes as faint as Kp~24 with a photometric precision of 10% over 6.5 hours, even in the densest regions. At the bright end, our photometric precision reaches ~30 parts-per-million. Our method leads to a considerable level of improvement at the faint magnitudes (Kp>15.5) with respect to the classical aperture photometry. This improvement is more significant in crowded regions. We also extracted raw light curves of ~60,000 stars and detrended them for systematic effects induced by spacecraft motion and other artifacts that harms K2 photometric precision. We present a list of 2133 variables.

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.

We present a new reduction pipeline for the VIRCAM@VISTA detector and describe the method developed to obtain high-precision astrometry with the VISTA Variables in the V\'ia L\'actea (VVV) data set. We derive an accurate geometric-distortion correction using as calibration field the globular cluster NGC 5139, and showed that we are able to reach a relative astrometric precision of about 8 mas per coordinate per exposure for well-measured stars over a field of view of more than 1 square degree. This geometric-distortion correction is made available to the community. As a test bed, we chose a field centered around the globular cluster NGC 6656 from the VVV archive and computed proper motions for the stars within. With 45 epochs spread over four years, we show that we are able to achieve a precision of 1.4 mas/yr and to isolate each population observed in the field (cluster, Bulge and Disk) using proper motions. We used proper-motion-selected field stars to measure the motion difference between Galactic disk and bulge stars. Our proper-motion measurements are consistent with UCAC4 and PPMXL, though our errors are much smaller. Models have still difficulties in reproducing the observations in this highly-reddened Galactic regions.

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.

Hercules is the prototype of the ultra-faint dwarf (UFD) galaxies. To date, there are still no firm constraints on its total luminosity, due to the difficulty of disentangling Hercules bona-fide stars from the severe Galactic field contamination. In order to better constrain Hercules properties we aim at removing foreground and background contaminants in the galaxy field using the proper motions of the Milky Way stars and the colour-colour diagram. We have obtained images of Hercules in the rSloan, BBessel and Uspec bands with the Large Binocular Telescope (LBT) and LBC-BIN mode capabilities. The rSloan new data-set combined with data from the LBT archive span a time baseline of about 5 yr, allowing us to measure for the first time proper motions of stars in the Hercules direction. The Uspec data along with existing LBT photometry allowed us to use colour-colour diagram to further remove the field contamination. Thanks to a highly-accurate procedure to derive the rSloan-filter geometric distortion solution for the LBC-red, we were able to measure stellar relative proper motions to a precision of better than 5 mas yr^-1 down to rSloan=22 mag and disentangle a significant fraction (\>90\%) of Milky Way contaminants. We ended up with a sample of 528 sources distributed over a large portion of the galaxy body (0.12 deg^2). Of these sources, 171 turned out to be background galaxies and additional foreground stars, from the analysis of the Uspec - BBessel vs. BBessel - rSloan colour-colour diagram. This leaves us with a sample of 357 likely members of the Hercules UFD. We compared the cleaned colour-magnitude diagram (CMD) with evolutionary models and synthetic CMDs, confirming the presence in Hercules of an old population (t=12\pm 2 Gyr), with a wide spread in metallicity (-3.3\<[Fe/H]\<-1.8).

High-precision astrometry requires accurate point-spread function modeling and accurate geometric-distortion corrections. This paper demonstrates that it is possible to achieve both requirements with data collected at the high acuity wide-field K-band imager (HAWK-I), a wide-field imager installed at the Nasmyth focus of UT4/VLT ESO 8m telescope. Our final astrometric precision reaches ~3 mas per coordinate for a well-exposed star in a single image with a systematic error less than 0.1 mas. We constructed calibrated astro-photometric catalogs and atlases of seven fields: the Baade's Window, NGC 6656, NGC 6121, NGC 6822, NGC 6388, NGC 104, and the James Webb Space Telescope calibration field in the Large Magellanic Cloud. We make these catalogs and images electronically available to the community. Furthermore, as a demonstration of the efficacy of our approach, we combined archival material taken with the optical wide-field imager at the MPI/ESO 2.2m with HAWK-I observations. We showed that we are able to achieve an excellent separation between cluster members and field objects for NGC 6656 and NGC 6121 with a time base-line of about 8 years. Using both HST and HAWK-I data, we also study the radial distribution of the SGB populations in NGC 6656 and conclude that the radial trend is flat within our uncertainty. We also provide membership probabilities for most of the stars in NGC 6656 and NGC 6121 catalogs and estimate membership for the published variable stars in these two fields.