The evolutionary status of ultraluminous X-ray sources

We analyze models for quasi-stationary, ultraluminous X-ray sources (ULXs) with luminosities 10³⁸–10⁴⁰ erg/s exceeding the Eddington limit for a ～1.4M_⊙ neutron star. With the exception of relatively rare stationary ULXs that are associated with supernova remnants or background quasars, most ULXs are close binary systems containing a massive stellar black hole (BH) that accretes matter donated by a stellar companion. To explain the observed luminosities of ～10⁴⁰ erg/s, the mass of the BH must be ～40M_⊙ if the accreted matter is helium and ～60M_⊙ if the accreted matter has the solar chemical composition. We consider donors in the form of main-sequence stars, red giants, red supergiants, degenerate helium dwarfs, heavy disks that are the remnants of disrupted degenerate dwarfs, helium nondegenerate stars, and Wolf-Rayet stars. The most common ULXs in galaxies with active star formation are BHs with Roche-lobe-filling main-sequence companions with masses ～7M_⊙ or close Wolf-Rayet companions, which support the required mass-exchange rate via their strong stellar winds. The most probable candidate ULXs in old galaxies are BHs surrounded by massive disks and close binaries containing a BH and degenerate helium-dwarf, red-giant, or red-supergiant donor.     

Wolf-Rayet stars and relativistic objects: Distinctions between the mass distributions in close binary systems

The observed properties of Wolf-Rayet stars and relativistic objects in close binary systems are analyzed. The final masses M _CO ^f for the carbon-oxygen cores of WR stars in WR + O binaries are calculated taking into account the radial loss of matter via stellar wind, which depends on the mass of the star. The analysis includes new data on the clumpy structure of WR winds, which appreciably decreases the required mass-loss rates $\dot M_{WR}$ for the WR stars. The masses M _CO ^f lie in the range (1–2)M _⊙–(20–44)M _⊙ and have a continuous distribution. The masses of the relativistic objects M _x are 1–20M _⊙ and have a bimodal distribution: the mean masses for neutron stars and black holes are 1.35 ± 0.15M _⊙ and 8–10M _⊙, respectively, with a gap from 2–4M _⊙ in which no neutron stars or black holes are observed in close binaries. The mean final CO-core mass is $\overline M _{CO}^f = 7.4 - 10.3M_ \odot$ , close to the mean mass for the black holes. This suggests that it is not only the mass of the progenitor that determines the nature of the relativistic object, but other parameters as well-rotation, magnetic field, etc. One SB1R Wolf-Rayet binary and 11 suspected WR + C binaries that may have low-mass companions (main-sequence or subgiant M-A stars) are identified; these could be the progenitors of low-mass X-ray binaries with neutron stars and black holes.     

Evolution of close binaries and gamma-ray bursts

We analyze the late stages of evolution of massive (M ₀ ? 8 M _⊙) close binaries, from the point of view of possible mechanisms for the generation of gamma-ray bursts. It is assumed that a gamma-ray burst requires the formation of a massive (～1 M _⊙), compact (R ? 10 km) accretion disk around a Kerr black hole or neutron star. Such Kerr black holes are produced by core collapses of Wolf-Rayet stars in very close binaries, as well as by mergers of neutron stars and black holes or two neutron stars in binaries. The required accretion disks can also form around neutron stars that were formed via the collapse of ONeMg white dwarfs. We estimate the Galactic rate of events resulting in the formation of rapidly rotating relativistic objects. The computations were carried out using the “Scenario Machine.”     

Nonlinear radial pulsations of Wolf-Rayet stars

The evolution of Population I stars with initial masses 60 M _⊙ ≤ M _ZAMS ≤ 120 M _⊙ is computed up to the Wolf-Rayet stage, when the central helium abundance decreases to Y _c ≈ 0.05. Several models from evolutionary sequences in the core helium-burning stage were used as initial conditions when solving the equations of radiative hydrodynamics for self-exciting stellar radial pulsations. The low-density envelope surrounding the compact core during the core helium burning is unstable against radial oscillations in a wide range of effective temperatures extending to T _eff ～ 10⁵ K. The e-folding time of the amplitude growth is comparable to the dynamical time scale of the star, and, when the instability ceases growing, the radial displacement of the outer layers is comparable to the stellar radius. Evolutionary changes of the stellar radius and luminosity are accompanied by a decrease in the amplitude of radial pulsations, but, at the effective temperature T _eff ≈ 10⁵ K, the stellar oscillations are still nonlinear, with a maximum expansion velocity of the outer layers of about one-third the local escape velocity. The period of the radial oscillations decreases from 9 hr to 4 min as stellar mass decreases from M = 28 M _⊙ to M = 6 M _⊙ in the course of evolution. The nonlinear oscillations lead to a substantial increase of the radii of the Lagrangian mass zones compared to their equilibrium radii throughout the instability region. The instability of Wolf-Rayet stars against radial oscillations is due to the action of the κ mechanism in the iron-group ionization zone, which has a temperature of T ～ 2 × 10⁵ K.     

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.     

Gamma-ray bursts—tracers of the history of star formation in the Universe

The rate of gamma-ray bursts (GRBs) in the Galaxy is estimated assuming that these events result from the formation of rapidly rotating Kerr black holes during the core collapse of massive, helium, Wolf-Rayet secondary components in very close binary systems. This process brings about rapid rotation of the cores of such Wolf-Rayet stars, inevitably resulting in the formation of Kerr black holes during type Ib,c supernovae. The current rate of formation of Kerr black holes (GRBs) in the Galaxy is about 3×10^?5/year. Collimation of the gamma-ray radiation into a small solid angle (about 0.1–0.01 sr) brings this rate into consistency with the observed rate of GRBs, estimated to be 10^?6–10^?7/year. Possible immediate progenitors of GRBs are massive X-ray binaries with X-ray luminosities of 10³⁸–10⁴⁰ erg/s. Due to the short lifetimes of the progenitors and the very high brightnesses of GRBs, the GRB rate can provide information about the history of star formation in the Universe on the Hubble time scale. A model in which the star-formation rate is determined by the conditions for ionization of the interstellar gas, whose density and volume are determined by supernovae, yields a Galactic star-formation history that can be viewed as representing the history of star formation in the Universe. The theoretical history of star formation is in satisfactory agreement with the history reconstructed from observations. The theoretical model for the history of star formation in the Galaxy can also be used to assess the influence of dust on optical observations of supernovae and GRBs in galaxies of various ages.     

Radial pulsations of Wolf-Rayet stars and the problem of mass loss in them

In order to study the problem of mass loss by Wolf-Rayet stars, we carried out numerical simulations of non-linear pulsations of these objects. Although our computations do not show direct dynamical mass loss, qualitative estimates show that certain as yet unobserved pulsations could in principle, create conditions facilitating mass outflows from Wolf-Rayet stars.     

Modeling of rapid variability in the spectral line profiles of Wolf-Rayet stars

Profiles of variable emission lines in the spectra of Wolf-Rayet stars are calculated using a stochastic cloud model for the inhomogeneous atmospheres of early-type stars. The model assumes that most of the line flux is formed in cold, dense condensations (clouds) that move through a rarified inter-cloud medium whose density monotonically decreases outwards. The formation of clouds is taken to be stochastic. Wavelet analysis is used to estimate the parameters of cloud ensembles. The model can reproduce the general pattern of line-profile variability observed in the spectra of Wolf-Rayet stars.     

The formation and evolution of the most massive stars

We consider the formation of massive stars under the assumption that a young star accretes material from the protostellar cloud through its accretion disk while losing gas in the polar directions via its stellar wind. The mass of the star reaches its maximum when the intensity of the gradually strengthening stellar wind of the young star becomes equal to the accretion rate. We show that the maximum mass of the forming stars increases with the temperature of gas in the protostellar cloud T ₀, since the rate at which the protostellar matter is accreted increases with T ₀. Numerical modeling indicates that the maximum mass of the forming stars increases to ～900 M _⊙ for T ₀ ～ 300 K. Such high temperatures of the protostellar gas can be reached either in dense star-formation regions or in the vicinity of bright active galactic nuclei. It is also shown that, the lower the abundance of heavy elements in the initial stellar material Z, the larger the maximum mass of the star, since the mass-loss rate due to the stellar wind decreases with decreasing Z. This suggests that supermassive stars with masses up to 10⁶ M _⊙ could be formed at early stages in the evolution of the Universe, in young galaxies that are almost devoid of heavy elements. Under the current conditions, for T ₀ = (30–100) K, the maximum mass of a star can reach ～100M _⊙, as is confirmed by observations. Another opportunity for the most massive stars to increase their masses emerges in connection with the formation and early stages of evolution of the most massive close binary systems: the most massive stars can be produced either by coalescence of the binary components or via mass transfer in such systems.     

The Spots and Activity of Stars in the Beehive Cluster Observed by the Kepler Space Telescope (K2)

The spottedness parameters S (the fraction of the visible surface of the star occupied by spots) characterizing the activity of 674 stars in the Beehive Cluster (age 650 Myr) are estimated, together with variations of this parameter as a function of the rotation period, Rossby number Ro and other characteristics of the stars. The activity of the stars in this cluster is lower than the activity of stars in the younger Pleiades (125 Myr). The average S value for the Beehive Cluster stars is 0.014, while Pleiades stars have the much higher average value 0.052. The activity parameters of 61 solar-type stars in the Beehive Cluster, similar Hyades stars (of about the same age), and stars in the younger Pleiades are compared. The average S value of such objects in the Beehive Cluster is 0.014± 0.008, nearly coincident with the estimate obtained for solar-type Hyades stars. The rotation periods of these objects are 9.1 ± 3.4 day, on average, in agreement with the average rotation period of the Hyades stars (8.6 ^d ). Stars with periods exceeding 3–4 ^d are more numerous in the Beehive Cluster than in the Pleiades, and their periods have a larger range, 3–30 ^d . The characteristic dependence with a kink at Ro (saturation) = 0.13 is not observed in the S–Rossby number diagram for the Beehive and Hyades stars, only a clump of objects with Rossby numbers Ro > 0.7. The spottedness data for the Beehive Cluster and Hyades stars are in good agreement with the S values for dwarfs with ages of 600–700 Myr. This provides evidence for the reliability of the results of gyrochronological calibrations. The data for the Beehive and Pleiades stars are used to analyze variations in the spot-forming activity for a large number of stars of the same age that are members of a single cluster. A joint consideration of the data for two clusters can be used to draw conclusions about the time evolution of the activity of stars of different masses (over a time interval of the order of 500 Myr).     

Stellar sources of the short-lived radionuclides in the early solar system

We discuss the possible stellar sources of short-lived radionuclides (SLRs) known to have been present in the early solar system (²⁶Al, ³⁶Cl, ⁴¹Ca, ⁵³Mn, ⁶⁰Fe, ¹⁰⁷Pd, ¹²⁹I, ¹⁸²Hf, ²⁴⁴Pu). SLRs produced primarily by irradiation (⁷Be, ¹⁰Be) are not discussed in this paper. We evaluate the role of the galactic background in explaining the inventory of SLRs in the early solar system. We review the nucleosynthetic processes that produce the different SLRs and place the processes in the context of stellar evolution of stars from 1 to 120 M_⊙. The ejection of newly synthesized SLRs from these stars is also discussed. We then examine the extent to which each stellar source can, by itself, explain the relative abundances of the different SLRs in the early solar system, and the probability that each source would have been in the right place at the right time to provide the SLRs. We conclude that intermediate-mass AGB stars and massive stars in the range from ∼20 to ∼60 M_⊙ are the most plausible sources. Low-mass AGB stars fail to produce enough ⁶⁰Fe. Core-collapse Type II supernovae from stars with initial masses of <20 M_⊙ produce too much ⁶⁰Fe and ⁵³Mn. Sources such as novae, Type Ia supernovae, and core-collapse supernovae of O-Ne-Mg white dwarfs do not appear to provide the SLRs in the correct proportions. However, intermediate-mass AGB stars cannot provide ⁵³Mn or the r-process elements, so if an AGB star provided the ⁴¹Ca, ³⁶Cl, ²⁶Al, ⁶⁰Fe, and ¹⁰⁷Pd, and if a late stellar source is required for ⁵³Mn and the r-process elements, then two types of sources would be required. A separate discussion of the production of r-process elements highlights the difficulties in modeling their production. There appear to be two sources of r-process elements, one that produces the heavy r-process elements, including the actinides, and one that produces the elements from N to Ge and the elements ∼110 < A < ∼130. These can be assigned to SNII explosions of stars of ?11 M_⊙ and stars of 12-25 M_⊙, respectively. More-massive stars, which leave black holes as supernova remnants, apparently do not produce r-process elements.     

Sources of efficient acceleration of solar flare particles: Observational aspects

This work continues earlier statistical analyses of catalogued proton events in cosmic rays. The spectra of proton enhancements identified from the logarithm of the ratio of the proton fluxes with energies exceeding 100 and 10 MeV, δ = log(F ₁₀₀/F ₁₀), are studied focusing on 172 powerful events with favorable conditions for escaping from the corona and subsequent propagation in the interplanetary space. The δ distribution for the flares is Gaussian, with a comparatively weak spread in δ. The distribution maximum corresponds to an excess of the 10 MeV flux over the 100 MeV flux by a factor of 30. The fact such a frequent spectrum does exist supports the idea that both soft and hard protons are efficiently accelerated, probably by a single mechanism that operates during the explosive phases of the flares. The sizes of the loops of M2-X4 flares observed by the Yohkoh Hard X-Ray Telescope at energies exceeding 50 keV indicate low heights for the main acceleration regions. There is some excess over the Gaussian distribution for “softer” events. Some post-eruptive phenomena can be clearly distinguished in these events, and they display a correlation between δ and the total flare duration. Thus, there are two sources of the particle acceleration operating at low and high altitudes during the explosive and post-eruptive phases of the solar flares, respectively. The second source, which is manifested through some prolonged flares and filament-ejection phenomena, accelerates protons only to energies of 10–30 MeV.     

Bismuth and 208Pb microdistributions in enstatite chondrites

Polished sections of 5 enstatite chondrites have been irradiated with 30 MeV ⁴He ions to produce the alpha-radioactive nuclei ²¹¹At and ²¹⁰Po from ²⁰⁹Bi and ²⁰⁸Pb, respectively. The distribution of alpha activity can be mapped, using cellulose nitrate as an alpha track detector, to give the corresponding Bi or Pb distributions in the meteorite. No strong localization of Bi or ²⁰⁸Pb was found; relatively uniform track distributions were observed. In particular, metal or sulfide grains are not enriched in Bi or Pb (relative to bulk), which is in agreement with the predictions of nebular condensation calculations. While the track distributions appear uniform, the results of detailed, track-by-track mappings of the Bi detectors indicate that the Bi is not totally randomly distributed; the statistical fluctuations in the observed track density are different for the cases where the Bi is totally randomly distributed and where the Bi is localized in point sources. Assuming that the Bi in a given sample is localized in identical point sources which are uniformly distributed throughout the sample, the observed relative population densities of clusters (‘stars’) of small numbers of tracks (2–5) corresponds to Bi being localized, with ~90% in grains with about 10^?16g-Bi (~3 × 10⁵Biatoms), and with ~10% in 4 × 10^?14 g-Bi sources. If these are elemental Bi, as predicted theoretically, they are ~ 10² Å and 10³ Å in size, respectively.     

The Wolf–Rayet star HD 192163 as a possible evolutionary progenitor of a low-mass X-ray binary

We have used 46 high-resolution echelle spectra of the Wolf-Rayet star HD 192163 taken in 2005–2009 at the Cassegrain focus of the 2-m Zeiss-2000 telescope of the Shamakha Astrophysical Observatory to study profiles of the five strongest emission lines (HeII 4859, HeII 5411, CIV 5808, HeI 5875, (HeII + Hα) 6560). We also obtained four echelle spectrograms of the Wolf-Rayet star HD 191765 for a comparative study of the NaI 5890 (D₂) and NaI 5896 (D₁) interstellar absorption lines. The echelle spectrograms were reduced using the DECH20 code. We determined the equivalent widths, radial velocities, central intensities, and half-widths of the emission lines. We detected variations in the violetwing of the (HeII + Hα) 6560 emission band (between λ ～ 6496 Å and λ ～ 6532 Å). Our statistical analysis of the radial velocities available for the (HeII + Hα) 6560 emission band reveals a peak at the ～1% level at the frequency ν = 0.195 d^?1, corresponding to the period P = 5.128^d. We also studied the NaI 5890 (D₂) and NaI 5896 (D₁) interstellar absorption lines, which are important for understanding the nature of the nebula NGC 6888, whose origin is related to HD 192163. Asymmetric profiles were found for the NaI 5890 and NaI 5896 interstellar absorption lines, interpreted as reflecting a contribution from NGC 6888 to these lines. We suggest that the detected profile variations of the (HeII + Hα) 6560 emission band and the periodic variations of its Doppler shifts (P = 5.128^d) are due to the existence of a low-mass K-M star companion to the Wolf-Rayet star. HD 192163 is a possible evolutionary progenitor of a low-mass X-ray binary.     

Evolution of the components of the close binary WR 20a

We have analyzed the evolution of the components of the unique massive binary system WR 20a, which consists of a Wolf-Rayet nitrogen star and an Of star with an extremely small separation. The estimated masses of the components are 83 and 82 M _⊙, which are among the highest stellar mass inferred. We have carried out numerical modeling of the evolution of the components, taking into account the mass loss due to the stellar wind inherent to massive stars. In a scenario in which the systemis detached from the time the components reach the main sequence until its present state, the initial component masses are inferred to be close to 110 M _⊙, if the initial masses of the stars were equal, or 120 and 100 M _⊙, if they were different. Currently, the components are evolved main-sequence stars, whose surfaces are relatively little enriched by helium. The further evolution of the system will result in one of the components filling its Roche lobe and evolution within a common envelope. As a result, the components may coalesce, leading to the formation of a single massive black hole the supernova explosion. Otherwise, depending on the masses of the resulting black holes, either a binary system with two black holes or two free black holes will be formed. In the latter case, gamma-ray bursts will be observed.     

A model for the population of helium stars in the Galaxy: Low-mass stars

We model the Galactic ensemble of helium stars using population synthesis techniques, assuming that all helium stars are formed in binaries. In this picture, single helium stars are produced by mergers of helium remnants of the components of close binaries (mainly, the merging of helium white dwarfs) or in the disruption of binaries with helium components during supernova explosions. The estimated total birthrate of helium stars in the Galaxy is 0.043 yr^?1; the total number is 4 × 10⁶; and the binarity rate is 76%. We construct a subsample of low-mass (M_He ? 2M_⊙) helium stars defined by observational selection effects: the limiting magnitude (V_He ≤ 16), ratio of the magnitudes of the components in binaries (V_He ≤ V_comp), and lower limit for the semiamplitude of the radial velocity required for detecting binarity (K_min = 30 km s^?1). The parameters of this subsample are in satisfactory agreement with observations of helium subdwarfs. In particular, the binarity rate in the selection-limited sample is 58%. We analyze the relations between the orbital periods and masses of helium subdwarfs and their companions in systems with various combinations of components. We predict that the overwhelming majority (～97%) of unobserved companions to helium stars will be white dwarfs, predominantly, carbon-oxygen white dwarfs.     

Spots and the Activity of Stars in the Hyades Cluster from Observations with the Kepler Space Telescope (K2)

Observations of the K2 mission (continuing the program of the Kepler Space Telescope) are used to estimate the spot coverage S (the fractional area of spots on the surface of an active star) for stars of the Hyades cluster. The analysis is based on data on the photometric variations of 47 confirmed single cluster members, together with their atmospheric parameters, masses, and rotation periods. The resulting values of S for these Hyades objects are lower than those stars of the Pleiades cluster (on average, by ΔS ~ 0.05?0.06). A comparison of the results of studies of cool, low-mass dwarfs in the Hyades and Pleiades clusters, as well as the results of a study of 1570 M stars from the main field observed in the Kepler SpaceMission, indicates that the Hyades stars are more evolved than the Pleiades stars, and demonstrate lower activity. The activity of seven solar-type Hyades stars (S = 0.013 ± 0.006) almost approaches the activity level of the present-day Sun, and is lower than the activity of solar-mass stars in the Pleiades (S = 0.031 ± 0.003). Solar-type stars in the Hyades rotate faster than the Sun (〈P〉 = 8.6 ^d ), but slower than similar Pleiades stars.     

The OSACA database and a kinematic analysis of stars in the solar neighborhood

We transformed radial velocities compiled from more than 1400 published sources, including the Geneva-Copenhagen survey of the solar neighborhood (CORAVEL-CfA), into a uniform system based on the radial velocities of 854 standard stars in our list. This enabled us to calculate the average weighted radial velocities for more than 25000 HIPPARCOS stars located in the local Galactic spiral arm (Orion arm) with a median error of ±1 km/s. We use these radial velocities together with the stars’ coordinates, parallaxes, and proper motions to determine their Galactic coordinates and space velocities. These quantities, along with other parameters of the stars, are available from the continuously updated Orion Spiral Arm Catalogue (OSACA) and the associated database. We perform a kinematic analysis of the stars by applying an Ogorodnikov-Milne model to the OSACA data. The kinematics of the nearest single and multiple main-sequence stars differ substantially. We used distant (-r ≈ 0.2 kpc) stars of mixed spectral composition to estimate the angular velocity of the Galactic rotation, ω_o = ?25.7 ± 1.2 kms^?1 kpc^?1, and the vertex deviation, l = 13° ± 2°, and detected a negative K effect. This negative K effect is most conspicuous in the motion of A0–A5 giants and is equal to K = ?13.1 ± 2.0 kms^?1 kpc^?1.