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

Stuart L. Shapiro

Stuart L. Shapiro 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-phastro-ph.SRastro-ph.COastro-ph.GAhep-phhep-th

Catalog footprint

What is connected

61works
8topics
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

61 published item(s)

preprint2022arXiv

Jet Launching from Binary Neutron Star Mergers: Incorporating Neutrino Transport and Magnetic Fields

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, \barν_e, ν_{\rm 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 $\sim 10^{53} \rm \, 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 $\sim 10^{52-53} \rm \,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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Accretion onto black holes inside neutron stars with piecewise-polytropic equations of state: analytic and numerical treatments

We consider spherically symmetric accretion onto a small, possibly primordial, black hole residing at the center of a neutron star governed by a cold nuclear equation of state (EOS). We generalize the relativistic Bondi solution for such EOSs, approximated by piecewise polytropes, and thereby obtain analytical expressions for the steady-state matter profiles and accretion rates. We compare these rates with those found by time-dependent, general relativistic hydrodynamical simulations upon relaxation and find excellent agreement. We consider several different candidate EOSs, neutron star masses and central densities and find that the accretion rates vary only little, resulting in an accretion rate that depends primarily on the black hole mass, and only weakly on the properties of the neutron star.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Black hole-neutron star coalescence: effects of the neutron star spin on jet launching and dynamical ejecta mass

Black hole-neutron star (BHNS) mergers are thought to be sources of gravitational waves (GWs) with coincident electromagnetic (EM) counterparts. To further probe whether these systems are viable progenitors of short gamma-ray bursts (sGRBs) and kilonovae, and how one may use (the lack of) EM counterparts associated with LIGO/Virgo candidate BHNS GW events to sharpen parameter estimation, we study the impact of neutron star spin in BHNS mergers. Using dynamical spacetime magnetohydrodynamic simulations of BHNSs initially on a quasicircular orbit, we survey configurations that differ in the BH spin ($a_{\rm BH}/M_{\rm BH}=0$ and $0.75$), the NS spin ($a_{\rm NS}/M_{\rm NS}=-0.17,\,0,\,0.23$ and $0.33$), and the binary mass ratio ($q\equiv M_{\rm BH}:M_{\rm NS}=3:1$ and $5:1$). The general trend we find is that increasing the NS prograde spin increases both the rest mass of the accretion disk onto the remnant black hole, and the rest mass of dynamically ejected matter. By a time~$Δt\sim 3500-5500M\sim 88-138(M_{\rm NS}/1.4M_\odot)\,\rm ms$ after the peak gravitational wave amplitude, a magnetically--driven jet is launched only for $q=3:1$ regardless of the initial NS spin. The lifetime of the jets [$Δt\sim 0.5-0.8(M_{\rm NS}/1.4 M_\odot)\,\rm s$] and their outgoing Poynting luminosity [$L_{\rm Poyn}\sim 10^{51.5\pm 0.5}\,\rm erg/s$] are consistent with typical sGRBs luminosities and expectations from the Blandford-Znajek mechanism. By the time we terminate our simulations, we do not observe either an outflow or a large-scale magnetic field collimation for the other systems we considered. The mass range of dynamically ejected matter is $10^{-4.5}-10^{-2}~(M_{\rm NS}/1.4M_\odot)M_\odot$, which can power kilonovae with peak bolometric luminosities $L_{\rm knova}\sim 10^{40}-10^{41.4}$ erg/s with rise times $\lesssim 6.5\,\rm h$ and potentially detectable by the LSST.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Magnetohydrodynamic simulations of self-consistent rotating neutron stars with mixed poloidal and toroidal magnetic fields

We perform the first magnetohydrodynamic simulations in full general relativity of self-consistent rotating neutron stars (NSs) with ultrastrong mixed poloidal and toroidal magnetic fields. The initial uniformly rotating NS models are computed assuming perfect conductivity, stationarity, and axisymmetry. Although the specific geometry of the mixed field configuration can delay or accelerate the development of various instabilities known from analytic perturbative studies, all our models finally succumb to them. Differential rotation is developed spontaneously in the cores of our magnetars which, after sufficient time, is converted back to uniform rotation. The rapidly rotating magnetars show a significant amount of ejecta, which can be responsible for transient kilonova signatures. However no highly collimated, helical magnetic fields or incipient jets, which are necessary for gamma-ray bursts, arise at the poles of these magnetars by the time our simulations are terminated.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Faint Dark Matter Annihilation Signals and the Milky Way's Supermassive Black Hole

A wide range of mechanisms predict present-day s-wave dark matter (DM) annihilation cross-sections that are orders of magnitude below current experimental sensitivity. We explore the capability of DM density spikes around the Milky Way's supermassive black hole to probe such faint signals of DM annihilations, considering a range of possible spike and halo distributions. As an exemplar of a theory with a suppressed s-wave annihilation cross-section, we consider a hidden sector axion portal model of DM. In this model, the leading contribution to the annihilation cross-section in the early universe is p-wave, while s-wave annihilations occur at higher order in the coupling constant. We provide a unified treatment of DM freezeout in this model including both s- and p-wave annihilations and analytically determine the photon spectrum for the dominant DM annihilation process in the universe today. We find that Fermi and H.E.S.S. observations of the Galactic Center offer excellent sensitivity to this model over a wide range of parameter space, with prospects depending sensitively on the properties of the DM spike as well as the central halo.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gravitational Waves from Disks Around Spinning Black Holes: Simulations in Full General Relativity

We present fully general-relativistic numerical evolutions of self-gravitating tori around spinning black holes with dimensionless spin $a/M = 0.7$ parallel or anti-parallel to the disk angular momentum. The initial disks are unstable to the hydrodynamic Papaloizou-Pringle Instability which causes them to grow persistent orbiting matter clumps. The effect of black hole spin on the growth and saturation of the instability is assessed. We find that the instability behaves similarly to prior simulations with non-spinning black holes, with a shift in frequency due to spin-induced changes in disk orbital period. Copious gravitational waves are generated by these systems, and we analyze their detectability by current and future gravitational wave observatories for large range of masses. We find that systems of $10 M_\odot$ - relevant for black hole-neutron star mergers - are detectable by Cosmic Explorer out to $\sim300$ Mpc, while DECIGO (LISA) will be able to detect systems of $1000 M_\odot$ ($10^5M_\odot$) - relevant for disks forming in collapsing supermassive stars - out to cosmological redshift of $z\sim5$ ($z\sim 1$). Computing the accretion rate of these systems we find that these systems may also be promising sources of coincident electromagnetic signals.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Great Impostors: Extremely Compact, Merging Binary Neutron Stars in the Mass Gap Posing as Binary Black Holes

Can one distinguish a binary black hole undergoing a merger from a binary neutron star if the individual compact companions have masses that fall inside the so-called mass gap of $3-5\ M_\odot$? For neutron stars, achieving such masses typically requires extreme compactness and in this work we present initial data and evolutions of binary neutron stars initially in quasiequilibrium circular orbits having a compactness $C=0.336$. These are the most compact, nonvacuum, quasiequilibrium binary objects that have been constructed and evolved to date, including boson stars. The compactness achieved is only slightly smaller than the maximum possible imposed by causality, $C_{\rm max}=0.355$, which requires the sound speed to be less than the speed of light. By comparing the emitted gravitational waveforms from the late inspiral to merger and postmerger phases between such a binary neutron star vs a binary black hole of the same total mass we identify concrete measurements that serve to distinguish them. With that level of compactness, the binary neutron stars exhibit no tidal disruption up until merger, whereupon a prompt collapse is initiated even before a common core forms. Within the accuracy of our simulations the black hole remnants from both binaries exhibit ringdown radiation that is not distinguishable from a perturbed Kerr spacetime. However, their inspiral leads to phase differences of the order of $\sim 5$ rad over an $\sim 81$ km separation (1.7 orbits) while typical neutron stars exhibit phase differences of $\geq 20$ rad. Although a difference of $\sim 5$ rad can be measured by current gravitational wave laser interferometers (e.g., aLIGO/Virgo), uncertainties in the individual masses and spins will likely prevent distinguishing such compact, massive neutron stars from black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Locating ergostar models in parameter space

Recently, we have shown that dynamically stable ergostar solutions (equilibrium neutron stars that contain an ergoregion) with a compressible and causal equation of state exist [A. Tsokaros, M. Ruiz, L. Sun, S. L. Shapiro, and K. Uryū, Phys. Rev. Lett. 123, 231103 (2019)]. These stars are hypermassive, differentially rotating, and highly compact. In this work, we make a systematic study of equilibrium models in order to locate the position of ergostars in parameter space. We adopt four equations of state that differ in the matching density of a maximally stiff core. By constructing a large number of models both with uniform and differential rotation of different degrees, we identify the parameters for which ergostars appear. We find that the most favorable conditions for the appearance of dynamically stable ergostars are a significant finite density close to the surface of the star (i.e., similar to self-bound quark stars) and a small degree of differential rotation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetohydrodynamic Simulations of Binary Neutron Star Mergers in General Relativity: Effects of Magnetic Field Orientation on Jet Launching

Binary neutron star (NSNS) mergers can be sources of gravitational waves coincident with electromagnetic counterpart emission. To solidify their role as multimessenger sources, we present fully 3D, general relativistic, magnetohydrodynamic simulations of spinning NSNSs initially on quasicircular orbits that merge and undergo delayed collapse to a black hole (BH). The NSNSs consist of two identical stars modeled as $Γ=2$ polytropes with spin $χ_{NS}= 0.36$ aligned along the direction of the total orbital angular momentum $L$. Each star is initially threaded by a dynamical unimportant interior dipole B-field. The field is extended into the exterior where a nearly force-free magnetosphere resembles that of a pulsar. The magnetic dipole moment $μ$ is either aligned or perpendicular to $L$ and has the same initial magnitude for each orientation. For comparison, we also impose symmetry across the orbital plane in one case where $μ$ in both stars is aligned along $L$. We find that the lifetime of the transient hypermassive neutron star remnant, the jet launching time, and the ejecta are very sensitive to the B-field orientation. By contrast, the physical properties of the BH + disk remnant, such as the mass and spin of the BH, the accretion rate, and the electromagnetic luminosity, are roughly independent of the initial B-field orientation. In addition, we find imposing symmetry across the orbital plane does not play a significant role in the final outcome of the mergers. Our results show that an incipient jet emerges only when the seed B-field has a sufficiently large-scale poloidal component aligned to $L$. The lifetime [$Δt\gtrsim 140(M_{NS}/1.625M_\odot)\rm ms$] and Poynting luminosities [$L_{EM}\simeq 10^{52}$erg/s] of the jet, when it forms, are consistent with typical short gamma-ray bursts, as well as with the Blandford--Znajek mechanism for launching jets.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Binary neutron star mergers: a jet engine for short gamma-ray bursts

We perform magnetohydrodynamic simulations in full general relativity (GRMHD) of quasi-circular, equal-mass, binary neutron stars that undergo merger. The initial stars are irrotational, $n=1$ polytropes and are magnetized. We explore two types of magnetic-field geometries: one where each star is endowed with a dipole magnetic field extending from the interior into the exterior, as in a pulsar, and the other where the dipole field is initially confined to the interior. In both cases the adopted magnetic fields are initially dynamically unimportant. The merger outcome is a hypermassive neutron star that undergoes delayed collapse to a black hole (spin parameter $a/M_{\rm BH} \sim 0.74$) immersed in a magnetized accretion disk. About $4000M \sim 60(M_{\rm NS}/1.625M_\odot)$ ms following merger, the region above the black hole poles becomes strongly magnetized, and a collimated, mildly relativistic outflow --- an incipient jet --- is launched. The lifetime of the accretion disk, which likely equals the lifetime of the jet, is $Δt \sim 0.1 (M_{\rm NS}/1.625M_\odot)$ s. In contrast to black hole--neutron star mergers, we find that incipient jets are launched even when the initial magnetic field is confined to the interior of the stars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Relativistic Simulations of Eccentric Binary Neutron Star Mergers: One-arm Spiral Instability and Effects of Neutron Star Spin

We perform general-relativistic hydrodynamical simulations of dynamical capture binary neutron star mergers, emphasizing the role played by the neutron star spin. Dynamical capture mergers may take place in globular clusters, as well as other dense stellar systems, where most neutron stars have large spins. We find significant variability in the merger outcome as a function of initial neutron star spin. For cases where the spin is aligned with the orbital angular momentum, the additional centrifugal support in the remnant hypermassive neutron star can prevent the prompt collapse to a black hole, while for antialigned cases the decreased total angular momentum can facilitate the collapse to a black hole. We show that even moderate spins can significantly increase the amount of ejected material, including the amount unbound with velocities greater than half the speed of light, leading to brighter electromagnetic signatures associated with kilonovae and interaction of the ejecta with the interstellar medium. Furthermore, we find that the initial neutron star spin can strongly affect the already rich phenomenology in the postmerger gravitational wave signatures that arise from the oscillation modes of the hypermassive neutron star. In several of our simulations, the resulting hypermassive neutron star develops the one-arm ($m=1$) spiral instability, the most pronounced cases being those with small but non-negligible neutron star spins. For long-lived hypermassive neutron stars, the presence of this instability leads to improved prospects for detecting these events through gravitational waves, and thus may give information about the neutron star equation of state.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Weak annihilation cusp inside the dark matter spike about a black hole

We reinvestigate the effect of annihilations on the distribution of collisionless dark matter (DM) in a spherical density spike around a massive black hole. We first construct a very simple, pedagogic, analytic model for an isotropic phase space distribution function that accounts for annihilation and reproduces the "weak cusp" found by Vasiliev for DM deep within the spike and away from its boundaries. The DM density in the cusp varies as $r^{-1/2}$ for $s$-wave annihilation, where $r$ is the distance from the central black hole, and is not a flat "plateau" profile. We then extend this model by incorporating a loss cone that accounts for the capture of DM particles by the hole. The loss cone is implemented by a boundary condition that removes capture orbits, resulting in an anisotropic distribution function. Finally, we evolve an initial spike distribution function by integrating the Boltzmann equation to show how the weak cusp grows and its density decreases with time. We treat two cases, one for $s$-wave and the other for $p$-wave DM annihilation, adopting parameters characteristic of the Milky Way nuclear core and typical WIMP models for DM. The cusp density profile for $p$-wave annihilation is weaker, varying like $\sim r^{-0.34}$, but is still not a flat plateau.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A black hole window into p-wave dark matter annihilation

We present a new method to measure or constrain p-wave-suppressed cross sections for dark matter (DM) annihilations inside the steep density spikes induced by supermassive black holes. We demonstrate that the high DM densities, together with the increased velocity dispersion, within such spikes combine to make thermal p-wave annihilation cross-sections potentially visible in gamma-ray observations of the Galactic center (GC). The resulting DM signal is a bright central point source with emission originating from DM annihilations in the absence of a detectable spatially-extended signal from the halo. We define two simple reference theories of DM with a thermal p-wave annihilation cross-section and establish new limits on the combined particle and astrophysical parameter space of these models, demonstrating that Fermi is currently sensitive to thermal p-wave DM over a wide range of possible scenarios for the DM distribution in the GC.

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
preprint2015arXiv

One-arm Spiral Instability in Hypermassive Neutron Stars Formed by Dynamical-Capture Binary Neutron Star Mergers

Using general-relativistic hydrodynamical simulations, we show that merging binary neutron stars can form hypermassive neutrons stars that undergo the one-arm spiral instability. We study the particular case of a dynamical capture merger where the stars have a small spin, as may arise in globular clusters, and focus on an equal-mass scenario where the spins are aligned with the orbital angular momentum. We find that this instability develops when post-merger fluid vortices lead to the generation of a toroidal remnant - a configuration whose maximum density occurs in a ring around the center-of-mass - with high vorticity along its rotation axis. The instability quickly saturates on a timescale of $\sim 10$ ms, with the $m=1$ azimuthal density multipole mode dominating over higher modes. The instability also leaves a characteristic imprint on the post-merger gravitational wave signal that could be detectable if the instability persists in long-lived remnants.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Relativistic simulations of black hole-neutron star coalescence: the jet emerges

We perform magnetohydrodynamic simulations in full general relativity (GRMHD) of a binary black hole-neutron star on a quasicircular orbit that undergoes merger. The binary mass ratio is 3:1, the black hole initial spin parameter $a/m=0.75$ ($m$ is the black hole Christodoulou mass) aligned with the orbital angular momentum, and the neutron star is an irrotational $Γ=2$ polytrope. About two orbits prior to merger (at time $t=t_B$), we seed the neutron star with a dynamically weak interior dipole magnetic field that extends into the stellar exterior. At $t=t_B$ the exterior has a low-density atmosphere with constant plasma parameter $β\equiv P_{\rm gas}/P_{\rm mag}$. Varying $β$ at $t_B$ in the exterior from $0.1$ to $0.01$, we find that at a time $\sim 4000{\rm M} \sim 100(M_{\rm NS}/1.4M_\odot)$ms [M is the total (ADM) mass] following the onset of accretion of tidally disrupted debris, magnetic winding above the remnant black hole poles builds up the magnetic field sufficiently to launch a mildly relativistic, collimated outflow - an incipient jet. The duration of the accretion and the lifetime of the jet is $Δt\sim 0.5(M_{\rm NS}/1.4M_\odot)$s. Our simulations furnish the first explicit examples in GRMHD which show that a jet can emerge following a black hole - neutron star merger.

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

Galactic Center Gamma-Ray Excess from Dark Matter Annihilation: Is There A Black Hole Spike?

If the supermassive black hole Sgr A* at the center of the Milky Way grew adiabatically from an initial seed embedded in an NFW dark matter (DM) halo, then the DM profile near the hole has steepened into a spike. We calculate the dramatic enhancement to the gamma ray flux from the Galactic center (GC) from such a spike if the 1-3 GeV excess observed in Fermi data is due to DM annihilations. We find that for the parameter values favored in recent fits, the point source-like flux from the spike is 35 times greater than the flux from the inner $1^\circ$ of the halo, far exceeding all Fermi point source detections near the GC. We consider the dependence of the spike signal on astrophysical and particle parameters and conclude that if the GC excess is due to DM, then a canonical adiabatic spike is disfavored by the data. We discuss alternative Galactic histories that predict different spike signals, including: (i) the nonadiabatic growth of the black hole, possibly associated with halo and/or black hole mergers, (ii) gravitational interaction of DM with baryons in the dense core, such as heating by stars, or (iii) DM self-interactions. We emphasize that the spike signal is sensitive to a different combination of particle parameters than the halo signal, and that the inclusion of a spike component to any DM signal in future analyses would provide novel information about both the history of the GC and the particle physics of DM annihilations.

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
preprint2014arXiv

Pulsar spin-down luminosity: Simulations in general relativity

Adopting our new method for matching general relativistic, ideal magnetohydrodynamics to its force-free limit, we perform the first systematic simulations of force-free pulsar magnetospheres in general relativity. We endow the neutron star with a general relativistic dipole magnetic field, model the interior with ideal magnetohydrodynamics, and adopt force-free electrodynamics in the exterior. Comparing the spin-down luminosity to its corresponding Minkowski value, we find that general relativistic effects give rise to a modest enhancement: the maximum enhancement for $n=1$ polytropes is $\sim 23\%$. Evolving a rapidly rotating $n=0.5$ polytrope we find an even greater enhancement of $\sim 35\%$. Using our simulation data, we derive fitting formulas for the pulsar spin-down luminosity as a function of the neutron star compaction, angular speed, and dipole magnetic moment. We expect stiffer equations of state and more rapidly spinning neutron stars to lead to even larger enhancements in the spin-down luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-interacting dark matter cusps around massive black holes

We adopt the conduction fluid approximation to model the steady-state distribution of matter around a massive black hole at the center of a weakly collisional cluster of particles. By "`weakly collisional" we mean a cluster in which the mean free time between particle collisions is much longer than the characteristic particle crossing (dynamical) time scale, but shorter than the cluster lifetime. When applied to a star cluster, we reproduce the familiar Bahcall-Wolf power-law cusp solution for the stars bound to the black hole. Here the star density scales with radius as $r^{-7/4}$ and the velocity dispersion as $r^{-1/2}$ throughout most of the gravitational well of the black hole. When applied to a relaxed, self-interacting dark matter (SIDM) halo with a velocity-dependent cross section $σ\sim v^{-a}$, the gas again forms a power-law cusp, but now the SIDM density scales as $r^{-β}$, where $β= (a+3)/4$, while its velocity dispersion again varies as $r^{-1/2}$. Results are obtained first in Newtonian theory and then in full general relativity. Although the conduction fluid model is a simplification, it provides a reasonable first approximation to the matter profiles and is much easier to implement than a full Fokker-Planck treatment or an $N$-body simulation of the Boltzmann equation with collisional perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A new scheme for matching general relativistic ideal magnetohydrodynamics to its force-free limit

We present a new computational method for smoothly matching general relativistic ideal magnetohydrodynamics (MHD) to its force-free limit. The method is based on a flux-conservative formalism for MHD and its force-free limit, and a vector potential formulation for the induction equation to maintain the zero divergence constraint for the magnetic field. The force-free formulation we adopt evolves the magnetic field and the Poynting vector, instead of the magnetic and electric fields. We show that our force-free code passes a robust suite of tests, performed both in 1D flat spacetime and in 3D black hole spacetimes. We also demonstrate that our matching technique successfully reproduces the aligned rotator force-free solution. Our new techniques are suitable for studying electromagnetic effects and predicting electromagnetic signals arising in many different curved spacetime scenarios. For example, we can treat spinning neutron stars, either in isolation or in compact binaries, that have MHD interiors and force-free magnetospheres.

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

Accretion Disks Around Binary Black Holes: A Simple GR-Hybrid Evolution Model

We consider a geometrically thin, Keplerian disk in the orbital plane of a binary black hole (BHBH) consisting of a spinning primary and low-mass secondary (mass ratio q < 1). To account for the principle effects of general relativity (GR), we propose a modification of the standard Newtonian evolution equation for the (orbit-averaged) time-varying disk surface density. In our modified equation the viscous torque in the disk is treated in full GR, while the tidal torque is handled in the Newtonian limit. Our GR-hybrid treatment is reasonable because the tidal torque is concentrated near the orbital radius of the secondary and is most important prior to binary-disk decoupling, when the orbital separation is large and resides in the weak-field regime. The tidal torque on the disk diminishes during late merger and vanishes altogether following merger. By contrast, the viscous torque drives the flow into the strong-field region and onto the primary during all epochs. Following binary coalescence, the viscous torque alone governs the time-dependent accretion onto the remnant, as well as the temporal behavior, strength and spectrum of the aftermath electromagnetic radiation from the disk. We solve our GR-hybrid equation for a representative BHBH-disk system, identify several observable EM signatures of the merger, and compare results obtained for the gas and EM radiation with those found with the Newtonian prescription.

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

Accretion Disks Around Binary Black Holes: A Quasistationary Model

Tidal torques acting on a gaseous accretion disk around a binary black hole can create a gap in the disk near the orbital radius. At late times, when the binary inspiral timescale due to gravitational wave emission becomes shorter than the viscous timescale in the disk, the binary decouples from the disk and eventually merges. Prior to decoupling the balance between tidal and viscous torques drives the disk to a quasistationary equilibrium state, perturbed slightly by small amplitude, spiral density waves emanating from the edges of the gap. We consider a black hole binary with a companion of smaller mass and construct a simple Newtonian model for a geometrically thin, Keplerian disk in the orbital plane of the binary. We solve the disk evolution equations in steady state to determine the quasistationary, (orbit-averaged) surface density profile prior to decoupling. We use our solution, which is analytic up to simple quadratures, to compute the electromagnetic flux and approximate radiation spectrum during this epoch. A single nondimensional parameter Td/Tvis, equal to the ratio of the tidal to viscous torque at the orbital radius, determines the disk structure, including the surface density profile, the extent of the gap, the existence of an inner disk, and the accretion rate. The solution reduces to the Shakura-Sunyaev profile for a stationary accretion disk around a single black hole in the limit of small Td/Tvis. Our solution may be useful for choosing physical parameters and setting up quasistationary disk initial data for detailed numerical simulations that begin prior to decoupling and track the subsequent evolution of a black hole binary-disk system.

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

Gravity darkening and brightening in binaries

We apply a von Zeipel gravity darkening model to corotating binaries to obtain a simple, analytical expression for the emergent radiative flux from a tidally distorted primary orbiting a point-mass secondary. We adopt a simple Roche model to determine the envelope structure of the primary, assumed massive and centrally condensed, and use the results to calculate the flux. As for single rotating stars, gravity darkening reduces the flux along the stellar equator of the primary, but, unlike for rotating stars, we find that gravity brightening enhances the flux in a region around the stellar poles. We identify a critical limiting separation beyond which hydrostatic equilibrium no longer is possible, whereby the flux vanishes at the point on the stellar equator of the primary facing the companion. For equal-mass binaries, the total luminosity is reduced by about 13 % when this limiting separation is reached.

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 merger of binary white dwarf-neutron stars: Simulations in full general relativity

We present fully general relativistic (GR) simulations of binary white dwarf-neutron star (WDNS) inspiral and merger. The initial binary is in a circular orbit at the Roche critical separation. The goal is to determine the ultimate fate of such systems. We focus on binaries whose total mass exceeds the maximum mass (Mmax) a cold, degenerate EOS can support against gravitational collapse. The time and length scales span many orders of magnitude, making fully general relativistic hydrodynamic (GRHD) simulations computationally prohibitive. For this reason, we model the WD as a "pseudo-white dwarf" (pWD) as in our binary WDNS head-on collisions study [PRD83:064002,2011]. Our GRHD simulations of a pWDNS system with a 0.98-solar-mass WD and a 1.4-solar-mass NS show that the merger remnant is a spinning Thorne-Zytkow-like Object (TZlO) surrounded by a massive disk. The final total rest mass exceeds Mmax, but the remnant does not collapse promptly. To assess whether the object will ultimately collapse after cooling, we introduce radiative thermal cooling. We first apply our cooling algorithm to TZlOs formed in WDNS head-on collisions, and show that these objects collapse and form black holes on the cooling time scale, as expected. However, when we cool the spinning TZlO formed in the merger of a circular-orbit WDNS binary, the remnant does not collapse, demonstrating that differential rotational support is sufficient to prevent collapse. Given that the final total mass exceeds Mmax, magnetic fields and/or viscosity may redistribute angular momentum and ultimately lead to delayed collapse to a BH. We infer that the merger of realistic massive WDNS binaries likely will lead to the formation of spinning TZlOs that undergo delayed collapse.

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
preprint2011arXiv

Binary black hole mergers in gaseous disks: Simulations in general relativity

Simultaneous gravitational and electromagnetic wave observations of merging black hole binaries (BHBHs) can provide unique opportunities to study gravitation physics, accretion and cosmology. Here we perform fully general relativistic, hydrodynamic simulations of equal-mass, nonspinning BHBHs coalescing in a circumbinary disk. We evolve the metric using the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation of Einstein's field equations with standard moving puncture gauge conditions. We handle the hydrodynamics via a high-resolution shock-capturing (HRSC) scheme. We track the inspiral starting from a binary separation of 10M, where M is the total binary mass. We take the disks to have an inner radius at R_in~15M to account for the hollow created by the binary torques. Our disks extend to R=65M and have an initial scale height of H/R=0.03-0.11. The gas is governed by a Gamma-law EOS, with Gamma equal to 5/3, 4/3, and 1.1. Disks are allowed to relax in the "early inspiral" epoch to provide quasistationary realistic initial data. We then evolve the metric and matter during the "late inspiral and merger" epoch. The later simulations are designed to track BHBH inspiral following disk-binary decoupling, through merger and ringdown, terminating before viscosity has time to fill the hollow about the remnant. We compute the gas flow and accretion rate and estimate the electromagnetic luminosity due to bremsstrahlung and synchrotron emission as a perturbation for optically thin disks. The synchrotron component of the luminosity peaks in the infrared band and should be detectable by WFIRST and possibly the LSST for a 10^8 M_sun binary embedded in a disk with a density n~10^12/cm^3 at z=1, beginning with a maximum value of $L~10^46 n_12^2 M_8^3 erg/s at decoupling, and decreasing steadily over a timescale of ~100 M_8 hours to a value of L~10^45 n_12^2 M_8^3 erg/s at merger.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Constraint propagation equations of the 3+1 decomposition of f(R) gravity

Theories of gravity other than general relativity (GR) can explain the observed cosmic acceleration without a cosmological constant. One such class of theories of gravity is f(R). Metric f(R) theories have been proven to be equivalent to Brans-Dicke (BD) scalar-tensor gravity without a kinetic term. Using this equivalence and a 3+1 decomposition of the theory it has been shown that metric f(R) gravity admits a well-posed initial value problem. However, it has not been proven that the 3+1 evolution equations of metric f(R) gravity preserve the (hamiltonian and momentum) constraints. In this paper we show that this is indeed the case. In addition, we show that the mathematical form of the constraint propagation equations in BD-equilavent f(R) gravity and in f(R) gravity in both the Jordan and Einstein frames, is exactly the same as in the standard ADM 3+1 decomposition of GR. Finally, we point out that current numerical relativity codes can incorporate the 3+1 evolution equations of metric f(R) gravity by modifying the stress-energy tensor and adding an additional scalar field evolution equation. We hope that this work will serve as a starting point for relativists to develop fully dynamical codes for valid f(R) models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Head-on collisions of binary white dwarf--neutron stars: Simulations in full general relativity

We simulate head-on collisions from rest at large separation of binary white dwarf -- neutron stars (WDNSs) in full general relativity. Our study serves as a prelude to our analysis of the circular binary WDNS problem. We focus on compact binaries whose total mass exceeds the maximum mass that a cold degenerate star can support, and our goal is to determine the fate of such systems. A fully general relativistic hydrodynamic computation of a realistic WDNS head-on collision is prohibitive due to the large range of dynamical time scales and length scales involved. For this reason, we construct an equation of state (EOS) which captures the main physical features of NSs while, at the same time, scales down the size of WDs. We call these scaled-down WD models "pseudo-WDs (pWDs)". Using pWDs, we can study these systems via a sequence of simulations where the size of the pWD gradually increases toward the realistic case. We perform two sets of simulations; One set studies the effects of the NS mass on the final outcome, when the pWD is kept fixed. The other set studies the effect of the pWD compaction on the final outcome, when the pWD mass and the NS are kept fixed. All simulations show that 14%-18% of the initial total rest mass escapes to infinity. All remnant masses still exceed the maximum rest mass that our cold EOS can support (1.92 solar masses), but no case leads to prompt collapse to a black hole. This outcome arises because the final configurations are hot. All cases settle into spherical, quasiequilibrium configurations consisting of a cold NS core surrounded by a hot mantle, resembling Thorne-Zytkow objects. Extrapolating our results to realistic WD compactions, we predict that the likely outcome of a head-on collision of a realistic, massive WDNS system will be the formation of a quasiequilibrium Thorne-Zytkow-like object.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Oppenheimer-Snyder Collapse in Moving-Puncture Coordinates

Moving-puncture coordinates are commonly used in numerical simulations of black holes. Their properties for vacuum Schwarzschild black holes have been analyzed in a number of studies. The behavior of moving-puncture coordinates in spacetimes containing matter, however, is less well understood. In this paper we explore the behavior of these coordinates for Oppenheimer-Snyder collapse, i.e., the collapse of a uniform density, pressureless sphere of dust initially at rest to a black hole. Oppenheimer-Snyder collapse provides a stringent test of the singularity-avoiding properties of moving-puncture coordinates, since the singularity can form more quickly than it would for matter with pressure. Our results include analytical expressions for the matter density, lapse function, and mean curvature at early times, as well as interesting limits for later times. We also carry out numerical simulations to obtain the full solution and these show that, even in the absence of pressure, moving-puncture coordinates are able to avoid the singularity. At late times the geometry settles down to a trumpet slice of a vacuum black hole.

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

Binary Black Hole Mergers in Gaseous Environments: "Binary Bondi" and "Binary Bondi-Hoyle-Lyttleton" Accretion

Merging supermassive black hole-black hole (BHBH) binaries produced in galaxy mergers are promising sources of detectable gravitational waves. If such a merger takes place in a gaseous environment, there is a possibility of a simultaneous detection of electromagnetic and gravitational radiation, as the stirring, shock heating and accretion of the gas may produce variability and enhancements in the electromagnetic flux. Such a simultaneous detection can provide a wealth of opportunities to study gravitational physics, accretion physics, and cosmology. We investigate this scenario by performing fully general relativistic, hydrodynamic simulations of merging, equal-mass, nonspinning BHBH binaries embedded in gas clouds. We evolve the metric using the BSSN formulation with standard moving puncture gauge conditions and handle the hydrodynamics via a high-resolution shock-capturing (HRSC) scheme. We consider both "binary Bondi accretion" in which the binary is at rest relative to the ambient gas cloud, as well as "binary Bondi-Hoyle-Lyttleton accretion" in which the binary moves relative to the gas cloud. The gas cloud is assumed to be homogeneous far from the binary and governed by a Γ-law equation of state. We vary Γbetween 4/3 and 5/3. For each simulation, we compute the gas flow and accretion rate and estimate the electromagnetic luminosity due to bremsstrahlung and synchrotron emission. We find evidence for significant enhancements in both the accretion rate and luminosity over values for a single black hole of the same mass as the binary. We estimate that this luminosity enhancement should be detectable by LSST for a 10^6 M_sun binary in a hot gas cloud of density n~10/cm^3 and temperature T~10^6 K at z=1, reaching a maximum of L~3x10^43 erg/s, with the emission peaking in the visible band.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Filling the disk hollow following binary black hole merger: The transient accretion afterglow

Tidal torques from a binary black hole (BHBH) empty out the central regions in any circumbinary gaseous accretion disk. The balance between tidal torques and viscosity maintain the inner edge of the disk at a radius r ~ 1.5a -- 2a, where a is the binary semimajor axis. Eventually, the inspiraling binary decouples from disk and merges, leaving behind a central hollow ("donut hole") in the disk orbiting the remnant black hole. We present a simple, time-dependent, Newtonian calculation that follows the secular (viscous) evolution of the disk as it fills up the hollow down to the black hole innermost stable circular orbit and then relaxes to stationary equilibrium. We use our model to calculate the electromagnetic radiation ("afterglow") spectrum emitted during this transient accretion epoch. Observing the temporal increase in the total electromagnetic flux and the hardening of the spectrum as the donut hole fills may help confirm a BHBH merger detected by a gravitational wave interferometer. We show how the very existence of the initial hollow can lead to super-Eddington accretion during this secular phase if the rate is not very far below Eddington prior to decoupling. Our model, though highly idealized, may be useful in establishing some of the key parameters, thermal emission features and scalings that characterize this transient. It can serve as a guide in the design and calibration of future radiation-magnetohydrodynamic simulations in general relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Merger of white dwarf-neutron star binaries: Prelude to hydrodynamic simulations in general relativity

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.

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
preprint2007arXiv

Magnetized Hypermassive Neutron Star Collapse: a candidate central engine for short-hard GRBs

Hypermassive neutron stars (HMNSs) are equilibrium configurations supported against collapse by rapid differential rotation and likely form as transient remnants of binary neutron star mergers. Though HMNSs are dynamically stable, secular effects such as viscosity or magnetic fields tend to bring HMNSs into uniform rotation and thus lead to collapse. We simulate the evolution of magnetized HMNSs in axisymmetry using codes which solve the Einstein-Maxwell-MHD system of equations. We find that magnetic braking and the magnetorotational instability (MRI) both contribute to the eventual collapse of HMNSs to rotating black holes surrounded by massive, hot accretion tori and collimated magnetic fields. Such hot tori radiate strongly in neutrinos, and the resulting neutrino-antineutrino annihilation could power short-hard GRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The Final Fate of Binary Neutron Stars: What Happens After the Merger?

The merger of two neutron stars usually produces a remnant with a mass significantly above the single (nonrotating) neutron star maximum mass. In some cases, the remnant will be stabilized against collapse by rapid, differential rotation. MHD-driven angular momentum transport eventually leads to the collapse of the remnant's core, resulting in a black hole surrounded by a massive accretion torus. Here we present simulations of this process. The plausibility of generating short duration gamma ray bursts through this scenario is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Quasi-equilibrium Binary Black Hole Initial Data for Dynamical Evolutions

We present a formalism for constructing quasi-equilibrium binary black hole initial data suitable for numerical evolution. We construct quasi-equilibrium models by imposing an approximate helical Killing symmetry appropriate for quasi-circular orbits. We use the sum of two Kerr-Schild metrics as our background metric, thereby improving on conformally flat backgrounds that do not accommodate rotating black holes and providing a horizon-penetrating lapse, convenient for implementing black hole excision. We set inner boundary conditions at an excision radius well inside the apparent horizon and construct these boundary conditions to incorporate the quasi-equilibrium condition and recover the solution for isolated black holes in the limit of large separation. We use our formalism both to generate initial data for binary black hole evolutions and to construct a crude quasi-equilibrium, inspiral sequence for binary black holes of fixed irreducible mass.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Black Hole Spin Evolution

We consider a subset of the physical processes that determine the spin j = a/M of astrophysical black holes. These include: (1) Initial conditions. Recent models suggest that the collapse of supermassive stars are likely to produce black holes with j ~ 0.7. (2) Major mergers. The outcome of a nearly equal mass black hole-black hole merger is not yet known, but we review the current best guesses and analytic bounds. (3) Minor mergers. We recover the result of Blandford & Hughes that accretion of small companions with isotropically distributed orbital angular momenta results in spindown, with j ~ M^{-7/3}. (4) Accretion. We present new results from fully relativistic magnetohydrodynamic accretion simulations. These show that, at least for one sequence of flow models, spin equilibrium (dj/dt = 0) is reached for j ~ 0.9, far less than the canonical value 0.998 of Thorne that was derived in the absence of MHD effects. This equilibrium value may not apply to all accretion flows, particularly thin disks. Nevertheless, it opens the possibility that black holes that have grown primarily through accretion are not maximally rotating.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Magnetic Braking in Differentially Rotating, Relativistic Stars

We study the magnetic braking and viscous damping of differential rotation in incompressible, uniform density stars in general relativity. Differentially rotating stars can support significantly more mass in equilibrium than nonrotating or uniformly rotating stars. The remnant of a binary neutron star merger or supernova core collapse may produce such a "hypermassive" neutron star. Although a hypermassive neutron star may be stable on a dynamical timescale, magnetic braking and viscous damping of differential rotation will ultimately alter the equilibrium structure, possibly leading to delayed catastrophic collapse. Here we consider the slow-rotation, weak-magnetic field limit in which E_rot << E_mag << W, where E_rot is the rotational kinetic energy, E_mag is the magnetic energy, and W is the gravitational binding energy of the star. We assume the system to be axisymmetric and solve the MHD equations in both Newtonian gravitation and general relativity. Toroidal magnetic fields are generated whenever the angular velocity varies along the initial poloidal field lines. We find that the toroidal fields and angular velocities oscillate independently along each poloidal field line, which enables us to transform the original 2+1 equations into 1+1 form and solve them along each field line independently. The incoherent oscillations on different field lines stir up turbulent-like motion in tens of Alfven timescales ("phase mixing"). In the presence of viscosity, the stars eventually are driven to uniform rotation, with the energy contained in the initial differential rotation going into heat. Our evolution calculations serve as qualitative guides and benchmarks for future, more realistic MHD simulations in full 3+1 general relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

One-Armed Spiral Instability in Differentially Rotating Stars

We investigate the dynamical instability of the one-armed spiral m=1 mode in differentially rotating stars by means of 3+1 hydrodynamical simulations in Newtonian gravitation. We find that both a soft equation of state and a high degree of differential rotation in the equilibrium star are necessary to excite a dynamical m=1 mode as the dominant instability at small values of the ratio of rotational kinetic to potential energy, T/|W|. We find that this spiral mode propagates outward from its point of origin near the maximum density at the center to the surface over several central orbital periods. An unstable m=1 mode triggers a secondary m=2 bar mode of smaller amplitude, and the bar mode can excite gravitational waves. As the spiral mode propagates to the surface it weakens, simultaneously damping the emitted gravitational wave signal. This behavior is in contrast to waves triggered by a dynamical m=2 bar instability, which persist for many rotation periods and decay only after a radiation-reaction damping timescale.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Effects of Differential Rotation on the Maximum Mass of Neutron Stars

The merger of binary neutron stars is likely to lead to differentially rotating remnants. In this paper we numerically construct models of differentially rotating neutron stars in general relativity and determine their maximum allowed mass. We model the stars adopting a polytropic equation of state and tabulate maximum allowed masses as a function of differential rotation and stiffness of the equation of state. We also provide a crude argument that yields a qualitative estimate of the effect of stiffness and differential rotation on the maximum allowed mass.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Black Hole Emergence in Supernovae

If a black hole formed in a core-collapse supernova is accreting material from the base of the envelope, the accretion luminosity could be observable in the supernova light curve. We present results of a fully relativistic numerical investigation of the fallback of matter onto a black hole in a supernova and examine conditions which would be favorable for detection of the black hole. In general, heating by radioactive decays is likely to prevent practical detection of the black hole, but we show that low energy explosions of more massive stars may provide an important exception. We emphasize the particular case of SN1997D in NGC1536, for which we predict that the presence of a black hole could be inferred observationally within the next year.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Coalescing Binary Neutron Stars

Coalescing compact binaries with neutron star or black hole components provide the most promising sources of gravitational radiation for detection by the LIGO/VIRGO/GEO/TAMA laser interferometers now under construction. This fact has motivated several different theoretical studies of the inspiral and hydrodynamic merging of compact binaries. Analytic analyses of the inspiral waveforms have been performed in the Post-Newtonian approximation. Analytic and numerical treatments of the coalescence waveforms from binary neutron stars have been performed using Newtonian hydrodynamics and the quadrupole radiation approximation. Numerical simulations of coalescing black hole and neutron star binaries are also underway in full general relativity. Recent results from each of these approaches will be described and their virtues and limitations summarized.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

On the Maximum Mass of Differentially Rotating Neutron Stars

We construct relativistic equilibrium models of differentially rotating neutron stars and show that they can support significantly more mass than their nonrotating or uniformly rotating counterparts. We dynamically evolve such ``hypermassive'' models in full general relativity and show that there do exist configurations which are dynamically stable against radial collapse and bar formation. Our results suggest that the remnant of binary neutron star coalescence may be temporarily stabilized by differential rotation, leading to delayed collapse and a delayed gravitational wave burst.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Numerical Evolution of Black Holes with a Hyperbolic Formulation of General Relativity

We describe a numerical code that solves Einstein's equations for a Schwarzschild black hole in spherical symmetry, using a hyperbolic formulation introduced by Choquet-Bruhat and York. This is the first time this formulation has been used to evolve a numerical spacetime containing a black hole. We excise the hole from the computational grid in order to avoid the central singularity. We describe in detail a causal differencing method that should allow one to stably evolve a hyperbolic system of equations in three spatial dimensions with an arbitrary shift vector, to second-order accuracy in both space and time. We demonstrate the success of this method in the spherically symmetric case.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

ON BLUE STRAGGLER FORMATION BY DIRECT COLLISIONS OF MAIN-SEQUENCE STARS

We report the results of new SPH calculations of parabolic collisions between main-sequence (MS) stars. The stars are assumed to be close to the MS turn-off point in a globular cluster and are therefore modeled as $n=3$, $Γ=5/3$ polytropes. We find that the high degree of central mass concentration in these stars has a profound effect on the hydrodynamics. In particular, very little hydrodynamic mixing occurs between the dense, helium-rich inner cores and the outer envelopes. As a result, and in contrast to what has been assumed in previous studies, blue stragglers formed by direct stellar collisions are not necessarily expected to have anomalously high helium abundances in their envelopes or to have their cores replenished with fresh hydrogen fuel.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Collapse to Black Holes in Brans-Dicke Theory: II. Comparison with General Relativity

We discuss a number of long-standing theoretical questions about collapse to black holes in the Brans-Dicke theory of gravitation. Using a new numerical code, we show that Oppenheimer-Snyder collapse in this theory produces black holes that are identical to those of general relativity in final equilibrium, but are quite different from those of general relativity during dynamical evolution. We find that there are epochs during which the apparent horizon of such a black hole passes {\it outside\/} the event horizon, and that the surface area of the event horizon {\it decreases\/} with time. This behavior is possible because theorems which prove otherwise assume $R_{ab}l^al^b \ge 0$ for all null vectors $l^a$. We show that dynamical spacetimes in Brans-Dicke theory can violate this inequality, even in vacuum, for any value of $ω$.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Disk collapse in general relativity

The radial collapse of a homogeneous disk of collisionless particles can be solved analytically in Newtonian gravitation. To solve the problem in general relativity, however, requires the full machinery of numerical relativity. The collapse of a disk is the simplest problem that exhibits the two most significant and challenging features of strong-field gravitation: black hole formation and gravitational wave generation. We carry out dynamical calculations of several different relativistic disk systems. We explore the growth of ring instabilities in equilibrium disks, and how they are suppressed by sufficient velocity dispersion. We calculate wave forms from oscillating disks, and from disks that undergo gravitational collapse to black holes. Studies of disk collapse to black holes should also be useful for developing new techniques for numerical relativity, such as apparent horizon boundary conditions for black hole spacetimes.

0 citations0 reviewsNo code linkedOpen access