Estimation of the accuracy of methods for determining component masses for low-mass X-ray binary systems

Modern modeling of the population of low-mass X-ray binary systems containing black holes applying standard assumptions leads to a lack of agreement between the modeled and observed mass distributions for the optical components, with the observed masses being lower. This makes the task of estimating the systematic errors in the derived component masses due to imperfect models relevant. To estimate the influence of systematic errors in the derived masses of stars in X-ray binary systems, we considered two approximations for the tidally deformed star in a Roche model. Approximating the star as a sphere with a volume equal to that of the Roche lobe leads to slight overestimation of the equatorial rotational velocity V _rot sin i, and hence to slight underestimation of the mass ratio q = M _x /M _v . Approximating the star as a flat, circular disk with constant local line profiles and a linear limb-darkening law (a classical rotational broadeningmodel) is an appreciably cruder approach, and leads to overestimation of V _rot sin i by about 20%. In the case of high values of q = M _x /M _v , this approximation leads to substantial underestimation of the mass ratio q, which can reach several tens of percent. The mass of the optical star is overestimated by a factor of 1.5 in this case, while the mass of the black hole is changed only slightly. Since most estimates of component mass ratios for X-ray binary systems are carried out using a classical rotational broadening model for the lines, this leads to the need for appreciable corrections to (reductions of) previously published masses for the optical stars, which enhances the contradiction with the standard evolutionary scenario for low-mass X-ray binaries containing black holes.     

Evolution of the masses of neutron stars in binary systems

We study the growth of the masses of neutron stars in binary systems due to the accumulation of mass from the optical donors accreted onto the neutron-star surface. Possible scenarios for this accretion are considered. The masses and magnetic-field strengths of radio pulsars derived using population-synthesis methods are compared to the observational data. The population-synthesis analysis indicates that a neutron star can increase its mass from the standard value of m _x ? 1.35M_⊙ to the Oppenheimer-Volkoff limit, m _x ? 2.5M_⊙, via accretion from a companion.     

Evolution of close stellar binaries with black holes under the action of gravitational radiation and magnetic and induced stellar winds of the donor

We show that semi-detached close binary systems with massive (4–25M_⊙) black holes are formed in the evolution of massive stellar binaries in which the initial mass of the primary exceeds ～25M_⊙. The mass exchange in such systems is maintained by the nuclear evolution of the donor and by its magnetic and induced stellar winds. The donor in such systems can be a main-sequence star, subgiant, non-degenerate helium star, or white dwarf. The evolution of corresponding systems with black-hole masses of 10M_⊙ is investigated.     

Logarithmic density range and the percentage of stars in the cores of various subsystems of the globular clusters M56, M12, NGC 6535, NGC 6171, NGC 5466, and M92

Lower limits for the percentages of stars with various luminosities in the cores of six globular clusters are derived using stellar spatial density distributions f(r) to deep limiting B magnitudes obtained earlier. For NGC 6535 and NGC 5466, the logarithmic density range and Kholopov parameters D_f and D_r are also determined. These two parameters are correlated with the mean masses of stars of various subsystems and the total mass (number) of stars in the cluster.     

The mass of the black hole in the X-ray binary M33 X-7 and the evolutionary status of M33 X-7 and IC 10 X-1

We have analyzed the observed radial-velocity curve for the X-ray binary M33 X-7 in a Roche model. We have analyzed the dependence between the component masses and the degree of filling of the optical star’s Roche lobe to obtain the ratio of the masses of the optical star and compact object. For the most probable mass of the optical star, m _v = 70 M⊙, the mass of the compact object is m _x = 15.55 ± 3.20 M⊙. It has been shown that black holes with masses of mx = 15 M⊙ and even higher can form in binaries. We present characteristic evolutionary tracks for binary systems passing through an evolutionary stage with properties similar to M33 X-7-type objects. According to population-synthesis analyses, such binaries should be present in galaxies with masses of at least 10¹¹ M⊙. The present number of such systems in M33 should be of the order of unity. We have also studied the evolutionary status of the X-ray binary IC 10 X-1 with a Wolf-Rayet component, which may contain a massive black hole. The final stages of the evolution of the M33 X-7 and IC 10 X-1 systems should be accompanied by the radiation of gravitational waves.     

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.     

Modeling the mass-metallicity relation for dSph dwarf galaxies of the local group

The relationship between the masses and metallicities of galaxies could be non-monotonic, due to the outflow of matter in these systems. It is shown using a simple, one-zone, chemical-dynamicalmodel that the metallicity should be a non-monotonic function of the mass for spheroidal dwarf galaxies with low masses of M ≤ 10⁸ M _⊙, and a monotonically growing function for galaxies with higher masses. This is consistent with observations.     

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.     

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.     

Masses of stellar black holes and testing theories of gravitation

The paper analyzes the mass distribution of stellar black holes derived from the light and radial-velocity curves of optical stars in close binary systems using dynamical methods. The systematic errors inherent in this approach are discussed. These are associated primarily with uncertainties in models for the contribution from gaseous structures to the optical brightness of the systems under consideration. The mass distribution is nearly flat in the range 4–15M _⊙. This is compared with the mass distribution for black holes in massive close binaries, which can be manifest as ultrabright X-ray sources (L _x >10³⁹ erg/s) observed in other galaxies. If the X-ray luminosities of these objects correspond to the Eddington limit, the black-hole mass distribution should be described by a power law, which is incompatible with the flat shape derived dynamically from observations of close binaries in our Galaxy. One possible explanation of this discrepancy is the rapid evaporation of stellar-mass black holes predicted in recent multi-dimensional models of gravity. This hypothesis can be verified by refining the stellar black-hole mass spectrum or finding isolated or binary black holes with masses below ～3M _⊙.     

A possible relation between the masses of black holes in galactic nuclei and the parameters of the host galaxies

Data on about forty virialized galaxy clusters with bright central galaxies, for which both the galactic velocity dispersion (?? _gal) and the stellar velocity dispersion in the brightest galaxies (??*) are measured, have been used to obtain several approximate relations between ?? _gal, ??*, the absolute B magnitude of the brightest central galaxyM _B ^BCG , and the mass of the central massive black holeM _BH: $\begin{gathered} \log \sigma _* = (0.12 \pm 0.14)\log \sigma _{gal} + (2.1 \pm 0.4), \hfill \\ \log \sigma _* = - (0.15 \pm 0.02)M_B^{BCG} + (0.85 \pm 0.5), \hfill \\ \log M_{BH} = 0.51\log \sigma _{gal} + 7.28. \hfill \\ \end{gathered} $ . These relations can be used to derive crude estimates ofMBH in the nuclei of the brightest galaxies using the parameters of the both host galaxies and the host galaxy clusters. The last relation above confirms earlier suggestions of a quadratic relation between the masses of the coronas of the host systems and the masses their central objects: M _hg ^halo ?? M _cent ² . The relations obtained are consistent with the common evolution of subsystems with different scales and masses formed in the process of hierarchical clustering.     

Evolution of the orbital-plane orientations in the two-protoplanet three-body problem with variable masses

The three body problem with variable masses with two of the bodies being protoplanets is analyzed. The protoplanetary masses are assumed to be much less than the protosolar mass: m ₁(t) ? m ₀(t), m ₂(t) ? m ₀(t). The variations of the body masses over time are assumed to be known. The masses vary isotropically with different rates: ? ₀/m ₀≠? ₁/m ₁, ? ₀/m ₀≠? ₂/m ₂, ? ₁/m ₁≠? ₂/m ₂. The bodies are treated like material points. The problem is described by analogy to the second system of Poincaré elements, based on the equations of motion in a Jacobian coordinate system. Individual aperiodic motions in a quasi-conical cross section are used as the initial, unperturbed, intermediate motions. The expression for the perturbing function does not include terms proportional to third and higher powers of the small masses m ₁ and m ₂. A new analytical expression for the perturbing function analogous to the second system of the Poincaré variables is obtained in the formulation considered using a classical scheme. The analogs of the eccentricities e ₁ and e ₂ and the orbital inclinations i ₁ and i ₂ are considered to be small. The perturbing function accurate to within terms of second order in the small quantities e ₁, e ₂, i ₁, and i ₂ is calculated in a symbolic form using Mathematica package. The equations for the secular perturbations in this protoplanetary three-body problem, with the bodies treated as points with masses varying isotropically with different rates, are obtained. General rigorous analytical solutions to these equations for the secular perturbations describing the evolution of the orbital planes are derived for oblique elements, for arbitrary mass-variation laws. An analog of the Laplace theorem is proved for the orbital inclinations. Analytical formulas are obtained for the temporal variation of the longitudes of the ascending nodes and the inclinations for arbitrary mass variations with different rates. It is shown that the Jacobian node theorem, which is valid in the classical three-body problem with constantmasses, is violated in this problem, unless special initial conditions apply.     

Spectroscopic study of the envelopes of the novae V2467 Cyg and V2491 Cyg

Spectrophotometric observations are used to study the envelopes of the FeII nova V2467 Cyg and the HeN nova V2491 Cyg. The abundances of several elements in the nova envelopes and the envelope masses are estimated. The nitrogen mass abundance in the V2467 Cyg envelope is higher than the solar value by a factor of 186 and the oxygen abundance by the factor of 10. The nitrogen abundance in the envelope of V2491 Cyg exceeds the solar value by a factor of 56, the oxygen abundance by a factor of 12, and the neon abundance by a factor of 8. The masses of the envelopes were estimated to be 8.5×10^?5 M _⊙ for V2467 Cyg and 1.5×10^?5 M _⊙ for V2491 Cyg. These envelope elemental abundances and masses are in good agreement with those of low-mass CO white dwarfs (0.8 M _⊙) and ONe white dwarfs (1.15 M _⊙).     

The mass of the compact object in the low-mass X-ray binary 2S 0921-630

We interpret the observed radial-velocity curve of the optical star in the low-mass X-ray binary 2S 0921-630 using a Roche model, taking into account the X-ray heating of the optical star and screening of X-rays coming from the relativistic object by the accretion disk. Consequences of possible anisotropy of the X-ray radiation are considered. We obtain relations between the masses of the optical and compact (X-ray) components, m _v and m _x , for orbital inclinations i = 60°, 75°, and 90°. Including X-ray heating enabled us to reduce the compact object’s mass by ～0.5–1 M _⊙, compared to the case with no heating. Based on the K0III spectral type of the optical component (with a probable mass of m _v ? 2.9 M _⊙), we concluded that m _x ? 2.45?2.55 M _⊙ (for i = 75°?90°). If the K0III star has lost a substantial part of its mass as a result of mass exchange, as in the V404 Cyg and GRS 1905+105 systems, and its mass is m _v ? 0.65?0.75 M _⊙, the compact object’s mass is close to the standard mass of a neutron star, m _x ? 1.4 M _⊙ (for i = 75°?90°). Thus, it is probable that the X-ray source in the 2S 0921-630 binary is an accreting neutron star.     

Rotational mass loss by Be stars

We have carried out a numerical study of rotational mass loss by rapidly rotating Be stars assuming preservation of rigid-body rotation during their main-sequence evolution. Evolutionary models are computed for stars with solar chemical composition and initial masses of 3, 10 and M_⊙. As a result of their rapid initial rotation, these stars can lose one to four percent of their initial mass during the main-sequence stage. The amount of mass lost increases with the initial mass of the star. The matter lost by Be stars can form gas-dust disks with masses comparable to the masses of planets, which, in principle, makes possible the formation of planetary systems around such stars.     

The Fragmentation in a Protoplanetary Disk and the Properties of the Resultant Self-Gravitating Clumps

We study the fragmentation properties in the protoplanetary disk and properties of the resultant self-gravitating clumps using our newly constructed disk model. Our disk model includes the mass inflall term from a molecular cloud core and the photoevaporation winds effect. We adopt the conventional fragmentation criterion to judge whether a protoplanetary disk can fragment. In this work, we follow our previous work to investigate the properties of the resultant self-gravitating clumps. In our calculation, the initial masses of the resultant self-gravitating clumps lie in the range of tens of M_J to more than one hundred of M_J, where M_J is the Jupiter mass. These initial masses can seemingly account for the masses of extrasolar planets in magnitude. We also calculate the subsequent gas accretion of clumps in 1.27 × 10⁴ yr after the formation of self-gravitating clumps. We find that the subsequent gas accretion of self-gravitating clumps is very efficient, and the clump masses grow to hundreds of M_J and the physical radii R_c of clumps increase to about 10 AU. Additionally, we also calculate the orbital migration of clumps. We find that most clumps have short migration timescale to be accreted onto the protostar, and only a small fraction of clumps have long migration timescale (>10⁶ yr) to successfully become gas giant planets. These results are consistent with previous studies.     

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.     

Minimum velocity dispersion in stable stellar disks. Numerical simulations

N-body dynamical simulations are used to analyze the conditions for the gravitational stability of a three-dimensional stellar disk in the gravitational field of two rigid spherical components—a bulge and halo whose central concentrations and relative masses vary over wide ranges. The number of point masses N in the simulations varies from 40 to 500 000 and the evolution of the simulated systems is followed over 10–20 rotation periods of the outer edge of the disk. The initially unstable disks are heated and, as a rule, reach a quasi-stationary equilibrium with a steady-state radial-velocity dispersion c_r over five to eight turns. The radial behavior of the Toomre stability parameter Q_T(r) for the final state of the disk is estimated. Simple models are used to analyze the dependence of the gravitational stability of the disk on the relative masses of the spherical components, disk thickness, degree of differential rotation, and initial state of the disk. Formal application of existing, analytical, local criteria for marginal stability of the disk can lead to errors in c_r of more than a factor of 1.5. It is suggested that the approximate constancy of Q_T?1.2–1.5 for r?(1–2)×L (where L is the radial scale of disk surface density), valid for a wide range of models, can be used to estimate upper limits for the mass and density of a disk based on the observed distributions of the rotational velocity of the gaseous component and of the stellar velocity dispersion.     

Physical and dynamical parameters of the multiple system HD 222326

We present the results of our study of the physical and dynamical parameters of the multiple system HD 222326. A new method for determining the individual radial velocities of components in wide binary and multiple systems in the case of small radial-velocity differences (δV _r ≤ the FWHMfor the line profiles) is suggested and tested for both model systems and the binary HD 10009. This testing yielded the component radial velocities V _{r 1,2} for HD 10009, enabling us to derive the center-of mass velocity, V _γ, for the first time. We determined the radial velocities of the components of HD 222326 from high-resolution spectra, and refined the orbital parameters of the subsystems using speckle-interferometric observations. A combined spectroscopic and speckle interferometric analysis enabled us to find the positions of the components in the spectral type—luminosity diagram and to estimate their masses. It is likely that the components are all in various evolutionary stages after leaving the main sequence. We analyzed the dynamical evolution of the system using numerical modeling in the gravitational three-body problem and the known stability criteria for triple systems. The system is probably stable on time scales of at least 10⁶ years. The presence of a fourth component in the system is also suggested.     

Supermassive black holes and central star clusters: Connection with the host galaxy kinematics and color

The relationship between the masses of the central, supermassive black holes (M _bh) and of the nuclear star clusters (M _nc) of disk galaxies with various parameters galaxies are considered: the rotational velocity at R = 2 kpc V ₍₂₎, the maximum rotational velocity V _max, the indicative dynamical mass M ₂₅, the integrated mass of the stellar populationM _*, and the integrated color index B-V. The rotational velocities andmasses of the central objects were taken from the literature. ThemassM _nc correlatesmore closely with the kinematic parameters and the disk mass than M _bh, including with the velocity V _max, which is closely related to the virial mass of the dark halo. On average, lenticular galaxies are characterized by higher massesM _bh compared to other types of galaxies with similar characteristics. The dependence of the blackhole mass on the color index is bimodal: galaxies of the red group (red-sequence) with B-V >0.6–0.7 which are mostly early-type galaxies with weak star formation, differ appreciably from blue galaxies, which have higher values of M _nc and M _bh. At the dependences we consider between the masses of the central objects and the parameters of the host galaxies (except for the dependence of M _bh on the central velocity dispersion), the red-group galaxies have systematically higher M _bh values, even when the host-galaxy parameters are similar. In contrast, in the case of nuclear star clusters, the blue and red galaxies form unified sequences. The results agree with scenarios in which most red-group galaxies form as a result of the partial or complete loss of interstellar gas in a stage of high nuclear activity in galaxies whose central black-hole masses exceed 10⁶?10⁷ M _⊙ (depending on the mass of the galaxy itself). The bulk of disk galaxies with M _bh > 10⁷ M _⊙ are lenticular galaxies (types S0, E/S0) whose disks are practically devoid of gas.