Search for: All records

Abstract The bulk abundances of exoplanetesimals can be measured when they are accreted by white dwarfs. Recently, lithium from the accretion of exoplanetesimals was detected in relatively high levels in multiple white dwarfs. There are presently three proposed hypotheses to explain the detection of excess lithium in white dwarf photospheres: Big Bang and Galactic nucleosynthesis, continental crust, and an exomoon formed from spalled ring material. We present new observations of three previously known lithium-polluted white dwarfs (WD J1824+1213, WD J2317+1830, and LHS 2534), and one with metal pollution without lithium (SDSS J1636+1619). We also present atmospheric model fits to these white dwarfs. We then evaluate the abundances of these white dwarfs and two additional lithium-polluted white dwarfs that were previously fit using the same atmospheric models (WD J1644-0449 and SDSS J1330+6435) in the context of the three extant hypotheses for explaining lithium excesses in polluted white dwarfs. We find Big Bang and Galactic nucleosynthesis to be the most plausible explanation of the abundances in WD J1644-0449, WD J1824+1213, and WD J2317+1830. SDSS J1330+6435 will require stricter abundances to determine its planetesimal’s origins, and LHS 2534, as presently modeled, defies all three hypotheses. We find the accretion of an exomoon formed from spalled ring material to be highly unlikely to be the explanation of the lithium excess in any of these cases. 

ABSTRACT A few per cent of red giants are enriched in lithium with $$A(\mathrm{Li}) \gt 1.5$$. Their evolutionary status has remained uncertain because these Li-rich giants can be placed both on the red giant branch (RGB) near the bump luminosity and in the red clump (RC) region. However, thanks to asteroseismology, it has been found that most of them are actually RC stars. Starting at the bump luminosity, RGB progenitors of the RC stars experience extra mixing in the radiative zone separating the H-burning shell from the convective envelope followed by a series of convective He-shell flashes at the RGB tip, known as the He-core flash. The He-core flash was proposed to cause fast extra mixing in the stars at the RGB tip that is needed for the Cameron–Fowler mechanism to produce Li. We propose that the RGB stars are getting enriched in Li by the RGB extra mixing that is getting enhanced and begins to produce Li, instead of destroying it, when the stars are approaching the RGB tip. After a discussion of several mechanisms of the RGB extra mixing, including the joint operation of rotation-driven meridional circulation and turbulent diffusion, the azimuthal magnetorotational instability (AMRI), thermohaline convection, buoyancy of magnetic flux tubes, and internal gravity waves, and based on results of (magneto-) hydrodynamics simulations and asteroseismology observations, we are inclined to conclude that it is the mechanism of the AMRI or magnetically enhanced thermohaline convection, that is most likely to support our hypothesis. 

Abstract Infrared-faint white dwarfs are cool white dwarfs exhibiting significant infrared flux deficits, most often attributed to collision-induced absorption (CIA) from H₂–He in mixed hydrogen–helium atmospheres. We present James Webb Space Telescope (JWST) near- and mid-infrared spectra of three such objects using Near-Infrared Spectrograph (0.6–5.3μm) and Mid-Infrared Instrument (5–14μm): LHS 3250, WD J1922+0233, and LHS 1126. Surprisingly, for LHS 3250, we detect no H₂–He CIA absorption at 2.4μm, instead observing an unexpected small flux bump at this wavelength. WD J1922+0233 exhibits the anticipated strong absorption feature centered at 2.4μm, but with an unexpected narrow emission-like feature inside this absorption band. LHS 1126 shows no CIA features and follows aλ⁻²power law in the mid-infrared. LHS 1126's lack of CIA features suggests a very low hydrogen abundance, with its infrared flux depletion likely caused by He–He–He CIA. For LHS 3250 and WD J1922+0233, the absence of a 1.2μm CIA feature in both stars argues against ultracool temperatures, supporting recent suggestions that infrared-faint (IR-faint) white dwarfs are warmer and more massive than previously thought. This conclusion is further solidified by Keck near-infrared spectroscopy of seven additional objects. We explore possible explanations for the unexpected emission-like features in both stars, and temperature inversions above the photosphere emerge as a promising hypothesis. Such inversions may be common among the IR-faint population, and since they significantly affect the infrared spectral energy distribution, this would impact their photometric fits. Further JWST observations are needed to confirm the prevalence of this phenomenon and guide the development of improved atmospheric models. 

Abstract We increase the spectroscopic completeness of the 100 pc white dwarf sample in the Sloan Digital Sky Survey footprint with 840 additional spectra. Our spectroscopy is 86% complete for white dwarfs hotter thanT_eff = 5000 K, where Hαremains visible and provides reliable constraints on the atmospheric composition. We identify 2108 DA white dwarfs with pure hydrogen atmospheres, and show that ultramassive DA white dwarfs withM≥ 1.1M_⊙are an order of magnitude less common below 10,000 K. This is consistent with a fraction of them getting stuck on the crystallization sequence due to²²Ne distillation. In addition, there are no ultramassive DA white dwarfs withM≥ 1.1M_⊙andT_eff≤ 6000 K in our sample, likely because Debye cooling makes them rapidly fade away. We detect a significant trend in the fraction of He atmosphere white dwarfs as a function of temperature; the fraction increases from 9% at 20,000 K to 32% at 6000 K. This provides direct evidence of convective mixing in cool DA white dwarfs. Finally, we detect a relatively tight sequence of low-mass DQ white dwarfs in color–magnitude diagrams for the first time. We discuss the implications of this tight DQ sequence, and conclude with a discussion of the future prospects from the upcoming Ultraviolet Transient Astronomy Satellite mission and the large-scale multi-fiber spectroscopic surveys. 

Abstract We present 3D hydrodynamical simulations of core convection with a stably stratified envelope of a 25M_⊙star in the early phase of the main sequence. We use the explicit gas-dynamics codePPMstar, which tracks two fluids and includes radiation pressure and radiative diffusion. Multiple series of simulations with different luminosities and radiative thermal conductivities are presented. The entrainment rate at the convective boundary, internal gravity waves in and above the boundary region, and the approach to dynamical equilibrium shortly after a few convective turnovers are investigated. We perform very long simulations on 896³grids accelerated by luminosity boost factors of 1000, 3162 and 10,000. In these simulations, the growing penetrative convection reduces the initially unrealistically large entrainment. This reduction is enabled by a spatial separation that develops between the entropy gradient and the composition gradient. The convective boundary moves outward much more slowly at the end of these simulations. Finally, we present a 1D method to predict the extent and character of penetrative convection beyond the Schwarzschild boundary. The 1D model is based on a spherically averaged reduced entropy equation that takes the turbulent dissipation as input from the 3D hydrodynamic simulation and takes buoyancy and all other energy sources and sinks into account. This 1D method is intended to be ultimately deployed in 1D stellar evolution calculations and is based on the properties of penetrative convection in our simulations carried forward through the local thermal timescale. 

Abstract We present follow-up spectroscopy and a detailed model atmosphere analysis of 29 wide double white dwarfs, including eight systems with a crystallized C/O core member. We use the state-of-the-art evolutionary models to constrain the physical parameters of each star, including the total age. Assuming that the members of wide binaries are coeval, any age difference between the binary members can be used to test the cooling physics for white dwarf stars, including potential delays due to crystallization and²²Ne distillation. We use our control sample of 14 wide binaries with noncrystallized members to show that this method works well; the control sample shows an age difference of only ΔAge = −0.03 ± 0.15 Gyr between its members. For the eight crystallized C/O core systems we find a cooling anomaly of ΔAge = ${1.13}_{-1.07}^{+1.20}$Gyr. Even though our results are consistent with a small additional cooling delay (∼1 Gyr) from²²Ne distillation and other neutron-rich impurities, the large uncertainties make this result not statistically significant. Nevertheless, we rule out cooling delays longer than 3.6 Gyr at the 99.7% (3σ) confidence level for 0.6–0.9M_⊙white dwarfs. Further progress requires larger samples of wide binaries with crystallized massive white dwarf members. We provide a list of subgiant + white dwarf binaries that could be used for this purpose in the future. 

Abstract We present a detailed model atmosphere analysis of massive white dwarfs withM> 0.9M_⊙andT_eff≥ 11,000 K in the Montreal White Dwarf Database 100 pc sample and the Pan-STARRS footprint. We obtained follow-up optical spectroscopy of 109 objects with no previous spectral classification in the literature. Our spectroscopic follow-up is now complete for all 204 objects in the sample. We find 118 normal DA white dwarfs, including 45 massive DAs near the ZZ Ceti instability strip. There are no normal massive DBs: the six DBs in the sample are strongly magnetic and/or rapidly rotating. There are 20 massive DQ white dwarfs in our sample, and all are found in the crystallization sequence. In addition, 66 targets are magnetic (32% of the sample). We use magnetic white dwarf atmosphere models to constrain the field strength and geometry using offset dipole models. We also use magnetism, kinematics, and rotation measurements to constrain the fraction of merger remnant candidates among this population. The merger fraction of this sample increases from 25% for 0.9–1M_⊙white dwarfs to 49% for 1.2–1.3M_⊙. However, this fraction is as high as ${78}_{-7}^{+4}$% for 1.1–1.2M_⊙white dwarfs. Previous works have demonstrated that 5%–9% of high-mass white dwarfs stop cooling for ∼8 Gyr due to the²²Ne distillation process, which leads to an overdensity of Q-branch stars in the solar neighborhood. We demonstrate that the overabundance of the merger remnant candidates in our sample is likely due to the same process. 

ABSTRACT Recent photometric observations of massive stars have identified a low-frequency power excess which appears as stochastic low-frequency variability in light-curve observations. We present the oscillation properties of high-resolution hydrodynamic simulations of a $$25\,\,{\rm{M}_\odot }$$ star performed with the PPMstar code. The model star has a convective core mass of $$\approx 12\,\,{\rm{M}_\odot }$$ and approximately half of the envelope simulated. From this simulation, we extract light curves from several directions, average them over each hemisphere, and process them as if they were real photometric observations. We show how core convection excites waves with a similar frequency as the convective time-scale in addition to significant power across a forest of low and high angular degree l modes. We find that the coherence of these modes is relatively low as a result of their stochastic excitation by core convection, with lifetimes of the order of 10s of days. Thanks to the still significant power at higher l and this relatively low coherence, we find that integrating over a hemisphere produces a power spectrum that still contains measurable power up to the Brunt–Väisälä frequency. These power spectra extracted from the stable envelope are qualitatively similar to observations, with the same order of magnitude yet lower characteristic frequency. This work further shows the potential of long-duration, high-resolution hydrodynamic simulations for connecting asteroseismic observations to the structure and dynamics of core convection and the convective boundary. 

Abstract Four years after the discovery of a unique DAQ white dwarf with a hydrogen-dominated and carbon-rich atmosphere, we report the discovery of four new DAQ white dwarfs, including two that were not recognized properly in the literature. We find all five DAQs in a relatively narrow mass and temperature range ofM= 1.14–1.19M_⊙andT_eff= 13,000–17,000 K. In addition, at least two show photometric variations due to rapid rotation with ≈10 minute periods. All five are also kinematically old, but appear photometrically young, with estimated cooling ages of about 1 Gyr based on standard cooling tracks, and their masses are roughly twice the mass of the most common white dwarfs in the solar neighborhood. These characteristics are smoking gun signatures of white dwarf merger remnants. Comparing the DAQ sample with warm DQ white dwarfs, we demonstrate that there is a range of hydrogen abundances among the warm DQ population and that the distinction between DAQ and warm DQ white dwarfs is superficial. We discuss the potential evolutionary channels for the emergence of the DAQ subclass, suggesting that DAQ white dwarfs are trapped on the crystallization sequence and may remain there for a significant fraction of the Hubble time. 

ABSTRACT The inner structure of core helium burning (CHeB) stars remains uncertain due to the yet unknown nature of mixing at the boundary of their cores. Large convective cores beyond a bare Schwarzschild model are favoured both from theoretical arguments and from asteroseismological constraints. However, the exact nature of this extra mixing, and in particular the possible presence of semiconvective layers, is still debated. In this work, we approach this problem through a new avenue by performing the first full-sphere 3D hydrodynamics simulations of the interiors of CHeB stars. We use the ppmstar explicit gas dynamics code to simulate the inner 0.45$$\, {\rm M}_{\odot }$$ of a 3 M⊙ CHeB star. Simulations are performed using different Cartesian grid resolutions (7683, 11523, and 17283) and heating rates. We use two different initial states, one based on mesas's predictive mixing scheme (which significantly extends the core beyond the Schwarzschild boundary) and one based on the convective premixing approach (which exhibits a semiconvective interface). The general behaviour of the flow in the convective core and in the stable envelope (where internal gravity waves are observed) is consistent with our recent simulations of core convection in massive main-sequence stars, and so are the various luminosity scaling relations. The semiconvective layers are dominated by strong internal gravity waves that do not produce measurable species mixing, but overshooting motions from the convective core gradually homogenize the semiconvective interface. This process can possibly completely erase the semiconvective layers, which would imply that CHeB stars do not harbour a semiconvection zone.