Search

Aims. We produce a clean and well-characterised catalogue of objects within 100 pc of the Sun from the Gaia Early Data Release 3. We characterise the catalogue through comparisons to the full data release, external catalogues, and simulations. We carry out a first analysis of the science that is possible with this sample to demonstrate its potential and best practices for its use., Methods. Theselection of objects within 100 pc from the full catalogue used selected training sets, machine-learning procedures, astrometric quantities, and solution quality indicators to determine a probability that the astrometric solution is reliable. The training set construction exploited the astrometric data, quality flags, and external photometry. For all candidates we calculated distance posterior probability densities using Bayesian procedures and mock catalogues to define priors. Any object with reliable astrometry and a non-zero probability of being within 100 pc is included in the catalogue., Results. We have produced a catalogue of 331 312 objects that we estimate contains at least 92% of stars of stellar type M9 within 100 pc of the Sun. We estimate that 9% of the stars in this catalogue probably lie outside 100 pc, but when the distance probability function is used, a correct treatment of this contamination is possible. We produced luminosity functions with a high signal-to-noise ratio for the main-sequence stars, giants, and white dwarfs. We examined in detail the Hyades cluster, the white dwarf population, and wide-binary systems and produced candidate lists for all three samples. We detected local manifestations of several streams, superclusters, and halo objects, in which we identified 12 members of Gaia Enceladus. We present the first direct parallaxes of five objects in multiple systems within 10 pc of the Sun., Conclusions. We provide the community with a large, well-characterised catalogue of objects in the solar neighbourhood. This is a primary benchmark for measuring and understanding fundamental parameters and descriptive functions in astronomy.

Estimating age distributions, or star formation histories, of stellar populations in the Milky Way is important in order to study the evolution of trends in elemental abundances and kinematics. We build on previous work to develop an algorithm for estimating the age-metallicity distribution which uses the full age-metallicity probability density functions (PDFs) of individual stars. No assumptions are made about the shape of the underlying distribution, and the only free parameter of the algorithm is used to ensure a smooth solution. In this work we use individual age-metallicity PDFs from isochrone fitting of stars with known metallicities. The method is tested with synthetic samples and is found to recover the input age-metallicity distribution more accurately than the distribution of individually estimated ages and metallicities. The recovered sample age distribution is always more accurate than the distribution of individual ages, also when restricted to the most precise individual ages. By applying the method to the stars in the Geneva-Copenhagen survey, we detect a possible minimum in the star formation history of the Solar neighbourhood at an age of 10 Gyr which is not seen in the distribution of individual ages. Although we apply the method only to age-metallicity distributions, the algorithm is described more generally and can in principle be applied in other parameter spaces. It is also not restricted to individual parameter distributions from isochrone fitting, meaning that a sample age distribution can be estimated based on individual age PDFs from other methods such as asteroseismology or gyrochronology., 20 pages, 16 figures. To be published in MNRAS. An implementation of the algorithm can be found at https://doi.org/10.5281/zenodo.3941733

Rowell, N., et al., [Context] The unprecedented astrometric precision of the Gaia mission relies on accurate estimates of the locations of sources in the Gaia data stream. This is ultimately performed by point spread function (PSF) fitting, which in turn requires an accurate reconstruction of the PSF, including calibrations of all the major dependences. These include a strong colour dependence due to Gaia's broad G band and a strong time dependence due to the evolving contamination levels and instrument focus. Accurate PSF reconstruction is also important for photometry. Gaia Early Data Release 3 (EDR3) will, for the first time, use a PSF calibration that models several of the strongest dependences, leading to signficantly reduced systematic errors. [Aims] We describe the PSF model and calibration pipeline implemented for Gaia EDR3, including an analysis of the calibration results over the 34 months of data. We include a discussion of the limitations of the current pipeline and directions for future releases. This will be of use both to users of Gaia data and as a reference for other precision astrometry missions. [Methods] We develop models of the 1D line spread function (LSF) and 2D PSF profiles based on a linear combination of basis components. These are designed for flexibility and performance, as well as to meet several mathematical criteria such as normalisation. We fit the models to selected primary sources in independent time ranges, using simple parameterisations for the colour and other dependences. Variation in time is smoothed by merging the independent calibrations in a square root information filter, with resets at certain mission events that induce a discontinuous change in the PSF. [Results] The PSF calibration shows strong time and colour dependences that accurately reproduce the varying state of the Gaia astrometric instrument. Analysis of the residuals reveals both the performance and the limitations of the current models and calibration pipeline, and indicates the directions for future development. [Conclusions] The PSF modelling and calibration carried out for Gaia EDR3 represents a major step forwards in the data processing and will lead to reduced systematic errors in the core mission data products. Further significant improvements are expected in the future data releases., The work described in this paper has been financially supported by the United Kingdom Particle Physics and Astronomy Research Council (PPARC), the United Kingdom Science and Technology Facilities Council (STFC), and the United Kingdom Space Agency (UKSA) through the following grants to the University of Bristol, the University of Cambridge, the University of Edinburgh, the University of Leicester, the Mullard Space Sciences Laboratory of University College London, and the United Kingdom Rutherford Appleton Laboratory (RAL): PP/D006511/1, PP/D006546/1, PP/D006570/1, ST/I000852/1, ST/J005045/1, ST/K00056X/1, ST/K000209/1, ST/K000756/1, ST/L006561/1, ST/N000595/1, ST/N000641/1, ST/N000978/1, ST/N001117/1, ST/S000089/1, ST/S000976/1, ST/S001123/1, ST/S001948/1, ST/S002103/1, and ST/V000969/1; the Spanish Ministry of Economy (MINECO/FEDER, UE) through grants ESP2016-80079-C2-1-R, RTI2018-095076-B-C21 and the Institute of Cosmos Sciences University of Barcelona (ICCUB, Unidad de Excelencia “María de Maeztu”) through grants MDM-2014-0369 and CEX2019-000918-M; the German Aerospace Agency (Deutsches Zentrum für Luft- und Raumfahrt e.V., DLR) through grants 50QG0501, 50QG0601, 50QG0602, 50QG0701, 50QG0901, 50QG1001, 50QG1101, 50QG1401, 50QG1402, 50QG1403, and 50QG1404 and the Centre for Information Services and High Performance Computing (ZIH) at the Technische Universität (TU) Dresden for generous allocations of computer time; the Agenzia Spaziale Italiana (ASI) through contracts I/037/08/0, I/058/10/0, 2014-025-R.0, 2014-025-R.1.2015 and 2018-24-HH.0 to the Italian Istituto Nazionale di Astrofisica; and the Swedish National Space Board (SNSB/Rymdstyrelsen). A.B. additionally acknowledges financial support from the Netherlands Research School for Astronomy (NOVA).

Aims. We aim to demonstrate the scientific potential of the Gaia Early Data Release 3 (EDR3) for the study of different aspects of the Milky Way structure and evolution and we provide, at the same time, a description of several practical aspects of the data and examples of their usage., Methods. We used astrometric positions, proper motions, parallaxes, and photometry from EDR3 to select different populations and components and to calculate the distances and velocities in the direction of the anticentre. In this direction, the Gaia astrometric data alone enable the calculation of the vertical and azimuthal velocities; also, the extinction is relatively low compared to other directions in the Galactic plane. We then explore the disturbances of the current disc, the spatial and kinematical distributions of early accreted versus in situ stars, the structures in the outer parts of the disc, and the orbits of open clusters Berkeley 29 and Saurer 1., Results. With the improved astrometry and photometry of EDR3, we find that: (i) the dynamics of the Galactic disc are very complex with oscillations in the median rotation and vertical velocities as a function of radius, vertical asymmetries, and new correlations, including a bimodality with disc stars with large angular momentum moving vertically upwards from below the plane, and disc stars with slightly lower angular momentum moving preferentially downwards; (ii) we resolve the kinematic substructure (diagonal ridges) in the outer parts of the disc for the first time; (iii) the red sequence that has been associated with the proto-Galactic disc that was present at the time of the merger with Gaia-Enceladus-Sausage is currently radially concentrated up to around 14 kpc, while the blue sequence that has been associated with debris of the satellite extends beyond that; (iv) there are density structures in the outer disc, both above and below the plane, most probably related to Monoceros, the Anticentre Stream, and TriAnd, for which the Gaia data allow an exhaustive selection of candidate member stars and dynamical study; and (v) the open clusters Berkeley 29 and Saurer 1, despite being located at large distances from the Galactic centre, are on nearly circular disc-like orbits., Conclusions. Even with our simple preliminary exploration of the Gaia EDR3, we demonstrate how, once again, these data from the European Space Agency are crucial for our understanding of the different pieces of our Galaxy and their connection to its global structure and history.

We present the early installment of the third Gaia data release, Gaia EDR3, consisting of astrometry and photometry for 1.8 billion sources brighter than magnitude 21, complemented with the list of radial velocities from Gaia DR2. Gaia EDR3 contains celestial positions and the apparent brightness in G for approximately 1.8 billion sources. For 1.5 billion of those sources, parallaxes, proper motions, and the (G_BP-G_RP) colour are also available. The passbands for G, G_BP, and G_RP are provided as part of the release. For ease of use, the 7 million radial velocities from Gaia DR2 are included in this release, after the removal of a small number of spurious values. New radial velocities will appear as part of Gaia DR3. Finally, Gaia EDR3 represents an updated materialisation of the celestial reference frame (CRF) in the optical, the Gaia-CRF3, which is based solely on extragalactic sources. The creation of the source list for Gaia EDR3 includes enhancements that make it more robust with respect to high proper motion stars, and the disturbing effects of spurious and partially resolved sources. The source list is largely the same as that for Gaia DR2, but it does feature new sources and there are some notable changes. The source list will not change for Gaia DR3. Gaia EDR3 represents a significant advance over Gaia DR2, with parallax precisions increased by 30 percent, proper motion precisions increased by a factor of 2, and the systematic errors in the astrometry suppressed by 30--40 percent for the parallaxes and by a factor ~2.5 for the proper motions. The photometry also features increased precision, but above all much better homogeneity across colour, magnitude, and celestial position. A single passband for G, G_BP, and G_RP is valid over the entire magnitude and colour range, with no systematics above the 1 percent level.

Gaia Early Data Release 3 (Gaia EDR3) gives trigonometric parallaxes for nearly 1.5 billion sources. Inspection of the EDR3 data for sources identified as quasars reveals that their parallaxes are biased, that is systematically offset from the expected distribution around zero, by a few tens of microarcsec. We attempt to map the main dependencies of the parallax bias in EDR3. In principle this could provide a recipe for correcting the EDR3 parallaxes. For faint sources the quasars provide the most direct way to estimate parallax bias. In order to extend this to brighter sources and a broader range of colours, we use differential methods based on physical pairs (binaries) and sources in the Large Magellanic Cloud. The functional forms of the dependencies are explored by mapping the systematic differences between EDR3 and DR2 parallaxes. The parallax bias is found to depend in a non-trivial way on (at least) the magnitude, colour, and ecliptic latitude of the source. Different dependencies apply to the five- and six-parameter solutions in EDR3. While it is not possible to derive a definitive recipe for the parallax correction, we give tentative expressions to be used at the researcher's discretion and point out some possible paths towards future improvements., Comment: 32 pages, 30 figures. Submitted to Astronomy & Astrophysics. Revised version where typos in (A.3) have been corrected. For associated code, see https://www.cosmos.esa.int/web/gaia/edr3-code (code was not affected by the typos)

Context. Gaia Early Data Release 3 (Gaia EDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions., Aims. The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar systembarycentre with respect to the rest frame of the Universe. Apart from being an important scientific result by itself, the acceleration measured in this way is a good quality indicator of the Gaia astrometric solution., Methods. Theeffect of the acceleration was obtained as a part of the general expansion of the vector field of proper motions in vector spherical harmonics (VSH). Various versions of the VSH fit and various subsets of the sources were tried and compared to get the most consistent result and a realistic estimate of its uncertainty. Additional tests with the Gaia astrometric solution were used to get a better idea of the possible systematic errors in the estimate., Results. Our best estimate of the acceleration based on Gaia EDR3 is (2.32 +/- 0.16) x 10(-10) m s(-2) (or 7.33 +/- 0.51 km s(-1) Myr-1) towards alpha = 269.1 degrees +/- 5.4 degrees, delta = -31.6 degrees +/- 4.1 degrees, corresponding to a proper motion amplitude of 5.05 +/- 0.35 mu as yr(-1). This is in good agreement with the acceleration expected from current models of the Galactic gravitational potential. We expect that future Gaia data releases will provide estimates of the acceleration with uncertainties substantially below 0.1 mu as yr(-1).

Under certain conditions, stellar radial velocities can be determined from astrometry, without any use of spectroscopy. This enables us to identify phenomena, other than the Doppler effect, that are displacing spectral lines. The change of stellar proper motions over time (perspective acceleration) is used to determine radial velocities from accurate astrometric data, which are now available from the Gaia and Hipparcos missions. Positions and proper motions at the epoch of Hipparcos are compared with values propagated back from the epoch of the Gaia Early Data Release 3. This propagation depends on the radial velocity, which obtains its value from an optimal fit assuming uniform space motion relative to the solar system barycentre. For 930 nearby stars we obtain astrometric radial velocities with formal uncertainties better than 100 km/s; for 55 stars the uncertainty is below 10 km/s, and for seven it is below 1 km/s. Most stars that are not components of double or multiple systems show good agreement with available spectroscopic radial velocities. Astrometry offers geometric methods to determine stellar radial velocity, irrespective of complexities in stellar spectra. This enables us to segregate wavelength displacements caused by the radial motion of the stellar centre-of-mass from those induced by other effects, such as gravitational redshifts in white dwarfs., Comment: 7 pages, 2 figures, 2 tables. Compared with the original submission, systematic corrections to the Gaia data have been improved as detailed in Section 3. Resulting changes in Tables 1 and 2 are small and mostly insignificant. Accepted for publication in Astronomy & Astrophysics

0, info:eu-repo/semantics/published

In the era of large-scale surveys of stars in the Milky Way, stellar ages are crucial for studying the evolution of the Galaxy. But determining ages of field stars is notoriously difficult; therefore, we attempt to determine benchmark ages for the extensively studied Gaia benchmark stars which can be used for validation purposes. By searching the literature for age estimates from different methods and deriving new ages based on Bayesian isochrone fitting, we are able to put reliable limits on the ages of 16 out of the 33 benchmark stars. The giants with well-defined ages are all young, and an expansion of the sample to include older giants with asteroseismic ages would be beneficial. Some of the stars have surface parameters inconsistent with isochrones younger than 16 Gyr. Including $\alpha$-enhancement in the models when relevant resolves some of these cases, but others clearly highlight discrepancies between the models and observations. We test the impact of atomic diffusion on the age estimates by fitting to the actual surface metallicity of the models instead of the initial value and find that the effect is negligible except for a single turn-off star. Finally, we show that our ability to determine isochrone-based ages for large spectroscopic surveys largely mirrors our ability to determine ages for these benchmark stars, except for stars with $\log g \gtrsim 4.4$ dex since their location in the HR diagram is almost age insensitive. Hence, isochrone fitting does not constrain their ages given the typical uncertainties of spectroscopic stellar parameters., Comment: Accepted in MNRAS. 69 pages (18 for main text, 11 for appendix, and 40 for extra figures)

We report the discovery of the first likely black hole in a non-interacting binary system with a field red giant. By combining radial velocity measurements from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) with photometric variability data from the All-Sky Automated Survey for Supernovae (ASAS-SN), we identified the bright rapidly-rotating giant 2MASS J05215658+4359220 as a binary system with a massive unseen companion. Subsequent radial velocity measurements reveal a system with an orbital period of 83 days and near-zero eccentricity. The photometric variability period of the giant is consistent with the orbital period, indicative of star spots and tidal synchronization. Constraints on the giant's mass and radius from its luminosity, surface gravity, and temperature imply an unseen companion with mass of $3.3^{+2.8}_{-0.7}$ M$_\odot$, indicating a low-mass black hole or an exceedingly massive neutron star. Measurement of the astrometric binary motion by {\it Gaia} will further characterize the system. This discovery demonstrates the potential of massive spectroscopic surveys like APOGEE and all-sky, high-cadence photometric surveys like ASAS-SN to revolutionize our understanding of the compact object mass function, and to test theories of binary star evolution and the supernova mechanism.

A deep understanding of the Milky Way galaxy, its formation and evolution requires observations of huge numbers of stars. Stellar photometry, therefore, provides an economical method to obtain intrinsic stellar parameters. With the addition of distance information - deriving reliable ages from photometry is a possibility. We have developed a Bayesian method that generates 2D probability maps of a star's age and metallicity from photometry and parallax using isochrones. Our synthetic tests show that including a near-UV passband enables us to break the degeneracy between a star's age and metallicity for certain evolutionary stages. It is possible to find well-constrained ages and metallicities for turn-off and sub-giant stars with colours including a U band and a parallax with uncertainty less than ~20%. Metallicities alone are possible for the main sequence and giant branch. We find good agreement with the literature when we apply our method to the Gaia benchmark stars, particularly for turn-off and young stars. Further tests on the old open cluster NGC 188, however, reveal significant limitations in the stellar isochrones. The ages derived for the cluster stars vary with evolutionary stage, such that turn-off ages disagree with those on the sub-giant branch, and metallicities vary significantly throughout. Furthermore, the parameters vary appreciably depending on which colour combinations are used in the derivation. We identify the causes of these mismatches and show that improvements are needed in the modelling of giant branch stars and in the creation and calibration of synthetic near-UV photometry. Our results warn against applying isochrone fitting indiscriminately. In particular, the uncertainty on the stellar models should be quantitatively taken into account. Further efforts to improve the models will result in significant advancements in our ability to study the Galaxy., Comment: Accepted for publication in A&A. 22 pages, 19 figures (+ 6 pages of appendix). Abstract abridged for the arXiv

(abridged) We develop Bayesian methods and detection criteria for orbital fitting, and revise the detectability of exoplanets in light of the in-flight properties of Gaia. Limiting ourselves to one-planet systems as a first step of the development, we simulate Gaia data for exoplanet systems over a grid of S/N, orbital period, and eccentricity. The simulations are then fit using Markov chain Monte Carlo methods. We investigate the detection rate according to three information criteria and the delta chi^2. For the delta chi^2, the effective number of degrees of freedom depends on the mission length. We find that the choice of the Markov chain starting point can affect the quality of the results; we therefore consider two limit possibilities: an ideal case, and a very simple method that finds the starting point assuming circular orbits. Using Jeffreys' scale of evidence, the fraction of false positives passing a strong evidence criterion is < ~0.2% (0.6%) when considering a 5 yr (10 yr) mission and using the Akaike information criterion or the Watanabe-Akaike information criterion, and, 18 pages, 11 figures. New version accepted by A&A for publication

Gaia Data Release 1 (Gaia DR1) contains astrometric results for more than 1 billion stars brighter than magnitude 20.7 based on observations collected by the Gaia satellite during the first 14 months of its operational phase. We give a brief overview of the astrometric content of the data release and of the model assumptions, data processing, and validation of the results. For stars in common with the Hipparcos and Tycho-2 catalogues, complete astrometric single-star solutions are obtained by incorporating positional information from the earlier catalogues. For other stars only their positions are obtained by neglecting their proper motions and parallaxes. The results are validated by an analysis of the residuals, through special validation runs, and by comparison with external data. Results. For about two million of the brighter stars (down to magnitude ~11.5) we obtain positions, parallaxes, and proper motions to Hipparcos-type precision or better. For these stars, systematic errors depending e.g. on position and colour are at a level of 0.3 milliarcsecond (mas). For the remaining stars we obtain positions at epoch J2015.0 accurate to ~10 mas. Positions and proper motions are given in a reference frame that is aligned with the International Celestial Reference Frame (ICRF) to better than 0.1 mas at epoch J2015.0, and non-rotating with respect to ICRF to within 0.03 mas/yr. The Hipparcos reference frame is found to rotate with respect to the Gaia DR1 frame at a rate of 0.24 mas/yr. Based on less than a quarter of the nominal mission length and on very provisional and incomplete calibrations, the quality and completeness of the astrometric data in Gaia DR1 are far from what is expected for the final mission products. The results nevertheless represent a huge improvement in the available fundamental stellar data and practical definition of the optical reference frame., Accepted for publication in Astronomy & Astrophysics

Context. As part of the data processing for Gaia Data Release 1 (Gaia DR1) a special astrometric solution was computed, the so-called auxiliary quasar solution. This gives positions for selected extragalactic objects, including radio sources in the second realisation of the International Celestial Reference Frame (ICRF2) that have optical counterparts bright enough to be observed with Gaia. A subset of these positions was used to align the positional reference frame of Gaia DR1 with the ICRF2. Although the auxiliary quasar solution was important for internal validation and calibration purposes, the resulting positions are in general not published in Gaia DR1. Aims. We describe the properties of the Gaia auxiliary quasar solution for a subset of sources matched to ICRF2, and compare their optical and radio positions at the sub-mas level. Methods. Descriptive statistics are used to characterise the optical data for the ICRF sources and the optical-radio differences. The most discrepant cases are examined using online resources to find possible alternative explanations than a physical optical-radio offset of the quasars. Results. In the auxiliary quasar solution 2191 sources have good optical positions matched to ICRF2 sources with high probability. Their formal standard errors are better than 0.76 milliarcsec (mas) for 50% of the sources and better than 3.35 mas for 90%. Optical magnitudes are obtained in Gaia’s unfiltered photometric G band. The Gaia results for these sources are given as a separate table in Gaia DR1. The comparison with the radio positions of the defining sources shows no systematic differences larger than a few tenths of a mas. The fraction of questionable solutions, not readily accounted for by the statistics, is less than 6%. Normalised differences have extended tails requiring case-by-case investigations for around 100 sources, but we have not seen any difference indisputably linked to an optical-radio offset in the sources. Conclusions. With less than a quarter of the data expected from the nominal mission it has been possible to obtain positions at the sub-mas level for most of the ICRF sources having an optical counterpart brighter than 20.5 mag.

Context. Gaia Data Release 2 (DR2) offers unparalleled precision on stars’ parallaxes and proper motions. This allows the prediction of microlensing events for which the lens stars (and any planets they possess) are nearby and may be well studied and characterised. Aims. We identify a number of potential microlensing events that will occur before the year 2035.5, 20 years from the Gaia DR2 reference epoch. Methods. We query Gaia DR2 for potential lenses within 100 pc, extract parallaxes and proper motions of the lenses and background sources, and identify potential lensing events. We estimate the lens masses from Priam effective temperatures and use these to calculate peak magnifications and the size of the Einstein radii relative to the lens stars’ habitable zones. Results. We identify seven future events with a probability >10% of an alignment within one Einstein radius. Of particular interest is DR2 5918299904067162240 (WISE J175839.20–583931.6), magnitude G = 14.9, which will lens a G = 13.9 background star in early 2030, with a median 23% net magnification. Other pairs are typically fainter, hampering characterisation of the lens (if the lens is faint) or the ability to accurately measure the magnification (if the source is much fainter than the lens). Of timely interest is DR2 4116504399886241792 (2MASS J17392440–2327071), which will lens a background star in July 2020, albeit with weak net magnification (0.03%). Median magnifications for the other five high-probability events range from 0.3% to 5.3%. The Einstein radii for these lenses are one to ten times the radius of the habitable zone, allowing these lensing events to pick out cold planets around the ice line, and filling a gap between transit and current microlensing detections of planets around very low-mass stars. Conclusions. We provide a catalogue of the predicted events to aid future characterisation efforts. Current limitations include a lack of many high-proper-motion objects in Gaia DR2 and often large uncertainties on the proper motions of the background sources (or only two-parameter solutions). Both of these deficiencies will be rectified with Gaia DR3 in 2020. Further characterisation of the lenses is also warranted to better constrain their masses and predict the photometric magnifications.

Accurate characterization of the errors in the global astrometric solution for Gaia is essential for making optimal use of the catalogue data. We investigate the structure of the covariance between the estimated astrometric parameters by studying the properties of the astrometric least squares solution. We find that astrometric errors can be separated in a star and an attitude part, due to the estimation of the star and attitude parameters respectively. Hence the covariances can be separated in a star, an attitude and a cross term. This is demonstrated using our scalable simulation tool AGISLab, where the covariances are estimated statistically using Monte Carlo techniques.

The ESA space astrometry mission Gaia, due for launch in early 2012, will in addition to its huge output of fundamental astrometric and astrophysical data also provide stringent tests of general relativity. In this paper we present an updated analysis of Gaia's capacity to measure the PPN parameter γ as part of its core astrometric solution. The analysis is based on small-scale astrometric solutions taking into account the simultaneous determination of stellar astrometric parameters and the satellite attitude. In particular, the statistical correlation between PPN γ and the stellar parallaxes is considered. Extrapolating the results to a full-scale solution using some 100 million stars, we find that PPN γ could be obtained to about 10−6, which is significantly better than today's best estimate from the Cassini mission of 2 × 10−5.

Accurate characterization of the astrometric errors in the forthcoming Gaia catalogue is essential for making optimal use of the data. Using small-scale numerical simulations of the astrometric solution, we investigate the expected spatial correlation between the astrometric errors of stars as function of their angular separation. Extrapolating to the full-scale solution for the final Gaia catalogue, we find that the expected correlations are generally very small, but could reach some fraction of a percent for angular separations smaller than about one degree. The spatial correlation length is related to the size of the field of view of Gaia, while the maximum correlation coefficient is related to the mean number of stars present in the field at any time. Our scalable simulation tool (AGISLab) makes it possible to characterize the astrometric errors and correlations, e.g., as functions of position and magnitude.

We investigate whether it is possible that viable microbes could have been transported to Earth from the planets in extra-solar systems by means of natural vehicles such as ejecta expelled by comet or asteroid impacts on such planets. The probabilities of close encounters with other solar systems are taken into account as well as the limitations of bacterial survival times inside ejecta in space, caused by radiation and DNA decay. The conclusion is that no potentially DNA/RNA life-carrying ejecta from another solar system in the general Galactic star field landed on Earth before life already existed on Earth, not even if microbial survival time in space is as long as tens of millions of years. However, if the Sun formed initially as a part of a star cluster, as is commonly assumed, we cannot rule out the possibility of transfer of life from one of the sister systems to us. Likewise, there is a possibility that some extra-solar planets carry life that originated in our solar system. It will be of great interest to identify the members of the Sun's birth cluster of stars and study them for evidence for planets and life on the planets. The former step may be accomplished by the GAIA mission, the latter step by the SIM and DARWIN missions. Therefore it may not be too long until we have experimental knowledge on the question whether the natural transfer of life from one solar system to another has actually taken place., Comment: 25 pages, 1 table, accepted to ApJ

The ESA space astrometry mission Gaia will measure the positions, parallaxes and proper motions of the 1 billion brightest stars on the sky. Expected accuracies are in the 7–25 μas range down to 15 mag and sub-mas accuracies at the faint limit (20 mag). The astrometric data are complemented by low-resolution spectrophotometric data in the 330–1000 nm wavelength range and, for the brighter stars, radial velocity measurements. The scientific case covers an extremely wide range of topics in galactic and stellar astrophysics, solar system and exoplanet science, as well as the establishment of a very accurate, dense and faint optical reference frame. With a planned launch around 2012 and an (extended) operational lifetime of 6 years, final results are expected around 2021. We give a brief overview of the science goals of Gaia, the overall project organisation, expected performance, and some key technical features and challenges.

Gaia's very accurate astrometric measurements will allow the International Celestial Reference Frame (ICRF) to be improved by a few orders of magnitude in the optical. Several sets of quasars are used to define a kinematically stable non-rotating reference frame with the barycentre of the Solar System as its origin. Gaia will also observe a large number of galaxies which could obtain accurate positions and proper motions although they are not point-like. The optical stability of the quasars is critical and we investigate how accurately the reference frame can be recovered. Various proper motion patterns are also present in the data, the best known is caused by the acceleration of the Solar System Barycentre, presumably, towards the Galactic centre. We review some other less-well-known effects that are not part of standard astrometric models. We model quasars and galaxies using realistic sky distributions, magnitudes and redshifts. Position variability is introduced using a Markov chain model. The reference frame is determined using the algorithm developed for the Gaia mission which also determines the acceleration of the Solar System. We also test a method to measure the velocity of the Solar System barycentre in a cosmological frame. We simulate the recovery of the reference frame and the acceleration of the Solar System and conclude that they are not significantly disturbed in the presence of quasar variability which is statistically averaged. However, the effect of a non-uniform sky distribution of the quasars can result in a correlation between the reference frame and acceleration which degrades the solution. Our results suggest that an attempt should be made to astrometrically determine the redshift dependent apparent drift of galaxies due to our velocity relative to the CMB, which in principle could allow the determination of the Hubble parameter., Comment: published in A&A, revised version (v2) (Abstract is same as v1 as the character limit is 1920, see the pdf for v2)

Context. The Gaia project will determine positions, proper motions, and parallaxes for more than one billion stars in our Galaxy. It is known that Gaia's two telescopes are affected by a small but significant variation of the basic angle between them. Unless this variation is taken into account during data processing, e.g. using on-board metrology, it causes systematic errors in the astrometric parameters, in particular a shift of the parallax zero-point. Previously, we suggested an early reduction of Gaia data for the subset of Tycho-2 stars (Tycho-Gaia Astrometric Solution; TGAS). Aims. We aim to investigate whether quasars can be used to independently verify the parallax zero-point already in early data reductions. This is not trivially possible as the observation interval is too short to disentangle parallax and proper motion for the quasar subset. Methods. We repeat TGAS simulations but additionally include simulated Gaia observations of quasars from ground-based surveys. All observations are simulated with basic angle variations. To obtain a full astrometric solution for the quasars in TGAS we explore the use of prior information for their proper motions. Results. It is possible to determine the parallax zero-point for the quasars with a few {\mu}as uncertainty, and it agrees to a similar precision with the zero-point for the Tycho-2 stars. The proposed strategy is robust even for quasars exhibiting significant fictitious proper motion due to a variable source structure, or when the quasar subset is contaminated with stars misidentified as quasars. Conclusions. Using prior information about quasar proper motions we could provide an independent verification of the parallax zero-point in early solutions based on less than one year of Gaia data., Comment: Astronomy & Astrophysics, accepted 25 October 2015, in press. Version 2 contains a few language improvements and a terminology change from 'fictitious proper motions' to 'spurious proper motions'

We present a new method, using Bayesian estimation, to determine stellar ages and their uncertainties from observational data and theoretical isochrones. The result for an individual star is obtained as the relative posterior probability density as function of the age ("G function"). From this can be derived the most probable age and confidence intervals. The convoluted morphology of isochrones and strong non-linearities make the age determination by any method difficult and susceptible to statistical biases, and as a result age uncertainties havee often been underestimated in the literature. From simulations we find that the G functions provide a general, robust and reliable way to quantify age information. Resulting age estimates are at least as accurate as those obtained with conventional isochrone fitting methods, and in some cases much better, especially when the observational uncertainties are large. We also find that undetected binaries, on the whole, have a surprisingly small effect on the age determinations. For a stellar sample, the individual G functions can be combined to derive the star formation history of the population; this will be developed in a forthcoming paper. For a coeval population the combination simplifies to computing the product of the individual G functions, and we apply that method to estimate the ages of the two open clusters IC 4651 and M67, using Padova isochrones and photometric data from the literature. For IC 4651 we find an estimated age of 1.56 +/- 0.03 Gyr, assuming a true distance modulus of 9.80. For M67 we find 4.05 +/- 0.05 Gyr for true distance modulus 9.48. The small formal errors of these age estimates do not include the ( much larger) uncertainties from calibration and model errors, but illustrate the statistical power of combining G functions. Our statistical approach to the age determination problem is well suited for the mass treatment of data resulting from large-scale surveys such as the Gaia mission. (Less)

Wavelength measurements in stellar spectra cannot readily be interpretedas true stellar motion on the sub-km s-1 accuracy level dueto the presence of many other effects, such as gravitational redshiftand stellar convection, which also produce line shifts. Following arecommendation by the IAU, the result of an accurate spectroscopicradial-velocity observation should therefore be given as the``barycentric radial-velocity measure'', i.e. the absolute spectralshift as measured by an observer at zero gravitational potential locatedat the solar-system barycentre. Standard procedures for reducingaccurate radial-velocity observations should be reviewed to take intoaccount this recommendation. We describe a procedure to determineaccurate barycentric radial-velocity measures of bright stars, based ondigital cross-correlation of spectra obtained with the ELODIEspectrometer (Observatoire de Haute-Provence) with a synthetic templateof Fe I lines. The absolute zero point of the radial-velocity measuresis linked to the wavelength scale of the Kurucz (1984) Solar Flux Atlasvia ELODIE observations of the Moon. Results are given for the Sun and42 stars, most of them members of the Hyades and Ursa Major clusters.The median internal standard error is 27 m s-1. The externalerror is estimated at around 120 m s-1, mainly reflecting theuncertainty in the wavelength scale of the Solar Flux Atlas. For the Sunwe find a radial-velocity measure of +257+/- 11 m s-1referring to the full-disk spectrum of the selected Fe I lines. Based onobservations made at Observatoire de Haute-Provence

We discuss the point spread function for the Spectro instrument based on realistic effects of the time delayed integration, transverse image motion and pixelization, but assuming rather small wave front errors. Stars brighter than 10 mag, which will saturate the CCD, may be observed photometrically using the image wings, but the central wavelengths will be slightly shifted towards the red.

While the precise relationship between the Milky Way disk and the symmetry planes of the dark matter halo remains somewhat uncertain, a time-varying disk orientation with respect to an inertial reference frame seems probable. Hierarchical structure formation models predict that the dark matter halo is triaxial and tumbles with a characteristic rate of ~2 rad/Hubble time (~30 muas/yr). These models also predict a time-dependent accretion of gas, such that the angular momentum vector of the disk should be misaligned with that of the halo. These effects, as well as tidal effects of the LMC, will result in the rotation of the angular momentum vector of the disk population with respect to the quasar reference frame. We assess the accuracy with which the positions and proper motions from Gaia can be referred to a kinematically non-rotating system, and show that the spin vector of the transformation from any rigid self-consistent catalog frame to the quasi-inertial system defined by quasars should be defined to better than 1 muas/yr. Determination of this inertial frame by Gaia will reveal any signature of the disk orientation varying with time, improve models of the potential and dynamics of the Milky Way, test theories of gravity, and provide new insights into the orbital evolution of the Sagittarius dwarf galaxy and the Magellanic Clouds., Comment: 16 pages; accepted for publication in ApJ

Hipparcos observations of double stars, mainly A and F primaries, are modelled in order to check the statistics of double-star solutions in the Hipparcos Catalogue and to set quantitative limits on binary distributions. The number of companions per primary is estimated to 0.25 ± 0.05 for semi-major axes from 1 to 10 AU, and to 0.43 ± 0.20 for the range 10–100 AU. The total multiplicity found is between 0.9 and 1.2.

Context. Kinematical data such as the mean velocities and velocity dispersions of stellar samples are useful tools to study galactic structure and evolution. However, observational data are often incomplete (e.g., lacking the radial component of the motion) and may have significant observational errors. For example, the majority of faint stars observed with Gaia will not have their radial velocities measured. Aims. Our aim is to formulate and test a new maximum likelihood approach to estimating the kinematical parameters for a local stellar sample when only the transverse velocities are known (from parallaxes and proper motions). Methods. Numerical simulations using synthetically generated data as well as real data (based on the Geneva-Copenhagen survey) are used to investigate the statistical properties (bias, precision) of the method, and to compare its performance with the much simpler "projection method" described by Dehnen & Binney (1998). Results. The maximum likelihood method gives more correct estimates of the dispersion when observational errors are important, and guarantees a positive-definite dispersion matrix, which is not always obtained with the projection method. Possible extensions and improvements of the method are discussed., 7 pages, 2 figures. Accepted for publication in Astronomy & Astrophysics

Context. Accurate characterization of the astrometric errors in the forthcoming Gaia Catalogue will be essential for making optimal use of the data. This includes the correlations among the estimated astrometric parameters of the stars as well as their standard uncertainties, i.e., the complete (variance-)covariance matrix of the relevant astrometric parameters. Aims. Because a direct computation of the covariance matrix is infeasible due to the large number of parameters, approximate methods must be used. The aim of this paper is to provide a mathematical basis for estimating the variance-covariance of any pair of astrometric parameters, and more generally the covariance matrix for multidimensional functions of the astrometric parameters. The validation of this model by means of numerical simulations will be considered in a forthcoming paper. Methods. Based on simplifying assumptions (in particular that calibration errors can be neglected), we derive and analyse a series expansion of the covariance matrix of the least-squares solution. A recursive relation for successive terms is derived and interpreted in terms of the propagation of errors from the stars to the attitude and back. We argue that the expansion should converge rapidly to useful precision. The recursion is vastly simplified by using a kinematographic (step-wise) approximation of the attitude model. Results. Low-order approximations of arbitrary elements from the covariance matrix can be computed efficiently in terms of a limited amount of pre-computed data representing compressed observations and the structural relationships among them. It is proposed that the user interface to the Gaia Catalogue should provide the tools necessary for such computations. (Less)

The HIPPARCOS satellite observes 118,000 stars from a predefined list. Some 10% of these are known doubles or multiples, and another 2—3% may ultimately be found to be non-single with separations above 0″.l and magnitude–differences below 4–5. This paper describes briefly the special reduction methods used by NDAC (Northern Data Analysis Consortium) for treating these non-singles. A key feature is the use of the main reduction results to calibrate and collect the data for individual doubles in ‘case history files’. This enables a global solution of both absolute and relative parameters for a double to be carried out in a single step. Using provisional data tapes covering small parts of a 14-month interval, tests have been made of major parts of the double star reductions. These early results are compared with independent HIPPARCOS data and with ground–based speckle interferometry.

We investigate the astrometric effects of stellar surface structures as a practical limitation to ultra-high-precision astrometry, e.g. in the context of exoplanet searches, and to quantify the expected effects in different regions of the HR-diagram. Stellar surface structures are likely to produce fluctuations in the integrated flux and radial velocity of the star, as well as a variation of the observed photocentre, i.e. astrometric jitter, and closure phase. We use theoretical considerations supported by Monte Carlo simulations to derive statistical relations between the corresponding astrometric, photometric, and radial-velocity effects. For most stellar types the astrometric jitter due to stellar surface structures is expected to be of order 10 micro-AU or greater. This is more than the astrometric displacement typically caused by an Earth-size exoplanet in the habitable zone, which is about 1-4 micro-AU for long-lived main-sequence stars. Only for stars with extremely low photometric variability (, 12 pages, 4 figures, submitted to A&A

The European Gaia astrometry mission is due for launch in 2011. Gaia will rely on the proven principles of ESA's Hipparcos mission to create an all-sky survey of about one billion stars throughout our Galaxy and beyond, by observing all objects down to 20th magnitude. Through its massive measurement of stellar distances, motions and multi-colour photometry it will provide fundamental data necessary for unravelling the structure, formation and evolution of the Galaxy. This paper presents the design and performance of the broad- and medium-band set of photometric filters adopted as the baseline for Gaia. The nineteen selected passbands (extending from the ultraviolet to the far-red), the criteria, and the methodology on which this choice has been based are discussed in detail. We analyse the photometric capabilities for characterizing the luminosity, temperature, gravity and chemical composition of stars. We also discuss the automatic determination of these physical parameters for the large number of observations involved, for objects located throughout the entire Hertzsprung-Russell diagram. Finally, the capability of the photometric system to deal with the main Gaia science case is outlined., MNRAS, accepted 1-Dec-2005; 26 pages, 15 figures

Spectral-line displacements away from the wavelengths naively expected from the Doppler shift caused by stellar radial motion may originate as convective shifts (correlated velocity and brightness patterns in the photosphere), as gravitational redshifts, or perhaps be induced by wave motions. Absolute lineshifts, in the past studied only for the Sun, are now accessible also for other stars thanks to astrometric determination of stellar radial motion, and spectrometers with accurate wavelength calibration. Comparisons between spectroscopic apparent radial velocities and astrometrically determined radial motions reveal greater spectral blueshifts in F-type stars than in the Sun (as theoretically expected from their more vigorous convection), further increasing in A-type stars (possibly due to atmospheric shockwaves). An important near-future development to enable a further analysis of stellar surface structure will be the study of wavelength variations across spatially resolved stellar disks, e.g., the center-to-limb wavelength changes along a stellar diameter, and their spatially resolved time variability., Comment: Proc. 13th Cool Stars Workshop, Hamburg, 5-9 July 2004, F.Favata et al., eds. 7 pages, 8 figures

ESA’s Gaia mission will provide accurate astrometric, photometric and spectroscopic data for millions of stars in the Galaxy and beyond. We discuss the prospects of using this vast amount of astrophysical information for understanding star formation histories in the Galaxy and its immediate neighbors within the Local Group.

We discuss the tracing of star formation histories with ESA's space astrometry mission GAIA, emphasizing the advantages of AGB stars for this purpose. GAIA's microarcsecond-level astrometry, multi-band photometry and spectroscopy will provide individual distances, motions, effective temperatures, gravities and metallicities for vast numbers of AGB stars in the Galaxy and beyond. Reliable ages of AGB stars can be determined to distances of \~200 kpc in a wide range of ages and metallicities, allowing star formation histories to be studied in a diversity of astrophysical environments.

Spectroscopic radial velocities for several nearby open clusters suggest that spectra of (especially earlier-type) rapidly rotating stars are systematically blueshifted by 3 km/s or more, relative to the spectra of slowly rotating ones. Comparisons with astrometrically determined radial motions in the Hyades suggests this to be an absolute blueshift, relative to wavelengths naively expected from stellar radial motion and gravitational redshift. Analogous trends are seen also in most other clusters studied (Pleiades, Coma Berenices, Praesepe, Alpha Persei, IC 2391, NGC 6475, IC 4665, NGC 1976 and NGC 2516). Possible mechanisms are discussed, including photospheric convection, stellar pulsation, meridional circulation, and shock-wave propagation, as well as effects caused by template mismatch in determining wavelength displacements. For early-type stars, a plausible mechanism is shock-wave propagation upward through the photospheric line-forming regions. Such wavelength shifts thus permit studies of certain types of stellar atmospheric dynamics and - irrespective of their cause - may influence deduced open-cluster membership (when selected from common velocity) and deduced cluster dynamics (some types of stars might show fortuitous velocity patterns)., Comment: Accepted by Astronomy & Astrophysics; 6 pages, 3 figures

This paper discusses the IAU Resolutions B1.3, B1.4, B1.5 and B1.9 (2000) that were adopted during the 24th General Assembly in Manchester, 2000 and provides details and explanations for these Resolutions. It is explained why they present significant progress over the corresponding IAU 1991 Resolutions and why they are necessary in the light of present accuracies in astrometry, celestial mechanics and metrology. In fact most of these Resolutions are consistent with astronomical models and software already in use. The metric tensors and gravitational potentials of both the Barycentric Celestial Reference System and Geocentric Celestial Reference System are defined and discussed. The necessity and relevance of the two celestial reference systems are explained. The transformations of coordinates and gravitational potentials are discussed. Potential coefficients parameterizing the post-Newtonian gravitational potentials are expounded. Simplified versions of the time transformations suitable for modern clock accuracies are elucidated. Various approximations used in the Resolutions are explicated and justified. Some models (e.g. for higher spin moments) that serve the purpose for estimating orders of magnitude have actually never been published before., 45 pages, 1 figure, submitted to Astronomical Journal

Radial motions of stars in nearby moving clusters are determined from accurate proper motions and trigonometric parallaxes, without any use of spectroscopy. Assuming that cluster members share the same velocity vector (apart from a random dispersion), we apply a maximum-likelihood method on astrometric data from Hipparcos to compute radial and space velocities (and their dispersions) in the Ursa Major, Hyades, Coma Berenices, Pleiades, and Praesepe clusters, and for the alpha Persei, Scorpius-Centaurus, and `HIP 98321' associations. The radial motion of the Hyades cluster is determined to within 0.47 km/s (standard error), and that of its individual stars to within 0.6 km/s. For other clusters, Hipparcos data yield astrometric radial velocities with typical accuracies of a few km/s. A comparison of these astrometric values with spectroscopic radial velocities in the literature shows a good general agreement and, in the case of the best-determined Hyades cluster, also permits searches for subtle astrophysical differences, such as evidence for enhanced convective blueshifts of F-dwarf spectra, and decreased gravitational redshifts in giants. Similar comparisons for the Scorpius OB2 complex indicate some expansion of its associations, albeit slower than expected from their ages. As a by-product from the radial-velocity solutions, kinematically improved parallaxes for individual stars are obtained, enabling Hertzsprung-Russell diagrams with unprecedented accuracy in luminosity. For the Hyades (parallax accuracy 0.3 mas), its main sequence resembles a thin line, possibly with wiggles in it., 18 pages, 7 figures, accepted for publication in A&A

The GAIA astrometric mission has recently been approved as one of the next two `cornerstones' of ESA's science programme, with a launch date target of not later than mid-2012. GAIA will provide positional and radial velocity measurements with the accuracies needed to produce a stereoscopic and kinematic census of about one billion stars throughout our Galaxy (and into the Local Group), amounting to about 1 per cent of the Galactic stellar population. GAIA's main scientific goal is to clarify the origin and history of our Galaxy, from a quantitative census of the stellar populations. It will advance questions such as when the stars in our Galaxy formed, when and how it was assembled, and its distribution of dark matter. The survey aims for completeness to V=20 mag, with accuracies of about 10 microarcsec at 15 mag. Combined with astrophysical information for each star, provided by on-board multi-colour photometry and (limited) spectroscopy, these data will have the precision necessary to quantify the early formation, and subsequent dynamical, chemical and star formation evolution of our Galaxy. Additional products include detection and orbital classification of tens of thousands of extra-Solar planetary systems, and a comprehensive survey of some 10^5-10^6 minor bodies in our Solar System, through galaxies in the nearby Universe, to some 500,000 distant quasars. It will provide a number of stringent new tests of general relativity and cosmology. The complete satellite system was evaluated as part of a detailed technology study, including a detailed payload design, corresponding accuracy assesments, and results from a prototype data reduction development., Accepted by A&A: 25 pages, 8 figures

High-accuracy astrometry permits the determination of not only stellar tangential motion, but also the component along the line-of-sight. Such non-spectroscopic (i.e. astrometric) radial velocities are independent of stellar atmospheric dynamics, spectral complexity and variability, as well as of gravitational redshift. Three methods are analysed: (1) changing annual parallax, (2) changing proper motion and (3) changing angular extent of a moving group of stars. All three have significant potential in planned astrometric projects. Current accuracies are still inadequate for the first method, while the second is marginally feasible and is here applied to 16 stars. The third method reaches high accuracy (, 13 pages, 2 figures. Accepted for publication in Astronomy & Astrophysics (main journal)

The Hipparcos Transit Data are a collection of partially reduced, fully calibrated observations of (mostly) double and multiple stars obtained with the ESA Hipparcos astrometry satellite. The data are publicly available, as part of the CD-ROM set distributed with the Hipparcos and Tycho Catalogues (ESA SP--1200, 1997), for about a third of the Hipparcos Catalogue entries including all confirmed or suspected non-single stars. The Transit Data consist of signal modulation parameters derived from the individual transits of the targets across the Hipparcos focal grid. The Transit Data permit re-reduction of the satellite data for individual objects, using arbitrarily complex object models in which time-variable photometric as well as geometric characteristics may be taken into account. We describe the structure and contents of the Transit Data files and give examples of how the data can be used. Some of the applications use standard astronomical software: Difmap or AIPS for aperture synthesis imaging, and GaussFit for detailed model fitting. Fortran code converting the data into formats suitable for these application programs has been made public in order to encourage and facilitate the use of Hipparcos Transit Data., Comment: A&AS, accepted for publication, 17 pages, 9 figures, 1 Table, Software available via http://www.astro.lu.se/~lennart/TD/index.html, Figures 4, 5, 6 and 7 need to copied separately, A complete postscript file can be found at http://www.astro.lu.se/~lennart/TD/ds1699.ps.gz

To assess the current status in the theoretical modeling of the spectral properties of late-type giants, we provide a comparison of synthetic photometric colors of late-type giants (calculated with PHOENIX, MARCS and ATLAS model atmospheres) with observations, at [M/H]=0.0 and -2.0. Overall, there is a good agreement between observed and synthetic colors, and synthetic colors and published Teff-color relations, both at [M/H]=0.0 and -2.0. Deviations from the observed trends in Teff-color planes are generally within \pm 150K (or less) in the effective temperature range Teff=3500-4800K. Synthetic colors calculated with different stellar atmosphere models typically agree to ~100K, within a large range of effective temperatures and gravities. Some discrepancies are seen in the Teff-(B-V) plane below Teff~3800K at [M/H]=0.0, due to difficulties in reproducing the 'turn-off' to the bluer colors which is seen in the observed data at Teff~3600K. Note that at [M/H]=-2.0 effective temperatures given by the scale of Alonso et al. (1999) are generally lower than those resulting from other Teff-color relations based both on observed and synthetic colors., Comment: 2 pages, 1 figure. Proceedings of the IAU Symposium 232 "The Scientific Requirements for Extremely Large Telescopes", eds. P. Whitelock, B. Leibundgut, and M. Dennefeld

A number of large spectroscopic surveys of stars in the Milky Way are under way or are being planned. In this context it is important to discuss the extent to which elemental abundances can be used as discriminators between different (known and unknown) stellar populations in the Milky Way. We aim to establish the requirements in terms of precision in elemental abundances, as derived from spectroscopic surveys of the Milky Way's stellar populations, in order to detect interesting substructures in elemental abundance space. We present a simple relation between the minimum number of stars needed to detect a given substructure and the precision of the measurements. The results are in agreement with recent small- and large-scale studies, with high and low precision, respectively. Large-number statistics cannot fully compensate for low precision in the abundance measurements and each survey should carefully evaluate what the main science drivers are for the survey and ensure that the chosen observational strategy will result in the precision necessary to answer the questions posed., Comment: 6 pages, 6 figures. Accepted for publication in Astronomy & Astrophysics

GAIA represents a preliminary concept for an astrometric mission being considered in the context of ESA's `Horizon 2000 Plus' long-term scientific program. It comprises three stacked Fizeau interferometers viewing different directions within an instantaneous scanning circle, each interferometer consisting of two 50 cm aperture mirrors with a baseline separation of 2.5 m. Equipped with a modulating grid, and using CCDs, at least as the baseline detector, repeated scanning of the celestial sphere over a period of five years is estimated to lead to positions, proper motions, and parallaxes of some 50 million objects, down to about V equals 15 mag, with an accuracy of better than 10 microarcsec, along with multi-color multi-epoch photometry of each object. The scientific case for such a mission is compelling: distances and kinematical motions for objects throughout our Galaxy would be obtained, along with valuable information on the space-time metric ((gamma) ), angular diameters of hundreds of stars, and a vast body of information on double and multiple systems. Screening of all 100,000 stars within 100 pc for periodic photocentric motions would provide the most powerful and systematic method of detecting possible planetary companions proposed to date.

Instrumentation and observing methods are being developed for a program in optical high-speed astrophysics, an exploratory project entering the domains of milli-, micro-, and nanosecond variability. Current studies include accretion flows around compact objects, stellar scintillation, and astronomical quantum optics. To study such rapid phenomena is not possible everywhere in the spectrum (e.g., X-ray studies are constrained by the photon count rates feasible with current spacecraft). The ground-based optical is a promising region, for which we have constructed a dedicated instrument, QVANTOS (‘Quantum-Optical Spectrometer’). Its first version was used on La Palma to study atmospheric scintillation on timescales between 100 milli-and 100 nanoseconds. For very high time resolution, light curves are of little use, and statistical functions of variability have to be measured. The noise in such functions decreases dramatically with increased light collecting power, making very large telescopes much more sensitive for the study of rapid variability than 2–4m class ones.

Context. To use the data in the future Gaia catalogue it is important to have accurate estimates of the statistical uncertainties and correlations of the errors in the astrometric data given in the catalogue. Aims. In a previous paper we derived a mathematical model for computing the covariances of the astrometric data based on series expansions and a simplified attitude description. The aim of the present paper is to determine to what extent this model provides an accurate representation of the expected random errors in the astrometric solution for Gaia. Methods. We simulate the astrometric core solution by making least-squares solutions of the astrometric parameters for one million stars and the attitude parameters for a five-year mission, using nearly one billion simulated elementary observations for a total of 26 million unknowns. Two cases are considered: one in which all stars have the same magnitude, and another with 30% brighter and 70% fainter stars. The resulting astrometric errors are statistically compared with the model predictions. Results. In all cases considered, and within the statistical uncertainties of the numerical experiments (typically below 0.4%), the theoretically calculated variances and covariances are consistent with the simulations. To achieve this it is however necessary to expand the covariances to at least third or fourth order, and to apply a (theoretically motivated and derived) "fudge factor" in the kinematographic model. Conclusions. The model provides a feasible method to estimate the covariance of arbitrary astrometric data, accurate enough for most applications, and as such it should be available as part of the user's interface to the Gaia catalogue. A main assumption in the current model is that the observational errors are uncorrelated (e.g., photon noise), and further studies are needed on how correlated modelling errors, in particular in the attitude, can be taken into account.

We describe a concept for an interferometric space mission dedicated to global (wide-angle) astrometry. The GAIA satellite contains two small (baseline APEQ 3 m) optical interferometers of the Fizeau type, mechanically set at a large and fixed angle to each other. Each interferometer has a field of view of about one degree. Continuous rotation of the whole satellite provides angular connections between the stars passing through the two fields of view. Positions, absolute parallaxes and annual proper motions can be determined with accuracies on the 20 micro-arcsec level. The observing programme may consist of all objects to a limiting magnitude around V = 15-16, including 50 million stars. The GAIA concept, which has been proposed for a Cornerstone Mission within the European Space Agency's long-term science programme, is based on the same general principles as the very successful ESA Hipparcos mission, but takes advantage of the much higher resolution and efficiency permitted by interferometry and modern detector techniques.

The Gaia space astrometry mission (to be launched in 2012) will use a continuously spinning spacecraft to construct a global system of positions, proper motions and absolute parallaxes from relative position measurements made in an astrometric focal plane. This astrometric reduction can be cast as a classical least-squares problem, and the adopted baseline method for its solution uses a simple iteration algorithm. A potential weakness of this approach, as opposed to a direct solution, is that any finite number of iterations results in truncation errors that are difficult to quantify. Thus it is of interest to investigate alternative approaches, in particular the feasibility of a direct (non-iterative) solution. A simplified version of the astrometric reduction problem is studied in which the only unknowns are the astrometric parameters for a subset of the stars and the continuous three-axis attitude, thus neglecting further calibration issues. The specific design of the Gaia spacecraft and scanning law leads to an extremely large and sparse normal equations matrix. Elimination of the star parameters leads to a much smaller but less sparse system. We try different reordering schemes and perform symbolic Cholesky decomposition of this reduced normal matrix to study the fill-in for successively longer time span of simulated observations. Extrapolating to the full mission length, we conclude that a direct solution is not feasible with today’s computational capabilities. Other schemes, e.g., eliminating the attitude parameters or orthogonalizing the observation equations, lead to similar or even worse problems. This negative result appears to be a consequence of the strong spatial and temporal connectivity among the unknowns achieved by two superposed fields of view and the scanning law, features that are in fact desirable and essential for minimizing large-scale systematic errors in the Gaia reference frame. We briefly consider also an approximate decomposition method a la Hipparcos, but conclude that it is either sub-optimal or effectively leads to an iterative solution.