The detailed chemical composition of stars is important in many astrophysical fields, among which the characterisation of exoplanetary systems. Previous studies seem to indicate an anomalous chemical pattern of the youngest stellar population in the solar vicinity with a sub-solar metal content. This can influence various observational relations linking the properties of exoplanets to the characteristics of the host stars, for example the giant planet-metallicity relation. In this framework, we aim to expand our knowledge of the chemical composition of intermediate-age stars and understand whether these peculiarities are real or related to spectroscopic analysis techniques. We analysed high-resolution optical and near-infrared GIARPS spectra of intermediate-age stars (< 700Myr). To overcome issues related to the young ages of the stars, we applied a new spectroscopic method that uses titanium lines to derive the atmospheric parameters, in particular surface gravities and microturbulence velocity parameter. We also derived abundances of 14 different atomic species. The lack of systematic trends between elemental abundances and effective temperatures validates our method. However, we observed that the coolest (<5400 K) stars in the sample, display higher abundances for the Cr II, and for high-excitation potential C I lines. We found a positive correlation between the higher abundances measured of C I and Cr II and the activity index logR$_{HK}$. Instead, we found no correlations between the C abundances obtained from CH molecular band at 4300\AA, and both effective temperatures and activity. Thus, we suggest that these are better estimates for C abundances in young and cool stars. Finally, we found an indication of an increasing abundance ratio [X/H] with the condensation temperature for HD167389, indicating possible episodes of planet engulfment.

Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate.

Young nearby stars are good candidates in the search for planets with both radial velocity (RV) and direct imaging techniques. This, in turn, allows for the computation of the giant planet occurrence rates at all separations. The RV search around young stars is a challenge as they are generally faster rotators than older stars of similar spectral types and they exhibit signatures of magnetic activity (spots) or pulsation in their RV time series. Specific analyses are necessary to characterize, and possibly correct for, this activity. Our aim is to search for planets around young nearby stars and to estimate the giant planet (GP) occurrence rates for periods up to 1000 days. We used the HARPS spectrograph on the 3.6m telescope at La Silla Observatory to observe 89 A-M young (< 600 Myr) stars. We used our SAFIR (Spectroscopic data via Analysis of the Fourier Interspectrum Radial velocities ) software to compute the RV and other spectroscopic observables. Then, we computed the companion occurrence rates on this sample. We confirm the binary nature of HD177171, HD181321 and HD186704. We report the detection of a close low mass stellar companion for HIP36985. No planetary companion was detected. We obtain upper limits on the GP (< 13 MJup) and BD (13-80 MJup) occurrence rates based on 83 young stars for periods less than 1000 days, which are set, 2_-2^+3 % and 1_-1^+3 %.

Context. HD72946 is a bright and nearby solar-type star hosting a low-mass companion at long period (P~16 yr) detected with the radial velocities (RV) method. The companion has a minimum mass of 60.4+/-2.2 MJ and might be a brown dwarf. Its expected semi-major axis of ~243 mas makes it a suitable target for further characterization with high-contrast imaging, in particular to measure its inclination, mass, and spectrum and thus definitely establish its substellar nature. Aims. We aim to further characterize the orbit, atmosphere, and physical nature of HD72946B. Methods. We present high-contrast imaging data in the near-infrared with the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. We also use proper motion measurements of the star from Hipparcos and Gaia. Results. The SPHERE data reveal a point source with a contrast of ~9 mag at a projected separation of ~235 mas. No other point sources are detected in the field of view. By jointly fitting the RV, imaging, and proper motion data, we constrain all the orbital parameters of HD72946B and assess a dynamical mass of 72.4+/-1.6 MJ and a semi-major axis of 6.45$^{+0.08}_{-0.07}$ au. Empirical comparison of its IFS spectrum to template dwarfs indicates a spectral type of L5.0+/-1.5. The J-H3 color is close to the expectations of the DUSTY models and suggests a cloudy atmosphere. Comparison with atmospheric models of the spectrophotometry suggests an effective temperature of ~1700 K. The bolometric luminosity (log(L/LS)=-4.11+/-0.10 dex) and dynamical mass of HD72946B are more compatible with evolutionary models for an age range of ~0.9-3 Gyr. The formation mechanism of the companion is currently unclear as it appears slightly away from the bulk of model predictions. HD72946B is currently the closest benchmark brown dwarf companion to a solar-type star with imaging, RV, and proper motion measurements.

R~CrA is the brightest member of the Coronet star forming region and it is the closest Herbig AeBe star with a spectrum dominated by emission lines. Its luminosity has been monitored since the end of the 19th century, but the origin of its variability, which shows a stable period of $65.767\pm 0.007$~days, is still unknown. We studied photometric and spectroscopic data for this star to investigate the nature of the variability of R~CrA. We exploited the fact that near infrared luminosity of the Herbig AeBe stars is roughly proportional to the total luminosity of the stars to derive the absorption, and then mass and age of R~CrA. In addition, we model the periodic modulation of the light curve as due to partial attenuation of a central binary by a circumbinary disk. This model reproduces very well the observations. We found that the central object in R~CrA is a very young ($1.5\pm 1.5$~Myr), highly absorbed ($A_V=5.47\pm 0.4$~mag) binary; we obtain masses of $M_A=3.02\pm 0.43$~M$_\odot$ and $M_B=2.32\pm 0.35$~M$_\odot$ for the two components. We propose that the secular decrease of the R~CrA apparent luminosity is due to a progressive increase of the disk absorption. This might be related to precession of a slightly inclined disk caused by the recently discovered M-dwarf companion. Thus, R~CrA may be a triple system hosting a disk.

Context. Direct imaging of debris discs gives important information about their nature, their global morphology, and allows us to identify specific structures possibly in connection with the presence of gravitational perturbers. It is the most straightforward technique to observe planetary systems as a whole. Aims. We present the first resolved images of the debris disc around the young F-type star HD 160305, detected in scattered light using the VLT/SPHERE instrument in the near infrared. Methods. We used a post-processing method based on angular differential imaging and synthetic images of debris discs produced with a disc modelling code (GRaTer) to constrain the main characteristics of the disc around HD 160305. All of the point sources in the field of the IRDIS camera were analysed with an astrometric tool to determine whether they are bound objects or background stars. Results. We detect a very inclined (~ 82{\deg}) ring-like debris disc located at a stellocentric distance of about 86au (deprojected width ~27 au). The disc displays a brightness asymmetry between the two sides of the major axis, as can be expected from scattering properties of dust grains. We derive an anisotropic scattering factor g>0.5. A second right-left asymmetry is also observed with respect to the minor axis. We measure a surface brightness ratio of 0.73 $\pm$ 0.18 between the bright and the faint sides. Because of the low signal-to-noise ratio (S/N) of the images we cannot easily discriminate between several possible explanations for this left-right asymmetry, such as perturbations by an unseen planet, the aftermath of the breakup of a massive planetesimal, or the pericenter glow effect due to an eccentric ring. Two epochs of observations allow us to reject the companionship hypothesis for the 15 point sources present in the field.

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.

Young low-mass stars of equal-mass exhibit a distribution of rotation periods. At the very early phases of stellar evolution, this distribution is set by the star-disc locking mechanism. The primordial disc lifetime and, consequently, the duration of the disc-locking mechanism, can be significantly shortened by the presence of a close companion, making the rotation period distribution of close binaries different from that of either single stars or wide binaries. We use new data to investigate and better constrain the range of ages, the components separation and the mass ratio dependence at which the rotation period distribution has been significantly affected by the disc dispersal that is enhanced by close companions. We select a sample of close binaries in the Upper Scorpius association (age $\sim$8 Myr) whose components have measured the separation and the rotation periods and compare their period distribution with that of coeval stars that are single stars. We find that components of close binaries have on average rotation periods shorter than single stars. More precisely, binaries with about equal-mass components (0.9 $\le$ M2/M1 $\le$ 1.0) have rotation periods on average by $\sim$0.4 d shorter than single stars; binaries with smaller mass ratios (0.8 $<$ M2/M1 $<$ 0.9) have rotation periods on average by $\sim$1.9 d the primary components, and by $\sim$1.0 d the secondary components shorter than single stars. A comparison with the older 25-Myr $\beta$ Pictoris association shows that, whereas in the latter all close binaries with projected separation $\rho$ $\le$ 80 AU all rotate faster than single stars, in the Upper Scorpius that has happened for about 70% stars, yet. We interpret the enhanced rotation in close binaries with respect to single stars as the consequence of an early disc dispersal induced by the presence of close companions.

Context. The 51 Eridani system harbors a complex architecture with its primary star forming a hierarchical system with the binary GJ 3305AB at a projected separation of 2000 au, a giant planet orbiting the primary star at 13 au, and a low-mass debris disk around the primary star with possibly a cold component and a warm component inferred from the spectral energy distribution. Aims. We aim to better constrain the orbital parameters of the known giant planet. Methods. We monitored the system over three years from 2015 to 2018 with the VLT/SPHERE exoplanet imaging instrument. Results. We measure an orbital motion for the planet of ~130 mas with a slightly decreasing separation (~10 mas) and find a hint of curvature. This potential curvature is further supported at 3$\sigma$ significance when including literature GPI astrometry corrected for calibration systematics. Fits of the SPHERE and GPI data using three complementary approaches provide broadly similar results. The data suggest an orbital period of 32$^{+17}_{-9}$ yr (i.e. 12$^{+4}_{-2}$ au in semi-major axis), an inclination of 133$^{+14}_{-7}$ deg, an eccentricity of 0.45$^{+0.10}_{-0.15}$, and an argument of periastron passage of 87$^{+34}_{-30}$ deg [mod 180 deg]. The time at periastron passage and the longitude of node exhibit bimodal distributions because we do not detect yet if the planet is accelerating or decelerating along its orbit. Given the inclinations of the planet's orbit and of the stellar rotation axis (134-144 deg), we infer alignment or misalignment within 18 deg for the star-planet spin-orbit. Further astrometric monitoring in the next 3-4 years is required to confirm at a higher significance the curvature in the planet's motion, determine if the planet is accelerating or decelerating on its orbit, and further constrain its orbital parameters and the star-planet spin-orbit.

M stars are ideal targets to search for Earth-like planets. However, they usually have high levels of magnetic activity, which could affect their habitability and make difficult the detection of exoplanets orbiting around them. Unfortunately, long-term variability of dM stars has not been extensively studied, due to their low intrinsic brightness. For this reason, in 1999 we started the HKalpha project, which systematically observes the spectra of a large number of stars, in particular dM stars, at the Complejo Astronomico El Leoncito (CASLEO). In this work, we study the long-term activity of the young active dM1 star AU Microscopii. We analyze the Mount Wilson index S derived from CASLEO spectra obtained between 2004 and 2016, which we complement with the S-index derived from HARPS, FEROS and UVES public spectra. We also analyze the simultaneous photometric counterpart provided by the ASAS public database for this star between 2000 and 2009, and our own photometry. In both totally independent time series, we detect a possible activity cycle of period around 5 years. We also derived a precise rotation period for this star P rot = 4.85 days, consistent with the literature. This activity cycle reflects that an alphaOmega dynamo could be operating in this star.

We detect the accreting low-mass companion HD142527B at a separation of 73 mas (11.4 au) from the star. No other companions with mass greater than 10 MJ are visible in the field of view of IFS (\sim 100 au centered on the star) or in the IRDIS field of view (\sim 400 au centered on the star). Measurements from IFS, SAM IFS, and IRDIS suggest an M6 spectral type for HD142527B, with an uncertainty of one spectral subtype, compatible with an object of M=0.11 \pm 0.06 MSun and R=0.15 \pm 0.07 RSun. The determination of the mass remains a challenge using contemporary evolutionary models, as they do not account for the energy input due to accretion from infalling material. We consider that the spectral type of the secondary may also be earlier than the type we derived from IFS spectra. From dynamical considerations, we further constrain the mass to 0.26^{+0.16}_{-0.14} MSun , which is consistent with both our spectroscopic analysis and the values reported in the literature. Following previous methods, the lower and upper dynamical mass values correspond to a spectral type between M2.5 and M5.5 for the companion. By fitting the astrometric points, we find the following orbital parameters: a period of P=35-137 yr; an inclination of i=121-130 deg.; , a value of Omega=124-135 deg for the longitude of node, and an 68% confidence interval of \sim 18 - 57 au for the separation at periapsis. Eccentricity and time at periapsis passage exhibit two groups of values: \sim0.2-0.45 and \sim0.45-0.7 for e, and \sim 2015-2020 and \sim2020-2022 for T_0. While these orbital parameters might at first suggest that HD142527B is not the companion responsible for the outer disk truncation, a previous hydrodynamical analysis of this system showed that they are compatible with a companion that is able to produce the large cavity and other observed features.

With an orbital distance comparable to that of Saturn in the solar system, \bpic b is the closest (semi-major axis $\simeq$\,9\,au) exoplanet that has been imaged to orbit a star. Thus it offers unique opportunities for detailed studies of its orbital, physical, and atmospheric properties, and of disk-planet interactions. With the exception of the discovery observations in 2003 with NaCo at the Very Large Telescope (VLT), all following astrometric measurements relative to \bpic have been obtained in the southwestern part of the orbit, which severely limits the determination of the planet's orbital parameters. We aimed at further constraining \bpic b orbital properties using more data, and, in particular, data taken in the northeastern part of the orbit. We used SPHERE at the VLT to precisely monitor the orbital motion of beta \bpic b since first light of the instrument in 2014. We were able to monitor the planet until November 2016, when its angular separation became too small (125 mas, i.e., 1.6\,au) and prevented further detection. We redetected \bpic b on the northeast side of the disk at a separation of 139\,mas and a PA of 30$^{\circ}$ in September 2018. The planetary orbit is now well constrained. With a semi-major axis (sma) of $a = 9.0 \pm 0.5$ au (1 $\sigma $), it definitely excludes previously reported possible long orbital periods, and excludes \bpic b as the origin of photometric variations that took place in 1981. We also refine the eccentricity and inclination of the planet. From an instrumental point of view, these data demonstrate that it is possible to detect, if they exist, young massive Jupiters that orbit at less than 2 au from a star that is 20 pc away.

Young stars exhibit short-term photometric variability caused by mass accretion events from circumstellar disks, the presence of dusty warps within the inner disks, starspots that rotate across the stellar surfaces, and flares. Long-term variability also occurs owing to starspot longevity and cycles, and from changes in stellar angular momenta and activity as the stars age. We propose to observe the Carina star-forming region in different bands with a cadence of 30 minutes every night for one week per year to clarify the nature of both the short-term and long-term variability of the thousands of young stars in this region. By obtaining well-sampled multicolor lightcurves of this dense young cluster, LSST would acquire the first statistically significant data on how these objects vary on both short and long timescales. This information will allow us to relate the observed variability to stellar properties such as mass, age, binarity, and to environmental properties such as location within or exterior to the H II region, and to the presence or absence of a circumstellar disk.

One of the aims of LSST is to perform a systematic survey of star clusters and star forming regions (SFRs) in our Galaxy. In particular, the observations obtained with LSST will make a big difference in Galactic regions that have been poorly studied in the past, such as the anticenter and the disk beyond the Galactic center, and they will have a strong impact in discovering new distant SFRs. These results can be achieved by exploiting the exquisite depth that will be attained if the wide-fast-deep (WFD) observing strategy of the main survey is also adopted for the Galactic plane, in the g, r, and i filters.

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.

The Gaia Data Release 2 (DR2): we summarise the processing and results of the identification of variable source candidates of RR Lyrae stars, Cepheids, long period variables (LPVs), rotation modulation (BY Dra-type) stars, delta Scuti & SX Phoenicis stars, and short-timescale variables. In this release we aim to provide useful but not necessarily complete samples of candidates. The processed Gaia data consist of the G, BP, and RP photometry during the first 22 months of operations as well as positions and parallaxes. Various methods from classical statistics, data mining and time series analysis were applied and tailored to the specific properties of Gaia data, as well as various visualisation tools. The DR2 variability release contains: 228'904 RR Lyrae stars, 11'438 Cepheids, 151'761 LPVs, 147'535 stars with rotation modulation, 8'882 delta Scuti & SX Phoenicis stars, and 3'018 short-timescale variables. These results are distributed over a classification and various Specific Object Studies (SOS) tables in the Gaia archive, along with the three-band time series and associated statistics for the underlying 550'737 unique sources. We estimate that about half of them are newly identified variables. The variability type completeness varies strongly as function of sky position due to the non-uniform sky coverage and intermediate calibration level of this data. The probabilistic and automated nature of this work implies certain completeness and contamination rates which are quantified so that users can anticipate their effects. This means that even well-known variable sources can be missed or misidentified in the published data. The DR2 variability release only represents a small subset of the processed data. Future releases will include more variable sources and data products; however, DR2 shows the (already) very high quality of the data and great promise for variability studies.

We present observations with the planet finder SPHERE of a selected sample of the most promising radial velocity (RV) companions for high-contrast imaging. Using a Monte Carlo simulation to explore all the possible inclinations of the orbit of wide RV companions, we identified the systems with companions that could potentially be detected with SPHERE. We found the most favorable RV systems to observe are : HD\,142, GJ\,676, HD\,39091, HIP\,70849, and HD\,30177 and carried out observations of these systems during SPHERE Guaranteed Time Observing (GTO). To reduce the intensity of the starlight and reveal faint companions, we used Principle Component Analysis (PCA) algorithms alongside angular and spectral differential imaging. We injected synthetic planets with known flux to evaluate the self-subtraction caused by our data reduction and to determine the 5$\sigma$ contrast in the J band $vs$ separation for our reduced images. We estimated the upper limit on detectable companion mass around the selected stars from the contrast plot obtained from our data reduction. Although our observations enabled contrasts larger than 15 mag at a few tenths of arcsec from the host stars, we detected no planets. However, we were able to set upper mass limits around the stars using AMES-COND evolutionary models. We can exclude the presence of companions more massive than 25-28 \MJup around these stars, confirming the substellar nature of these RV companions.

We present the methods devised to identify the BY Dra variables candidates in Gaia DR2 and infer their variability parameters. BY Dra candidates are pre-selected from their position in the HR diagram, built from Gaia parallaxes, $G$ magnitudes, and $(G_{BP} - G_{RP})$ colours. Since the time evolution of the stellar active region can disrupt the coherence of the signal, segments not much longer than their expected evolution timescale are extracted from the entire photometric time-series and period search algorithms are applied to each segment. For the Gaia DR2, we select sources having similar period in at least two segments as candidates BY Dra. Results are further filtered considering the time series phase coverage and the expected approximate light curve shape. Gaia DR2 includes rotational periods and modulation amplitudes of 147 535 BY Dra candidates. The data unveil the existence of two populations with distinctive period and amplitude distributions. The sample covers 38% of the whole sky when divided in bins (HEALPix) of $\approx$0.84 square degrees and we estimate that represents 0.7 -- 5 % of all BY Dra stars potentially detectable by Gaia. The preliminary data contained in Gaia DR2 illustrate the vast and unique information that the mission is going to provide on stellar rotation and magnetic activity. This information, complemented by Gaia exquisite parallaxes, proper motions, and astrophysical parameter, is opening new and unique perspectives for our understanding of the evolution of stellar angular momentum and dynamo action.

Evolutionary models are widely used to infer the mass of stars, brown dwarfs, and giant planets. Their predictions are thought to be less reliable at young ages ($<$ 200 Myr) and in the low-mass regime ($\mathrm{<1~M_{\odot}}$). GJ 2060 AB and TWA 22 AB are two rare astrometric M-dwarf binaries respectively members of the AB Doradus and Beta Pictoris moving groups. As their dynamical mass can be measured within a few years, they can be used to calibrate the evolutionary tracks and set new constraints on the age of young moving groups. We find a total mass of $\mathrm{0.18\pm 0.02~M_\odot}$ for TWA 22. That mass is in good agreement with model predictions at the age of the Beta Pic moving group. We obtain a total mass of $\mathrm{1.09 \pm 0.10~M_{\odot}}$ for GJ 2060. We estimate a spectral type of M$1\pm0.5$, $\mathrm{L/L_{\odot}=-1.20\pm0.05}$ dex, and $\mathrm{T_{eff}=3700\pm100}$ K for GJ 2060 A. The B component is a M$3\pm0.5$ dwarf with $\mathrm{L/L_{\odot}=-1.63\pm0.05}$ dex and $\mathrm{T_{eff}=3400\pm100}$ K. The dynamical mass of GJ 2060 AB is inconsistent with the most recent models predictions (BCAH15, PARSEC) for an ABDor age in the range 50-150 Myr. It is 10 to 20\% (1-2 sigma, depending on the assumed age) above the models predictions, corresponding to an underestimation of $0.10$ to $0.20~\mathrm{M_\odot}$. Coevality suggests a young age for the system ($\sim$ 50 Myr) according to most evolutionary models. TWA 22 validates the predictions of recent evolutionary tracks at $\sim$20 Myr. On the other hand, we evidence a 1-2 sigma mismatch between the predicted and observed mass of GJ 2060 AB. This slight departure may indicate that one of the star hosts a tight companion. Alternatively, this would confirm the models tendency to underestimate the mass of young low-mass stars.

Context. Stellar activity is currently challenging the detection of young planets via the radial velocity (RV) technique. Aims. We attempt to definitively discriminate the nature of the RV variations for the young active K5 star BD+20 1790, for which visible (VIS) RV measurements show divergent results on the existence of a substellar companion. Methods. We compare VIS data with high precision RVs in the near infrared (NIR) range by using the GIANO - B and IGRINS spectrographs. In addition, we present for the first time simultaneous VIS-NIR observations obtained with GIARPS (GIANO - B and HARPS - N) at Telescopio Nazionale Galileo (TNG). Orbital RVs are achromatic, so the RV amplitude does not change at different wavelengths, while stellar activity induces wavelength-dependent RV variations, which are significantly reduced in the NIR range with respect to the VIS. Results. The NIR radial velocity measurements from GIANO - B and IGRINS show an average amplitude of about one quarter with respect to previously published VIS data, as expected when the RV jitter is due to stellar activity. Coeval multi-band photometry surprisingly shows larger amplitudes in the NIR range, explainable with a mixture of cool and hot spots in the same active region. Conclusions. In this work, the claimed massive planet around BD+20 1790 is ruled out by our data. We exploited the crucial role of multi- wavelength spectroscopy when observing young active stars: thanks to facilities like GIARPS that provide simultaneous observations, this method can reach its maximum potential.

Debris disks are usually detected through the infrared excess over the photospheric level of their host star. The most favorable stars for disk detection are those with spectral types between A and K, while the statistics for debris disks detected around low-mass M-type stars is very low, either because they are rare or because they are more difficult to detect. Terrestrial planets, on the other hand, may be common around M-type stars. Here, we report on the discovery of an extended (likely) debris disk around the M-dwarf GSC 07396-00759. The star is a wide companion of the close accreting binary V4046 Sgr. The system probably is a member of the $\beta$ Pictoris Moving Group. We resolve the disk in scattered light, exploiting high-contrast, high-resolution imagery with the two near-infrared subsystems of the VLT/SPHERE instrument, operating in the YJ bands and the H2H3 doublet. The disk is clearly detected up to 1.5" ($\sim110$ au) from the star and appears as a ring, with an inclination $i\sim83$ degree, and a peak density position at $\sim 70$ au. The spatial extension of the disk suggests that the dust dynamics is affected by a strong stellar wind, showing similarities with the AU Mic system that has also been resolved with SPHERE. The images show faint asymmetric structures at the widest separation in the northwest side. We also set an upper limit for the presence of giant planets to $2 M_J$. Finally, we note that the 2 resolved disks around M-type stars of 30 such stars observed with SPHERE are viewed close to edge-on, suggesting that a significant population of debris disks around M dwarfs could remain undetected because of an unfavorable orientation.

Stellar activity is the ultimate source of radial-velocity (RV) noise in the search for Earth-mass planets orbiting late-type main-sequence stars. We analyse the performance of four different indicators and the chromospheric index $\log R'_{\rm HK}$ in detecting RV variations induced by stellar activity in 15 slowly rotating ($v\sin i \leq 5$ km/s), weakly active ($\log R'_{\rm HK} \leq -4.95$) solar-like stars observed with the high-resolution spectrograph HARPS-N. We consider indicators of the asymmetry of the cross-correlation function (CCF) between the stellar spectrum and the binary weighted line mask used to compute the RV, that is the bisector inverse span (BIS), $\Delta V$, and a new indicator $V_{\rm asy(mod)}$ together with the full width at half maximum (FWHM) of the CCF. We present methods to evaluate the uncertainties of the CCF indicators and apply a kernel regression (KR) between the RV, the time, and each of the indicators to study their capability of reproducing the RV variations induced by stellar activity. The considered indicators together with the KR prove to be useful to detect activity-induced RV variations in $47 \pm 18$ percent of the stars over a two-year time span when a significance (two-sided p-value) threshold of one percent is adopted. In those cases, KR reduces the standard deviation of the RV time series by a factor of approximately two. The BIS, the FWHM, and the newly introduced $V_{\rm asy(mod)}$ are the best indicators, being useful in $27 \pm 13$, $13 \pm 9$, and $13 \pm 9$ percent of the cases, respectively. The relatively limited performances of the activity indicators are related to the very low activity level and $v\sin i$ of the considered stars. For the application of our approach to sun-like stars, a spectral resolution of at least $10^5$ and highly stabilized spectrographs are recommended.

Transition disks are considered sites of ongoing planet formation, and their dust and gas distributions could be signposts of embedded planets. The transition disk around the T Tauri star RY Lup has an inner dust cavity and displays a strong silicate emission feature. Using high-resolution imaging we study the disk geometry, including non-axisymmetric features, and its surface dust grain, to gain a better understanding of the disk evolutionary process. Moreover, we search for companion candidates, possibly connected to the disk. We obtained high-contrast and high angular resolution data in the near-infrared with the VLT/SPHERE extreme adaptive optics instrument whose goal is to study the planet formation by detecting and characterizing these planets and their formation environments through direct imaging. We performed polarimetric imaging of the RY~Lup disk with IRDIS (at 1.6 microns), and obtained intensity images with the IRDIS dual-band imaging camera simultaneously with the IFS spectro-imager (0.9-1.3 microns). We resolved for the first time the scattered light from the nearly edge-on circumstellar disk around RY~Lup, at projected separations in the 100 \,au range. The shape of the disk and its sharp features are clearly detectable at wavelengths ranging from 0.9 to 1.6 microns. We show that the observed morphology can be interpreted as spiral arms in the disk. This interpretation is supported by in-depth numerical simulations. We also demonstrate that these features can be produced by one planet interacting with the disk. We also detect several point sources which are classified as probable background objects.

HD95086 (A8V, 17Myr) hosts a rare planetary system for which a multi-belt debris disk and a giant planet of 4-5MJup have been directly imaged. Our study aims to characterize the physical and orbital properties of HD95086b, search for additional planets at short and wide orbits and image the cold outer debris belt in scattered light. We used HARPS at the ESO 3.6m telescope to monitor the radial velocity of HD95086 over 2 years and investigate the existence of giant planets at less than 3au orbital distance. With the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE at VLT, we imaged the faint circumstellar environment beyond 10au at six epochs between 2015 and 2017. We do not detect additional giant planets around HD95086. We identified the nature (bound companion or background contaminant) of all point-like sources detected in the IRDIS field of view. None of them correspond to the ones recently discovered near the edge of the cold outer belt by ALMA. HD95086b is resolved for the first time in J-band with IFS. Its near-infrared spectral energy distribution is well fitted by a few dusty and/or young L7-L9 dwarf spectral templates. The extremely red 1-4um spectral distribution is typical of low-gravity objects at the L/T spectral type transition. The planet's orbital motion is resolved between January 2015 and May 2017. Together with past NaCo measurements properly re-calibrated, our orbital fitting solutions favor a retrograde low to moderate-eccentricity orbit e=0.2 (0.0 to 0.5), with a semi-major axis 52au corresponding to orbital periods of 288$ yrs and an inclination that peaks at i = 141deg, which is compatible with a planet-disk coplanar configuration. Finally, we report the detection in polarimetric differential imaging of the cold outer debris belt between 100 and 300au, consistent in radial extent with recent ALMA 1.3mm resolved observations.

The SHINE program is a large high-contrast near-infrared survey of 600 young, nearby stars. It is aimed at searching for and characterizing new planetary systems using VLT/SPHERE's unprecedented high-contrast and high-angular resolution imaging capabilities. It also intends at placing statistical constraints on the occurrence and orbital properties of the giant planet population at large orbits as a function of the stellar host mass and age to test planet formation theories. We use the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images and spectra of the young A2 star HIP65426. It is a member of the ~17 Myr old Lower Centaurus-Crux association. At a separation of 830 mas (92 au projected) from the star, we detect a faint red companion. Multi-epoch observations confirm that it shares common proper motion with HIP65426. Spectro-photometric measurements extracted with IFS and IRDIS between 0.95 and 2.2um indicate a warm, dusty atmosphere characteristic of young low surface-gravity L5-L7 dwarfs. Hot-start evolutionary models predict a luminosity consistent with a 6-12 MJup, Teff=1300-1600 K and R=1.5 RJup giant planet. Finally, the comparison with Exo-REM and PHOENIX BT-Settl synthetic atmosphere models gives consistent effective temperatures but with slightly higher surface gravity solutions of log(g)=4.0-5.0 with smaller radii (1.0-1.3 RJup). Given its physical and spectral properties, HIP65426b occupies a rather unique placement in terms of age, mass and spectral-type among the currently known imaged planets. It represents a particularly interesting case to study the presence of clouds as a function of particle size, composition, and location in the atmosphere, to search for signatures of non-equilibrium chemistry, and finally to test the theory of planet formation and evolution.

AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory designed to extend its energy range of detection and to directly measure the muon content of the cosmic ray primary particle showers. The array will be formed by an infill of surface water-Cherenkov detectors associated with buried scintillation counters employed for muon counting. Each counter is composed of three scintillation modules, with a 10 m$^2$ detection area per module. In this paper, a new generation of detectors, replacing the current multi-pixel photomultiplier tube (PMT) with silicon photo sensors (aka. SiPMs), is proposed. The selection of the new device and its front-end electronics is explained. A method to calibrate the counting system that ensures the performance of the detector is detailed. This method has the advantage of being able to be carried out in a remote place such as the one where the detectors are deployed. High efficiency results, i.e. 98 % efficiency for the highest tested overvoltage, combined with a low probability of accidental counting ($\sim$2 %), show a promising performance for this new system.

The substellar companion HD206893b has recently been discovered by direct imaging of its disc-bearing host star with the SPHERE instrument. We investigate the atypical properties of the companion, which has the reddest near-infrared colours among all known substellar objects, either orbiting a star or isolated, and we provide a comprehensive characterisation of the host star-disc-companion system. We conducted a follow-up of the companion with adaptive optics imaging and spectro-imaging with SPHERE, and a multiinstrument follow-up of its host star. We obtain a R=30 spectrum from 0.95 to 1.64 micron of the companion and additional photometry at 2.11 and 2.25 micron. We carried out extensive atmosphere model fitting for the companions and the host star in order to derive their age, mass, and metallicity. We found no additional companion in the system in spite of exquisite observing conditions resulting in sensitivity to 6MJup (2MJup) at 0.5" for an age of 300 Myr (50 Myr). We detect orbital motion over more than one year and characterise the possible Keplerian orbits. We constrain the age of the system to a minimum of 50 Myr and a maximum of 700 Myr, and determine that the host-star metallicity is nearly solar. The comparison of the companion spectrum and photometry to model atmospheres indicates that the companion is an extremely dusty late L dwarf, with an intermediate gravity (log g 4.5-5.0) which is compatible with the independent age estimate of the system. Though our best fit corresponds to a brown dwarf of 15-30 MJup aged 100-300 Myr, our analysis is also compatible with a range of masses and ages going from a 50 Myr 12MJup planetary-mass object to a 50 MJup Hyades-age brown dwarf...

In this paper we present the results of the SPHERE observation of the HD 284149 system, aimed at a more detailed characterisation of both the primary and its brown dwarf companion. We observed HD 284149 in the near-infrared with SPHERE, using the imaging mode (IRDIS+IFS) and the long-slit spectroscopy mode (IRDIS-LSS). The data were reduced using the dedicated SPHERE pipeline, and algorithms such as PCA and TLOCI were applied to reduce the speckle pattern. The IFS images revealed a previously unknown low-mass (~0.16$M_{\odot}$) stellar companion (HD 294149 B) at ~0.1$^{\prime\prime}$, compatible with previously observed radial velocity differences, as well as proper motion differences between Gaia and Tycho-2 measurements. The known brown dwarf companion (HD 284149 b) is clearly visible in the IRDIS images. This allowed us to refine both its photometry and astrometry. The analysis of the medium resolution IRDIS long slit spectra also allowed a refinement of temperature and spectral type estimates. A full reassessment of the age and distance of the system was also performed, leading to more precise values of both mass and semi-major axis. As a result of this study, HD 284149 ABb therefore becomes the latest addition to the (short) list of brown dwarfs on wide circumbinary orbits, providing new evidence to support recent claims that object in such configuration occur with a similar frequency to wide companions to single stars.

Low-mass members of young stellar associations exhibit a wide spread of rotation periods. Such a spread originates from distributions of masses and initial rotation periods. However, multiplicity can also play a significant role. We investigate the role played by physical companions in shortening the primordial disc lifetime. We have compiled the most extensive list of low-mass members of the young 25-Myr beta Pictoris association. We have measured the rotation periods of about all members and used updated UVWXYZ components to assess their membership. We built the rotation period distribution distinguishing between bona fide members and candidate members and according to their multiplicity status. We found that single stars and components of multiple systems in wide orbits (>80 AU) have rotation periods that exhibit a well defined sequence arising from mass distribution. All components of multiple systems in close orbits (<80 AU) have rotation periods significantly shorter than their equal-mass single counterparts. A comparison with the younger 13 Myr h Per cluster and with the older 40-Myr open clusters/stellar associations NGC2547, IC2391, Argus, and IC2602 and the 130-Myr Pleiades shows that whereas the evolution of F-G stars is well reproduced by angular momentum evolution models, this is not the case for the slow K and early-M stars. Finally, we found that the amplitude of their light curves is correlated neither with rotation nor with mass. Once single stars and wide components of multiple systems are separated from close components of multiple systems, the rotation period distributions exhibit a well defined dependence on mass that allows to make a meaningful comparison with similar distributions of either younger or older associations/clusters. Such cleaned distributions allow to use the stellar rotation period as age indicator, meaningfully for F and G type stars.

Magnetic cycles have been detected in tens of solar-like stars. The relationship between the cycle properties and global stellar parameters is not fully understood yet. We searched for activity cycles in 90 solar-like stars with ages between 4 and 95 Myr aiming to investigate the properties of activity cycles in this age range. We measured the length $P_{ cyc}$ of a given cycle by analyzing the long-term time-series of three activity indexes. For each star, we computed also the global magnetic activity index <IQR> that is proportional to the amplitude of the rotational modulation and is a proxy of the mean level of the surface magnetic activity. We detected activity cycles in 67 stars. Secondary cycles were also detected in 32 stars. The lack of correlation between $P_{ cyc}$ and $P_{ rot}$ suggest that these stars belong to the Transitional Branch and that the dynamo acting in these stars is different from the solar one. This statement is also supported by the analysis of the butterfly diagrams. We computed the Spearman correlation coefficient $r_{ S}$ between $P_{ cyc}$, <IQR> and different stellar parameters. We found that $P_{ cyc}$ is uncorrelated with all the investigated parameters. The <IQR> index is positively correlated with the convective turn-over time-scale, the magnetic diffusivity time-scale $\tau_{ diff}$, and the dynamo number $D_{ N}$, whereas it is anti-correlated with the effective temperature $T_{ eff}$, the photometric shear $\Delta\Omega_{\rm phot}$ and the radius $R_{ C}$ at which the convective zone is located. We found that $P_{ cyc}$ is about constant and that <IQR> decreases with the stellare age in the range 4-95 Myr. We investigated the magnetic activity of AB Dor A by merging ASAS time-series with previous long-term photometric data. We estimated the length of the AB Dor A primary cycle as $P_{ cyc} = 16.78 \pm 2 \rm yr$.

We report a multi-resolution search for anisotropies in the arrival directions of cosmic rays detected at the Pierre Auger Observatory with local zenith angles up to $80^\circ$ and energies in excess of 4 EeV ($4 \times 10^{18}$ eV). This search is conducted by measuring the angular power spectrum and performing a needlet wavelet analysis in two independent energy ranges. Both analyses are complementary since the angular power spectrum achieves a better performance in identifying large-scale patterns while the needlet wavelet analysis, considering the parameters used in this work, presents a higher efficiency in detecting smaller-scale anisotropies, potentially providing directional information on any observed anisotropies. No deviation from isotropy is observed on any angular scale in the energy range between 4 and 8 EeV. Above 8 EeV, an indication for a dipole moment is captured; while no other deviation from isotropy is observed for moments beyond the dipole one. The corresponding $p$-values obtained after accounting for searches blindly performed at several angular scales, are $1.3 \times 10^{-5}$ in the case of the angular power spectrum, and $2.5 \times 10^{-3}$ in the case of the needlet analysis. While these results are consistent with previous reports making use of the same data set, they provide extensions of the previous works through the thorough scans of the angular scales.

Context. The SAO 206462 (HD 135344B) disk is one of the few known transitional disks showing asymmetric features in scattered light and thermal emission. Near-infrared scattered-light images revealed two bright outer spiral arms and an inner cavity depleted in dust. Giant protoplanets have been proposed to account for the disk morphology. Aims. We aim to search for giant planets responsible for the disk features and, in the case of non-detection, to constrain recent planet predictions using the data detection limits. Methods. We obtained new high-contrast and high-resolution total intensity images of the target spanning the Y to the K bands (0.95-2.3 mic) using the VLT/SPHERE near-infrared camera and integral field spectrometer. Results. The spiral arms and the outer cavity edge are revealed at high resolutions and sensitivities without the need for image post-processing techniques, which introduce photometric biases. We do not detect any close-in companions. For the derivation of the detection limits on putative giant planets embedded in the disk, we show that the knowledge of the disk aspect ratio and viscosity is critical for the estimation of the attenuation of a planet signal by the protoplanetary dust because of the gaps that these putative planets may open. Given assumptions on these parameters, the mass limits can vary from ~2-5 to ~4-7 Jupiter masses at separations beyond the disk spiral arms. The SPHERE detection limits are more stringent than those derived from archival NaCo/L' data and provide new constraints on a few recent predictions of massive planets (4-15 MJ) based on the spiral density wave theory. The SPHERE and ALMA data do not favor the hypotheses on massive giant planets in the outer disk (beyond 0.6). There could still be low-mass planets in the outer disk and/or planets inside the cavity.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Thanks to the detections of more than 3000 exoplanets these last 20 years, statistical studies have already highlighted some properties in the distribution of the planet parameters. Nevertheless, few studies have yet investigated the planet populations from short to large separations around the same star since this requires the use of different detection techniques that usually target different types of stars. We wish to develop a tool that combines direct and indirect methods so as to correctly investigate the giant planet populations at all separations. We developed the MESS2 code, a Monte Carlo simulation code combining radial velocity and direct imaging data obtained at different epochs for a given star to estimate the detection probability of giant planets spanning a wide range of physical separations. It is based on the generation of synthetic planet populations. We apply MESS2 on a young M1-type, the nearby star AUMic observed with HARPS and NACO/ESO. We show that giant planet detection limits are significantly improved at intermediate separations (~20au in the case of AUMic). We show that the traditional approach of analysing independently the RV and DI detection limits systematically overestimates the planet detection limits and hence planet occurrence rates. The use of MESS2 allows to obtain correct planet occurrence rates in statistical studies, making use of multi-epoch DI data and/or RV measurements. We also show that MESS2 can optimise the schedule of future DI observations.

Understanding the formation and evolution of giant planets ($\ge$1 $M_{Jup}$) at wide orbital separation ($\ge$5 AU) is one of the goals of direct imaging. Over the past 15 years, many surveys have placed strong constraints on the occurrence rate of wide-orbit giants, mostly based on non-detections, but very few have tried to make a direct link with planet formation theories. In the present work, we combine the results of our previously published VLT/NaCo large program with the results of 12 past imaging surveys to constitute a statistical sample of 199 FGK stars within 100 pc, including 3 stars with sub-stellar companions. Using Monte Carlo simulations and assuming linear flat distributions for the mass and semi-major axis of planets, we estimate the sub-stellar companion frequency to be within 0.75-5.7% at the 68% confidence level (CL) within 20-300 AU and 0.5-75 $M_{Jup}$, which is compatible with previously published results. We also compare our results with the predictions of state-of-the-art population synthesis models based on the gravitational instability (GI) formation scenario by Forgan & Rice (2013), with and without scattering. We estimate that in both the scattered and non-scattered populations, we would be able to detect more than 30% of companions in the 1-75 $M_{Jup}$ range (95% CL). With the 3 sub-stellar detections in our sample, we estimate the fraction of stars that host a planetary system formed by GI to be within 1.0-8.6% (95% CL). We also conclude that even though GI is not common, it predicts a mass distribution of wide-orbit massive companions that is much closer to what is observed than what the core accretion scenario predicts. Finally, we associate the present paper with the release of the Direct Imaging Virtual Archive (DIVA, http://cesam.lam.fr/diva/), a public database that aims at gathering the results of past, present, and future direct imaging surveys.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

The ESA Gaia mission provides a unique time-domain survey for more than one billion sources brighter than G=20.7 mag. Gaia offers the unprecedented opportunity to study variability phenomena in the Universe thanks to multi-epoch G-magnitude photometry in addition to astrometry, blue and red spectro-photometry, and spectroscopy. Within the Gaia Consortium, Coordination Unit 7 has the responsibility to detect variable objects, classify them, derive characteristic parameters for specific variability classes, and provide global descriptions of variable phenomena. We describe the variability processing and analysis that we plan to apply to the successive data releases, and we present its application to the G-band photometry results of the first 14 months of Gaia operations that comprises 28 days of Ecliptic Pole Scanning Law and 13 months of Nominal Scanning Law. Out of the 694 million, all-sky, sources that have calibrated G-band photometry in this first stage of the mission, about 2.3 million sources that have at least 20 observations are located within 38 degrees from the South Ecliptic Pole. We detect about 14% of them as variable candidates, among which the automated classification identified 9347 Cepheid and RR Lyrae candidates. Additional visual inspections and selection criteria led to the publication of 3194 Cepheid and RR Lyrae stars, described in Clementini et al. (2016). Under the restrictive conditions for DR1, the completenesses of Cepheids and RR Lyrae stars are estimated at 67% and 58%, respectively, numbers that will significantly increase with subsequent Gaia data releases. Data processing within the Gaia Consortium is iterative, the quality of the data and the results being improved at each iteration. The results presented in this article show a glimpse of the exceptional harvest that is to be expected from the Gaia mission for variability phenomena. [abridged]

In case of pre-main sequence objects, the only way to determine age and mass is by fitting theoretical isochrones on color-magnitude (alternatively luminosity-temperature) diagrams. Since young stellar objects exhibit photometric variability over wide range in magnitude and colors, the age and mass determined by fitting isochrones is expected to be inaccurate, if not erroneous. These in turn will badly affect any study carried out on age spread and process of star formation. Since we have carried out very extensive photometric observations of the Orion Nebula Cluster (ONC), we decided to use our multi-band data to explore the influence of variability in determining mass and age of cluster members. In this study, we get the amplitudes of the photometric variability in V, R, and I optical bands of a sample of 346 ONC members and use it to investigate how the variability affects the inferred masses and ages and if it alone can take account for the age spread among the ONC members reported by earlier studies. We find that members that show periodic and smooth photometric rotational modulation have their masses and ages unaffected by variability. On other hand, we found that members with periodic but very scattered photometric rotational modulation and members with irregular variability have their masses and ages significantly affected. Moreover, using Hertzsprung-Russell (HR) diagrams we find that the observed I band photometric variability can take account of only a fraction (about 50%) of the inferred age spread, whereas the V band photometric variability is large enough to mask any age spread.

On September 14, 2015 the Advanced LIGO detectors observed their first gravitational-wave (GW) transient GW150914. This was followed by a second GW event observed on December 26, 2015. Both events were inferred to have arisen from the merger of black holes in binary systems. Such a system may emit neutrinos if there are magnetic fields and disk debris remaining from the formation of the two black holes. With the surface detector array of the Pierre Auger Observatory we can search for neutrinos with energy above 100 PeV from point-like sources across the sky with equatorial declination from about -65 deg. to +60 deg., and in particular from a fraction of the 90% confidence-level (CL) inferred positions in the sky of GW150914 and GW151226. A targeted search for highly-inclined extensive air showers, produced either by interactions of downward-going neutrinos of all flavors in the atmosphere or by the decays of tau leptons originating from tau-neutrino interactions in the Earth's crust (Earth-skimming neutrinos), yielded no candidates in the Auger data collected within $\pm 500$ s around or 1 day after the coordinated universal time (UTC) of GW150914 and GW151226, as well as in the same search periods relative to the UTC time of the GW candidate event LVT151012. From the non-observation we constrain the amount of energy radiated in ultrahigh-energy neutrinos from such remarkable events.

There are a variety of different techniques available to estimate the ages of pre-main-sequence stars. Components of physical pairs, thanks to their strict coevality and the mass difference, such as the binary system analysed in this paper, are best suited to test the effectiveness of these different techniques. We consider the system WW Psa + TX Psa whose membership of the 25-Myr beta Pictoris association has been well established by earlier works. We investigate which age dating technique provides the best agreement between the age of the system and that of the association. We have photometrically monitored WW Psa and TX Psa and measured their rotation periods as P = 2.37d and P = 1.086d, respectively. We have retrieved from the literature their Li equivalent widths and measured their effective temperatures and luminosities. We investigate whether the ages of these stars derived using three independent techniques are consistent with the age of the beta Pictoris association. We find that the rotation periods and the Li contents of both stars are consistent with the distribution of other bona fide members of the cluster. On the contrary, the isochronal fitting provides similar ages for both stars, but a factor of about four younger than the quoted age of the association, or about 30% younger when the effects of magnetic fields are included. We explore the origin of the discrepant age inferred from isochronal fitting, including the possibilities that either the two components may be unresolved binaries or that the basic stellar parameters of both components are altered by enhanced magnetic activity. The latter is found to be the more reasonable cause, suggesting that age estimates based on the Li content is more reliable than isochronal fitting for pre-main-sequence stars with pronounced magnetic activity.

We intended to compile the most complete catalog of bona fide members and candidate members of the beta Pictoris association, and to measure their rotation periods and basic properties from our own observations, public archives, and exploring the literature. We carried out a multi-observatories campaign to get our own photometric time series and collected all archived public photometric data time series for the stars in our catalog. Each time series was analyzed with the Lomb-Scargle and CLEAN periodograms to search for the stellar rotation periods. We complemented the measured rotational properties with detailed information on multiplicity, membership, and projected rotational velocity available in the literature and discussed star by star. We measured the rotation periods of 112 out of 117 among bona fide members and candidate members of the beta Pictoris association and, whenever possible, we also measured the luminosity, radius, and inclination of the stellar rotation axis. This represents to date the largest catalog of rotation periods of any young loose stellar association. We provided an extensive catalog of rotation periods together with other relevant basic properties useful to explore a number of open issues, such as the causes of spread of rotation periods among coeval stars, evolution of angular momentum, and lithium-rotation connection.

In order to understand the observed physical and orbital diversity of extrasolar planetary systems, a full investigation of these objects and of their host stars is necessary. Within this field, one of the main purposes of the GAPS observing project with HARPS-N@TNG is to provide a more detailed characterisation of already known systems. In this framework we monitored the star, hosting two giant planets, HD108874, with HARPS-N for three years in order to refine the orbits, to improve the dynamical study and to search for additional low-mass planets in close orbits. We subtracted the radial velocity (RV) signal due to the known outer planets, finding a clear modulation of 40.2 d period. We analysed the correlation between RV residuals and the activity indicators and modelled the magnetic activity with a dedicated code. Our analysis suggests that the 40.2 d periodicity is a signature of the rotation period of the star. A refined orbital solution is provided, revealing that the system is close to a mean motion resonance of about 9:2, in a stable configuration over 1 Gyr. Stable orbits for low-mass planets are limited to regions very close to the star or far from it. Our data exclude super-Earths with Msin i \gtrsim 5 M_Earth within 0.4 AU and objects with Msin i \gtrsim 2 M_Earth with orbital periods of a few days. Finally we put constraints on the habitable zone of the system, assuming the presence of an exomoon orbiting the inner giant planet.

We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the electromagnetic and muonic components. It allows an accurate determination of the spread of primary masses in the cosmic-ray flux. Up till now, constraints on the spread of primary masses have been dominated by systematic uncertainties. The present correlation measurement is not affected by systematics in the measurement of the depth of shower maximum or the signal in the water Cherenkov stations. The analysis relies on general characteristics of air showers and is thus robust also with respect to uncertainties in hadronic event generators. The observed correlation in the energy range around the `ankle' at $\lg(E/{\rm eV})=18.5-19.0$ differs significantly from expectations for pure primary cosmic-ray compositions. A light composition made up of proton and helium only is equally inconsistent with observations. The data are explained well by a mixed composition including nuclei with mass $A > 4$. Scenarios such as the proton dip model, with almost pure compositions, are thus disfavoured as the sole explanation of the ultrahigh-energy cosmic-ray flux at Earth.

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

(Abridged) Understanding stellar activity in M dwarfs is crucial for the physics of stellar atmospheres as well as for ongoing radial velocity exoplanet programmes. Despite the increasing interest in M dwarfs, our knowledge of the chromospheres of these stars is far from being complete. We aim to test whether the relations between activity, rotation, and stellar parameters and flux-flux relationships also hold for early-M dwarfs on the main-sequence. We analyse in an homogeneous and coherent way a well defined sample of 71 late-K/early-M dwarfs that are currently being observed in the framework of the HArps-n red Dwarf Exoplanet Survey (HADES). Rotational velocities are derived using the cross-correlation technique while emission flux excesses in the Ca II H & K and Balmer lines from Halpha up to Hepsilon are obtained by using the spectral subtraction technique. The relationships between the emission excesses and the stellar parameters are studied. Relations between pairs of fluxes of different chromospheric lines are also studied. We find that the strength of the chromospheric emission in the Ca II H & K and Balmer lines is roughly constant for stars in the M0-M3 spectral range. Our data suggest that a moderate but significant correlation between activity and rotation might be present as well as a hint of kinematically selected young stars showing higher levels of emission. We find our sample of M dwarfs to be complementary in terms of chromospheric and X-ray fluxes with those of the literature, extending the analysis of the flux-flux relationships to the very low flux domain. Our results agree with previous works suggesting that the activity-rotation-age relationship known to hold for solar-type stars also applies to early-M dwarfs. We also confirm previous findings that the field stars which deviate from the bulk of the empirical flux-flux relationships show evidence of youth.

M dwarfs are prime targets for planet search programs, particularly of those focused on the detection and characterization of rocky planets in the habitable zone. Understanding their magnetic activity is important because it affects our ability to detect small planets, and it plays a key role in the characterization of the stellar environment. We analyze observations of the Ca II H&K and H{\alpha} lines as diagnostics of chromospheric activity for low-activity early-type M dwarfs. We analyze the time series of spectra of 71 early-type M dwarfs collected for the HADES project for planet search purposes. The HARPS-N spectra provide simultaneously the H&K doublet and the H{\alpha} line. We develop a reduction scheme able to correct the HARPS-N spectra for instrumental and atmospheric effects, and to provide flux-calibrated spectra in units of flux at the stellar surface. The H&K and H{\alpha} fluxes are compared with each other, and their variability is analyzed. We find that the H and K flux excesses are strongly correlated with each other, while the H{\alpha} flux excess is generally less correlated with the H&K doublet. We also find that H{\alpha} emission does not increase monotonically with the H&K line flux, showing some absorption before being filled in by chromospheric emission when H&K activity increases. Analyzing the time variability of the emission fluxes, we derive a tentative estimate of the rotation period (of the order of a few tens of days) for some of the program stars, and the typical lifetime of chromospheric active regions (a few stellar rotations). Our results are in good agreement with previous studies. In particular, we find evidence that the chromospheres of early-type M dwarfs could be characterized by different filaments coverage, affecting the formation mechanism of the H{\alpha} line. We also show that chromospheric structure is likely related to spectral type.

We search for signs of ongoing planet-disk interaction and study the distribution of small grains at the surface of the transition disk around RXJ1615.3-3255 (RX J1615). We observed RXJ1615 with VLT/SPHERE. We image the disk for the first time in scattered light and detect two arcs, two rings, a gap and an inner disk with marginal evidence for an inner cavity. The shapes of the arcs suggest that they probably are segments of full rings. Ellipse fitting for the two rings and inner disk yield a disk inclination i = 47 \pm 2 degrees and find semi-major axes of 1.50 \pm 0.01" (278 au), 1.06 \pm 0.01" (196 au) and 0.30 \pm 0.01" (56 au), respectively. We determine the scattering surface height above the midplane, based on the projected ring center offsets. Nine point sources are detected between 2.1" and 8.0" separation and considered as companion candidates. With NACO data we recover four of the nine point sources, which we determine not to be co-moving, and therefore unbound to the system. We present the first detection of the transition disk of RXJ1615 in scattered light. The height of the rings indicate limited flaring of the disk surface, which enables partial self-shadowing in the disk. The outermost arc either traces the bottom of the disk or it is another ring with semi-major axis > 2.35" (435 au). We explore both scenarios, extrapolating the complete shape of the feature, which will allow to distinguish between the two in future observations. The most interesting scenario, where the arc traces the bottom of the outer ring, requires the disk truncated at r ~ 360 au. The closest companion candidate, if indeed orbiting the disk at 540 au, would then be the most likely cause for such truncation. This companion candidate, as well as the remaining four, require follow up observations to determine if they are bound to the system.

We present a search for ultra-relativistic magnetic monopoles with the Pierre Auger Observatory. Such particles, possibly a relic of phase transitions in the early universe, would deposit a large amount of energy along their path through the atmosphere, comparable to that of ultrahigh-energy cosmic rays (UHECRs). The air shower profile of a magnetic monopole can be effectively distinguished by the fluorescence detector from that of standard UHECRs. No candidate was found in the data collected between 2004 and 2012, with an expected background of less than 0.1 event from UHECRs. The corresponding 90% confidence level (C.L.) upper limits on the flux of ultra-relativistic magnetic monopoles range from $10^{-19}$ (cm$^{2}$ sr s)$^{-1}$ for a Lorentz factor $\gamma=10^9$ to $2.5 \times10^{-21}$ (cm$^{2}$ sr s)$^{-1}$ for $\gamma=10^{12}$. These results - the first obtained with a UHECR detector - improve previously published limits by up to an order of magnitude.

The recent development of brand new observational techniques and theoretical models have greatly advanced the exoplanet research field. Despite significant achievements, which have allowed the detection of thousands extrasolar systems, a comprehensive understanding of planetary formation and evolution mechanisms is still desired. One relevant limitation is given by the accuracy in the measurements of planet-host star ages. The star GJ 504 has been found to host a substellar companion whose nature is strongly debated. There has been a recent difference of opinion in the literature owing to the uncertainty on the age of the system: a young age of $\sim$ 160 Myr would imply a giant planet as a companion, but a recent revision pointing to a solar age ($\sim$ 4 Gyr) instead suggests a brown dwarf. With the aim of shedding light on this debated topic, we have carried out a high-resolution spectroscopic study of GJ 504 to derive stellar parameters, metallicity, and abundances of both light and heavy elements, providing a full chemical characterisation. The main objective is to infer clues on the evolutionary stage (hence the age) of this system. We performed a strictly differential (line-by-line) analysis of GJ 504 with respect to two reference stars, that is the planet-host dwarf $\iota$ Hor and the subgiant HIP 84827. The former is crucial in this context because its stellar parameters (hence the evolutionary stage) is well constrained from asteroseismic observations. Regardless of the zero point offsets, our differential approach allows us to put tight constraints on the age of GJ 504 with respect to $\iota$ Hor, thereby minimising the internal uncertainties. We found that the surface gravity of GJ 504 is 0.2 $\pm$ 0.07 dex lower than that of the main-sequence star $\iota$ Hor, suggesting a past turn-off evolution[...]

Equal-mass stars in young open clusters and loose associations exhibit a wide spread of rotation periods, which likely originates from differences in the initial rotation periods and in the primordial disc lifetimes. We want to explore if the gravitational effects by nearby companions may play an additional role in producing the observed rotation period spread. We measure the photometric rotation periods of components of multiple stellar systems and look for correlations of the period differences among the components to their reciprocal distances. In this paper, we analysed the triple system AU Mic + AT Mic A&B in the 25-Myr beta Pictoris Association. We have retrieved from the literature the rotation period of AU Mic (P = 4.85d) and measured from photometric archival data the rotation periods of both components of AT Mic (P = 1.19d and P = 0.78d) for the first time. Moreover, we detected a high rate of flare events from AT Mic. Whereas the distant component AU Mic has evolved rotationally as a single star, the A and B components of AT Mic, separated by about 27 AU, exhibit a rotation rate a factor 5 larger than AU Mic. Moreover, the A and B components, despite have about equal mass, show a significant difference (about 40%) between their rotation periods. A possible explanation is that the gravitational forces between the A and B components of AT Mic (that are a factor about 7.3 x 10^6 more intense than those between AU Mic and AT Mic) have enhanced the dispersal of the AT Mic primordial disc, shortening its lifetime and the disc-locking phase duration, making the component A and B of AT Mic to rotate faster than the more distant AU Mic. We suspect that a different level of magnetic activity between the A and B components of AT Mic may be the additional parameter responsible for the difference between their rotation periods.

There is evidence in the 125-Myr Pleiades cluster, and more recently in the 5-Myr NGC 2264 cluster, that rotation plays a key role in the Lithium (Li) depletion processes among low-mass stars. Fast rotators appear to be less Li-depleted than equal-mass slow rotators. We intend to explore the existence of a Li depletion - rotation connection among the beta Pictoris members at an age of about 24 Myr, and to use such correlation either to confirm or to improve the age estimate based on the Lithium Depletion Boundary (LDB) modeling. We have photometrically monitored all the known members of the beta Pictoris association with at least one Lithium equivalent width (Li EW) measurement from the literature. We measured the rotation periods of 30 members for the first time and retrieved from the literature the rotation periods for other 36 members, building a catalogue of 66 members with measured rotation period and Li EW. We find that in the 0.3 < M < 0.8 Msun range, there is a strong correlation between rotation and Li EW. For higher mass stars, no significant correlation is found. For very low mass stars in the Li depletion onset, at about 0.1 Msun, data are too few to infer a significant correlation. The observed Li EWs are compared with those predicted by the Dartmouth stellar evolutionary models that incorporate the effects of magnetic fields. After decorrelating the Li EW from the rotation period, we find that the hot side of the LDB is fitted well by Li EW values corresponding to an age of 25$\pm$3 Myr in good agreement with independent estimates from the literature.

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

Aims. In this work, we characterize the low mass companion of the A0 field star HR3549. Methods. We observed HR3549AB in imaging mode with the the NIR branch (IFS and IRDIS) of SPHERE@VLT, with IFS in YJ mode and IRDIS in the H band. We also acquired a medium resolution spectrum with the IRDIS long slit spectroscopy mode. The data were reduced using the dedicated SPHERE GTO pipeline, purposely designed for this instrument. We employed algorithms such as PCA and TLOCI to reduce the speckle noise. Results. The companion was clearly visible both with IRDIS and IFS.We obtained photometry in four different bands as well as the astrometric position for the companion. Based on our astrometry, we confirm that it is a bound object and put constraints on its orbit. Although several uncertainties are still present, we estimate an age of ~100-150 Myr for this system, yielding a most probable mass for the companion of 40-50MJup and T_eff ~300-2400 K. Comparing with template spectra points to a spectral type between M9 and L0 for the companion, commensurate with its position on the color-magnitude diagram.

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

The space telescope CoRoT has provided light curves of T Tauri stars belonging to the star-forming region of NGC 2264 with unprecedented continuity and precision in the framework of a coordinated multi-wavelength observational project. We perform spot modelling of the optical light curves of five weak-line T Tauri stars whose variability is dominated by starspots. We apply two-spot and evolving single spot models in the framework of a Bayesian Monte Carlo Markov Chain approach to derive the a posteriori distribution of the starspot parameters and the inclination of the star rotation axis. We focus on the rotation periods of the spots that can provide evidence for differential rotation in those stars. We find meaningful results in the case of three stars with an inclination higher than 50 deg and a slow variation of the light curve amplitude. The relative difference of the spot rotation periods ranges from 0.02 to 0.05 that is 3-10 times larger than the amplitude of the differential rotation found in similar stars with Doppler imaging techniques. We conclude that the intrinsic starspot evolution, although very slow, has a significant impact on the determination of the differential rotation by means of our spot modelling approach. We estimate typical timescales for the evolution of the starspot pattern between about 20 and 50 rotation periods in our stars.

We investigate here the lithium content of young low-mass stars in the 5~Myr-old star forming region NGC~2264 and its relationship with rotation. We combine lithium equivalent width measurements, EW(Li), from the Gaia-ESO Survey with the determination of rotational periods from the CSI~2264 survey. We consider only bona fide non accreting cluster members in order to minimize uncertainties on EW(Li). We report the existence of a relationship between lithium content and rotation in NGC~2264 at an age of 5~Myr. The Li-rotation connection is seen over a restricted temperature range, Teff=3800--4400K, where fast rotators are Li-rich compared to slow ones. This correlation is similar to, albeit of lower amplitude than, the Li-rotation connection previously reported for K dwarfs in the 125 Myr-old Pleiades cluster. We investigate whether the non-standard pre-main sequence models developed so far to explain the Pleiades results, which are based on episodic accretion, pre-main sequence core-envelope decoupling, and/or radius inflation due to enhanced magnetic activity, can account for an early development of the Li-rotation connection. While radius inflation appears to be the most promising possibility, each of these models has issues. We therefore also discuss external causes that might operate during the first few Myr of pre-main sequence evolution, such as planet engulfment and/or steady disk accretion, as possible candidates for the common origin for Li-excess and fast rotation in young low-mass pre-main sequence stars. The emergence of a connection between lithium content and rotation rate at such an early age as 5~Myr suggests a complex link between accretion processes, early angular momentum evolution, and possibly planet formation, which likely impacts early stellar evolution and still is to be fully deciphered.

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec \theta)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec \theta)_\mathrm{max}$.