Black hole tidal problem in the Fermi normal coordinates

We derive a tidal potential for a self-gravitating fluid star orbiting Kerr black hole along a timelike geodesic extending previous works by Fishbone and Marck. In this paper, the tidal potential is calculated up to the third and fourth-order terms in $R/r$, where $R$ is the stellar radius and $r$ the orbital separation, in the Fermi-normal coordinate system following the framework developed by Manasse and Misner. The new formulation is applied for determining the tidal disruption limit (Roche limit) of corotating Newtonian stars in circular orbits moving on the equatorial plane of Kerr black holes. It is demonstrated that the third and fourth-order terms quantitatively play an important role in the Roche limit for close orbits with R/r \agt 0.1. It is also indicated that the Roche limit of neutron stars orbiting a stellar-mass black hole near the innermost stable circular orbit may depend sensitively on the equation of state of the neutron star.Comment: Correct typo

SHORTCUT METHOD OF SOLUTION OF GEODESIC EQUATIONS FOR SCHWARZSCHILD BLACK HOLE

It is shown how the use of the Kerr-Schild coordinate system can greatly simplify the formulation of the geodesic equation of the Schwarzschild solution. An application of this formulation to the numerical computation of the aspect of a non-rotating black hole is presented. The generalization to the case of the Kerr solution is presented too.Comment: 11 pages, 2 PostScript figures (available as uuencoded compressed tar file), uses epsfig.tex). Accepted on February 1995 for publication in Classical and Quantum Gravit

Tidal Interaction between a Fluid Star and a Kerr Black Hole in Circular Orbit

We present a semi-analytic study of the equilibrium models of close binary systems containing a fluid star (mass $m$ and radius $R_0$) and a Kerr black hole (mass $M$) in circular orbit. We consider the limit $M\gg m$ where spacetime is described by the Kerr metric. The tidally deformed star is approximated by an ellipsoid, and satisfies the polytropic equation of state. The models also include fluid motion in the stellar interior, allowing binary models with nonsynchronized stellar spin (as expected for coalescing neutron star-black hole binaries) to be constructed. Tidal disruption occurs at orbital radius $r_{\rm tide}\sim R_0(M/m)^{1/3}$, but the dimensionless ratio $\hat r_{\rm tide}=r_{\rm tide}/[R_0(M/m)^{1/3}]$ depends on the spin parameter of the black hole as well as on the equation of state and the internal rotation of the star. We find that the general relativistic tidal field disrupts the star at a larger $\hat r_{\rm tide}$ than the Newtonian tide; the difference is particularly prominent if the disruption occurs in the vicinity of the black hole's horizon. In general, $\hat r_{\rm tide}$ is smaller for a (prograde rotating) Kerr black hole than for a Schwarzschild black hole. We apply our results to coalescing black hole-neutron star and black hole-white dwarf binaries. The tidal disruption limit is important for characterizing the expected gravitational wave signals and is relevant for determining the energetics of gamma ray bursts which may result from such disruption.Comment: 29 pages including 8 figures. Minor changes and update. To appear in ApJ, March 20, 2000 (Vol.532, #1

An approximate binary-black-hole metric

An approximate solution to Einstein's equations representing two widely-separated non-rotating black holes in a circular orbit is constructed by matching a post-Newtonian metric to two perturbed Schwarzschild metrics. The spacetime metric is presented in a single coordinate system valid up to the apparent horizons of the black holes. This metric could be useful in numerical simulations of binary black holes. Initial data extracted from this metric have the advantages of being linked to the early inspiral phase of the binary system, and of not containing spurious gravitational waves.Comment: 20 pages, 1 figure; some changes in Sec. IV B,C and Sec.

Burst dynamics during drainage displacements in porous media: Simulations and experiments

We investigate the burst dynamics during drainage going from low to high injection rate at various fluid viscosities. The bursts are identified as pressure drops in the pressure signal across the system. We find that the statistical distribution of pressure drops scales according to other systems exhibiting self-organized criticality. The pressure signal was calculated by a network model that properly simulates drainage displacements. We compare our results with corresponding experiments.Comment: 7 pages, 4 figures. Submitted to Europhys. Let

Numerical approach for high precision 3-D relativistic star models

A multi-domain spectral method for computing very high precision 3-D stellar models is presented. The boundary of each domain is chosen in order to coincide with a physical discontinuity (e.g. the star's surface). In addition, a regularization procedure is introduced to deal with the infinite derivatives on the boundary that may appear in the density field when stiff equations of state are used. Consequently all the physical fields are smooth functions on each domain and the spectral method is absolutely free of any Gibbs phenomenon, which yields to a very high precision. The power of this method is demonstrated by direct comparison with analytical solutions such as MacLaurin spheroids and Roche ellipsoids. The relative numerical error reveals to be of the order of $10^{-10}$. This approach has been developed for the study of relativistic inspiralling binaries. It may be applied to a wider class of astrophysical problems such as the study of relativistic rotating stars too.Comment: Minor changes, Phys. Rev. D in pres

Retarded coordinates based at a world line, and the motion of a small black hole in an external universe

In the first part of this article I present a system of retarded coordinates based at an arbitrary world line of an arbitrary curved spacetime. The retarded-time coordinate labels forward light cones that are centered on the world line, the radial coordinate is an affine parameter on the null generators of these light cones, and the angular coordinates are constant on each of these generators. The spacetime metric in the retarded coordinates is displayed as an expansion in powers of the radial coordinate and expressed in terms of the world line's acceleration vector and the spacetime's Riemann tensor evaluated at the world line. The formalism is illustrated in two examples, the first involving a comoving world line of a spatially-flat cosmology, the other featuring an observer in circular motion in the Schwarzschild spacetime. The main application of the formalism is presented in the second part of the article, in which I consider the motion of a small black hole in an empty external universe. I use the retarded coordinates to construct the metric of the small black hole perturbed by the tidal field of the external universe, and the metric of the external universe perturbed by the presence of the black hole. Matching these metrics produces the MiSaTaQuWa equations of motion for the small black hole.Comment: 20 pages, revtex4, 2 figure

Numerical models of irrotational binary neutron stars in general relativity

We report on general relativistic calculations of quasiequilibrium configurations of binary neutron stars in circular orbits with zero vorticity. These configurations are expected to represent realistic situations as opposed to corotating configurations. The Einstein equations are solved under the assumption of a conformally flat spatial 3-metric (Wilson-Mathews approximation). The velocity field inside the stars is computed by solving an elliptical equation for the velocity scalar potential. Results are presented for sequences of constant baryon number (evolutionary sequences). Although the central density decreases much less with the binary separation than in the corotating case, it still decreases. Thus, no tendency is found for the stars to individually collapse to black hole prior to merger.Comment: Minor corrections, improved figure, 5 pages, REVTeX, Phys. Rev. Lett. in pres