Source author record

A. C. Robin

A. C. Robin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2022arXiv

3D extinction mapping of the Milky Way using Convolutional Neural Networks: Presentation of the method and demonstration in the Carina Arm region

Context. Several methods have been proposed to build 3D extinction maps of the Milky Way (MW), most often based on Bayesian approaches. Although some studies employed machine learning (ML) methods in part of their procedure, or to specific targets, no 3D extinction map of a large volume of the MW solely based on a Neural Network method has been reported so far. Aims. We aim to apply deep learning as a solution to build 3D extinction maps of the MW. Methods. We built a convolutional neural network (CNN) using the CIANNA framework, and trained it with synthetic 2MASS data. We used the Besançon Galaxy model to generate mock star catalogs, and 1D Gaussian random fields to simulate the extinction profiles. From these data we computed color-magnitude diagrams (CMDs) to train the network, using the corresponding extinction profiles as targets. A forward pass with observed 2MASS CMDs provided extinction profile estimates for a grid of lines of sight. Results. We trained our network with data simulating lines of sight in the area of the Carina spiral arm tangent and obtained a 3D extinction map for a large sector in this region ($l = 257 - 303$ deg, $|b| \le 5$ deg), with distance and angular resolutions of $100$ pc and $30$ arcmin, respectively, and reaching up to $\sim 10$ kpc. Although each sightline is computed independently in the forward phase, the so-called fingers-of-God artifacts are weaker than in many other 3D extinction maps. We found that our CNN was efficient in taking advantage of redundancy across lines of sight, enabling us to train it with only 9 sightlines simultaneously to build the whole map. Conclusions. We found deep learning to be a reliable approach to produce 3D extinction maps from large surveys. With this methodology, we expect to easily combine heterogeneous surveys without cross-matching, and therefore to exploit several surveys in a complementary fashion.

preprint2022arXiv

Earth through the looking glass: how frequently are we detected by other civilisations through photometric microlensing?

Microlensing is proving to be one of the best techniques to detect distant, low-mass planets around the most common stars in the Galaxy. In principle, Earth's microlensing signal could offer the chance for other technological civilisations to find the Earth across Galactic distances. We consider the photometric microlensing signal of Earth to other potential technological civilisations and dub the regions of our Galaxy from which Earth's photometric microlensing signal is most readily observable as the "Earth Microlensing Zone" (EMZ). The EMZ can be thought of as the microlensing analogue of the Earth Transit Zone (ETZ) from where observers see Earth transit the Sun. Just as for the ETZ, the EMZ could represent a game-theoretic Schelling point for targeted searches for extra-terrestrial intelligence (SETI). To compute the EMZ, we use the Gaia DR2 catalogue with magnitude G<20 to generate Earth microlensing probability and detection rate maps to other observers. Whilst our Solar system is a multi-planet system, we show that Earth's photometric microlensing signature is almost always well approximated by a binary lens assumption. We then show that the Earth is in fact well-hidden to observers with technology comparable to our own. Specifically, even if observers are located around every Gaia DR2 star with G<20, we expect photometric microlensing signatures from the Earth to be observable on average only tens per year by any of them. In addition, the EMZs overlap with the ETZ near the Galactic centre which could be the main areas for future SETI searches.

preprint2016arXiv

Chemical tagging with APOGEE: Discovery of a large population of N-rich stars in the inner Galaxy

Formation of globular clusters (GCs), the Galactic bulge, or galaxy bulges in general, are important unsolved problems in Galactic astronomy. Homogeneous infrared observations of large samples of stars belonging to GCs and the Galactic bulge field are one of the best ways to study these problems. We report the discovery by APOGEE of a population of field stars in the inner Galaxy with abundances of N, C, and Al that are typically found in GC stars. The newly discovered stars have high [N/Fe], which is correlated with [Al/Fe] and anti-correlated with [C/Fe]. They are homogeneously distributed across, and kinematically indistinguishable from, other field stars in the same volume. Their metallicity distribution is seemingly unimodal, peaking at [Fe/H]~-1, thus being in disagreement with that of the Galactic GC system. Our results can be understood in terms of different scenarios. N-rich stars could be former members of dissolved GCs, in which case the mass in destroyed GCs exceeds that of the surviving GC system by a factor of ~8. In that scenario, the total mass contained in so-called "first-generation" stars cannot be larger than that in "second-generation" stars by more than a factor of ~9 and was certainly smaller. Conversely, our results may imply the absence of a mandatory genetic link between "second generation" stars and GCs. Last, but not least, N-rich stars could be the oldest stars in the Galaxy, the by-products of chemical enrichment by the first stellar generations formed in the heart of the Galaxy.

preprint2016arXiv

Close encounters involving RAVE stars beyond the 47 Tucanae tidal radius

The most accurate 6D phase-space information from the Radial Velocity Experiment (RAVE) was used to integrate the orbits of 105 stars around the galactic globular cluster 47 Tucanae, to look for close encounters between them in the past, with a minimum distance approach less than the cluster tidal radius. The stars are currently over the distance range 3.0 kpc $<$ d $<$ 5.5 kpc. Using the uncertainties in the current position and velocity vector for both, star and cluster, 105 pairs of star-cluster orbits were generated in a Monte Carlo numerical scheme, integrated over 2 Gyr and considering an axisymmetric and non-axisymmetric Milky-Way-like Galactic potential, respectively. In this scheme, we identified 20 potential cluster members that had close encounters with the globular cluster 47 Tucanae, all of which have a relative velocity distribution (V$_{rel}$) less than 200 km s$^{-1}$ at the minimum distance approach. Among these potential members, 9 had close encounters with the cluster with velocities less than the escape velocity of 47 Tucanae, therefore a scenario of tidal stripping seems likely. These stars have been classified with a 93\% confidence level, leading to the identification of extratidal cluster stars. For the other 11 stars, V$_{rel}$ exceeds the escape velocity of the cluster, therefore they were likely ejected or are unassociated interlopers.

preprint2016arXiv

Discovery of a Metal-Poor Field Giant with a Globular Cluster Second-Generation Abundance Pattern

We report on detection, from observations obtained with the APOGEE spectroscopic survey, of a metal-poor ([Fe/H] $= -1.3$ dex) field giant star with an extreme Mg-Al abundance ratio ([Mg/Fe] $= -0.31$ dex; [Al/Fe] $= 1.49$ dex). Such low Mg/Al ratios are seen only among the second-generation population of globular clusters, and are not present among Galactic disk field stars. The light element abundances of this star, 2M16011638-1201525, suggest that it could have been born in a globular cluster. We explore several origin scenarios, in particular studying the orbit of the star to check the probability of it being kinematically related to known globular clusters. We performed simple orbital integrations assuming the estimated distance of 2M16011638-1201525 and the available six-dimensional phase-space coordinates of 63 globular clusters, looking for close encounters in the past with a minimum distance approach within the tidal radius of each cluster. We found a very low probability that 2M16011638-1201525 was ejected from most globular clusters; however, we note that the best progenitor candidate to host this star is globular cluster $ω$ Centauri (NGC 5139). Our dynamical investigation demonstrates that 2M16011638-1201525 reaches a distance $|Z_{max}| < 3 $ kpc from the Galactic plane and a minimum and maximum approach to the Galactic center of $R_{min}<0.62$ kpc and $R_{max}<7.26$ kpc in an eccentric ($e\sim0.53$) and retrograde orbit. Since the extreme chemical anomaly of 2M16011638-1201525 has also been observed in halo field stars, this object could also be considered a halo contaminant, likely been ejected into the Milky Way disk from the halo. We conclude that, 2M16011638-20152 is also kinematically consistent with the disk but chemically consistent with halo field stars.

preprint2016arXiv

The microlensing rate and distribution of free-floating planets towards the Galactic bulge

Ground-based optical microlensing surveys have provided tantalising, if inconclusive, evidence for a significant population of free-floating planets (FFPs). Both ground and space-based facilities are being used and developed which will be able to probe the distrubution of FFPs with much better sensitivity. It is vital also to develop a high-precision microlensing simulation framework to evaluate the completeness of such surveys. We present the first signal-to-noise limited calculations of the FFP microlensing rate using the Besancon Galactic model. The microlensing distribution towards the Galactic centre is simulated for wide-area ground-based optical surveys such as OGLE or MOA, a wide-area ground-based near-IR survey, and a targeted space-based near-IR survey which could be undertaken with Euclid or WFIRST. We present a calculation framework for the computation of the optical and near-infrared microlensing rate and optical depth for simulated stellar catalogues which are signal-to-noise limited, and take account of extinction, unresolved stellar background light and finite source size effects, which can be significant for FFPs. We find that the global ground-based I-band yield over a central 200 deg^2 region covering the Galactic centre ranges from 20 Earth-mass FFPs year^-1 up to 3,500 year^-1 for Jupiter FFPs in the limit of 100% detection efficiency, and almost an order of magnitude larger for a K-band survey. For ground-based surveys we find that the inclusion of finite source and the unresolved background reveals a mass-dependent variation in the spatial distribution of FFPs. For a space-based H-band covering 2 deg^2, the yield depends on the target field but maximizes close to the Galactic centre with around 76 Earth through to 1,700 Jupiter FFPs year^-1. For near-IR space-based surveys the spatial distribution of FFPs is found to be largely insensitive to the FFP mass scale.

preprint2015arXiv

Deep SDSS optical spectroscopy of distant halo stars II. Iron, calcium, and magnesium abundances

We analyze a sample of 3,944 low-resolution (R ~ 2000) optical spectra from the Sloan Digital Sky Survey (SDSS), focusing on stars with effective temperatures 5800 < Teff < 6300 K, and distances from the Milky Way plane in excess of 5 kpc, and determine their abundances of Fe, Ca, and Mg. We followed the same methodology as in the previous paper in this series, deriving atmospheric parameters by chi2 minimization, but this time we obtained the abundances of individual elements by fitting their associated spectral lines. Distances were calculated from absolute magnitudes obtained by a statistical comparison of our stellar parameters with stellar-evolution models. The observations reveal a decrease in the abundances of iron, calcium, and magnesium at large distances from the Galactic center. The median abundances for the halo stars analyzed are fairly constant up to a Galactocentric distance r ~ 20 kpc, rapidly decrease between r ~ 20 and r ~ 40 kpc, and flatten out to significantly lower values at larger distances, consistent with previous studies. In addition, we examine the [Ca/Fe] and [Mg/Fe] as a function of Fe/H and Galactocentric distance. Our results show that the most distant parts of the halo show a steeper variation of the [Ca/Fe] and [Mg/Fe] with iron. We found that at the range -1.6 < [Fe/H] < -0.4 [Ca/Fe] decreases with distance, in agreement with earlier results based on local stars. However, the opposite trend is apparent for [Mg/Fe]. Our conclusion that the outer regions of the halo are more metal-poor than the inner regions, based on in situ observations of distant stars, agrees with recent results based on inferences from the kinematics of more local stars, and with predictions of recent galaxy formation simulations for galaxies similar to the Milky Way.

preprint2015arXiv

Mapping the inner stellar halo of the Milky Way from 2MASS and SDSS-III/APOGEE survey

The Besançon Galaxy model was used to compare the infrared colour distribution of synthetic stars with those from 2MASS observations taking the selection function of the data into account, in order to study the shape of the stellar halo of the Milky Way, with complemetary spectroscopic data from SDSS-III/APOGEE survey. Furthermore, we compared the generated mock metallicity distribution of the Besançon Galaxy model, to the intrinsic metallicity distribution with reliable stellar parameters from the APOGEE Stellar Parameters and Chemical Abundances Pipeline (ASPCAP). The comparison was carried accross a large volume of the inner part of the Galaxy, revealing that a metal-poor population, [M/H]$<-1.2$ dex, could fill an extended component of the inner galactic halo. With this data set, we are able to model a more realistic mass density distribution of the stellar halo component of the Milky Way, assuming a six-parameters double power-law model, and reconstruct the behaviour of the rotation curve in the inner part of the Galaxy.

preprint2015arXiv

Quasi integral of motion for axisymmetric potentials

We present an estimate of the third integral of motion for axisymmetric three-dimensional potentials. This estimate is based on a Staeckel approximation and is explicitly written as a function of the potential. We tested this scheme for the Besancon Galactic model and two other disc-halo models and find that orbits of disc stars have an accurately conserved third quasi integral. The accuracy ranges from of 0.1% to 1% for heights varying from z = 0~kpc to z= 6 kpc and Galactocentric radii R from 5 to 15kpc. We also tested the usefulness of this quasi integral in analytic distribution functions of disc stellar populations: we show that the distribution function remains approximately stationary and that it allows to recover the potential and forces by applying Jeans equations to its moments.

preprint2015arXiv

RAVE stars tidally stripped/ejected from $ω$ Centauri globular cluster

Using six-dimesional phase-space information from the Fourth Data release of the Radial Velocity Experiment (RAVE) over the range of Galactic longitude 240$^{\circ}< l <$ 360$^{\circ}$ and $V_{LSR} < -239$ kms$^{-1}$, we have computed orbits for 329 RAVE stars that were originally selected as chemically and kinematically related to $ω$ Centauri. The orbits were integrated in a Milky-Way-like axisymmetric Galactic potential, ignoring the effects of the dynamical evolution of $ω$ Centauri due to the tidal effects of the Galaxy disk on the cluster along time. We also ignored secular changes in the Milky Way potential over time. In a Monte Carlo scheme, and under the assumption that the stars may have been ejected with velocities greater than the escape velocity ($V_{rel}>V_{esc,0}$) from the cluster, we identified 15 stars as having close encounters with $ω$ Centauri: (\textit{i}) 8 stars with relative velocities $V_{rel}< 200 $ kms$^{-1}$ may have been ejected $\sim$ 200 Myr ago from $ω$ Centauri; (\textit{ii}) other group of 7 stars were identified with high relative velocity $V_{rel}> 200 $ kms$^{-1}$ during close encounters, and seems unlikely that they have been ejected from $ω$ Centauri. We also confirm the link between J131340.4-484714 as potential member of $ω$ Centauri, and probably ejected $\sim$ 2.0 Myr ago, with a relative velocity $V_{rel}\sim80$ kms$^{-1}$.

preprint2014arXiv

A Stellar Population Synthesis Model for the Study of Ultraviolet Star Counts of the Galaxy

GALEX, the first all sky imaging UV satellite, has imaged a large part of the sky providing an excellent opportunity for studying UV star counts. The aim of our study is to investigate in detail the observed UV star counts obtained by GALEX vis-a-vis the model simulated catalogs produced by the Besancon model of stellar population synthesis in various Galactic directions, and to explore the potential for studying the structure of our Galaxy from images in multiple NUV and FUV filters of the forthcoming Ultraviolet Imaging Telescope (UVIT) to be flown onboard ASTROSAT. We have upgraded the Besancon model of stellar population synthesis to include the UV bands of GALEX and UVIT. Depending on the availability of contiguous GALEX, SDSS, WISE and 2MASS overlapping regions, we have chosen a set of 19 GALEX fields which spread over a range of Galactic directions. We cross-matched GALEX sources with the WISE+2MASS and SDSS catalogs and UV stars in the GALEX catalog are identified by choosing a suitable IR colour, J - W1 (W1 is a WISE band at 3.4 microns). The IR colour cut method, which is used for the first time for separation of stars, is discussed in comparison with the GALEX+SDSS star counts method. We present the results of the UV star counts analysis carried out using the data from GALEX. We find that the Besancon model simulations represent the observed star counts of both the GALEX AIS and MIS well within the error bars in various Galactic directions. Based on the model analysis, we separated out white dwarfs of the disc and blue horizontal branch stars of the halo from the observed sample by selecting a suitable FUV - NUV colour. The Besancon model is now ready for further comparisons in the UV domain and will be used for prospective studies for the UVIT instrument to be flown onboard ASTROSAT.

preprint2014arXiv

Chemodynamics of the Milky Way. I. The first year of APOGEE data

We investigate the chemo-kinematic properties of the Milky Way disc by exploring the first year of data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE), and compare our results to smaller optical high-resolution samples in the literature, as well as results from lower resolution surveys such as GCS, SEGUE and RAVE. We start by selecting a high-quality sample in terms of chemistry ($\sim$ 20.000 stars) and, after computing distances and orbital parameters for this sample, we employ a number of useful subsets to formulate constraints on Galactic chemical and chemodynamical evolution processes in the Solar neighbourhood and beyond (e.g., metallicity distributions -- MDFs, [$α$/Fe] vs. [Fe/H] diagrams, and abundance gradients). Our red giant sample spans distances as large as 10 kpc from the Sun. We find remarkable agreement between the recently published local (d $<$ 100 pc) high-resolution high-S/N HARPS sample and our local HQ sample (d $<$ 1 kpc). The local MDF peaks slightly below solar metallicity, and exhibits an extended tail towards [Fe/H] $= -$1, whereas a sharper cut-off is seen at larger metallicities. The APOGEE data also confirm the existence of a gap in the [$α$/Fe] vs. [Fe/H] abundance diagram. When expanding our sample to cover three different Galactocentric distance bins, we find the high-[$α$/Fe] stars to be rare towards the outer zones, as previously suggested in the literature. For the gradients in [Fe/H] and [$α$/Fe], measured over a range of 6 $ < $ R $ <$ 11 kpc in Galactocentric distance, we find a good agreement with the gradients traced by the GCS and RAVE dwarf samples. For stars with 1.5 $<$ z $<$ 3 kpc, we find a positive metallicity gradient and a negative gradient in [$α$/Fe].

preprint2014arXiv

Constraining the thick disc formation scenario of the Milky Way

We study the shape of the thick disc using photometric data at high and intermediate latitudes from SDSS and 2MASS surveys. We use the population synthesis approach using an ABC-MCMC method to characterize the thick disc shape, scale height, scale length, local density and flare, and we investigate the extend of the thick disc formation period by simulating several formation episodes. We find that the vertical variation in density is not exponential, but much closer to an hyperbolic secant squared. Assuming a single formation epoch, the thick disc is better fitted with a sech2 scale height of 470 pc and a scale length of 2.3 kpc. However if one simulates two successive formation episodes, mimicking an extended formation period, the older episode has a higher scale height and a larger scale length than the younger episode, indicating a contraction during the collapse phase. The scale height decreases from 800 pc to 340 pc, and the scale length from 3.2 kpc to 2 kpc. The star formation increases from the old episode to the young one. During the fitting process, the halo parameters are also determined. The constraint on the halo shows that a transition between a inner and outer halo, if exists, cannot be at a distance of less than about 30 kpc, which is the limit of our investigation using turnoff halo stars. Finally, we show that extrapolating the thick disc towards the bulge region explains well the stellar populations observed there, that there is no longer need to invoke a classical bulge. To explain these results, the most probable scenario for the thick disc is that it formed while the Galaxy was gravitationally collapsing from well mixed gas-rich giant clumps sustained by high turbulence which awhile prevent a thin disc to form, as proposed by Bournaud et al. (2009). This scenario explains the observations in the thick disc region as well as in the bulge region. (abridged)

preprint2014arXiv

Deep SDSS optical spectroscopy of distant halo stars I. Atmospheric parameters and stellar metallicity distribution

We analyze a sample of tens of thousands of spectra of halo turnoff stars, obtained with the optical spectrographs of the Sloan Digital Sky Survey (SDSS), to characterize the stellar halo population "in situ" out to a distance of a few tens of kpc from the Sun. In this paper we describe the derivation of atmospheric parameters. We also derive the overall stellar metallicity distribution based on F-type stars observed as flux calibrators for the Baryonic Oscillations Spectroscopic Survey (BOSS). Our analysis is based on an automated method that determines the set of parameters of a model atmosphere that reproduces each observed spectrum best. We used an optimization algorithm and evaluate model fluxes by means of interpolation in a precomputed grid. In our analysis, we account for the spectrograph's varying resolution as a function of fiber and wavelength. Our results for early SDSS (pre-BOSS upgrade) data compare well with those from the SEGUE Stellar Parameter Pipeline (SSPP), except for stars with logg (cgs units) lower than 2.5. An analysis of stars in the globular cluster M13 reveals a dependence of the inferred metallicity on surface gravity for stars with logg < 2.5, confirming the systematics identified in the comparison with the SSPP. We find that our metallicity estimates are significantly more precise than the SSPP results. We obtain a halo metallicity distribution that is narrower and more asymmetric than in previous studies. The lowest gravity stars in our sample, at tens of kpc from the Sun, indicate a shift of the metallicity distribution to lower abundances, consistent with what is expected from a dual halo system in the Milky Way.

preprint2014arXiv

ExELS: an exoplanet legacy science proposal for the ESA Euclid mission. II. Hot exoplanets and sub-stellar systems

The Exoplanet Euclid Legacy Survey (ExELS) proposes to determine the frequency of cold exoplanets down to Earth mass from host separations of ~1 AU out to the free-floating regime by detecting microlensing events in Galactic Bulge. We show that ExELS can also detect large numbers of hot, transiting exoplanets in the same population. The combined microlensing+transit survey would allow the first self-consistent estimate of the relative frequencies of hot and cold sub-stellar companions, reducing biases in comparing "near-field" radial velocity and transiting exoplanets with "far-field" microlensing exoplanets. The age of the Bulge and its spread in metallicity further allows ExELS to better constrain both the variation of companion frequency with metallicity and statistically explore the strength of star-planet tides. We conservatively estimate that ExELS will detect ~4100 sub-stellar objects, with sensitivity typically reaching down to Neptune-mass planets. Of these, ~600 will be detectable in both Euclid's VIS (optical) channel and NISP H-band imager, with ~90% of detections being hot Jupiters. Likely scenarios predict a range of 2900-7000 for VIS and 400-1600 for H-band. Twice as many can be expected in VIS if the cadence can be increased to match the 20-minute H-band cadence. The separation of planets from brown dwarfs via Doppler boosting or ellipsoidal variability will be possible in a handful of cases. Radial velocity confirmation should be possible in some cases, using 30-metre-class telescopes. We expect secondary eclipses, and reflection and emission from planets to be detectable in up to ~100 systems in both VIS and NISP-H. Transits of ~500 planetary-radius companions will be characterised with two-colour photometry and ~40 with four-colour photometry (VIS,YJH), and the albedo of (and emission from) a large sample of hot Jupiters in the H-band can be explored statistically.

preprint2014arXiv

Gaia photometry for white dwarfs

Context. White dwarfs can be used to study the structure and evolution of the Galaxy by analysing their luminosity function and initial mass function. Among them, the very cool white dwarfs provide the information for the early ages of each population. Because white dwarfs are intrinsically faint only the nearby (about 20 pc) sample is reasonably complete. The Gaia space mission will drastically increase the sample of known white dwarfs through its 5-6 years survey of the whole sky up to magnitude V = 20-25. Aims. We provide a characterisation of Gaia photometry for white dwarfs to better prepare for the analysis of the scientific output of the mission. Transformations between some of the most common photometric systems and Gaia passbands are derived. We also give estimates of the number of white dwarfs of the different galactic populations that will be observed. Methods. Using synthetic spectral energy distributions and the most recent Gaia transmission curves, we computed colours of three different types of white dwarfs (pure hydrogen, pure helium, and mixed composition with H/He= 0.1). With these colours we derived transformations to other common photometric systems (Johnson-Cousins, Sloan Digital Sky Survey, and 2MASS). We also present numbers of white dwarfs predicted to be observed by Gaia. Results. We provide relationships and colour-colour diagrams among different photometric systems to allow the prediction and/or study of the Gaia white dwarf colours. We also include estimates of the number of sources expected in every galactic population and with a maximum parallax error. Gaia will increase the sample of known white dwarfs tenfold to about 200 000. Gaia will be able to observe thousands of very cool white dwarfs for the first time, which will greatly improve our understanding of these stars and early phases of star formation in our Galaxy.

preprint2014arXiv

Metallicity and kinematics of the bar in-situ

Constraints on the Galactic bulge/bar structure and formation history from stellar kinematics and metallicities mainly come from relatively high-latitude fields (|b|>4) where a complex mix of stellar population is seen. We aim here to constrain the formation history of the Galactic bar by studying the radial velocity and metallicity distributions of stars in-situ (|b|<1). We observed red clump stars in four fields along the bar's major axis (l=10,6,-6 and b=0 plus a field at l=0,b=1) with low-resolution spectroscopy from VLT/FLAMES, observing around the CaII triplet. We developed robust methods for extracting radial velocity and metallicity estimates from these low signal-to-noise spectra. We derived distance probability distributions using Bayesian methods rigorously handling the extinction law. We present radial velocities and metallicity distributions, as well as radial velocity trends with distance. We observe an increase in the radial velocity dispersion near the Galactic plane. We detect the streaming motion of the stars induced by the bar in fields at l=+/-6, the highest velocity components of this bar stream being metal-rich ([Fe/H]~0.2 dex). Our data is consistent with a bar inclined at 26+/-3 from the Sun-Galactic centre line. We observe a significant fraction of metal-poor stars, in particular in the field at l=0,b=1. We confirm the flattening of the metallicity gradient along the minor axis when getting closer to the plane, with a hint that it could actually be inverted. Our stellar kinematics corresponds to the expected behaviour of a bar issued from the secular evolution of the Galactic disc. The mix of several populations, seen further away from the plane, is also seen in the bar in-situ since our metallicity distributions highlight a different spatial distribution between metal-poor and metal-rich stars, the more metal-poor stars being more centrally concentrated.

preprint2014arXiv

Overview and stellar statistics of the expected Gaia Catalogue using the Gaia Object Generator

Aims: An effort has been undertaken to simulate the expected Gaia Catalogue, including the effect of observational errors. A statistical analysis of this simulated Gaia data is performed in order to better understand what can be obtained from the Gaia astrometric mission. This catalogue is used in order to investigate the potential yield in astrometric, photometric and spectroscopic information, and the extent and effect of observational errors on the true Gaia Catalogue. This article is a follow-up to Robin et. al. (2012), where the expected Gaia Catalogue content was reviewed but without the simulation of observational errors. Methods: The Gaia Object Generator (GOG) catalogue is analysed using the Gaia Analysis Tool (GAT), producing a number of statistics on the catalogue. Results: A simulated catalogue of one billion objects is presented, with detailed information on the 523 million individual single stars it contains. Detailed information is provided for the expected errors in parallax, position, proper motion, radial velocity, photometry in the four Gaia bands, and physical parameter determination including temperature, metallicity and line of sight extinction.

preprint2014arXiv

Searching for tidal tails around $ω$ Centauri using RR Lyrae Stars

We present a survey for RR Lyrae stars in an area of 50 deg$^2$ around the globular cluster $ω$ Centauri, aimed to detect debris material from the alleged progenitor galaxy of the cluster. We detected 48 RR Lyrae stars of which only 11 have been previously reported. Ten among the eleven previously known stars were found inside the tidal radius of the cluster. The rest were located outside the tidal radius up to distances of $\sim 6$ degrees from the center of the cluster. Several of those stars are located at distances similar to that of $ω$ Centauri. We investigated the probability that those stars may have been stripped off the cluster by studying their properties (mean periods), calculating the expected halo/thick disk population of RR Lyrae stars in this part of the sky, analyzing the radial velocity of a sub-sample of the RR Lyrae stars, and finally, studying the probable orbits of this sub-sample around the Galaxy. None of these investigations support the scenario that there is significant tidal debris around $ω$ Centauri, confirming previous studies in the region. It is puzzling that tidal debris have been found elsewhere but not near the cluster itself.

preprint2014arXiv

The Besancon Galaxy model renewed I. Constraints on the local star formation history from Tycho data

The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allows us to test hypotheses on the star formation history, star evolution, and chemical and dynamical evolution of the Galaxy. The new version of the Besancon Galaxy model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the initial mass function (IMF) and star formation rate (SFR) of the Galactic disc. We present a new strategy for the generation of thin disc stars, which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. The local dynamical self-consistency is maintained in this new scheme. We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. For the first time, a full sky comparison is performed between BGM and data. This strategy allows us to constrain the IMF slope at high masses, which is found to be close to 3.0 and excludes a shallower slope such as Salpeter's one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 Mo/pc^3 implying that there is no compelling evidence for the significant amount of dark matter in the disc. While the model is fitted to Tycho-2 data, which is a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

preprint2013arXiv

3D kinematics through the X-shaped Milky Way bulge

It has recently been discovered that the Galactic bulge is X-shaped, with the two southern arms of the X both crossing the lines of sight at l=0 and |b|>4, hence producing a double red clump (RC) in the bulge CMD. Dynamical models predict the formation of X-shaped bulges, as extreme cases of boxy-peanut bulges. However, since X-shaped bulges were known to be present only in external galaxies, models have never been compared to 3D kinematical data for individual stars. We study the orbital motion of Galactic bulge stars, in the two arms of the X in the southern hemisphere. The goal is to provide observational constraints to bulge formation models that predict the formation of X-shapes through bar dynamical instabilities. Radial velocities have been obtained for a sample of 454 bulge giants, roughly equally distributed between the bright and faint RC, in a field at (l,b)=(0,-6). Proper motions were derived for all RC stars in the same field by combining images from two epochs obtained 11 years apart. The proper motions for the spectroscopic subsample are analyzed taking into account the radial velocities and metallicities measured from near-IR Calcium triplet lines, for a subsample of 352 stars. The radial velocity distribution of stars in the bright RC, tracing the closer overdensity of bulge stars, shows an excess of stars moving towards the Sun. Similarly, an excess of stars receding from the Sun is seen in the far overdensity, traced by faint RC stars. This is explained with the presence of stars on elongated orbits, most likely streaming along the arms of the X-shaped bulge. Proper motions for these stars are consistent with qualitative predictions of dynamical models of peanut-shaped bulges. Surprisingly, stars on elongated orbits have preferentially metal poor metallicities, while the metal rich ones, in both overdensities, are preferentially found in more axysimmetric orbits.

preprint2013arXiv

A Sub-Earth-Mass Moon Orbiting a Gas Giant Primary or a High Velocity Planetary System in the Galactic Bulge

We present the first microlensing candidate for a free-floating exoplanet-exomoon system, MOA-2011-BLG-262, with a primary lens mass of M_host ~ 4 Jupiter masses hosting a sub-Earth mass moon. The data are well fit by this exomoon model, but an alternate star+planet model fits the data almost as well. Nevertheless, these results indicate the potential of microlensing to detect exomoons, albeit ones that are different from the giant planet moons in our solar system. The argument for an exomoon hinges on the system being relatively close to the Sun. The data constrain the product M pi_rel, where M is the lens system mass and pi_rel is the lens-source relative parallax. If the lens system is nearby (large pi_rel), then M is small (a few Jupiter masses) and the companion is a sub-Earth-mass exomoon. The best-fit solution has a large lens-source relative proper motion, mu_rel = 19.6 +- 1.6 mas/yr, which would rule out a distant lens system unless the source star has an unusually high proper motion. However, data from the OGLE collaboration nearly rule out a high source proper motion, so the exoplanet+exomoon model is the favored interpretation for the best fit model. However, the alternate solution has a lower proper motion, which is compatible with a distant (so stellar) host. A Bayesian analysis does not favor the exoplanet+exomoon interpretation, so Occam's razor favors a lens system in the bulge with host and companion masses of M_host = 0.12 (+0.19 -0.06) M_solar and m_comp = 18 (+28 -100 M_earth, at a projected separation of a_perp ~ 0.84 AU. The existence of this degeneracy is an unlucky accident, so current microlensing experiments are in principle sensitive to exomoons. In some circumstances, it will be possible to definitively establish the low mass of such lens systems through the microlensing parallax effect. Future experiments will be sensitive to less extreme exomoons.

preprint2013arXiv

CFBDS J111807-064016: A new L/T transition brown dwarf in a binary system

Stellar-substellar binary systems are quite rare, and provide interesting benchmarks. They constrain the complex physics of substellar atmospheres, because several physical parameters of the substellar secondary can be fixed from the much better characterized main sequence primary. We report the discovery of CFBDS J111807-064016, a T2 brown dwarf companion to 2MASS J111806.99-064007.8, a low-mass M4.5-M5 star. The brown-dwarf was identified from the Canada France Brown Dwarf Survey. At a distance of 50-120 pc, the 7.7 arcsec angular separation corresponds to projected separations of 390-900 AU. The primary displays no Halpha emission, placing a lower limit on the age of the system of about 6 Gyr. The kinematics is also consistent with membership in the old thin disc. We obtained near-infrared spectra, which together with recent atmosphere models allow us determine the effective temperature and gravity of both components. From these parameters and the age constraint, evolutionary models estimate masses of 0.10 to 0.15 Msol for the M dwarf, and 0.06 to 0.07 Msol for the T dwarf. This system is a particularly valuable benchmark because the brown dwarf is an early T: the cloud-clearing that occurs at the L/T transition is very sensitive to gravity, metallicity, and detailed dust properties, and produces a large scatter in the colours. This T2 dwarf, with its metallicity measured from the primary and its mass and gravity much better constrained than those of younger early-Ts, will anchor our understanding of the colours of L/T transition brown dwarfs. It is also one of the most massive T dwarfs, just below the hydrogen-burning limit, and all this makes it a prime probe a brown dwarf atmosphere and evolution models.

preprint2013arXiv

ExELS: an exoplanet legacy science proposal for the ESA Euclid mission I. Cold exoplanets

Euclid is the second M-class mission of the ESA Cosmic Vision programme, with the principal science goal of studying dark energy. Euclid is also expected to undertake additional Legacy Science programmes. One proposal is the Exoplanet Euclid Legacy Survey (ExELS) which will be the first survey able to measure the abundance of exoplanets down to Earth mass for host separations from ~1AU out to the free-floating (unbound) regime. The cold and free-floating exoplanet regimes represent a crucial discovery space for testing planet formation theories. ExELS will use the gravitational microlensing technique and will detect 1000 microlensing events per month over 1.6 deg^2 of the Galactic bulge. We assess how many of these events will have detectable planetary signatures using a detailed multi-wavelength microlensing simulator (MABuLS) which incorporates the Besancon Galactic model with 3D extinction. MABuLS is the first theoretical simulation of microlensing to treat the effects of point spread function (PSF) blending self-consistently with the underlying Galactic model. We use MABuLS, together with current numerical models for the Euclid PSFs, to explore a number of designs and de-scope options for ExELS, including the exoplanet yield as a function of filter choice and slewing time, and the effect of systematic photometry errors. Using conservative extrapolations of current empirical exoplanet mass functions determined from ground-based microlensing and radial velocity surveys, ExELS can expect to detect a few hundred cold exoplanets around mainly G, K and M-type stellar hosts, including ~45 Earth-mass planets and ~6 Mars-mass planets for an observing programme totalling 10 months. ExELS will be capable of measuring the cold exoplanet mass function down to Earth mass or below, with orbital separations from ~1AU to the free-floating regime. (Abridged)

preprint2013arXiv

Kinematics and chemical properties of the Galactic stellar populations. The HARPS FGK dwarfs sample

(Abridged) We analyze chemical and kinematical properties of about 850 FGK solar neighborhood long-lived dwarfs observed with the HARPS high-resolution spectrograph. The stars in the sample have logg > 4 dex, 5000 < Teff < 6500 K, and -1.39 < [Fe/H] < 0.55 dex. We apply a purely chemical analysis approach based on the [alpha/Fe] vs. [Fe/H] plot to separate Galactic stellar populations into the thin disk, thick disk and high-alpha metal-rich (hamr). Our analysis shows a negative gradient of the rotational velocity of the thin disk stars with [Fe/H] (-17 km s^-1 dex^-1), and a steep positive gradient for both the thick disk and hamr stars with the same magnitude of about +42 km s^-1 dex^-1. For the thin disk stars we observed no correlation between orbital eccentricities and metallicity, but observed a steep negative gradient for the thick disk and hamr stars with practically the same magnitude (about -0.18 dex^-1). Our results suggest that radial migration played an important role in the formation and evolution of the thin disk. For the thick disk stars it is not possible to reach a firm conclusion about their origin. Based on the eccentricity distribution of the thick disk stars only their accretion origin can be ruled out, and the heating and migration scenario could explain the positive steep gradient of V_phi with [Fe/H]. Analyzing the hamr stellar population we found that they share properties of both the thin and thick disk population. A comparison of the properties of the hamr stars with that of the subsample of stars from the N-body/SPH simulation using radial migration suggest that they may have originated from the inner Galaxy. Further detailed investigations would help to clarify their exact nature and origin.

preprint2013arXiv

The Effective Temperature Scale of M dwarfs

Despite their large number in the Galaxy, M dwarfs remain elusive objects and the modeling of their photospheres has long remained a challenge (molecular opacities, dust cloud formation). Our objectives are to validate the BT-Settl model atmospheres, update the M dwarf T eff -spectral type relation, and find the atmospheric parameters of the stars in our sample. We compare two samples of optical spectra covering the whole M dwarf sequence with the most recent BT-Settl synthetic spectra and use a \c{hi}2 minimization technique to determine Teff . The first sample consists of 97 low-resolution spectra obtained with NTT at La Silla Observatory. The second sample contains 55 mid-resolution spectra obtained at the Siding Spring Observatory (SSO). The spectral typing is realized by comparison with already classified M dwarfs. We show that the BT-Settl synthetic spectra reproduce the slope of the spectral energy distribution and most of its features. Only the CaOH band at 5570Å and AlH and NaH hydrides in the blue part of the spectra are still missing in the models. The Teff-scale obtained with the higher resolved SSO 2.3 m spectra is consistent with that obtained with the NTT spectra. We compare our Teff scale with those of other authors and to published isochrones using the BT-Settl colors. We also present relations between effective temperature, spectral type and colors of the M dwarfs.

preprint2012arXiv

Gaia Universe Model Snapshot : A statistical analysis of the expected contents of the Gaia catalogue

Context. This study has been developed in the framework of the computational simulations executed for the preparation of the ESA Gaia astrometric mission. Aims. We focus on describing the objects and characteristics that Gaia will potentially observe without taking into consideration instrumental effects (detection efficiency, observing errors). Methods. The theoretical Universe Model prepared for the Gaia simulation has been statistically analyzed at a given time. Ingredients of the model are described, giving most attention to the stellar content, the double and multiple stars, and variability. Results. In this simulation the errors have not been included yet. Hence we estimate the number of objects and their theoretical photometric, astrometric and spectroscopic characteristics in the case that they are perfectly detected. We show that Gaia will be able to potentially observe 1.1 billion of stars (single or part of multiple star systems) of which about 2% are variable stars, 3% have one or two exoplanets. At the extragalactic level, observations will be potentially composed by several millions of galaxies, half million to 1 million of quasars and about 50,000 supernovas that will occur during the 5 years of mission. The simulated catalogue will be made publicly available by the DPAC on the Gaia portal of the ESA web site http://www.rssd.esa.int/gaia/.

preprint2012arXiv

Three dimensional interstellar extinction map towards the Galactic Bulge

Studies of the properties of the inner Galactic Bulge depend strongly on the assumptions about the interstellar extinction. Most of the extinction maps available in the literature lack the information about the distance. We combine the observations with the Besancon model of the Galaxy to investigate the variations of extinction along different lines of sight towards the inner Galactic bulge as a function of distance. In addition we study the variations in the extinction law in the Bulge. We construct color-magnitude diagrams with the following sets of colors: H-Ks and J-Ks from the VVV catalogue as well as Ks-[3.6], Ks-[4.5], Ks-[5.8] and Ks-[8.0] from GLIMPSE-II catalogue matched with 2MASS. Using the newly derived temperature-color relation for M giants that match better the observed color-magnitude diagrams we then use the distance-color relations to derive the extinction as a function of distance. The observed colors are shifted to match the intrinsic colors in the Besançon model, as a function of distance, iteratively thereby creating an extinction map with three dimensions: two spatial and one distance dimension along each line of sight towards the bulge. Colour excess maps are presented at a resolution of 15' x 15' for 6 different combinations of colors, in distance bins of 1 kpc. The high resolution and depth of the photometry allows to derive extinction maps to 10 kpc distance and up to 35 magnitudes of extinction in Av (3.5 mag in Aks). Integrated maps show the same dust features and consistent values with the other 2D maps. Starting from the color excess between the observations and the model we investigate the extinction law in near-infrared and its variation along different lines of sight.

preprint2011arXiv

Selection constraints on high redshift quasar searches in the VISTA kilo-degree infrared galaxy survey

The European Southern Observatory's (ESO) Visible and Infrared Survey Telescope for Astronomy (VISTA) is a 4-m class survey telescope for wide-field near-infrared imaging. VISTA is currently running a suite of six public surveys, which will shortly deliver their first Europe wide public data releases to ESO. The VISTA Kilo-degree Infrared Galaxy Survey (VIKING) forms a natural intermediate between current wide shallow, and deeper more concentrated surveys, by targeting two patches totalling 1500 sq.deg in the northern and southern hemispheres with measured 5-sigma limiting depths of Z ~ 22.4, Y ~ 21.4, J ~ 20.9, H ~ 19.9 and Ks ~19.3 (Vega). This architecture forms an ideal working parameter space for the discovery of a significant sample of 6.5 <= z <= 7.5 quasars. In the first data release priority has been placed on small areas encompassing a number of fields well sampled at many wavelengths, thereby optimising science gains and synergy whilst ensuring a timely release of the first products. For rare object searches e.g. high-z quasars, this policy is not ideal since photometric selection strategies generally evolve considerably with the acquisition of data. Without a reasonably representative data set sampling many directions on the sky it is not clear how a rare object search can be conducted in a highly complete and efficient manner. In this paper, we alleviate this problem by supplementing initial data with a realistic model of the spatial, luminosity and colour distributions of sources known to heavily contaminate photometric quasar selection spaces, namely dwarf stars of spectral type M, L and T. We use this model along with a subset of available data to investigate contamination of quasar selection space by cool stars and galaxies and lay down a set of benchmark selection constraints that limit contamination to reasonable levels whilst maintaining high completeness...

preprint2010arXiv

GYES, a multifibre spectrograph for the CFHT

We have chosen the name of GYES, one of the mythological giants with one hundred arms, offspring of Gaia and Uranus, for our instrument study of a multifibre spectrograph for the prime focus of the Canada-France-Hawaii Telescope. Such an instrument could provide an excellent ground-based complement for the Gaia mission and a northern complement to the HERMES project on the AAT. The CFHT is well known for providing a stable prime focus environment, with a large field of view, which has hosted several imaging instruments, but has never hosted a multifibre spectrograph. Building upon the experience gained at GEPI with FLAMES-Giraffe and X-Shooter, we are investigating the feasibility of a high multiplex spectrograph (about 500 fibres) over a field of view 1 degree in diameter. We are investigating an instrument with resolution in the range 15000 to 30000, which should provide accurate chemical abundances for stars down to 16th magnitude and radial velocities, accurate to 1 km/s for fainter stars. The study is led by GEPI-Observatoire de Paris with a contribution from Oxford for the study of the positioner. The financing for the study comes from INSU CSAA and Observatoire de Paris. The conceptual study will be delivered to CFHT for review by October 1st 2010.

preprint2008arXiv

High precision microlensing maps of the Galactic bulge

We present detailed maps of the microlensing optical depth and event density over an area of 195 sq. deg towards the Galactic bulge. The maps are computed from synthetic stellar catalogues generated from the Besancon Galaxy Model, which comprises four stellar populations and a three-dimensional extinction map calibrated against the Two-Micron All-Sky Survey. The optical depth maps have a resolution of 15 arcminutes, corresponding to the angular resolution of the extinction map. We compute optical depth and event density maps for all resolved sources above I=19, for unresolved (difference image) sources magnified above this limit, and for bright standard candle sources in the bulge. We show that the resulting optical depth contours are dominated by extinction effects, exhibiting fine structure in stark contrast to previous theoretical optical depth maps. Optical depth comparisons between Galactic models and optical microlensing survey measurements cannot safely ignore extinction or assume it to be smooth. We show how the event distribution for hypothetical J and K-band microlensing surveys, using existing ground-based facilities such as VISTA, UKIRT or CFHT, would be much less affected by extinction, especially in the K band. The near infrared provides a substantial sensitivity increase over current I-band surveys and a more faithful tracer of the underlying stellar distribution, something which upcoming variability surveys such as VVV will be able to exploit. Synthetic population models offer a promising way forward to fully exploit large microlensing datasets for Galactic structure studies.

preprint2006arXiv

Modelling the Galactic Interstellar Extinction Distribution in Three Dimensions

The Two Micron All Sky Survey, along with the Stellar Population Synthesis Model of the Galaxy, developed in Besancon, is used to calculate the extinction distribution along different lines of sight. By combining many lines of sight, the large scale distribution of interstellar material can be deduced. The Galaxy model is used to provide the intrinsic colour of stars and their probable distances, so that the near infrared colour excess, and hence the extinction, may be calculated and its distance evaluated. Such a technique is dependent on the model used, however we are able to show that moderate changes in the model parameters result in insignificant changes in the predicted extinction. This technique has now been applied to over 64000 lines of sight, each separated by 15 arcmin, in the inner Galaxy (|l|<=100 deg, |b|<=10 deg). Using our extinction map, we have derived the main characteristics of the large scale structure of the dust distribution: scale height and warp of the ISM disc as well as the angle of the dust in the Galactic Bar. This resulting extinction map will be useful for studies of the inner Galaxy and its stellar populations.

preprint2006arXiv

The stellar content of the COSMOS field as derived from morphological and SED based gtar/galaxy Separation

We report on the stellar content of the COSMOS two degree field, as derived from a rigorous star/galaxy separation approach developed for using stellar sources to define the point spread function variation map used in a study of weak galaxy lensing. The catalog obtained in one filter from the ACS (Advanced Camera for Survey on the Hubble Space Telescope) is cross-identified with ground based multi-wavelength catalogs. The classification is reliable to magnitude $F_{814W}=24$ and the sample is complete even fainter. We construct a color-magnitude diagram and color histograms and compare them with predictions of a standard model of population synthesis. We find features corresponding to the halo subdwarf main sequence turnoff, the thick disk, and the thin disk. This data set provides constraints on the thick disk and spheroid density laws and on the IMF at low mass. We find no evidence of a sharp spheroid edge out to this distance. We identify a blue population of white dwarfs with counts that agree with model predictions. We find a hint for a possible slight stellar overdensity at about 22-34 kpc but the data are not strong enough at present to claim detection of a stream feature in the halo (abridged).

preprint2005arXiv

Stellar populations in the CFHTLS. I. New constraints on the IMF at low mass

We present a stellar populations analysis of the first release of the CFHTLS (Canada-France-Hawai Telescope Legacy Survey) data. A detailed comparison between the Besançon model of the Galaxy and the first data release of the CFHTLS-Deep survey is performed by implementing the MEGACAM photometric system in this model using stellar atmosphere model libraries. The reliability of the theoretical libraries to reproduce the observed colours in the MEGACAM system is investigated. The locations of various stellar species like subdwarfs, white dwarfs, late-type and brown dwarfs, binary systems are identified. The contamination of the stellar sample by quasars and compact galaxies is quantified using spectroscopic data from the VIMOS-VLT Deep Survey (VVDS) as a function of $i'$ magnitude and $r'-i'$ colour. A comparison between simulated counts using the standard IMF at low masses show that the number of very low mass dwarfs may have been underestimated in previous studies. These observations favour a power law IMF following $d(n)/dm \propto m^{-α}$ with $α=2.5$ for $m < $ 0.25 \Msun or $α=3.0$ for $m < $ 0.2 \Msun for single stars. The resulting LF is in agreement with the local LF as measured from the 5 or 25 pc samples. It is in strong disagreement with the Zheng et al (2001) LF measured from deep HST data. We show that this discrepancy can be understood as an indication of a different IMF at low masses at early epochs of the Galaxy compared to the local thin disc IMF.