Refine
Year of publication
Language
- English (35)
Is part of the Bibliography
- yes (35)
Keywords
- stars: mass-loss (35) (remove)
Institute
We investigate the connections between the magnetic fields and the X-ray emission from massive stars. Our study shows that the X-ray properties of known strongly magnetic stars are diverse: while some comply to the predictions of the magnetically confined wind model, others do not. We conclude that strong, hard, and variable X-ray emission may be a sufficient attribute of magnetic massive stars, but it is not a necessary one. We address the general properties of X-ray emission from "normal" massive stars, especially the long standing mystery about the correlations between the parameters of X-ray emission and fundamental stellar properties. The recent development in stellar structure modeling shows that small-scale surface magnetic fields may be common. We suggest a "hybrid" scenario which could explain the X-ray emission from massive stars by a combination of magnetic mechanisms on the surface and shocks in the stellar wind. The magnetic mechanisms and the wind shocks are triggered by convective motions in sub-photospheric layers. This scenario opens the door for a natural explanation of the well established correlation between bolometric and X-ray luminosities.
The star zeta Ophiuchi is one of the brightest massive stars in the northern hemisphere and was intensively studied in various wavelength domains. The currently available observational material suggests that certain observed phenomena are related to the presence of a magnetic field. We acquired spectropolarimetric observations of zeta Oph with FORS 1 mounted on the 8-m Kueyen telescope of the VLT to investigate if a magnetic field is indeed present in this star. Using all available absorption lines, we detect a mean longitudinal magnetic field < B(z)>(all) = 141 +/- 45 G, confirming the magnetic nature of this star. We review the X-ray properties of zeta Oph with the aim to understand whether the X-ray emission of zeta Oph is dominated by magnetic or by wind instability processes.
The Galactic WC stars Stellar parameters from spectral analyses indicate a new evolutionary sequence
(2012)
Context. The life cycles of massive stars from the main sequence to their explosion as supernovae or gamma ray bursts are not yet fully clear, and the empirical results from spectral analyses are partly in conflict with current evolutionary models. The spectral analysis of Wolf-Rayet stars requires the detailed modeling of expanding stellar atmospheres in non-LTE. The Galactic WN stars have been comprehensively analyzed with such models of the latest stage of sophistication, while a similarly comprehensive study of the Galactic WC sample remains undone.
Aims. We aim to establish the stellar parameters and mass-loss rates of the Galactic WC stars. These data provide the empirical basis of studies of (i) the role of WC stars in the evolution of massive stars, (ii) the wind-driving mechanisms, and (iii) the feedback of WC stars as input to models of the chemical and dynamical evolution of galaxies.
Methods. We analyze the nearly complete sample of un-obscured Galactic WC stars, using optical spectra as well as ultraviolet spectra when available. The observations are fitted with theoretical spectra, using the Potsdam Wolf-Rayet (PoWR) model atmosphere code. A large grid of line-blanked models has been established for the range of WC subtypes WC4 - WC8, and smaller grids for the WC9 parameter domain. Both WO stars and WN/WC transit types are also analyzed using special models.
Results. Stellar and atmospheric parameters are derived for more than 50 Galactic WC and two WO stars, covering almost the whole Galactic WC population as far as the stars are single, and un-obscured in the visual. In the Hertzsprung-Russell diagram, the WC stars reside between the hydrogen and the helium zero-age main sequences, having luminosities L from 10(4.9) to 10(5.6) L-circle dot. The mass-loss rates scale very tightly with L-0.8. The two WO stars in our sample turn out to be outstandingly hot (approximate to 200 kK) and do not fit into the WC scheme.
Conclusions. By comparing the empirical WC positions in the Hertzsprung-Russell diagram with evolutionary models, and from recent supernova statistics, we conclude that WC stars have evolved from initial masses between 20 solar masses and 45 M-circle dot. In contrast to previous assumptions, it seems that WC stars in general do not descend from the most massive stars. Only the WO stars might stem from progenitors that have been initially more massive than 45 M-circle dot.
On the weak-wind problem in massive stars X-ray spectra reveal a massive hot wind in mu columbaea
(2012)
mu Columbae is a prototypical weak-wind O star for which we have obtained a high-resolution X-ray spectrum with the Chandra LETG/ACIS instrument and a low-resolution spectrum with Suzaku. This allows us, for the first time, to investigate the role of X-rays on the wind structure in a bona fide weak-wind system and to determine whether there actually is a massive hot wind. The X-ray emission measure indicates that the outflow is an order of magnitude greater than that derived from UV lines and is commensurate with the nominal wind-luminosity relationship for O stars. Therefore, the "weak-wind problem"-identified from cool wind UV/optical spectra-is largely resolved by accounting for the hot wind seen in X-rays. From X-ray line profiles, Doppler shifts, and relative strengths, we find that this weak-wind star is typical of other late O dwarfs. The X-ray spectra do not suggest a magnetically confined plasma-the spectrum is soft and lines are broadened; Suzaku spectra confirm the lack of emission above 2 keV. Nor do the relative line shifts and widths suggest any wind decoupling by ions. The He-like triplets indicate that the bulk of the X-ray emission is formed rather close to the star, within five stellar radii. Our results challenge the idea that some OB stars are "weak-wind" stars that deviate from the standard wind-luminosity relationship. The wind is not weak, but it is hot and its bulk is only detectable in X-rays.
Recent studies have claimed the existence of very massive stars (VMS) up to 300 M⊙ in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 M⊙, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.
Context. The true mass-loss rates from massive stars are important for many branches of astrophysics. For the correct modeling of the resonance lines, which are among the key diagnostics of stellar mass-loss, the stellar wind clumping has been found to be very important. To incorporate clumping into a radiative transfer calculation, three-dimensional (3D) models are required. Various properties of the clumps may have a strong impact on the resonance line formation and, therefore, on the determination of empirical mass-loss rates.
Aims. We incorporate the 3D nature of the stellar wind clumping into radiative transfer calculations and investigate how different model parameters influence the resonance line formation.
Methods. We develop a full 3D Monte Carlo radiative transfer code for inhomogeneous expanding stellar winds. The number density of clumps follows the mass conservation. For the first time, we use realistic 3D models that describe the dense as well as the tenuous wind components to model the formation of resonance lines in a clumped stellar wind. At the same time, we account for non-monotonic velocity fields.
Results. The 3D density and velocity wind inhomogeneities show that there is a very strong impact on the resonance line formation. The different parameters describing the clumping and the velocity field results in different line strengths and profiles. We present a set of representative models for various sets of model parameters and investigate how the resonance lines are affected. Our 3D models show that the line opacity is lower for a larger clump separation and shallower velocity gradients within the clumps.
Conclusions. Our model demonstrates that to obtain empirically correct mass-loss rates from the UV resonance lines, the wind clumping and its 3D nature must be taken into account.
A considerable fraction of the central stars of planetary nebulae (CSPNe) are hydrogen-deficient. Almost all of these H-deficient central stars (CSs) display spectra with strong carbon and helium lines. Most of them exhibit emission-line spectra resembling those of massive WC stars. Therefore these stars are classed as CSPNe of spectral type [WC]. Recently, quantitative spectral analysis of two emission-line CSs, PB 8 and IC 4663, revealed that these stars do not belong to the [WC] class. Instead PB 8 has been classified as [WN/WC] type and IC 4663 as [WN] type. In this work we report the spectroscopic identification of another rare [WN] star, the CS of Abell 48. We performed a spectral analysis of Abell 48 with the Potsdam Wolf-Rayet (PoWR) models for expanding atmospheres. We find that the expanding atmosphere of Abell 48 is mainly composed of helium (85 per cent by mass), hydrogen (10 per cent) and nitrogen (5 per cent). The residual hydrogen and the enhanced nitrogen abundance make this object different from the other [WN] star IC 4663. We discuss the possible origin of this atmospheric composition.
Context. Recent studies of O-type stars have demonstrated that discrepant mass-loss rates are obtained when different diagnostic methods are employed. Fitting the unsaturated UV resonance lines (e.g., P v) gives drastically lower values than obtained from the Ha emission. Wind inhomogeneity (so-called "clumping") may be the main cause of this discrepancy. Aims. In a previous paper, we presented 3D Monte-Carlo calculations for the formation of scattering lines in a clumped stellar wind. In the present paper we select five O-type supergiants (from 04 to 07) and test whether the reported discrepancies can be resolved this way. Methods. In the first step, the analyses started with simulating the observed spectra with Potsdam Wolf-Rayet (PoWR) non-LTE model atmospheres. The mass-loss rates are adjusted to fit to the observed Ha emission lines best. For the unsaturated UV resonance lines (i.e., P v) we then applied our 3D Monte-Carlo code, which can account for wind clumps of any optical depths ("macroclumping"), a non-void interclump medium, and a velocity dispersion inside the clumps. The ionization stratifications and underlying photospheric spectra were adopted from the PoWR models. The properties of the wind clumps were constrained by fitting the observed resonance line profiles. Results. Our results show that with the mass-loss rates that fit Ha (and other Balmer and He II lines), the UV resonance lines (especially the unsaturated doublet of P v) can also be reproduced with no problem when macroclumping is taken into account. There is no need to artificially reduce the mass-loss rates or to assume a subsolar phosphorus abundance or an extremely high clumping factor, unlike what was claimed by other authors. These consistent mass-loss rates are lower by a factor of 1.3 to 2.6, compared to the mass-loss rate recipe from Vink et al. Conclusions. Macroclumping resolves the previously reported discrepancy between Ha and P v mass-loss diagnostics.
Context. During their evolution massive stars can reach the phase of critical rotation when a further increase in rotational speed is no longer possible. Direct centrifugal ejection from a critically or near-critically rotating surface forms a gaseous equatorial decretion disk. Anomalous viscosity provides the efficient mechanism for transporting the angular momentum outwards. The outer part of the disk can extend up to a very large distance from the parent star.
Aims. We study the evolution of density, radial and azimuthal velocity, and angular momentum loss rate of equatorial decretion disks out to very distant regions. We investigate how the physical characteristics of the disk depend on the distribution of temperature and viscosity.
Methods. We calculated stationary models using the Newton-Raphson method. For time-dependent hydrodynamic modeling we developed the numerical code based on an explicit finite difference scheme on an Eulerian grid including full Navier-Stokes shear viscosity.
Results. The sonic point distance and the maximum angular momentum loss rate strongly depend on the temperature profile and are almost independent of viscosity. The rotational velocity at large radii rapidly drops accordingly to temperature and viscosity distribution. The total amount of disk mass and the disk angular momentum increase with decreasing temperature and viscosity.
Conclusions. The time-dependent one-dimensional models basically confirm the results obtained in the stationary models as well as the assumptions of the analytical approximations. Including full Navier-Stokes viscosity we systematically avoid the rotational velocity sign change at large radii. The unphysical drop of the rotational velocity and angular momentum loss at large radii (present in some models) can be avoided in the models with decreasing temperature and viscosity.
Context. Comprehensive studies of Wolf-Rayet stars were performed in the past for the Galactic and the LMC population. The results revealed significant differences, but also unexpected similarities between the WR populations of these different galaxies. Analyzing the WR stars in M 31 will extend our understanding of these objects in different galactic environments.
Aims. The present study aims at the late-type WN stars in M 31. The stellar and wind parameters will tell about the formation of WR stars in other galaxies with different metallicity and star formation histories. The obtained parameters will provide constraints to the evolution of massive stars in the environment of M 31.
Methods. We used the latest version of the Potsdam Wolf-Rayet model atmosphere code to analyze the stars via fitting optical spectra and photometric data. To account for the relatively low temperatures of the late WN I 0 and WN I I subtypes, our WN models have been extended into this temperature regime.
Results. Stellar and atmospheric parameters are derived for all known late-type WN stars in M 31 with available spectra. All of these stars still have hydrogen in their outer envelopes, some of them up to 50% by mass. The stars are located on the cool side of the zero age main sequence in the Hertzsprung-Russell diagram, while their luminosities range from 105 to 1064). It is remarkable that no star exceeds 106 L.
Conclusions. If formed via single-star evolution, the late-type WN stars in M 31 stem from an initial mass range between 20 and 60 M-circle dot. From the very late-type WN9-11 stars, only one star is located in the S Doradus instability strip. We do not find any late-type WN stars with the high luminosities known in the Milky Way.
The Wolf-Rayet stars in the Large Magellanic Cloud - A comprehensive analysis of the WN class
(2014)
Context. Massive stars, although being important building blocks of galaxies, are still not fully understood. This especially holds true for Wolf-Rayet (WR) stars with their strong mass loss, whose spectral analysis requires adequate model atmospheres. Aims. Following our comprehensive studies of the WR stars in the Milky Way, we now present spectroscopic analyses of almost all known WN stars in the LMC.
Methods. For the quantitative analysis of the wind-dominated emission-line spectra, we employ the Potsdam Wolf-Rayet (PoWR) model atmosphere code. By fitting synthetic spectra to the observed spectral energy distribution and the available spectra (ultraviolet and optical), we obtain the physical properties of 107 stars.
Results. We present the fundamental stellar and wind parameters for an almost complete sample of WN stars in the LMC. Among those stars that are putatively single, two different groups can be clearly distinguished. While 12% of our sample are more luminous than 10(6) L-circle dot and contain a significant amount of hydrogen, 88% of the WN stars, with little or no hydrogen, populate the luminosity range between log (L/L-circle dot) = 5.3 ... 5.8.
Conclusions. While the few extremely luminous stars (log (L/L-circle dot) > 6), if indeed single stars, descended directly from the main sequence at very high initial masses, the bulk of WN stars have gone through the red-supergiant phase. According to their luminosities in the range of log (L/L-circle dot) = 5.3 ... 5.8, these stars originate from initial masses between 20 and 40 M-circle dot. This mass range is similar to the one found in the Galaxy, i.e. the expected metallicity dependence of the evolution is not seen. Current stellar evolution tracks, even when accounting for rotationally induced mixing, still partly fail to reproduce the observed ranges of luminosities and initial masses. Moreover, stellar radii are generally larger and effective temperatures correspondingly lower than predicted from stellar evolution models, probably due to subphotospheric inflation.
We present an overview of four deep phase-constrained Chandra HETGS X-ray observations of delta Ori A. Delta Ori A is actually a triple system that includes the nearest massive eclipsing spectroscopic binary, delta Ori Aa, the only such object that can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, delta Ori Aa2, has a much lower X-ray luminosity than the brighter primary (delta Ori Aa1), delta Ori Aa provides a unique system with which to test the spatial distribution of the X-ray emitting gas around delta Ori Aa1 via occultation by the photosphere of, and wind cavity around, the X-ray dark secondary. Here we discuss the X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time of nearly 500 ks and covering nearly the entire binary orbit. The companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities obtained simultaneously with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectra. We find that the X-ray emission is dominated by embedded wind shock emission from star Aa1, with little contribution from the tertiary star Ab or the shocked gas produced by the collision of the wind of Aa1 against the surface of Aa2. We find a similar temperature distribution to previous X-ray spectrum analyses. We also show that the line half-widths are about 0.3-0.5 times the terminal velocity of the wind of star Aa1. We find a strong anti-correlation between line widths and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the strong and weak lines of Fe XVII and Ne X are inconsistent with model predictions, which may be an effect of resonance scattering.
Context. Wolf-Rayet (WR) stars have a severe impact on their environments owing to their strong ionizing radiation fields and powerful stellar winds. Since these winds are considered to be driven by radiation pressure, it is theoretically expected that the degree of the wind mass-loss depends on the initial metallicity of WR stars.
Aims. Following our comprehensive studies of WR stars in the Milky Way, M31, and the LMC, we derive stellar parameters and mass-loss rates for all seven putatively single WN stars known in the SMC. Based on these data, we discuss the impact of a low-metallicity environment on the mass loss and evolution of WR stars.
Methods. The quantitative analysis of the WN stars is performed with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The physical properties of our program stars are obtained from fitting synthetic spectra to multi-band observations.
Results. In all SMC WN stars, a considerable surface hydrogen abundance is detectable. The majority of these objects have stellar temperatures exceeding 75 kK, while their luminosities range from 10(5.5) to 10(6.1) L-circle dot. The WN stars in the SMC exhibit on average lower mass-loss rates and weaker winds than their counterparts in the Milky Way, M31, and the LMC.
Conclusions. By comparing the mass-loss rates derived for WN stars in different Local Group galaxies, we conclude that a clear dependence of the wind mass-loss on the initial metallicity is evident, supporting the current paradigm that WR winds are driven by radiation. A metallicity effect on the evolution of massive stars is obvious from the HRD positions of the SMC WN stars at high temperatures and high luminosities. Standard evolution tracks are not able to reproduce these parameters and the observed surface hydrogen abundances. Homogeneous evolution might provide a better explanation for their evolutionary past.
Context. Spectroscopic analysis remains the most common method to derive masses of massive stars, the most fundamental stellar parameter. While binary orbits and stellar pulsations can provide much sharper constraints on the stellar mass, these methods are only rarely applicable to massive stars. Unfortunately, spectroscopic masses of massive stars heavily depend on the detailed physics of model atmospheres.
Aims. We demonstrate the impact of a consistent treatment of the radiative pressure on inferred gravities and spectroscopic masses of massive stars. Specifically, we investigate the contribution of line and continuum transitions to the photospheric radiative pressure. We further explore the effect of model parameters, e.g., abundances, on the deduced spectroscopic mass. Lastly, we compare our results with the plane-parallel TLUSTY code, commonly used for the analysis of massive stars with photospheric spectra.
Methods. We calculate a small set of O-star models with the Potsdam Wolf-Rayet (PoWR) code using different approaches for the quasi-hydrostatic part. These models allow us to quantify the effect of accounting for the radiative pressure consistently. We further use PoWR models to show how the Doppler widths of line profiles and abundances of elements such as iron affect the radiative pressure, and, as a consequence, the derived spectroscopic masses.
Results. Our study implies that errors on the order of a factor of two in the inferred spectroscopic mass are to be expected when neglecting the contribution of line and continuum transitions to the radiative acceleration in the photosphere. Usage of implausible microturbulent velocities, or the neglect of important opacity sources such as Fe, may result in errors of approximately 50% in the spectroscopic mass. A comparison with TLUSTY model atmospheres reveals a very good agreement with PoWR at the limit of low mass-loss rates.
We present new grids of Potsdam Wolf-Rayet (PoWR) model atmospheres for Wolf-Rayet stars of the nitrogen sequence (WN stars). The models have been calculated with the latest version of the PoWR stellar atmosphere code for spherical stellar winds. The WN model atmospheres include the non-LTE solutions of the statistical equations for complex model atoms, as well as the radiative transfer equation in the co-moving frame. Iron-line blanketing is treated with the help of the superlevel approach, while wind inhomogeneities are taken into account via optically thin clumps. Three of our model grids are appropriate for Galactic metallicity. The hydrogen mass fraction of these grids is 50%, 20%, and 0%, thus also covering the hydrogen-rich late-type WR stars that have been discovered in recent years. Three grids are adequate for LMC WN stars and have hydrogen fractions of 40%, 20%, and 0%. Recently, additional grids with SMC metallicity and with 60%, 40%, 20%, and 0% hydrogen have been added. We provide contour plots of the equivalent widths of spectral lines that are usually used for classification and diagnostics.
A rather strong mean longitudinal magnetic field of the order of a few hundred gauss was detected a few years ago in the Of?p star CPD -28 degrees 2561 using FORS2 (FOcal Reducer low dispersion Spectrograph 2) low-resolution spectropolarimetric observations. In this work, we present additional low-resolution spectropolarimetric observations obtained during several weeks in 2013 December using FORS 2 mounted at the 8-m Antu telescope of the Very Large Telescope (VLT). These observations cover a little less than half of the stellar rotation period of 73.41 d mentioned in the literature. The behaviour of the mean longitudinal magnetic field is consistent with the assumption of a single-wave variation during the stellar rotation cycle, indicating a dominant dipolar contribution to the magnetic field topology. The estimated polar strength of the surface dipole B-d is larger than 1.15 kG. Further, we compared the behaviour of the line profiles of various elements at different rotation phases associated with different magnetic field strengths. The strongest contribution of the emission component is observed at the phases when the magnetic field shows a negative or positive extremum. The comparison of the spectral behaviour of CPD -28 degrees 2561 with that of another Of?p star, HD 148937 of similar spectral type, reveals remarkable differences in the degree of variability between both stars. Finally, we present new X-ray observations obtained with the Suzaku X-ray Observatory. We report that the star is X-ray bright with log L-X/L-bol approximate to -5.7. The low-resolution X-ray spectra reveal the presence of a plasma heated up to 24 MK. We associate the 24 MK plasma in CPD -28 degrees 2561 with the presence of a kG strong magnetic field capable to confine stellar wind.
We report on both high-precision photometry from the Microvariability and Oscillations of Stars (MOST) space telescope and ground-based spectroscopy of the triple system delta Ori A, consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7 days, and a more distant tertiary (O9 IV P > 400 years). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for three weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy, we have a well-constrained radial velocity (RV) curve of Aa1. While we are unable to recover RV variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries that show evidence for tidally induced pulsations.
We report on interferometric observations with the CHARAArray of two classical Wolf-Rayet (WR) stars in suspected binary systems, namely WR 137 and WR 138. In both cases, we resolve the component stars to be separated by a few milliarcseconds. The data were collected in the H band, and provide a measure of the fractional flux for both stars in each system. We find that the WR star is the dominant H-band light source in both systems (fWR, 137 = 0.59 +/- 0.04; fWR, 138 = 0.67 +/- 0.01), which is confirmed through both comparisons with estimated fundamental parameters for WR stars and O dwarfs, as well as through spectral modelling of each system. Our spectral modelling also provides fundamental parameters for the stars and winds in these systems. The results on WR 138 provide evidence that it is a binary system which may have gone through a previous mass-transfer episode to create the WR star. The separation and position of the stars in the WR 137 system together with previous results from the IOTA interferometer provides evidence that the binary is seen nearly edgeon. The possible edge-on orbit of WR 137 aligns well with the dust production site imaged by the Hubble Space Telescope during a previous periastron passage, showing that the dust production may be concentrated in the orbital plane.
WR 148 (HD 197406) is an extreme runaway system considered to be a potential candidate for a short-period (4.3173 d) rare WR + compact object binary. Provided with new high-resolution, high signal-to-noise spectra from the Keck observatory, we determine the orbital parameters for both the primary WR and the secondary, yielding respective projected orbital velocity amplitudes of 88.1 ± 3.8 km s−1 and 79.2 ± 3.1 km s−1 and implying a mass ratio of 1.1 ± 0.1. We then apply the shift-and-add technique to disentangle the spectra and obtain spectra compatible with a WN7ha and an O4-6 star. Considering an orbital inclination of ∼67°, derived from previous polarimetry observations, the system's total mass would be a mere 2–3M⊙⁠, an unprecedented result for a putative massive binary system. However, a system comprising a 37M⊙ secondary (typical mass of an O5V star) and a 33M⊙ primary (given the mass ratio) would infer an inclination of ∼18°. We therefore reconsider the previous methods of deriving the orbital inclination based on time-dependent polarimetry and photometry. While the polarimetric results are inconclusive requiring better data, the photometric results favour low inclinations. Finally, we compute WR 148’s space velocity and retrace the runaway's trajectory back to the Galactic plane (GP). With an ejection velocity of 198 ± 27 km s−1 and a travel time of 4.7 ± 0.8 Myr to reach its current location, WR 148 was most likely ejected via dynamical interactions in a young cluster.
Context. For more than two decades, stellar atmosphere codes have been used to derive the stellar and wind parameters of massive stars. Although they have become a powerful tool and sufficiently reproduce the observed spectral appearance, they can hardly be used for more than measuring parameters. One major obstacle is their inconsistency between the calculated radiation field and the wind stratification due to the usage of prescribed mass-loss rates and wind-velocity fields. Aims. We present the concepts for a new generation of hydrodynamically consistent non-local thermodynamical equilibrium (nonLTE) stellar atmosphere models that allow for detailed studies of radiation-driven stellar winds. As a first demonstration, this new kind of model is applied to a massive O star. Methods. Based on earlier works, the PoWR code has been extended with the option to consistently solve the hydrodynamic equation together with the statistical equations and the radiative transfer in order to obtain a hydrodynamically consistent atmosphere stratification. In these models, the whole velocity field is iteratively updated together with an adjustment of the mass-loss rate. Results. The concepts for obtaining hydrodynamically consistent models using a comoving-frame radiative transfer are outlined. To provide a useful benchmark, we present a demonstration model, which was motivated to describe the well-studied O4 supergiant zeta Pup. The obtained stellar and wind parameters are within the current range of literature values. Conclusions. For the first time, the PoWR code has been used to obtain a hydrodynamically consistent model for a massive O star. This has been achieved by a profound revision of earlier concepts used for Wolf-Rayet stars. The velocity field is shaped by various elements contributing to the radiative acceleration, especially in the outer wind. The results further indicate that for more dense winds deviations from a standard beta-law occur.