We assess the partition function and ionization degree of magnetized hydrogen atoms at thermodynamic equilibrium for a wide range of field intensities, $B\approx 10^5$-$10^{12}$~G. Evaluations include fitting formulae for an arbitrary number of binding energies, the coupling between the internal atomic structure and the center-of-mass motion across the magnetic field, and the formation of the so-called decentered states (bound states with the electron shifted from the Coulomb well). Non-ideal gas effects are treated within the occupational probability method. We also present general mathematical expressions for the bound state correspondence between the limits of zero-field and high-field. This let us evaluate the atomic partition function in a continuous way from the Zeeman perturbative regime to very strong fields. Results are shown for conditions found in atmospheres of magnetic white dwarf stars (MWDs), with temperatures $T\approx 5000$-$80000$~K and densities $\rho\approx 10^{-12}$-$10^{-3}$~g~cm$^3$. Our evaluations show a marked reduction of the gas ionization due to the magnetic field in the atmospheres of strong MWDs. We also found that decentered states could be present in the atmospheres of currently known hot MWDs, giving a significant contribution to the partition function in the strongest magnetized atmospheres.

Binary systems composed of a recycled millisecond pulsar and a stellar companion in close orbit could be excellent sites to diagnose pulsar winds. In such systems, the pulsar outflow irradiates and heats up the companion atmosphere, which can lead to the observation of strong day/night modulations in temperature. We demonstrate with particle shower simulations that the particle energy of the wind affects the heating depth in the atmosphere: the wind heat can be deposited above or below the photosphere, leading to different signatures in the observed spectra. We apply our method to four specific systems: We find that systems with cool night side companions showing strong temperature variations can give interesting lower limits on the particle energy in the winds. In particular, if the companion night side of PSR B1957+20 were to be suddenly irradiated, deep heating would only take place if particles with energy > 100 TeV were present. Observational evidence of deep heating in this system thus suggests that i) such particles exist in the pulsar wind and/or ii) binary evolution non-trivially takes the companion to the observed temperature asymmetry. Besides, the observed temperature difference can be maintained only with particle energies of the order of 100 MeV.

Ultra-massive white dwarfs are powerful tools to study various physical processes in the Asymptotic Giant Branch (AGB), type Ia supernova explosions and the theory of crystallization through white dwarf asteroseismology. Despite the interest in these white dwarfs, there are few evolutionary studies in the literature devoted to them. Here, we present new ultra-massive white dwarf evolutionary sequences that constitute an improvement over previous ones. In these new sequences, we take into account for the first time the process of phase separation expected during the crystallization stage of these white dwarfs, by relying on the most up-to-date phase diagram of dense oxygen/neon mixtures. Realistic chemical profiles resulting from the full computation of progenitor evolution during the semidegenerate carbon burning along the super-AGB phase are also considered in our sequences. Outer boundary conditions for our evolving models are provided by detailed non-gray white dwarf model atmospheres for hydrogen and helium composition. We assessed the impact of all these improvements on the evolutionary properties of ultra-massive white dwarfs, providing up-dated evolutionary sequences for these stars. We conclude that crystallization is expected to affect the majority of the massive white dwarfs observed with effective temperatures below $40\,000\, \rm K$. Moreover, the calculation of the phase separation process induced by crystallization is necessary to accurately determine the cooling age and the mass-radius relation of massive white dwarfs. We also provide colors in the GAIA photometric bands for our H-rich white dwarf evolutionary sequences on the basis of new models atmospheres. Finally, these new white dwarf sequences provide a new theoretical frame to perform asteroseismological studies on the recently detected ultra-massive pulsating white dwarfs.

A selection of 17 stars in the Sloan Digital Sky Survey, previously identified as DC class white dwarfs (WDs), has been reported to show spectra very close to blackbody radiation in the wavelength range from ultraviolet to infrared. Due to the absence of lines and other details in their spectra, the surface gravity of these objects has not been previously well constrained and their effective temperatures have been determined by fits to the continuum spectrum using pure helium atmosphere models. We compute model atmospheres with pure helium and H/He mixtures and use Gaia DR2 parallaxes available for 16 out of the 17 selected stars to analyze their physical properties. We find that the atmospheres of the selected stars are very probably contaminated with a trace amount of hydrogen as $-6 \leq \log{N_H/N_{He}} \leq -5.4$. For the 16 stars with Gaia parallaxes, we calculate a mean stellar mass $0.606\pm 0.076\,M_\odot$, which represents typical mass values and surface gravities ($7.8<\log g<8.3$) for WDs.

Iglesias et al. (2002) showed that the Rayleigh scattering from helium atoms decreases by collective effects in the atmospheres of cool white dwarf stars. Their study is here extended to consider an accurate evaluation of the atomic polarizability and the density effects involved in the Rayleigh cross section over a wide density-temperature region. The dynamic dipole polarizability of helium atoms in the ground state is determinated with the oscillator-strength distribution approach. The spectral density of oscillator strength considered includes most significant single and doubly excited transitions to discrete and continuum energies. Static and dynamic polarizability results are confronted with experiments and other theoretical evaluations shown a very good agreement. In addition, the refractive index of helium is evaluated with the Lorentz-Lorenz equation and shows a satisfactory agreement with the most recent experiments. The effect of spatial correlation of atoms on the Rayleigh scattering is calculated with Monte Carlo simulations and effective energy potentials that represent the particle interactions, covering fluid densities between 0.005 and a few g/cm$^3$ and temperatures between $1000$ K and $15000$ K. We provide analytical fits from which the Rayleigh cross section of fluid helium can be easily calculated at wavelength $\lambda>505.35$ \AA. Collision-induced light scattering was estimated to be the dominant scattering process at densities greater than 1-2 g/cm$^3$ depending on the temperature.

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.

The development of surveys which will be able to cover a large region of the sky several times per year will allow the massive detection of transient events taking place in timescales of years. In addition, the projected full digitalization of the Harvard plate collection will open a new window to identify slow transients taking place in timescales of centuries. In particular, these projects will allow the detection of stars undergoing slow eruptions as those expected during late helium flashes in the post-AGB evolution. In order to identify those transients which correspond with late helium flashes the development of synthetic light curves of those events is mandatory. In this connection we present preliminary results of a project aimed at computing grids of theoretical light curves of born again stars.

White dwarf evolution is essentially a gravothermal cooling process, which,for cool white dwarfs, sensitively depends on the treatment of the outer boundary conditions. We provide detailed outer boundary conditions appropriate for computing the evolution of cool white dwarfs employing detailed non-gray model atmospheres for pure H composition. We also explore the impact on the white dwarf cooling times of different assumptions for energy transfer in the atmosphere of cool white dwarfs. Detailed non-gray model atmospheres are computed taken into account non-ideal effects in the gas equation of state and chemical equilibrium, collision-induced absorption from molecules, and the Lyman alpha quasi-molecular opacity. Our results show that the use of detailed outer boundary conditions becomes relevant for effective temperatures lower than 5800 and 6100K for sequences with 0.60 and 0.90 M_sun, respectively. Detailed model atmospheres predict ages that are up to approx 10% shorter at log L/L_sun=-4 when compared with the ages derived using Eddington-like approximations at tau_Ross=2/3. We also analyze the effects of various assumptions and physical processes of relevance in the calculation of outer boundary conditions. In particular, we find that the Ly_alpha red wing absorption does not affect substantially the evolution of white dwarfs. White dwarf cooling timescales are sensitive to the surface boundary conditions for T_eff < 6000K. Interestingly enough, non-gray effects have little consequences on these cooling times at observable luminosities. In fact, collision-induced absorption processes, which significantly affect the spectra and colors of old white dwarfs with hydrogen-rich atmospheres, have not noticeable effects in their cooling rates, except throughout the Rosseland mean opacity.

We introduce a conditional pair distribution function (CPDF) which characterizes the probability density of finding an object (e.g., a particle in a fluid) to certain distance of other, with each of these two having a nearest neighbor to a fixed but otherwise arbitrary distance. This function describes special four-body configurations, but also contains contributions due to the so-called mutual nearest neighbor (two-body) and shared neighbor (three-body) configurations. The CPDF is introduced to improve a Helmholtz free energy method based on space partitions. We derive exact expressions of the CPDF and various associated quantities for randomly distributed, non-interacting points at Euclidean spaces of one, two and three dimensions. Results may be of interest in many diverse scientific fields, from fluid physics to social and biological sciences.

We argue that NSV 11749, an eruption observed in the early twentieth century, was a rare event known as "very late thermal pulse" (VLTP). To support our argument we compare the lightcurve of NSV 11749 with those of the two bonafide VLTP objects known to date, V4334 Sgr and V605 Aql, and with those predicted by state of the art stellar evolution models. Next, we explore the IPHAS and 2MASS catalogues for possible counterparts of the eruption. Our analysis shows that the VLTP scenario outperforms all other proposed scenarios as an explanation of NSV 11749. We identify an IPHAS/2MASS source at the eruption location of NSV 11749. The derived colors suggest that the object is not enshrouded in a thick dust shell as V605 Aql and V4334 Sgr. Also the absence of an apparent planetary nebula (PN) at the eruption location suggests differences with known VLTP objects which might be linked to the intensity of the eruption and the mass of the object. Further exploration of this source and scenario seems desirable. If NSV 11749 was a born again star, it would be the third event of its kind to have been observed and will strongly help us to increase our understanding on the later stages of stellar evolution and violent reactive convective burning.

Kowalski & Saumon (2006) identified the missing absorption mechanism in the observed spectra of cool white dwarf stars as the Ly-alpha red wing formed by the collisions between atomic and molecular hydrogen and successfully explained entire spectra of many cool DA-type white dwarfs. Owing to the important astrophysical implications of this issue, we present here an independent assessment of the process. For this purpose, we compute free-free quasi-molecular absorption in Lyman-alpha due to collisions with H and H2 within the one-perturber, quasi-static approximation. Line cross-sections are obtained using theoretical molecular potentials to describe the interaction between the radiating atom and the perturber. The variation of the electric-dipole transition moment with the interparticle distance is also considered. Six and two allowed electric dipole transitions due to H-H and H-H2 collisions, respectively, are taken into account. The new theoretical Lyman-alpha line profiles are then incorporated in our stellar atmosphere program for the computation of synthetic spectra and colours of DA-type white dwarfs. Illustrative model atmospheres and spectral energy distributions are computed, which show that Ly-alpha broadening by atoms and molecules has a significant effect on the white dwarf atmosphere models. The inclusion of this collision-induced opacity significantly reddens spectral energy distributions and affects the broadband colour indices for model atmospheres with Teff<5000 K. These results confirm those previously obtained by Kowalski & Saumon (2006). Our study points out the need for reliable evaluations of H3 potential energy surfaces covering a large region of nuclear configurations, in order to obtain a better description of H-H2 collisions and a more accurate evaluation of their influence on the spectrum of cool white dwarfs.

Motivated by the strong discrepancy between the main sequence turn-off age and the white dwarf cooling age in the metal-rich open cluster NGC 6791, we compute a grid of white dwarf evolutionary sequences that incorporates for the first time the energy released by the processes of 22Ne sedimentation and of carbon/oxygen phase separation upon crystallization. The grid covers the mass range from 0.52 to 1.0 Msun, and it is appropriate for the study of white dwarfs in metal-rich clusters. The evolutionary calculations are based on a detailed and self-consistent treatment of the energy released from these two processes, as well as on the employment of realistic carbon/oxygen profiles, of relevance for an accurate evaluation of the energy released by carbon/oxygen phase separation. We find that 22Ne sedimentation strongly delays the cooling rate of white dwarfs stemming from progenitors with high metallicities at moderate luminosities, whilst carbon/oxygen phase separation adds considerable delays at low luminosities. Cooling times are sensitive to possible uncertainties in the actual value of the diffusion coefficient of 22Ne. Changing the diffusion coefficient by a factor of 2, leads to maximum age differences of approx. 8-20% depending on the stellar mass. We find that the magnitude of the delays resulting from chemical changes in the core is consistent with the slow down in the white dwarf cooling rate that is required to solve the age discrepancy in NGC 6791.

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.

We present full evolutionary calculations appropriate for the study of hydrogen-rich DA white dwarfs. This is done by evolving white dwarf progenitors from the zero age main sequence, through the core hydrogen burning phase, the helium burning phase and the thermally pulsing asymptotic giant branch phase to the white dwarf stage. Complete evolutionary sequences are computed for a wide range of stellar masses and for two different metallicities: Z=0.01, which is representative of the solar neighborhood, and Z=0.001, which is appropriate for the study of old stellar systems, like globular clusters. During the white dwarf cooling stage we compute self-consistently the phase in which nuclear reactions are still important, the diffusive evolution of the elements in the outer layers and, finally, we also take into account all the relevant energy sources in the deep interior of the white dwarf, like the release of latent heat and the release of gravitational energy due to carbon-oxygen phase separation upon crystallization. We also provide colors and magnitudes for these sequences, based on a new set of improved non-gray white dwarf model atmospheres, which include the most up-to-date physical inputs like the Lyman alpha quasi-molecular opacity. The calculations are extended down to an effective temperature of 2,500 K. Our calculations provide a homogeneous set of evolutionary cooling tracks appropriate for mass and age determinations of old DA white dwarfs and for white dwarf cosmochronology of the different Galactic populations.

We present full evolutionary calculations appropriate for the study of hot hydrogen-deficent DO white dwarfs, PG 1159 stars, and DB white dwarfs. White dwarf sequences are computed for a wide range of stellar masses and helium envelopes on the basis of a complete treatment of the evolutionary history of progenitors stars, including the core hydrogen and helium burning phases, the thermally-pulsing AGB phase, and the born-again episode that is responsible for the hydrogen deficiency. We also provide colors and magnitudes for the new sequences for $T_{\rm eff} < 40 000$ K, where the NLTE effects are not dominant. These new calculations provide an homogeneous set of evolutionary tracks appropriate for mass and age determinations for both PG 1159 stars and DO white dwarfs. The calculations are extended down to an effective temperature of 7 000 K. We applied these new tracks to redetermine stellar masses and ages of all known DO white dwarfs with spectroscopically-determined effective temperatures and gravities, and compare them with previous results. We also compare for the first time consistent mass determinations for both DO and PG 1159 stars, and find a considerably higher mean mass for the DO white dwarfs. We discuss as well the chemical profile expected in the envelope of variable DB white dwarfs from the consideration of the evolutionary history of progenitor stars. Finally, we present tentative evidence for a different evolutionary channel, other than that involving the PG 1159 stars, for the formation of hot, hydrogen-deficient white dwarfs.

Motivated by the recent detection of single and binary He-core white dwarfs in metal-rich clusters, we present a full set of evolutionary calculations and colors appropriate for the study of such white dwarfs. The paper is also aimed at investigating whether stable hydrogen burning may constitute a main source of energy for massive He-core white dwarfs resulting from high-metallicity progenitors. White dwarf sequences are derived by taking into account the evolutionary history of progenitor stars with supersolar metallicities. We also incorporate a self-consistent, time-dependent treatment of gravitational settling and chemical diffusion, as well as of the residual nuclear burning. We find that the influence of residual nuclear burning during the late stages of white dwarf evolution is strongly dependent on the occurrence of chemical diffusion at the base of the hydrogen-rich envelope. When no diffusion is considered, residual hydrogen burning strongly influences the advanced stages of white dwarf cooling, introducing evolutionary delays of several Gyr. By contrast, when diffusion is taken into account the role of residual nuclear burning is strongly mitigated, and the evolution is dictated only by the thermal content stored in the ions. In addition, for all of our sequences, we provide accurate color and magnitudes on the basis of new and improved non gray model atmospheres which explicitly include Ly$\alpha$ quasi-molecular opacity.

We present evolutionary calculations and colors for massive white dwarfs with oxygen-neon cores for masses between 1.06 and 1.28 Mo. The evolutionary stages computed cover the luminosity range from log(L/Lo) approx. 0.5 down to -5.2. Our cooling sequences are based on evolutionary calculations that take into account the chemical composition expected from massive white dwarf progenitors that burned carbon in partially degenerate conditions. The use of detailed non-gray model atmospheres provides us with accurate outer boundary conditions for our evolving models at low effective temperatures. We examine the cooling age, colors and magnitudes of our sequences. We find that massive white dwarfs are characterized by very short ages to such an extent that they reach the turn-off in their colors and become blue at ages well below 10 Gyr. Extensive tabulations for massive white dwarfs, accessible from our web site, are also presented.

Splaver and coworkers have measured the masses of the white dwarf and the neutron star components of the PSR J1713+0747 binary system pair by Shapiro Delay. We attempt to find the original configuration of this system performing a set of binary evolution calculations to simultaneously account for the masses of both stars and the orbital period. We considered initial masses of 1.5 and 1.4 \msun for the normal (donor) and the neutron star, respectively. We assumed two metallicity values (Z = 0.010 and 0.020), and an initial orbital period near 3 days. We assume that the neutron star is only able to retain \lesssim 0.10 of the matter transferred by the donor star. Calculations were performed employing our binary hydro code that handles the mass transfer rate in a fully implicit way together with state-of-the-art physical ingredients, diffusion and a non-grey atmospheres. We compare the structure of the resulting white dwarfs with the characteristic age of PSR J1713+0747 finding a nice agreement with observations by Lundgren et al. especially for the case of a donor star with Z= 0.010. This result indicates that the evolution of this kind of binary system is well understood. The models predict that, due to diffusion, the atmosphere of the white dwarf is an almost hydrogen-pure one. We find that such structures are unable to account for the colours measured by Lundgren et al. within their error bars. Thus, some discrepances in the white dwarf emergent radiation remain to be explained.

The new mathematical framework based on the free energy of pure classical fluids presented in [R. D. Rohrmann, Physica A 347, 221 (2005)] is extended to multi-component systems to determine thermodynamic and structural properties of chemically complex fluids. Presently, the theory focuses on $D$-dimensional mixtures in the low-density limit (packing factor $\eta < 0.01$). The formalism combines the free-energy minimization technique with space partitions that assign an available volume $v$ to each particle. $v$ is related to the closeness of the nearest neighbor and provides an useful tool to evaluate the perturbations experimented by particles in a fluid. The theory shows a close relationship between statistical geometry and statistical mechanics. New, unconventional thermodynamic variables and mathematical identities are derived as a result of the space division. Thermodynamic potentials $\mu_{il}$, conjugate variable of the populations $N_{il}$ of particles class $i$ with the nearest neighbors of class $l$ are defined and their relationships with the usual chemical potentials $\mu_i$ are established. Systems of hard spheres are treated as illustrative examples and their thermodynamics functions are derived analytically. The low-density expressions obtained agree nicely with those of scaled-particle theory and Percus-Yevick approximation. Several pair distribution functions are introduced and evaluated. Analytical expressions are also presented for hard spheres with attractive forces due to K\^ac-tails and square-well potentials. Finally, we derive general chemical equilibrium conditions.

The present paper focuses on the evolution of hydrogen-deficient white dwarfs with the aim of exploring the consequences of different initial envelope structures on the carbon abundances expected in helium-rich, carbon-contaminated DQ white dwarfs. In particular, the evolutionary link between the DQs with low detected carbon abundances and the PG1159, extreme horizontal branch, and helium-rich R Coronae Borealis (RCrB) stars is explored. We present full evolutionary calculations that take a self-consistent treatment of element diffusion into account as well as expectations for the outer layer chemical stratification of progenitor stars upon entering the white dwarf regime. We find that PG1159 stars cannot be related to any DQ white dwarfs with low C abundances. Instead, we suggest that the latter could constitute the progeny of the giant, helium-rich RCrB stars.

An equilibrium theory of classical fluids based on the space distribution among the particles is derived in the framework of the energy minimization method. This study is motivated by current difficulties of evaluation of optical properties in atmospheres of degenerate stars. Present paper focuses on diluted one-component systems, where the interaction energy is calculated as a sum of binary contributions. The spatial configuration of the gas is described in terms of a particle-state variable v which roughly measures the free space surrounding each particle. The formalism offers a unified treatment of both thermodynamics and structure of fluids, since it not only provides the state function (the Helmholtz free energy) of a fluid, but also automatically gives information on the microstructure of the system (e.g. the nearest-neighbor distribution function). Equations of state and nearest-neighbor distribution functions of perfect and hard-body fluids are obtained in straightforward way in one, two and three dimensions. The formalism is applied to describe the atomic population in a partially ionized hydrogen gas. Combined and self-consistent evaluations of atomic populations and internal effects on bound states are performed in a detailed form. The present theory allow us to resolve the atomic density at each internal level into groups of atoms which experiment different perturbation intensities according to the size of their spaces v.

The present work is designed to explore the evolution of helium-core white dwarf (HeWD) stars for the case of metallicities much lower than the solar one (Z=0.001 and Z=0.0002). Evolution is followed in a self-consistent way with the predictions of detalied and new non-grey atmospheres, time-dependent element diffusion and the history of the white dwarf progenitor. Reliable initial models for low mass HeWDs are obtained by applying mass loss rates to a 1msun stellar model. The loss of angular momentum caused by gravitational wave emission and magnetic stellar wind braking are considered. Model atmospheres, based on a detailed treatment of the microphysics entering the WD atmosphere enable us to provide accurate colours and magnitudes at both early and advanced evolutionary stages. We find that most of our evolutionary sequences experience several episodes of hydrogen thermonuclear flashes. In particular, the lower the metallicity, the larger the minimum stellar mass for the occurrence fo flashes induced by CNO cycle reactions. The existence of a mass-threshold for the occurrence of diffusion-induced CNO flashes leadss to a marked dichotomy in the age of our models. Another finding of this study is that our HeWD models experience unstable hydrogen burning via PP nuclear reactions at late cooling stages as a result of hydrogen chemically diffusing inwards. Such PP flashes take place in models with very low metal content. We also find that models experiencing CNO flashes exhibit a pronouncede turn-off in most of their colours at M_V=16 approximately. Finally, colour-magnitude diagrams for our models are presented and compared with recent observational data of HeWD candidates in the globular clusters NGC 6397 and 47 Tucanae.

We examine the emergent fluxes from helium-core white dwarfs following their evolution from the end of pre-white dwarf stages down to advanced cooling stages. For this purpose, we include a detailed treatment of the physics of the atmosphere, particularly an improved representation of the state of the gas by taking into account non-ideal effects according to the so-called occupation probability formalism. The present calculations also incorporate hydrogen line opacity from Lyman, Balmer and Paschen series, pseudo-continuum absorptions and new updated induced-dipole absorption from H$_2$-H$_2$, H$_2$-He and H-He pairs. We find that the non-ideal effects and line absorption alter the appearance of the stellar spectrum and have a significant influence upon the photometric colours in the UBVRI-JHKL system. This occurs specially for hot models $T_{\rm eff}\ga 8000$ due to line and pseudo-continuum opacities, and for cool models $T_{\rm eff}\la 4000$ where the perturbation of atoms and molecules by neighbour particles affects the chemical equilibrium of the gas. In the present study, we also include new cooling sequences for helium-core white dwarfs of very low mass (0.160 and 0.148 M$_\odot$) with metallicity $Z=0.02$. These computations provide theoretical support to search for and identify white dwarfs of very low mass, specially useful for recent and future observational studies of globular cluster, where these objects have began to be detected.

The purpose of this work is to explore the evolution of helium-core white dwarf stars in a self-consistent way with the predictions of detailed non-gray model atmospheres and element diffusion. To this end, we consider helium-core white dwarf models with stellar masses of 0.406, 0.360, 0.327, 0.292, 0.242, 0.196 and 0.169 solar masses and follow their evolution from the end of mass loss episodes during their pre-white dwarf evolution down to very low surface luminosities. We find that when the effective temperature decreases below 4000K, the emergent spectrum of these stars becomes bluer within time-scales of astrophysical interest. In particular, we analyse the evolution of our models in the colour-colour and colour-magnitude diagrams and we find that helium-core white dwarfs with masses ranging from approx. 0.18 to 0.3 solar masses can reach the turn-off in their colours and become blue again within cooling times much less than 15 Gyr and then remain brighter than M_V approx. 16.5. In view of these results, many low-mass helium white dwarfs could have had time enough to evolve to the domain of collision-induced absorption from molecular hydrogen, showing blue colours.

We present a detailed calculation of model atmospheres for DA white dwarfs. Our atmosphere code solves the atmosphere structure in local thermodynamic equilibrium with a standard partial linearization technique, which takes into account the energy transfer by radiation and convection. This code incorporates recent improved and extended data base of collision induced absorption by molecular hydrogen. We analyse the thermodynamic structure and emergent flux of atmospheres in a range 2500 < Teff < 60000$ K and 6.5 < log g < 9.0. Bolometric correction and colour indices are provided for a subsample of the model grid. Comparison of the colours is made with published observational material and results of other recent model calculations. Motivated by the increasing interest on helium core white dwarfs, we analyse the photometric characteristics of these stars during their cooling, using evolutionary models recently available. Effective temperatures, surface gravities, masses and ages have been determined for some helium core white dwarf candidates, and their possible binary nature is briefly discussed.