Modeling contact binaries, II. The effect of energy transfer

Matthias Fabry,Pablo Marchant,Norbert Langer,Hugues Sana
DOI: https://doi.org/10.1051/0004-6361/202346277
2023-03-16
Abstract:Context. It is common for massive stars to engage in binary interaction. In close binaries, the components can enter a contact phase, where both stars overflow their respective Roche lobes simultaneously. While there exist observational constraints on the stellar properties of such systems, the most detailed stellar evolution models that feature a contact phase are not fully reconcilable with those measurements. Aims. We aim to consistently model contact phases of binary stars in a 1D stellar evolution code. To this end, we develop the methodology to account for energy transfer in the common contact layers. Methods. We implement an approximative model for energy transfer between the components of a contact binary based on the von Zeipel theorem in the stellar evolution code MESA. We compare structure and evolution models with and without this transfer and analyze the implications for the observable properties of the contact phase. Results. Implementing energy transfer helps eliminating baroclinicity in the common envelope between the components of a contact binary, which, if present, would drive strong thermal flows. We find that accounting for energy transfer in massive contact binaries significantly alters the mass ratio evolution and can extend the lifetime of an unequal mass ratio contact system.
Solar and Stellar Astrophysics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the inconsistency between theory and observation in the modeling process of contact binaries in binary star systems. Specifically, the paper focuses on the influence of energy transfer (ET) in contact binaries, especially when the mass ratios of the two stars are not equal. Existing stellar evolution models fail to fully explain the observed stellar characteristics when dealing with contact binaries, especially when these stars are in the contact phase. Therefore, the authors developed a new method, introducing an approximate model of energy transfer in the one - dimensional stellar evolution code, in order to more accurately simulate the structure and evolution process of contact binaries. ### Paper Background - **Contact Binaries**: In a close binary star system, when two stars simultaneously overflow their respective Roche lobes, they will enter a contact configuration, forming so - called contact binaries. - **Inconsistency between Observation and Model**: Although there are some observational constraints, the most detailed stellar evolution models still cannot fully match the observational data when dealing with contact binaries. - **Research Objective**: The paper aims to consistently simulate the contact phase of contact binaries by introducing an energy transfer model in the one - dimensional stellar evolution code, and analyze the influence of this energy transfer on the observable properties of contact binaries. ### Main Research Contents 1. **Implementation of the Energy Transfer Model**: - The authors implemented an energy transfer model based on the von Zeipel theorem in the stellar evolution code MESA. - The model takes into account the energy transfer in the contact layer, which helps to eliminate baroclinicity in the common envelope, thereby reducing the strong heat flux. 2. **Model Comparison**: - The authors compared the structural and evolution models with and without energy transfer. - The results show that considering energy transfer significantly changes the evolution of the mass ratio and can extend the lifetime of contact systems with unequal mass ratios. 3. **Comparison between Theory and Observation**: - By introducing the energy transfer model, the authors found that the properties of observed contact binaries can be better explained, especially when the mass ratio is not equal. ### Key Conclusions - **Importance of Energy Transfer**: Introducing an energy transfer model in contact binaries can eliminate baroclinicity in the common envelope, reduce heat flux, and thus more accurately simulate the evolution of contact binaries. - **Mass Ratio Evolution**: Considering energy transfer significantly changes the evolution of the mass ratio, especially for contact systems with unequal mass ratios, and can extend their lifetimes. - **Observation Consistency**: The new model better explains the observational data, especially in contact binaries with unequal mass ratios. ### Formula Summary - **Mass - Radius Relationship**: \[ \frac{R_2}{R_1} \bigg|_{\text{ZAMS}} = \left(\frac{M_2}{M_1}\right)^{0.57} = q^{0.57} \] \[ \frac{R_2}{R_1} \bigg|_{\text{Roche}} \approx \left(\frac{M_2}{M_1}\right)^{0.46} = q^{0.46} \] - **Luminosity Ratio**: \[ \frac{L_2}{L_1} \approx \frac{M_2}{M_1} = q \] - **Amount of Energy Transfer**: \[ f \approx \frac{q - q^\alpha}{1 + q^\alpha} \] \[ L_{\text{trans}} = L_1 \frac{S_2 \langle g \rangle_2 - L_2 S_1 \langle g \rangle_1}{S_1 \langle g \rangle_1 + S_2 \langle g \rangle_2} \] Through these models and analyses, the paper provides new perspectives and tools for understanding the evolution of contact binaries.