ESA recently called for new "Science Ideas" to be investigated in terms of feasibility and technological developments -- for technologies not yet sufficiently mature. These ideas may in the future become candidates for M or L class missions within the ESA Science Program. With the launch of Gaia in December 2013, Europe entered a new era of space astrometry following in the footsteps of the very successful Hipparcos mission from the early 1990s. Gaia is the successor to Hipparcos, both of which operated in optical wavelengths, and Gaia is two orders of magnitude more accurate in the five astrometric parameters and is surveying four orders of magnitude more stars in a vast volume of the Milky Way. The combination of the Hipparcos/Tycho-2 catalogues with the first early Gaia data release will give improved proper motions over a long ~25 year baseline. The final Gaia solution will also establish a new optical reference frame by means of quasars, by linking the optical counterparts of radio (VLBI) sources defining the orientation of the reference frame, and by using the zero proper motion of quasars to determine a non-rotating frame. A weakness of Gaia is that it only operates at optical wavelengths. However, much of the Galactic centre and the spiral arm regions, important for certain studies, are obscured by interstellar extinction and this makes it difficult for Gaia to deeply probe. Traditionally, this problem is overcome by switching to the infra-red but this was not possible with Gaia's CCDs. Additionally, to scan the entire sky and make global absolute parallax measurements the spacecraft must have a constant rotation and this requires that the CCDs operate in TDI mode, increasing their complexity.

From the approximately $\sim$3,500 planetary nebulae (PNe) discovered in our Galaxy, only 14 are known to be members of the Galactic halo. Nevertheless, a systematic search for halo PNe has never been performed. In this study, we present new photometric diagnostic tools to identify compact PNe in the Galactic halo by making use of the novel 12-filter system projects, J-PLUS (Javalambre Photometric Local Universe Survey) and S-PLUS (Southern-Photometric Local Universe Survey). We reconstructed the IPHAS (Isaac Newton Telescope (INT) Photometric H${\alpha}$ Survey of the Northern Galactic Plane) diagnostic diagram and propose four new ones using i) the J-PLUS and S-PLUS synthetic photometry for a grid of photo-ionisation models of halo PNe, ii) several observed halo PNe, as well as iii) a number of other emission-line objects that resemble PNe. All colour-colour diagnostic diagrams are validated using two known halo PNe observed by J-PLUS during the scientific verification phase and the first data release (DR1) of S-PLUS and the DR1 of J-PLUS. By applying our criteria to the DR1s ($\sim$1,190 deg$^2$), we identified one PN candidate. However, optical follow-up spectroscopy proved it to be a H II region belonging to the UGC 5272 galaxy. Here, we also discuss the PN and two H II galaxies recovered by these selection criteria. Finally, the cross-matching with the most updated PNe catalogue (HASH) helped us to highlight the potential of these surveys, since we recover all the known PNe in the observed area. The tools here proposed to identify PNe and separate them from their emission-line contaminants proved to be very efficient thanks to the combination of many colours, even when applied -like in the present work- to an automatic photometric search that is limited to compact PNe.

We present the first release of GTC OSIRIS Broad Band data archive. This is an effort conducted in the framework of the Spanish Virtual Observatory to help optimize science from the Gran Telescopio Canarias Archive. Data Release 1 includes 6 788 broad-band images in the Sloan griz filters obtained between April 2009 and January 2014 and the associated catalogue with roughly 6.23 million detections of more than 630 000 unique sources. The catalogue contains standard PSF and Kron aperture photometry with a mean accuracy better than 0.09 and 0.15 mag, respectively. The relative astrometric residuals are always better than 30 mas and better than 15 mas in most cases. The absolute astrometric uncertainty of the catalogue is of 0.12 arcsec. In this paper we describe the procedure followed to build the image archive and the associated catalogue, as well as the quality tests carried out for validation. To illustrate some of the scientific potential of the catalogue, we also provide two examples of its scientific exploitation: discovery and identification of asteroids and cool dwarfs.

The nature and physical properties of asteroids, in particular those orbiting in the near-Earth space, are of scientific interest and practical importance. Exoplanet surveys can be excellent resources to detect asteroids, both already known and new objects. This is due their similar observing requirements: large fields of view, long sequences, and short cadence. If the targeted fields are not located far from the ecliptic, many asteroids will cross occasionally the field of view. We present two complementary methodologies to identify asteroids serendipitously observed in large-area astronomical surveys. One methodology focuses on detecting already known asteroids using the Virtual Observatory tool SkyBoT, which predicts their positions and motions in the sky at a specific epoch. The other methodology applies the ssos pipeline, which is able to identify known and new asteroids based on their apparent motion. The application of these methods to the 6.4 deg 2 of the sky covered by the Wide-Field CAMera Transit Survey in the J-band is described. We identified 15 661 positions of 1 821 different asteroids. Of them, 182 are potential new discoveries. A publicly accessible online, Virtual Observatory compliant catalogue was created. We obtained the shapes and periods for five of our asteroids from their light-curves built with additional photometry taken from external archives. We demonstrated that our methodologies are robust and reliable approaches to find, at zero cost of observing time, asteroids observed by chance in astronomical surveys. Our future goal is to apply them to other surveys with adequate temporal coverage.

We analyse the 100pc Gaia white dwarf volume-limited sample by means of VOSA (Virtual Observatory SED Analyser) with the aim of identifying candidates for displaying infrared excesses. Our search focuses on the study of the spectral energy distribution (SED) of 3,733 white dwarfs with reliable infrared photometry and GBP-GRP colours below 0.8 mag, a sample which seems to be nearly representative of the overall white dwarf population. Our search results in 77 selected candidates, 52 of which are new identifications. For each target we apply a two-component SED fitting implemented in VOSA to derive the effective temperatures of both the white dwarf and the object causing the excess. We calculate a fraction of infrared-excess white dwarfs due to the presence of a circumstellar disk of 1.6+-0.2%, a value which increases to 2.6+-0.3% if we take into account incompleteness issues. Our results are in agreement with the drop in the percentage of infrared excess detections for cool (<8,000K) and hot (>20,000K) white dwarfs obtained in previous analyses. The fraction of white dwarfs with brown dwarf companions we derive is ~0.1-0.2%.

We present the photometric calibration of the twelve optical passbands observed by the Javalambre Photometric Local Universe Survey (J-PLUS). The proposed calibration method has four steps: (i) definition of a high-quality set of calibration stars using Gaia information and available 3D dust maps; (ii) anchoring of the J-PLUS gri passbands to the Pan-STARRS photometric solution, accounting for the variation of the calibration with the position of the sources on the CCD; (iii) homogenization of the photometry in the other nine J-PLUS filters using the dust de-reddened instrumental stellar locus in (X - r) versus (g - i) colours, where X is the filter to calibrate. The zero point variation along the CCD in these filters was estimated with the distance to the stellar locus. Finally, (iv) the absolute colour calibration was obtained with the white dwarf locus. We performed a joint Bayesian modelling of eleven J-PLUS colour-colour diagrams using the theoretical white dwarf locus as reference. This provides the needed offsets to transform instrumental magnitudes to calibrated magnitudes outside the atmosphere. The uncertainty of the J-PLUS photometric calibration, estimated from duplicated objects observed in adjacent pointings and accounting for the absolute colour and flux calibration errors, are ~19 mmag in u, J0378 and J0395, ~11 mmag in J0410 and J0430, and ~8 mmag in g, J0515, r, J0660, i, J0861, and z. We present an optimized calibration method for the large area multi-filter J-PLUS project, reaching 1-2% accuracy within an area of 1 022 square degrees without the need for long observing calibration campaigns or constant atmospheric monitoring. The proposed method will be adapted for the photometric calibration of J-PAS, that will observe several thousand square degrees with 56 narrow optical filters.

A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost 2 billion common stars the combination of Gaia and a new all-sky NIR astrometry mission would provide much improved proper motions, answering key science questions -- from the solar system and stellar systems, including exoplanet systems, to compact galaxies, quasars, neutron stars, binaries and dark matter substructures. The addition of NIR will result in up to 8 billion newly measured stars in some of the most obscured parts of our Galaxy, and crucially reveal the very heart of the Galactic bulge region. In this white paper we argue that rather than improving on the accuracy, a greater overall science return can be achieved by going deeper than Gaia and by expanding the wavelength range to the NIR.

Ultracool dwarfs (UCDs) are objects with spectral types equal or later than M7. Most of them have been discovered using wide-field imaging surveys. The Virtual Observatory (VO) has proven to be of great utility to efficiently exploit these astronomical resources. We aim to validate a VO methodology designed to discover and characterize UCDs in the J-PLUS photometric survey. J-PLUS is a multiband survey carried out with the wide angle T80Cam optical camera mounted on the 0.83-m telescope JAST/T80 in the Observatorio Astrof\'isico de Javalambre. In this work we make use of the Internal Data Release (IDR) covering 528 deg$^2$. We complement J-PLUS photometry with other catalogues in the optical and IR using VOSA, a VO tool that estimates physical parameters from the spectral energy distribution fitting to collections of theoretical models. Objects identified as UCDs are distinguished from background M giants and highly reddened stars using parallaxes and proper motions from Gaia DR2. We identify 559 UCDs, ranging from i=16.2 to 22.4 mag, of which 187 are candidate UCDs not previously reported in the literature. This represents an increase in the number of known UCDs of about 50% in the studied region of the sky, particularly at the faint end of our sensitivity, which is interesting as reference for future wide and deep surveys such as Euclid. Three candidates constitute interesting targets for exoplanet surveys because of their proximity (<40 pc). We also analyze the kinematics of UCDs in our catalogue and find evidence that it is consistent with a Galactic thin-disk population, except for 6 objects that might be members of the thick disk. The results obtained validate the proposed methodology, which will be used in future J-PLUS and J-PAS releases. Considering the region of the sky covered by the IDR used, we foresee to discover 3,000-3,500 new UCDs at the end of the J-PLUS project.

Gaia-DR2 has provided an unprecedented number of white dwarf candidates of our Galaxy. In particular, it is estimated that Gaia-DR2 has observed nearly 400,000 of these objects and close to 18,000 up to 100 pc from the Sun. This large quantity of data requires a thorough analysis in order to uncover their main Galactic population properties, in particular the thin and thick disk and halo components. Taking advantage of recent developments in artificial intelligence techniques, we make use of a detailed Random Forest algorithm to analyse an 8-dimensional space (equatorial coordinates, parallax, proper motion components and photometric magnitudes) of accurate data provided by Gaia-DR2 within 100 pc from the Sun. With the aid of a thorough and robust population synthesis code we simulated the different components of the Galactic white dwarf population to optimize the information extracted from the algorithm for disentangling the different population components. The algorithm is first tested in a known simulated sample achieving an accuracy of 85.3%. Our methodology is thoroughly compared to standard methods based on kinematic criteria demonstrating that our algorithm substantially improves previous approaches. Once trained, the algorithm is then applied to the Gaia-DR2 100 pc white dwarf sample, identifying 12,227 thin disk, 1,410 thick disk and 95 halo white dwarf candidates, which represent a proportion of 74:25:1, respectively. Hence, the numerical spatial densities are $(3.6\pm0.4)\times10^{-3}\,{\rm pc^{-3}}$, $(1.2\pm0.4)\times10^{-3}\,{\rm pc^{-3}}$ and $(4.8\pm0.4)\times10^{-5}\,{\rm pc^{-3}}$ for the thin disk, thick disk and halo components, respectively. The populations thus obtained represent the most complete and volume-limited samples to date of the different components of the Galactic white dwarf population.

Binary and multiple stars have long provided an effective empirical method of testing stellar formation and evolution theories. In particular, the existence of wide binary systems (separations $>$20,000 au) is particularly challenging to binary formation models as their physical separations are beyond the typical size of a collapsing cloud core ($\sim$5,000$-$10,000 au). We mined the recently published Gaia-DR2 catalog to identify bright comoving systems in the five-dimensional space (sky position, parallax, and proper motion). We identified 3,741 comoving binary and multiple stellar candidate systems, out of which 575 have compatible RVs for all the members of the system. The candidate systems have separations between $\sim$400 and 500,000 au. We used the analysis tools of the Virtual Observatory to characterize the comoving system members and to assess their reliability. The comparison with previous comoving systems catalogs obtained from TGAS showed that these catalogs contain a large number of false systems. In addition, we were not able to confirm the ultra-wide binary population presented in these catalogs. The robustness of our methodology is demonstrated by the identification of well known comoving star clusters and by the low contamination rate for comoving binary systems with projected physical separations $<$50,000 au. These last constitute a reliable sample for further studies. The catalog is available online at the Spanish Virtual Observatory portal (http://svo2.cab.inta-csic.es/vocats/v2/comovingGaiaDR2/).

In very young clusters, stellar age distribution is the empirical proof of the duration of star formation (SF) and of the physical mechanisms involved in the process. We derived accurate stellar ages for the cluster NGC6530, associated with the Lagoon Nebula to infer its SF history. We use the Gaia-ESO survey observations and Gaia DR2 data, to derive cluster membership and fundamental stellar parameters. We identified 652 confirmed and 9 probable members. The reddening inferred for members and non-members allows us to distinguish MS stars and giants, in agreement with the distances inferred from Gaia DR2 data. The foreground and background stars show a spatial pattern that traces the 3D structure of the nebular dust component. We derive stellar ages for 382 confirmed cluster members and we find that the gravity-sensitive gamma index distribution for M stars is correlated with stellar age. For all members with Teff<5500 K, the mean logarithmic age is 5.84 (units of years) with a dispersion of 0.36 dex. The age distribution of stars with accretion and/or disk (CTTSe) is similar to that of stars without accretion and without disk (WTTSp). We interpret this dispersion as evidence of a real age spread since the total uncertainties on age determinations, derived from Monte Carlo simulations, are significantly smaller than the observed spread. This conclusion is supported by the evidence of a decreasing of the gravity-sensitive gamma index as a function of stellar ages. The presence of the age spread is also supported by the spatial distribution and the kinematics of old and young members. In particular, members with accretion and/or disk, formed in the last 1 Myr, show evidence of subclustering around the cluster center, in the Hourglass Nebula and in the M8-E region, suggesting a possible triggering of star formation events by the O-type star ionization fronts.

The GES survey using FLAMES at the VLT has obtained high-resolution UVES spectra for a large number of giant stars, allowing a determination of the abundances of the key chemical elements C and N at their surface. The surface abundances of these chemical species are well-known to change in stars during their evolution on the red giant branch after the first dredge-up episod, as a result of extra-mixing phenomena. We investigate the effects of thermohaline mixing on C and N abundances using the first comparison between the GES [C/N] determinations with simulations of the observed fields using a model of stellar population synthesis. We explore the effects of thermohaline mixing on the chemical properties of giants through stellar evolutionary models computed with the stellar evolution code STAREVOL. We include these stellar evolution models in the Besan\c{c}on Galaxy model to simulate the [C/N] distributions determined from the UVES spectra of the GES and compare them with the observations. Theoretical predictions including the effect of thermohaline mixing are in good agreement with the observations. However, the field stars in the GES with C and N-abundance measurements have a metallicity close to solar, where the efficiency of thermohaline mixing is not very large. The C and N abundances derived by the GES in open and globular clusters clearly show the impact of thermohaline mixing at low-metallicity, allowing to explain the [C/N] ratio observed in lower-mass and older giant stars. Using independent observations of carbon isotopic ratio in clump field stars and open clusters, we also confirm that thermohaline mixing should be taken into account to explain the behavior of 12C/13C ratio as a function of stellar age. Overall the current model including thermohaline mixing is able to reproduce very well the C- and N-abundances over the whole metallicity range investigated by the GES data.

We present a catalogue of 73,221 white dwarf candidates extracted from the astrometric and photometric data of the recently published Gaia DR2 catalogue. White dwarfs were selected from the Gaia Hertzsprung-Russell diagram with the aid of the most updated population synthesis simulator. Our analysis shows that Gaia has virtually identified all white dwarfs within 100 pc from the Sun. Hence, our sub-population of 8,555 white dwarfs within this distance limit and the colour range considered, $-\,0.52<(G_{\rm BP}-G_{\rm RP})<0.80$, is the largest and most complete volume-limited sample of such objects to date. From this sub-sample we identified 8,343 CO-core and 212 ONe-core white dwarf candidates and derived a white dwarf space density of $4.9\pm0.4\times10^{-3}\,{\rm pc^{-3}}$. A bifurcation in the Hertzsprung-Russell diagram for these sources, which our models do not predict, is clearly visible. We used the Virtual Observatory tool VOSA to derive effective temperatures and luminosities for our sources by fitting their spectral energy distributions, that we built from the UV to the NIR using publicly available photometry through the Virtual Observatory. From these parameters, we derived the white dwarf radii. Interpolating the radii and effective temperatures in hydrogen-rich white dwarf cooling sequences, we derived the surface gravities and masses. The Gaia 100 pc white dwarf population is clearly dominated by cool ($\sim$ 8,000 K) objects and reveals a significant population of massive ($M \sim 0.8 M_{\odot}$) white dwarfs, of which no more than $\sim$ $30-40 \%$ can be attributed to hydrogen-deficient atmospheres, and whose origin remains uncertain.

Context. The origin and dynamical evolution of star clusters is an important topic in stellar astrophysics. Several models have been proposed to understand the formation of bound and unbound clusters and their evolution, and these can be tested by examining the kinematical and dynamical properties of clusters over a wide range of ages and masses. Aims. We use the Gaia-ESO Survey products to study four open clusters (IC 2602, IC 2391, IC 4665, and NGC 2547) that lie in the age range between 20 and 50 Myr. Methods. We employ the gravity index $\gamma$ and the equivalent width of the lithium line at 6708 $\AA$, together with effective temperature $\rm{T_{eff}}$, and the metallicity of the stars in order to discard observed contaminant stars. Then, we derive the cluster radial velocity dispersions $\sigma_c$, the total cluster mass $\rm{M}_{tot}$, and the half mass radius $r_{hm}$. Using the $Gaia$-DR1 TGAS catalogue, we independently derive the intrinsic velocity dispersion of the clusters from the astrometric parameters of cluster members. Results. The intrinsic radial velocity dispersions derived by the spectroscopic data are larger than those derived from the TGAS data, possibly due to the different masses of the considered stars. Using $\rm{M}_{tot}$ and $r_{hm}$ we derive the virial velocity dispersion $\sigma_{vir}$ and we find that three out of four clusters are supervirial. This result is in agreement with the hypothesis that these clusters are dispersing, as predicted by the "residual gas expulsion" scenario. However, recent simulations show that the virial ratio of young star clusters may be overestimated if it is determined using the global velocity dispersion, since the clusters are not fully relaxed.

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.

Milky Way open clusters are very diverse in terms of age, chemical composition, and kinematic properties. Intermediate-age and old open clusters are less common, and it is even harder to find them inside the solar Galactocentric radius, due to the high mortality rate and strong extinction inside this region. NGC 6802 is one of the inner disk open clusters (IOCs) observed by the $Gaia$-ESO survey (GES). This cluster is an important target for calibrating the abundances derived in the survey due to the kinematic and chemical homogeneity of the members in open clusters. Using the measurements from $Gaia$-ESO internal data release 4 (iDR4), we identify 95 main-sequence dwarfs as cluster members from the GIRAFFE target list, and eight giants as cluster members from the UVES target list. The dwarf cluster members have a median radial velocity of $13.6\pm1.9$ km s$^{-1}$, while the giant cluster members have a median radial velocity of $12.0\pm0.9$ km s$^{-1}$ and a median [Fe/H] of $0.10\pm0.02$ dex. The color-magnitude diagram of these cluster members suggests an age of $0.9\pm0.1$ Gyr, with $(m-M)_0=11.4$ and $E(B-V)=0.86$. We perform the first detailed chemical abundance analysis of NGC 6802, including 27 elemental species. To gain a more general picture about IOCs, the measurements of NGC 6802 are compared with those of other IOCs previously studied by GES, that is, NGC 4815, Trumpler 20, NGC 6705, and Berkeley 81. NGC 6802 shows similar C, N, Na, and Al abundances as other IOCs. These elements are compared with nucleosynthetic models as a function of cluster turn-off mass. The $\alpha$, iron-peak, and neutron-capture elements are also explored in a self-consistent way.

The Gaia ESO Public Spectroscopic Survey (GES) is providing the astronomical community with high-precision measurements of many stellar parameters including radial velocities (RVs) of stars belonging to several young clusters and star-forming regions. One of the main goals of the young cluster observations is to study of their dynamical evolution and provide insight into their future, revealing if they will eventually disperse to populate the field, rather than evolve into bound open clusters. In this paper we report the analysis of the dynamical state of L1688 in the $\rho$~Ophiuchi molecular cloud using the dataset provided by the GES consortium. We performed the membership selection of the more than 300 objects observed. Using the presence of the lithium absorption and the location in the Hertzspung-Russell diagram, we identify 45 already known members and two new association members. We provide accurate RVs for all 47 confirmed members.A dynamical analysis, after accounting for unresolved binaries and errors, shows that the stellar surface population of L1688 has a velocity dispersion $\sigma \sim$1.14$\pm$0.35 km s$^{-1}$ that is consistent with being in virial equilibrium and is bound with a $\sim$80% probability. We also find a velocity gradient in the stellar surface population of $\sim$1.0 km s$^{-1}$pc$^{-1}$ in the northwest/southeast direction, which is consistent with that found for the pre-stellar dense cores, and we discuss the possibility of sequential and triggered star formation in L1688.

Context. Extremely reddened AGB stars lose mass at high rates of >10^-5 Msun/yr. This is the very last stage of AGB evolution, in which stars in the mass range 2.0--4.0 Msun (for solar metallicity) should have been converted to C stars already. The extremely reddened AGB stars in the Galactic bulge are however predominantly O-rich, implying that they might be either low-mass stars or stars at the upper end of the AGB mass range. Aims. To determine the mass range of the most reddened AGB stars in the Galactic bulge. Methods. Using Virtual Observatory tools, we constructed spectral energy distributions of a sample of 37 evolved stars in the Galactic bulge with extremely red IRAS colours. We fitted DUSTY models to the observational data to infer the bolometric fluxes. Applying individual corrections for interstellar extinction and adopting a common distance, we determined luminosities and mass-loss rates, and inferred the progenitor mass range from comparisons with AGB evolutionary models. Results. The observed spectral energy distributions are consistent with a classification as reddened AGB stars, except for two stars, which are proto-planetary nebula candidates. For the AGB stars, we found luminosities in the range 3000--30,000 Lsun and mass-loss rates 10^-5--3x10^-4 Msun/yr. The corresponding mass range is 1.1--6.0 Msun assuming solar metallicity. Conclusions. Contrary to the predictions of the evolutionary models, the luminosity distribution is continuous, with many O-rich AGB stars in the mass range in which they should have been converted into C stars already. We suspect that bulge AGB stars have higher than solar metallicity and therefore may avoid the conversion to C-rich. The presence of low-mass stars in the sample shows that their termination of the AGB evolution also occurs during a final phase of very high mass-loss rate, leading to optically thick circumstellar shells.

Context. Stars are born together from giant molecular clouds and, if we assume that the priors were chemically homogeneous and well-mixed, we expect them to share the same chemical composition. Most of the stellar aggregates are disrupted while orbiting the Galaxy and most of the dynamic information is lost, thus the only possibility of reconstructing the stellar formation history is to analyze the chemical abundances that we observe today. Aims. The chemical tagging technique aims to recover disrupted stellar clusters based merely on their chemical composition. We evaluate the viability of this technique to recover co-natal stars that are no longer gravitationally bound. Methods. Open clusters are co-natal aggregates that have managed to survive together. We compiled stellar spectra from 31 old and intermediate-age open clusters, homogeneously derived atmospheric parameters, and 17 abundance species, and applied machine learning algorithms to group the stars based on their chemical composition. This approach allows us to evaluate the viability and efficiency of the chemical tagging technique. Results. We found that stars at different evolutionary stages have distinct chemical patterns that may be due to NLTE effects, atomic diffusion, mixing, and biases. When separating stars into dwarfs and giants, we observed that a few open clusters show distinct chemical signatures while the majority show a high degree of overlap. This limits the recovery of co-natal aggregates by applying the chemical tagging technique. Nevertheless, there is room for improvement if more elements are included and models are improved.

Aim of this work is to determine C, N, and O abundances in stars of Galactic open clusters of the Gaia-ESO survey and to compare the observed abundances with those predicted by current stellar and Galactic evolution models. In this pilot paper, we investigate the first three intermediate-age open clusters. High-resolution spectra, observed with the FLAMES-UVES spectrograph on the ESO VLT telescope, were analysed using a differential model atmosphere method. Abundances of carbon were derived using the C2 band heads at 5135 and 5635.5 {\AA}. The wavelength interval 6470- 6490 {\AA}, with CN features, was analysed to determine nitrogen abundances. Oxygen abundances were determined from the [O i] line at 6300 {\AA}. The mean values of the elemental abundances in Trumpler 20 as determined from 42 stars are: [Fe/H] = 0.10 +- 0.08 (s.d.), [C/H] = -0.10 +- 0.07, [N/H] = 0.50 +- 0.07, and consequently C/N = 0.98 +- 0.12. We measure from five giants in NGC 4815: [Fe/H] = -0.01 +- 0.04, [C/H] = -0.17 +- 0.08, [N/H] = 0.53 +- 0.07, [O/H] = 0.12 +- 0.09, and C/N = 0.79 +- 0.08. We obtain from 27 giants in NGC 6705: [Fe/H] = 0.0 +- 0.05, [C/H] = -0.08 +- 0.06, [N/H] = 0.61 +- 0.07, [O/H] = 0.13 +- 0.05, and C/N = 0.83 +- 0.19. The C/N ratios of stars in the investigated open clusters were compared with the ratios predicted by stellar evolutionary models. For the corresponding stellar turn-off masses from 1.9 to 3.3 solar masses, the observed C/N ratio values are very close to the predictions of standard first dredge-up models as well as to models of thermohaline extra-mixing. The average [O/H] abundance ratios of NGC 4815 and NGC 6705 are compared with the predictions of two Galactic chemical evolution models. The data are consistent with the evolution at the solar radius within the errors.

(Abridged) AVOCADO aims at establishing firm conclusions on the formation and evolution of dwarf galaxies by constructing and analysing a homogeneous, multiwavelength dataset for a statistically significant sample of approximately 6500 nearby dwarfs (M_i > -18 mag). We present the sample selection criteria and describe the suite of analysis tools, some of them developed in the framework of the Virtual Observatory. We make use of optical spectra and UV-to-NIR imaging of the dwarf sample to derive SFRs, stellar masses, ages and metallicities - which are further supplemented with structural parameters that are used to classify them morphologically. This dataset, coupled with a detailed characterisation of each dwarf's environment, allows for a fully comprehensive investigation of their origins and to track the potential evolutionary paths between the different dwarf types. We characterise the local environment of all dwarfs in our sample, paying special attention to trends with current star formation activity. We find that virtually all quiescent dwarfs are located in the vicinity (projected distances < 1.5 Mpc) of L* companions, consistent with recent results. While star-forming dwarfs are preferentially found at separations of the order of 1 Mpc, there appears to be a tail towards low separations (<100 kpc) in the distribution of projected distances. We speculate that, modulo projection effects, this probably represents a genuine population of late-type dwarfs caught upon first infall about their host and before environmental quenching has fully operated. In this context, these results suggest that internal mechanisms -such as gas exhaustion via star formation or feedback effects- are not sufficient to completely cease the star formation activity in dwarf galaxies, and that becoming the satellite of a massive central galaxy appears to be a necessary condition to create a quiescent dwarf.

(abridged) We want to provide further evidence of the origin of the proposed stellar members of Chamaeleon and to identify interlopers from the foreground \epsilon Cha and \eta Cha associations. To this aim, we compile lists of spectroscopically confirmed members of Chamaeleon I and II, \epsilon Cha and \eta Cha, and of background objects in the same line of sight. Using Virtual Observatory tools, we cross-match these lists with the UCAC3 catalogue to get the proper motions of the objects. In the vector point diagram, we identify the different moving groups, and use this information to study the membership of proposed candidate members of the associations from the literature. For those objects with available radial velocities, we compute their Galactic space velocities. We look for correlations between the known properties of the objects and their proper motions. The members of the dark clouds exhibit clearly different proper motions from those of the foreground associations and of the background stars. The data suggest that Chamaeleon II could have different dynamical properties from Chamaeleon I. Although the two foreground clusters \epsilon and \eta Chamaeleontis constitute two different proper motion groups, they have similar spatial motions, which are different from the spatial motion of Chamaeleon I. On the other hand, the space motions of the Chamaeleon II stars look more similar to those of the foreground clusters than to the Chamaeleon I stars, but the numbers are low. Hence, with the available data it is unclear to what extent the stellar populations in both clouds are physically connected to each other. We find no correlations between the proper motions and the properties of the objects in either of the clouds.

The observation of blank fields, defined as regions of the sky that are devoid of stars down to a given threshold magnitude, constitutes one of the most relevant calibration procedures required for the proper reduction of astronomical data obtained following typical observing strategies. In this work, we have used the Delaunay triangulation to search for deep blank fields throughout the whole sky, with a minimum size of 10 arcmin in diameter and an increasing threshold magnitude from 15 to 18 in the R band of the USNO-B Catalog of the United States Naval Observatory. The result is a catalogue with the deepest blank fields known so far. A short sample of these regions has been tested with the 10.4m Gran Telescopio Canarias, and it has been shown to be extremely useful for medium and large size telescopes. Because some of the regions found could also be suitable for new extragalactic studies, we have estimated the galactic extinction in the direction of each deep blank field. This catalogue is accessible through the Virtual Observatory tool TESELA, and the user can retrieve - and visualize using Aladin - the deep blank fields available near a given position in the sky.

Several clusters of red supergiants have been discovered in a small region of the Milky Way close to the base of the Scutum-Crux Arm and the tip of the Long Bar. Population synthesis models indicate that they must be very massive to harbour so many supergiants. Among them, Stephenson 2, with a core grouping of 26 RSGs, is a strong candidate to be the most massive cluster in the Galaxy. It is located close to a region where a strong over-density of RSGs had been found. We explore the actual cluster size and its possible connection to this over-density. We have performed a cross-match between DENIS, USNO-B1 and 2MASS to identify candidate obscured luminous red stars around Ste 2, and in a control nearby region, finding >600 candidates. More than 400 are sufficiently bright in I to allow observation with a 4-m class telescope. We have observed a subsample of ~250 stars, using AF2 on the WHT telescope in La Palma, obtaining intermediate-resolution spectroscopy in the 7500--9000A range. We derive spectral types and luminosity classes for all these objects and measure their radial velocities. Our targets are G and K supergiants, late (>=M4) M giants, and M-type bright giants (luminosity class II) and supergiants. We find ~35 RSGs with radial velocities similar to Ste 2 members, spread over the two areas surveyed. In addition, we find ~40 RSGs with radial velocities incompatible in principle with a physical association. Our results show that Ste 2 is not an isolated cluster, but part of a huge structure likely containing hundreds of RSGs, with radial velocities compatible with the terminal velocity at this Galactic longitude (and a distance ~6kpc). In addition, we find evidence of several populations of massive stars at different distances along this line of sight [ABRIDGED].

Aims: With available Virtual Observatory tools, we looked for new M dwarfs in the solar neighbourhood and M giants with high tangential velocities. Methods: From an all-sky cross-match between the optical Tycho-2 and the near-infrared 2MASS catalogues, we selected objects with proper motions >50mas/yr and very red V-Ks colours. For the most interesting targets, we collected multi-wavelength photometry, constructed spectral energy distributions, estimated effective temperatures and surface gravities from fits to atmospheric models, performed time-series analysis of ASAS V-band light curves, and assigned spectral types from low-resolution spectroscopy obtained with CAFOS at the 2.2m Calar Alto telescope. Results: We got a sample of 59 bright red high proper-motion objects, including fifty red giants, four red dwarfs, and five objects reported in this work for the first time. The five new stars have magnitudes V~10.8-11.3mag, reduced proper motions midway between known dwarfs and giants, near-infrared colours typical of giants, and effective temperatures Teff~2900-3400K. From our time-series analysis, we discovered a long secondary period in Ruber 4 and an extremely long primary period in Ruber 6. With the CAFOS spectra, we confirmed the red giant nature of Ruber 7 and 8, the last of which seems to be one of the brightest metal-poor M giants ever identified.

The observation of blank fields, regions of the sky devoid of stars down to a given threshold magnitude, constitutes one of the typical important calibration procedures required for the proper reduction of astronomical data obtained in imaging mode. This work describes a method, based on the use of the Delaunay triangulation on the surface of a sphere, that allows the easy generation of blank fields catalogues. In addition to that, a new tool named TESELA, accessible through the WEB, has been created to facilitate the user to retrieve, and visualise using the VO-tool Aladin, the blank fields available near a given position in the sky.

With available Virtual Observatory tools, we looked for new bright blue high proper motion objects in the entire sky: white dwarfs, hot subdwarfs, runaway OB stars, and early-type stars in nearby young moving groups. We performed an all-sky cross-match between the optical Tycho-2 and near-infrared 2MASS catalogues with Aladin, and selected objects with proper motions >50mas/yr and colours Vt-Ks<-0.5mag with TOPCAT. We also collected multi-wavelength photometry, constructed the spectral energy distributions and estimated effective temperatures from fits to atmospheric models with VOSA for the most interesting targets. We assembled a sample of 32 bright blue high proper motion objects, including ten sdO/B subdwarfs, nine DA white dwarfs, five young early-type stars (two of which are runaway stars), two blue horizontal branch stars, one star with poor information, and five objects reported for the first time in this work. These last five objects have magnitudes Bt~11.0-11.6mag, effective temperatures ~24,000-30,000K, and are located in the region of known white dwarfs and hot subdwarfs in a reduced proper motion-colour diagram. We confirmed the hot subdwarf nature of one of the new objects, Albus 5, with public far-ultraviolet spectroscopic data obtained with FUSE.

Aims: By using kinematical information, we test the membership of the new Lupus candidate members proposed by the Cores to Disks (c2d) Spitzer Legacy Program program and by a complementary optical survey. We also investigate the relationship between the proper motions (pm) and other properties, in order to get some clues about their formation and early evolution. Methods: We compiled a list of members and possible members of Lupus 1, 3, and 4, together with all information on their spectral types, disks, and physical parameters. Using VO-tools, we cross-matched this list with the astrometric catalogues to get pm. Our final sample contains sources with magnitudes I<16mag and estimated masses >~0.1Msun. Results: According to the kinematic information, our sources can be divided into two main groups. The first one contains sources with higher pm in agreement with other Gould Belt populations and with spatial distribution, optical and near-infrared colours, and disk composition consistent with these objects belonging to the Lupus clouds. In the second group, sources have lower pm with random orientations, and they are mostly located outside the cloud cores, making their association with the Lupus complex more doubtful. We investigate the properties of the higher pm group, but cannot find any correlations with spatial location, binarity, the presence of a circumstellar disk, or with physical properties such as effective temperature, luminosity, mass, or age. Conclusions: We conclude that the lower pm group probably represents a background population or mixture of populations unrelated to the Lupus clouds. The higher pm group, on the other hand, has properties consistent with it being a genuine population of the Lupus star-forming region. More accurate pm and/or radial velocity information are required for a more detailed study of the kinematic properties of the Lupus stellar members.

Aiming to perform a study of the warm dust and gas in the luminous blue variable star G79.29+0.46 and its associated nebula, we present infrared Spitzer imaging and spectroscopy, and new CO J=2-->1 and 4-->3 maps obtained with the IRAM 30m radio telescope and with the Submillimeter Telescope, respectively. We have analyzed the nebula detecting multiple shells of dust and gas connected to the star. Using Infrared Spectrograph-Spitzer spectra, we have compared the properties of the central object, the nebula, and their surroundings. These spectra show a rich variety of solid-state features (amorphous silicates, polycyclic aromatic hydrocarbons, and CO2 ices) and narrow emission lines, superimposed on a thermal continuum. We have also analyzed the physical conditions of the nebula, which point to the existence of a photo-dissociation region.

High spatial resolution images of PNe have shown their extremely complex morphology. However, the circumstellar envelopes of their progenitors, the AGB stars, are strikingly spherical. In order to understand the carving processes leading to axisymmetric nebulae, we are carrying out a study of a large sample of pre-PNe. Our emission model of the nebular molecular gas (12CO & 13CO) will allow us to determine important physical parameters (mass, linear momentum, kinetic energy) of the fast bipolar and slow spherical nebular components separately. We will study in an innovative way the properties for each source individually, and put our results in an evolutionary context with the help of the data obtained by us and collected from the literature.

We performed a near-IR/optical monitoring programme from 1999 to 2005 in order to study the variability properties of the 'Arecibo sample of OH/IR stars' (periods, amplitudes, and colour variations). Here we describe this multi-wavelength long-term monitoring programme. Data analysis is still in process. Our ultimate goal is to study in particular the oxygen-rich AGB stars with masses larger than 2 Solar Mass, which are probably rare among AGB stars discovered optically and/or in the near-IR, but are common in samples discovered in the mid- to far-IR.

We report the detection of nebular emission lines in the optical and mid-infrared spectra of IRAS 17347-3139, a heavily obscured OH/IR star which may be rapidly evolving from the AGB to the PN stage. The presence of emission lines is interpreted as a clear indication that the ionization of its circumstellar envelope has already started. This source belongs to the rare class of objects known as `OHPNe' displaying both OH maser and radio continuum emission. However, unlike the rest of stars in this class, prominent C-rich dust features are detected in its mid-infrared spectrum, which makes the analysis of this star particularly interesting.

We present the results of a near-infrared monitoring programme of a selected sample of stars, initially suspected to be Mira variables and OH/IR stars, covering more than a decade of observations. The objects monitored cover the typical range of IRAS colours shown by O-rich stars on the Asymptotic Giant Branch and show a surprisingly large diversity of variability properties. 16 objects are confirmed as large-amplitude variables. Periods between 360 and 1800 days and typical amplitudes from 1 to 2 magnitudes could be determined for nine of them. In three light curves we find a systematic decrease of the mean brightness, two light curves show pronounced asymmetry. One source, IRAS 07222-2005, shows infrared colours typical of Mira variables but pulsates with a much longer period (approx. 1200 days) than a normal Mira. Two objects are ither close to (IRAS 03293+6010) or probably in (IRAS 18299-1705) the post-AGB phase. In IRAS 16029-3041 we found a systematic increase of the H-K colour of approximately 1 magnitude, which we interpret as evidence of a recent episode of enhanced mass loss. IRAS 18576+0341, a heavily obscured Luminous Blue Variable was also monitored. The star showed a continued decrease of brightness over a period of 7 years (1995 - 2002).

We present optical and near-infrared finding charts taken from the DSS and 2MASS surveys of 94 IRAS sources selected from the GLMP catalogue (Garc\'{\i}a-Lario 1992), and accurate astrometry (~0.2") for most of them. Selection criteria were very red IRAS colours representative for OH/IR stars with optically thick circumstellar shells and the presence of variability according to the IRAS variability index (VAR>50). The main photometric properties of the stars in this `GLMP sample' are presented, discussed and compared with the correspondent properties of the `Arecibo sample' of OH/IR stars studied in Jim\'enez-Esteban et al. (2005a). We find that 37% of the sample (N=34) have no counterpart in the 2MASS, implying extremely high optical depths of their shells. Most of the sources identified in the 2MASS are faint (K>~8) and of very red colour in the near-infrared, as expected. The brightest 2MASS counterpart (K=5.3mag) was found for IRAS 18299--1705. Its blue colour H--K=1.3 suggests that IRAS 8299--1705 is a post-AGB star. A couple of GLMP sources have faint but relatively blue counterparts. They might be misidentifed field stars or stars hich experienced recently a drop of their mass loss rates. The `GLMP sample' in general is made of oxygen-rich AGB stars, which are highly obscured by their circumstellar shells. They belong to the same population as the reddest OH/IR stars in the `Arecibo sample'.

In this paper we present optical and near-infrared finding charts, accurate astrometry (~1") and single-epoch near-infrared photometry for 371 IRAS sources, 96% of those included in the so-called Arecibo sample of OH/IR stars (Eder et al. 1988; Lewis et al. 1990a; Chengalur et al. 1993). The main photometric properties of the stars in the sample are presented and discussed as well as the problems found during the process of identification of the optical/near-infrared counterparts. In addition, we also identify suitable reference stars in each field to be used for differential photometry purposes in the future. We find that 39% of the sources (144 in number) have no optical counterpart, 8 of them being invisible even at near infrared wavelengths. The relative distribution of sources with and without optical counterpart in the IRAS two-colour diagram and their characteristic near infrared colours are interpreted as the consequence of the increasing thickness of their circumstellar shells. Among the objects not detected at near infrared wavelengths four non-variable sources are proposed to be heavily obscured post-AGB stars which have just very recently left the AGB. Eight additional objects with unusually bright and/or blue near-infrared colours are identified as candidate post-AGB stars and/or proto-planetary nebulae.

The `O-rich AGB sequence' is a sequence of colours describing the location of O-rich AGB stars in the IRAS two-colour diagram [12]--[25] vs [25]--[60]. We propose an evolutionary scenario for this sequence in which all stars, independent of their progenitor mass, start the AGB phase in the blue part of the `O-rich AGB sequence' and then evolve toward redder colors, although only the more massive stars would reach the very end of the `O-rich AGB sequence'. The sources located in the blue part of the sequence are mainly Mira variables, whose mean period is increasing with the IRAS colours. Most of them will evolve into O-rich Type II (and III) Planetary Nebulae. Part of the stars located in the red part of the sequence will change their chemical composition from O-rich to C-rich during their evolution in the AGB phase, and might evolve into C-rich Type II Planetary Nebulae. Hot bottom burning may prevent the conversion to carbon stars of the rest of sources located in the red part of the sequence and they will end up as N-rich Type I Planetary Nebulae.