No Arabic abstract
We analyze the volume-limited nearly complete 100 pc sample of 95 halo white dwarf candidates identified by the second data release of Gaia. Based on a detailed population synthesis model, we apply a method that relies on Gaia astrometry and photometry to accurately derive the individual white dwarf parameters (mass, radius, effective temperature, bolometric luminosity and age). This method is tested with 25 white dwarfs of our sample for which we took optical spectra and performed spectroscopic analysis. We build and analyse the halo white dwarf luminosity function, for which we find for the first time possible evidences of the cut-off at its faintest end, leading to an age estimate of $simeq12pm0.5 $Gyr. The mass distribution of the sample peaks at $0.589,M_{odot}$, with $71%$ of the white dwarf masses below $0.6,M_{odot}$ and just two massive white dwarfs of more than $0.8,M_{odot}$. From the age distribution we find three white dwarfs with total ages above 12 Gyr, of which J1312-4728 is the oldest white dwarf known with an age of $12.41pm0.22 $Gyr. We prove that the star formation history is mainly characterised by a burst of star formation that occurred from 10 to 12 Gyr in the past, but extended up to 8 Gyr. We also find that the peak of the star formation history is centered at around 11 Gyr, which is compatible with the current age of the Gaia-Enceladus encounter. Finally, $13%$ of our halo sample is contaminated by high-speed young objects (total age<7 Gyr). The origin of these white dwarfs is unclear but their age distribution may be compatible with the encounter with the Sagittarius galaxy.
As the remnants of stars with initial masses $lesssim$ 8 M$_{odot}$, white dwarfs contain valuable information on the formation histories of stellar populations. In this paper, we use deep, high-quality, u-band photometry from the Canada France Imaging Survey (CFIS), griz photometry from Pan-STARRS 1 (PS1), as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over $sim$4500 deg$^2$ using a reduced proper motion diagram. We develop a new white dwarf population synthesis code that returns mock observations of the Galactic field white dwarf population for a given star formation history, while simultaneously taking into account the geometry of the Milky Way, survey parameters, and selection effects. We use this model to derive the star formation histories of the thin disk, thick disk, and stellar halo. Our results show that the Milky Way disk began forming stars (11.3 $pm$ 0.5) Gyr ago, with a peak rate of (8.8 $pm$ 1.4) M$_{odot}$yr$^{-1}$ at (9.8 $pm$ 0.4) Gyr, before a slow decline to a constant rate until the present day --- consistent with recent results suggesting a merging event with a satellite galaxy. Studying the residuals between the data and best-fit model shows evidence for a slight increase in star formation over the past 3 Gyr. We fit the local fraction of helium-atmosphere white dwarfs to be (21 $pm$ 3) %. Incorporating this methodology with data from future wide-field surveys such as LSST, Euclid, CASTOR, and WFIRST should provide an unprecedented view into the formation of the Milky Way at its earliest epoch through its white dwarfs.
White dwarfs are the fossils left by the evolution of low-and intermediate-mass stars, and have very long evolutionary timescales. This allows us to use them to explore the properties of old populations, like the Galactic halo. We present a population synthesis study of the luminosity function of halo white dwarfs, aimed at investigating which information can be derived from the currently available observed data. We employ an up-to-date population synthesis code based on Monte Carlo techniques, that incorporates the most recent and reliable cooling sequences for metal poor progenitors as well as an accurate modeling of the observational biases. We find that because the observed sample of halo white dwarfs is restricted to the brightest stars only the hot branch of the white dwarf luminosity function can be used for such purposes, and that its shape function is almost insensitive to the most relevant inputs, like the adopted cooling sequences, the initial mass function, the density profile of the stellar spheroid, or the adopted fraction of unresolved binaries. Moreover, since the cut-off of the observed luminosity has not been yet determined only lower limits to the age of the halo population can be placed. We conclude that the current observed sample of the halo white dwarf population is still too small to obtain definite conclusions about the properties of the stellar halo, and the recently computed white dwarf cooling sequences which incorporate residual hydrogen burning should be assessed using metal-poor globular clusters.
Gaia-DR2 has provided an unprecedented number of white dwarf candidates of our Galaxy. In particular, it is estimated that Gaia-DR2 has observed nearly 400,000 of these objects and close to 18,000 up to 100 pc from the Sun. This large quantity of data requires a thorough analysis in order to uncover their main Galactic population properties, in particular the thin and thick disk and halo components. Taking advantage of recent developments in artificial intelligence techniques, we make use of a detailed Random Forest algorithm to analyse an 8-dimensional space (equatorial coordinates, parallax, proper motion components and photometric magnitudes) of accurate data provided by Gaia-DR2 within 100 pc from the Sun. With the aid of a thorough and robust population synthesis code we simulated the different components of the Galactic white dwarf population to optimize the information extracted from the algorithm for disentangling the different population components. The algorithm is first tested in a known simulated sample achieving an accuracy of 85.3%. Our methodology is thoroughly compared to standard methods based on kinematic criteria demonstrating that our algorithm substantially improves previous approaches. Once trained, the algorithm is then applied to the Gaia-DR2 100 pc white dwarf sample, identifying 12,227 thin disk, 1,410 thick disk and 95 halo white dwarf candidates, which represent a proportion of 74:25:1, respectively. Hence, the numerical spatial densities are $(3.6pm0.4)times10^{-3},{rm pc^{-3}}$, $(1.2pm0.4)times10^{-3},{rm pc^{-3}}$ and $(4.8pm0.4)times10^{-5},{rm pc^{-3}}$ for the thin disk, thick disk and halo components, respectively. The populations thus obtained represent the most complete and volume-limited samples to date of the different components of the Galactic white dwarf population.
We develop a novel method to simultaneously determine the vertical potential, force and stellar $z-v_z$ phase space distribution function (DF) in our local patch of the Galaxy. We assume that the Solar Neighborhood can be treated as a one-dimensional system in dynamical equilibrium and directly fit the number density in the $z-v_z$ plane to what we call the Rational Linear DF (RLDF) model. This model can be regarded as a continuous sum of isothermal DFs though it has only one more parameter than the isothermal model. We apply our method to a sample of giant stars from Gaia Data Release 2 and show that the RLDF provides an excellent fit to the data. The well-known phase space spiral emerges in the residual map of the $z-v_z$ plane. We use the best-fit potential to plot the residuals in terms of the frequency and angle of vertical oscillations and show that the spiral maps into a straight line. From its slope, we estimate that the phase spirals were generated by a perturbation $sim540$ Myr years ago. We also determine the differential surface density as a function of vertical velocity dispersion, a.k.a. the vertical temperature distribution. The result is qualitatively similar to what was previously found for SDSS/SEGUE G dwarfs. Finally, we address parameter degeneracies and the validity of the 1D approximation. Particularly, the mid-plane density derived from a cold subsample, where the 1D approximation is more secure, is closer to literature values than that derived from the sample as a whole.
Very precise observational data are needed for studying the stellar cluster parameters (distance, reddening, age, metallicity) and cluster internal kinematics. In turn, these give us an insight into the properties of our Galaxy, for example, by giving us the ability to trace Galactic spiral structure, star formation rates and metallicity gradients. We investigated the available Gaia DR2 catalogue of 1229 open clusters and studied cluster distances, sizes and membership distributions in the 3D space. An appropriate analysis of the parallaxto-distance transformation problem is presented in the context of getting distances toward open clusters and estimating their sizes. Based on our investigation of the Gaia DR2 data we argue that, within 2 kpc, the inverse-parallax method gives comparable results (distances and sizes) as the Bayesian approach based on the exponentially decreasing volume density prior. Both of these methods show very similar dependence of the line-of-sight elongation of clusters (needle-like shapes resulting from the parallax uncertainties) on the distance. We also looked at a measure of elongations of the studied clusters and find the maximum distance of 665 pc at which a spherical fit still contains about half of the stellar population of a cluster. It follows from these results that the 3D structure of an open cluster cannot be properly studied beyond about 500 pc when using any of mentioned standard transformations of parallaxes to distances.