Search

We used the Javalambre Photometric Local Universe Survey (J-PLUS) DR2 photometry in twelve optical bands over 2176 deg2 to estimate the fraction of white dwarfs with presence of CaII H+K absorption along the cooling sequence. We compared the J-PLUS photometry against metal-free theoretical models to estimate the equivalent width in the J0395 passband of 10 nm centered at 395 nm (EW_J0395), a proxy to detect calcium absorption. A total of 4399 white dwarfs within 30000 > Teff > 5500 K and mass M > 0.45 Msun were analyzed. Their EW_J0395 distribution was modeled using two populations, corresponding to polluted and non-polluted systems, to estimate the fraction of calcium white dwarfs (f_Ca) as a function of Teff. The probability for each individual white dwarf of presenting calcium absorption, pca, was also computed. The comparison with both the measured Ca/He abundance and the metal pollution from spectroscopy shows that EW_J0395 correlates with the presence of calcium. The fraction of calcium white dwarfs increases from f_Ca = 0 at Teff = 13500 K to f_Ca = 0.15 at Teff = 5500 K. We compare our results with the fractions derived from the 40 pc spectroscopic sample and from SDSS spectra. The trend found in J-PLUS observations is also present in the 40 pc sample, however SDSS shows a deficit of metal-polluted objects at Teff < 12000 K. Finally, we found 39 white dwarfs with pca > 0.99. Twenty of them have spectra presented in previous studies, whereas we observed six additional targets. These 26 objects were all confirmed as metal-polluted systems. The J-PLUS optical data provide a robust statistical measurement for the presence of CaII H+K absorption in white dwarfs. We find a 15 +- 3 % increase in the fraction of calcium white dwarfs from Teff = 13500 K to 5500 K, which reflects their selection function in the optical from the total population of metal-polluted systems., Comment: Submitted to Astronomy and Astrophysics. Comments are welcome

With its 12 optical filters, the Javalambre-Photometric Local Universe Survey (J-PLUS) provides an unprecedented multicolor view of the local Universe. The third data release (DR3) covers 3,192 deg$^2$ and contains 47.4 million objects. However, the classification algorithms currently implemented in its pipeline are deterministic and based solely on the sources morphology. Our goal is classify the sources identified in the J-PLUS DR3 images into stars, quasi-stellar objects (QSOs), and galaxies. For this task, we present BANNJOS, a machine learning pipeline that uses Bayesian neural networks to provide the probability distribution function (PDF) of the classification. BANNJOS is trained on photometric, astrometric, and morphological data from J-PLUS DR3, Gaia DR3, and CatWISE2020, using over 1.2 million objects with spectroscopic classification from SDSS DR18, LAMOST DR9, DESI EDR, and Gaia DR3. Results are validated using $1.4 10^5$ objects and cross-checked against theoretical model predictions. BANNJOS outperforms all previous classifiers in terms of accuracy, precision, and completeness across the entire magnitude range. It delivers over 95% accuracy for objects brighter than $r = 21.5$ mag, and ~90% accuracy for those up to $r = 22$ mag, where J-PLUS completeness is < 25%. BANNJOS is also the first object classifier to provide the full probability distribution function (PDF) of the classification, enabling precise object selection for high purity or completeness, and for identifying objects with complex features, like active galactic nuclei with resolved host galaxies. BANNJOS has effectively classified J-PLUS sources into around 20 million galaxies, 1 million QSOs, and 26 million stars, with full PDFs for each, which allow for later refinement of the sample. The upcoming J-PAS survey, with its 56 color bands, will further enhance BANNJOS's ability to detail each source's nature., Comment: Submitted to Astronomy and Astrophysics. 29 pages, 23 figures, 4 tables

The Gaia-ESO Survey is an European Southern Observatory (ESO) public spectroscopic survey that targeted $10^5$ stars in the Milky Way covering the major populations of the disk, bulge and halo. The observations were made using FLAMES on the VLT obtaining both UVES high ($R\sim47,000$) and GIRAFFE medium ($R\sim20,000$) resolution spectra. The analysis of the Gaia-ESO spectra was the work of multiple analysis teams (nodes) within five working groups (WG). The homogenisation of the stellar parameters within WG11 (high resolution observations of FGK stars) and the homogenisation of the stellar parameters within WG10 (medium resolution observations of FGK stars) is described here. In both cases, the homogenisation was carried out using a bayesian Inference method developed specifically for the Gaia-ESO Survey by WG11. The WG10 homogenisation primarily used the cross-match of stars with WG11 as the reference set in both the stellar parameter and chemical abundance homogenisation. In this way the WG10 homogenised results have been placed directly onto the WG11 stellar parameter and chemical abundance scales. The reference set for the metal-poor end was sparse which limited the effectiveness of the homogenisation in that regime. For WG11, the total number of stars for which stellar parameters were derived was 6,231 with typical uncertainties for Teff, log g and [Fe/H] of 32~K, 0.05 and 0.05 respectively. One or more chemical abundances out of a possible 39 elements were derived for 6,188 of the stars. For WG10, the total number of stars for which stellar parameters were derived was 76,675 with typical uncertainties for Teff, log g and [Fe/H] of 64~K, 0.15 and 0.07 respectively. One or more chemical abundances out of a possible 30 elements were derived for 64,177 of the stars., Comment: 31 pages, 19 figures

The Gaia-ESO Survey (GES) is a large public spectroscopic survey which acquired spectra for more than 100000 stars across all major components of the Milky Way. In addition to atmospheric parameters and stellar abundances that have been derived in previous papers of this series, the GES spectra allow us to detect spectroscopic binaries with one (SB1), two (SB2) or more (SBn $\ge$ 3) components. Cross-correlation functions (CCFs) have been re-computed thanks to a dozen spectral masks probing a range of effective temperatures, surface gravities and metallicities. By optimising the mask choice for a given spectrum, the new computed so-called Nacre (Narrow cross-correlation experiment) CCFs are narrower and allow to unblend more stellar components than standard masks. The Doe (Detection of Extrema) extremum-finding code then selects the individual components and provides their radial velocities. From the sample of HR10 and HR21 spectra corresponding to 37565 objects, the present study leads to the detection of 322 SB2, ten (three of them being tentative) SB3, and two tentative SB4. In particular, compared to our previous study, the Nacre CCFs allow us to multiply the number of SB2 candidates by $\approx$ 1.5. The colour-magnitude diagram reveals, as expected, the shifted location of the SB2 main sequence. A comparison between the SB identified in Gaia DR3 and the ones detected in the present work is performed and the complementarity of the two censuses is discussed. An application to mass-ratio determination is presented, and the mass-ratio distribution of the GES SB2 is discussed. When accounting for the SB2 detection rate, an SB2 frequency of $\approx$ 1.4% is derived within the present stellar sample of mainly FGK-type stars. As primary outliers identified within the GES data, SBn spectra produce a wealth of information and useful constraints for the binary population synthesis studies.

The characterization of white dwarf atmospheres is crucial for accurately deriving stellar parameters such as effective temperature, mass, and age. We aim to classify the population of white dwarfs up to 500 pc into hydrogen-rich or hydrogen-deficient atmospheres based on Gaia spectra and to derive an accurate spectral type-temperature distribution of white dwarfs as a function of the effective temperature for the largest observed unbiased sample of these objects. We took advantage of the recent Gaia low-resolution spectra available for 76,657 white dwarfs up to 500 pc. We calculated synthetic J-PAS narrow-band photometry and fitted the spectral energy distribution of each object with up-to-date models for hydrogen-rich and helium-rich white dwarf atmospheres. We estimated the probability for a white dwarf to have a hydrogen-rich atmosphere and validated the results using the Montreal White Dwarf Database. Finally, precise effective temperature values were derived for each object using La Plata evolutionary models. We have successfully classified a total of 65,310 white into DAs and non-DAs with an accuracy of 94%. An unbiased subsample of nearly 34,000 objects was built, from which we computed a precise spectral distribution spanning an effective temperature range from 5,500 to 40,000 K, while accounting for potential selection effects. Some characteristic features of the spectral evolution, such as the deficit of helium-rich stars at T_eff $\approx$35,000-40,000 K and in the range 22,000 < T_eff < 25,000 K, as well as a gradual increase from 18,000K to T_eff $\approx$7,000K, where the non-DA stars percentage reaches its maximum of 41%, followed by a decrease for cooler temperatures, are statistically significant. These findings will provide precise constraints for the proposed models of spectral evolution., Comment: 7 pages, 11 pages. Accepted for publication in A&A

The Gaia-ESO Survey is a public spectroscopic survey that has targeted $\gtrsim10^5$ stars covering all major components of the Milky Way from the end of 2011 to 2018, delivering its public final release in May 2022. Unlike other spectroscopic surveys, Gaia-ESO is the only survey that observed stars across all spectral types with dedicated, specialised analyses: from O ($T_\mathrm{eff} \sim 30,000-52,000$~K) all the way to K-M ($\gtrsim$3,500~K). The physics throughout these stellar regimes varies significantly, which has previously prohibited any detailed comparisons between stars of significantly different type. In the final data release (internal data release 6) of the Gaia-ESO Survey, we provide the final database containing a large number of products such as radial velocities, stellar parameters and elemental abundances, rotational velocity, and also, e.g., activity and accretion indicators in young stars and membership probability in star clusters for more than 114,000 stars. The spectral analysis is coordinated by a number of Working Groups (WGs) within the Survey, which specialise in the various stellar samples. Common targets are analysed across WGs to allow for comparisons (and calibrations) amongst instrumental setups and spectral types. Here we describe the procedures employed to ensure all Survey results are placed on a common scale to arrive at a single set of recommended results for all Survey collaborators to use. We also present some general quality and consistency checks performed over all Survey results., Comment: A&A accepted, minor revision, 36 pages, 38 figures

We present the photometric calibration of the twelve optical passbands for the Javalambre Photometric Local Universe Survey (J-PLUS) third data release (DR3), comprising 1642 pointings of two square degrees each. We selected nearly 1.5 million main sequence stars with a signal-to-noise ratio larger than ten in the twelve J-PLUS passbands and available low-resolution (R = 20-80) spectrum from the blue and red photometers (BP/RP) in Gaia DR3. We compared the synthetic photometry from BP/RP spectra with the J-PLUS instrumental magnitudes, after correcting for the magnitude and color terms between both systems, to obtain an homogeneous photometric solution for J-PLUS. To circumvent the current limitations in the absolute calibration of the BP/RP spectra, the absolute color scale was derived using the locus of 109 white dwarfs closer than 100 pc with a negligible interstellar extinction. Finally, the absolute flux scale was anchored to the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) photometry in the r band. The precision of the J-PLUS photometric calibration, estimated from duplicated objects observed in adjacent pointings and by comparison with the spectro-photometric standard star GD 153, is ~12 mmag in u, J0378, and J0395; and ~7 mmag in J0410, J0430, g, J0515, r, J0660, i, J0861, and z. The estimated accuracy in the calibration along the surveyed area is better than 1% for all the passbands. The Gaia BP/RP spectra provide a high-quality, homogeneous photometric reference in the optical range across the full-sky, in spite of their current limitations as an absolute reference. The calibration method for J-PLUS DR3 reaches an absolute precision and accuracy of 1% in the twelve optical filters within an area of 3284 square degrees., Comment: Submitted to Astronomy and Astrophysics. 17 pages, 14 figures, 2 appendix. Comments are welcome. arXiv admin note: text overlap with arXiv:2101.12407

The third data release of Gaia has provided low resolution spectra for ~100,000 white dwarfs (WDs) that, together with the excellent photometry and astrometry, represent an unrivalled benchmark for the study of this population. In this work, we first built a highly-complete volume-limited sample consisting in 12,718 WDs within 100 pc from the Sun. The use of VOSA tool allowed us to perform an automated fitting of their spectral energy distributions to different atmospheric models. In particular, the use of spectrally derived J-PAS photometry from Gaia spectra led to the classification of DA and non-DA WDs with an accuracy >90%, tested in already spectroscopically labelled objects. The excellent performance achieved was extended to practically the whole population of WDs with effective temperatures above 5500 K. Our results show that, while the A branch of the Gaia WD Hertzsprung-Russell diagram is practically populated by DA WDs, the B branch is largely formed by non-DAs (65%). The remaining 35% of DAs within the B branch implies a second peak at ~0.8 Mo in the DA-mass distribution. Additionally, the Q branch and its extension to lower temperatures can be observed for both DA and non-DA objects due to core crystallisation. Finally, we derived a detailed spectral evolution function, which confirms a slow increase of the fraction of non-DAs as the effective temperature decreases down to 10,500 K, where it reaches a maximum of 36% and then decreases for lower temperatures down to ~31%., Comment: Accepted for publication in MNRAS

With a unique set of 54 overlapping narrow-band and two broader filters covering the entire optical range, the incoming Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will provide a great opportunity for stellar physics and near-field cosmology. In this work, we use the miniJPAS data in 56 J-PAS filters and 4 complementary SDSS-like filters to explore and prove the potential of the J-PAS filter system in characterizing stars and deriving their atmospheric parameters. We obtain estimates for the effective temperature with a good precision (<150 K) from spectral energy distribution fitting. We have constructed the metallicity-dependent stellar loci in 59 colours for the miniJPAS FGK dwarf stars, after correcting certain systematic errors in flat-fielding. The very blue colours, including uJAVA-r, J0378-r, J0390-r, uJPAS-r, show the strongest metallicity dependence, around 0.25 mag/dex. The sensitivities decrease to about 0.1 mag/dex for the J0400-r, J0410-r, and J0420-r colours. The locus fitting residuals show peaks at the J0390, J0430, J0510, and J0520 filters, suggesting that individual elemental abundances such as [Ca/Fe], [C/Fe], and [Mg/Fe] can also be determined from the J-PAS photometry. Via stellar loci, we have achieved a typical metallicity precision of 0.1 dex. The miniJPAS filters also demonstrate strong potential in discriminating dwarfs and giants, particularly the J0520 and J0510 filters. Our results demonstrate the power of the J-PAS filter system in stellar parameter determinations and the huge potential of the coming J-PAS survey in stellar and Galactic studies., Comment: 16 pages, 16 figures, 3 tables, accepted by MNRAS

Ultracool dwarfs (UCDs) comprise the lowest mass members of the stellar population and brown dwarfs, from M7 V to cooler objects with L, T, and Y spectral types. Most of them have been discovered using wide-field imaging surveys, for which the Virtual Observatory (VO) has proven to be of great utility. We aim to perform a search for UCDs in the entire Javalambre Photometric Local Universe Survey (J-PLUS) second data release (2176 deg$^2$) following a VO methodology. We also explore the ability to reproduce this search with a purely machine learning (ML)-based methodology that relies solely on J-PLUS photometry. We followed three different approaches based on parallaxes, proper motions, and colours, respectively, using the VOSA tool to estimate the effective temperatures. For the ML methodology, we built a two-step method based on principal component analysis and support vector machine algorithms. We identified a total of 7827 new candidate UCDs, which represents an increase of about 135% in the number of UCDs reported in the sky coverage of the J-PLUS second data release. Among the candidate UCDs, we found 122 possible unresolved binary systems, 78 wide multiple systems, and 48 objects with a high Bayesian probability of belonging to a young association. We also identified four objects with strong excess in the filter corresponding to the Ca II H and K emission lines and four other objects with excess emission in the H$\alpha$ filter. With the ML approach, we obtained a recall score of 92% and 91% in the test and blind test, respectively. We consolidated the proposed search methodology for UCDs, which will be used in deeper and larger upcoming surveys such as J-PAS and Euclid. We concluded that the ML methodology is more efficient in the sense that it allows for a larger number of true negatives to be discarded prior to analysis with VOSA, although it is more photometrically restrictive., Comment: Accepted in A&A

With its origin coming from several sources (Big Bang, stars, cosmic rays) and given its strong depletion during its stellar lifetime, the lithium element is of great interest as its chemical evolution in the Milky Way is not well understood at present. To help constrain stellar and galactic chemical evolution models, numerous and precise lithium abundances are necessary for a large range of evolutionary stages, metallicities, and Galactic volume. In the age of stellar parametrization on industrial scales, spectroscopic surveys such as APOGEE, GALAH, RAVE, and LAMOST have used data-driven methods to rapidly and precisely infer stellar labels (atmospheric parameters and abundances). To prepare the ground for future spectroscopic surveys such as 4MOST and WEAVE, we aim to apply machine learning techniques to lithium measurements and analyses. We trained a convolution neural network (CNN), coupling Gaia-ESO Survey iDR6 stellar labels (Teff, log(g), [Fe/H], and A(Li)) and GIRAFFE HR15N spectra, to infer the atm parameters and lithium abundances for ~40,000 stars. We show that the CNN properly learns the physics of the stellar labels, from relevant spectral features through a broad range of evolutionary stages and stellar parameters. The Li feature at 6707.8 A is successfully singled out by our CNN, among the thousands of lines. Rare objects such as Li-rich giants are found in our sample. This level of performance is achieved thanks to a meticulously built, high-quality, and homogeneous training sample. The CNN approach is very well adapted for the next generations of spectroscopic surveys aimed at studying (among other elements) lithium, such as the 4MIDABLE-LR/HR (4MOST Milky Way disk and bulge low- and high-resolution) surveys. In this context, the caveats of ML applications should be appropriately investigated, along with the realistic label uncertainties and upper limits for abundances., Comment: A&A, accepted 05 December 2022, 23 pages, 23 figures. Codes and trained models available at https://github.com/SamirNepal/Li\_CNN\_2022

We aim to prepare the machine-learning ground for the next generation of spectroscopic surveys, such as 4MOST and WEAVE. Our goal is to show that convolutional neural networks can predict accurate stellar labels from relevant spectral features in a physically meaningful way. We built a neural network and trained it on GIRAFFE spectra with associated stellar labels from the sixth internal Gaia-ESO data release. Our neural network predicts the atmospheric parameters Teff and log(g) as well as the chemical abundances [Mg/Fe], [Al/Fe], and [Fe/H] for 30115 stellar spectra. The scatter of predictions from eight slightly different network models shows a high internal precision of the network results: 24 K for Teff, 0.03 for log(g), 0.02 dex for [Mg/Fe], 0.03 dex for [Al/Fe], and 0.02 dex for [Fe/H]. The network gradients reveal that the network is inferring the labels in a physically meaningful way from spectral features. Validation with benchmark stars and several scientific applications confirm that our network predictions are accurate for individual stars and recover the properties of different stellar populations in the Milky Way galaxy. Such a study provides very good insights into the application of machine-learning for the spectral analysis of large-scale spectroscopic surveys, such as WEAVE and 4MIDABLE-LR and -HR (4MOST Milky Way disk and bulge low- and high-resolution). The community will have to put a substantial effort into building proactive training sets for machine-learning methods to minimize the possible systematics., Comment: Submitted to A&A

The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper (arXiv:2206.02901) introduces the survey results. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. The Gaia-ESO Survey obtained 202,000 spectra of 115,000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022., Comment: 38 pages. A&A in press

A renewed interest about the origin of \emph{r}-process elements has been stimulated by the multi-messenger observation of the gravitational event GW170817, with the detection of both gravitational waves and electromagnetic waves corresponding to the merger of two neutron stars. Such phenomenon has been proposed as one of the main sources of the \emph{r}-process. However, the origin of the \emph{r}-process elements at different metallicities is still under debate. We aim at investigating the origin of the \emph{r}-process elements in the Galactic thin disc population. From the sixth internal data release of the \emph{Gaia}-ESO we have collected a large sample of Milky Way thin- and thick-disc stars for which abundances of Eu, O, and Mg are available. The sample consists of members of 62 open clusters, located at a Galactocentric radius from $\sim 5$ kpc to $\sim 20$ kpc in the disc, in the metallicity range $[-0.5, 0.4]$ and covering an age interval from 0.1 to 7 Gy, and about 1300 Milky Way disc field stars in the metallicity range $[-1.5, 0.5]$. We compare the observations with the results of a chemical evolution model, in which we varied the nucleosynthesis sources for the three considered elements. Our main result is that Eu in the thin disc is predominantly produced by sources with short lifetimes, such as magneto-rotationally driven SNe. There is no strong evidence for additional sources at delayed times. Our findings do not imply that there cannot be a contribution from mergers of neutron stars in other environments, as in the halo or in dwarf spheroidal galaxies, but such a contribution is not needed to explain Eu abundances at thin disc metallicities., Comment: 20 pages, accepted for publication in A&A

In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars using FLAMES on the ESO VLT (both Giraffe and UVES spectrographs), covering all the Milky Way populations, with a special focus on open star clusters. This article provides an overview of the survey implementation (observations, data quality, analysis and its success, data products, and releases), of the open cluster survey, of the science results and potential, and of the survey legacy. A companion article (Gilmore et al.) reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. The GES has determined homogeneous good-quality radial velocities and stellar parameters for a large fraction of its more than 110,000 unique target stars. Elemental abundances were derived for up to 31 elements for targets observed with UVES. Lithium abundances are delivered for about 1/3 of the sample. The analysis and homogenisation strategies have proven to be successful; several science topics have been addressed by the Gaia-ESO consortium and the community, with many highlight results achieved. The final catalogue has been released through the ESO archive at the end of May 2022, including the complete set of advanced data products. In addition to these results, the Gaia-ESO Survey will leave a very important legacy, for several aspects and for many years to come., Comment: Accepted for publication in Astronomy and Astrophysics. 30 pages, 30 figures, 4 tables

Context. Stellar parameters are among the most important characteristics in studies of stars, which are based on atmosphere models in traditional methods. However, time cost and brightness limits restrain the efficiency of spectral observations. The J-PLUS is an observational campaign that aims to obtain photometry in 12 bands. Owing to its characteristics, J-PLUS data have become a valuable resource for studies of stars. Machine learning provides powerful tools to efficiently analyse large data sets, such as the one from J-PLUS, and enable us to expand the research domain to stellar parameters. Aims. The main goal of this study is to construct a SVR algorithm to estimate stellar parameters of the stars in the first data release of the J-PLUS observational campaign. Methods. The training data for the parameters regressions is featured with 12-waveband photometry from J-PLUS, and is cross-identified with spectrum-based catalogs. These catalogs are from the LAMOST, the APOGEE, and the SEGUE. We then label them with the stellar effective temperature, the surface gravity and the metallicity. Ten percent of the sample is held out to apply a blind test. We develop a new method, a multi-model approach in order to fully take into account the uncertainties of both the magnitudes and stellar parameters. The method utilizes more than two hundred models to apply the uncertainty analysis. Results. We present a catalog of 2,493,424 stars with the Root Mean Square Error of 160K in the effective temperature regression, 0.35 in the surface gravity regression and 0.25 in the metallicity regression. We also discuss the advantages of this multi-model approach and compare it to other machine-learning methods., Comment: 18 pages, 10 figures, accepted by AA

The chemical makeup of a star provides the fossil information of the environment where it formed. Under this premise, it should be possible to use chemical abundances to tag stars that formed within the same stellar association. This idea - known as chemical tagging - has not produced the expected results, especially within the thin disk where open stellar clusters have chemical patterns that are difficult to disentangle. The ultimate goal of this study is to probe the feasibility of chemical tagging within the thin disk population using high-quality data from a controlled sample of stars. We also aim at improving the existing techniques of chemical tagging and giving guidance on different strategies of clustering analysis in the elemental abundance space. Here we develop the first blind search of open clusters' members through clustering analysis in the elemental abundance space using the OPTICS algorithm applied to data from the Gaia-ESO survey. First, we evaluate different strategies of analysis, determining which ones are more performing. Second, we apply these methods to a data set including both field stars and open clusters attempting a blind recover of as many open clusters as possible. We show how specific strategies of data analysis can improve the final results. Specifically, we demonstrate that open clusters can be more efficaciously recovered with the Manhattan metric and on a space whose dimensions are carefully selected. Using these (and other) prescriptions we are able to recover open clusters hidden in our data set and find new members of these stellar associations. Our results indicate that there are chances of recovering open clusters' members via clustering analysis in the elemental abundance space. Presumably, the performances of chemical tagging will further increase with higher quality data and more sophisticated clustering algorithms., Comment: Accepted for publication by A&A. 21 pages, 11 figures

We analyze the white dwarf population in miniJPAS, the first square degree observed with 56 medium-band, 145 A in width optical filters by the Javalambre Physics of the accelerating Universe Astrophysical Survey (J-PAS), to provide a data-based forecast for the white dwarf science with low-resolution (R ~ 50) photo-spectra. We define the sample of the bluest point-like sources in miniJPAS with r < 21.5 mag, point-like probability larger than 0.5, (u-r) < 0.80 mag, and (g-i) < 0.25 mag. This sample comprises 33 sources with spectroscopic information, 11 white dwarfs and 22 QSOs. We estimate the effective temperature (Teff), the surface gravity, and the composition of the white dwarf population by a Bayesian fitting to the observed photo-spectra. The miniJPAS data permit the classification of the observed white dwarfs into H-dominated and He-dominated with 99% confidence, and the detection of calcium absorption and polluting metals down to r ~ 21.5 mag at least for sources with 7000 < Teff < 22000 K, the temperature range covered by the white dwarfs in miniJPAS. The effective temperature is estimated with a 2% uncertainty, close to the 1% from spectroscopy. A precise estimation of the surface gravity depends on the available parallax information. In addition, the white dwarf population at Teff > 7000 K can be segregated from the bluest extragalactic QSOs, providing a clean sample based on optical photometry alone. The J-PAS low-resolution photo-spectra provide precise and accurate effective temperatures and atmospheric compositions for white dwarfs, complementing the data from Gaia. J-PAS will also detect and characterize new white dwarfs beyond the Gaia magnitude limit, providing faint candidates for spectroscopic follow up., Comment: Submitted to Astronomy and Astrophyscis. 12 pages, 9 figures, 4 tables. Comments are welcome

The Gaia-ESO Survey (GES) is a large public spectroscopic survey that has collected, over a period of 6 years, spectra of ~ 10^5 stars. This survey provides not only the reduced spectra, but also the stellar parameters and abundances resulting from the analysis of the spectra. The GES dataflow is organised in 19 working groups. Working group 13 (WG13) is responsible for the spectral analysis of the hottest stars (O, B and A type, with a formal cut-off of Teff > 7000 K) that were observed as part of GES. We present the procedures and techniques that have been applied to the reduced spectra, in order to determine the stellar parameters and abundances of these stars. The procedure used is similar to that of other working groups in GES. A number of groups (called `Nodes') each independently analyse the spectra, using their state-of-the-art techniques and codes. Specific for the analysis in WG13 is the large temperature range that is covered (Teff = 7000 - 50,000 K), requiring the use of different analysis codes. Most Nodes can therefore only handle part of the data. Quality checks are applied to the results of these Nodes by comparing them to benchmark stars, and by comparing them one to another. For each star the Node values are then homogenised into a single result: the recommended parameters and abundances. Eight Nodes each analysed (part of) the data. In total 17,693 spectra of 6462 stars were analysed, most of them in 37 open star clusters. The homogenisation led to stellar parameters for 5584 stars. Abundances were determined for a more limited number of stars. Elements studied are He, C, N, O, Ne, Mg, Al, Si and Sc. Abundances for at least one of those elements were determined for 292 stars. The hot-star data analysed here, as well as the Gaia-ESO Survey data in general, will be of considerable use in future studies of stellar evolution and open clusters., Comment: 26 pages, 14 figures, accepted for publication in Astronomy & Astrophysics; language-edited version; two appendices merged

The last decade has seen a revolution in our knowledge of the Galaxy thanks to the Gaia and asteroseismic space missions and the ground-based spectroscopic surveys. To complete this picture, it is necessary to map the ages of its stellar populations. During recent years, the dependence on time of abundance ratios involving slow (s) neutron-capture and $\alpha$ elements (called chemical-clocks) has been used to provide estimates of stellar ages, usually in a limited volume close to the Sun. We aim to analyse the relations of chemical clocks in the Galactic disc extending the range to R$_{\rm GC}\sim$6-20~kpc. Using the sixth internal data release of the Gaia-ESO survey, we calibrated several relations between stellar ages and abundance ratios [s/$\alpha$] using a sample of open clusters, the largest one so far used with this aim. Thanks to their wide galactocentric coverage, we investigated the radial variations of the shape of these relations, confirming their non-universality. We estimated our accuracy and precision in recovering the global ages of open clusters, and the ages of their individual members. We applied the multi-variate relations with the highest correlation coefficients to the field star population. We confirm that there is no single age-chemical clock relationship valid for the whole disc, but that there is a dependence on the galactocentric position, which is related to the radial variation of the star formation history combined with the non-monotonic dependence on metallicity of the yields of the s-process elements from low- and intermediate-mass stars. Finally, the abundance ratios [Ba/$\alpha$] are more sensitive to age than those with [Y/$\alpha$] for young disc stars, and their slopes vary less with galactocentric distance., Comment: 14 pages, 10 figures + Appendix (3 tables and 2 figures)

The $Gaia$ mission has provided an unprecedented wealth of information about the white dwarf population of our Galaxy. In particular, our studies show that the sample up to 100\,pc from the Sun can be considered as practically complete. This fact allows us to estimate a precise fraction of double-degenerate ($1.18\pm 0.10$\%) and white dwarf plus main-sequence stars ($6.31\pm0.23$\%) among all white dwarfs through comoving pairs identification. With the aid of a detailed population synthesis code we are able to reproduce synthetic white dwarf populations with nearly identical fractions as the observed ones, thus obtaining valuable information about the binary fraction, $f_{\rm b}$, initial mass ratio distribution, $n(q)$, and initial separation distribution, $f(a)$, among other parameters. Our best-fit model is achieved within a $1\sigma$ confidence level for $f(a)\propto a^{-1}$, $n(q)\propto q^{n_q}$, with $n_q=-1.13^{+0.12}_{-0.10}$ and $f_{\rm b}=0.32\pm 0.02$. The fraction of white dwarf mergers generated by this model is $9\sim16\%$, depending on the common-envelope treatment. As sub-products of our modelling we find that around $1\sim3\%$ of the white dwarf population are unresolved double-degenerates and that only $\sim1\%$ of all white dwarfs contain a He-core. Finally, only a mild kick during white dwarf formation seems to be necessary for fitting the observed sky separation of double-degenerate systems., Comment: 14 pages, 8 figures, 5 tables, accepted for publication in MNRAS

The Gaia-ESO Survey (GES) is a public, high-resolution spectroscopic survey with FLAMES@VLT. GES targeted in particular a large sample of open clusters (OCs) of all ages. The different kinds of OCs are useful to reach the main science goals, which are the study of the OC structure and dynamics, the use of OCs to constrain and improve stellar evolution models, and the definition of Galactic disc properties (e.g. metallicity distribution). GES is organised in 19 working groups (WGs). We describe here the work of three of them, WG4 in charge of the selection of the targets within each cluster), WG1 responsible for defining the most probable candidate members, and WG6 in charge of the preparation of the observations. As GES has been conducted before Gaia DR2, we could not make use of the Gaia astrometry to define cluster members. We made use of public and private photometry to select the stars to be observed with FLAMES. Candidate target selection was based on ground-based proper motions, radial velocities, and X-ray properties when appropriate, and it was mostly used to define the position of the clusters' evolutionary sequences in the colour-magnitude diagrams. Targets for GIRAFFE were selected near the sequences in an unbiased way. We used available information on membership only for the few UVES stars. We collected spectra for 62 confirmed OCs (a few more were taken from the ESO archive). Among them are very young clusters, where the main targets are pre-main sequence stars, clusters with very hot and massive stars currently on the main sequence, intermediate-age and old clusters where evolved stars are the main targets. The selection of targets was as inclusive and unbiased as possible and we observed a representative fraction of all possible targets, thus collecting the largest, most accurate, and most homogeneous spectroscopic data set on ever achieved. [abridged], Comment: Accepted for publication on A&A. Abstract abridged from the accepted version

Pre-main sequence models with inflated radii are needed to simultaneously reproduce the colour-magnitude diagram and the lithium depletion pattern in young open clusters. We tested a new set of PMS models including radius inflation due to starspots or magnetic inhibition of convection, using five clusters observed by the Gaia-ESO Survey, spanning the age range ~10-100 Myr where such effects could be important. Gaia-ESO radial velocities were combined with Gaia EDR3 astrometry to obtain clean lists of high-probability members for the five clusters. A Bayesian maximum likelihood method was adopted to derive the best model parameters and the cluster reddening and age. Models were calculated for different values of the mixing length parameter ($\alpha_{ML}=2.0$, 1.5 and 1.0), without spots or with effective spot coverage $\beta_{spot}=0.2$ and 0.4. To reproduce the CMD and the Li depletion pattern in Gamma Vel A and B and in 25 Ori we need both a reduced convection efficiency $\alpha_{ML}=1.0$ and an effective spot coverage of ~20%. We obtained ages of 18 Myr and 21 Myr for Gamma Vel A and B, respectively, and 19 Myr for 25 Ori. However, a single isochrone is not sufficient to account for the Li dispersion, and an increasing level of spot coverage as mass decreases seems to be required. The older clusters (NGC2451B at 30 Myr, NGC2547 at 35 Myr, and NGC2516 at 138 Myr) are consistent with standard models ($\alpha_{ML}=2.0$ and no spots) except at low masses: a 20% spot coverage seems to reproduce the sequence of M-type stars better and might explain the observed abundance spread. The quality of Gaia-ESO data combined with Gaia allows us to gain important insights on PMS evolution. Models including starspots can provide a consistent explanation of the cluster sequences and Li abundances of young clusters, although a range of starspot coverage is required to fully reproduce the data., Comment: Accepted for publication in A&A, 18 pages, 10 figures. Revised version after language editing, revised acknowledgements

We estimated the spectral evolution of white dwarfs with effective temperature using the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2), that provides twelve photometric optical passbands over 2176 deg2. We analysed 5926 white dwarfs with r <= 19.5 mag in common between a white dwarf catalog defined from Gaia EDR3 and J-PLUS DR2. We performed a Bayesian analysis by comparing the observed J-PLUS photometry with theoretical models of hydrogen (H) and helium (He) dominated atmospheres. We estimated the PDF for effective temperature (Teff), surface gravity, parallax, and spectral type; and the probability of having a H-dominated atmosphere (pH) for each source. We applied a prior in parallax, using Gaia EDR3 measurements as reference, and derived a self-consistent prior for the atmospheric composition as a function of Teff. We described the fraction of He-dominated atmosphere white dwarfs (fHe) with a linear function of Teff at 5000 < Teff < 30000 K. We found fHe = 0.24 +- 0.01 at Teff = 10000 K, a change rate along the cooling sequence of 0.14 +- 0.02 per 10 kK, and a minimum He-dominated fraction of 0.08 +- 0.02 at the high-temperature end. We tested the obtained pH by comparison with spectroscopic classifications, finding that it is reliable. We estimated the mass distribution for the 351 sources with distance d < 100 pc, mass M > 0.45 Msun, and Teff > 6000 K. The result for H-dominated white dwarfs agrees with previous work, with a dominant M = 0.59 Msun peak and the presence of an excess at M ~ 0.8 Msun. This high-mass excess is absent in the He-dominated distribution, which presents a single peak. The J-PLUS optical data provides a reliable statistical classification of white dwarfs into H- and He-dominated atmospheres. We find a 21 +- 3 % increase in the fraction of He-dominated white dwarfs from Teff = 20000 K to Teff = 5000 K., Comment: Published in Astronomy and Astrophysics. 18 pages, 15 figures, 4 tables

Spectroscopy from the final internal data release of the Gaia-ESO Survey (GES) has been combined with Gaia EDR3 to assign membership probabilities to targets observed towards 63 Galactic open clusters and 7 globular clusters. The membership probabilities are based chiefly on maximum likelihood modelling of the 3D kinematics of the targets, separating them into cluster and field populations. From 43211 observed targets, 13985 are identified as highly probable cluster members ($P>0.9$), with an average membership probability of 0.993. The addition of GES radial velocities successfully drives down the fraction of false positives and we achieve better levels of discrimination in most clusters over the use of astrometric data alone, especially those at larger distances. Since the membership selection is almost purely kinematic, the union of this catalogue with GES and Gaia is ideal for investigating the photometric and chemical properties of clusters as a function of stellar mass, age and Galactic position., Comment: Accepted version for publication in MNRAS. 16 pages + 38 pages of Appendices

We explore the stellar content of the Javalambre Photometric Local Universe Survey (J-PLUS) Data Release 2 and show its potential to identify low-metallicity stars using the Stellar Parameters Estimation based on Ensemble Methods (SPEEM) pipeline. SPEEM is a tool to provide determinations of atmospheric parameters for stars and separate stellar sources from quasars, using the unique J-PLUS photometric system. The adoption of adequate selection criteria allows the identification of metal-poor star candidates suitable for spectroscopic follow-up. SPEEM consists of a series of machine learning models which uses a training sample observed by both J-PLUS and the SEGUE spectroscopic survey. The training sample has temperatures Teff between 4\,800 K and 9\,000 K; $\log g$ between 1.0 and 4.5, and $-3.1<[Fe/H]<+0.5$. The performance of the pipeline has been tested with a sample of stars observed by the LAMOST survey within the same parameter range. The average differences between the parameters of a sample of stars observed with SEGUE and J-PLUS, which were obtained with the SEGUE Stellar Parameter Pipeline and SPEEM, respectively, are $\Delta Teff\sim 41$ K, $\Delta \log g\sim 0.11$ dex, and $\Delta [Fe/H]\sim 0.09$ dex. A sample of 177 stars have been identified as new candidates with $[Fe/H]<-2.5$ and 11 of them have been observed with the ISIS spectrograph at the William Herschel Telescope. The spectroscopic analysis confirms that $64\%$ of stars have $[Fe/H]<-2.5$, including one new star with $[Fe/H]<-3.0$. SPEEM in combination with the J-PLUS filter system has shown the potential to estimate the stellar atmospheric parameters (Teff, $\log g$, and [Fe/H]). The spectroscopic validation of the candidates shows that SPEEM yields a success rate of $64\%$ on the identification of very metal-poor star candidates with $[Fe/H]<-2.5$., Comment: Accepted for publication in the Astronomy & Astrophysics Journal

It has recently been suggested that all giant stars with mass below 2 $M_{\odot}$ suffer an episode of surface lithium enrichment between the tip of the red giant branch (RGB) and the red clump (RC). We test if the above result can be confirmed in a sample of RC and RGB stars that are members of open clusters. We discuss Li abundances in six open clusters with ages between 1.5 and 4.9 Gyr (turn-off masses between 1.1 and 1.7 $M_{\odot}$). These observations are compared with the predictions of different models that include rotation-induced mixing, thermohaline instability, mixing induced by the first He flash, and energy losses by neutrino magnetic moment. In six clusters, we find about 35\% RC stars with Li abundances that are similar or higher than those of upper RGB stars. This can be a sign of fresh Li production. Because of the extra-mixing episode connected to the luminosity bump, the expectation was for RC stars to have systematically lower surface Li abundances. However, we cannot confirm that the possible Li production is ubiquitous. For about 65\% RC giants we can only determine abundance upper limits that could be hiding very low Li abundances. Our results indicate a possible production of Li during the RC, at levels that would not classify the stars as Li rich. Determination of their carbon isotopic ratio would help to confirm that the RC giants have suffered extra mixing followed by Li enrichment. The Li abundances of the RC stars can be qualitatively explained by the models with an additional mixing episode close to the He flash., Comment: Accepted for publication in A&A, 8 pages, 5 figures, tables available online and under request

Young open clusters (t<200 Myr) have been observed to exhibit several peculiarities in their chemical compositions, from a slightly sub-solar iron content, super-solar abundances of some atomic species (e.g. ionised chromium), and atypical enhancements of [Ba/Fe], with values up to +0.7 dex. Regarding the behaviour of the other $s$-process elements like yttrium, zirconium, lanthanum, and cerium, there is general disagreement in the literature. In this work we expand upon our previous analysis of a sample of five young open clusters (IC2391, IC2602, IC4665, NGC2516, and NGC2547) and one star-forming region (NGC2264), with the aim of determining abundances of different neutron-capture elements, mainly CuI, SrI, SrII, YII, ZrII, BaII, LaII, and CeII. We analysed high-resolution, high signal-to-noise spectra of 23 solar-type stars observed within the \textit{Gaia}-ESO survey. We find that our clusters have solar [Cu/Fe] within the uncertainties, while we confirm the super-solar [Ba/Fe] values (from +0.22 to +0.64 dex). Our analysis also points to mildly enhanced [Y/Fe] values (from 0 and +0.3 dex). For the other $s$-process elements we find that [X/Fe] ratios are solar at all ages. It is not possible to reconcile the anomalous behaviour of Ba and Y at young ages with standard stellar yields and Galactic chemical evolution model predictions. Thus, we explore different possible scenarios related to the behaviour of spectral lines, from the sensitivity to the presence of magnetic fields to the first ionisation potential effect. We also investigate the possibility that they may arise from alterations of the structure of the stellar photosphere due to higher levels of activity in such young stars. We are still unable to explain these enhancements, but we suggest that other elements (i.e. La) might be more reliable tracer of the $s$-process at young ages and encourage further observations., Comment: 26 pages, 14 figures, accepted for publication in A&A

We use the data provided by the Gaia Early Data Release 3 to search for a highly-complete volume-limited sample of unresolved binaries consisting of a white dwarf and a main sequence companion (i.e. WDMS binaries) within 100pc. We select 112 objects based on their location within the Hertzsprung-Russell diagram, of which 97 are new identifications. We fit their spectral energy distributions (SED) with a two-body fitting algorithm implemented in VOSA (Virtual Observatory SED Analyser) to derive the effective temperatures, luminosities and radii (hence surface gravities and masses) of both components. The stellar parameters are compared to those from the currently largest catalogue of close WDMS binaries, from the Sloan Digital Sky Survey (SDSS). We find important differences between the properties of the Gaia and SDSS samples. In particular, the Gaia sample contains WDMS binaries with considerably cooler white dwarfs and main sequence companions (some expected to be brown dwarfs). The Gaia sample also shows an important population of systems consisting of cool and extremely low-mass white dwarfs, not present in the SDSS sample. Finally, using a Monte Carlo population synthesis code, we find that the volume-limited sample of systems identified here seems to be highly complete (~80+-9 per cent), however it only represents ~9 per cent of the total underlying population. The missing ~91 per cent includes systems in which the main sequence companions entirely dominate the SEDs. We also estimate an upper limit to the total space density of close WDMS binaries of ~(3.7+-1.9)x10^{-4} pc{-3}., Comment: Accepted for publication by MNRAS

We exploit the unique characteristics of a sample of open clusters (OCs) and field stars for which high-precision 7Li abundances and stellar parameters are homogeneously derived by the Gaia-ESO Survey (GES). We derive possibly undepleted 7Li abundances for 26 OCs and star forming regions with ages from young to old spanning a large range of Galactocentric distances, which allows us to reconstruct the local late Galactic evolution of lithium as well as its current abundance gradient along the disc. Field stars are added to look further back in time and to constrain 7Li evolution in other Galactic components. The data are then compared to theoretical tracks from chemical evolution models that implement different 7Li forges. We find that the upper envelope of the 7Li abundances measured in field stars of nearly solar metallicities traces very well the level of lithium enrichment attained by the ISM as inferred from observations of cluster stars. We confirm previous findings that the abundance of 7Li in the solar neighbourhood does not decrease at supersolar metallicity. The comparison of the data with the chemical evolution model predictions favours a scenario in which the majority of the 7Li abundance in meteorites comes from novae. Current data also seem to suggest that the nova rate flattens out at later times. This requirement might have implications for the masses of the white dwarf nova progenitors and deserves further investigation. Neutrino-induced reactions taking place in core-collapse supernovae also produce some fresh lithium. This likely makes a negligible contribution to the meteoritic abundance, but could be responsible for a mild increase of the 7Li abundance in the ISM of low-metallicity systems that would counterbalance the astration processes., Comment: 17 pages, 10 figures, 3 tables (full tables 1 and 3 only available at the CDS), accepted for publication in Astronomy and Astrophysics. Version 2 includes corrections from language editor

We aim to constrain the mixing processes in low-mass stars by investigating the behaviour of the Li surface abundance after the main sequence. We take advantage of the data from the sixth internal data release of Gaia-ESO, idr6, and from the Gaia Early Data Release 3, edr3. We select a sample of main sequence, sub-giant, and giant stars in which Li abundance is measured by the Gaia-ESO survey, belonging to 57 open clusters with ages from 120~Myr to about 7 Gyr and to Milky Way fields, covering a range in [Fe/H] between -1.0 and +0.5dex. We study the behaviour of the Li abundances as a function of stellar parameters. We compare the observed Li behaviour in field giant stars and in giant stars belonging to individual clusters with the predictions of a set of classical models and of models with mixing induced by rotation and thermohaline instability. The comparison with stellar evolution models confirms that classical models cannot reproduce the lithium abundances observed in the metallicity and mass regimes covered by the data. The models that include the effects of both rotation-induced mixing and thermohaline instability account for the Li abundance trends observed in our sample, in all metallicity and mass ranges. The differences between the results of the classical models and of the rotation models largely differ (up to ~2 dex), making lithium the best element to constrain stellar mixing processes in low-mass stars. For stars with well-determined masses, we find a better agreement between observed surface abundances and models with rotation-induced and thermohaline mixings, the former dominating during the main sequence and the first phases of the post-main sequence evolution and the latter after the bump in the luminosity function., Comment: Accepted for publication in A&A, 19 pages, 16 figures

We present the results of a near-infrared (NIR) monitoring program carried out between 1999 and 2005 to determine the variability properties of the Arecibo sample of OH/IR stars. The sample consists of 385 IRAS-selected Asymptotic Giant Branch (AGB) candidates, for which their O-rich chemistry has been proven by the detection of 1612 MHz OH maser emission. The monitoring data was complemented by data collected from public optical and NIR surveys. We fitted the light curves obtained in the optical and NIR bands with a model using an asymmetric cosine function, and derived a period for 345 sources (90% of the sample). Based on their variability properties, most of the Arecibo sources are classified as long-period large-amplitude variable stars (LPLAV), 4% as (candidate) post-AGB stars, and 3% remain unclassified although they are likely post-AGB stars or highly obscured AGB stars. The period distribution of the LPLAVs peaks at 400 days, with periods between 300 and 800 days for most of the sources, and has a long tail up to 2100 days. Typically, the amplitudes are between 1 and 3 mag in the NIR and between 2 and 6 mag in the optical. We find correlations between periods and amplitudes, with larger amplitudes associated to longer periods, as well as between the period and the infrared colours, with the longer periods linked to the redder sources. Among the post-AGB stars, the light curve of IRAS 19566+3423 was exceptional, showing a large systematic increase in K-band brightness over 7 years., Comment: Published in MNRAS

We present the updated photometric calibration of the twelve optical passbands for the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2), comprising 1088 pointings of two square degrees, and study the systematic impact of metallicity in the stellar locus technique. The [Fe/H] metallicity from LAMOST DR5 for 146184 high-quality calibration stars, defined with S/N > 10 in J-PLUS passbands and S/N > 3 in Gaia parallax, was used to compute the metallicity-dependent stellar locus (ZSL). The initial homogenization of J-PLUS photometry, performed with a unique stellar locus, was refined by including the metallicity effect in colours via the ZSL. The variation of the average metallicity along the Milky Way produces a systematic offset in J-PLUS calibration. This effect is well above 1% for the bluer passbands and amounts 0.07, 0.07, 0.05, 0.03, and 0.02 mag in u, J0378, J0395, J0410, and J0430, respectively. We modelled this effect with the Milky Way location of the J-PLUS pointing, providing also an updated calibration for those observations without LAMOST information. The estimated accuracy in the calibration after including the metallicity effect is at 1% level for the bluer J-PLUS passbands and below for the rest. We conclude that photometric calibration with the stellar locus technique is prone to significant systematic bias along the Milky Way location for passbands bluer than lambda = 4500 A. The updated calibration method for J-PLUS DR2 reaches 1-2% precision and 1% accuracy for twelve optical filters within an area of 2176 square degrees., Comment: Published in Astronomy & Astrophysics. 18 pages, 12 figures, 3 tables, 2 appendices

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., Comment: 12 pages, 10 figures, accepted for publication at A&A

Clusterix 2.0 is a web-based, Virtual Observatory-compliant, interactive tool for the determination of membership probabilities in stellar clusters based on proper motion data using a fully non-parametric method. In the area occupied by the cluster, the frequency function is made up of two contributions: cluster and field stars. The tool performs an empirical determination of the frequency functions from the Vector-Point Diagram without relying in any previous assumption about their profiles. Clusterix 2.0 allows to search in an interactive way the appropriate spatial areas until an optimal separation of the two populations is obtained. Several parameters can be adjusted to make the calculation computationally feasible without interfering in the quality of the results. The system offers the possibility to query different catalogues, such as Gaia, or upload the user own data. The results of the membership determination can be sent via SAMP to VO tools like TopCat. We apply Clusterix 2.0 to several open clusters with different properties and environments to show the capabilities of the tool: an area of five degrees around NGC 2682 (M 67), an old, well known cluster; a young cluster NGC 2516 with a striking elongate structure extended up to four degrees; NGC 1750 & NGC 1758, a pair of partly overlapping clusters; in the area of NGC 1817 we confirm a little-known cluster, Juchert 23; and in an area with many clusters we disentangle the existence of two overlapping clusters where only one was previously known: Ruprecht 26 and the new, Clusterix 1., Comment: 17 pages, 22 figures + Appendix with 20 figures. Submitted to MNRAS

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%., Comment: Accepted for publication by MNRAS

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., Comment: Submitted to Astronomy and Astrophysics. 14 figures, 5 tables. Comments are welcome. Extra information about the photometry update can be found at www.j-plus.es/datareleases/dr1_swdl_calibration

In the era of large high-resolution spectroscopic surveys, high-quality spectra can contribute to our understanding of the Galactic chemical evolution, providing chemical abundances belonging to the different nucleosynthesis channels, and also providing constraints to stellar age. Some abundance ratios have been proven to be excellent indicators of stellar ages. We aim at providing an empirical relationship between stellar ages and [C/N] using, as calibrators, open star clusters observed by both the Gaia-ESO and APOGEE surveys. We use stellar parameters and abundances from the Gaia-ESO and APOGEE of the Galactic field and open cluster stars. Ages of star clusters are retrieved from the literature sources and validated using a common set of isochrones. We use the same isochrones to determine, for each age and metallicity, the surface gravity at which the first dredge-up and red giant branch bump occur. We study the effect of extra-mixing processes in our sample of giant stars, and we derive the mean [C/N] in evolved stars, including only stars without evidence of extra-mixing. Combining the Gaia-ESO and APOGEE samples of open clusters, we derive a linear relationship between [C/N] and logarithmic cluster ages. We apply our relationship to selected giant field stars in both Gaia-ESO and APOGEE. We find an age separation between thin and thick disc stars and age trends within their populations, with an increasing age towards lower metallicity populations. With such empirical relationship, we are able to provide an age estimate for giant stars in which C and N abundances are measured. Isochrone fitting is less sensitive for giant than dwarf stars at the turn off. The present method can be thus considered as an additional tool to give an independent estimate of the age of giant stars, with uncertainties in their ages comparable to those obtained using isochrone fitting for dwarf stars., Comment: 28 pages, 24 figures. Accepted by A&A

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

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., Comment: 18 pages, 11 figures and 3 tables. Accepted for publication in MNRAS

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., Comment: Accepted for publication in A&A; 23 pages, 17 figures, 3 tables

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/)., Comment: Accepted in AJ

The abundance ratio N/O is a useful tool to study the interplay of galactic processes, e.g. star formation efficiency, time-scale of infall and outflow loading factor We aim to trace log(N/O) versus [Fe/H] in the Milky Way and to compare it with a set of chemical evolution models to understand the role of infall, outflow and star formation efficiency in the building-up of the Galactic disc. We use the abundances from idr2-3, idr4, idr5 data releases of the Gaia-ESO Survey both for Galactic field and open cluster stars.We determine membership and average composition of open clusters and we separate thin and thick disc field stars.We consider the effect of mixing in the abundance of N in giant stars. We compute a grid of chemical evolution models, suited to reproduce the main features of our Galaxy, exploring the effects of the star formation efficiency, the infall time-scale and the differential outflow. With our samples, we map the metallicity range -0.6<[Fe/H]<0.3 with a corresponding -1.2

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., Comment: Accepted in MNRAS (2018 July 16)

Several works have found an increase of the abundances of the s-process neutron-capture elements in the youngest Galactic stellar populations, giving important constraints to stellar and Galactic evolution. We aim to trace the abundance patterns and the time-evolution of five s-process elements in the first peak, Y and Zr, and in the second peak, Ba, La and Ce using the Gaia-ESO idr5 results. From the UVES spectra of cluster member stars, we determined the average composition of clusters with ages >0.1 Gyr. We derived statistical ages and distances of field stars, and we separated them in thin and thick disc populations. We studied the time evolution and dependence on metallicity of abundance ratios using open clusters and field stars. Using our large and homogeneous sample of open clusters, thin and thick disc stars, spanning an age range larger than 10 Gyr, we confirm an increase towards young ages of s-process abundances in the Solar neighbourhood. These trends are well defined for open clusters and stars located nearby the solar position and they may be explained by a late enrichment due to significant contribution to the production of these elements from long-living low-mass stars. At the same time, we found a strong dependence of the s-process abundance ratios with the Galactocentric distance and with the metallicity of the clusters and field stars. Our results, derived from the largest and homogeneous sample of s-process abundances in the literature, confirm the growth with decreasing stellar ages of the s-process abundances in both field and open cluster stars. At the same time, taking advantage of the abundances of open clusters located in a wide Galactocentric range, they open a new view on the dependence of the s-process evolution on the metallicity and star formation history, pointing to different behaviours at various Galactocentric distances., Comment: 16 pages, 11 figures, accepted for publication in A&A

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., Comment: Accepted for publication in A&A ; 12 pages, 12 figures

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., Comment: 14 pages, 6 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

We present a spectroscopic study of the dynamics of the ionized and neutral gas throughout the Lagoon nebula (M8), using VLT/FLAMES data from the Gaia-ESO Survey. We explore the connections between the nebular gas and the stellar population of the associated star cluster NGC6530. We characterize through spectral fitting emission lines of H-alpha, [N II] and [S II] doublets, [O III], and absorption lines of sodium D doublet, using data from the FLAMES/Giraffe and UVES spectrographs, on more than 1000 sightlines towards the entire face of the Lagoon nebula. Gas temperatures are derived from line-width comparisons, densities from the [S II] doublet ratio, and ionization parameter from H-alpha/[N II] ratio. Although doubly-peaked emission profiles are rarely found, line asymmetries often imply multiple velocity components along the line of sight. This is especially true for the sodium absorption, and for the [O III] lines. Spatial maps for density and ionization are derived, and compared to other known properties of the nebula and of its massive stars 9 Sgr, Herschel 36 and HD 165052 which are confirmed to provide most of the ionizing flux. The detailed velocity fields across the nebula show several expanding shells, related to the cluster NGC6530, the O stars 9 Sgr and Herschel 36, and the massive protostar M8East-IR. The origins of kinematical expansion and ionization of the NGC6530 shell appear to be different. We are able to put constrains on the line-of-sight (relative or absolute) distances between some of these objects and the molecular cloud. The large obscuring band running through the middle of the nebula is being compressed by both sides, which might explain its enhanced density. We also find an unexplained large-scale velocity gradient across the entire nebula. At larger distances, the transition from ionized to neutral gas is studied using the sodium lines., Comment: 26 pages, 31 figures, accepted on Astronomy and Astrophysics journal