-
SPHINCS_BSSN: Numerical Relativity with Particles
Authors:
S. Rosswog,
P. Diener
Abstract:
In this book chapter we describe the {\em Lagrangian} numerical relativity code \sphi. This code evolves spacetimes in full General Relativity by integrating the BSSN equations on structured meshes with a simple dynamical mesh refinement strategy. The fluid is evolved by means of freely moving Lagrangian particles, that are evolved using a modern Smooth Particle Hydrodynamics (SPH) formulation. To…
▽ More
In this book chapter we describe the {\em Lagrangian} numerical relativity code \sphi. This code evolves spacetimes in full General Relativity by integrating the BSSN equations on structured meshes with a simple dynamical mesh refinement strategy. The fluid is evolved by means of freely moving Lagrangian particles, that are evolved using a modern Smooth Particle Hydrodynamics (SPH) formulation. To robustly and accurately capture shocks, our code uses artificial dissipation terms, but, similar to Finite Volume schemes, we apply a slope-limited reconstruction within the dissipative terms and we use in addition time-dependent dissipation parameters, so that dissipation is only applied where needed. The technically most complicated, but absolutely crucial part of the methodology, is the coupling between the particles and the mesh. For the mapping of the energy-momentum tensor $T_{μν}$ from the particles to the mesh, we use a sophisticated combination of "Local Regression Estimate" (LRE) method and a "multi-dimensional optimal order detection" (MOOD) approach which we describe in some detail. The mapping of the metric quantities from the grid to the particles is achieved by a quintic Hermite interpolation. Apart from giving an introduction to our numerical methods, we demonstrate the accurate working of our code by presenting a set of representative relativistic hydrodynamics tests. We begin with a relativistic shock tube test, then compare the frequencies of a fully relativistic neutron star with reference values from the literature and, finally, we present full-blown merger simulations of irrotational binary systems, one case where a central remnant survives and another where a black hole forms, and of a binary where only one of the stars is rapidly spinning.
△ Less
Submitted 21 May, 2024; v1 submitted 24 April, 2024;
originally announced April 2024.
-
Mergers of double neutron stars with one high-spin component: brighter kilonovae and fallback accretion, weaker gravitational waves
Authors:
S. Rosswog,
P. Diener,
F. Torsello,
T. M. Tauris,
N. Sarin
Abstract:
Neutron star (NS) mergers where both stars have negligible spins are commonly considered as the most likely ``standard'' case. In globular clusters, however, the majority of NSs have been spun up to millisecond (ms) periods and, based on observed systems, we estimate that a non-negligible fraction of all double NS mergers ($\sim 4\pm2\;\%$) contains one component with a spin of a (few) ms. We use…
▽ More
Neutron star (NS) mergers where both stars have negligible spins are commonly considered as the most likely ``standard'' case. In globular clusters, however, the majority of NSs have been spun up to millisecond (ms) periods and, based on observed systems, we estimate that a non-negligible fraction of all double NS mergers ($\sim 4\pm2\;\%$) contains one component with a spin of a (few) ms. We use the Lagrangian numerical relativity code \SpB to simulate mergers where one star has no spin and the other has a dimensionless spin parameter of $χ=0.5$. Such mergers exhibit several distinct signatures compared to irrotational cases. They form only one, very pronounced spiral arm and they dynamically eject an order of magnitude more mass of unshocked material at the original, very low electron fraction. One can therefore expect particularly bright, red kilonovae. Overall, the spinning case collisions are substantially less violent and they eject smaller amounts of shock-generated semi-relativistic material. Therefore, the ejecta produce a weaker blue/UV kilonova {\em precursor} signal, but -- since the total amount is larger -- brighter kilonova {\em afterglows} months after the merger. The spinning cases also have significantly more fallback accretion and thus could power late-time X-ray flares. Since the post-merger remnant loses energy and angular momentum significantly less efficiently to gravitational waves, such systems can delay a potential collapse to a black hole and are therefore candidates for merger-triggered gamma-ray bursts with longer emission time scales.
△ Less
Submitted 10 February, 2024; v1 submitted 24 October, 2023;
originally announced October 2023.
-
Thermodynamics of magnetized dense neutron-rich matter
Authors:
J. P. W. Diener
Abstract:
A neutron star is one of the possible end states of a massive star. It is compressed by gravity and stabilized by the nuclear degeneracy pressure. Despite its name, the composition of these objects is not exactly known. However, from the inferred densities, neutrons will most likely compose a significant fraction of the star's interior. While all neutron stars are expected to have a magnetic field…
▽ More
A neutron star is one of the possible end states of a massive star. It is compressed by gravity and stabilized by the nuclear degeneracy pressure. Despite its name, the composition of these objects is not exactly known. However, from the inferred densities, neutrons will most likely compose a significant fraction of the star's interior. While all neutron stars are expected to have a magnetic field, some neutron stars (''magnetars'') are much more highly magnetized than others with inferred magnetar surface magnetic field is between $10^{14}$ to $10^{15}$ gauss. While neutron stars are macroscopic objects, due to the extreme value of the stars' energy, pressure, and magnetic field the thermodynamics on the microscopic scale can be imprinted on the star's large scale behaviour. This contribution focusses on describing the thermodynamics of magnetized dense neutron matter, its equation of state and to explore conditions of a possible ferromagnetic state, contributions from the magnetized vacuum, as well as possible observational implications.
△ Less
Submitted 9 June, 2022;
originally announced June 2022.
-
Thinking outside the box: Numerical Relativity with particles
Authors:
S. Rosswog,
P. Diener,
F. Torsello
Abstract:
To date, essentially all simulation codes that solve the full set of Einstein's equations are performed in the framework of Eulerian hydrodynamics. The exception is our recently developed Numerical Relativity code SPHINCS_BSSN which solves the commonly used BSSN formulation of the Einstein equations on a structured mesh and the matter equations via Lagrangian particles. We show here, for the first…
▽ More
To date, essentially all simulation codes that solve the full set of Einstein's equations are performed in the framework of Eulerian hydrodynamics. The exception is our recently developed Numerical Relativity code SPHINCS_BSSN which solves the commonly used BSSN formulation of the Einstein equations on a structured mesh and the matter equations via Lagrangian particles. We show here, for the first time, SPHINCS_BSSN neutron star merger simulations with piecewise polytropic approximations to four nuclear matter equations of state. In this set of neutron star merger simulations we focus on perfectly symmetric binary systems that are irrotational and have 1.3 $M_\odot$ masses. We introduce some further methodological refinements (a new way of steering dissipation, an improved particle-mesh mapping) and we explore the impact of the exponent that enters in the calculation of the thermal pressure contribution. We find that it leaves a noticeable imprint on the gravitational wave amplitude (calculated via both quadrupole approximation and the $Ψ_4$-formalism) and has a noticeable impact on the amount of dynamic ejecta. Consistent with earlier findings, we only find a few times $10^{-3}$ \Msun as dynamic ejecta in the studied equal mass binary systems, with softer equations of state (which are more prone to shock formation) ejecting larger amounts of matter. In all of the cases, we see a credible high-velocity ($\sim0.5 .. 0.7c$) ejecta component of $\sim 10^{-4}$ \Msun that is launched at contact from the interface between the two neutron stars. Such a high-velocity component has been suggested to produce an early, blue precursor to the main kilonova emission and it could also potentially cause a kilonova afterglow.
△ Less
Submitted 15 June, 2022; v1 submitted 17 May, 2022;
originally announced May 2022.
-
Simulating neutron star mergers with the Lagrangian Numerical Relativity code SPHINCS_BSSN
Authors:
P. Diener,
S. Rosswog,
F. Torsello
Abstract:
We present the first neutron star merger simulations performed with the newly developed Numerical Relativity code SPHINCS_BSSN. This code evolves the spacetime on a mesh using the BSSN formulation, but matter is evolved via Lagrangian particles according to a high-accuracy version of general-relativistic Smooth Particle Hydrodynamics (SPH). Our code contains a number of new methodological elements…
▽ More
We present the first neutron star merger simulations performed with the newly developed Numerical Relativity code SPHINCS_BSSN. This code evolves the spacetime on a mesh using the BSSN formulation, but matter is evolved via Lagrangian particles according to a high-accuracy version of general-relativistic Smooth Particle Hydrodynamics (SPH). Our code contains a number of new methodological elements compared to other Numerical Relativity codes. The main focus here is on the new elements that were introduced to model neutron star mergers. These include a) a refinement (fixed in time) of the spacetime-mesh, b) corresponding changes in the particle--mesh mapping algorithm and c) a novel way to construct SPH initial data for binary systems via the recently developed "Artificial Pressure Method." This latter method makes use of the spectral initial data produced by the library LORENE, and is implemented in a new code called SPHINCS_ID. While our main focus is on introducing these new methodological elements and documenting the current status of SPHINCS_BSSN, we also show as a first application a set of neutron star merger simulations employing "soft" ($Γ=2.00$) and "stiff" ($Γ=2.75$) polytropic equations of state.
△ Less
Submitted 4 April, 2022; v1 submitted 12 March, 2022;
originally announced March 2022.
-
SPHINCS_BSSN: A general relativistic Smooth Particle Hydrodynamics code for dynamical spacetimes
Authors:
S. Rosswog,
P. Diener
Abstract:
We present a new methodology for simulating self-gravitating general-relativistic fluids. In our approach the fluid is modelled by means of Lagrangian particles in the framework of a general-relativistic (GR) Smooth Particle Hydrodynamics (SPH) formulation, while the spacetime is evolved on a mesh according to the BSSN formulation that is also frequently used in Eulerian GR-hydrodynamics. To the b…
▽ More
We present a new methodology for simulating self-gravitating general-relativistic fluids. In our approach the fluid is modelled by means of Lagrangian particles in the framework of a general-relativistic (GR) Smooth Particle Hydrodynamics (SPH) formulation, while the spacetime is evolved on a mesh according to the BSSN formulation that is also frequently used in Eulerian GR-hydrodynamics. To the best of our knowledge this is the first Lagrangian fully general relativistic hydrodynamics code (all previous SPH approaches used approximations to GR-gravity). A core ingredient of our particle-mesh approach is the coupling between the gas (represented by particles) and the spacetime (represented by a mesh) for which we have developed a set of sophisticated interpolation tools that are inspired by other particle-mesh approaches, in particular by vortex-particle methods. One advantage of splitting the methodology between matter and spacetime is that it gives us more freedom in choosing the resolution, so that -- if the spacetime is smooth enough -- we obtain good results already with a moderate number of grid cells and can focus the computational effort on the simulation of the matter. Further advantages of our approach are the ease with which ejecta can be tracked and the fact that the neutron star surface remains well-behaved and does not need any particular treatment. In the hydrodynamics part of the code we use a number of techniques that are new to SPH, such as reconstruction, slope limiting and steering dissipation by monitoring entropy conservation. We describe here in detail the employed numerical methods and demonstrate the code performance in a number of benchmark problems ranging from shock tube tests, over Cowling approximations to the fully dynamical evolution of neutron stars in self-consistently evolved spacetimes.
△ Less
Submitted 11 March, 2021; v1 submitted 27 December, 2020;
originally announced December 2020.
-
The spin-polarized ferromagnetic state of a cold Fermi gas
Authors:
J. P. W. Diener,
F. G. Scholtz
Abstract:
The spin-polarized ferromagnetic state of a cold Fermi gas is investigated for interacting and non-interacting charge-neutral and $β$-equilibrated gases. The standard minimal couplings between the magnetic field and the fermions' charges and magnetic dipole moments define the fermions' interaction with the magnetic field. Assuming a variable coupling strength between the magnetic field and the fer…
▽ More
The spin-polarized ferromagnetic state of a cold Fermi gas is investigated for interacting and non-interacting charge-neutral and $β$-equilibrated gases. The standard minimal couplings between the magnetic field and the fermions' charges and magnetic dipole moments define the fermions' interaction with the magnetic field. Assuming a variable coupling strength between the magnetic field and the fermion (baryon) dipole moments, it is shown that a ferromagnetized state can be achieved that corresponds to a lower energy spin-polarized state with a magnetic field entirely due to the gas's magnetic response. We find that, depending on the density, a very large increase in the baryon dipole moments is needed to achieve this ferromagnetized state. While the required increase seems unlikely, the induced magnetic field is of the order $\sim10^{17}$ gauss. Furthermore, while externally magnetized Fermi gases have an anisotropic pressure, the pressure of the ferromagnetized gas is completely isotropic and the thermodynamically preferred magnetized state.
△ Less
Submitted 19 November, 2020; v1 submitted 9 June, 2020;
originally announced June 2020.
-
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.
-
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.
-
Design and Testing of Kinetic Inductance Detectors Made of Titanium Nitride
Authors:
P. Diener,
H. G. Leduc,
S. J. C. Yates,
Y. J. Y. Lankwarden,
J. J. A. Baselmans
Abstract:
To use highly resistive material for Kinetic Inductance Detectors (KID), new designs have to be done, in part due to the impedance match needed between the KID chip and the whole 50 ohms readout circuit. Chips from two new hybrid designs, with an aluminum throughline coupled to titanium nitride microresonators, have been measured and compared to a TiN only chip. In the hybrid chips, parasitic temp…
▽ More
To use highly resistive material for Kinetic Inductance Detectors (KID), new designs have to be done, in part due to the impedance match needed between the KID chip and the whole 50 ohms readout circuit. Chips from two new hybrid designs, with an aluminum throughline coupled to titanium nitride microresonators, have been measured and compared to a TiN only chip. In the hybrid chips, parasitic temperature dependent box resonances are absent. The dark KID properties have been measured in a large set of resonators. A surprisingly long lifetime, up to 5.6 ms is observed in a few KIDs. For the other more reproducible devices, the mean electrical Noise Equivalent Power is 5.4 10-19 W.Hz1/2.
△ Less
Submitted 17 September, 2013;
originally announced September 2013.
-
Ferromagnetism in neutron and charge neutral beta-equilibrated nuclear matter
Authors:
J. P. W. Diener,
F. G. Scholtz
Abstract:
Ferromagnetism in infinite neutron matter as well as beta equilibrated, charge neutral, dense, and infinite nuclear matter is investigated using a model of interacting baryons and mesons. The standard minimal couplings between the magnetic field and the particle charges as well as the baryon dipole moments are included in the Lagrangian density. Minimizing the energy density with respect to the ma…
▽ More
Ferromagnetism in infinite neutron matter as well as beta equilibrated, charge neutral, dense, and infinite nuclear matter is investigated using a model of interacting baryons and mesons. The standard minimal couplings between the magnetic field and the particle charges as well as the baryon dipole moments are included in the Lagrangian density. Minimizing the energy density with respect to the magnetic field yields a self-consistent expression for the ferromagnetic field. We calculate the phase boundary at a given density by increasing the strength of the baryon dipole moments till the energy density of magnetized matter is lower than that of unmagnetized matter. We find that, depending on the density, it is crossed when the baryon dipole moments are increased by a factor of 35. It is also sensitive to the details of the nuclear matter parameterizations and crossing it induces a magnetic field of $\sim 10^{17}$ gauss.
△ Less
Submitted 11 June, 2020; v1 submitted 10 June, 2013;
originally announced June 2013.
-
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.
-
Photon noise limited radiation detection with lens-antenna coupled Microwave Kinetic Inductance Detectors
Authors:
S. J. C. Yates,
J. J. A. Baselmans,
A. Endo,
R. M. J. Janssen,
L. Ferrari,
P. Diener,
A. M. Baryshev
Abstract:
Microwave Kinetic Inductance Detectors (MKIDs) have shown great potential for sub-mm instrumentation because of the high scalability of the technology. Here we demonstrate for the first time in the sub-mm band (0.1...2 mm) a photon noise limited performance of a small antenna coupled MKID detector array and we describe the relation between photon noise and MKID intrinsic generation-recombination n…
▽ More
Microwave Kinetic Inductance Detectors (MKIDs) have shown great potential for sub-mm instrumentation because of the high scalability of the technology. Here we demonstrate for the first time in the sub-mm band (0.1...2 mm) a photon noise limited performance of a small antenna coupled MKID detector array and we describe the relation between photon noise and MKID intrinsic generation-recombination noise. Additionally we use the observed photon noise to measure the optical efficiency of detectors to be 0.8+-0.2.
△ Less
Submitted 21 July, 2011;
originally announced July 2011.
-
Number fluctuations of sparse quasiparticles in a superconductor
Authors:
P. J. de Visser,
J. J. A. Baselmans,
P. Diener,
S. J. C. Yates,
A. Endo,
T. M. Klapwijk
Abstract:
We have directly measured quasiparticle number fluctuations in a thin film superconducting Al resonator in thermal equilibrium. The spectrum of these fluctuations provides a measure of both the density and the lifetime of the quasiparticles. We observe that the quasiparticle density decreases exponentially with decreasing temperature, as theoretically predicted, but saturates below 160 mK to 25-55…
▽ More
We have directly measured quasiparticle number fluctuations in a thin film superconducting Al resonator in thermal equilibrium. The spectrum of these fluctuations provides a measure of both the density and the lifetime of the quasiparticles. We observe that the quasiparticle density decreases exponentially with decreasing temperature, as theoretically predicted, but saturates below 160 mK to 25-55 per cubic micron. We show that this saturation is consistent with the measured saturation in the quasiparticle lifetime, which also explains similar observations in qubit decoherence times.
△ Less
Submitted 25 March, 2011; v1 submitted 3 March, 2011;
originally announced March 2011.
-
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.
-
Relativistic mean-field theory applied to the study of neutron star properties
Authors:
J. P. W. Diener
Abstract:
Nuclear physics can be applied in various ways to the study of neutron stars. This thesis reports on one such application, where the relativistic mean-field approximation has been employed to calculate the equations of state of matter in the neutron star interior. In particular the equations of state of nuclear and neutron star matter of the NL3, PK1 and FSUGold parameter sets were derived. A surv…
▽ More
Nuclear physics can be applied in various ways to the study of neutron stars. This thesis reports on one such application, where the relativistic mean-field approximation has been employed to calculate the equations of state of matter in the neutron star interior. In particular the equations of state of nuclear and neutron star matter of the NL3, PK1 and FSUGold parameter sets were derived. A survey of available literature on neutron stars is presented and we use the derived equations of state to reproduce the properties of saturated nuclear matter as well as the mass-radius relationship of a static, spherical symmetric neutron star. Results are compared to published values of the properties of saturated nuclear matter and to available observational data of the mass-radius relationship of neutron stars.
△ Less
Submitted 24 June, 2010; v1 submitted 4 June, 2008;
originally announced June 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.
-
A template bank for gravitational waveforms from coalescing binary black holes: non-spinning binaries
Authors:
P. Ajith,
S. Babak,
Y. Chen,
M. Hewitson,
B. Krishnan,
A. M. Sintes,
J. T. Whelan,
B. Bruegmann,
P. Diener,
N. Dorband,
J. Gonzalez,
M. Hannam,
S. Husa,
D. Pollney,
L. Rezzolla,
L. Santamaria,
U. Sperhake,
J. Thornburg
Abstract:
Gravitational waveforms from the inspiral and ring-down stages of the binary black hole coalescences can be modelled accurately by approximation/perturbation techniques in general relativity. Recent progress in numerical relativity has enabled us to model also the non-perturbative merger phase of the binary black-hole coalescence problem. This enables us to \emph{coherently} search for all three…
▽ More
Gravitational waveforms from the inspiral and ring-down stages of the binary black hole coalescences can be modelled accurately by approximation/perturbation techniques in general relativity. Recent progress in numerical relativity has enabled us to model also the non-perturbative merger phase of the binary black-hole coalescence problem. This enables us to \emph{coherently} search for all three stages of the coalescence of non-spinning binary black holes using a single template bank. Taking our motivation from these results, we propose a family of template waveforms which can model the inspiral, merger, and ring-down stages of the coalescence of non-spinning binary black holes that follow quasi-circular inspiral. This two-dimensional template family is explicitly parametrized by the physical parameters of the binary. We show that the template family is not only \emph{effectual} in detecting the signals from black hole coalescences, but also \emph{faithful} in estimating the parameters of the binary. We compare the sensitivity of a search (in the context of different ground-based interferometers) using all three stages of the black hole coalescence with other template-based searches which look for individual stages separately. We find that the proposed search is significantly more sensitive than other template-based searches for a substantial mass-range, potentially bringing about remarkable improvement in the event-rate of ground-based interferometers. As part of this work, we also prescribe a general procedure to construct interpolated template banks using non-spinning black hole waveforms produced by numerical relativity.
△ Less
Submitted 11 May, 2009; v1 submitted 11 October, 2007;
originally announced October 2007.
-
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.
-
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.
-
Phenomenological template family for black-hole coalescence waveforms
Authors:
P. Ajith,
S. Babak,
Y. Chen,
M. Hewitson,
B. Krishnan,
J. T. Whelan,
B. Bruegmann,
P. Diener,
J. Gonzalez,
M. Hannam,
S. Husa,
M. Koppitz,
D. Pollney,
L. Rezzolla,
L. Santamaria,
A. M. Sintes,
U. Sperhake,
J. Thornburg
Abstract:
Recent progress in numerical relativity has enabled us to model the non-perturbative merger phase of the binary black-hole coalescence problem. Based on these results, we propose a phenomenological family of waveforms which can model the inspiral, merger, and ring-down stages of black hole coalescence. We also construct a template bank using this family of waveforms and discuss its implementatio…
▽ More
Recent progress in numerical relativity has enabled us to model the non-perturbative merger phase of the binary black-hole coalescence problem. Based on these results, we propose a phenomenological family of waveforms which can model the inspiral, merger, and ring-down stages of black hole coalescence. We also construct a template bank using this family of waveforms and discuss its implementation in the search for signatures of gravitational waves produced by black-hole coalescences in the data of ground-based interferometers. This template bank might enable us to extend the present inspiral searches to higher-mass binary black-hole systems, i.e., systems with total mass greater than about 80 solar masses, thereby increasing the reach of the current generation of ground-based detectors.
△ Less
Submitted 11 July, 2007; v1 submitted 30 April, 2007;
originally announced April 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.