Showing 1–42 of 42 results for author: Schnetter, E

AsterX: a new open-source GPU-accelerated GRMHD code for dynamical spacetimes

Authors: Jay V. Kalinani, Liwei Ji, Lorenzo Ennoggi, Federico G. Lopez Armengol, Lucas Timotheo Sanches, Bing-Jyun Tsao, Steven R Brandt, Manuela Campanelli, Riccardo Ciolfi, Bruno Giacomazzo, Roland Haas, Erik Schnetter, Yosef Zlochower

Abstract: We present AsterX, a novel open-source, GPU-accelerated, fully general relativistic magnetohydrodynamic (GRMHD) code designed for dynamic spacetimes in 3D Cartesian coordinates, and tailored for exascale computing. We utilize block-structured adaptive mesh refinement (AMR) through CarpetX, the new driver for the Einstein Toolkit, which is built on AMReX, a software framework for massively parallel… ▽ More We present AsterX, a novel open-source, GPU-accelerated, fully general relativistic magnetohydrodynamic (GRMHD) code designed for dynamic spacetimes in 3D Cartesian coordinates, and tailored for exascale computing. We utilize block-structured adaptive mesh refinement (AMR) through CarpetX, the new driver for the Einstein Toolkit, which is built on AMReX, a software framework for massively parallel applications. AsterX employs the Valencia formulation for GRMHD, coupled with the 'Z4c' formalism for spacetime evolution, while incorporating high resolution shock capturing schemes to accurately handle the hydrodynamics. AsterX has undergone rigorous testing in both static and dynamic spacetime, demonstrating remarkable accuracy and agreement with other codes in literature. Benchmarking the code through scaling tests on OLCF's Frontier supercomputer, we demonstrate a weak scaling efficiency of about 67%-77% on 4096 nodes compared to an 8-node performance. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 41 pages, 12 figures, to be submitted to Classical and Quantum Gravity

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

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

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)

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

The Progenitor Dependence of Three-Dimensional Core-Collapse Supernovae

Authors: C. D. Ott, L. F. Roberts, A. da Silva Schneider, J. M. Fedrow, R. Haas, E. Schnetter

Abstract: We present a first study of the progenitor star dependence of the three-dimensional (3D) neutrino mechanism of core-collapse supernovae. We employ full 3D general-relativistic multi-group neutrino radiation-hydrodynamics and simulate the post-bounce evolutions of progenitors with zero-age main sequence masses of $12$, $15$, $20$, $27$, and $40\,M_\odot$. All progenitors, with the exception of the… ▽ More We present a first study of the progenitor star dependence of the three-dimensional (3D) neutrino mechanism of core-collapse supernovae. We employ full 3D general-relativistic multi-group neutrino radiation-hydrodynamics and simulate the post-bounce evolutions of progenitors with zero-age main sequence masses of $12$, $15$, $20$, $27$, and $40\,M_\odot$. All progenitors, with the exception of the $12\,M_\odot$ star, experience shock runaway by the end of their simulations. In most cases, a strongly asymmetric explosion will result. We find three qualitatively distinct evolutions that suggest a complex dependence of explosion dynamics on progenitor density structure, neutrino heating, and 3D flow. (1) Progenitors with massive cores, shallow density profiles, and high post-core-bounce accretion rates experience very strong neutrino heating and neutrino-driven turbulent convection, leading to early shock runaway. Accretion continues at a high rate, likely leading to black hole formation. (2) Intermediate progenitors experience neutrino-driven, turbulence-aided explosions triggered by the arrival of density discontinuities at the shock. These occur typically at the silicon/silicon-oxygen shell boundary. (3) Progenitors with small cores and density profiles without strong discontinuities experience shock recession and develop the 3D standing-accretion shock instability (SASI). Shock runaway ensues late, once declining accretion rate, SASI, and neutrino-driven convection create favorable conditions. These differences in explosion times and dynamics result in a non-monotonic relationship between progenitor and compact remnant mass. △ Less

Submitted 4 December, 2017; originally announced December 2017.

Comments: 8 pages, 6 figures, submitted

Report number: YITP-17-122

A Survey of High Level Frameworks in Block-Structured Adaptive Mesh Refinement Packages

Authors: Anshu Dubey, Ann Almgren, John Bell, Martin Berzins, Steve Brandt, Greg Bryan, Phillip Colella, Daniel Graves, Michael Lijewski, Frank Löffler, Brian O'Shea, Erik Schnetter, Brian Van Straalen, Klaus Weide

Abstract: Over the last decade block-structured adaptive mesh refinement (SAMR) has found increasing use in large, publicly available codes and frameworks. SAMR frameworks have evolved along different paths. Some have stayed focused on specific domain areas, others have pursued a more general functionality, providing the building blocks for a larger variety of applications. In this survey paper we examine a… ▽ More Over the last decade block-structured adaptive mesh refinement (SAMR) has found increasing use in large, publicly available codes and frameworks. SAMR frameworks have evolved along different paths. Some have stayed focused on specific domain areas, others have pursued a more general functionality, providing the building blocks for a larger variety of applications. In this survey paper we examine a representative set of SAMR packages and SAMR-based codes that have been in existence for half a decade or more, have a reasonably sized and active user base outside of their home institutions, and are publicly available. The set consists of a mix of SAMR packages and application codes that cover a broad range of scientific domains. We look at their high-level frameworks, and their approach to dealing with the advent of radical changes in hardware architecture. The codes included in this survey are BoxLib, Cactus, Chombo, Enzo, FLASH, and Uintah. △ Less

Submitted 27 October, 2016; originally announced October 2016.

Journal ref: Journal of Parallel and Distributed Computing, Volume 74, Issue 12, December 2014, Pages 3217-3227

SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics

Authors: Lawrence E. Kidder, Scott E. Field, Francois Foucart, Erik Schnetter, Saul A. Teukolsky, Andy Bohn, Nils Deppe, Peter Diener, François Hébert, Jonas Lippuner, Jonah Miller, Christian D. Ott, Mark A. Scheel, Trevor Vincent

Abstract: We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics problems such as core-collapse supernovae and binary neutron star mergers. The robustness of the discontinuous Galerkin method allows for the use of high-resoluti… ▽ More We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics problems such as core-collapse supernovae and binary neutron star mergers. The robustness of the discontinuous Galerkin method allows for the use of high-resolution shock capturing methods in regions where (relativistic) shocks are found, while exploiting high-order accuracy in smooth regions. A task-based parallelism model allows efficient use of the largest supercomputers for problems with a heterogeneous workload over disparate spatial and temporal scales. We argue that the locality and algorithmic structure of discontinuous Galerkin methods will exhibit good scalability within a task-based parallelism framework. We demonstrate the code on a wide variety of challenging benchmark problems in (non)-relativistic (magneto)-hydrodynamics. We demonstrate the code's scalability including its strong scaling on the NCSA Blue Waters supercomputer up to the machine's full capacity of 22,380 nodes using 671,400 threads. △ Less

Submitted 21 July, 2017; v1 submitted 31 August, 2016; originally announced September 2016.

Comments: 41 pages, 13 figures, and 7 tables. Ancillary data contains simulation input files

Journal ref: Journal of Computational Physics, Volume 335, 2017, Pages 84-114

General Relativistic Three-Dimensional Multi-Group Neutrino Radiation-Hydrodynamics Simulations of Core-Collapse Supernovae

Authors: Luke F. Roberts, Christian D. Ott, Roland Haas, Evan P. O'Connor, Peter Diener, Erik Schnetter

Abstract: We report on a set of long-term general-relativistic three-dimensional (3D) multi-group (energy-dependent) neutrino-radiation hydrodynamics simulations of core-collapse supernovae. We employ a full 3D two-moment scheme with the local M1 closure, three neutrino species, and 12 energy groups per species. With this, we follow the post-core-bounce evolution of the core of a nonrotating $27$-$M_\odot$… ▽ More We report on a set of long-term general-relativistic three-dimensional (3D) multi-group (energy-dependent) neutrino-radiation hydrodynamics simulations of core-collapse supernovae. We employ a full 3D two-moment scheme with the local M1 closure, three neutrino species, and 12 energy groups per species. With this, we follow the post-core-bounce evolution of the core of a nonrotating $27$-$M_\odot$ progenitor in full unconstrained 3D and in octant symmetry for $\gtrsim$$ 380\,\mathrm{ms}$. We find the development of an asymmetric runaway explosion in our unconstrained simulation. We test the resolution dependence of our results and, in agreement with previous work, find that low resolution artificially aids explosion and leads to an earlier runaway expansion of the shock. At low resolution, the octant and full 3D dynamics are qualitatively very similar, but at high resolution, only the full 3D simulation exhibits the onset of explosion. △ Less

Submitted 6 April, 2017; v1 submitted 26 April, 2016; originally announced April 2016.

Comments: Accepted to ApJ

Report number: YITP-16-54

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)

Neutrino-Driven Convection in Core-Collapse Supernovae: High-Resolution Simulations

Authors: David Radice, Christian D. Ott, Ernazar Abdikamalov, Sean M. Couch, Roland Haas, Erik Schnetter

Abstract: We present results from high-resolution semi-global simulations of neutrino-driven convection in core-collapse supernovae. We employ an idealized setup with parametrized neutrino heating/cooling and nuclear dissociation at the shock front. We study the internal dynamics of neutrino-driven convection and its role in re-distributing energy and momentum through the gain region. We find that even if b… ▽ More We present results from high-resolution semi-global simulations of neutrino-driven convection in core-collapse supernovae. We employ an idealized setup with parametrized neutrino heating/cooling and nuclear dissociation at the shock front. We study the internal dynamics of neutrino-driven convection and its role in re-distributing energy and momentum through the gain region. We find that even if buoyant plumes are able to locally transfer heat up to the shock, convection is not able to create a net positive energy flux and overcome the downwards transport of energy from the accretion flow. Turbulent convection does, however, provide a significant effective pressure support to the accretion flow as it favors the accumulation of energy, mass and momentum in the gain region. We derive an approximate equation that is able to explain and predict the shock evolution in terms of integrals of quantities such as the turbulent pressure in the gain region or the effects of non-radial motion of the fluid. We use this relation as a way to quantify the role of turbulence in the dynamics of the accretion shock. Finally, we investigate the effects of grid resolution, which we change by a factor 20 between the lowest and highest resolution. Our results show that the shallow slopes of the turbulent kinetic energy spectra reported in previous studies are a numerical artefact. Kolmogorov scaling is progressively recovered as the resolution is increased. △ Less

Submitted 22 March, 2016; v1 submitted 16 October, 2015; originally announced October 2015.

Comments: 19 pages, 14 figures. Version accepted on ApJ. An animated version of Fig. 3 is available at: https://youtu.be/yLub83WP3WA

Journal ref: The Astrophysical Journal, 820, 76 (2016)

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

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

The Transient Gravitational-Wave Sky

Authors: Nils Andersson, John Baker, Kris Belczynski, Sebastiano Bernuzzi, Emanuele Berti, Laura Cadonati, Pablo Cerda-Duran, James Clark, Marc Favata, Lee Samuel Finn, Chris Fryer, Bruno Giacomazzo, Jose Antonio Gonzalez, Martin Hendry, Ik Siong Heng, Stefan Hild, Nathan Johnson-McDaniel, Peter Kalmus, Sergei Klimenko, Shiho Kobayashi, Kostas Kokkotas, Pablo Laguna, Luis Lehner, Janna Levin, Steve Liebling , et al. (22 additional authors not shown)

Abstract: Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techn… ▽ More Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients. △ Less

Submitted 3 May, 2013; originally announced May 2013.

Comments: 33 pages, 2 figures

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

The runaway instability in general relativistic accretion disks

Authors: O. Korobkin, E. Abdikamalov, N. Stergioulas, E. Schnetter, B. Zink, S. Rosswog, C. D. Ott

Abstract: When an accretion disk falls prey to the runaway instability, a large portion of its mass is devoured by the black hole within a few dynamical times. Despite decades of effort, it is still unclear under what conditions such an instability can occur. The technically most advanced relativistic simulations to date were unable to find a clear sign for the onset of the instability. In this work, we pre… ▽ More When an accretion disk falls prey to the runaway instability, a large portion of its mass is devoured by the black hole within a few dynamical times. Despite decades of effort, it is still unclear under what conditions such an instability can occur. The technically most advanced relativistic simulations to date were unable to find a clear sign for the onset of the instability. In this work, we present three-dimensional relativistic hydrodynamics simulations of accretion disks around black holes in dynamical space-time. We focus on the configurations that are expected to be particularly prone to the development of this instability. We demonstrate, for the first time, that the fully self-consistent general relativistic evolution does indeed produce a runaway instability. △ Less

Submitted 27 April, 2013; v1 submitted 3 October, 2012; originally announced October 2012.

Comments: 5 pages, 3 figures, minor corrections to match published version in MNRAS, +link to animation

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 431, Issue 1, pp.349-354 (2013)

Correlated Gravitational Wave and Neutrino Signals from General-Relativistic Rapidly Rotating Iron Core Collapse

Authors: C. D. Ott, E. Abdikamalov, E. O'Connor, C. Reisswig, R. Haas, P. Kalmus, S. Drasco, A. Burrows, E. Schnetter

Abstract: We present results from a new set of 3D general-relativistic hydrodynamic simulations of rotating iron core collapse. We assume octant symmetry and focus on axisymmetric collapse, bounce, the early postbounce evolution, and the associated gravitational wave (GW) and neutrino signals. We employ a finite-temperature nuclear equation of state, parameterized electron capture in the collapse phase, and… ▽ More We present results from a new set of 3D general-relativistic hydrodynamic simulations of rotating iron core collapse. We assume octant symmetry and focus on axisymmetric collapse, bounce, the early postbounce evolution, and the associated gravitational wave (GW) and neutrino signals. We employ a finite-temperature nuclear equation of state, parameterized electron capture in the collapse phase, and a multi-species neutrino leakage scheme after bounce. The latter captures the important effects of deleptonization, neutrino cooling and heating and enables approximate predictions for the neutrino luminosities in the early evolution after core bounce. We consider 12-solar-mass and 40-solar-mass presupernova models and systematically study the effects of (i) rotation, (ii) progenitor structure, and (iii) postbounce neutrino leakage on dynamics, GW, and, neutrino signals. We demonstrate, that the GW signal of rapidly rotating core collapse is practically independent of progenitor mass and precollapse structure. Moreover, we show that the effects of neutrino leakage on the GW signal are strong only in nonrotating or slowly rotating models in which GW emission is not dominated by inner core dynamics. In rapidly rotating cores, core bounce of the centrifugally-deformed inner core excites the fundamental quadrupole pulsation mode of the nascent protoneutron star. The ensuing global oscillations (f~700-800 Hz) lead to pronounced oscillations in the GW signal and correlated strong variations in the rising luminosities of antineutrino and heavy-lepton neutrinos. We find these features in cores that collapse to protoneutron stars with spin periods <~ 2.5 ms and rotational energies sufficient to drive hyper-energetic core-collapse supernova explosions. Hence, joint GW + neutrino observations of a core collapse event could deliver strong evidence for or against rapid core rotation. [abridged] △ Less

Submitted 25 March, 2013; v1 submitted 2 April, 2012; originally announced April 2012.

Comments: 29 pages, 14 figures. Replaced with version matching published version

Journal ref: Phys. Rev. D. 86, 024026 (2012)

A New Monte Carlo Method for Time-Dependent Neutrino Radiation Transport

Authors: Ernazar Abdikamalov, Adam Burrows, Christian D. Ott, Frank Löffler, Evan O'Connor, Joshua C. Dolence, Erik Schnetter

Abstract: Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy to extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such a… ▽ More Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy to extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such as core-collapse supernovae. The aim of this paper is to explore Monte Carlo methods for modeling neutrino transport in core-collapse supernovae. We generalize the Implicit Monte Carlo photon transport scheme of Fleck & Cummings and gray discrete-diffusion scheme of Densmore et al. to energy-, time-, and velocity-dependent neutrino transport. Using our 1D spherically-symmetric implementation, we show that, similar to the photon transport case, the implicit scheme enables significantly larger timesteps compared with explicit time discretization, without sacrificing accuracy, while the discrete-diffusion method leads to significant speed-ups at high optical depth. Our results suggest that a combination of spectral, velocity-dependent, Implicit Monte Carlo and discrete-diffusion Monte Carlo methods represents a robust approach for use in neutrino transport calculations in core-collapse supernovae. Our velocity-dependent scheme can easily be adapted to photon transport. △ Less

Submitted 11 July, 2012; v1 submitted 13 March, 2012; originally announced March 2012.

Comments: 29 pages, 13 figures, Accepted for publication in ApJ

The Einstein Toolkit: A Community Computational Infrastructure for Relativistic Astrophysics

Authors: Frank Löffler, Joshua Faber, Eloisa Bentivegna, Tanja Bode, Peter Diener, Roland Haas, Ian Hinder, Bruno C. Mundim, Christian D. Ott, Erik Schnetter, Gabrielle Allen, Manuela Campanelli, Pablo Laguna

Abstract: We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The Toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, c… ▽ More We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The Toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, compact objects, and collapsing stars, as well as a full suite of analysis tools. The Einstein Toolkit is currently based on the Cactus Framework for high-performance computing and the Carpet adaptive mesh refinement driver. It implements spacetime evolution via the BSSN evolution system and general-relativistic hydrodynamics in a finite-volume discretization. The toolkit is under continuous development and contains many new code components that have been publicly released for the first time and are described in this article. We discuss the motivation behind the release of the toolkit, the philosophy underlying its development, and the goals of the project. A summary of the implemented numerical techniques is included, as are results of numerical test covering a variety of sample astrophysical problems. △ Less

Submitted 14 November, 2011; originally announced November 2011.

Comments: 62 pages, 20 figures

Extreme black hole simulations: collisions of unequal mass black holes and the point particle limit

Authors: Ulrich Sperhake, Vitor Cardoso, Christian D. Ott, Erik Schnetter, Helvi Witek

Abstract: Numerical relativity has seen incredible progress in the last years, and is being applied with success to a variety of physical phenomena, from gravitational-wave research and relativistic astrophysics to cosmology and high-energy physics. Here we probe the limits of current numerical setups, by studying collisions of unequal mass, non-rotating black holes of mass-ratios up to 1:100 and making con… ▽ More Numerical relativity has seen incredible progress in the last years, and is being applied with success to a variety of physical phenomena, from gravitational-wave research and relativistic astrophysics to cosmology and high-energy physics. Here we probe the limits of current numerical setups, by studying collisions of unequal mass, non-rotating black holes of mass-ratios up to 1:100 and making contact with a classical calculation in General Relativity: the infall of a point-like particle into a massive black hole. Our results agree well with the predictions coming from linearized calculations of the infall of point-like particles into non-rotating black holes. In particular, in the limit that one hole is much smaller than the other, and the infall starts from an infinite initial separation, we recover the point-particle limit. Thus, numerical relativity is able to bridge the gap between fully non-linear dynamics and linearized approximations, which may have important applications. Finally, we also comment on the "spurious" radiation content in the initial data and the linearized predictions. △ Less

Submitted 26 May, 2011; originally announced May 2011.

Comments: 7 pages, 2 pages

Dynamics and Gravitational Wave Signature of Collapsar Formation

Authors: C. D. Ott, C. Reisswig, E. Schnetter, E. O'Connor, U. Sperhake, F. Loeffler, P. Diener, E. Abdikamalov, I. Hawke, A. Burrows

Abstract: We perform 3+1 general relativistic simulations of rotating core collapse in the context of the collapsar model for long gamma-ray bursts. We employ a realistic progenitor, rotation based on results of stellar evolution calculations, and a simplified equation of state. Our simulations track self-consistently collapse, bounce, the postbounce phase, black hole formation, and the subsequent early hyp… ▽ More We perform 3+1 general relativistic simulations of rotating core collapse in the context of the collapsar model for long gamma-ray bursts. We employ a realistic progenitor, rotation based on results of stellar evolution calculations, and a simplified equation of state. Our simulations track self-consistently collapse, bounce, the postbounce phase, black hole formation, and the subsequent early hyperaccretion phase. We extract gravitational waves from the spacetime curvature and identify a unique gravitational wave signature associated with the early phase of collapsar formation. △ Less

Submitted 8 December, 2010; originally announced December 2010.

Comments: 4 pages, 4 figures, submitted to PRL

Journal ref: Phys.Rev.Lett.106:161103,2011

Gravitational Wave Extraction in Simulations of Rotating Stellar Core Collapse

Authors: Christian Reisswig, Christian D. Ott, Ulrich Sperhake, Erik Schnetter

Abstract: We perform simulations of general relativistic rotating stellar core collapse and compute the gravitational waves (GWs) emitted in the core bounce phase of three representative models via multiple techniques. The simplest technique, the quadrupole formula (QF), estimates the GW content in the spacetime from the mass quadrupole tensor. It is strictly valid only in the weak-field and slow-motion app… ▽ More We perform simulations of general relativistic rotating stellar core collapse and compute the gravitational waves (GWs) emitted in the core bounce phase of three representative models via multiple techniques. The simplest technique, the quadrupole formula (QF), estimates the GW content in the spacetime from the mass quadrupole tensor. It is strictly valid only in the weak-field and slow-motion approximation. For the first time, we apply GW extraction methods in core collapse that are fully curvature-based and valid for strongly radiating and highly relativistic sources. We employ three extraction methods computing (i) the Newman-Penrose (NP) scalar Psi_4, (ii) Regge-Wheeler-Zerilli-Moncrief (RWZM) master functions, and (iii) Cauchy-Characteristic Extraction (CCE) allowing for the extraction of GWs at future null infinity, where the spacetime is asymptotically flat and the GW content is unambiguously defined. The latter technique is the only one not suffering from residual gauge and finite-radius effects. All curvature-based methods suffer from strong non-linear drifts. We employ the fixed-frequency integration technique as a high-pass waveform filter. Using the CCE results as a benchmark, we find that finite-radius NP extraction yields results that agree nearly perfectly in phase, but differ in amplitude by ~1-7% at core bounce, depending on the model. RWZM waveforms, while in general agreeing in phase, contain spurious high-frequency noise of comparable amplitudes to those of the relatively weak GWs emitted in core collapse. We also find remarkably good agreement of the waveforms obtained from the QF with those obtained from CCE. They agree very well in phase but systematically underpredict peak amplitudes by ~5-11% which is comparable to the NP results and is within the uncertainties associated with core collapse physics. (abridged) △ Less

Submitted 8 March, 2011; v1 submitted 2 December, 2010; originally announced December 2010.

Comments: 26 pages, 10 figures, 5 tables, matches published version

Journal ref: Phys.Rev.D83:064008,2011

Stability of general-relativistic accretion disks

Authors: Oleg Korobkin, Ernazar B. Abdikamalov, Erik Schnetter, Nikolaos Stergioulas, Burkhard Zink

Abstract: Self-gravitating relativistic disks around black holes can form as transient structures in a number of astrophysical scenarios such as binary neutron star and black hole-neutron star coalescences, as well as the core-collapse of massive stars. We explore the stability of such disks against runaway and non-axisymmetric instabilities using three-dimensional hydrodynamics simulations in full general… ▽ More Self-gravitating relativistic disks around black holes can form as transient structures in a number of astrophysical scenarios such as binary neutron star and black hole-neutron star coalescences, as well as the core-collapse of massive stars. We explore the stability of such disks against runaway and non-axisymmetric instabilities using three-dimensional hydrodynamics simulations in full general relativity using the THOR code. We model the disk matter using the ideal fluid approximation with a $Γ$-law equation of state with $Γ=4/3$. We explore three disk models around non-rotating black holes with disk-to-black hole mass ratios of 0.24, 0.17 and 0.11. Due to metric blending in our initial data, all of our initial models contain an initial axisymmetric perturbation which induces radial disk oscillations. Despite these oscillations, our models do not develop the runaway instability during the first several orbital periods. Instead, all of the models develop unstable non-axisymmetric modes on a dynamical timescale. We observe two distinct types of instabilities: the Papaloizou-Pringle and the so-called intermediate type instabilities. The development of the non-axisymmetric mode with azimuthal number m = 1 is accompanied by an outspiraling motion of the black hole, which significantly amplifies the growth rate of the m = 1 mode in some cases. Overall, our simulations show that the properties of the unstable non-axisymmetric modes in our disk models are qualitatively similar to those in Newtonian theory. △ Less

Submitted 12 November, 2010; originally announced November 2010.

Comments: 30 pages, 21 figures

Journal ref: Phys.Rev.D83:043007,2011

On the frequency band of the f-mode CFS instability

Authors: Burkhard Zink, Oleg Korobkin, Erik Schnetter, Nikolaos Stergioulas

Abstract: Rapidly rotating neutron stars can be unstable to the gravitational-wave-driven CFS mechanism if they have a neutral point in the spectrum of nonaxisymmetric f-modes. We investigate the frequencies of these modes in two sequences of uniformly rotating polytropes using nonlinear simulations in full general relativity, determine the approximate locations of the neutral points, and derive limits on… ▽ More Rapidly rotating neutron stars can be unstable to the gravitational-wave-driven CFS mechanism if they have a neutral point in the spectrum of nonaxisymmetric f-modes. We investigate the frequencies of these modes in two sequences of uniformly rotating polytropes using nonlinear simulations in full general relativity, determine the approximate locations of the neutral points, and derive limits on the observable frequency band available to the instability in these sequences. We find that general relativity enhances the detectability of a CFS-unstable neutron star substantially, both by widening the instability window and enlarging the band into the optimal range for interferometric detectors like LIGO, VIRGO, and GEO-600. △ Less

Submitted 3 March, 2010; originally announced March 2010.

Comments: 13 pages, 9 figures

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

Computational Models of Stellar Collapse and Core-Collapse Supernovae

Authors: C. D. Ott, E. Schnetter, A. Burrows, E. Livne, E. O'Connor, F. Loeffler

Abstract: Core-collapse supernovae are among Nature's most energetic events. They mark the end of massive star evolution and pollute the interstellar medium with the life-enabling ashes of thermonuclear burning. Despite their importance for the evolution of galaxies and life in the universe, the details of the core-collapse supernova explosion mechanism remain in the dark and pose a daunting computational… ▽ More Core-collapse supernovae are among Nature's most energetic events. They mark the end of massive star evolution and pollute the interstellar medium with the life-enabling ashes of thermonuclear burning. Despite their importance for the evolution of galaxies and life in the universe, the details of the core-collapse supernova explosion mechanism remain in the dark and pose a daunting computational challenge. We outline the multi-dimensional, multi-scale, and multi-physics nature of the core-collapse supernova problem and discuss computational strategies and requirements for its solution. Specifically, we highlight the axisymmetric (2D) radiation-MHD code VULCAN/2D and present results obtained from the first full-2D angle-dependent neutrino radiation-hydrodynamics simulations of the post-core-bounce supernova evolution. We then go on to discuss the new code Zelmani which is based on the open-source HPC Cactus framework and provides a scalable AMR approach for 3D fully general-relativistic modeling of stellar collapse, core-collapse supernovae and black hole formation on current and future massively-parallel HPC systems. We show Zelmani's scaling properties to more than 16,000 compute cores and discuss first 3D general-relativistic core-collapse results. △ Less

Submitted 23 July, 2009; originally announced July 2009.

Comments: 16 pages, 5 figures, to appear in the proceedings of the DOE/SciDAC 2009 conference. A version with high-resolution figures is available from http://stellarcollapse.org/papers/Ott_SciDAC2009.pdf

Multi-patch methods in general relativistic astrophysics - I. Hydrodynamical flows on fixed backgrounds

Authors: Burkhard Zink, Erik Schnetter, Manuel Tiglio

Abstract: Many systems of interest in general relativistic astrophysics, including neutron stars, accreting compact objects in X-ray binaries and active galactic nuclei, core collapse, and collapsars, are assumed to be approximately spherically symmetric or axisymmetric. In Newtonian or fixed-background relativistic approximations it is common practice to use spherical polar coordinates for computational… ▽ More Many systems of interest in general relativistic astrophysics, including neutron stars, accreting compact objects in X-ray binaries and active galactic nuclei, core collapse, and collapsars, are assumed to be approximately spherically symmetric or axisymmetric. In Newtonian or fixed-background relativistic approximations it is common practice to use spherical polar coordinates for computational grids; however, these coordinates have singularities and are difficult to use in fully relativistic models. We present, in this series of papers, a numerical technique which is able to use effectively spherical grids by employing multiple patches. We provide detailed instructions on how to implement such a scheme, and present a number of code tests for the fixed background case, including an accretion torus around a black hole. △ Less

Submitted 3 December, 2007; originally announced December 2007.

Comments: 26 pages, 20 figures. A high-resolution version is available at http://www.cct.lsu.edu/~bzink/papers/multipatch_1.pdf

Journal ref: Phys.Rev.D77:103015,2008

The final spin from the coalescence of aligned-spin black-hole binaries

Authors: Luciano Rezzolla, Peter Diener, Ernst Nils Dorband, Denis Pollney, Christian Reisswig, Erik Schnetter, Jennifer Seiler

Abstract: Determining the final spin of a black-hole (BH) binary is a question of key importance in astrophysics. Modelling this quantity in general is made difficult by the fact that it depends on the 7-dimensional space of parameters characterizing the two initial black holes. However, in special cases, when symmetries can be exploited, the description can become simpler. For black-hole binaries with un… ▽ More Determining the final spin of a black-hole (BH) binary is a question of key importance in astrophysics. Modelling this quantity in general is made difficult by the fact that it depends on the 7-dimensional space of parameters characterizing the two initial black holes. However, in special cases, when symmetries can be exploited, the description can become simpler. For black-hole binaries with unequal masses but with equal spins which are aligned with the orbital angular momentum, we show that the use of recent simulations and basic but exact constraints derived from the extreme mass-ratio limit allow to model this quantity with a simple analytic expression. Despite the simple dependence, the expression models very accurately all of the available estimates, with errors of a couple of percent at most. We also discuss how to use the fit to predict when a Schwarzschild BH is produced by the merger of two spinning BHs, when the total angular momentum of the spacetime ``flips'' sign, or under what conditions the final BH is ``spun-up'' by the merger. Finally, suggest an extension of the fit to include unequal-spin binaries, thus potentially providing a complete description of the final spin from the coalescence of generic black-hole binaries with spins aligned to the orbital angular momentum. △ Less

Submitted 2 February, 2008; v1 submitted 17 October, 2007; originally announced October 2007.

Comments: Version matching the published one; small changes throughout to fit space constraints; corrects error in vii) about spin-up/down

Journal ref: Astrophys.J.674:L29-L32,2008

Spin Diagrams for Equal-Mass Black-Hole Binaries with Aligned Spins

Authors: Luciano Rezzolla, Ernst Nils Dorband, Christian Reisswig, Peter Diener, Denis Pollney, Erik Schnetter, Bela Szilagyi

Abstract: Binary black-hole systems with spins aligned with the orbital angular momentum are of special interest as they may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with spins aligned with the orbital angular momentum but otherwise arbit… ▽ More Binary black-hole systems with spins aligned with the orbital angular momentum are of special interest as they may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with spins aligned with the orbital angular momentum but otherwise arbitrary. By least-square fitting the results of these simulations we have constructed two "spin diagrams" which provide straightforward information about the recoil velocity |v_kick| and the final black-hole spin a_fin in terms of the dimensionless spins a_1 and a_2 of the two initial black holes. Overall they suggest a maximum recoil velocity of |v_kick|=441.94 km/s, and minimum and maximum final spins a_fin=0.3471 and a_fin=0.9591, respectively. △ Less

Submitted 23 February, 2008; v1 submitted 29 August, 2007; originally announced August 2007.

Comments: 4 pages, 3 figs; small changes matching published version

Journal ref: Astrophys.J. 679:1422-1426,2008

Recoil velocities from equal-mass binary black-hole mergers: a systematic investigation of spin-orbit aligned configurations

Authors: Denis Pollney, Christian Reisswig, Luciano Rezzolla, Bela Szilagyi, Marcus Ansorg, Barrett Deris, Peter Diener, Ernst Nils Dorband, Michael Koppitz, Alessandro Nagar, Erik Schnetter

Abstract: Binary black-hole systems with spins aligned with the orbital angular momentum are of special interest, as studies indicate that this configuration is preferred in nature. If the spins of the two bodies differ, there can be a prominent beaming of the gravitational radiation during the late plunge, causing a recoil of the final merged black hole. We perform an accurate and systematic study of rec… ▽ More Binary black-hole systems with spins aligned with the orbital angular momentum are of special interest, as studies indicate that this configuration is preferred in nature. If the spins of the two bodies differ, there can be a prominent beaming of the gravitational radiation during the late plunge, causing a recoil of the final merged black hole. We perform an accurate and systematic study of recoil velocities from a sequence of equal-mass black holes whose spins are aligned with the orbital angular momentum, and whose individual spins range from a = +0.584 to -0.584. In this way we extend and refine the results of a previous study and arrive at a consistent maximum recoil of 448 +- 5 km/s for anti-aligned models as well as to a phenomenological expression for the recoil velocity as a function of spin ratio. This relation highlights a nonlinear behavior, not predicted by the PN estimates, and can be readily employed in astrophysical studies on the evolution of binary black holes in massive galaxies. An essential result of our analysis is the identification of different stages in the waveform, including a transient due to lack of an initial linear momentum in the initial data. Furthermore we are able to identify a pair of terms which are largely responsible for the kick, indicating that an accurate computation can be obtained from modes up to l=3. Finally, we provide accurate measures of the radiated energy and angular momentum, finding these to increase linearly with the spin ratio, and derive simple expressions for the final spin and the radiated angular momentum which can be easily implemented in N-body simulations of compact stellar systems. Our code is calibrated with strict convergence tests and we verify the correctness of our measurements by using multiple independent methods whenever possible. △ Less

Submitted 17 July, 2007; originally announced July 2007.

Comments: 24 pages, 15 figures, 5 tables

Journal ref: Phys.Rev.D76:124002,2007

Fragmentation of general relativistic quasi-toroidal polytropes

Authors: Burkhard Zink, Nikolaos Stergioulas, Ian Hawke, Christian D. Ott, Erik Schnetter, Ewald Mueller

Abstract: We investigate the role of rotational instabilities in the context of black hole formation in relativistic stars. In addition to the standard scenario - an axially symmetric dynamical instability forming a horizon at the star's center - the recently found low-$T/|W|$ instabilities are shown to lead to fragmentation and off-center horizon formation in differentially rotating stars. This process m… ▽ More We investigate the role of rotational instabilities in the context of black hole formation in relativistic stars. In addition to the standard scenario - an axially symmetric dynamical instability forming a horizon at the star's center - the recently found low-$T/|W|$ instabilities are shown to lead to fragmentation and off-center horizon formation in differentially rotating stars. This process might be an alternative pathway to produce SMBHs from supermassive stars with inefficient angular momentum transport. △ Less

Submitted 3 April, 2007; originally announced April 2007.

Comments: 2 pages, 1 figure. Submitted to the Proceedings of the 11th Marcel Grossmann Meeting (MG11) in Berlin, Germany, July 23-29, 2006

Rotating Collapse of Stellar Iron Cores in General Relativity

Authors: Christian D. Ott, Harald Dimmelmeier, Andreas Marek, Hans-Thomas Janka, Burkhard Zink, Ian Hawke, Erik Schnetter

Abstract: We present results from the first 2+1 and 3+1 simulations of the collapse of rotating stellar iron cores in general relativity employing a finite-temperature equation of state and an approximate treatment of deleptonization during collapse. We compare full 3+1 and conformally-flat spacetime evolution methods and find that the conformally-flat treatment is sufficiently accurate for the core-colla… ▽ More We present results from the first 2+1 and 3+1 simulations of the collapse of rotating stellar iron cores in general relativity employing a finite-temperature equation of state and an approximate treatment of deleptonization during collapse. We compare full 3+1 and conformally-flat spacetime evolution methods and find that the conformally-flat treatment is sufficiently accurate for the core-collapse supernova problem. We focus on the gravitational wave (GW) emission from rotating collapse, core bounce, and early postbounce phases. Our results indicate that the GW signature of these phases is much more generic than previously estimated. In addition, we track the growth of a nonaxisymmetric instability of dominant m = 1 character in two of our models that leads to prolonged narrow-band GW emission at ~ 920-930 Hz over several tens of milliseconds. △ Less

Submitted 3 April, 2007; v1 submitted 21 December, 2006; originally announced December 2006.

Comments: 19 pages, 6 figures, minor corrections and additions; accepted for publication in the Proceedings of the New Frontiers in Numerical Relativity Conference, AEI Golm, Germany, 2006

Journal ref: Class.Quant.Grav.24:S139-S154,2007

Non-axisymmetric instability and fragmentation of general relativistic quasi-toroidal stars

Authors: Burkhard Zink, Nikolaos Stergioulas, Ian Hawke, Christian D. Ott, Erik Schnetter, Ewald Mueller

Abstract: In a recent publication, we have demonstrated that differentially rotating stars admit new channels of black hole formation via fragmentation instabilities. Since a higher order instability of this kind could potentially transform a differentially rotating supermassive star into a multiple black hole system embedded in a massive accretion disk, we investigate the dependence of the instability on… ▽ More In a recent publication, we have demonstrated that differentially rotating stars admit new channels of black hole formation via fragmentation instabilities. Since a higher order instability of this kind could potentially transform a differentially rotating supermassive star into a multiple black hole system embedded in a massive accretion disk, we investigate the dependence of the instability on parameters of the equilibrium model. We find that many of the models constructed exhibit non-axisymmetric instabilities with corotation points, even for low values of T/|W|, which lead to a fission of the stars into one, two or three fragments, depending on the initial perturbation. At least in the models selected here, an m=1 mode becomes unstable at lower values of T/|W|, which would seem to favor a scenario where one black hole with a massive accretion disk forms. In this case, we have gained evidence that low values of compactness of the initial model can lead to a stabilization of the resulting fragment, thus preventing black hole formation in this scenario. △ Less

Submitted 19 November, 2006; originally announced November 2006.

Comments: 30 pages, 36 figures. A version with high-resolution color figures can be found at http://www.cct.lsu.edu/~bzink/material/qts_instability.pdf

Journal ref: Phys.Rev.D76:024019,2007

3D Collapse of Rotating Stellar Iron Cores in General Relativity including Deleptonization and a Nuclear Equation of State

Authors: C. D. Ott, H. Dimmelmeier, A. Marek, H. -T. Janka, I. Hawke, B. Zink, E. Schnetter

Abstract: We present results from the first 2D and 3D simulations of the collapse of rotating stellar iron cores in general relativity employing a finite-temperature equation of state and an approximate treatment of deleptonization during collapse. We compare fully nonlinear and conformally flat spacetime evolution methods and find that the conformally flat treatment is sufficiently accurate for the core-… ▽ More We present results from the first 2D and 3D simulations of the collapse of rotating stellar iron cores in general relativity employing a finite-temperature equation of state and an approximate treatment of deleptonization during collapse. We compare fully nonlinear and conformally flat spacetime evolution methods and find that the conformally flat treatment is sufficiently accurate for the core-collapse supernova problem. We focus on the gravitational wave (GW) emission from rotating collapse, core bounce, and early postbounce phases. Our results indicate that the GW signature of these phases is much more generic than previously estimated. In addition, we track the growth of a nonaxisymmetric instability of dominant m = 1 character in one of our models that leads to prolonged narrow-band GW emission at ~930 Hz over several tens of milliseconds. △ Less

Submitted 1 May, 2007; v1 submitted 29 September, 2006; originally announced September 2006.

Comments: 4 pages, 4 figures, accepted for publication in Phys. Rev. Lett

Report number: AEI-2006-073

Journal ref: Phys.Rev.Lett.98:261101,2007

A numerical study of the quasinormal mode excitation of Kerr black holes

Authors: Ernst Nils Dorband, Emanuele Berti, Peter Diener, Erik Schnetter, Manuel Tiglio

Abstract: We present numerical results from three-dimensional evolutions of scalar perturbations of Kerr black holes. Our simulations make use of a high-order accurate multi-block code which naturally allows for fixed adaptivity and smooth inner (excision) and outer boundaries. We focus on the quasinormal ringing phase, presenting a systematic method for extraction of the quasinormal mode frequencies and… ▽ More We present numerical results from three-dimensional evolutions of scalar perturbations of Kerr black holes. Our simulations make use of a high-order accurate multi-block code which naturally allows for fixed adaptivity and smooth inner (excision) and outer boundaries. We focus on the quasinormal ringing phase, presenting a systematic method for extraction of the quasinormal mode frequencies and amplitudes and comparing our results against perturbation theory. The amplitude of each mode depends exponentially on the starting time of the quasinormal regime, which is not defined unambiguously. We show that this time-shift problem can be circumvented by looking at appropriately chosen relative mode amplitudes. From our simulations we extract the quasinormal frequencies and the relative and absolute amplitudes of corotating and counterrotating modes (including overtones in the corotating case). We study the dependence of these amplitudes on the shape of the initial perturbation, the angular dependence of the mode and the black hole spin, comparing against results from perturbation theory in the so-called asymptotic approximation. We also compare the quasinormal frequencies from our numerical simulations with predictions from perturbation theory, finding excellent agreement. Finally we study under what conditions the relative amplitude between given pairs of modes gets maximally excited and present a quantitative analysis of rotational mode-mode coupling. The main conclusions and techniques of our analysis are quite general and, as such, should be of interest in the study of ringdown gravitational waves produced by astrophysical gravitational wave sources. △ Less

Submitted 18 August, 2006; originally announced August 2006.

Comments: 21 pages, 9 figures

Report number: LSU-REL-081806

Journal ref: Phys.Rev.D74:084028,2006

Black hole head-on collisions and gravitational waves with fixed mesh-refinement and dynamic singularity excision

Authors: U. Sperhake, B. Kelly, P. Laguna, K. L. Smith, E. Schnetter

Abstract: We present long-term-stable and convergent evolutions of head-on black hole collisions and extraction of gravitational waves generated during the merger and subsequent ring-down. The new ingredients in this work are the use of fixed mesh-refinement and dynamical singularity excision techniques. We are able to carry out head-on collisions with large initial separations and demonstrate that our ex… ▽ More We present long-term-stable and convergent evolutions of head-on black hole collisions and extraction of gravitational waves generated during the merger and subsequent ring-down. The new ingredients in this work are the use of fixed mesh-refinement and dynamical singularity excision techniques. We are able to carry out head-on collisions with large initial separations and demonstrate that our excision infrastructure is capable of accommodating the motion of the individual black holes across the computational domain as well as their their merger. We extract gravitational waves from these simulations using the Zerilli-Moncrief formalism and find the ring-down radiation to be, as expected, dominated by the l=2, m=0 quasi-normal mode. The total radiated energy is about 0.1 % of the total ADM mass of the system. △ Less

Submitted 20 June, 2005; v1 submitted 16 March, 2005; originally announced March 2005.

Comments: Revised version, 1 figure added, accepted for publication in Phys.Rev.D, 15 pages, 10 figures, revtex 4.0

Journal ref: Phys.Rev.D71:124042,2005

Gravitational-Wave Emission from Rotating Gravitational Collapse in three Dimensions

Authors: L. Baiotti, I. Hawke, L. Rezzolla, E. Schnetter

Abstract: We present the first calculation of gravitational wave emission produced in the gravitational collapse of uniformly rotating neutron stars to black holes in fully three-dimensional simulations. The initial stellar models are relativistic polytropes which are dynamically unstable and with angular velocities ranging from slow rotation to the mass-shedding limit. An essential aspect of these simula… ▽ More We present the first calculation of gravitational wave emission produced in the gravitational collapse of uniformly rotating neutron stars to black holes in fully three-dimensional simulations. The initial stellar models are relativistic polytropes which are dynamically unstable and with angular velocities ranging from slow rotation to the mass-shedding limit. An essential aspect of these simulations is the use of progressive mesh-refinement techniques which allow to move the outer boundaries of the computational domain to regions where gravitational radiation attains its asymptotic form. The waveforms have been extracted using a gauge-invariant approach in which the numerical spacetime is matched with the non-spherical perturbations of a Schwarzschild spacetime. Overall, the results indicate that the waveforms have features related to the properties of the initial stellar models (in terms of their w-mode oscillations) and of the newly produced rotating black holes (in terms of their quasi-normal modes). While our waveforms are in good qualitative agreement with those computed by Stark and Piran in two-dimensional simulations, our amplitudes are about one order of magnitude smaller and this difference is mostly likely due to our less severe pressure reduction. For a neutron star rotating uniformly near mass-shedding and collapsing at 10 kpc, the signal-to-noise ratio computed uniquely from the burst is S/N ~ 0.25, but this grows to be S/N <~ 4 in the case of LIGO II. △ Less

Submitted 4 March, 2005; originally announced March 2005.

Comments: 4 pages, Phys. Rev. Lett., in press

Journal ref: Phys.Rev.Lett. 94 (2005) 131101

Black hole formation through fragmentation of toroidal polytropes

Authors: Burkhard Zink, Nikolaos Stergioulas, Ian Hawke, Christian D. Ott, Erik Schnetter, Ewald Mueller

Abstract: We investigate new paths to black hole formation by considering the general relativistic evolution of a differentially rotating polytrope with toroidal shape. We find that this polytrope is unstable to nonaxisymmetric modes, which leads to a fragmentation into self-gravitating, collapsing components. In the case of one such fragment, we apply a simplified adaptive mesh refinement technique to fo… ▽ More We investigate new paths to black hole formation by considering the general relativistic evolution of a differentially rotating polytrope with toroidal shape. We find that this polytrope is unstable to nonaxisymmetric modes, which leads to a fragmentation into self-gravitating, collapsing components. In the case of one such fragment, we apply a simplified adaptive mesh refinement technique to follow the evolution to the formation of an apparent horizon centered on the fragment. This is the first study of the one-armed instability in full general relativity. △ Less

Submitted 10 February, 2005; v1 submitted 26 January, 2005; originally announced January 2005.

Comments: 4 pages, 3 figures

Journal ref: Phys.Rev.Lett.96:161101,2006

Computational Relativistic Astrophysics With Adaptive Mesh Refinement: Testbeds

Authors: Edwin Evans, Sai Iyer, Erik Schnetter, Wai-Mo Suen, Jian Tao, Randy Wolfmeyer, Hui-Min Zhang

Abstract: We have carried out numerical simulations of strongly gravitating systems based on the Einstein equations coupled to the relativistic hydrodynamic equations using adaptive mesh refinement (AMR) techniques. We show AMR simulations of NS binary inspiral and coalescence carried out on a workstation having an accuracy equivalent to that of a $1025^3$ regular unigrid simulation, which is, to the best… ▽ More We have carried out numerical simulations of strongly gravitating systems based on the Einstein equations coupled to the relativistic hydrodynamic equations using adaptive mesh refinement (AMR) techniques. We show AMR simulations of NS binary inspiral and coalescence carried out on a workstation having an accuracy equivalent to that of a $1025^3$ regular unigrid simulation, which is, to the best of our knowledge, larger than all previous simulations of similar NS systems on supercomputers. We believe the capability opens new possibilities in general relativistic simulations. △ Less

Submitted 3 March, 2005; v1 submitted 20 January, 2005; originally announced January 2005.

Comments: 7 pages, 16 figures

Report number: WUGRAV Preprint 05-01

Journal ref: Phys.Rev. D71 (2005) 081301

A modified SPH approach for fluids with large density differences

Authors: Frank Ott, Erik Schnetter

Abstract: We introduce a modified SPH approach that is based on discretising the particle density instead of the mass density. This approach makes it possible to use SPH particles with very different masses to simulate multi-phase flows with large differences in mass density between the phases. We test our formulation with a simple advection problem, with sound waves encountering a density discontinuity,… ▽ More We introduce a modified SPH approach that is based on discretising the particle density instead of the mass density. This approach makes it possible to use SPH particles with very different masses to simulate multi-phase flows with large differences in mass density between the phases. We test our formulation with a simple advection problem, with sound waves encountering a density discontinuity, and with shock tubes containing a contact discontinuity between air and Diesel oil. For all examined problems where particles have different masses, the new formulation yields better results than standard SPH. This is also the case for problems in which different spatial resolutions are needed while the mass density does not change. △ Less

Submitted 20 October, 2003; v1 submitted 27 March, 2003; originally announced March 2003.

Comments: 11 pages, 4 figures. Use more stringent test cases, clarify motivation

Grazing Collisions of Black Holes via the Excision of Singularities

Authors: Steve Brandt, Randall Correll, Roberto Gomez, Mijan Huq, Pablo Laguna, Luis Lehner, Pedro Marronetti, Richard A. Matzner, David Neilsen, Jorge Pullin, Erik Schnetter, Deirdre Shoemaker, Jeffrey Winicour

Abstract: We present the first simulations of non-headon (grazing) collisions of binary black holes in which the black hole singularities have been excised from the computational domain. Initially two equal mass black holes $m$ are separated a distance $\approx10m$ and with impact parameter $\approx2m$. Initial data are based on superposed, boosted (velocity $\approx0.5c$) solutions of single black holes… ▽ More We present the first simulations of non-headon (grazing) collisions of binary black holes in which the black hole singularities have been excised from the computational domain. Initially two equal mass black holes $m$ are separated a distance $\approx10m$ and with impact parameter $\approx2m$. Initial data are based on superposed, boosted (velocity $\approx0.5c$) solutions of single black holes in Kerr-Schild coordinates. Both rotating and non-rotating black holes are considered. The excised regions containing the singularities are specified by following the dynamics of apparent horizons. Evolutions of up to $t \approx 35m$ are obtained in which two initially separate apparent horizons are present for $t\approx3.8m$. At that time a single enveloping apparent horizon forms, indicating that the holes have merged. Apparent horizon area estimates suggest gravitational radiation of about 2.6% of the total mass. The evolutions end after a moderate amount of time because of instabilities. △ Less

Submitted 14 September, 2000; v1 submitted 13 September, 2000; originally announced September 2000.

Comments: 2 References corrected, reference to figure updated

Journal ref: Phys.Rev.Lett.85:5496-5499,2000