Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://lnfm1.sai.msu.ru/~marat/articles/eng/Vr_forma_eng.pdf Äàòà èçìåíåíèÿ: Tue Nov 22 14:42:28 2005 Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 23:28:19 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: supernova

Astronomy Reports, Vol. 49, No. 10, 2005, pp. 801­813. Translated from Astronomicheski Zhurnal, Vol. 82, No. 10, 2005, pp. 900­913. i Original Russian Text Copyright c 2005 by Abubekerov, Antokhina, Cherepashchuk.

Dependence of the Absorption-Line Profiles and Radial-Velocity Curve of the Optical Star in an X-ray Binary on the Orbital Inclination and Component-Mass Ratio
´ M. K. Abubekerov, E. A. Antokhina, and A. M. Cherepashchuk
Sternberg Astronomical Institute, Universitetski i pr. 13, Moscow, 119992 Russia
Received December 18, 2004; in final form, May 18, 2005

Abstract--Theoretical absorption-line profiles and radial-velocity curves for tidally deformed optical stars in X-ray binary systems are calculated assuming LTE. The variations in the profile shapes and radialvelocity curve of the optical star are analyzed as a function of the orbital inclination of the X-ray binary system. The dependence of the shape of the radial-velocity curve on the orbital inclination i increases with decreasing component-mass ratio q = mx /mv . The integrated line profiles and radial-velocity curves of the optical star are calculated for the Cyg X-1 binary, which are then used to estimate the orbital inclination and mass of the relativistic object: i < 43 and mx = 8.2-12.8 M . These estimates are in good agreement with earlier results of fitting the radial-velocity curve of Cyg X-1 using a simpler model (i < 45 , mx = 9.0-13.2 M ). c 2005 Pleiades Publishing, Inc.

1. INTRODUCTION The optical component in an X-ray binary system is tidally deformed and has a complex temperature distribution on its surface due to the effects of gravitational darkening and X-ray heating. These effects of the interaction between the components give rise to orbital variability of the absorption-line profiles of an optical star. The orbital variability of the absorption profiles leads to a dependence of both the half-amplitude and the shape of the radial-velocity curve on the orbital inclination i and the componentmass ratio q = mx /mv . The dependence of the shape of the radial-velocity curve of a tidally deformed optical star on the parameters of a close binary system was first pointed out by Sofia and Wilson [1]. Antokhina and Cherepashchuk [2] and Shabaz [3] proposed a new method for determining the component-mass ratio q = mx /mv and orbital inclination i of an X-ray binary based on the orbital variability of the absorptionline profiles in the spectrum of the optical star. Abubekerov et al. [4] estimated the orbital inclination of the X-ray binary Cyg X-1 using a high-accuracy observational radial-velocity curve. Here, we present the results of theoretical modeling of the radial-velocity curve of the optical star in a Roche model assuming LTE for various orbital inclinations, for X-ray systems with low-, moderate-, and high-mass optical stars. We analyze variations in the radial-velocity curve with variations of i for various values of q , together with the corresponding variations in the H absorption-line profile.

2. SYNTHESIS OF THE RADIAL-VELOCITY CURVES FOR THE OPTICAL STAR The synthesis of the theoretical absorption-line profiles and radial-velocity curves for the optical star in an X-ray binary system was carried out using the algorithm described in detail by Antokhina et al. [5, 6]. We will briefly summarize the basis of this method here. In the Roche model, the X-ray binary system consists of an optical star and a pointlike X-ray source. The star is tidally deformed and has a nonuniform surface-temperature distribution due to the effects of gravitational darkening and heating of the stellar surface by the X-ray emission of the relativistic object. The surface of the optical star was divided into 2600 area elements, for each of which we calculated the emergent local radiation assuming LTE. Each area element corresponds to a local temperature Tloc , local gravitational acceleration gloc , and local value lo of the parameter kx c , which is equal to the ratio of the incident X-ray flux and the outgoing radiation flux without allowance for external irradiation of the atmosphere. Using these parameter values at a given point of the surface, a model for the atmosphere is calculated by solving the equation of line radiative transfer in the presence of incident external X-ray radiation. In this way, we can compute the intensity of the outgoing radiation in the line and continuum for each local area element. At different phases of the orbital period, the contributions of the areas to the total radiation are summed taking into account Doppler

1063-7729/05/4910-0801$26.00 c 2005 Pleiades Publishing, Inc.


802

ABUBEKEROV et al.

Table 1. Numerical values of parameters used to synthesize the radial-velocity curves of the optical component in the Roche model P ,day mx , M mv , M e i,deg µ f Teff ,K kx A u


5.0 10 var


Period Mass of the compact object Mass of the optical star Eccentricity Orbital inclination Roche lobe filling coefficient for the optical star Asynchronicity coefficient for the rotation of the optical star Effective temperature of the optical star Gravitational-darkening coefficient Ratio of the X-ray luminosity of the relativistic component to the bolometric luminosity of the optical component Lx /Lv Reprocessing coefficient for the X-ray radiation Limb-darkening coefficient

0.0 30, 60, 90 1.0 1.0 var 0.25




0.1 0.5 0.3

Parameter of the X-ray binary was varied during the model computations. When the mass of the optical star was mv = 1 M , the gravitational-darkening coefficient was taken to be = 0.08.

effects and the conditions for visibility of the areas by an observer. This yields the total radiation flux from the star in the direction toward the observer in the continuum and a rotationally broadened spectral-line profile that can be used to derive the radial velocity of the star. Since the observed radial-velocity curves for OB stars have been derived primarily from hydrogen Balmer absorption lines, we synthesized theoretical radial-velocity curves for the H line of the optical star. The radial velocity at a given orbital phase was calculated using the mean wavelength at the 1/3, 1/2, and 2/3 levels of the maximum depth of the integrated absorption-line profile. In addition, in order to estimate the uncertainties in the modeling, we carried out computations for the same situations using our previous algorithm for synthesizing theoretical radial-velocity curves, used earlier to analyze the radial-velocity curves of OB supergiants in X-ray binaries with neutron stars [7] and in the Cyg X-1 system [4]. This algorithm was proposed by Antokhina and Cherepashchuk [8] in 1994. The main difference from our more modern algorithm [5, 6] is that the local profile of an area element is found, not by constructing a model atmosphere and calculating the intensity of the outgoing radiation in the line and continuum (taking into account reprocessing of the external X-ray radiation), but instead using

computed and tabulated Kurucz Balmer absorptionline profiles for various effective temperatures Teff and gravitational accelerations g. In addition, the effect of heating of the stellar atmosphere by X-ray radiation from its companion was taken into account only in a simple way, by adding the outgoing and incident flux units without taking into account radiative transfer in the stellar atmosphere. In addition to its simplified treatment of the reflection effect, this method for calculating the hydrogen absorption-line profiles is not entirely correct, since the tables of Kurucz [9] present the theoretical line profiles in relative flux units and not intensities. However, since we are using the theoretical line profiles to derive radial velocities, and not for comparisons with observed spectral lines, we consider this approximation to be acceptable. In addition, calculating the theoretical radial-velocity curves using the algorithm described based on the tables of Kurucz requires comparatively little computer time. The use of the more modern algorithm [5, 6], which calculates a model atmosphere for each local area element, requires appreciably more computer time and became possible only relatively recently with the appearance of computers with processing rates of 1 GHz or higher. One drawback of the new algorithm is the absence of a contribution to the integrated absorption-line profile of the optical star from the local profiles of areas near the limb of the stellar disk (the requirement of
ASTRONOMY REPORTS Vol. 49 No. 10 2005


OPTICAL STAR IN AN X-RAY BINARY Table 2. Mass and corresponding effective temperature of the optical star Te ,K ff
eff eff

803

mv , M

T

,K

T

,K

Range of atmospheric parameter values on the stellar surface in the Roche model log g Teff ,K 2800­5930 6570­12 820 9190­17 060 14 020­24 540 18 390­31 010

S, %

1 5 10 20 30
eff eff

­ 14 000 17 000 26 000 29 000

5500 12 000 15 000 23 000 28 000

5500 12 000 16 000 23 000 29 000

2.02­3.15 2.02­3.16 2.13­3.22 2.23­3.25 2.28­3.26

89 93 94 94 96
eff

Note: T and T are the temperatures of the optical star according to the mass­luminosity relations of [11] and [12]; T is the temperature of the optical star used to compute the radial-velocity curves in the current study; and S the fraction of the area of the stellar surface where the temperature differs from the mean effective temperature Te by no more than 10%. ff

the boundary conditions [5, 6]). At the same time, for real optical stars, the local profiles of areas near the disk limb are nonzero and will make a significant contribution to an integrated absorption-line profile. Note that the old algorithm was free of this effect of "zeroing" the contributions of the local profiles for areas near the disk limb. We used both these algorithms to synthesize the theoretical radial-velocity curves for all the model problems considered. The results of these computations using these two methods are compared below. Table 1 presents the parameters of the X-ray binary for which the modeling was carried out. The numerical values of the parameters of the modeled close binary were adopted based on the catalog [10], as being the most characteristic values. We synthesized radial-velocity curves for optical stars with masses mv = 1, 5, 10, 20, 30 M (the remaining parameters of the modeled binary system are presented in Table 1). In order to investigate the dependence of the shape of the radial-velocity curve on the orbital inclination i, we synthesized radialvelocity curves for i = 30 , 60 ,and 90 . When modeling the radial-velocity curves of an optical star with mass mv = 30 M using the tables of Kurucz [9], the local gravitational acceleration gloc and local temperature Tloc at the surface of the optical star fell outside the range of tabulated values, so that there were no tabulated profiles for some of the local areas. The number of such areas was modest (about 10­20% of the total) and they were all located on the "nose" of the filled Roche lobe of the optical component. In this case, the profiles of all the areas were taken to be the same. We used the profile of the H line for the mean effective temperature and the mean gravitational acceleration at the stellar surface
ASTRONOMY REPORTS Vol. 49 No. 10 2005

as the local profile. The accuracy of this approximation was tested in test computation 1 presented below. We should also qualify the situation with regard to the mean effective temperature of the optical star. The mean temperature of the Roche lobe filling star is given by the expression Teff = Tloc dS dS , where the integration is carried out over the entire surface of the tidally deformed star. Table 2 presents estimates of the effective temperature based on the mass­ luminosity relations of Herrero [11] and Straizhis [12]. We can see that the effective temperatures of the star obtained in different ways are appreciably different. Therefore, we carried out a test computation (test computation 2) in order to quantitatively estimate the influence of the effective temperature of the optical star on the theoretical radial-velocity curve. We took the mean effective temperature of the tidally deformed star to be the mean effective temperature of a spherical star with the same volume. The results of the test computations are presented below. Table 2 also presents the range of effective temperatures and gravitational accelerations at the surface of the optical star in the Roche model. Note that the number of local area elements whose temperature Tloc differs from the adopted effective temperature Te by ff more than 10% is very small.

Test Computation 1: Influence of the Local-Profile Approximation on the Theoretical Radial-Velocity Curve
As we indicated above, the test computations check the difference in the radial velocities calculated using the two methods and the tables of Kurucz [9]. We synthesized the radial-velocity curve for an optical


804
Vr, km/s 130 70 10 ­50 ­ 110 ­ 170 , km/s 2 1 0 ­1 ­2 0 (a)

ABUBEKEROV et al.

Vr, km/s 110 70 30 ­10 ­50 ­90 ­ 130

(a)

(b)

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

, km/s (b) 3 2 1 0 ­1 ­2 ­3 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

Fig. 1. Test computation 1. (a) Model radial-velocity curve for an optical star with mv = 10 M and Teff = 15 000 Kfor i = 90 (solid) together with the same curve obtained assuming that the shape of the local profile is constant (dashed). The two radial-velocity curves virtually coincide on the scale shown. These curves were calculated using the tabulated H line profiles (in flux units) of Kurucz. (b) Difference between the absolute values of the radial velocities. See text for more detail.

Fig. 2. Same as Fig. 1 for the parameters of the optical star mv = 20 M and Teff = 23 000 K.

mv = 20 M . The resulting radial-velocity curves are shown in Fig. 2a. In this case, the discrepancy between the radial-velocity curves does not exceed 2 km/s, or 1.6% of the radial-velocity half-amplitude (Fig. 2b).

star with mass mv = 10 M and mean effective surface temperature Teff = 15 000 K (the remaining binary parameters can be found in Table 1); each local area is specified in the LTE approximation in accordance with its local H absorption profile from the tables of Kurucz (solid curve in Fig. 1a). We then constructed the radial-velocity curve for the optical star specifying the shapes of the local profiles for all the area elements to be the same, but taking into account normalization to the continuum over the stellar surface when summing the area profiles. For the constant profile shape, we used the H absorption profile for the mean effective temperature and gravitational acceleration of the optical star. In the case of a star with mv = 10 M , we used the tabulated profile of Kurucz [9] corresponding to Tloc = 15 000 K and log gloc = 3.2. The resulting radial-velocity curve is shown by the dotted curve in Fig. 1a. Figure 1b shows the difference between the absolute values of the radial velocities obtained using the "intrinsic" local profiles and assuming identical local profiles for each area element. We can see that the difference between the curves does not exceed 1.7 km/s, or 1% of the half-amplitude of the radial-velocity curve. An analogous computation was carried out for a close X-ray binary with an optical star with mass

Test Computation 2: Influence of the Effective Temperature of the Optical Star on the Theoretical Radial-Velocity Curve
As we noted above, the effective temperature of the optical star is usually not known exactly. We can see from Table 2 that the observed effective temperatures derived from the mass­luminosity relations of [11, 12] are somewhat different. Therefore, we carried out a test computation to estimate the influence of the effective temperature of the optical star on the theoretical radial-velocity curve. We synthesized a radial-velocity curve for a close binary with an optical star with mass mv = 10 M and orbital inclination i = 90 for Teff = 10 000 K and Teff = 17 000 K (having especially chosen a wide range of variation for the effective temperature). The resulting radial-velocity curves are shown in Fig. 3a (see Table 1 for the remaining binary parameters). Since the discrepancy between the curves is insignificant, Fig. 3b presents the difference between their absolute values, = |Vr (Teff = 17 000 K)|-|Vr (Teff = 10 000 K)|, where |Vr (Teff = 17 000 K)| and |Vr (Teff = 10 000 K)| are the absolute values of the radial velocity of the optical star
ASTRONOMY REPORTS Vol. 49 No. 10 2005


OPTICAL STAR IN AN X-RAY BINARY

805

Vr, km/s 180 120 60 0 ­60 ­ 120 ­ 180 , km/s 4 2 0 ­2 ­4 0

(a)

I 1.0

0.9 (b)
i i i i = = = = 30° 30° 60° 90° = = = = 0 0.35 0.35 0.35

0.8 6437 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

6438

6439

6440

6441 , å

Fig. 3. Test computation 2. (a) Model radial-velocity curves of an optical star with mv = 10 M , i = 90 ,and Teff = 10 000 K (solid) and Teff = 17 000 K (dashed). The curves were calculated using the tabulated H profiles (in flux units) of Kurucz. (b) Difference between the absolute values of the radial velocities = |Vr (Teff = 17 000 K)|- |Vr (Teff = 10 000 K)|. See text for more detail.

° Fig. 4. Theoretical CaI 6439 A absorption profiles without allowance for an instrumental profile (calculated with the new algorithm [5, 6]) at orbital phase = 0.0 for orbital inclination i = 30 (solid) and at orbital phase = 0.35 for i = 30 (dotted), 60 (dash-dotted), and 90 (dashed). The model profiles were obtained assuming LTE and a mass and effective temperature for the optical star mv = 1 M and Teff = 5500 K (seeTable 1for ° the remaining parameters). The CaI 6439 A absorption profiles for phase 0.35 have been corrected for the Doppler shifts due to the orbital motion.

in the Roche model for mean effective surface temperatures of Teff = 17 000 K and Teff = 10 000 K. We can see from Fig. 3b that the maximum discrepancy between these two values occurs at orbital phases 0.35­0.45, and reaches 5 km/s, or 2.6% of the radialvelocity half-amplitude. We can see from the test computation that the uncertainty in the mean effective temperature of the optical star (Teff 5000-7000 K) has an appreciable effect on the shape of the radial-velocity curve. The variations in the radial-velocity curve are different at different orbital phases (Fig. 3b). These variations are maximum at phases 0.35­0.45, where they reach 3% of the radial-velocity half-amplitude. As we noted above, the uncertainty in the effective temperature was artificially increased in the test computation. According to Table 2, the maximum uncertainty in the effective temperature does not exceed 3000 K, so that the corresponding variations in the radialvelocity curve should not exceed 1%. Our computations indicate that the mean effective temperature of the optical star, Teff , should be known as accurately as possible when using radialvelocity curves to determine the parameters of a binary system.
ASTRONOMY REPORTS Vol. 49 No. 10 2005

3. DEPENDENCE OF THE SHAPE OF THE RADIAL-VELOCITY CURVE ON THE COMPONENT-MASS RATIO q AND ORBITAL INCLINATION i We also calculated radial-velocity curves with and without allowing for the influence of the instrumental function on the theoretical integrated profile.

Calculated Radial-Velocity Curves without Account of the Instrumental Function on the Theoretical Absorption Profile
We synthesized radial-velocity curves in the Roche model assuming LTE for optical stars with masses mv = 1, 5, 10, 20, and 30 M (see Table 1 for the remaining parameters of the X-ray binary). For the optical star with mass mv = 1 M , the synthesis ° was carried out for the CaI 6439 A absorption line. The radial-velocity curve syntheses for the stars with masses mv = 5, 10, 20, and 30 M were performed for the H absorption line using the two methods described above (calculating the intensity of the local profile of each area element based on a constructed model atmosphere [5, 6] and based on the tabulated line profiles in flux units of Kurucz [8, 13]). The ° theoretical integrated CaI 6439 Aand H absorption profiles are presented in Figs. 4 and 5. For each value of mv , we synthesized radial-velocity curves for


806

ABUBEKEROV et al.
Vr, km/s 220 140 60 ­20 i = 30° ­ 100 i = 60° ­ 180 i = 90° ­ 260 Vnorm 0.9 0.5 0.1 ­ 0.3 i = 30° i = 60° ­ 0.7 i = 90° ­ 1.1 I 0.012 0.010 0.008 0.006 0.004 0.002 0 ­ 0.002 ­ 0.004 0 0.1 0.2 0.3 0.4 (a)

I 1.0

0.9

0.8 0.7 0.6 4334
i i i i = = = = 30° 30° 60° 90° = = = = 0 0.35 0.35 0.35

(b)

4336

4338

4340

4342

4344

4346 , å

Fig. 5. Same as Fig. 4 for the H absorption line with mv = 20 M and Teff = 23 000 K.

(c)

orbital inclinations i = 30 , i = 60 ,and i = 90 .The resulting theoretical radial-velocity curves for optical stars with masses mv = 1 M , 20 M , and 30 M are presented in Figs. 6a, 7a, and 8a. The half-amplitude of the radial-velocity curve Kv increases with increasing orbital inclination i (Figs. 6a, 7a, 8a). Therefore, each theoretical radialvelocity curve was normalized to its value of Kv , which was taken to be the maximum radial velocity in the orbital-phase interval 0.0­0.5. The normalized radial-velocity curves for optical stars with masses mv = 1, 20, 30 M are presented in Figs. 6b, 7b, 8b. Since the effect of the differences in the shapes of the curves is small (relative to the half-amplitude Kv ), we calculated the differences of the absolute values of the normalized radial velocities = |Vnorm (i = 90 )|-|Vnorm (i = 30 )| and = |Vnorm (i = 90 )|- |Vnorm (i = 60 )|,where |Vnorm (i = 90 )|, |Vnorm (i = 60 )|, and |Vnorm (i = 30 )| are the absolute values of the normalized radial velocities for orbital inclinations of i = 90 , i = 60 , and i = 30 . The maximum differences between the normalized radial-velocity curves obtained using the two methods for optical stars with masses mv = 1, 20, 30 M are reached at orbital phases 0.35­0.45 and 0.55­0.65 (Tables 3 and 4). We denote I and F to be the maximum variations in the shape of the radialvelocity curve for variations in the orbital inclination obtained using the new algorithm of [5, 6] and the theoretical line profiles of Kurucz [8, 13], respectively. The difference in the absolute values of the normalized radial velocities for optical stars with masses mv = 1, 20, 30 M are presented in Figs. 6c, 7c, 8c. We can see from these figures that the maximum difference between the radial-velocity curves obtained

0.5 0.6 0.7 0.8 0.9 1.0

Fig. 6. (a) Theoretical radial-velocity curve of an optical star with mv = 1 M , Teff = 5500 K calculated without allowance for the effect of the instrumental func° tion on the integrated CaI 6439 A absorption profile for orbital inclinations i = 30 (dotted), 60 (dash-dotted), and 90 (solid) (see Table 1 for remaining parameters). The radial-velocity curves were calculated using the new algorithm [5, 6]. (b) The same curves normalized to their half-amplitudes for i = 30 (dotted), 60 (dashdotted), and 90 (solid). The curves nearly coincide on this scale. (c) Difference between the absolute values of the normalized radial velocities I = |Vnorm (i = 90 )|- |Vnorm (i = 60 )| (dashed) and I = |Vnorm (i = 90 )|- |Vnorm (i = 30 )| (solid). The quantity I is given in units of the half-amplitude of the radial-velocity curve Kv (for more detail, see text).

for different values of i occurs at orbital phases 0.35­0.45 and 0.55­0.65. It follows from the computations that the variations in the shape of the radial-velocity curve depend on both the orbital inclination i and the componentmass ratio q = mx /mv . For example, the maximum variations in this curve (in units of the radial-velocity half-amplitude) when the inclination is changed from i = 30 to i = 90 for optical stars with masses mv = 1, 5, 10, 20, and 30 M are F = 0.9, 3.2, 4.6, 5.7, and 8%, respectively. Thus, the variations in the shape of the radialASTRONOMY REPORTS Vol. 49 No. 10 2005


OPTICAL STAR IN AN X-RAY BINARY
Vr, km/s 110 70 30 ­10 ­50 ­90 ­ 130 Vnorm 0.9 0.5 0.1 ­ 0.3 ­ 0.7 ­ 1.1 F 0.07 0.06 0.05 0.04 0.03 0.02 0.01 0 ­ 0.01
i = 30° i = 60° i = 90°

807
(a)

(a)

i = 30° i = 60° i = 90°

(b)

Vr, km/s 120 80 40 0 ­40 ­80 ­ 120 Vnorm 0.9 0.5 0.1 ­ 0.3 ­ 0.7 ­ 1.1 I 0.07 0.06 0.05 0.04 0.03 0.02 0.01 0 ­ 0.01 ­ 0.02

i = 30° i = 60° i = 90°

(b)

i = 30° i = 60° i = 90°

(c)

(c)

0

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

0

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

Fig. 7. () Theoretical radial-velocity curve of an optical star with mv = 20 M , Teff = 23 000 K calculated without allowance for the effect of the instrumental function on the integrated H absorption profile for orbital inclinations i = 30 (dotted), 60 (dash-dotted), and 90 (solid) (see Table 1 for remaining parameters). The radial-velocity curves were calculated using the old algorithm [8, 13] (using the tabulated H absorption profiles of Kurucz in flux units). (b) The same curves normalized to their half-amplitudes for i = 30 (dotted), 60 (dash-dotted), and 90 (solid). (c) Difference between the absolute values of the normalized radial velocities F = |Vnorm (i = 90 )|- |Vnorm (i = 60 )| (dashed) and F = |Vnorm (i = 90 )|- |Vnorm (i = 30 )| (solid). The quantity F is given in units of the halfamplitude of the radial-velocity curve Kv (for more detail see text).

Fig. 8. Same as Fig. 6 for the H absorption line and with mv = 30 M and Teff = 29 000 K.

Computation of Radial-Velocity Curves Including the Effect of the Instrumental Function on the Theoretical Absorption Profiles
As in the previous case, the radial-velocity curve synthesis was carried out for optical stars with masses mv = 1, 5, 10, 20, and 30 M (see Table 1 for the remaining binary parameters). We used a Gaussian profile for the instrumental function of the ° spectrograph. The CaI 6439 A line was convolved with an instrumental function with a full width ° at half maximum intensity FWHM = 1 A. When synthesizing the curves for stars with masses mv = 5, 10, 20, and 30 M , the theoretical H profile was convolved with an instrumental function with ° ° FWHM = 7 A. The convolved theoretical CaI 6439 A and H absorption profiles are presented in Figs. 9 and 10, respectively. We investigated the variations of the radialvelocity curves with variations of the orbital inclination analogous to those described above. As earlier, the maximum changes in the curves occurred at phases 0.35­0.45. The variations in the radialvelocity curve obtained when the inclination was changed from i = 60 to i = 90 are given in Table 5.

velocity curve with changes in the orbital inclination are most clearly expressed in systems with low component-mass ratios q = mx /mv . This is due to the fact that, when q < 1, the center of mass of the binary system lies inside the body of the optical star, so that, during the orbital motion, the part of the stellar surface lying inside the inner Lagrangian point moves in the same direction as the relativistic object. This leads to strong distortion of the integrated absorption-line profile that depends on both i and q .
ASTRONOMY REPORTS Vol. 49 No. 10 2005


808

ABUBEKEROV et al. Table 4. Maximum variation of the shape of the radialvelocity curve in the Roche model for variations of the orbital inclination from 30 to 90 mv ,M 1 5 10 20 30 I , % 0.9 1.5 2.5 4.2 6.5 F , % ­ 3.2 4.6 5.7 8.0

Table 3. Maximum variation of the shape of the radialvelocity curve in the Roche model for variations of the orbital inclination from 60 to 90 mv ,M 1 5 10 20 30 I , % 0.3 0.4 0.6 1.1 1.9 F , % ­ 1.1 1.5 2.3 2.9

Note: The quantities I and F are expressed in units of the half-amplitude Kv . See text for more detail.

Note: The quantities I and F are expressed in units of the half-amplitude Kv . See text for more detail.

The analogous results for the case when the inclination is changed from i = 30 to i = 90 are given in Table 6. A comparison of Tables 3, 4 and 5, 6 shows that the results obtained with and without the effect of the instrumental profile coincide within reasonable errors for the computed integrated profiles. Thus, the variations in the shape of the radialvelocity curve with orbital inclination are not "washed out" by the instrumental function of the spectrograph, leaving intact the possibility of estimating the orbital inclination of a binary system using accurate observations of the radial-velocity curve. Estimating the orbital inclination of a system with a given mass for its optical star reduces to a required accuracy for the radial-velocity observations. For example, for binary systems with optical-star masses mv = 20-30 M (close to filling their Roche lobes), the accuracy in the radial velocities must be better than 7­8% of the halfamplitude of the radial-velocity curve, Kv . For sysI 1.0

tems with optical-star masses mv = 1-10 M (close to filling their Roche lobes), this accuracy must be better than 1­4% of Kv . The orbital inclination of the Cyg X-1 system estimated using an accurate radial-velocity curve is presented in [4]; the accuracy in the observed radial velocities Vr was 3% of the half-amplitude Kv . For an optical-star mass of mv 20 M , this made it possible to place a limit on the orbital inclination of i < 45 based purely on the observed radial-velocity curve. 4. FITTING THE OBSERVED RADIAL-VELOCITY CURVE OF THE Cyg X-1 SYSTEM The fitting of the mean radial-velocity curve for Cyg X-1 presented in [4] was carried out with the earlier algorithm, which used the tabulated H profiles of Kurucz [9] in fluxes for the local profiles of the area elements. However, as we can see from Tables 3­6, the effects of varying the shapes of the radial-velocity curves calculated using the old and new algorithms are somewhat different. For this reason, we decided to fit the accurate mean radial-velocity curve considered in [4] using the new algorithm, in which the integrated H absorption profiles are calculated in intensity units based on model atmospheres for the local area elements [5, 6]. The calculations based on the old algorithm yielded a mass for the compact object in the Cyg X-1 binary of mx = 9.0-13.2 M [4]. No model for the binary system with an orbital inclination i < 45 is consistent with the accurate observational radial-velocity curve [4]. We used the Fisher statistical criterion to test the adequacy of the model description of the observational data. The method used is described in detail in [4], and we do not present it here. As in our earlier work [4], the test of the model description of the observational data wasmadefor the = 5% significance level.
ASTRONOMY REPORTS Vol. 49 No. 10 2005

FWHM = 1 å

0.9
i i i i = = = = 30° 30° 60° 90° = = = = 0 0.35 0.35 0.35

0.8 6437

6438

6439

6440

6441 , å

Fig. 9. Same as Fig. 4 but including the effect of an ° instrumental profile with FWHM = 1 A.


OPTICAL STAR IN AN X-RAY BINARY Table 5. Maximum variation of the shape of the radialvelocity curve in the Roche model for variations of the