Issue |
A&A
Volume 566, June 2014
|
|
---|---|---|
Article Number | A119 | |
Number of page(s) | 15 | |
Section | Catalogs and data | |
DOI | https://doi.org/10.1051/0004-6361/201423636 | |
Published online | 24 June 2014 |
Overview and stellar statistics of the expected Gaia Catalogue using the Gaia Object Generator
1
Dept. d’Astronomia i Meteorologia, Institut de Ciències del Cosmos,
Universitat de Barcelona (IEEC-UB),
Martí Franquès 1,
08028
Barcelona,
Spain
e-mail:
mpalmer@am.ub.es
2
GEPI, Observatoire de Paris, CNRS, Université Paris
Diderot, 5 place Jules
Janssen, 92190
Meudon,
France
3
Scientific Support Office of the Directorate of Science and
Robotic Exploration of the European Space Agency, European Space Research and
Technology Centre, Keplerlaan 1, 2201 AZ
Noordwijk, The
Netherlands
4
Institut UTINAM, CNRS UMR 6213, Observatoire des Sciences de
l’Univers THETA Franche-Comté-Bourgogne, Université de Franche-Comté, Observatoire de
Besançon, BP 1615,
25010
Besançon Cedex,
France
Received: 13 February 2014
Accepted: 2 April 2014
Aims. An effort has been made to simulate the expected Gaia Catalogue, including the effect of observational errors. We statistically analyse this simulated Gaia data to better understand what can be obtained from the Gaia astrometric mission. This catalogue is used 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 previous work, where the expected Gaia Catalogue content was reviewed but without the simulation of observational errors.
Methods. We analysed the Gaia Object Generator (GOG) catalogue using the Gaia Analysis Tool (GAT), thereby producing a number of statistics about 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, and the photometry in the four Gaia bands. Information is also given on the expected performance of physical parameter determination, including temperature, metallicity, and line-of-sight extinction.
Key words: stars: statistics / galaxies: statistics / Galaxy: stellar content / methods: data analysis / astrometry / catalogs
© ESO, 2014
1. Introduction
Gaia is a cornerstone ESA mission, launched in December 2013, and will produce the fullest 3D Galactic census to date. The mission is expected to yield a huge advancement in our understanding of the composition, structure, and evolution of the Galaxy (Perryman et al. 2001). Through Gaia’s photometric instruments, object detection up to G = 20 mag will be possible (see Jordi et al. (2010) for a definition of G magnitude), including measurements of positions, proper motions, and parallaxes up to micro arcsecond accuracy. The on-board radial velocity spectrometer will provide radial velocity measurements for stars down to a limit of GRVS = 17 mag. With low-resolution spectra providing information on effective temperature, line-of-sight extinction, surface gravity, and chemical composition, Gaia will yield a detailed catalogue that contains roughly 1% of the entire Galactic stellar population.
Gaia will represent a huge advance on its predecessor, Hipparcos (Perryman & ESA 1997), both in terms of massive increases in precision and in the numbers of objects observed. Thanks to accurate observations of large numbers of stars of all kinds, including rare objects, large numbers of Cepheids and other variable stars, and direct parallax measurements for stars in all Galactic populations (thin disk, thick disk, halo, and bulge), Gaia data is expected to have a strong impact on luminosity calibration and improvement of the distance scale. This, along with applications to studies of Galactic dynamics and evolution and of fields ranging from exoplanets to general relativity and cosmology, Gaia’s impact is expected to be significant and far reaching.
During its five years of data collection, Gaia is expected to transmit some 150 terabytes of raw data to Earth, leading to production of a catalogue of 109 individual objects. After on-ground processing, the full database is expected to be in the range of one to two petabytes of data. Preparation for acquiring this huge amount of data is essential. Work has begun to model the expected output of Gaia in order to predict the content of the Gaia Catalogue, to facilitate the production of tools required to effectively validate the real data before publication, and to analyse the real data set at the end of the mission.
To this end, the Gaia Data Processing and Analysis Consortium (DPAC) has been preparing a set of simulators, including a simulator called the Gaia Object Generator (GOG), which simulates the end-of-mission catalogue, including observational errors. Here a full description of GOG is provided, including the models assumed for the performance of the Gaia satellite and an overview of its simulated end-of-mission catalogue. A selection of statistics from this catalogue is provided to give an idea of the performance and output of Gaia.
In Sect. 2, a brief description of the Gaia instrument and an overview of the current simulation effort is given, followed by definitions of the error models assumed for the performance of the Gaia satellite in Sect. 3. In Sect. 4, the methods used for searching the simulated catalogue and generating statistics are described. In Sect. 5, we present the results of the full sky simulation, broken up into sections for each parameter in the catalogue and specific object types of interest. Finally, in Sect. 6 we provide a summery and conclusions.
2. The Gaia satellite
The Gaia space astrometry mission will map the entire sky in the visible G band over the course of its five-year mission. Located at Lagrangian point L2, Gaia will be constantly and smoothly spinning. It has two telescopes separated by a basic angle of 106.5°. Light from stars that are observed in either telescope is collected and reflected to transit across the Gaia focal plane.
The Gaia focal plane can be split into several main components. The majority of the area is taken up by CCDs for the broad band G magnitude measurements in white light, used in taking the astrometry measurements. Second, there is a pair of low-resolution spectral photometers, one red and one blue, producing low-resolution spectra with integrated magnitude GRP and GBP, respectively. Finally, there is a radial velocity spectrometer observing at near-infrared, with integrated magnitude GRVS. The magnitudes G, GRP, and GBP will be measured for all Gaia sources (G ≤ 20), whereas GRVS will be measured for objects up to GRVS ≤ 17 magnitude. For an exact definition of the Gaia focal plane and the four Gaia bands, see Figs. 1 and 3 of Jordi et al. (2010).
The motion of Gaia is complex, with rotations on its own spin axis occurring every six hours. This spin axis is itself precessing, and is held at a constant 45° degrees from the Sun. From its position at L2, Gaia will orbit the Sun over the course of a year. Thanks to the combination of these rotations, the entire sky will be observed repeatedly. The Gaia scanning law gives the number of times a region will be re-observed by Gaia over its five-year mission, and comes from this spinning motion of the satellite and its orbit around the Sun. Objects in regions with more observations have greater precision, while regions with fewer repeat observations have lower precision. The average number of observations per object is 70, although it can be as low as a few tens or as high as 200.
3. Simulations
Simulation of many aspects of the Gaia mission has been carried out in order to test and improve instrument design, data reduction algorithms, and tools for using the final Gaia Catalogue data. The Gaia Simulator is a collection of three data generators designed for this task: the Gaia Instrument and Basic Image Simulator (GIBIS, Babusiaux 2005), the Gaia System Simulator (GASS, Masana et al. 2010), and the Gaia Object Generator (GOG, described here). Through these three packages, the production of the simulated Gaia telemetry stream, observation images down to pixel level and intermediate or final catalogue data is possible.
3.1. Gaia Universe Model snapshot and the Besançon galaxy model
One basic component of the Gaia Simulator is its Universe Model (UM), which is used to create object catalogues down to a particular limiting magnitude (in our case G = 20 mag for Gaia). For stellar sources, the UM is based on the Besançon galaxy model (Robin et al. 2003). This model simulates the stellar content of the Galaxy, including stellar distribution and a number of object properties. It produces stellar objects based on the four main stellar populations (thin disk, thick disk, halo, and bulge), each population with its own star formation history and stellar evolutionary models. Additionally, a number of object-specific properties are also assigned to each object, dependent on its type. Possible objects are stars (single and multiple), nebulae, stellar clusters, diffuse light, planets, satellites, asteroids, comets, resolved galaxies, unresolved extended galaxies, quasars, AGN, and supernovae. Therefore, the UM is capable of simulating almost every object type that Gaia can potentially observe. It can therefore construct simulated object catalogues down to Gaia’s limiting magnitude.
Building on this, the UM creates for any time, over any section of the sky (or the whole sky), a set of objects with positions and assigns each a set of observational properties (Robin et al. 2012). These properties include distances, apparent magnitudes, spectral characteristics, and kinematics.
Clearly the models and probability distributions used in order to create the objects with their positions and properties are highly important in producing a realistic catalogue. The UM has been designed so that the objects it creates fit as closely as possible to observed statistics and to the latest theoretical formation and evolution models (Robin et al. 2003). For a statistical analysis of the underlying potentially observable population (with G ≤ 20 mag) using the Gaia UM without satellite instrument specifications and error models, see Robin et al. (2012).
3.2. The Gaia Object Generator (GOG)
The GOG is capable of transforming this UM catalogue into Gaia’s simulated intermediate and final catalogue data. This is achieved through the use of analytical and numerical error models to create realistic observational errors in astrometric, photometric, and spectroscopic parameters (Isasi et al. 2010). In this way, GOG transforms “true” object properties from the UM into “observed” quantities that have an associated error that depends on the object’s properties, Gaia’s instrument capabilities, and the type and number of observations made.
3.3. Error models
DPAC is divided into a number of coordination units (CUs), each of which specialises in a specific area. In GOG we have taken the recommendations from the various CUs in order to include the most complete picture of Gaia performance as possible. The CUs are divided into the following areas: CU1, system architecture; CU2, simulations; CU3, core processing (astrometry); CU4, object processing (multiple stars, exoplanets, solar system objects, extended objects); CU5, photometric processing; CU6, spectroscopic processing; CU7, variability processing; CU8, astrophysical parameters; and CU9, catalogue access.
Models for specific parameters have been provided by the various CUs, and only an outline is given here. In the following description, true refers to UM data (without errors), epoch to simulated individual observations (including errors), and observed to the simulated observed data for the end of the Gaia mission (including standard errors). Throughout, error refers to the formal standard error on a measurement.
3.3.1. Astrometric error model
The formal error on the parallax, σϖ, is calculated following the expression: (1)
-
ση is the CCD centroid positioning error. It uses the Cramer Rao (CR) lower bound in its discrete form, which defines the best possible precision of the maximum likelihood centroid location estimator. The CR lower bound requires the line spread function (LSF) derivative for each sample1, the background, the readout noise, and the source integrated signal.
-
m is the contingency margin that is used to take scientific and environmental effects into account, for example: uncertainties in the on-ground processing, such as uncertainties in relativistic corrections and solar system ephemeris; effects such as having an imperfect calibrating LSF; errors in estimating the sky background; and other effects when dealing with real stars. The default value assumed for the Gaia mission has been set by ESA as 1.2 and is used in GOG.
-
gϖ = 1.47 / sin ξ is a geometrical factor where ξ is known as the solar aspect angle, with a value of 45.
-
Neff is the number of elementary CCD transits (Nstrip × Ntransit).
-
Ntransit is the number of field of view transits.
-
Nstrip is the number of CCDs in a row on the Gaia focal plane. It has a value of 9, except for the row that includes the wavefront sensor, which has 8 CCDs.
-
σcal is the calibration noise. A constant value of 5.7 μas has been applied. This takes into account that the end-of-mission precision on the astrometric parameters not only depends on the error due to the location estimation with each CCD. There are calibration errors from the CCD calibrations, the uncertainty of the attitude of the satellite and the uncertainty on the basic angle.
We are enabling the activation of gates, as described in the Gaia Parameter Database. The Gaia satellite will be smoothly rotating and will constantly image the sky by collecting the photons from each source as they pass along the focal plane. The total time for a source to pass along the focal plane will be 107 seconds, and the electrons accumulated in the a CCD pixel will be passed along the CCD at the same rate as the source. To avoid saturation for brighter sources, and the resulting loss of astrometric precision, gates can be activated that limit the exposure time. Here we are using the default gating system, which could change during the mission.
Following de Bruijne (2012)2, we have assumed that the errors on the positional coordinates at the mean epoch (middle of the mission), and the error in the proper motion coordinates can be given respectively by
-
σα = 0.787σϖ
-
σδ = 0.699σϖ
-
σμα = 0.556σϖ
-
σμδ = 0.496σϖ.
3.3.2. Photometric error model
GOG uses the single CCD transit photometry error σp,j (Jordi et al. 2010) defined as (2)to compute either the epoch errors or the end-of-mission errors.
We assume, following an “aperture photometry” approach, that the object flux sj is measured in a rectangular “aperture” (window) of ns along-scan object samples. The sky background bj is assumed to be measured in nb background samples, and r is CCD readout noise. The faperture·sj is expressed in units of photo-electrons (e−), and denotes the object flux in photometric band j contained in the “aperture” (window) of ns samples, after a CCD crossing. The number faperture thus represents the fraction of the object flux measured in the aperture window.
For the epoch and end-of-mission data, the same expression is used for the standard deviation calculation (3)where Neff is the number of elementary CCD transits (Nstrip × Ntransit), with (Nstrip = 1 and Ntransit = 1) for epoch photometry and equal to the number of transits for the end-of-mission photometry. The calibration noise σcal has a fixed value of σcal = 30 mmag. Use of a fixed value of the centroiding error is possible because this error is only relevant for brighter stars, because their centroiding errors are smaller than the calibration error.
3.3.3. Radial velocity spectrometer errors
CU6 tables (see Table 1) using the stars’ physical parameters and apparent magnitude are used to obtain σVr. They were computed following the prescriptions of Katz et al. (2004) and later updates. Those tables have been provided for one and 40 field-of-view transits, therefore the value for 40 transits is used here to calculate the average end-of-mission errors in RVS.
Given the information on the apparent Johnson V magnitude and the atmospheric parameters of each star (from the UM), we select from Table 1 the closest spectral type and return the corresponding radial velocity error. Since Table 1 is given for [ Fe / H ] = 0 alone, a correction is made on the apparent magnitude in order to take different metallicities into account: for each variation in metallicity of Δ [ Fe / H ] = 1.5 dex, the magnitude is increased by V = 0.5 mag.
We have set a lower limit on the wavelength calibration error, giving a lower limit on the radial velocity error of 1 km s-1. For the faintest stars the spectra will be of poor quality and will not contain enough information to enable accurate estimation of the radial velocity. Owing to the limited bandwidth in the downlink of Gaia data to Earth, poor quality spectra will not be transmitted. We therefore set an upper limit on the radial velocity error of 20 km s-1, beyond which we assume that there will be no data. The exact point at which the data will be assumed to have too low a quality is still unknown.
Average end-of-mission formal error in radial velocity with an assumed average of 40 field-of-view transits, in km s-1, for each spectral type.
3.3.4. Physical parameters
GOG uses the stellar parametrisation performance given by CU8 to calculate error estimations for effective temperature, line-of-sight extinction, metallicity, and surface gravity. The colour-independent extinction parameter A0 is used in preference to the band specific extinctions AV or AG, because A0 is a property of the interstellar medium alone (Bailer-Jones 2011). CU8 use three different algorithms to calculate physical parameters using spectrophotometry (see Liu et al. 2012).
It should be noted that the errors calculated here are calculated only as a function of apparent magnitude. However, as described in Liu et al. (2012), there are clear dependencies on the spectral type of the star, because some star types may or may not exhibit spectral features required for parameter determination. Additionally, Liu et al. (2012) report a strong correlation between the estimation of effective temperature and extinction. This correlation is not simulated in GOG. Following the recommendation of CU8, calculating errors of physical parameters depends on apparent magnitude and is split into two cases, objects with A0< 1 mag and A0 ≥ 1 mag.
In GOG, σTeff, σA0, σFe / H and σlog g are calculated from a Gamma distribution, with shape parameter α and scale parameter θ : (4)where α and θ are obtained from the following expressions, which have been calculated to give each σ a close approximation to the CU8 algorithm results. A gamma distribution was selected for ease of implementation and for its ability to statistically recreate the CU8 results to a reasonable approximation. A gamma distribution is also only non-zero for positive values of sigma. This is essential when modelling errors because, of course, it is impossible to have a negative error.
-
For stars with A0< 1 mag:
αA0 = 0.204−0.032G + 0.001G2
αlog g = 0.151−0.019G + 0.001G2
αFe / H = 0.295−0.047G + 0.002G2
αTeff = 78.2−10.3G + 0.46G2
θA0 = 0.084
θlog g = 0.160
θFe / H = 0.121
θTeff = 28.2.
-
For stars with A0 ≥ 1 mag:
αA0 = 0.178−0.026G + 0.001G2
αlog g = 0.319−0.044G + 0.002G2
αFe / H = 0.717−0.115G + 0.005G2
αTeff = 67.3−7.85G + 0.35G2
θA0 = 0.096
θlog g = 0.179
θFe / H = 0.353
θTeff = 33.5.
The Gamma distributions thus obtained for each parameter are used to generate a formal error for each parameter for each individual star, aiming to statistically (but not individually) reproduce the results that will be obtained by the application of the CU8 algorithms and then included in the Gaia Catalogue.
It should be noted that in the stellar parametrisation algorithms used in Liu et al. (2012), a degeneracy is reported between extinction and effective temperature owing to the lack of resolved spectral lines only sensitive to effective temperature. In GOG, this degeneracy has not been taken into account, and the precisions of each of the four stellar parameters is simulated independently.
Additionally, the results of Liu et al. (2012) have recently been updated, and Bailer-Jones et al. (2013) gives the latest results regarding the capabilities of physical parameter determination. This latest paper has not been included in the current version of GOG.
3.4. Limitations
In the present paper, only the results for single stars are given in detail. All of the figures and tables represent the numbers and statistics of only individual single stars, excluding all binary and multiple systems. Since the performance of the Gaia satellite is largely unknown for binary and multiple systems, the implementation into GOG of realistic error models has not yet been possible. While the results presented in Sect. 5 are expected to be reliable under current assumptions for the performance of Gaia, the real Gaia Catalogue will differ from these results thanks to the presence of binary and multiple systems. By removing binaries from the latter, direct comparison of the results presented here with the real Gaia Catalogue will not be possible because of the presence of unresolved binaries, which are difficult to detect. As a simulator, GOG relies heavily on all inputs and assumptions supplied both from the UM or from the Gaia predicted performance and error models.
In our simulations we used an exact cut at G = 20 mag, beyond which no stars are observed. In reality, in regions of low density observations of stars up to 20.5 mag could be possible. Inversely, very crowded regions may not be complete up to 20 mag, or the numbers of observations per star over the five-year mission may be reduced in these regions.
There is no simulation of the impact of crowding on object detection or the detection of components in binary and multiple systems. This can lead to unrealistic quality in all observed data in the most crowded regions of the plane of the Galaxy, to overestimates for star counts in the bulge, and to a lack of features related to the disk and bulge in Figs. 3 and 19.
Additionally, GOG uses the nominal Gaia scanning law to calculate the number of field of view transits per object over the five years of the mission while Gaia is operating in normal mode. There will be an additional one-month period at the start of the mission using an ecliptic pole scanning law, and this has not been taken into account. It may lead to a slight underestimation of the number of transits, and therefore a slight overestimation of errors, for some stars near the ecliptic poles.
There is the possibility that the Gaia mission will be extended above the nominal five-year mission. Since this idea is under discussion and has not yet been confirmed or discarded, we only present results for the Gaia mission as originally planned.
If the length of the mission is increased, the number of field-of-view transits will increase, and the precision per object will improve. If the proposal is accepted, the GOG simulator could be used to provide updated statistics for the expected catalogue without extensive modification.
4. Methods and statistics
Considering current computing capabilities, it is not straightforward to make statistics and visualisations when dealing with catalogues of such a large size. A specific tool has been created which is capable of extracting information and visualising results, with excellent scalability allowing its use for huge datasets and distributed computing systems.
The Gaia Analysis Tool (GAT) is a data analysis package that allows, through three distinct frameworks, generation of statistics, validation of data, and generation of catalogues. It currently handles both UM- and GOG-generated data, and could be adapted to handle other data types.
Every statistical analysis is performed by a Statistical Analysis Module (SAM), with several grouped into a single XML file as an input to GAT. Each SAM can contain a set of filters, enabling analysis of specific subsets of the data. This allows the production of a wide range of statistics for objects satisfying any number of specific user-defined criteria or for the catalogue as a whole.
GAT creates a number of different statistics outputs including histograms, sky density maps and HR diagrams. After the GAT execution, statistics output are stored to either generate a report or to be analysed using the GAT Displaying tool.
Because we have information from not only observations of a population but also of the observed population itself, comparison is possible between the simulated Gaia Catalogue and the simulated “true” population, allowing large scope for investigating the precision3 of the observations and differences between the two catalogues. Clearly this is only possible with simulated data and cannot be attempted with the true Gaia Catalogue, so it is an effective way to investigate the possible extent and effect of observational errors and selection bias on the real Gaia Catalogue, where this kind of comparison is not possible.
GOG can be used in the preparations for validating the true Gaia Catalogue, by testing the GOG catalogue for accuracy3 and precision. In special cases, observational biases could even be implemented into the code to allow thorough testing of validation methods.
Fig. 1 Skymap of total integrated flux over the Milky Way, in the G band. The colour bar represents a relative scale, from maximum flux in white to minimum flux in black. The figure is plotted in Galactic coordinates with the Galactic-longitude orientation swapped left to right. |
5. Results
The GOG simulator has been used to generate the simulated final mission catalogue for Gaia, down to magnitude G = 20. The simulation was performed on the MareNostrum super computer at the Barcelona Supercomputing Centre (Centre Nacional de Supercomputació), and it took 400 thousand CPU hours. An extensive set of validations and statistics have been produced using GAT to validate performance of the simulator. Below we include a subset of these statistics for the most interesting cases to give an overview of the expected Gaia Catalogue.
5.1. General
In total, GOG has produced a catalogue of about one billion objects, consisting of 523 million individual single stars and 484 million binary or multiple systems. The total number of stars, including the components of binary and multiple systems is 1.6 billion. The skymap of the total flux detected over the entire five-year mission is given in Fig. 1. Although GOG can produce extragalactic sources, none have been simulated here.
The following discussion is split into sections for different types of objects of interest. Section 5.2 covers all Galactic stellar sources. Section 5.3 covers all variable objects. Section 5.4 is a discussion of physical parameters estimated by Gaia. All objects in these sections are within the Milky Way.
To make the presentation of performance as clear as possible, all binary and multiple systems and their components have been removed from the following statistics. This is due to complicating effects that arise when dealing with binary systems, some of which GOG is not yet capable of correctly simulating; for example, GOG does not yet contain an orbital solution in its astrometric error models, and the effects of unresolved systems on photometry and astrophysical parameter determination have not yet been well determined. Therefore, the numbers presented are only for individual single stars and do not include the full one billion objects simulated. Of the single stars presented in this paper, 74 million are within the radial velocity spectrometer magnitude range.
Table 2 gives the mean and median error for each of the observed parameters discussed in this paper, along with the upper and lower 25% quartile.
In Fig. 2, the mean error for parallax, position, proper motion, and photometry in the four Gaia bands are given as a function of G magnitude. Also, the mean error in radial velocity is given as a function of GRVS magnitude. The sharp jumps in the mean error in astrometric parameters between 8 and 12 mag are due to the activation of gates for the brighter sources in an attempt to prevent CCD saturation (see Sect. 3.3.1).
Mean and median value of the end-of-mission error in each observable, along with the upper (UQ) and lower (LQ) 25% quartile.
Fig. 2 Mean end-of-mission error as a function of G magnitude for parallax, position, proper motion, and photometry in the four Gaia passbands. Additionally the mean end-of-mission error in radial velocity as a function of GRVS magnitude. |
5.2. Stars
5.2.1. Parallax
The distribution of parallax measurements for all single stars is given in Fig. 4. The mean parallax error for all single stars is 147.0 μas. The number of single stars falling below three relative parallax error limits is given in Table 3 for each spectral type and in Table 4 for each luminosity class. For those interested in a specific type of star, Table 5 gives the full breakdown of the number of single stars falling below three relative parallax error limits for every spectral type and luminosity class. The distribution of parallax errors is given for each stellar population in Fig. 5 and for each spectral type in Fig. 7. The relative parallax error σϖ/ϖ is given in Fig. 6 for stars split by spectral type.
Fig. 3 Sky map (healpix) of mean parallax error for all single stars in equatorial coordinates. Colour scale is mean parallax error in μas. The red area is the location of the bulge. |
Fig. 4 Histogram of parallax for all single stars. The histogram contains 99.5% of all data. |
The error in parallax measurements for Gaia depends on the magnitude of the source, the number of observations made, and the true value of the parallax. Figure 3 shows the mean parallax error over the sky. Its shape clearly follows that of the Gaia scanning law. The red area corresponding to the region of worst precision is due to the bulge population, which suffers from high levels of reddening. The faint ring around the centre of the figure corresponds to the disk of the Galaxy, remembering that the plot is given in equatorial coordinates. The blue areas corresponding to regions of improved mean precision are areas with a higher number of observations. The characteristic shape of this plot is due to the Gaia scanning law. The error in parallax as a function of measured G magnitude is given in Fig. 8 and as a function of the real parallax in Fig. 9.
Total number of single stars for each spectral type, along with the percentage of those that fall below each relative parallax error limit; e.g., 68% of M-type stars have a relative parallax error better than 20%.
Total number of single stars for each luminosity class, along with the percentage that fall below each relative parallax error limit.
Total number of single stars for each stellar classification, along with the percentage that fall below each relative parallax error limit.
Fig. 5 Histogram of end-of-mission parallax error for all single stars, split by stellar population. |
Fig. 6 Cumulative histogram of relative parallax error for all single stars, split by spectral type. The histogram range displays 74% of all data. |
Fig. 7 Histogram of end-of-mission parallax error for all single stars, split by spectral type. |
Fig. 8 End-of-mission parallax error against G magnitude for all single stars. The colour scale represents the log of density of objects in a bin size of 80 mmag by 2.5 μas. White area represents zero stars. |
Fig. 9 End-of-mission parallax error against parallax for all single stars. The colour scale represents the log of density of objects in a bin size of 10 by 2.5 μas. White area represents zero stars. |
5.2.2. Position
Gaia will be capable of measuring the position of each observed star at an unprecedented accuracy, producing the most precise full sky position catalogue to date.
The mean error is 90 μas for right ascension and 103 μas for declination. The distribution of error in right ascension and declination as a function of the true value, along with a histogram of the error, are given in Figs. 10 and 11. The overdensities are due to the bulge of the Galaxy.
Fig. 10 Right ascension error against real right ascension. The colour scale is linear, with a factor of 105. Histograms are computed for both right ascension and right ascension error. The colour scale represents log density of objects in a bin size of 2 degrees by 7.5 μas. White area represents zero stars. |
Fig. 11 Declination error against real declination. The colour scale is linear, with a factor of 105. Histograms are computed for both declination and declination error. The colour scale represents log density of objects in a bin size of 1 degrees by 5 μas. White area represents zero stars. |
5.2.3. Proper motion and radial velocity
In addition to parallax measurements, Gaia will also measure proper motions for all stars it detects. The proper motion in right ascension and declination is labelled as μα and μδ, respectively. The mean error in μα is 81.7 μas yr-1, and in μδ is 72.9 μas yr-1. The distribution of errors in both components of proper motion is given in Fig. 12.
The radial velocity is measured by the on-board radial velocity spectrometer. This instrument is only sensitive to stars down to GRVS = 17 mag. We assume an upper limit on the error in radial velocity of 20 km s-1, and assume that stars with a precision worse than this will not be given any radial velocity information.
Of the 523 million measured individual Milky Way stars, 74 million have a radial velocity measurement. The mean error in the radial velocity measurement is 8.0 km s-1. The distribution of radial velocity error is given for each G magnitude in Fig. 13, and in Fig. 14 split by spectral type. The radial velocity error is given as a function of GRVS magnitude in Fig. 16.
Fig. 12 Error in proper motion for alpha and delta for all single stars. |
Fig. 13 Histogram of radial velocity error split by G magnitude range. The histogram contains 100% of all data that have radial velocity information. |
Fig. 14 Histogram of radial velocity error split by spectral type. The histogram contains 100% of all data that have radial velocity information. |
Fig. 15 2D histograms showing error in proper motion against G magnitude. The colour scale represents the log density of objects in a bin size of 80 mmag by 2 μas yr-1. Left is proper motion in right ascension, and right is proper motion in declination. White area represents zero stars. |
5.2.4. Photometry
The end-of-mission error in each measurement as a function of G magnitude is given in Fig. 17.
Gaia will produce low-resolution spectra, in addition to measuring the magnitude of each source in the Gaia bands G, GBP, GRP, and GRVS. Whilst GOG is capable of simulating these spectra, they have not been included in the present simulations owing to the long computation time and the large storage space requirement of a catalogue of spectra for one billion sources.
Figure 18 shows the distribution in the error of each photometric measurement. As can be seen in this figure, the error in G is much lower than for the other instruments, and for all stars it is less than 8 mmag. The mean error in G is 3.0 mmag. The mean error in GBP and GRP is 14.6 mmag and 7.7 mmag, respectively. The mean error in GRVS is 13.2 mmag, although it must be remembered that the radial velocity spectroscopy instrument is limited to brighter than GRVS = 17.
Figure 19 shows the mean photometric error as a function of position on the sky for the four Gaia photometric passbands. The structure seen in all four maps is derived from the Gaia scanning law.
It is interesting to point out the ring in the four plots of Fig. 19 caused by the disk of the Galaxy. Owing to significant levels of interstellar dust in the disk of the Galaxy, visible objects are generally much redder. This reddening causes objects to lose flux at the bluer end of the spectrum, making them appear fainter to the GBP photometer. Therefore the plane of the Galaxy can be seen as an increase in the mean photometric error in the GBP error map.
Conversely, the disk of the Galaxy shows as a ring of decreased mean photometric error in the GRP and GRVS maps, since the sensitivity of their spectra is skewed more towards the redder end of the spectrum. It is important to note, however, that the effect of crowding on photometry is not accounted for in GOG.
Fig. 16 End-of-mission error in radial velocity against GRVS magnitude. The colour scale represents the log density in a bin size of 50 mmag by 1 km s-1. White area represents zero stars. |
Fig. 17 End-of-mission errors in photometry as a function of G magnitude. The colour scale represents the log density of objects in a bin size of 80 mmag by 0.4 mmag. Top left, G magnitude; top right, GBP; bottom left, GRP; bottom right, GRVS. White area represents zero stars. |
5.3. Variables
Gaia will be continuously imaging the sky over its full five-year mission, and each individual object will be observed 70 times on average. The scanning law means that the time between repeated observations varies, and Gaia will be incredibly useful for detecting many types of variable stars. GOG produces a total of 10.8 million single variable objects. This number comes from the UM (Robin et al. 2012) and assumes 100% variability detection. The exact detection rates and the classification accuracy for each variability type are still unknown. In fact, the numbers of variable objects in the catalogue is expected to be higher than 10.8 million because some variable star types have not yet been implemented (see Robin et al. 2012 for a more detailed description).
The distribution of relative parallax error is given for each type of variable star in Fig. 20. The numbers of each type of variable produced by GOG are given in Table 6, along with the number of each type that falls below each relative parallax error limit.
In general, the numbers of variables presented in this paper are lower than in Robin et al. (2012) by a factor of two or three. This is expected, because in the present paper we are excluding all variables that are part of binary or multiple systems, and presenting the number of single variable stars alone.
However, the number of emission variables is higher in the present paper. This is due to implementation of new types of emission stars: Oe, Ae, dMe, and WR stars. These are now included as emission variables but were not simulated in Robin et al. (2012). Additionally, the number of Mira variable stars is higher in the present paper. This is from an implementation error in the version of the UM used in Robin et al. (2012), which has been fixed in the version used in the present paper.
5.3.1. Cepheids and RR-Lyrae
Cepheids and RR-Lyrae are types of pulsating variable stars. Their regular pulsation and a tight period–luminosity relation make them excellent standard candles, and therefore of particular interest in studies of Galactic structure and the distance scale. Figure 21 shows the histogram of error in parallax specifically for Cepheid and RR-Lyrae variable stars, while Fig. 22 shows the errors in proper motions for Cepheids and RR-Lyrae.
5.4. Physical parameters
Adding low-resolution spectral photometers on-board Gaia will make it capable of providing information on several object parameters including an estimate of line-of-sight extinction, effective temperature, metallicity, and surface gravity. Discussion of each individual physical parameter is given below.
Provided here are results for an approximation of the results of Liu et al. (2012), which reproduces CU8 results statistically but not individually for each star. Therefore for detailed analysis of specific object types, care should be taken. Again, due to the very long tails of the error distributions caused by large numbers of extremely faint stars, the mean values given below should be taken with caution.
Fig. 18 Histogram of error in G, GRVS, GRP, and GBP for all single stars. |
5.4.1. Extinction
Across many fields of astronomy, the effects of extinction on the apparent magnitude and colour of stars can play a major role in contributing to uncertainty. An accurate estimation of extinction will prove highly useful for many applications of the Gaia Catalogue.
Figure 23 shows the comparison between true extinction and the simulated Gaia estimate. For the vast majority of stars, the Gaia estimated extinction lies very close to the true value. This could prove very useful when, for example, using parallax and apparent magnitude data from Gaia because accurate extinction estimates are required to constrain the absolute magnitude of an object.
Fig. 19 HealPixMap in equatorial coordinates of the mean error in: top left: G; top right: GBP; lower left: GRP; lower right: GRVS. The colour scale gives the mean photometric error in mmag. The colour scales are different due to differences in the maximum mean magnitude. |
Additionally, these results show that the Gaia data will be highly useful in terms of mapping Galactic extinction in three dimensions, thanks to the combination of a large number of accurate parallax and extinction measurements. The negative extinction values in Fig. 23 are of course non-physical and are simply the result of applying a Gaussian random error to stars with near zero extinction.
Fig. 20 Cumulative histogram of the relative parallax error for all single stars, split by variability type. The histogram range displays 85% of all data. |
Fig. 21 Histogram of parallax error for Cepheid and RR-Lyrae variable stars. RR-Lyrae is a combination of the two sub-populations RR-ab and RR-c. |
Fig. 22 Histogram of proper motion error μα and μδ for Cepheid and RR-Lyrae variable stars. RR-Lyrae is a combination of the two sub-populations RR-ab and RR-c. |
The discontinuity at A0 = 1 in the top left-hand panel of Fig. 23 comes from the distinction made between high and low extinction stars in the presentation of the results in Liu et al. (2012). Our algorithm is based on results given in that paper, where the dependence on the extinction has been simplified to two cases, stars with A0< 1 and those with A0> 1. This distinction was made only for presentation of the results, and the real results from the DPAC algorithms will not show this discontinuity. Liu et al. (2012) report a degeneracy between extinction and effective temperature due to the lack of resolved spectral lines sensitive only to effective temperature.
Total number of single stars of each variability type, and the percentage of each that falls below each relative parallax error limit: 500%, 100%, 50%, 20%, 5%, and 1%.
5.4.2. Effective temperature
For all objects in the GOG catalogue, the measured effective temperature ranges between 850 and 102 000 K. The error in effective temperature is less than 640 K for all stars, with a mean value of 388 K. Figure 23 shows the comparison between true object effective temperature and the Gaia estimation. The thin lines visible in Fig. 23 are an artefact from the UM, which uses a Hess diagram to produce stars, leading to some quantisation in the effective temperature of simulated stars.
5.4.3. Metallicity
Metallicity can be estimated by Gaia in the form of [Fe/H]. Measured values range from −6.5 to +4.6. The mean error in metallicity estimate is 0.57 dex. The relatively high error in metallicity estimate can lead a large difference between real and observed values, as seen in Fig. 23.
5.4.4. Surface gravity
The mean error in surface gravity is 0.45 dex. The comparison between real and observed surface gravity can be seen in Fig. 23. As with metallicity, the lines at regular intervals at high gravity in this plot are due to the UM (Robin et al. 2012).
Fig. 23 Comparison of the true values of physical parameters with the GOG “observed” values for: top left, extinction; top right, metallicity; bottom left, surface gravity; and bottom right, effective temperature. The colour scales represent log density of objects in a bin size of: top left, 50 by 50 mmag; top right, 0.4 by 0.4 dex; bottom left, 0.5 by 0.5 dex; and bottom right, 100 by 100 K. White area represents zero stars. |
6. Conclusions
The Gaia Object Generator provides the most complete picture to date of what can be expected from the Gaia astrometric mission. Its simulated catalogue provides useful insight into how various types of objects will be observed and how each of their observables will appear after including observational errors and instrument effects. The simulated catalogue includes directly observed quantities, such as sky position and parallax, as well as derived quantities, such as interstellar extinction and metallicity.
Additionally, the full sky simulation described here is useful for gaining an idea of the size and format of the eventual Gaia Catalogue, for preparing tools and hardware for hosting and distribution of the data, and for becomeing familiar with working with such a large and rich dataset.
In addition to the stellar simulation described in this paper, there are plans to generate other simulated catalogues of interest, such as open clusters, Magallanic Clouds, supernovae, and other types of extragalactic objects, so that a more complete version of the simulated Gaia Catalogue can be compiled.
Here we have focused on the simulated catalogue from the inbuilt Gaia Universe Model, based on the Besançon Galaxy model. However, GOG can alternatively be supplied with an input catalogue generated by the user. This way, simulated data from any other model can be processed with GOG to obtain simulated Gaia observations of specific interest to the individual user. The input can be either synthetic data on a specific star or catalogue, or an entire simulated survey such as those generated using Galaxia (Sharma et al. 2011), provided a minimum of input information is supplied (e.g. position, distance, apparent magnitude, and colour).
With GOG, the capabilities of the instrument can be explored, and it is possible to gain insight into the expected performance for specific types of objects. While only a subset of the available statistics have been reproduced here, it is possible to obtain the full set of available statistics at request.
We are working to make the full simulated catalogue publicly available, so that interested individuals can begin working with data similar to the forthcoming Gaia Catalogue.
Acknowledgments
GOG is the product of many years of work from a number of people involved in DPAC and specifically CU2. The authors would like the thank the various CUs for contributing predicted error models for Gaia. The processing of the GOG data made significant use of the Barcelona Supercomputing Center (MareNostrum), and the authors would specifically like to thank Javi Castañeda, Marcial Clotet, and Aidan Fries for handling our computation and data handling needs. Additionally, thanks go to Sergi Blanco-Cuaresmo for his help with Matplotlib. This work was supported by the Marie Curie Initial Training Networks grant number PITN-GA-2010-264895 ITN “Gaia Research for European Astronomy Training”, and MINECO (Spanish Ministry of Economy) - FEDER through grants AYA2009-14648-C02-01, AYA2010-12176-E, AYA2012-39551-C02-01, and CONSOLIDER CSD2007-00050.
References
- Babusiaux, C. 2005, in The Three-Dimensional Universe with Gaia, eds. C. Turon, K. S. O’Flaherty, & M. A. C. Perryman, ESA SP, 576, 417 [Google Scholar]
- Bailer-Jones, C. A. L. 2011, MNRAS, 411, 435 [NASA ADS] [CrossRef] [Google Scholar]
- Bailer-Jones, C. A. L., Andrae, R., Arcay, B., et al. 2013, A&A, 559, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- de Bruijne, J. H. J. 2012, Ap&SS, 341, 31 [NASA ADS] [CrossRef] [Google Scholar]
- Isasi, Y., Figueras, F., Luri, X., & Robin, A. C. 2010, in Highlights of Spanish Astrophysics V, eds. J. M. Diego, L. J. Goicoechea, J. I. González-Serrano, & J. Gorgas (Heidelberg, Berlin: Springer-Verlag), 415 [Google Scholar]
- Jordi, C., Gebran, M., Carrasco, J. M., et al. 2010, A&A, 523, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Katz, D., Munari, U., Cropper, M., et al. 2004, MNRAS, 354, 1223 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, C., Bailer-Jones, C. A. L., Sordo, R., et al. 2012, MNRAS, 426, 2463 [NASA ADS] [CrossRef] [Google Scholar]
- Masana, E., Isasi, Y., Luri, X., & Peralta, J. 2010, in Highlights of Spanish Astrophysics V, eds. J. M. Diego, L. J. Goicoechea, J. I. González-Serrano, & J. Gorgas (Heidelberg, Berlin: Springer-Verlag), 515 [Google Scholar]
- Perryman, M. A. C., & ESA 1997, The Hipparcos and Tycho catalogues. Astrometric and photometric star catalogues derived from the ESA Hipparcos Space Astrometry Mission ESA SP (Noordwijk: ESA), 1200 [Google Scholar]
- Perryman, M. A. C., de Boer, K. S., Gilmore, G., et al. 2001, A&A, 369, 339 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Robin, A. C., Reylé, C., Derrière, S., & Picaud, S. 2003, A&A, 409, 523 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Robin, A. C., Luri, X., Reylé, C., et al. 2012, A&A, 543, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sharma, S., Bland-Hawthorn, J., Johnston, K. V., & Binney, J. 2011, ApJ, 730, 3 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Average end-of-mission formal error in radial velocity with an assumed average of 40 field-of-view transits, in km s-1, for each spectral type.
Mean and median value of the end-of-mission error in each observable, along with the upper (UQ) and lower (LQ) 25% quartile.
Total number of single stars for each spectral type, along with the percentage of those that fall below each relative parallax error limit; e.g., 68% of M-type stars have a relative parallax error better than 20%.
Total number of single stars for each luminosity class, along with the percentage that fall below each relative parallax error limit.
Total number of single stars for each stellar classification, along with the percentage that fall below each relative parallax error limit.
Total number of single stars of each variability type, and the percentage of each that falls below each relative parallax error limit: 500%, 100%, 50%, 20%, 5%, and 1%.
All Figures
Fig. 1 Skymap of total integrated flux over the Milky Way, in the G band. The colour bar represents a relative scale, from maximum flux in white to minimum flux in black. The figure is plotted in Galactic coordinates with the Galactic-longitude orientation swapped left to right. |
|
In the text |
Fig. 2 Mean end-of-mission error as a function of G magnitude for parallax, position, proper motion, and photometry in the four Gaia passbands. Additionally the mean end-of-mission error in radial velocity as a function of GRVS magnitude. |
|
In the text |
Fig. 3 Sky map (healpix) of mean parallax error for all single stars in equatorial coordinates. Colour scale is mean parallax error in μas. The red area is the location of the bulge. |
|
In the text |
Fig. 4 Histogram of parallax for all single stars. The histogram contains 99.5% of all data. |
|
In the text |
Fig. 5 Histogram of end-of-mission parallax error for all single stars, split by stellar population. |
|
In the text |
Fig. 6 Cumulative histogram of relative parallax error for all single stars, split by spectral type. The histogram range displays 74% of all data. |
|
In the text |
Fig. 7 Histogram of end-of-mission parallax error for all single stars, split by spectral type. |
|
In the text |
Fig. 8 End-of-mission parallax error against G magnitude for all single stars. The colour scale represents the log of density of objects in a bin size of 80 mmag by 2.5 μas. White area represents zero stars. |
|
In the text |
Fig. 9 End-of-mission parallax error against parallax for all single stars. The colour scale represents the log of density of objects in a bin size of 10 by 2.5 μas. White area represents zero stars. |
|
In the text |
Fig. 10 Right ascension error against real right ascension. The colour scale is linear, with a factor of 105. Histograms are computed for both right ascension and right ascension error. The colour scale represents log density of objects in a bin size of 2 degrees by 7.5 μas. White area represents zero stars. |
|
In the text |
Fig. 11 Declination error against real declination. The colour scale is linear, with a factor of 105. Histograms are computed for both declination and declination error. The colour scale represents log density of objects in a bin size of 1 degrees by 5 μas. White area represents zero stars. |
|
In the text |
Fig. 12 Error in proper motion for alpha and delta for all single stars. |
|
In the text |
Fig. 13 Histogram of radial velocity error split by G magnitude range. The histogram contains 100% of all data that have radial velocity information. |
|
In the text |
Fig. 14 Histogram of radial velocity error split by spectral type. The histogram contains 100% of all data that have radial velocity information. |
|
In the text |
Fig. 15 2D histograms showing error in proper motion against G magnitude. The colour scale represents the log density of objects in a bin size of 80 mmag by 2 μas yr-1. Left is proper motion in right ascension, and right is proper motion in declination. White area represents zero stars. |
|
In the text |
Fig. 16 End-of-mission error in radial velocity against GRVS magnitude. The colour scale represents the log density in a bin size of 50 mmag by 1 km s-1. White area represents zero stars. |
|
In the text |
Fig. 17 End-of-mission errors in photometry as a function of G magnitude. The colour scale represents the log density of objects in a bin size of 80 mmag by 0.4 mmag. Top left, G magnitude; top right, GBP; bottom left, GRP; bottom right, GRVS. White area represents zero stars. |
|
In the text |
Fig. 18 Histogram of error in G, GRVS, GRP, and GBP for all single stars. |
|
In the text |
Fig. 19 HealPixMap in equatorial coordinates of the mean error in: top left: G; top right: GBP; lower left: GRP; lower right: GRVS. The colour scale gives the mean photometric error in mmag. The colour scales are different due to differences in the maximum mean magnitude. |
|
In the text |
Fig. 20 Cumulative histogram of the relative parallax error for all single stars, split by variability type. The histogram range displays 85% of all data. |
|
In the text |
Fig. 21 Histogram of parallax error for Cepheid and RR-Lyrae variable stars. RR-Lyrae is a combination of the two sub-populations RR-ab and RR-c. |
|
In the text |
Fig. 22 Histogram of proper motion error μα and μδ for Cepheid and RR-Lyrae variable stars. RR-Lyrae is a combination of the two sub-populations RR-ab and RR-c. |
|
In the text |
Fig. 23 Comparison of the true values of physical parameters with the GOG “observed” values for: top left, extinction; top right, metallicity; bottom left, surface gravity; and bottom right, effective temperature. The colour scales represent log density of objects in a bin size of: top left, 50 by 50 mmag; top right, 0.4 by 0.4 dex; bottom left, 0.5 by 0.5 dex; and bottom right, 100 by 100 K. White area represents zero stars. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.