Search

We investigate the ability of the Euclid telescope to detect galaxy-scale gravitational lenses. To do so, we perform a systematic visual inspection of the $0.7\,\rm{deg}^2$ Euclid ERO data towards the Perseus cluster using both the high-resolution VIS $I_{\scriptscriptstyle\rm E}$ band, and the lower resolution NISP bands. We inspect every extended source brighter than magnitude $23$ in $I_{\scriptscriptstyle\rm E}$ with $41$ expert human classifiers. This amounts to $12\,086$ stamps of $10^{\prime\prime}\,\times\,10^{\prime\prime}$. We find $3$ grade A and $13$ grade B candidates. We assess the validity of these $16$ candidates by modelling them and checking that they are consistent with a single source lensed by a plausible mass distribution. Five of the candidates pass this check, five others are rejected by the modelling and six are inconclusive. Extrapolating from the five successfully modelled candidates, we infer that the full $14\,000\,{\rm deg}^2$ of the Euclid Wide Survey should contain $100\,000^{+70\,000}_{-30\,000}$ galaxy-galaxy lenses that are both discoverable through visual inspection and have valid lens models. This is consistent with theoretical forecasts of $170\,000$ discoverable galaxy-galaxy lenses in Euclid. Our five modelled lenses have Einstein radii in the range $0.\!\!^{\prime\prime}68\,<\,\theta_\mathrm{E}\,<1.\!\!^{\prime\prime}24$, but their Einstein radius distribution is on the higher side when compared to theoretical forecasts. This suggests that our methodology is likely missing small Einstein radius systems. Whilst it is implausible to visually inspect the full Euclid data set, our results corroborate the promise that Euclid will ultimately deliver a sample of around $10^5$ galaxy-scale lenses., Comment: 21 pages, 20 figures, submitted to A&A

Octans is one of the most distant ($d\sim150$pc) young stellar associations of the solar neighbourhood. Its age is still poorly constrained in the literature and requires further investigation. We take advantage of the state-of-the-art astrometry delivered by the third data release of the Gaia space mission combined with radial velocity measurements obtained from high-resolution spectroscopy to compute the 3D positions and 3D spatial velocities of the stars and derive the dynamical traceback age of the association. We performed an extensive traceback analysis using different subsamples of stars, different metrics to define the size of the association, and different models for the Galactic potential to integrate the stellar orbits in the past. We derive a dynamical age of $34^{+2}_{-2}$Myr that is independent from stellar models and represents the most precise age estimate currently available for the Octans association. After correcting the radial velocity of the stars for the effect of gravitational redshift, we obtain a dynamical age of $33^{+3}_{-1}$Myr, which is in very good agreement with our first solution. This shows that the effect of gravitational redshift is small for such a distant young stellar association. Our result is also consistent with the less accurate age estimates obtained in previous studies from lithium depletion (30-40Myr) and isochrones (20-30Myr). By integrating the stellar orbits in time, we show that the members of Octans and Octans-Near had different locations in the past, which indicates that the two associations are unrelated despite the close proximity in the sky. Our results confirm that it is possible to derive precise dynamical ages via the traceback method for $\sim30$Myr old stellar clusters at about $\sim150$pc with the same precision level that has been achieved in other studies for young stellar groups within 50pc of the Sun., Comment: 12 pages, 8 figures, A&A (accepted)

The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exceptional depth and spatial resolution of this wide-field multi-band data enable the simultaneous detection and characterisation of both bright and low surface brightness galaxies, along with their globular cluster systems, from the optical to the NIR. This study advances beyond previous analyses of the cluster and enables a range of scientific investigations summarised here. We derive the luminosity and stellar mass functions (LF and SMF) of the Perseus cluster in the Euclid IE band, thanks to supplementary u,g,r,i,z and Halpha data from the CFHT. We adopt a catalogue of 1100 dwarf galaxies, detailed in the corresponding ERO paper. We identify all other sources in the Euclid images and obtain accurate photometric measurements using AutoProf or AstroPhot for 138 bright cluster galaxies, and SourceExtractor for half a million compact sources. Cluster membership for the bright sample is determined by calculating photometric redshifts with Phosphoros. Our LF and SMF are the deepest recorded for the Perseus cluster, highlighting the groundbreaking capabilities of the Euclid telescope. Both the LF and SMF fit a Schechter plus Gaussian model. The LF features a dip at M(IE)=-19 and a faint-end slope of alpha_S = -1.2 to -1.3. The SMF displays a low-mass-end slope of alpha_S = -1.2 to -1.35. These observed slopes are flatter than those predicted for dark matter halos in cosmological simulations, offering significant insights for models of galaxy formation and evolution., Comment: Submitted to A&A, 44 pages, 35 figures, Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

We provide an early assessment of the imaging capabilities of the Euclid space mission to probe deeply into nearby star-forming regions and associated very young open clusters, and in particular to check to what extent it can shed light on the new-born free-floating planet population. This paper focuses on a low-reddening region observed in just one Euclid pointing where the dust and gas has been cleared out by the hot sigma Orionis star. One late-M and six known spectroscopically confirmed L-type substellar members in the sigma Orionis cluster are used as benchmarks to provide a high-purity procedure to select new candidate members with Euclid. The exquisite angular resolution and depth delivered by the Euclid instruments allow us to focus on bona-fide point sources. A cleaned sample of sigma Orionis cluster substellar members has been produced and the initial mass function (IMF) has been estimated by combining Euclid and Gaia data. Our sigma Orionis substellar IMF is consistent with a power-law distribution with no significant steepening at the planetary-mass end. No evidence of a low-mass cutoff is found down to about 4 Jupiter masses at the young age (3 Myr) of the sigma Orionis open cluster., Comment: 13 pages, 10 figures, submitted to the A&A special issue on "Euclid on the sky"

The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline to create visually compelling images while simultaneously meeting the scientific demands within months of launch, leveraging a pragmatic, data-driven development strategy. The pipeline's key requirements are to preserve the image quality and to provide flux calibration and photometry for compact and extended sources. The pipeline's five pillars are: removal of instrumental signatures; astrometric calibration; photometric calibration; image stacking; and the production of science-ready catalogues for both the VIS and NISP instruments. We report a PSF with a full width at half maximum of 0.16" in the optical and 0.49" in the three NIR bands. Our VIS mean absolute flux calibration is accurate to about 1%, and 10% for NISP due to a limited calibration set; both instruments have considerable colour terms. The median depth is 25.3 and 23.2 AB mag with a SNR of 10 for galaxies, and 27.1 and 24.5 AB mag at an SNR of 5 for point sources for VIS and NISP, respectively. Euclid's ability to observe diffuse emission is exceptional due to its extended PSF nearly matching a pure diffraction halo, the best ever achieved by a wide-field, high-resolution imaging telescope. Euclid offers unparalleled capabilities for exploring the LSB Universe across all scales, also opening a new observational window in the NIR. Median surface-brightness levels of 29.9 and 28.3 AB mag per square arcsec are achieved for VIS and NISP, respectively, for detecting a 10 arcsec x 10 arcsec extended feature at the 1 sigma level., Comment: Submitted to A&A, 44 pages, 36 figures - Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance., Comment: Accepted for publication in the A&A special issue`Euclid on Sky'

Lithium (Li) abundance is an age indicator for G, K, and M stellar types, as its abundance decreases over time for these spectral types. However, despite the observational efforts made over the past few decades, the role of rotation, activity, and metallicity in the depletion of Li is still unclear. We have investigated how Li depletion is affected by rotation and metallicity in G and K members of the roughly Pleiades-aged open cluster M35. To do so, we have collected a sample of 165 candidate members observed with the WIYN/Hydra spectrograph. In addition, we have taken advantage of three previous spectroscopic studies of Li in M35. As a result, we have collected a final sample of 396 stars which we have classified as members and non-members of the cluster. We have measured iron abundances, Li equivalent widths, and Li abundances for the 110 M35 members added to the existing sample by this study. Finally, rotation periods for cluster members have been obtained from the literature or derived from Zwicky Transient Facility light curves. As a result, we have confirmed that fast G and K rotators are Li-rich in comparison with slow rotators of similar effective temperature. Furthermore, while we derived subsolar metallicity for M35 from our spectra, the distribution of Li in this cluster is similar to those observed for the Pleiades and M34, which have solar metallicity and slightly different ages. In addition, we have shown that an empirical relationship proposed to remove the contribution of the Fe I line at 670.75 nm to the blended feature at 670.78 nm overestimates the contribution of this iron line for M35 members. We conclude that a 0.2-0.3 dex difference in metallicity makes little difference in the Li distributions of open clusters with ages between 100 and 250 Myr., Comment: 14 pages, 11 figures, accepted for publication in A\&A

We aim at characterizing the low-mass (sub)stellar population of the central portion (2.4 pc$^2$) of the $\sim$2 Myr old cluster NGC 2244 using near infrared spectroscopy. By studying this cluster, characterized by a low stellar density and numerous OB stars, we aim at exploring the effect that OB stars may have on the production of BDs. We obtain near infrared HK spectroscopy of 85 faint candidate members of NGC 2244. We derive the spectral type and extinction by comparison with spectral templates. We evaluate cluster membership using three gravity-sensitive spectral indices based on the shape of the $H$-band. Furthermore, we evaluate the infrared excess from Spitzer of all the candidate members of the cluster. Finally, we estimate the mass of all the candidate members of the cluster and derive the initial mass function, star-to-BD number ratio and disk fraction. The initial mass function is well represented by a power law ($dN/dM \propto M^{-\alpha}$) below 0.4 $M_\odot$, with a slope $\alpha$ = 0.7-1.1 depending on the fitted mass range. We calculate a star-to-BD number ratio of 2.2-2.8. We find the low-mass population of NGC 2244 to be consistent with nearby star-forming regions, although it is at the high-end of BD production. We find BDs in NGC 2244 to be on average closer to OB stars than to low-mass stars, which could potentially be the first evidence of OB stars affecting the formation of BDs. We find a disk fraction of all the members with spectral type later than K0 of 39$\pm$9% which is lower than typical values found in nearby star-forming regions of similar ages., Comment: Accepted by A&A, KMOS spectroscopic data will be made public on Vizier upon publication

Context. Star-forming regions are excellent benchmarks for testing and validating theories of star formation and stellar evolution. The Perseus star-forming region being one of the youngest (<10 Myr), closest (280-320 pc), and most studied in the literature, is a fundamental benchmark. Aims. We aim to study the membership, phase-space structure, mass, and energy (kinetic plus potential) distribution of the Perseus star-forming region using public catalogues (Gaia, APOGEE, 2MASS, PanSTARRS). Methods. We use Bayesian methodologies accounting for extinction to identify the Perseus physical groups in the phase-space, retrieve their candidate members, derive their properties (age, mass, 3D positions, 3D velocities, and energy), and attempt to reconstruct their origin. Results. We identify 1052 candidate members in seven physical groups (one of them new) with ages between 3 and 10 Myr, dynamical super-virial states, and large fractions of energetically unbound stars. Their mass distributions are broadly compatible with that of Chabrier for masses >0.1 $M_\odot$ and do not show hints of over-abundance of low-mass stars in NGC1333 with respect to IC348. These groups' ages, spatial structure, and kinematics are compatible with at least three generations of stars. Future work is still needed to clarify if the formation of the youngest was triggered by the oldest. Conclusions. The exquisite Gaia data complemented with public archives and mined with comprehensive Bayesian methodologies allow us to identify 31% more members than in previous studies, discover a new physical group (Gorgophone: 7 Myr, 191 members, and 145 $M_\odot$), and confirm that the spatial, kinematic, and energy distributions of these groups support the hierarchical star-formation scenario., Comment: Accepted for publication in Astronomy & Astrophysics, 29 pages, 22 Figures, 8 tables one of them electronic

(Pre-)transitional disks show gaps and cavities that can be related with on-going planet formation. According to theory, young embedded planets can accrete material from the circumplanetary and circumstellar disks, so that they could be detected in accretion tracers, like the H$_{\alpha}$ emission line. In this work, we present spectral angular differential imaging AO-assisted observations of five (pre-)transitional disks obtained with SPHERE/ZIMPOL at the Very Large Telescope (VLT). They were obtained in the H$_{\alpha}$ line and the adjacent continuum. We have combined spectral and angular differential imaging techniques to increase the contrast in the innermost regions close to the star, and search for the signature of young accreting protoplanets. As a result, the reduced images do not show any clear H$_{\alpha}$ point source around any of the targets. We report faint H$_{\alpha}$ emissions around TW Hya and HD163296: while the former is most probably an artifact related with a spike, the nature of the latter remains unclear. The spectral and angular differential images yield contrasts of 6--8 magnitudes at separations of $\sim$ 100 mas from the central stars, except in the case of LkCa15, with values of $\sim$3 mag. We have estimated upper limits to the accretion luminosity of potential protoplanets, obtaining that planetary models provide an average value of $L_{\rm acc} \sim 10^{-4}$ $L_{\odot}$ at 200 mas, which is $\sim$2 orders of magnitude higher than the $L_{\rm acc}$ estimated from the extrapolation of the $L_{H_{\alpha}}$ - $L_{acc}$ stellar relationship. We explain the lack of protoplanet detections as a combination of different factors, like e.g. episodic accretion, extinction from the circumstellar and circumplanetray disks, and/or a majority of low-mass, low-accreting planets., Comment: 17 pages, Accepted for publication in A&A

We aim to study spatial and kinematic substructures of the region encompassed by Upper Scorpius and Ophiuchus star-forming regions, determining dynamical traceback ages and studying the star formation history of the complex. We identified seven different groups in this region. Four groups (nu Sco, beta Sco, sigma Sco and delta Sco) are part of Upper Scorpius, two groups (rho Oph and alpha Sco) are in Ophiuchus, and another group (pi Sco) is a nearby young population. We found an age gradient from the rho Oph group (the youngest) to the delta Sco group (<~5 Myr), showing that star formation was a sequential process for the past 5 Myr. Our traceback analysis shows that Upper Scorpius and rho Oph groups share a common origin. The closer group of pi Sco is probably older, and the traceback analysis suggests that this group and the alpha Sco group have a different origin, likely related to other associations in the Sco-Cen complex. Our study shows that this region has a complex star formation history that goes beyond the current formation scenario, likely a result of stellar feedback from massive stars, supernovae explosions, and dynamic interactions between stellar groups and the molecular gas. In particular, we speculate that photo-ionisation from the massive delta Sco star could have triggered star formation first in the beta Sco group and then in the nu Sco group. The perturbations of stellar orbits due to stellar feedback and dynamical interactions could also be responsible for the 1-3 Myr difference that we found between dynamical traceback ages and isochronal ages., Comment: Accepted for publication in A&A

Context: A rich population of low-mass brown dwarfs and isolated planetary mass objects has been reported recently in the Upper Scorpius and Ophiuchus star forming complex. Aims: We investigate the membership, nature and properties of 17 of these isolated planetary mass candidates using low-resolution near-infrared spectra. Methods: We investigate the membership by looking for evidences of youth using four diagnostics: the slope of the continuum between the J and Ks band, the Hcont and TLI-g gravity sensitive indices, and by comparing the spectra to young and field (old) M and L-dwarf standards. Results: All the targets but one are confirmed as young ultracool objects, with spectral types between L0 and L6 and masses in the range 0.004-0.013 M according to evolutionary models. The status of the last target is unclear at this point. Conclusions: Only one possible contaminant has been identified among the 17 targets, suggesting that the contamination level of the original sample must be lower than 6%, Comment: Accepted on 2022 June 2. Pending language revision

Stellar ages are key to improving our understanding of different astrophysical phenomena. However, many techniques to estimate stellar ages are highly model-dependent. The lithium depletion boundary (LDB), based on the presence or absence of lithium in low-mass stars, can be used to derive ages in stellar associations of between 20 and 500~Ma. The purpose of this work is to revise former LDB ages in stellar associations in a consistent way, taking advantage of the homogeneous \textit{Gaia} parallaxes as well as bolometric luminosity estimations that do not rely on monochromatic bolometric corrections. We studied nine open clusters and three moving groups characterised by a previous determination of the LDB age. We gathered all the available information from our data and the literature: membership, distances, photometric data, reddening, metallicity, and surface gravity. We re-assigned membership and calculated bolometric luminosities and effective temperatures using distances derived from Gaia DR2 and multi-wavelength photometry for individual objects around the former LDB. We located the LDB using a homogeneous method for all the stellar associations. Finally, we estimated the age by comparing it with different evolutionary models. We located the LDB for the twelve stellar associations and derived their ages using several theoretical evolutionary models. We compared the LDB ages among them, along with data obtained with other techniques, such as isochrone fitting, ultimately finding some discrepancies among the various approaches. Finally, we remark that the 32 Ori MG is likely to be composed of at least two populations of different ages., Comment: 78 pages, accepted by A&A

Context: The $\chi^{1}$ Fornacis cluster (Alessi 13) is one of a few open clusters of its age and distance in the Solar neighbourhood that ought to benefit from more attention as it can serve as a cornerstone for numerous future studies related to star and planet formation. Aims: We take advantage of the early installment of the third data release of the Gaia space mission in combination with archival data and our own observations, to expand the census of cluster members and revisit some properties of the cluster. Methods: We applied a probabilistic method to infer membership probabilities over a field of more than 1000 deg${^2}$ to select the most likely cluster members and derive the distances, spatial velocities, and physical properties of the stars in this sample. Results: We identify 164 high-probability cluster members (including 61 new members) covering the magnitude range from 5.1 to 19.6 mag in the G-band. Our sample of cluster members is complete down to 0.04 M$_{\odot}$. We derive the distance of $108.4\pm0.3$ pc from Bayesian inference and confirm that the cluster is comoving with the Tucana-Horologium, Columba, and Carina young stellar associations. We investigate the kinematics of the cluster from a subsample of stars with measured radial velocities and we do not detect any significant expansion or rotation effects in the cluster. Our results suggest that the cluster is somewhat younger (about 30 Myr) than previously thought. Based on spectroscopic observations, we argue that the cluster is mass-segregated and that the distribution of spectral types shows little variation compared to other young stellar groups. Conclusions: In this study, we deliver the most complete census of cluster members that can be done with Gaia data alone and we use this new sample to provide an updated picture on the 6D structure of the cluster., Comment: 16 pages, 12 Figures, online material, A&A (in press)

Exoplanet searches through space-based photometric time series have shown to be very efficient in recent years. However, follow-up efforts on the detected planet candidates have been demonstrated to be critical to uncover the true nature of the transiting objects. In this paper we show a detailed analysis of one of those false positives hidden as planetary signals. In this case, the candidate KOI-3886.01 showed clear evidence of a planetary nature from various techniques. Indeed, the properties of the fake planet set it among the most interesting and promising for the study of planetary evolution as the star leaves the main sequence. To unveil the true nature of this system, we present a complete set of observational techniques including high-spatial resolution imaging, high-precision photometric time series (showing eclipses, phase curve variations and asteroseismology signals), high-resolution spectroscopy and derived radial velocities, to unveil the true nature of this planet candidate. We find that KOI-3886.01 is an interesting false positive case: a hierarchical triple system composed by a $\sim$K2III giant star (KOI-3886A) accompanied by a close-in eclipsing binary formed by a subgiant $\sim$G4IV star (KOI-3886B) and a brown dwarf (KOI-3886C). In particular, KOI-3886C is one of the most irradiated brown dwarfs known to date, showing the largest radius in this substellar regime. It is also the first eclipsing brown dwarf known around an evolved star. In this paper we highlight the relevance of complete sets of follow-up observations to extrasolar planets detected by the transit technique using large-pixel photometers such as Kepler and TESS, and in the future, PLATO. In particular, multi-color high-spatial resolution imaging was the first hint toward ruling out the planet scenario in this system., Comment: Accepted for publication in A&A, 21 pages, 18 figures, 5 tables (language corrected version)

Context. The analysis of luminosity and mass distributions of young stellar clusters is essential to understanding the star-formation process. However, the gas and dust left over by this process extinct the light of the newborn stars and can severely bias both the census of cluster members and its luminosity distribution. Aims. We aim to develop a Bayesian methodology to infer, with minimal biases due to photometric extinction, the candidate members and magnitude distributions of embedded young stellar clusters. Methods. We improve a previously published methodology and extend its application to embedded stellar clusters. We validate the method using synthetically extincted data sets of the Pleiades cluster with varying degrees of extinction. Results. Our methodology can recover members from data sets extincted up to Av ~ 6 mag with accuracies, true positive, and contamination rates that are better than 99%, 80%, and 9%, respectively. Missing values hamper our methodology by introducing contaminants and artifacts into the magnitude distributions. Nonetheless, these artifacts vanish through the use of informative priors in the distribution of the proper motions. Conclusions. The methodology presented here recovers, with minimal biases, the members and distributions of embedded stellar clusters from data sets with a high percentage of sources with missing values (>96%)., Comment: 21 pages, 24 figures

Context: Chamaeleon is the southernmost low-mass star-forming complex within 200 pc from the Sun. Its stellar population has been extensively studied in the past, but the current census of the stellar content is not complete yet and deserves further investigation. Aims: We take advantage of the second data release of the \textit{Gaia} space mission to expand the census of stars in Chamaeleon and to revisit the properties of the stellar populations associated to the Chamaeleon I (Cha I) and Chamaeleon II (Cha II) dark clouds. Methods: We perform a membership analysis of the sources in the \textit{Gaia} catalogue over a field of 100 deg$^{2}$ encompassing the Chamaeleon clouds, and use this new census of cluster members to investigate the 6D structure of the complex. Results: We identify 188 and 41 high-probability members of the stellar populations in Cha I and Cha II, respectively, including 19 and 7 new members. Our sample covers the magnitude range from $G=6$ to $G=20$ mag in Cha I, and from $G=12$ to $G=18$ mag in Cha II. We confirm that the northern and southern subgroups of Cha I are located at different distances ($191.4^{+0.8}_{-0.8}$ pc and $186.7^{+1.0}_{-1.0}$ pc), but they exhibit the same space motion within the reported uncertainties. Cha II is located at a distance of $197.5^{+1.0}_{-0.9}$ pc and exhibits a space motion that is consistent with Cha I within the admittedly large uncertainties on the spatial velocities of the stars that come from radial velocity data. The median age of the stars derived from the Hertzsprung-Russell diagram (HRD) and stellar models is about 1-2 Myr, suggesting that they are somewhat younger than previously thought. We do not detect significant age differences between the Chamaeleon subgroups, but we show that Cha II exhibits a higher fraction of disc-bearing stars compared to Cha I., Comment: 22 pages, 13 figures, online material (Tables A.1 - A.7). Accepted for publication in A&A

Ultracool dwarfs encompass a wide variety of compact stellar-like objects with spectra classified as late-M, L, T and Y. Most of them have been discovered using wide-field imaging surveys. The Virtual Observatory has proven to be of great utility to efficiently exploit these astronomical resources. We aim to validate a Virtual Observatory methodology designed to discover and characterize ultracool dwarfs in deep extragalactic surveys like ALHAMBRA and COSMOS. Three complementary searches based on parallaxes, proper motions and colours, respectively were carried out. A total of 897 candidate ultracool dwarfs were found, with only 16 previously reported in SIMBAD. Most of the new UCDs reported here are likely late-M and L dwarfs because of the limitations imposed by the utilization of optical ($Gaia$ DR2 and r-band) data. We complement ALHAMBRA and COSMOS photometry with other catalogues in the optical and infrared using VOSA, a Virtual Observatory tool that estimates effective temperatures from the spectral energy distribution fitting to collections of theoretical models. The agreement between the number of UCDs found in the COSMOS field and theoretical estimations together with the low false negative rate (known UCDs not discovered in our search) validates the methodology proposed in this work, which will be used in the forthcoming wide and deep surveys provided by the Euclid space mission. Simulations of Euclid number counts for UCDs detectable in different photometric passbands are presented for a wide survey area of 15,000 square degrees, and the limitations of applicability of Euclid data to detect UCDs using the methods employed in this paper are discussed.

Context: Lupus is recognised as one of the closest star-forming regions, but the lack of trigonometric parallaxes in the pre-Gaia era hampered many studies on the kinematic properties of this region and led to incomplete censuses of its stellar population. Aims: We use the second data release of the Gaia space mission combined with published ancillary radial velocity data to revise the census of stars and investigate the 6D structure of the Lupus complex. Methods: We performed a new membership analysis of the Lupus association based on astrometric and photometric data over a field of 160 deg2 around the main molecular clouds of the complex and compared the properties of the various subgroups in this region. Results: We identified 137 high-probability members of the Lupus association of young stars, including 47 stars that had never been reported as members before. Many of the historically known stars associated with the Lupus region identified in previous studies are more likely to be field stars or members of the adjacent Scorpius-Centaurus association. Our new sample of members covers the magnitude and mass range from G=8 to G=18 mag and from 0.03 to 2.4Msun, respectively. We compared the kinematic properties of the stars projected towards the molecular clouds Lupus 1 to 6 and showed that these subgroups are located at roughly the same distance (about 160~pc) and move with the same spatial velocity. Our age estimates inferred from stellar models show that the Lupus subgroups are coeval (with median ages ranging from about 1 to 3 Myr). The Lupus association appears to be younger than the population of young stars in the Corona-Australis star-forming region recently investigated by our team using a similar methodology. The initial mass function of the Lupus association inferred from the distribution of spectral types shows little variation compared to other star-forming regions., Comment: 20 pages, 14 figures, A&A in press

Context: Stellar clusters are benchmarks for theories of star formation and evolution. The high precision parallax data of the Gaia mission allows significant improvements in the distance determination to stellar clusters and its stars. In order to have accurate and precise distance determinations, systematics like the parallax spatial correlations need to be accounted for, especially for stars in small sky regions. Aims: Provide the astrophysical community with a free and open code designed to simultaneously infer cluster parameters (i.e. distance and size) and the distances to its stars using Gaia parallax measurements. It includes cluster oriented prior families and is specifically designed to deal with the Gaia parallax spatial correlations. Methods: A Bayesian hierarchical model is created to allow the inference of both the cluster parameters and distances to its stars. Results: Using synthetic data that mimics Gaia parallax uncertainties and spatial correlations, we observe that our cluster oriented prior families result in distance estimates with smaller errors than those obtained with an exponentially decreasing space density prior. In addition, the treatment of the parallax spatial correlations minimizes errors in the estimated cluster size and stellar distances and avoids the underestimation of uncertainties. Although neglecting the parallax spatial correlations has no impact on the accuracy of cluster distance determinations, it underestimates the uncertainties and may result in measurements that are incompatible with the true value. Conclusions: The combination of prior knowledge with the treatment of Gaia parallax spatial correlations produces accurate (error <10%) and trustworthy estimates (i.e. true values contained within the 2$\sigma$ uncertainties) of clusters distances for clusters up to ~5 kpc, and cluster sizes for clusters up to ~1 kpc., Comment: 14 pages, 10 figures, accepted in A&A, for the associated code visit https://github.com/olivares-j/Kalkayotl

Context: The $\beta$ Pictoris moving group is one of the most well-known young associations in the solar neighbourhood and several members are known to host circumstellar discs, planets, and comets. Measuring its age with precision is basic to study several astrophysical processes such as planet formation and disc evolution which are strongly age dependent. Aims: We aim to determine a precise and accurate dynamical traceback age for the $\beta$ Pictoris moving group. Methods: Our sample combines the extremely precise Gaia DR2 astrometry with ground-based radial velocities measured in an homogeneous manner. We use an updated version of our algorithm to determine dynamical ages. The new approach takes into account a robust estimate of the spatial and kinematic covariance matrices of the association to improve the sample selection process and to perform the traceback analysis. Results: We estimate a dynamical age of $18.5_{-2.4}^{+2.0}$ Myr for the $\beta$ Pictoris moving group. We investigated the spatial substructure of the association at birth time and we propose the existence of a core of stars more concentrated. We also provide precise radial velocity measurements for 81 members of $\beta$ Pic, including ten stars with the first determination of their radial velocities. Conclusions: Our dynamical traceback age is three times more precise than previous traceback age estimates and, more important, for the first time, reconciles the traceback age with the most recent estimates of other dynamical, lithium depletion boundary, and isochronal ages. This has been possible thanks to the excellent astrometric and spectroscopic precisions, the homogeneity of our sample, and the detailed analysis of binaries and membership.

Near-infrared spectroscopic surveys have uncovered a population of short-period, blended-light spectral binaries composed of low-mass stars and brown dwarfs. These systems are amenable to orbit determination and individual mass measurements via astrometric monitoring. Here we present first results of a multi-year campaign to obtain high-precision absolute astrometry for spectral binaries using the Gemini-South and Gemini-North GMOS imagers. We measure the complete astrometric orbits for two systems: 2M0805$+$48 and 2M1059$-$21. Our astrometric orbit of 2M0805$+$48 is consistent with its 2-year radial velocity orbit determined previously and we find a mass of $66^{+5}_{-14} M_\mathrm{Jup}$ for its T5.5 companion. For 2M1059$-$21 we find a 1.9 year orbital period and a mass of $67^{+4}_{-5} M_\mathrm{Jup}$ for its T3.5 companion. We demonstrate that sub-milliarcsecond absolute astrometry can be obtained with both GMOS imagers and that this is an efficient avenue for confirming and characterising ultracool binary systems., Comment: 14 pages, 13 figures, accepted for publication in MNRAS

Context. As a fragile element, lithium is a sensitive probe of physical processes occurring in stellar interiors. Aims. We aim at investigating the relationship between lithium abundance and rotation rate in low-mass members of the newly discovered 125~Myr-old Psc-Eri stellar stream. Methods. We obtained high resolution optical spectra and measure the equivalent width of the 607.8 nm LiI line for 40 members of the Psc-Eri stream, whose rotational periods have been derived by arXiv:1905.10588. Results. We show that a tight correlation exists between lithium content and rotation rate among the late-G to early K-type stars of the Psc-Eri stream. Fast rotators are systematically Li-rich, while slow rotators are Li-depleted. This trend mimics the one previously reported for the similar age Pleiades cluster. Conclusions. The lithium-rotation connection thus seems to be universal over a restricted effective temperature range for low-mass stars at or close to the zero-age main sequence, and does not depend on environmental conditions., Comment: 5 pages, 2 figures, 1 table. Accepted to A&A letters

Corona-Australis is one of the nearest regions to the Sun with recent and ongoing star formation, but the current picture of its stellar (and substellar) content is not complete yet. We take advantage of the second data release of the Gaia space mission to revisit the stellar census and search for additional members of the young stellar association in Corona-Australis. We applied a probabilistic method to infer membership probabilities based on a multidimensional astrometric and photometric data set over a field of 128 deg$^{2}$ around the dark clouds of the region. We identify 313 high-probability candidate members to the Corona-Australis association, 262 of which had never been reported as members before. Our sample of members covers the magnitude range between $G\gtrsim5$ mag and $G\lesssim20$ mag, and it reveals the existence of two kinematically and spatially distinct subgroups. There is a distributed `off-cloud' population of stars located in the north of the dark clouds that is twice as numerous as the historically known `on-cloud' population that is concentrated around the densest cores. By comparing the location of the stars in the HR-diagram with evolutionary models, we show that these two populations are younger than 10 Myr. Based on their infrared excess emission, we identify 28 Class II and 215 Class III stars among the sources with available infrared photometry, and we conclude that the frequency of Class II stars (i.e. `disc-bearing' stars) in the on-cloud region is twice as large as compared to the off-cloud population. The distance derived for the Corona-Australis region based on this updated census is $d=149.4^{+0.4}_{-0.4}$ pc, which exceeds previous estimates by about 20 pc.In this paper we provide the most complete census of stars in Corona-Australis available to date that can be confirmed with Gaia data., Comment: 17 pages, 11 figures, online data. Abridged abstract. Accepted for publication in A&A (in press)

Aims:We take advantage of the second data release of the Gaia space mission and the state-of-the-art astrometry delivered from very long baseline interferometry observations to revisit the structure and kinematics of the nearby Taurus star-forming region. Methods: We apply a hierarchical clustering algorithm for partitioning the stars in our sample into groups (i.e., clusters) that are associated with the various molecular clouds of the complex, and derive the distance and spatial velocity of individual stars and their corresponding molecular clouds. Results: We show that the molecular clouds are located at different distances and confirm the existence of important depth effects in this region reported in previous studies. For example, we find that the L 1495 molecular cloud is located at $d=129.9^{+0.4}_{-0.3}$ pc, while the filamentary structure connected to it (in the plane of the sky) is at $d=160.0^{+1.2}_{-1.2}$ pc. We report B 215 and L 1558 as the closest ($d=128.5^{+1.6}_{-1.6}$ pc) and most remote ($d=198.1^{+2.5}_{-2.5}$ pc) substructures of the complex, respectively. The median inter-cloud distance is 25 pc and the relative motion of the subgroups is on the order of a few km/s. We find no clear evidence for expansion (or contraction) of the Taurus complex, but signs of the potential effects of a global rotation. Finally, we compare the radial velocity of the stars with the velocity of the underlying $^{13}$CO molecular gas and report a mean difference of $0.04\pm0.12$ km/s (with r.m.s. of 0.63 km/s) confirming that the stars and the gas are tightly coupled., Comment: Accepted for publication in A&A

Context. The study of star formation is extremely challenging due to the lack of complete and clean samples of young, nearby clusters, and star forming regions. The recent Gaia DR2 catalogue complemented with the deep, ground based COSMIC DANCe catalogue offers a new database of unprecedented accuracy to revisit the membership of clusters and star forming regions. The 30 Myr open cluster IC 4665 is one of the few well-known clusters of this age and it is an excellent target where to test evolutionary models and study planetary formation. Aims. We aim to provide a comprehensive membership analysis of IC 4665 and to study the following properties: empirical isochrones, distance, magnitude distribution, present-day system mass function, and spatial distribution. Methods. We use the Gaia DR2 catalogue together with the DANCe catalogue to look for members using a probabilistic model of the distribution of the observable quantities in both the cluster and background populations. Results. We obtain a final list of 819 candidate members which cover a 12.4 magnitude range (7 < J < 19.4). We find that 50% are new candidates, and we estimate a conservative contamination rate of 20%. This unique sample of members allows us to obtain a present-day system mass function in the range of 0.02-6 Msun, which reveals a number of details not seen in previous studies. In addition, they favour a spherically symmetric spatial distribution for this young open cluster. Conclusions. Our membership analysis represents a significant increase in the quantity and quality (low-contamination) with respect to previous studies. As such, it offers an excellent opportunity to revisit other fundamental parameters such as the age.

Context. Ruprecht 147 is the oldest (2.5 Gyr) open cluster in the solar vicinity (< 300 pc), making it an important target for stellar evolution studies and exoplanet searches. Aims. Derive a census of members and the luminosity, mass, and spatial distributions of the cluster. Methods. We use an astro-photometric data set including all available information from the literature together with our own observations. We process the data with an updated version of an existent membership selection methodology. Results. We identify 259 high-probability candidate members, including 58 previously unreported. All these candidates cover the luminosity interval between G > 6 mag to i< 21 mag. The cluster luminosity and mass distributions are derived with an unprecedented level of details allowing us to recognize, among other features, the Wielen dip. The mass distribution in the low-mass regime drops sharply at 0.4 $M_{\odot}$ even though our data are sensitive to stellar masses down to 0.1 $M_{\odot}$, suggesting that most very-low-mass members left the cluster as the result of its dynamical evolution. In addition, the cluster is highly elongated (ellipticity $\sim$ 0.5) towards the galactic plane, and mass segregated. Conclusions. Our combined Gaia+DANCe data set allows us to obtain an extended list of cluster candidate members, and to derive luminosity, mass and projected spatial distributions in the oldest open cluster of the solar vicinity., Comment: 25 pages

There has been increasing progress toward detailed characterization of exoplanetary atmospheres, in both observations and theoretical methods. Improvements in observational facilities and data reduction and analysis techniques are enabling increasingly higher quality spectra, especially from ground-based facilities. The high data quality also necessitates concomitant improvements in models required to interpret such data. In particular, the detection of trace species such as metal oxides has been challenging. Extremely irradiated exoplanets (~3000 K) are expected to show oxides with strong absorption signals in the optical. However, there are only a few hot Jupiters where such signatures have been reported. Here we aim to characterize the atmosphere of the ultra-hot Jupiter WASP-33b using two primary transits taken 18 orbits apart. Our atmospheric retrieval, performed on the combined data sets, provides initial constraints on the atmospheric composition of WASP-33b. We report a possible indication of aluminum oxide (AlO) at 3.3-sigma significance. The data were obtained with the long slit OSIRIS spectrograph mounted at the 10-meter Gran Telescopio Canarias. We cleaned the brightness variations from the light curves produced by stellar pulsations, and we determined the wavelength-dependent variability of the planetary radius caused by the atmospheric absorption of stellar light. A simultaneous fit to the two transit light curves allowed us to refine the transit parameters, and the common wavelength coverage between the two transits served to contrast our results. Future observations with HST as well as other large ground-based facilities will be able to further constrain the atmospheric chemical composition of the planet., Comment: 15 pages, 11 figures

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., Comment: 30 pages, 21 figures, 1 online table

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., Comment: 5 pages, accepted for publication in A&A letters

The evolution of lithium abundance over a star's lifetime is indicative of transport processes operating in the stellar interior. We revisit the relationship between lithium content and rotation rate previously reported for cool dwarfs in the Pleiades cluster. We derive new LiI 670.8 nm equivalent width measurements from high-resolution spectra obtained for low-mass Pleiades members. We combine these new measurements with previously published ones, and use the Kepler/K2 rotational periods recently derived for Pleiades cool dwarfs to investigate the lithium-rotation connection in this 125 Myr-old cluster. The new data confirm the correlation between lithium equivalent width and stellar spin rate for a sample of 51 early K-type members of the cluster, where fast rotating stars are systematically lithium-rich compared to slowly rotating ones. The correlation is valid for all stars over the (J-Ks) color range 0.50-0.70 mag, corresponding to a mass range from about 0.75 to 0.90 solar mass, and may extend down to lower masses. We argue that the dispersion in lithium equivalent widths observed for cool dwarfs in the Pleiades cluster reflects an intrinsic scatter in lithium abundances, and suggest that the physical origin of the lithium dispersion pattern is to be found in the pre-main sequence rotational history of solar-type stars., Comment: 10 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

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)., Comment: Accepted for publication in A&A

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., Comment: 39 pages, 31 figures

We describe the methodologies that, taking advantage of Gaia-DR1 and the Gaia-ESO Survey data, enable the comparison of observed open star cluster sequences with stellar evolutionary models. The final, long-term goal is the exploitation of open clusters as age calibrators. We perform a homogeneous analysis of eight open clusters using the Gaia-DR1 TGAS catalogue for bright members, and information from the Gaia-ESO Survey for fainter stars. Cluster membership probabilities for the Gaia-ESO Survey targets are derived based on several spectroscopic tracers. The Gaia-ESO Survey also provides the cluster chemical composition. We obtain cluster parallaxes using two methods. The first one relies on the astrometric selection of a sample of bona fide members, while the other one fits the parallax distribution of a larger sample of TGAS sources. Ages and reddening values are recovered through a Bayesian analysis using the 2MASS magnitudes and three sets of standard models. Lithium depletion boundary (LDB) ages are also determined using literature observations and the same models employed for the Bayesian analysis. For all but one cluster, parallaxes derived by us agree with those presented in Gaia Collaboration et al. (2017), while a discrepancy is found for NGC 2516; we provide evidence supporting our own determination. Inferred cluster ages are robust against models and are generally consistent with literature values. The systematic parallax errors inherent in the Gaia DR1 data presently limit the precision of our results. Nevertheless, we have been able to place these eight clusters onto the same age scale for the first time, with good agreement between isochronal and LDB ages where there is overlap. Our approach appears promising and demonstrates the potential of combining Gaia and ground-based spectroscopic datasets., Comment: Accepted for publication by A&A. 22 pages main text plus 8 pages appendices, 12 figures main text, 7 figures appendices

The search for minor bodies in the solar system promises insights into its formation history. Wide imaging surveys offer the opportunity to serendipitously discover and identify these traces of planetary formation and evolution. We aim to present a method to acquire position, photometry, and proper motion measurements of solar system objects in surveys using dithered image sequences. The application of this method on the Kilo-Degree Survey is demonstrated. Optical images of 346 square degree fields of the sky are searched in up to four filters using the AstrOmatic software suite to reduce the pixel to catalog data. The solar system objects within the acquired sources are selected based on a set of criteria depending on their number of observation, motion, and size. The Virtual Observatory SkyBoT tool is used to identify known objects. We observed 20,221 SSO candidates, with an estimated false-positive content of less than 0.05%. Of these SSO candidates, 53.4% are identified by SkyBoT. KiDS can detect previously unknown SSOs because of its depth and coverage at high ecliptic latitude, including parts of the Southern Hemisphere. Thus we expect the large fraction of the 46.6% of unidentified objects to be truly new SSOs. Our method is applicable to a variety of dithered surveys such as DES, LSST, and Euclid. It offers a quick and easy-to-implement search for solar system objects. SkyBoT can then be used to estimate the completeness of the recovered sample., Comment: Final version to be published in A&A. Supplementary data will be published at CDS. 11 pages, 15 figures

One-opposition near-Earth asteroids (NEAs) are growing in number, and they must be recovered to prevent loss and mismatch risk, and to improve their orbits, as they are likely to be too faint for detection in shallow surveys at future apparitions. We aimed to recover more than half of the one-opposition NEAs recommended for observations by the Minor Planet Center (MPC) using the Isaac Newton Telescope (INT) in soft-override mode and some fractions of available D-nights. During about 130 hours in total between 2013 and 2016, we targeted 368 NEAs, among which 56 potentially hazardous asteroids (PHAs), observing 437 INT Wide Field Camera (WFC) fields and recovering 280 NEAs (76% of all targets). Engaging a core team of about ten students and amateurs, we used the THELI, Astrometrica, and the Find_Orb software to identify all moving objects using the blink and track-and-stack method for the faintest targets and plotting the positional uncertainty ellipse from NEODyS. Most targets and recovered objects had apparent magnitudes centered around V~22.8 mag, with some becoming as faint as V~24 mag. One hundred and three objects (representing 28% of all targets) were recovered by EURONEAR alone by Aug 2017. Orbital arcs were prolonged typically from a few weeks to a few years; our oldest recoveries reach 16 years. The O-C residuals for our 1,854 NEA astrometric positions show that most measurements cluster closely around the origin. In addition to the recovered NEAs, 22,000 positions of about 3,500 known minor planets and another 10,000 observations of about 1,500 unknown objects (mostly main-belt objects) were promptly reported to the MPC by our team. Four new NEAs were discovered serendipitously in the analyzed fields, increasing the counting to nine NEAs discovered by the EURONEAR in 2014 and 2015., Comment: Accepted for publication in Astronomy and Astrophysics (11 Oct 2017). Version 2 adding two co-authors and fixing the affiliation page overflow

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., Comment: 29 pages, 25 figures. Accepted for publication by A&A

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., Comment: Accepted for publication by A&A. 21 pages main text plus 46 pages appendices. 34 figures main text, 38 figures appendices. 8 table in main text, 19 tables in appendices

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., Comment: Accepted for publication on Astronomy and Astrophysics

Aims: We provide a new distance estimate for the Pleiades based on the moving cluster method, which will be useful to further discuss the so-called Pleiades distance controversy and compare it with the very precise parallaxes from the Gaia space mission. Methods: We apply a refurbished implementation of the convergent point search method to an updated census of Pleiades stars to calculate the convergent point position of the cluster from stellar proper motions. Then, we derive individual parallaxes for 64 cluster members using radial velocities compiled from the literature, and approximate parallaxes for another 1146 stars based on the spatial velocity of the cluster. This represents the largest sample of Pleiades stars with individual distances to date. Results: The parallaxes derived in this work are in good agreement with previous results obtained in different studies (excluding Hipparcos) for individual stars in the cluster. We report a mean parallax of $7.44\pm 0.08$~mas and distance of $134.4^{+2.9}_{-2.8}$pc that is consistent with the weighted mean of $135.0\pm 0.6$pc obtained from the non-Hipparcos results in the literature. Conclusions: Our result for the distance to the Pleiades open cluster is not consistent with the Hipparcos catalog, but favors the recent and more precise distance determination of $136.2\pm 1.2$pc obtained from Very Long Baseline Interferometry observations. It is also in good agreement with the mean distance of $133\pm 5$pc obtained from the first trigonometric parallaxes delivered by the Gaia satellite for the brightest cluster members in common with our sample., Comment: 22 pages, 20 figures. Accepted for publication in A&A (in press)

We aim at performing a kinematic census of young stellar objects (YSOs) in the Rho Ophiuchi F core and partially in the E core of the L1688 dark cloud. We run a proper motion program at the ESO New Technology Telescope (NTT) with the Son of ISAAC (SOFI) instrument over nine years in the near-infrared. We complemented these observations with various public image databases to enlarge the time base of observations and the field of investigation to 0.5 deg X 0.5 deg. We derived positions and proper motions for 2213 objects. From these, 607 proper motions were derived from SOFI observations with a ~1.8 mas/yr accuracy while the remaining objects were measured only from auxiliary data with a mean precision of about ~3 mas/yr. We performed a kinematic analysis of the most accurate proper motions derived in this work, which allowed us to separate cluster members from field stars and to derive the mean properties of the cluster. From the kinematic analysis we derived a list of 68 members and 14 candidate members, comprising 26 new objects with a high membership probability. These new members are generally fainter than the known ones. We measured a mean proper motion of (mu_RA*, mu_DEC)=(-8.2, -24.3)+/-0.8 mas/yr for the L1688 dark cloud. A supervised classification was applied to photometric data of members to allocate a spectral energy distribution (SED) classification to the unclassified members.} We kinematically confirmed that the 56 members that were known from previous studies of the Rho Ophiuchi F cluster and that were also part of our survey are members of the cluster, and we added 26 new members. We defined the evolutionary status of the unclassified members of the cluster. We showed that a large part (23) of these new members are probably brown dwarfs, which multiplies the number of known substellar objects in the cluster by a factor of 3.3., Comment: 16 pages, 13 figures, 3 tables in electronic form

This paper continues our study of the foreground population to the Orion molecular clouds. The goal is to characterize the foreground population north of NGC 1981 and to investigate the star formation history in the large Orion star-forming region. We focus on a region covering about 25 square degrees, centered on the $\epsilon$ Orionis supergiant (HD 37128, B0\,Ia) and covering the Orion Belt asterism. We used a combination of optical (SDSS) and near-infrared (2MASS) data, informed by X-ray (\textit{XMM-Newton}) and mid-infrared (WISE) data, to construct a suite of color-color and color-magnitude diagrams for all available sources. We then applied a new statistical multiband technique to isolate a previously unknown stellar population in this region. We identify a rich and well-defined stellar population in the surveyed region that has about 2\,000 objects that are mostly M stars. We infer the age for this new population to be at least 5\, Myr and likely $\sim10$\,Myr and estimate a total of about 2\,500 members, assuming a normal IMF. This new population, which we call the Orion Belt population, is essentially extinction-free, disk-free, and its spatial distribution is roughly centered near $\epsilon$ Ori, although substructure is clearly present. The Orion Belt population is likely the low-mass counterpart to the Ori OB Ib subgroup. Although our results do not rule out Blaauw's sequential star formation scenario for Orion, we argue that the recently proposed blue streams scenario provides a better framework on which one can explain the Orion star formation region as a whole. We speculate that the Orion Belt population could represent the evolved counterpart of an Orion nebula-like cluster., Comment: 13 pages, 10 figures. Accepted for publication to A&A

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., Comment: Accepted Astronomy & Astrophysics, three big tables submitted to CDS

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., Comment: Accepted for publication in A&A on May 18th 2016, 9 pages, 5 figures, 4 tables. The first two authors contributed equally to the work contained in this paper

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., Comment: 26 pages, 15 figures; electronic-only tables and figures available upon request to the author. Accepted by AJ

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., Comment: 30 pages, 26 figures, accepted to AJ. Paper I is arXiv:1606.00052 and Paper II is arXiv:1606.00055

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., Comment: 28 pages, 15 figures. Accepted by AJ

Transitional disks are protoplanetary disks with opacity gaps/cavities in their dust distribution, a feature that may be linked to planet formation. We perform Bayesian modeling of the three transitional disks SZ Cha, CS Cha and T25 including photometry from the Herschel Space Observatory to quantify the improvements added by these new data. We find disk dust masses between 2x10^-5 and 4x10^-4 Msun, and gap radii in the range of 7-18 AU, with uncertainties of ~ one order of magnitude and ~ 4 AU, respectively. Our results show that adding Herschel data can significantly improve these estimates with respect to mid-infrared data alone, which have roughly twice as large uncertainties on both disk mass and gap radius. We also find weak evidence for different density profiles with respect to full disks. These results open exciting new possibilities to study the distribution of disk masses for large samples of disks., Comment: 10 pages, 3 figures. Accepted for publication in MNRAS

Context. Chamaeleon II molecular cloud is an active star forming region that offers an excellent opportunity for studying the formation of brown dwarfs in the southern hemisphere. Aims. Our aims are to identify a population of pre- and proto- brown dwarfs (5 sigma mass limit threshold of ~0.015 Msun) and provide information on the formation mechanisms of substellar objects. Methods. We performed high sensitivity observations at 870 microns using the LABOCA bolometer at the APEX telescope towards an active star forming region in Chamaeleon II. The data are complemented with an extensive multiwavelength catalogue of sources from the optical to the far-infrared to study the nature of the LABOCA detections. Results. We detect fifteen cores at 870 microns, and eleven of them show masses in the substellar regime. The most intense objects in the surveyed field correspond to the submillimeter counterparts of the well known young stellar objects DK Cha and IRAS 12500-7658. We identify a possible proto-brown dwarf candidate (ChaII-APEX-L) with IRAC emission at 3.6 and 4.5 microns. Conclusions. Our analysis indicates that most of the spatially resolved cores are transient, and that the point-like starless cores in the sub-stellar regime (with masses between 0.016 Msun and 0.066 Msun) could be pre-brown dwarfs cores gravitationally unstable if they have radii smaller than 220 AU to 907 AU (1.2" to 5" at 178 pc) respectively for different masses. ALMA observations will be the key to reveal the energetic state of these pre-brown dwarfs candidates., Comment: 11 pages, 6 figures

We present optical depth and temperature maps of the Perseus molecular cloud, obtained combining dust emission data from the Herschel and Planck satellites and 2MASS/NIR dust extinction maps. The maps have a resolution of 36 arcsec in the Herschel regions, and of 5 arcmin elsewhere. The dynamic range of the optical depth map ranges from $1\times10^{-2}\, \mathrm{mag}$ up to $20 \,\mathrm{mag}$ in the equivalent K band extinction. We also evaluate the ratio between the $2.2 \,\mathrm{\mu m}$ extinction coefficient and the $850 \,\mathrm{\mu m}$ opacity. The value we obtain is close to the one found in the Orion B molecular cloud. We show that the cumulative and the differential area function of the data (which is proportional to the probability distribution function of the cloud column density) follow power laws with index respectively $\simeq -2$, and $\simeq -3$. We use WISE data to improve current YSO catalogues based mostly on \emph{Spitzer} data and we build an up-to-date selection of Class~I/0 objects. Using this selection, we evaluate the local Schmidt law, $\Sigma_{\mathrm{YSO}} \propto \Sigma_{\mathrm{gas}}^{\beta}$, showing that $\beta=2.4 \pm 0.6$. Finally, we show that the area-extinction relation is important for determining the star formation rate in the cloud, which is in agreement with other recent works., Comment: 15 pages, 19 figures. Accepted for publication in Astronomy & Astrophysics