non-LTE Model Atmospheres Research Articles

Spectral analysis of the He-enriched sdO-star HD 127493

Abstract The bright sdO star HD127493 is known to be of mixed H/He composition and excellent archival spectra covering both optical and ultraviolet ranges are available. UV spectra play a key role as they give access to many chemical species that do not show spectral lines in the optical, such as iron and nickel. This encouraged the quantitative spectral analysis of this prototypical mixed H/He composition sdO star. We determined atmospheric parameters for HD127493 in addition to the abundance of C, N, O, Si, S, Fe, and Ni in the atmosphere using non-LTE model atmospheres calculated with TLUSTY/SYNSPEC. A comparison between the parallax distance measured by Hipparcos and the derived spectroscopic distance indicate that the derived atmospheric parameters are realistic. From our metal abundance analysis, we find a strong CNO signature and enrichment in iron and nickel.

Chemical fingerprints of He-sdO stars

Abstract The chemical composition of helium-rich hot subluminous O stars plays an important role to understand and model their formation history. We present a spectroscopic analysis of four He-sdO stars,CD-31° 4800, [CW83] 0904- 02, LSS 1274 and LS IV +10° 9. The analysis is based on archival optical and UV high-resolution spectra. We used Tlusty200/Synspec48 to compute line blanketed non-LTE model atmospheres and their corresponding synthetic spectra and derive the atmospheric parameters as well as the abundances of the most prominent elements. All stars have helium-dominated atmospheres with hardly any hydrogen and temperatures between 42000 K and 47000 K while their surface gravity spans between log g = 5.4 and 5.7. CD-31° 4800 shows an enrichment of nitrogen and the characteristic pattern of hydrogen burning via the CNO-cycle, while the rest of the elements have about the solar abundance. This points to the slow merger of two helium white dwarfs as the most likely origin for this system. The other three stars are enriched in carbon, nitrogen and neon while their intermediate mass element’s abundance scatters around the solar value. They were possibly formed in the deep mixing late hot flasher scenario.

Bi-Abundance Ionisation Structure of the Wolf–Rayet Planetary Nebula PB 8

AbstractThe planetary nebula PB 8 around a [WN/WC]-hybrid central star is one of planetary nebulae with moderate abundance discrepancy factors (ADFs ~ 2–3), which could be an indication of a tiny fraction of metal-rich inclusions embedded in the nebula (bi-abundance). In this work, we have constructed photoionisation models to reproduce the optical and infrared observations of the planetary nebula PB 8 using a non-LTE stellar model atmosphere ionising source. A chemically homogeneous model initially used cannot predict the optical recombination lines. However, a bi-abundance model provides a better fit to most of the observed optical recombination lines from N and O ions. The metal-rich inclusions in the bi-abundance model occupy 5.6% of the total volume of the nebula, and are roughly 1.7 times cooler and denser than the mean values of the surrounding nebula. The N/H and O/H abundance ratios in the metal-rich inclusions are ~1.0 and 1.7 dex larger than the diffuse warm nebula, respectively. To reproduce the Spitzer spectral energy distribution of PB 8, dust grains with a dust-to-gas ratio of 0.01 (by mass) were also included. It is found that the presence of metal-rich inclusions can explain the heavy element optical recombination lines, while a dual-dust chemistry with different grain species and discrete grain sizes likely produces the infrared continuum of this planetary nebula. This study demonstrates that the bi-abundance hypothesis, which was examined in a few planetary nebulae with large abundance discrepancies (ADFs > 10), could also be applied to those typical planetary nebulae with moderate abundance discrepancies.

Candidate exoplanet host HD 131399A: a nascent Am star

Direct imaging suggests that there is a Jovian exoplanet around the primary A-star in the triple-star system HD131399. We investigate a high-quality spectrum of the primary component HD131399A obtained with FEROS on the ESO/MPG 2.2m telescope, aiming to characterise the star's atmospheric and fundamental parameters, and to determine elemental abundances at high precision and accuracy. The aim is to constrain the chemical composition of the birth cloud of the system and therefore the bulk composition of the putative planet. A hybrid non-local thermal equilibrium (non-LTE) model atmosphere technique is adopted for the quantitative spectral analysis. Comparison with the most recent stellar evolution models yields the fundamental parameters. The atmospheric and fundamental stellar parameters of HD131399A are constrained to Teff=9200+-100 K, log g=4.37+-0.10, M=1.95+0.08-0.06 Msun, R=1.51+0.13-0.10 Rsun, and log L/Lsun=1.17+-0.07, locating the star on the zero-age main sequence. Non-LTE effects on the derived metal abundances are often smaller than 0.1dex, but can reach up to ~0.8dex for individual lines. The observed lighter elements up to calcium are overall consistent with present-day cosmic abundances, with a C/O ratio of 0.45$\pm$0.07 by number, while the heavier elements show mild overabundances. We conclude that the birth cloud of the system had a standard chemical composition, but we witness the onset of the Am phenomenon in the slowly rotating star. We furthermore show that non-LTE analyses have the potential to solve the remaining discrepancies between observed abundances and predictions by diffusion models for Am stars. Moreover, the present case allows mass loss, not turbulent mixing, to be identified as the main transport process competing with diffusion in very young Am stars.

The Origin of B-type Runaway Stars: Non-LTE Abundances as a Diagnostic

Abstract There are two accepted mechanisms to explain the origin of runaway OB-type stars: the binary supernova (SN) scenario and the cluster ejection scenario. In the former, an SN explosion within a close binary ejects the secondary star, while in the latter close multibody interactions in a dense cluster cause one or more of the stars to be ejected from the region at high velocity. Both mechanisms have the potential to affect the surface composition of the runaway star. tlusty non-LTE model atmosphere calculations have been used to determine the atmospheric parameters and the C, N, Mg, and Si abundances for a sample of B-type runaways. These same analytical tools were used by Hunter et al. for their analysis of 50 B-type open-cluster Galactic stars (i.e., nonrunaways). Effective temperatures were deduced using the Si-ionization balance technique, surface gravities from Balmer line profiles, and microturbulent velocities derived using the Si spectrum. The runaways show no obvious abundance anomalies when compared with stars in the open clusters. The runaways do show a spread in composition that almost certainly reflects the Galactic abundance gradient and a range in the birthplaces of the runaways in the Galactic disk. Since the observed Galactic abundance gradients of C, N, Mg, and Si are of a similar magnitude, the abundance ratios (e.g., N/Mg) are as obtained essentially uniform across the sample.

HOT SUBDWARF STARS OBSERVED IN LAMOST DR1—ATMOSPHERIC PARAMETERS FROM SINGLE-LINED SPECTRA

ABSTRACT We present a catalog of 166 spectroscopically identified hot subdwarf stars from LAMOST DR1, 44 of which show the characteristics of cool companions in their optical spectra. Atmospheric parameters of 122 subdwarf stars with non-composite spectra were measured by fitting the profiles of hydrogen (H) and helium (He) lines with synthetic spectra from non-LTE model atmospheres. Most of the sdB stars scatter near the Extreme Horizontal Branch in the T eff– log g ?> diagram and two well defined groups can be outlined. A clustering of He-enriched sdO stars appears near T eff = 45,000 K and log g = 5.8. The sdB population separates into several nearly parallel sequences in the T eff–He abundance diagram with clumps corresponding to those in the T eff– log g ?> diagram. Over 38,000 K (sdO) stars show abundance extremes; they are either He-rich or He-deficient and we observe only a few stars in the abundance range − 1 < log y < 0 ?> . With increasing temperature these extremes become less prominent and the He abundance approaches log y ∼ −0.5. A unique property of our sample is that it covers a large range in apparent magnitude and galactic latitude, therefore it contains a mix of stars from different populations and galactic environments. Our results are consistent with the findings of Hirsch and we conclude that He-rich and He-deficient sdB stars (log y < 1) probably originate from different populations. We also find that most sdO and sdB stars lie in a narrow strip in the plane of luminosity and helium abundance, which suggests that these atmospheric parameters are correlated.

Two bi-stability jumps in theoretical wind models for massive stars and the implications for luminous blue variable supernovae

Luminous Blue Variables have been suggested to be the direct progenitors of supernova types IIb and IIn, with enhanced mass loss prior to explosion. However, the mechanism of this mass loss is not yet known. Here, we investigate the qualitative behaviour of theoretical stellar wind mass-loss as a function of Teff across two bi-stability jumps in blue supergiant regime and also in proximity to the Eddington limit, relevant for LBVs. To investigate the physical ingredients that play a role in the radiative acceleration we calculate blue supergiant wind models with the CMFGEN non-LTE model atmosphere code over an effective temperature range between 30000 and 8800 K. Although our aim is not to provide new mass-loss rates for BA supergiants, we study and confirm the existence of two bi-stability jumps in mass-loss rates predicted by Vink, de Koter, & Lamers (1999). However, they are found to occur at somewhat lower Teff (20000 and 9000 K, respectively) than found previously, which would imply that stars may evolve towards lower Teff before strong mass-loss is induced by the bi-stability jumps. When the combined effects of the second bi-stability jump and the proximity to Eddington limit are accounted for, we find a dramatic increase in the mass-loss rate by up to a factor of 30. Further investigation of both bi-stability jumps is expected to lead to a better understanding of discrepancies between empirical modelling and theoretical mass-loss rates reported in the literature, and to provide key inputs for the evolution of both normal AB supergiants and LBVs, as well as their subsequent supernova type II explosions.

B fields in OB stars (BOB): Detection of a magnetic field in the He-strong star CPD −57° 3509

We report the detection of a magnetic field in the helium-strong star CPD-57 3509 (B2 IV), a member of the Galactic open cluster NGC3293, and characterise the star's atmospheric and fundamental parameters. Spectropolarimetric observations with FORS2 and HARPSpol are analysed using two independent approaches to quantify the magnetic field strength. A high-S/N FLAMES/GIRAFFE spectrum is analysed using a hybrid non-LTE model atmosphere technique. Comparison with stellar evolution models constrains the fundamental parameters of the star. We obtain a firm detection of a surface averaged longitudinal magnetic field with a maximum amplitude of about 1 kG. Assuming a dipolar configuration of the magnetic field, this implies a dipolar field strength larger than 3.3 kG. Moreover, the large amplitude and fast variation (within about 1 day) of the longitudinal magnetic field implies that CPD-57 3509 is spinning very fast despite its apparently slow projected rotational velocity. The star should be able to support a centrifugal magnetosphere, yet the spectrum shows no sign of magnetically confined material; in particular, emission in H{\alpha} is not observed. Apparently, the wind is either not strong enough for enough material to accumulate in the magnetosphere to become observable or, alternatively, some leakage process leads to loss of material from the magnetosphere. The quantitative spectroscopic analysis of the star yields an effective temperature and a logarithmic surface gravity of 23750+-250 K and 4.05+-0.10, respectively, and a surface helium fraction of 0.28+-0.02 by number. The surface abundances of C, N, O, Ne, S, and Ar are compatible with the cosmic abundance standard, whereas Mg, Al, Si, and Fe are depleted by about a factor of 2. This abundance pattern can be understood as the consequence of a fractionated stellar wind. CPD-57 3509 is one of the most evolved He-strong stars known.

Quantitative spectral analysis of the sdB star HD 188112: A helium-core white dwarf progenitor

HD 188112 is a bright (V = 10.2 mag) hot subdwarf B (sdB) star with a mass too low to ignite core helium burning and is therefore considered as a pre-extremely low mass (ELM) white dwarf (WD). ELM WDs (M $\le$ 0.3 Msun) are He-core objects produced by the evolution of compact binary systems. We present in this paper a detailed abundance analysis of HD 188112 based on high-resolution Hubble Space Telescope (HST) near and far-ultraviolet spectroscopy. We also constrain the mass of the star's companion. We use hybrid non-LTE model atmospheres to fit the observed spectral lines and derive the abundances of more than a dozen elements as well as the rotational broadening of metallic lines. We confirm the previous binary system parameters by combining radial velocities measured in our UV spectra with the already published ones. The system has a period of 0.60658584 days and a WD companion with M $\geq$ 0.70 Msun. By assuming a tidally locked rotation, combined with the projected rotational velocity (v sin i = 7.9 $\pm$ 0.3 km s$^{-1}$) we constrain the companion mass to be between 0.9 and 1.3 Msun. We further discuss the future evolution of the system as a potential progenitor of a (underluminous) type Ia supernova. We measure abundances for Mg, Al, Si, P, S, Ca, Ti, Cr, Mn, Fe, Ni, and Zn, as well as for the trans-iron elements Ga, Sn, and Pb. In addition, we derive upper limits for the C, N, O elements and find HD 188112 to be strongly depleted in carbon. We find evidence of non-LTE effects on the line strength of some ionic species such as Si II and Ni II. The metallic abundances indicate that the star is metal-poor, with an abundance pattern most likely produced by diffusion effects.

Limb-darkening coefficients from line-blanketed non-LTE hot-star model atmospheres

We present grids of limb-darkening coefficients computed from non-LTE, line-blanketed TLUSTY model atmospheres, covering effective-temperature and surface-gravity ranges of 15--55kK and 4.75 dex (cgs) down to the effective Eddington limit, at 1x, 1x, 0.5x (LMC), 0.2x (SMC), and 0.1x solar. Results are given for the Bessell UBVRIJKHL, Sloan ugriz, Stromgren ubvy, WFCAM ZYJHK, Hipparcos, Kepler, and Tycho passbands, in each case characterized by several different limb-darkening `laws'. We examine the sensitivity of limb darkening to temperature, gravity, metallicity, microturbulent velocity, and wavelength, and make a comparison with LTE models. The dependence on metallicity is very weak, but limb darkening is a moderately strong function of log(g) in this temperature regime.

Light-travel-time diagnostics in early supernova spectra: substantial mass-loss of the IIb progenitor of SN 2013cu through a superwind

The progenitors of type-IIb supernovae (SNe) are believed to have lost their H-rich envelopes almost completely in the direct pre-SN phase. Recently the first 'flash spectrum' of a SN IIb (SN2013cu) has been presented, taken early enough to study its immediate circumstellar medium (CSM). Similar to a previous study by Groh (2014) we analyse the structure and chemical composition of the optically-thick CSM using non-LTE model atmospheres. For the first time we take light-travel time (LTT) effects on the spectrum formation into account, which affect the shapes and strengths of the observable emission lines, as well as the inferred SN luminosity. Based on the new CSM parameters we estimate a lower limit of ~0.3Msun for the CSM mass, which is a factor 10-100 higher than previous estimates. The spectral fit implies a CSM in the form of a homogeneous and spherically symmetric superwind whose mass-loss rate exceeds common expectations by up to two orders of magnitude. The derived chemical composition is in agreement with a progenitor that has just left, or is just about to leave the Red-Supergiant (RSG) stage, confirming the standard picture for the origin of SNe IIb. Due to its extreme mass loss the SN progenitor will likely appear as extreme RSG, Luminous Blue Variable (LBV), or Yellow Hypergiant (YHG). The direct detection of a superwind, and the high inferred CSM mass suggest that stellar wind mass loss may play an important role in the formation of SNe IIb.

A non-LTE analysis of the hot subdwarf O star BD+28°4211

We present the second part of our detailed analysis of the hot sdO and spectroscopic standard star BD+28 4211, in which we focus on the optical spectrum. In the first part of our study, we determined the abundances of some 11 metals detected in the atmosphere of BD+28 4211 using UV spectra of the star and corroborated the fundamental parameters estimated in past studies (Teff $\sim$ 82,000 K, log g $\sim$ 6.2, and solar N(He)/N(H)). In this work, we aim at rederiving these secured parameters on the sole basis of high-quality optical spectra. A first grid of non-LTE line-blanketed model atmospheres, including metals with the abundances derived from the UV spectrum, does not give satisfactory results when we apply a standard simultaneous fitting procedure to the observed H and He lines of our optical spectra. The line profiles are not finely reproduced and the resulting effective temperatures, in particular, are too low by $\sim$10,000 K. We next investigate the probable cause of this failure, that is, the importance of missing opacity sources on the atmospheric stratification. We compare line profiles computed from models with artificially boosted metallicities, from solar abundances to 15$\times$ these values. We find that the structural effects saturate for a metallicity of $\sim$10x solar, and use this to compute a second full grid of models and synthetic spectra. This metal-enriched grid allows us to achieve significantly improved spectral fits with models having the expected parameters. Our test case thus reveals that there is still a need for models with enhanced metallicity for better estimating the atmospheric parameters of objects such as hot subdwarfs and hot white dwarfs if only optical spectra are available.

The VLT-FLAMES Tarantula Survey

TLUSTY non-LTE model atmosphere calculations have been used to determine atmospheric parameters and nitrogen (N) abundances for 34 single and 18 binary B-type supergiants (BSGs). The effects of flux contribution from an unseen secondary were considered for the binary sample. We present the first systematic study of the incidence of binarity for a sample of BSGs across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the BSGs it may be necessary to extend the TAMS to lower temperatures. This is consistent with the derived distribution of mass discrepancies, projected rotational velocities (vsini) and N abundances, provided that stars cooler than this temperature are post RSG objects. For the BSGs in the Tarantula and previous FLAMES surveys, most have small vsini. About 10% have larger vsini (>100 km/s) but surprisingly these show little or no N enhancement. All the cooler BSGs have low vsini of <70km/s and high N abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. A lack of cool binaries, possibly reflects the small sample size. Single star evolutionary models, which include rotation, can account for the N enhancement in both the single and binary samples. The detailed distribution of N abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. The first comparative study of single and binary BSGs has revealed that the main sequence may be significantly wider than previously assumed, extending to Teff=20000K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolution of BSGs. Full abstract in paper

Analysis of cool DO-type white dwarfs from the Sloan Digital Sky Survey data release 10

We report on the identification of 22 new cool DO-type white dwarfs (WD) detected in Data Release 10 (DR10) of the Sloan Digital Sky Survey (SDSS). Among them, we found one more member of the so-called hot-wind DO WDs, which show ultrahigh excitation absorption lines. Our non-LTE model atmosphere analyses of these objects and two not previously analyzed hot-wind DO WDs, revealed effective temperatures and gravities in the ranges Teff=45-80kK and log g= 7.50-8.75. In eight of the spectra we found traces of C (0.001-0.01, by mass). Two of these are the coolest DO WDs ever discovered that still show a considerable amount of C in their atmospheres. This is in strong contradiction with diffusion calculations, and probably, similar to what is proposed for DB WDs, a weak mass-loss is present in DO WDs. One object is the most massive DO WD discovered so far with a mass of 1.07 M_sun if it is an ONe-WD or 1.09 M_sun if it is a CO-WD. We furthermore present the mass distribution of all known hot non-DA (pre-) WDs and derive the hot DA to non-DA ratio for the SDSS DR7 spectroscopic sample. The mass distribution of DO WDs beyond the wind limit strongly deviates from the mass distribution of the objects before the wind limit. We address this phenomenon by applying different evolutionary input channels. We argue that the DO WD channel may be fed by about 13% by post-extreme-horizontal branch stars and that PG1159 stars and O(He) stars may contribute in a similar extent to the non-DA WD channel.

Weak metal lines in optical high-resolution Very Large Telescope and Keck spectra of “cool” PG 1159 stars

PG 1159 stars are very hot (effective temperatures Teff = 75 000–200 000 K), hydrogen-deficient (pre-) white dwarfs. They probably are the result of a late helium-shell flash that laid bare the He, C, and O rich intershell matter of the progenitor Asymptotic Giant Branch (AGB) star. Their chemical surface composition thus allows to conclude on details of AGB-star nucleosynthesis. Due to their very high effective temperatures, detailed spectral analyses are usually completely reliant on ultraviolet observations, except for some species in the hottest PG 1159 stars (Teff > 130 000 K), which do exhibit highly excited lines from the CNO elements and neon (C iv ,N v ,O vi ,N evii-viii) in optical spectra. Particularly problematic are, however, the coolest members of the PG 1159 class that exclusively show C iv lines in the optical. Access to the nitrogen abundance is important to decide which of the late-thermal pulse evolutionary scenarios was experienced by a particular star, while a high oxygen abundance is an important marker that the star could pulsate. In the present paper, we investigate high-resolution high signal-to-noise optical spectra of three “cool” PG 1159 stars (PG 0122+200, PG 2131+066, MCT 0130−1937, Teff = 80 000–95 000 K). With the help of non-LTE model atmospheres and synthetic spectra, we are able to identify a large number of weak CNO lines (C iii ,N iv ,O iii-v) that were not detected before in these stars. They allow abundance determinations and enable us to constrain the effective temperature to high precision through ionization equilibria without the requirement to access the ultraviolet spectral range.

The rapid evolution of the exciting star of the Stingray nebula

SAO244567, the exciting star of the Stingray nebula, is rapidly evolving. Previous analyses suggested that it has heated up from an effective temperature of about 21kK in 1971 to over 50kK in the 1990s. Canonical post-asymptotic giant branch evolution suggests a relatively high mass while previous analyses indicate a low-mass star. Fitting line profiles from static and expanding non-LTE model atmospheres to the observed UV and optical spectra, taken during 1988-2013, allowed us to study the temporal change of effective temperature, surface gravity, mass-loss rate, and terminal wind velocity. In addition, we determined the chemical composition of the atmosphere. We find that the central star has steadily increased its effective temperature from 38kK in 1988 to a peak value of 60kK in 2002. During the same time, the star was contracting, as concluded from an increase in surface gravity from log g = 4.8 to 6.0 and a drop in luminosity. Simultaneously, the mass-loss rate declined from log (dM/dt/Msun/yr)=-9.0 to -11.6 and the terminal wind velocity increased from 1800km/s to 2800km/s. Since around 2002, the star stopped heating and has cooled down again to 55kK by 2006. It has a largely solar surface composition with the exception of slightly subsolar carbon, phosphorus, and sulfur. By comparison with stellar-evolution calculations, we confirm that SAO244567 must be a low-mass star (M < 0.55 Msun). However, the slow evolution of the respective stellar evolutionary models is in strong contrast to the observed fast evolution and the young planetary nebula with a kinematical age of only about 1000 years. We speculate that the star could be a late He-shell flash object. Alternatively, it could be the outcome of close-binary evolution. Then SAO244567 would be a low-mass (0.354 Msun) helium prewhite dwarf after the common-envelope phase, during which the planetary nebula was ejected.

The impact of rotation on the line profiles of Wolf-Rayet stars

Massive Wolf-Rayet stars are recognized today to be in a very common, but short, evolutionary phase of massive stars. While our understanding of Wolf-Rayet stars has increased dramatically over the past decades, it remains unclear whether rapid rotators are among them. There are various indications that rapidly rotating Wolf-Rayet stars should exist. Unfortunately, due to their expanding atmospheres, rotational velocities of Wolf-Rayet stars are very difficult to measure. However, recently observed spectra of several Wolf-Rayet stars reveal peculiarly broad and round emission lines. Could these spectra imply rapid rotation? In this work, we model the effects of rotation on the atmospheres of Wolf-Rayet stars. We further investigate whether the peculiar spectra of five Wolf-Rayet stars may be explained with the help of stellar rotation, infer appropriate rotation parameters, and discuss the implications of our results. We make use of the Potsdam Wolf-Rayet (PoWR) non-LTE model atmosphere code. Since the observed spectra of Wolf-Rayet stars are mainly formed in their expanding atmospheres, rotation must be accounted for with a 3D integration scheme of the formal integral. For this purpose, we assume a rotational velocity field consisting of an inner co-rotating domain and an outer domain, where the angular momentum is conserved. We find that rotation can reproduce the unique spectra analyzed here. However, the inferred rotational velocities at the stellar surface are large (~200 km/s), and the inferred co-rotation radii (~10 stellar radii) suggest the existence of very strong photospheric magnetic fields (~20 kG).

Macroclumping as solution of the discrepancy between Hαand P v mass loss diagnostics for O-type stars

Recent studies of O-type stars 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 H{\alpha} emission. Wind clumping may be the main cause for this discrepancy. In a previous paper, we have presented 3-D 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 O4 to O7) and test whether the reported discrepancies can be resolved this way. In the first step, the analyses start with simulating the observed spectra with Potsdam Wolf-Rayet (PoWR) non-LTE model atmospheres. The mass-loss rates are adjusted to fit best to the observed H{\alpha} emission lines. For the unsaturated UV resonance lines (i.e. P v) we then apply our 3-D Monte-Carlo code, which can account for wind clumps of any optical depths, a non-void inter-clump medium, and a velocity dispersion inside the clumps. The ionization stratifications and underlying photospheric spectra are adopted from the PoWR models. From fitting the observed resonance line profiles, the properties of the wind clumps are constrained. Our results show that with the mass-loss rates that fit H{\alpha} (and other Balmer and He II lines), the UV resonance lines (especially the unsaturated doublet of P v) can also be reproduced without problem when macroclumping is taken into account. There is no need to artificially reduce the mass-loss rates, nor to assume a sub-solar phosphorus abundance or an extremely high clumping factor, contrary to 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. Macroclumping resolves the previously reported discrepancy between H{\alpha} and P v mass-loss diagnostics.

Orbital and physical parameters of the spectroscopic binary HD37737

We report the physical and orbital parameters of the visible component of the spectroscopic binary HD37737 (m V = 8.03). The observations were performed with the 1.2-m telescope of the Kourovka Astronomical Observatory of the Ural Federal University in 2012 and the 6-m BTA telescope of the SAO RAS in 2007 and 2009. Radial velocities were measured separately from each spectral line of the list by the cross-correlation method with a synthetic spectrum. The latter was calculated using the grids of non-LTE model atmospheres with solar chemical compositions. A significant difference in the epochs of observations (2005–2012) allowed to refine the orbital period of the star (7 · d 84705) and the orbital elements of the binary system. We obtained an estimate of the mass function f(m) = 0.23 ± 0.02M ⊙. The best agreement between the synthetic and observed spectra is achieved at T eff = 30 000 K and log g = 3.50 according to the observations on both instruments. The obtained parameters correspond to a star of spectral type O9.5 III, with mass estimated at 26 ± 2M ⊙. The minimum mass estimate of the secondary component of the binary is 6.2 ± 0.5M ⊙. We have discovered a fact that the velocities, obtained from different spectral lines, differ, which is typical for giant stars. Engaging additional spectra, obtained in 2005 with the 2.1-m KPNO telescope, we investigated the effect of this fact on the estimate of the speed of the system’s center of mass. The difference in the velocities of various lines is approximately the same in the spectra, obtained at all the three instruments. The obtained ratios suggest that the deeper layers of the atmosphere of the star are moving with a greater velocity than the outer layers. Depending on the line, the estimate of the heliocentric velocity of the binary’s center of mass varies in the range from −11 to 1 km/s.

Testing the opacity and equation of state of LTE and non-LTE model atmospheres with OPAL and OP data for early-type stars

Context. Complex investigations of stars are made by studying their atmospheres, evolutionary states and oscillation properties, but the opacities and equation of state are still uncertain. Aims. Compatibility of model atmospheres with internal structure of early-type stars are investigated at photospheric and subphotospheric regions. Methods. The problem was studied quantitatively by means of diagrams that involve the Rosseland-mean opacities, temperatures, densities, and gradients of radiation pressure as functions of depth using stellar atmosphere and envelope models. Results. Two new grids of radiative-equilibrium models of atmospheres were calculated assuming the local thermodynamic equilibrium (LTE). The first one is based on Kurucz’s ATLAS12 computer code implemented with the occupation probability formalism, which accounts for the destruction of loosely bound states of atoms by interactions with particles in the plasma. The second grid of models is based on monochromatic opacities and equations of state taken from the well-known international opacity project (OP). Non-LTE model atmospheres are also calculated and some aspects of the modeling are pointed out to improve agreement with ATLAS12, OP data, as well as with results of the Rogers & Iglesias OPAL opacity computing code (OPAL). The model atmospheres are also discussed by analyzing emerging fluxes of radiation. Conclusions. Commonly used LTE and non-LTE models of stellar atmospheres of early-type stars differ markedly from each other and do not fit OPAL or OP stellar envelope models at great optical depths. The results of the OP project distributed as the OPCD v 3.3 base are useful for calculating not only stellar envelopes but also fully line-blanketed LTE model atmospheres. These models have diagnostic values for studying the atomic physics used for modeling of photospheric and sub-photospheric regions of stars. The OPAL and OP opacities are markedly underestimated in comparison with the Rosseland-mean opacities taken from the Castelli & Kurucz (2003, IAU Symp., 210, 10) atmosphere models in which a new opacity bump appears at lgT ≈ 5.06. This additional opacity bump affects the OPAL- and OP-driving zone for stellar pulsations and therefore the new envelope models may markedly change predicted spectra of unstable oscillation modes. Differences between current non-LTE and LTE models based on the TLUSTY200 and ATLAS codes, respectively, cannot be assigned entirely to non-LTE effects, but for B stars they reflect mainly opacity effects.