A well-behaved class of charged analogue of Durgapal solution

We present a well behaved class of charged analogue of M.C. Durgapal (J. Phys. A, Math. Gen. 15:2637, 1982) solution. This solution describes charged fluid balls with positively finite central pressure, positively finite central density; their ratio is less than one and causality condition is obeyed at the centre. The outmarch of pressure, density, pressure-density ratio and the adiabatic speed of sound is monotonically decreasing, however, the electric intensity is monotonically increasing in nature. This solution gives us wide range of parameter for every positive value of n for which the solution is well behaved hence, suitable for modeling of super dense stars. Keeping in view of well behaved nature of this solution, one new class of solution is being studied extensively. Moreover, this class of solution gives us wide range of constant K (0≤K≤2.2) for which the solution is well behaved hence, suitable for modeling of super dense stars like strange quark stars, neutron stars and pulsars. For this class of solution the mass of a star is maximized with all degree of suitability, compatible with quark stars, neutron stars and pulsars. By assuming the surface density ρ _b =2×10¹⁴ g/cm³ (like, Brecher and Capocaso, Nature 259:377, 1976), corresponding to K=0 with X=0..235, the resulting well behaved model has the mass M=4.03M _Θ , radius r _b =19.53 km and moment of inertia I=1.213×10⁴⁶ g?cm²; for K=1.5 with X=0.235, the resulting well behaved model has the mass M=4.43M _Θ , radius r _b =18.04 km and moment of inertia I=1.136×10⁴⁶ g?cm²; for K=2.2 with X=0.235, the resulting well behaved model has the mass M=4.56M _Θ , radius r _b =17.30 km and moment of inertia I=1.076×10⁴⁶ g?cm². These values of masses and moment of inertia are found to be consistent with the crab pulsars.     

Mass—luminosity relation and pulsational properties of Wolf—Rayet stars

The evolution of Population I stars with initial masses 70M _⊙ ≤ M _ZAMS ≤ 130M _⊙ is considered. The computations were performed under various assumptions about the mass loss rate and were terminated at the phase of gravitational contraction after core helium exhaustion. The mass loss rate at the helium burning phase, ?_3α , is shown to be the main parameter that determines the coefficients of the mass—luminosity relation for Wolf—Rayet stars. Several more accurate mass—luminosity relations for mass loss rates ? = f _3α ?_3α , where 0.5 ≤ f _3α ≤ 3, are suggested, along with the mass—luminosity relation that combines all of the evolutionary sequences considered. The results of the stellar evolution computations were used as initial conditions in solving the hydrodynamic equations describing the spherically symmetric motions of a self-gravitating gas. The outer layers of massive Population I stars are unstable against radial oscillations throughout the helium burning phase. The oscillation amplitude is largest at enhanced carbon and oxygen abundances in the outer stellar layers, i.e., at a lower initial stellar mass M _ZAMS or a lower mass loss rate during the entire preceding evolution. In the course of evolution, the radial oscillation amplitude decreases and the small nonlinearity of the oscillations at M < 10M _⊙ allow the integral of mechanical work W done by an elementary spherical layer of gas in a closed thermodynamic cycle to be calculated with the necessary accuracy. The maximum of the radial dependence of W is shown to be located in layers with a gas temperature T ～ 2 × 10⁵ K, where the oscillations are excited by the iron Z-bump κ-mechanism. Comparison of the radial dependences of the integral of mechanical work W and the amplitude of the radiative flux variations suggests that the nonlinear radial oscillations of more massive Wolf—Rayet stars are also excited by the κ-mechanism.     

Neutron stars in globular clusters: Formation and observational manifestations

Population synthesis is used to model the number of neutron stars in globular clusters that are observed as low-mass X-ray sources and millisecond radio pulsars. The dynamical interactions between binary and single stars in a cluster are assumed to take place only with a continuously replenished “background” of single stars whose properties keep track of the variations in parameters of the cluster as a whole and the evolution of single stars. We use the hypothesis that the neutron stars forming in binary systems from components with initial masses of ～8–12 M _⊙ during the collapse of degenerate O-Ne-Mg cores through electron captures do not acquire a high space velocity. The remaining neutron stars (from single stars with masses >8 M _⊙ or from binary components with masses >12 M _⊙) are assumed to be born with high space velocities. According to this hypothesis, a sizeable fraction of the forming neutron stars remain in globular clusters (about 1000 stars in a cluster with a mass of 5 × 10⁵ M _⊙). The number of millisecond radio pulsars forming in such a cluster in the case of accretion-driven spinup in binary systems is found to be ～10, in agreement with observations. Our modeling also reproduces the observed shape of the X-ray luminosity function for accreting neutron stars in binary systems with normal and degenerate components and the distribution of spin periods for millisecond pulsars.     

Maximal red shifts of neutron stars

It has been recently established that there exists a maximal red shiftz _max for a homogeneous star of given massM. The relationshipz _max(M) is obtained for neutron stars in the mass range 0.71M/M 12.06.     

Calculating Method and Characteristics of the Distribution of Neutron Exposures in the Radiative S-process Nucleosynthesis Model for AGB Stars

A study on the distribution of neutron exposures in a low-mass asymptotic giant branch (AGB) star is presented, according to the s-process nucleosynthesis model with the ¹²C(α, n)¹⁶O reaction occurred under radiative conditions in the interpulse phases. The model parameters, such as the over- lap factor r of two successive convective thermal pulses, the mass ratio q of the ¹³C shell with respect to the He intershell, and the effective mass of ¹³C in the ¹³C shell, vary with the pulse number. Considering these factors, a calculating method for the distribution of neutron exposures in the He intershell has been presented. This method has the features of simplicity and universality. Using this method, the exposure distribution for the stellar model of a star with the mass of 3 M_? and the solar metallicity has been calculated. The results suggest that under the reasonable assumption that the number density of neutrons is uniform in the ¹³C shell, the ?nal exposure distribution approaches to an exponential distribution. For a stellar model with the de?nite initial mass and metallicity, there is a de?nite relation between the mean neutron exposure τ₀ and the neutron exposure Δτ of each pulse, namely τ₀ = 0.434λ(q1, q2, …, qm_max +1, …, r₁, r₂, …, rm_max +1)Δτ, where m_max is the total number of thermal pulses with the third dredge-up episode, and the proportional coeffcient λ(q₁, q₂, …, qm_max +1, …, r₁, r₂, …, rm_max +1) can be determined by an exponential curve ?tting to the ?nal exposure distribution. This new formula quantitatively uni?es the classical model with the s-process nu- cleosynthesis model by means of neutron exposure distribution, and makes the classical model continue to offer guidance and constraints to the s-process nu- merical calculations of stellar models.     

Stellar dynamics in young clusters: the formation of massive runaways and very massive runaway mergers

In the present paper we combine an N-body code that simulates the dynamics of young dense stellar systems with a massive star evolution handler that accounts in a realistic way for the effects of stellar wind mass loss. We discuss two topics.

1. The formation and the evolution of very massive stars (with masses >120 M_⊙) is followed in detail. These very massive stars are formed in the cluster core as a consequence of the successive (physical) collisions of the 10–20 most massive stars in the cluster (this process is known as ‘runaway merging’). The further evolution is governed by stellar wind mass loss during core hydrogen and core helium burning (the WR phase of very massive stars). Our simulations reveal that, as a consequence of runaway merging in clusters with solar and supersolar values, massive black holes can be formed, but with a maximum mass ≈70 M_⊙. In low-metallicity clusters, however, it cannot be excluded that the runaway-merging process is responsible for pair-instability supernovae or for the formation of intermediate-mass black holes with a mass of several 100 M_⊙.
2. Massive runaways can be formed via the supernova explosion of one of the components in a binary system (the Blaauw scenario), or via dynamical interaction of a single star and a binary or between two binaries in a star cluster. We explore the possibility that the most massive runaways (e.g. ζ Pup, λ Cep, BD+43°3654) are the product of the collision and merger of two or three massive stars.

     

A family of exact solutions of Einstein-Maxwell field equations in isotropic coordinates: an application to optimization of quark star mass

In the present paper, we have obtained a class of charged super dense star models, starting with a static spherically symmetric metric in isotropic coordinates for perfect fluid by considering Hajj-Boutros (in J. Math. Phys. 27:1363, 1986) type metric potential and a specific choice of electrical intensity which involves a parameter K. The resulting solutions represent charged fluid spheres joining smoothly with the Reissner-Nordstrom metric at the pressure free interface. The solutions so obtained are utilized to construct the models for super-dense star like neutron stars (ρ _b =2 and 2.7×10¹⁴ g/cm³) and Quark stars (ρ _b =4.6888×10¹⁴ g/cm³). Our solution is well behaved for all values of n satisfying the inequalities \(4 < n \le4(4 + \sqrt{2} )\) and K satisfying the inequalities 0≤K≤0.24988, depending upon the value of n. Corresponding to n=4.001 and K=0.24988, we observe that the maximum mass of quark star M=2.335M _⊙ and radius R=10.04 km. Further, this maximum mass limit of quark star is in the order of maximum mass of stable Strange Quark Star established by Dong et al. (in arXiv:1207.0429v3, 2013). The robustness of our results is that the models are alike with the recent discoveries.     

Unsteady mass outflow from Wolf-Rayet stars

We show that hydrostatically equilibrium models for the thin photospheres of helium stars based on new opacities κ_R (OPAL and OP) can be constructed only for masses M<5M _⊙. The parameter Г=κL/4πGMc, defined as the ratio of light pressure to gravity, exceeds a critical value of 1.0 for larger masses, which must result in mass outflow under light pressure. This mass limit matches the observed lower limit for the masses of Wolf-Rayet stars (M _WR>5M _⊙)), which is an additional argument that the Wolf-Rayet stellar cores are actually helium stars. By solving the equation of radiative transfer in extended atmospheres, we construct a semiempirical model for a WN5 star (M _WN5=10M _⊙)) with a helium core and an expanding envelope, whose physical and geometric parameters are known mainly from light-curve solution for the eclipsing binary V444 Cyg (WN5+06): outflow rate $\dot M \approx 1.0 \times 10^{ - 5} M_ \odot yr^{ - 1} $ , terminal velocity V _∞≈2000 km s^?1, and expanding-envelope optical depth τ_env≈25. The temperature at the outer boundary of the photosphere of a helium star surrounded by such an envelope is approximately 130 kK higher than that in the absence of an envelope, being T _ph≈240 kK. Because of the high temperatures, the absorption coefficients at the corresponding photospheric levels are smaller than those in models with no envelope; therefore, the photosphere turns out to be in hydrostatic equilibrium and stable against light pressure (Г_max≈0.9). As a way out of this conflicting situation (an expanding envelope together with a hydrostatically equilibrium photosphere), we propose a model of discrete mass outflow, which is also supported by the observed cloudy structure of the envelopes in this type of stars. To quantitatively estimate parameters of the nonuniform outflow model requires detailed gas-dynamical calculations.     

Nonlinear pulsations of red supergiants

Excitation of radial oscillations in population I (X = 0.7, Z = 0.02) red supergiants is investigated using the solution of the equations of radiation hydrodynamics and turbulent convection. The core helium burning stars with masses 8M _⊙ ≤ M ≤ 20M _⊙ and effective temperatures T _eff < 4000 K are shown to be unstable against radial pulsations in the fundamental mode. The oscillation periods range between 45 and 1180 days. The pulsational instability is due to the κ-mechanism in the hydrogen and heliumionization zones. Radial pulsations of stars with mass M < 15M _⊙ are strictly periodic with the light amplitude ΔM _bol ≤ 0?5. The pulsation amplitude increases with increasing stellar mass and for M > 15M _⊙ the maximum expansion velocity of outer layers is as high as one third of the escape velocity. The mean radii of outer Lagrangean mass zones increase due to nonlinear oscillations by ≤30% in comparison with the initial equilibrium. The approximate method (with uncertainty of a factor of 1.5) to evaluate the mass of the pulsating red supergiant with the known period of radial oscillations is proposed. The approximation of the pulsation constant Q as a function of the mass-to-radius ratio is given. Masses of seven galactic red supergiants are evaluated using the period-mean density relation.     

Main stages of evolution of matter in the universe. II

A possible scenario for the evolution of the universe following the big bang at t > 10^-5 sec is considered. The necessary conditions that must be present for the formation of stars and stellar systems to be possible are formulated. As a condition for the formation of stars we take kT_s≤ GM_sm_p(3R), and for stellar systems HR ? (GM/R)^1/2, where T_s is the temperature of the cosmic plasma, m_p is the mass of a proton, M_s is the mass of a star, M is the mass of a stellar cluster, R is the radius of these celestial bodies, and H is the bubble parameter for the corresponding time. In accordance with these criteria, we assume that in the course of cosmological expansion, neutron stars should have been formed first (times 2.10^-4 ? t ? 1 sec, densities 0.07 ? ρ_B? 2.10⁴ g/cm³) and then, in chronological order, appeared white dwarfs (t ≈ 10² sec, ρ_B ? 5.10^-3 g/cm³), ordinary stars (t ≈ 4.10⁶ sec, ?_B ≈ 10^-11 g/cm³), galactic nuclei (t ≈ 3.10¹¹ sec, ?_B ≈ 5.10^-19 g/cm³, globular clusters (t ≈ 10¹³ sec, ?_B ≈ 4.10^-21 g/cm³), and galaxies (t ≈ 10¹⁵ sec, ?_B ≈ 10^-24 g/cm³), where ?_B is the average density of ordinary (baryon) matter in the universe. It is shown that a galactic nucleus is a stellar system in statistical equilibrium and consists mainly of neutron stars and white dwarfs. The formation of some pulsars (neutron stars with angular rotation rates 1 < Ω < 200 sec^-1) may occur in a galactic nucleus. Observed pulsars should therefore contain some fraction of neutron stars from the nucleus of the Galaxy that were able to escape it over the relaxation time (the tail of the Maxwell distribution, with star velocities v > v₀, where v₀ is the velocity corresponding to the work function 2GMM_s/R, M being the mass and R the radius of the Galaxy’s nucleus.     

Properties of Wolf-Rayet stars in eclipsing binary systems

The results of investigations of a number of eclipsing Wolf-Rayet binaries are presented. The ‘core’ radiuses, the ‘core’ temperatures and masses of WR stars in the eclipsing WR+OB binary systems V 444 Cyg, CX Cep, CQ Cep, and CV Ser are obtained (see Table I). The results obtained from the light curves analysis of the V 444 Cyg in the range λλ2460 Å-3.5μ give strong evidence for the Beals (1944) model of WR phenomenon. The chromospheric-coronal effects in the WN5 extended atmosphere are not observed up to a distance ofr?20R _⊙. In the Hertzsprung—Russell diagram all the WR stars lie on the left side from the main sequence between the main sequence and the sequence of uniform helium stars (see Figure 9). Their locations are close to those of the helium remnants formed as a result of mass exchange in massive close binary systems. The period variations in the systems V 444 Cyg and CQ Cep have been discovered and a reliable value of the mass loss rateM=10^?5 M _⊙yr^?1 is obtained, for the two WR stars. The results of the photometric and spectroscopic investigations of the WR stars with low mass companions (post X-ray binary stage?) are presented too (see Table II). The masses of the companions are (1–2)M _⊙, their optical luminosity is ～10³⁶, erg s^?1 which implies that these companions cannot be the normal stars. It is possible that these companions are neutron stars accreting from the stellar wind of the WR stars. Low values of the X-ray luminosities of such WR stars with low mass companions imply that the accretion of matter in such systems is distinct from the accretion process in classical X-ray binary systems. It is noted also that the parameters of low massive companions coupled with WR stars are close to those of helium stars.     

The masses of sdB and sdO stars

Mass is a fundamental parameter, but the masses are not well known for most hot subdwarfs. We propose a method of determining the masses of hot subdwarfs. Using this method, we studied the masses of hot subdwarfs from the ESO supernova Ia progenitor survey and the Hamburg quasar survey. The study shows that most of the subdwarf B stars have masses between 0.42 and 0.54 M _⊙, whilst most sdO stars are in the range 0.40～0.55 M _⊙. Comparing our study to the theoretical mass distributions of Han et al. (Mon. Not. R. Astron. Soc. 341:669, 2003), we found that sdO stars with mass less than ～0.5 M _⊙ may evolve from sdB stars, whilst most high-mass (>0.5 M _⊙) sdO stars result from mergers directly.     

Pulsations of intermediate–mass stars on the asymptotic giant branch

Evolutionary tracks from the zero age main sequence to the asymptotic giant branch were computed for stars with initial masses 2 M _⊙ ≤ M _ZAMS ≤ 5 M _⊙ and metallicity Z = 0.02. Some models of evolutionary sequences were used as initial conditions for equations of radiation hydrodynamics and turbulent convection describing radial stellar pulsations. The early asymptotic giant branch stars are shown to pulsate in the fundamental mode with periods 30 day ? Π ? 400day. The rate of period change gradually increases as the star evolves but is too small to be detected (Π?/Π < 10^?5 yr^?1). Pulsation properties of thermally pulsing AGB stars are investigated on time intervals comprising 17 thermal pulses for evolutionary sequences with initial masses M _ZAMS = 2 M _⊙ and 3 M _⊙ and 6 thermal pulses for M _ZAMS = 4 M _⊙ and 5 M _⊙. Stars with initial masses M _ZAMS ≤ 3 M _⊙ pulsate either in the fundamental mode or in the first overtone, whereas more massive red giants (M _ZAMS ≥ 4 M _⊙) pulsate in the fundamental mode with periods Π ? 10³ day. Most rapid pulsation period change with rate ?0.02 yr^?1 ? Π?/Π ? ?0.01 yr^?1 occurs during decrease of the surface luminosity after the maximum of the luminosity in the helium shell source. The rate of subsequent increase of the period is Π?/Π ? 5 × 10^?3 yr^?1.     

On the mass-luminosity relation

The results of a least-squares study of the mass-luminosity relation for eclipsing and visual binary stars consisting of main sequence components are presented. Two methods are discussed. In Part A, the values of the coefficientsA andB in the relation logM=A+BM _Bol are determined. Part B presents a technique which permits the determination of α and β in the relationM=αL ^β, when only the sum of the masses, and not the individual masses of each component, is known. The results and a comparison of the two methods are discussed. It is found that the following massluminosity relation represents the observational data satisfactorily: $$log M = 0.504 - 0.103M_{BOL,} {\text{ }} - \leqslant M_{BOL} \leqslant + 10.5$$ . A discussion of the data and of the possibility that separate mass-luminosity relations may exist for visual and eclipsing binaries is given. The possiblity that more than one mass-luminosity relation is required in the range ?8≤M _Bol ≤+13 is also discussed.     

A rare reflection effect eclipsing sdB+dM binary: 2M 1533+3759

In this paper, we report a rare reflection effect eclipsing sdB+dM binary, 2M?1533+3759. It is the seventh eclipsing sdB+dM binary that has been discovered to date. This system has an orbital period of 0.16177042 day and a velocity semi-amplitude of 71.1 km?s^?1. Using a grid of zero-metallicity NLTE model atmospheres, we derived T _eff=29250 K, log?g=5.58 and [He/H]=?2.37 from spectra taken near the reflection effection minimum. Lightcurve modeling resulted in a system mass ratio of 0.301 and an orbital inclination angle of 86.6°. The derived primary mass for 2M?1533+3759, 0.376±0.055 M _⊙, is significantly lower than the canonical mass (0.48 M _⊙) found for most previously investigated sdB stars. This implies an initial progenitor mass >1.8 M _⊙, at least a main sequence A star and perhaps even one massive enough to undergo non-degenerate helium ignition.     

Mass loss and shell masses of close binary stars

From the values of period changes for 6 close binary stars the mass transfer rate was calculated. Comparing these values M_t with the values of shell masses M_sh, the expression $$lg \dot M_t = \begin{array}{*{20}c} {4.24} \\ { \pm 24} \\ \end{array} + \begin{array}{*{20}c} {0.63} \\ { \pm 6} \\ \end{array} lg M_{sh} $$ Was derived. The analysis of this expression points out the initial character of the outflow of matter, and one may determine the time interval of the substitution of the shell matter. So one may conclude that for a certain mass transfer rate, a certain amount of matter accumulates in the nearby regions of the system. The study of orbital period changes of close binary stellar systems led to the idea that these secular and irregular changes are due to the mass loss and to the redistribution of masses in a close binary. Secular changes of orbital periods are known for approximately 400 eclipsing binary stars. For many stars, including cataclysmic binaries, irregular period changes are known. Thus, the mass loss and the matter redistribution in close binaries are often observed phenomena.     

The Crust Properties and Cooling of Strange Stars

We investigate the influence of the following parameters on the crust properties of strange stars: the strange quark mass (m _s), the strong coupling constant (α_c) and the vacuum energy density (B). It is found that the mass density at the crust base of strange stars cannot reach the neutron drip density. For a conventional parameter set of m _s=200 MeV, B ^1/4 = 145 MeV and α_c = 0.3, the maximum density at the crust base of a typical strange star is only 5.5 × 10¹⁰ gcm^-3, and correspondingly the maximum crust mass is 1.4 ×10^-6 M_⊙. Subsequently, we present the thermal structure and the cooling behavior of strange stars with crusts of different thickness, and under different diquark pairing gaps. Our work might provide important clues for distinguishing strange stars from neutron stars.     

Late stages of the evolution of close binaries

Close binaries can evolve through various ways of interaction into compact objects (white dwarfs, neutron stars, black holes). Massive binary systems (mass of the primaryM ₁ larger than 14 to 15M ₀) are expected to leave, after the first stage of mass transfer a compact component orbiting a massive star. These systems evolve during subsequent stages into massive X-ray binaries. Systems with initial large periode evolve into Be X-ray binaries.Low mass X-ray sources are probably descendants of lower mass stars, and various channels for their production are indicated. The evolution of massive close binaries is examined in detail and different X-ray stages are discussed. It is argued that a first X-ray stage is followed by a reverse extensive mass transfer, leading to systems like SS 433, Cir X1. During further evolution these systems would become Wolf-Rayet runaways. Due to spiral in these system would then further evolve into ultra short X-ray binaries like Cyg X-3.Finally the explosion of the secondary will in most cases disrupt the system. In an exceptional case the system remains bound, leading to binary pulsars like PSR 1913+16. In such systems the orbit will shrink due to gravitational radiation and finally the two neutron stars will coalesce. It is argued that the millisecond pulsar PSR 1937+214 could be formed in this way.A complete scheme starting from two massive ZAMS stars, ending with a millisecond pulsar is presented.Paper presented at the Lembang-Bamberg IAU Colloquium No. 80 on Double Stars: Physical Properties and Generic Relations, held at Bandung, Indonesia 3–7 June, 1983.     

Study of physical properties of spectroscopic binary stars

The main results of a study of a catalogue of physical parameters of 1041 spectroscopic binaries are presented. The distribution of spectroscopic binaries over all main parametersM ₁, a, e, M₁/M₂, P, and certain dependencies between some of them have been found.

1. It appears that among bright (m _v?3 ^m –5 ^m ) stars withM?1M _⊙, about 40% are apparently spectroscopic binaries with comparable masses of components.
2. The majority of spectroscopic binaries with the ratio of the large semiaxis of the orbit to the radius of the primarya/R ₁?20, have eccentricities close to zero. This is probably a consequence of the tidal circularization of orbits of close binaries by viscous friction.
3. The discovery of duplicity of double-line spectroscopic binaries is possible only if the semiamplitude of radial velocityK ₁ is almost 10 times higher than the semiamplitude of the radial velocity of a single-line spectroscopic binary of the same mass.
4. Double-line spectroscopic binaries witha/R _⊙?6(M ₁/M _⊙)^1/3,M ₁≈M ₂?1.5M _⊙ are almost almost absent, and the number of stars witha/R _⊙?6(M ₁/M _⊙)^1/3,M ₁≈1.5M _⊙ is relatively low.
5. The distribution of unevolved SB stars over the large semiaxis may be described by the expression d(N _d/N_t)≈0.2 d loga for 6(M ₁/M _⊙)^1/3?a/R _⊙?100.
6. The intial mass-function for primaries of spectroscopic binaries is the same Salpeter function dN _d≈M ₁ ^?2.35 dM ₁ for 1?M ₁/M _⊙?30.
7. It is possible to explain the observed ratio of the number of single-line spectroscopic binaries to the number of double-line binaries if one assumes that the average initial mass ratio is close to 1 and that the mass of the postmass-exchange remnant of the primary exceeds the theoretical one and/or that half of the angular momentum of the system is lost during mass-exchange.
8. The above-mentioned distributions ofM ₁ anda and assumptions on the mass of remnant and/or momentum loss also allow us to explain the observed shapes of dN/dM, dN/dq, and dN/da distributions after some selection effects are taken into account.

     

The evolution of massive close binaries

The final state of the primaries of binary systems with initial massesM _1i=10M to 15M is derived from the mass of their C/O-cores. The possibility of a second stage of mass transfer towards the secondary is considered. It turns out that the critical mass for the bifurcation is about 14M : stars with larger masses in this range are the progenitors of neutron stars, while the lower mass stars are the ancestors of white dwarfs.Research supported by the National Foundation of Collective Fundamental Research of Belgium (F.K.F.O.) under No. 10303.