No Arabic abstract
The standard, or fast, solutions of m-CAK line-driven wind theory cannot account for slowly outflowing disks like the ones that surround Be stars. It has been previously shown that there exists another family of solutions --- the $Omega$-slow solutions --- that is characterized by much slower terminal velocities and higher mass-loss rates. We have solved the one-dimensional m-CAK hydrodynamical equation of rotating radiation-driven winds for this latter solution, starting from standard values of the line force parameters ($alpha$, $k$, and $delta$), and then systematically varying the values of $alpha$ and $k$. Terminal velocities and mass-loss rates that are in good agreement with those found in Be stars are obtained from the solutions with lower $alpha$ and higher $k$ values. Furthermore, the equatorial densities of such solutions are comparable to those that are typically assumed in ad hoc models. For very high values of $k$, we find that the wind solutions exhibit a new kind of behavior.
As the disk formation mechanism(s) in Be stars is(are) as yet unknown, we investigate the role of rapidly rotating radiation-driven winds in this process. We implemented the effects of high stellar rotation on m-CAK models accounting for: the shape of the star, the oblate finite disk correction factor, and gravity darkening. For a fast rotating star, we obtain a two-component wind model, i.e., a fast, thin wind in the polar latitudes and an $Omega$-slow, dense wind in the equatorial regions. We use the equatorial mass densities to explore H$alpha$ emission profiles for the following scenarios: 1) a spherically symmetric star, 2) an oblate shaped star with constant temperature, and 3) an oblate star with gravity darkening. One result of this work is that we have developed a novel method for solving the gravity darkened, oblated m-CAK equation of motion. Furthermore, from our modeling we find a) the oblate finite disk correction factor, for the scenario considering the gravity darkening, can vary by at least a factor of two between the equatorial and polar directions, influencing the velocity profile and mass-loss rate accordingly, b) the H$alpha$ profiles predicted by our model are in agreement with those predicted by a standard power-law model for following values of the line-force parameters: $1.5 lesssim k lesssim 3$, $ , alpha sim 0.6$ and $, delta gtrsim 0.1$, and c) the contribution of the fast wind component to the H$alpha$ emission line profile is negligible; therefore, the line profiles arise mainly from the equatorial disks of Be stars.
Accurate mass-loss rates and terminal velocities from massive stars winds are essential to obtain synthetic spectra from radiative transfer calculations and to determine the evolutionary path of massive stars. From a theoretical point of view, analytical expressions for the wind parameters and velocity profile would have many advantages over numerical calculations that solve the complex non-linear set of hydrodynamic equations. In a previous work, we obtained an analytical description for the fast wind regime. Now, we propose an approximate expression for the line-force in terms of new parameters and obtain a velocity profile closed-form solution (in terms of the Lambert $W$ function) for the $delta$-slow regime. Using this analytical velocity profile, we were able to obtain the mass-loss rates based on the m-CAK theory. Moreover, we established a relation between this new set of line-force parameters with the known stellar and m-CAK line-force parameters. To this purpose, we calculated a grid of numerical hydrodynamical models and performed a multivariate multiple regression. The numerical and our descriptions lead to good agreement between their values.
Massive stars present strong stellar that which are described by the radiation driven wind theory. Accurate mass-loss rates are necessary to properly describe the stellar evolution across the Hertzsprung--Russel Diagram. We present a self-consistent procedure that coupled the hydrodynamics with calculations of the line-force, giving as results the line-force parameters, the velocity field, and the mass-loss rate. Our calculations contemplate the contribution to the line-force multiplier from more than $sim 900,000$ atomic transitions, an NLTE radiation flux from the photosphere and a quasi-LTE approximation for the occupational numbers. A full set of line-force parameters for $T_text{eff}ge 32,000$ K, surface gravities higher than 3.4 dex for two different metallicities are presented, with their corresponding wind parameters (terminal velocities and mass-loss rates). The already known dependence of line-force parameters on effective temperature is enhanced by the dependence on $log g$. The terminal velocities present a stepper scaling relation with respect to the escape velocity, this might explain the scatter values observed in the hot side of the bistability jump. Moreover, a comparison of self-consistent mass-loss rates with empirical values shows a good agreement. Self-consistent wind solutions are used as input in FASTWIND to calculate synthetic spectra. We show, comparing with the observed spectra for three stars, that varying the clumping factor, the synthetic spectra rapidly converge into the neighbourhood region of the solution. It is important to stress that our self-consistent procedure significantly reduces the number of free parameters needed to obtain a synthetic spectrum.
Hot massive stars present strong stellar winds that are driven by absorption, scattering and re-emission of photons by the ions of the atmosphere (textit{line-driven winds}). A better comprehension of this phenomenon, and a more accurate calculation of hydrodynamics and radiative acceleration is required to reduce the number of free parameters in spectral fitting, to determine accurate wind parameters such as mass-loss rates and velocity profiles. We use the non-LTE model-atmosphere code CMFGEN to numerically solve the radiative transfer equation in the stellar atmosphere and to calculate the radiative acceleration $g_text{rad}(r)$. Under the assumption that the radiative acceleration depends only on the radial coordinate, we solve analytically the equation of motion by means of the Lambert $W$-function. An iterative procedure between the solution of the radiative transfer and the equation of motion is executed in order to obtain a final self-consistent velocity field that is no longer based on any $beta$-law. We apply the Lambert-procedure to three O supergiant stars ($zeta$-Puppis, HD~165763 and $alpha$-Cam) and discuss the Lambert-solutions for the velocity profiles. It is found that, even without recalculation of the mass-loss rate, the Lambert-procedure allows the calculation of consistent velocity profiles that reduce the number of free parameters when a spectral fitting using CMFGEN is performed. Synthetic spectra calculated from our Lambert-solutions show significant differences compared to the initial $beta$-law CMFGEN models. The results indicate the importance of consistent velocity profile calculation in the CMFGEN code and its usage in a fitting procedure and interpretation of observed spectra.
We present two self-consistent procedures that couple the hydrodynamics with calculations of the line-force in the frame of radiation wind theory. These procedures give us the line-force parameters, the velocity field, and the mass-loss rate. The first one is based on the so-called m-CAK theory. A full set of line-force parameters for $T_text{eff}ge 32,000$ K and surface gravities higher than 3.4 dex for two different metallicities are presented, along with their corresponding wind parameters. We find that the dependence of line-force parameters on effective temperature is enhanced by the dependence on $log g$. For the case of homogeneous winds (without clumping) comparison of self-consistent mass-loss rates shows a good agreement with empirical values. We also consider self-consistent wind solutions that are used as input in FASTWIND to calculate synthetic spectra. By comparison with the observed spectra for three stars with clumped winds, we found that varying the clumping factor the synthetic spectra rapidly converge into the neighbourhood region of the solution. Therefore, this self-consistent m-CAK procedure significantly reduces the number of free parameters needed to obtain a synthetic spectrum. The second procedure (called Lambert-procedure) provides a self-consistent solution beyond m-CAK theory, and line-acceleration is calculated by the full NLTE radiative transfer code CMFGEN. Both the mass-loss rate and the clumping factor are set as free parameters, hence their values are obtained by spectral fitting after the respective self-consistent hydrodynamics is calculated. Since performing the Lambert-procedure requires significant computational power, the analysis is made only for the star z-Puppis. The promising results gives a positive balance about the future applications for the self-consistent solutions presented on this thesis.