The Ophiuchus cloud complex is one of the best laboratories to study the earlier stages of the stellar and protoplanetary disc evolution. The wealth of accurate astrometric measurements contained in the Gaia Data Release 2 can be used to update the census of Ophiuchus member candidates. We seek to find potential new members of Ophiuchus and identify those surrounded by a circumstellar disc. We constructed a control sample composed of 188 bona fide Ophiuchus members. Using this sample as a reference we applied three different density-based machine learning clustering algorithms (DBSCAN, OPTICS, and HDBSCAN) to a sample drawn from the Gaia catalogue centred on the Ophiuchus cloud. The clustering analysis was applied in the five astrometric dimensions defined by the three-dimensional Cartesian space and the proper motions in right ascension and declination. The three clustering algorithms systematically identify a similar set of candidate members in a main cluster with astrometric properties consistent with those of the control sample. The increased flexibility of the OPTICS and HDBSCAN algorithms enable these methods to identify a secondary cluster. We constructed a common sample containing 391 member candidates including 166 new objects, which have not yet been discussed in the literature. By combining the Gaia data with 2MASS and WISE photometry, we built the spectral energy distributions from 0.5 to $22\microm$ for a subset of 48 objects and found a total of 41 discs, including 11 Class II and 1 Class III new discs. Density-based clustering algorithms are a promising tool to identify candidate members of star forming regions in large astrometric databases. If confirmed, the candidate members discussed in this work would represent an increment of roughly 40% of the current census of Ophiuchus.

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.

We want to study whether the astrometric and photometric accuracies obtained for the Carte du Ciel plates digitized with a commercial digital camera are high enough for scientific exploitation of the plates. We use a digital camera Canon EOS~5Ds, with a 100mm macrolens for digitizing. We analyze six single-exposure plates and four triple-exposure plates from the Helsinki zone of Carte du Ciel (+39 degr < delta < 47 degr). Each plate is digitized using four images, with a significant central area being covered twice for quality control purposes. The astrometric calibration of the digitized images is done with the data from the Gaia TGAS (Tycho-Gaia Astrometric Solution) of the first Gaia data release (Gaia DR1), Tycho-2, HSOY (Hot Stuff for One Year), UCAC5 (USNO CCD Astrograph Catalog), and PMA catalogs. The best astrometric accuracy is obtained with the UCAC5 reference stars. The astrometric accuracy for single-exposure plates is sigma(R.A.)=0.16" and sigma(Dec.)=0.15" expressed as a Gaussian deviation of the astrometric residuals. For triple-exposure plates the astrometric accuracy is sigma(R.A.)=0.12" and sigma(Dec.)=0.13". The 1-sigma uncertainty of photometric calibration is about 0.28 mag and 0.24 mag for single- and triple-exposure plates, respectively. We detect the photographic adjacency (Kostinsky) effect in the triple-exposure plates. We show that accuracies at least of the level of scanning machines can be achieved with a digital camera, without any corrections for possible distortions caused by our instrumental setup. This method can be used to rapidly and inexpensively digitize and calibrate old photographic plates enabling their scientific exploitation.

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.

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.

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. 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. 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 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.

The European Space Agency's Gaia satellite was launched into orbit around L2 in December 2013 with a payload containing 106 large-format scientific CCDs. The primary goal of the mission is to repeatedly obtain high-precision astrometric and photometric measurements of one thousand million stars over the course of five years. The scientific value of the down-linked data, and the operation of the onboard autonomous detection chain, relies on the high performance of the detectors. As Gaia slowly rotates and scans the sky, the CCDs are continuously operated in a mode where the line clock rate and the satellite rotation spin-rate are in synchronisation. Nominal mission operations began in July 2014 and the first data release is being prepared for release at the end of Summer 2016. In this paper we present an overview of the focal plane, the detector system, and strategies for on-orbit performance monitoring of the system. This is followed by a presentation of the performance results based on analysis of data acquired during a two-year window beginning at payload switch-on. Results for parameters such as readout noise and electronic offset behaviour are presented and we pay particular attention to the effects of the L2 radiation environment on the devices. The radiation-induced degradation in the charge transfer efficiency (CTE) in the (parallel) scan direction is clearly diagnosed; however, an extrapolation shows that charge transfer inefficiency (CTI) effects at end of mission will be approximately an order of magnitude less than predicted pre-flight. It is shown that the CTI in the serial register (horizontal direction) is still dominated by the traps inherent to the manufacturing process and that the radiation-induced degradation so far is only a few per cent. Finally, we summarise some of the detector effects discovered on-orbit which are still being investigated.

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.

The first release of astrometric data from Gaia is expected in 2016. It will contain the mean stellar positions and magnitudes from the first year of observations. For more than 100 000 stars in common with the Hipparcos Catalogue it will be possible to compute very accurate proper motions due to the time difference of about 24 years between the two missions. This Hundred Thousand Proper Motions (HTPM) project will be part of the first release. Our aim is to investigate how early Gaia data can be optimally combined with information from the Hipparcos Catalogue in order to provide the most accurate and reliable results for HTPM. The Astrometric Global Iterative Solution (AGIS) was developed to compute the astrometric core solution based on the Gaia observations and will be used for all releases of astrometric data from Gaia. We adapt AGIS to process Hipparcos data in addition to Gaia observations, and use simulations to verify and study the joint solution method. For the HTPM stars we predict proper motion accuracies between 14 and 134 muas/yr, depending on stellar magnitude and amount of Gaia data available. Perspective effects will be important for a significant number of HTPM stars, and in order to treat these effects accurately we introduce a scaled model of kinematics. We define a goodness-of-fit statistic which is sensitive to deviations from uniform space motion, caused for example by binaries with periods of 10-50 years. HTPM will significantly improve the proper motions of the Hipparcos Catalogue well before highly accurate Gaia- only results become available. Also, HTPM will allow us to detect long period binary and exoplanetary candidates which would be impossible to detect from Gaia data alone. The full sensitivity will not be reached with the first Gaia release but with subsequent data releases. Therefore HTPM should be repeated when more Gaia data become available.

Starting in 2013, Gaia will deliver highly accurate astrometric data, which eventually will supersede most other stellar catalogues in accuracy and completeness. It is, however, lim- ited to observations from magnitude 6 to 20 and will therefore not include the brightest stars. Nano-JASMINE, an ultrasmall Japanese astrometry satellite, will observe these bright stars, but with much lower accuracy. Hence, the Hipparcos catalogue from 1997 will likely remain the main source of accurate distances to bright nearby stars. We are investigating how this might be improved by optimally combining data from all three missions in a joint astrometric solu- tion. This would take advantage of the unique features of each mission: the historic bright-star measurements of Hipparcos, the updated bright-star observations of Nano-JASMINE, and the very accurate reference frame of Gaia. The long temporal baseline between the missions pro- vides additional benefits for the determination of proper motions and binary detection, which indirectly improve the parallax determination further. We present a quantitative analysis of the expected gains based on simulated data for all three missions.

Reliable timing calibration is essential for the accurate comparison of XMM-Newton light curves with those from other observatories, to ultimately use them to derive precise physical quantities. The XMM-Newton timing calibration is based on pulsar analysis. However, as pulsars show both timing noise and glitches, it is essential to monitor these calibration sources regularly. To this end, the XMM-Newton observatory performs observations twice a year of the Crab pulsar to monitor the absolute timing accuracy of the EPIC-pn camera in the fast Timing and Burst modes. We present the results of this monitoring campaign, comparing XMM-Newton data from the Crab pulsar (PSR B0531+21) with radio measurements. In addition, we use five pulsars (PSR J0537-69, PSR B0540-69, PSR B0833-45, PSR B1509-58 and PSR B1055-52) with periods ranging from 16 ms to 197 ms to verify the relative timing accuracy. We analysed 38 XMM-Newton observations (0.2-12.0 keV) of the Crab taken over the first ten years of the mission and 13 observations from the five complementary pulsars. All the data were processed with the SAS, the XMM-Newton Scientific Analysis Software, version 9.0. Epoch folding techniques coupled with \chi^{2} tests were used to derive relative timing accuracies. The absolute timing accuracy was determined using the Crab data and comparing the time shift between the main X-ray and radio peaks in the phase folded light curves. The relative timing accuracy of XMM-Newton is found to be better than 10^{-8}. The strongest X-ray pulse peak precedes the corresponding radio peak by 306\pm9 \mus, which is in agreement with other high energy observatories such as Chandra, INTEGRAL and RXTE. The derived absolute timing accuracy from our analysis is \pm48 \mus.

The Hipparcos mission (1989-1993) resulted in the first space-based stellar catalogue including measurements of positions, parallaxes and annual proper motions accurate to about one milli-arcsecond. More space astrometry missions will follow in the near future. The ultra-small Japanese mission Nano-JASMINE (launch in late 2013) will determine positions and annual proper motions with some milli-arcsecond accuracy. In mid 2013 the next-generation ESA mission Gaia will deliver some tens of micro-arcsecond accurate astrometric parameters. Until the final Gaia catalogue is published in early 2020 the best way of improving proper motion values is the combination of positions from different missions separated by long time intervals. Rather than comparing positions from separately reduced catalogues, we propose an optimal method to combine the information from the different data sets by making a joint astrometric solution. This allows to obtain good results even when each data set alone is insufficient for an accurate reduction. We demonstrate our method by combining Hipparcos and simulated Nano-JASMINE data in a joint solution. We show a significant improvement over the conventional catalogue combination.

The ESA space astrometry mission Gaia, planned to be launched in 2013, has been designed to make angular measurements on a global scale with micro-arcsecond accuracy. A key component of the data processing for Gaia is the astrometric core solution, which must implement an efficient and accurate numerical algorithm to solve the resulting, extremely large least-squares problem. The Astrometric Global Iterative Solution (AGIS) is a framework that allows to implement a range of different iterative solution schemes suitable for a scanning astrometric satellite. In order to find a computationally efficient and numerically accurate iteration scheme for the astrometric solution, compatible with the AGIS framework, we study an adaptation of the classical conjugate gradient (CG) algorithm, and compare it to the so-called simple iteration (SI) scheme that was previously known to converge for this problem, although very slowly. The different schemes are implemented within a software test bed for AGIS known as AGISLab, which allows to define, simulate and study scaled astrometric core solutions. After successful testing in AGISLab, the CG scheme has been implemented also in AGIS. The two algorithms CG and SI eventually converge to identical solutions, to within the numerical noise (of the order of 0.00001 micro-arcsec). These solutions are independent of the starting values (initial star catalogue), and we conclude that they are equivalent to a rigorous least-squares estimation of the astrometric parameters. The CG scheme converges up to a factor four faster than SI in the tested cases, and in particular spatially correlated truncation errors are much more efficiently damped out with the CG scheme.

The Gaia satellite will observe about one billion stars and other point-like sources. The astrometric core solution will determine the astrometric parameters (position, parallax, and proper motion) for a subset of these sources, using a global solution approach which must also include a large number of parameters for the satellite attitude and optical instrument. The accurate and efficient implementation of this solution is an extremely demanding task, but crucial for the outcome of the mission. We provide a comprehensive overview of the mathematical and physical models applicable to this solution, as well as its numerical and algorithmic framework. The astrometric core solution is a simultaneous least-squares estimation of about half a billion parameters, including the astrometric parameters for some 100 million well-behaved so-called primary sources. The global nature of the solution requires an iterative approach, which can be broken down into a small number of distinct processing blocks (source, attitude, calibration and global updating) and auxiliary processes (including the frame rotator and selection of primary sources). We describe each of these processes in some detail, formulate the underlying models, from which the observation equations are derived, and outline the adopted numerical solution methods with due consideration of robustness and the structure of the resulting system of equations. Appendices provide brief introductions to some important mathematical tools (quaternions and B-splines for the attitude representation, and a modified Cholesky algorithm for positive semidefinite problems) and discuss some complications expected in the real mission data.

This paper provides a description of the Java software framework which has been constructed to run the Astrometric Global Iterative Solution for the Gaia mission. This is the mathematical framework to provide the rigid reference frame for Gaia observations from the Gaia data itself. This process makes Gaia a self calibrated, and input catalogue independent, mission. The framework is highly distributed typically running on a cluster of machines with a database back end. All code is written in the Java language. We describe the overall architecture and some of the details of the implementation.

In recent years Java has matured to a stable easy-to-use language with the flexibility of an interpreter (for reflection etc.) but the performance and type checking of a compiled language. When we started using Java for astronomical applications around 1999 they were the first of their kind in astronomy. Now a great deal of astronomy software is written in Java as are many business applications. We discuss the current environment and trends concerning the language and present an actual example of scientific use of Java for high-performance distributed computing: ESA's mission Gaia. The Gaia scanning satellite will perform a galactic census of about 1000 million objects in our galaxy. The Gaia community has chosen to write its processing software in Java. We explore the manifold reasons for choosing Java for this large science collaboration. Gaia processing is numerically complex but highly distributable, some parts being embarrassingly parallel. We describe the Gaia processing architecture and its realisation in Java. We delve into the astrometric solution which is the most advanced and most complex part of the processing. The Gaia simulator is also written in Java and is the most mature code in the system. This has been successfully running since about 2005 on the supercomputer "Marenostrum" in Barcelona. We relate experiences of using Java on a large shared machine. Finally we discuss Java, including some of its problems, for scientific computing.

Brief outline of Science Operations Centre activities for Gaia.

The ESAC Gaia team engages in a form of eXtreme programming while the DPAC will follow a series of six month development cycles modeled on this approach. As a project within the European Space Agency the European Committee for Space Standardization (ECSS) standards are required. We present the bringing together of these realms.

Gaia is ESA's ambitious space astrometry mission the main objective of which is to astrometrically and spectro-photometrically map 1000 Million celestial objects (mostly in our galaxy) with unprecedented accuracy. The announcement of opportunity for the data processing will be issued by ESA late in 2006. The Gaia Data Processing and Analysis Consortium (DPAC) has been formed recently and is preparing an answer. The satellite will downlink close to 100 TB of raw telemetry data over 5 years. To achieve its required accuracy of a few 10s of Microarcsecond astrometry, a highly involved processing of this data is required. In addition to the main astrometric instrument Gaia will host a Radial Velocity instrument, two low-resolution dispersers for multi-color photometry and two Star Mappers. Gaia is a flying Giga Pixel camera. The various instruments each require relatively complex processing while at the same time being interdependent. We describe the overall composition of the DPAC and the envisaged overall architecture of the Gaia data processing system. We shall delve further into the core processing - one of the nine, so-called, coordination units comprising the Gaia processing system.

The EPIC pn CCD camera on board of XMM-Newton is designed to perform high throughput imaging and spectroscopy as well as high resolution timing observations in the energy range of 0.1-15 keV. A temporal resolution of milliseconds or microseconds, depending on the instrument mode and detector, is outstanding for CCD based X-ray cameras. In order to calibrate the different observing modes of the EPIC pn CCD, XMM-Newton observations of the pulsars PSR B1509-58, PSR B0540-69 and the Crab were performed during the calibration and performance verification phase. To determine the accuracy of the on board clock against Coordinated Universal Time (UTC), PSR B1509-58 was observed simultaneously with XMM-Newton and RXTE in addition. The paper summarizes the current status of the clock calibration.

We present 0.1-10 keV observations of Her X-1 obtained with the Low-Energy Spectrometer Concentrator instrument onboard the BeppoSAX satellite during the main on-state of the 35 day cycle. We confirm the existence of an intense 0.093 keV blackbody component and a broad emission feature at 0.94 keV. The pulse phase dependence of these components is similar, suggesting a common origin. This is most likely fluorescent excitation of moderately ionized ($\xi\sim$10--100) material located at the inner edge of the accretion disk.

We present a 0.2-10 keV spectrum of the low-mass X-ray binary Cygnus X-2 obtained using the Low Energy Concentrator Spectrometer on-board Bepposax. The spectrum can be described by a cut-off power-law model with absorption of (2.28 +/- 0.07) 10^{21} atoms/cm^2, a power-law index of 0.78 +/- 0.02 and a cut-off energy of 4.30 +/- 0.08 keV (68% confidence errors), except at energies near ~1 keV where excess emission is present. This can be modeled by a broad Gaussian line feature with an energy of 1.02 +/- 0.04 keV, a full width half-maximum of 0.47 +/- 0.07 keV and an equivalent width of 74 +/-^{25}_{11} eV. This result confirms earlier reports of line emission near 1 keV and shows the intensity and structure of the feature to be variable.