Showing 1–36 of 36 results for author: Moesta, P

Outflows from Short-Lived Neutron-Star Merger Remnants Can Produce a Blue Kilonova

Authors: Sanjana Curtis, Pablo Bosch, Philipp Mösta, David Radice, Sebastiano Bernuzzi, Albino Perego, Roland Haas, Erik Schnetter

Abstract: We present a 3D general-relativistic magnetohydrodynamic simulation of a short-lived neutron star remnant formed in the aftermath of a binary neutron star merger. The simulation uses an M1 neutrino transport scheme to track neutrino-matter interactions and is well-suited to studying the resulting nucleosynthesis and kilonova emission. We find that the ejecta in our simulations under-produce $r$-pr… ▽ More We present a 3D general-relativistic magnetohydrodynamic simulation of a short-lived neutron star remnant formed in the aftermath of a binary neutron star merger. The simulation uses an M1 neutrino transport scheme to track neutrino-matter interactions and is well-suited to studying the resulting nucleosynthesis and kilonova emission. We find that the ejecta in our simulations under-produce $r$-process abundances beyond the second $r$-process peak. For sufficiently long-lived remnants, these outflows \textit{alone} can produce blue kilonovae, including the blue kilonova component observed for AT2017gfo. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: 9 pages, 5 figures, submitted to ApJL

Density Profiles of Collapsed Rotating Massive Stars Favor Long Gamma-Ray Bursts

Authors: Goni Halevi, Belinda Wu, Philipp Moesta, Ore Gottlieb, Alexander Tchekhovskoy, David R. Aguilera-Dena

Abstract: Long-duration gamma-ray bursts (lGRBs) originate in relativistic collimated outflows -- jets -- that drill their way out of collapsing massive stars. Accurately modeling this process requires realistic stellar profiles for the jets to propagate through and break out of. Most previous studies have used simple power laws or pre-collapse models for massive stars. However, the relevant stellar profile… ▽ More Long-duration gamma-ray bursts (lGRBs) originate in relativistic collimated outflows -- jets -- that drill their way out of collapsing massive stars. Accurately modeling this process requires realistic stellar profiles for the jets to propagate through and break out of. Most previous studies have used simple power laws or pre-collapse models for massive stars. However, the relevant stellar profile for lGRB models is in fact that of a star after its core has collapsed to form a compact object. To self-consistently compute such a stellar profile, we use the open-source code GR1D to simulate the core-collapse process for a suite of low-metallicity, rotating, massive stellar progenitors that have undergone chemically homogeneous evolution. Our models span a range of zero-age main sequence (ZAMS) masses: $M_\mathrm{ZAMS} = 13, 18, 21, 25, 35, 40$, and $45 M_\odot$. All of these models, at the onset of core-collapse, feature steep density profiles, $ρ\propto r^{-α}$ with $α\approx 2.5$, which would result in jets that are inconsistent with lGRB observables. We follow the collapse of four out of our seven models until they form BHs and the other three proto-neutron stars (PNSs). We find, across all models, that the density profile outside of the newly-formed BH or PNS is well-represented by a flatter power law with $α\approx 1.35{-}1.55$. Such flat density profiles are conducive to successful formation and breakout of BH-powered jets and, in fact, required to reproduce observable properties of lGRBs. Future models of lGRBs should be initialized with shallower \textit{post-collapse} stellar profiles like those presented here instead of the much steeper pre-collapse profiles that are typically used. △ Less

Submitted 1 February, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 10 pages, 4 figures+1 table, accepted to ApJL

GRaM-X: A new GPU-accelerated dynamical spacetime GRMHD code for Exascale computing with the Einstein Toolkit

Authors: Swapnil Shankar, Philipp Mösta, Steven R. Brandt, Roland Haas, Erik Schnetter, Yannick de Graaf

Abstract: We present GRaM-X (General Relativistic accelerated Magnetohydrodynamics on AMReX), a new GPU-accelerated dynamical-spacetime general relativistic magnetohydrodynamics (GRMHD) code which extends the GRMHD capability of Einstein Toolkit to GPU-based exascale systems. GRaM-X supports 3D adaptive mesh refinement (AMR) on GPUs via a new AMR driver for the Einstein Toolkit called CarpetX which in turn… ▽ More We present GRaM-X (General Relativistic accelerated Magnetohydrodynamics on AMReX), a new GPU-accelerated dynamical-spacetime general relativistic magnetohydrodynamics (GRMHD) code which extends the GRMHD capability of Einstein Toolkit to GPU-based exascale systems. GRaM-X supports 3D adaptive mesh refinement (AMR) on GPUs via a new AMR driver for the Einstein Toolkit called CarpetX which in turn leverages AMReX, an AMR library developed for use by the United States DOE's Exascale Computing Project (ECP). We use the Z4c formalism to evolve the equations of GR and the Valencia formulation to evolve the equations of GRMHD. GRaM-X supports both analytic as well as tabulated equations of state. We implement TVD and WENO reconstruction methods as well as the HLLE Riemann solver. We test the accuracy of the code using a range of tests on static spacetime, e.g. 1D MHD shocktubes, the 2D magnetic rotor and a cylindrical explosion, as well as on dynamical spacetimes, i.e. the oscillations of a 3D TOV star. We find excellent agreement with analytic results and results of other codes reported in literature. We also perform scaling tests and find that GRaM-X shows a weak scaling efficiency of $\sim 40-50\%$ on 2304 nodes (13824 NVIDIA V100 GPUs) with respect to single-node performance on OLCF's supercomputer Summit. △ Less

Submitted 21 November, 2022; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 22 pages, 8 figures, to be submitted to Classical and Quantum Gravity

Propagation of a realistic magnetar jet through binary neutron star merger medium and implications for short gamma-ray bursts

Authors: Gustavo Soares, Pablo Bosch, Davide Lazzati, Philipp Mösta

Abstract: The origin of short gamma-ray bursts (sGRBs) is associated with outflows powered by the remnant of a binary neutron star merger. This remnant can be either a black hole or a highly magnetized, fastly spinning neutron star, also known as a magnetar. Here, we present the results of two relativistic magnetohydrodynamical (RMHD) simulations aimed at investigating the large-scale dynamics and propagati… ▽ More The origin of short gamma-ray bursts (sGRBs) is associated with outflows powered by the remnant of a binary neutron star merger. This remnant can be either a black hole or a highly magnetized, fastly spinning neutron star, also known as a magnetar. Here, we present the results of two relativistic magnetohydrodynamical (RMHD) simulations aimed at investigating the large-scale dynamics and propagation of magnetar collimated outflows through the medium surrounding the remnant. The first simulation evolves a realistic jet by injecting external simulation data, while the second evolves an analytical model jet with similar properties for comparison. We find that both outflows remain collimated and successfully emerge through the static medium surrounding the remnant. However, they fail to attain relativistic velocities and only reach a mean maximum speed of ~0.7c for the realistic jet, and ~0.6c for the analytical jet. We also find that the realistic jet has a much more complex structure. The lack of highly relativistic speeds, that makes these jets unsuitable as short GRB sources, is due to numerical limitations and not general to all possible magnetar outflows. A jet like the one we study, however, could give rise to or augment a blue kilonova component. In addition, it would make the propagation of a relativistic jet easier, should one be launched after the neutron star collapses into a black hole. △ Less

Submitted 29 June, 2023; v1 submitted 15 September, 2022; originally announced September 2022.

Comments: 11 pages, 7 figures, accepted to ApJ

Magnetic field effects on nucleosynthesis and kilonovae from neutron star merger remnants

Authors: Sebastiaan de Haas, Pablo Bosch, Philipp Mösta, Sanjana Curtis, Nathanyel Schut

Abstract: We investigate the influence of parametric magnetic field configurations of a hypermassive neutron star (HMNS) on electromagnetic (EM) observables, specifically the kilonova lightcurves and nucleosynthesis yields. We perform three-dimensional (3D) dynamical-spacetime general-relativistic magnetohydrodynamic (GRMHD) simulations, including a neutrino leakage scheme, microphysical finite-temperature… ▽ More We investigate the influence of parametric magnetic field configurations of a hypermassive neutron star (HMNS) on electromagnetic (EM) observables, specifically the kilonova lightcurves and nucleosynthesis yields. We perform three-dimensional (3D) dynamical-spacetime general-relativistic magnetohydrodynamic (GRMHD) simulations, including a neutrino leakage scheme, microphysical finite-temperature equation of state (EOS), and an initial poloidal magnetic field. We find that varying the magnetic field strength and falloff impacts the formation of magnetized winds or mildy-relativistic jets, which in turn has profound effects on the outflow properties. All of the evolved configurations collapse to a black hole (BH) $\sim 21-23$ ms after the onset of the simulations, however, the ones forming jets may be considerably more effective at transporting angular momentum out of the system, resulting in earlier collapse times. Larger mass ejecta rates and radial velocities of unbound material characterise the systems that form jets. The bolometric light curves of the kilonovae and $r$-process yields change considerably with different magnetic field parameters. We conclude that the magnetic field strength and falloff have robust effects on the outflow properties and electromagnetic observables. This can be particularly important as the total ejecta mass from our simulations ($\simeq 10^{-3}\;M_{\odot}$) makes the ejecta from HMNS a compelling source to power kilonova through radioactive decay of $r$-process elements. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Comments: 11 pages, 8 figures

Black hole to photosphere: 3D GRMHD simulations of collapsars reveal wobbling and hybrid composition jets

Authors: Ore Gottlieb, Matthew Liska, Alexander Tchekhovskoy, Omer Bromberg, Aretaios Lalakos, Dimitrios Giannios, Philipp Mösta

Abstract: Long-duration $γ$-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth by a black-hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over 6 orders of magnitude in space and time. The collapsing stellar… ▽ More Long-duration $γ$-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth by a black-hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over 6 orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important and launches bipolar jets. The jets reach the photosphere at $\sim10^{12}$ cm with an opening angle $θ_j\sim6^\circ$ and a Lorentz factor $Γ_j\lesssim 30$, unbinding $\gtrsim90\%$ of the star. We find that (i) the disk-jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle $θ_t\sim12^\circ$, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection $θ_j+θ_t$ suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs may be rarer than currently thought and emerge in only $\sim 0.1\%$ of supernovae Ib/c, implying that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where $\sim 10\%$ of the gas maintains magnetization $σ\gtrsim 0.1$. This indicates that both a photospheric component and reconnection may play a role in the prompt emission. △ Less

Submitted 20 June, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: For movies of the simulation, see https://oregottlieb.com/collapsar.html

$r$-process Nucleosynthesis and Kilonovae from Hypermassive Neutron Star Remnants

Authors: Sanjana Curtis, Philipp Mösta, Zhenyu Wu, David Radice, Luke Roberts, Giacomo Ricigliano, Albino Perego

Abstract: We investigate $r$-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage sch… ▽ More We investigate $r$-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network $\texttt{SkyNet}$. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction $Y_e$, peaking between $\sim$0.25-0.4 depending on the neutrino luminosity employed. We find that the resulting $r$-process abundance patterns differ from solar, with no significant production of material beyond the second $r$-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code $\texttt{SNEC}$ and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the $r$-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission. △ Less

Submitted 1 December, 2021; originally announced December 2021.

Comments: 10 pages, 13 figures

EuCAPT White Paper: Opportunities and Challenges for Theoretical Astroparticle Physics in the Next Decade

Authors: R. Alves Batista, M. A. Amin, G. Barenboim, N. Bartolo, D. Baumann, A. Bauswein, E. Bellini, D. Benisty, G. Bertone, P. Blasi, C. G. Böhmer, Ž. Bošnjak, T. Bringmann, C. Burrage, M. Bustamante, J. Calderón Bustillo, C. T. Byrnes, F. Calore, R. Catena, D. G. Cerdeño, S. S. Cerri, M. Chianese, K. Clough, A. Cole, P. Coloma , et al. (112 additional authors not shown)

Abstract: Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, und… ▽ More Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, under the coordination of the European Consortium for Astroparticle Theory (EuCAPT). Addressed to the whole astroparticle physics community, it explores upcoming theoretical opportunities and challenges for our field of research, with particular emphasis on the possible synergies among different subfields, and the prospects for solving the most fundamental open questions with multi-messenger observations. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: White paper of the European Consortium for Astroparticle Theory (EuCAPT). 135 authors, 400 endorsers, 133 pages, 1382 references

Estimating outflow masses and velocities in merger simulations: impact of r-process heating and neutrino cooling

Authors: Francois Foucart, Philipp Moesta, Teresita Ramirez, Alex James Wright, Siva Darbha, Daniel Kasen

Abstract: The determination of the mass, composition, and geometry of matter outflows in black hole-neutron star and neutron star-neutron star binaries is crucial to current efforts to model kilonovae, and to understand the role of neutron star merger in r-process nucleosynthesis. In this manuscript, we review the simple criteria currently used in merger simulations to determine whether matter is unbound an… ▽ More The determination of the mass, composition, and geometry of matter outflows in black hole-neutron star and neutron star-neutron star binaries is crucial to current efforts to model kilonovae, and to understand the role of neutron star merger in r-process nucleosynthesis. In this manuscript, we review the simple criteria currently used in merger simulations to determine whether matter is unbound and what the asymptotic velocity of ejected material will be. We then show that properly accounting for both heating and cooling during r-process nucleosynthesis is important to accurately predict the mass and kinetic energy of the outflows. These processes are also likely to be crucial to predict the fallback timescale of any bound ejecta. We derive a model for the asymptotic veloicity of unbound matter and binding energy of bound matter that accounts for both of these effects and that can easily be implemented in merger simulations. We show, however, that the detailed velocity distribution and geometry of the outflows can currently only be captured by full 3D fluid simulations of the outflows, as non-local effect ignored by the simple criteria used in merger simulations cannot be safely neglected when modeling these effects. Finally, we propose the introduction of simple source terms in the fluid equations to approximately account for heating/cooling from r-process nucleosynthesis in future seconds-long 3D simulations of merger remnants, without the explicit inclusion of out-of-nuclear statistical equilibrium reactions in the simulations. △ Less

Submitted 17 November, 2021; v1 submitted 1 September, 2021; originally announced September 2021.

Comments: Accepted by Phys.Rev.D

Proto-magnetar jets as central engines for broad-lined type Ic supernovae

Authors: Swapnil Shankar, Philipp Mösta, Jennifer Barnes, Paul C. Duffell, Daniel Kasen

Abstract: A subset of type Ic supernovae (SNe Ic), broad-lined SNe Ic (SNe Ic-bl), show unusually high kinetic energies ($\sim 10^{52}$ erg) which cannot be explained by the energy supplied by neutrinos alone. Many SNe Ic-bl have been observed in coincidence with long gamma-ray bursts (GRBs) which suggests a connection between SNe and GRBs. A small fraction of core-collapse supernovae (CCSNe) form a rapidly… ▽ More A subset of type Ic supernovae (SNe Ic), broad-lined SNe Ic (SNe Ic-bl), show unusually high kinetic energies ($\sim 10^{52}$ erg) which cannot be explained by the energy supplied by neutrinos alone. Many SNe Ic-bl have been observed in coincidence with long gamma-ray bursts (GRBs) which suggests a connection between SNe and GRBs. A small fraction of core-collapse supernovae (CCSNe) form a rapidly-rotating and strongly-magnetized protoneutron star (PNS), a proto-magnetar. Jets from such magnetars can provide the high kinetic energies observed in SNe Ic-bl and also provide the connection to GRBs. In this work we use the jetted outflow produced in a 3D CCSN simulation from a consistently formed proto-magnetar as the central engine for full-star explosion simulations. We extract a range of central engine parameters and find that the extracted engine energy is in the range of $6.231 \times 10^{51}-1.725 \times 10^{52}$ erg, the engine time-scale in the range of $0.479-1.159$ s and the engine half-opening angle in the range of $\sim 9-19^{\circ}$. Using these as central engines, we perform 2D special-relativistic (SR) hydrodynamic (HD) and radiation transfer simulations to calculate the corresponding light curves and spectra. We find that these central engine parameters successfully produce SNe Ic-bl which demonstrates that jets from proto-magnetars can be viable engines for SNe Ic-bl. We also find that only the central engines with smaller opening angles ($\sim 10^{\circ}$) form a GRB implying that GRB formation is likely associated with narrower jet outflows and Ic-bl's without GRBs may be associated with wider outflows. △ Less

Submitted 17 May, 2021; originally announced May 2021.

Comments: 13 pages, 12 figures

A magnetar engine for short GRBs and kilonovae

Authors: Philipp Mösta, David Radice, Roland Haas, Erik Schnetter, Sebastiano Bernuzzi

Abstract: We investigate the influence of magnetic fields on the evolution of binary neutron-star (BNS) merger remnants via three-dimensional (3D) dynamical-spacetime general-relativistic (GR) magnetohydrodynamic (MHD) simulations. We evolve a postmerger remnant with an initial poloidal magnetic field, resolve the magnetoturbulence driven by shear flows, and include a microphysical finite-temperature equati… ▽ More We investigate the influence of magnetic fields on the evolution of binary neutron-star (BNS) merger remnants via three-dimensional (3D) dynamical-spacetime general-relativistic (GR) magnetohydrodynamic (MHD) simulations. We evolve a postmerger remnant with an initial poloidal magnetic field, resolve the magnetoturbulence driven by shear flows, and include a microphysical finite-temperature equation of state (EOS). A neutrino leakage scheme that captures the overall energetics and lepton number exchange is also included. We find that turbulence induced by the magnetorotational instability (MRI) in the hypermassive neutron star (HMNS) amplifies magnetic field to beyond magnetar-strength ($10^{15}\, \mathrm{G}$). The ultra-strong toroidal field is able to launch a relativistic jet from the HMNS. We also find a magnetized wind that ejects neutron-rich material with a rate of $\dot{M}_{\mathrm{ej}} \simeq 1 \times10^{-1}\, \mathrm{M_{\odot}\, s^{-1}}$. The total ejecta mass in our simulation is $5\times 10^{-3}\, \mathrm{M_{\odot}}$. This makes the ejecta from the HMNS an important component in BNS mergers and a promising source of $r$-process elements that can power a kilonova. The jet from the HMNS reaches a terminal Lorentz factor of $\sim 5$ in our highest-resolution simulation. The formation of this jet is aided by neutrino-cooling preventing the accretion disk from protruding into the polar region. As neutrino pair-annihilation and radiative processes in the jet (which were not included in the simulations) will boost the Lorentz factor in the jet further, our simulations demonstrate that magnetars formed in BNS mergers are a viable engine for short gamma-ray bursts (sGRBs). △ Less

Submitted 28 October, 2020; v1 submitted 12 March, 2020; originally announced March 2020.

Comments: ApJL published version

Journal ref: ApJL, 901, L37 (2020)

Catching Element Formation In The Act

Authors: Chris L. Fryer, Frank Timmes, Aimee L. Hungerford, Aaron Couture, Fred Adams, Wako Aoki, Almudena Arcones, David Arnett, Katie Auchettl, Melina Avila, Carles Badenes, Eddie Baron, Andreas Bauswein, John Beacom, Jeff Blackmon, Stephane Blondin, Peter Bloser, Steve Boggs, Alan Boss, Terri Brandt, Eduardo Bravo, Ed Brown, Peter Brown, Steve Bruenn. Carl Budtz-Jorgensen, Eric Burns , et al. (194 additional authors not shown)

Abstract: Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-ray… ▽ More Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions. △ Less

Submitted 7 February, 2019; originally announced February 2019.

Comments: 14 pages including 3 figures

Report number: LA-UR-18-29748

Gas Flows Within Cavities of Circumbinary Discs in Eccentric Binary Protostellar Systems

Authors: Philipp Mösta, Ronald E. Taam, Paul C. Duffell

Abstract: The structure and evolution of gas flows within the cavity of a circumbinary disk (CBD) surrounding the stellar components in eccentric binaries are examined via two-dimensional hydrodynamical simulations. The degree to which gas fills the cavity between the circumstellar disks (CSDs) and the CBD is found to be greater for highly eccentric systems, in comparison to low-eccentricity systems, reflec… ▽ More The structure and evolution of gas flows within the cavity of a circumbinary disk (CBD) surrounding the stellar components in eccentric binaries are examined via two-dimensional hydrodynamical simulations. The degree to which gas fills the cavity between the circumstellar disks (CSDs) and the CBD is found to be greater for highly eccentric systems, in comparison to low-eccentricity systems, reflecting the spatial extent over which mass enters into the cavity throughout the orbit. The pattern of the gas flow in the cavity differs for eccentric binaries from that of binaries in a circular orbit. In particular, the former reveals tightly wound gas streams and figure-eight-like structures for systems characterized by eccentricities, $e \ge 0.4$, whereas the latter only reveal relatively loosely bent streams from the CBD to the CSDs. Hence, the description of the stream structures can be a probe of sufficient non-circularity of the binary orbital motion. Given that the inner edge of the CBD is not very well defined for highly eccentric systems due to the complex gas structures, it is suggested that the area of the cavity for high-sensitivity imaging observations may prove to be a more useful diagnostic for probing the effectiveness of CBD clearing in the future. △ Less

Submitted 22 April, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 8 pages, 5 figures

Journal ref: ApJL, 875, L21 (2019)

$r$-Process Nucleosynthesis from Three-dimensional Jet-driven Core-Collapse Supernovae with Magnetic Misalignments

Authors: Goni Halevi, Philipp Mösta

Abstract: We investigate $r$-process nucleosynthesis in three-dimensional general relativistic magnetohydrodynamic simulations of jet-driven supernovae resulting from rapidly rotating, strongly magnetized core-collapse. We explore the effect of misaligning the pre-collapse magnetic field with respect to the rotation axis by performing four simulations: one aligned model and models with 15, 30, and 45 degree… ▽ More We investigate $r$-process nucleosynthesis in three-dimensional general relativistic magnetohydrodynamic simulations of jet-driven supernovae resulting from rapidly rotating, strongly magnetized core-collapse. We explore the effect of misaligning the pre-collapse magnetic field with respect to the rotation axis by performing four simulations: one aligned model and models with 15, 30, and 45 degree misalignments. The simulations we present employ a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to post-bounce neutrino emission and absorption. We track the thermodynamic properties of the ejected material with Lagrangian tracer particles and analyse its composition with the nuclear reaction network SkyNet. By using different neutrino luminosities in post-processing the tracer data with SkyNet, we constrain the impact of uncertainties in neutrino luminosities. We find that, for the aligned model considered here, the use of an approximate leakage scheme results in neutrino luminosity uncertainties corresponding to a factor of 100-1000 uncertainty in the abundance of third peak $r$-process elements. Our results show that for misalignments of 30 degrees or less, $r$-process elements are robustly produced as long as neutrino luminosities are reasonably low ($\lesssim 5 \times 10^{52}$ erg s$^{-1}$). For a more extreme misalignment of 45 degrees, we find the production of $r$-process elements beyond the second peak significantly reduced. We conclude that robust $r$-process nucleosynthesis in magnetorotational supernovae requires a progenitor stellar core with a large poloidal magnetic field component that is at least moderately (within $\sim 30$ degrees) aligned with the rotation axis. △ Less

Submitted 20 March, 2018; v1 submitted 26 January, 2018; originally announced January 2018.

Comments: 10 pages, 9 figures, 1 table, submitted to MNRAS

R-process Nucleosynthesis from Three-Dimensional Magnetorotational Core-Collapse Supernovae

Authors: Philipp Mösta, Luke F. Roberts, Goni Halevi, Christian D. Ott, Jonas Lippuner, Roland Haas, Erik Schnetter

Abstract: We investigate r-process nucleosynthesis in three-dimensional (3D) general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. The simulations include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission and absorption. We track th… ▽ More We investigate r-process nucleosynthesis in three-dimensional (3D) general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. The simulations include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission and absorption. We track the composition of the ejected material using the nuclear reaction network SkyNet. Our results show that the 3D dynamics of magnetorotational core-collapse supernovae (CCSN) are important for their nucleosynthetic signature. We find that production of r-process material beyond the second peak is reduced by a factor of 100 when the magnetorotational jets produced by the rapidly rotating core undergo a kink instability. Our results indicate that 3D magnetorotationally powered CCSNe are a robust r-process source only if they are obtained by the collapse of cores with unrealistically large precollapse magnetic fields of order $10^{13}$G. Additionally, a comparison simulation that we restrict to axisymmetry, results in overly optimistic r-process production for lower magnetic field strengths. △ Less

Submitted 26 December, 2017; originally announced December 2017.

Comments: 10 pages, 9 figures, 2 tables. submitted to ApJ

Report number: YITP-17-129; LANL LA-UR-17-31278

Recovery schemes for primitive variables in general-relativistic magnetohydrodynamics

Authors: Daniel M. Siegel, Philipp Mösta, Dhruv Desai, Samantha Wu

Abstract: General-relativistic magnetohydrodynamic (GRMHD) simulations are an important tool to study a variety of astrophysical systems such as neutron star mergers, core-collapse supernovae, and accretion onto compact objects. A conservative GRMHD scheme numerically evolves a set of conservation equations for 'conserved' quantities and requires the computation of certain primitive variables at every time… ▽ More General-relativistic magnetohydrodynamic (GRMHD) simulations are an important tool to study a variety of astrophysical systems such as neutron star mergers, core-collapse supernovae, and accretion onto compact objects. A conservative GRMHD scheme numerically evolves a set of conservation equations for 'conserved' quantities and requires the computation of certain primitive variables at every time step. This recovery procedure constitutes a core part of any conservative GRMHD scheme and it is closely tied to the equation of state (EOS) of the fluid. In the quest to include nuclear physics, weak interactions, and neutrino physics, state-of-the-art GRMHD simulations employ finite-temperature, composition-dependent EOSs. While different schemes have individually been proposed, the recovery problem still remains a major source of error, failure, and inefficiency in GRMHD simulations with advanced microphysics. The strengths and weaknesses of the different schemes when compared to each other remain unclear. Here we present the first systematic comparison of various recovery schemes used in different dynamical spacetime GRMHD codes for both analytic and tabulated microphysical EOSs. We assess the schemes in terms of (i) speed, (ii) accuracy, and (iii) robustness. We find large variations among the different schemes and that there is not a single ideal scheme. While the computationally most efficient schemes are less robust, the most robust schemes are computationally less efficient. More robust schemes may require an order of magnitude more calls to the EOS, which are computationally expensive. We propose an optimal strategy of an efficient three-dimensional Newton-Raphson scheme and a slower but more robust one-dimensional scheme as a fall-back. △ Less

Submitted 24 June, 2018; v1 submitted 20 December, 2017; originally announced December 2017.

Comments: 13 pages, 3 figures, 2 tables; matches published version

Journal ref: The Astrophysical Journal 859:71 (2018)

Turbulence in Core-Collapse Supernovae

Authors: David Radice, Ernazar Abdikamalov, Christian D. Ott, Philipp Moesta, Sean M. Couch, Luke F. Roberts

Abstract: Multidimensional simulations show that non-radial, turbulent, fluid motion is a fundamental component of the core-collapse supernova (CCSN) explosion mechanism. Neutrino-driven convection, the standing accretion shock instability, and relic-perturbations from advanced stages of nuclear burning can all impact the outcome of core collapse in a qualitative and quantitative way. Here, we review the cu… ▽ More Multidimensional simulations show that non-radial, turbulent, fluid motion is a fundamental component of the core-collapse supernova (CCSN) explosion mechanism. Neutrino-driven convection, the standing accretion shock instability, and relic-perturbations from advanced stages of nuclear burning can all impact the outcome of core collapse in a qualitative and quantitative way. Here, we review the current understanding of these phenomena and their role in the explosion of massive stars. We also discuss the role of protoneutron star convection and of magnetic fields in the context of the delayed neutrino mechanism. △ Less

Submitted 23 March, 2018; v1 submitted 3 October, 2017; originally announced October 2017.

Comments: 27 pages, 8 figures. Invited review for J. Phys. G special issue: "Focus on microphysics in core-collapse supernovae: 30 years since SN1987A". Accepted version

Journal ref: J. Phys. G: Nucl. Part. Phys. 45 (2018) 053003

How loud are neutron star mergers?

Authors: Sebastiano Bernuzzi, David Radice, Christian D. Ott, Luke F. Roberts, Philipp Moesta, Filippo Galeazzi

Abstract: We present results from the first large parameter study of neutron star mergers using fully general relativistic simulations with finite-temperature microphysical equations of state and neutrino cooling. We consider equal and unequal-mass binaries drawn from the galactic population and simulate each binary with three different equations of state. Our focus is on the emission of energy and angular… ▽ More We present results from the first large parameter study of neutron star mergers using fully general relativistic simulations with finite-temperature microphysical equations of state and neutrino cooling. We consider equal and unequal-mass binaries drawn from the galactic population and simulate each binary with three different equations of state. Our focus is on the emission of energy and angular momentum in gravitational waves in the postmerger phase. We find that the emitted gravitational-wave energy in the first $\sim$$10\,\mathrm{ms}$ of the life of the resulting hypermassive neutron star (HMNS) is about twice the energy emitted over the entire inspiral history of the binary. The total radiated energy per binary mass is comparable to or larger than that of nonspinning black hole inspiral-mergers. About $0.8-2.5\%$ of the binary mass-energy is emitted at kHz frequencies in the early HMNS evolution. We find a clear dependence of the postmerger GW emission on binary configuration and equation of state and show that it can be encoded as a broad function of the binary tidal coupling constant $κ^T_2$. Our results also demonstrate that the dimensionless spin of black holes resulting from subsequent HMNS collapse are limited to $\lesssim0.7-0.8$. This may significantly impact the neutrino pair annihilation mechanism for powering short gamma-ray bursts (sGRB). △ Less

Submitted 13 July, 2016; v1 submitted 20 December, 2015; originally announced December 2015.

Comments: 6 pages, 3 figures. Waveforms available at https://zenodo.org/record/57844

Journal ref: Phys. Rev. D 94, 024023 (2016)

A large scale dynamo and magnetoturbulence in rapidly rotating core-collapse supernovae

Authors: Philipp Mösta, Christian D. Ott, David Radice, Luke F. Roberts, Erik Schnetter, Roland Haas

Abstract: Magnetohydrodynamic (MHD) turbulence is of key importance in many high-energy astrophysical systems, where MHD instabilities can amplify local magnetic field over very short time scales. Specifically, the magnetorotational instability (MRI) and dynamo action have been suggested as a mechanism to grow magnetar-strength magnetic field ($\ge10^{15} G$) and magnetorotationally power the explosion of a… ▽ More Magnetohydrodynamic (MHD) turbulence is of key importance in many high-energy astrophysical systems, where MHD instabilities can amplify local magnetic field over very short time scales. Specifically, the magnetorotational instability (MRI) and dynamo action have been suggested as a mechanism to grow magnetar-strength magnetic field ($\ge10^{15} G$) and magnetorotationally power the explosion of a rotating massive star. Such stars are progenitor candidates for type Ic-bl hypernova explosions and make up all supernovae connected to long gamma-ray bursts (GRBs). The MRI has been studied with local high-resolution shearing box simulations in 3D or with global 2D simulations, but it is an open question whether MRI-driven turbulence can result in the creation of a large-scale ordered and dynamically relevant field. Here we report results from global 3D general-relativistic magnetohydrodynamic (GRMHD) turbulence simulations and show that MRI-driven MHD turbulence in rapidly rotating protoneutron stars produces an inverse cascade of energy. We find a large-scale ordered toroidal field that is consistent with the formation of bipolar magnetorotationally driven outflows. Our results demonstrate that rapidly rotating massive stars are plausible progenitors for both type Ic-bl supernovae and long GRBs, present a viable formation scenario for magnetars, and may account for potentially magnetar-powered superluminous supernovae. △ Less

Submitted 28 March, 2018; v1 submitted 2 December, 2015; originally announced December 2015.

Comments: Author's accepted version of Nature letter DOI:10.1038/nature15755, published online, Nature, Nov 30, 2015

Journal ref: Volume 528, pp. 376-379 (2015)

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

Authors: Zachariah B. Etienne, Vasileios Paschalidis, Roland Haas, Philipp Moesta, Stuart L. Shapiro

Abstract: 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 relativ… ▽ More 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. △ Less

Submitted 18 August, 2015; v1 submitted 28 January, 2015; originally announced January 2015.

Comments: 37 pages, 6 figures, single column. Matches published version

Journal ref: Class. Quantum Grav. 32 (2015) 175009

Neutrino-driven Turbulent Convection and Standing Accretion Shock Instability in Three-Dimensional Core-Collapse Supernovae

Authors: E. Abdikamalov, C. D. Ott, D. Radice, L. F. Roberts, R. Haas, C. Reisswig, P. Moesta, H. Klion, E. Schnetter

Abstract: We conduct a series of numerical experiments into the nature of three-dimensional (3D) hydrodynamics in the postbounce stalled-shock phase of core-collapse supernovae using 3D general-relativistic hydrodynamic simulations of a $27$-$M_\odot$ progenitor star with a neutrino leakage/heating scheme. We vary the strength of neutrino heating and find three cases of 3D dynamics: (1) neutrino-driven conv… ▽ More We conduct a series of numerical experiments into the nature of three-dimensional (3D) hydrodynamics in the postbounce stalled-shock phase of core-collapse supernovae using 3D general-relativistic hydrodynamic simulations of a $27$-$M_\odot$ progenitor star with a neutrino leakage/heating scheme. We vary the strength of neutrino heating and find three cases of 3D dynamics: (1) neutrino-driven convection, (2) initially neutrino-driven convection and subsequent development of the standing accretion shock instability (SASI), (3) SASI dominated evolution. This confirms previous 3D results of Hanke et al. 2013, ApJ 770, 66 and Couch & Connor 2014, ApJ 785, 123. We carry out simulations with resolutions differing by up to a factor of $\sim$4 and demonstrate that low resolution is artificially favorable for explosion in the 3D convection-dominated case, since it decreases the efficiency of energy transport to small scales. Low resolution results in higher radial convective fluxes of energy and enthalpy, more fully buoyant mass, and stronger neutrino heating. In the SASI-dominated case, lower resolution damps SASI oscillations. In the convection-dominated case, a quasi-stationary angular kinetic energy spectrum $E(\ell)$ develops in the heating layer. Like other 3D studies, we find $E(\ell) \propto \ell^{-1}$ in the "inertial range," while theory and local simulations argue for $E(\ell) \propto \ell^{-5/3}$. We argue that current 3D simulations do not resolve the inertial range of turbulence and are affected by numerical viscosity up to the energy containing scale, creating a "bottleneck" that prevents an efficient turbulent cascade. △ Less

Submitted 6 June, 2015; v1 submitted 24 September, 2014; originally announced September 2014.

Comments: 24 pages, 15 figures. Accepted for publication in The Astrophysical Journal. Added one figure and made minor modifications to text according to suggestions from the referee

Improved Upper Limits on the Stochastic Gravitational-Wave Background from 2009-2010 LIGO and Virgo Data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, J. Aasi, B. P. Abbott, R. Abbott, T. Abbott, M. R. Abernathy, T. Accadia, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, C. Affeldt, M. Agathos, N. Aggarwal, O. D. Aguiar, A. Ain, P. Ajith, A. Alemic, B. Allen, A. Allocca, D. Amariutei, M. Andersen , et al. (824 additional authors not shown)

Abstract: Gravitational waves from a variety of sources are predicted to superpose to create a stochastic background. This background is expected to contain unique information from throughout the history of the universe that is unavailable through standard electromagnetic observations, making its study of fundamental importance to understanding the evolution of the universe. We carry out a search for the st… ▽ More Gravitational waves from a variety of sources are predicted to superpose to create a stochastic background. This background is expected to contain unique information from throughout the history of the universe that is unavailable through standard electromagnetic observations, making its study of fundamental importance to understanding the evolution of the universe. We carry out a search for the stochastic background with the latest data from LIGO and Virgo. Consistent with predictions from most stochastic gravitational-wave background models, the data display no evidence of a stochastic gravitational-wave signal. Assuming a gravitational-wave spectrum of Omega_GW(f)=Omega_alpha*(f/f_ref)^alpha, we place 95% confidence level upper limits on the energy density of the background in each of four frequency bands spanning 41.5-1726 Hz. In the frequency band of 41.5-169.25 Hz for a spectral index of alpha=0, we constrain the energy density of the stochastic background to be Omega_GW(f)<5.6x10^-6. For the 600-1000 Hz band, Omega_GW(f)<0.14*(f/900 Hz)^3, a factor of 2.5 lower than the best previously reported upper limits. We find Omega_GW(f)<1.8x10^-4 using a spectral index of zero for 170-600 Hz and Omega_GW(f)<1.0*(f/1300 Hz)^3 for 1000-1726 Hz, bands in which no previous direct limits have been placed. The limits in these four bands are the lowest direct measurements to date on the stochastic background. We discuss the implications of these results in light of the recent claim by the BICEP2 experiment of the possible evidence for inflationary gravitational waves. △ Less

Submitted 12 August, 2014; v1 submitted 17 June, 2014; originally announced June 2014.

Comments: 11 pages, 2 figures, 1 table

Journal ref: Phys. Rev. Lett. 113, 231101 (2014)

First all-sky search for continuous gravitational waves from unknown sources in binary systems

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, J. Aasi, B. P. Abbott, R. Abbott, T. Abbott, M. R. Abernathy, T. Accadia, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, C. Affeldt, M. Agathos, N. Aggarwal, O. D. Aguiar, A. Ain, P. Ajith, A. Alemic, B. Allen, A. Allocca, D. Amariutei, M. Andersen , et al. (827 additional authors not shown)

Abstract: We present the first results of an all-sky search for continuous gravitational waves from unknown spinning neutron stars in binary systems using LIGO and Virgo data. Using a specially developed analysis program, the TwoSpect algorithm, the search was carried out on data from the sixth LIGO Science Run and the second and third Virgo Science Runs. The search covers a range of frequencies from 20 Hz… ▽ More We present the first results of an all-sky search for continuous gravitational waves from unknown spinning neutron stars in binary systems using LIGO and Virgo data. Using a specially developed analysis program, the TwoSpect algorithm, the search was carried out on data from the sixth LIGO Science Run and the second and third Virgo Science Runs. The search covers a range of frequencies from 20 Hz to 520 Hz, a range of orbital periods from 2 to ~2,254 h and a frequency- and period-dependent range of frequency modulation depths from 0.277 to 100 mHz. This corresponds to a range of projected semi-major axes of the orbit from ~0.6e-3 ls to ~6,500 ls assuming the orbit of the binary is circular. While no plausible candidate gravitational wave events survive the pipeline, upper limits are set on the analyzed data. The most sensitive 95% confidence upper limit obtained on gravitational wave strain is 2.3e-24 at 217 Hz, assuming the source waves are circularly polarized. Although this search has been optimized for circular binary orbits, the upper limits obtained remain valid for orbital eccentricities as large as 0.9. In addition, upper limits are placed on continuous gravitational wave emission from the low-mass x-ray binary Scorpius X-1 between 20 Hz and 57.25 Hz. △ Less

Submitted 17 September, 2014; v1 submitted 30 May, 2014; originally announced May 2014.

Comments: 16 pages, 6 figures

Journal ref: Phys. Rev. D 90, 062010 (2014)

Search for gravitational radiation from intermediate mass black hole binaries in data from the second LIGO-Virgo joint science run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, J. Aasi, B. P. Abbott, R. Abbott, T. Abbott, M. R. Abernathy, T. Accadia, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, C. Affeldt, M. Agathos, N. Aggarwal, O. D. Aguiar, A. Ain, P. Ajith, A. Alemic, B. Allen, A. Allocca, D. Amariutei, M. Andersen , et al. (825 additional authors not shown)

Abstract: This paper reports on an unmodeled, all-sky search for gravitational waves from merging intermediate mass black hole binaries (IMBHB). The search was performed on data from the second joint science run of the LIGO and Virgo detectors (July 2009 - October 2010) and was sensitive to IMBHBs with a range up to $\sim 200$ Mpc, averaged over the possible sky positions and inclinations of the binaries wi… ▽ More This paper reports on an unmodeled, all-sky search for gravitational waves from merging intermediate mass black hole binaries (IMBHB). The search was performed on data from the second joint science run of the LIGO and Virgo detectors (July 2009 - October 2010) and was sensitive to IMBHBs with a range up to $\sim 200$ Mpc, averaged over the possible sky positions and inclinations of the binaries with respect to the line of sight. No significant candidate was found. Upper limits on the coalescence-rate density of nonspinning IMBHBs with total masses between 100 and $450 \ \mbox{M}_{\odot}$ and mass ratios between $0.25$ and $1\,$ were placed by combining this analysis with an analogous search performed on data from the first LIGO-Virgo joint science run (November 2005 - October 2007). The most stringent limit was set for systems consisting of two $88 \ \mbox{M}_{\odot}$ black holes and is equal to $0.12 \ \mbox{Mpc}^{-3} \ \mbox{Myr}^{-1}$ at the $90\%$ confidence level. This paper also presents the first estimate, for the case of an unmodeled analysis, of the impact on the search range of IMBHB spin configurations: the visible volume for IMBHBs with nonspinning components is roughly doubled for a population of IMBHBs with spins aligned with the binary's orbital angular momentum and uniformly distributed in the dimensionless spin parameter up to 0.8, whereas an analogous population with antialigned spins decreases the visible volume by $\sim 20\%\,$. △ Less

Submitted 19 June, 2014; v1 submitted 8 April, 2014; originally announced April 2014.

Comments: 15 pages, 10 figures, data for plots and archived public version at https://dcc.ligo.org/LIGO-P1300158/public, see also the public announcement at http://www.ligo.org/science/Publication-S6IMBH/index.php, updated bibliography

Report number: LIGO-P1300158

Journal ref: Phys. Rev. D 89, 122003 (2014)

Magnetorotational Core-Collapse Supernovae in Three Dimensions

Authors: P. Mösta, S. Richers, C. D. Ott, R. Haas, A. L. Piro, K. Boydstun, E. Abdikamalov, C. Reisswig, E. Schnetter

Abstract: We present results of new three-dimensional (3D) general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. These simulations are the first of their kind and include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission. Our result… ▽ More We present results of new three-dimensional (3D) general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. These simulations are the first of their kind and include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission. Our results show that the 3D dynamics of magnetorotational core-collapse supernovae are fundamentally different from what was anticipated on the basis of previous simulations in axisymmetry (2D). A strong bipolar jet that develops in a simulation constrained to 2D is crippled by a spiral instability and fizzles in full 3D. While multiple (magneto-)hydrodynamic instabilities may be present, our analysis suggests that the jet is disrupted by an m=1 kink instability of the ultra-strong toroidal field near the rotation axis. Instead of an axially symmetric jet, a completely new, previously unreported flow structure develops. Highly magnetized spiral plasma funnels expelled from the core push out the shock in polar regions, creating wide secularly expanding lobes. We observe no runaway explosion by the end of the full 3D simulation at 185 ms after bounce. At this time, the lobes have reached maximum radii of 900 km. △ Less

Submitted 5 March, 2014; originally announced March 2014.

Comments: 6 pages, 4 figures. Submitted to ApJL. Additional material and animations available at http://stellarcollapse.org/cc3dgrmhd

Journal ref: Astrophys.J. 785 (2014) L29

Implementation of an F-statistic all-sky search for continuous gravitational waves in Virgo VSR1 data

Authors: J. Aasi, B. P. Abbott, R. Abbott, T. Abbott, M. R. Abernathy, T. Accadia, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, C. Affeldt, M. Agathos, N. Aggarwal, O. D. Aguiar, A. Ain, P. Ajith, A. Alemic, B. Allen, A. Allocca, D. Amariutei, M. Andersen, R. Anderson, S. B. Anderson , et al. (826 additional authors not shown)

Abstract: We present an implementation of the $\mathcal{F}$-statistic to carry out the first search in data from the Virgo laser interferometric gravitational wave detector for periodic gravitational waves from a priori unknown, isolated rotating neutron stars. We searched a frequency $f_0$ range from 100 Hz to 1 kHz and the frequency dependent spindown $f_1$ range from… ▽ More We present an implementation of the $\mathcal{F}$-statistic to carry out the first search in data from the Virgo laser interferometric gravitational wave detector for periodic gravitational waves from a priori unknown, isolated rotating neutron stars. We searched a frequency $f_0$ range from 100 Hz to 1 kHz and the frequency dependent spindown $f_1$ range from $-1.6\,(f_0/100\,{\rm Hz}) \times 10^{-9}\,$ Hz/s to zero. A large part of this frequency - spindown space was unexplored by any of the all-sky searches published so far. Our method consisted of a coherent search over two-day periods using the $\mathcal{F}$-statistic, followed by a search for coincidences among the candidates from the two-day segments. We have introduced a number of novel techniques and algorithms that allow the use of the Fast Fourier Transform (FFT) algorithm in the coherent part of the search resulting in a fifty-fold speed-up in computation of the $\mathcal{F}$-statistic with respect to the algorithm used in the other pipelines. No significant gravitational wave signal was found. The sensitivity of the search was estimated by injecting signals into the data. In the most sensitive parts of the detector band more than 90% of signals would have been detected with dimensionless gravitational-wave amplitude greater than $5 \times 10^{-24}$. △ Less

Submitted 10 April, 2014; v1 submitted 20 February, 2014; originally announced February 2014.

Comments: 27 pages, 10 figures, submitted to CQG, corrected affiliations and autors list; Science summary of results available at http://www.ligo.org/science/Publication-VSR1CWAllSkyFstat/index.php

Report number: LIGO Document No. LIGO-P1300133

The NINJA-2 project: Detecting and characterizing gravitational waveforms modelled using numerical binary black hole simulations

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the NINJA-2 Collaboration, :, J. Aasi, B. P. Abbott, R. Abbott, T. Abbott, M. R. Abernathy, T. Accadia, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, C. Affeldt, M. Agathos, N. Aggarwal, O. D. Aguiar, A. Ain, P. Ajith, A. Alemic, B. Allen, A. Allocca , et al. (867 additional authors not shown)

Abstract: The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of… ▽ More The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\odot}+10M_{\odot}$ ($50M_{\odot}+50M_{\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF] △ Less

Submitted 5 January, 2014; originally announced January 2014.

Comments: 52 pages, 11 figures, 7 tables

Report number: LIGO-P1300199

Journal ref: Class. Quantum Grav. 31 115004, 2014

Formation and Coalescence of Cosmological Supermassive Black Hole Binaries in Supermassive Star Collapse

Authors: C. Reisswig, C. D. Ott, E. Abdikamalov, R. Haas, P. Moesta, E. Schnetter

Abstract: We study the collapse of rapidly rotating supermassive stars that may have formed in the early Universe. By self-consistently simulating the dynamics from the onset of collapse using three-dimensional general-relativistic hydrodynamics with fully dynamical spacetime evolution, we show that seed perturbations in the progenitor can lead to the formation of a system of two high-spin supermassive blac… ▽ More We study the collapse of rapidly rotating supermassive stars that may have formed in the early Universe. By self-consistently simulating the dynamics from the onset of collapse using three-dimensional general-relativistic hydrodynamics with fully dynamical spacetime evolution, we show that seed perturbations in the progenitor can lead to the formation of a system of two high-spin supermassive black holes, which inspiral and merge under the emission of powerful gravitational radiation that could be observed at redshifts z>10 with the DECIGO or Big Bang Observer gravitational-wave observatories, assuming supermassive stars in the mass range 10^4-10^6 Msol. The remnant is rapidly spinning with dimensionless spin a^*=0.9. The surrounding accretion disk contains ~10% of the initial mass. △ Less

Submitted 20 September, 2013; v1 submitted 29 April, 2013; originally announced April 2013.

Comments: 6 pages, 3 figures, 2 tables, version accepted by PRL

Three-Dimensional General-Relativistic Hydrodynamic Simulations of Binary Neutron Star Coalescence and Stellar Collapse with Multipatch Grids

Authors: C. Reisswig, R. Haas, C. D. Ott, E. Abdikamalov, P. Moesta, D. Pollney, E. Schnetter

Abstract: We present a new three-dimensional general-relativistic hydrodynamic evolution scheme coupled to dynamical spacetime evolutions which is capable of efficiently simulating stellar collapse, isolated neutron stars, black hole formation, and binary neutron star coalescence. We make use of a set of adapted curvi-linear grids (multipatches) coupled with flux-conservative cell-centered adaptive mesh ref… ▽ More We present a new three-dimensional general-relativistic hydrodynamic evolution scheme coupled to dynamical spacetime evolutions which is capable of efficiently simulating stellar collapse, isolated neutron stars, black hole formation, and binary neutron star coalescence. We make use of a set of adapted curvi-linear grids (multipatches) coupled with flux-conservative cell-centered adaptive mesh refinement. This allows us to significantly enlarge our computational domains while still maintaining high resolution in the gravitational-wave extraction zone, the exterior layers of a star, or the region of mass ejection in merging neutron stars. The fluid is evolved with a high-resolution shock capturing finite volume scheme, while the spacetime geometry is evolved using fourth-order finite differences. We employ a multi-rate Runge-Kutta time integration scheme for efficiency, evolving the fluid with second-order and the spacetime geometry with fourth-order integration, respectively. We validate our code by a number of benchmark problems: a rotating stellar collapse model, an excited neutron star, neutron star collapse to a black hole, and binary neutron star coalescence. The test problems, especially the latter, greatly benefit from higher resolution in the gravitational-wave extraction zone, causally disconnected outer boundaries, and application of Cauchy-characteristic gravitational-wave extraction. We show that we are able to extract convergent gravitational-wave modes up to (l,m)=(6,6). This study paves the way for more realistic and detailed studies of compact objects and stellar collapse in full three dimensions and in large computational domains. The multipatch infrastructure and the improvements to mesh refinement and hydrodynamics codes discussed in this paper will be made available as part of the open-source Einstein Toolkit. △ Less

Submitted 5 December, 2012; originally announced December 2012.

Comments: 36 pages, 26 figures

Journal ref: PhysRevD 87 064023 2013

General-Relativistic Simulations of Three-Dimensional Core-Collapse Supernovae

Authors: C. D. Ott, E. Abdikamalov, P. Moesta, R. Haas, S. Drasco, E. O'Connor, C. Reisswig, C. Meakin, E. Schnetter

Abstract: We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27-solar-mass star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a 3-species neutrino leakage scheme. The leakage scheme captures the essential aspects… ▽ More We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27-solar-mass star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a 3-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27-solar-mass progenitor was studied in 2D by B. Mueller et al. (ApJ 761:72, 2012), who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability. △ Less

Submitted 25 March, 2013; v1 submitted 24 October, 2012; originally announced October 2012.

Comments: 24 emulateapj pages, 15 figures. Accepted for publication in the Astrophysical Journal

Implications For The Origin Of GRB 051103 From LIGO Observations

Authors: The LIGO Scientific Collaboration, J. Abadie, B. P. Abbott, T. D. Abbott, R. Abbott, M. Abernathy, C. Adams, R. Adhikari, C. Affeldt, P. Ajith, B. Allen, G. S. Allen, E. Amador Ceron, D. Amariutei, R. S. Amin, S. B. Anderson, W. G. Anderson, K. Arai, M. A. Arain, M. C. Araya, S. M. Aston, D. Atkinson, P. Aufmuth, C. Aulbert, B. E. Aylott , et al. (546 additional authors not shown)

Abstract: We present the results of a LIGO search for gravitational waves (GWs) associated with GRB 051103, a short-duration hard-spectrum gamma-ray burst (GRB) whose electromagnetically determined sky position is coincident with the spiral galaxy M81, which is 3.6 Mpc from Earth. Possible progenitors for short-hard GRBs include compact object mergers and soft gamma repeater (SGR) giant flares. A merger pro… ▽ More We present the results of a LIGO search for gravitational waves (GWs) associated with GRB 051103, a short-duration hard-spectrum gamma-ray burst (GRB) whose electromagnetically determined sky position is coincident with the spiral galaxy M81, which is 3.6 Mpc from Earth. Possible progenitors for short-hard GRBs include compact object mergers and soft gamma repeater (SGR) giant flares. A merger progenitor would produce a characteristic GW signal that should be detectable at the distance of M81, while GW emission from an SGR is not expected to be detectable at that distance. We found no evidence of a GW signal associated with GRB 051103. Assuming weakly beamed gamma-ray emission with a jet semi-angle of 30 deg we exclude a binary neutron star merger in M81 as the progenitor with a confidence of 98%. Neutron star-black hole mergers are excluded with > 99% confidence. If the event occurred in M81 our findings support the the hypothesis that GRB 051103 was due to an SGR giant flare, making it the most distant extragalactic magnetar observed to date. △ Less

Submitted 17 April, 2012; v1 submitted 20 January, 2012; originally announced January 2012.

Comments: 8 pages, 3 figures. For a repository of data used in the publication, go to: https://dcc.ligo.org/cgi-bin/DocDB/ShowDocument?docid=15166 . Also see the announcement for this paper on ligo.org at: http://www.ligo.org/science/Publication-GRB051103/index.php

Implementation and testing of the first prompt search for gravitational wave transients with electromagnetic counterparts

Authors: The LIGO Scientific Collaboration, Virgo Collaboration, J. Abadie, B. P. Abbott, R. Abbott, T. D. Abbott, M. Abernathy, T. Accadia, F. Acernese, C. Adams, R. Adhikari, C. Affeldt, P. Ajith, B. Allen, G. S. Allen, E. Amador Ceron, D. Amariutei, R. S. Amin, S. B. Anderson, W. G. Anderson, K. Arai, M. A. Arain, M. C. Araya, S. M. Aston, P. Astone , et al. (794 additional authors not shown)

Abstract: Aims. A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly ident… ▽ More Aims. A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize GW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify GW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with ~50% or better probability with a few pointings of wide-field telescopes. △ Less

Submitted 12 January, 2012; v1 submitted 15 September, 2011; originally announced September 2011.

Comments: 17 pages. This version (v2) includes two tables and 1 section not included in v1. Accepted for publication in Astronomy & Astrophysics

Report number: LIGO-P1000061

On the detectability of dual jets from binary black holes

Authors: Philipp Moesta, Daniela Alic, Luciano Rezzolla, Olindo Zanotti, Carlos Palenzuela

Abstract: We revisit the suggestion that dual jets can be produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. By performing independent calculations and by computing the electromagnetic emission in a way which is consistent with estimates using the Poynting flux, we show that a dual-jet structure is present… ▽ More We revisit the suggestion that dual jets can be produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. By performing independent calculations and by computing the electromagnetic emission in a way which is consistent with estimates using the Poynting flux, we show that a dual-jet structure is present but energetically subdominant with respect to a non-collimated and predominantly quadrupolar emission, which is similar to the one computed when the binary is in electrovacuum. While our findings set serious restrictions on the detectability of dual jets from coalescing binaries, they also increase the chances of detecting an EM counterpart from these systems. △ Less

Submitted 3 April, 2012; v1 submitted 6 September, 2011; originally announced September 2011.

Comments: 4 pages, 3 figures

Journal ref: ApJ, 749, L32, 2012

Beating the spin-down limit on gravitational wave emission from the Vela pulsar

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, J. Abadie, B. P. Abbott, R. Abbott, M. Abernathy, T. Accadia, F. Acernese, C. Adams, R. Adhikari, C. Affeldt, B. Allen, G. S. Allen, E. Amador Ceron, D. Amariutei, R. S. Amin, S. B. Anderson, W. G. Anderson, F. Antonucci, K. Arai, M. A. Arain, M. C. Araya, S. M. Aston, P. Astone, D. Atkinson , et al. (725 additional authors not shown)

Abstract: We present direct upper limits on continuous gravitational wave emission from the Vela pulsar using data from the Virgo detector's second science run. These upper limits have been obtained using three independent methods that assume the gravitational wave emission follows the radio timing. Two of the methods produce frequentist upper limits for an assumed known orientation of the star's spin axis… ▽ More We present direct upper limits on continuous gravitational wave emission from the Vela pulsar using data from the Virgo detector's second science run. These upper limits have been obtained using three independent methods that assume the gravitational wave emission follows the radio timing. Two of the methods produce frequentist upper limits for an assumed known orientation of the star's spin axis and value of the wave polarization angle of, respectively, $1.9\ee{-24}$ and $2.2\ee{-24}$, with 95% confidence. The third method, under the same hypothesis, produces a Bayesian upper limit of $2.1\ee{-24}$, with 95% degree of belief. These limits are below the indirect {\it spin-down limit} of $3.3\ee{-24}$ for the Vela pulsar, defined by the energy loss rate inferred from observed decrease in Vela's spin frequency, and correspond to a limit on the star ellipticity of $\sim 10^{-3}$. Slightly less stringent results, but still well below the spin-down limit, are obtained assuming the star's spin axis inclination and the wave polarization angles are unknown. △ Less

Submitted 15 April, 2011; v1 submitted 14 April, 2011; originally announced April 2011.

Comments: 37 pages, 6 figures corrected typo in the Authors field

Journal ref: Astrophys.J.737:93,2011

Search for Gravitational Wave Bursts from Six Magnetars

Authors: J. Abadie, B. P. Abbott, R. Abbott, M. Abernathy, T. Accadia, F. Acerneseac, C. Adams, R. Adhikari, C. Affeldt, B. Allen, G. S. Allen, E. Amador Ceron, D. Amariutei, R. S. Amin, S. B. Anderson, W. G. Anderson, F. Antonuccia, K. Arai, M. A. Arain, M. C. Araya, S. M. Aston, P. Astonea, D. Atkinson, P. Aufmuth, C. Aulbert , et al. (743 additional authors not shown)

Abstract: Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are thought to be magnetars: neutron stars powered by extreme magnetic fields. These rare objects are characterized by repeated and sometimes spectacular gamma-ray bursts. The burst mechanism might involve crustal fractures and excitation of non-radial modes which would emit gravitational waves (GWs). We present the results of a search… ▽ More Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are thought to be magnetars: neutron stars powered by extreme magnetic fields. These rare objects are characterized by repeated and sometimes spectacular gamma-ray bursts. The burst mechanism might involve crustal fractures and excitation of non-radial modes which would emit gravitational waves (GWs). We present the results of a search for GW bursts from six galactic magnetars that is sensitive to neutron star f-modes, thought to be the most efficient GW emitting oscillatory modes in compact stars. One of them, SGR 0501+4516, is likely ~1 kpc from Earth, an order of magnitude closer than magnetars targeted in previous GW searches. A second, AXP 1E 1547.0-5408, gave a burst with an estimated isotropic energy >10^{44} erg which is comparable to the giant flares. We find no evidence of GWs associated with a sample of 1279 electromagnetic triggers from six magnetars occurring between November 2006 and June 2009, in GW data from the LIGO, Virgo, and GEO600 detectors. Our lowest model-dependent GW emission energy upper limits for band- and time-limited white noise bursts in the detector sensitive band, and for f-mode ringdowns (at 1090 Hz), are 3.0x10^{44} d_1^2 erg and 1.4x10^{47} d_1^2 erg respectively, where d_1 = d_{0501} / 1 kpc and d_{0501} is the distance to SGR 0501+4516. These limits on GW emission from f-modes are an order of magnitude lower than any previous, and approach the range of electromagnetic energies seen in SGR giant flares for the first time. △ Less

Submitted 15 April, 2011; v1 submitted 17 November, 2010; originally announced November 2010.

Comments: 9 pages, 3 figures; v2 minor clarifications and new references

Journal ref: Astrophys.J.734:L35,2011

Vacuum Electromagnetic Counterparts of Binary Black-Hole Mergers

Authors: Philipp Mösta, Carlos Palenzuela, Luciano Rezzolla, Luis Lehner, Shin'ichirou Yoshida, Denis Pollney

Abstract: As one step towards a systematic modeling of the electromagnetic (EM) emission from an inspiralling black hole binary we consider a simple scenario in which the binary moves in a uniform magnetic field anchored to a distant circumbinary disc. We study this system by solving the Einstein-Maxwell equations in which the EM fields are chosen with astrophysically consistent strengths. We consider bin… ▽ More As one step towards a systematic modeling of the electromagnetic (EM) emission from an inspiralling black hole binary we consider a simple scenario in which the binary moves in a uniform magnetic field anchored to a distant circumbinary disc. We study this system by solving the Einstein-Maxwell equations in which the EM fields are chosen with astrophysically consistent strengths. We consider binaries with spins aligned or anti-aligned with the orbital angular momentum and study the dependence of gravitational and EM signals with these spin configurations. Overall we find that the EM radiation in the lowest l=2, m=2 multipole accurately reflects the gravitational one, with identical phase evolutions and amplitudes that differ only by a scaling factor. We also compute the efficiency of the energy emission in EM waves and find that it is given by E^rad_EM/M ~ 10^-15 (M/10^8 M_Sun)^2 (B/10^4 G)^2, hence 13 orders of magnitude smaller than the gravitational energy for realistic magnetic fields. The corresponding luminosity is much smaller than the accretion luminosity if the system is accreting at near the Eddington rate. Most importantly, this EM emission is at frequencies of 10^-4 (10^8 M_Sun/M) Hz, well outside those accessible to astronomical radio observations. As a result, it is unlikely that the EM emission discussed here can be detected directly and simultaneously with the gravitational-wave one. However, indirect processes, driven by changes in the EM fields behavior could yield observable events. In particular if the accretion rate of the circumbinary disc is small and sufficiently stable over the timescale of the final inspiral, then the EM emission may be observable indirectly as it will alter the accretion rate through the magnetic torques exerted by the distorted magnetic field lines. △ Less

Submitted 11 December, 2009; originally announced December 2009.

Journal ref: Phys.Rev.D81:064017,2010