White Dwarf Research Articles

Main sequence dynamo magnetic fields emerging in the white dwarf phase

Main sequence dynamo magnetic fields emerging in the white dwarf phase

Finite-field Cholesky decomposed coupled-cluster techniques (ff-CD-CC): theory and application to pressure broadening of Mg by a He atmosphere and a strong magnetic field.

For the interpretation of spectra of magnetic stellar objects such as magnetic white dwarfs (WDs), highly accurate quantum chemical predictions for atoms and molecules in finite magnetic field are required. Especially the accurate description of electronically excited states and their properties requires established methods such as those from coupled-cluster (CC) theory. However, respective calculations are computationally challenging even for medium-sized systems. Cholesky decomposition (CD) techniques may be used to alleviate memory bottlenecks. In finite magnetic field computations, the latter are increased due to the reduction of permutational symmetry within the electron-repulsion-integrals (ERIs) as well as the need for complex-valued data types. CD enables a memory-efficient, approximate description of the ERIs with rigorous error control and thus the treatment of larger systems at the CC level becomes feasible. In order to treat molecules in a finite magnetic field, we present in this work the working equations of the left and right-hand side equations for finite field (ff)-EOM-CD-CCSD for various EOM flavours as well as for the approximate ff-EOM-CD-CC2 method. The methods are applied to the study of the modification of the spectral lines of a magnesium transition by a helium atmosphere that can be found on magnetic WD stars.

Characterising high and low accretion states in VY Scl CVs using ZTF and TESS data

Abstract VY Scl binaries are a sub-class of cataclysmic variable (CV) which show extended low states, but do not show outbursts which are seen in other classes of CV. To better determine how often these systems spend in low states and to resolve the state transitions we have analysed ZTF data on eight systems and TESS data on six systems. Half of the sample spent most of the time in a high state; three show a broad range and one spends roughly half the time transitioning between high and low states. Using the ZTF data we explore the colour variation as a function of brightness. In KR Aur, we identify a series of repeating outburst events whose brightness appears to increase over time. Using TESS data we searched for periods other than the orbital. In LN UMa we find evidence for a peak whose period varies between 3–6 d. We outline the current models which aim to explain the observed properties of VY Scl systems which includes disc irradiation and a white dwarf having a significant magnetic field.

Nuclear burning in an accretion flow around a stellar-mass black hole embedded within an AGN disk

Abstract A stellar-mass black hole, embedded within the accretion disk of an active galactic nuclei (AGN), has the potential to accrete gas at a rate that can reach approximately ∼109 times the Eddington limit. This study explores the potential for nuclear burning in the rapidly accreting flow towards this black hole and studies how nucleosynthesis affects metal production. Using numerical methods, we have obtained the disk structure while considering nuclear burning and assessed the stability of the disk. In contrast to gas accretion onto the surface of a neutron star or white dwarf, the disk remains stable against the thermal and secular instabilities because advection cooling offsets the nuclear heating effects. The absence of a solid surface for a black hole prevents excessive mass accumulation in the inner disk region. Notably, nuclear fusion predominantly takes place in the inner disk region, resulting in substantial burning of $\rm ^{12}C$ and $\rm ^{3}He$, particularly for black holes around M = 10 M⊙ with accretion rates exceeding approximately ∼107 times the Eddington rate. The ejection of carbon-depleted gas through outflows can lead to an increase in the mass ratio of oxygen or nitrogen to carbon, which may be reflected in observed line ratios such as $\rm N\, V/C\, IV$ and $\rm O\, IV/C\, IV$. Consequently, these elevated spectral line ratios could be interpreted as indications of super-solar metallicity in the broad line region.

A catalog of binary stars from phase modulation in the first four years of TESS Mission photometry

Abstract We present a catalog of binary companions to δ Scuti stars, detected through phase modulations of their pulsations in TESS data. Pulsation timing has provided orbits for hundreds of pulsating stars in binaries from space-based photometry. We have applied this technique to δ Sct stars observed in the first four years of TESS Mission photometry. We searched the 2-min cadence light curves of 1161 short-period instability strip pulsators for variations in pulsation phase caused by the dynamical influence of an unseen companion. We discovered 53 new binaries and we present orbital parameters and mass functions for the 24 systems with solvable orbits. For the brightest star in our sample α Pictoris, we perform a joint fit of the pulsation timing and Hipparcos astrometry. We present the first orbit for the α Pictoris system, obtaining an orbital period of 1316±2 days and a mass for α Pic B of 1.05±0.05 M⊙. We revisit pulsation timing binaries from Kepler with Gaia kinematics, finding four systems that are members of the Galatic thick disk or halo. This suggests that they have been rejuvenated by mass transfer, and their companions are now white dwarfs. Further follow up of these systems may yield valuable constraints of the galactic blue straggler population.

Thermal Conductivity and Thermal Hall Effect in Dense Electron-Ion Plasma

In this study, we examine thermal conductivity and the thermal Hall effect in electron-ion plasmas relevant to hot neutron stars, white dwarfs, and binary neutron star mergers, focusing on densities found in the outer crusts of neutron stars and the interiors of white dwarfs. We consider plasma consisting of single species of ions, which could be either iron [56]Fe or carbon [12]C nuclei. The temperature range explored is from the melting temperature of the solid T∼109 K up to 1011 K. This covers both degenerate and non-degenerate electron regimes. We find that thermal conductivity increases with density and temperature for which we provide analytical scaling relations valid in different regimes. The impact of magnetic fields on thermal conductivity is also analyzed, showing anisotropy in low-density regions and the presence of the thermal Hall effect characterized by the Righi–Leduc coefficient. The transition froma degenerate to non-degenerate regime is characterized by a minimum ratio of thermal conductivity to temperature, which is analogous to the minimum observed already in the case of electrical conductivity. We provide also formulas fit to our numerical results, which can be used in dissipative magneto-hydrodynamics simulations of warm compact stars.

Cataclysmic variables from Sloan Digital Sky Survey - V (2020-2023) identified using machine learning.

Abstract SDSS-V is carrying out a dedicated survey for white dwarfs, single and in binaries, and we report the analysis of the spectroscopy of 504 cataclysmic variables (CVs) and CV candidates obtained during the first 34 months of observations of SDSS-V. We developed a convolutional neural network (CNN) to aid with the identification of CV candidates among the over 2 million SDSS-V spectra obtained with the BOSS spectrograph. The CNN reduced the number of spectra that required visual inspection to ≃ 2 per cent of the total. We identified 776 CV spectra among the CNN-selected candidates, plus an additional 27 CV spectra that the CNN misclassified, but that were found serendipitously by human inspection of the data. Analysing the SDSS-V spectroscopy and ancillary data of the 504 CVs in our sample, we report 61 new CVs, spectroscopically confirm 248 and refute 13 published CV candidates, and we report 82 new or improved orbital periods. We discuss the completeness and possible selection biases of the machine learning methodology, as well as the effectiveness of targeting CV candidates within SDSS-V. Finally, we re-assess the space density of CVs, and find 1.2 × 10−5 pc−3.

WISE J152614.95-111326.4, an Unusual Variable Star

ABSTRACT New time series photometry of WISE J152614.95-111326.4, an eclipsing binary candidate, has been obtained. Full cycles of variation were covered in five filters, ranging from B to z. Archival time series photometry is also available from several sources. The phased light curve shape changes from a double wave form in the red, to a single wave at shorter wavelengths. Analysis of the spectral energy distribution and SALT spectra shows the presence of a cool (∼ 7250−7900 K) white dwarf and an M6 star. The light curves can be explained by a hot spot on the opposing hemisphere of the white dwarf. The star may be in a pre-cataclysmic variable phase with a very low rate of mass flow from the red dwarf to the white dwarf, such that no flickering is evident. Evidence in favour of this hypothesis is that the period of the system (2.25 h) is in the cataclysmic variable period gap. It is speculated that a weak magnetic field associated with the white dwarf funnels accreted material onto a magnetic pole. Amplitudes of the W1 and W2 WISE light curves are anomalously large. The possibility is discussed that variability in this spectral region is primarily driven by electron cyclotron radiation.

Semi-supervised spectral classification of DESI white dwarfs by dimensionality reduction

Abstract As a new generation of large-sky spectroscopic surveys comes online, the enormous data volume poses unprecedented challenges in classifying spectra. Modern unsupervised techniques have the power to group spectra based on their dominant features, circumventing the complete reliance on training data suffered by supervised methods. We outline the use of dimensionality reduction to generate a 2D map of the structure of an intermediate-resolution spectroscopic dataset. This technique efficiently separates white dwarfs of different spectral classes in the Dark Energy Spectroscopic Instrument’s Early Data Release (DESI EDR), identifying spectral features that had been missed even by visual classification. By focusing the method on particular spectral regions, we identify white dwarfs with helium features at 90percnt recall, and cataclysmic variables at 100percnt recall, illustrating rapid selection of low-contamination samples from spectroscopic surveys. We also demonstrate the use of dimensionality reduction in a supervised manner, outlining a procedure to classify any white dwarf spectrum in comparison with those in the DESI EDR. With upcoming surveys promising tens of millions of spectra, our work highlights the potential for semi-supervised techniques as an efficient means of classification and dataset visualisation.

Earth May Survive the Sun’s Demise

A distant white dwarf hosts an Earth-like planet in an orbit that might be similar to Earth’s if it survives the Sun’s red giant phase.

J-PLUS: The fraction of calcium white dwarfs along the cooling sequence

Aims. We used the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2) photometry in 12 optical bands over 2176 deg2 to estimate the fraction of white dwarfs with the presence of Ca II H+K absorption along the cooling sequence. Methods. We compared the J-PLUS photometry against metal-free theoretical models to estimate the equivalent width in the J0395 passband of 10 nm centered at 395 nm (EWJ0395), a proxy to detect calcium absorption. A total of 4399 white dwarfs with effective temperatures within 30 000 > Teff > 5500 K and masses of M > 0.45 M⊙ were analyzed. Their EWJ0395 distribution was modeled using two populations, corresponding to polluted and non-polluted systems, to estimate the fraction of calcium white dwarfs (fCa) as a function of Teff. The probability of each individual white dwarf presenting calcium absorption, pCa, was also computed. Results. The comparison of EWJ0395 with both the measured Ca/He abundance and the identification of metal pollution from spectroscopy shows that EWJ0395 correlates with the presence of Ca II H+K absorption. The fraction of calcium white dwarfs changes along the cooling sequence, increasing from fCa ≈ 0 at Teff = 13 500 K to fCa ≈ 0.15 at Teff = 5500 K. This trend reflects the selection function of calcium white dwarfs in the optical. We compare our results with the fractions derived from the 40 pc spectroscopic sample and from Sloan Digital Sky Survey (SDSS) spectra. The trend found in J-PLUS observations is also present in the 40 pc sample; however, SDSS shows a deficit of metal-polluted objects at Teff < 12 000 K. Finally, we found 39 white dwarfs with pCa > 0.99. Twenty of them have spectra presented in previous studies, whereas we obtained follow-up spectroscopic observations for six additional targets. These 26 objects were all confirmed as metal-polluted systems. Conclusions. The J-PLUS optical data provide a robust statistical measurement for the presence of Ca II H+K absorption in white dwarfs. We find a 15 ± 3% increase in the fraction of calcium white dwarfs from Teff = 13 500 K to 5500 K, which reflects their selection function in the optical from the total population of metal-polluted systems.

He-accreting Oxygen–Neon White Dwarfs and Accretion-induced Collapse Events

It has been widely accepted that mass-accreting white dwarfs (WDs) are the progenitors of Type Ia supernovae (SNe) or electron-capture supernovae. Previous work has shown that the accretion rate could affect the elemental abundance on the outer layers of CO WDs, and therefore affect the observational characteristics after they exploded as SNe Ia. However, it has not been well studied how elemental abundance changes on the outer layers of He-accreting ONe WDs as they approach the Chandrasekhar mass limit. In this paper, we investigated the evolution of He-accreting ONe WDs with MESA. We found that a CO-rich mantle will accumulate beneath the He layers resulting from the He burning, after which the ignition of the CO-rich mantle could transform carbon into silicon (Si). The amount of Si produced by carbon burning is strongly anticorrelated with the accretion rate. As the ONe WD nearly approaches the Chandrasekhar mass limit (M ch) through accretion, it is likely to undergo accretion-induced collapse, resulting in the formation of a neutron star.

An analysis of spectroscopic, seismological, astrometric, and photometric masses of pulsating white dwarf stars

Context. A central challenge in the field of stellar astrophysics lies in accurately determining the mass of stars, particularly when dealing with isolated ones. However, for pulsating white dwarf stars, the task becomes more tractable due to the availability of multiple approaches such as spectroscopy, asteroseismology, astrometry, and photometry, each providing valuable insights into the mass properties of white dwarf stars. Aims. Numerous asteroseismological studies of white dwarfs have been published, focusing on determining stellar mass using pulsational spectra and comparing it with spectroscopic mass, which uses surface temperature and gravity. The objective of this work is to compare these mass values in detail and, in turn, to compare them with the mass values derived using astrometric parallaxes or distances and photometry data from Gaia, employing astrometric and photometric methods. Methods. Our analysis involves a selection of pulsating white dwarfs with different surface chemical abundances that define the main classes of variable white dwarfs. We calculated their spectroscopic masses, compiled seismological masses, and determined astrometric masses. We also derived photometric masses, when possible. Subsequently, we compared all the sets of stellar masses obtained through these different methods. To ensure consistency and robustness in our comparisons, we used identical white dwarf models and evolutionary tracks across all four methods. Results. The analysis suggests a general consensus among the four methods regarding the masses of pulsating white dwarfs with hydrogen-rich atmospheres, known as DAV or ZZ Ceti stars, especially for objects with masses below approximately 0.75 M⊙, although notable disparities emerge for certain massive stars. For pulsating white dwarf stars with helium-rich atmospheres, called DBV or V777 Her stars, we find that astrometric masses generally exceed seismological, spectroscopic, and photometric masses. Finally, while there is agreement among the sets of stellar masses for pulsating white dwarfs with carbon-, oxygen-, and helium-rich atmospheres (designated as GW Vir stars), outliers exist, where mass determinations by various methods show significant discrepancies. Conclusions. Although a general agreement exists among different methodologies for estimating the mass of pulsating white dwarfs, significant discrepancies are prevalent in many instances. This shows the need to redo the determination of spectroscopic parameters and the parallax and/or improve asteroseismological models for many stars.

A JWST Medium-resolution MIRI Spectrum and Models of the Type Ia Supernova 2021aefx at +415 days

We present a JWST MIRI medium-resolution spectrometer spectrum (5–27 μm) of the Type Ia supernova (SN Ia) SN 2021aefx at +415 days past B-band maximum. The spectrum, which was obtained during the iron-dominated nebular phase, has been analyzed in combination with previous JWST observations of SN 2021aefx to provide the first JWST time series analysis of an SN Ia. We find that the temporal evolution of the [Co iii] 11.888 μm feature directly traces the decay of 56Co. The spectra, line profiles, and their evolution are analyzed with off-center delayed-detonation models. Best fits were obtained with white dwarf (WD) central densities of ρ c = 0.9−1.1 × 109 g cm−3, a WD mass of M WD = 1.33–1.35 M ⊙, a WD magnetic field of ≈106 G, and an off-center deflagration-to-detonation transition at ≈0.5 M ⊙ seen opposite to the line of sight of the observer (−30°). The inner electron capture core is dominated by energy deposition from γ-rays, whereas a broader region is dominated by positron deposition, placing SN 2021aefx at +415 days in the transitional phase of the evolution to the positron-dominated regime. The formerly “flat-tilted” profile at 9 μm now has a significant contribution from [Ni iv], [Fe ii], and [Fe iii] and less from [Ar iii], which alters the shape of the feature as positrons mostly excite the low-velocity Ar. Overall, the strength of the stable Ni features in the spectrum is dominated by positron transport rather than the Ni mass. Based on multidimensional models, our analysis is consistent with a single-spot, close-to-central ignition with an indication of a preexisting turbulent velocity field and excludes a multiple-spot, off-center ignition.

Impact of Current Uncertainties in the 12C+12C Nuclear Reaction Rate on Intermediate-mass Stars and Massive White Dwarfs

Recent determinations of the total rate of the 12C+12C nuclear reaction show non-negligible differences with the reference reaction rate commonly used in previous stellar simulations. In addition, the current uncertainties in determining each exit channel constitute one of the main uncertainties in shaping the inner structure of super asymptotic giant branch stars that could have a measurable impact on the properties of pulsating ultramassive white dwarfs (WDs). We explore how new determinations of the nuclear reaction rate and its branching ratios affect the evolution of WD progenitors. We show that the current uncertainties in the branching ratios constitute the main uncertainty factor in determining the inner composition of ultramassive WDs and their progenitors. We found that the use of extreme branching ratios leads to differences in the central abundances of 20Ne of at most 17%, which are translated into differences of at most 1.3% and 0.8% in the cooling times and size of the crystallized core, respectively. However, the impact on the pulsation properties is small, less than 1 s for the asymptotic period spacing. We found that the carbon burns partially in the interior of ultramassive WD progenitors within a particular range of masses, leaving a hybrid CONe-core composition in their cores. The evolution of these new kinds of predicted objects differs substantially from the evolution of objects with pure CO cores. Differences in the size of the crystallized core and cooling times of up to 15% and 6%, respectively, lead to distinct patterns in the period spacing distribution.

A Broadband X-Ray Investigation of Fast-spinning Intermediate Polar CTCV J2056–3014

Abstract We report on XMM-Newton, NuSTAR, and NICER X-ray observations of CTCV J2056–3014, a cataclysmic variable (CV) with one of the fastest-spinning white dwarfs (WDs) at P = 29.6 s. While previously classified as an intermediate polar, CJ2056 also exhibits the properties of WZ Sge–type CVs, such as dwarf novae and superoutbursts. With XMM-Newton and NICER, we detected the spin period up to ∼2 keV with 7σ significance. We constrained its derivative to | P ̇ | < 1.8 × 10 − 12 s s−1 after correcting for binary orbital motion. The pulse profile is characterized by a single broad peak with ∼25% modulation. NuSTAR detected a fourfold increase in unabsorbed X-ray flux coincident with an optical flare, in 2022 November. The XMM-Newton and NICER X-ray spectra at 0.310 keV are best characterized by an absorbed, optically thin three-temperature thermal plasma model (kT = 0.3, 1.0, and 4.9 keV), while the NuSTAR spectra at 3–30 keV are best fit by a single-temperature thermal plasma model (kT = 8.4 keV), both with Fe abundance Z Fe/Z ⊙ = 0.3. CJ2056 exhibits similarities to other fast-spinning CVs, such as low plasma temperatures and no significant X-ray absorption at low energies. As the WD’s magnetic field strength is unknown, we applied both nonmagnetic and magnetic CV spectral models (MKCFLOW and MCVSPEC) to determine the WD mass. The derived WD mass range (M = 0.7–1.0 M ⊙) is above the centrifugal breakup mass limit of 0.56 M ⊙ and consistent with the mean WD mass of local CVs (M ≈ 0.8–0.9 M ⊙).

The First Catalog of Candidate White Dwarf–Main-sequence Binaries in Open Star Clusters: A New Window into Common Envelope Evolution

Close binary systems are the progenitors to both Type Ia supernovae and the compact object mergers that can be detected via gravitational waves. To achieve a binary with a small radial separation, it is believed that the system likely undergoes common envelope (CE) evolution. Despite its importance, CE evolution may be one of the largest uncertainties in binary evolution due to a combination of computational challenges and a lack of observed benchmarks where both the post-CE and pre-CE conditions are known. Identifying post-CE systems in star clusters can partially circumvent this second issue by providing an independent age constraint on the system. For the first time, we conduct a systematic search for white dwarf and main-sequence binary systems in 299 Milky Way open star clusters. Coupling Gaia DR3 photometry and kinematics with multiband photometry from Pan-STARRS1 and the Two Micron All Sky Survey, we apply a machine learning-based approach and find 52 high-probability candidates in 38 open clusters. For a subset of our systems, we present follow-up spectroscopy from the Gemini and Lick Observatories and archival light curves from the Transiting Exoplanet Survey Satellite, Kepler/K2, and the Zwicky Transient Facility. Examples of M dwarfs with hot companions are spectroscopically observed, along with regular system variability. While the kinematics of our candidates are consistent with their host clusters, some systems have spatial positions offset relative to their hosts, potentially indicative of natal kicks. Ultimately, this catalog is a first step to obtaining a set of observational benchmarks to better link post-CE systems to their pre-CE progenitors.

Asteroseismology of One of the Most Rapidly Rotating DBV Stars: EPIC 248705247

We present a comprehensive analysis of oscillations of the helium-rich white dwarf DBV star EPIC 248705247 (SDSS J102106.69+082724.8). We extract the light curves from 79.7 days of K2 Campaign 14 data and identify 21 independent frequencies, including three complete and three incomplete dipole modes, with an average frequency splitting of 20.37 ± 0.40 μHz. Our analyses reveal a rotation period of EPIC 248705247 of approximately 6.82 ± 0.07 hr, making it one of the most rapidly rotating DBV stars currently known. We utilize the White Dwarf Evolution Code (WDEC) to compute models to look for the seismic best-fitted model, yielding the following stellar parameters: M/M ⊙ = 0.650 ± 0.003, T eff = 23,660 ± 68 K, log(Menv/M*)=−2.01±0.01 and log(MHe/M*)=−4.90±0.05 . We also explore the phenomenon of mode trapping, as well as the potential influence of a weak magnetic field. We find that this star is located near the red edge of the theoretical instability strip of DBVs. Hence, our asteroseismological analysis of EPIC 248705247 provides the properties and behaviors of a rapidly rotating DBV.

Common Envelopes, Gamma Rays, and Sudden Spectral Changes of Novae

A common envelope (CE) is proposed as the origin of the early post-outburst spectra of many novae. A simple model is proposed to explain the properties of the CE based on the emission line strengths and an assumed density distribution. Rapid changes in the spectrum during post-outburst decline are suggested as possible evidence for a CE. Time-resolved spectra from the ARAS group show sudden spectral shifts that are correlated with detected γ-ray emission, suggestive of its possible origin on the white dwarf that produces a change in condition within the CE. Episodic mass loss, formation of transient heavy element absorption systems, and dissipation of the CE may possibly be triggered by γ-ray emission.

The HST Large Programme on ω Centauri

We present a study of the white dwarf (WD) cooling sequence (CS) in the globular cluster (GC) Omega Centauri (or NGC 5139; hereafter, ω Cen), the primary goal of a dedicated Hubble Space Telescope (HST) programme. We find that the peak at the termination of the WD CS is located at mF606W = 30.1 ± 0.2 (equivalent to V ∼ 31). The brighter part of ω Cen’s WD CS is consistent with the presence of massive He-core WDs, in agreement with previous HST analyses with ultraviolet and blue filters. Comparative analyses of the WD luminosity function (LF) and theoretical counterparts show that a single-age population for the cluster is compatible with the data. However, an analysis of only the WD LF cannot entirely exclude the possibility of an age range, due to uncertainties in the present-day WD mass function, with a star formation history potentially spanning up to 5 billion years, predominantly comprising stars about 13 Gyr old, with a minority potentially as young as 8 Gyr. This underscores the need for global spectroscopic and photometric investigations that simultaneously include both the WD populations and the previous evolutionary phases, in order to fully understand the cluster’s diverse chemical compositions and ages.