We analyze a large brightening event that lasted for roughly half a day in the pulsating hydrogen-atmosphere white dwarf GD 1212 during K2 Campaign 12 of the extended Kepler mission. For the other 80 days of K2 observations, GD 1212 exhibited a rich spectrum of long-period ( 1100 s) pulsations that underwent rapid variations in frequency and amplitude but did not exhibit any additional outbursts. We refine previous attempts at mode identification and find a likely sequence of dipole and quadrupole splittings that reveal an overall rotation rate of roughly 17.0 hr. The outburst at Day 61 is fully resolved by the 60-second-cadence K2 data, with the entire white dwarf becoming up to 17.5
Theoretical spectral line shapes describe the distribution of opacity due to bound electronic transitions in hot dense plasmas and are used to fit emergent spectra from many astrophysical sources. Deficiencies in line broadening theory have been proposed as a possible explanation for unresolved discrepancies between theoretical and observed spectra in white dwarf star atmospheres and laboratory experiments at white dwarf star photosphere conditions. One possible source of these discrepancies is the formation of quasi-molecules. Quasi-molecules are close (unbound) collisions between atoms, which broaden line shapes and create additional satellite lines. Quasi-molecules are challenging to implement into traditional line shape codes and have historically required a number of physical approximations beyond what is used in standard Stark broadening models. Here we present a new approach to calculating line shapes with quasi-molecular resonances, using a novel multiple-basis method that considers both atomic and molecular states. We implement this approach into a simulation line shape code, present hydrogen Lyman-series line shapes with quasi-H_2^+ resonances, and demonstrate the impact our new line shapes have on hydrogen-atmosphere white dwarf star model spectra. We find that our new approach leads to broader quasi-molecular features that agree well with observed spectra in initial comparisons.
The structure and emergent flux of hydrogen atmosphere white dwarfs depend on the opacity of the Ly α and Ly β spectral lines. The opacity here is set by the strength and broadening of these lines; the latter is dictated by the far line wing, which is in the “quasi-static” limit of electron broadening, placing it in the incomplete collision regime, and describes the transient parts of electron and ion collisions. These transient stages of the collision form resonances: In the case of ions, they manifest as molecular resonances, while for electrons they are H ^− resonances, both of which can only be captured quantum-mechanically. Quantum-mechanical calculations have historically preserved only a handful of broadening terms that are most important near the center of the line. However, in the wings of the line, the previously neglected terms that describe the transient stages of the collision need to be included. This requirement arises because, in the line wings, the broadening from the 1 s ground state, which is generally assumed to be extremely small compared to the broadening of the upper state, is no longer negligible within a quantum-mechanical model that takes into account exchange interactions. The inclusion of all the transient terms results in asymmetries and extra broadening. The increased broadening of Ly α increases the opacity at the energy where most of the flux leaves the star. The broader Ly α lines also impact the visible flux, raising it by an amount that exceeds previously estimated errors.
White dwarfs (WDs) polluted by exoplanetary material provide the unprecedented opportunity to directly observe the interiors of exoplanets. However, spectroscopic surveys are often limited by brightness constraints, and WDs tend to be very faint, making detections of large populations of polluted WDs difficult. In this paper, we aim to increase considerably the number of WDs with multiple metals in their atmospheres. Using 96,134 WDs with Gaia DR3 BP/RP (XP) spectra, we constructed a 2D map using an unsupervised machine-learning technique called Uniform Manifold Approximation and Projection (UMAP) to organize the WDs into identifiable spectral regions. The polluted WDs are among the distinct spectral groups identified in our map. We have shown that this selection method could potentially increase the number of known WDs with five or more metal species in their atmospheres by an order of magnitude. Such systems are essential for characterizing exoplanet diversity and geology.
Accurate models of cooling white dwarfs must treat the energy released as their cores crystallize. This phase transition slows the cooling by releasing latent heat and also gravitational energy, which results from phase separation: liquid C is released from the solid C/O core, driving an outward carbon flux. The Gaia color-magnitude diagram provides striking confirmation of this theory by revealing a mass-dependent overdensity of white dwarfs, indicating slowed cooling at the expected location. However, the observed overdensity is enhanced relative to the models. Additionally, it is associated with increased magnetism, suggesting a link between crystallization and magnetic field generation. Recent works aimed at explaining an enhanced cooling delay and magnetic field generation employ a uniform mixing prescription that assumes large-scale turbulent motions; we show here that these calculations are not self-consistent. We also show that thermohaline mixing is most likely efficient enough to provide the required chemical redistribution during C/O phase separation, and that the resulting velocities and mixing lengths are much smaller than previous estimates. These reduced fluid motions cannot generate measurable magnetic fields, suggesting any link with crystallization needs to invoke a separate mechanism. Finally, this mixing alters the chemical profiles, which in turn affects the frequencies of the pulsation modes.
White dwarfs (WDs) are useful across a wide range of astrophysical contexts. The appropriate interpretation of their spectra relies on the accuracy of WD atmosphere models. One essential ingredient of atmosphere models is the theory used for the broadening of spectral lines. To date, the models have relied on Vidal et al., known as the unified theory of line broadening (VCS). There have since been advancements in the theory; however, the calculations used in model atmosphere codes have only received minor updates. Meanwhile, advances in instrumentation and data have uncovered indications of inaccuracies: spectroscopic temperatures are roughly 10% higher and spectroscopic masses are roughly 0.1 M ⊙ higher than their photometric counterparts. The evidence suggests that VCS-based treatments of line profiles may be at least partly responsible. Gomez et al. developed a simulation-based line-profile code Xenomorph using an improved theoretical treatment that can be used to inform questions around the discrepancy. However, the code required revisions to sufficiently decrease noise for use in model spectra and to make it computationally tractable and physically realistic. In particular, we investigate three additional physical effects that are not captured in the VCS calculations: ion dynamics, higher-order multipole expansion, and an expanded basis set. We also implement a simulation-based approach to occupation probability. The present study limits the scope to the first three hydrogen Balmer transitions (Hα, Hβ, and Hγ). We find that screening effects and occupation probability have the largest effects on the line shapes and will likely have important consequences in stellar synthetic spectra.
For isolated white dwarf (WD) stars, fits to their observed spectra provide the most precise estimates of their effective temperatures and surface gravities. Even so, recent studies have shown that systematic offsets exist between such spectroscopic parameter determinations and those based on broadband photometry. These large discrepancies (10% in T eff , 0.1 M ⊙ in mass) provide scientific motivation for reconsidering the atomic physics employed in the model atmospheres of these stars. Recent simulation work of ours suggests that the most important remaining uncertainties in simulation-based calculations of line shapes are the treatment of 1) the electric field distribution and 2) the occupation probability (OP) prescription. We review the work that has been done in these areas and outline possible avenues for progress.
PG 1159-035 is the prototype of the DOV hot pre-white dwarf pulsators. It was observed during the Kepler satellite K2 mission for 69 days in 59 s cadence mode and by the TESS satellite for 25 days in 20 s cadence mode. We present a detailed asteroseismic analysis of those data. We identify a total of 107 frequencies representing 32 l=1 modes, 27 frequencies representing 12 l=2 modes, and 8 combination frequencies. The combination frequencies and the modes with very high k values represent new detections. The multiplet structure reveals an average splitting of 4.0+/-0.4 muHz for l=1 and 6.8+/-0.2 muHz for l=2, indicating a rotation period of 1.4+/-0.1 days in the region of period formation. In the Fourier transform of the light curve, we find a significant peak at 8.904+/-0.003 muHz suggesting a surface rotation period of 1.299+/-0.002 days. We also present evidence that the observed periods change on timescales shorter than those predicted by current evolutionary models. Our asteroseismic analysis finds an average period spacing for l=1 of 21.28+/-0.02 s. The l=2 modes have a mean spacing of 12.97+/-0.4 s. We performed a detailed asteroseismic fit by comparing the observed periods with those of evolutionary models. The best fit model has Teff=129600+/- 11100 K, mass M*=0.565+/-0.024 Msun, and log g=7.41+0.38-0.54, within the uncertainties of the spectroscopic determinations. We argue for future improvements in the current models, e.g., on the overshooting in the He-burning stage, as the best-fit model does not predict excitation for all the pulsations detected in PG 1159-03.
Accurate helium White Dwarf (DB) masses are critical for understanding the star’s evolution. DB masses derived from the spectroscopic and photometric methods are inconsistent. Photometric masses agree better with currently accepted DB evolutionary theories and are mostly consistent across a large range of surface temperatures. Spectroscopic masses rely on untested He i Stark line-shape and Van der Waals broadening predictions, show unexpected surface temperature trends, and are thus viewed as less reliable. To test this conclusion, we present in this paper detailed He i Stark line-shape measurements at conditions relevant to DB atmospheres ( T electron ≈12,000–17,000 K, n electron ≈ 10 17 cm −3 ). We use X-rays from Sandia National Laboratories’ Z -machine to create a uniform ≈120 mm long hydrogen–helium mixture plasma. Van der Waals broadening is negligible at our experimental conditions, allowing us to measure He i Stark profiles only. H β , which has been well-studied in our platform and elsewhere, serves as the n e diagnostic. We find that He i Stark broadening models used in DB analyses are accurate within errors at tested conditions. It therefore seems unlikely that line-shape models are solely responsible for the observed spectroscopic mass trends. Our results should motivate the WD community to further scrutinize the validity of other spectroscopic and photometric input parameters, like atmospheric structure assumptions and convection corrections. These parameters can significantly change the derived DB mass. Identifying potential weaknesses in any input parameters could further our understanding of DBs, help elucidate their evolutionary origins, and strengthen confidence in both spectroscopic and photometric masses.
We present a novel method to detect variable astrophysical objects and transient phenomena using anomalous excess scatter in repeated measurements from public catalogs of Gaia DR2 and Zwicky Transient Facility (ZTF) DR3 photometry. We first provide a generalized, all-sky proxy for variability using only Gaia DR2 photometry, calibrated to white dwarf stars. To ensure more robust candidate detection, we further employ a method combining Gaia with ZTF photometry and alerts. To demonstrate its efficacy, we apply this latter technique to a sample of roughly 12,100 white dwarfs within 200 pc centered on the ZZ Ceti instability strip, where hydrogen-atmosphere white dwarfs are known to pulsate. By inspecting the top 1% of the samples ranked by these methods, we demonstrate that both the Gaia-only and ZTF-informed techniques are highly effective at identifying known and new variable white dwarfs, which we verify using follow-up, high-speed photometry. We confirm variability in all 33 out of 33 (100%) observed white dwarfs within our top 1% highest-ranked candidates, both inside and outside the ZZ Ceti instability strip. In addition to dozens of new pulsating white dwarfs, we also identify five white dwarfs highly likely to show transiting planetary debris; if confirmed, these systems would more than triple the number of white dwarfs known to host transiting debris.
Understanding how atoms interact with hot dense matter is essential for astrophysical and laboratory plasmas. Interactions in high-density plasmas broaden spectral lines, providing a rare window into interactions that govern, for example, radiation transport in stars. However, up to now, spectral line-shape theories employed at least one of three common approximations: second-order Taylor treatment of broadening operator, dipole-only interactions between atom and plasma, and classical treatment of perturbing electrons. In this Letter, we remove all three approximations simultaneously for the first time and test the importance for two applications: neutral hydrogen and highly ionized magnesium and oxygen. We found 15%-50% change in the spectral line widths, which are sufficient to impact applications including white-dwarf mass determination, stellar-opacity research, and laboratory plasma diagnostics.
White Dwarf (WD) stars are the most common stellar remnant in the universe. WDs usually have a hydrogen or helium atmosphere, and helium WD (called DB) spectra can be used to solve outstanding problems in stellar and galactic evolution. DB origins, which are still a mystery, must be known to solve these problems. DB masses are crucial for discriminating between different proposed DB evolutionary hypotheses. Current DB mass determination methods deliver conflicting results. The spectroscopic mass determination method relies on line broadening models that have not been validated at DB atmosphere conditions. We performed helium benchmark experiments using the White Dwarf Photosphere Experiment (WDPE) platform at Sandia National Laboratories' Z-machine that aims to study He line broadening at DB conditions. Using hydrogen/helium mixture plasmas allows investigating the importance of He Stark and van der Waals broadening simultaneously. Accurate experimental data reduction methods are essential to test these line-broadening theories. In this paper, we present data calibration methods for these benchmark He line shape experiments. We give a detailed account of data processing, spectral power calibrations, and instrument broadening measurements. Uncertainties for each data calibration step are also derived. We demonstrate that our experiments meet all benchmark experiment accuracy requirements: WDPE wavelength uncertainties are <1 Å, spectral powers can be determined to within 15%, densities are accurate at the 20% level, and instrumental broadening can be measured with 20% accuracy. Fulfilling these stringent requirements enables WDPE experimental data to provide physically meaningful conclusions about line broadening at DB conditions.
and 2) the occupation probability (OP) prescription. We review the work that has been done in these areas and outline possible avenues for progress.