Search

Magnetic activity and rotation are related to the age of low-mass main-sequence stars. To further constrain these relations, we study a sample of 574 main-sequence stars members of common proper motion pairs with white dwarfs, identified thanks to Gaia astrometry. We use the white dwarfs as age indicators, while the activity indexes and rotational velocities are obtained from the main-sequence companions using standard procedures. We find that stars older than 5 Gyr do not display H a nor Ca II H&K emission unless they are fast rotators due to tidal locking from the presence of unseen companions and that the rotational velocities tend to decrease over time, thus supporting the so-called gyrochronology. However, we also find moderately old stars (2–6 Gyr) that are active presumably because they rotate faster than they should for their given ages. This indicates that they may be suffering from weakened magnetic braking or that they possibly evolved through wind accretion processes in the past. The activity fractions that we measure for all stars younger than 5 Gyr range between 10 and 40 per cent. This is line with the expectations, since our sample is composed of F, G, K, and early M stars, which are thought to have short (<2 Gyr) activity lifetimes. Finally, we observe that the H a fractional luminosities and the RHK indexes for our sample of (slowly rotating) stars show a spread (-4 >log(LH a/Lbol); log(RHK) > -5) typically found in inactive M stars or weakly active/inactive F, G, K stars., This work has been partially supported by the (Spanish Ministerio de Asuntos Economicos ´ y Transformacion´ Digital) MINECO grant PID2020-117252GB-I00 and by the AGAUR/Generalitat de Catalunya grant SGR-386/2021. RR acknowledges support from Grant RYC2021-030837-I funded by MCIN/(Agencia Estatal de Investigacion) ´ AEI/ 10.13039/501100011033 and by ‘European Union NextGenerationEU/PRTR’. MEC acknowledges grant RYC2021-032721-I, funded by MCIN/AEI/10.13039/501100011033 and by the ‘European Union NextGenerationEU/PRTR’. MAH was supported by grant ST/V000853/1 from the Science and Technology Facilities Council (STFC). This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 101020057)., Peer Reviewed, Postprint (published version)

One of the key research themes identified by the Astro2020 decadal survey is Worlds and Suns in Context. The Advanced X-ray Imaging Satellite (AXIS) is a proposed NASA APEX mission that will become the prime high-energy instrument for studying star-planet connections from birth to death. This work explores the major advances in this broad domain of research that will be enabled by the AXIS mission, through X-ray observations of stars in clusters spanning a broad range of ages, flaring M-dwarf stars known to host exoplanets, and young stars exhibiting accretion interactions with their protoplanetary disks. In addition, we explore the ability of AXIS to use planetary nebulae, white dwarfs, and the Solar System to constrain important physical processes from the microscopic (e.g., charge exchange) to the macroscopic (e.g., stellar wind interactions with the surrounding interstellar medium)., Comment: This White Paper is part of a series commissioned for the AXIS Probe Concept Mission

We present an initial-final mass relation derived from the spectroscopically-complete volume-limited 40 pc sample of white dwarfs. The relation is modelled using population synthesis methods to derive an initial stellar population which can be fit to the observed mass distribution of white dwarfs. The population synthesis accounts for binary evolution, where higher-mass white dwarfs are more likely to be merger products than their lower-mass counterparts. Uncertainties are accounted from the initial mass function, stellar metallicity and age of the Galactic disc. We also consider biases induced by the spectral type of the white dwarf where pure-hydrogen atmosphere white dwarfs are likely to have more accurate masses, whilst the full white dwarf sample will have fewer biases arising from spectral evolution. We provide a four-piece segmented linear regression using Monte Carlo methods to sample the 1-$\sigma$ range of uncertainty on the initial stellar population. The derived initial-final mass relation provides a self-consistent determination of the progenitor mass for white dwarfs in the Solar neighbourhood which will be useful to study the local stellar formation history., Comment: Eleven pages. Accepted for publication in MNRAS

The atmospheres of a large proportion of white dwarf stars are polluted by heavy elements that are expected to sink out of visible layers on short timescales. This has been interpreted as a signature of ongoing accretion of debris from asteroids, comets, and giant planets. This scenario is supported by the detection of debris discs and transits of planetary fragments around some white dwarfs. However, photospheric metals are only indirect evidence for ongoing accretion, and the inferred accretion rates and parent body compositions heavily depend on models of diffusion and mixing processes within the white dwarf atmosphere. Here we report a 4.4$\sigma$ detection of X-rays from a polluted white dwarf, G29$-$38, using a 106 ks exposure with the Chandra X-ray Observatory, demonstrating directly that the star is currently accreting. From the measured X-ray luminosity, we find an instantaneous accretion rate of $\dot{M_{\rm X}}=1.63^{+1.29}_{-0.40}\times 10^{9}\mathrm{\,g\,s^{-1}}$. This is the first direct measurement of the accretion rate onto the white dwarf, which is independent of stellar atmosphere models. This rate exceeds estimates based on past studies of the photospheric abundances by more than a factor two, and implies that convective overshoot has to be accounted for in modelling the spectra of debris-accreting white dwarfs. We measure a low plasma temperature of $kT=0.5\pm0.2\,\mathrm{keV}$, corroborating the predicted bombardment solution for white dwarfs accreting at low accretion rates. Offering a new method for studying evolved planetary systems, these observations provide the opportunity to independently measure the instantaneous accretion rate of planetary material, and therefore investigate the timescale of accretion onto white dwarfs, and the evolution and replenishment of debris disks., Comment: Author submitted version. Published in Nature, Feb 10 2022 issue

]Context] Kinematic and chemical tagging of stellar populations have both revealed much information on the past and recent history of the Milky Way, including its formation history, merger events, and mixing of populations across the Galactic disk and halo. [Aims] We present the first detailed 3D kinematic analysis of a sample of 3133 white dwarfs that used Gaia astrometry plus radial velocities, which were measured either by Gaia or by ground-based spectroscopic observations. The sample includes either isolated white dwarfs that have direct radial velocity measurements, or white dwarfs that belong to common proper motion pairs that contain nondegenerate companions with available radial velocities. A subset of common proper motion pairs also have metal abundances that have been measured by large-scale spectroscopic surveys or by our own follow-up observations. [Methods] We used the white dwarfs as astrophysical clocks by determining their masses and total ages through interpolation with dedicated evolutionary models. We also used the nondegenerate companions in common proper motions to chemically tag the population. Combining accurate radial velocities with Gaia astrometry and proper motions, we derived the velocity components of our sample in the Galactic rest frame and their Galactic orbital parameters. [Results] The sample is mostly located within ∼300 pc from the Sun. It predominantly contains (90-95%) thin-disk stars with almost circular Galactic orbits, while the remaining 5-10% of stars have more eccentric trajectories and belong to the thick disk. We identified seven isolated white dwarfs and two common proper motion pairs as halo members. We determined the age-velocity dispersion relation for the thin-disk members, which agrees with previous results that were achieved from different white dwarf samples without published radial velocities. The age-velocity dispersion relation shows signatures of dynamical heating and saturation after 4-6 Gyr. We observed a m

We identify two ultra-cool ($T_\mathrm{eff} < 4000$ K) metal-polluted (DZ) white dwarfs WDJ2147$-$4035 and WDJ1922$+$0233 as the coolest and second coolest DZ stars known to date with $T_\mathrm{eff} \approx 3050$ K and $T_\mathrm{eff} \approx 3340$ K, respectively. Strong atmospheric collision-induced absorption (CIA) causes the suppression of red optical and infra-red flux in WDJ1922$+$0233, resulting in an unusually blue colour given its low temperature. WDJ2147$-$4035 has moderate infra-red CIA yet has the reddest optical colours known for a DZ white dwarf. Microphysics improvements to the non-ideal effects and CIA opacities in our model atmosphere code yields reasonable solutions to observations of these ultra-cool stars. WDJ2147$-$4035 has a cooling age of over 10 Gyr which is the largest known for a DZ white dwarf, whereas WDJ1922$+$0233 is slightly younger with a cooling age of 9 Gyr. Galactic kinematics calculations from precise Gaia EDR3 astrometry reveal these ultra-cool DZ stars as likely members of the Galactic disc thus they could be pivotal objects in future studies constraining an upper age limit for the disc of the Milky Way. We present intermediate-resolution spectroscopy for both objects, which provides the first spectroscopic observations of WDJ2147$-$4035. Detections of sodium and potassium are made in both white dwarfs, in addition to calcium in WDJ1922$+$0233 and lithium in WDJ2147$-$4035. We identify the magnetic nature of WDJ2147$-$4035 from Zeeman splitting in the lithium line and also make a tentative detection of carbon, so we classify this star as DZQH. WDJ1922$+$0233 likely accreted planetary crust debris, while the debris composition that polluted WDJ2147$-$4035 remains unconstrained., Comment: 18 pages, 8 figures. Accepted by MNRAS

The atmospheres of a large proportion of white dwarf stars are polluted by heavy elements that are expected to sink out of visible layers on short timescales. This has been interpreted as a signature of ongoing accretion of debris from asteroids, comets, and giant planets. This scenario is supported by the detection of debris discs and transits of planetary fragments around some white dwarfs. However, photospheric metals are only indirect evidence for ongoing accretion, and the inferred accretion rates and parent body compositions heavily depend on models of diffusion and mixing processes within the white dwarf atmosphere. Here we report a 4.4$\sigma$ detection of X-rays from a polluted white dwarf, G29$-$38, using a 106 ks exposure with the Chandra X-ray Observatory, demonstrating directly that the star is currently accreting. From the measured X-ray luminosity, we find an instantaneous accretion rate of $\dot{M_{\rm X}}=1.63^{+1.29}_{-0.40}\times 10^{9}\mathrm{\,g\,s^{-1}}$. This is the first direct measurement of the accretion rate onto the white dwarf, which is independent of stellar atmosphere models. This rate exceeds estimates based on past studies of the photospheric abundances by more than a factor two, and implies that convective overshoot has to be accounted for in modelling the spectra of debris-accreting white dwarfs. We measure a low plasma temperature of $kT=0.5\pm0.2\,\mathrm{keV}$, corroborating the predicted bombardment solution for white dwarfs accreting at low accretion rates. Offering a new method for studying evolved planetary systems, these observations provide the opportunity to independently measure the instantaneous accretion rate of planetary material, and therefore investigate the timescale of accretion onto white dwarfs, and the evolution and replenishment of debris disks., Comment: Author submitted version. Published in Nature, Feb 10 2022 issue

Our knowledge of white dwarf planetary systems predominately arises from the region within a few Solar radii of the white dwarfs, where minor planets break up, form rings and discs, and accrete onto the star. The entry location, angle and speed into this Roche sphere has rarely been explored but crucially determines the initial geometry of the debris, accretion rates onto the photosphere, and ultimately the composition of the minor planet. Here we evolve a total of over 10^5 asteroids with single-planet N-body simulations across the giant branch and white dwarf stellar evolution phases to quantify the geometry of asteroid injection into the white dwarf Roche sphere as a function of planetary mass and eccentricity. We find that lower planetary masses increase the extent of anisotropic injection and decrease the probability of head-on (normal to the Roche sphere) encounters. Our results suggest that one can use dynamical activity within the Roche sphere to make inferences about the hidden architectures of these planetary systems., Comment: Accepted for publication in MNRAS

White dwarfs with metal-polluted atmospheres have been studied widely in the context of the accretion of rocky debris from evolved planetary systems. One open question is the geometry of accretion and how material arrives and mixes in the white dwarf surface layers. Using the three-dimensional (3D) radiation hydrodynamics code (COBOLD)-B-5, we present the first transport coefficients in degenerate star atmospheres that describe the advection-diffusion of a passive scalar across the surface plane. We couple newly derived horizontal diffusion coefficients with previously published vertical diffusion coefficients to provide theoretical constraints on surface spreading of metals in white dwarfs. Our grid of 3D simulations probes the vast majority of the parameter space of convective white dwarfs, with pure-hydrogen atmospheres in the effective temperature range of 6000-18 000K and pure-helium atmospheres in the range of 12 000-34 000 K. Our results suggest that warm hydrogen-rich atmospheres (DA; greater than or similar to 13 000 K) and helium-rich atmospheres (DB and DBA; greater than or similar to 30 000 K) are unable to efficiently spread the accreted metals across their surface, regardless of the time dependence of accretion. This result may be at odds with the current non-detection of surface abundance variations in white dwarfs with debris discs. For cooler hydrogen- and helium-rich atmospheres, we predict a largely homogeneous distribution of metals across the surface within a vertical diffusion time-scale. This is typically less than 0.1 per cent of disc lifetime estimates, a quantity that is revisited in this paper using the overshoot results. These results have relevance for studies of the bulk composition of evolved planetary systems and models of accretion disc physics.

Our knowledge of white dwarf planetary systems predominately arises from the region within a few Solar radii of the white dwarfs, where minor planets break up, form rings and discs, and accrete onto the star. The entry location, angle and speed into this Roche sphere has rarely been explored but crucially determines the initial geometry of the debris, accretion rates onto the photosphere, and ultimately the composition of the minor planet. Here we evolve a total of over 10^5 asteroids with single-planet N-body simulations across the giant branch and white dwarf stellar evolution phases to quantify the geometry of asteroid injection into the white dwarf Roche sphere as a function of planetary mass and eccentricity. We find that lower planetary masses increase the extent of anisotropic injection and decrease the probability of head-on (normal to the Roche sphere) encounters. Our results suggest that one can use dynamical activity within the Roche sphere to make inferences about the hidden architectures of these planetary systems., Comment: Accepted for publication in MNRAS

White dwarfs with metal-polluted atmospheres have been studied widely in the context of the accretion of rocky debris from evolved planetary systems. One open question is the geometry of accretion and how material arrives and mixes in the white dwarf surface layers. Using the 3D radiation-hydrodynamics code CO$^5$BOLD, we present the first transport coefficients in degenerate star atmospheres which describe the advection-diffusion of a passive scalar across the surface-plane. We couple newly derived horizontal diffusion coefficients with previously published vertical diffusion coefficients to provide theoretical constraints on surface spreading of metals in white dwarfs. Our grid of 3D simulations probes the vast majority of the parameter space of convective white dwarfs, with pure-hydrogen atmospheres in the effective temperature range 6000-18000 K and pure-helium atmospheres in the range 12000-34000 K. Our results suggest that warm hydrogen-rich atmospheres (DA; $\gtrsim$13000 K) and helium-rich atmospheres (DB, DBA; $\gtrsim$30000 K) are unable to efficiently spread the accreted metals across their surface, regardless of the time dependence of accretion. This result may be at odds with the current non-detection of surface abundance variations at white dwarfs with debris discs. For cooler hydrogen- and helium-rich atmospheres, we predict a largely homogeneous distribution of metals across the surface within a vertical diffusion timescale. This is typically less than 0.1 per cent of disc lifetime estimates, a quantity which is revisited in this paper using the overshoot results. These results have relevance for studies of the bulk composition of evolved planetary systems and models of accretion disc physics., Comment: Accepted for publication in MNRAS. 22 pages, 17 figures

We present an overview of the sample of northern hemisphere white dwarfs within 40 pc of the Sun detected from $Gaia$ Data Release 2 (DR2). We find that 521 sources are spectroscopically confirmed degenerate stars, 111 of which were first identified as white dwarf candidates from $Gaia$ DR2 and followed-up recently with the William Herschel Telescope and Gran Telescopio Canarias. Three additional white dwarf candidates remain spectroscopically unobserved and six unresolved binaries are known to include a white dwarf but were not in our initial selection of white dwarfs in the $Gaia$ DR2 Hertzsprung-Russell diagram (HRD). Atmospheric parameters are calculated from $Gaia$ and Pan-STARRS photometry for all objects in the sample, confirming most of the trends previously observed in the much smaller 20 pc sample. Local white dwarfs are overwhelmingly consistent with Galactic disc kinematics, with only four halo candidates. We find that DAZ white dwarfs are significantly less massive than the overall DA population ($\overline{M}_\mathrm{DAZ} = 0.59\,\mathrm{M}_\odot$, $\overline{M}_\mathrm{DA} = 0.66\,\mathrm{M}_\odot$). It may suggest that planet formation is less efficient at higher mass stars, producing more massive white dwarfs. We detect a sequence of crystallised white dwarfs in the mass range from $0.6\ \lesssim M/\mathrm{M}_\odot \lesssim\ 1.0\,$ and find that the vast majority of objects on the sequence have standard kinematic properties that correspond to the average of the sample, suggesting that their nature can be explained by crystallisation alone. We also detect 56 wide binaries including a white dwarf and 26 double degenerates., Comment: 20 pages, 14 figures. Accepted for publication in MNRAS. Online tables available at http://deneb.astro.warwick.ac.uk/phrgwr/40pcTables/index.html

Ultracompact accreting binary systems each consist of a stellar remnant accreting helium-enriched material from a compact donor star. Such binaries include two related sub-classes, AM CVn-type binaries and helium cataclysmic variables, in both of which the central star is a white dwarf. We present a spectroscopic and photometric study of six accreting binaries with orbital periods in the range of 40--70 min, including phase-resolved VLT spectroscopy and high-speed ULTRACAM photometry. Four of these are AM CVn systems and two are helium cataclysmic variables. For four of these binaries we are able to identify orbital periods (of which three are spectroscopic). SDSS J1505+0659 has an orbital period of 67.8 min, significantly longer than previously believed, and longer than any other known AM CVn binary. We identify a WISE infrared excess in SDSS J1505+0659 that we believe to be the first direct detection of an AM CVn donor star in a non-direct impacting binary. The mass ratio of SDSS J1505+0659 is consistent with a white dwarf donor. CRTS J1028-0819 has an orbital period of 52.1 min, the shortest period of any helium cataclysmic variable. MOA 2010-BLG-087 is co-aligned with a K-class star that dominates its spectrum. ASASSN-14ei and ASASSN-14mv both show a remarkable number of echo outbursts following superoutbursts (13 and 10 echo outbursts respectively). ASASSN-14ei shows an increased outburst rate over the years following its superoutburst, perhaps resulting from an increased accretion rate., Comment: 22 pages, 21 figures. Accepted for publication by MNRAS

We present an overview of the sample of northern hemisphere white dwarfs within 40 pc of the Sun detected from $Gaia$ Data Release 2 (DR2). We find that 521 sources are spectroscopically confirmed degenerate stars, 111 of which were first identified as white dwarf candidates from $Gaia$ DR2 and followed-up recently with the William Herschel Telescope and Gran Telescopio Canarias. Three additional white dwarf candidates remain spectroscopically unobserved and six unresolved binaries are known to include a white dwarf but were not in our initial selection of white dwarfs in the $Gaia$ DR2 Hertzsprung-Russell diagram (HRD). Atmospheric parameters are calculated from $Gaia$ and Pan-STARRS photometry for all objects in the sample, confirming most of the trends previously observed in the much smaller 20 pc sample. Local white dwarfs are overwhelmingly consistent with Galactic disc kinematics, with only four halo candidates. We find that DAZ white dwarfs are significantly less massive than the overall DA population ($\overline{M}_\mathrm{DAZ} = 0.59\,\mathrm{M}_\odot$, $\overline{M}_\mathrm{DA} = 0.66\,\mathrm{M}_\odot$). It may suggest that planet formation is less efficient at higher mass stars, producing more massive white dwarfs. We detect a sequence of crystallised white dwarfs in the mass range from $0.6\ \lesssim M/\mathrm{M}_\odot \lesssim\ 1.0\,$ and find that the vast majority of objects on the sequence have standard kinematic properties that correspond to the average of the sample, suggesting that their nature can be explained by crystallisation alone. We also detect 56 wide binaries including a white dwarf and 26 double degenerates., Comment: 20 pages, 14 figures. Accepted for publication in MNRAS. Online tables available at http://deneb.astro.warwick.ac.uk/phrgwr/40pcTables/index.html

Ultracompact accreting binary systems each consist of a stellar remnant accreting helium-enriched material from a compact donor star. Such binaries include two related sub-classes, AM CVn-type binaries and helium cataclysmic variables, in both of which the central star is a white dwarf. We present a spectroscopic and photometric study of six accreting binaries with orbital periods in the range of 40--70 min, including phase-resolved VLT spectroscopy and high-speed ULTRACAM photometry. Four of these are AM CVn systems and two are helium cataclysmic variables. For four of these binaries we are able to identify orbital periods (of which three are spectroscopic). SDSS J1505+0659 has an orbital period of 67.8 min, significantly longer than previously believed, and longer than any other known AM CVn binary. We identify a WISE infrared excess in SDSS J1505+0659 that we believe to be the first direct detection of an AM CVn donor star in a non-direct impacting binary. The mass ratio of SDSS J1505+0659 is consistent with a white dwarf donor. CRTS J1028-0819 has an orbital period of 52.1 min, the shortest period of any helium cataclysmic variable. MOA 2010-BLG-087 is co-aligned with a K-class star that dominates its spectrum. ASASSN-14ei and ASASSN-14mv both show a remarkable number of echo outbursts following superoutbursts (13 and 10 echo outbursts respectively). ASASSN-14ei shows an increased outburst rate over the years following its superoutburst, perhaps resulting from an increased accretion rate., Comment: 22 pages, 21 figures. Accepted for publication by MNRAS

We present an overview of the sample of Northern hemisphere white dwarfs within 40 pc of the Sun detected from Gaia Data Release 2 (DR2). We find that 521 sources are spectroscopically confirmed degenerate stars, 111 of which were first identified as white dwarf candidates from Gaia DR2 and followed up recently with the William Herschel Telescope and Gran Telescopio Canarias. Three additional white dwarf candidates remain spectroscopically unobserved and six unresolved binaries are known to include a white dwarf but were not in our initial selection in the Gaia DR2 Hertzsprung-Russell diagram. Atmospheric parameters are calculated from Gaia and Pan-STARRS photometry for all objects in the sample, confirming most of the trends previously observed in the much smaller 20 pc sample. Local white dwarfs are overwhelmingly consistent with Galactic disc kinematics, with only four halo candidates. We find that DAZ white dwarfs are significantly less massive than the overall DA population (M¯DAZ = 0.59 M⊙, M¯DA = 0.66 M⊙). It may suggest that planet formation is less efficient at higher mass stars, producing more massive white dwarfs. We detect a sequence of crystallized white dwarfs in the mass range from 0.6 ≤M/M⊙ ≤ 1.0 and find that the vast majority of objects on the sequence have standard kinematic properties that correspond to the average of the sample, suggesting that their nature can be explained by crystallization alone. We also detect 26 double degenerates and white dwarf components in 56 wide binary systems.

We present a theoretical description of macroscopic diffusion caused by convective overshoot in pure-hydrogen DA white dwarfs using 3D, closed-bottom, radiation hydrodynamics (COBOLD)-B-5 simulations. We rely on a new grid of deep 3D white dwarf models in the temperature range 11 400 <= T-eff <= 18 000 K where tracer particles and a tracer density are used to derive macroscopic diffusion coefficients driven by convective overshoot. These diffusion coefficients are compared to microscopic diffusion coefficients from 1D structures. We find that the mass of the fully mixed region is likely to increase by up to 2.5 orders of magnitude while inferred accretion rates increase by a more moderate order of magnitude. We present evidence that an increase in settling time of up to 2 orders of magnitude is to be expected, which is of significance for time-variability studies of polluted white dwarfs. Our grid also provides the most robust constraint on the onset of convective instabilities in DA white dwarfs to be in the effective temperature range from 18 000 to 18 250 K.

We present a theoretical description of macroscopic diffusion caused by convective overshoot in pure-hydrogen DA white dwarfs using 3D, closed-bottom, radiation hydrodynamics (COBOLD)-B-5 simulations. We rely on a new grid of deep 3D white dwarf models in the temperature range 11 400 <= T-eff <= 18 000 K where tracer particles and a tracer density are used to derive macroscopic diffusion coefficients driven by convective overshoot. These diffusion coefficients are compared to microscopic diffusion coefficients from 1D structures. We find that the mass of the fully mixed region is likely to increase by up to 2.5 orders of magnitude while inferred accretion rates increase by a more moderate order of magnitude. We present evidence that an increase in settling time of up to 2 orders of magnitude is to be expected, which is of significance for time-variability studies of polluted white dwarfs. Our grid also provides the most robust constraint on the onset of convective instabilities in DA white dwarfs to be in the effective temperature range from 18 000 to 18 250 K.

We present a theoretical description of macroscopic diffusion caused by convective overshoot in pure-hydrogen DA white dwarfs using 3D, closed-bottom, radiation hydrodynamics (COBOLD)-B-5 simulations. We rely on a new grid of deep 3D white dwarf models in the temperature range 11 400 <= T-eff <= 18 000 K where tracer particles and a tracer density are used to derive macroscopic diffusion coefficients driven by convective overshoot. These diffusion coefficients are compared to microscopic diffusion coefficients from 1D structures. We find that the mass of the fully mixed region is likely to increase by up to 2.5 orders of magnitude while inferred accretion rates increase by a more moderate order of magnitude. We present evidence that an increase in settling time of up to 2 orders of magnitude is to be expected, which is of significance for time-variability studies of polluted white dwarfs. Our grid also provides the most robust constraint on the onset of convective instabilities in DA white dwarfs to be in the effective temperature range from 18 000 to 18 250 K.

We present a study of the hypothesis that white dwarfs undergo a spectral change from hydrogen- to helium-dominated atmospheres using a volume-limited photometric sample drawn from the $Gaia$ DR2 catalogue, the Sloan Digital Sky Survey (SDSS) and the $Galaxy~Evolution~Explorer$ (GALEX). We exploit the strength of the Balmer jump in hydrogen-atmosphere DA white dwarfs to separate them from helium-dominated objects in SDSS colour space. Across the effective temperature range from 20000K to 9000K we find that 22% of white dwarfs will undergo a spectral change, with no spectral evolution being ruled out at 5$\sigma$. The most likely explanation is that the increase in He-rich objects is caused by the convective mixing of DA stars with thin hydrogen layers, in which helium is dredged up from deeper layers by a surface hydrogen convection zone. The rate of change in the fraction of He-rich objects as a function of temperature, coupled with a recent grid of 3D radiation-hydrodynamic simulations of convective DA white dwarfs - which include the full overshoot region - lead to a discussion on the distribution of total hydrogen mass in white dwarfs. We find that 60% of white dwarfs must have a hydrogen mass larger than $M_{\rm H}/M_{\rm WD} = 10^{-10}$, another 25% have masses in the range $M_{\rm H}/M_{\rm WD} = 10^{-14}-10^{-10}$, and 15% have less hydrogen than $M_{\rm H}/M_{\rm WD} = 10^{-14}$. These results have implications for white dwarf asteroseismology, stellar evolution through the asymptotic giant branch (AGB) and accretion of planetesimals onto white dwarfs., Comment: 14 pages, 13 figures, accepted for publication in MNRAS

White dwarfs are stellar embers depleted of nuclear energy sources that cool over billions of years. These stars, which are supported by electron degeneracy pressure, reach densities of 1e7 grams per cubic centimetre in their cores. It has been predicted that a first-order phase transition occurs during white-dwarf cooling, leading to the crystallization of the non-degenerate carbon and oxygen ions in the core, which releases a considerable amount of latent heat and delays the cooling process by about one billion years. However, no direct observational evidence of this effect has been reported so far. Here we report the presence of a pile-up in the cooling sequence of evolving white dwarfs within 100 parsecs of the Sun, determined using photometry and parallax data from the Gaia satellite. Using modelling, we infer that this pile-up arises from the release of latent heat as the cores of the white dwarfs crystallize. In addition to the release of latent heat, we find strong evidence that cooling is further slowed by the liberation of gravitational energy from element sedimentation in the crystallizing cores. Our results describe the energy released by crystallization in strongly coupled Coulomb plasmas, and the measured cooling delays could help to improve the accuracy of methods used to determine the age of stellar populations from white dwarfs., Comment: 17 pages, author accepted version. Published in Nature 9 January 2019

We present a theoretical description of macroscopic diffusion caused by convective overshoot in pure-hydrogen DA white dwarfs using three-dimensional (3D), closed-bottom, radiation hydrodynamics CO$^5$BOLD simulations. We rely on a new grid of deep 3D white dwarf models in the temperature range 11400 K $\leq T_{\mathrm{eff}} \leq$ 18000 K where tracer particles and a tracer density are used to derive macroscopic diffusion coefficients driven by convective overshoot. These diffusion coefficients are compared to microscopic diffusion coefficients from one-dimensional structures. We find that the mass of the fully mixed region is likely to increase by up to 2.5 orders of magnitude while inferred accretion rates increase by a more moderate order of magnitude. We present evidence that an increase in settling time of up to 2 orders of magnitude is to be expected which is of significance for time-variability studies of polluted white dwarfs. Our grid also provides the most robust constraint on the onset of convective instabilities in DA white dwarfs to be in the effective temperature range from 18000 to 18250 K., Comment: 21 pages, 24 figures, accepted for publication in MNRAS

Practically all known planet hosts will evolve into white dwarfs, and large parts of their planetary systems will survive this transition - the same is true for the solar system beyond the orbit of Mars. Spectroscopy of white dwarfs accreting planetary debris provides the most accurate insight into the bulk composition of exo-planets. Ground-based spectroscopic surveys of ~260, 000 white dwarfs detected with Gaia will identify >1000 evolved planetary systems, and high-throughput high-resolution space-based ultraviolet spectroscopy is essential to measure in detail their abundances. So far, evidence for two planetesimals orbiting closely around white dwarfs has been obtained, and their study provides important constraints on the composition and internal structure of these bodies. Major photometric and spectroscopic efforts will be necessary to assemble a sample of such close-in planetesimals that is sufficiently large to establish their properties as a population, and to deduce the architectures of the outer planetary systems from where they originated. Mid-infrared spectroscopy of the dusty disks will provide detailed mineralogical information of the debris, which, in combination with the elemental abundances measured from the white dwarf spectroscopy, will enable detailed physical modelling of the chemical, thermodynamic, and physical history of the accreted material. Flexible multi-epoch infrared observations are essential to determine the physical nature, and origin of the variability observed in many of the dusty disks. Finally, the direct detection of the outer reservoirs feeding material to the white dwarfs will require sensitive mid- and far-infrared capabilities., Comment: 5 pages, 3 figures Science White Paper submitted to the Astro2020 Decadal Survey

We present a theoretical description of macroscopic diffusion caused by convective overshoot in pure-hydrogen DA white dwarfs using 3D, closed-bottom, radiation hydrodynamics (COBOLD)-B-5 simulations. We rely on a new grid of deep 3D white dwarf models in the temperature range 11 400 <= T-eff <= 18 000 K where tracer particles and a tracer density are used to derive macroscopic diffusion coefficients driven by convective overshoot. These diffusion coefficients are compared to microscopic diffusion coefficients from 1D structures. We find that the mass of the fully mixed region is likely to increase by up to 2.5 orders of magnitude while inferred accretion rates increase by a more moderate order of magnitude. We present evidence that an increase in settling time of up to 2 orders of magnitude is to be expected, which is of significance for time-variability studies of polluted white dwarfs. Our grid also provides the most robust constraint on the onset of convective instabilities in DA white dwarfs to be in the effective temperature range from 18 000 to 18 250 K.

We present a theoretical description of macroscopic diffusion caused by convective overshoot in pure-hydrogen DA white dwarfs using three-dimensional (3D), closed-bottom, radiation hydrodynamics CO$^5$BOLD simulations. We rely on a new grid of deep 3D white dwarf models in the temperature range 11400 K $\leq T_{\mathrm{eff}} \leq$ 18000 K where tracer particles and a tracer density are used to derive macroscopic diffusion coefficients driven by convective overshoot. These diffusion coefficients are compared to microscopic diffusion coefficients from one-dimensional structures. We find that the mass of the fully mixed region is likely to increase by up to 2.5 orders of magnitude while inferred accretion rates increase by a more moderate order of magnitude. We present evidence that an increase in settling time of up to 2 orders of magnitude is to be expected which is of significance for time-variability studies of polluted white dwarfs. Our grid also provides the most robust constraint on the onset of convective instabilities in DA white dwarfs to be in the effective temperature range from 18000 to 18250 K., Comment: 21 pages, 24 figures, accepted for publication in MNRAS

Practically all known planet hosts will evolve into white dwarfs, and large parts of their planetary systems will survive this transition - the same is true for the solar system beyond the orbit of Mars. Spectroscopy of white dwarfs accreting planetary debris provides the most accurate insight into the bulk composition of exo-planets. Ground-based spectroscopic surveys of ~260, 000 white dwarfs detected with Gaia will identify >1000 evolved planetary systems, and high-throughput high-resolution space-based ultraviolet spectroscopy is essential to measure in detail their abundances. So far, evidence for two planetesimals orbiting closely around white dwarfs has been obtained, and their study provides important constraints on the composition and internal structure of these bodies. Major photometric and spectroscopic efforts will be necessary to assemble a sample of such close-in planetesimals that is sufficiently large to establish their properties as a population, and to deduce the architectures of the outer planetary systems from where they originated. Mid-infrared spectroscopy of the dusty disks will provide detailed mineralogical information of the debris, which, in combination with the elemental abundances measured from the white dwarf spectroscopy, will enable detailed physical modelling of the chemical, thermodynamic, and physical history of the accreted material. Flexible multi-epoch infrared observations are essential to determine the physical nature, and origin of the variability observed in many of the dusty disks. Finally, the direct detection of the outer reservoirs feeding material to the white dwarfs will require sensitive mid- and far-infrared capabilities., Comment: 5 pages, 3 figures Science White Paper submitted to the Astro2020 Decadal Survey

We present a study of the hypothesis that white dwarfs undergo a spectral change from hydrogen- to helium-dominated atmospheres using a volume-limited photometric sample drawn from the $Gaia$ DR2 catalogue, the Sloan Digital Sky Survey (SDSS) and the $Galaxy~Evolution~Explorer$ (GALEX). We exploit the strength of the Balmer jump in hydrogen-atmosphere DA white dwarfs to separate them from helium-dominated objects in SDSS colour space. Across the effective temperature range from 20000K to 9000K we find that 22% of white dwarfs will undergo a spectral change, with no spectral evolution being ruled out at 5$\sigma$. The most likely explanation is that the increase in He-rich objects is caused by the convective mixing of DA stars with thin hydrogen layers, in which helium is dredged up from deeper layers by a surface hydrogen convection zone. The rate of change in the fraction of He-rich objects as a function of temperature, coupled with a recent grid of 3D radiation-hydrodynamic simulations of convective DA white dwarfs - which include the full overshoot region - lead to a discussion on the distribution of total hydrogen mass in white dwarfs. We find that 60% of white dwarfs must have a hydrogen mass larger than $M_{\rm H}/M_{\rm WD} = 10^{-10}$, another 25% have masses in the range $M_{\rm H}/M_{\rm WD} = 10^{-14}-10^{-10}$, and 15% have less hydrogen than $M_{\rm H}/M_{\rm WD} = 10^{-14}$. These results have implications for white dwarf asteroseismology, stellar evolution through the asymptotic giant branch (AGB) and accretion of planetesimals onto white dwarfs., Comment: 14 pages, 13 figures, accepted for publication in MNRAS

White dwarfs are stellar embers depleted of nuclear energy sources that cool over billions of years. These stars, which are supported by electron degeneracy pressure, reach densities of 1e7 grams per cubic centimetre in their cores. It has been predicted that a first-order phase transition occurs during white-dwarf cooling, leading to the crystallization of the non-degenerate carbon and oxygen ions in the core, which releases a considerable amount of latent heat and delays the cooling process by about one billion years. However, no direct observational evidence of this effect has been reported so far. Here we report the presence of a pile-up in the cooling sequence of evolving white dwarfs within 100 parsecs of the Sun, determined using photometry and parallax data from the Gaia satellite. Using modelling, we infer that this pile-up arises from the release of latent heat as the cores of the white dwarfs crystallize. In addition to the release of latent heat, we find strong evidence that cooling is further slowed by the liberation of gravitational energy from element sedimentation in the crystallizing cores. Our results describe the energy released by crystallization in strongly coupled Coulomb plasmas, and the measured cooling delays could help to improve the accuracy of methods used to determine the age of stellar populations from white dwarfs., Comment: 17 pages, author accepted version. Published in Nature 9 January 2019

We present a theoretical description of macroscopic diffusion caused by convective overshoot in pure-hydrogen DA white dwarfs using 3D, closed-bottom, radiation hydrodynamics (COBOLD)-B-5 simulations. We rely on a new grid of deep 3D white dwarf models in the temperature range 11 400 <= T-eff <= 18 000 K where tracer particles and a tracer density are used to derive macroscopic diffusion coefficients driven by convective overshoot. These diffusion coefficients are compared to microscopic diffusion coefficients from 1D structures. We find that the mass of the fully mixed region is likely to increase by up to 2.5 orders of magnitude while inferred accretion rates increase by a more moderate order of magnitude. We present evidence that an increase in settling time of up to 2 orders of magnitude is to be expected, which is of significance for time-variability studies of polluted white dwarfs. Our grid also provides the most robust constraint on the onset of convective instabilities in DA white dwarfs to be in the effective temperature range from 18 000 to 18 250 K.

We present a catalogue of white dwarf candidates selected from the second data release of $Gaia$ (DR2). We used a sample of spectroscopically confirmed white dwarfs from the Sloan Digital Sky Survey (SDSS) to map the entire space spanned by these objects in the $Gaia$ Hertzsprung-Russell diagram. We then defined a set of cuts in absolute magnitude, colour, and a number of $Gaia$ quality flags to remove the majority of contaminating objects. Finally, we adopt a method analogous to the one presented in our earlier SDSS photometric catalogues to calculate a probability of being a white dwarf ($P_{\mathrm{WD}}$) for all $Gaia$ sources which passed the initial selection. The final catalogue is composed of $486\,641$ stars with calculated $P_{\mathrm{WD}}$ from which it is possible to select a sample of $\simeq 260\,000$ high-confidence white dwarf candidates in the magnitude range $8 7000$ K, at high Galactic latitudes ($|b|>20^{\circ}$). However, the completeness drops at low Galactic latitudes, and the magnitude limit of the catalogue varies significantly across the sky as a function of $Gaia$'s scanning law. We also provide the list of objects within our sample with available SDSS spectroscopy. We use this spectroscopic sample to characterise the observed structure of the white dwarf distribution in the H-R diagram., Comment: 23 pages, 21 figures, 4 tables. Accepted for publication in MNRAS. Version 3 accepted. Gaia DR2 catalogue of white dwarf candidates available . SDSS spectra sample available

We present a catalogue of white dwarf candidates selected from the second data release of $Gaia$ (DR2). We used a sample of spectroscopically confirmed white dwarfs from the Sloan Digital Sky Survey (SDSS) to map the entire space spanned by these objects in the $Gaia$ Hertzsprung-Russell diagram. We then defined a set of cuts in absolute magnitude, colour, and a number of $Gaia$ quality flags to remove the majority of contaminating objects. Finally, we adopt a method analogous to the one presented in our earlier SDSS photometric catalogues to calculate a probability of being a white dwarf ($P_{\mathrm{WD}}$) for all $Gaia$ sources which passed the initial selection. The final catalogue is composed of $486\,641$ stars with calculated $P_{\mathrm{WD}}$ from which it is possible to select a sample of $\simeq 260\,000$ high-confidence white dwarf candidates in the magnitude range $8 7000$ K, at high Galactic latitudes ($|b|>20^{\circ}$). However, the completeness drops at low Galactic latitudes, and the magnitude limit of the catalogue varies significantly across the sky as a function of $Gaia$'s scanning law. We also provide the list of objects within our sample with available SDSS spectroscopy. We use this spectroscopic sample to characterise the observed structure of the white dwarf distribution in the H-R diagram., Comment: 23 pages, 21 figures, 4 tables. Accepted for publication in MNRAS. Version 3 accepted. Gaia DR2 catalogue of white dwarf candidates available . SDSS spectra sample available

The Palace and Landscape Project at Palaikastro excavations, conducted under the auspices of the British School at Athens from 2013 to 2015, uncovered remains of three buildings in a heretofore unexplored area of the site. These buildings, two of which were completely revealed, show significant differences in construction, plan, and situation in what is essentially a different neighborhood of the Bronze Age town. With preserved floor deposits from Late Minoan IB and Late Minoan IIIB and contexts going back to Middle Minoan II, these newly revealed structures advance our understanding of the layout of the site and have implications for social organisation and chronological developments at Palaikastro and elsewhere on Crete. At the settlement level, we discuss the integration of this new sector in the overall layout of the town. At the neighbourhood scale, this paper focuses on a description of the three buildings, both in terms of architectural techniques and spatial organization, while drawing comparisons with other buildings from the Roussolakkos area. Finally, we address the functional development of the sector which witnessed important changes when, in Postpalatial times, the domestic na-ture of the buildings seemed to have dwindled in favor of a clear burst of ritual activities, attested both by various deposits (e.g., cultic material, feasting refuse) and architectural alterations of the preexisting built environment.

This volume comprises the proceedings of a workshop with the same title that took place in October 2015.

Presentation on system producibility, DARPA/USAF Wright Laboratory (ManTech)

info:eu-repo/semantics/published