Search

We present a new magnetic-atmosphere model code for obtaining synthetic spectral fluxes of hydrogen-rich magnetic white dwarfs. To date, observed spectra have been analyzed with models that neglect the magnetic field's effects on the atomic populations. In this work, we incorporate state-of-art theory in the evaluation of numerical densities of atoms, free electrons, and ions in local thermodynamical equilibrium under the action of a magnetic field. The energy distribution of atoms is rigorously evaluated for arbitrary field strength. This energy pattern includes going from tightly bound states to metastable or truly bound, highly excited states embedded in the continuum, that is, over the first Landau level. Finite nuclear mass effects and the coupling between the internal atomic structure and the motion of the atom across the magnetic field are also considered. Synthetic fluxes are generated with integrations of numerical solutions of polarized radiative transfer over the visible stellar disk using a spherical t-design method. The atmosphere code is tested with observations from the Sloan Digital Sky Survey for a group of known magnetic white dwarfs. Physical stellar parameters are obtained from least-squares fits to the observed energy distribution and compared with results of previous works. We show that the use of zerofield ionization equilibrium in spectral analyses can lead to underestimated effective temperatures for highly magnetic white dwarfs., Comment: 14 pages, 11 figures. Accepted for publication in Astronomy & Astrophysics

We report a peculiar variable blue star in the globular cluster NGC 6397, using Hubble Space Telescope optical imaging. Its position in the colour-magnitude diagrams, and its spectrum, are consistent with this star being a helium core white dwarf (He WD) in a binary system. The optical light curve shows a periodicity at 18.5 hours. We argue that this periodicity is due to the rotation of the WD and possibly due to magnetic spots on the surface of the WD. This would make this object the first candidate magnetic He WD in any globular cluster (GC), and the first candidate magnetic WD in a detached binary system in any GC and one of the few He WDs with a known rotation period and of magnetic nature. Another possibility is that this system is a He WD in a binary system with another WD or another degenerate object, which would make this object one of the few candidate non-accreting double degenerate binaries in any GC., Comment: Accepted for publication in MNRAS

The non-relativistic cross section from Rayleigh scattering by hydrogen atoms in the ground state is calculated over a wide range of photon energies ($< 0.8$ keV). Evaluations are performed in terms of the real and imaginary components of the atomic polarizability. The sum over intermediate states that characterizes this second-order radiative process is performed using exact analytic expressions for oscillator strengths of bound and continuum states. Damping terms associated with the finite lifetimes of excited states and their splitting into two fine-structure levels ($p_{1/2}$ and $p_{3/2}$) are taken into account in resonance cross sections. Fitting formulas required for cross-section evaluation are presented for incident photon energy i) redward of the first resonance (Lyman-$\alpha_{1/2}$), ii) in the spectral region corresponding to resonances (for an arbitrary number of them), and iii) above the ionization threshold., Comment: 8 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics

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., Comment: 15 pages, 10 figures. Accepted for publication in Astronomy & Astrophysics

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., Comment: 13 pages, 15 figures, 4 tables. Accepted for publication in Astronomy & Astrophysics

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., Comment: Accepted for publication in MNRAS; 14 pages, 14 figures

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., Comment: 40 pages, 13 figures. To be published in ApJ

We use the Percus-Yevick approach in the chemical-potential route to evaluate the equation of state of hard hyperspheres in five dimensions. The evaluation requires the derivation of an analytical expression for the contact value of the pair distribution function between particles of the bulk fluid and a solute particle with arbitrary size. The equation of state is compared with those obtained from the conventional virial and compressibility thermodynamic routes and the associated virial coefficients are computed. The pressure calculated from all routes is exact up to third density order, but it deviates with respect to simulation data as density increases, the compressibility and the chemical-potential routes exhibiting smaller deviations than the virial route. Accurate linear interpolations between the compressibility route and either the chemical-potential route or the virial one are constructed., Comment: 9 pages, 6 figures; v2: Change in one reference

The coupling-parameter method, whereby an extra particle is progressively coupled to the rest of the particles, is applied to the sticky-hard-sphere fluid to obtain its equation of state in the so-called chemical-potential route ($\mu$ route). As a consistency test, the results for one-dimensional sticky particles are shown to be exact. Results corresponding to the three-dimensional case (Baxter's model) are derived within the Percus-Yevick approximation by using different prescriptions for the dependence of the interaction potential of the extra particle on the coupling parameter. The critical point and the coexistence curve of the gas-liquid phase transition are obtained in the $\mu$ route and compared with predictions from other thermodynamics routes and from computer simulations. The results show that the $\mu$ route yields a general better description than the virial, energy, compressibility, and zero-separation routes., Comment: 13 pages, 7 figures; v2: Results from the zero-separation route have been added

The chemical potentials of multicomponent fluids are derived in terms of the pair correlation functions for arbitrary number of components, interaction potentials, and dimensionality. The formally exact result is particularized to hard-sphere mixtures with zero or positive nonadditivity. As a simple application, the chemical potentials of three-dimensional additive hard-sphere mixtures are derived from the Percus-Yevick theory and the associated equation of state is obtained. This Percus-Yevick chemical-route equation of state is shown to be more accurate than the virial equation of state. An interpolation between the chemical-potential and compressibility routes exhibits a better performance than the well-known Boubl\'ik-Mansoori-Carnahan-Starling-Leland equation of state., Comment: 9 pages, 1 figure; v2: minor changes

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., Comment: To be published in the proceedings of the 18th European White Dwarf Workshop. 4 pages, 1 figure

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., Comment: 12 pages, 11 figures, v2: new section, appendix and references, plus some other minor changes; to be published in Phys. Rev. E

Structural and thermodynamic properties of multicomponent hard-sphere fluids at odd dimensions have recently been derived in the framework of the rational function approximation (RFA) [Rohrmann and Santos, Phys. Rev. E \textbf{83}, 011201 (2011)]. It is demonstrated here that the RFA technique yields the exact solution of the Percus-Yevick (PY) closure to the Ornstein-Zernike (OZ) equation for binary mixtures at arbitrary odd dimensions. The proof relies mainly on the Fourier transforms $\hat{c}_{ij}(k)$ of the direct correlation functions defined by the OZ relation. From the analysis of the poles of $\hat{c}_{ij}(k)$ we show that the direct correlation functions evaluated by the RFA method vanish outside the hard core, as required by the PY theory., Comment: 6 pages

Mixtures of hard hyperspheres in odd space dimensionalities are studied with an analytical approximation method. This technique is based on the so-called Rational Function Approximation and provides a procedure for evaluating equations of state, structure factors, radial distribution functions, and direct correlations functions of additive mixtures of hard hyperspheres with any number of components and in arbitrary odd-dimension space. The method gives the exact solution of the Ornstein--Zernike equation coupled with the Percus--Yevick closure, thus extending to arbitrary odd dimension the solution for hard-sphere mixtures [J. L. Lebowitz, Phys.\ Rev.\ \textbf{133}, 895 (1964)]. Explicit evaluations for binary mixtures in five dimensions are performed. The results are compared with computer simulations and a good agreement is found., Comment: 16 pages, 8 figures; v2: slight change of notation

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., Comment: 10 pages, 6 figures, to be published in The Astrophysical Journal

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., Comment: 3 pages, 2 figures, published in Nature

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., Comment: 36 pages, 9 figures. Accepted for publication in The Astrophysical Journal

A recently derived method [R. D. Rohrmann and A. Santos, Phys. Rev. E. {\bf 76}, 051202 (2007)] to obtain the exact solution of the Percus-Yevick equation for a fluid of hard spheres in (odd) $d$ dimensions is used to investigate the convergence properties of the resulting virial series. This is done both for the virial and compressibility routes, in which the virial coefficients $B_j$ are expressed in terms of the solution of a set of $(d-1)/2$ coupled algebraic equations which become nonlinear for $d \geq 5$. Results have been derived up to $d=13$. A confirmation of the alternating character of the series for $d\geq 5$, due to the existence of a branch point on the negative real axis, is found and the radius of convergence is explicitly determined for each dimension. The resulting scaled density per dimension $2 \eta^{1/d}$, where $\eta$ is the packing fraction, is wholly consistent with the limiting value of 1 for $d \to \infty$. Finally, the values for $B_j$ predicted by the virial and compressibility routes in the Percus-Yevick approximation are compared with the known exact values [N. Clisby and B. M. McCoy, J. Stat. Phys. {\bf 122}, 15 (2006)], Comment: 9 pages, 6 figures

The structural properties of single component fluids of hard hyperspheres in odd space dimensionalities $d$ are studied with an analytical approximation method that generalizes the Rational Function Approximation earlier introduced in the study of hard-sphere fluids [S. B. Yuste and A. Santos, Phys. Rev. A {\bf 43}, 5418 (1991)]. The theory makes use of the exact form of the radial distribution function to first order in density and extends it to finite density by assuming a rational form for a function defined in Laplace space, the coefficients being determined by simple physical requirements. Fourier transform in terms of reverse Bessel polynomials constitute the mathematical framework of this approximation, from which an analytical expression for the static structure factor is obtained. In its most elementary form, the method recovers the solution of the Percus-Yevick closure to the Ornstein-Zernike equation for hyperspheres at odd dimension. The present formalism allows one to go beyond by yielding solutions with thermodynamic consistency between the virial and compressibility routes to any desired equation of state. Excellent agreement with available computer simulation data at $d=5$ and $d=7$ is obtained. As a byproduct of this study, an exact and explicit polynomial expression for the intersection volume of two identical hyperspheres in arbitrary odd dimensions is given., Comment: 18 pages, 7 figures; v2: new references added plus minor changes; to be published in PRE