Simulating binary neutron stars: dynamics and gravitational waves

We model two mergers of orbiting binary neutron stars, the first forming a black hole and the second a differentially rotating neutron star. We extract gravitational waveforms in the wave zone. Comparisons to a post-Newtonian analysis allow us to compute the orbital kinematics, including trajectories and orbital eccentricities. We verify our code by evolving single stars and extracting radial perturbative modes, which compare very well to results from perturbation theory. The Einstein equations are solved in a first order reduction of the generalized harmonic formulation, and the fluid equations are solved using a modified convex essentially non-oscillatory method. All calculations are done in three spatial dimensions without symmetry assumptions. We use the \had computational infrastructure for distributed adaptive mesh refinement.Comment: 14 pages, 16 figures. Added one figure from previous version; corrected typo

Relativistic MHD with Adaptive Mesh Refinement

This paper presents a new computer code to solve the general relativistic magnetohydrodynamics (GRMHD) equations using distributed parallel adaptive mesh refinement (AMR). The fluid equations are solved using a finite difference Convex ENO method (CENO) in 3+1 dimensions, and the AMR is Berger-Oliger. Hyperbolic divergence cleaning is used to control the $\nabla\cdot {\bf B}=0$ constraint. We present results from three flat space tests, and examine the accretion of a fluid onto a Schwarzschild black hole, reproducing the Michel solution. The AMR simulations substantially improve performance while reproducing the resolution equivalent unigrid simulation results. Finally, we discuss strong scaling results for parallel unigrid and AMR runs.Comment: 24 pages, 14 figures, 3 table

Energy and Momentum Distributions of the Magnetic Solution to (2+1) Einstein-Maxwell Gravity

We use Moeller's energy-momentum complex in order to explicitly evaluate the energy and momentum density distributions associated with the three-dimensional magnetic solution to the Einstein-Maxwell equations. The magnetic spacetime under consideration is a one-parametric solution describing the distribution of a radial magnetic field in a three-dimensional AdS background, and representing the superposition of the magnetic field with a 2+1 Einstein static gravitational field.Comment: LaTex, 13 pages; v2 clarifying comments and references added, Conclusions improved, to appear in Mod. Phys. Lett.

Relativistic MHD and black hole excision: Formulation and initial tests

A new algorithm for solving the general relativistic MHD equations is described in this paper. We design our scheme to incorporate black hole excision with smooth boundaries, and to simplify solving the combined Einstein and MHD equations with AMR. The fluid equations are solved using a finite difference Convex ENO method. Excision is implemented using overlapping grids. Elliptic and hyperbolic divergence cleaning techniques allow for maximum flexibility in choosing coordinate systems, and we compare both methods for a standard problem. Numerical results of standard test problems are presented in two-dimensional flat space using excision, overlapping grids, and elliptic and hyperbolic divergence cleaning.Comment: 22 pages, 8 figure

An Axisymmetric Gravitational Collapse Code

We present a new numerical code designed to solve the Einstein field equations for axisymmetric spacetimes. The long term goal of this project is to construct a code that will be capable of studying many problems of interest in axisymmetry, including gravitational collapse, critical phenomena, investigations of cosmic censorship, and head-on black hole collisions. Our objective here is to detail the (2+1)+1 formalism we use to arrive at the corresponding system of equations and the numerical methods we use to solve them. We are able to obtain stable evolution, despite the singular nature of the coordinate system on the axis, by enforcing appropriate regularity conditions on all variables and by adding numerical dissipation to hyperbolic equations.Comment: 19 pages, 9 figure

Positron annihilation analysis of nanopores and growth mechanism of oblique angle evaporated TiO2 and SiO2 thin films and multilayers

The nano-porosity embedded into the tilted and separated nanocolumns characteristic of the microstructure of evaporated thin films at oblique angles has been critically assessed by various variants of the positron annihilation spectroscopy. This technique represents a powerful tool for the analysis of porosity, defects and internal interfaces of materials, and has been applied to different as-deposited SiO and TiO thin films as well as SiO/TiO multilayers prepared by electron beam evaporation at 70° and 85° zenithal angles. It is shown that, under same deposition conditions, the concentration of internal nano-pores in SiO is higher than in TiO nanocolumns, while the situation is closer to this latter in TiO/SiO multilayers. These features have been compared with the predictions of a Monte Carlo simulation of the film growth and explained by considering the influence of the chemical composition on the growth mechanism and, ultimately, on the structure of the films

Gravitational collapse of massless scalar field and radiation fluid

Several classes of conformally-flat and spherically symmetric exact solutions to the Einstein field equations coupled with either a massless scalar field or a radiation fluid are given, and their main properties are studied. It is found that some represent the formation of black holes due to the gravitational collapse of the matter fields. When the spacetimes have continuous self-similarity (CSS), the masses of black holes take a scaling form $M_{BH} \propto (P - P^{*})^{\gamma}$, where $\gamma = 0.5$ for massless scalar field and $\gamma = 1$ for radiation fluid. The reasons for the difference between the values of $\gamma$ obtained here and those obtained previously are discussed. When the spacetimes have neither CSS nor DSS (Discrete self-similarity), the masses of black holes always turn on with finite non-zero values.Comment: Two figures have been removed, and the text has been re-written. To appear in Phys. Rev.

Tips for implementing multigrid methods on domains containing holes

As part of our development of a computer code to perform 3D `constrained evolution' of Einstein's equations in 3+1 form, we discuss issues regarding the efficient solution of elliptic equations on domains containing holes (i.e., excised regions), via the multigrid method. We consider as a test case the Poisson equation with a nonlinear term added, as a means of illustrating the principles involved, and move to a "real world" 3-dimensional problem which is the solution of the conformally flat Hamiltonian constraint with Dirichlet and Robin boundary conditions. Using our vertex-centered multigrid code, we demonstrate globally second-order-accurate solutions of elliptic equations over domains containing holes, in two and three spatial dimensions. Keys to the success of this method are the choice of the restriction operator near the holes and definition of the location of the inner boundary. In some cases (e.g. two holes in two dimensions), more and more smoothing may be required as the mesh spacing decreases to zero; however for the resolutions currently of interest to many numerical relativists, it is feasible to maintain second order convergence by concentrating smoothing (spatially) where it is needed most. This paper, and our publicly available source code, are intended to serve as semi-pedagogical guides for those who may wish to implement similar schemes.Comment: 18 pages, 11 figures, LaTeX. Added clarifications and references re. scope of paper, mathematical foundations, relevance of work. Accepted for publication in Classical & Quantum Gravit

Spherical and planar three-dimensional anti-de Sitter black holes

The technique of dimensional reduction was used in a recent paper (Zanchin et al, Phys. Rev. D66, 064022,(2002)) where a three-dimensional (3D) Einstein-Maxwell-Dilaton theory was built from the usual four-dimensional (4D) Einstein-Maxwell-Hilbert action for general relativity. Starting from a class of 4D toroidal black holes in asymptotically anti-de Sitter (AdS) spacetimes several 3D black holes were obtained and studied in such a context. In the present work we choose a particular case of the 3D action which presents Maxwell field, dilaton field and an extra scalar field, besides gravity field and a negative cosmological constant, and obtain new 3D static black hole solutions whose horizons may have spherical or planar topology. We show that there is a 3D static spherically symmetric solution analogous to the 4D Reissner-Nordstr\"om-AdS black hole, and obtain other new 3D black holes with planar topology. From the static spherical solutions, new rotating 3D black holes are also obtained and analyzed in some detail.Comment: 27 pages, uses "iopclass" files (Latex2e