Flat rotation curves in disk galaxies represent the main evidence for large amounts of surrounding dark matter. Despite of the difficulty in identifying the dark matter contribution to the total mass density in our Galaxy, stellar kinematics, as tracer of gravitational potential, is the most reliable observable for gauging different matter components. This work tests the flatness of the MW rotation curve with a simple general relativistic model suitable to represent the geometry of a disk as a stationary axisymmetric dust metric at a sufficiently large distance from a central body. Circular velocities of unprecedented accuracy were derived from the Gaia DR2 data for a carefully selected sample of disk stars. We then fit these velocities to both the classical, i.e. including a dark matter halo, rotation curve model and a relativistic analogue, as derived form the solution of Einstein's equation. The GR-compliant MW rotational curve model results statistically indistinguishable from its state-of-the-art DM analogue. This supports our ansatz that a stationary and axisymmetric galaxy-scale metric could "fill the gap" in a baryons-only Milky Way, suggestive of star orbits dragged along the background geometry. We confirmed that geometry is a manifestation of gravity according to the Einstein theory, in particular the weak gravitational effect due to the off-diagonal term of the metric could mimic for a "DM-like" effect in the observed flatness of the MW rotation curve. In the context of Local Cosmology, our findings are suggestive of a Galaxy phase-space as the exterior gravitational field of a Kerr-like source (inner rotating bulge) without the need of extra-matter.

We present new results from the Parallaxes of Southern Extremely Cool dwarfs program to measure parallaxes, proper motions and multi-epoch photometry of L and early T dwarfs. The observations were made on 108 nights over the course of 8 years using the Wide Field Imager on the ESO 2.2m telescope. We present 118 new parallaxes of L \& T dwarfs of which 52 have no published values and 24 of the 66 published values are preliminary estimates from this program. The parallax precision varies from 1.0 to 15.5 mas with a median of 3.8 mas. We find evidence for 2 objects with long term photometric variation and 24 new moving group candidates. We cross-match our sample to published photometric catalogues and find standard magnitudes in up to 16 pass-bands from which we build spectral energy distributions and H-R diagrams. This allows us to confirm the theoretically anticipated minimum in radius between stars and brown dwarfs across the hydrogen burning minimum mass. We find the minimum occurs between L2 and L6 and verify the predicted steep dependence of radius in the hydrogen burning regime and the gentle rise into the degenerate brown dwarf regime. We find a relatively young age of 2 Gyr from the kinematics of our sample.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Gaia Data Release 2 (Gaia DR2) contains results for 1693 million sources in the magnitude range 3 to 21 based on observations collected by the European Space Agency Gaia satellite during the first 22 months of its operational phase. We describe the input data, models, and processing used for the astrometric content of Gaia DR2, and the validation of these results performed within the astrometry task. Some 320 billion centroid positions from the pre-processed astrometric CCD observations were used to estimate the five astrometric parameters (positions, parallaxes, and proper motions) for 1332 million sources, and approximate positions at the reference epoch J2015.5 for an additional 361 million mostly faint sources. Special validation solutions were used to characterise the random and systematic errors in parallax and proper motion. For the sources with five-parameter astrometric solutions, the median uncertainty in parallax and position at the reference epoch J2015.5 is about 0.04 mas for bright (G<14 mag) sources, 0.1 mas at G=17 mag, and 0.7 mas at G=20 mag. In the proper motion components the corresponding uncertainties are 0.05, 0.2, and 1.2 mas/yr, respectively. The optical reference frame defined by Gaia DR2 is aligned with ICRS and is non-rotating with respect to the quasars to within 0.15 mas/yr. From the quasars and validation solutions we estimate that systematics in the parallaxes depending on position, magnitude, and colour are generally below 0.1 mas, but the parallaxes are on the whole too small by about 0.03 mas. Significant spatial correlations of up to 0.04 mas in parallax and 0.07 mas/yr in proper motion are seen on small (<1 deg) and intermediate (20 deg) angular scales. Important statistics and information for the users of the Gaia DR2 astrometry are given in the appendices.

A key objective of the ESA Gaia satellite is the realization of a quasi-inertial reference frame at visual wavelengths by means of global astrometric techniques. This requires an accurate mathematical and numerical modeling of relativistic light propagation, as well as double-blind-like procedures for the internal validation of the results, before they are released to the scientific community at large. Aim of this work is to specialize the Time Transfer Functions (TTF) formalism to the case of the Gaia observer and prove its applicability to the task of Global Sphere Reconstruction (GSR), in anticipation of its inclusion in the GSR system, already featuring the suite of RAMOD models, as an additional semi-external validation of the forthcoming Gaia baseline astrometric solutions. We extend the current GSR framework and software infrastructure (GSR2) to include TTF relativistic observation equations compatible with Gaia's operations. We use simulated data generated by the Gaia Data Reduction and Analysis Consortium (DPAC) to obtain different least-squares estimations of the full stellar spheres and gauge results. These are compared to analogous solutions obtained with the current RAMOD model in GSR2 and to the catalog generated with GREM, the model baselined for Gaia and used to generate the DPAC synthetic data. Linearized least-squares TTF solutions are based on spheres of about 132,000 primary stars uniformly distributed on the sky and simulated observations spanning the entire 5-yr range of Gaia's nominal operational lifetime. The statistical properties of the results compare well with those of GREM. Finally, comparisons to RAMOD@GSR2 solutions confirmed the known lower accuracy of that model and allowed us to establish firm limits on the quality of the linearization point outside of which an iteration for non-linearity is required for its proper convergence.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

We use methods of differential astrometry to construct a small field inertial reference frame stable at the micro-arcsecond level. Such a high level of astrometric precision can be expected with the end-of-mission standard errors to be achieved with the Gaia space satellite using global astrometry. We harness Gaia measurements of field angles and look at the influence of the number of reference stars and the star's magnitude as well as astrometric systematics on the total error budget with the help of Gaia-like simulations around the Ecliptic Pole in a differential astrometric scenario. We find that the systematic errors are modeled and reliably estimated to the $\mu$as level even in fields with a modest number of 37 stars with G $<$13 mag over a 0.24 sq.degs. field of view for short time scales of the order of a day with high-cadence observations such as those around the North Ecliptic Pole during the EPSL scanning mode of Gaia for a perfect instrument. The inclusion of the geometric instrument model over such short time scales accounting for large-scale calibrations requires fainter stars down to G = 14 mag without diminishing the accuracy of the reference frame. We discuss several future perspectives of utilizing this methodology over different and longer timescales.

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.

Given the extreme accuracy of modern space astrometry, a precise relativistic modeling of observations is required. Moreover, the availability of several models formulated in different and independent ways is a security against the presence of systematic errors in the analysis of future experimental results, like in the case of the Gaia mission. In this work, we simulate a series of observations using the two models to be used for the data analysis of Gaia, the Gaia RElativistic Model (GREM) and the Relativistic Astrometric MODel (RAMOD), and we compare them with the results of our astrometric model based on the Time Transfer Functions.

We review the mathematical models available for relativistic astrometry, discussing the different approaches and their accuracies in the context of the modern experiments from space like Gaia and GAME, and we show how these models can be applied to the real world, and their consequences from the mathematical and numerical point of view, with specific reference to the case of Gaia, whose launch is due before the end of the year.

The discovery and subsequent detailed study of T dwarfs has provided many surprises and pushed the physics and modeling of cool atmospheres in unpredicted directions. Distance is a critical parameter for studies of these objects to determine intrinsic luminosities, test binarity and measure their motion in the Galaxy. We describe a new observational program to determine distances across the full range of T dwarf sub-types using the NTT/SOFI telescope/instrument combination. We present preliminary results for ten objects, five of which represent new distances.

Advancement in astronomical observations and technical instrumentation implies taking into account the general relativistic effects due the gravitational fields encountered by the light while propagating from the star to the observer. Therefore, data exploitation for Gaia-like space astrometric mission (ESA, launch 2013) requires a fully relativistic interpretation of the inverse ray-tracing problem, namely the development of a highly accurate astrometric models in accordance with the geometrical environment affecting light propagation itself and the precepts of the theory of measurement. This could open a new rendition of the stellar distances and proper motions, or even an alternative detection perspective of many subtle relativistic effects suffered by light while it is propagating and subsequently recorded in the physical measurements.

Advancement in astronomical observations and technical instrumentation requires coding light propagation at high level of precision; this could open a new detection window of many subtle relativistic effects suffered by light while it is propagating and entangled in the physical measurements. Light propagation and its subsequent detection should indeed be conceived in a fully relativistic context, in order to interpret the results of the observations in accordance with the geometrical environment affecting light propagation itself, as an unicum surrounding universe. One of the most intriguing aspects is the boost towards the development of highly accurate models able to recon- struct the light path consistently with General Relativity and the precepts of measurements. This paper deals with the complexity of such a topic by showing how the geometrical framework of models like RAMOD, initially developed for astrometric observations, constitutes an appropriate physical environment for back tracing a light ray conforming to the intrinsic accuracy of space-time. This article discusses the reasons why RAMOD stands out among the existent approaches applied to the light propagation problem and provides a proof of its capability in recasting recent literature cases.

We present results from the PARallaxes of Southern Extremely Cool objects (PARSEC) program, an observational program begun in April 2007 to determine parallaxes for 122 L and 28 T southern hemisphere dwarfs using the Wide Field Imager on the ESO 2.2m telescope. The results presented here include parallaxes of 10 targets from observations over 18 months and a first version proper motion catalog. The proper motions were obtained by combining PARSEC observations astrometrically reduced with respect to the UCAC2 Catalog, and the 2MASS Catalog. The resulting median proper motion precision is 5mas/yr for 195,700 sources. The 140 0.3deg2 fields sample the southern hemisphere in an unbiased fashion with the exception of the galactic plane due to the small number of targets in that region. We present preliminary parallaxes with a 4.2 mas median precision for 10 brown dwarfs, 2 of which are within 10pc. These increase by 20% the present number of L dwarfs with published parallaxes. Of the 10 targets, 7 have been previously discussed in the literature: two were thought to be binary but the PARSEC observations show them to be single, one has been confirmed as a binary companion and another has been found to be part of a binary system, both of which will make good benchmark systems. Observations for the PARSEC program will end in early 2011 providing 3-4 years of coverage for all targets. The main expected outputs are: more than a 100% increase of the number of L dwarfs with parallaxes; to increment - in conjuction with published results - to at least 10 the number of objects per spectral subclass up to L9, and; to put sensible limits on the general binary fraction of brown dwarfs. We aim to contribute significantly to the understanding of the faint end of the H-R diagram and of the L/T transition region.

Modern astrometry is based on angular measurements at the micro-arcsecond level. At this accuracy a fully general relativistic treatment of the data reduction is required. This paper concludes a series of articles dedicated to the problem of relativistic light propagation, presenting the final microarcsecond version of a relativistic astrometric model which enable us to trace back the light path to its emitting source throughout the non-stationary gravity field of the moving bodies in the Solar System. The previous model is used as test-bed for numerical comparisons to the present one. Here we also test different versions of the computer code implementing the model at different levels of complexity to start exploring the best trade-off between numerical efficiency and the micro-arcsecond accuracy needed to be reached.

The detectors designed for Gaia, the next ESA space astrometry mission to be launched in 2011, will allow to observe repeatedly stars very close to Jupiter's limb. This will open a unique opportunity to test General Relativity by performing many Eddington-like experiments through the comparison between the pattern of a starfield observed with or without Jupiter. We have derived the main formulas relevant for the monopole and quadrupole light deflection by an oblate planet and developed a simulator to investigate the processing of the Gaia astrometric observation in the vicinity of the planet. The results show that such an experiment carried out with the Gaia data will provide a new fully independent determination of the PPN parameter gamma by means of differential astrometric measurements and, more importantly, for the first time will evidence the bending effect due to the quadrupole moment with a 3-sigma confidence level. Given the accuracy of the experiment for the monopole deflection, this will permit to test alternative modelling of the light bending by moving masses.

We develop here a new approach for the relativistic modeling of the photons moving into a quasi-Minkowskian space-time, where the metric is generated by an arbitrary n-body distribution within an isolated Solar System. Our model is built on the prescriptions of the theory of General Relativity and leaves the choices of the metric, as well as that of the motion of the observer, arbitrary. Adopting a quasi-Minkowskian expression of the metric accurate to order (v/c)^2, a thorough numerical test campaign is conducted to verify correctness and reliability of the model equations. The test results show that the model behaves according to predictions. Specifically, comparisons to true (simulated) data demonstrate that stellar distances are reconstructed up to the specified level of accuracy. Although the (v/c)^2 approximation is not always sufficient for its application to future astrometric experiments, which require modeling to (v/c)^3, this work serves also as a natural test-ground for the higher order model, whose formulation is now close to completion, and will be presented in a forthcoming paper.

In this article we outline the structure of a general relativistic astrometric model which has been developed to deduce the position and proper motion of stars from 1-microarcsecond optical observations made by an astrometric satellite orbiting around the Sun. The basic assumption of our model is that the Solar System is the only source of gravity, hence we show how we modeled the satellite observations in a many-body perturbative approach limiting ourselves to the order of accuracy of $(v/c)^2$. The microarcsecond observing scenario outlined is that for the GAIA astrometric mission.

The global astrometric observations of a GAIA-like satellite were modeled within the PPN formulation of Post-Newtonian gravitation. An extensive experimental campaign based on realistic end-to-end simulations was conducted to establish the sensitivity of global astrometry to the PPN parameter \gamma, which measures the amount of space curvature produced by unit rest mass. The results show that, with just a few thousands of relatively bright, photometrically stable, and astrometrically well behaved single stars, among the ~10^9 objects that will be observed by GAIA, \gamma can be estimated after 1 year of continuous observations with an accuracy of ~10^{-5} at the 3\sigma level. Extrapolation to the full 5-year mission of these results based on the scaling properties of the adjustment procedure utilized suggests that the accuracy of \simeq 2x10^{-7}, at the same 3\sigma level, can be reached with \~10^6 single stars, again chosen as the most astrometrically stable among the millions available in the magnitude range V=12-13. These accuracies compare quite favorably with recent findings of scalar-tensor cosmological models, which predict for \gamma a present-time deviation, |1-\gamma|, from the General Relativity value between 10^{-5} and 10^{-7}.