We present an analysis of the kinematics of Galactic disc white dwarf stars in the Solar neighbourhood using data from Gaia Data Release 2. Selection of white dwarfs based on parallax provides the first large, kinematically unbiased sample of Solar neighbourhood white dwarfs to date. Various classical properties of the Solar neighbourhood kinematics have been detected for the first time in the WD population. The disc white dwarf population exhibits a correlation between absolute magnitude and mean age, which we exploit to obtain an independent estimate of the Solar motion with respect to the Local Standard of Rest. This is found to be $(U,V,W)_{\odot} = (9.5\pm1.2, 7.5\pm1.2, 8.2\pm1.2)$ kms$^{-1}$. The $UW$ components agree with studies based on main sequence stars, however the $V$ component differs and may be affected by systematics arising from metallicity gradients in the disc. The velocity ellipsoid is shown to vary strongly with magnitude, and exhibits a significant vertex deviation in the $UV$ plane of around 15 degrees, due to the non-axisymmetric Galactic potential. The results of this study provide an important input to proper motion surveys for white dwarfs, which require knowledge of the velocity distribution in order to correct for missing low velocity stars that are culled from the sample to reduce subdwarf contamination.

We use 156 044 white dwarf candidates with $\geq5\sigma$ significant parallax measurements from the Gaia mission to measure the velocity dispersion of the Galactic disc; $(\sigma_U,\sigma_V,\sigma_W) = (30.8, 23.9, 20.0)$ km s$^{-1}$. We identify 142 objects that are inconsistent with disc membership at the $>5\sigma$ level. This is the largest sample of field halo white dwarfs identified to date. We perform a detailed model atmosphere analysis using optical and near-infrared photometry and parallaxes to constrain the mass and cooling age of each white dwarf. The white dwarf cooling ages of our targets range from 7 Myr for J1657+2056 to 10.3 Gyr for J1049-7400. The latter provides a firm lower limit of 10.3 Gyr for the age of the inner halo based on the well-understood physics of white dwarfs. Including the pre-white dwarf evolutionary lifetimes, and limiting our sample to the recently formed white dwarfs with cooling ages of $<500$ Myr, we estimate an age of $10.9 \pm 0.4$ Gyr (internal errors only) for the Galactic inner halo. The coolest white dwarfs in our sample also give similar results. For example, J1049-7400 has a total age of 10.9-11.1 Gyr. Our age measurements are consistent with other measurements of the age of the inner halo, including the white dwarf based measurements of the globular clusters M4, NGC 6397, and 47 Tuc.

A large sample of white dwarfs is selected by both proper motion and colours from the Pan-STARRS 1 3{\pi} Steradian Survey Processing Version 2 to construct the White Dwarf Luminosity Functions of the discs and halo in the solar neighbourhood. Four-parameter astrometric solutions were recomputed from the epoch data. The generalised maximum volume method is then used to calculate the density of the populations. After removal of crowded areas near the Galactic plane and centre, the final sky area used by this work is 7.833 sr, which is 83% of the 3{\pi} sky and 62% of the whole sky. By dividing the sky using Voronoi tessellation, photometric and astrometric uncertainties are recomputed at each step of the integration to improve the accuracy of the maximum volume. Interstellar reddening is considered throughout the work. We find a disc-to-halo white dwarf ratio of about 100.

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 use Gaia Data Release 2 to identify 13,928 white dwarfs within 100 pc of the Sun. The exquisite astrometry from Gaia reveals for the first time a bifurcation in the observed white dwarf sequence in both Gaia and the Sloan Digital Sky Survey (SDSS) passbands. The latter is easily explained by a helium atmosphere white dwarf fraction of 36%. However, the bifurcation in the Gaia colour-magnitude diagram depends on both the atmospheric composition and the mass distribution. We simulate theoretical colour-magnitude diagrams for single and binary white dwarfs using a population synthesis approach and demonstrate that there is a significant contribution from relatively massive white dwarfs that likely formed through mergers. These include white dwarf remnants of main-sequence (blue stragglers) and post-main sequence mergers. The mass distribution of the SDSS subsample, including the spectroscopically confirmed white dwarfs, also shows this massive bump. This is the first direct detection of such a population in a volume-limited sample.

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.

The European Space Agency Gaia satellite was launched into orbit around L2 in December 2013. This ambitious mission has strict requirements on residual systematic errors resulting from instrumental corrections in order to meet a design goal of sub-10 microarcsecond astrometry. During the design and build phase of the science instruments, various critical calibrations were studied in detail to ensure that this goal could be met in orbit. In particular, it was determined that the video-chain offsets on the analogue side of the analogue-to-digital conversion electronics exhibited instabilities that could not be mitigated fully by modifications to the flight hardware. We provide a detailed description of the behaviour of the electronic offset levels on microsecond timescales, identifying various systematic effects that are known collectively as offset non-uniformities. The effects manifest themselves as transient perturbations on the gross zero-point electronic offset level that is routinely monitored as part of the overall calibration process. Using in-orbit special calibration sequences along with simple parametric models, we show how the effects can be calibrated, and how these calibrations are applied to the science data. While the calibration part of the process is relatively straightforward, the application of the calibrations during science data processing requires a detailed on-ground reconstruction of the readout timing of each charge-coupled device (CCD) sample on each device in order to predict correctly the highly time-dependent nature of the corrections. We demonstrate the effectiveness of our offset non-uniformity models in mitigating the effects in Gaia data. We demonstrate for all CCDs and operating instrument and modes on board Gaia that the video-chain noise-limited performance is recovered in the vast majority of science samples.

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 discuss the Gaia Data Release 1 (September 2016) and preliminary work on maximising the benefit for cool white dwarf (WD) science in advance of the full parallax catalogue which will appear around one year later in DR2. The Tycho catalogue is used in conjunction with the all-sky ground based astrometric/photometric SuperCOSMOS Sky Survey in order to identify candidate faint common proper motion objects to the Tycho stars. Gaia DR1 is supplemented by the Tycho-Gaia Astrometric Solution catalogue containing some 2 million parallaxes with Hipparcos-like precision for Tycho stars. While hotter, brighter WDs are present in Tycho, cooler examples are much rarer (if present at all) and CPM offers one method to infer precision distances for a statistically useful sample of these very faint WDs.

Context. This paper presents an overview of the photometric data that are part of the first Gaia data release. Aims. The principles of the processing and the main characteristics of the Gaia photometric data are presented. Methods. The calibration strategy is outlined briefly and the main properties of the resulting photometry are presented. Results. Relations with other broadband photometric systems are provided. The overall precision for the Gaia photometry is shown to be at the milli-magnitude level and has a clear potential to improve further in future releases.

The traditional Schmidt density estimator has been proven to be unbiased and effective in a magnitude-limited sample. Previously, efforts have been made to generalize it for populations with non-uniform density and proper motion-limited cases. This work shows that the then-good assumptions for a proper motion-limited sample are no longer sufficient to cope with modern data. Populations with larger differences in the kinematics as compared to the local standard of rest are most severely affected. We show that this systematic bias can be removed by treating the discovery fraction inseparable from the generalized maximum volume integrand. The treatment can be applied to any proper motion-limited sample with good knowledge of the kinematics. This work demonstrates the method through application to a mock catalogue of a white dwarf-only solar neighbourhood for various scenarios and compared against the traditional treatment using a survey with Pan-STARRS-like characteristics.

The termination in the white dwarf luminosity function is a standard diagnostic tool for measuring the total age of nearby stellar populations. In this paper, an algorithm is presented for inverting the full white dwarf luminosity function to obtain a maximum likelihood estimate of the time varying star formation rate of the host stellar population. Tests with synthetic data demonstrate that the algorithm converges over a wide class of underlying star formation rate forms. The algorithm successfully estimates the moving average star formation rate as a function of lookback time in the presence of realistic measurement noise, though suffers from degeneracies around discontinuities in the underlying star formation rate. The inversion results are most sensitive to the choice of white dwarf cooling models, with the models produced by different groups giving quite different results. The results are relatively insensitive to the progenitor metallicity, initial mass function, initial-final mass relation and ratio of H/He atmosphere white dwarfs. Application to two independent determinations of the Solar neighbourhood white dwarf luminosity function gives similar results. The star formation rate has a bimodal form, with broad peaks at 2-3 Gyr and 7-9 Gyr in the past, separated by a significant lull of magnitude 30-90% depending on choice of cooling models. The onset of star formation occurs around 8-10 Gyr ago. The total integrated star formation rate is ~0.014 stars/pc3 in the Solar neighbourhood, for stars more massive than 0.6M_{solar}.

We discuss the prospects for new deep, wide-angle surveys of the Galactic cool white dwarf populations using data from Pan-STARRS: the Panoramic Survey Telescope & Rapid Response System.

I present an algorithm for inverting the luminosity function for white dwarfs to obtain a maximum likelihood estimate of the star formation rate of the host stellar population. The algorithm is of the general class of Expectation Maximization, and involves iteratively improving an initial guess of the star formation rate. Tests show that the inversion results are quite sensitive to the assumed metallicity and initial mass function, but relatively insensitive to the initial-final mass relation and ratio of H/He atmosphere white dwarfs. Application to two independent determinations of the Solar neighbourhood white dwarf luminosity function gives similar results: the star formation rate is characterised by an early burst, and more recent peak at 2-3 Gyr in the past.

We present a magnitude and proper motion limited catalogue of ~10,000 white dwarf candidates, obtained from the SuperCOSMOS Sky Survey by means of reduced proper motion selection. This catalogue extends to magnitudes r~19.75 and proper motions as low as mu ~0.05"/yr, and covers nearly three quarters of the sky. Photometric parallaxes provide distance estimates accurate to ~50%. This catalogue is used to measure the luminosity functions for disk and spheroid white dwarfs, using strict velocity cuts to isolate subsamples belonging to each population. Disk luminosity functions measured in this manner are really a conglomerate of thin and thick disk objects, due to the significant velocity overlap between these populations. We introduce a new statistical approach to the stellar luminosity function for nearby objects that succesfully untangles the contributions from the different kinematic populations, without the need for stringent velocity cuts. This improves the statistical power by allowing all stars to contribute to the luminosity function, even at tangential velocities where the populations are indistinguishable. This method is particularly suited to white dwarfs, for which population discrimination by chemical tagging is not possible. We use this technique to obtain the first measurement of the thick disk white dwarf luminosity function, while also improving constraint on both the thin disk and spheroid luminosity functions. We find that the thin disk, thick disk and spheroid populations contribute to the local white dwarf density in roughly 79%/16%/5% proportions.

We present the techniques and early results of our program to measure the luminosity function for White Dwarfs in the SuperCOSMOS Sky Survey. Our survey covers over three quarters of the sky to a mean depth of I~19.2, and finds ~9,500 Galactic disk WD candidates on applying a conservative lower tangential velocity cut of 30kms^-1. Novel techniques introduced in this survey include allowing the lower proper motion limit to vary according to apparent magnitude, fully exploiting the accuracy of proper motion measurements to increase the sample size. Our luminosity function shows good agreement with that measured in similar works. We find a pronounced drop in the local number density of WDs at a M_bol~15.75, and an inflexion in the luminosity function at M_bol~12.

We present photographic B, R and I photometry, and optical and near-infrared spectroscopy, of a new ultracool white dwarf (UCWD) discovered in the SuperCOSMOS Sky Survey. The spectrum of SSSJ1556-0806 shows strong flux suppression due to the presence of collisionally induced absorption by molecular hydrogen (H2CIA), a feature characteristic of the cool, high density environments found in the atmospheres of ultracool white dwarfs. SSSJ1556-0806 therefore joins a list of <10 ultracool white dwarfs displaying extreme flux suppression. Synthetic model fitting suggests an effective temperature <3000K, which if true would make this one of the coolest white dwarfs currently known. We also exploit the similarity between the SEDs of SSSJ1556-0806 and the well-studied UCWD LHS 3250 to aid in the determination of the atmospheric parameters in a regime where models consistently fail to reproduce observations. SSSJ1556-0806 is relatively bright (R ~ 17.8), making it particularly amenable to follow up observations to obtain trigonometric parallax and IR photometry.