New ephemeris of the ADC source 2A 1822371: a stable orbitalperiod derivative over 30 years
L. Burderi^{1}  T. Di Salvo^{2}  A. Riggio^{1,3}  A. Papitto^{1,3}  R. Iaria^{2}  A. D'Aì^{2}  M. T. Menna^{4}
1  Dipartimento di Fisica, Università degli Studi di Cagliari,
SP MonserratoSestu, KM 0.7, Monserrato 09042, Italy
2 
Dipartimento di Scienze Fisiche ed Astronomiche, Università di Palermo,
via Archirafi 36, Palermo 90123, Italy
3 
INAF Osservatorio Astronomico di Cagliari, Loc. Poggio dei Pini,
Strada 54, 09012 Capoterra (CA), Italy
4 
INAF Osservatorio Astronomico di Roma, via Frascati 33, Monteporzio
Catone 00040, Italy
Received 13 July 2009 / Accepted 4 March 2010
Abstract
We report on a timing of the eclipse arrival times of the low mass
Xray binary and Xray pulsar 2A 1822371 performed using all
available observations of the Proportional Counter Array on
board the Rossi Xray Timing Explorer, XMMNewton pn, and Chandra.
These observations span the years from 1996 to 2008.
Combining these eclipse arrival time measurements with those already
available covering the period from 1977 to 1996,
we obtain an orbital solution valid for more than thirty
years.
The time delays calculated with respect to a constant orbital
period model show a clear parabolic trend, implying that the orbital
period in this source constantly increases with time at a rate
s/s. This is 3 orders of
magnitude larger than what is expected from conservative mass transfer
driven by magnetic braking and gravitational radiation.
From the conservation
of the angular momentum of the system we find that to explain
the high and positive value of the orbital period derivative
the mass transfer rate must not be less than 3 times the Eddington
limit for a neutron star, suggesting that the mass transfer has to be
partially nonconservative.
With the hypothesis that the neutron star accretes at the Eddington limit
we find a consistent solution in which at least 70% of the transferred
mass has to be expelled from the system.
Key words: eclipses  ephemerides  pulsars: individual: 2A 1822371  Xrays: binaries  stars: winds, outflows  accretion, accretion disks
1 Introduction
The source 2A 1822371 is a wellknown low mass Xray binary (hereafter LMXB) seen almost edgeon with an inclination angle of (Hellier & Mason 1989). The observed average unabsorbed flux of the source in the 0.1100 keV energy range is erg cm^{2} s^{1} (Iaria et al. 2001). This corresponds to an unabsorbed luminosity of erg s^{1}, adopting a distance of 2.5 kpc (Mason & Cordova 1982). However, it has been noted (Parmar et al. 2000) that the mean ratio of the Xray over optical luminosity, , for 2A 1822371 is about 20, while the average value for LMXBs is about 500 (van Paradijs & McClintock 1994). This would imply an unobscured Xray luminosity as high as erg/s for the assumed distance of 2.5 kpc. The apparent low luminosity of the source has therefore to be ascribed to the high inclination of the system with respect to the line of sight. Indeed, the light curve of 2A 1822371 shows both dips and eclipses of the Xray source by the companion star. The partial nature of the eclipse indicates that the Xray emitting region is extended and that the observed Xrays are scattered in an accretion disk corona (ADC, White et al. 1981). The Xray light curve shows clear signs of orbital modulation with a binary orbital period of 5.57 h. This Xray modulation is probably caused by the obscuration of the ADC by the thick rim of an accretion disk. The orbital period has been measured from eclipse timing to increase gradually (Hellier et al. 1990). Parmar et al. (2000) gave the best ephemeris of this source before this work. In particular they found a significant positive orbital period derivative of s/s.
Table 1: New Xray eclipse times for 2A 1822371.
Jonker & van der Klis (2001) reported on the discovery of 0.59 s Xray pulsations in this source in an RXTE observation performed in 1998. The timing analysis of the pulse arrival times indicates a circular orbit with an eccentricity e < 0.03 (95% c.l.) and an for the neutron star of 1.006(5) lts, implying a mass function of . The comparison between the pulse period measured by RXTE in 1996 and 1998 also indicates that the neutron star in this system is spinning up at a rate of s/s. Jonker & van der Klis (2001) inferred a bolometric Xray luminosity of about erg/s assuming a magnetic field of Gauss. From spectroscopic measurements of the radial velocity curve of the companion, Jonker et al. (2003) derived a lower limit to the mass of the neutron star and to that of the companion star of and , respectively (1, including uncertainties in the inclination), and an accurate estimate of the system inclination angle, .
In this paper we report on the analysis of Xray observations of 2A 1822371 performed from 1996 to 2008 by RXTE, XMMNewton, and Chandra with the aim to derive eclipse arrival times and to improve the orbital ephemeris. We confirm with higher precision and over a much larger time span (about 31 years) the ephemeris found by Parmar et al. (2000). In particular we find that the orbital period derivative has remained constant during the last 30 years. Finally we discuss the implications of a high and positive value of the orbital period derivative on the mass transfer rate and secular evolution of this source.
2 Timing analysis and results
We analysed all available Xray observations of 2A 1822371 performed over the period from 1996 to 2008. In particular we used observations from the PCA on board RXTE performed in 1996 (P10115), 1998 (P30060), 2001 (P50048), 2002 (P70036), 20022003 (P70037), one observation from XMMNewton performed in 2001 (Obs ID: 0111230101 and 0111230201), and two Chandra observations performed in 2000 (Obs ID: 671) and in 2008 (Obs ID: 9076 and 9858), respectively. The arrival times of all events were referred to the solar system barycenter, using as the best estimate for the source coordinates those derived from the 2008 Chandra observations (RA: 18 25 46.81, Dec: 37 06 18.5, uncertainty: 0.6'').
The typical eclipse duration is around 2.2 ks, which corresponds to of the binary orbital period. In order to improve the statistics for the measure of the eclipse epochs and to have the possibility of fitting a complete orbital light curve we decided to perform a folding of these data using the known binary orbital period of the source, after verifying that this folding does not affect the results reported here in any case. Folding the data is not an important issue for the two Chandra observations and the XMM observation, where just one or two consecutive eclipses are observed. But it is important for the RXTE observations, because these are short and sparse, and also because the RXTE observations are continuously interrupted by the Earth occultation at every RXTE orbit (lasting approximately 1.5 h). In this case the folding is required to sample a complete orbital light curve from the source, because this is important for a meaningful fitting of the eclipse. For each of these observations we hence folded the data using the local orbital period as derived from the ephemeris published by Parmar et al. (2000). The 20022003 RXTE dataset (P70036 and P70037) was long enough and we decided to divide it into the following four periods: i) 2002 June 710; ii) 2002 August 218; iii) 2002 September 230; and iv) 2003 August 31September 3. In this way we obtained a total of 10 orbital light curves in which the eclipses were clearly visible (see Table 1 for details on the used observations).
We then fitted these orbital light curves to derive eclipse arrival times with the procedure described below. Because the eclipses are asymmetrical and partial, the exact eclipse centroid times crucially depend on the model adopted to describe their shape as well as the variable continuum they are superimposed on. In order to be conservative in our estimates, we then decided to fit the folded light curves using 10 different models. The first model is that used by Parmar et al. (2000) consisting of a Gaussian and a constant fitted on a phase interval of 0.1 around the eclipse. The second and third models consist again of a Gaussian and a constant plus a linear term (second model) and a linear and quadratic term (third model) fitted on a phase interval of 0.3 around the eclipse. The fourth model is as the third model plus a cubic term fitted on a phase interval of 0.4 around the eclipse. The fifth model consists of a Gaussian and a constant plus a sinusoid of period fixed to the orbital period fitted on the whole 01 phase interval. The models from the sixth one to the tenth one are as the fifth model, plus from 2 to 6 sinusoids with periods fixed to 1/2 up to 1/6 of the orbital period, respectively. The addition of higher harmonic components was required to better describe the overall orbital light curve shape, which differs from a pure sinusoid. We restricted our fitting to the first six harmonics because the addition of higher harmonic components was not statistically significant based on an Ftest.
Thus we obtained 10 eclipse arrival times (each corresponding to one of the models described above) for each orbital light curve. The final eclipse arrival time for each orbital light curve was chosen to be the average of these 10 values, and the associated uncertainty was chosen to be half of the maximum range spanned by these values (1 error included). The uncertainty derived in this way fully takes into account significant discrepancies among the different eclipse arrival times found with a particular model to describe the eclipse and the orbital modulation. The obtained values of the eclipse epochs for each of the 10 orbital light curves and the relative uncertainties are reported in Table 1.
Figure 1: Eclipse time delays with respect to a constant orbital period model plotted vs. the orbital cycle for all the available eclipse time measures spanning the period from 1977 to 2008 together with the bestfit parabola ( top panel), and residuals in units of with respect to the bestfit parabola ( bottom panel). Different symbols indicate different datasets: black squares are from Hellier et al. (1990), red stars are from Parmar et al. (2000), cyan squares are from RXTE data, the magenta dot is from XMMNewton data, and the green triangles are from Chandra data. 

Open with DEXTER 
We then computed the eclipse time delays by subtracting from our measures the eclipse arrival times predicted by a constant orbital period model adopting the orbital period, , and the reference time, , given by Parmar et al. (2000). These time delays were plotted versus the orbital cycle number N. The integer N is the exact number of orbital cycles elapsed since ; i.e., N is the closest integer to under the assumption that that we have verified a posteriori. These results are shown in Fig. 1 together with all delays computed from previously available eclipse times, namely those given by Hellier et al. (1990) and by Parmar et al. (2000), respectively.
These points show a clear parabolic trend that we fitted to the equation
where the correction to the adopted value of the eclipse time, , and to the adopted value of the orbital period, , and the orbital period derivative, , are the fit parameters. We get a very good fit with a . In agreement with previous results, we find a highly significant derivative of the orbital period, which indicates that the orbital period in this system is increasing at a rate of s/s. The bestfit values for the orbital parameters, calculated with the corrections we found from the fit of the parabolic trend of the eclipse epochs with Eq. (1), are shown in Table 2. Note that a similar orbital period derivative was recently found with new measures of optical eclipses by Bayless et al. (2009).
Table 2: Bestfit orbital solution for 2A 1822371 derived from the analysis of the eclipse arrival times from 1977 to 2008.
3 Orbital evolution of 2A 1822371
Apart from mass transfer between the companion and the neutron star,
the orbital evolution of this binary system is expected to be driven
by the emission of gravitational waves and by magnetic braking. Under
the further assumption of conservative mass transfer, orbital
evolution calculations show that the orbital period derivative should
be
=  
(2) 
(see Verbunt 1993; Di Salvo et al. 2008; see also Rappaport et al. 1987), where m_{1} and m are the mass of the primary, M_{1}, and the total mass, M_{1} + M_{2}, in units of respectively, is the mass of the secondary in units of , is the orbital period in units of 5 h (that is appropriate for 2A 1822371 because h), n is the index of the massradius relation of the secondary , and where the term takes into account the effect of the magnetic braking.
In line with Verbunt & Zwaan (1981), Verbunt (1993), and
King (1988; see Tauris 2001, for a review)
we can parametrise this term as
where f is a dimensionless parameter of order of unity: preferred values are f = 0.73 (Skumanich 1972) or f = 1.78 (Smith 1979), and k_{0.277}is the radius of gyration of the star k in units of 0.277, which is the appropriate value for a ZAMS star (Claret & Gimenez 1989). Note that the expression for the MB term given in Verbunt (1993) is recovered from the above adopting f = 1 and k for a ZAMS star. Actually, Tauris (2001) discussed three different expression for the angular momentum losses due to MB, , namely that proposed by Skumanich (Verbunt & Zwaan 1981), that proposed by Stepien (1995), and that proposed by Rappaport et al. (1983), respectively. However with somewhat between them. Because we found that to describe 2A 1822371 a quite large is required, we decided to adopt which resulted in the term adopted above.
The orbital period derivative we measured cannot be explained by a conservative scenario however. A positive orbital period derivative certainly indicates a massradius index n < 1/3; this is indeed a quite general result, which does not depend on the details of the angular momentum losses (see also Eq. (4) below). However, the orbital period derivative we measured, s/s, is about three orders of magnitude larger than what is expected even including the (strongest) MB term! This discrepancy is embarrassingly large suggesting that the conservative evolutionary scenario cannot be applied in this case. A similar conclusion was reached by Bayless et al. (2009), who give an improved ephemeris for this source based on new optical eclipse measures; these authors also note that an extremely high mass accretion rate onto the neutron star, corresponding to about four times the Eddington limit, would be required to explain the observed large orbital period derivative, and conclude that much of the transferred mass must be lost from the system. Below we show how the orbital period derivative we measured can be used to constrain the mass transfer in the system, and how this strongly indicates that a large fraction of the mass which the companion tries to transfer to the neutron star is lost by the system.
The massloss rate from the secondary can be easily calculated as a
function of the orbital period of the system and the measured orbital
period derivative
combining the third Kepler law, which must be always satisfied by the
orbital parameters of the system, with
the condition that in this persistent system the neutron star is accreting mass
through Roche Lobe overflow. This means that the radius of the secondary follows
the evolution of the secondary Roche Lobe radius:
,
where for the secondary
we adopted a massradius relation
and for the radius
of the secondary Roche Lobe we adopted the Paczynski (1971)
approximation
R_{L2} = 2/3^{4/3} [q/(1+q)]^{1/3} a, where a is the
orbital separation, which is valid for small mass ratios,
.
From these conditions it is possible to derive a relation between
the massloss rate from the secondary and the orbital period derivative
where is the secondary mass derivative (negative since the secondary star looses mass) in units of , and is the orbital period derivative in units of 10^{10}. We stress that in Eq. (4) an expression for the angular momentum losses mechanism that drives the evolution of the system (e.g. MB or GR) does not explicitly appear. This is quite relevant because at present there is no general consensus on the absolute strength of the MB term nor on its functional dependence on the other orbital parameters. Indeed the effects of the driving mechanism are implicitly considered through the orbital period derivative, which is a measured quantity in our case.
Equation (4) can be inverted to derive the masstransfer timescale
On this short timescale the response of the secondary star must be adiabatic. For (see below) the envelope is convective and the appropriate index is n = 1/3, in agreement with the condition n < 1/3 discussed above. With this value we find for 2A 1822371. We note that the Eddington limit (in units of ) for accretion onto a neutron star is , where R_{6}(m_{1}) is the neutron star radius in units of 10^{6} cm, which slightly depends on the neutron star mass once an equation of state (EoS) for the ultradense matter is adopted.
Thus, adopting (Jonker et al. 2003) and n = 1/3 in Eq. (4), we have to conclude that the secondary mass loss rate in 2A 1822371 is superEddington. We are therefore forced to conclude that the evolution of the system is highly nonconservative.
In order to search for a possible evolutionary scenario for 2A 1822371 we make the assumption that the neutron star is accreting at the maximum possible rate, i.e. the Eddington limit. It has to be noted that the Eddington limit strictly holds for a spherical geometry, and may not be a constraint for highly magnetised neutron stars for which the accreting matter is channeled onto the magnetic polar caps and the geometry of the matter distribution over the Alfvén surface may not be symmetric (see e.g. Basko & Sunyaev 1976). However, our assumption is justified because the luminosity function for highly magnetized neutron stars (usually found in High Mass Xray Binaries) does not disagree with this assumption (see e.g. Grimm et al. 2002). In particular, no highly magnetized neutron star is known to accrete at a rate much higher than the Eddington limit, and the most luminous high mass Xray binaries containing a neutron star in our Galaxy reach luminosities of the order of the Eddington limit. Moreover, the extrapolated Xray luminosity of 2A 1822371 does not indicate an extremely high Xray luminosity. Hence we do not have any evidence that the limiting mass accretion rate in this source is very different from the Eddington limit.
This results in the following condition:
where is the fraction of the mass lost by the secondary, which is accreted by the neutron star, namely , where and are the mass derivatives of the primary and the secondary, respectively. We consider two EoS, namely the moderately soft FPS and the stiffer L (Cook et al. 1994), which give the relation R_{6}(m_{1}) for the neutron star radius. With this and considering the mass function of the system derived from the timing analysis of the neutron star spin (Jonker & van der Klis 2001), Eqs. (4) and (6) can be solved to derive and for any value of m_{1} between and , which depends on the particular EoS adopted, which is 1.8 and 2.5for EoS FPS and L, respectively. In Fig. 2 , in units of the Eddington mass transfer rate, and are plotted for the appropriate range of neutron star masses for the FPS case (for the L EoS the value of is 30% higher, while in units of the Eddington limit is 30% lower).
We now consider Eq. (3) of Di Salvo et al. (2008), which expresses
the conservation of the angular momentum of the system
where, in this case,
represents all the possible losses of angular momentum from the system caused by MB and GR, where is given by Eq. (3) and
takes into account the effects of angular momentum losses because of mass loss from the system. is the specific angular momentum of the mass leaving the system, , in units of the specific angular momentum of the secondary, that is: , where r_{2} is the distance of the secondary star from the center of mass of the system, and a is the orbital separation.
Figure 2: Secondary mass loss rate in units of the Eddington limit for the FPS EoS ( top panel, solid line), fraction of the mass lost by the secondary star which is accreted onto the neutron star ( top panel, dashed line), and specific angular momentum of the mass leaving the system, , in units of the specific angular momentum of the secondary computed for two different prescription of the magnetic braking strength, namely f = 0.73 (Skumanich 1972) and f = 1.78 (Smith 1979) ( bottom panel). The vertical line at represents the lower limit of the neutron star mass derived by Jonker et al. (2003). 

Open with DEXTER 
Adopting the two values of f discussed above, namely f = 0.73 (Skumanich 1972) or f = 1.78 (Smith 1979), and k_{0.277}=1, Eqs. (3), (7), (8), and (9) can be solved to derive as a function of m_{1}. In Fig. 2 is plotted for the appropriate range of neutron star masses for the FPS case (for the L EoS the values of are 7% higher).
The values of we obtain are in between the specific angular momentum at the inner Lagrangian point, for 2A 1822371, and the specific angular momentum of the secondary, , and actually quite close to . This is expected if the mass lost by the secondary star is blown away because of the radiation pressure exerted by the Eddington luminosity generated by the accretion onto the neutron star.
For both the adopted EoS and all the possible values of the neutron star mass,
the values of
are in the range 0.130.29, which means that the
mass transfer in 2A 1822371 is not conservative, at least, at 70% level, which,
as we already noted, is true independently of any assumption on the
particular angular momentum losses.
Interestingly, this is the key that opens the possibility
of constructing a consistent secular evolution for this system. Indeed the
contact condition,
,
can be solved to
derive a theoretical prediction for the massloss rate once a prescription
is given for the possible losses of angular momentum from the system caused
by MB and GR as in Eq. (8):
where
The function is very sensitive to the scenario adopted: for a conservative scenario () while for in the range 0.130.22 (which is appropriate for the FPS EoS and all the possible values of the neutron star mass) . This means that the term acting to shrink the secondary Roche Lobe  because of the extra angular momentum losses caused by the mass expelled from the system  determines an amplification of the massloss rate through the function with respect to a fully conservative case. This amplified massloss rate in turn determines the high value of the orbital period derivative observed in this system. Consequential angular momentum losses (CAML), i.e. angular momentum losses that are themselves the result of mass transfer, have been proposed in the context of cataclysmic variables (CVs) evolution (see e.g. King & Kolb 1995).
Inserting the values determined in this paper for 2A 1822371 in Eq. (5) we find . This means that the system as it is observed now will probably end on this timescale, possibly with the tidal disruption of the companion star. Indeed, King & Kolb (1995) argued that mass transfer could be unstable when CAML are present. This timescale is extremely short, which indicates that it is possible that some short orbital periodLMXBs can last much shorter than what was previously thought. This evolutionary phase, characterised by a superEddington mass transfer rate, may be a common phase in the evolution of LMXBs, albeit shortliving. Because this phase should not last more than a few million years, there may be very few observed systems in this phase (e.g. the socalled Zsources, which are persistently bright LMXBs). This could have profound implications for the estimate of the actual number of LMXBs produced in the Galaxy as inferred from the observed ones, and also for the predicted number of millisecond binary pulsar. We note that this would help to bring the number of LMXBs in line with the estimated number of millisecond binary pulsars. But a detailed analysis of this delicate and longstanding problem needs a dedicated study of this almost unstable phase of the orbital evolution, which is beyond the scope of this paper and will be discussed in a forthcoming paper.
AcknowledgementsThis work is supported by the Italian Space Agency, ASIINAF I/088/06/0 contract for High Energy Astrophysics. We thank the referee for useful comments on the manuscript.
References
 Basko, M. M., & Sunyaev, R. A. 1976, MNRAS, 175, 395 [NASA ADS] [CrossRef] [Google Scholar]
 Bayless, A. J., Robinson, E. L., Hynes, R. I., Ashcraft, T. A., & Cornell, M. E. 2010, ApJ, 709, 251 [NASA ADS] [CrossRef] [Google Scholar]
 Claret, A., & Gimenez, A. 1989, A&AS, 81, 37 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Cook, G. B., Shapiro, S. L., & Teukolsky, S. A. 1994, ApJ, 424, 823 [NASA ADS] [CrossRef] [Google Scholar]
 Di Salvo, T., Burderi, L., Riggio, A., Papitto, A., & Menna, M. T. 2008, MNRAS, 389, 1851 [NASA ADS] [CrossRef] [Google Scholar]
 Grimm, H. J., Gilfanov, M., & Sunyaev, R. 2002, A&A, 391, 923 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hellier, C., & Mason, K. O. 1989, MNRAS, 239, 715 [NASA ADS] [CrossRef] [Google Scholar]
 Hellier, C., Mason, K. O., Smale, A. P., & Kilkenny, D. 1990, MNRAS, 244, 39 [NASA ADS] [Google Scholar]
 Iaria, R., Di Salvo, T., Burderi, L., & Robba, N. R. 2001, ApJ, 557, 24 [NASA ADS] [CrossRef] [Google Scholar]
 Jonker, P. G., & van der Klis, M. 2001, ApJ, 553, L43 [NASA ADS] [CrossRef] [Google Scholar]
 Jonker, P. G., van der Klis, M., & Groot, P. 2003, MNRAS, 339, 663 [NASA ADS] [CrossRef] [Google Scholar]
 King, A. R. 1988, QJRAS, 29, 1 [NASA ADS] [Google Scholar]
 King, A. R., & Kolb, U. 1995, ApJ, 439, 330 [NASA ADS] [CrossRef] [Google Scholar]
 Mason, K. O., & Cordova, F. A. 1982, ApJ, 262, 253 [NASA ADS] [CrossRef] [Google Scholar]
 Paczynski, B. 1971, ARA&A, 9, 183 [NASA ADS] [CrossRef] [Google Scholar]
 Parmar, A. N., Oosterbroek, T., Del Sordo, S., et al. 2000, A&A, 356, 175 [NASA ADS] [Google Scholar]
 Rappaport, S., Ma, C. P., Joss, P. C., & Nelson, L. A. 1987, ApJ, 322, 842 [NASA ADS] [CrossRef] [Google Scholar]
 Rappaport, S., Verbunt, F., & Joss, P. C. 1983, ApJ, 275, 713 [NASA ADS] [CrossRef] [Google Scholar]
 Skumanich, A. 1972, ApJ, 171, 565 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, M. A. 1979, PASP, 91, 737 [NASA ADS] [CrossRef] [Google Scholar]
 Stepien, K. 1995, MNRAS, 274, 1019 [NASA ADS] [Google Scholar]
 Tauris, T. M. 2001, ASP Conf. Ser., 229, 145 [NASA ADS] [Google Scholar]
 van Paradijs, J., & McClintock, J. E. 1994, A&A, 290, 133 [NASA ADS] [Google Scholar]
 Verbunt, F. 1993, ARA&A, 31, 93 [NASA ADS] [CrossRef] [Google Scholar]
 Verbunt, F., & Zwaan, C. 1981, A&A, 100, L7 [NASA ADS] [Google Scholar]
 White, N. E., Becker, R. H., Boldt, E. A., et al. 1981, ApJ, 247, 994 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Table 1: New Xray eclipse times for 2A 1822371.
Table 2: Bestfit orbital solution for 2A 1822371 derived from the analysis of the eclipse arrival times from 1977 to 2008.
All Figures
Figure 1: Eclipse time delays with respect to a constant orbital period model plotted vs. the orbital cycle for all the available eclipse time measures spanning the period from 1977 to 2008 together with the bestfit parabola ( top panel), and residuals in units of with respect to the bestfit parabola ( bottom panel). Different symbols indicate different datasets: black squares are from Hellier et al. (1990), red stars are from Parmar et al. (2000), cyan squares are from RXTE data, the magenta dot is from XMMNewton data, and the green triangles are from Chandra data. 

Open with DEXTER  
In the text 
Figure 2: Secondary mass loss rate in units of the Eddington limit for the FPS EoS ( top panel, solid line), fraction of the mass lost by the secondary star which is accreted onto the neutron star ( top panel, dashed line), and specific angular momentum of the mass leaving the system, , in units of the specific angular momentum of the secondary computed for two different prescription of the magnetic braking strength, namely f = 0.73 (Skumanich 1972) and f = 1.78 (Smith 1979) ( bottom panel). The vertical line at represents the lower limit of the neutron star mass derived by Jonker et al. (2003). 

Open with DEXTER  
In the text 
Copyright ESO 2010