The [N II] Kinematics of R Aquarii

We report a kinematic study of the symbiotic star system R Aqr derived from [N H]lambda 6584 emission observations with a Fabry-Perot imaging spectrometer. The [N II] spatial structure of the R Aqr jet, first observed circa 1977, and surrounding hourglass-shaped nebulosity, due to an explosion approximately 660 years ago, are derived from 41 velocity planes spaced at approximately 12 km/s intervals. Fabry-Perot imagery shows the elliptical nebulosity comprising the waist of the hourglass shell is consistent with a circular ring expanding radially at 55 km/s as seen at an inclination angle, i approximately 70 deg. Fabry-Perot imagery shows the two-sided R Aqr jet is collimated flow in opposite directions. The intensity-velocity structure of the strong NE jet component is shown in contrast to the amorphous SW jet component. We offer a idealized schematic model for the R Aqr jet motion which results in a small-scale helical structure forming around a larger-scale helical path. The implications of such a jet model are discussed. We present a movie showing a side-by-side comparison of the spatial structure of the model and the data as a function of the 41 velocity planes

New, efficient, and accurate high order derivative and dissipation operators satisfying summation by parts, and applications in three-dimensional multi-block evolutions

We construct new, efficient, and accurate high-order finite differencing operators which satisfy summation by parts. Since these operators are not uniquely defined, we consider several optimization criteria: minimizing the bandwidth, the truncation error on the boundary points, the spectral radius, or a combination of these. We examine in detail a set of operators that are up to tenth order accurate in the interior, and we surprisingly find that a combination of these optimizations can improve the operators' spectral radius and accuracy by orders of magnitude in certain cases. We also construct high-order dissipation operators that are compatible with these new finite difference operators and which are semi-definite with respect to the appropriate summation by parts scalar product. We test the stability and accuracy of these new difference and dissipation operators by evolving a three-dimensional scalar wave equation on a spherical domain consisting of seven blocks, each discretized with a structured grid, and connected through penalty boundary conditions.Comment: 16 pages, 9 figures. The files with the coefficients for the derivative and dissipation operators can be accessed by downloading the source code for the document. The files are located in the "coeffs" subdirector

Characteristic extraction in numerical relativity: binary black hole merger waveforms at null infinity

The accurate modeling of gravitational radiation is a key issue for gravitational wave astronomy. As simulation codes reach higher accuracy, systematic errors inherent in current numerical relativity wave-extraction methods become evident, and may lead to a wrong astrophysical interpretation of the data. In this paper, we give a detailed description of the Cauchy-characteristic extraction technique applied to binary black hole inspiral and merger evolutions to obtain gravitational waveforms that are defined unambiguously, that is, at future null infinity. By this method we remove finite-radius approximations and the need to extrapolate data from the near zone. Further, we demonstrate that the method is free of gauge effects and thus is affected only by numerical error. Various consistency checks reveal that energy and angular momentum are conserved to high precision and agree very well with extrapolated data. In addition, we revisit the computation of the gravitational recoil and find that finite radius extrapolation very well approximates the result at \scri. However, the (non-convergent) systematic differences to extrapolated data are of the same order of magnitude as the (convergent) discretisation error of the Cauchy evolution hence highlighting the need for correct wave-extraction.Comment: 41 pages, 8 figures, 2 tables, added references, fixed typos. Version matches published version

Matching post-Newtonian and numerical relativity waveforms: systematic errors and a new phenomenological model for non-precessing black hole binaries

We present a new phenomenological gravitational waveform model for the inspiral and coalescence of non-precessing spinning black hole binaries. Our approach is based on a frequency domain matching of post-Newtonian inspiral waveforms with numerical relativity based binary black hole coalescence waveforms. We quantify the various possible sources of systematic errors that arise in matching post-Newtonian and numerical relativity waveforms, and we use a matching criteria based on minimizing these errors; we find that the dominant source of errors are those in the post-Newtonian waveforms near the merger. An analytical formula for the dominant mode of the gravitational radiation of non-precessing black hole binaries is presented that captures the phenomenology of the hybrid waveforms. Its implementation in the current searches for gravitational waves should allow cross-checks of other inspiral-merger-ringdown waveform families and improve the reach of gravitational wave searches.Comment: 22 pages, 11 figure

Numerical relativity with characteristic evolution, using six angular patches

The characteristic approach to numerical relativity is a useful tool in evolving gravitational systems. In the past this has been implemented using two patches of stereographic angular coordinates. In other applications, a six-patch angular coordinate system has proved effective. Here we investigate the use of a six-patch system in characteristic numerical relativity, by comparing an existing two-patch implementation (using second-order finite differencing throughout) with a new six-patch implementation (using either second- or fourth-order finite differencing for the angular derivatives). We compare these different codes by monitoring the Einstein constraint equations, numerically evaluated independently from the evolution. We find that, compared to the (second-order) two-patch code at equivalent resolutions, the errors of the second-order six-patch code are smaller by a factor of about 2, and the errors of the fourth-order six-patch code are smaller by a factor of nearly 50.Comment: 12 pages, 5 figures, submitted to CQG (special NFNR issue

Initial data transients in binary black hole evolutions

We describe a method for initializing characteristic evolutions of the Einstein equations using a linearized solution corresponding to purely outgoing radiation. This allows for a more consistent application of the characteristic (null cone) techniques for invariantly determining the gravitational radiation content of numerical simulations. In addition, we are able to identify the {\em ingoing} radiation contained in the characteristic initial data, as well as in the initial data of the 3+1 simulation. We find that each component leads to a small but long lasting (several hundred mass scales) transient in the measured outgoing gravitational waves.Comment: 18 pages, 4 figure

Exact boundary conditions in numerical relativity using multiple grids: scalar field tests

Cauchy-Characteristic Matching (CCM), the combination of a central 3+1 Cauchy code with an exterior characteristic code connected across a time-like interface, is a promising technique for the generation and extraction of gravitational waves. While it provides a tool for the exact specification of boundary conditions for the Cauchy evolution, it also allows to follow gravitational radiation all the way to infinity, where it is unambiguously defined. We present a new fourth order accurate finite difference CCM scheme for a first order reduction of the wave equation around a Schwarzschild black hole in axisymmetry. The matching at the interface between the Cauchy and the characteristic regions is done by transfering appropriate characteristic/null variables. Numerical experiments indicate that the algorithm is fourth order convergent. As an application we reproduce the expected late-time tail decay for the scalar field.Comment: 14 pages, 5 figures. Included changes suggested by referee

AMR, stability and higher accuracy

Efforts to achieve better accuracy in numerical relativity have so far focused either on implementing second order accurate adaptive mesh refinement or on defining higher order accurate differences and update schemes. Here, we argue for the combination, that is a higher order accurate adaptive scheme. This combines the power that adaptive gridding techniques provide to resolve fine scales (in addition to a more efficient use of resources) together with the higher accuracy furnished by higher order schemes when the solution is adequately resolved. To define a convenient higher order adaptive mesh refinement scheme, we discuss a few different modifications of the standard, second order accurate approach of Berger and Oliger. Applying each of these methods to a simple model problem, we find these options have unstable modes. However, a novel approach to dealing with the grid boundaries introduced by the adaptivity appears stable and quite promising for the use of high order operators within an adaptive framework

The Samurai Project: verifying the consistency of black-hole-binary waveforms for gravitational-wave detection

We quantify the consistency of numerical-relativity black-hole-binary waveforms for use in gravitational-wave (GW) searches with current and planned ground-based detectors. We compare previously published results for the $(\ell=2,| m | =2)$ mode of the gravitational waves from an equal-mass nonspinning binary, calculated by five numerical codes. We focus on the 1000M (about six orbits, or 12 GW cycles) before the peak of the GW amplitude and the subsequent ringdown. We find that the phase and amplitude agree within each code's uncertainty estimates. The mismatch between the $(\ell=2,| m| =2)$ modes is better than $10^{-3}$ for binary masses above $60 M_{\odot}$ with respect to the Enhanced LIGO detector noise curve, and for masses above $180 M_{\odot}$ with respect to Advanced LIGO, Virgo and Advanced Virgo. Between the waveforms with the best agreement, the mismatch is below $2 \times 10^{-4}$. We find that the waveforms would be indistinguishable in all ground-based detectors (and for the masses we consider) if detected with a signal-to-noise ratio of less than $\approx14$, or less than $\approx25$ in the best cases.Comment: 17 pages, 9 figures. Version accepted by PR

A pilgrimage to gravity on GPUs

In this short review we present the developments over the last 5 decades that have led to the use of Graphics Processing Units (GPUs) for astrophysical simulations. Since the introduction of NVIDIA's Compute Unified Device Architecture (CUDA) in 2007 the GPU has become a valuable tool for N-body simulations and is so popular these days that almost all papers about high precision N-body simulations use methods that are accelerated by GPUs. With the GPU hardware becoming more advanced and being used for more advanced algorithms like gravitational tree-codes we see a bright future for GPU like hardware in computational astrophysics.Comment: To appear in: European Physical Journal "Special Topics" : "Computer Simulations on Graphics Processing Units" . 18 pages, 8 figure