# Frame dragging and other precessional effects in black hole-pulsar binaries

## 1 Introduction

So far the best arguments for the existence of stellar mass black holes (BHs) are based on dynamical mass estimates in X-ray binaries. The measurement of absorption-line velocities of the secondary star allows us to determine a lower limit to the mass of the compact companion. If the mass of the companion exceeds the calculated maximum mass of a neutron star ($`3M_{}`$) we call it a BH candidate (see Wijers (1996) for a list of BH candidates). But to argue convincingly that the companion is indeed a BH one must rule out possible alternatives to stellar mass BHs. In case of cold, non-rotating neutron stars (NSs) the permissible maximum mass depends on a rather complicated equation of state (EOS) for condensed matter at and above nuclear density. Since our knowledge of the EOS for matter above nuclear density is yet imperfect the maximum masses of NSs range widely, from $`1.5M_{}`$ to about $`2.5M_{}`$ (Friedman, Ipser, & Parker ; Komatsu, Erguchi, & Hachisu (1989); Lattimer et al. (1990); Weber & Glendenning (1992)). This is well below the dynamical mass estimations for many of the suspected BH binaries. However, if we are very conservative and assuming complete ignorance about the EOS above a density of $`10^{14}`$ g cm<sup>-3</sup>, dropping the causality condition (i.e. allowing a dispersive medium), and allowing rotation, then the maximum mass goes up to about $`14M_{}`$ (Friedman & Ipser (1987)). A compact object with a maximum mass of $`14M_{}`$ can explain all the dynamical mass estimations in BH candidates. There are other ways of exceeding the dynamical mass limits, such as rather exotic compact stars like Q-stars (Bahcall, Lynn, & Selipsky ; Miller, Shahbaz, & Nolan (1998)) or boson stars made from a strongly self-interacting scalar field (Colpi, Shapiro, & Wasserman (1986)), the latter having in principle no maximum-mass limit. Finally, one might even think of abandoning general relativity as the correct theory for the description of the strong gravitational fields inside a neutron star (although pulsar timing tests of the strong equivalence principle support the validity of general relativity in these strong field regime (Damour & Schäfer (1991); Wex (1997)).

Recent advances in the theory of accretion physics (e.g., the advection-dominated accretion flow model) lead to the conclusion that the matter accreting onto the compact companion in some of the BH candidates does not hit a hard surface. This may indirectly indicate the presence of an event horizon through which matter and radiation can fall in but from which nothing escapes, a fundamental property of BHs (Narayan, Garcia, & McClintock (1997)). At present the concept of an event horizon in these BH candidates has rather the character of a “most natural explanation” than of an inevitable necessity. Certainly the absence of a hard surface does not prove the presence of a BH as predicted by Einstein’s theory of gravity, i.e. a Kerr BH.

Astrophysical BHs are expected to rotate. The spin of the BH gives rise to a so called gravitomagnetic field (Thorne & Hartle (1985)) which causes the relativistic dragging of inertial frames in the vicinity of the BH. The dragging of inertial frames has two important consequences for an accretion disc around the BH. First, the radius of the last (marginally) stable orbit of the accretion disk is a function of the BH spin (Bardeen, Press, & Teukolsky (1972)). Secondly, any deviation of particle motion from the equatorial plane will cause a precession of the particle orbit (Lense & Thirring (1918); Bardeen & Petterson (1975); Shakura (1987)). In principle, these effects allow the determination of the spin of the accreting BH using observations of low-energy X-ray spectra (Zhang, Cui & Chen (1997)) or high frequency quasi-periodic oscillations (Cui, Zhang, & Chen (1998)). While it is still unclear whether the observed quasi-periodic oscillations are a result of the Lense-Thirring precession of inner parts of a warped accretion disk, spectral analysis seems to show the presence of the gravitomagnetic field of the BH (Marković & Lamb (1998); McClintock (1998)). Nevertheless, the arguments are still require certain plausible assumptions (e.g., the accretion disk reaches all the way down to the last stable orbit, etc.) and the precision with which the spin is determined is rather poor, even in systems where all the input parameters are comparably well known.

The discovery (Hulse & Taylor (1975)) and continuous observation of binary pulsars (i.e. pulsars in orbit with a gravitationally bound companion) has opened an entirely new testing ground for gravity theories. So far binary pulsars provide the only laboratories where one can probe the gravitational radiation properties and strong field aspects of relativistic gravity (Taylor & Weisberg (1989); Damour & Taylor (1992); Damour & Esposito-Farèse 1996a ; Damour & Esposito-Farèse 1996b ). A useful review on pulsar timing and relativistic gravity can be found in Taylor (1993), where the reader can find the basic ideas of this subject.

A pulsar orbiting a BH in a close orbit would certainly be a high precision laboratory for BH physics (Paczynski & Trimble (1979)). After a short observation period it would be possible to measure the relativistic advance of periastron, $`\dot{\omega }`$, and thus obtain the total mass of the binary system to high precision. At this stage the mass of the BH would be fairly well constrained since the mass of the pulsar is unlikely to exceed $`2M_{}`$. The measurement of a second post-Keplerian (PK) parameter would then yield the mass of the pulsar, $`m_p`$ ($`_pM_{}`$), and the mass of the companion, $`m_c`$ ($`_cM_{}`$). It was even suggested (Narayan, Piran, & Shemi (1991), Laguna & Wolszczan (1997)) that a pulsar in a close, nearly edge-on orbit, with a spinning BH might be used to probe the rotation of the BH via the timing delay produced by the dragging of inertial frames caused by the spin of the BH (frame-dragging propagation effect). Unfortunately so far none of the binary pulsars seems to orbit a stellar mass BH of typically $`10M_{}`$ (see Fig. 1).

However, it is possible that pulsars in close orbit with BHs have already been missed in previous pulsar surveys due to the Doppler smearing of the signal during the integration. Whilst it is true that this effect can be largely removed by the application of an “acceleration code” (Johnston & Kulkarni (1991)), only a number of globular-cluster searches have utilized this technique due to previous limitations in computing power (Anderson (1992)). The continual improvement and availability of more powerful computers, however, means that many more surveys are now able to incorporate acceleration searches (Lorimer (1998)).

In this paper we discuss what one can learn about a stellar mass BH by timing its pulsar companion. We focus, in particular, on the observational determination of the BH’s rotation (angular velocity and orientation) and assume throughout the paper that in addition to the five Keplerian parameters, at least two PK parameters are measured with high precision, so that we know $`_p`$, $`_c`$, and therefore $`\mathrm{sin}i`$, i.e. the angle of orbital inclination, $`i`$, up to the ambiguity $`i\pi i`$. We begin with a brief introduction to the definition of the spin and quadrupole moments of a rotating body in general relativity and relations between them (§2). In §3 we give an analytic treatment of the frame-dragging propagation effect and show that it is practically impossible to observe this effect in stellar mass BH-pulsar binaries. In §4 we concentrate on the orbital dynamics of a pulsar orbiting a rotating BH and give a detailed investigation of the precession of the pulsar orbit caused by the spin (relativistic spin-orbit coupling) and quadrupole moment (classical spin-orbit coupling) of the BH. In §5 we calculate the secular changes in observable quantities caused by the relativistic spin-orbit coupling. We show that, in principle, one can extract the spin of the BH from nonlinear-in-time changes of the observables which can be approximated by polynomials in time with sufficient accuracy. We discuss additional effects that can cause secular changes in the observed parameters reducing the measurement accuracy for the BH spin. In §6 we outline the prospects of BH spin determination once a BH-pulsar binary is discovered. In §7 we present a method of extracting the quadrupole moment of the rotating BH from the orbital dynamics of the pulsar and evaluate its actual measurability. We summarize in §8.

## 2 Spin and quadrupole moment of a compact body

The external metric of a stationary, axially symmetric body can be written in the standard form (Bardeen & Wagoner (1971))

$$ds^2=e^{2\nu }c^2dt^2+R^2\mathrm{sin}^2\theta B^2e^{2\nu }(d\varphi \omega dt)^2+e^{2\alpha }(dR^2+R^2d\theta ^2),$$
(1)   
where $`c`$ is the speed of light; $`R`$, $`\varphi `$, $`\theta `$ are ordinary spherical coordinates, and the potentials $`\nu `$, $`B`$, $`\omega `$, and $`\alpha `$ are functions of $`R`$ and $`\theta `$. Butterworth and Ipser (1976) have calculated the asymptotic behavior of these potentials for large $`R`$. In particular, the asymptotic behavior of the potentials $`\nu `$ and $`\omega `$ is (using $`cG1`$)

$`\nu `$ $`=`$ $`{\displaystyle \frac{M}{R}}{\displaystyle \frac{1}{12}}\left({\displaystyle \frac{M}{R}}\right)^3q\left({\displaystyle \frac{M}{R}}\right)^3P_2(\mathrm{cos}\theta )+O\left({\displaystyle \frac{1}{R^4}}\right),`$ (2)   
$`\omega `$ $`=`$ $`2\chi \left({\displaystyle \frac{M}{R}}\right)^3+O\left({\displaystyle \frac{1}{R^4}}\right),`$ (3)   
(Butterworth & Ipser (1976), Laarakkers & Poisson (1997)); $`M`$ is the body’s mass and $`\chi `$ is a dimensionless measure of the body’s angular momentum (spin), $`S`$:

$$\chi \frac{c}{G}\frac{S}{M^2}.$$
(4)   
The dimensionless parameter $`q`$ in equation (2) is related to the quadrupole moment, $`Q`$, of the rotating body:

$$q=\frac{c^4}{G^2}\frac{Q}{M^3}.$$
(5)   
In the Newtonian limit one has

$$Q=\varrho (R^{},\theta ^{})R^2P_2(\mathrm{cos}\theta ^{})𝑑V^{},$$
(6)   
where $`\varrho `$ is the mass density of matter inside the star. The Newtonian theory of self-gravitating, rotating bodies predicts a certain relationship between the angular velocity of the body’s rotation and its quadrupole moment due to the oblateness in the mass distribution (Chandrasekhar (1969)). For the same reason, one now may expect a complicated relation between the rotational parameter of a star, $`\chi `$, and its quadrupole moment, $`q`$, in relativistic gravity theories. For a rotating Kerr BH in general relativity $`\chi 1`$ for if one had $`\chi >1`$ it would give an unacceptable naked singularity (see Hawking & Ellis (1973)). Furthermore, for a Kerr BH the quadrupole parameter, $`q`$, is uniquely determined by the rotational parameter $`\chi `$ (Thorne (1980), Thorne, Price & Macdonald (1986)):

$$q=\chi ^2.$$
(7)   
Similarly for solid extended bodies, Laarakkers and Poisson (1997) have computed the quadrupole moment of rotating neutron stars in the interval between 1.0 and 1.8 solar masses using four different EOS. They have found the maximum values for $`\chi `$ lying in between 0.62 (softest EOS) and 0.73 (stiffest EOS). For the quadrupole moment of the mass distribution, $`q`$, they have derived a dependence on $`\chi `$ which is well reproduced by the relation

$$qC\chi ^2,$$
(8)   
where the constant $`C`$ takes a value between 2.0 and 12.1 depending on the mass of the NS and the EOS.

For spinning boson stars with large self-interaction the maximum value for the parameter $`\chi `$ seems to be between 3 and 4, and the relation between $`\chi `$ and $`q`$ is rather complicated once $`\chi `$ exceeds 0.2 (Ryan (1997)). The following inequality, however, always holds:

$$q10\chi ^2.$$
(9)   
Comparing Eqs. (7) – (9) one sees that simultaneous measurements of the mass, $`M`$, spin parameter, $`\chi `$, and the quadrupole moment, $`q`$, of the pulsar’s compact companion can lead to a unequivocal identification of its physical nature.<sup>1</sup><sup>1</sup>1We are not aware of any calculations done for Q-stars. Such measurements seem to be a straightforward way to observationally verify the existence of a rotating Kerr BH.

## 3 The frame-dragging propagation effect in highly inclined binary systems

In Newtonian gravity the rotation of a star contributes only indirectly to its gravitational field through the rotationally induced oblateness of the mass distribution. The centrifugal flattening of the surfaces of equal mass density caused by the rotation of the star gives rise to a gravitational quadrupole field. In general relativity, however, not only mass itself contributes to the gravitational field. The rotation of a star gives rise to a mass current which contributes in a specific way to its gravitational field. The rotation is the source of the gravitomagnetic field of the rotating body (Thorne, Price & Macdonald (1986)). This gravitomagnetic field influences the motion of test particles and the propagation of light in the vicinity of the rotating body (see, for instance Ciufolini & Wheeler (1995), Mashhoon (1997) and references therein). It is well-known that BHs, although purely geometrical in their nature (i.e. a solution of Einstein’s vacuum-field equations), are physical objects and as such also can carry angular momentum which gives rise to a gravitomagnetic field. In fact, the strength and configuration of the gravitomagnetic field of the BH is used to define its rotation. The rotation of the BH will influence the propagation of photons in a specific way (Bardeen (1973)). For pulsar-timing experiments its influence on the propagation time of radio signals are of particular interest. We call this contribution to the propagation time of the pulsar signals the frame-dragging propagation effect.

Recently, Laguna and Wolszczan (1997) have suggested on the basis of numerical ray-tracing calculations that pulsar-timing experiments could measure the gravitomagnetic field of a rotating BH companion by its influence on the propagation of the pulsar signals emitted from different points of the pulsar’s orbit while the pulsar is passing through the point of the superior conjunction with the companion. In this section we give a more detailed investigation of this frame-dragging propagation effect based on analytic calculations and show the problems preventing its direct measurement.

For a discussion of the effect in question it is convenient to use the metric of Kerr spacetime, written in terms of Boyer-Lindquist time $`t`$ and asymptotically Cartesian spatial coordinates $`X^j`$. One has (Thorne, Price & Macdonald (1986)):

$$ds^2=\alpha ^2c^2dt^2+g_{jk}(dX^j+\beta ^jdt)(dX^k+\beta ^kdt),$$
(10)   
where at large distances from the BH one finds for the lapse, $`\alpha `$, and shift, $`\beta ^i`$, functions the following expressions (Thorne, Price & Macdonald (1986)):

$$\alpha ^2=1+𝒪\left(\frac{1}{R}\right),𝜷=\frac{2G}{c^2}\frac{𝐒\times 𝐗}{R^3}+𝒪\left(\frac{1}{R^3}\right),g_{jk}=\delta _{jk}+𝒪\left(\frac{1}{R}\right),$$
(11)   
where $`R=|𝐗|`$ is the distance from the BH to the current point in space. The shift function $`𝜷`$ represents the vector potential of the gravitomagnetic vector field which drags local inertial frames into rotation (frame dragging effect). Since for photons $`ds^2=0`$ the contribution of the frame dragging effect to the propagation time of a photon from the pulsar to the observer is

$$dt_{FD}\frac{𝜷d𝐗}{c^2}=\frac{2G}{c^4}\frac{(𝐒\times 𝐗)d𝐗}{R^3}.$$
(12)   
After substituting into the right-hand side of equation (12) the unperturbed trajectory of the photon, which is a straight line, the equation is easily integrated in a coordinate system where $`𝐗=(x,0,b)`$. The constant $`b`$ is the minimum distance of the unperturbed trajectory of the photon to the BH.<sup>2</sup><sup>2</sup>2Details about the propagation of light rays in the stationary field of relativistic gravitational multipoles can be found in (Kopeikin 1997a ) Consequently, the timing delay caused by the frame-dragging effect assumes the form

$$\mathrm{\Delta }_{FD}=\frac{2G}{c^4}S_yb_{x_0}^{\mathrm{}}\frac{dx}{(x^2+b^2)^{3/2}}=\frac{2G}{c^4}\frac{S_y}{b}\left(1\frac{x_0}{r_0}\right),$$
(13)   
where $`r_0=|𝐗_0|`$ is the distance between the BH and the point of emission of the photon which is located at the pulsar’s orbital position. We have assumed that the observer is far away from the rotating BH; for this reason, the upper limit of the integral in (13) was taken to be infinity. Indeed, from equation (12) it is already clear that most of the contribution to this frame-dragging (FD) propagation effect comes from a region around the BH with a few impact parameters, $`b`$, in radius.

Let $`𝐑`$ denote the vector pointing from the BH to the pulsar and $`𝐊_0`$ denote the unit vector pointing from the observer to the pulsar (see Fig. 2). Then, $`R|r_0|`$, and one has (cf. eq. (43) for $`l=1`$ in Kopeikin 1997a )

$$\mathrm{\Delta }_{FD}=\frac{2G}{c^4}\frac{𝐒(𝐊_0\times 𝐑)}{R(R𝐊_0𝐑)}.$$
(14)   
Using the angles as defined in Fig. 2 we find for the FD propagation effect in a binary pulsar system

$$\mathrm{\Delta }_{FD}=\mathrm{\Lambda }^1\left\{𝒜_{FD}\mathrm{cos}i\mathrm{sin}\left[\omega +A_e(u)\right]+_{FD}\mathrm{cos}\left[\omega +A_e(u)\right]\right\},$$
(15)   
where the function $`\mathrm{\Lambda }`$ is defined by

$$\mathrm{\Lambda }=1e\mathrm{cos}us[\mathrm{sin}\omega (\mathrm{cos}ue)+(1e^2)^{1/2}\mathrm{cos}\omega \mathrm{sin}u],$$
(16)   
the eccentric anomaly angle is

$$A_e(u)=2\mathrm{arctan}\left[\left(\frac{1+e}{1e}\right)^{1/2}\mathrm{tan}\left(\frac{u}{2}\right)\right],$$
(17)   
and

$`𝒜_{FD}`$ $`=`$ $`+4T_{}^{5/3}\left({\displaystyle \frac{2\pi }{P_b}}\right)^{2/3}{\displaystyle \frac{_{}^2}{(_p+_{})^{1/3}}}\chi \mathrm{sin}\lambda _{}\mathrm{cos}\eta _{},`$ (18)   
$`_{FD}`$ $`=`$ $`4T_{}^{5/3}\left({\displaystyle \frac{2\pi }{P_b}}\right)^{2/3}{\displaystyle \frac{_{}^2}{(_p+_{})^{1/3}}}\chi \mathrm{sin}\lambda _{}\mathrm{sin}\eta _{}.`$ (19)   
When $`_p_{}`$ we find the following numerical estimate for the constant factor

$$T_{}^{5/3}\left(\frac{2\pi }{P_b}\right)^{2/3}\frac{_{}^2}{(_p+_{})^{1/3}}(0.0001\mu \mathrm{s})\left(\frac{P_b}{1\mathrm{day}}\right)^{2/3}\left(\frac{_{}}{10}\right)^{5/3}.$$
(20)   
To illustrate the expected strength of the FD propagation effect we assume a pulsar in a circular orbit where $`\eta _{}=0`$ and the inclination angle $`i`$ is close to $`90^{}`$, i.e. $`e=0`$, $`\omega =0`$, $`A_e=u`$, $`s=\mathrm{sin}i1`$, and $`\mathrm{cos}i1`$. Then $`𝒜_{FD}0`$, $`_{FD}=0`$, and

$$\frac{\mathrm{sin}(\omega +A_e)}{\mathrm{\Lambda }}=\frac{\mathrm{sin}u}{1s\mathrm{sin}u}.$$
(21)   
This function has maximum at $`u=\pi /2`$ and minimum at $`u=3\pi /2`$. Thus, for the maximal magnitude of the FD effect we have:

$$\mathrm{max}\left(\mathrm{\Delta }_{FD}\right)\frac{0.0008\mu \mathrm{s}}{|\mathrm{cos}i|}\left(\frac{P_b}{1\mathrm{day}}\right)^{2/3}\left(\frac{_{}}{10}\right)^{5/3}\chi \mathrm{sin}\lambda _{}.$$
(22)   
Hence, for BH binary pulsars with an orbital inclination $`i`$ very close to $`90\mathrm{°}`$, the FD propagation effect will have a measurable influence on the TOAs, provided that the pulsar is a millisecond pulsar where one can expect (with present-day technology) a timing accuracy of better than one microsecond. This fact was already pointed out by Laguna & Wolszczan (1997) using numerical methods to study the propagation of the pulsar signals in the spacetime of the Kerr-BH companion. In particular, for a binary pulsar with an orbital period $`P_b`$ about 1 day and an extreme $`10M_{}`$ Kerr BH companion with $`\chi =1`$ the strength of the FD propagation effect will be one microsecond if $`\mathrm{cos}i0.001`$ or $`|i90\mathrm{°}|0.05\mathrm{°}`$.

A measurement of at least one of the parameters $`𝒜_{FD}`$ or $`_{FD}`$ will give values for either $`\chi \mathrm{sin}\lambda _{}\mathrm{cos}\eta _{}`$ or $`\chi \mathrm{sin}\lambda _{}\mathrm{sin}\eta _{}`$. Since we have a priori no idea about the orientation of the BH spin it is not possible to extract $`\chi `$ from these measurements. But we get a lower limit on $`|\chi |`$, which has to be less than one, otherwise it is not a Kerr BH. In the lucky case of having good values for both parameters, $`𝒜_{FD}`$ and $`_{FD}`$, we can calculate $`\eta _{}`$ and consequently $`\chi \mathrm{sin}\lambda _{}`$ which gives a definite lower limit for the $`\chi `$ of the companion. An upper limit to $`\chi `$ can be set only with a certain probability excluding values of $`\lambda _{}`$ being close to either $`0\mathrm{°}`$ or $`180\mathrm{°}`$.

Unfortunately, the measurement of the FD propagation effect is complicated (and may be impossible) by a competing effect which occurs during the superior conjunction, the bending delay (Doroshenko & Kopeikin (1995)) The delay originates because of the modulation of the pulsar’s rotational phase by the effect of gravitational deflection of light in the field of the pulsar’s companion. The expression for the bending delay in the framework of general relativity is given by the formula (Doroshenko & Kopeikin (1995))

$$\mathrm{\Delta }_B=\mathrm{\Lambda }^1[𝒜_B\mathrm{cos}i\mathrm{sin}(\omega +A_e)+_B\mathrm{cos}(\omega +A_e)],$$
(23)   
where

$`𝒜_B`$ $`=`$ $`+{\displaystyle \frac{T_{}^{2/3}}{\pi \nu _p}}\left({\displaystyle \frac{2\pi }{P_b}}\right)^{2/3}{\displaystyle \frac{_{}}{(_p+_{})^{1/3}}}{\displaystyle \frac{\mathrm{cos}\eta _p}{\mathrm{sin}\lambda _p}},`$ (24)   
$`_B`$ $`=`$ $`{\displaystyle \frac{T_{}^{2/3}}{\pi \nu _p}}\left({\displaystyle \frac{2\pi }{P_b}}\right)^{2/3}{\displaystyle \frac{_{}}{(_p+_{})^{1/3}}}{\displaystyle \frac{\mathrm{sin}\eta _p}{\mathrm{sin}\lambda _p}}.`$ (25)   
Here, $`\nu _p`$ is the rotational frequency of the pulsar, and the angles $`\lambda _p,\eta _p`$ have the same meaning as $`\lambda _{},\eta _{}`$ but define now the orientation of the pulsar’s angular velocity vector in space. For $`_p_{}`$ we find

$$\frac{T_{}^{2/3}}{\pi \nu _p}\left(\frac{2\pi }{P_b}\right)^{2/3}\frac{_{}}{(_p+_{})^{1/3}}(0.0007\mu \mathrm{s})\left(\frac{P}{1\mathrm{ms}}\right)\left(\frac{P_b}{1\mathrm{day}}\right)^{2/3}\left(\frac{_{}}{10}\right)^{2/3},$$
(26)   
where $`P`$ denotes the rotational period of the pulsar.

The bending delay has the same functional dependence on the orbital motion as the FD propagation effect and, for this reason, one can observe only a linear combination of the FD effect and the bending delay. Hence, in timing observations the parameters of the FD effect are not separated from those of the bending delay. For pulsars with rotational periods close to 1 ms the amplitudes of these two effects are comparable if $`_{}20`$. For smaller BH masses and/or longer pulsar periods the bending delay is clearly the dominating effect. Only for BH companions with more than $`100M_{}`$ the FD propagation effect will dominate over the bending delay.

In principle, one may extract $`\eta _p`$ and $`\lambda _p`$ from additional pulse structure analysis (Damour & Taylor (1992)). In this case the bending delay can be predicted and subtracted from the total effect. In practice, pulse structure observations and subsequent analysis are connected with rather large uncertainties in the corresponding angles. Therefore, while using results of such an analysis one can still hope to get a comparably good measurement for $`\lambda _p`$, the angle $`\eta _p`$ affects the polarization pattern of the pulsar only at the negligible small $`x/P_b`$ level (Damour & Taylor (1992)) and, thus, will remain unobservable.

## 4 Spin-orbit coupling and precession of the binary orbit

Eighty years ago Lense & Thirring (1918) pointed out that the gravitomagnetic field of a central rotating body will cause a precession of the orbit of a test particle. In the same way, the rotation of one or both components of a binary system will cause a precession of the binary orbit (Brumberg (1991)). In this and the following sections we will show that the observation of such a precession can lead to the direct determination of the spin of a BH companion.

The typical moment of inertia for a pulsar is of order $`10^{45}\mathrm{cm}^2\mathrm{g}`$ (Arnett & Browns (1977)). Thus, assuming the rotational period of the pulsar to be just 1 ms we obtain a spin of order $`6\times 10^{48}\mathrm{cm}^2\mathrm{g}/\mathrm{s}`$. Using equation (4) we find that the spin of a $`10M_{}`$ extreme Kerr BH ($`\chi =1`$) has a value more than 100 times bigger. Since for even a very soft EOS the rotational period of a pulsar should lie above $``$0.5 ms (Weber & Glendenning (1992)), the spin of the BH will dominate the orbital precession. For this reason, we neglect the spin of the pulsar in our subsequent calculations. The secular precession of the orbit is then given by two vector differential equations of first order (Brumberg (1991))

$$\dot{𝐋}=𝛀_{\mathrm{prec}}\times 𝐋,\dot{𝐀}=𝛀_{\mathrm{prec}}\times 𝐀,$$
(27)   
where $`𝐋`$ is the orbital angular momentum and $`𝐀`$ is the Laplace-Runge-Lenz vector which points to the instantaneous position of periastron of the precessing orbit. The angular velocity vector $`𝛀_{\mathrm{prec}}`$ is the sum of the well-known relativistic (post-Newtonian) periastron advance, $`𝛀_{PN}`$, the gravitomagnetic Lense-Thirring precession, $`𝛀_S`$, caused by the coupling of the orbital angular momentum vector to the spin of the BH, and the classical precession, $`𝛀_Q`$, due to the Newtonian coupling of the orbital angular momentum vector to the quadrupole moment of the rotating BH

$$𝛀_{\mathrm{prec}}=𝛀_{PN}+𝛀_S+𝛀_Q,$$
(28)   
where

$`𝛀_{PN}`$ $`=`$ $`\mathrm{\Omega }_{PN}^{}\widehat{𝐋},`$ (29)   
$`𝛀_S`$ $`=`$ $`\mathrm{\Omega }_S^{}[3(\widehat{𝐋}\widehat{𝐒})\widehat{𝐋}\widehat{𝐒}],`$ (30)   
$`𝛀_Q`$ $`=`$ $`\mathrm{\Omega }_Q^{}[\{5(\widehat{𝐋}\widehat{𝐒})^21\}\widehat{𝐋}2(\widehat{𝐋}\widehat{𝐒})\widehat{𝐒}],`$ (31)   
and a hat on a vector indicates the unit vector in the same direction as that of the vector itself. It is important to stress that $`\mathrm{\Omega }_S^{}`$ and $`\mathrm{\Omega }_Q^{}`$ are not equal to the absolute values of the corresponding vectors. (We have deliberately used the asterisk to avoid such a possible confusion.) One finds (see Barker & O’Connell (1975)):

$`\mathrm{\Omega }_{PN}^{}`$ $`=`$ $`{\displaystyle \frac{T_{}^{2/3}}{1e^2}}\left({\displaystyle \frac{2\pi }{P_b}}\right)^{5/3}3(_p+_{})^{2/3},`$ (32)   
$`\mathrm{\Omega }_S^{}`$ $`=`$ $`{\displaystyle \frac{\chi T_{}}{(1e^2)^{3/2}}}\left({\displaystyle \frac{2\pi }{P_b}}\right)^2{\displaystyle \frac{_{}(3_p+4_{})}{2(_p+_{})}},`$ (33)   
$`\mathrm{\Omega }_Q^{}`$ $`=`$ $`{\displaystyle \frac{qT_{}^{4/3}}{(1e^2)^2}}\left({\displaystyle \frac{2\pi }{P_b}}\right)^{7/3}{\displaystyle \frac{3_{}^2}{4(_p+_{})^{2/3}}}.`$ (34)   
It is worth emphasizing that both vectors $`𝐋`$ and $`𝐒`$ precess around the conserved total angular momentum of the binary system $`𝐉𝐋+𝐒`$, while their absolute values, $`L|𝐋|`$ and $`S|𝐒|`$, are conserved quantities (averaged over one orbital revolution).

The precession of the orbital plane and the longitude of periastron are best described by the angles $`\mathrm{\Phi }`$ and $`\mathrm{\Psi }`$ defined with respect to the invariable plane, i.e. the plane perpendicular to $`𝐉`$ (see Fig. 3 for details). From the identities

$`𝛀_{\mathrm{prec}}\times 𝐋=\dot{𝐋}=\dot{\mathrm{\Phi }}\widehat{𝐉}\times 𝐋,`$ (35)   
$`𝛀_{\mathrm{prec}}\times 𝐀=\dot{𝐀}=(\dot{\mathrm{\Phi }}\widehat{𝐉}+\dot{\mathrm{\Psi }}\widehat{𝐋})\times 𝐀`$ (36)   
one obtains the precession in terms of the angles $`\mathrm{\Phi }`$ and $`\mathrm{\Psi }`$:

$$\dot{\mathrm{\Phi }}=\dot{\mathrm{\Phi }}_S+\dot{\mathrm{\Phi }}_Q=\mathrm{const}.,$$
(37)   
where

$`\dot{\mathrm{\Phi }}_S`$ $`=`$ $`\mathrm{\Omega }_S^{}\left({\displaystyle \frac{\mathrm{sin}\theta }{\mathrm{sin}\theta _J}}\right),`$ (38)   
$`\dot{\mathrm{\Phi }}_Q`$ $`=`$ $`\mathrm{\Omega }_Q^{}\left({\displaystyle \frac{\mathrm{sin}2\theta }{\mathrm{sin}\theta _J}}\right),`$ (39)   
and

$$\dot{\mathrm{\Psi }}=\dot{\mathrm{\Psi }}_{PN}+\dot{\mathrm{\Psi }}_S+\dot{\mathrm{\Psi }}_Q=\mathrm{const}.,$$
(40)   
where

$`\dot{\mathrm{\Psi }}_{PN}`$ $`=`$ $`\mathrm{\Omega }_{PN}^{}=\dot{\omega }_{GR0},`$ (41)   
$`\dot{\mathrm{\Psi }}_S`$ $`=`$ $`\mathrm{\Omega }_S^{}(2\mathrm{cos}\theta +\mathrm{sin}\theta \mathrm{cot}\theta _J),`$ (42)   
$`\dot{\mathrm{\Psi }}_Q`$ $`=`$ $`\mathrm{\Omega }_Q^{}\left(\frac{1}{2}+\frac{3}{2}\mathrm{cos}2\theta +\mathrm{sin}2\theta \mathrm{cot}\theta _J\right).`$ (43)   
In particular, eqs. (39), (43) can be found in Smarr & Blandford (1976) (see also Lai, Bildsten, & Kaspi (1995)).

It is worth noting that the rate of the precession depends on two angles $`\theta `$ and $`\theta _J`$. However, given the spin of the BH, $`S`$, one can express $`\theta _J`$ as a function of $`\theta `$, since $`L\mathrm{sin}\theta _J=S\mathrm{sin}(\theta \theta _J)`$. We find

$`{\displaystyle \frac{1}{\mathrm{sin}\theta _J}}`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{sin}\theta }}\left(1+{\displaystyle \frac{L^2}{S^2}}+2{\displaystyle \frac{L}{S}}\mathrm{cos}\theta \right)^{1/2},`$ (44)   
$`\mathrm{cot}\theta _J`$ $`=`$ $`\mathrm{cot}\theta +{\displaystyle \frac{L}{S\mathrm{sin}\theta }},`$ (45)   
where

$$\frac{S}{L}=\frac{\chi T_{}^{1/3}}{(1e^2)^{1/2}}\left(\frac{2\pi }{P_b}\right)^{1/3}\frac{_{}(_p+_{})^{1/3}}{_p}.$$
(46)   
Hence, the rate of the orbital precession may be expressed as a function of only one angle.

We conclude this section by giving numerical values for the various precessional effects. Assuming a pulsar mass of $`1.4M_{}`$ and a Kerr BH mass $`_{}10`$ allows us to approximate equations (32) to (34) by the expressions

$`\mathrm{\Omega }_{PN}^{}`$ $``$ $`{\displaystyle \frac{1}{1e^2}}\left({\displaystyle \frac{P_b}{1\mathrm{day}}}\right)^{5/3}\left({\displaystyle \frac{_{}}{10}}\right)^{2/3}\mathrm{deg}/\mathrm{yr},`$ (47)   
$`\mathrm{\Omega }_S^{}`$ $``$ $`0.9\times 10^3{\displaystyle \frac{\chi }{(1e^2)^{3/2}}}\left({\displaystyle \frac{P_b}{1\mathrm{day}}}\right)^2\left({\displaystyle \frac{_{}}{10}}\right)\mathrm{deg}/\mathrm{yr},`$ (48)   
$`\mathrm{\Omega }_Q^{}`$ $``$ $`0.5\times 10^6{\displaystyle \frac{\chi ^2}{(1e^2)^2}}\left({\displaystyle \frac{P_b}{1\mathrm{day}}}\right)^{7/3}\left({\displaystyle \frac{_{}}{10}}\right)^{4/3}\mathrm{deg}/\mathrm{yr}.`$ (49)   
Now it is easy to see that the precession in $`\mathrm{\Psi }`$ is clearly dominated by $`\mathrm{\Omega }_{PN}^{}`$ while that in $`\mathrm{\Phi }`$ is dominated by the relativistic spin-orbit coupling $`\mathrm{\Omega }_S^{}`$. The precession $`\mathrm{\Omega }_Q^{}`$ caused by the quadrupole moment of the BH is three orders of magnitude smaller than the relativistic spin-orbit precession and, therefore, will be omitted from the following discussion.

## 5 Spin-orbit coupling and observable quantities

The angles $`\mathrm{\Phi }`$ and $`\mathrm{\Psi }`$, which change linearly in time, are not directly observable in analyzing pulse arrival times. Instead, one can extract from the timing observations the projected semi-major axis of the pulsar orbit, $`x=a_p\mathrm{sin}i/c`$, and the longitude of periastron, $`\omega `$, which are connected with $`\mathrm{\Phi }`$ and $`\mathrm{\Psi }`$ through the trigonometric relationships (see Fig. 3)

$$\mathrm{cos}i=\mathrm{cos}i_J\mathrm{cos}\theta _J\mathrm{sin}\theta _J\mathrm{sin}i_J\mathrm{cos}\mathrm{\Phi },$$
(50)   
and

$`\mathrm{sin}i\mathrm{sin}\omega `$ $`=`$ $`(\mathrm{sin}\theta _J\mathrm{cos}i_J+\mathrm{cos}\theta _J\mathrm{sin}i_J\mathrm{cos}\mathrm{\Phi })\mathrm{sin}\mathrm{\Psi }+\mathrm{sin}i_J\mathrm{sin}\mathrm{\Phi }\mathrm{cos}\mathrm{\Psi },`$ (51)   
$`\mathrm{sin}i\mathrm{cos}\omega `$ $`=`$ $`(\mathrm{sin}\theta _J\mathrm{cos}i_J+\mathrm{cos}\theta _J\mathrm{sin}i_J\mathrm{cos}\mathrm{\Phi })\mathrm{cos}\mathrm{\Psi }\mathrm{sin}i_J\mathrm{sin}\mathrm{\Phi }\mathrm{sin}\mathrm{\Psi }.`$ (52)   
Consequently, due to the non-linear character of these relationships, the linear-in-time precession of $`\mathrm{\Phi }`$ and $`\mathrm{\Psi }`$ will cause a non-linear-in-time precessional evolution of the observed parameters $`x`$ and $`\omega `$. Although the changes in $`\omega `$ and $`x`$ can be expressed in a closed analytic form, a Taylor expansion in powers of time $`tT_0`$, where $`T_0`$ is the time of periastron passage, is more suitable for observational purposes, as will become clear in the next section. Thus, one has:

$`x^{\mathrm{prec}}(t)`$ $`=`$ $`x_0+\dot{x}^{\mathrm{prec}}(tt_0)+{\displaystyle \frac{1}{2}}\ddot{x}^{\mathrm{prec}}(tt_0)^2+\mathrm{},`$ (53)   
$`\omega ^{\mathrm{prec}}(t)`$ $`=`$ $`\omega _0+\dot{\omega }^{\mathrm{prec}}(tt_0)+{\displaystyle \frac{1}{2}}\ddot{\omega }^{\mathrm{prec}}(tt_0)^2+\mathrm{},`$ (54)   
where $`x_0`$, $`\omega _0`$ are initial values of $`x`$, $`\omega `$ at the initial epoch $`T_0`$,

$`\dot{x}^{\mathrm{prec}}`$ $`=`$ $`\dot{x}_S+\dot{x}_Q`$ (55)   
$`\dot{\omega }^{\mathrm{prec}}`$ $`=`$ $`\dot{\omega }_{GR0}+\dot{\omega }_S+\dot{\omega }_Q,`$ (56)   
and subindeces $`S`$ and $`Q`$ mean that the corresponding quantity is caused by the influence of the BH spin or the BH quadrupole moment, respectively. The quantity $`\dot{\omega }_{GR0}`$ describes the standard relativistic periastron advance.

As mentioned previously, the contribution of the quadrupole moment of the BH companion to any observable parameter is extremely small. For this reason, we again omit, in what follows, all quantities with the subindex $`Q`$. Further simplification is achieved if one makes use of the fact that, in general, $`SL`$ and, thus, as a consequence of equation (44), $`\theta _JS\mathrm{sin}\theta /L1`$. Indeed, when $`_{}10`$,

$$\frac{S}{L}0.01\frac{\chi }{(1e^2)^{1/2}}\left(\frac{P_b}{1\mathrm{day}}\right)^{1/3}\left(\frac{_{}}{10}\right)^{4/3}.$$
(57)   
Therefore, even for a binary system with an orbital period $`P_b`$ of 0.5 days, an eccentricity $`e`$ of 0.8, and a 20$`M_{}`$ BH companion $`\theta _J`$ will be smaller than $`3\mathrm{°}`$. Taking this into account equations (50) and (51) assume the form:

$$\mathrm{sin}i=\mathrm{sin}i_J+\theta _J\mathrm{cos}i_J\mathrm{cos}\mathrm{\Phi }+𝒪(\theta _J^2)$$
(58)   
and

$$\omega =\mathrm{\Psi }+\mathrm{\Phi }\theta _J\mathrm{cot}i_J\mathrm{sin}\mathrm{\Phi }+𝒪(\theta _J^2),$$
(59)   
(cf. Wex (1998)). Thus, one obtains, to leading order in $`\theta _J`$,

$`\dot{\omega }_S`$ $``$ $`B_S^{}\chi \mathrm{sin}\theta \mathrm{cos}\mathrm{\Phi }_02A_S^{}\chi \mathrm{cos}\theta ,`$ (60)   
$`{\displaystyle \frac{\dot{x}_S}{x}}`$ $``$ $`B_S^{}\chi \mathrm{sin}\theta \mathrm{sin}\mathrm{\Phi }_0,`$ (61)   
$`\ddot{\omega }_S`$ $``$ $`+C_S^{}\chi \mathrm{sin}\theta \mathrm{sin}\mathrm{\Phi }_0,`$ (62)   
$`{\displaystyle \frac{\ddot{x}_S}{x}}`$ $``$ $`C_S^{}\chi \mathrm{sin}\theta \mathrm{cos}\mathrm{\Phi }_0,`$ (63)   
where $`\mathrm{\Phi }_0`$ is the numerical value of the angle $`\mathrm{\Phi }`$ at the initial epoch $`T_0`$ and the constant coefficients

$$A_S^{}\mathrm{\Omega }_S^{}/\chi ,B_S^{}A_S^{}\mathrm{cot}i,C_S^{}B_S^{}\mathrm{\Omega }_S^{}L/S$$
(64)   
are independent of the angles $`\theta `$, $`\mathrm{\Phi }_0`$, and the spin parameter $`\chi `$. The values for these coefficients are fixed once we know the masses of the binary system.

At the beginning of the timing project one will just be able to see the linear trend in the precession, $`\dot{\omega }`$, and some time later $`\dot{x}`$. But after a while, depending on the compactness of the binary system, one will start to see $`\ddot{\omega }`$ and $`\ddot{x}`$. Since $`\dot{\omega }`$ is clearly dominated by $`\dot{\omega }_{GR0}`$ the first sign for the existence of spin-orbit coupling in the binary system will come from the measurement of $`\dot{x}`$. Under favorable circumstances (see equation (65) and related discussion) this allows us to calculate $`\chi \mathrm{sin}\theta \mathrm{sin}\mathrm{\Phi }_0`$ if we assume that the masses of pulsar and BH and, thus, the angle of the orbital inclination, $`i`$, are known from the measurement of at least two PK parameters. Since $`|\chi \mathrm{sin}\theta \mathrm{sin}\mathrm{\Phi }_0||\chi |`$ we can get a lower limit for the spin of the BH or, if we find $`|\chi |>1`$, a Kerr BH is ruled out. In fact, it is interesting to note that if one assumes that the precession of PSR J0045–7319 (see Kaspi et al. 1996) is caused by a BH companion, the analysis presented here would give $`|\chi |100`$ which is certainly not a Kerr BH.<sup>3</sup><sup>3</sup>3There is clear observational evidence that the companion of PSR J0045–7319 is a B star (Bell et al. 1995). The precession of the orbit is caused by the large quadrupole moment of the B star (Kaspi et al. 1996) The measurement of $`\ddot{x}`$ will then allow the separate determination of the angle $`\mathrm{\Phi }_0`$ and $`\chi \mathrm{sin}\theta `$. If we can determine the masses of the bodies comprising the binary system without measuring $`\dot{\omega }`$ (e.g., using $`\gamma `$ and $`\dot{P}_b`$) then, in principle, we can extract $`\dot{\omega }_S`$ from the total observed value of $`\dot{\omega }`$ and finally determine $`\chi `$. It is important to see that neither $`\ddot{\omega }`$ nor even higher derivatives of $`x`$ and $`\omega `$ contain a term $`\chi \mathrm{cos}\theta `$. Thus the only way to fully determine the spin of the BH is the separate measurement of $`\dot{\omega }_S`$.

Unfortunately, the precession of the orbit due to spin-orbit coupling is not the only effect that gives rise to the parameter $`\dot{x}`$ in the timing model. Following Damour & Taylor (1992) and Kopeikin (1994, 1996) one finds, in fact, four alternative effects causing a secular change of the observed value, $`x^{\mathrm{obs}}`$, of the parameter $`x`$:

$$\left(\frac{\dot{x}}{x}\right)^{\mathrm{obs}}=\left(\frac{\dot{x}}{x}\right)^{\mathrm{prec}}+\left(\frac{\dot{a}_p}{a_p}\right)^{\mathrm{gw}}+\left(\frac{\dot{x}}{x}\right)^{\mathrm{pm}}+\frac{d\epsilon _A}{dt}\frac{\dot{D}}{D}.$$
(65)   
The first contribution to the right-hand side of equation (65) is due to the precession of the binary orbit discussed at the beginning of this section. The contribution to the right-hand side of equation (65) from the second (gw) term is caused by the shrinking of the pulsar orbit due to the emission of gravitational waves by the binary system. Given the masses of pulsar and companion this effect can be calculated (Peters (1964)) and subtracted. The third (pm) contribution is due to the proper motion of the binary system, which causes a gradual change in the apparent geometrical orientation of the orbital plane (Kopeikin 1994, 1996; see also Arzoumanian et al. (1996)).<sup>4</sup><sup>4</sup>4Apart from a contribution to $`\dot{x}`$ this effect changes also the observed value for the advance of periastron, $`\dot{\omega }`$ (Kopeikin 1994, 1996). If the pulsar is close to the solar system and/or has a high proper motion this contribution can become quite significant, in which case, additional astrometric information will be needed for measurement and subtraction of $`\dot{x}^{\mathrm{pm}}`$. The forth contribution to the right-hand side of the equation (65) is due to a varying aberration caused by the deSitter precession of the pulsar’s spin while the pulsar orbits its companion. The pulsar spin precesses around the direction of the orbital angular momentum, $`𝐋`$, with angular velocity (Barker & O’Connell (1975))

$$\mathrm{\Omega }_p^{\mathrm{spin}}=\frac{T_{}^{2/3}}{1e^2}\left(\frac{2\pi }{P_b}\right)^{5/3}\frac{_{}(4_p+3_{})}{2(_p+_{})^{4/3}}$$
(66)   
and, consequently, (Damour & Taylor (1992))

$$\frac{d\epsilon _A}{dt}=\frac{P}{P_b}\frac{\mathrm{\Omega }_p^{\mathrm{spin}}}{(1e^2)^{1/2}}\frac{\mathrm{cot}\lambda _p\mathrm{sin}2\eta _p+\mathrm{cot}i\mathrm{cos}\eta _p}{\mathrm{sin}\lambda _p},$$
(67)   
where $`P`$ is the pulsar’s rotational period. A comparison with equation (61) leads to (excluding values of $`\lambda _p`$ close to 0°and 180°)

$$\left|\frac{\dot{\epsilon }_A}{\dot{x}_S/x}\right|\frac{P}{P_b}\frac{\mathrm{\Omega }_p^{\mathrm{spin}}}{\mathrm{\Omega }_S^{}}0.006\left(\frac{P}{1\mathrm{s}}\right)\left(\frac{P_b}{1\mathrm{day}}\right)^{2/3}\left(\frac{_{}}{10}\right)^{1/3}.$$
(68)   
Thus, the varying aberration effect will be of the same order as the precession of the orbit only for slow pulsars ($`P1`$ s) in very close orbits with a BH companion. But in this case the varying aberration will also show up as a change of the eccentricity (Damour & Taylor (1992))

$$\left(\frac{\dot{e}}{e}\right)^{\mathrm{obs}}=\left(\frac{\dot{e}}{e}\right)^{\mathrm{gw}}+\frac{d\epsilon _A}{dt},$$
(69)   
and can therefore be observed and subtracted. (The change due to gravitational wave damping (gw) can be calculated if the masses are known (see Peters (1964)).)

The last term in equation (65) is due to a varying Doppler shift, $`D`$, caused by a secular change of the distance between the solar-system barycenter and the binary-pulsar system. The effect under discussion is non-linear in time and can be caused both by the acceleration of gravitational fields in our galaxy and an apparent acceleration due to the proper motion of the binary system in the sky (Shklovskii (1970); Damour & Taylor (1991)):

$$\frac{\dot{D}}{D}=\frac{1}{c}𝐊_0(𝐚_{\mathrm{binary}\mathrm{system}}𝐚_{\mathrm{solar}\mathrm{system}})+\frac{V_T^2}{cd}.$$
(70)   
Here, $`V_T`$ denotes the velocity of the pulsar being transverse to the line of sight, $`d`$ is the distance between the pulsar and the solar system, and $`𝐚_{\mathrm{binary}\mathrm{system}}`$ and $`𝐚_{\mathrm{solar}\mathrm{system}}`$ are galactic accelerations at the locations of the solar and binary systems respectively. The galactic acceleration at the solar location is of order $`2\times 10^8\mathrm{cm}\mathrm{s}^2`$ which is a typical value for the acceleration fields in the galactic disk (Carlberg & Innanen (1987); Kuijken & Gilmore (1989) for a model of the Galactic acceleration field). Hence, unless the pulsar is located in the core of a globular cluster where accelerations are much higher, the contribution of galactic gravitational fields to $`\dot{x}`$ is negligible. The same is true for the Shklovskii term if we assume typical distances ($`1`$ kpc) and typical pulsar velocities (few 100 km/s). Thus, we conclude that a varying Doppler shift should not complicate the interpretation of timing observations of BH-pulsar binaries under ordinary circumstances.

As pointed out above, for a full determination of the rotational parameter $`\chi `$ of the pulsar companion, one cannot make use of $`\dot{\omega }`$ to calculate the masses of pulsar and BH. Thus, one needs a third PK parameter for the separate mass determination which will be most likely $`\dot{P}_b`$ (caused by gravitational wave damping). Since (Damour & Taylor (1991))

$$\left(\frac{\dot{P}_b}{P_b}\right)^{\mathrm{obs}}=\left(\frac{\dot{P}_b}{P_b}\right)^{\mathrm{gw}}\frac{\dot{D}}{D},$$
(71)   
and

$$\left(\frac{\dot{P}_b}{P_b}\right)^{\mathrm{gw}}\left(1.6\times 10^{18}\mathrm{s}^1\right)\left(\frac{P_b}{1\mathrm{day}}\right)^{2/3}\left(\frac{_{}}{10}\right)^{2/3}$$
(72)   
one has to be cautious about the contribution of the galactic acceleration and Shklovskii’s term. Even if we understand these effects to about 10% for BH-pulsar binaries with orbital periods greater than a few days, a reliable determination of $`\dot{\omega }_S`$ will be impossible. Exceptions are cases where the system is seen sufficiently edge-on so that one can measure the two PK parameters related with the Shapiro delay, $`r`$ and $`s`$, with necessary precision. For PSR B1913+16, the limit on the determination of $`\dot{P}_b^{\mathrm{gw}}`$ due to galactic accelerations has been clearly demonstrated by Damour & Taylor (1991).

For a BH-pulsar binary the deSitter precession of the pulsar spin, equation (66), is approximately $`\dot{\omega }_{PN}/2`$. For short orbital periods this can be a few degrees per year or even more. Then, if the pulsar spin is tilted with respect to the orbital angular momentum the deSitter precession can lead to significant changes in the pulse profile (cf. Weisberg, Romani, & Taylor (1989); Kramer (1998)) which might, in principle, limit the timing accuracy. In extreme cases the pulsar could even cease to beam in our direction after a couple of years.

Only if the influence of all these additional effects listed above is well understood we can hope to perform a spin determination for the BH companion of the pulsar.

## 6 Parameter estimations

So far we have shown that in principle, one can use the pulsar timing observations of a BH-pulsar binary system to determine the mass and spin of the rotating BH companion. In the previous section, we discussed in detail the extraction of the BH mass and spin from the observed values of Keplerian and PK parameters. However, the number of measured parameters one hopes to determine from timing observations depends critically on both the accuracy of the measured TOAs and the compactness of the binary system’s orbit. For example, we show that an ordinary one-second period pulsar in a one-year orbit with a 10 $`M_{}`$ BH will allow at most the determination of the mass function and sets only a lower limit on the mass of the BH companion. However, for a millisecond pulsar in an eccentric one-day orbit with a $`10M_{}`$ BH, the masses of the pulsar and companion can be determined with high precision after a month of regular timing observations from the measurement of two PK parameters (most likely $`\dot{\omega }`$ and $`\gamma `$). After a one-year observing campaign one can expect to see additional relativistic effects, such as gravitational radiation damping and relativistic spin-orbit precession. In this section, we explore the potential measurability of various relativistic effects in BH-pulsar binaries. We use a standard tool, the information matrix, to asses the errors of our parameter estimation method.

Let $`N`$ denote the total number of measured TOAs, $`\tau _i`$ ($`i=1,2,\mathrm{},N`$), and $`M`$ the number of fitting parameters, $`\xi _a`$ ($`a=1,2,\mathrm{},M`$), in the timing model used for the a priori estimate of $`𝒩(\tau ,𝝃)`$ the number of the pulse arriving at the time $`\tau `$. The parameters are determined using a least-square minimization method where the goodness of fit parameter, $`\chi ^2`$, is given by

$$\chi ^2=\underset{i=1}{\overset{N}{}}\left[\frac{n_i𝒩(\tau _i,𝝃)}{\nu _p\sigma _i}\right]^2.$$
(73)   
Here $`\nu _p`$ is the pulsar’s rotational frequency, $`n_i`$ is the closest integer to $`𝒩(\tau _i,𝝃)`$, and $`\sigma _i`$ is the estimated uncertainty of the $`i`$-th TOA. The information matrix is defined as

$$L_{ab}=\underset{i=1}{\overset{N}{}}\frac{𝒩(\tau _i)}{\xi _a}\frac{𝒩(\tau _i)}{\xi _b}.$$
(74)   
If only white noise is present in the TOA residuals, the inverse of the information matrix, $`M_{ab}L_{ab}^1`$, is the correlation matrix of the fitting parameters. The elements of the main diagonal of $`M_{ab}`$ give the variations of the measured parameters, $`\xi _a^2=M_{aa}`$, and the off-diagonal elements of $`M_{ab}`$ represent the correlations between them (Bard (1974)).

Two kinds of BH-pulsar binaries seem likely to be present in our Galaxy. First, young binary pulsars located in the Galactic plane, where the pulsar is the result of the second supernova explosion. Numerical simulations indicate, that such a pulsar is expected to be in eccentric wide ($`P_b10`$ days) orbits (Lipunov et al. (1994)). Secondly, a millisecond pulsar captured by a BH, most likely located in a globular cluster (Kulkarni, Hut, & McMillan (1993)).<sup>5</sup><sup>5</sup>5The ratio of orbital binding energy to potential energy in the gravitational field of the cluster is only a very weak function of the companion mass, once $`m_c1.4M_{}`$. Thus, BH-pulsar binaries should have roughly the same probability to remain in or near the cluster after formation as NS-NS binaries, like PSR B2127+11C. The aim of our numerical estimates is to understand which parameters one can expect to measure in these two different types of BH-pulsar binaries within a reasonable time span of observations ($`20`$ years). Here we assume that all measured TOAs, $`\tau _i`$, are of a similar quality, i.e. all $`\sigma _i`$ are the same ($`\sigma _i=\sigma _{\mathrm{TOA}}`$). Furthermore, there is at least one timing observation every month after an initial set of 10 timing observations along one full orbit which give the Keplerian parameters required for subsequent exploration of the binary system. We further assume that the pulsar is in a highly eccentric orbit ($`e=0.8`$) with a $`10M_{}`$ BH companion.

Figure 4 shows the estimated time span of regular timing observations, $`T_{\mathrm{obs}}`$, needed for the measurement of certain astrophysical parameters or effects, as a function of the orbital period, $`P_b`$. As astrophysical parameters and effects we consider:

* Total mass of the binary system — measurement of $`\dot{\omega }`$.
* Mass of the BH — measurement of $`\dot{\omega }`$ and $`\gamma `$.
* Testing the emission of gravitational waves — measurement of $`\dot{\omega }`$, $`\gamma `$ and $`\dot{P}_b`$ (cf. Taylor & Weisberg (1989)).
* Frame-dragging effect — measurement of $`\dot{x}`$.
* Spin of the BH (magnitude and orientation) — measurement of $`\dot{\omega }`$, $`\gamma `$, $`\dot{P}_b`$, $`\dot{x}`$, and $`\ddot{x}`$. The mass determination has to be done using $`\gamma `$ and $`\dot{P}_b`$. The accuracy for the total mass has to better than $`\epsilon \dot{\omega }_S/\dot{\omega }`$, where $`\epsilon `$ denotes the aspired accuracy for the BH spin measurement.

In figure 4a we have assumed a timing accuracy of $`\sigma _{\mathrm{TOA}}100\mu `$s which represents the optimum one can expect for timing observations of a young pulsar, with a typical rotational period of 100 ms to a few seconds. In figure 4b we have assumed the presence of a millisecond pulsar, i.e. $`\sigma _{\mathrm{TOA}}1\mu `$s.

We conclude from figure 4a that for a young pulsar only in systems with orbital periods below 0.2 days the measurement of the spin of the BH (magnitude and orientation) is possible within a reasonable observing time span. For systems with orbital periods below 3 days one can expect to see the frame-dragging caused by the extreme Kerr BH, i.e. measure $`\dot{x}`$. Systems with orbital periods below 10 days will allow a mass determination for the BH with an accuracy of better than 5% after less than 20 years of observing. The periastron advance will be a prominent effect even for orbital periods exceeding one year and thus the total mass of the binary system should be a measureable quantity for all orbital periods $``$1000 days.

In the case of a millisecond pulsar in orbit with a 10 solar mass extreme Kerr BH, where one expects timing accuracies two orders of magnitude better than in the previous case, the corresponding time spans for observations are clearly shorter (see figure 4b). Therefore the spin of the BH will be measurable for systems with orbital periods up to one day, and the frame-dragging will be seen even in systems with 10 days orbital period. For systems with orbital periods below a hundred days a good mass estimates for the BH companion will be possible, and the total mass of the binary system will be measurable even in very wide systems with orbital periods of a few years.

The previous calculations were done for a BH mass of $`10M_{}`$, which we assume to be typical for stellar-mass BHs. In a $`30M_{}`$ BH-pulsar binary spin-related effects are expected to be an order of magnitude larger, since the spin of a BH scales with its mass squared (see equation (4)). On the other hand, there might be a fairly large population of BHs in our Galaxy with masses slightly above $`1.5M_{}`$ (Brown & Bethe (1994)). Our calculations show that only for very short orbital periods ($`P_b0.2`$ days) and high orbital eccentricities, ($`e0.8`$) one can expect to see the relativistic spin-orbit precession ($`\dot{x}_S`$) within less than 20 years of observing, given an uncertainty in the TOAs which is not worse than $`1\mu `$s. In this case the measurement of $`\dot{x}_S`$, which gives a lower limit for $`\chi `$, might have the potential to distinguish between a NS and a BH, since it seems that $`\chi 0.7`$ for NSs of comparable mass (Laarakkers & Poisson (1997)). It would require that the angle $`\theta `$ is close to $`90\mathrm{°}`$ (the spin of the BH lies almost in the orbital plane) and $`\mathrm{\Phi }_0`$ is either close to $`90\mathrm{°}`$ or close to $`270\mathrm{°}`$. Then, if $`\chi `$ is close to one for the BH companion, the minimum value for $`\chi `$ derived from the $`\dot{x}`$ measurement can exceed the maximum value for NSs.

Let us point out that on time scales of several years, timing observations are usually contaminated by the increasing amount of red (low frequency) timing noise (Cordes & Helfand (1980), Taylor (1991), Lyne (1996), Kopeikin 1997b ). If the red noise dominates in TOAs residuals, it makes the problem of measuring relativistic parameters more difficult. In particular for young pulsars in wide orbits, the separation of orbital effects and timing noise may cause a severe problem (cf. Wex et al. (1998)).

## 7 Measuring the quadrupole of the BH companion

Equation (7) implies that with the determination of the mass and the spin of a Kerr BH we also know its quadrupole moment. Therefore, if we are able to extract independently the quadrupole moment of the companion from our timing observations we could actually test whether the observed pulsar is orbiting a Kerr BH or another compact relativistic object like neutron or boson star. As discussed in Section 4 the quadrupole moment of the BH companion will lead to an additional precession in the angles $`\mathrm{\Phi }`$ and $`\mathrm{\Psi }`$ and, thus, to a secular time evolution in the parameters $`\omega `$ and $`x`$. Unfortunately these secular changes in the orientation of the orbit caused by the quadrupole moment of a BH companion are typically three orders of magnitude smaller than the changes caused by the relativistic spin-orbit coupling (see equation (49)). Thus, the chance of extracting the quadrupole moment of the BH from the measurement of the parameters $`\dot{\omega }`$ and $`\dot{x}`$ is unlikely. On the other hand, the anisotropic nature of the quadrupole component of the external gravitational field (2) will lead to characteristic short-term periodic effects every time the pulsar gets close to the oblate companion. A detailed investigation of these short-term effects can be found in Wex (1998). These short-term periodic effects lead to a unique signature, which will show up in the post-fit TOA residuals if one uses a timing model that accounts only for the secular changes of the precession and the short-term periodic effects caused by the mass-monopole (Damour & Deruelle (1986)) and spin (Wex (1995)) of the BH. Figure 5a shows the result of numerical simulations of TOA residuals for a pulsar orbiting a 30 $`M_{}`$ BH in an 0.1 day orbit with an eccentricity of 0.9. This represents certainly an unrealistic case from the observational point of view since the lifetime of such a system due to gravitational radiation damping is just $`10^4`$ years which makes a discovery of such a system extremely unlikely. In addition, the accuracy needed to measure this effect would require a nano-second timing precision. Even for the best millisecond pulsars the present timing precision is still more than two orders of magnitude worse and, thus, we conclude that the quadrupole of a stellar-mass BH companion is not measurable in timing observations. However, since the quadrupole moment of a BH scales with the mass cubed (see equation 5), these short-term periodic effects under discussion gain importance for very massive BHs which could sit in the center of globular clusters as a result of a collapse of the cluster core. Figure 5b gives numerical simulations of TOA residuals for a pulsar orbiting a $`10^4`$ $`M_{}`$ BH in a 10 day orbit with an eccentricity of 0.9. A timing accuracy of a few hundred nanoseconds would allow the measurement of the quadrupole moment. Of course, even higher BH masses and/or smaller pulsar orbits will make the quadrupole a prominent feature in the TOAs.

## 8 Summary

In this paper we have examined how timing observations of a pulsar orbiting a rotating stellar mass BH can be used to study the physical properties of the BH companion. We pointed out that the measurement of two post-Keplerian parameters can lead to a mass determination of the black hole with very high accuracy, far better than it is possible for present black hole candidates. It was shown that the frame dragging propagation effect discussed by Laguna & Wolszczan (1997) will not be separable from the bending delay for those BH companions less than a hundred solar masses, since it would require determining the orientation of the pulsar’s spin with an accuracy that is not likely possible with pulse structure analysis. We gave detailed calculations about the (secular) influence of the BH rotation on the orbital motion of the pulsar. In particular the long term behavior of this relativistic spin-orbit coupling can be described by two angles, defined with respect to the invariable plane, which change linearly in time. For an observer this linear-in-time precession converts into a non-linear in time evolution of the orbital inclination and the longitude of periastron. In a timing model for such a binary system, this can be taken into account by fitting for the parameters $`\dot{x}`$, $`\ddot{x}`$, $`\dot{\omega }`$, and $`\ddot{\omega }`$ as is done for pulsars with main-sequence star companions where the classical spin-orbit coupling is the cause of the precession (see Wex (1998)). The measurement of $`\dot{x}`$, $`\ddot{x}`$, and $`\dot{\omega }`$ could, in principle, lead to the determination of the BH spin, if the masses of the BH and the pulsar can be determined without making use of the advance of periastron, using, for example, the measurement of the Einstein delay ($`\gamma `$) and the orbital period decay ($`\dot{P}_b`$). In order to do so, a fractional measurement precision of $`10^4`$ is required for $`\dot{\omega }`$, $`\gamma `$, and $`\dot{P}_b`$. If this measurement precision is not achieved, or if there are Galactic and kinematic effects present which are not sufficiently understood, then only the projection of the BH spin onto the orbital angular momentum will be measurable. If the timing data allow just the determination of $`\dot{x}`$ then we can calculate only a lower limit for the projected BH spin, which, of course, is also a lower limit for the BH spin. However if we have just the lower limit for the spin, the presence of the $`\dot{x}`$ would already indicate frame dragging if an oblate non-compact companion can be ruled out (for example by optical observations, the absence of eclipses of the pulsar, absence of tidal effects, etc.). On the other hand, the method under discussion can be used to rule out a Kerr BH companion as has been demonstrated here in the case of PSR J0045–7319.

We have studied the measurability of the relativistic spin orbit coupling for two different cases. First, a young pulsar orbiting a $`10M_{}`$ extreme Kerr BH and second, a millisecond pulsar orbiting a $`10M_{}`$ extreme Kerr BH. In the first case only for binary pulsars with orbital periods less than three days one can expect to measure $`\dot{x}`$, i.e. the frame dragging caused by the rotating BH, with sufficient accuracy, but in the second case a measurement of $`\dot{x}`$ seems very likely. For BH-millisecond pulsar binaries with an orbital period below one day one can expect a measurement of $`\ddot{x}`$ (and $`\ddot{\omega }`$) within a reasonable time span of observations, and one will be able to fully determine the BH spin parameters (magnitude and orientation). However, if the pulsar is located close to the core of a globular cluster one might not have the necessary understanding for the gravitational cluster potential in order to take properly into account its contribution to the measured binary parameters. For a $`1.5M_{}`$ BH-millisecond pulsar binary with $`P_b0.2`$ days and $`e0.8`$ one might be able to measure $`\dot{x}_S`$. In principle, this could help to exclude a neutron star, if the magnitude and orientation of the BH companion is such that $`\chi \mathrm{sin}\theta \mathrm{sin}\mathrm{\Phi }_00.7`$, since NSs seem to have $`\chi 0.7`$.

Finally we pointed out that the quadrupole field of the rotating BH will lead to a distinctive signature in the post-fit TOA residuals during each periastron passage. However, our numerical simulations suggest that this would be observable only for very massive BHs, typically more than $`10^4M_{}`$. Whether BHs of this size do exist, e.g. in the center of globular clusters, is still unclear. The discovery of a pulsar orbiting the super-massive BH in the center of our Galaxy could provide a possible setting for these measurements.

###### Acknowledgements.

We thank B. Paczynski for attracting our attention to the problem, P. Laguna, D. Lorimer, S. Thorsett, and A. Wolszczan for valuable discussions, and R. Dewey for carefully reading the manuscript. Norbert Wex acknowledges the hospitality of the Department of Physics at Princeton University during tenure of the Otto-Hahn Prize. S. M. Kopeikin is pleasured to acknowledge the hospitality of G. Neugebauer and G. Schäfer and other members of the Institute for Theoretical Physics of the Friedrich Schiller University of Jena. This work has been partially supported by the Thüringer Ministerium für Wissenschaft, Forschung und Kultur grant No B501-96060.