We present a spectroscopic survey of 230 white dwarf candidates within 40 pc of the Sun from the William Herschel Telescope and Gran Telescopio Canarias. All candidates were selected from $\textit{Gaia}$ Data Release 2 (DR2) and in almost all cases had no prior spectroscopic classifications. We find a total of 191 confirmed white dwarfs and 39 main-sequence star contaminants. The majority of stellar remnants in the sample are relatively cool ($\langle T_{\rm eff} \rangle$ = 6200 K), showing either hydrogen Balmer lines or a featureless spectrum, corresponding to 89 DA and 76 DC white dwarfs, respectively. We also recover two DBA white dwarfs and 9--10 magnetic remnants. We find two carbon-bearing DQ stars and 14 new metal-rich white dwarfs. This includes the possible detection of the first ultra-cool white dwarf with metal lines. We describe three DZ stars for which we find at least four different metal species, including one which is strongly Fe- and Ni-rich, indicative of the accretion of a planetesimal with core-Earth composition. We find one extremely massive (1.31 $\pm$ 0.01 M$_{\odot}$) DA white dwarf showing weak Balmer lines, possibly indicating stellar magnetism. Another white dwarf shows strong Balmer line emission but no infrared excess, suggesting a low-mass sub-stellar companion. High spectroscopic completeness ($>$99%) has now been reached for $\textit{Gaia}$ DR2 sources within 40 pc sample, in the northern hemisphere ($\delta >$ 0 deg) and located on the white dwarf cooling track in the Hertzsprung-Russell diagram. A statistical study of the full northern sample is presented in a companion paper.

We perform an analysis of the single white dwarf and the double degenerate binary populations in the solar neighbourhood following a population synthesis approach to investigate the effects of unresolved double degenerates in the white dwarf luminosity function. We consider all unresolved synthetic binaries to be associated with fictitious effective temperatures and surface gravities that are obtained in the same way as if these objects were observed as single point sources. We evaluate the effects of unresolved double white dwarfs assuming that the synthetic samples are observed both by the magnitude-limited SDSS and the volume-limited Gaia surveys, the latter limited to a distance of no more than 100pc. We find that, for our standard model, the impact of unresolved double degenerates in the white dwarf luminosity function derived from the Gaia sample is nearly negligible. Unresolved double degenerates are hence expected to have no effect on the age of the Galactic disc, nor on the star formation history from this population. However, for the SDSS sample, the effect of unresolved double degenerates is significant at the brighter bins (Mbol<6.5 mag), with the fraction of such systems reaching ~40% of the total white dwarf population at Mbol=6 mag. This indicates unresolved double degenerates may influence the constraints on the star formation history derived from the SDSS white dwarf sample.

We analyse the 100pc Gaia white dwarf volume-limited sample by means of VOSA (Virtual Observatory SED Analyser) with the aim of identifying candidates for displaying infrared excesses. Our search focuses on the study of the spectral energy distribution (SED) of 3,733 white dwarfs with reliable infrared photometry and GBP-GRP colours below 0.8 mag, a sample which seems to be nearly representative of the overall white dwarf population. Our search results in 77 selected candidates, 52 of which are new identifications. For each target we apply a two-component SED fitting implemented in VOSA to derive the effective temperatures of both the white dwarf and the object causing the excess. We calculate a fraction of infrared-excess white dwarfs due to the presence of a circumstellar disk of 1.6+-0.2%, a value which increases to 2.6+-0.3% if we take into account incompleteness issues. Our results are in agreement with the drop in the percentage of infrared excess detections for cool (<8,000K) and hot (>20,000K) white dwarfs obtained in previous analyses. The fraction of white dwarfs with brown dwarf companions we derive is ~0.1-0.2%.

Comets in the Oort cloud evolve under the influence of internal and external perturbations, such as giant planets, stellar passages, and the galactic tidal field. We aim to study the dynamical evolution of the comets in the Oort cloud, accounting for external perturbations (passing stars and the galactic tide). We first construct an analytical model of stellar encounters. We find that individual perturbations do not modify the dynamics of the comets in the cloud unless very close (< 0.5pc) encounters occur. Using proper motions, parallaxes, and radial velocities from Gaia DR2, we construct an astrometric catalogue of 14,659 stars that are within 50pc from the Sun. For all these stars we calculate the time and the closest distance to the Sun. We find that the cumulative effect of relatively distant ($\leq1$ pc) passing stars can perturb the comets in the Oort cloud. Finally, we study the dynamical evolution of the comets in the Oort cloud under the influence of multiple stellar encounters within 2.5pc from the Sun and the galactic tidal field over $\pm10$Myr. We considered two models for the Oort cloud, compact (a $\leq$0.25 pc) and extended (a$ \leq0.5$ pc). We find that the cumulative effect of stellar encounters is the major perturber of the Oort cloud for a compact configuration while for the extended, the galactic tide is the major perturber. In both cases, the effect of passing stars and the galactic tide raises the semi-major axis of $\sim1.1$\% of the comets at the edge of the cloud up to interstellar regions ($a >0.5$pc). This leads to the creation of transitional interstellar comets, which might become interstellar objects due to external perturbations. This raises the question about the existence of a cloud of objects in the interstellar space which might overlap with our Oort cloud if we consider that other planetary systems face similar processes for the ejection of comets.

We analyzed the velocity space of the thin and thick-disk Gaia white dwarf population within 100 pc looking for signatures of the Hercules stellar stream. We aimed to identify those objects belonging to the Hercules stream and, by taking advantage of white dwarf stars as reliable cosmochronometers, to derive a first age distribution. We applied a kernel density estimation to the $UV$ velocity space of white dwarfs. For the region where a clear overdensity of stars was found, we created a 5-D space of dynamic variables. We applied a hierarchichal clustering method, HDBSCAN, to this 5-D space, identifying those white dwarfs that share similar kinematic characteristics. Finally, under general assumptions and from their photometric properties, we derived an age estimate for each object. The Hercules stream was firstly revealed as an overdensity in the $UV$ velocity space of the thick-disk white dwarf population. Three substreams were then found: Hercules $a$ and Hercules $b$, formed by thick-disk stars with an age distribution peaked $4\,$Gyr in the past and extended to very old ages; and Hercules $c$, with a ratio of 65:35 thin:thick stars and a more uniform age distribution younger than 10 Gyr.

We analyze the effect of the sedimentation of $^{22}$Ne on the local white dwarf luminosity function by studying scenarios under different Galactic metallicity models. We make use of an up-to-date population synthesis code based on Monte Carlo techniques to derive the synthetic luminosity function. Constant solar metallicity models are not able to simultaneously reproduce the peak and cut-off of the white dwarf luminosity function. The extra release of energy due to $^{22}$Ne sedimentation piles up more objects in brighter bins of the faint end of the luminosity function. The contribution of a single burst thick disk population increases the number of stars in the magnitude interval centered around $M_{\rm bol}=15.75$. Among the metallicity models studied, the one following a Twarog's profile is disposable. Our best fit model was obtained when a dispersion in metallicities around the solar metallicity value is considered along with a $^{22}$Ne sedimentation model, a thick disk contribution and an age of the thin disk of $8.8\pm0.2$ Gyr. Our population synthesis model is able to reproduce the local white dwarf luminosity function with a high degree of precision when a dispersion in metallicities around the solar value model is adopted. Although the effects of $^{22}$Ne sedimentation are only marginal and the contribution of a thick disk population is minor, both of them help in better fitting the peak and the cut-off regions of the white dwarf luminosity function.

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.6\pm0.4)\times10^{-3}\,{\rm pc^{-3}}$, $(1.2\pm0.4)\times10^{-3}\,{\rm pc^{-3}}$ and $(4.8\pm0.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.

Cool subdwarfs are metal-poor low-mass stars that formed during the early stages of the evolution of our Galaxy. Because they are relatively rare in the vicinity of the Sun, we know of few cool subdwarfs in the solar neighbourhood, and none with both the mass and the radius accurately determined. This hampers our understanding of stars at the low-mass end of the main-sequence. Here we report the discovery of SDSSJ235524.29+044855.7 as an eclipsing binary containing a cool subdwarf star, with a white dwarf companion. From the light-curve and the radial-velocity curve of the binary we determine the mass and the radius of the cool subdwarf and we derive its effective temperature and luminosity by analysing its spectral energy distribution. Our results validate the theoretical mass-radius-effective temperature-luminosity relations for low-mass low-metallicity stars.

The Carnegie Supernova Project-II (CSP-II) was an NSF-funded, four-year program to obtain optical and near-infrared observations of a "Cosmology" sample of $\sim100$ Type Ia supernovae located in the smooth Hubble flow ($0.03 \lesssim z \lesssim 0.10$). Light curves were also obtained of a "Physics" sample composed of 90 nearby Type Ia supernovae at $z \leq 0.04$ selected for near-infrared spectroscopic time-series observations. The primary emphasis of the CSP-II is to use the combination of optical and near-infrared photometry to achieve a distance precision of better than 5%. In this paper, details of the supernova sample, the observational strategy, and the characteristics of the photometric data are provided. In a companion paper, the near-infrared spectroscopy component of the project is presented.

We present an analysis of the final data release of the Carnegie Supernova Project I, focusing on the absolute calibration of the luminosity-decline-rate relation for Type Ia supernovae (SNeIa) using new intrinsic color relations with respect to the color-stretch parameter, $s_{BV}$, enabling improved dust extinction corrections. We investigate to what degree the so-called fast-declining SNeIa can be used to determine accurate extragalactic distances. We estimate the intrinsic scatter in the luminosity-decline-rate relation, and find it ranges from $\pm 0.13$ mag to $\pm 0.18$ mag with no obvious dependence on wavelength. Using the Cepheid variable star data from the SH0ES project (Riess et al., 2016), the SNIa distance scale is calibrated and the Hubble constant is estimated using our optical and near-infrared sample, and these results are compared to those determined exclusively from a near-infrared sub-sample. The systematic effect of the supernova's host galaxy mass is investigated as a function of wavelength and is found to decrease toward redder wavelengths, suggesting this effect may be due to dust properties of the host. Using estimates of the dust extinction derived from optical and NIR wavelengths, and applying these to H band, we derive a Hubble constant $H_0 = 73.2 \pm 2.3$ km/s/Mpc, whereas using a simple $B-V$ color-correction applied to B band yields $H_0 = 72.7 \pm 2.1$ km/s/Mpc. Photometry of two calibrating SNeIa from the CSP-II sample, SN2012ht and SN2015F, is presented and used to improve the calibration of the SNIa distance ladder.

Double white dwarf binaries with merger timescales smaller than the Hubble time and with a total mass near the Chandrasekhar limit (i.e. classical Chandrasekhar population) or with high-mass primaries (i.e. sub-Chandrasekhar population) are potential supernova type Ia (SNIa) progenitors. However, we have not yet unambiguously confirmed the existence of these objects observationally, a fact that has been often used to criticise the relevance of double white dwarfs for producing SNIa. We analyse whether this lack of detections is due to observational effects. To that end we simulate the double white dwarf binary population in the Galaxy and obtain synthetic spectra for the SNIa progenitors. We demonstrate that their identification, based on the detection of Halpha double-lined profiles arising from the two white dwarfs in the synthetic spectra, is extremely challenging due to their intrinsic faintness. This translates into an observational probability of finding double white dwarf SNIa progenitors in the Galaxy of (2.1+-1.0)x10^{-5} and (0.8+-0.4)x10^{-5} for the classical Chandrasekhar and the sub-Chandrasekhar progenitor populations, respectively. Eclipsing double white dwarf SNIa progenitors are found to suffer from the same observational effect. The next generation of large-aperture telescopes are expected to help in increasing the probability for detection by ~1 order of magnitude. However, it is only with forthcoming observations such as those provided by LISA that we expect to unambiguously confirm or disprove the existence of double white dwarf SNIa progenitors and to test their importance for producing SNIa.

We present a catalogue of 73,221 white dwarf candidates extracted from the astrometric and photometric data of the recently published Gaia DR2 catalogue. White dwarfs were selected from the Gaia Hertzsprung-Russell diagram with the aid of the most updated population synthesis simulator. Our analysis shows that Gaia has virtually identified all white dwarfs within 100 pc from the Sun. Hence, our sub-population of 8,555 white dwarfs within this distance limit and the colour range considered, $-\,0.52<(G_{\rm BP}-G_{\rm RP})<0.80$, is the largest and most complete volume-limited sample of such objects to date. From this sub-sample we identified 8,343 CO-core and 212 ONe-core white dwarf candidates and derived a white dwarf space density of $4.9\pm0.4\times10^{-3}\,{\rm pc^{-3}}$. A bifurcation in the Hertzsprung-Russell diagram for these sources, which our models do not predict, is clearly visible. We used the Virtual Observatory tool VOSA to derive effective temperatures and luminosities for our sources by fitting their spectral energy distributions, that we built from the UV to the NIR using publicly available photometry through the Virtual Observatory. From these parameters, we derived the white dwarf radii. Interpolating the radii and effective temperatures in hydrogen-rich white dwarf cooling sequences, we derived the surface gravities and masses. The Gaia 100 pc white dwarf population is clearly dominated by cool ($\sim$ 8,000 K) objects and reveals a significant population of massive ($M \sim 0.8 M_{\odot}$) white dwarfs, of which no more than $\sim$ $30-40 \%$ can be attributed to hydrogen-deficient atmospheres, and whose origin remains uncertain.

We conduct a population synthesis study of the common envelope evolution (CEE) of a white dwarf (WD) and an asymptotic giant branch (AGB) star and find that the potential number of type Ia supernovae (SNe Ia) from the core degenerate (CD) and from the double degenerate (DD) scenarios are of the same order of magnitude. For the CD scenario we consider cases where a carbon oxygen rich (CO) AGB core and a CO WD merge during the CEE and leave a WD remnant with a mass of >1.35Mo, and for the DD scenario we count surviving CO WD binary systems that merge within a time of 10 Gyr. When either the AGB core or the WD are oxygen neon rich (ONe) we assume that the outcome might be a peculiar SN Ia. We find that the number of potential peculiar SNe Ia in the channels we study, that do not include peculiar SNe Ia that involve helium WDs or helium donors, is non-negligible, but less than the number of normal SNe Ia. If a SN Ia or a peculiar SN Ia explosion takes place within about million year after CEE, whether in the CD scenario or in the DD scenario, a massive circumstellar matter is present at explosion time. Our results are compatible with the suggestion that Chandrasekhar-mass SNe Ia mostly come from the CD scenario, and sub-Chandrasekhar SNe Ia mostly come from the DD scenario.

We study the origin of the interstellar object 1I/2017 U1 'Oumuamua by juxtaposing estimates based on the observations with simulations. We speculate that objects like 'Oumuamua are formed in the debris disc as left over from the star and planet formation process, and subsequently liberated. The liberation process is mediated either by interaction with other stars in the parental star-cluster, by resonant interactions within the planetesimal disc or by the relatively sudden mass loss when the host star becomes a compact object. Integrating backward in time in the Galactic potential together with stars from the Gaia-TGAS catalogue we find that about 1.3Myr ago 'Oumuamua passed the nearby star HIP 17288 within a mean distance of $1.3$pc. By comparing nearby observed L-dwarfs with simulations of the Galaxy we conclude that the kinematics of 'Oumuamua is consistent with relatively young objects of $1.1$--$1.7$Gyr. We just met 'Oumuamua by chance, and with a derived mean Galactic density of $\sim 3\times 10^{5}$ similarly sized objects within 100\,au from the Sun or $\sim 10^{14}$ per cubic parsec we expect about 2 to 12 such visitors per year within 1au from the Sun.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

The evolution of white dwarfs is a simple gravothermal process of cooling. Since the shape of their luminosity function is sensitive to the characteristic cooling time, it is possible to use its slope to test the existence of additional sources or sinks of energy, such as those predicted by alternative physical theories. The aim of this paper is to study if the changes in the slope of the white dwarf luminosity function around bolometric magnitudes ranging from 8 to 10 and previously attributed to axion emission are, effectively, a consequence of the existence of axions and not an artifact introduced by the star formatio rate. We compute theoretical luminosity functions of the thin and thick disk, and of the stellar halo including axion emission and we compare them with the existing observed luminosity functions. Since these stellar populations have different star formation histories, the slope change should be present in all of them at the same place if it is due to axions or any other intrinsic cooling mechanism. The signature of an unexpected cooling seems to be present in the luminosity functions of the thin and thick disks, as well as in the halo luminosity function. This additional cooling is compatible with axion emission, thus supporting the idea that DFSZ axions, with a mass in the range of 4 to 10 meV, could exist. If this were the case, these axions could be detected by the future solar axioscope IAXO.

Recent Hubble Space Telescope observations have unveiled the white dwarf cooling sequence of the Galactic bulge. Although the degenerate sequence can be well fitted employing the most up-to-date theoretical cooling sequences, observations show a systematic excess of red objects that cannot be explained by the theoretical models of single carbon-oxygen white dwarfs of the appropriate masses. Here we present a population synthesis study of the white dwarf cooling sequence of the Galactic bulge that takes into account the populations of both single white dwarfs and binary systems containing at least one white dwarf. These calculations incorporate state-of-the-art cooling sequences for white dwarfs with hydrogen-rich and hydrogen-deficient atmospheres, for both white dwarfs with carbon-oxygen and helium cores, and also take into account detailed prescriptions of the evolutionary history of binary systems. Our Monte Carlo simulator also incorporates all the known observational biases. This allows us to model with a high degree of realism the white dwarf population of the Galactic bulge. We find that the observed excess of red stars can be partially attributed to white dwarf plus main sequence binaries, and to cataclysmic variables or dwarf novae. Our best fit is obtained with a higher binary fraction and an initial mass function slope steeper than standard values, as well as with the inclusion of differential reddening and blending. Our results also show that the possible contribution of double degenerate systems or young and thick-disk bulge stars is negligible.

We present final natural system optical (ugriBV) and near-infrared (YJH) photometry of 134 supernovae (SNe) with probable white dwarf progenitors that were observed in 2004-2009 as part of the first stage of the Carnegie Supernova Project (CSP-I). The sample consists of 123 Type Ia SNe, 5 Type Iax SNe, 2 super-Chandrasekhar SN candidates, 2 Type Ia SNe interacting with circumstellar matter, and 2 SN 2006bt-like events. The redshifts of the objects range from z = 0.0037 to 0.0835; the median redshift is 0.0241. For 120 (90%) of these SNe, near-infrared photometry was obtained. Average optical extinction coefficients and color terms are derived and demonstrated to be stable during the five CSP-I observing campaigns. Measurements of the CSP-I near-infrared bandpasses are also described, and near-infrared color terms are estimated through synthetic photometry of stellar atmosphere models. Optical and near-infrared magnitudes of local sequences of tertiary standard stars for each supernova are given, and a new calibration of Y-band magnitudes of the Persson et al. (1998) standards in the CSP-I natural system is presented.

Exoplanet discoveries over recent years have shown that terrestrial planets are exceptionally common. Many of these planets are in compact systems that result in complex orbital dynamics. A key step toward determining the surface conditions of these planets is understanding the latitudinally dependent flux incident at the top of the atmosphere as a function of orbital phase. The two main properties of a planet that influence the time-dependent nature of the flux are the obliquity and orbital eccentricity of the planet. We derive the criterion for which the flux variation due to obliquity is equivalent to the flux variation due to orbital eccentricity. This equivalence is computed for both the maximum and average flux scenarios, the latter of which includes the effects of the diurnal cycle. We apply these calculations to four known multi-planet systems (GJ 163, K2-3, Kepler-186, and Proxima Centauri), where we constrain the eccentricity of terrestrial planets using orbital dynamics considerations and model the effect of obliquity on incident flux. We discuss the implications of these simulations on climate models for terrestrial planets and outline detectable signatures of planetary obliquity.

We present a Monte Carlo population synthesis study of white dwarf-main sequence (WD+MS) binaries in the Galactic disk aimed at reproducing the ensemble properties of the entire population observed by the Sloan Digital Sky Survey (SDSS) Data Release 12. Our simulations take into account all known observational biases and use the most up-to-date stellar evolutionary models. This allows us to perform a sound comparison between the simulations and the observational data. We find that the properties of the simulated and observed parameter distributions agree best when assuming low values of the common envelope efficiency (0.2-0.3), a result that is in agreement with previous findings obtained by observational and population synthesis studies of close SDSS WD+MS binaries. We also show that all synthetic populations that result from adopting an initial mass ratio distribution with a positive slope are excluded by observations. Finally, we confirm that the properties of the simulated WD+MS binary populations are nearly independent of the age adopted for the thin disk, on the contribution of WD+MS binaries from the thick disk (0-17 per cent of the total population) and on the assumed fraction of the internal energy that is used to eject the envelope during the common envelope phase (0.1-0.5).

We use the Sloan Digital Sky Survey Data Release 12, which is the largest available white dwarf catalog to date, to study the evolution of the kinematical properties of the population of white dwarfs in the Galactic disc. We derive masses, ages, photometric distances and radial velocities for all white dwarfs with hydrogen-rich atmospheres. For those stars for which proper motions from the USNO-B1 catalog are available the true three-dimensional components of the stellar space velocity are obtained. This subset of the original sample comprises 20,247 objects, making it the largest sample of white dwarfs with measured three-dimensional velocities. Furthermore, the volume probed by our sample is large, allowing us to obtain relevant kinematical information. In particular, our sample extends from a Galactocentric radial distance $R_{\rm G}=7.8$~kpc to 9.3~kpc, and vertical distances from the Galactic plane ranging from $Z=-0.5$~kpc to 0.5~kpc. We examine the mean components of the stellar three-dimensional velocities, as well as their dispersions with respect to the Galactocentric and vertical distances. We confirm the existence of a mean Galactocentric radial velocity gradient, $\partial\langle V_{\rm R}\rangle/\partial R_{\rm G}=-3\pm5$~km~s$^{-1}$~kpc$^{-1}$. We also confirm North-South differences in $\langle V_{\rm z}\rangle$. Specifically, we find that white dwarfs with $Z>0$ (in the North Galactic hemisphere) have $\langle V_{\rm z}\rangle<0$, while the reverse is true for white dwarfs with $Z<0$. The age-velocity dispersion relation derived from the present sample indicates that the Galactic population of white dwarfs may have experienced an additional source of heating, which adds to the secular evolution of the Galactic disc.

We present a set of full evolutionary sequences for white dwarfs with hydrogen-deficient atmospheres. We take into account the evolutionary history of the progenitor stars, all the relevant energy sources involved in the cooling, element diffusion in the very outer layers, and outer boundary conditions provided by new and detailed non-gray white dwarf model atmospheres for pure helium composition. These model atmospheres are based on the most up-to-date physical inputs. Our calculations extend down to very low effective temperatures, of $\sim 2\,500$~K, provide a homogeneous set of evolutionary cooling tracks that are appropriate for mass and age determinations of old hydrogen-deficient white dwarfs, and represent a clear improvement over previous efforts, which were computed using gray atmospheres.

Some globular clusters host multiple stellar populations with different chemical abundance patterns. This is particularly true for $\omega$ Centauri, which shows clear evidence of a helium- enriched sub-population characterized by a helium abundance as high as $Y= 0.4$. We present a whole and consistent set of evolutionary tracks from the ZAMS to the white dwarf stage appropriate for the study of the formation and evolution of white dwarfs resulting from the evolution of helium-rich progenitors. Different issues of the white dwarf evolution and their helium-rich progenitors have been explored. In particular, the final mass of the remnants, the role of overshooting during the thermally-pulsing phase, and the cooling of the resulting white dwarfs differ markedly from the evolutionary predictions of progenitor stars with standard initial helium abundance. Finally, the pulsational properties of the resulting white dwarfs are also explored. We find that, for the range of initial masses explored in this paper, the final mass of the helium-rich progenitors is markedly larger than the final mass expected from progenitors with the usual helium abundance. We also find that progenitors with initial mass smaller than $M_\star \simeq 0.65\,M_\odot$ evolve directly into helium-core white dwarfs in less than 14~Gyr, and that for larger progenitor masses the evolution of the resulting low-mass carbon-oxygen white dwarfs is dominated by residual nuclear burning. For helium-core white dwarfs, we find that they evolve markedly faster than their counterparts coming from standard progenitors. Also, in contrast with what occurs for white dwarfs resulting from progenitors with the standard helium abundance, the impact of residual burning on the cooling time of white dwarfs is not affected by the occurrence of overshooting during the thermally-pulsing phase of progenitor stars.

The white dwarf luminosity function is an important tool to understand the properties of the Solar neighborhood, like its star formation history, and its age. Here we present a population synthesis study of the white dwarf population within 40~pc from the Sun, and compare the results of this study with the properties of the observed sample. We use a state-of-the-art population synthesis code based on Monte Carlo techniques, that incorporates the most recent and reliable white dwarf cooling sequences, an accurate description of the Galactic neighborhood, and a realistic treatment of all the known observational biases and selection procedures. We find a good agreement between our theoretical models and the observed data. In particular, our simulations reproduce a previously unexplained feature of the bright branch of the white dwarf luminosity function, which we argue is due to a recent episode of star formation. We also derive the age of the Solar neighborhood employing the position of the observed cut-off of the white dwarf luminosity function, obtaining ~8.9+-0.2 Gyr. We conclude that a detailed description of the ensemble properties of the population of white dwarfs within 40pc of the Sun allows us to obtain interesting constraints on the history of the Solar neighborhood.

NGC 6397 is one of the most interesting, well observed and theoretically studied globular clusters. The existing wealth of observations allows us to study the reliability of the theoretical white dwarf cooling sequences of low metallicity progenitors,to determine its age and the percentage of unresolved binaries, and to assess other important characteristics of the cluster, like the slope of the initial mass function, or the fraction of white dwarfs with hydrogen deficient atmospheres. We present a population synthesis study of the white dwarf population of NGC 6397. In particular, we study the shape of the color-magnitude diagram, and the corresponding magnitude and color distributions. We do this using an up-to-date Monte Carlo code that incorporates the most recent and reliable cooling sequences and an accurate modeling of the observational biases. We find a good agreement between our theoretical models and the observed data. In particular, we find that this agreement is best for those cooling sequences that take into account residual hydrogen burning. This result has important consequences for the evolution of progenitor stars during the thermally-pulsing asymptotic giant branch phase, since it implies that appreciable third dredge-up in low-mass, low-metallicity progenitors is not expected to occur. Using a standard burst duration of 1.0 Gyr, we obtain that the age of the cluster is 12.8+0.50-0.75 Gyr. Larger ages are also compatible with the observed data, but then realistic longer durations of the initial burst of star formation are needed to fit the luminosity function. We conclude that a correct modeling of the white dwarf opulation of globular clusters, used in combination with the number counts of main sequence stars provides an unique tool to model the properties of globular clusters.

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.

Modern large-scale surveys have allowed the identification of large numbers of white dwarfs. However, these surveys are subject to complicated target selection algorithms, which make it almost impossible to quantify to what extent the observational biases affect the observed populations. The LAMOST (Large Sky Area Multi-Object Fiber Spectroscopic Telescope) Spectroscopic Survey of the Galactic anti-center (LSS-GAC) follows a well-defined set of criteria for selecting targets for observations. This advantage over previous surveys has been fully exploited here to identify a small yet well-characterised magnitude-limited sample of hydrogen-rich (DA) white dwarfs. We derive preliminary LSS-GAC DA white dwarf luminosity and mass functions. The space density and average formation rate of DA white dwarfs we derive are 0.83+/-0.16 x 10^{-3} pc^{-3} and 5.42 +/- 0.08 x 10^{-13} pc^{-3} yr^{-1}, respectively. Additionally, using an existing Monte Carlo population synthesis code we simulate the population of single DA white dwarfs in the Galactic anti-center, under various assumptions. The synthetic populations are passed through the LSS-GAC selection criteria, taking into account all possible observational biases. This allows us to perform a meaningful comparison of the observed and simulated distributions. We find that the LSS-GAC set of criteria is highly efficient in selecting white dwarfs for spectroscopic observations (80-85 per cent) and that, overall, our simulations reproduce well the observed luminosity function. However, they fail at reproducing an excess of massive white dwarfs present in the observed mass function. A plausible explanation for this is that a sizable fraction of massive white dwarfs in the Galaxy are the product of white dwarf-white dwarf mergers.

47 Tucanae is one of the most interesting and well observed and theoretically studied globular clusters. This allows us to study the reliability of our understanding of white dwarf cooling sequences, to confront different methods to determine its age, and to assess other important characteristics, like its star formation history. Here we present a population synthesis study of the cooling sequence of the globular cluster 47 Tucanae. In particular, we study the distribution of effective temperatures, the shape of the color-magnitude diagram, and the corresponding magnitude and color distributions. We do so using an up-to-date population synthesis code based on Monte Carlo techniques, that incorporates the most recent and reliable cooling sequences and an accurate modeling of the observational biases. We find a good agreement between our theoretical models and the observed data. Thus, our study, rules out previous claims that there are still missing physics in the white dwarf cooling models at moderately high effective temperatures. We also derive the age of the cluster using the termination of the cooling sequence, obtaining a good agreement with the age determinations using the main-sequence turn-off. Finally, we find that the star formation history of the cluster is compatible with that btained using main sequence stars, which predict the existence of two distinct populations. We conclude that a correct modeling of the white dwarf population of globular clusters, used in combination with the number counts of main sequence stars provides an unique tool to model the properties of globular clusters.

This is the product of a Gaia GREAT workshop on "Gaia and the end states of stellar evolution" held at the University of Leicester in April 2011.

Detached WD+MS PCEBs are perhaps the most suitable objects for testing predictions of close-compact binary-star evolution theories, in particular, CE evolution. The population of WD+MS PCEBs has been simulated by several authors in the past and compared with observations. However, most of those predictions did not take the possible contributions to the envelope ejection from additional sources of energy (mostly recombination energy) into account. Here we update existing binary population models of WD+MS PCEBs by assuming that a fraction of the recombination energy available within the envelope contributes to ejecting the envelope. We performed Monte Carlo simulations of 10^7 MS+MS binaries for 9 different models using standard assumptions for the initial primary mass function, binary separations, and initial-mass-ratio distribution and evolved these systems using the publicly available BSE code. Including a fraction of recombination energy leads to a clear prediction of a large number of long orbital period (>~10 days) systems mostly containing high-mass WDs. The fraction of systems with He-core WD primaries increases with the CE efficiency and the existence of very low-mass He WDs is only predicted for high values of the CE efficiency (>~0.5). All models predict on average longer orbital periods for PCEBs containing C/O-core WDs than for PCEBs containing He WDs. This effect increases with increasing values of both efficiencies. Longer periods after the CE phase are also predicted for systems containing more massive secondary stars. The initial-mass-ratio distribution affects the distribution of orbital periods, especially the distribution of secondary star masses. Our simulations, in combination with a large and homogeneous observational sample, can provide constraints on the values of the CE efficiencies, as well as on the initial-mass-ratio distribution for MS+MS binary stars.

It has been known for a long time that stellar metallicity plays a significant role in the determination of the ages of the different Galactic stellar populations, when main sequence evolutionary tracks are employed. Here we analyze the role that metallicity plays on the white dwarf luminosity function of the Galactic disk, which is often used to determine its age. We employ a Monte Carlo population synthesis code that accounts for the properties of the population of Galactic disk white dwarfs. Our code incorporates the most up-to-date evolutionary cooling sequences for white dwarfs with hydrogen-rich and hydrogen-deficient atmospheres for both carbon-oxygen and oxygen-neon cores. We use two different models to assess the evolution of the metallicity, one in which the adopted metallicity is constant with time, but with a moderate dispersion, and a second one in which the metallicity increases with time. We found that our theoretical results are in a very satisfactory agreement with the observational luminosity functions obtained from the Sloan Digital Sky Survey (SDSS) and from the SuperCOSMOS Sky Survey (SSS), independently of the adopted age-metallicity law. In particular, we found that the age-metallicity law has no noticeable impact in shaping the bright branch of the white dwarf luminosity function, and that the position of its cut-off is almost insensitive to the adopoted age-metallicity relationship. Because the shape of the bright branch of the white dwarf luminosity function is insensitive to the age-metallicity law, it can be safely employed to test the theoretical evolutionary sequences, while due to the limited sensitivity of the position of the drop-off to the distribution of metallicities, its location provides a robust indicator of the age of the Galactic disk.

Detached white dwarf + main sequence (WD+MS) systems represent the simplest population of post-common envelope binaries (PCEBs). Since the ensemble properties of this population carries important information about the characteristics of the common-envelope (CE) phase, it deserves close scrutiny. However, most population synthesis studies do not fully take into account the effects of the observational selection biases of the samples used to compare with the theoretical simulations. Here we present the results of a set of detailed Monte Carlo simulations of the population of WD+MS binaries in the Sloan Digital Sky Survey (SDSS) Data Release 7. We used up-to-date stellar evolutionary models, a complete treatment of the Roche lobe overflow episode, and a full implementation of the orbital evolution of the binary systems. Moreover, in our treatment we took into account the selection criteria and all the known observational biases. Our population synthesis study allowed us to make a meaningful comparison with the available observational data. In particular, we examined the CE efficiency, the possible contribution of internal energy, and the initial mass ratio distribution (IMRD) of the binary systems. We found that our simulations correctly reproduce the properties of the observed distribution of WD+MS PCEBs. In particular, we found that once the observational biases are carefully taken into account, the distribution of orbital periods and of masses of the WD and MS stars can be correctly reproduced for several choices of the free parameters and different IMRDs, although models in which a moderate fraction (<=10%) of the internal energy is used to eject the CE and in which a low value of CE efficiency is used (<=0.3) seem to fit better the observational data. We also found that systems with He-core WDs are over-represented in the observed sample, due to selection effects.

We present a coherent and detailed Monte Carlo simulation of the population of hot white dwarfs. We assess the statistical significance of the hot end of the white dwarf luminosity function and the role played by the bolometric corrections of hydrogen-rich white dwarfs at high effective temperatures. We use the most up-to-date stellar evolutionary models and implement a full description of the observational selection biases to obtain realistic simulations of the observed white dwarf population. Our theoretical results are compared with the luminosity function of hot white dwarfs obtained from the Sloan Digital Sky Survey (SDSS), for both DA and non-DA white dwarfs. We find that the theoretical results are in excellent agreement with the observational data for the population of white dwarfs with hydrogen deficient atmospheres (non-DA white dwarfs). For the population of white dwarfs with hydrogen-rich atmospheres (white dwarfs of the DA class), our simulations show some discrepancies with the observations for the brightest luminosity bins. These discrepancies can be attributed to the way in which the masses of the white dwarfs contributing to this luminosity bin have been computed, as most of them have masses smaller than the theoretical lower limit for carbon-oxygen white dwarfs. We conclude that the way in which the observational luminosity function of hot white dwarfs is obtained is very sensitive to the particular implementation of the method used to derive the masses of the sample. We also provide a revised luminosity function for hot white dwarfs with hydrogen-rich atmospheres.

A secular variation of $G$ modifies the structure and evolutionary time scales of white dwarfs. Using an state-of-the-art stellar evolutionary code, an up-to-date pulsational code, and a detailed population synthesis code we demonstrate that the effects of a running $G$ are obvious both in the properties of individual white dwarfs, and in those of the white dwarf populations in clusters. Specifically, we show that the white dwarf evolutionary sequences depend on both the value of $\dot G/G$, and on the value of $G$ when the white dwarf was born. We show as well that the pulsational properties of variable white dwarfs can be used to constrain $\dot G/G$. Finally, we also show that the ensemble properties of of white dwarfs in clusters can also be used to set upper bounds to $\dot G/G$. Precisely, the tightest bound --- $\dot G/G \sim -1.8 10^{-12}$ yr$^{-1}$ --- is obtained studying the population of the old, metal-rich, well populated, open cluster NGC 6791. Less stringent upper limits can be obtained comparing the theoretical results obtained taking into account the effects of a running $G$ with the measured rates of change of the periods of two well studied pulsating white dwarfs, G117--B15A and R548. Using these white dwarfs we obtain $\dot G/G\sim -1.8\times 10^{-10}$ yr$^{-1}$, and $\dot G/G\sim -1.3\times 10^{-10}$ yr$^{-1}$, respectively, which although less restrictive than the previous bound, can be improved measuring the rate of change of the period of massive white dwarfs.

We present a partial-sky survey of the radio continuum at 2.3 GHz within the scope of the Galactic Emission Mapping (GEM) project, an observational program conceived and developed to reveal the large-scale properties of Galactic synchrotron radiation through a set of self-consistent surveys of the radio continuum between 408 MHz and 10 GHz. The GEM experiment uses a portable and double-shielded 5.5-m radiotelescope in altazimuthal configuration to map 60-degree-wide declination bands from different observational sites by circularly scanning the sky at zenithal angles of 30 deg from a constantly rotating platform. The observations were accomplished with a total power receiver, whose front-end High Electron Mobility Transistor (HEMT) amplifier was matched directly to a cylindrical horn at the prime focus of the parabolic reflector. The Moon was used to calibrate the antenna temperature scale and the preparation of the map required direct subtraction and destriping algorithms to remove ground contamination as the most significant source of systematic error. We used 484 hours of total intensity observations from two locations in Colombia and Brazil to yield 66% sky coverage from DEC = -51.73 deg to DEC = +34.78 deg. The zero-level uncertainty of the combined survey is 103 mK with a temperature scale error of 5% after direct correlation with the Rhodes/HartRAO survey at 2326 MHz on a T-T plot. The sky brightness distribution into regions of low and high emission in the GEM survey is consistent with the appearance of a transition region as seen in the Haslam 408 MHz and WMAP K-band surveys. Preliminary results also show that the temperature spectral index between 408 MHz and the 2.3 GHz band of the GEM survey has a weak spatial correlation with these regions; but it steepens significantly from high to low emission regions with respect to the WMAP K-band survey.

The evolution of white dwarfs is a simple gravothermal process. This process can be tested in two ways, through the luminosity function of these stars and through the secular variation of the period of pulsation of those stars that are variable. Here we show how the mass of the axion can be constrained using the white dwarf luminosity function.

We argue that the multiple shells of circumstellar material (CSM) and the supernovae (SN) ejecta interaction with the CSM starting 59 days after the explosion of the Type Ia SN (SN Ia) PTF 11kx, are best described by a violent prompt merger. In this prompt merger scenario the common envelope (CE) phase is terminated by a merger of a WD companion with the hot core of a massive asymptotic giant (AGB) star. In most cases the WD is disrupted and accreted onto the more massive core. However, in the rare cases where the merger takes place when the WD is denser than the core, the core will be disrupted and accreted onto the cooler WD. In such cases the explosion might occur with no appreciable delay, i.e., months to years after the termination of the CE phase. This, we propose, might be the evolutionary route that could lead to the explosion of PTF 11kx. This scenario can account for the very massive CSM within ~1000 AU of the exploding PTF 11kx star, for the presence of hydrogen, and for the presence of shells in the CSM.

We present a detailed Monte Carlo simulation of the population of the hot branch of the white dwarf luminosity function. We used the most up-to-date stellar evolutionary models and we implemented a full description of the observational selection biases. Our theoretical results are compared with the luminosity function of hot white dwarfs obtained from the Sloan Digital Sky Survey (SDSS), for both DA and non-DA white dwarfs. For non-DA white dwarfs we find an excellent agreement with the observational data, while for DA white dwarfs our simulations show some discrepancies with the observations for the brightest luminosity bins, those corresponding to L>= 10 L_sun.

High-field magnetic white dwarfs have been long suspected to be the result of stellar mergers. However, the nature of the coalescing stars and the precise mechanism that produces the magnetic field are still unknown. Here we show that the hot, convective, differentially rotating corona present in the outer layers of the remnant of the merger of two degenerate cores is able to produce magnetic fields of the required strength that do not decay for long timescales. We also show, using an state-of-the-art Monte Carlo simulator, that the expected number of high-field magnetic white dwarfs produced in this way is consistent with that found in the solar neighborhood.

White dwarfs are almost completely degenerate objects that cannot obtain energy from thermonuclear sources, so their evolution is just a gravothermal cooling process. Recent improvements in the accuracy and precision of the luminosity function and in pulsational data of variable white dwarfs suggest that they are cooling faster than expected from conventional theory. In this contribution we show that the inclusion of an additional cooling term due to axions able to interact with electrons with a coupling constant g_ae ~(2-7)x10^{-13} allows to fit better the observations.

High-field magnetic white dwarfs have been long suspected to be the result of stellar mergers. However, the nature of the coalescing stars and the precise mechanism that produces the magnetic field are still unknown. Here we show that the hot, convective, differentially rotating corona present in the outer layers of the remnant of the merger of two degenerate cores is able to produce magnetic fields of the required strength that do not decay for long timescales. We also show, using an state-of-the-art Monte Carlo simulator, that the expected number of high-field magnetic white dwarfs produced in this way is consistent with that found in the Solar neighborhood.

NGC 6791 is a well-studied, metal-rich open cluster that is so close to us that can be imaged down to luminosities fainter than that of the termination of its white dwarf cooling sequence, thus allowing for an in-depth study of its white dwarf population. We use a Monte Carlo simulator that employs up-to-date evolutionary cooling sequences for white dwarfs with hydrogen-rich and hydrogen-deficient atmospheres, with carbon-oxygen and helium cores. The cooling sequences for carbon-oxygen cores account for the delays introduced by both Ne^22 sedimentation in the liquid phase and by carbon-oxygen phase separation upon crystallization. We do not find evidence for a substantial fraction of helium-core white dwarfs, and hence our results support the suggestion that the origin of the bright peak of the white dwarf luminosity function can only be attributed to a population of unresolved binary white dwarfs. Moreover, our results indicate that the number distribution of secondary masses of the population of unresolved binaries has to increase with increasing mass ratio between the secondary and primary components of the progenitor system. We also find that the observed cooling sequence appears to be able to constrain the presence of progenitor sub-populations with different chemical compositions and the fraction of non-DA white dwarfs. Our simulations place interesting constraints on important characteristics of the stellar populations of NGC 6791. In particular, we find that the fraction of single helium-core white dwarfs must be smaller than 5%, that a sub-population of stars with zero metallicity must be <12%, while if the adopted metallicity of the sub-population is solar the upper limit is ~8%. Finally, we also find that the fraction of non-DA white dwarfs in this particular cluster is surprinsingly small <6%.

A variation of the gravitational constant over cosmological ages modifies the main sequence lifetimes and white dwarf cooling ages. Using an state-of-the-art stellar evolutionary code we compute the effects of a secularly varying G on the main sequence ages and, employing white dwarf cooling ages computed taking into account the effects of a running G, we place constraints on the rate of variation of Newton's constant. This is done using the white dwarf luminosity function and the distance of the well studied open Galactic cluster NGC 6791. We derive an upper bound G'/G ~ -1.8 10^{-12} 1/yr. This upper limit for the secular variation of the gravitational constant compares favorably with those obtained using other stellar evolutionary properties, and can be easily improved if deep images of the cluster allow to obtain an improved white dwarf luminosity function.

Within the theoretical framework of some modern unification theories the constants of nature are functions of cosmological time. White dwarfs offer the possibility of testing a possible variation of G and, thus, to place constraints to these theories. We present full white dwarf evolutionary calculations in the case that G decreases with time. White dwarf evolution is computed in a self-consistent way, including the most up-to-date physical inputs, non-gray model atmospheres and a detailed core chemical composition that results from the calculation of the full evolution of progenitor stars. We find that the mechanical structure and the energy balance of white dwarfs are strongly modified by the presence of a varying G. In particular, for certain values of the rate of change of G, the evolution of cool white dwarfs is markedly affected. The impact of a varying G is more notorious in the case of more massive white dwarfs. In view of the recent results reporting that a very accurate white dwarf cooling age can be derived for the old and metal-rich open cluster NGC 6791, our study suggests that this cluster could be a potential target to constrain or detect a ypothetical secular variation of G.

White dwarfs are almost completely degenerate objects that cannot obtain energy from the thermonuclear sources and their evolution is just a gravothermal process of cooling. The simplicity of these objects, the fact that the physical inputs necessary to understand them are well identified, although not always well understood, and the impressive observational background about white dwarfs make them the most well studied Galactic population. These characteristics allow to use them as laboratories to test new ideas of physics. In this contribution we discuss the robustness of the method and its application to the axion case.

NGC 6791 is a well studied open cluster1 that it is so close to us that can be imaged down to very faint luminosities. The main sequence turn-off age (~8 Gyr) and the age derived from the termination of the white dwarf cooling sequence (~6 Gyr) are significantly different. One possible explanation is that as white dwarfs cool, one of the ashes of helium burning, 22Ne, sinks in the deep interior of these stars. At lower temperatures, white dwarfs are expected to crystallise and phase separation of the main constituents of the core of a typical white dwarf, 12C and 16O, is expected to occur. This sequence of events is expected to introduce significant delays in the cooling times, but has not hitherto been proven. Here we report that, as theoretically anticipated, physical separation processes occur in the cores of white dwarfs, solving the age discrepancy for NGC 6791.

The nature of the several microlensing events observed by the MACHO team towards the Large Magellanic Cloud (LMC) is still a subject of debate. Low-mass substellar objects and stars with masses larger than ~M_{sun} have been ruled out as major components of a Massive Astrophysical Compact Halo Object (MACHO) Galactic halo, while stars of half a solar mass seem to be viable candidates. Main sequence stars have been already discarded, and there are tight restrictions on the role played by white dwarfs with hydrogen-dominated atmospheres. In this paper we evaluate the contribution to the dark matter content of the Galaxy of white dwarfs with hydrogen-deficient atmospheres. For this purpose we use a Monte Carlo simulator which incorporates up-to-date evolutionary sequences of white dwarfs with hydrogen-rich and hydrogen-deficient atmospheres. We also take into account detailed descriptions of the thick disk and the halo of our Galaxy as well as of a reliable model of the LMC. We find that the contribution of white dwarfs with hydrogen-deficient atmospheres moderately increases the theoretical estimate of the optical depth with respect to the value obtained when only hydrogen-rich white dwarfs are considered. We also find that the contribuiton of the thick disk population of white dwarfs is comparable to the halo contribution. However, the contributions of both the halo and the thick disk white-dwarf populations are still insufficient to explain the number of events observed by the MACHO team. Finally, we find that the contribution to the halo dark matter of the entire population under study is less than 10% at the 95% conficence level.

We present final results on the angular power spectrum of total intensity anisotropies in the CMB from the CBI. Our analysis includes all primordial anisotropy data collected between January 2000 and April 2005, and benefits significantly from an improved maximum likelihood analysis pipeline. It also includes results from a 30 GHz foreground survey conducted with the Green Bank Telescope (GBT) which places significant constraints on the possible contamination due to foreground point sources. We improve on previous CBI results by about a factor of two in the damping tail. These data confirm, at ~3-sigma, the existence of an excess of power over intrinsic CMB anisotropy on small angular scales (l > 1800). Using the GBT survey, we find currently known radio source populations are not capable of generating the power; a new population of faint sources with steeply rising spectral indices would be required to explain the excess with sources... We also present a full cosmological parameter analysis of the new CBI power spectrum... With CBI alone, the full parameter analysis finds the excess is 1.6-sigma above the level expected for a sigma_8=0.8 universe. We find the addition of high-l CMB data substantially improves constraints on cosmic string contributions to the TT power spectrum as well as the running of the scalar spectral index... We also present forecasts for what other experiments should see at different frequencies and angular resolutions given the excess power observed by CBI. We find that the reported high-l bandpowers from current high resolution CMB bolometer experiments are consistent with each other and CBI if the excess power is due to the SZE at the CBI-level of 2.5 +/- 1 times the sigma_8=0.8 standard SZ template. <Abridged>

White dwarfs are the end-product of the lifes of intermediate- and low-mass stars and their evolution is described as a simple cooling process. Recently, it has been possible to determine with an unprecedented precision their luminosity function, that is, the number of stars per unit volume and luminosity interval. We show here that the shape of the bright branch of this function is only sensitive to the averaged cooling rate of white dwarfs and we propose to use this property to check the possible existence of axions, a proposed but not yet detected weakly interacting particle. Our results indicate that the inclusion of the emission of axions in the evolutionary models of white dwarfs noticeably improves the agreement between the theoretical calculations and the observational white dwarf luminosity function. The best fit is obtained for m_a cos^2 \beta ~ 5 meV, where m_a is the mass of the axion and cos^2 \beta is a free parameter. We also show that values larger than 10 meV are clearly excluded. The existing theoretical and observational uncertainties do not allow yet to confirm the existence of axions, but our results clearly show that if their mass is of the order of few meV, the white dwarf luminosity function is sensitive enough to detect their existence.

The evolution of white dwarfs can be described as a simple cooling process. Recently, it has been possible to determine with an unprecedented precision their luminosity function, that is, the number of stars per unit volume and luminosity interval. Since the shape of the bright branch of this function is only sensitive to the average cooling rate, we use this property to check the possible existence of axions, a proposed but not yet detected weakly interacting particle. We show here that the inclusion of the axion emissivity in the evolutionary models of white dwarfs noticeably improves the agreement between the theoretical calculations and the observational white dwarf luminosity function, thus providing the first positive indication that axions could exist. Our results indicate that the best fit is obtained for m_a cos^2beta ~ 2-6 meV, where m_a is the mass of the axion and cos^2beta is a free parameter, and that values larger than 10 meV are clearly excluded.

We explore different evolutionary scenarios to explain the helium deficiency observed in H1504+65, the most massive known PG1159 star. We concentrate mainly on the possibility that this star could be the result of mass loss shortly after the born-again and during the subsequent evolution through the [WCL] stage. This possibility is sustained by recent observational evidence of extensive mass-loss events in Sakurai's object and is in line with the recent finding that such mass losses give rise to PG1159 models with thin helium-rich envelopes and large rates of period change, as demanded by the pulsating star PG1159-035. We compute the post born again evolution of massive sequences by taking into account different mass-loss rate histories. Our results show that stationary winds during the post-born-again evolution fail to remove completely the helium-rich envelope so as to explain the helium deficiency observed in H1504+65. Stationary winds during the Sakurai and [WCL] stages only remove at most half of the envelope surviving the violent hydrogen burning during the born-again phase. In view of our results, the recently suggested evolutionary connection born-again stars --> H1504+65 --> white dwarfs with carbon-rich atmospheres is difficult to sustain unless the whole helium-rich envelope could be ejected by non-stationary mass-loss episodes during the Sakurai stage.

The nature of the several microlensing events observed by the MACHO team towards the LMC still remains controversial. Low-mass substellar objects and stars with masses larger than ~1 M_{sun} have been ruled out as major components of a MACHO galactic halo, while stars of half solar masses are the most probable candidates. In this paper we assess jointly the relative contributions of both red dwarfs and white dwarfs to the mass budget of the galactic halo. In doing so we use a Monte Carlo simulator which incorporates up-to-date evolutionary sequences of both red dwarfs and white dwarfs as well as detailed descriptions of our Galaxy and of the LMC. We explore the complete mass range between 0.08 and 1 M_{sun} as possible microlensing candidates and we compare the synthetic populations obtained with our simulator with the results obtained by the MACHO and EROS experiments. Our results indicate that, despite that the contribution of the red dwarf population increases by a factor of 2 the value of the optical depth obtained when taking into account the white dwarf population alone, it is still insufficient to explain the number of events claimed by the MACHO team. Finally, we find that the contribution to the halo dark matter of the whole population under study is smaller than 10% at the 95% conficence level.

The interpretation of microlensing results towards the Large Magellanic Cloud (LMC) still remains controversial. White dwarfs have been proposed to explain these results and, hence, to contribute significantly to the mass budget of our Galaxy. However, several constraints on the role played by regular carbon-oxygen white dwarfs exist. Massivewhite dwarfs are thought to be made of a mixture of oxygen and neon. Correspondingly, their cooling rate is larger than those of typical carbon-oxygen white dwarfs and they fade to invisibility in short timescales. Consequently, they constitute a good candidate for explaining the microlensing results. Here, we examine in detail this hypothesis by using the most recent and up-to-date cooling tracks for massive white dwarfs and a Monte Carlo simulator which takes into account the most relevant Galactic inputs. We find that oxygen-neon white dwarfs cannot account for a substantial fraction of the microlensing depth towards the LMC, independently of the adopted initial mass function, although some microlensing events could be due to oxygen--neon white dwarfs. The white dwarf population contributes at most a 5% to the mass of the Galactic halo.

Only the 1/V_max method has been employed so far for observationally determining the white dwarf luminosity function, whereas for other kind of luminosity functions several other methods have been frequently used. Moreover, the procedures to determine the white dwarf luminosity function are not free of biases. These biases have two different origins: they can either be of statistical nature or a consequence of the measurement errors. In a previous paper we carried out an in-depth study of the first category of biases for several luminosity function estimators. In this paper we focus on the biases introduced by the measurement errors and on the effects of the degree of contamination of the input sample used to build the disc white dwarf luminosity function by different kinematical populations. To assess the extent of these biases we use a Monte Carlo simulator to generate a controlled synthetic population and analyse the behaviour of the disc white dwarf luminosity function for several assumptions about the magnitude of the measurement errors and for several degrees of contamination, comparing the performances of the most robust luminosity function estimators under such conditions.

We present new measurements of the power spectra of the E-mode of CMB polarization, the temperature T, the cross-correlation of E and T, and upper limits on the B-mode from 2.5 years of dedicated Cosmic Background Imager (CBI) observations. Both raw maps and optimal signal images in the uv-plane and real space show strong detections of the E-mode (11.7 sigma for the EE power spectrum overall) and no detection of the B-mode. The power spectra are used to constrain parameters of the flat tilted adiabatic Lambda-CDM models: those determined from EE and TE bandpowers agree with those from TT, a powerful consistency check. There is little tolerance for shifting polarization peaks from the TT-forecast locations, as measured by the angular sound crossing scale theta = 100 ell_s = 1.03 +/- 0.02 from EE and TE cf. 1.044 +/- 0.005 with the TT data included. The scope for extra out-of-phase peaks from subdominant isocurvature modes is also curtailed. The EE and TE measurements of CBI, DASI and BOOMERANG are mutually consistent, and, taken together rather than singly, give enhanced leverage for these tests.

Over the years, several methods have been proposed to compute galaxy luminosity functions, from the most simple ones -counting sample objects inside a given volume- to very sophisticated ones -like the C- method, the STY method or the Choloniewski method, among others. However, only the V/Vmax method is usually employed in computing the white dwarf luminosity function and other methods have not been applied so far to the observational sample of spectroscopically identified white dwarfs. Moreover, the statistical significance of the white dwarf luminosity function has also received little attention and a thorough study still remains to be done. In this paper we study, using a controlled synthetic sample of white dwarfs generated using a Monte Carlo simulator, which is the statistical significance of the white dwarf luminosity function and which are the expected biases. We also present a comparison between different estimators for computing the white dwarf luminosity function. We find that for sample sizes large enough the V/Vmax method provides a reliable characterization of the white dwarf luminosity function, provided that the input sample is selected carefully. Particularly, the V/Vmax method recovers well the position of the cut-off of the white dwarf luminosity function. However, this method turns out to be less robust than the Choloniewski method when the possible incompletenesses of the sample are taken into account. We also find that the Choloniewski method performs better than the V/Vmax method in estimating the overall density of white dwarfs, but misses the exact location of the cut-off of the white dwarf luminosity function.

One of the most promising space missions of ESA is the astrometric satellite Gaia, which will provide very precise astrometry and multicolour photometry, for all 1.3 billion objects to V~20, and radial velocities with accuracies of a few km/s for most stars brighter than V ~ 17. Consequently, full homogeneous six-dimensional phase-space information for a huge number of stars will become available. Our Monte Carlo simulator has been used to estimate the number of white dwarfs potentially observable by Gaia. From this we assess which would be the white dwarf luminosity functions which Gaia will obtain and discuss in depth the scientific returns of Gaia in the specific field of white dwarf populations. Scientific attainable goals include, among others, a reliable determination of the age of the Galactic disk, a better knowledge of the halo of the Milky Way and the reconstruction of the star formation history of the Galactic disk. Our results also demonstrate the potential impact of a mission like Gaia in the current understanding of the white dwarf cooling theory.

Polarization observations of the cosmic microwave background with the Cosmic Background Imager from September 2002 to May 2004 provide a significant detection of the E-mode polarization and reveal an angular power spectrum of polarized emission showing peaks and valleys that are shifted in phase by half a cycle relative to those of the total intensity spectrum. This key agreement between the phase of the observed polarization spectrum and that predicted based on the total intensity spectrum provides support for the standard model of cosmology, in which dark matter and dark energy are the dominant constituents, the geometry is close to flat, and primordial density fluctuations are predominantly adiabatic with a matter power spectrum commensurate with inflationary cosmological models.

Two years of microwave background observations with the Cosmic Background Imager (CBI) have been combined to give a sensitive, high resolution angular power spectrum over the range 400 < l < 3500. This power spectrum has been referenced to a more accurate overall calibration derived from WMAP. The data cover 90 deg^2 including three pointings targeted for deep observations. The uncertainty on the l > 2000 power previously seen with the CBI is reduced. Under the assumption that any signal in excess of the primary anisotropy is due to a secondary Sunyaev-Zeldovich anisotropy in distant galaxy clusters we use CBI, ACBAR, and BIMA data to place a constraint on the present-day rms mass fluctuation sigma_8. We present the results of a cosmological parameter analysis on the l < 2000 primary anisotropy data which show significant improvements in the parameters as compared to WMAP alone, and we explore the role of the small-scale cosmic microwave background data in breaking parameter degeneracies.

The interpretation of microlensing results towards the Large Magellanic Cloud (LMC) still remains controversial. Whereas white dwarfs have been proposed to explain these results and, hence, to contribute significantly to the mass budget of our Galaxy, there are as well several constraints on the role played by white dwarfs. In this paper we analyze self-consistently and simultaneously four different results, namely, the local halo white dwarf luminosity function, the microlensing results reported by the MACHO team towards the LMC, the results of Hubble Deep Field (HDF) and the results of the EROS experiment, for several initial mass functions and halo ages. We find that the proposed log-normal initial mass functions do not contribute to solve the problem posed by the observed microlensing events and, moreover, they overproduce white dwarfs when compared to the results of the HDF and of the EROS survey. We also find that the contribution of hydrogen-rich white dwarfs to the dynamical mass of the halo of the Galaxy cannot be more than $\sim 4%$.

The interpretation of the old, cool white dwarfs recently found by Oppenheimer et al. (2001) is still controversial. Whereas these authors claim that they have finally found the elusive ancient halo white dwarf population that contributes significantly to the mass budget of the galactic halo, there have been several other contributions that argue that these white dwarfs are not genuine halo members but, instead, thick disk stars. We show here that the interpretation of this sample is based on the adopted distances, which are obtained from a color--magnitude calibration, and we demonstrate that when the correct distances are used a sizeable fraction of these putative halo white dwarfs belong indeed to the disk population. We also perform a maximum likelihood analysis of the remaining set of white dwarfs and we find that they most likely belong to the thick disk population. However, another possible explanation is that this sample of white dwarfs has been drawn from a 1:1 mixture of the halo and disk white dwarf populations.

We analyze in detail various possible sources of systematic errors on the distances of globular clusters derived by fitting a local template DA white dwarf sequence to the cluster counterpart (the so-called WD-fitting technique). We find that the unknown thickness of the hydrogen layer of white dwarfs in clusters plays a non negligible role. For reasonable assumptions - supported by the few sparse available observational constraints - about the unknown mass and thickness of the hydrogen layer for the cluster white dwarfs, a realistic estimate of the systematic error on the distance is within +-0.10 mag. However, particular combinations of white dwarf masses and envelope thicknesses - which at present cannot be excluded a priori - could produce larger errors. Contamination of the cluster DA sequence by non-DA white dwarfs introduces a very small systematic error of about -0.03 mag in the Mv/(V-I) plane, but in the Mv/(B-V) plane the systematic error amounts to ~ +0.20 mag. Contamination by white dwarfs with helium cores should not influence appreciably the WD-fitting distances. Finally, we obtain a derivative D((m-M)v)/D(E(B-V))~ -5.5 for the WD-fitting distances, which is very similar to the dependence found when using the Main Sequence fitting technique.

We have mapped the total sky brightness at 1465 MHz in two adjacent 60-degree declination bands with the portable 5.5-m parabolic reflector of the Galactic Emission Mapping (GEM) project, an on-going international collaboration to survey the radio continuum of the sky in decimeter and centimeter wavelengths. The observations were conducted from two locations, one in the USA and the other in Brazil, using a novel instrumental approach to overcome the well-known shortcomings of survey experiments. Our strategy consists of a 1-rpm rotating dish to circularly scan the sky at 30 degrees from zenith. The dish uses a rim-halo to re-direct the spillover sidelobes of its backfire helical feed toward the sky and the entire assembly has been enclosed inside a wire mesh ground shield in order to minimize and level out the contamination from the ground. The diffraction characteristics of this set-up have been succesfully modelled and undesired systematic striping across the observed bands has been carefully removed by a baseline propagation method which exploits the time-forward and time-backward intersections of the circular scans. The map displays nearly 300 hours of our best quality data taken with a HPBW of 5.4 degrees at a sensitivity of 20 mK.

We have analyzed observations from a radioastronomical experiment to survey the sky at decimetric wavelengths along with feed pattern measurements in order to account for the level of ground contamination entering the sidelobes. A major asset of the experiment is the use of a wire mesh fence around the rim-halo shielded antenna with the purpose of levelling out and reducing this source of stray radiation for zenith-centered 1-rpm circular scans. We investigate the shielding performance of the experiment by means of a geometric diffraction model in order to predict the level of the spillover and diffraction sidelobes in the direction of the ground. Using 408 MHz and 1465 MHz feed measurements, the model shows how a weakly-diffracting and unshielded antenna configuration becomes strongly-diffracting and double-shielded as far-field diffraction effects give way to near-field ones. Due to the asymmetric response of the feeds, the orientation of their radiation fields with respect to the secondary must be known a priori before comparing model predictions with observational data. By adjusting the attenuation coefficient of the wire mesh the model is able to reproduce the amount of differential ground pick-up observed during test measurements at 1465 MHz.

The white dwarf luminosity function has proven to be an excellent tool to study some properties of the galactic disk such as its age and the past history of the local star formation rate. The existence of an observational luminosity function for halo white dwarfs could provide valuable information about its age, the time that the star formation rate lasted, and could also constrain the shape of the allowed Initial Mass Functions (IMF). However, the main problem is the scarce number of white dwarfs already identified as halo stars. In this Letter we show how an artificial intelligence algorithm can be succesfully used to classify the population of spectroscopically identified white dwarfs allowing us to identify several potential halo white dwarfs and to improve the significance of its luminosity function.

In order to understand the dynamical and chemical evolution of our Galaxy it is of fundamental importance to study the local neighborhood. White dwarf stars are ideal candidates to probe the history of the solar neighborhood, since these ``fossil'' stars have very long evolutionary time-scales and, at the same time, their evolution is relatively well understood. In fact, the white dwarf luminosity function has been used for this purpose by several authors. However, a long standing problem arises from the relatively poor statistics of the samples, especially at low luminosities. In this paper we assess the statistical reliability of the white dwarf luminosity function by using a Monte Carlo approach.

Halo white dwarfs can provide important information about the properties and evolution of the galactic halo. In this paper we compute, assuming a standard IMF and updated models of white dwarf cooling, the expected luminosity function, both in luminosity and in visual magnitude, for different star formation rates. We show that a deep enough survey (limiting magnitude > 20) could provide important information about the halo age and the duration of the formation stage. We also show that the number of white dwarfs produced using the recently proposed biased IMFs cannot represent a large fraction of the halo dark matter if they are constrained by the presently observed luminosity function. Furthermore, we show that a robust determination of the bright portion of the luminosity function can provide strong constraints on the allowable IMF shapes.

The geodesics followed by cosmic microwave background (CMB) photons show different behaviours depending on the geometry of space. Namely, the effect of `mixing geodesics' predicts a distinct signature in CMB maps: threshold-independent elongated anisotropy spots in negatively curved geometries. We have found statistically significant sign for spot elongation in the COBE four year maps. This can be a direct indication for the negative curvature of the Universe.

When considering the statistical properties of a bundle of cosmic microwave background (CMB) photons propagating through space, the effect of `mixing geodesics' appears with a distinct signature that depends on the geometry of space. In a Universe with negative curvature this effect is expected to produce elongated anisotropy spots on CMB maps. We used COBE-DMR data to look for such effect. Based on the analysis of a measure of eccentricity of hot spots it appears that there is a clear indication of an excess eccentricity of hot spots with respect to that expected from noise alone. This result must be interpreted with caution as this effect can be due in part to galactic emission. If the detected eccentricity of anisotropy spots can be attributed to the effect of mixing it implies the negative curvature of the Universe and a value of $\Omega < 1$.

The GEM (Galactic Emission Mapping) project is an international collaboration established with the aim of surveying the full sky at long wavelengths with a multi-frequency radio telescope. A total of 745 hours of observation at 408 MHz were completed from an Equatorial site in Colombia. The observations cover the celestial band $0^h < \alpha < 24^h$, and $-24^{\circ} \ 22^{\prime} < \delta < +35^{\circ} \ 37^{\prime}$. Preliminary results of this partial survey will be discussed. A review of the instrumental setup and a $\sim 10^{\circ}$ resolution sky map at 408 MHz is presented.

A precise determination of the primordial spectrum of matter density fluctuations at super-horizon scales is essential in understanding large scale structure in the universe. Attempts to constrain or obtain the primordial spectrum using data on cosmic microwave background (CMB) anisotropies has relied on statistical and correlation analyses that assume a power-law spectrum. We propose a method to derive $P(k)$ directly from the CMB angular power spectrum which does not presuppose the need to know anything about its functional form. The method consists of a direct inversion of the Sachs-Wolfe formula. Using this new analysis technique and COBE data we obtain an empirical $P(k)$ which 1) supports a power-law parameterization and 2) has an amplitude and spectral index consistent with previous analyses of the same data. We obtained for the spectral index, $n=1.52 \pm 0.4$ when the 2nd year COBE data is used and $n=1.22 \pm 0.3$ using the 4 year data set.

We perform the stastistics of temperature maxima and minima in COBE-DMR 2-year maps. For power-law spectra the surface distribution of peaks implies an amplitude consistent with more conventional analyses of COBE data (for instance, we get $Q_{{\rm rms-PS}}=17\pm 3\ \mu $K for a spectral index $% n=1),$ but not with the measured quadrupole $Q_{{\rm rms}}=6\pm 3\ \mu $K. This provides further support for the existence an infrared cutoff in the cosmic spectrum.

We have made a likelihood statistical analysis of the angular correlations in the {\it COBE}-DMR two-year sky maps by Monte Carlo simulation of the temperature fluctuations. We assume an open universe and consider as primordial power spectrum the Harrison-Zeldovich one, $P(k)=Ak$. We find that the flatness of the universe is not implied by the data. The quadrupole normalization amplitude, $Q_{rms-PS}$, is related to the density parameter, $\Omega$, by $Q_{rms-PS} = 10.67 + 55.81 \Omega - 128.59 \Omega^2 + 81.26 \Omega^3\ \mu$K. We have determined the p.d.f. of $\Omega$ due to cosmic plus sampling (i.e. $20^\circ$ galactic cut) variance which generically shows a bimodal shape. The uncertainty as given by the r.m.s. is $\approx 0.35$, therefore to better constrain $\Omega$ experiments sensitive to higher multipoles ($l>20$) should be considered.

A statistical analysis of texture on the {\it COBE}-DMR first year sky maps based on the genus and spot number is presented. A generalized $\chi^2$ statistic is defined in terms of ``observable'' quantities: the genus and spot density that would be measured by different cosmic observers. This strategy together with the use of Monte Carlo simulations of the temperature fluctuations, including all the relevant experimental parameters, represent the main difference with previous analyses. Based on the genus analysis we find a strong anticorrelation between the quadrupole amplitude $Q_{rms-PS}$ and the spectral index $n$ of the density fluctuation power spectrum at recombination of the form $Q_{rms-PS}= 22.2 \pm 1.7 - (4.7 \pm 1.3) \times n\ \mu$K for fixed $n$, consistent with previous works. The result obtained based on the spot density is consistent with this $Q_{rms-PS} (n)$ relation. In addition to the previous results we have determined, using Monte Carlo simulations, the minimum uncertainty due to cosmic variance for the determination of the spectral index with the genus analysis. This uncertainty is $\delta n\approx 0.2$.

After the detection of structure in the microwave sky, the characterization of the observed features may be a useful guide to current and future experiments. The first year COBE-DMR sky maps were analyzed in order to identify the most significant hot and cold spots. The statistical significance of potential spots of cosmic origin was evaluated with Monte Carlo simulations. The area, eccentricity, and location of the most significant spots are given.

Geometric characteristics of random fields are exploited to test the consistency of density perturbation model spectra with COBE data. These CMB maps are analyzed using the number of anisotropy hot spots and their boundary curvature. CMB maps which account for instrumental effects and sky coverage are Monte Carlo generated. These simulations show that a scale invariant Harrison-Zeldovich primordial Gaussian density fluctuation spectrum is consistent with the data. The CMB fluctuation coherence angle, based on boundary curvature, gives a spectral index n = 1.2 +- 0.3.

A phenomenological power spectrum of primordial density perturbations has been constructed by using both COBE data to probe the large wavelength region, and a double power law, locally deduced from galaxy catalogs, which describes the matter correlation function up to tens of Megaparsec. The shape of the spectrum P(k) of density fluctuations exhibits a peak that singles out a characteristic wavelength proportional to the cutoff radius R0 in the matter distribution (comparable to the distance at which matter becomes anticorrelated). {}From a least squares fit to COBE's angular correlation function we obtain R0 = 35 +- 12 Mpc/h for the correlation length, and n = 0.76 +- 0.3 for the spectral index of P(k) in the large wavelength region. The inferred scale in the spectrum is 51 +- 18 Mpc/h. This number agrees with that derived from the analysis of the correlation function of matter and with a preferred scale identified in IRAS PSC.