Search

We performed a series of 739 full numerical simulations of high energy collision of black holes to search for the maximum gravitational energy emitted $E_{rad}$, during their merger. We consider equal mass binaries with spins pointing along their orbital angular momentum $\vec{L}$ and perform a search over impact parameters $b$ and initial linear momenta $p/m=\gamma v$ to find the maximum $E_{rad}$ for a given spin $\vec{S}$. The total radiated energy proves to have a weak dependence on the intrinsic spin $s$ of the holes, for the sequence $s=+0.8, 0.0, -0.8$ studied here. We thus estimate the maximum $E_{rad}^{max}/M_{ADM}\approx25\%\pm2\%$ for these direct merger encounters. We also explore the radiated angular momentum and the maximum spin of the merger remnant (within these configurations), finding $\alpha_f^{max}=0.987$. We then use the zero frequency limit expansion to analytically model the radiated energy in the small impact parameter and large initial linear momentum regime., Comment: 11 pages, 7 figures, 4 tables

We performed a series of 1381 full numerical simulations of high energy collision of black holes to search for the maximum recoil velocity after their merger. We consider equal mass binaries with opposite spins pointing along their orbital plane and perform a search of spin orientations, impact parameters, and initial linear momenta to find the maximum recoil for a given spin magnitude $s$. This spin sequence for $s=0.4, 0.7, 0.8, 0.85, 0.9$ is then extrapolated to the extreme case, $s=1$, to obtain an estimated maximum recoil velocity of $28,562\pm 342$ km/s, thus approximately bounded by $10\%$ the speed of light., Comment: 8 figures, 1 Table

We revisit the definition of transverse frames and tetrad choices with regards to its application to numerically generated spacetimes, in particular those from the merger of binary black holes. We introduce the concept of local and approximate algebraic Petrov types in the strong field regime. We define an index $\mathcal{D}=\sqrt{12/I}\left(\Psi_2 - \Psi_3^2/\Psi_4\right)$ able to discriminate between Petrov types II and D and define regions of spacetime of those approximate types when used in conjunction with the speciality invariant $S=27J^2/I^3$. We provide an explicit example applying this method to Brill-Lindquist initial data corresponding to two nonspinning black holes from rest at a given initial separation. We find a doughnut-like region that is approximately of Petrov type II surrounded by an approximately Petrov type D region. We complete the study by proposing a totally symmetric tetrad fixing of the transverse frame that can be simply implemented in numerically generated spacetimes through the computation of spin coefficients ratios. We provide an application by explicitly deriving the Kerr-perturbative equations in this tetrad., Comment: 18 pages, 11 figures

We give a unified description of the flip-flop effect in spinning binary black holes and the anti-alignment instability in terms of real and imaginary flip-flop frequencies. We find that this instability is only effective for mass ratios $0.5

We use fully nonlinear numerical relativity techniques to study high energy head-on collision of nonspinning, equal-mass black holes to estimate the maximum gravitational radiation emitted by these systems. Our simulations include improvements in the construction of initial data, subsequent full numerical evolutions, and the computation of waveforms at infinity. The new initial data significantly reduces the spurious radiation content, allowing for initial speeds much closer to the speed of light, i.e. $v\sim0.99c$. Using these new techniques, We estimate the maximum radiated energy from head-on collisions to be $E_{\text{max}}/M_{\text{ADM}}=0.13\pm0.01$. This value differs from the second-order perturbative $(0.164)$ and zero-frequency-limit $(0.17)$ analytic computations, but is close to those obtained by thermodynamic arguments $(0.134)$ and by previous numerical estimates $(0.14\pm0.03)$., Comment: 11 pages, 10 figures

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011., Comment: 68 pages, 4 Figures

We study the statistical distributions of the spins of generic black-hole binaries during the inspiral and merger, as well as the distributions of the remnant mass, spin, and recoil velocity. For the inspiral regime, we start with a random uniform distribution of spin directions S1 and S2 and magnitudes S1=S2=0.97 for different mass ratios. Starting from a fiducial initial separation of ri=50m, we perform 3.5PN evolutions down to rf=5m. At this final fiducial separation, we compute the angular distribution of the spins with respect to the final orbital angular momentum, L. We perform 16^4 simulations for six mass ratios between q=1 and q=1/16 and compute the distribution of the angles between L and Delta and L and S, directly related to recoil velocities and total angular momentum. We find a small but statistically significant bias of the distribution towards counter-alignment of both scalar products. To study the merger of black-hole binaries, we turn to full numerical techniques. We introduce empirical formulae to describe the final remnant black hole mass, spin, and recoil velocity for merging black-hole binaries with arbitrary mass ratios and spins. We then evaluate those formulae for randomly chosen directions of the individual spins and magnitudes as well as the binary's mass ratio. We found that the magnitude of the recoil velocity distribution decays as P(v) \exp(-v/2500km/s), =630km/s, and sqrt{ - ^2}= 534km/s, leading to a 23% probability of recoils larger than 1000km/s, and a highly peaked angular distribution along the final orbital axis. The final black-hole spin magnitude show a universal distribution highly peaked at Sf/mf^2=0.73 and a 25 degrees misalignment with respect to the final orbital angular momentum., Comment: 22 pages, 28 figures, abridged abstract. Major addenda in section IV and discussion

We solve the problem of expressing the Weyl scalars $\psi $ that describe gravitational perturbations of a Kerr black hole in terms of Cauchy data. To do so we use geometrical identities (like the Gauss-Codazzi relations) as well as Einstein equations. We are able to explicitly express $\psi $ and $\partial _t\psi $ as functions only of the extrinsic curvature and the three-metric (and geometrical objects built out of it) of a generic spacelike slice of the spacetime. These results provide the link between initial data and $\psi $ to be evolved by the Teukolsky equation, and can be used to compute the gravitational radiation generated by two orbiting black holes in the close limit approximation. They can also be used to extract waveforms from spacetimes completely generated by numerical methods., Comment: 5 pages, REVTEX, no figures

We revisit the question of the imposition of initial data representing astrophysical gravitational perturbations of black holes. We study their dynamics for the case of nonrotating black holes by numerically evolving the Teukolsky equation in the time domain. In order to express the Teukolsky function Psi explicitly in terms of hypersurface quantities, we relate it to the Moncrief waveform phi_M through a Chandrasekhar transformation in the case of a nonrotating black hole. This relation between Psi and phi_M holds for any constant time hypersurface and allows us to compare the computation of the evolution of Schwarzschild perturbations by the Teukolsky and by the Zerilli and Regge-Wheeler equations. We explicitly perform this comparison for the Misner initial data in the close limit approach. We evolve numerically both, the Teukolsky (with the recent code of Ref. [1]) and the Zerilli equations, finding complete agreement in resulting waveforms within numerical error. The consistency of these results further supports the correctness of the numerical code for evolving the Teukolsky equation as well as the analytic expressions for Psi in terms only of the three-metric and the extrinsic curvature., Comment: 14 pages, 7 Postscript figures, to appear in Phys. Rev. D

Gravitational perturbations about a Kerr black hole in the Newman-Penrose formalism are concisely described by the Teukolsky equation. New numerical methods for studying the evolution of such perturbations require not only the construction of appropriate initial data to describe the collision of two orbiting black holes, but also to know how such new data must be imposed into the Teukolsky equation. In this paper we show how Cauchy data can be incorporated explicitly into the Teukolsky equation for non-rotating black holes. The Teukolsky function $% \Psi $ and its first time derivative $\partial_t \Psi $ can be written in terms of only the 3-geometry and the extrinsic curvature in a gauge invariant way. Taking a Laplace transform of the Teukolsky equation incorporates initial data as a source term. We show that for astrophysical data the straightforward Green function method leads to divergent integrals that can be regularized like for the case of a source generated by a particle coming from infinity., Comment: 9 pages, REVTEX. Misprints corrected in formulas (2.4)-(2.7). Final version to appear in PRD

Numerical relativity codes now being developed will evolve initial data representing colliding black holes at a relatively late stage in the collision. The choice of initial data used for code development has been made on the basis of mathematical definitiveness and usefulness for computational implementation. By using the ``particle limit'' (the limit of an extreme ratio of masses of colliding holes) we recently showed that the standard choice is not a good representation of astrophysically generated initial data. Here we show that, for the particle limit, there is a very simple alternative choice that appears to give excellent results. That choice, ``convective'' initial data is, roughly speaking, equivalent to the start of a time sequence of parameterized solutions of the Hamiltonian constraint; for a particle in circular orbit, it is the initial data of the steady state solution on any hypersurface. The implementation of related schemes for equal mass holes is discussed., Comment: 10 pages, REVTEX, 14 PS figs

We show that the radial Teukolsky equation (in the frequency domain) with sources that extend to infinity has well-behaved solutions. To prove that, we follow Poisson approach to regularize the non-rotating hole, and extend it to the rotating case. To do so we use the Chandrasekhar transformation among the Teukolsky and Regge-Wheeler-like equations, and express the integrals over the source in terms of solutions to the homogeneous Regge-Wheeler-like equation, to finally regularize the resulting integral. We then discuss the applicability of these results., Comment: 14 pages, 1 Table, REVTEX

Numerical relativity, applied to collisions of black holes, starts with initial data for black holes already in each other's strong field. The initial hypersurface data typically used for computation is based on mathematical simplifying prescriptions, such as conformal flatness of the 3-geometry and longitudinality of the extrinsic curvature. In the case of head on collisions of equal mass holes, there is evidence that such prescriptions work reasonably well, but it is not clear why, or whether this success is more generally valid. Here we study these questions by considering the ``particle limit'' for head on collisions of nonspinning holes. Einstein's equations are linearized in the mass of the small hole, and described by a single gauge invariant spacetime function psi, for each multipole. The resulting equations have been solved by numerical evolution for collisions starting from various initial separations, and the evolution is studied on a sequence of hypersurfaces. In particular, we extract hypersurface data, that is psi and its time derivative, on surfaces of constant background Schwarzschild time. These evolved data can then be compared with ``prescribed'' data, evolved data can be replaced by prescribed data on any hypersurface, and evolved further forward in time, a gauge invariant measure of deviation from conformal flatness can be evaluated, etc. The main findings of this study are: (i) For holes of unequal mass the use of prescribed data on late hypersurfaces is not successful. (ii) The failure is likely due to the inability of the prescribed data to represent the near field of the smaller hole. (iii) The discrepancy in the extrinsic curvature is more important than in the 3-geometry. (iv) The use of the more general conformally flat longitudinal data does not notably improve this picture., Comment: 20 pages, REVTEX, 26 PS figures included

I present here the results of the study of the gravitational radiation generated by the infall (from rest at radius $r_0$) of a point particle of mass $m_0$ into a Schwarzschild black hole of mass $M$. We use Laplace's transform methods and find that the spectra of radiation for $\sim5M

We find a self-consistent pp-gravitational shock wave solution to the semiclassical Einstein equations resulting from the $1/N$ approach to the effective action. We model the renormalized matter stress-energy-momentum tensor by $N$ massless scalar fields in the Minkowski vacuum plus a classical particle. We show that quantum effects generate a milder singularity at the position of the particle than the classical solution, but the singularity does not disappear. At large distances from the particle, the quantum correction decreases slowly, as $1/\rho^2$ ($\rho$ being the distance to the particle in the shock wave plane). We argue that this large distance correction is a necessary consequence of quantum gravity., Comment: 13 pages, REVTEX, 4 PS figures. Revised version contains a derivation of the solution compatible with its distributional character. The final results, though, are the same

We compute gravitational radiation waveforms, spectra and energies for a point particle of mass $m_0$ falling from rest at radius $r_0$ into a Schwarzschild hole of mass $M$. This radiation is found to lowest order in $(m_0/M)$ with the use of a Laplace transform. In contrast with numerical relativity results for head-on collisions of equal-mass holes, the radiated energy is found not to be a monotonically increasing function of initial separation; there is a local radiated-energy maximum at $r_0\approx4.5M$. The present results, along with results for infall from infinity, provide a complete catalog of waveforms and spectra for particle infall. We give a representative sample from that catalog and an interesting observation: Unlike the simple spectra for other head-on collisions (either of particle and hole, or of equal mass holes) the spectra for $\infty>r_0>\sim5M$ show a series of evenly spaced bumps. A simple explanation is given for this. Lastly, our energy vs. $r_0$ results are compared with approximation methods used elsewhere, for small and for large initial separation., Comment: 15 pages, REVTeX, 25 figures

In this paper we show how quantum corrections, although perturbatively small, may play an important role in the analysis of the existence of some classical models. This, in fact, appears to be the case of static, uniform--density models of the interior metric of cosmic strings and neutron stars. We consider the fourth order semiclassical equations and first look for perturbative solutions in the coupling constants $\alpha$ and $\beta$ of the quadratic curvature terms in the effective gravitational Lagrangian. We find that there is not a consistent solution; neither for strings nor for spherical stars. We then look for non--perturbative solutions and find an explicit approximate metric for the case of straight cosmic strings. We finally analyse the contribution of the non--local terms to the renormalized energy--momentum tensor and the possibility of this terms to allow for a perturbative solution. We explicitly build up a particular renormalized energy--momentum tensor to fulfill that end. These state--dependent corrections are found by simple considerations of symmetry, conservation law and trace anomaly, and are chosen to compensate for the local terms. However, they are not only ad hoc, but have to depend on $\alpha$ and $\beta$, what is not expected to first perturbative order. We then conclude that non--perturbative solutions are valuable for describing certain physical situations., Comment: 19 pages, REVTEX, no figures

We study the classical dynamics of a bosonic string in the $D$--dimensional flat Friedmann--Robertson--Walker and Schwarzschild backgrounds. We make a perturbative development in the string coordinates around a {\it null} string configuration; the background geometry is taken into account exactly. In the cosmological case we uncouple and solve the first order fluctuations; the string time evolution with the conformal gauge world-sheet $\tau$--coordinate is given by $X^0(\sigma, \tau)=q(\sigma)\tau^{1\over1+2\beta}+c^2B^0(\sigma, \tau)+\cdots$, $B^0(\sigma,\tau)=\sum_k b_k(\sigma)\tau^k$ where $b_k(\sigma)$ are given by Eqs.\ (3.15), and $\beta$ is the exponent of the conformal factor in the Friedmann--Robertson--Walker metric, i.e. $R\sim\eta^\beta$. The string proper size, at first order in the fluctuations, grows like the conformal factor $R(\eta)$ and the string energy--momentum tensor corresponds to that of a null fluid. For a string in the black hole background, we study the planar case, but keep the dimensionality of the spacetime $D$ generic. In the null string expansion, the radial, azimuthal, and time coordinates $(r,\phi,t)$ are $r=\sum_n A^1_{n}(\sigma)(-\tau)^{2n/(D+1)}~,$ $\phi=\sum_n A^3_{n}(\sigma)(-\tau)^{(D-5+2n)/(D+1)}~,$ and $t=\sum_n A^0_{n} (\sigma)(-\tau)^{1+2n(D-3)/(D+1)}~.$ The first terms of the series represent a {\it generic} approach to the Schwarzschild singularity at $r=0$. First and higher order string perturbations contribute with higher powers of $\tau$. The integrated string energy-momentum tensor corresponds to that of a null fluid in $D-1$ dimensions. As the string approaches the $r=0$ singularity its proper size grows indefinitely like $\sim(-\tau)^{-(D-3)/(D+1)}$. We end the paper giving three particular exact string solutions inside the black hole., Comment: 17 pages, REVTEX, no figures

We find an {\it exact} pp--gravitational wave solution of the fourth order gravity field equations. Outside the (delta--like) source this {\it not} a vacuum solution of General Relativity. It represents the contribution of the massive, $m=(-\beta)^{-1/2}$, spin--two field associated to the Ricci squared term in the gravitational Lagrangian. The fourth order terms tend to make milder the singularity of the curvature at the point where the particle is located. We generalize this analysis to $D$--dimensions, extended sources, and higher than fourth order theories. We also briefly discuss the scattering of fields by this kind of plane gravitational waves., Comment: 12 pages, REVTEX, Fully revised version. Amplitude of the wave computed. Discussion section added. Figure added. To appear in Phys. Rev. D

We consider the computation of the entanglement entropy in curved backgrounds with event horizons. We use a Hamiltonian approach to the problem and perform numerical computations on a spherical lattice of spacing $a$. We study the cosmological case and make explicit computations for the Friedmann-Robertson-Walker universe. Our results for a massless, minimally coupled scalar field can be summarized by $S_{ent}=0.30 r_H^2/a^2$,which resembles the flat space formula, although here the horizon radius, $r_H$, is time-dependent., Comment: 12 pages, RevTex 3.0, 2 figures as uuencoded compressed Postscript files

We study first order fluctuations of a relativistic membrane in the curved background of a black hole. The zeroth-order solution corresponds to a spherical membrane tightly covering the event horizon. We obtain a massive Klein-Gordon equation for the fluctuations of the membrane's radial coordinate on the 2+1 dimensional world-volume. We finally suggest that quantization of the fluctuations can be related to black hole's mass quantization and the corresponding entropy is computed. This entropy is proportional to the membrane area and is related to the one-loop correction to the thermodynamical entropy $A_H/4$. With regards to the membrane model for describing effectively a quantum black hole, we connect these results with previous work on critical phenomena in black hole thermodynamics., Comment: 13 pages, RevTEX. Paper rewritten; energy levels and entropy recomputed

We deal with the black hole information loss paradox by showing that the stimulated emission component of the black hole radiation contains information about the initial state of the system. The nonlocal behaviour that allows the recovery of information about the matter that falls behind the horizon appears in a natural way. We calculate the expectation value and probability distribution of particles at ${\cal J}^+$ for a non-vacuum initial state. The entropy of the final state is compared to that of a thermal state with the same energy per mode. We find that the information recovered about the initial state increases with the number $r$ of the initially incoming particles, reaching for example over 30\% for $r = 1000$. We point out that recovering information about the initial state, nevertheless, does not automatically imply the purity of the final state., Comment: Latex, 9 pages, 2 figures appended as postscript files