BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Zachariah B. Etienne

Zachariah B. Etienne appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.SRastro-ph.COastro-ph.GAphysics.ins-detphysics.app-phphysics.comp-phphysics.opticsSymbolic Computation

Catalog footprint

What is connected

24works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

NRPyElliptic: A Fast Hyperbolic Relaxation Elliptic Solver for Numerical Relativity, I: Conformally Flat, Binary Puncture Initial Data

We introduce NRPyElliptic, an elliptic solver for numerical relativity (NR) built within the NRPy+ framework. As its first application, NRPyElliptic sets up conformally flat, binary black hole (BBH) puncture initial data (ID) on a single numerical domain, similar to the widely used TwoPunctures code. Unlike TwoPunctures, NRPyElliptic employs a hyperbolic relaxation scheme, whereby arbitrary elliptic PDEs are trivially transformed into a hyperbolic system of PDEs. As consumers of NR ID generally already possess expertise in solving hyperbolic PDEs, they will generally find NRPyElliptic easier to tweak and extend than other NR elliptic solvers. When evolved forward in (pseudo)time, the hyperbolic system exponentially reaches a steady state that solves the elliptic PDEs. Notably NRPyElliptic accelerates the relaxation waves, which makes it many orders of magnitude faster than the usual constant-wavespeed approach. While it is still ${\sim}12$x slower than TwoPunctures at setting up full-3D BBH ID, NRPyElliptic requires only ${\approx}0.3\%$ of the runtime for a full BBH simulation in the Einstein Toolkit. Future work will focus on improving performance and generating other types of ID, such as binary neutron star.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

NRPyLaTeX: A LaTeX interface to computer algebra systems for general relativity

While each computer algebra system (CAS) contains its own unique syntax for inputting mathematical expressions, LaTeX is perhaps the most widespread language for typesetting mathematics. NRPyLaTeX (NL) enables direct LaTeX input of complex tensorial expressions (written in Einstein notation) relevant to general relativity and differential geometry into the SymPy CAS. As SymPy also supports output compatible with the Mathematica and Maple CASs, NL lowers the learning curve for inputting and manipulating tensorial expressions in three widely used CASs. LaTeX however is a typesetting language, and as such is not designed to resolve ambiguities in mathematical expressions. To address this, NL implements a convenient configuration interface that, e.g., defines variables with certain attributes. Configuration commands appear as LaTeX comments, so that entire NL workflows can fit seamlessly into the LaTeX source code of scientific papers without interfering with the rendered mathematical expressions. Further, NL adopts NRPy+'s rigid syntax for indexed symbols (e.g., tensors), which enables NL output to be directly converted into highly optimized C/C++-code kernels using NRPy+. Finally NL has robust and user-friendly error-handling, which catches common tensor indexing errors and reports unresolved ambiguities, further expediting the input and validation of LaTeX expressions into a CAS.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Active Optical Table Tilt Stabilization

We show that a simple modification to an optical table with pneumatic vibration isolation can be used to actively reduce the long term drift in the tilt of the table by nearly a factor of 1000. Without active stabilization, we measure a root-mean-square (RMS) tilt variation of \SI{270}{\upmu rad} over three days. The active stabilization can be used to limit the tilt to \SI{0.35}{\upmu rad} RMS over the same time period. This technique can be used to minimize drift in tilt-sensitive experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

High sensitivity accelerometry with a feedback-cooled magnetically levitated microsphere

We show that a magnetically levitated microsphere in high vacuum can be used as an accelerometer by comparing its response to that of a commercially available geophone. This system shows great promise for ultrahigh acceleration sensitivities without the need for large masses or cryogenics. With feedback cooling, the transient decay time is reduced and the center-of-mass motion is cooled to \SI{9}{K} or less. Remarkably, the levitated particle accelerometer has a sensitivity down to \SI{3.6e-8}{\mathit{g}/\sqrt{Hz}} and gives measurements similar to those of the commercial geophone at frequencies up to \SI{14}{Hz} despite a test mass that is four billion times smaller. With no free parameters in the calibration, the responses of the accelerometers match within \num{3}\% at \SI{5}{Hz}. The system reaches this sensitivity due to a relatively large particle mass of \SI{0.25}{\upmu g}, a low center of mass oscillation frequency of \SI{1.75}{Hz}, and a novel image analysis method that can measure the displacement with an uncertainty of \SI{1.6}{nm} in a single image.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical relativity in spherical coordinates: A new dynamical spacetime and general relativistic MHD evolution framework for the Einstein Toolkit

We present SphericalNR, a new framework for the publicly available Einstein Toolkit that numerically solves the Einstein field equations coupled to the equations of general relativistic magnetohydrodynamics (GRMHD) in a 3+1 split of spacetime in spherical coordinates without symmetry assumptions. The spacetime evolution is performed using reference-metric versions of either the Baumgarte-Shapiro-Shibata-Nakamura equations or the fully covariant and conformal Z4 system with constraint damping. We have developed a reference-metric version of the Valencia formulation of GRMHD with a vector potential method, guaranteeing the absence of magnetic monopoles during the evolution. In our framework, every dynamical field (both spacetime and matter) is evolved using its components in an orthonormal basis with respect to the spherical reference-metric. Furthermore, all geometric information about the spherical coordinate system is encoded in source terms appearing in the evolution equations. This allows for the straightforward extension of Cartesian high-resolution shock-capturing finite volume codes to use spherical coordinates with our framework. To this end, we have adapted GRHydro, a Cartesian finite volume GRMHD code already available in the Einstein Toolkit, to use spherical coordinates. We present the full evolution equations of the framework, as well as details of its implementation in the Einstein Toolkit. We validate SphericalNR by demonstrating it passes a variety of challenging code tests in static and dynamical spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimizing spinning time-domain gravitational waveforms for Advanced LIGO data analysis

The Spinning Effective One Body-Numerical Relativity (SEOBNR) series of gravitational wave approximants are among the best available for Advanced LIGO data analysis. Unfortunately, SEOBNR codes as they currently exist within LALSuite are generally too slow to be directly useful for standard Markov-Chain Monte Carlo-based parameter estimation (PE). Reduced-Order Models (ROMs) of SEOBNR have been developed for this purpose, but there is no known way to make ROMs of the full eight-dimensional intrinsic parameter space more efficient for PE than the SEOBNR codes directly. So as a proof of principle, we have sped up the original LALSuite SEOBNRv2 approximant code, which models waveforms from aligned-spin systems, by nearly 300x. Our optimized code shortens the timescale for conducting PE with this approximant to months, assuming a purely serial analysis, so that even modest parallelization combined with our optimized code will make running the full PE pipeline with SEOBNR codes directly a realistic possibility. A number of our SEOBNRv2 optimizations have already been applied to SEOBNRv3, a new approximant capable of modeling sources with all eight (precessing) intrinsic degrees of freedom. We anticipate that once all of our optimizations have been applied to SEOBNRv3, a similar speed-up may be achieved.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

IllinoisGRMHD: An Open-Source, User-Friendly GRMHD Code for Dynamical Spacetimes

In the extreme violence of merger and mass accretion, compact objects like black holes and neutron stars are thought to launch some of the most luminous outbursts of electromagnetic and gravitational wave energy in the Universe. Modeling these systems realistically is a central problem in theoretical astrophysics, but has proven extremely challenging, requiring the development of numerical relativity codes that solve Einstein's equations for the spacetime, coupled to the equations of general relativistic (ideal) magnetohydrodynamics (GRMHD) for the magnetized fluids. Over the past decade, the Illinois Numerical Relativity (ILNR) Group's dynamical spacetime GRMHD code has proven itself as a robust and reliable tool for theoretical modeling of such GRMHD phenomena. However, the code was written "by experts and for experts" of the code, with a steep learning curve that would severely hinder community adoption if it were open-sourced. Here we present IllinoisGRMHD, which is an open-source, highly-extensible rewrite of the original closed-source GRMHD code of the ILNR Group. Reducing the learning curve was the primary focus of this rewrite, with the goal of facilitating community involvement in the code's use and development, as well as the minimization of human effort in generating new science. IllinoisGRMHD also saves computer time, generating roundoff-precision identical output to the original code on adaptive-mesh grids, but nearly twice as fast at scales of hundreds to thousands of cores.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Accretion disks around binary black holes of unequal mass: GRMHD simulations of postdecoupling and merger

We report results from simulations in general relativity of magnetized disks accreting onto merging black hole binaries, starting from relaxed disk initial data. The simulations feature an effective, rapid radiative cooling scheme as a limiting case of future treatments with radiative transfer. Here we evolve the systems after binary-disk decoupling through inspiral and merger, and analyze the dependence on the binary mass ratio with $q\equiv m_{\rm bh}/M_{\rm BH}=1,1/2,$ and $1/4$. We find that the luminosity associated with local cooling is larger than the luminosity associated with matter kinetic outflows, while the electromagnetic (Poynting) luminosity associated with bulk transport of magnetic field energy is the smallest. The cooling luminosity around merger is only marginally smaller than that of a single, non-spinning black hole. Incipient jets are launched independently of the mass ratio, while the same initial disk accreting on a single non-spinning black hole does not lead to a jet, as expected. For all mass ratios we see a transient behavior in the collimated, magnetized outflows lasting $2-5 ( M/10^8M_\odot ) \rm days$ after merger: the outflows become increasingly magnetically dominated and accelerated to higher velocities, boosting the Poynting luminosity. These sudden changes can alter the electromagnetic emission across the jet and potentially help distinguish mergers of black holes in AGNs from single accreting black holes based on jet morphology alone.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Improved Moving Puncture Gauge Conditions for Compact Binary Evolutions

Robust gauge conditions are critically important to the stability and accuracy of numerical relativity (NR) simulations involving compact objects. Most of the NR community use the highly robust---though decade-old---moving-puncture (MP) gauge conditions for such simulations. It has been argued that in binary black hole (BBH) evolutions adopting this gauge, noise generated near adaptive-mesh-refinement (AMR) boundaries does not converge away cleanly with increasing resolution, severely limiting gravitational waveform accuracy at computationally feasible resolutions. We link this noise to a sharp (short-wavelength), initial outgoing gauge wave crossing into progressively lower resolution AMR grids, and present improvements to the standard MP gauge conditions that focus on stretching, smoothing, and more rapidly settling this outgoing wave. Our best gauge choice greatly reduces gravitational waveform noise during inspiral, yielding less fluctuation in convergence order and $\sim 40%$ lower waveform phase and amplitude errors at typical resolutions. Noise in other physical quantities of interest is also reduced, and constraint violations drop by more than an order of magnitude. We expect these improvements will carry over to simulations of all types of compact binary systems, as well as other $N$+1 formulations of gravity for which MP-like gauge conditions can be chosen.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Accretion disks around binary black holes of unequal mass: GRMHD simulations near decoupling

We report on simulations in general relativity of magnetized disks onto black hole binaries. We vary the binary mass ratio from 1:1 to 1:10 and evolve the systems when they orbit near the binary-disk decoupling radius. We compare (surface) density profiles, accretion rates (relative to a single, non-spinning black hole), variability, effective $α$-stress levels and luminosities as functions of the mass ratio. We treat the disks in two limiting regimes: rapid radiative cooling and no radiative cooling. The magnetic field lines clearly reveal jets emerging from both black hole horizons and merging into one common jet at large distances. The magnetic fields give rise to much stronger shock heating than the pure hydrodynamic flows, completely alter the disk structure, and boost accretion rates and luminosities. Accretion streams near the horizons are among the densest structures; in fact, the 1:10 no-cooling evolution results in a refilling of the cavity. The typical effective temperature in the bulk of the disk is $\sim 10^5 (M/10^8 M_\odot)^{-1/4} (L/L_{\rm edd})^{1/4} {\rm K}$ yielding characteristic thermal frequencies $\sim 10^{15} (M/10^8 M_\odot)^{-1/4} (L/L_{\rm edd})^{1/4}(1+z)^{-1}{\rm Hz} $. These systems are thus promising targets for many extragalactic optical surveys, such as LSST, WFIRST, and PanSTARRS.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration

The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios <= 4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

General relativistic simulations of binary black hole-neutron stars: Precursor electromagnetic signals

We perform the first general relativistic force-free simulations of neutron star (NS) magnetospheres in orbit about spinning and non-spinning black holes. We find promising precursor electromagnetic emission: typical Poynting luminosities at, e.g., an orbital separation of 6.6 times the NS radius are L ~ 6 x 10^{42} erg/s for a 1.4 solar-mass NS with a 10^{13}G polar magnetic field. The Poynting flux peaks within a broad beam of ~40 degrees in the azimuthal direction and within ~60 degrees from the orbital plane, establishing a possible lighthouse effect. Our calculations, though preliminary, preview more detailed simulations of these systems that we plan to perform in the future.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical Relativity Simulations of Magnetized Black Hole-Neutron Star Mergers

We present new numerical techniques we developed for launching the first parameter study of {\it magnetized} black hole--neutron star (BHNS) mergers, varying the magnetic fields seeded in the initial neutron star. We found that magnetic fields have a negligible impact on the gravitational waveforms and bulk dynamics of the system during merger, regardless of magnetic field strength or BH spin. In a recent simulation, we seeded the remnant disk from an unmagnetized BHNS merger simulation with large-scale, purely poloidal magnetic fields, which are otherwise absent in the full simulation. The outcome appears to be a viable sGRB central engine.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Binary Black-Hole Mergers in Magnetized Disks: Simulations in Full General Relativity

We present results from the first fully general relativistic, magnetohydrodynamic (GRMHD) simulations of an equal-mass black hole binary (BHBH) in a magnetized, circumbinary accretion disk. We simulate both the pre and post-decoupling phases of a BHBH-disk system and both "cooling" and "no-cooling" gas flows. Prior to decoupling, the competition between the binary tidal torques and the effective viscous torques due to MHD turbulence depletes the disk interior to the binary orbit. However, it also induces a two-stream accretion flow and mildly relativistic polar outflows from the BHs. Following decoupling, but before gas fills the low-density "hollow" surrounding the remnant, the accretion rate is reduced, while there is a prompt electromagnetic (EM) luminosity enhancement following merger due to shock heating and accretion onto the spinning BH remnant. This investigation, though preliminary, previews more detailed GRMHD simulations we plan to perform in anticipation of future, simultaneous detections of gravitational and EM radiation from a merging BHBH-disk system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

General relativistic simulations of black hole-neutron star mergers: Effects of magnetic fields

As a neutron star (NS) is tidally disrupted by a black hole (BH) companion at the end of a BH-NS binary inspiral, its magnetic fields will be stretched and amplified. If sufficiently strong, these magnetic fields may impact the gravitational waveforms, merger evolution and mass of the remnant disk. Formation of highly-collimated magnetic field lines in the disk+spinning BH remnant may launch relativistic jets, providing the engine for a short-hard GRB. We analyze this scenario through fully general relativistic, magnetohydrodynamic (GRMHD) BHNS simulations from inspiral through merger and disk formation. Different initial magnetic field configurations and strengths are chosen for the NS interior for both nonspinning and moderately spinning (a/M=0.75) BHs aligned with the orbital angular momentum. Only strong interior (Bmax~10^17 G) initial magnetic fields in the NS significantly influence merger dynamics, enhancing the remnant disk mass by 100% and 40% in the nonspinning and spinning BH cases, respectively. However, detecting the imprint of even a strong magnetic field may be challenging for Advanced LIGO. Though there is no evidence of mass outflows or magnetic field collimation during the preliminary simulations we have performed, higher resolution, coupled with longer disk evolutions and different initial magnetic field configurations, may be required to definitively assess the possibility of BHNS binaries as short-hard GRB progenitors.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

General relativistic simulations of black hole-neutron star mergers: Effects of tilted magnetic fields

Black hole--neutron star (BHNS) binary mergers can form disks in which magnetorotational instability (MRI)-induced turbulence may drive accretion onto the remnant BH, supporting relativistic jets and providing the engine for a short-hard gamma-ray burst (SGRB). Our earlier study of magnetized BHNSs showed that NS tidal disruption winds the magnetic field into a toroidal configuration, with poloidal fields so weak that capturing MRI with full-disk simulations would require $\sim 10^8$ CPU-hours. In that study we imposed equatorial symmetry, suppressing poloidal magnetic fields that might be generated from plasma crossing the orbital plane. Here we show that initial conditions that break this symmetry (i.e., {\it tilted} poloidal magnetic fields in the NS) generate much stronger poloidal fields in the disk, indicating that asymmetric initial conditions may be necessary for establishing BHNS mergers as SGRB progenitors via fully general relativistic MHD simulations. We demonstrate that BHNS mergers may form an SGRB engine under the right conditions by seeding the remnant disk from an unmagnetized BHNS simulation with purely poloidal fields dynamically unimportant initially, but strong enough to resolve MRI. Magnetic turbulence occurs in the disk, driving accretion and supporting Poynting-dominated jet outflows sufficient to power an SGRB.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Importance of cooling in triggering the collapse of hypermassive neutron stars

The inspiral and merger of a binary neutron star (NSNS) can lead to the formation of a hypermassive neutron star (HMNS). As the HMNS loses thermal pressure due to neutrino cooling and/or centrifugal support due to gravitational wave (GW) emission, and/or magnetic breaking of differential rotation it will collapse to a black hole. To assess the importance of shock-induced thermal pressure and cooling, we adopt an idealized equation of state and perform NSNS simulations in full GR through late inspiral, merger, and HMNS formation, accounting for cooling. We show that thermal pressure contributes significantly to the support of the HMNS against collapse and that thermal cooling accelerates its "delayed" collapse. Our simulations demonstrate explicitly that cooling can induce the catastrophic collapse of a hot hypermassive neutron star formed following the merger of binary neutron stars. Thus, cooling physics is important to include in NSNS merger calculations to accurately determine the lifetime of the HMNS remnant and to extract information about the NS equation of state, cooling mechanisms, bar instabilities and B-fields from the GWs emitted during the transient phase prior to BH formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Relativistic MHD in dynamical spacetimes: Improved EM gauge condition for AMR grids

We recently developed a new general relativistic magnetohydrodynamic code with adaptive mesh refinement that evolves the electromagnetic (EM) vector potential (A) instead of the magnetic fields directly. Evolving A enables one to use any interpolation scheme on refinement level boundaries and still guarantee that the magnetic field remains divergenceless. As in classical EM, a gauge choice must be made when evolving A, and we chose a straightforward "algebraic" gauge condition to simplify the A evolution equation. However, magnetized black hole-neutron star (BHNS) simulations in this gauge exhibit unphysical behavior, including the spurious appearance of strong magnetic fields on refinement level boundaries. This spurious behavior is exacerbated when matter crosses refinement boundaries during tidal disruption of the NS. Applying Kreiss-Oliger dissipation to the evolution of the magnetic vector potential A slightly weakens this spurious magnetic effect, but with undesired consequences. We demonstrate via an eigenvalue analysis and a numerical study that zero-speed modes in the algebraic gauge, coupled with the frequency filtering that occurs on refinement level boundaries, are responsible for the creation of spurious magnetic fields. We show that the EM Lorenz gauge exhibits no zero-speed modes, and as a consequence, spurious magnetic effects are quickly propagated away, allowing for long-term, stable magnetized BHNS evolutions. Our study demonstrates how the EM gauge degree of freedom can be chosen to one's advantage, and that for magnetized BHNS simulations the Lorenz gauge constitutes a major improvement over the algebraic gauge.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search and parameter-estimation algorithms using numerically generated waveforms, and to foster closer collaboration between the numerical relativity and data analysis communities. The first NINJA project used only a small number of injections of short numerical-relativity waveforms, which limited its ability to draw quantitative conclusions. The goal of the NINJA-2 project is to overcome these limitations with long post-Newtonian - numerical relativity hybrid waveforms, large numbers of injections, and the use of real detector data. We report on the submission requirements for the NINJA-2 project and the construction of the waveform catalog. Eight numerical relativity groups have contributed 63 hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. We summarize the techniques used by each group in constructing their submissions. We also report on the procedures used to validate these submissions, including examination in the time and frequency domains and comparisons of waveforms from different groups against each other. These procedures have so far considered only the $(\ell,m)=(2,2)$ mode. Based on these studies we judge that the hybrid waveforms are suitable for NINJA-2 studies. We note some of the plans for these investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evolution of near-extremal-spin black holes using the moving puncture technique

We propose a new radial coordinate to write the Kerr metric in puncture form. Unlike the quasi-radial coordinate introduced previously, the horizon radius remains finite in our radial coordinate in the extreme Kerr limit a/M -> 1. This significantly improves the accuracy of the evolution of black holes with spins close to the extreme Kerr limit. We are able to evolve accurately both stationary and boosted black holes with spins as high as a/M=0.99 using initial data constructed in these new puncture coordinates. Initial data of compact binaries with rapidly spinning black holes can be constructed using our proposed new puncture metric for the background conformal metric. Our simulations for single black holes suggest that such initial data can be evolved successfully by the moving puncture technique.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Relativistic magnetohydrodynamics in dynamical spacetimes: A new AMR implementation

We have written and tested a new general relativistic magnetohydrodynamics (GRMHD) code, capable of evolving MHD fluids in dynamical spacetimes with adaptive-mesh refinement (AMR). Our code solves the Einstein-Maxwell-MHD system of coupled equations in full 3+1 dimensions, evolving the metric via the Baumgarte-Shapiro Shibata-Nakamura (BSSN) formalism and the MHD and magnetic induction equations via a conservative, high-resolution shock-capturing scheme. The induction equations are recast as an evolution equation for the magnetic vector potential, which exists on a grid that is staggered with respect to the hydrodynamic and metric variables. The divergenceless constraint div(B)=0 is enforced by the curl of the vector potential. Our MHD scheme is fully compatible with AMR, so that fluids at AMR refinement boundaries maintain div(B)=0. In simulations with uniform grid spacing, our MHD scheme is numerically equivalent to a commonly used, staggered-mesh constrained-transport scheme. We present code validation test results, both in Minkowski and curved spacetimes. They include magnetized shocks, nonlinear Alfvén waves, cylindrical explosions, cylindrical rotating disks, magnetized Bondi tests, and the collapse of a magnetized rotating star. Some of the more stringent tests involve black holes. We find good agreement between analytic and numerical solutions in these tests, and achieve convergence at the expected order.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Status of NINJA: the Numerical INJection Analysis project

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.

0 citations0 reviewsNo code linkedOpen access