WR + OB Binary Systems: Observational Evidence of Their Formation as a Result of Mass Exchange

An analysis of observational data shows that, in most cases,Wolf–Rayet (WR) stars in known WR+ OB binary systems were formed as a result of mass transfer in initial OB + OB systems, rather than through radial mass loss by the more massive OB star via its stellar wind.     

Wolf-Rayet stars,black holes,and gamma-ray bursters in close binaries

We consider the evolutionary status of observed close binary systems containing black holes and Wolf-Rayet (WR) stars. When the component masses and the orbital period of a system are known, the reason for the formation of a WR star in an initial massive system of two main-sequence stars can be established. Such WR stars can form due to the action of the stellar wind from a massive OB star (M_OB≥50M_⊙), conservative mass transfer between components with close initial masses, or the loss of the common envelope in a system with a large (up to ～25) initial component mass ratio. The strong impact of observational selection effects on the creation of samples of close binaries with black holes and WR stars is demonstrated. We estimate theoretical mass-loss rates for WR stars, which are essential for our understanding the observed ratio of the numbers of carbon and nitrogen WR stars in the Galaxy \(\dot M_{WR} (M_ \odot yr^{ - 1} ) = 5 \times 10^{ - 7} (M_{WR} /M_ \odot )^{1.3} \). We also estimate the minimum initial masses of the components in close binaries producing black holes and WR stars to be ～25M_⊙. The spatial velocities of systems with black holes indicate that, during the formation of a black hole from a WR star, the mass loss reaches at least several solar masses. The rate of formation of rapidly rotating Kerr black holes in close binaries in the Galaxy is ～3×10^?6 yr^?1. Their formation may be accompanied by a burst of gamma radiation, possibly providing clues to the nature of gamma-ray bursts. The initial distribution of the component mass ratios for close binaries is dN～dq=dM₂/M₁ in the interval 0.04?q₀≤1, suggesting a single mechanism for their formation.     

Circumstellar nebulae of WR and OB stars with strong winds: Simultaneous excitation by stellar radiation and shocks

We use a two-phase model for the structure of the circumstellar nebulae of hot stars to analyze the radiative cooling of a dense, compact cloud behind the shock produced by the compression of the cloud by hot gas from the stellar wind, taking into account ionization and heating by radiation from the central star. We can distinguish three stages of the evolution of the cloud during its compression. In the first stage, relevant for the entire cloud before compression and the gas ahead of the shock front, the state of the gas is determined purely by ionization by the stellar radiation. The next stage is characterized by the simultaneous action of two gas excitation mechanisms—photoionization by the stellar radiation and shock heating. In this stage, the gas intensively radiates thermal energy received at the shock front. After radiative cooling, in the final stage, ionization and heating of the gas are again determined mainly by the star. To compute the spectrum of the cloud radiation, we solved for the propagation of a plane-parallel, homogeneous flux through the shock front in the radiation field of the hot star. The computations show that a combination of two excitation mechanisms considerably enriches the theoretical spectrum. The relative intensities of emission lines of a single cloud may resemble either those for an HII region or of a supernova remnant.     

Spots and activity cycles of HD 199178

We present our analysis of photometry for the FK Com star HD 199178 (V1794 Cyg). The V-band light curves are used to restore the distribution of temperature inhomogeneities on the stellar surface. The spots on the surface of HD 199178 are concentrated at two preferred longitudes separated by 0.5 in phase (180° in longitude). In addition to the quasi-periodic switching of the most active area between these two longitudes, which occurs in cycles of 2.1–2.4 or 4.1 years, we suspect that the two active areas moved toward each other across the stellar surface, possibly merging into a single formation. The detected cycle in the star’s brightness variations of about 8.0 years is also clearly visible in variations of the star’s spottedness. Themagnetic activity of the FK Com star HD199178 is in many ways similar to that observed for the prototype of this group.     

Massive close binary stars and gamma-ray bursts

We analyze the observed parameters of massive extremely close binaries containing Wolf-Rayet stars and black holes, and identify those systems whose supernova outbursts lead to the formation of rapidly rotating Kerr black holes. It is proposed that the formation of such a black hole is accompanied by a strong gamma-ray burst. Several types of observed systems satisfy the conditions necessary for the formation of a Kerr black hole: BH+WR, BH+OB, WR+O, and BH+K,M.     

Evolution of Wolf-Rayet stars in binary systems: An analysis of the mass and orbital-eccentricity distributions

We have undertaken a statistical study of the component mass ratios and the orbital eccentricities of WR + O close binary, detached main-sequence (DMS), contact early-type (CE), and semidetached (SD) systems. A comparison of the characteristics of WR + O systems and of DMS, CE, and SD systems has enabled us to draw certain conclusions about the evolutionary paths of WR + O binaries and to demonstrate that up to 90% of all known WR + O binaries formed as a result of mass transfer in massive close O + O binary systems. Since there is a clear correlation between the component masses in SD systems with subgiants, the absence of an anticorrelation between the masses of the WR stars and O stars in WR + O binaries cannot be considered evidence against the formation of WR + O binaries via mass transfer. The spectroscopic transitional orbital period P _tr ^sp corresponding to the transition from nearly circular orbits (e _sp<0.1) to elliptical orbits (e _sp≥0.1) is ～14^d for WR + O systems and ～2^d–3^d for OB + OB systems. The period range in which all WR + O orbits are circular \((1\mathop d\limits_. 6 \leqslant P \leqslant 14^d )\) is close to the range for SD systems with subgiants, \(0\mathop d\limits_. 7 \leqslant P \leqslant 15^d \). The large difference between the P _tr ^sp values for WR + O and OB + OB systems suggests that a mechanism of orbit circularization additional to that for OB + OB systems at the DMS stage (tidal dissipation of the orbital energy due to radiative damping of the dynamical tides) acts in WR + O binaries. It is natural to suggest mass transfer in the parent O + O binaries as this supplementary orbit-circularization mechanism. Since the transitional period between circular and elliptical orbits for close binaries with convective envelopes and ages of 5×10⁹ years is \(P_{tr} = 12\mathop d\limits_. 4\), the orbits of most known SD systems with subgiants had enough time to circularize during the DMS stage, prior to the mass transfer. Thus, for most SD systems, mass transfer plays a secondary role in circularization of their orbits.In many cases, the initial orbital eccentricities of the O + O binary progenitors of WR + O systems are preserved, due to the low viscosity of the O-star envelopes and the short timescale for their nuclear evolution until the primary O star fills its Roche lobe and the mass transfer begins. The mass transfer in the parent O + O systems is short-lived, and the number of orbital cycles during the early mass-transfer stage is relatively low (lower than for the progenitors of SD systems by three or four orders of magnitude). The continued transfer of mass from the less massive to the more massive star after the component masses have become equal leads to the formation of a WR + O system, and the orbit's residual eccentricity increases to the observed value. The increase of the orbital eccentricity is also facilitated by variable radial mass loss via the wind from the WR star in the WR + O system during its motion in the elliptical orbit. The result is that WR + O binaries can have considerable orbital eccentricities, despite their intense mass transfer. For this reason, the presence of appreciable eccentricities among WR + O binaries with large orbital periods cannot be considered firm evidence against mass transfer in the parent O + O binary systems. Only for the WR + O binaries with the longest orbital periods (4 of 35 known systems, or 11 %) can the evolution of the parent O + O binaries occur without filling of the Roche lobe by the primary O star, being governed by radial outflow in the form of the stellar wind and possibly by the LBV phenomenon, as in the case of HD 5980.     

A pre-outburst of the classical symbiotic star BF Cyg in 1987–1989

We analyze previously unpublished spectrophotometric data for the classical symbiotic star BF Cyg obtained in 1982–1989. These data show that BF Cyg already started a new period of activity in 1987–1989, earlier believed to correspond to quiescence. An optical spectrum taken on September 29, 1987 displays several features characteristic of activity in classical symbiotic stars. During the same time interval, excess near-UV radiation was detected, possibly from the accretion disk, and evidence for Rayleigh scattering of the hot component’s radiation on atomic hydrogen was detected in the far-UV, at phases previously believed not to be characteristic of this effect. This latter finding may testify to an enhanced (compared to 1979–1986) rate of mass loss from the cool component of BF Cyg. We discuss the interrelation between an increased mass-loss rate of the cool component and the hot component’s activity in the context of classical symbiotic stars. Original Russian Text ￠ A.A. Tatarnikova, E.A. Kolotilov, A.M. Tatarnikov, 2008, published in Astronomicheskiĭ Zhurnal, 2008, Vol. 85, No. 2, pp. 151–162.     

Disk wind in the radiation of young intermediate-mass stars

We consider the role of the disk wind in the formation of hydrogen emission spectra of young intermediate-mass Herbig Ae stars. The assumed parameters for the kinematic characteristics of the wind region are typical of the theory of magnetocentrifugal disk winds. Our computations of the excitation and ionization states for atoms in the emitting region are based on the Sobolev approximation for a medium with a large velocity gradient. The calculations show that hydrogen-line emission profiles can have a complex structure that depends on the disk-wind parameters, as well as on the system’s orientation relative to the line of sight. The model proposed is able to explain most types of Hα profiles observed for Herbig Ae stars.     

Thermodynamic properties and equation of state of fcc aluminum and bcc iron, derived from a lattice vibrational method

We use a lattice vibrational technique to derive thermophysical and thermochemical properties of the pure elements aluminum and iron in pressure–temperature space. This semi-empirical technique is based on either the Mie–Grüneisen–Debye (MGD) approach or an extension of Kieffer’s model to incorporate details of the phonon spectrum. It includes treatment of intrinsic anharmonicity, electronic effects based on the free electron gas model, and magnetic effects based on the Calphad approach. We show that Keane’s equation of state for the static lattice is better suitable to represent thermodynamic data for aluminum from 1 bar to pressures in the multi-megabar region relative to Vinet’s universal and the Birch–Murnaghan equation of state. It appears that the MGD and Mie–Grüneisen–Kieffer approach produce similar results, but that the last one better represents heat capacity below room temperature. For iron we show that the high temperature behavior of thermal expansivity can be explained within the Calphad approach by a pressure-dependent Curie temperature with a slope between –1 and 0 K/GPa.     

The role of metals in cooling gas behind shock fronts in the atmospheres of cool stars

The influence of various chemical elements on radiative cooling of the gas flowing from a viscous jump is investigated in a model with a stationary shock in the atmosphere of a cool star. A closed system of equations is written for the thermal energy per heavy particle, the electron temperature, and the relative concentrations of elements in all ionization states. In addition to hydrogen and helium, atomic, singly ionized, and doubly ionized carbon, nitrogen, oxygen, sodium, magnesium, aluminum, silicon, sulfur, potassium, calcium, and iron are included, assuming they have their normal cosmic abundances. The high optical depth in Lyman-series lines leads to a return of the thermal energy to electrons via secondary collisions. As a result, the contribution of hydrogen to the cooling rate falls to the level of the contribution of metals, mainly carbon, magnesium, and iron. Thus, such shock models are able to explain the presence of bright metal lines in the spectra of cool and solar-type stars.     

The origin of intergalactic stars in galaxy clusters

The paper analyzes possible origins of stars located in intergalactic space that are not bound to specific galaxies, which comprise 15–50% of all stars in galaxy clusters. Some such stars can form in streams of intergalactic gas flowing around gas-rich disk galaxies moving in the cluster. Others may be the products of the decay of young, low-mass, spheroidal galaxies after the loss of their gaseous components during an initial burst of star formation. The decay of low-mass disk galaxies moving at high speeds after they have lost their gaseous components due to the pressure of the incident flow of dense intergalactic gas is possible in the cluster core. The largest fraction of intergalactic stars are probably produced by the partial disruption of galaxies as a result of close passages, collisions, or mergers. Collisions of low-mass, gas-rich galaxies are especially good suppliers of intergalactic stars. Both stars from decaying stellar components of galaxies and stars arising in the gaseous components of colliding galaxies can be supplied to the intergalactic medium. The merger of galaxies harboring supermassive black holes in their nuclei could lead to the partial or total disruption of these galaxies during the deceleration of the binary black hole that is formed during the merger. An enhanced density of intergalactic stars is observed in the cores of galaxy clusters, underscoring the role of galaxy collisions in the formation of the intergalactic stellar population, since the frequency of galaxy collisions grows with their density.     

Instabilities in the ionization zones around the first stars

We consider the evolution of the ionization zone around Population III stars with M _* ?? 25?C200M _?? in protogalaxies with M ?? 10⁷ M _?? at redshifts z = 12, assuming that the dark-energy profile is a modified isothermal sphere. We study the conditions for the growth of instabilities in the ionization zones. The Rayleigh-Taylor and thermal instabilities develop efficiently in the ionization zones around 25?C40M _?? stars, while this efficiency is lower for stars withM _* ?? 120M _??. For more massive stars (??200M _??), the flux of ionizing photons is strong enough to considerably reduce the gas density in the ionization zone, and the typical lifetimes of stars (??2 Myr) are insufficient for the growth of instabilities. The gas in a protogalaxy with M ?? 10⁷ M _?? with a 200M _?? central star is completely ionized by the end of the star??s lifetime; in the case of a 120M _?? central star, only one-third of the total mass of gas is ionized. Thus, ionizing photons from stars with M _* ? 120M _?? cannot leave protogalaxies with M ? 10⁷ M _??. If the masses of the central stars are 25 and 40M _??, the gas in protogalaxies of this mass remains essentially neutral. We discuss the consequences of the evolution of the ionization zones for the propagation of the envelope after the supernova explosions of the strs and the efficiency of enrichment of the intergalactic medium in heavy elements.     

Dust formation in binaries with OB and WR components in a two-phase stellar-wind model

We describe the formation of carbon dust in binary systems with hot components as a result of the collisions of clouds in a two-phase stellar-wind model. Calculations are made for the well studied system WR 140. The collisions lead to the formation of composite clouds and shock waves, with the temperature at the shock front equal to about 3×10⁸ K along both sides of the interface boundary. During isobaric deexcitation to (0.5–0.7)×10⁴ K, the cloud density increases by a factor of several thousand; its thickness in the direction of the shock decreases by the same factor. After deexcitation, the hydrogen inside the composite cloud is in its atomic state, while the carbon remains ionized. The deexcitation is followed by expansion of the cloud, which moves away from both stars. During the first 10⁶ s, its thickness remains relatively small, so that the expansion is one-dimensional. The radiation field inside the cloud decays, resulting in the recombination of the carbon. Further expansion of the cloud leads to adiabatic cooling, and the formation of dust particles becomes possible. After the dimensions of the cloud have become roughly the same in all directions, its expansion is isotropic, so that it becomes transparent within approximately 10⁶ s, and the dust is heated to (1.0–1.4)×10³ K, observed as an IR “lare.” The time required for the cloud to move from the exciting star and heat the dust is comparable to the observed delay in the increased IR emission relative to the time of periastron.     

Periodic variations in the emission spectrum of T Tauri

We have analyzed a many-year homogeneous series of spectroscopic data for the star T Tau N. Our optical observations of hydrogen lines and the Hand K lines of CaII have revealed periodic variations of the emission spectrum with a period of 33 ± 1.5 days. The detected period is also observed in IUE spectrograms in the ultraviolet. Variations of the polarization of the star’s light also agree with the derived period. Our results suggest the possible presence of a comet-like extended object orbiting T Tau N.     

Line-profile microvariability in OB-star spectra: the Supergiant λ Cep (O6If(n))

We observed the bright O6If(n) supergiant λ Cep in 1997 with the 6-m optical telescope of the Special Astrophysical Observatory and in 2007 with the 1.8-m telescope of the Bohyunsan Optical Astronomy Observatory (South Korea). A total of 90 spectra of the star were acquired, with good time resolution (10 minutes), signal-to-noise ratios 150–300, and spectral resolutions of 45 000–60 000. We detected line-profile variations of H, HeI, and HeII lines. It is suggested that the detected variations are due to non-radial photospheric pulsations and the star’s rotation (rotational profile modulation).     

Modeling of Spectroscopic and Interferometric Observations of the Herbig Star VV Ser with Hybrid Models

Modeling of hydrogen emission lines is a powerful tool to study physical processes in the nearest vicinity of young stars because spectral lines carry information on the kinematics and physical conditions of the gas. One of the lines that probe emitting regions closest to the star is the Br$$\gamma $$ line. We consider different types of hybrid models to reproduce both interferometric VLTI-AMBER observations and LBT-LUCIFER spectroscopic observations of the single-peak profile of the Br$$\gamma $$ line of the Herbig AeBe star (HAEBE) VV Ser, a member of the UX Ori type subclass. We computed models of a magneto-centrifugal disk wind, a magnetospheric accretion region (magnetosphere), Cranmer’s polar wind, and scattered light from circumstellar polar dust. Furthermore, we calculated hybrid two-component models consisting of a disk wind and one of the aforementioned models. We computed visibilities and line profiles for all types of models and compared them with the available interferometric observations to constrain model parameters. We conclude that for the inclinations reported for this star (60°–70°), the disk wind alone cannot explain the Br$$\gamma $$ line profile although it may be a dominant contributor to the hydrogen line radiation. However, magneto-centrifugal disk wind in combination with aforementioned emitting regions (magnetosphere, polar wind, or scattered light from polar dust) may be able to reproduce the observations.     

New variable stars on digitized Moscow collection plates. The field of 66 Ophiuchi

In the course of a program to digitize the astronomical plates of the Sternberg Astronomical Institute’s plate stacks, we are developing algorithms for searching for new variable stars and studying them using digitized photographic plates. We have discovered and studied 480 new variable stars in a 10° × 10° field of view centered on 66 Ophiuchi. The digitized plate negatives used are from the 40-cm astrograph, and are 30 × 30 cm in size. These stars include three new Cepheids of the Galaxy’s spherical component, 157 eclipsing binaries, 11 high-amplitude δ Scuti stars (HADSs), 144 RR Lyrae stars, 110 irregular variables (109 LB and one white star), and 55 semi-regular red variables. New important information has been obtained for 43 known variables, which we have classified and derived or improved their light elements; an erroneous identification of the Mira V404 Oph has been corrected. We have also identified more than 50 suspect brightness variables; a program of CCD observations of these suspected variables has been initiated. Our discoveries of new variable stars were performed in a star field with a large number of known variables, detected earlier photographically or using CCD techniques. The discovery of hundreds of new variables in a well-studied region of sky demonstrates that archive photographs possess a large information potential that has remained unrealized.     

Circumstellar envelopes of hot stars: Formation of extended regions of weakly ionized gas

We propose a model explaining the presence of vast regions of partially ionized gas in the interstellar medium. The circumstellar envelope of a hot star absorbs soft ionizing radiation, but transmits an appreciable fraction of the hard photons, which are absorbed much more weakly than photons with energies close to the ionization limit. For this reason, the radiation attenuated by the envelope becomes harder, and can penetrate to larger distances. For stars hotter than 50 000 K, the transition zone between the ionized and neutral gas can extend to tens or hundreds of parsecs. Thus, a region of partially ionized hydrogen, with a small gradient of the degree of ionization without a well-defined inner HII zone, can form in the interstellar medium.     

The role of external factors in the evolution of galaxies

We consider the evolution of galaxies in dense galactic clusters. Observations and theoretical estimates indicate that this evolution may be specified to a large extent by collisions between galaxies, as well as interactions between the gaseous components of disk galaxies and intergalactic gas. We analyze collisions between disk galaxies with gaseous components using a simple model based on a comparison of the duration of a collision and the characteristic cooling time for the gas heated by the collision, and also of the relative masses of stars and gas in the colliding disk galaxies. This model is used to analyze scenarios for collisions between disk galaxies with various masses as a function of their relative velocities. Our analysis indicates that galaxies can merge, lose one or both of their gaseous components, or totally disintegrate as a result of a collision; ultimately, a new galaxy may form from the gas lost by the colliding galaxies. Disk galaxies with mass M _G and velocities exceeding ～300 (M _G/10¹⁰ M _⊙)^1/2 km/s in intergalactic gas in clusters with densities ～10^?27 g/cm³ can lose their gas due to the pressure of inflowing intergalactic gas, thereby developing into E(SO) galaxies.     

The nature hypervelocity stars

We list and analyze the main currently known mechanisms for accelerating the space motions of stars. A high space velocity of a star can be a consequence of its formation in the early stages of the evolution of a massive galaxy, when it was spheroidal and non-stationary, so that stars were born with velocities close to the escape velocity for the galaxy. Another possibility is that the star arrived from another galaxy with a velocity that is high for our Galaxy. The decay of unstable close multiple stars or supernova explosions in close binaries can also provide velocities of up to several hundreds of km/s to main-sequence stars and velocities of up to ∼1000 km/s to degenerate stars, neutron stars, and stellar-mass black holes. The merger of components of a binary system containing two neutron stars or a neutron star and a black hole due to gravitational-wave radiation can accelerate the nascent black hole to a velocity∼1000 km/s. Hypervelocity relativistic stars can be born due to asymmetric neutrino ejection during a supernova explosion. Stars can be efficiently accelerated by single and binary supermassive black holes (with masses from several millions to several billions of solar masses) in the nuclei of galaxies. Thanks to their gravitational field and fast orbital motion (in the case of binary objects), supermassive black holes are able to accelerate even main-sequence stars to relativistic velocities.