Massive stars

We describe the present state of massive star research seen from the viewpoint of stellar evolution, with special emphasis on close binaries. Statistics of massive close binaries are reasonably complete for the Solar neighbourhood. We defend the thesis that within our knowledge, many scientific results where the effects of binaries are not included, have an academic value, but may be far from reality. In chapter I, we summarize general observations of massive stars where we focus on the HR diagram, stellar wind mass loss rates, the stellar surface chemistry, rotation, circumstellar environments, supernovae. Close binaries can not be studied separately from single stars and vice versa. First, the evolution of single stars is discussed (chapter I). We refer to new calculations with updated stellar wind mass loss rate formalisms and conclusions are proposed resulting from a comparison with representative observations. Massive binaries are considered in chapter II. Basic processes are briefly described, i.e. the Roche lobe overflow and mass transfer, the common envelope process, the spiral-in process in binaries with extreme mass ratio, the effects of mass accretion and the merging process, the implications of the (asymmetric) supernova explosion of one of the components on the orbital parameters of the binary. Evolutionary computations of interacting close binaries are discussed and general conclusions are drawn. The enormous amount of observational data of massive binaries is summarized. We separately consider the non-evolved and evolved systems. The latter class includes the semi-detached and contact binaries, the WR binaries, the X-ray binaries, the runaways, the single and binary pulsars. A general comparison between theoretical evolution and observations is combined with a discussion of specially interesting binaries: the evolved binaries HD 163181, HD 12323, HD 14633, HD 193516, HD 25638, HD 209481, Per and Sgr; the WR+OB binary V444 Cyg; the high mass X-ray binaries Vela X-1, Wray 977, Cyg X-1; the low mass X-ray binaries Her X-1 and those with a black hole candidate; the runaway Pup, the WR+compact companion candidates Cyg X-3, HD 50896 and HD 197406. We finally propose an overall evolutionary model of massive close binaries as a function of primary mass, mass ratio and orbital period. Chapter III deals with massive star population synthesis with a realistic population of binaries. We discuss the massive close binary frequency, mass ratio and period distribution, the observations that allow to constrain possible asymmetries during the supernova explosion of a massive star. We focuss on the comparison between observed star numbers (as a function of metallicity) and theoretically predicted numbers of stellar populations in regions of continuous star formation and in starburst regions. Special attention is given to the O-type star/WR star/red supergiant star population, the pulsar and binary pulsar population, the supernova rates. Received 17 July 1998     

共生星双星的星风物质损失率的测定

共生星双星是一颗有强大星风物质损失的红巨星与一颗早型热星组成的特殊双星系统。由于早型热星在充满红巨星的星风物质的空间中环绕运行，可以产生Ｐ—Ｃｙｇｎｉ型谱线，通过对Ｐ—Ｃｙｇｎｉ型谱线的理论分析可以精确测定共生星双星的星风物质损失率。文章介绍了国际上比较有代表意义的一些共生星双星的工作，其中的方法几乎都是近十年中发展起来的。同时，文章还介绍了作者的一些工作。在谱线形成计算中考虑了较多的因素，如氢—氦混合气体的多能级跃迁问题、轨道运动引起的密度非径向分布问题等，并在轨道形状方面做了一些简化。反映了共生星双星谱线形成关键的周期性相位变化的特征，取得了比较满意的结果，对这个方法存在的问题和改进方向进行了一些简要讨论，此外，还介绍了作者在线性化分离法求解Ｎｏｎ－ＬＴＥ大气模型中所做的工作。     

钡星系统轨道根数分布及丰度的Monte-Carlo模拟计算

采用星风质量吸积的角动量守恒模型，用Monte—Carlo方法研究了普通红巨星双星系统和钡星的轨道根数的变化规律，由于钡星系统是由普通红巨星双星系统演化而来，因此钡星系统的轨道偏心率及周期的分布显示了经过质量吸积后双星系统的最终轨道特征。计算结果表明，随着星风吸积过程的进行，在星风质量损失阶段系统轨道半长轴将增大，导致轨道周期增大，而偏心率变化不大，由此可以解释普通红巨星双星系统和钡星系统的轨道根数的分布规律和变化情况以及钡星重元素丰度分布特征。     

Variable stars in the open cluster NGC 2141

We report the results of a search for variable stars in the open cluster NGC 2141. Ten variable stars are detected, among which nine are new variable stars and they are classified as three short-period W UMa-type eclipsing binaries, two EAtype eclipsing binaries, one EB-type eclipsing binary, one very short-period RS CVntype eclipsing binary, one d-type RR Lyrae variable star, and one unknown type of variable star. The membership and physical properties are discussed, based on their light curves, positions in the color magnitude diagrams, spatial locations and periods.A known EB-type eclipsing binary is also identified as a blue straggler candidate in the cluster. Furthermore, we find that all eclipsing contact binaries have prominent asymmetric eclipses and display the O'Connell effect, which increases with a decrease in orbital periods. This suggests that the O'Connell effect is probably related to the evolution of the orbital period in short period eclipsing binary systems.     

X-ray binaries

Summary The various types and classes of X-ray binary are reviewed high-lighting recent results. The high mass X-ray binaries (HMXRBs) can be used to probe the nature of the mass loss from the OB star in these systems. Absorption measurements through one orbital cycle of the supergiant system X1700-37 are well modelled by a radiation driven wind and also require a gas stream trailing behind the X-ray source. In Cen X-3 the gas stream is accreted by the X-ray source via an accretion disk. Changes in the gas stream can cause the disk to thicken and the disk to obscure the X-ray source. How close the supergiant is to corotation seems to be as much a critical factor in these systems as how close it is to filling its Roche lobe. In the Be star X-ray binaries a strong correlation between the neutron stars rotation period and its orbital period has been explained as due to the neutron star being immersed in a dense, slow moving equatorial wind from the Be star. For the X-ray pulsars in the transient Be X-ray binaries a centrifugal barrier to accretion is important in determining the X-ray lightcurve and the spin evolution. The X-ray orbital modulations from the low mass X-ray binaries, LMXRBs, include eclipses by the companion and/or periodic dipping behaviour from structure at the edge of the disk. The corresponding optical modulations show a smooth sinusoidal like component and in some cases a sharp eclipse by the companion. The orbital period of the LMXRB XB1916-05 is 1% longer in the optical compared to that given by the X-ray dip period. The optical period has been interpreted as the orbital period, but this seems inconsistent with the well established view of the origin of the X-ray modulations in LMXRB. A new model is presented that assumes the X-ray dip period is the true orbital period. The 5.2 h eclipsing LMXRB XB2129+47 recently entered a low state and optical observations unexpectedly reveal an F star which is too big to fit into the binary. This is probably the first direct evidence that an X-ray binary is part of a hierarchical triple. Finally the class of X-ray binaries containing black hole candidates is reviewed focusing on the value of using X-ray signatures to identify new candidates.     

Modelling the spectra of colliding winds in the Wolf–Rayet WC7+O binaries WR 42 and WR 79

We have obtained complete phase coverage of the WC7+O binaries WR 42 = HD 97152 and WR 79 = HD 152270 with high signal-to-noise ratio (S/N), moderate-resolution spectra. Remarkable orbital phase-locked profile variations of the C  iii λ 5696 line are observed and interpreted as arising from colliding wind effects. Within this scenario, we have modelled the spectra using a purely geometrical model that assumes a cone-shaped wind–wind interaction region which partially wraps around the O star. Such modelling holds the exciting promise of revealing a number of interesting parameters for WR+O binaries, such as the orbital inclination, the streaming velocity of material in the interaction region and the ratio of wind momentum flux. Knowledge of these parameters in turn leads to the possibility of a better understanding of WR star masses, mass-loss rates and wind region characteristics.     

Orbital period modulation and magnetic cycles in close binaries

We discuss the observed orbital period modulations in close binaries, and focus on the mechanism proposed by Applegate relating the changes of the stellar internal rotation associated with a magnetic activity cycle with the variation of the gravitational quadrupole moment of the active component; the variation of this quadrupole moment in turn forces the orbital motion of the binary stars to follow the activity level of the active star. We generalize this approach by considering the details of this interaction, and develop some illustrative examples in which the problem can be easily solved in analytical form. Starting from such results, we consider the interplay between rotation and magnetic field generation in the framework of different types of dynamo models, which have been proposed to explain solar and stellar activity. We show how the observed orbital period modulation in active binaries may provide new constraints for discriminating between such models. In particular, we study the case of the prototype active binary RS Canum Venaticorum, and suggest that torsional oscillations — driven by a stellar magnetic dynamo — may account for the observed behaviour of this star. Further possible applications of the relationship between magnetic activity and orbital period modulation, related to the recent discovery of binary systems containing a radio pulsar and a convecting upper main-sequence or a late-type low-mass companion, are discussed.     

A possible detection mechanism for pre-low-mass X-ray binaries

This work presents a possible detection mechanism for close, detached, neutron star–red dwarf binaries, which are expected to be the evolutionary precursors of low-mass X-ray binaries (LMXBs). Although this pre-low-mass X-ray binary (pre-LMXB) phase of evolution is predicted theoretically, as yet no such systems have been identified observationally. The calculations presented here suggest that the X-ray luminosity of neutron star wind accretion in a pre-LMXB system can be expected to exceed the intrinsic X-ray luminosity of the red dwarf secondary star. Furthermore, the temperature of the radiation emitted from the neutron star wind accretion process is expected, within the confines of a reasonable set of conditions, to lie within the detection range of X-ray satellites. Sources with X-ray luminosities greater than that expected for a red dwarf star, but the positions of which coincide with that of a red dwarf star, are then candidate pre-LMXB systems. These candidate systems should be surveyed for the radial velocity shifts that would occur as a result of the orbital motion of a red dwarf star within a close binary system containing a high-mass compact object.     

Are the W Ursae Majoris-type systems EK Comae Berenices and UX Eridani surrounded by circumstellar matter?

The variations of the orbital periods of two nearly neglected W UMa-type eclipsing binaries, EK Comae Berenices and UX Eridani, are presented through a detailed analysis of the O–C diagrams. It is found that the orbital period of EK Com is decreasing and the period of UX Eridani is increasing, and several sudden jumps have occurred in the orbital periods of both binaries. We analyze the mechanism(s), which might underlie the changes of the orbital periods of both systems, and obtain some new results. The long-term decrease of the orbital period of EK Comae Berenices might be caused by the decrease of the orbital angular momentum due to a magnetic stellar wind (MSW) or by mass transfer from the more massive to the less massive component. The secular increase in the orbital period of UX Eridani might be caused by mass transfer from the less massive to the more massive star. The possible mechanisms, which underlie the sudden changes in the orbital periods of the close binary systems are as the followings: (1) the variations of the structure due to the variation of the magnetic field; (2) the rapid mass exchange between the close binaries and their circumstellar matter. Finally, the evolutionary status of the systems EK Comae Berenices and UX Eridani is discussed.     

The orbital period distribution of subdwarf-B binaries

We present the results of a 3.5 year long campaign to measure orbital periods of subdwarf-B (sdB) star binaries. We directly compare our observed orbital period distribution with that predicted by using binary population synthesis. Up to now, most of our systems seem to have been formed through two of the formation channels discussed by Han et al. (2003, MNRAS 341, p. 669), i.e. the first and the second common envelope ejection (CE) channels. At this point, thanks to the long baseline of our observations, we are starting to detect also very long orbital period systems. These have probably come from a complete different formation path, the first stable Roche Lobe overflow (RLOF) channel in which the first mass transfer phase is stable. This channel is expected to lead to the formation of very wide binaries with typical orbital periods ranging from 1 month to 1 year.     

Evolution of massive binary stars in the LMC and its implications for radio pulsar population

The Hertzsprung-Russell diagram of the Large Magellanic Cloud compiled recently by Fitzpatrick & Garmany (1990) shows that there are a number of supergiant stars immediately redward of the main sequence although theoretical models of massive stars with normal hydrogen abundance predict that the region 4.5 ≤ logT _eff ≤ 4.3 should be un-populated (“gap”). Supergiants having surface enrichment of helium acquired for example from a previous phase of accretion from a binary companion, however, evolve in a way so that the evolved models and observed data are consistent — an observation first made by Tuchman & Wheeler (1990). We compare the available optical data on OB supergiants with computed evolutionary tracks of massive stars of metallicity relevant to the LMC with and without helium-enriched envelopes and conclude that a large fraction (≈ 60 per cent) of supergiant stars may occur in binaries. As these less evolved binaries will later evolve into massive X-ray binaries, the observed number and orbital period distribution of the latter can constrain the evolutionary scenarios of the supergiant binaries. The distributions of post main sequence binaries and closely related systems like WR + O stars are bimodal-consisting of close and wide binaries in which the latter type is numerically dominating. When the primary star explodes as a supernova leaving behind a neutron star, the system receives a kick and in some cases can lead to runaway O-stars. We calculate the expected space velocity distribution for these systems. After the second supernova explosion, the binaries in most cases, will be disrupted leading to two runaway neutron stars. In between the two explosions, the first born neutron star’s spin evolution will be affected by accretion of mass from the companion star. We determine the steady-state spin and radio luminosity distributions of single pulsars born from the massive stars under some simple assumptions. Due to their great distance, only the brightest radio pulsars may be detected in a flux-limited survey of the LMC. A small but significant number of observable single radio pulsars arising out of the disrupted massive binaries may appear in the short spin period range. Most pulsars will have a low velocity of ejection and therefore may cluster around the OB associations in the LMC.     

Post-common-envelope binaries from SDSS – I. 101 white dwarf main-sequence binaries with multiple Sloan Digital Sky Survey spectroscopy

We present a detailed analysis of 101 white dwarf main-sequence binaries (WDMS) from the Sloan Digital Sky Survey (SDSS) for which multiple SDSS spectra are available. We detect significant radial velocity variations in 18 WDMS, identifying them as post-common-envelope binaries (PCEBs) or strong PCEB candidates. Strict upper limits to the orbital periods are calculated, ranging from 0.43 to 7880 d. Given the sparse temporal sampling and relatively low spectral resolution of the SDSS spectra, our results imply a PCEB fraction of ≳15 per cent among the WDMS in the SDSS data base. Using a spectral decomposition/fitting technique we determined the white dwarf effective temperatures and surface gravities, masses and secondary star spectral types for all WDMS in our sample. Two independent distance estimates are obtained from the flux-scaling factors between the WDMS spectra, and the white dwarf models and main-sequence star templates, respectively. Approximately one-third of the systems in our sample show a significant discrepancy between the two distance estimates. In the majority of discrepant cases, the distance estimate based on the secondary star is too large. A possible explanation for this behaviour is that the secondary star spectral types that we determined from the SDSS spectra are systematically too early by one to two spectral classes. This behaviour could be explained by stellar activity, if covering a significant fraction of the star by cool dark spots will raise the temperature of the interspot regions. Finally, we discuss the selection effects of the WDMS sample provided by the SDSS project.     

Speckle interferometry of nearby multiple stars. IV. Measurements in 2004 and new orbits

The results of speckle interferometric observations of 104 binary and 6 triple stars performed at the BTA 6 m telescope in 2004 October are presented. Nearby low-mass stars are mostly observed for the program, among which 59 there are new binaries recently discovered by the Hipparcos astrometric satellite. Concurrently with the diffraction-limited position measurements we obtained 154 brightness ratio measurements of binary and multiple star components in different bands of the visible spectrum. New, first-resolved binaries are the symbiotic star CH Cyg with a weak companion at 0.043″ separation and the pair of red dwarfs, GJ 913 = HIP 118212. In addition, we derived the orbital parameters for two interferometric systems: the CN-giant pair HD 210211 = HIP 109281 (P = 10.7 yr) and the G2V-K2V G2V-K2V binary GJ 9830 = HIP 116259 (P = 15.7 yr).     

A New Algorithm in Computations of Binary Stars’ Observational Quantities and Its Application

The determination of the orbital parameters of binaries is an important constituent of fundamental astronomy. It is the only reliable way to determine the stellar mass, which is crucial for the studies of stellar dynamics and evolution. And the orbital parameters can be used to compute the motions of component stars, which could improve the population of reference stars in the bright star catalogue reference frame. However, the previous studies on the determination of binaries’ orbital parameters adopted some approximations in the computation of observational quantities. According to the statistical results of hundreds of binaries with bright component stars, it could be found that the approximations in the variation of the tangent plane and the prospective effect of the reference plane have made the deviations of computed observational quantities become larger than the present observational precision (1 mas). In this paper, a new algorithm based on the two-body model and the rigorous relation of solid geometry for computing observational quantities is presented. This new algorithm has been verified by fitting the simulation and actually-observed data. And the fitting results show that the fitted parameters solved by the new algorithm usually have a higher confidence than that solved by the approximate algorithm. Therefore, the new algorithm has certain improvements in the determination of the orbital parameters of binaries, especially, for the wide binaries nearby the solar system.     

The period gap and the AM Herculis type systems

We argue that the period distribution of AM Herculis binaries in the enlarged sample incorporating results from the recent ROSAT X-ray survey differs significantly from that of other cataclysmic variables. In particular, there is no evidence for a pronounced period gap at 2—3 hr and the significance of the period spike at about 114 min is brought into question. We present an alternative evolutionary scenario for the AM Herculis binaries based on the hypothesis that magnetic braking by the stellar wind of the M star secondary either ceases or is severely curtailed when the rotation of the magnetic white dwarf becomes synchronised with the orbital motion, The orbital evolution of the AM Herculis binaries is thereafter driven mainly by angular momentum loss due to gravitational radiation. This scenario not only explains the higher proportion of AM Herculis binaries in the period gap when compared with other cataclysmic variables but also provides a natural explanation for the low mass transfer rates in these binaries and for the existence of an apparent upper limit for the surface magnetic fields of their white dwarfs.     

Population Synthesis of Binary Be Stars with Degenerate Companions in the Galaxy

Using the numerical code (`Scenario Machine') we study of number and physical properties of binary Be stars. Evolutionary tracks leading to a formation of the observational binary systems are presented. We conclude that synchronization must be taken into account when calculating binary Be star evolution and calculate the minimal orbital period for Be/evolved companion binary. The obtained distributions over orbital parameters are in good agreement with the observational lack of short-period Be/X-ray binaries. According to our calculations 70% of all Be stars must have a white dwarf. The white dwarfs in these systems should be hot enough with the surface temperature distribution peaking at 10000–20000 K. Their detection is possible during the period of the lack of Be star envelope by the detection of white dwarf extremely UV and soft X-ray emission. This method of registration appears to be particularly promising for `single' early-type Be stars because in these systems the white dwarfs must have a very high surface temperature. However, the loss of the Be disc-like envelope does not often occur and it is a rather rare event for many Be stars. The best possibility of white dwarf detection is given by the study of helium spectral lines found in emission from several Be stars. The ultraviolet continuum energy of these Be stars is found to be not enough to produce the observed helium emission. Besides, we also discuss the orbital properties of binary Be star systems with other evolved companions such as helium stars and neutron stars and give a possible explanation for the lack of Be/black hole binaries. This revised version was published online in July 2006 with corrections to the Cover Date.     

Do the physical properties of Ap binaries depend on their orbital elements?

We reveal sufficient evidence that the physical characteristics of Ap stars are related to binarity. The Ap star peculiarity [represented by the  Δ( V 1- G )  value and magnetic field strength] diminishes with eccentricity, and it may also increase with orbital period ( P _orb). This pattern, however, does not hold for large orbital periods. A striking gap that occurs in the orbital period distribution of Ap binaries at 160–600 d might well mark a discontinuity in the above-mentioned behaviour. There is also an interesting indication that the Ap star eccentricities are relatively lower than those of corresponding B9–A2 normal binaries for   P _orb>10 d  . All this gives serious support to the pioneering idea of Abt & Snowden concerning a possible interplay between the magnetism of Ap stars and their binarity. Nevertheless, we argue instead in favour of another mechanism, namely that it is binarity that affects magnetism and not the opposite, and suggest the presence of a new magnetohydrodynamical mechanism induced by the stellar companion and stretching to surprisingly large P _orb.     

Orbital period changes of contact binary systems: direct evidence for thermal relaxation oscillation theory

Orbital period changes of ten contact binary systems (S Ant, ε CrA, EF Dra, UZ Leo, XZ Leo, TY Men, V566 Oph, TY Pup, RZ Tau and AG Vir) are studied based on the analysis of their     curves. It is discovered that the periods of the six systems, S Ant, ε CrA, EF Dra, XZ Leo, TY Men and TY Pup, show secular increases. For UZ Leo, its secular period increase rate is revised. For the three systems, V566 Oph, RZ Tau and AG Vir, weak evidence is presented that a periodic oscillation (with periods of 20.4, 28.5 and 40.9 yr respectively) is superimposed on a secular period increase. The cyclic period changes can be explained by the presence of an unseen third body in the three systems. All the sample stars studied are contact binaries with     .
Furthermore, orbital period changes of 27 hot contact binaries have been checked. It is found that, apart from AW UMa with the lowest mass ratio     , none shows an orbital period decrease. The relatively weak magnetic activity in the hotter contact binaries means little angular momentum loss (AML) from the systems via magnetic stellar winds. The period increases of these W UMa binaries can be explained by mass transfer from the secondary to the primary components, which is in agreement with the prediction of the thermal relaxation oscillation (TRO) models. This suggests that the evolution of a hotter W UMa star is mainly controlled by TRO. On the other hand, for a cooler W UMa star     , its evolution may be TRO plus AML, which coincides with the recent results of Qian.     

Companions to sdB binaries−degenerate or non-degenerate?

We report the results of an on-going programme to collect lightcurves for EHB binaries with orbital periods P ? 1 d. Degenerate and non-degenerate companions are clearly distinguished by the presence or absence of a reflection effect in these binaries. The amplitude of the reflection effect can be combined with other data to provide an estimate of the mass of the companion star. We find that the fraction of sdB binaries in our sample which have non-degenerate companions is 0.08±0.06 (1?σ error). These non-degenerate companions have very low masses (? 0.1M⊙). This property is not predicted by existing population synthesis models.     

A photometric study of two recently discovered,contact binaries-II: UCAC4 479-113658 and UCAC4 479-113711

In a previous paper, we have presented photometric studies of four recently discovered contact binaries (Djurašević, et al., 2016; hereafter Paper I) that were observed by the 1.88 m telescope of the Kottamia Astronomical Observatory (KAO) in Egypt. In this paper we continue studying two more recently discovered, by some of us, contact binaries UCAC4 479-113658 and UCAC4 479-113711, and analyze their light curves to determine orbital and physical parameters using the modeling program of G. Djurašević as in Paper I. The analysis shows that the two systems are W UMa contact binaries of late type main-sequence stars with a bright spot on the secondary less massive cool star of the first system; and spot on each component of the later system. We also conclude that both systems are most likely to be of the A-type class.