Refine
Has Fulltext
- no (48) (remove)
Year of publication
Language
- English (48)
Is part of the Bibliography
- yes (48)
Keywords
- stars: early-type (48) (remove)
Institute
We perform a statistical clustering analysis of upper main-sequence stars in the Large Magellanic Cloud (LMC) using data from the Visible and Infrared Survey Telescope for Astronomy survey of the Magellanic Clouds. We map over 2500 young stellar structures at 15 significance levels across similar to 120 square degrees centred on the LMC. The structures have sizes ranging from a few parsecs to over 1 kpc. We find that the young structures follow power-law size and mass distributions. From the perimeter-area relation, we derive a perimeter-area dimension of 1.44 +/- 0.20. From the mass-size relation and the size distribution, we derive two-dimensional fractal dimensions of 1.50 +/- 0.10 and 1.61 +/- 0.20, respectively. We find that the surface density distribution is well represented by a lognormal distribution. We apply the Larson relation to estimate the velocity dispersions and crossing times of these structures. Our results indicate that the fractal nature of the young stellar structures has been inherited from the gas clouds from which they form and that this architecture is generated by supersonic turbulence. Our results also suggest that star formation in the LMC is scale-free from 10 to 700 pc.
Context
Line-driven wind instability is expected to cause small-scale wind inhomogeneities, X-ray emission, and wind line profile variability. The instability can already develop around the sonic point if it is initiated close to the photosphere due to stochastic turbulent motions. In such cases, it may leave its imprint on the light curve as a result of wind blanketing.
Aims
We study the photometric signatures of the line-driven wind instability.
Methods
We used line-driven wind instability simulations to determine the wind variability close to the star. We applied two types of boundary perturbations: a sinusoidal one that enables us to study in detail the development of the instability and a stochastic one given by a Langevin process that provides a more realistic boundary perturbation. We estimated the photometric variability from the resulting mass-flux variations. The variability was simulated assuming that the wind consists of a large number of independent conical wind sectors. We compared the simulated light curves with TESS light curves of OB stars that show stochastic variability.
Results
We find two typical signatures of line-driven wind instability in photometric data: a knee in the power spectrum of magnitude fluctuations, which appears due to engulfment of small-scale structure by larger structures, and a negative skewness of the distribution of fluctuations, which is the result of spatial dominance of rarefied regions. These features endure even when combining the light curves from independent wind sectors.
Conclusions
The stochastic photometric variability of OB stars bears certain signatures of the line-driven wind instability. The distribution function of observed photometric data shows negative skewness and the power spectra of a fraction of light curves exhibit a knee. This can be explained as a result of the line-driven wind instability triggered by stochastic base perturbations.
Classical Wolf-Rayet (cWR) stars are at a crucial evolutionary stage for constraining the fates of massive stars. The feedback of these hot, hydrogen-depleted stars dominates their surrounding by tremendous injections of ionizing radiation and kinetic energy. The strength of a Wolf-Rayet (WR) wind decides the eventual mass of its remnant, likely a massive black hole. However, despite their major influence and importance for gravitational wave detection statistics, WR winds are particularly poorly understood. In this paper, we introduce the first set of hydrodynamically consistent stellar atmosphere models for cWR stars of both the carbon (C) and the nitrogen (N) sequence, i.e. WC and WN stars, as a function of stellar luminosity-to-mass ratio (or Eddington Gamma) and metallicity. We demonstrate the inapplicability of the CAK wind theory for cWR stars and confirm earlier findings that their winds are launched at the (hot) iron (Fe) opacity peak. For log Z/Z(circle dot) > -2, Fe is also the main accelerator throughout the wind. Contrasting previous claims of a sharp lower mass-loss limit forWR stars, we obtain a smooth transition to optically thin winds. Furthermore, we find a strong dependence of the mass-loss rates on Eddington Gamma, both at solar and subsolar metallicity. Increases inWCcarbon and oxygen abundances turn out to slightly reduce the predicted mass-loss rates. Calculations at subsolar metallicities indicate that below the metallicity of the Small Magellanic Cloud, WR mass-loss rates decrease much faster than previously assumed, potentially allowing for high black hole masses even in the local Universe.
Context
The recently claimed discovery of a massive (M-BH = 68(-13)(+11) M-circle dot) black hole in the Galactic solar neighborhood has led to controversial discussions because it severely challenges our current view of stellar evolution.
Aims
A crucial aspect for the determination of the mass of the unseen black hole is the precise nature of its visible companion, the B-type star LSV +22 25. Because stars of different mass can exhibit B-type spectra during the course of their evolution, it is essential to obtain a comprehensive picture of the star to unravel its nature and, thus, its mass.
Methods
To this end, we study the spectral energy distribution of LSV +22 25 and perform a quantitative spectroscopic analysis that includes the determination of chemical abundances for He, C, N, O, Ne, Mg, Al, Si, S, Ar, and Fe.
Results
Our analysis clearly shows that LSV +22 25 is not an ordinary main sequence B-type star. The derived abundance pattern exhibits heavy imprints of the CNO bi-cycle of hydrogen burning, that is, He and N are strongly enriched at the expense of C and O. Moreover, the elements Mg, Al, Si, S, Ar, and Fe are systematically underabundant when compared to normal main-sequence B-type stars. We suggest that LSV +22 25 is a stripped helium star and discuss two possible formation scenarios. Combining our photometric and spectroscopic results with the Gaia parallax, we infer a stellar mass of 1.1 +/- 0.5 M-circle dot. Based on the binary system's mass function, this yields a minimum mass of 2-3 M-circle dot for the compact companion, which implies that it may not necessarily be a black hole but a massive neutron- or main sequence star.
Conclusions
The star LSV +22 25 has become famous for possibly having a very massive black hole companion. However, a closer look reveals that the star itself is a very intriguing object. Further investigations are necessary for complete characterization of this object.
Context:
Massive binaries play a crucial role in the Universe. Knowing the distributions of their orbital parameters is important for a wide range of topics from stellar feedback to binary evolution channels and from the distribution of supernova types to gravitational wave progenitors, yet no direct measurements exist outside the Milky Way.
Aims:
The Tarantula Massive Binary Monitoring project was designed to help fill this gap by obtaining multi-epoch radial velocity (RV) monitoring of 102 massive binaries in the 30 Doradus region.
Methods:
In this paper we analyze 32 FLAMES/GIRAFFE observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single-lined (SB1) and 31 double-lined (SB2) spectroscopic binaries.
Results:
Overall, the binary fraction and orbital properties across the 30 Doradus region are found to be similar to existing Galactic samples. This indicates that within these domains environmental effects are of second order in shaping the properties of massive binary systems. A small difference is found in the distribution of orbital periods, which is slightly flatter (in log space) in 30 Doradus than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, orbital periods in 30 Doradus can be as short as 1.1 d, somewhat shorter than seen in Galactic samples. Equal mass binaries (q> 0.95) in 30 Doradus are all found outside NGC 2070, the central association that surrounds R136a, the very young and massive cluster at 30 Doradus’s core. Most of the differences, albeit small, are compatible with expectations from binary evolution. One outstanding exception, however, is the fact that earlier spectral types (O2–O7) tend to have shorter orbital periods than later spectral types (O9.2–O9.7).
Conclusions:
Our results point to a relative universality of the incidence rate of massive binaries and their orbital properties in the metallicity range from solar (Z⊙) to about half solar. This provides the first direct constraints on massive binary properties in massive star-forming galaxies at the Universe’s peak of star formation at redshifts z ~ 1 to 2 which are estimated to have Z ~ 0.5 Z⊙.
B fields in OB stars (BOB)
(2017)
Aims. The B fields in OB stars (BOB) Collaboration is based on an ESO Large Programme to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. Methods. In the framework of this program, we carried out low-resolution spectropolarimetric observations of a large sample of massive stars using FORS2 installed at the ESO VLT 8m telescope. Results. We determined the magnetic field values with two completely independent reduction and analysis pipelines. Our in-depth study of the magnetic field measurements shows that differences between our two pipelines are usually well within 3 sigma errors. From the 32 observations of 28 OB stars, we were able to monitor the magnetic fields in CPD -57 degrees 3509 and HD164492C, confirm the magnetic field in HD54879, and detect a magnetic field in CPD -62 degrees 2124. We obtain a magnetic field detection rate of 6 +/- 3% for the full sample of 69 OB stars observed with FORS 2 within the BOB program. For the preselected objects with a nu sin i below 60 km s(-1), we obtain a magnetic field detection rate of 5 +/- 5%. We also discuss X-ray properties and multiplicity of the objects in our FORS2 sample with respect to the magnetic field detections.
Context. Massive stars severely influence their environment by their strong ionizing radiation and by the momentum and kinetic energy input provided by their stellar winds and supernovae. Quantitative analyses of massive stars are required to understand how their feedback creates and shapes large scale structures of the interstellar medium. The giant H II region N206 in the Large Magellanic Cloud contains an OB association that powers a superbubble filled with hot X-ray emitting gas, serving as an ideal laboratory in this context. Aims. We aim to estimate stellar and wind parameters of all OB stars in N206 by means of quantitative spectroscopic analyses. In this first paper, we focus on the nine Of-type stars located in this region. We determine their ionizing flux and wind mechanical energy. The analysis of nitrogen abundances in our sample probes rotational mixing. Methods. We obtained optical spectra with the multi-object spectrograph FLAMES at the ESO-VLT. When possible, the optical spectroscopy was complemented by UV spectra from the HST, IUE, and FUSE archives. Detailed spectral classifications are presented for our sample Of-type stars. For the quantitative spectroscopic analysis we used the Potsdam Wolf-Rayet model atmosphere code. We determined the physical parameters and nitrogen abundances of our sample stars by fitting synthetic spectra to the observations. Results. The stellar and wind parameters of nine Of-type stars, which are largely derived from spectral analysis are used to construct wind momentum luminosity relationship. We find that our sample follows a relation close to the theoretical prediction, assuming clumped winds. The most massive star in the N206 association is an Of supergiant that has a very high mass-loss rate. Two objects in our sample reveal composite spectra, showing that the Of primaries have companions of late O subtype. All stars in our sample have an evolutionary age of less than 4 million yr, with the O2-type star being the youngest. All these stars show a systematic discrepancy between evolutionary and spectroscopic masses. All stars in our sample are nitrogen enriched. Nitrogen enrichment shows a clear correlation with increasing projected rotational velocities. Conclusions. The mechanical energy input from the Of stars alone is comparable to the energy stored in the N206 superbubble as measured from the observed X-ray and H alpha emission.
Context:
The first directly detected gravitational waves (GW 150914) were emitted by two coalescing black holes (BHs) with masses of ≈ 36 M⊙ and ≈ 29 M⊙. Several scenarios have been proposed to put this detection into an astrophysical context. The evolution of an isolated massive binary system is among commonly considered models.
Aims:
Various groups have performed detailed binary-evolution calculations that lead to BH merger events. However, the question remains open as to whether binary systems with the predicted properties really exist. The aim of this paper is to help observers to close this gap by providing spectral characteristics of massive binary BH progenitors during a phase where at least one of the companions is still non-degenerate.
Methods:
Stellar evolution models predict fundamental stellar parameters. Using these as input for our stellar atmosphere code (Potsdam Wolf-Rayet), we compute a set of models for selected evolutionary stages of massive merging BH progenitors at different metallicities.
Results:
The synthetic spectra obtained from our atmosphere calculations reveal that progenitors of massive BH merger events start their lives as O2-3V stars that evolve to early-type blue supergiants before they undergo core-collapse during the Wolf-Rayet phase. When the primary has collapsed, the remaining system will appear as a wind-fed high-mass X-ray binary. Based on our atmosphere models, we provide feedback parameters, broad band magnitudes, and spectral templates that should help to identify such binaries in the future.
Conclusions:
While the predicted parameter space for massive BH binary progenitors is partly realized in nature, none of the known massive binaries match our synthetic spectra of massive BH binary progenitors exactly. Comparisons of empirically determined mass-loss rates with those assumed by evolution calculations reveal significant differences. The consideration of the empirical mass-loss rates in evolution calculations will possibly entail a shift of the maximum in the predicted binary-BH merger rate to higher metallicities, that is, more candidates should be expected in our cosmic neighborhood than previously assumed.
Circumstantial evidence suggests that magnetism and enhanced X-ray emission are likely correlated in early B-type stars: similar fractions of them (similar to 10%) are strong and hard X-ray sources and possess strong magnetic fields. It is also known that some B-type stars have spots on their surface. Yet up to now no X-ray activity associated with spots on early-type stars was detected. In this Letter we report the detection of a magnetic field on the B2V star rho Oph A. Previously, we assessed that the X-ray activity of this star is associated with a surface spot, herewith we establish its magnetic origin. We analyze spectra of rho Oph A obtained with the FORS2 spectrograph at ESO Very Large Telescope (VLT) at two epochs, and detect a longitudinal component of the magnetic field of the order of similar to 500 G in one of the datasets. The detection of the magnetic field only at one epoch can be explained by stellar rotation which is also invoked to explain observed periodic X-ray activity. From archival HARPS ESO VLT high resolution spectra we derived the fundamental stellar parameters of rho Oph A and further constrained its age. We conclude that rho Oph A provides strong evidence for the presence of active X-ray emitting regions on young magnetized early type stars.
Aims. To gain a realistic picture of the donor star in Vela X-1, we constructed a hydrodynamically consistent atmosphere model describing the wind stratification while properly reproducing the observed donor spectrum. To investigate how X-ray illumination affects the stellar wind, we calculated additional models for different X-ray luminosity regimes. Methods. We used the recently updated version of the Potsdam Wolf-Rayet code to consistently solve the hydrodynamic equation together with the statistical equations and the radiative transfer. Results. The wind flow in Vela X-1 is driven by ions from various elements, with Fe III and S III leading in the outer wind. The model-predicted mass-loss rate is in line with earlier empirical studies. The mass-loss rate is almost unaffected by the presence of the accreting NS in the wind. The terminal wind velocity is confirmed at u(infinity) approximate to 600 km s(-1). On the other hand, the wind velocity in the inner region where the NS is located is only approximate to 100 km s(-1), which is not expected on the basis of a standard beta-velocity law. In models with an enhanced level of X-rays, the velocity field in the outer wind can be altered. If the X-ray flux is too high, the acceleration breaks down because the ionization increases. Conclusions. Accounting for radiation hydrodynamics, our Vela X-1 donor atmosphere model reveals a low wind speed at the NS location, and it provides quantitative information on wind driving in this important HMXB.
Context. Runaway stars form bow shocks by ploughing through the interstellar medium at supersonic speeds and are promising sources of non-thermal emission of photons. One of these objects has been found to emit non-thermal radiation in the radio band. This triggered the development of theoretical models predicting non-thermal photons from radio up to very-high-energy (VHE, E >= 0.1 TeV) gamma rays. Subsequently, one bow shock was also detected in X-ray observations. However, the data did not allow discrimination between a hot thermal and a non-thermal origin. Further observations of different candidates at X-ray energies showed no evidence for emission at the position of the bow shocks either. A systematic search in the Fermi-LAT energy regime resulted in flux upper limits for 27 candidates listed in the E-BOSS catalogue. Aims. Here we perform the first systematic search for VHE gamma-ray emission from bow shocks of runaway stars. Methods. Using all available archival H.E.S.S. data we search for very-high-energy gamma-ray emission at the positions of bow shock candidates listed in the second E-BOSS catalogue release. Out of the 73 bow shock candidates in this catalogue, 32 have been observed with H.E.S.S. Results. None of the observed 32 bow shock candidates in this population study show significant emission in the H.E.S.S. energy range. Therefore, flux upper limits are calculated in five energy bins and the fraction of the kinetic wind power that is converted into VHE gamma rays is constrained. Conclusions. Emission from stellar bow shocks is not detected in the energy range between 0.14 and 18 TeV. The resulting upper limits constrain the level of VHE gamma-ray emission from these objects down to 0.1-1% of the kinetic wind energy.
We present new radio/millimeter measurements of the hot magnetic star HR5907 obtained with the VLA and ALMA interferometers. We find that HR5907 is the most radio luminous early type star in the cm-mm band among those presently known. Its multi-wavelength radio light curves are strongly variable with an amplitude that increases with radio frequency. The radio emission can be explained by the populations of the non-thermal electrons accelerated in the current sheets on the outer border of the magnetosphere of this fast-rotating magnetic star. We classify HR5907 as another member of the growing class of strongly magnetic fast-rotating hot stars where the gyro-synchrotron emission mechanism efficiently operates in their magnetospheres. The new radio observations of HR5907 are combined with archival X-ray data to study the physical condition of its magnetosphere. The X-ray spectra of HR5907 show tentative evidence for the presence of non-thermal spectral component. We suggest that non-thermal X-rays originate a stellar X-ray aurora due to streams of non-thermal electrons impacting on the stellar surface. Taking advantage of the relation between the spectral indices of the X-ray power-law spectrum and the non-thermal electron energy distributions, we perform 3-D modelling of the radio emission for HR5907. The wavelength-dependent radio light curves probe magnetospheric layers at different heights above the stellar surface. A detailed comparison between simulated and observed radio light curves leads us to conclude that the stellar magnetic field of HR 5907 is likely non-dipolar, providing further indirect evidence of the complex magnetic field topology of HR5907.
A small fraction of the radiative flux emitted by hot stars is absorbed by their winds and redistributed towards longer wavelengths. This effect, which leads also to the heating of the stellar photosphere, is termed wind blanketing. For stars with variable winds, the effect of wind blanketing may lead to the photometric variability. We have studied the consequences of line driven wind instability and wind blanketing for the light variability of O stars. We combined the results of wind hydrodynamic simulations and of global wind models to predict the light variability of hot stars due to the wind blanketing and instability. The wind instability causes stochastic light variability with amplitude of the order of tens of millimagnitudes and a typical timescale of the order of hours for spatially coherent wind structure. The amplitude is of the order of millimagnitudes when assuming that the wind consists of large number of independent concentric cones. The variability with such amplitude is observable using present space borne photometers. We show that the simulated light curve is similar to the light curves of O stars obtained using BRITE and CoRoT satellites.
We report the detection of the auroral radio emission from the early-type magnetic star HD142301. New VLA observations of HD142301 detected highly polarized amplified emission occurring at fixed stellar orientations. The coherent emission mechanism responsible for the stellar auroral radio emission amplifies the radiation within a narrow beam, making the star where this phenomenon occurs similar to a radio lighthouse. The elementary emission process responsible for the auroral radiation mainly amplifies one of the two magneto-ionic modes of the electromagnetic wave. This explains why the auroral pulses are highly circularly polarized. The auroral radio emission of HD142301 is characterized by a reversal of the sense of polarization as the star rotates. The effective magnetic field curve of HD142301 is also available making it possible to correlate the transition from the left to the right-hand circular polarization sense ( and vice versa) of the auroral pulses with the known orientation of the stellar magnetic field. The results presented in this letter have implications for the estimation of the dominant magneto-ionic mode amplified within the HD142301 magnetosphere.
The study of massive stars in different metallicity environments is a central topic of current stellar research. The spectral analysis of massive stars requires adequate model atmospheres. The computation of such models is difficult and time-consuming. Therefore, spectral analyses are greatly facilitated if they can refer to existing grids of models. Here we provide grids of model atmospheres for OB-type stars at metallicities corresponding to the Small and Large Magellanic Clouds, as well as to solar metallicity. In total, the grids comprise 785 individual models. The models were calculated using the state-of-the-art Potsdam Wolf-Rayet (PoWR) model atmosphere code. The parameter domain of the grids was set up using stellar evolution tracks. For all these models, we provide normalized and flux-calibrated spectra, spectral energy distributions, feedback parameters such as ionizing photons, Zanstra temperatures, and photometric magnitudes. The atmospheric structures (the density and temperature stratification) are available as well. All these data are publicly accessible through the PoWR website.
We present new, contemporaneous Hubble Space Telescope STIS and XMM-Newton observations of the O7. III(n) ((f)) star xi Per. We supplement the new data with archival IUE spectra, to analyze the variability of the wind lines and X-ray flux of xi Per. The variable wind of this star is known to have a 2.086-day periodicity. We use a simple, heuristic spot model that fits the low-velocity (near-surface) IUE wind line variability very well, to demonstrate that the low-velocity absorption in the new STIS spectra of N IV lambda 1718 and Si IV lambda 1402 vary with the same 2.086-day period. It is remarkable that the period and amplitude of the STIS data agree with those of the IUE spectra obtained 22 yr earlier. We also show that the time variability of the new XMM-Newton fluxes is also consistent with the 2.086-day period. Thus, our new, multiwavelength coordinated observations demonstrate that the mechanism that causes the UV wind line variability is also responsible for a significant fraction of the X-rays in single O stars. The sequence of events for the multiwavelength light-curve minima is Si IV lambda 1402, N IV lambda 1718, and X-ray flux, each separated by a phase of about 0.06 relative to the 2.086-day period. Analysis of the X-ray fluxes shows that they become softer as they weaken. This is contrary to expectations if the variability is caused by periodic excess absorption. Furthermore, the high-resolution X-ray spectra suggest that the individual emission lines at maximum are more strongly blueshifted. If we interpret the low-velocity wind line light curves in terms of our model, it implies that there are two bright regions, i.e., regions with less absorption, separated by 180 degrees, on the surface of the star. We note that the presence and persistence of two spots separated by 180 degrees suggest that a weak dipole magnetic field is responsible for the variability of the UV wind line absorption and X-ray flux in xi Per.
Hypervelocity stars are rare objects, mostly main-sequence (MS) B stars, traveling so fast that they will eventually escape from the Milky Way. Recently, it has been shown that the popular Hills mechanism, in which a binary system is disrupted via a close encounter with the supermassive black hole at the Galactic center, may not be their only ejection mechanism. The analyses of Gaia data ruled out a Galactic center origin for some of them, and instead indicated that they are extreme disk runaway stars ejected at velocities exceeding the predicted limits of classical scenarios (dynamical ejection from star clusters or binary supernova ejection). We present the discovery of a new extreme disk runaway star, PG 1610+062, which is a slowly pulsating B star bright enough to be studied in detail. A quantitative analysis of spectra taken with ESI at the Keck Observatory revealed that PG 1610+062 is a late B-type MS star of 4–5 M⊙ with low projected rotational velocity. Abundances (C, N, O, Ne, Mg, Al, Si, S, Ar, and Fe) were derived differentially with respect to the normal B star HD 137366 and indicate that PG 1610+062 is somewhat metal rich. A kinematic analysis, based on our spectrophotometric distance (17.3 kpc) and on proper motions from Gaia’s second data release, shows that PG 1610+062 was probably ejected from the Carina-Sagittarius spiral arm at a velocity of 550 ± 40 km s−1, which is beyond the classical limits. Accordingly, the star is in the top five of the most extreme MS disk runaway stars and is only the second among the five for which the chemical composition is known.
In this paper, we investigate the multiwavelength properties of the magnetic early B-type star HR 7355. We present its radio light curves at several frequencies, taken with the Jansky Very Large Array, and X-ray spectra, taken with the XMM-Newton X-ray telescope. Modelling of the radio light curves for the Stokes I and V provides a quantitative analysis of the HR 7355 magnetosphere. A comparison between HR 7355 and a similar analysis for the Ap star CU Vir allows us to study how the different physical parameters of the two stars affect the structure of the respective magnetospheres where the non-thermal electrons originate. Our model includes a cold thermal plasma component that accumulates at high magnetic latitudes that influences the radio regime, but does not give rise to X-ray emission. Instead, the thermal X-ray emission arises from shocks generated by wind stream collisions close to the magnetic equatorial plane. The analysis of the X-ray spectrum of HR 7355 also suggests the presence of a non-thermal radiation. Comparison between the spectral index of the power-law X-ray energy distribution with the non-thermal electron energy distribution indicates that the non-thermal X-ray component could be the auroral signature of the non-thermal electrons that impact the stellar surface, the same non-thermal electrons that are responsible for the observed radio emission. On the basis of our analysis, we suggest a novel model that simultaneously explains the X-ray and the radio features of HR 7355 and is likely relevant for magnetospheres of other magnetic early-type stars.
We present the results of a 140 ks XMM-Newton observation of the B2 star rho Oph A. The star has exhibited strong X-ray variability: a cusp-shaped increase of rate, similar to that which we partially observed in 2013, and a bright flare. These events are separated in time by about 104 ks, which likely correspond to the rotational period of the star (1.2 days). Time resolved spectroscopy of the X-ray spectra shows that the first event is caused by an increase of the plasma emission measure, while the second increase of rate is a major flare with temperatures in excess of 60 MK (kT similar to 5 keV). From the analysis of its rise, we infer a magnetic field of >= 300 G and a size of the flaring region of similar to 1.4-1.9 x 10(11) cm, which corresponds to similar to 25%-30% of the stellar radius. We speculate that either an intrinsic magnetism that produces a hot spot on its surface or an unknown low mass companion are the source of such X-rays and variability. A hot spot of magnetic origin should be a stable structure over a time span of >= 2.5 yr, and suggests an overall large scale dipolar magnetic field that produces an extended feature on the stellar surface. In the second scenario, a low mass unknown companion is the emitter of X-rays and it should orbit extremely close to the surface of the primary in a locked spin-orbit configuration, almost on the verge of collapsing onto the primary. As such, the X-ray activity of the secondary star would be enhanced by its young age, and the tight orbit as in RS Cvn systems. In both cases rho Oph would constitute an extreme system that is worthy of further investigation.
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.