New self-consistent theoretical descriptions for mass-loss rates of O-type stars

TL;DR

This paper presents a self-consistent theoretical framework for calculating mass-loss rates of O-type stars, incorporating detailed line acceleration and stellar atmosphere models, leading to improved predictions aligned with observations.

Contribution

It introduces a novel self-consistent approach coupling hydrodynamics with detailed line-force calculations for O-type star winds, considering multiple atomic configurations.

Findings

Mass-loss rates decrease with more elements in the radiation field.

Theoretical prescriptions for mass-loss rates show high correlation with models (r > 0.92).

Results align with observed wind momentum-luminosity relationships.

Abstract

Massive O-type stars lose a significant fraction of their mass through radiation-driven winds, a process that critically shapes their evolution and feedback into the interstellar medium. Accurate predictions of mass-loss rates are essential for models of stellar structure and population synthesis. We computed wind parameters for O-type stars using a self-consistent approach that couples the hydrodynamics of the wind with detailed calculations of the line acceleration. This approach follows the theory of radiation-driven stellar winds and allows us to derive mass-loss rate distributions for different atomic configurations of the stellar flux. We used the TLUSTY code for stellar atmosphere models to compute non-local thermodynamic equilibrium models; these models served as input radiation fields for the calculation of the line-force parameters, for which we used the LOCUS code. These line-force parameters were then iteratively coupled with the HYDWIND code to solve the wind hydrodynamics. The procedure was applied across a grid of stellar parameters for three chemical configurations. We obtain self-consistent wind parameters for a broad set of O-type stellar models. The results show a systematic decrease in mass-loss rates with the inclusion of more elements in the radiation field, which is attributed to a strong effect on the UV region of the spectral energy distribution. As more elements are included, resulting in a larger number of spectral lines, the contribution from the UV diminishes, leading to lower mass-loss rates. We fitted three theoretical prescriptions for $\dot{M}$ using a Bayesian approach; this yielded Pearson correlation values greater than 0.92 for all three model grids. It also allowed for the estimation of the wind momentum-luminosity relationships for each of the grids, yielding results similar to those based on observations of O-type stars.