We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Aims. We develop, test and characterise of a new statistical tool (intelligent system) for the sifting and analysis of nearby young open cluster (NYOC) populations. Methods. Using a Bayesian formalism, this statistical tool is able to obtain the posterior distributions of parameters governing the cluster model. It also uses hierarchical bayesian models to establish weakly informative priors, and incorporates the treatment of missing values and non-homogeneous (heteroscedastic) observational uncertainties. Results. From simulations, we estimate that this statistical tool renders kinematic (proper motion) and photometric (luminosity) distributions of the cluster population with a contamination rate of $5.8 \pm 0.2$ %. The luminosity distributions and present day mass function agree with the ones found by Bouy et al. (2015b) on the completeness interval of the survey. At the probability threshold of maximum accuracy, the classifier recovers $\sim$ 90% of Bouy et al. (2015b) candidate members and finds 10% of new ones. Conclusions. A new statistical tool for the analysis of NYOC is introduced, tested and characterised. Its comprehensive modelling of the data properties allows it to get rid of the biases present in previous works. In particular, those resulting from the use of only completely observed (non-missing) data and the assumption of homoskedastic uncertainties. Also, its Bayesian framework allows it to properly propagate observational uncertainties into membership probabilities and cluster velocity and luminosity distributions. Our results are in a general agreement with those from the literature, although we provide the most up-to-date and extended list of candidate members of the Pleiades cluster.

Methods. We compute Bayesian evidences and Bayes Factors for a set of variations of the classical radial models by King (1962), Elson et al. (1987) and Lauer et al. (1995). The variations incorporate different degrees of model freedom and complexity, amongst which we include biaxial (elliptical) symmetry, and luminosity segregation. As a by-product of the model comparison, we obtain posterior distributions and maximum a posteriori estimates for each set of model parameters. Results. We find that the model comparison results depend on the spatial extent of the region used for the analysis. For a circle of 11.5 parsecs around the cluster centre (the most homogeneous and complete region), we find no compelling reason to abandon Kings model, although the Generalised King model, introduced in this work, has slightly better fitting properties. Furthermore, we find strong evidence against radially symmetric models when compared to the elliptic extensions. Finally, we find that including mass segregation in the form of luminosity segregation in the J band, is strongly supported in all our models. Conclusions. We have put the question of the projected spatial distribution of the Pleiades cluster on a solid probabilistic framework, and inferred its properties using the most exhaustive and least contaminated list of Pleiades candidate members available to date. Our results suggest however that this sample may still lack about 20% of the expected number of cluster members. Therefore, this study should be revised when the completeness and homogeneity of the data can be extended beyond the 11.5 parsecs limit. Such study will allow a more precise determination of the Pleiades spatial distribution, its tidal radius, ellipticity, number of objects and total mass.

MWC758 is a young star surrounded by a transitional disk. Recently, a protoplanet candidate has been detected around MWC758 through high-resolution $L'$-band observations. The candidate is located inside the disk cavity at a separation of $\sim$111 mas from the central star, and at an average position angle of $\sim$165.5 degrees. We have performed simultaneous adaptive optics observations of MWC758 in the H$_{\alpha}$ line and the adjacent continuum using SPHERE/ZIMPOL at the Very Large Telescope (VLT). We aim at detecting accreting protoplanet candidates through spectral angular differential imaging observations. The data analysis does not reveal any H$_{\alpha}$ signal around the target. The derived contrast curve in the B_Ha filter allows us to derive a 5$\sigma$ upper limit of $\sim$7.6 mag at 111 mas, the separation of the previously detected planet candidate. This contrast translates into a H$_{\alpha}$ line luminosity of $L_{\rm H_{\alpha}}\lesssim$ 5$\times$10$^{-5}$ $L_{\odot}$ at 111 mas, and an accretion luminosity of $L_{acc} <$3.7$\times$10$^{-4}\,L_{\odot}$. For the predicted mass range of MWC758b, 0.5-5 $M_{\rm Jup}$, this implies accretion rates of $\dot M \lesssim$ 3.4$\times$(10$^{-8}$-10$^{-9})\,M_{\odot}/yr$, for an average planet radius of 1.1 $R_{\rm Jup}$. Therefore, our estimates are consistent with the predictions of accreting circumplanetary accretion models for $R_{\rm in} = 1 R_{\rm Jup}$. In any case, the non-detection of any H$_{\alpha}$ emitting source in the ZIMPOL images does not allow us to unveil the true nature of the $L'$ detected source.

The CARMENES radial velocity (RV) survey is observing 324 M dwarfs to search for any orbiting planets. In this paper, we present the survey sample by publishing one CARMENES spectrum for each M dwarf. These spectra cover the wavelength range 520--1710nm at a resolution of at least $R > 80,000$, and we measure its RV, H$\alpha$ emission, and projected rotation velocity. We present an atlas of high-resolution M-dwarf spectra and compare the spectra to atmospheric models. To quantify the RV precision that can be achieved in low-mass stars over the CARMENES wavelength range, we analyze our empirical information on the RV precision from more than 6500 observations. We compare our high-resolution M-dwarf spectra to atmospheric models where we determine the spectroscopic RV information content, $Q$, and signal-to-noise ratio. We find that for all M-type dwarfs, the highest RV precision can be reached in the wavelength range 700--900nm. Observations at longer wavelengths are equally precise only at the very latest spectral types (M8 and M9). We demonstrate that in this spectroscopic range, the large amount of absorption features compensates for the intrinsic faintness of an M7 star. To reach an RV precision of 1ms$^{-1}$ in very low mass M dwarfs at longer wavelengths likely requires the use of a 10m class telescope. For spectral types M6 and earlier, the combination of a red visual and a near-infrared spectrograph is ideal to search for low-mass planets and to distinguish between planets and stellar variability. At a 4m class telescope, an instrument like CARMENES has the potential to push the RV precision well below the typical jitter level of 3-4ms$^{-1}$.

Context: The main goal of the CARMENES survey is to find Earth-mass planets around nearby M-dwarf stars. Seven M-dwarfs included in the CARMENES sample had been observed before with HIRES and HARPS and either were reported to have one short period planetary companion (GJ15A, GJ176, GJ436, GJ536 and GJ1148) or are multiple planetary systems (GJ581 and GJ876). Aims: We aim to report new precise optical radial velocity measurements for these planet hosts and test the overall capabilities of CARMENES. Methods: We combined our CARMENES precise Doppler measurements with those available from HIRES and HARPS and derived new orbital parameters for the systems. Bona-fide single planet systems are fitted with a Keplerian model. The multiple planet systems were analyzed using a self-consistent dynamical model and their best fit orbits were tested for long-term stability. Results: We confirm or provide supportive arguments for planets around all the investigated stars except for GJ15A, for which we find that the post-discovery HIRES data and our CARMENES data do not show a signal at 11.4 days. Although we cannot confirm the super-Earth planet GJ15Ab, we show evidence for a possible long-period ($P_{\rm c}$ = 7025$_{-629}^{+972}$ d) Saturn-mass ($m_{\rm c} \sin i$ = 51.8$_{-5.8}^{+5.5}M_\oplus$) planet around GJ15A. In addition, based on our CARMENES and HIRES data we discover a second planet around GJ1148, for which we estimate a period $P_{\rm c}$ = 532.6$_{-2.5}^{+4.1}$ d, eccentricity $e_{\rm c}$ = 0.34$_{-0.06}^{+0.05}$ and minimum mass $m_{\rm c} \sin i$ = 68.1$_{-2.2}^{+4.9}M_\oplus$. Conclusions: The CARMENES optical radial velocities have similar precision and overall scatter when compared to the Doppler measurements conducted with HARPS and HIRES. We conclude that CARMENES is an instrument that is up to the challenge of discovering rocky planets around low-mass stars.

In this work we present ALMA continuum observations at 880 $\mu$m of 30 sub-mm cores previously identified with APEX/LABOCA at 870$\mu$m in the Barnard 30 cloud. The main goal is to characterize the youngest and lowest mass population in the cloud. As a result, we report the detection of five (out of 30) spatially unresolved sources with ALMA, with estimated masses between 0.9 and 67 M$_{\rm Jup}$. From these five sources, only two show gas emission. The analysis of multi-wavelength photometry from these two objects, namely B30-LB14 and B30-LB19, is consistent with one Class II- and one Class I low-mass stellar object, respectively. The gas emission is consistent with a rotating disk in the case of B30-LB14, and with an oblate rotating envelope with infall signatures in the case of LB19. The remaining three ALMA detections do not have infrared counterparts and can be classified as either deeply embedded objects or as starless cores if B30 members. In the former case, two of them (LB08 and LB31) show internal luminosity upper limits consistent with Very Low Luminosity objects, while we do not have enough information for LB10. In the starless core scenario, and taking into account the estimated masses from ALMA and the APEX/LABOCA cores, we estimate final masses for the central objects in the substellar domain, so they could be classified as pre-BD core candidates.

The early evolutionary stage of brown dwarfs are not very well characterized, specially during the embedded phase. To gain insight into the dominant formation mechanism of very low-mass objects and brown dwarfs, we conducted deep observations at 870$\mu$m with the LABOCA bolometer at the APEX telescope. Our goal was to identify young sub-mm sources in the Barnard 30 dark cloud. We complemented these data with multi-wavelength observations from the optical to the far-IR and. As a result, we have identified 34 submm sources and a substantial number of possible and probable Barnard 30 members within each individual APEX/LABOCA beam. They can be classified in three distinct groups. First, 15 out of these 34 have a clear optical or IR counterpart to the submm peak and nine of them are potential proto-BDs candidates. Moreover, a substantial number of them could be multiple systems. A second group of 13 sources comprises candidate members with significant infrared excesses located away from the central submm emission. All of them include brown dwarf candidates, some displaying IR excess, but their association with submm emission is unclear. In addition, we have found six starless cores and, based on the total dust mass estimate, three might be pre-substellar (or pre-BDs) cores. Finally, the complete characterization of our APEX/LABOCA sources, focusing on those detected at 24 and/or 70 $\mu$m, indicates that in our sample of 34 submm sources there are, at least: two WTTs, four CTTs, five young stellar objects (YSOs), eight proto-BD candidates (with another three dubious cases), and one Very Low Luminosity object (VeLLO).

The detection of Earth-like planets, exocomets or Kuiper belts show that the different components found in the solar system should also be present in other planetary systems. Trojans are one of these components and can be considered fossils of the first stages in the life of planetary systems. Their detection in extrasolar systems would open a new scientific window to investigate formation and migration processes. In this context, the main goal of the TROY project is to detect exotrojans for the first time and to measure their occurrence rate (eta-Trojan). In this first paper, we describe the goals and methodology of the project. Additionally, we used archival radial velocity data of 46 planetary systems to place upper limits on the mass of possible trojans and investigate the presence of co-orbital planets down to several tens of Earth masses. We used archival radial velocity data of 46 close-in (P<5 days) transiting planets (without detected companions) with information from high-precision radial velocity instruments. We took advantage of the time of mid-transit and secondary eclipses (when available) to constrain the possible presence of additional objects co-orbiting the star along with the planet. This, together with a good phase coverage, breaks the degeneracy between a trojan planet signature and signals coming from additional planets or underestimated eccentricity. We identify nine systems for which the archival data provide 1-sigma evidence for a mass imbalance between L4 and L5. Two of these systems provide 2-sigma detection, but no significant detection is found among our sample. We also report upper limits to the masses at L4/L5 in all studied systems and discuss the results in the context of previous findings.

The multi-planetary system HD 106315 was recently found in K2 data . The planets have periods of $P_b \sim9.55$ and $P_c \sim 21.06\,$days, and radii of $ r_b = 2.44 \pm 0.17\, $ and $r_c = 4.35 \pm 0.23\, $ $R_{\oplus}$. The brightness of the host star (V=9.0 mag) makes it an excellent target for transmission spectroscopy. However, to interpret transmission spectra it is crucial to measure the planetary masses. We obtained high precision radial velocities for HD~106315 to determine the mass of the two transiting planets discovered with Kepler K2. Our successful observation strategy was carefully tailored to mitigate the effect of stellar variability. We modelled the new radial velocity data together with the K2 transit photometry and a new ground-based partial transit of HD 106315c to derive system parameters. We estimate the mass of HD 106315b to be 12.6 $\pm$ 3.2 $M_{\oplus}$ and the density to be $4.7 \pm 1.7\, g\,cm^{-3}$, while for HD 106315c we estimate a mass of 15.2 $\pm$ 3.7 $M_{\oplus}$ and a density of $1.01 \pm 0.29\, $g\,cm$^{-3}$. Hence, despite planet c having a radius almost twice as large as planet b, their masses are consistent with one another. We conclude that HD 106315c has a thick hydrogen-helium gaseous envelope. A detailed investigation of HD 106315b using a planetary interior model constrains the core mass fraction to be 5-29\%, and the water mass fraction to be 10-50\%. An alternative, not considered by our model, is that HD 106315b is composed of a large rocky core with a thick H-He envelope. Transmission spectroscopy of these planets will give insight into their atmospheric compositions and also help constrain their core compositions.

We report the discovery of the exoplanet K2-110 b (previously EPIC212521166b) from K2 photometry orbiting in a 13.8637d period around an old, metal-poor K3 dwarf star. With a V-band magnitude of 11.9, K2-110 is particularly amenable to RV follow-up. A joint analysis of K2 photometry and high-precision RVs from 28 HARPS and HARPS-N spectra reveal it to have a radius of 2.6$\pm 0.1 R_{\oplus}$ and a mass of 16.7$\pm 3.2$~M$_{\oplus}$, hence a density of $5.2\pm1.2$ g.cm$^{-3}$, making it one of the most massive planets yet to be found with a sub-Neptune radius. When accounting for compression, the resulting Earth-like density is best fitted by a 0.2 M$_{\oplus}$ hydrogen atmosphere over an 16.5 M$_{\oplus}$ Earth-like interior, although the planet could also have significant water content. At 0.1~AU, even taking into account the old stellar age of $8 \pm 3$ Gyr, the planet is unlikely to have been significantly affected by EUV evaporation. However the planet likely disc-migrated to its current position making the lack of a thick H$_2$ atmosphere puzzling. This analysis has made K2-110 b one of the best-characterised mini-Neptunes with density constrained to less than 30%.

We identify and investigate known ultracool stars and brown dwarfs that are being observed or indirectly constrained by the Gaia mission. These objects will be the core of the Gaia ultracool dwarf sample composed of all dwarfs later than M7 that Gaia will provide direct or indirect information on. We match known L and T dwarfs to the Gaia first data release, the Two Micron All Sky Survey and the Wide-field Infrared Survey Explorer AllWISE survey and examine the Gaia and infrared colours, along with proper motions, to improve spectral typing, identify outliers and find mismatches. There are 321 L and T dwarfs observed directly in the Gaia first data release, of which 10 are later than L7. This represents 45 % of all the known LT dwarfs with estimated Gaia G magnitudes brighter than 20.3 mag. We determine proper motions for the 321 objects from Gaia and the Two Micron All Sky Survey positions. Combining the Gaia and infrared magnitudes provides useful diagnostic diagrams for the determination of L and T dwarf physical parameters. We then search the Tycho-Gaia astrometric solution Gaia first data release subset to find any objects with common proper motions to known L and T dwarfs and a high probability of being related. We find 15 new candidate common proper motion systems.

We present results from high-resolution, optical to near-IR imaging of host stars of Kepler Objects of Interest (KOIs), identified in the original Kepler field. Part of the data were obtained under the Kepler imaging follow-up observation program over seven years (2009 - 2015). Almost 90% of stars that are hosts to planet candidates or confirmed planets were observed. We combine measurements of companions to KOI host stars from different bands to create a comprehensive catalog of projected separations, position angles, and magnitude differences for all detected companion stars (some of which may not be bound). Our compilation includes 2297 companions around 1903 primary stars. From high-resolution imaging, we find that ~10% (~30%) of the observed stars have at least one companion detected within 1" (4"). The true fraction of systems with close (< ~4") companions is larger than the observed one due to the limited sensitivities of the imaging data. We derive correction factors for planet radii caused by the dilution of the transit depth: assuming that planets orbit the primary stars or the brightest companion stars, the average correction factors are 1.06 and 3.09, respectively. The true effect of transit dilution lies in between these two cases and varies with each system. Applying these factors to planet radii decreases the number of KOI planets with radii smaller than 2 R_Earth by ~2-23% and thus affects planet occurrence rates. This effect will also be important for the yield of small planets from future transit missions such as TESS.

The young (~2 Myr) cluster Chamaeleon I is one of the closest laboratories to study the early stages of star cluster dynamics in a low-density environment. We studied its structural and kinematical properties combining parameters from the high-resolution spectroscopic survey Gaia-ESO with data from the literature. Our main result is the evidence of a large discrepancy between the velocity dispersion (sigma = 1.14 \pm 0.35 km s^{-1}) of the stellar population and the dispersion of the pre-stellar cores (~0.3 km s^{-1}) derived from submillimeter observations. The origin of this discrepancy, which has been observed in other young star clusters is not clear. It may be due to either the effect of the magnetic field on the protostars and the filaments, or to the dynamical evolution of stars driven by two-body interactions. Furthermore, the analysis of the kinematic properties of the stellar population put in evidence a significant velocity shift (~1 km s^{-1}) between the two sub-clusters located around the North and South main clouds. This result further supports a scenario, where clusters form from the evolution of multiple substructures rather than from a monolithic collapse. Using three independent spectroscopic indicators (the gravity indicator $\gamma$, the equivalent width of the Li line, and the H_alpha 10\% width), we performed a new membership selection. We found six new cluster members located in the outer region of the cluster. Starting from the positions and masses of the cluster members, we derived the level of substructure Q, the surface density \Sigma and the level of mass segregation $\Lambda_{MSR}$ of the cluster. The comparison between these structural properties and the results of N-body simulations suggests that the cluster formed in a low density environment, in virial equilibrium or supervirial, and highly substructured.

Although mid-to-late type M dwarfs are the most common stars in our stellar neighborhood, our knowledge of these objects is still limited. Open questions include the evolution of their angular momentum, internal structures, dust settling in their atmospheres, age dispersion within populations. In addition, at young ages, late-type Ms have masses below the hydrogen burning limit and therefore are key objects in the debate on the brown dwarf mechanism of formation. In this work we determine and study in detail the physical parameters of two samples of young, late M-type sources belonging to either the Chamaeleon I Dark Cloud or the TW Hydrae Association and compare them with the results obtained in the literature for other young clusters and also for older, field, dwarfs. We used multi-wavelength photometry to construct and analyze SEDs to determine general properties of the photosphere and disk presence. We also used low resolution optical and near-infrared spectroscopy to study activity, accretion, gravity and effective temperature sensitive indicators. We propose a VO-based spectral index that is both temperature and age sensitive. We derived physical parameters using independent techniques confirming the already common feature/problem of the age/luminosity spread. In particular, we highlight two brown dwarfs showing very similar temperatures but clearly different surface gravity (explained invoking extreme early accretion). We also show how, despite large improvement in the dust treatment in theoretical models, there is still room for further progress in the simultaneous reproduction of the optical and near-infrared features of these cold young objects.

We show that non-magnetic models for the evolution of pre-main-sequence (PMS) stars *cannot* simultaneously describe the colour-magnitude diagram (CMD) and the pattern of lithium depletion seen in the cluster of young, low-mass stars surrounding $\gamma^2$ Velorum. The age of 7.5+/-1 Myr inferred from the CMD is much younger than that implied by the strong Li depletion seen in the cluster M-dwarfs and the Li depletion occurs at much redder colours than predicted. The epoch at which a star of a given mass depletes its Li and the surface temperature of that star are both dependent on its radius. We demonstrate that if the low-mass stars have radii ~10 per cent larger at a given mass and age, then both the CMD and Li depletion pattern of the Gamma Vel cluster are explained at a common age of 18-21 Myr. This radius inflation could be produced by some combination of magnetic suppression of convection and extensive cool starspots. Models that incorporate radius inflation suggest that PMS stars similar to those in the Gamma Vel cluster, in the range 0.2<M/Msun<0.7, are at least a factor of two older and ~7 per cent cooler than previously thought and that their masses are much larger (by >30 per cent) than inferred from conventional, non-magnetic models in the Hertzsprung-Russell diagram. Systematic changes of this size may be of great importance in understanding the evolution of young stars, disc lifetimes and the formation of planetary systems.

Stellar clusters are open windows to understand stellar evolution. Specifically, the change with time and the dependence on mass of different stellar properties. As such, they are our laboratories where different theories can be tested. We try to understand the origin of the connection between lithium depletion in F, G and K stars, rotation and activity, in particular in the Pleiades open cluster. We have collected all the relevant data in the literature, including information regarding rotation period, binarity and activity, and cross-matched with proper motions, multi-wavelength photometry and membership probability from the DANCe database. In order to avoid biases, only Pleiades single members with probabilities larger than p=0.75 have been included in the discussion. Results. The analysis confirms that there is a strong link between activity, rotation and the lithium equivalent width excess, specially for the range Lum(bol) = 0.5-0.2 Lsun (about K2-K7 spectral types or 0.75-0.95 Msun). It is not possible to disentangle these effects but we cannot exclude that the observed lithium overabundance is partially an observational effect due to enhanced activity, due to a large coverage by stellar spots induced by high rotation rates. Since a bona fide lithium enhancement is present in young, fast rotators, both activity and rotation should play a role in the lithium problem.

We report the discovery of EPIC201702477b, a transiting brown dwarf in a long period (40.73691 +/- 0.00037 day) and eccentric (e=0.2281 +/- 0.0026) orbit. This system was initially reported as a planetary candidate based on two transit events seen in K2 Campaign 1 photometry and later validated as an exoplanet. We confirm the transit and refine the ephemeris with two subsequent ground-based detections of the transit using the LCOGT 1m telescope network. We rule out any transit timing variations above the level of 30s. Using high precision radial velocity measurements from HARPS and SOPHIE we identify the transiting companion as a brown dwarf with a mass, radius, and bulk density of 66.9 +/- 1.7 M$_J$, 0.757 +/- 0.065 R$_J$, and 191+/-51 g.cm$^{-3}$ respectively. EPIC201702477b is the smallest radius brown dwarf yet discovered, with a mass just below the H-burning limit. It has the highest density of any planet, substellar mass object or main-sequence star discovered so far. We find evidence in the set of known transiting brown dwarfs for two populations of objects - high mass brown dwarfs and low mass brown dwarfs. The higher-mass population have radii in very close agreement to theoretical models, and show a lower-mass limit around 60 M$_J$. This may be the signature of mass-dependent ejection of systems during the formation process.

Light curve modulations in the sample of Kepler planet candidates allows the disentangling of the nature of the transiting object by photometrically measuring its mass. This is possible by detecting the effects of the gravitational pull of the companion (ellipsoidal modulations) and in some cases, the photometric imprints of the Doppler effect when observing in a broad band (Doppler beaming). We aim to photometrically unveil the nature of some transiting objects showing clear modulations in the phase-folded Kepler light curve. We selected a subsample among the large crop of Kepler objects of interest (KOIs) based on their chances to show detectable light curve modulations, i.e., close ($a<12~R_{\star}$) and large (in terms of radius) candidates. We modeled their phase-folded light curves with consistent equations for the three effects, namely, reflection, ellipsoidal and beaming (known as REB modulations). We provide detailed general equations for the fit of the REB modulations for the case of eccentric orbits. These equations are accurate to the photometric precisions achievable by current and forthcoming instruments and space missions. By using this mathematical apparatus, we find three close-in {very low-mass companions (two of them in the brown dwarf mass domain)} orbiting main-sequence stars (KOI-554, KOI-1074, and KOI-3728), and reject the planetary nature of the transiting objects (thus classifying them as false positives). In contrast, the detection of the REB modulations and transit/eclipse signal allows the measurement of their mass and radius that can provide important constraints for modeling their interiors since just a few cases of low-mass eclipsing binaries are known. Additionally, these new systems can help to constrain the similarities in the formation process of the more massive and close-in planets (hot Jupiters), brown dwarfs, and very low-mass companions.

We use high quality K2 light curves for hundreds of stars in the Pleiades to understand better the angular momentum evolution and magnetic dynamos of young, low mass stars. The K2 light curves provide not only rotational periods but also detailed information from the shape of the phased light curve not available in previous studies. A slowly rotating sequence begins at $(V-K_{\rm s})_0\sim$1.1 (spectral type F5) and ends at $(V-K_{\rm s})_0\sim$ 3.7 (spectral type K8), with periods rising from $\sim$2 to $\sim$11 days in that interval. Fifty-two percent of the Pleiades members in that color interval have periods within 30\% of a curve defining the slow sequence; the slowly rotating fraction decreases significantly redward of $(V-K_{\rm s})_0$=2.6. Nearly all of the slow-sequence stars show light curves that evolve significantly on timescales less than the K2 campaign duration. The majority of the FGK Pleiades members identified as photometric binaries are relatively rapidly rotating, perhaps because binarity inhibits star-disk angular momentum loss mechanisms during pre-main sequence evolution. The fully convective, late M dwarf Pleiades members (5.0 $<(V-K_{\rm s})_0<$ 6.0) nearly always show stable light curves, with little spot evolution or evidence of differential rotation. During PMS evolution from $\sim$3 Myr (NGC2264 age) to $\sim$125 Myr (Pleiades age), stars of 0.3 $M_{\odot}$ shed about half their angular momentum, with the fractional change in period between 3 and 125 Myr being nearly independent of mass for fully convective stars. Our data also suggest that very low mass binaries form with rotation periods more similar to each other and faster than would be true if drawn at random from the parent population of single stars.

Young (125 Myr), populous ($>$1000 members), and relatively nearby, the Pleiades has provided an anchor for stellar angular momentum models for both younger and older stars. We used K2 to explore the distribution of rotation periods in the Pleiades. With more than 500 new periods for Pleiades members, we are vastly expanding the number of Pleiads with periods, particularly at the low mass end. About 92\% of the members in our sample have at least one measured spot-modulated rotation period. For the $\sim$8\% of the members without periods, non-astrophysical effects often dominate (saturation, etc.), such that periodic signals might have been detectable, all other things being equal. We now have an unusually complete view of the rotation distribution in the Pleiades. The relationship between $P$ and $(V-K_{\rm s})_0$ follows the overall trends found in other Pleiades studies. There is a slowly rotating sequence for $1.1\lesssim(V-K_{\rm s})_0\lesssim 3.7$, and a primarily rapidly rotating population for $(V-K_{\rm s})_0\gtrsim 5.0$. There is a region in which there seems to be a disorganized relationship between $P$ and $(V-K_{\rm s})_0$ for $3.7 \lesssim(V-K_{\rm s})_0\lesssim 5.0$. Paper II continues the discussion, focusing on multi-period structures, and Paper III speculates about the origin and evolution of the period distribution in the Pleiades.

We use K2 to continue the exploration of the distribution of rotation periods in Pleiades that we began in Paper I. We have discovered complicated multi-period behavior in Pleiades stars using these K2 data, and we have grouped them into categories, which are the focal part of this paper. About 24% of the sample has multiple, real frequencies in the periodogram, sometimes manifesting as obvious beating in the light curves. Those having complex and/or structured periodogram peaks, unresolved multiple periods, and resolved close multiple periods are likely due to spot/spot group evolution and/or latitudinal differential rotation; these largely compose the slowly rotating sequence in $P$ vs.~$(V-K_{\rm s})_0$ identified in Paper I. The fast sequence in $P$ vs.~$(V-K_{\rm s})_0$ is dominated by single-period stars; these are likely to be rotating as solid bodies. Paper III continues the discussion, speculating about the origin and evolution of the period distribution in the Pleiades.

We report the discovery of the two hot-Jupiters K2-30b and K2-34b. The two planets were detected transiting their main-sequence star with periods $\sim$ 4.099 and $\sim$ 2.996 days, in campaigns 4 and 5 of the extension of the Kepler mission, K2. Subsequent ground-based radial velocity follow-up with SOPHIE, HARPS-N and CAFE, established the planetary nature of the transiting objects. We analyzed the transit signal, radial velocity and spectral energy distributions of the two systems to characterize their properties. Both planets (K2-30b and K2-34b) are bloated hot-Jupiters (1.25 $R_{\rm Jup}$ and 1.33 $R_{\rm Jup}$) around relatively bright (V =13.5 and V=11.5), slow rotating main-sequence (G8 and F9) stars. Thus, these systems are good candidates for detecting the Rossiter-MacLaughlin effect to measure their obliquity and for atmospheric studies.

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.

Extrasolar planets abound in almost any possible configuration. However, until five years ago, there was a lack of planets orbiting closer than 0.5 au to giant or subgiant stars. Since then, recent detections have started to populated this regime by confirming 13 planetary systems. We discuss the properties of these systems in terms of their formation and evolution off the main sequence. Interestingly, we find that $70.0\pm6.6$ % of the planets in this regime are inner components of multiplanetary systems. This value is 4.2$\sigma$ higher than for main-sequence hosts, which we find to be $42.4\pm0.1$ %. The properties of the known planets seem to indicate that the closest-in planets (a < 0.06 au) to main-sequence stars are massive (i.e., hot Jupiters) and isolated and that they are subsequently engulfed by their host as it evolves to the red giant branch, leaving only the predominant population of multiplanetary systems in orbits 0.06 < a < 0.5 au. We discuss the implications of this emerging observational trend in the context of formation and evolution of hot Jupiters.