Gloria Koenigsberger
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Research Interests
Gloria Koenigsberger has a total of 99 co-authors. Here are the top co-authors by number of shared publications:
- Werner Schmutz
Physikalisch-Meteorologisches Observatorium Davos, World Radiation Center
7 shared publications - D. M. Harrington
Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI, 96822, USA
6 shared publications - Diana Estrella-Trujillo
Instituto de Ciencias Físicas Universidad Nacional Autónoma de México Ave. Universidad S/N Chamilpa Cuernavaca México
6 shared publications - Edmundo Moreno
Instituto de Astronomía Universidad Nacional Autonoma de Mexico Mexico D.F Mexico
5 shared publications - S. Jane Arthur
Universidad Nacional Autónoma de México
5 shared publications - Catherine A. Pilachowski
Indiana University Bloomington
4 shared publications
Latest Publications
- Open Access
- 4 September 2026
- Published by EDP Sciences in Astronomy & Astrophysics
- Vol. 713, A71
Revised orbital parameters of the γ 2 Velorum system
Context.γ 2 Velorum is the closest and visually brightest Wolf-Rayet (WR) binary system. Its eccentric orbit modulates the X-rays observed from the wind-wind interaction, and its large separation enables spatial resolution of both components. Aims. We aim to strengthen the constraints on the properties of γ 2 Velorum and, in particular, to resolve the discrepancy between the eccentricity determined from the emission lines and that from the absorption lines. Methods. We obtained VLT/GRAVITY observations and combined them with earlier spatially resolved data covering different orbital phases. Results. We determine strong constraints on all orbital parameters and find that e = 0.322, close to the value derived from the emission lines. The X-ray light curve declines as s −3 after periastron, where s is the separation of the two stars, but its modulation is likely affected by absorption and occultation of the X-ray emitting region at other orbital phases. We find that previous discrepancies in the reddening value can be traced to a brighter K -band magnitude than predicted by the WR wind models. We derive E ( B − V ) = 0.02 ± 0.02 mag. Our more precise mass and radius values, combined with previously determined effective temperatures, provide strong constraints on evolutionary models. The closest match for the O star comes from a rotationally mixed model with initial mass M init = 28.7 M ⊙ and an M init = 32 M ⊙ model for the WR star, with negligible accretion onto the O star during the WR progenitor’s Roche lobe overflow phase. However, the temperature of the O star is higher and the mass of the WR star is smaller than predicted by the evolutionary tracks at the current epoch, consistent with the well-known “mass-discrepancy problem” in massive stars. Conclusions.γ 2 Velorum’s strongly constrained parameters make it ideal for testing structure, evolution, and stellar wind models.
- 4 September 2026
- Published by EDP Sciences in Astronomy & Astrophysics
- Vol. 713, A71
Open AccessRevised orbital parameters of the γ 2 Velorum system
Context.γ 2 Velorum is the closest and visually brightest Wolf-Rayet (WR) binary system. Its eccentric orbit modulates the X-rays observed from the wind-wind interaction, and its large separation enables spatial resolution of both components. Aims. We aim to strengthen the constraints on the properties of γ 2 Velorum and, in particular, to resolve the discrepancy between the eccentricity determined from the emission lines and that from the absorption lines. Methods. We obtained VLT/GRAVITY observations and combined them with earlier spatially resolved data covering different orbital phases. Results. We determine strong constraints on all orbital parameters and find that e = 0.322, close to the value derived from the emission lines. The X-ray light curve declines as s −3 after periastron, where s is the separation of the two stars, but its modulation is likely affected by absorption and occultation of the X-ray emitting region at other orbital phases. We find that previous discrepancies in the reddening value can be traced to a brighter K -band magnitude than predicted by the WR wind models. We derive E ( B − V ) = 0.02 ± 0.02 mag. Our more precise mass and radius values, combined with previously determined effective temperatures, provide strong constraints on evolutionary models. The closest match for the O star comes from a rotationally mixed model with initial mass M init = 28.7 M ⊙ and an M init = 32 M ⊙ model for the WR star, with negligible accretion onto the O star during the WR progenitor’s Roche lobe overflow phase. However, the temperature of the O star is higher and the mass of the WR star is smaller than predicted by the evolutionary tracks at the current epoch, consistent with the well-known “mass-discrepancy problem” in massive stars. Conclusions.γ 2 Velorum’s strongly constrained parameters make it ideal for testing structure, evolution, and stellar wind models.
- Open Access
- 2 June 2026
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 172 (1)
Elemental Abundances in the Binary Star V505 Per
We present a detailed chemical abundance analysis of the eclipsing binary system V505 Per. High-resolution spectra were analyzed using the MOOG spectrum analysis code, we determined abundances not only for iron and lithium but also for Si, Na, Ca, Mn, and Ni—elements that have not previously been analyzed in detail for this system. Abundances were computed across 15 temperature points using model atmospheres, with stellar parameters refined by minimizing abundance trends with excitation potential. We determined effective temperatures of T eff = 6650 ± 50 K for the primary and T eff = 6550 ± 50 K, with iron abundances of [Fe/H] = −0.10 ± 0.06 and [Fe/H] = −0.19 ± 0.07, respectively. Most [X/Fe] ratios are consistent with solar values, though manganese is deficient. Our analysis of the effective temperatures shows that both stars lie on the hot edge of the lithium dip, consistent with G. Koenigsberger et al., which may help resolve the inconsistency noted of the stars lithium abundance within the dip by Baugh et al.
- 2 June 2026
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 172 (1)
Open AccessElemental Abundances in the Binary Star V505 Per
We present a detailed chemical abundance analysis of the eclipsing binary system V505 Per. High-resolution spectra were analyzed using the MOOG spectrum analysis code, we determined abundances not only for iron and lithium but also for Si, Na, Ca, Mn, and Ni—elements that have not previously been analyzed in detail for this system. Abundances were computed across 15 temperature points using model atmospheres, with stellar parameters refined by minimizing abundance trends with excitation potential. We determined effective temperatures of T eff = 6650 ± 50 K for the primary and T eff = 6550 ± 50 K, with iron abundances of [Fe/H] = −0.10 ± 0.06 and [Fe/H] = −0.19 ± 0.07, respectively. Most [X/Fe] ratios are consistent with solar values, though manganese is deficient. Our analysis of the effective temperatures shows that both stars lie on the hot edge of the lithium dip, consistent with G. Koenigsberger et al., which may help resolve the inconsistency noted of the stars lithium abundance within the dip by Baugh et al.
- Open Access
Lithium in V505 Per
We determine the surface lithium abundance of the eclipsing binary components in V 505 Per (HIP 10961), A(Li)=2.65+/-0.07 and 2.35+/-0.1, which supports the rather unexpected conclusion that their surface Li abundances differ. We find effective temperatures 6600 K + 6550 K (~150 K higher than previously reported), which place the stars at the hot limit of the Lithium Dip, thus aleviating the previously suggested discrepancy with cluster stars of similar ages and temperatures. These temperatures are also more consistent with the system's Gaia spectral energy distribution. Our iron abundances, [Fe/H]=-0.15+\-0.07 and -0.25+\-0.1, agree with predictions of the higher temperatures deduced from our spectra and from evolutionary tracks. The surface rotation rate implied by our line profiles, 12.5+\-1 km/s, is smaller than the synchronous rotation rate, a curious result given the circular orbit and the age of the system.
Open AccessLithium in V505 Per
We determine the surface lithium abundance of the eclipsing binary components in V 505 Per (HIP 10961), A(Li)=2.65+/-0.07 and 2.35+/-0.1, which supports the rather unexpected conclusion that their surface Li abundances differ. We find effective temperatures 6600 K + 6550 K (~150 K higher than previously reported), which place the stars at the hot limit of the Lithium Dip, thus aleviating the previously suggested discrepancy with cluster stars of similar ages and temperatures. These temperatures are also more consistent with the system's Gaia spectral energy distribution. Our iron abundances, [Fe/H]=-0.15+\-0.07 and -0.25+\-0.1, agree with predictions of the higher temperatures deduced from our spectra and from evolutionary tracks. The surface rotation rate implied by our line profiles, 12.5+\-1 km/s, is smaller than the synchronous rotation rate, a curious result given the circular orbit and the age of the system.
- Open Access
- 7 January 2025
- Published by EDP Sciences in Astronomy & Astrophysics
- Vol. 693, A103
A Maia-type candidate was misclassified: V424 Cep is an eclipsing β Cep-type pulsator in a triple system
Context. V424 Cephei is an eclipsing binary system that has been classified as a Maia variable candidate. These objects are pulsators that apparently lie outside the theoretical instability strips. Aims. We determine the properties of V424 Cep and identify the nature of the pulsating variable. Methods. We analyzed photometric data obtained over the past three decades from TESS, Gaia, Hipparcos, and ground-based observations. Times of minimum light were determined, and the light curves and spectral energy distribution were analyzed. We analyzed the radial velocity curves of the double-lined system obtained from spectroscopy at the Observatorio Astronómico Nacional on the Sierra San Pedro Mártir. Results. The two eclipses observed in the light curves yield a refined orbital period P = 4.93 d, an eccentricity e = 0.02, and an apsidal period U = 730 yr. The eclipse O–C curve is not linear. Combined with the apsidal period, this indicates the presence of a third component in the system. This conclusion was confirmed by a comparison of predicted and observed absolute fluxes. The masses of the binary pair are 8.7 M⊙ and 6.3 M⊙, their radii are 6.0 R⊙ and 3.5 R⊙, and the luminosities are log(L/L⊙) 3.80 and 3.11. The TESS light curve shows oscillations with a dominant period of 0.17 d. They are coherent throughout the orbital cycles with a stable amplitude, except during the (partial) eclipse of the primary star. This indicates that this is the pulsator. Its temperature is 21 000 K. The third-light component contributes no more than 15% to the total light of the system in the TESS wavelength band and 20% in the K band. Conclusions. V424 Cep is not a Maia variable, but rather a β Cep star. This result highlights the importance of combining photometric light-curve solutions with absolute flux models and spectroscopic observations.
- 7 January 2025
- Published by EDP Sciences in Astronomy & Astrophysics
- Vol. 693, A103
Open AccessA Maia-type candidate was misclassified: V424 Cep is an eclipsing β Cep-type pulsator in a triple system
Context. V424 Cephei is an eclipsing binary system that has been classified as a Maia variable candidate. These objects are pulsators that apparently lie outside the theoretical instability strips. Aims. We determine the properties of V424 Cep and identify the nature of the pulsating variable. Methods. We analyzed photometric data obtained over the past three decades from TESS, Gaia, Hipparcos, and ground-based observations. Times of minimum light were determined, and the light curves and spectral energy distribution were analyzed. We analyzed the radial velocity curves of the double-lined system obtained from spectroscopy at the Observatorio Astronómico Nacional on the Sierra San Pedro Mártir. Results. The two eclipses observed in the light curves yield a refined orbital period P = 4.93 d, an eccentricity e = 0.02, and an apsidal period U = 730 yr. The eclipse O–C curve is not linear. Combined with the apsidal period, this indicates the presence of a third component in the system. This conclusion was confirmed by a comparison of predicted and observed absolute fluxes. The masses of the binary pair are 8.7 M⊙ and 6.3 M⊙, their radii are 6.0 R⊙ and 3.5 R⊙, and the luminosities are log(L/L⊙) 3.80 and 3.11. The TESS light curve shows oscillations with a dominant period of 0.17 d. They are coherent throughout the orbital cycles with a stable amplitude, except during the (partial) eclipse of the primary star. This indicates that this is the pulsator. Its temperature is 21 000 K. The third-light component contributes no more than 15% to the total light of the system in the TESS wavelength band and 20% in the K band. Conclusions. V424 Cep is not a Maia variable, but rather a β Cep star. This result highlights the importance of combining photometric light-curve solutions with absolute flux models and spectroscopic observations.
- Open Access
- 1 January 2025
- Published by Cambridge University Press (CUP) in Proceedings of the International Astronomical Union
- Vol. 21 (S402) , 224-230
Concluding remarks
The International Astronomical Union Symposium Number 402 focused on massive stars across redshifts in the era of JWST and large-scale surveys, and brought together specialists from a broad range of topics. I briefly summarize a few of the highlights from the excellent presentations and rich discussions and, because recent observations suggest that binary interactions dominated star formation and evolution during the early universe, I venture a few personal reflections on this topic.
- 1 January 2025
- Published by Cambridge University Press (CUP) in Proceedings of the International Astronomical Union
- Vol. 21 (S402) , 224-230
Open AccessConcluding remarks
The International Astronomical Union Symposium Number 402 focused on massive stars across redshifts in the era of JWST and large-scale surveys, and brought together specialists from a broad range of topics. I briefly summarize a few of the highlights from the excellent presentations and rich discussions and, because recent observations suggest that binary interactions dominated star formation and evolution during the early universe, I venture a few personal reflections on this topic.
- Open Access
- 1 January 2025
- Published by Cambridge University Press (CUP) in Proceedings of the International Astronomical Union
- Vol. 21 (S402) , 290-292
Tidal shear energy dissipation and stellar structure
Asynchronous rotation in binary stars induces a localized differential rotation structure which, in the presence of viscous turbulence, can inject energy into the outer stellar layers. This paper explores the effect of such energy injection in a 162 upper M Subscript circled dot $M_\odot$ star early in its evolutionary lifetime. We find that in response to the injected energy, the radius, luminosity and mass-loss rate increase, compared to a standard model. We also find that the computed turbulent viscosity values imply synchronization timescales for the inner layers that exceed the main sequence lifetime. Hence, binary stars on the main sequence may not be in true corotation even if their outer layers appear to be rotating synchronously, and they may appear more luminous and have larger mass-loss rates than their analogous single-star counterparts.- 1 January 2025
- Published by Cambridge University Press (CUP) in Proceedings of the International Astronomical Union
- Vol. 21 (S402) , 290-292
Open AccessTidal shear energy dissipation and stellar structure
Asynchronous rotation in binary stars induces a localized differential rotation structure which, in the presence of viscous turbulence, can inject energy into the outer stellar layers. This paper explores the effect of such energy injection in a 162 upper M Subscript circled dot $M_\odot$ star early in its evolutionary lifetime. We find that in response to the injected energy, the radius, luminosity and mass-loss rate increase, compared to a standard model. We also find that the computed turbulent viscosity values imply synchronization timescales for the inner layers that exceed the main sequence lifetime. Hence, binary stars on the main sequence may not be in true corotation even if their outer layers appear to be rotating synchronously, and they may appear more luminous and have larger mass-loss rates than their analogous single-star counterparts. - Open Access
- 1 January 2025
- Published by Cambridge University Press (CUP) in Proceedings of the International Astronomical Union
- Vol. 21 (S402) , 286-289
The low-metallicity binary black hole progenitor prototype HD 5980
The multiple Small Magellanic Cloud system HD 5980 contains two very massive stars displaying a Wolf-Rayet type spectrum and which appear to be following chemically homogeneous evolutionary tracks. This means that they will remain very massive and compact prior to ending their lives as supernovae, and are thus excellent binary black hole progenitor prototypes formed in low-metallicty environments. In this paper we summarize some of HD 5980’s extraordinary properties and refine its orbital period for the 2018–2025 timespan, based on the times of secondary minima in TESS observations.
- 1 January 2025
- Published by Cambridge University Press (CUP) in Proceedings of the International Astronomical Union
- Vol. 21 (S402) , 286-289
Open AccessThe low-metallicity binary black hole progenitor prototype HD 5980
The multiple Small Magellanic Cloud system HD 5980 contains two very massive stars displaying a Wolf-Rayet type spectrum and which appear to be following chemically homogeneous evolutionary tracks. This means that they will remain very massive and compact prior to ending their lives as supernovae, and are thus excellent binary black hole progenitor prototypes formed in low-metallicty environments. In this paper we summarize some of HD 5980’s extraordinary properties and refine its orbital period for the 2018–2025 timespan, based on the times of secondary minima in TESS observations.
- Open Access
- 17 May 2024
- Published by EDP Sciences in Astronomy & Astrophysics
- Vol. 685, A145
Eccentric binaries: Periastron events and tidal heating
Context. Binary stars cannot be fully understood without assessing the interaction effects between the two components and the impact of these effects on observational diagnostics. Periastron brightening events, also known as the heartbeat phenomenon, are a clear manifestation of this type of interaction. Aims. We aim to explore the role of tidal shear energy dissipation in stars undergoing periastron brightening events. Methods. We performed a computation from first principles that uses a quasi-hydrodynamic Lagrangian scheme to simultaneously solve the orbital motion and the equations of motion of a 3D grid of volume elements covering the inner, rigidly rotating “core” of a tidally perturbed star. The equations of motion include the gravitational acceleration of both stars, the centrifugal, Coriolis, gas pressure accelerations, and viscous coupling between volume elements. The method is illustrated for a grid of model binary systems with a 10 M⊙ primary that is perturbed by a 6.97 M⊙ companion in eccentric orbits (e = 0 − 0.7). The model is then applied to the heartbeat star MACHO 80.7443.1718. Results. We find an increase by factors ∼10−6–10−3 in tidal shear energy dissipation at periastron, consistent with the majority of observed heartbeat stars. The magnitude of the periastron effect correlates with the degree of departure from synchronicity: stars rotating much faster or much slower than the synchronous rate at periastron present the strongest effect. We confirm that for eccentricities ≤0.3, pseudo-synchronization occurs for 0.8 < ω/Ωave < 1, where Ωave is the average orbital angular velocity. The minimum energy rotation rate (pseudo-synchronism) for e = 0.5 and 0.7 occurs for 1.0 < ω/Ωave < 1.15. The tidal shear energy dissipation model reproduces from first principles the ∼23% maximum brightness enhancement at periastron of MACHO 80.7443.1718. Conclusions. Our results suggest that the magnitude and shape of the heartbeat signal may serve as diagnostics for the internal stellar rotation and turbulent viscosity values.
- 17 May 2024
- Published by EDP Sciences in Astronomy & Astrophysics
- Vol. 685, A145
Open AccessEccentric binaries: Periastron events and tidal heating
Context. Binary stars cannot be fully understood without assessing the interaction effects between the two components and the impact of these effects on observational diagnostics. Periastron brightening events, also known as the heartbeat phenomenon, are a clear manifestation of this type of interaction. Aims. We aim to explore the role of tidal shear energy dissipation in stars undergoing periastron brightening events. Methods. We performed a computation from first principles that uses a quasi-hydrodynamic Lagrangian scheme to simultaneously solve the orbital motion and the equations of motion of a 3D grid of volume elements covering the inner, rigidly rotating “core” of a tidally perturbed star. The equations of motion include the gravitational acceleration of both stars, the centrifugal, Coriolis, gas pressure accelerations, and viscous coupling between volume elements. The method is illustrated for a grid of model binary systems with a 10 M⊙ primary that is perturbed by a 6.97 M⊙ companion in eccentric orbits (e = 0 − 0.7). The model is then applied to the heartbeat star MACHO 80.7443.1718. Results. We find an increase by factors ∼10−6–10−3 in tidal shear energy dissipation at periastron, consistent with the majority of observed heartbeat stars. The magnitude of the periastron effect correlates with the degree of departure from synchronicity: stars rotating much faster or much slower than the synchronous rate at periastron present the strongest effect. We confirm that for eccentricities ≤0.3, pseudo-synchronization occurs for 0.8 < ω/Ωave < 1, where Ωave is the average orbital angular velocity. The minimum energy rotation rate (pseudo-synchronism) for e = 0.5 and 0.7 occurs for 1.0 < ω/Ωave < 1.15. The tidal shear energy dissipation model reproduces from first principles the ∼23% maximum brightness enhancement at periastron of MACHO 80.7443.1718. Conclusions. Our results suggest that the magnitude and shape of the heartbeat signal may serve as diagnostics for the internal stellar rotation and turbulent viscosity values.
- Open Access
- 1 May 2024
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 967 (1) , 69
Transitory Tidal Heating and Its Impact on Cluster Isochrones
The kinetic energy in tidal flows, when converted into heat, can affect the internal structure of a star and shift its location on a color–magnitude diagram from that of standard models. In this paper we explore the impact of injecting heat into stars with masses near the main-sequence turnoff mass (1.26 M ⊙) of the open cluster M67. The heating rate is obtained from the tidal shear energy dissipation rate, which is calculated from first principles by simultaneously solving the equations that describe orbital motion and the response of a star’s layers to the gravitational, Coriolis, centrifugal, gas pressure, and viscous forces. The stellar structure models are computed with MESA. We focus on the effects of injecting heat in pulses lasting 0.01 Gyr, a time frame consistent with the synchronization timescale in binary systems. We find that the location of the tidally perturbed stars in the M67 color–magnitude diagram is shifted to significantly higher luminosities and effective temperatures than predicted by the standard model isochrone and includes locations corresponding to some of the blue straggler stars. Because tidal heating takes energy from the orbit, causing it to shrink, blue straggler stars could be merger or mass transfer progenitors, as well as products of these processes.
- 1 May 2024
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 967 (1) , 69
Open AccessTransitory Tidal Heating and Its Impact on Cluster Isochrones
The kinetic energy in tidal flows, when converted into heat, can affect the internal structure of a star and shift its location on a color–magnitude diagram from that of standard models. In this paper we explore the impact of injecting heat into stars with masses near the main-sequence turnoff mass (1.26 M ⊙) of the open cluster M67. The heating rate is obtained from the tidal shear energy dissipation rate, which is calculated from first principles by simultaneously solving the equations that describe orbital motion and the response of a star’s layers to the gravitational, Coriolis, centrifugal, gas pressure, and viscous forces. The stellar structure models are computed with MESA. We focus on the effects of injecting heat in pulses lasting 0.01 Gyr, a time frame consistent with the synchronization timescale in binary systems. We find that the location of the tidally perturbed stars in the M67 color–magnitude diagram is shifted to significantly higher luminosities and effective temperatures than predicted by the standard model isochrone and includes locations corresponding to some of the blue straggler stars. Because tidal heating takes energy from the orbit, causing it to shrink, blue straggler stars could be merger or mass transfer progenitors, as well as products of these processes.
- Open Access
- 22 April 2024
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 966 (1) , 9
Investigating the Chemically Homogeneous Evolution Channel and Its Role in the Formation of the Enigmatic Binary Black Hole Progenitor Candidate HD 5980
Chemically homogeneous evolution (CHE) is a promising channel for forming massive binary black holes. The enigmatic, massive Wolf–Rayet binary HD 5980 A&B has been proposed to have formed through this channel. We investigate this claim by comparing its observed parameters with CHE models. Using MESA, we simulate grids of close massive binaries, then use a Bayesian approach to compare them with the stars’ observed orbital period, masses, luminosities, and hydrogen surface abundances. The most probable models, given the observational data, have initial periods ∼3 days, widening to the present-day ∼20 days orbit as a result of mass loss—correspondingly, they have very high initial stellar masses (≳150 M ⊙). We explore variations in stellar-wind mass loss and internal mixing efficiency, and find that models assuming enhanced mass loss are greatly favored to explain HD 5980, while enhanced mixing is only slightly favored over our fiducial assumptions. Our most probable models slightly underpredict the hydrogen surface abundances. Regardless of its prior history, this system is a likely binary black hole progenitor. We model its further evolution under our fiducial and enhanced wind assumptions, finding that both stars produce black holes with masses ∼19–37 M ⊙. The projected final orbit is too wide to merge within a Hubble time through gravitational waves alone. However, the system is thought to be part of a 2+2 hierarchical multiple. We speculate that secular effects with the (possible) third and fourth companions may drive the system to promptly become a gravitational-wave source.
- 22 April 2024
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 966 (1) , 9
Open AccessInvestigating the Chemically Homogeneous Evolution Channel and Its Role in the Formation of the Enigmatic Binary Black Hole Progenitor Candidate HD 5980
Chemically homogeneous evolution (CHE) is a promising channel for forming massive binary black holes. The enigmatic, massive Wolf–Rayet binary HD 5980 A&B has been proposed to have formed through this channel. We investigate this claim by comparing its observed parameters with CHE models. Using MESA, we simulate grids of close massive binaries, then use a Bayesian approach to compare them with the stars’ observed orbital period, masses, luminosities, and hydrogen surface abundances. The most probable models, given the observational data, have initial periods ∼3 days, widening to the present-day ∼20 days orbit as a result of mass loss—correspondingly, they have very high initial stellar masses (≳150 M ⊙). We explore variations in stellar-wind mass loss and internal mixing efficiency, and find that models assuming enhanced mass loss are greatly favored to explain HD 5980, while enhanced mixing is only slightly favored over our fiducial assumptions. Our most probable models slightly underpredict the hydrogen surface abundances. Regardless of its prior history, this system is a likely binary black hole progenitor. We model its further evolution under our fiducial and enhanced wind assumptions, finding that both stars produce black holes with masses ∼19–37 M ⊙. The projected final orbit is too wide to merge within a Hubble time through gravitational waves alone. However, the system is thought to be part of a 2+2 hierarchical multiple. We speculate that secular effects with the (possible) third and fourth companions may drive the system to promptly become a gravitational-wave source.