Numerical Relativity

The tetrad-based equations for vacuum gravity published by Estabrook, Robinson, and Wahlquist are simplified and adapted for numerical relativity. We show that the evolution equations as partial differential equations for the Ricci rotation coefficients constitute a rather simple first-order symmetrizable hyperbolic system, not only for the Nester gauge condition on the acceleration and angular velocity of the tetrad frames considered by Estabrook et al., but also for the Lorentz gauge condition of van Putten and Eardley and for a fixed gauge condition. We introduce a lapse function and a shift vector to allow general coordinate evolution relative to the timelike congruence defined by the tetrad vector field.

We derive the Teukolsky equation for perturbations of a Kerr spacetime when the spacetime metric is written in either ingoing or outgoing Kerr-Schild form. We also write explicit formulae for setting up the initial data for the Teukolsky equation in the time domain in terms of a three metric and an extrinsic curvature. The motivation of this work is to have in place a formalism to study the evolution in the "close limit" of two recently proposed solutions to the initial value problem in general relativity that are based on Kerr-Schild slicings. A perturbative formalism in horizon penetrating coordinates is also very desirable in connection with numerical relativity simulations using black hole "excision".

Black objects in higher dimensional space-times have a remarkably richer structure than their four dimensional counterparts. They appear in a variety of configurations (e.g. black holes, black branes, black rings, black Saturns), and display complex stability phase diagrams. They might also play a key role in high energy physics: for energies above the fundamental Planck scale, gravity is the dominant interaction which, together with the hoop-conjecture, implies that the trans-Planckian scattering of point particles should be well described by black hole scattering. Higher dimensional scenarios with a fundamental Planck scale of the order of TeV predict, therefore, black hole production at the LHC, as well as in future colliders with yet higher energies. In this setting, accurate predictions for the production cross-section and energy loss (through gravitational radiation) in the formation of black holes in parton-parton collisions is crucial for accurate phenomenological modelling in Monte Carlo event generators.

La primera detección directa de ondas gravitacionales, obtenida recientemente en los detectores LIGO, abrió una nueva ventana al Universo. Adicionalmente, dado el rol que juega la relatividad numérica en la elaboración de catálogos de perfiles de ondas, este hallazgo propicia un empuje importante al modelado computacional de la radiación gravitacional. En este sentido, teniendo presente la importancia de simular sistemas aislados en esta área, recientemente ha habido un renovado interés en resolver numéricamente el problema de valores iniciales hiperboloidal para las ecuaciones de campo de Einstein, en el que la frontera exterior de la malla numérica se ubica en el infinito nulo futuro (scri+). En este trabajo implementamos numéricamente el enfoque tetradial BSB presentado en [J.M. Bardeen, O. Sarbach y L.T. Buchman, Phys. Rev. D 83, 104045 (2011)], para el caso de un campo escalar autogravitante, minimamente acoplado, esféricamente simétrico. Debido a la dificultad de esta tarea, introduciremos muchos de los ingredientes que necesitamos con tres pruebas numéricas previas. Partimos resolviendo la ecuación de onda en 1+1 dimensiones, luego resolvemos esta ecuación en el contexto de la relatividad espacial sobre un espacio-tiempo de Minkowski, para posteriormente resolver la ecuación de onda esféricamente simétrica de un campo escalar no masivo sobre un fondo de Schwarzschild. En las dos últimas pruebas implementamos métodos conformes, utilizando foliaciones en curvas (o bien hipersuperficies) tipo espacio, de curvatura extrínseca media constante (CMC) positiva que intersectan scri+, junto con compactificaciones que nos permiten mapear nuestro espacio-tiempo infinito a una variedad finita. Únicamente despues de realizar estas pruebas, es que nos volcamos al objetivo principal de este trabajo. Estudiaremos en detalle el enfoque BSB, para el caso general sin simetrías, y luego reduciremos el escenario a simetría esférica y trabajaremos en la implementación numérica incluyendo el ya mencionado campo escalar autogravitante. En particular, el sistema de evolución se reduce a un sistema de ecuaciones simétrico-hiperbólico de primer orden, regular, para el campo escalar conformemente reescalado, el cual se acopla a un conjunto de constricciones elípticas singulares para los coeficientes métricos. Aquí mostraremos como resolver este sistema, basándonos en aproximaciones de diferencias finitas, obteniendo evoluciones numéricas estables para configuraciones iniciales de agujeros negros rodeados por un cascarón esférico de campo escalar, resultando que una parte de este último se dispersa hacia scri+ y la otra parte es acretada por el agujero negro. Como una prueba no trivial, estudiamos el decaimiento de cola del campo escalar a lo largo de diferentes líneas de mundo: una situada a lo largo del horizonte (aparente), otra correspondiente a un observador tipo tiempo a un radio finito fuera del horizonte, y una última a lo largo de scri+. Nuestros resultados coinciden con la ley usual de decaimientos polinomiales predicha por la teoría linealizada, así como por trabajos previos.

We describe the burgeoning field of numerical relativity, which aims to solve Einstein's equations of general relativity numerically. The field presents many questions that may interest numerical analysts, especially problems related to nonlinear partial differential equations: elliptic systems, hyperbolic systems, and mixed systems. There are many novel features, such as dealing with boundaries when black holes are excised from the computational domain, or how to even pose the problem computationally when the coordinates must be determined during the evolution from initial data. The most important unsolved problem is that there is no known general 3-dimensional algorithm that can evolve Einstein's equations with black holes that is stable. This review is meant to be an introduction that will enable numerical analysts and other computational scientists to enter the field. No previous knowledge of special or general relativity is assumed.

We consider the Universe deep inside the cell of uniformity. At these scales, the Universe is filled with inhomogeneously distributed discrete structures (galaxies, groups and clusters of galaxies), which perturb the background Friedmann model. Here, the mechanical approach (Eingorn & Zhuk, 2012) is the most appropriate to describe the dynamics of the inhomogeneities which is defined, on the one hand, by gravitational potentials of inhomogeneities and, on the other hand, by the cosmological expansion of the Universe. In this paper, we present additional arguments in favor of this approach. First, we estimate the size of the cell of uniformity. With the help of the standard methods of statistical physics and for the galaxies of the type of the Milky Way and Andromeda, we get that it is of the order of 190 Mpc which is rather close to observations. Then, we show that the nonrelativistic approximation (with respect to the peculiar velocities) is valid for $z \lesssim 10$, i.e. approximately for 13 billion years from the present moment. We consider scalar perturbations and, within the $\Lambda$CDM model, justify the main equations. Moreover, we demonstrate that radiation can be naturally incorporated into our scheme. This emphasizes the viability of our approach. This approach gives a possibility to analyze different cosmological models and compare them with the observable Universe. For example, we indicate some problematic aspects of the spatially flat models. Such models require a rather specific distribution of the inhomogeneities to get a finite potential at any points outside gravitating masses. We also criticize the application of the Schwarzschild-de Sitter solution to the description of the motion of test bodies on the cosmological background.

We give an introduction to the Einstein Toolkit, a mature, open-source computational infrastructure for numerical relativity based on the Cactus Framework, for the target group of new users. This toolkit is composed of several different modules, is developed by researchers from different institutions throughout the world and is in active continuous development. Documentation for the toolkit and its several modules is often scattered across different locations, a difficulty new users may at times have to struggle with. Scientific papers exist describing the toolkit and its methods in detail, but they might be overwhelming at first. With these lecture notes we hope to provide an initial overview for new users. We cover how to obtain, compile and run the toolkit, and give an overview of some of the tools and modules provided with it.

Assessing the quality of groundwater is important to ensure sustainable safe use of these resources. However, describing the overall water quality condition is difficult due to the spatial variability of multiple contaminants and the wide range of indicators (chemical, physical and biological) that could be measured. This contribution proposes a GIS-based groundwater quality index (GQI) which synthesizes different available water quality data (e.g., Cl − , Na + , Ca 2+ ) by indexing them numerically relative to the World Health Organization (WHO) standards. Also, introduces an objective procedure to select the optimum parameters to compute the GQI, incorporates the aspect of temporal variation to address the degree of water use sustainability and tests the sensitivity of the proposed model.

The isolated horizon formalism recently introduced by Ashtekar et al. aims at providing a quasi-local concept of a black hole in equilibrium in an otherwise possibly dynamical spacetime. In this formalism, a hierarchy of geometrical structures is constructed on a null hypersurface. On the other side, the 3+13+1 formulation of general relativity provides a powerful setting for studying the spacetime dynamics, in particular gravitational radiation from black hole systems. Here we revisit the kinematics and dynamics of null hypersurfaces by making use of some 3+13+1 slicing of spacetime. In particular, the additional structures induced on null hypersurfaces by the 3+13+1 slicing permit a natural extension to the full spacetime of geometrical quantities defined on the null hypersurface. This four-dimensional point of view facilitates the link between the null and spatial geometries. We proceed by reformulating the isolated horizon structure in this framework. We also reformulate previous works, such as Damour's black hole mechanics, and make the link with a previous 3+13+1 approach of black hole horizon, namely the membrane paradigm  . We explicit all geometrical objects in terms of 3+13+1 quantities, putting a special emphasis on the conformal 3+13+1 formulation. This is in particular relevant for the initial data problem of black hole spacetimes for numerical relativity. Illustrative examples are provided by considering various slicings of Schwarzschild and Kerr spacetimes.

Grid Portals, based on standard web technologies, are emerging as important and useful user interfaces to computational and data Grids. Grid Portals enable Virtual Organizations, comprised of distributed researchers to collaborate and access resources more efficiently and seamlessly. The Astrophysics Simulation Collaboratory (ASC) Grid Portal provides a framework to enable researchers in the field of numerical relativity to study astrophysical phenomenon by making use of the Cactus computational toolkit. We examine user requirements and describe the design and implementation of the ASC Grid Portal.

The ‘anomalous perihelion precession’ of Mercury, announced by Le Verrier in 1859, was a highly controversial topic for more than half a century and invoked many alternative theories until 1916, when Einstein presented his theory of general relativity as an alternative theory of gravitation and showed perihelion precession to be one of its potential manifestations. As perihelion precession was a directly derived result of the full General Theory and not just the Equivalence Principle, Einstein viewed it as the most critical test of his theory. This paper presents the computed value of the anomalous perihelion precession of Mercury's orbit using a new relativistic simulation model that employs a simple transformation factor for mass and time, proposed in an earlier paper. This computed value compares well with the prediction of general relativity and is, also, in complete agreement with the observed value within its range of uncertainty. No general relativistic equations have bee...

Perturbations of stars and black holes have been one of the main topics of relativistic astrophysics for the last few decades. They are of particular importance today, because of their relevance to gravitational wave astronomy. In this review we present the theory of quasi-normal modes of compact objects from both the mathematical and astrophysical points of view. The discussion includes perturbations of black holes (Schwarzschild, Reissner-Nordström, Kerr and Kerr-Newman) and relativistic stars (nonrotating and slowly-rotating). The properties of the various families of quasi-normal modes are described, and numerical techniques for calculating quasi-normal modes reviewed. The successes, as well as the limits, of perturbation theory are presented, and its role in the emerging era of numerical relativity and supercomputers is discussed.

Improvements in the performance of processors and networks make it both feasible and interesting to treat collections of workstations, servers, clusters, and supercomputers as integrated computational resources, or Grids. However, the highly heterogeneous and dynamic nature of such Grids can make application development difficult. Here we describe an architecture and prototype implementation for a Grid-enabled computational framework based on Cactus, the MPICH-G2 Grid-enabled message-passing library, and a variety of specialized features to support efficient execution in Grid environments. We have used this framework to perform record-setting computations in numerical relativity, running across four supercomputers and achieving scaling of 88% (1140 CPU's) and 63% (1500 CPUs). The problem size we were able to compute was about five times larger than any other previous run. Further, we introduce and demonstrate adaptive methods that automatically adjust computational parameters during run time, to increase dramatically the efficiency of a distributed Grid simulation, without modification of the application and without any knowledge of the underlying network connecting the distributed computers.

We developed a new method for time-domain signal processing based on deep (dilated) convolutional neural networks which can rapidly identify and extract signals much weaker than the background noise. We applied this method for gravitational wave analysis and demonstrated that it significantly outperforms existing matched-filtering and conventional machine learning techniques and also extends the range of gravitational waves that can be detected by advanced LIGO, allowing real-time processing of raw big data with minimal resources. This initiates a new paradigm for research which uses massively-parallel simulations to train artificial intelligence algorithms that exploit emerging hardware architectures such as deep-learning-optimized GPUs. This approach provides a natural framework to enable coincident detection campaigns of gravitational waves sources and their electromagnetic counterparts.

This article reviews some aspects in the current relationship between mathematical and numerical General Relativity. Focus is placed on the description of isolated systems, with a particular emphasis on recent developments in the study of black holes. Ideas concerning asymptotic flatness, the initial value problem, the constraint equations, evolution formalisms, geometric inequalities and quasi-local black hole horizons are discussed on the light of the interaction between numerical and mathematical relativists.

This is the second in a series of papers on the construction and validation of a three-dimensional code for the solution of the coupled system of the Einstein equations and of the general relativistic hydrodynamic equations, and on the application of this code to problems in general relativistic astrophysics. In particular, we report on the accuracy of our code in the long-term dynamical evolution of relativistic stars and on some new physics results obtained in the process of code testing. The following aspects of our code have been validated: the generation of initial data representing perturbed general relativistic polytropic models (both rotating and non-rotating), the long-term evolution of relativistic stellar models, and the coupling of our evolution code to analysis modules providing, for instance, the detection of apparent horizons or the extraction of gravitational waveforms. The tests involve single non-rotating stars in stable equilibrium, non-rotating stars undergoing radial and quadrupolar oscillations, non-rotating stars on the unstable branch of the equilibrium configurations migrating to the stable branch, non-rotating stars undergoing gravitational collapse to a black hole, and rapidly rotating stars in stable equilibrium and undergoing quasi-radial oscillations. We have carried out evolutions in full general relativity and compared the results to those obtained either with perturbation techniques, or with lower dimensional numerical codes, or in the Cowling approximation (in which all the perturbations of the spacetime are neglected). In all cases an excellent agreement has been found. The numerical evolutions have been carried out using different types of polytropic equations of state using either the rest-mass density only, or the rest-mass density and the internal energy as independent variables. New variants of the spacetime evolution and new high resolution shock capturing (HRSC) treatments based on Riemann solvers and slope limiters have been implemented and the results compared with those obtained from previous methods. In particular, we have found the "monotonized central differencing" (MC) limiter to be particularly effective in evolving the relativistic stellar models considered. Finally, we have obtained the first eigenfrequencies of rotating stars in full general relativity and rapid rotation. A long standing problem, such frequencies have not been obtained by other methods. Overall, and to the best of our knowledge, the results presented in this paper represent the most accurate long-term three-dimensional evolutions of relativistic stars available to date.

Initial data are the starting point for any numerical simulation. In the case of numerical relativity, Einstein's equations constrain our choices of these initial data. We will examine several of the formalisms used for specifying Cauchy initial data in the 3 + 1 decomposition of Einstein's equations. We will then explore how these formalisms have been used in constructing initial data for spacetimes containing black holes and neutron stars. In the topics discussed, emphasis is placed on those issues that are important for obtaining astrophysically realistic initial data for compact binary coalescence.

We describe a simple implementation of black hole excision in 3+1 numerical relativity. We apply this technique to a Schwarzschild black hole with octant symmetry in Eddington-Finkelstein coordinates and show how one can obtain accurate, long-term stable numerical evolutions. 04.25.Dm, 04.30.Db, 95.30.Sf, 97.60.Lf

The periodic standing wave method studies circular orbits of compact objects coupled to helically symmetric standing wave gravitational fields. From this solution an approximation is extracted for the strong field, slowly inspiralling motion of black holes and binary stars. Previous work on this model has dealt with nonlinear scalar models, and with linearized general relativity. Here we present the results of the method for the post-Minkowski (PM) approximation to general relativity, the first step beyond linearized gravity. We compute the PM approximation in two ways: first, via the standard approach of computing linearized gravitational fields and constructing from them quadratic driving sources for second-order fields, and second, by solving the second-order equations as an "exact" nonlinear system. The results of these computations have two distinct applications: (i) The computational infrastructure for the "exact" PM solution will be directly applicable to full general relativity. (ii) The results will allow us to begin supplying initial data to collaborators running general relativistic evolution codes.

We apply recently developed techniques for simulations of moving black holes to study dynamics and radiation generation in the last few orbits and merger of a binary black hole system. Our analysis produces a consistent picture from the gravitational wave forms and dynamical black hole trajectories for a set of simulations with black holes beginning on circular-orbit trajectories at a variety of initial separations. We find profound agreement at the level of 1% among the simulations for the last orbit, merger and ringdown, resulting in a final black hole with spin parameter a/m = 0.69. Consequently, we are confident that this part of our waveform result accurately represents the predictions from Einstein's General Relativity for the final burst of gravitational radiation resulting from the merger of an astrophysical system of equal-mass non-spinning black holes. We also find good agreement at a level of roughly 10% for the radiation generated in the preceding few orbits.

We perform general relativistic, magnetohydrodynamic (GRMHD) simulations of merging binary neutron stars incorporating neutrino transport and magnetic fields. Our new radiative transport module for neutrinos adopts a general relativistic, truncated-moment (M1) formalism. The binaries consist of two identical, irrotational stars modeled by the SLy nuclear equation of state (EOS). They are initially in quasicircular orbit and threaded with a poloidal magnetic field that extends from the stellar interior into the exterior, as in typical pulsars. We insert neutrino processes shortly after the merger and focus on the role of neutrinos in launching a jet following the collapse of the hypermassive neutron star (HMNS) remnant to a spinning black hole (BH). We treat two microphysical versions: one (a "warm-up") evolving a single neutrino species and considering only charged-current processes, and the other evolving three species (νe,νe, νx) and related processes. We trace the evolution until the system reaches a quasiequilibrium state after BH formation. We find that the BH + disk remnant eventually launches an incipient jet. The electromagnetic Poynting luminosity is ∼ 10 53 erg s −1 , consistent with that of typical short gamma-ray bursts (sGRBs). The effect of neutrino cooling shortens the lifetime of the HMNS, and lowers the amplitude of the major peak of the gravitational wave (GW) power spectrum somewhat. After BH formation, neutrinos help clear out the matter near the BH poles, resulting in lower baryon-loaded surrounding debris. The neutrino luminosity resides in the range ∼ 10 52−53 erg s −1 once quasiequilibrium is achieved. Comparing with the neutrino-free models, we observe that the inclusion of neutrinos yields similar ejecta masses and is inefficient in carrying off additional angular momentum.

Black hole-neutron star mergers resulting in the disruption of the neutron star and the formation of an accretion disk and/or the ejection of unbound material are prime candidates for the joint detection of gravitational-wave and electromagnetic signals when the next generation of gravitational-wave detectors comes online. However, the disruption of the neutron star and the properties of the post-merger remnant are very sensitive to the parameters of the binary. In this paper, we study the impact of the radius of the neutron star and the alignment of the black hole spin for systems within the range of mass ratio currently deemed most likely for field binaries (M_BH ~ 7 M_NS) and for black hole spins large enough for the neutron star to disrupt (J/M^2=0.9). We find that: (i) In this regime, the merger is particularly sensitive to the radius of the neutron star, with remnant masses varying from 0.3M_NS to 0.1M_NS for changes of only 2 km in the NS radius; (ii) 0.01-0.05M_sun of unbound material can be ejected with kinetic energy >10^51 ergs, a significant increase compared to low mass ratio, low spin binaries. This ejecta could power detectable optical and radio afterglows. (iii) Only a small fraction (<3%) of the Advanced LIGO events in this parameter range have gravitational-wave signals which could offer constraints on the equation of state of the neutron star. (iv) A misaligned black hole spin works against disk formation, with less neutron star material remaining outside of the black hole after merger, and a larger fraction of that material remaining in the tidal tail instead of the forming accretion disk. (v) Large kicks (v>300 km/s) can be given to the final black hole as a result of a precessing BHNS merger, when the disruption of the neutron star occurs just outside or within the innermost stable spherical orbit.

Collapsing supermassive stars (SMSs) with masses M 10 4−6 M have long been speculated to be the seeds that can grow and become supermassive black holes (SMBHs). We previously performed general rel-ativistic magnetohydrodynamic (GRMHD) simulations of marginally stable Γ = 4/3 polytropes uniformly rotating at the mass-shedding limit and endowed initially with a dynamically unimportant dipole magnetic field to model the direct collapse of SMSs. These configurations are supported entirely by thermal radiation pressure and reliably model SMSs with M 10 6 M. We found that around 90% of the initial stellar mass forms a spinning black hole (BH) remnant surrounded by a massive, hot, magnetized torus, which eventually launches a magnetically-driven jet. SMSs could be therefore sources of ultra-long gamma-ray bursts (ULGRBs). Here we perform GRMHD simulations of Γ 4/3, polytropes to account for the perturbative role of gas pressure in SMSs with M 10 6 M. We also consider different initial stellar rotation profiles. The stars are initially seeded with a dynamically weak dipole magnetic field that is either confined to the stellar interior or extended from its interior into the stellar exterior. We calculate the gravitational wave burst signal for the different cases. We find that the mass of the black hole remnant is 90%−99% of the initial stellar mass, depending sharply on Γ−4/3 as well as on the initial stellar rotation profile. After t ∼ 250 − 550M ≈ 1 − 2 × 10 3 (M/10 6 M)s following the appearance of the BH horizon, an incipient jet is launched and it lasts for ∼ 10 4 − 10 5 (M/10 6 M)s, consistent with the duration of long gamma-ray bursts. Our numerical results suggest that the Blandford-Znajek mechanism powers the incipient jet. They are also in rough agreement with our recently proposed universal model that estimates accretion rates and electromagnetic (Poynting) luminosities that characterize magnetized BH-disk remnant systems that launch a jet. This model helps explain why the outgoing electromagnetic luminosities computed for vastly different BH-disk formation scenarios all reside within a narrow range (∼ 10 52±1 erg s −1), roughly independent of M .

We apply the postquasistatic approximation, an iterative method for the evolution of selfgravitating spheres of matter, to study the evolution of anisotropic non-adiabatic radiating and dissipative distributions in General Relativity. Dissipation is described by viscosity and free-streaming radiation, assuming an equation of state to model anisotropy induced by the shear viscosity. We match the interior solution, in non-comoving coordinates, with the Vaidya exterior solution. Two simple models are presented, based on the Schwarzschild and Tolman VI solutions, in the nonadiabatic and adiabatic limit. In both cases the eventual collapse or expansion of the distribution is mainly controlled by the anisotropy induced by the viscosity. PACS numbers: 04.25.-g,04.25.D-,0.40.-b

Currently there is considerable interest in making use of many-core processor architectures, such as Nvidia and AMD graphics processing units (GPUs) for scientific computing. In this work we explore the use of the Open Computing Language (OpenCL) for a typical Numerical Relativity application: a time-domain Teukolsky equation solver (a linear, hyperbolic, partial differential equation solver using finite-differencing). OpenCL is the only vendor-agnostic and multi-platform parallel computing framework that has been adopted by all major processor vendors. Therefore, it allows us to write portable source-code and run it on a wide variety of compute hardware and perform meaningful comparisons. The outcome of our experimentation suggests that it is relatively straightforward to obtain order-of-magnitude gains in overall application performance by making use of many-core GPUs over multi-core CPUs and this fact is largely independent of the specific hardware architecture and vendor. We also observe that a single high-end GPU can match the performance of a small-sized, message-passing based CPU cluster.

In this paper we study scalar perturbations of the metric for nonlinear $f(R)$ models. We consider the Universe at the late stage of its evolution and deep inside the cell of uniformity. We investigate the astrophysical approach in the case of Minkowski spacetime background and two cases in the cosmological approach, the large scalaron mass approximation and the quasi-static approximation, getting explicit expressions for scalar perturbations for both these cases. In the most interesting quasi-static approximation, the scalar perturbation functions depend on both the nonlinearity function $f(R)$ and the scale factor $a$. Hence, we can study the dynamical behavior of the inhomogeneities (e.g., galaxies and dwarf galaxies) including into consideration their gravitational attraction and the cosmological expansion, and also taking into account the effects of nonlinearity. Our investigation is valid for functions $f(R)$ which have stable de Sitter points in future with respect to the present time, that is typical for the most popular $f(R)$ models.

Numerical relativity has come a long way in the last three decades and is now reaching a state of maturity. We are gaining a deeper understanding of the fundamental theoretical issues related to the field, from the well posedness of the Cauchy problem, to better gauge conditions, improved boundary treatment, and more realistic initial data. There has also been important work both in numerical methods and software engineering. All these developments have come together to allow the construction of several advanced fully three-dimensional codes capable of dealing with both matter and black holes. In this manuscript I make a brief review the current status of the field.

Recent demonstrations of unexcised black holes traversing across computational grids represent a significant advance in numerical relativity. Stable and accurate simulations of multiple orbits, and their radiated waves, result. This capability is critically undergirded by a careful choice of gauge. Here we present analytic considerations which suggest certain gauge choices, and numerically demonstrate their efficacy in evolving a single moving puncture black hole.

This is the second in a series of papers on the construction and validation of a three-dimensional code for the solution of the coupled system of the Einstein equations and of the general relativistic hydrodynamic equations, and on the application of this code to problems in general relativistic astrophysics. In particular, we report on the accuracy of our code in the long-term dynamical evolution of relativistic stars and on some new physics results obtained in the process of code testing. The tests involve single non-rotating stars in stable equilibrium, non-rotating stars undergoing radial and quadrupolar oscillations, non-rotating stars on the unstable branch of the equilibrium configurations migrating to the stable branch, non-rotating stars undergoing gravitational collapse to a black hole, and rapidly rotating stars in stable equilibrium and undergoing quasi-radial oscillations. The numerical evolutions have been carried out in full general relativity using different types of polytropic equations of state using either the rest-mass density only, or the rest-mass density and the internal energy as independent variables. New variants of the spacetime evolution and new high resolution shock capturing (HRSC) treatments based on Riemann solvers and slope limiters have been implemented and the results compared with those obtained from previous methods. Finally, we have obtained the first eigenfrequencies of rotating stars in full general relativity and rapid rotation. A long standing problem, such frequencies have not been obtained by other methods. Overall, and to the best of our knowledge, the results presented in this paper represent the most accurate long-term three-dimensional evolutions of relativistic stars available to date.

We present the first numerical simulations of an initially nonspinning black-hole binary with mass ratio as large as 10∶1 in full general relativity. The binary completes approximately three orbits prior to merger and radiates (0.415±0.017)% of the total energy and (12.48±0.62)% of the initial angular momentum in the form of gravitational waves. The single black hole resulting from the merger acquires a kick of (66.7±3.3)km/s relative to the original center of mass frame. The resulting gravitational waveforms are used to validate existing formulas for the recoil, final spin, and radiated energy over a wider range of the symmetric mass ratio parameter η=M1M2/(M1+M2)2 than previously possible. The contributions of ℓ>2 multipoles are found to visibly influence the gravitational wave signal obtained at fixed inclination angles.

White dwarf-neutron star binaries generate detectable gravitational radiation. We construct Newtonian equilibrium models of corotational white dwarf-neutron star (WDNS) binaries in circular orbit and find that these models terminate at the Roche limit. At this point the binary will undergo either stable mass transfer (SMT) and evolve on a secular time scale, or unstable mass transfer (UMT), which results in the tidal disruption of the WD. The path a given binary will follow depends primarily on its mass ratio. We analyze the fate of known WDNS binaries and use population synthesis results to estimate the number of LISA-resolved galactic binaries that will undergo either SMT or UMT. We model the quasistationary SMT epoch by solving a set of simple ordinary differential equations and compute the corresponding gravitational waveforms. Finally, we discuss in general terms the possible fate of binaries that undergo UMT and construct approximate Newtonian equilibrium configurations of merged WDNS remnants. We use these configurations to assess plausible outcomes of our future, fully relativistic simulations of these systems. If sufficient WD debris lands on the NS, the remnant may collapse, whereby the gravitational waves from the inspiral, merger, and collapse phases will sweep from LISA through LIGO frequency bands. If the debris forms a disk about the NS, it may fragment and form planets.

Einstein's field equations in general relativity admit a variety of solutions with spacetime singularities. Numerical relativity has recently revealed the properties of somewhat generic spacetime singularities. It has been found that in a variety of systems self-similar solutions can describe asymptotic or intermediate behaviour of more general solutions in an approach to singularities. The typical example is the convergence to an attractor self-similar solution in gravitational collapse. This is closely related to the cosmic censorship violation in the spherically symmetric collapse of a perfect fluid. The self-similar solution also plays an important role in critical phenomena in gravitational collapse. The critical phenomena are understood as the intermediate behaviour around a critical self-similar solution. We see that the convergence and critical phenomena are understood in a unified manner in terms of attractors of codimension zero and one, respectively, in renormalisation group flow.

This letter describes the implementation of second-order one-way wave-equation absorbing boundary conditions (ABCs) in two unconditionally stable finite-difference time-domain (FDTD) methods-namely the locally one-dimensional (LOD)-and the alternating-direction implicit (ADI)-FDTD methods. The Higdon second-order absorbing operator is discretized in the same way as when it is used with the conventional FDTD method. The resulting discrete expression is directly applied to the electric field in each time substep of the LOD-and the ADI-FDTD methods. Numerical examples are given to illustrate the validity of the proposed approach. Index Terms-Alternating-direction implicit finite-difference time-domain (ADI-FDTD) method, locally one-dimensional FDTD (LOD-FDTD) method, one-way wave-equation absorbing boundary condition.

The present work shows how black hole information, such as mass and angular momentum, can be extracted via a new theoretical and computational analysis of the hole’s so called ringdown radiation, which is a gravitational wave emitted by a freely oscillating black hole.

This analysis fits within the broader field of numerical relativity whose treatment requires fast and robust computational algorithms. A new spectral algorithm is presented in this work that solves the spin weight spheroidal equation, describing angular perturbations of a black hole. We then use it together with Leaver’s continued fraction method to calculate quasinormal modes (QNM) of a Kerr black hole for multipole number l = 2 through 12 and first 8 overtones with  \epsilon= 10^(−12) accuracy.

The computed QNM spectrum is then used to represent the ringdown signal of a black hole, which is a signature radiation that corresponds to the hole’s parameters. According to the no-hair theorem these parameters are limited to the mass, angular momentum and charge. We assume the ringing black hole is electrically neutral so only the information about its mass and angular momentum is present in the ringdown signal. The multimode nonlinear fitting formalism, that we developed, allows us to extract that information from the black hole.

We present a new evolution system for Einstein's field equations which is based on tetrad fields and conformally compactified hyperboloidal spatial hypersurfaces which reach future null infinity. The boost freedom in the choice of the tetrad is fixed by requiring that its timelike leg be orthogonal to the foliation, which consists of constant mean curvature slices. The rotational freedom in the tetrad is fixed by the 3D Nester gauge. With these conditions, the field equations reduce naturally to a first-order constrained symmetric hyperbolic evolution system which is coupled to elliptic equations for the gauge variables. The conformally rescaled equations are given explicitly, and their regularity at future null infinity is discussed. Our formulation is potentially useful for high accuracy numerical modeling of gravitational radiation emitted by inspiraling and merging black hole binaries and other highly relativistic isolated systems.

Numerical simulations of 15 orbits of an equal-mass binary black-hole system are presented. Gravitational waveforms from these simulations, covering more than 30 cycles and ending about 1.5 cycles before merger, are compared with those from quasicircular zero-spin post-Newtonian (PN) formulae. The cumulative phase uncertainty of these comparisons is about 0.05 radians, dominated by effects arising from the small residual spins of the black holes and the small residual orbital eccentricity in the simulations. Matching numerical results to PN waveforms early in the run yields excellent agreement (within 0.05 radians) over the first 15 cycles, thus validating the numerical simulation and establishing a regime where PN theory is accurate. In the last 15 cycles to merger, however, generic time-domain Taylor approximants build up phase differences of several radians. But, apparently by coincidence, one specific post-Newtonian approximant, TaylorT4 at 3.5PN order, agrees much better with the numerical simulations, with accumulated phase differences of less than 0.05 radians over the 30-cycle waveform. Gravitational-wave amplitude comparisons are also done between numerical simulations and post-Newtonian, and the agreement depends on the post-Newtonian order of the amplitude expansion: the amplitude difference is about 6%-7% for zeroth order and becomes smaller for increasing order. A newly derived 3.0PN amplitude correction improves agreement significantly ( < 1% amplitude difference throughout most of the run, increasing to 4% near merger) over the previously known 2.5PN amplitude terms.

Assessing the quality of groundwater is important to ensure sustainable safe use of these resources. However, describing the overall water quality condition is difficult due to the spatial variability of multiple contaminants and the wide range of indicators (chemical, physical and biological) that could be measured. This contribution proposes a GIS-based groundwater quality index (GQI) which synthesizes different available water quality data (e.g., Cl−, Na+, Ca2+) by indexing them numerically relative to the World Health Organization (WHO) standards. Also, introduces an objective procedure to select the optimum parameters to compute the GQI, incorporates the aspect of temporal variation to address the degree of water use sustainability and tests the sensitivity of the proposed model. The GQI indicated that the groundwater quality in the Nasuno basin, Tochigi Prefecture, Japan, is generally high (GQI V, 15–30%) compared to the middle part (V,

We present a new method for the calculation of black hole perturbations induced by extended sources in which the solution of the nonlinear hydrodynamics equations is coupled to a perturbative method based on Regge-Wheeler/Zerilli and Bardeen-Press-Teukolsky equations when these are solved in the frequency domain. In contrast to alternative methods in the time domain which may be unstable for rotating black-hole spacetimes, this approach is expected to be stable as long as an accurate evolution of the matter sources is possible. Hence, it could be used under generic conditions and also with sources coming from three-dimensional numerical relativity codes. As an application of this method we compute the gravitational radiation from an oscillating high-density torus orbiting around a Schwarzschild black hole and show that our method is remarkably accurate, capturing both the basic quadrupolar emission of the torus and the excited emission of the black hole.

We present improved post-Newtonian-inspired initial data for non-spinning black-hole binaries, suitable for numerical evolution with punctures. We revisit the work of Tichy et al. [W. Tichy, B. Brügmann, M. Campanelli, and P. Diener, Phys. Rev. D 67, 064008 (2003)], explicitly calculating the remaining integral terms. These terms improve accuracy in the far zone and, for the first time, include realistic gravitational waves in the initial data. We investigate the behavior of these data both at the center of mass and in the far zone, demonstrating agreement of the transversetraceless parts of the new metric with quadrupole-approximation waveforms. These data can be used for numerical evolutions, enabling a direct connection between the merger waveforms and the post-Newtonian inspiral waveforms.