🔭 astrophysics

Physics of Eclipsing Binaries. VI. Hot, compact stars

This paper introduces the PHOEBE 2.5 release, which enhances the modeling of eclipsing binaries containing hot, compact stars by incorporating specialized model atmospheres (TMAP and Montreal/Tremblay), enabling atmosphere blending for grid extrapolation, and providing new limb-darkening coefficient tables.

Original authors: David Jones, Nicole Reindl, Kyle Conroy, James Munday, Pier-Emmanuel Tremblay, Michael Abdul-Masih, Matthias Fabry, Joseph Giammarco, Kelly Hambleton, Herbert Pablo, Marcin Wrona, Andrej Prsa

Published 2026-07-10

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: David Jones, Nicole Reindl, Kyle Conroy, James Munday, Pier-Emmanuel Tremblay, Michael Abdul-Masih, Matthias Fabry, Joseph Giammarco, Kelly Hambleton, Herbert Pablo, Marcin Wrona, Andrej Prsa

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine trying to take a perfect photo of a cosmic dance between two stars. To get the picture right, you need to know exactly how bright each star is and how that brightness changes as you look at different parts of their surfaces. For a long time, the computer program scientists use to model these dances—called PHOEBE—had a bit of a blind spot. It was great at modeling "normal" stars, like our Sun, but when it came to the super-hot, super-dense stars (like white dwarfs and subdwarf O stars), the program had to guess. It used a "blackbody" approximation, which is like trying to describe a complex, colorful painting using only a single shade of gray.

In this paper, the team behind PHOEBE (led by David Jones and Nicole Reindl) announces a major upgrade: version 2.5. They've swapped out that boring gray guess for a high-definition, multi-colored map specifically designed for hot, compact stars.

The New "Weather Maps" for Stars
Think of the old "blackbody" method as a generic weather forecast that just says "it's hot." It's okay for a rough idea, but it fails to capture the real details. The new update brings in two specific sets of "weather maps" (model atmospheres) created by the TMAP and Montréal/Tremblay teams. These maps are like detailed, 3D topographical maps of the star's surface, showing exactly how light behaves at different temperatures and gravity levels.

The authors found that using the old "gray" guess was leading to some pretty big mistakes. In their simulations, when they modeled a hot, compact star using the old blackbody method, the program thought the star was about 25% cooler than it actually was. To make the simulation match the real light curve, they had to fudge the numbers, giving the companion star an "albedo" (reflectivity) of 2.0. That's like saying a mirror reflects twice as much light as it receives—physically impossible! By using the new TMAP maps, the program gets the temperature right without needing to invent impossible physics.

The "Blending" Safety Net
Here's a tricky problem: sometimes, a star in a binary system gets stretched or heated up by its partner so much that a few tiny spots on its surface end up in a "forbidden zone" where the new maps don't have data yet. In the past, the computer would just crash and say, "Error! I can't do this."

The new PHOEBE 2.5 has a clever safety net called "blending." Imagine you are walking from a well-lit room (where the data is perfect) into a foggy hallway (where the data is missing). Instead of tripping in the dark, the program smoothly blends the detailed map with the old "gray" guess as you step into the fog. It doesn't just stop; it gently transitions, so the simulation keeps running even if a few tiny parts of the star are outside the known data. This means scientists can now model weird, distorted stars without the whole program breaking.

New Tricks for the Program
The update also turns on a feature called "Doppler boosting" (or beaming). This is a relativistic effect where a star moving toward us looks slightly brighter, and moving away looks slightly dimmer. The program had turned this off in a previous version because the math was tricky, but now it's back on, letting users plug in their own precise numbers to see these subtle effects.

Finally, the team gave users a new "Lego kit" called "User-Defined Features." Before, if you wanted to model a star with a weirdly moving spot or a changing rotation speed, you were stuck. Now, you can write your own little code "hooks" to tell the program exactly how to move those spots or change the star's shape. It's like giving the user a remote control to tweak the simulation in real-time, allowing for things like spots that migrate across the star's surface as it spins.

The Bottom Line
The paper doesn't claim to have solved every mystery in the universe, but it has definitely sharpened the tools. By swapping out the rough "blackbody" guesses for detailed TMAP and Tremblay maps, and by adding safety nets for when data runs out, PHOEBE 2.5 allows astronomers to measure the true temperatures and sizes of hot, compact stars much more accurately. The authors show through their simulations that sticking with the old methods could lead to underestimating a star's temperature by a significant margin, but with these new tools, the cosmic dance is finally being captured in high definition.

Technical Summary: Physics of Eclipsing Binaries. VI. Hot, Compact Stars

Problem Statement
Accurate modeling of eclipsing binary light curves requires the assignment of appropriate emergent intensities to the surface elements of the binary components. While codes like PHOEBE utilize model atmospheres to derive these intensities based on effective temperature ( $T_{\text{eff}}$ ), surface gravity ( $\log g$ ), and chemical composition, existing implementations have historically been limited to specific stellar types. Previous versions relied on Castelli & Kurucz (CK2004) and PHOENIX models, which are best suited for main-sequence stars and low-temperature giants. Consequently, for hot, compact stars (such as hot white dwarfs and subdwarf O stars), models often default to blackbody approximations. This approximation is demonstrably poor for these objects, leading to significant inaccuracies in derived physical parameters. Furthermore, previous versions of the code would fail entirely if specific surface elements (due to irradiation, tidal distortion, or rapid rotation) fell outside the parameter bounds of the available model atmosphere grids.

Methodology
This paper details the implementation of new capabilities in the PHOEBE 2.5 release, focusing on the integration of non-LTE model atmospheres for hot, compact stars and the development of robust handling for out-of-bounds parameters.

Integration of New Atmospheres:
- TMAP Grids: The authors incorporated plane-parallel model atmospheres from the Tübingen Model Atmosphere Package (TMAP). These cover four specific grids: sdO, DA, DAO, and DO, spanning $T_{\text{eff}}$ from 20 kK to 200 kK and $\log g$ from 4.0 to 9.5. A new parameter, loghefrac (logarithmic Helium to Hydrogen number fraction), was introduced to define abundances for these grids, distinguishing them from the metallicity-based definitions of previous models.
- Montréal/Tremblay Grids: The release also includes LTE Montréal/Tremblay models for DA white dwarfs, covering $T_{\text{eff}}$ from 3.75 kK to 60 kK. The authors explicitly excluded the non-LTE TLUSTY models from Claret et al. (2020a) due to identified calculation errors in the intensities provided in that work.
Blending and Extrapolation:
To address the issue of surface elements falling outside grid boundaries, PHOEBE 2.5 implements two new operational modes:
- Extrapolation: Users can select 'nearest' (using the nearest defined node) or 'linear' (linear extrapolation from the nearest fully defined hypercube) to estimate values for out-of-bounds parameters.
- Blending: This mode creates a smooth transition between the extrapolated model atmosphere values and blackbody values. The blending ratio depends on the distance from the model atmosphere grid, preventing the code from failing when only a few surface elements are out of bounds.
Additional Functionalities:
- Doppler Boosting: Relativistic boosting (beaming) functionality has been reactivated. Unlike previous versions, users must now provide boosting factors calculated externally (e.g., from literature) rather than relying on internal interpolation.
- User-Defined Features: A new architecture allows users to inject custom code into the PHOEBE logic. This enables modifications to the mesh, local values (radial velocities, temperatures), or the final synthetic model, facilitating features like migrating spots, differential rotation, and variable third light.

Key Results and Validation

Atmosphere Consistency: The authors verified the consistency between the TMAP and Tremblay grids. Comparisons of normal intensities and limb-darkening coefficients in the SDSS $g'$ band show agreement within one standard deviation of unity (ratio $1.01 \pm 0.10$ ). The transition between the four TMAP grids (sdO, DA, DAO, DO) was shown to be continuous regarding limb-darkening coefficients and normal intensity as a function of abundance.
Impact of Blackbody Approximation: A synthetic binary system (a 127 kK pre-white dwarf with a K5V companion) was modeled. The study found that using a blackbody approximation for the hot component results in an underestimation of the effective temperature by approximately 25% (requiring a reduction to ~94 kK to match the TMAP light curve). Furthermore, reproducing the irradiation effect amplitude with a blackbody required an unphysical bolometric albedo of 2.0 for the companion, compared to a physical value of 0.6 in the TMAP model.
Limb-Darkening Coefficients: The authors generated and provided tables of limb-darkening coefficients derived from the new TMAP and Tremblay models for standard filter sets (SDSS and Johnson), suitable for use in other binary modeling codes.

Significance
The paper claims that the inclusion of TMAP and Tremblay atmospheres, coupled with the blending/extrapolation capabilities, significantly improves the fidelity of PHOEBE for hot and compact stars. The primary significance lies in eliminating the reliance on the blackbody approximation for these objects, which the authors demonstrate can lead to substantial errors in derived stellar parameters (specifically effective temperature and luminosity). By allowing the code to gracefully handle surface elements that fall outside grid limits, the new release enables the modeling of complex systems involving strong irradiation, tidal distortion, and rapid rotation without resorting to error states or inaccurate approximations. Additionally, the reactivation of Doppler boosting and the introduction of user-defined features expand the code's utility for analyzing compact binaries and dynamic stellar phenomena.

More like this