Evolution of massive close binaries and formation of neutron stars and black holes

Main results of computations of evolution for massive close binaries (10M +9.4M , 16M +15M , 32M +30M , 64M +60M ) up to oxygen exhaustion in the core are described. Mass exchange starting in core hydrogen, shell hydrogen and core helium burning stages was studied. Computations were performed assuming both the Ledoux and Schwarzschild stability criteria for semiconvection. The influence of UFI-neutrino emission on evolution of close binaries was investigated. The results obtained allow to outline the following evolutionary chain: two detached Main-Sequence stars — mass exchange — Wolf-Rayet star or blue supergiant plus main sequence star — explosion of the initially more massive star appearing as a supernova event — collapsed or neutron star plus Main-Sequence star, that may be observed as a runaway star — mass exchange leading to X-rays emission — collapsed or neutron star plus WR-star or blue supergiant — second explosion of supernova that preferentially disrupts the system and gives birth to two single high spatial velocity pulsars.Numerical estimates concerning the number and properties of WR-stars, pulsars and X-ray sources are presented. The results are in favour of the existence of UFI-neutrino and of the Ledoux criterion for describing semiconvection. Properties of several well-known X-ray sources and the binary pulsar are discussed on base of evolutionary chain of close binaries.     

Supernova SN 1961v: An explosion of a very massive star

An investigation of the outburst of the unique supernova SN 1961v in the galaxy NGC 1058 is carried out. An analysis of hydrodynamic models of supernova outbursts and a comparison with a considerable body of observational data on SN 1961v clearly show that the SN 1961v phenomenon is an explosion of a very massive star (VMS) of a mass of 2000M and a radius of about 100R which results in expelling the envelope with a kinetic energy of 1.8×10⁵² erg. The light curve of SN 1961v (Figures 1, 7a, and 7b) furnishes direct evidence for the heterogeneity of the presupernova interior. The chemical composition profile produced during the evolution of the VMS and in the final explosion must have a number of the essential features (Figure 11). In particular, hydrogen has to be underabundant relative to the solar content and distributed in a specific manner throughout the star. At the late stages from February, 1963 to February, 1967, the light curve of SN 1961v (Figure 1) may be accounted for by the interaction of the expelled envelope with the stellar wind of the presupernova. The latest observations of SN 1961v in 1968 and 1970 are virtually those of a giantHii region created by the VMS before the explosion. Two astrophysical phenomena-the peculiar outburst of SN 1961v and the most luminous object R136a in the Large Magellacnic Cloud (LMC) which reveals a striking similarity with the presupernova-are evidence for the existence of VMSs. The evolution of VMSs similar to the object R136a may be terminated by explosions like the outburst of SN 1961v. Such explosions give rise to the formation of energetic supernova remnants whose examples may be the Cygnus superbubble and the supergiant shells in the LMC. A comparison of the internal structure of the presupernova with the available evolutionary calculations allows one to conclude that the influence of mass loss on the evolution of VMSs is negligible.     

The neutrino radiation for a hot neutron star formation and the envelope outburst problem

With the equations of neutrino heat conductivity being used, the neutrino light curve is calculated for the spherically symmetrical collapse of an iron-oxygen 2M star (Figure 1) up to the formation of a hot hydrostatically equilibrium neutron star. The total energy, radiated in the form of muon and electron neutrinos, is 5.8×10⁵³ erg (0.16Mc ²). The mean neutrino particle energy is 12 MeV for all the time the collapse proceeds. The maximum neutrino luminosity value is equal to 3×10⁵³ erg s^–1. For a 10M star collapse, the luminosity maximum 3×10⁵⁴ erg s^–1 takes place just at the moment of the formation of a black hole inside the collapsing star. The total radiated energy in this case is about 0.08Mc ². The set of calculations, allowing for the deposition of momentum by means of neutrino-nuclear coherent scattering, brings us to a conclusion that the envelope outburst is only possible if the scattering cross-section is 50 times larger than the value experimentally accepted (inequality 20)).     

Low-Mass Quark Stars or Quark White Dwarfs

An equation of state is considered that, in superdense nuclear matter, results in a phase transition of the first kind from the nucleon state to the quark state with a transition parameter > 3/2 ( = _Q /( _N + P ₀/c ²)). A calculation of the integrated parameters of superdense stars on the basis of this equation of state shows that on the stable branch of the dependence of stellar mass on central pressure (dM/dP _c > 0), in the low-mass range, following the formation of a tooth-shaped break (M = 0.08 M , R = 200 km) due to quark formation, a new local maximum with M _max = 0.082 M and R = 1251 km is also formed. The mass and radius of the quark core of such a star turn out to be M _core = 0.005 M and R _core = 1.7 km, respectively. Mass accretion in this model can result in two successive transitions to a neutron star with a quark core, with energy release like supernova outbursts.     

A massive circumstellar envelope around the type-IIn supernova SN 1995G

We model the interaction of the supernova SN 1995G with a dense circumstellar (CS) gas in a thin-shell approximation. A model fit to the observed bolometric light curve combined with data on the supernova expansion velocity gives estimates for the density, mass (≈1 M_⊙), and age (≈8 yr) of the CS envelope. The determined CS-envelope density is shown to be virtually independent of the assumed mass of the supernova envelope because of the high CS-gas density at which the forward shock wave is essentially radiative. The derived CS-envelope density is consistent with the Hα luminosity and with the presence of distinct Thomson scattering in the red wing of this line. The mass of the CS envelope together with its expansion velocity and age indicate that the CS envelope was ejected by the presupernova eight years before the supernova explosion through violent energy release (～6×10⁴⁸ erg).     

Supernova remnants and content in high mass stars of the giant Hii region NGC 604

An optical survey of NGC 604, the brightestHII complex in M33, has lead to the discovery of a single supernova remnant, in agreement with the predominantly thermal radio spectrum of the region. For the same NGC 604, we have derived from the H flux and the evolutionary tracks of the stars, the population and the birth-rates of Main Sequence, high-mass stars. The computation is in agreement with the observations of the upper part of the stellar luminosity function. From the stellar population one predicts in the region a number of remnants significantly larger than 1. To explain this discrepancy, it is suggested the presence in NGC 604 of at least one overluminous, hot star (M _B9), which contributes to a large fraction of the ionizations.On leave from the Asiago Observatory.On leave from Istituto di Fisica Generale, Università di Torino.     

The bolometric light curve of SN 1993J and the nature of its progenitor

We have constructed the bolometric light curve of SN 1993J based on UBVRI(JHK) photometric data obtained from various sources and assumingA _V = 0 and a distance modulus of 27.6. Effective temperatures and photosphere radius at various times have been obtained from detailed blackbody fits. The bolometric light curve shows two maxima. The short rise time to the second maximum, and the luminosities at the minimum and the second maximum are used to constrain the properties of the progenitor star. The total mass of the hydrogen envelope M_H, in the star is found to be ≲ 0.2 M_⊙ at the time of explosion, and the explosion ejected about 0.05 M_⊙ of Ni⁵⁶. Thin hydrogen envelope combined with a sufficient presupernova luminosity suggest that the exploding star was in a binary with a probable period range of 5yr ≤P _orb ≤ 11yr.     

The accretion of matter by a collapsing star in the presence of a magnetic field. II. Selfconsistent stationary picture

The stationary two-dimensional magnetohydrodynamic solution for the accretion of the matter without pressure into a gravitating centre of a black hole is obtained. It is assumed that the magnetic field far from the collapsed star is homogeneous and its influence on the flow is negligible. Around the star, at the plane perpendicular to the direction of the magnetic field, the dense quasistationary disc is formed, the structure of which in a large extent is determined by dissipation processes. The structure is calculated for (a) a laminar disc with the Coulomb mechanism of dissipation; and (b) a turbulent disc.The estimations of the parameters of the shock which result from the infall of the matter onto the disc are given. In the last section the numerical estimation and approximate character of the radiation spectrum of the disc and the shock are obtained for two cases of 10M and 10⁵ M . The luminosity of collapsed objects withM=10M appears to be about solar, thus its observation is possibly only at the distances less than 300–1000 pc. The collapsed objects in the Galaxy withM=10⁵ M could constitute very bright sources in spectral regions from optical up to X-ray. The spectra of a laminar and a turbulent disc for 10M black hole are very different. The laminar disc radiates primarily in the ultraviolet. The turbulent disc radiates a large part of its flux in the infrared. Therefore, one cannot exclude the possibility that some of the galactic infrared star-like sources are individual black holes in the accretion state.     

Parameters of the classical type-IIP supernova SN 1999em

Based on observations of SN 1999em, we determined the physical parameters of this supernova using hydrodynamic calculations including nonequilibrium radiative transfer. Taking the distance to SN 1999em estimated by the expanding photosphere method (EPM) to be D = 7.5 Mpc, we found the parameters of the presupernova: radius R = 450R_⊙, mass M = 15M_⊙, and explosion energy E = 7 × 10⁵⁰ erg. For the distance D = 12 Mpc determined from Cepheids, R, M, and E must be increased to the following values: R = 1000R_⊙, M = 18M_⊙, and E = 10⁵¹ erg. We show that one cannot restrict oneself to using the simple analytical formulas relating the supernova and presupernova parameters to obtain reliable parameters for type-IIP presupernovae.     

Supernova models with slow energy pumping and galactic supernova remnants

The study of supernova (SN) models with slow energy pumping is continued. At maximum luminosity the main characteristics of a SN are shown to be independent of the initial structure of the model (Table I, Figure 1). However, they depend on the massM _e of the envelope, and on the intensity of energy pumpingL , with an increase ofM _e leading qualitatively to the same changes in the SN parameters as a decrease inL (Table I, Figures 2 and 3). A simple relationship connecting the important SN parameters is obtained (Equation (6)). From the inflection of the color indexB-V curve, the possibility of deriving the characteristic time of energy pumping with intensityL 10⁴⁴ erg s^–1 is pointed out. The comparison of the extragalactic type I SN observations with the results of calculations leads to the estimate ofM _e 0.3–0.7M.An investigation of the galactic type I SN remnants is carried out (Table III). The estimate ofM _e 0.2–0.3M is obtained for the remnants of supernovae SN 1006, SN 1572, and SN 1604. It completely fits the results for the extragalactic type I SNs. The total initial mass of SN 1604 presupernova was shown to be at least about 7M .It was established that the Crab nebula resulted from the outburst of a peculiar SN. The unique properties of such SNs, including SN 1054, are due to the low intensity of energy pumping (L 10⁴² erg s^–1). The mass of the envelope of the Crab nebula is evaluated to beM _e 0.7M . SN 1054 was shown to have m _max ^v =–4 _m ^. 0 at maximum luminosity.     

Neutrino mass,luminosity variation,and spectrum of SN1987A

In this paper, the neutrino mass has been determined from SN1987a observation in a manner that the simultaneity of neutrino emission is not regarded as the starting point, but is itself defined through the analysis by Monte-Carlo simulation. The result is that the neutrino mass lies in 3–4 eV, possiblym _v 3.6 eV. Neutrino luminosity variation and neutrino spectrum are also obtained. Comparison with theories gives further support to the mass determination, and also predicts the mass of progenitor star of SN1987a to be in the range of 12–25M .The project supported by the National Natural Science Foundation of China.     

On the pre-supernova parameters of the Hercules X-1 system

On the assumption that, at the moment of the supernova explosion, the presently nondegenerate companion of Her X-1 was a Main-Sequence star with polytropic indexn=3, it is found that the effects of impact and ablation cannot have removed more than about 10% of its pre-supernova mass. A remnant mass for a neutron star of 1.4M was adopted. The effects of impact and ablation were calculated in the manner given by Wheeleret al. (1975). Depending on the ejection velocity of the SN shell, it is found that with a pre-SN mass of 2.2M (i.e., the maximum possible value) for the non-degenerate component, the initial binary period was in the range 2–4 days.     

The parameters of planetary nebulae and their central stars derived from observations

On the basis of data on planetary nebula (PN) central star temperatures obtained by measurements in the ultraviolet (UV) range, the empirical calibration dependence between the number of Lyman photons emitted by a central starS and PN diameterD, is constructed. The temperatures of 118 PN central stars are estimated with this dependence. It is shown that the central star masses are distributed in a wide interval from 0.5 to 1.2M . About 60% of all stars have masses <0.6M , about 25% have masses >0.6M and the remainder have masses 0.6M . The averaged empirical tracks of evolution of low-mass (<0.6M ) and massive (>0.6M ) central stars differing considerably from each other are constructed. It is shown that the majority of central stars may possess hot chromospheres (T>2×10⁵ K) which spread for several tens of radii of the central star. The PN originates as a result of ionization of the matter ejected by a red giant at the superwind stage. The cause for this ionization is the UV radiation of the PN central star.     

Electromagnetic bursting of a rapidly rotating magnetized neutron star in a close binary system

We examine a possible manifestation of the electromagnetic activity of a magnetized, rotating neutron star in a binary system. Accreting matter from the companion is initially accumulated at the magnetosphere. When the accumulated mass is such that the inflow can start, together with the accretion flare there will be a burst due to the closure of electric currents. The luminosity associated to the latter effect may be as large as 10⁴² erg/s, if a neutron star possesses the following characteristics: massM =M , period of rotationP = 5 ms, magnetic fieldB ₀ = 10¹² G, and radiusr ₀ = 10⁶ cm. The electromagnetic activity might be relevant for understanding soft gamma ray repeaters.     

Models of pulsating low-massive yellow supergiants

The nonlinear self-excited oscillations of the envelopes of low-massive highly luminous stars are described. The parameters for these models wereM=0.8M ,M _bol=–5.5, –5.84 mag,T _eff=4500, 5000, 5500 K. The oscillations have been found to consist of the standing wave pulsation near the envelope bottom and running waves in outer layers. The ratio of the standing wave frequency _s to the average frequency of the running waves _r increases with the stellar luminosity: _s / _r =1.7 whenM _bol=–5.5 mag and _s / _r =2.4 whenM _bol=–5.84 mag. The frequency of oscillations near the photosphere is found to be in close agreement with the critical frequency for running waves. Mass loss from these stars is caused by shocks. It has been shown that agreement between FG Sge's period change observed during the last decade and the period-luminosity relation for double shell stars takes place when FG Sge's luminosity isM _bol=–5.96 mag.     

Mass accretion onto massive white dwarfs

We have studied the thermonuclear runaways which develop on white dwarfs of 1.205M and 1.358M accreting hydrogen rich material at 10^–10 M yr^–1. It is found that ignition of this material occurs at densities in excess of about 10⁴ gm cm^–3 and that the critical accumulated mass required to initiate the runaway is 0.7(1.5)×10^–4 M for a 1.358(1.205)M white dwarf.     

Hidden mass in the solar neighborhood

It is suggested that the minimum mass of a star at the time of its formation is approximately 0.01M . Making use of this fact and the stellar mass functionF(M) M ^–, it is found that the hidden mass (or the missing mass) in the solar neighborhood may be explained by the presence of a large number of invisible stars of very low mass (0.01M M<0.07M ).     

Evolutionary models of nucleosynthesis in the galaxy

A model of the galaxy is constructed and evolved in which the integrated influence of stellar and supernova nucleosynthesis on the composition of the interstellar gas is traced numerically. Our detailed assumptions concerning the character of the matter released from evolving stars and supernovae are guided by the results of recent stellar evolutionary calculations and hydrodynamic studies of supernova events. Stars of main sequence mass in the range 4M8M are assumed to give rise to supernova events, leaving remnants we identify with neutron stars and pulsars and forming both the carbon-to-iron nuclei and ther-process heavy elements in the explosive ejection of the core material. For more massive stars, we assume the core implosion will result in the formation of a Schwarzschild singularity, that is, a black hole or collapsar. The straightforward assumptions (1) that the gas content of the galaxy decreases exponentially with time to its present level of 5% and (2) that the luminosity function characteristic of young clusters and the solar neighborhood is appropriate throughout galactic history, lead to the prediction that 20% of the unevolved stars of approximately one solar mass (M ) in the galaxy today should have metal compositionsZ0.1Z . As Schmidt has argued from similar reasoning, this is quite inconsistent with current observations; an early generation dominated by more massive stars—which would by now have evolved—is suggested by this difficulty. Many of these massive stars, according to our assumptions, will end their lives as collapsed black hole remnants. It is difficult to visualize an epoch of massive star formation in the collapsing gas cloud which formed our galaxy which would enrich the gas rapidly enough to account for the level of heavy element abundances in halo population stars; we have therefore proposed a stage of star formation which is entirely pregalactic in character. We suggest that the Jeans' length-sized initial condensations in the expanding universe discussed by Peebles and Dicke may provide the appropriate setting for this first generation of stars. Guided by these considerations, and by the need for a substantial quantity of unseen mass to bind our local group of galaxies, we have constructed a model of the galaxy in which this violent early phase of massive star formation produces both (1) approximately 25% of the level of heavy elements observed in the solar system and (2) an enormous unseen mass in the form of black holes. The implications of our model for other features of the galaxy, including supernova nucleosynthesis, the cosmic ray production of the light elements, and cosmochronology, are discussed in detail.     

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.