Spectral methods in general relativistic astrophysics

We present spectral methods developed in our group to solve three-dimensional partial differential equations. The emphasis is put on equations arising from astrophysical problems in the framework of general relativity.Comment: 51 pages, elsart (Elsevier Preprint), 19 PostScript figures, submitted to Journal of Computational & Applied Mathematic

Merger of black hole-neutron star binaries: nonspinning black hole case

We perform a simulation for merger of a black hole (BH)-neutron star (NS) binary in full general relativity preparing a quasicircular state as initial condition. The BH is modeled by a moving puncture with no spin and the NS by the $\Gamma$-law equation of state with $\Gamma=2$. Corotating velocity field is assumed for the NS. The mass of the BH and the rest-mass of the NS are chosen to be $\approx 3.2 M_{\odot}$ and $\approx 1.4 M_{\odot}$ with relatively large radius of the NS $\approx 14$ km. The NS is tidally disrupted near the innermost stable orbit but $\sim 80$ of the material is swallowed into the BH with small disk mass $\sim 0.3M_{\odot}$ even for such small BH mass $\sim 3M_{\odot}$. The result indicates that the system of a BH and a massive disk of $\sim M_{\odot}$ is not formed from nonspinning BH-NS binaries, although a disk of mass $\sim 0.1M_{\odot}$ is a possible outcome.Comment: 5 pages. Phys. Rev. D 74, 121503 (R) (2006

Darwin-Riemann problems in general relativity

A review is given of recent results about the computation of irrotational Darwin-Riemann configurations in general relativity. Such configurations are expected to represent fairly well the late stages of inspiralling binary neutron stars.Comment: 20 pages, 11 PostScript figures, uses PTPTeX, to appear in the Proceedings of Yukawa International Seminar 99 "Black Holes and Gravitational Waves", edited by T. Nakamura & H. Kodama, Prog. Theor. Phys. Supp

Biodegradability and tissue reaction of random copolymers of L-leucine, L-aspartic acid, and L-aspartic acid esters

A series of copoly(α-amino acids) with varying percentages of hydrophilic (l-aspartic acid) and hydrophobic monomers (l-leucine, ß-methyl-l-aspartate, and ß-benzyl-l-aspartate) were implanted subcutaneously in rats and the macroscopic degradation behavior was studied. Three groups of materials (A, B, C) with different ranges of hydrophilicity were distinguished: A) hydrophobic materials showed no degradation after 12 weeks; B) more hydrophilic materials revealed a gradual reduction in size of the samples, but were still present after 12 weeks; and C) hydrophilic copolymers disappeared within 24 hr. \ud The tissue reactions caused by the materials of group A resembled that of silicone rubber, whereas those of group B showed a more cellular reaction

A Dynamical Systems Approach to Schwarzschild Null Geodesics

The null geodesics of a Schwarzschild black hole are studied from a dynamical systems perspective. Written in terms of Kerr-Schild coordinates, the null geodesic equation takes on the simple form of a particle moving under the influence of a Newtonian central force with an inverse-cubic potential. We apply a McGehee transformation to these equations, which clearly elucidates the full phase space of solutions. All the null geodesics belong to one of four families of invariant manifolds and their limiting cases, further characterized by the angular momentum L of the orbit: for |L|>|L_c|, (1) the set that flow outward from the white hole, turn around, then fall into the black hole, (2) the set that fall inward from past null infinity, turn around outside the black hole to continue to future null infinity, and for |L|<|L_c|, (3) the set that flow outward from the white hole and continue to future null infinity, (4) the set that flow inward from past null infinity and into the black hole. The critical angular momentum Lc corresponds to the unstable circular orbit at r=3M, and the homoclinic orbits associated with it. There are two additional critical points of the flow at the singularity at r=0. Though the solutions of geodesic motion and Hamiltonian flow we describe here are well known, what we believe is a novel aspect of this work is the mapping between the two equivalent descriptions, and the different insights each approach can give to the problem. For example, the McGehee picture points to a particularly interesting limiting case of the class (1) that move from the white to black hole: in the limit as L goes to infinity, as described in Schwarzschild coordinates, these geodesics begin at r=0, flow along t=constant lines, turn around at r=2M, then continue to r=0. During this motion they circle in azimuth exactly once, and complete the journey in zero affine time.Comment: 14 pages, 3 Figure

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