White dwarf masses in magnetic cataclysmic variables: multi-temperature fits to Ginga data

One method of obtaining the mass of the white dwarf in magnetic cataclysmic variables (mCVs) is through their hard X-ray spectra. However, previous mass estimates using this method give lower limits because the temperature of the plasma in the post-shock region (where the hard X-rays are emitted) is lower than the temperature of the shock itself. In AM Her systems, the additional cooling of the post-shock plasma by cyclotron emission will further lower the derived mass. Here we present estimates of the masses of the white dwarf in 13 mCVs derived using Ginga data and a model in which X-rays are emitted from a multi-temperature emission region with the appropriate temperature and density profile. We include in the model reflection from the surface of the white dwarf and a partially ionized absorber. We are able to achieve good fits to the data. We compare the derived masses with previous estimates and the masses for larger samples of isolated white dwarfs and those in CVs.     

X-ray emissions from two-temperature accretion flows within a dipole magnetic funnel

We investigate the hydrodynamics of accretion channelled by a dipolar magnetic field (funnel flows). We consider situations in which the electrons and ions in the flow cannot maintain thermal equilibrium [two-temperature (2T) effects] due to strong radiative loss, and determine the effects on the keV X-ray properties of the systems. We apply this model to investigate the accretion shocks of white dwarfs in magnetic cataclysmic variables (mCVs). We have found that the incorporation of 2T effects could harden the keV X-rays. Also, the dipolar model yields harder X-ray spectra than the standard planar model if white dwarf is sufficiently massive  (≳1 M_⊙)  . When fitting observed keV X-ray spectra of mCVs, the inclusion of 2T hydrodynamics and a dipolar accretion geometry lowers estimates for white dwarf masses when compared with masses inferred from models excluding these effects. We find mass reductions ≲9 per cent in the most massive cases.     

A source of high-velocity white dwarfs

We investigate whether the recently observed population of high-velocity white dwarfs can be derived from a population of binaries residing initially within the thin disc of the Galaxy. In particular, we consider binaries where the primary is sufficiently massive to explode as a Type II supernova. A large fraction of such binaries are broken up when the primary then explodes as a supernova, owing to the combined effects of the mass loss from the primary and the kick received by the neutron star on its formation. For binaries where the primary evolves to fill its Roche lobe, mass transfer from the primary leads to the onset of a common envelope phase during which the secondary and the core of the primary spiral together as the envelope is ejected. Such binaries are the progenitors of X-ray binaries if they are not broken up when the primary explodes. For those systems that are broken up, a large number of the secondaries receive kick velocities ∼100–200 km s⁻¹ and subsequently evolve into white dwarfs. We compute trajectories within the Galactic potential for this population of stars and relate the birth rate of these stars over the entire Galaxy to those seen locally with high velocities relative to the local standard of rest (LSR) . We show that for a reasonable set of assumptions concerning the Galactic supernova rate and the binary population, our model produces a local number density of high-velocity white dwarfs compatible with that inferred from observations. We therefore propose that a population of white dwarfs originating in the thin disc may make a significant contribution to the observed population of high-velocity white dwarfs.     

Two-temperature accretion flows in magnetic cataclysmic variables: structures of post-shock emission regions and X-ray spectroscopy

We use a two-temperature hydrodynamical formulation to determine the temperature and density structures of the post-shock accretion flows in magnetic cataclysmic variables (mCVs) and calculate the corresponding X-ray spectra. The effects of two-temperature flows are significant for systems with a massive white dwarf and a strong white-dwarf magnetic field. Our calculations show that two-temperature flows predict harder keV spectra than one-temperature flows for the same white-dwarf mass and magnetic field. This result is insensitive to whether the electrons and ions have equal temperature at the shock, but depends on the electron–ion exchange rate, relative to the rate of radiative loss along the flow. White-dwarf masses obtained by fitting the X-ray spectra of mCVs using hydrodynamic models including the two-temperature effects will be lower than those obtained using single-temperature models. The bias is more severe for systems with a massive white dwarf.     

The origin of the initial mass function and its dependence on the mean Jeans mass in molecular clouds

We investigate the dependence of stellar properties on the mean thermal Jeans mass in molecular clouds. We compare the results from the two largest hydrodynamical simulations of star formation to resolve the fragmentation process down to the opacity limit, the first of which was reported by Bate, Bonnell & Bromm. The initial conditions of the two calculations are identical except for the radii of the clouds, which are chosen so that the mean densities and mean thermal Jeans masses of the clouds differ by factors of 9 and 3, respectively.
We find that the denser cloud, with the lower mean thermal Jeans mass, produces a higher proportion of brown dwarfs and has a lower characteristic (median) mass of the stars and brown dwarfs. This dependence of the initial mass function (IMF) on the density of the cloud may explain the observation that the Taurus star-forming region appears to be deficient in brown dwarfs when compared with the Orion Trapezium cluster. The new calculation also produces wide binaries (separations >20 au), one of which is a wide binary brown dwarf system.
Based on the hydrodynamical calculations, we develop a simple accretion/ejection model for the origin of the IMF. In the model, all stars and brown dwarfs begin with the same mass (set by the opacity limit for fragmentation) and grow in mass until their accretion is terminated stochastically by their ejection from the cloud through dynamically interactions. The model predicts that the main variation of the IMF in different star-forming environments should be in the location of the peak (due to variations in the mean thermal Jeans mass of the cloud) and in the substellar regime. However, the slope of the IMF at high masses may depend on the dispersion in the accretion rates of protostars.     

The brown dwarf desert as a consequence of orbital migration

We show that the dearth of brown dwarfs in short-period orbits around Solar-mass stars – the brown dwarf desert – can be understood as a consequence of inward migration within an evolving protoplanetary disc. Brown dwarf secondaries forming at the same time as the primary star have masses which are comparable to the initial mass of the protoplanetary disc. Subsequent disc evolution leads to inward migration, and destruction of the brown dwarf, via merger with the star. This is in contrast with massive planets, which avoid this fate by forming at a later epoch when the disc is close to being dispersed. Within this model, a brown dwarf desert arises because the mass at the hydrogen-burning limit is coincidentally comparable to the initial disc mass for a Solar mass star. Brown dwarfs should be found in close binaries around very low mass stars, around other brown dwarfs, and around Solar-type stars during the earliest phases of star formation.     

Ionization structure in accretion shocks with a composite cooling function

We have investigated the ionization structure of the post-shock regions of magnetic cataclysmic variables, using an analytic density and temperature structure model in which effects caused by bremsstrahlung and cyclotron cooling are considered. We find that in the majority of the shock-heated region where H- and He-like lines of the heavy elements are emitted, the collisional-ionization and corona-condition approximations are justified. We have calculated the line emissivity and ionization profiles for iron as a function of height within the post-shock flow. For low-mass white dwarfs, line emission takes place near the shock. For high-mass white dwarfs, most of the line emission takes place in regions well below the shock and hence it is less sensitive to the shock temperature. Thus, the line ratios are useful to determine the white dwarf masses for the low-mass white dwarfs, but the method is less reliable when the white dwarfs are massive. Line spectra can, however, be used to map the hydrodynamic structure of the post-shock accretion flow.     

The lower boundary of the accretion column in magnetic cataclysmic variables

Using a parametrized function for the mass loss at the base of the post-shock region, we have constructed a formulation for magnetically confined accretion flows which avoids singularities, such as the infinity in density, at the base associated with all previous formulations. With the further inclusion of a term allowing for the heat input into the base from the accreting white dwarf, we are also able to obtain the hydrodynamic variables to match the conditions in the stellar atmosphere. (We do not, however, carry out a mutually consistent analysis for the match.) Changes to the emitted X-ray spectra are negligible unless the thickness of mass leakage region at the base approaches or exceeds one per cent of the height of the post-shock region. In this case the predicted spectra from higher-mass white dwarfs will be harder, and fits to X-ray data will predict lower white dwarf masses than previous formulations.     

Accretion and dynamical interactions in small-N star-forming clusters: N= 5

We present results from high-resolution hydrodynamical simulations that explore the effects of small-scale clustering in star-forming regions. A large ensemble of small- N clusters with five stellar seeds have been modelled and the resulting properties of stars and brown dwarfs statistically derived and compared with observational data.
Close dynamical interactions between the protostars and competitive accretion driven by the cloud collapse are shown to produce a distribution of final masses that is bimodal, with most of the mass residing in the binary components. When convolved with a suitable core mass function, the final distribution of masses resembles the observed initial mass function, in both the stellar and substellar regimes. Binaries and single stars are found to constitute two kinematically distinct populations, with about half of the singles attaining velocities ≥2 km s⁻¹, which might deprive low-mass star-forming regions of their lightest members in a few crossing times. The eccentricity distribution of binaries and multiples is found to follow a distribution similar to that of observed long-period (uncircularized) binaries.
The results obtained support a mechanism in which a significant fraction of brown dwarfs form under similar circumstances as those of normal stars but are ejected from the common envelope of unstable multiple systems before their masses exceed the hydrogen burning limit. We predict that many close binary stars should have wide brown dwarf companions. Brown dwarfs, and, in general, very low-mass stars, would be rare as pure binary companions. The binary fraction should be a decreasing function of primary mass, with low-mass or substellar primaries being scarce. Where such binaries exist, they are expected either to be close enough (semimajor axis ∼10 au) to survive strong interactions with more massive binaries or to be born in very small molecular cloud cores.     

Relation of Velocity Distribution and Mass for DA White Dwarfstwo

White dwarfs are the evolutionary endpoint of the low-and-medium mass stars. In the studies of white dwarfs, the mass of white dwarf is an important physical parameter. In this paper, we give an analysis about the velocity distribution of DA white dwarfs in the Sloan Digital Sky Survey (SDSS), and hope to find the relation between mass and velocity distribution of white dwarfs. We get the radial velocity and tangential velocity of every DA white dwarf according to their proper motion and spectral shift. Through analyzing the velocity distribution of DA white dwarfs, we find that the small-mass white dwarfs, which are produced from the single-star evolution channel, have a relatively large velocity dispersion.     

Magnetic cataclysmic variables: optical counterparts of X-ray and EUV sources

The magnetic cataclysmic variables (mCVs) are a small but astrophysically important class of X-ray source. This paper discusses some of the work carried out by the author and his collaborators on these objects, which were discovered as part of various optical identification programs of both X-ray and extreme ultraviolet sources observed by the HEAO-1 and ROSAT satellites. Definitions are given of the two classes of mCVs; the Polars and the Intermediate Polars. The reasons why mCVs are significant sources of X-ray and EUV radiation are also discussed. Individual systems which have been the subject of intensive follow-up programmes at SAAO are then described, and the implications arising from these observations are discussed.     

Evolutionary calculations of carbon dredge-up in helium envelope white dwarfs

We investigate the evolution of cooling helium atmosphere white dwarfs using a full evolutionary code, specifically developed to follow the effects of element diffusion and gravitational settling on white dwarf cooling. The major difference between this work and previous work is that we use more recent opacity data from the OPAL project. Since, in general, these opacities are higher than those available 10 years ago, at a given effective temperature, convection zones go deeper than in models with older opacity data. Thus convective dredge-up of observationally detectable carbon in helium atmosphere white dwarfs can occur for thicker helium layers than found by Pelletier et al. We find that the range of observed C to He ratios in different DQ white dwarfs of similar effective temperature is well explained by a range of initial helium layer mass between 10⁻³ and 10⁻² M⊙, in good agreement with stellar evolution theory, assuming a typical white dwarf mass of 0.6 M⊙. We also predict that oxygen will be present in DQ white dwarf atmospheres in detectable amounts if the helium layer mass is near the lower limit compatible with stellar evolution theory. Determination of the oxygen abundance has the potential of providing information on the profile of oxygen in the core and hence on the important ¹²C(α,γ)¹⁶O reaction rate.     

Low mass companions to white dwarfs

This paper summarizes the results of over 17 years of work searching for low mass stellar and substellar companions to more than 370 nearby white dwarfs. Roughly 60 low mass, unevolved companions were found and studied all together, with over 20 discovered in the last few years, including the first unambiguous brown dwarf companion to a white dwarf, GD 1400B. The resulting spectral type distributions for companions to white dwarfs and nearby cool field dwarfs are compared, and the implications for binary star formation are discussed. A brief analysis of GD 1400B, including new data, is also presented. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The origin of the magnetic fields in white dwarfs

Magnetic white dwarfs with fields in excess of ∼10⁶ G (the high field magnetic white dwarfs; HFMWDs) constitute about ∼10 per cent of all white dwarfs and show a mass distribution with a mean mass of  ∼0.93 M_⊙  compared to  ∼0.56 M_⊙  for all white dwarfs. We investigate two possible explanations for these observations. First, that the initial–final mass relationship (IFMR) is influenced by the presence of a magnetic field and that the observed HFMWDs originate from stars on the main sequence that are recognized as magnetic (the chemically peculiar A and B stars). Secondly, that the IFMR is essentially unaffected by the presence of a magnetic field, and that the observed HFMWDs have progenitors that are not restricted to these groups of stars. Our calculations argue against the former hypothesis and support the latter. The HFMWDs have a higher than average mass because on the average they have more massive progenitors and not because the IFMR is significantly affected by the magnetic field. A requirement of our model is that ∼40 per cent of main-sequence stars more massive than  ∼4.5 M_⊙  must either have magnetic fields in the range of ∼10–100 G, which is below the current level of detection, or generate fields during subsequent stellar evolution towards the white dwarf phase. In the former case, the magnetic fields of the HFMWDs could be fossil remnants from the main-sequence phase consistent with the approximate magnetic flux conservation.     

The impact of element diffusion on the formation and evolution of helium white dwarf stars

The aim of this work is to investigate the effect of element diffusion on the evolution of helium white dwarfs. To this end, we couple the multicomponent flow equations that describe gravitational settling, chemical and thermal diffusion to an evolutionary code. We compute the evolution of a set of helium white dwarf models with masses ranging from 0.169 to 0.406 M_⊙. In particular, several low-mass white dwarfs have been found in binary systems as companion to millisecond pulsars. In these systems, pulsar emission is activated by mass transfer episodes so that, if we place the zero-age point at the end of such mass transfer, then the pulsar and the white dwarf ages should be equal. Interestingly enough, available models of helium white dwarfs neglect element diffusion. Using such models, good agreement has been found between the ages of the components of the PSR J1012+5307 system. However, recent observations of the PSR B1855+09 system cast doubts on the correctness of such models, which predict a white dwarf age twice as long as the spin-down age of the pulsar. In this work, we find that element diffusion induces thermonuclear hydrogen shell flashes for models in the mass interval 0.18≲ M /M_⊙ ≲ 0.41 . We show, in particular, that the occurrence of these diffusion-induced flashes eventually leads to white dwarf models with hydrogen envelope masses too small to support any further nuclear burning, thus implying much shorter cooling ages than in the case when diffusion is neglected. In particular, excellent agreement is found between the ages of PSR B1855+09 system components, solving the age discrepancy from first principles.     

A new magnetic white dwarf: PG 2329+267

We have discovered that the white dwarf PG 2329+267 is magnetic, and, assuming a centred dipole structure, has a dipole magnetic field strength of approximately 2.3 MG. This makes it one of only approximately 4 per cent of isolated white dwarfs with a detectable magnetic field. Linear Zeeman splitting, as well as quadratic Zeeman shifts, is evident in the hydrogen Balmer sequence and circular spectropolarimetry reveals ∼10 per cent circular polarization in the two displaced σ components of Hα. We suggest from comparison with spectra of white dwarfs of known mass that PG 2329+267 is more massive than typical isolated white dwarfs, in agreement with the hypothesis that magnetic white dwarfs evolve from magnetic chemically peculiar Ap and Bp type main-sequence stars.     

The possiblity of detection of ultracool dwarfs with the UKIRT Infrared Deep Sky Survey

We present predictions for the numbers of ultracool dwarfs in the Galactic disc population that could be detected by the WFCAM/UKIDSS Large Area Survey and Ultra Deep Survey. Simulated samples of objects are created with masses and ages drawn from different mass functions and birthrates. Each object is then given absolute magnitudes in different passbands based on empirically derived bolometric correction versus effective temperature relationships (or model predictions for Y dwarfs). These are then combined with simulated space positions, velocities and photometric errors to yield observables such as apparent magnitudes and proper motions. Such observables are then passed through the survey selection mechanism to yield histograms in colour. This technique also produces predictions for the proper motion histograms for ultracool dwarfs and estimated numbers for the as yet undetected Y dwarfs. Finally, it is shown that these techniques could be used to constrain the ultra-low-mass mass function and birthrate of the Galactic disc population.     

The effect of photospheric heavy elements on the hot DA white dwarf temperature scale

Using the latest non-local thermodynamic equilibrium (non-LTE) synthetic spectra and stellar model calculations, we have evaluated the potential effect of the presence of heavy elements in the photospheres of hot H-rich DA white dwarfs. In particular, we have examined their influence on the effective temperature and surface gravity perceived from analysis of the Balmer line profiles. It is apparent that both the inclusion of non-LTE effects in the models and significant quantities of heavy elements act independently to lower the value of T _eff determined from a particular spectrum. Hence, the true effective temperatures of the heavy element-rich DA white dwarfs, currently estimated to be above 55 000 K, are apparently lower than previously reported from pure-H LTE analyses, by some 4000–7000 K. We do not see any similar influence on measurements of log g . This work concentrates on a group of relatively bright well-studied objects, for which heavy element abundances are known. As a consequence of this, establishment of correct temperatures for all other hot white dwarfs will require a programme of far-UV spectroscopy in order to obtain the essential compositional information. Since only stars with effective temperatures lying notionally in the range from ≈ 55 000 to 70 000 K (52 000–62 000 K when the non-LTE effects and heavy elements are taken into account) have been considered here, important questions remain regarding the magnitude of any similar effects in even hotter white dwarfs and pre-white dwarfs. The resulting implications for the plausibility of the evolutionary link between the main hot DA population and their proposed precursors, the H-rich central stars of planetary nebulae, need to be investigated.     

The initial–final mass relationship of white dwarfs revisited: effect on the luminosity function and mass distribution

The initial–final mass relationship connects the mass of a white dwarf with the mass of its progenitor in the main sequence. Although this function is of fundamental importance to several fields in modern astrophysics, it is not well constrained either from the theoretical or from the observational points of view. In this work, we revise the present semi-empirical initial–final mass relationship by re-evaluating the available data. The distribution obtained from grouping all our results presents a considerable dispersion, which is larger than the uncertainties. We have carried out a weighted least-squares linear fit of these data and a careful analysis to give some clues on the dependence of this relationship on some parameters such as metallicity or rotation. The semi-empirical initial–final mass relationship arising from our study covers the range of initial masses from 1.0 to  6.5 M_⊙  , including in this way the low-mass domain, poorly studied until recently. Finally, we have also performed a test of the initial–final mass relationship by studying its effect on the luminosity function and on the mass distribution of white dwarfs. This was done by using different initial–final mass relationships from the literature, including the expression derived in this work, and comparing the results obtained with the observational data from the Palomar Green Survey and the Sloan Digital Sky Survey. We find that the semi-empirical initial–final mass relationship derived here gives results in good agreement with the observational data, especially in the case of the white dwarf mass distribution.     

The evolution of a rapidly accreting helium white dwarf to become a low-luminosity helium star

We have examined the evolution of merged low-mass double white dwarfs which become low-luminosity (or high-gravity) extreme helium stars. We have approximated the merging process by the rapid accretion of matter, consisting mostly of helium, on to a helium white dwarf. After a certain mass is accumulated, a helium shell flash occurs, the radius and luminosity increase and the star becomes a yellow giant. Mass accretion is stopped artificially when the total mass reaches a pre-determined value. As the helium-burning shell moves inwards with repeating shell flashes, the effective temperature gradually increases as the star evolves towards the helium main sequence. When the mass interior to the helium‐burning shell is approximately 0.25 M_⊙, the star enters a regime where it is pulsationally unstable. We have obtained radial pulsation periods for these models.
These models have properties very similar to those of the pulsating helium star V652 Her. We have compared the rate of period change of the theoretical models with that observed in V652 Her, as well as with its position on the Hertzsprung–Russell diagram. We conclude that the merger between two helium white dwarfs can produce a star with properties remarkably similar to those observed in at least one extreme helium star, and is a viable model for their evolutionary origin. Such helium stars will evolve to become hot subdwarfs close to the helium main sequence. We also discuss the number of low-luminosity helium stars in the Galaxy expected for our evolution scenario.