Spectral modelling of massive binary systems

TL;DR

This paper models the spectra of massive binary star systems across their orbital phases, incorporating gravitational, radiative, and dynamical effects to better understand spectral variations and phenomena like the Struve-Sahade effect.

Contribution

It introduces a comprehensive spectral modeling approach that accounts for radiation pressure, gravity darkening, and mutual irradiation, improving the realism of binary star spectra simulations.

Findings

Models predict line strength variations over orbital cycles.

Including radiation pressure aligns surface temperature distributions with observations.

Classical von Zeipel theorem suffices for spectral calculations.

Abstract

Aims: We simulate the spectra of massive binaries at different phases of the orbital cycle, accounting for the gravitational influence of the companion star on the shape and physical properties of the stellar surface. Methods: We used the Roche potential modified to account for radiation pressure to compute the stellar surface of close circular systems and we used the TIDES code for surface computation of eccentric systems. In both cases, we accounted for gravity darkening and mutual heating generated by irradiation to compute the surface temperature. We then interpolated NLTE plane-parallel atmosphere model spectra in a grid to obtain the local spectrum at each surface point. We finally summed all contributions, accounting for the Doppler shift, limb-darkening, and visibility to obtain the total synthetic spectrum. We computed different orbital phases and sets of physical and orbital parameters. Results: Our models predict line strength variations through the orbital cycle, but fail to completely reproduce the Struve-Sahade effect. Including radiation pressure allows us to reproduce a surface temperature distribution that is consistent with observations of semi-detached binary systems. Conclusions: Radiation pressure effects on the stellar surface are weak in (over)contact binaries and well-detached systems but can become very significant in semi-detached systems. The classical von Zeipel theorem is sufficient for the spectral computation. Broad-band light curves derived from the spectral computation are different from those computed with a model in which the stellar surfaces are equipotentials of the Roche potential scaled by the instantaneous orbital separation. In many cases, the fit of two Gaussian/Lorentzian profiles fails to properly measure the equivalent width of the lines and leads to apparent variations that could explain some of the effects reported in the literature.