next up previous contents index
Next: Star Formation Rate and Up: No Title Previous: Calculation of Evolutionary Tracks

Initial Distributions of Binary Parameters


To start calculations, we choose the distributions of the initial binary parameters: the mass of the primary ZAMS component, tex2html_wrap_inline10264 , binary mass ratio, tex2html_wrap_inline10266 , orbital separation a. Zero initial eccentricity is assumed.

The distribution by orbital separations is deduced from observations (Abt 1983[3]),



whereas the mass ratio distribution is not. Especially important from the evolutionary point of view is how different are the initial masses of the components (see e.g. Trimble, 1983[193]). The zero assumption often made is that the mass ratio distribution has a flat shape, that is binaries with a high mass ratio occur as frequently as those with equal masses. As we cannot reliably determine this distribution directly from observations (due to a number of selection effects), we have parametrized it by a power-law shape, assuming the primary mass to obey Salpeter's power law: 



We should note that some apparently reasonable distributions - such as both the primary and secondary mass obeying Salpeter's law, or ``hierarchical'' distributions involving the assumption that the total binary mass and primary mass are distributed according to the Salpeter mass function  - all yield essentially flat distributions by mass ratio (i.e. with our parameter tex2html_wrap_inline10270 ).

We assume that a NS with a mass of tex2html_wrap_inline10272 is formed as a result of the collapse  of a star with the core mass before the collapse tex2html_wrap_inline10274 - tex2html_wrap_inline10276 corresponding to the initial mass range tex2html_wrap_inline894510- tex2html_wrap_inline10280 (we take into account that massive star can loose more than tex2html_wrap_inline8945 10-20 percent of their initial mass during evolution through a strong stellar wind (de Jager, 1980[41])).

As discussed in the previous section, magnetic fields of rotating compact objects (NS and WD) largely define the evolutionary stage of the compact object in a binary system. The initial magnetic dipole moment of the newborn NS is taken to have a flat distribution for a wide range of magnetic moments


and the initial rotational period of the NS is assumed to be 1 ms. Arguments pointing to the existence of a super-Schwinger ( tex2html_wrap_inline10284 - tex2html_wrap_inline10286  G) magnetic field strength in some NS follows from different studies: first, as early as 1982, the spin-up  and spin-down  analysis performed for a number of X-ray pulsar periods pointed indirectly to such high magnetic field strengths (Lipunov, 1982c[100]); presently, there are some indications that the soft gamma-ray repeaters (SGR) are young NS with a surface magnetic field  strength as high as tex2html_wrap_inline10286 G (Kulkarni 1995[90]).

It is not definitely clear as yet whether the magnetic field of NS decays  or not (for a comprehensive review see, Chanmugam, 1992[29]). Here we assume that the magnetic fields of NS decay exponentially on a times-scale of tex2html_wrap_inline9162 yr. A radiopulsar is assumed to be ``on'' until its period P (in seconds) has reached a ``death-line'' value defined from the relation tex2html_wrap_inline10294 , where tex2html_wrap_inline10296 is the dipole magnetic moment in units of tex2html_wrap_inline9563  G cm tex2html_wrap_inline8875 .

In Figure 15 we plot the computed and observed p-L diagram for tex2html_wrap_inline8873  G cm tex2html_wrap_inline8875 , which shows that the characteristics of the simulated X-ray pulsars with accepted parameters are close to the observed values.


Figure: Computed and observed p-L diagram for X-ray pulsars with tex2html_wrap_inline8873  G cm tex2html_wrap_inline8875 (Kornilov and Lipunov, 1983b). 

These calculations show, for example, that the choice of the quantity tex2html_wrap_inline10314  G cm tex2html_wrap_inline8875 would contradict the existence of a large number of the long-period X-ray pulsars (with periods >100 s). On the other hand, the choice tex2html_wrap_inline10320  G cm tex2html_wrap_inline8875 has no impact on the distribution of pulsars on the diagram since they normally fall into class G (georotators). 

The statistical calculations performed by Lipunov and Postnov (1988)[116] showed that the flat spectrum of WD magnetic moments does not contradict the observed properties of WD in CVs, polars  and intermediate polars.  The initial magnetic moments of WD are taken to be distributed as


If in the course of evolution a WD with magnetic moment tex2html_wrap_inline10324 collapses to form a NS (accretion induced collapse, AIC),  its magnetic moment obviously scales as tex2html_wrap_inline10326 . Rotational periods of WD are calculated from angular momentum conservation of a degenerate stellar core producing WD.

The mass limit for NS (the Oppenheimer-Volkoff limit)  is tex2html_wrap_inline10328 , which corresponds to a hard equation of state of NS matter.

The most massive stars are assumed to leave BH behind after the collapse,  provided that the progenitor mass before the collapse was tex2html_wrap_inline10330 . The mass of the BH are calculated as tex2html_wrap_inline10332 , where the parameter tex2html_wrap_inline10334 (see, Lipunov et al., (1994b)[124] for justification of this choice).

Another important point is that the collapse  of a massive star into a NS can be asymmetrical, so that the newborn NS can acquire an additional, presumably randomly oriented in space, kick velocity w.  In view of the importance of this parameter, we study the influence of this asymmetry on the evolutionary calculations in detail in a separate section below.

The location of the binary in the Galaxy needed for calculation of the gravitational wave stochastic background produced by binary systems is also calculated after the old stellar galactic distribution in the form


where tex2html_wrap_inline10340 kpc and tex2html_wrap_inline10342 kpc.

next up previous contents index
Next: Star Formation Rate and Up: No Title Previous: Calculation of Evolutionary Tracks

Mike E. Prokhorov
Sat Feb 22 18:38:13 MSK 1997