Showing 1–50 of 57 results for author: Palenzuela, C

Large Eddy Simulations of Magnetized Mergers of Black Holes and Neutron Stars

Authors: Manuel R. Izquierdo, Miguel Bezares, Steven Liebling, Carlos Palenzuela

Abstract: The LIGO-Virgo-Kagra collaboration has observed gravitational waves consistent with the mergers of a black hole and a neutron star, namely GW200105 and GW200115, providing evidence for such cataclysmic events. Although no electromagnetic counterpart was reported for either of these two events, under certain conditions black hole--neutron star mergers are expected to form a significant accretion di… ▽ More The LIGO-Virgo-Kagra collaboration has observed gravitational waves consistent with the mergers of a black hole and a neutron star, namely GW200105 and GW200115, providing evidence for such cataclysmic events. Although no electromagnetic counterpart was reported for either of these two events, under certain conditions black hole--neutron star mergers are expected to form a significant accretion disk and to produce both a short gamma ray burst and a kilonova, much as observed in the binary neutron star merger GW170817. Here, we extend our publicly available code \mhduet to study numerically the merger of a magnetized neutron star with a black hole. \mhduet employs Large Eddy Simulation techniques to help capture the magnetic field amplification resulting from turbulence and other sub-grid scale dynamics in the post-merger stage. In particular, we simulate a merger with parameters favorable to producing an accretion disk, focusing on the formation and dynamics of the turbulent disk and the resulting magnetic field amplification. Following the tidal disruption and during the formation of the accretion disk, the magnetic field undergoes significant amplification driven by the Kelvin-Helmholtz instability, reaching strengths of more than $10^{14}\,\rm{G}$ from a realistic initial strength of $10^{11}\,\rm{G}$ in short timescales of approximately $20\,\rm{ms}$. At later times in our simulations the magnetic field growth occurs at larger scales and is dominated mainly by magnetic winding. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 16 pages, 11 figures. Comments are welcome!

The Guided Moments formalism: a new efficient full-neutrino treatment for astrophysical simulations

Authors: Manuel R. Izquierdo, J. Fernando Abalos, Carlos Palenzuela

Abstract: We present the new Guided Moments ($\texttt{GM}$) formalism for neutrino modeling in astrophysical scenarios like core-collapse supernovae and neutron star mergers. The truncated moments approximation ($\texttt{M1}$) and Monte-Carlo ($\texttt{MC}$) schemes have been proven to be robust and accurate in solving the Boltzmann's equation for neutrino transport. However, it is well-known that each meth… ▽ More We present the new Guided Moments ($\texttt{GM}$) formalism for neutrino modeling in astrophysical scenarios like core-collapse supernovae and neutron star mergers. The truncated moments approximation ($\texttt{M1}$) and Monte-Carlo ($\texttt{MC}$) schemes have been proven to be robust and accurate in solving the Boltzmann's equation for neutrino transport. However, it is well-known that each method exhibits specific strengths and weaknesses in various physical scenarios. The $\texttt{GM}$ formalism effectively solves these problems, providing a comprehensive scheme capable of accurately capturing the optically thick limit through the exact $\texttt{M1}$ closure and the optically thin limit through a $\texttt{MC}$ based approach. In addition, the $\texttt{GM}$ method also approximates the neutrino distribution function with a reasonable computational cost, which is crucial for the correct estimation of the different neutrino-fluid interactions. Our work provides a comprehensive discussion of the formulation and application of the $\texttt{GM}$ method, concluding with a thorough comparison across several test problems involving the three schemes ($\texttt{M1}$, $\texttt{MC}$, $\texttt{GM}$) under consideration. △ Less

Submitted 31 January, 2024; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: 17+3 pages, 14 figures. Accepted in PRD

Journal ref: Phys. Rev. D 109, 043044 (2024)

The role of turbulence and winding in the development of large-scale, strong magnetic fields in long-lived remnants of binary neutron star mergers

Authors: Ricard Aguilera-Miret, Carlos Palenzuela, Federico Carrasco, Daniele Viganò

Abstract: We perform a long and accurate Large-Eddy Simulation of a binary neutron star merger, following the newly formed remnant up to 110 milliseconds. The combination of high-order schemes, high-resolution and the gradient subgrid-scale model allow us to have among the highest effective resolutions ever achieved. Our results show that, although the magnetic fields are strongly amplified by the Kelvin-He… ▽ More We perform a long and accurate Large-Eddy Simulation of a binary neutron star merger, following the newly formed remnant up to 110 milliseconds. The combination of high-order schemes, high-resolution and the gradient subgrid-scale model allow us to have among the highest effective resolutions ever achieved. Our results show that, although the magnetic fields are strongly amplified by the Kelvin-Helmholtz instability, they are coherent only over very short spatial scales until t \gtrsim 30 ms. Around that time, magnetic winding becomes more efficient leading to a linear growth of the toroidal component and slowly ordering the field to more axisymmetric, large scales. The poloidal component only starts to grow at small scales at much later times t \gtrsim 90 ms, in a way compatible with the magneto-rotational instability. No strong large-scale poloidal field or jet is produced in the timescales spanned by our simulation, although there is an helicoidal structure gradually developing at late times. We highlight that soon after the merger the topology is always strongly dominated by toroidal structures, with a complex distribution in the meridional plane and highly turbulent perturbations. Thus, starting with strong purely dipolar fields before the merger is largely inconsistent with the outcomes of a realistic evolution. Finally, we confirm the universality of the evolved topology, even when starting with very different magnetic fields confined to the outermost layers of each neutron star. △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 18 pages, 14 figures. arXiv admin note: text overlap with arXiv:2112.08413

Magnetic winding and turbulence in ultra-hot Jupiters

Authors: Clàudia Soriano-Guerrero, Daniele Viganò, Rosalba Perna, Taner Akgün, Carlos Palenzuela

Abstract: While magnetism in exoplanets remains largely unknown, Hot Jupiters have been considered as natural candidates to harbour intense magnetic fields, both due to their large masses and their high energy budgets coming from irradiation as a consequence of their vicinity to their host stars. In this work we perform MHD simulations of a narrow day-side atmospheric column of ultra-hot Jupiters, suitable… ▽ More While magnetism in exoplanets remains largely unknown, Hot Jupiters have been considered as natural candidates to harbour intense magnetic fields, both due to their large masses and their high energy budgets coming from irradiation as a consequence of their vicinity to their host stars. In this work we perform MHD simulations of a narrow day-side atmospheric column of ultra-hot Jupiters, suitable for very high local temperatures (T > 3000 K). Since the conductivity in this regime is very high, the dominant effect is winding due to the intense zonal winds. By including a forcing that mimics the wind profiles obtained in global circulation models, the shear layer induces a strong toroidal magnetic field (locally reaching hundreds of gauss), supported by meridional currents. Such fields and the sustaining currents don$'$t depend on the internally generated field, but are all confined in the thin (less than a scale-height) shear layer around 1 bar. Additionally, we add random perturbations that induce turbulent motions, which lead to further (but much smaller) magnetic field generation to a broader range of depths. These results allow an evaluation of the currents induced by the atmospheric dynamo. Although here we use ideal MHD and the only resistivity comes from the numerical scheme, we estimate a-posteriori the amount of Ohmic heat deposited in the outer layers, which could be employed in evolutionary models for Hot Jupiters' inflated radii. △ Less

Submitted 27 July, 2023; v1 submitted 14 April, 2023; originally announced April 2023.

Comments: Accepted for publication in MNRAS

Global high-order numerical schemes for the time evolution of the general relativistic radiation magneto-hydrodynamics equations

Authors: Manuel R. Izquierdo, Lorenzo Pareschi, Borja Miñano, Joan Massó, Carlos Palenzuela

Abstract: Modeling correctly the transport of neutrinos is crucial in some astrophysical scenarios such as core-collapse supernovae and binary neutron star mergers. In this paper, we focus on the truncated-moment formalism, considering only the first two moments (M1 scheme) within the grey approximation, which reduces Boltzmann seven-dimensional equation to a system of $3+1$ equations closely resembling the… ▽ More Modeling correctly the transport of neutrinos is crucial in some astrophysical scenarios such as core-collapse supernovae and binary neutron star mergers. In this paper, we focus on the truncated-moment formalism, considering only the first two moments (M1 scheme) within the grey approximation, which reduces Boltzmann seven-dimensional equation to a system of $3+1$ equations closely resembling the hydrodynamic ones. Solving the M1 scheme is still mathematically challenging, since it is necessary to model the radiation-matter interaction in regimes where the evolution equations become stiff and behave as an advection-diffusion problem. Here, we present different global, high-order time integration schemes based on Implicit-Explicit Runge-Kutta (IMEX) methods designed to overcome the time-step restriction caused by such behavior while allowing us to use the explicit RK commonly employed for the MHD and Einstein equations. Finally, we analyze their performance in several numerical tests. △ Less

Submitted 29 May, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: 18 + 6. Updated manuscript matching published version + additional appendix "Comparing the convergence order of IMEX and semi-implicit schemes"

Relativistic force-free models of the thermal X-ray emission in millisecond pulsars observed by NICER

Authors: Federico Carrasco, Joaquin Pelle, Oscar Reula, Daniele Viganò, Carlos Palenzuela

Abstract: Several important properties of rotation-powered millisecond pulsars (MSPs), such as their mass-radius ratio, equation of state and magnetic field topology, can be inferred from precise observations and modelling of their X-ray light curves. In the present study, we model the thermal X-ray signals originated in MSPs, all the way from numerically solving the surrounding magnetospheres up to the ray… ▽ More Several important properties of rotation-powered millisecond pulsars (MSPs), such as their mass-radius ratio, equation of state and magnetic field topology, can be inferred from precise observations and modelling of their X-ray light curves. In the present study, we model the thermal X-ray signals originated in MSPs, all the way from numerically solving the surrounding magnetospheres up to the ray tracing of the emitted photons and the final computation of their light curves and spectra. Our modelled X-ray signals are then compared against a set of very accurate NICER observations of four target pulsars: PSR J0437-4715, PSR J1231-1411, PSR J2124-3358 and PSR J0030+0451. We find very good simultaneous fits for the light curve and spectral distribution in all these pulsars. The magnetosphere is solved by performing general relativistic force-free simulations of a rotating neutron star (NS) endowed with a simple centered dipolar magnetic field, for many different stellar compactness and pulsar misalignments. From these solutions, we derive an emissivity map over the surface of the star, which is based on the electric currents in the magnetosphere. In particular, the emission regions (ERs) are determined in this model by spacelike four-currents that reach the NS. We show that this assumption, together with the inclusion of the gravitational curvature on the force-free simulations, lead to non-standard ERs facing the closed-zone of the pulsar, in addition to other ERs within the polar caps. The combined X-ray signals from these two kinds of ERs (both antipodal) allow to approximate the non-trivial interpulses found in all the target MSPs light curves. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: 16 pages, 9 figures

Large Eddy Simulations of Magnetized Mergers of Neutron Stars with Neutrinos

Authors: C. Palenzuela, S. L. Liebling, B. Miñano

Abstract: Neutron star mergers are very violent events involving extreme physical processes: dynamic, strong-field gravity, large magnetic field, very hot, dense matter, and the copious production of neutrinos. Accurate modeling of such a system and its associated multi-messenger signals, such as gravitational waves, short gamma ray burst, and kilonova, requires the inclusion of all these processes, and is… ▽ More Neutron star mergers are very violent events involving extreme physical processes: dynamic, strong-field gravity, large magnetic field, very hot, dense matter, and the copious production of neutrinos. Accurate modeling of such a system and its associated multi-messenger signals, such as gravitational waves, short gamma ray burst, and kilonova, requires the inclusion of all these processes, and is increasingly important in light of advancements in multi-messenger astronomy generally, and in gravitational wave astronomy in particular (such as the development of third-generation detectors). Several general relativistic codes have been incorporating some of these elements with different levels of realism. Here, we extend our code MHDuet, which can perform large eddy simulations of magnetohydrodynamics to help capture the magnetic field amplification during the merger, and to allow for realistic equations of state and neutrino cooling via a leakage scheme. We perform several tests involving isolated and binary neutron stars demonstrating the accuracy of the code. △ Less

Submitted 6 May, 2022; v1 submitted 6 April, 2022; originally announced April 2022.

Comments: 20 pages, 11 figures. Modified to match the accepted version of the paper

Turbulent magnetic field amplification in binary neutron star mergers

Authors: C. Palenzuela, R. Aguilera-Miret, F. Carrasco, R. Ciolfi, J. V. Kalinani, W. Kastaun, B. Miñano, D. Viganò

Abstract: Magnetic fields are expected to play a key role in the dynamics and the ejection mechanisms that accompany the merger of two neutron stars. General relativistic magnetohydrodynamic (MHD) simulations offer a unique opportunity to unravel the details of the ongoing physical processes. Nevertheless, current numerical studies are severely limited by the fact that any affordable resolution remains insu… ▽ More Magnetic fields are expected to play a key role in the dynamics and the ejection mechanisms that accompany the merger of two neutron stars. General relativistic magnetohydrodynamic (MHD) simulations offer a unique opportunity to unravel the details of the ongoing physical processes. Nevertheless, current numerical studies are severely limited by the fact that any affordable resolution remains insufficient to fully capture the small-scale dynamo, initially triggered by the Kelvin-Helmholtz instability, and later sourced by several MHD processes involving differential rotation. Here, we alleviate this limitation by using explicit large-eddy simulations, a technique where the unresolved dynamics occurring at the sub-grid scales (SGS) is modeled by extra terms, which are functions of the resolved fields and their derivatives. The combination of high-order numerical schemes, high resolutions, and the gradient SGS model allow us to capture the small-scale dynamos produced during the binary neutron star mergers. Here we follow the first 50 milliseconds after the merger and, for the first time, we find numerical convergence on the magnetic field amplification, in terms of integrated energy and spectral distribution over spatial scales. We also find that the average intensity of the magnetic field in the remnant saturates at $\sim 10^{16}$~G around $5$~ms after the merger. After $20-30$~ms, both toroidal and poloidal magnetic field components grow continuously, fed by the winding mechanism that provides a slow inverse cascade. We find no clear hints for magneto-rotational instabilities, and no significant impact of the magnetic field on the redistribution of angular momentum in the remnant in our simulations, probably due to the very turbulent and dynamical topology of the magnetic field at all stages, with small-scale components largely dominating over the large-scale ones. △ Less

Submitted 3 July, 2022; v1 submitted 15 December, 2021; originally announced December 2021.

Comments: 22 pages, 17 figures, version matching the published article

Universality of the turbulent magnetic field in hypermassive neutron stars produced by binary mergers

Authors: Ricard Aguilera-Miret, Daniele Viganò, Carlos Palenzuela

Abstract: The detection of a binary neutron star merger in 2017 through both gravitational waves and electromagnetic emission opened a new era of multimessenger astronomy. The understanding of the magnetic field amplification triggered by the Kelvin-Helmholtz instability during the merger is still a numerically unresolved problem because of the relevant small scales involved. One of the uncertainties comes… ▽ More The detection of a binary neutron star merger in 2017 through both gravitational waves and electromagnetic emission opened a new era of multimessenger astronomy. The understanding of the magnetic field amplification triggered by the Kelvin-Helmholtz instability during the merger is still a numerically unresolved problem because of the relevant small scales involved. One of the uncertainties comes from the simplifications usually assumed in the initial magnetic topology of merging neutron stars. We perform high-resolution, convergent large-eddy simulations of binary neutron star mergers, following the newly formed remnant for up to $30$ milliseconds. Here we specifically focus on the comparison between simulations with different initial magnetic configurations, going beyond the widespread-used aligned dipole confined within each star. The results obtained show that the initial topology is quickly forgotten, in a timescale of few miliseconds after the merger. Moreover, at the end of the simulations, the average intensity ($B\sim 10^{16}$ G) and the spectral distribution of magnetic energy over spatial scales barely depend on the initial configuration. This is expected due to the small-scale efficient dynamo involved, and thus it holds as long as: (i) the initial large-scale magnetic field is not unrealistically high (as often imposed in mergers studies); (ii) the turbulent instability is numerically (at least partially) resolved, so that the amplified magnetic energy is distributed across a wide range of scales and becomes orders-of-magnitude larger than the initial one. △ Less

Submitted 28 March, 2022; v1 submitted 15 December, 2021; originally announced December 2021.

Comments: 8 pages, 5 figures

Kinetic screening in nonlinear stellar oscillations and gravitational collapse

Authors: Miguel Bezares, Lotte ter Haar, Marco Crisostomi, Enrico Barausse, Carlos Palenzuela

Abstract: We consider k-essence, a scalar-tensor theory with first-order derivative self-interactions that can screen local scales from scalar fifth forces, while allowing for sizeable deviations from General Relativity on cosmological scales. We construct fully nonlinear static stellar solutions that show the presence of this screening mechanism, and we use them as initial data for simulations of stellar o… ▽ More We consider k-essence, a scalar-tensor theory with first-order derivative self-interactions that can screen local scales from scalar fifth forces, while allowing for sizeable deviations from General Relativity on cosmological scales. We construct fully nonlinear static stellar solutions that show the presence of this screening mechanism, and we use them as initial data for simulations of stellar oscillations and gravitational collapse in spherical symmetry. We find that for k-essence theories of relevance for cosmology, the screening mechanism works in the case of stellar oscillation and suppresses the monopole scalar emission to undetectable levels. In collapsing stars, we find that the Cauchy problem, although locally well posed, can lead to diverging characteristic speeds for the scalar field. By introducing a ''fixing equation'' in the spirit of J. Cayuso, N. Ortiz, and L. Lehner [Phys. Rev. D 96, 084043 (2017)], inspired in turn by dissipative relativistic hydrodynamics, we manage to evolve collapsing neutron stars past the divergence of the characteristic speeds. We show that, in these systems, the screening mechanism is less efficient than for oscillating and static stars, because the collapsing star must shed away all of its scalar hair before forming a black hole. For k-essence theories of relevance for cosmology, the characteristic frequency of the resulting scalar monopole signal is too low for terrestrial detectors, but we conjecture that space-borne interferometers such as LISA might detect it if a supernova explodes in the Galaxy. △ Less

Submitted 12 August, 2021; v1 submitted 28 May, 2021; originally announced May 2021.

Comments: 17 pages, 11 figures. Minor changes, matches version published in PRD

Journal ref: Phys. Rev. D 104, 044022 (2021)

Multimessenger signals from black hole-neutron star mergers without significant tidal disruption

Authors: William E. East, Luis Lehner, Steven L. Liebling, Carlos Palenzuela

Abstract: We study the multimessenger signals from the merger of a black hole with a magnetized neutron star using resistive magnetohydrodynamics simulations coupled to full general relativity. We focus on a case with a 5:1 mass ratio, where only a small amount of the neutron star matter remains post-merger, but we nevertheless find that significant electromagnetic radiation can be powered by the interactio… ▽ More We study the multimessenger signals from the merger of a black hole with a magnetized neutron star using resistive magnetohydrodynamics simulations coupled to full general relativity. We focus on a case with a 5:1 mass ratio, where only a small amount of the neutron star matter remains post-merger, but we nevertheless find that significant electromagnetic radiation can be powered by the interaction of the neutron star's magnetosphere with the black hole. In the lead-up to merger, strong twisting of magnetic field lines from the inspiral leads to plasmoid emission and results in a luminosity in excess of that expected from unipolar induction. We find that the strongest emission occurs shortly after merger during a transitory period in which magnetic loops form and escape the central region. The remaining magnetic field collimates around the spin axis of the remnant black hole before dissipating, an indication that, in more favorable scenarios (higher black hole spin/lower mass ratio) with larger accretion disks, a jet would form. △ Less

Submitted 23 April, 2021; v1 submitted 28 January, 2021; originally announced January 2021.

Comments: 9 pages, 5 figures; revised to match ApJL version

Journal ref: ApJL 912 L18 (2021)

Turbulent magnetic-field amplification in the first 10 milliseconds after a binary neutron star merger: comparing high-resolution and large eddy simulations

Authors: Ricard Aguilera-Miret, Daniele Viganò, Federico Carrasco, Borja Miñano, Carlos Palenzuela

Abstract: The detection of binary neutron star mergers represents one of the most important and complex astrophysical discoveries of the recent years. One of the unclear aspects of the problem is the turbulent magnetic field amplification, initially triggered by the Kelvin-Helmholtz instability at much smaller scales than any reachable numerical resolution nowadays. Here we present numerical simulations of… ▽ More The detection of binary neutron star mergers represents one of the most important and complex astrophysical discoveries of the recent years. One of the unclear aspects of the problem is the turbulent magnetic field amplification, initially triggered by the Kelvin-Helmholtz instability at much smaller scales than any reachable numerical resolution nowadays. Here we present numerical simulations of the first ten milliseconds of a binary neutron star merger. First, we confirm in detail how the simulated amplification depends on the numerical resolution and is distributed on a broad range of scales, as expected from turbulent MHD theory. We find that an initial large-scale magnetic field of $10^{11}\,$G inside each star is amplified in the remnant to root-mean-square values above $10^{16}\,$G within the first $5$ milliseconds for our highest-resolution run. Then, we run large eddy simulations, exploring the performance of the subgrid-scale gradient model, already tested successfully in previous turbulent box simulations. We show that the addition of this model is especially important in the induction equation, since it leads to an amplification of the magnetic field comparable to a higher-resolution run, but with a greatly reduced computational cost. In the first 10 milliseconds, there is no clear hint for an ordered, large-scale magnetic field, which should indeed occur in longer timescales through magnetic winding and the magneto-rotational instability. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 14 pages, 7 figures

Journal ref: Phys. Rev. D 102, 103006 (2020)

Dynamics of Screening in Modified Gravity

Authors: Lotte ter Haar, Miguel Bezares, Marco Crisostomi, Enrico Barausse, Carlos Palenzuela

Abstract: Gravitational theories differing from General Relativity may explain the accelerated expansion of the Universe without a cosmological constant. However, to pass local gravitational tests, a "screening mechanism" is needed to suppress, on small scales, the fifth force driving the cosmological acceleration. We consider the simplest of these theories, i.e. a scalar-tensor theory with first-order deri… ▽ More Gravitational theories differing from General Relativity may explain the accelerated expansion of the Universe without a cosmological constant. However, to pass local gravitational tests, a "screening mechanism" is needed to suppress, on small scales, the fifth force driving the cosmological acceleration. We consider the simplest of these theories, i.e. a scalar-tensor theory with first-order derivative self-interactions, and study isolated (static and spherically symmetric) non-relativistic and relativistic stars. We produce screened solutions and use them as initial data for non-linear numerical evolutions in spherical symmetry. We find that these solutions are stable under large initial perturbations, as long as they do not cause gravitational collapse. When gravitational collapse is triggered, the characteristic speeds of the scalar evolution equation diverge, even before apparent black-hole or sound horizons form. This casts doubts on whether the dynamical evolution of screened stars may be predicted in these effective field theories. △ Less

Submitted 1 March, 2021; v1 submitted 7 September, 2020; originally announced September 2020.

Comments: 6 pages, 3 figures; PRL version

Journal ref: Phys. Rev. Lett. 126, 091102 (2021)

Introduction to Numerical Relativity

Authors: Carlos Palenzuela

Abstract: Numerical Relativity is a multidisciplinary field including relativity, magneto-hydrodynamics, astrophysics and computational methods, among others, with the aim of solving numerically highly-dynamical, strong-gravity scenarios where no other approximations are available. Here we describe some of the foundations of the field, starting from the covariant Einstein equations and how to write them as… ▽ More Numerical Relativity is a multidisciplinary field including relativity, magneto-hydrodynamics, astrophysics and computational methods, among others, with the aim of solving numerically highly-dynamical, strong-gravity scenarios where no other approximations are available. Here we describe some of the foundations of the field, starting from the covariant Einstein equations and how to write them as a well-posed system of evolution equations, discussing the different formalisms, coordinate conditions and numerical methods commonly employed nowadays for the modeling of gravitational wave sources. △ Less

Submitted 29 August, 2020; originally announced August 2020.

Comments: Accepted by Frontiers Astronomy and Space Sciences, invited review for the Research Topic "Gravitational Waves: A New Window to the Universe"

K-dynamics: well-posed 1+1 evolutions in K-essence

Authors: Miguel Bezares, Marco Crisostomi, Carlos Palenzuela, Enrico Barausse

Abstract: We study the vacuum Cauchy problem for K-essence, i.e. cosmologically relevant scalar-tensor theories that involve first-order derivative self-interactions, and which pass all existing gravitational wave bounds. We restrict to spherical symmetry and show that there exists a large class of theories for which no breakdown of the Cauchy problem occurs outside apparent black hole horizons, even in the… ▽ More We study the vacuum Cauchy problem for K-essence, i.e. cosmologically relevant scalar-tensor theories that involve first-order derivative self-interactions, and which pass all existing gravitational wave bounds. We restrict to spherical symmetry and show that there exists a large class of theories for which no breakdown of the Cauchy problem occurs outside apparent black hole horizons, even in the presence of scalar shocks/caustics, except for a small set of initial data sufficiently close to critical black hole collapse. We characterise these problematic initial data, and show that they lead to large or even diverging (coordinate) speeds for the characteristic curves. We discuss the physical relevance of this problem and propose ways to overcome it. △ Less

Submitted 22 March, 2021; v1 submitted 17 August, 2020; originally announced August 2020.

Comments: 13 pages, 3 figures, references added, published version

Journal ref: JCAP03(2021)072

GRMHD large eddy simulations with gradient subgrid-scale model

Authors: Daniele Viganò, Ricard Aguilera-Miret, Federico Carrasco, Borja Miñano, Carlos Palenzuela

Abstract: The detection of binary neutron star mergers represents one of the most important astrophysical discoveries of the recent years. Due to the extreme matter and gravity conditions and the rich dynamics developed, it becomes a tremendous challenge to accurately simulate numerically all the scales present during the collision. Here we present how to study such systems by using large eddy simulations w… ▽ More The detection of binary neutron star mergers represents one of the most important astrophysical discoveries of the recent years. Due to the extreme matter and gravity conditions and the rich dynamics developed, it becomes a tremendous challenge to accurately simulate numerically all the scales present during the collision. Here we present how to study such systems by using large eddy simulations with a self-consistent subgrid-scale gradient model, that we generalized to the special relativistic case in a previous work and now extend to the general relativistic case. Adapted from nonrelativistic scenarios, the so-called gradient model allows to capture part of the effects of the hidden dynamics on the resolved scales, by means of a physically-agnostic, mathematically-based Taylor expansion of the nonlinear terms in the conservative evolution equations' fluxes. We assess the validity of this approach in bounding-box simulations of the magnetic Kelvin-Helmholtz instability. Several resolutions and a broad range of scenarios are considered in order to carefully test the performance of the model under three crucial aspects: (i) highly curved backgrounds, (ii) jumps on the fluid density profiles and (iii) strong shocks. The results suggest our extension of the gradient subgrid-scale model to general relativistic magnetohydrodynamics is a promising approach for studying binary neutron stars mergers, and potentially to other relevant astrophysical scenarios. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 15 pages, 7 figures Submitted to PRD

Journal ref: Phys. Rev. D 101, 123019 (2020)

Gradient sub-grid-scale model for relativistic MHD Large Eddy Simulations

Authors: Federico Carrasco, Daniele Viganò, Carlos Palenzuela

Abstract: MHD turbulence is likely to play an important role in several astrophysical scenarios where the magnetic Reynolds is very large. Numerically, these cases can be studied efficiently by means of Large Eddy Simulations, in which the computational resources are used to evolve the system only up to a finite grid size. The resolution is not fine enough to capture all the relevant small-scale physics at… ▽ More MHD turbulence is likely to play an important role in several astrophysical scenarios where the magnetic Reynolds is very large. Numerically, these cases can be studied efficiently by means of Large Eddy Simulations, in which the computational resources are used to evolve the system only up to a finite grid size. The resolution is not fine enough to capture all the relevant small-scale physics at play, which is instead effectively modeled by a set of additional terms in the evolution equations, dubbed as sub-grid-scale model. Here we extend such approach, commonly used in non-relativistic/non-magnetic/incompressible fluid dynamics, applying the so-called gradient model to a general set of balance-law equations, that includes the relevant case in which a non-trivial inversion of conserved to primitive fields is needed. In particular, we focus on the relativistic compressible ideal MHD scenario, providing for the first time (and for any equation of state) all the additional sub-grid-scale terms. As an application, we consider box simulations of the relativistic Kelvin-Helmholtz instability, which is also the first mechanism responsible for the magnetic field amplification in binary neutron star mergers and cannot yet be fully captured by the finest-grid and longest simulations available. The performance of our model is numerically assessed by comparing it to the residuals arising from the filtering of high-resolution simulations. We find that the model can fit very well those residuals from resolutions a few times higher. Although the application shown here explicitly considers the Minkowski metric, it can be directly extended to general relativity, thus settling the basis to implement the gradient sub-grid model in a GRMHD binary merger. Our results suggest that this approach will be potentially able to unveil much better the small-scale dynamics achievable in full GRMHD simulations. △ Less

Submitted 23 March, 2020; v1 submitted 4 August, 2019; originally announced August 2019.

Comments: 15 pages, 6 figures. Final version. arXiv admin note: text overlap with arXiv:1904.04099

Journal ref: Phys. Rev. D 101, 063003 (2020)

Extension of the sub-grid-scale gradient model for compressible magnetohydrodynamics turbulent instabilities

Authors: Daniele Viganò, Ricard Aguilera-Miret, Carlos Palenzuela

Abstract: Performing accurate large eddy simulations in compressible, turbulent magnetohydrodynamics is more challenging than in non-magnetized fluids due to the complex interplay between kinetic, magnetic and internal energy at different scales. Here we extend the sub-grid-scale gradient model, so far used in the momentum and induction equations, to account also for the unresolved scales in the energy evol… ▽ More Performing accurate large eddy simulations in compressible, turbulent magnetohydrodynamics is more challenging than in non-magnetized fluids due to the complex interplay between kinetic, magnetic and internal energy at different scales. Here we extend the sub-grid-scale gradient model, so far used in the momentum and induction equations, to account also for the unresolved scales in the energy evolution equation of a compressible ideal MHD fluid with a generic equation of state. We assess the model by considering box simulations of the turbulence triggered across a shear layer by the Kelvin-Helmholtz instability, testing cases where the small-scale dynamics cannot be fully captured by the resolution considered, such that the efficiency of the simulated dynamo effect depends on the resolution employed. This lack of numerical convergence is actually a currently common issue in several astrophysical problems, where the integral and fastest-growing-instability scales are too far apart to be fully covered numerically. We perform a-priori and a-posteriori tests of the extended gradient model. In the former, we find that, for many different initial conditions and resolutions, the gradient model outperforms other commonly used models in terms of correlation with the residuals coming from the filtering of a high-resolution run. In the second test, we show how a low-resolution run with the gradient model is able to quantitatively reproduce the evolution of the magnetic energy (the integrated value and the spectral distribution) coming from higher-resolution runs. This extension is the first step towards the implementation in relativistic magnetohydrodynamics. △ Less

Submitted 12 September, 2019; v1 submitted 8 April, 2019; originally announced April 2019.

Comments: 25 pages, 23 figures, 2 tables. Submitted as a first revision

Triggering magnetar outbursts in 3D force-free simulations

Authors: Federico Carrasco, Daniele Viganò, Carlos Palenzuela, Jose A. Pons

Abstract: In this letter, we present the first 3D force-free general relativity simulations of the magnetosphere dynamics related to the magnetar outburst/flare phenomenology. Starting from an initial dipole configuration, we adiabatically increase the helicity by twisting the footprints of a spot on the stellar surface and follow the succession of quasi-equilibrium states until a critical twist is reached.… ▽ More In this letter, we present the first 3D force-free general relativity simulations of the magnetosphere dynamics related to the magnetar outburst/flare phenomenology. Starting from an initial dipole configuration, we adiabatically increase the helicity by twisting the footprints of a spot on the stellar surface and follow the succession of quasi-equilibrium states until a critical twist is reached. Twisting beyond that point triggers instabilities that results in the rapid expansion of magnetic field lines, followed by reconnection, as observed in previous axi-symmetric simulations. If the injection of magnetic helicity goes on, the process is recurrent, periodically releasing a similar amount of energy, of the order of a few % of the total magnetic energy. From our current distribution, we estimate the local temperature assuming that dissipation occurs mainly in the highly resistive outermost layer of the neutron star. We find that the temperature smoothly increases with injected twist, being larger for spots located in the tropical regions than in polar regions, and rather independent of their sizes. After the injection of helicity ceases, the magnetosphere relaxes to a new stable state, in which the persistent currents maintain the footprints area slightly hotter than before the onset of the instability. △ Less

Submitted 25 January, 2019; originally announced January 2019.

Comments: 6 pages, 5 figures

A Simflowny-based high-performance 3D code for the generalized induction equation

Authors: Daniele Viganò, David Martínez-Gómez, José A. Pons, Carlos Palenzuela, Federico Carrasco, Borja Miñano, Antoni Arbona, Carles Bona, Joan Massó

Abstract: In the interior of neutron stars, the induction equation regulates the long-term evolution of the magnetic fields by means of resistivity, Hall dynamics and ambipolar diffusion. Despite the apparent simplicity and compactness of the equation, the dynamics it describes is not trivial and its understanding relies on accurate numerical simulations. While a few works in 2D have reached a mature stage… ▽ More In the interior of neutron stars, the induction equation regulates the long-term evolution of the magnetic fields by means of resistivity, Hall dynamics and ambipolar diffusion. Despite the apparent simplicity and compactness of the equation, the dynamics it describes is not trivial and its understanding relies on accurate numerical simulations. While a few works in 2D have reached a mature stage and a consensus on the general dynamics at least for some simple initial data, only few attempts have been performed in 3D, due to the computational costs and the need for a proper numerical treatment of the intrinsic non-linearity of the equation. Here, we carefully analyze the general induction equation, studying its characteristic structure, and we present a new Cartesian 3D code, generated by the user-friendly, publicly available {\em Simflowny} platform. The code uses high-order numerical schemes for the time and spatial discretization, and relies on the highly-scalable {\em SAMRAI} architecture for the adaptive mesh refinement. We present the application of the code to several benchmark tests, showing the high order of convergence and accuracy achieved and the capabilities in terms of magnetic shock resolution and three-dimensionality. This paper paves the way for the applications to a realistic, 3D long-term evolution of neutron stars interior and, possibly, of other astrophysical sources. △ Less

Submitted 11 December, 2018; v1 submitted 20 November, 2018; originally announced November 2018.

Comments: 23 pages, 13 figures. In press

Journal ref: Computer Physics Communications, 2018

Anisotropic stars as ultracompact objects in General Relativity

Authors: Guilherme Raposo, Paolo Pani, Miguel Bezares, Carlos Palenzuela, Vitor Cardoso

Abstract: Anisotropic stresses are ubiquitous in nature, but their modeling in General Relativity is poorly understood and frame dependent. We introduce the first study on the dynamical properties of anisotropic self-gravitating fluids in a covariant framework. Our description is particularly useful in the context of tests of the black hole paradigm, wherein ultracompact objects are used as black hole mimic… ▽ More Anisotropic stresses are ubiquitous in nature, but their modeling in General Relativity is poorly understood and frame dependent. We introduce the first study on the dynamical properties of anisotropic self-gravitating fluids in a covariant framework. Our description is particularly useful in the context of tests of the black hole paradigm, wherein ultracompact objects are used as black hole mimickers but otherwise lack a proper theoretical framework. We show that: (i) anisotropic stars can be as compact and as massive as black holes, even for very small anisotropy parameters; (ii) the nonlinear dynamics of the 1+1 system is in good agreement with linearized calculations, and shows that configurations below the maximum mass are nonlinearly stable; (iii) strongly anisotropic stars have vanishing tidal Love numbers in the black-hole limit; (iv) their formation will usually be accompanied by gravitational-wave echoes at late times. △ Less

Submitted 15 May, 2019; v1 submitted 19 November, 2018; originally announced November 2018.

Comments: 7+2 pages, 6 figures; v2: include extra material (general covariant framework for anisotropic fluids in General Relativity without symmetries and code validation); to appear in PRD

Journal ref: Phys. Rev. D 99, 104072 (2019)

Black holes, gravitational waves and fundamental physics: a roadmap

Authors: Leor Barack, Vitor Cardoso, Samaya Nissanke, Thomas P. Sotiriou, Abbas Askar, Krzysztof Belczynski, Gianfranco Bertone, Edi Bon, Diego Blas, Richard Brito, Tomasz Bulik, Clare Burrage, Christian T. Byrnes, Chiara Caprini, Masha Chernyakova, Piotr Chrusciel, Monica Colpi, Valeria Ferrari, Daniele Gaggero, Jonathan Gair, Juan Garcia-Bellido, S. F. Hassan, Lavinia Heisenberg, Martin Hendry, Ik Siong Heng , et al. (181 additional authors not shown)

Abstract: The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horiz… ▽ More The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress. △ Less

Submitted 1 February, 2019; v1 submitted 13 June, 2018; originally announced June 2018.

Comments: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

Pulsar magnetospheres in General Relativity

Authors: Federico Carrasco, Carlos Palenzuela, Oscar Reula

Abstract: The main contribution to the pulsar power can be calculated by assuming a rotating magnetically-dominated magnetosphere described by the force-free approximation. Although this simple model has been used thoroughly to study pulsar magnetospheres in the flat spacetime regime, only few works have considered the relativistic corrections introduced by the curvature and frame-dragging effects induced b… ▽ More The main contribution to the pulsar power can be calculated by assuming a rotating magnetically-dominated magnetosphere described by the force-free approximation. Although this simple model has been used thoroughly to study pulsar magnetospheres in the flat spacetime regime, only few works have considered the relativistic corrections introduced by the curvature and frame-dragging effects induced by a rotating neutron star. Here we revisit the problem and describe pulsar magnetospheres within full General Relativity, quantifying the corrections as a function of the angular velocity, the compactness of the star and the misalignment angle between the spin and the magnetic dipole. We provide analytical expressions for the pulsar luminosity by fitting our numerical results. Finally, we also analyze the effect of the relativistic corrections on the braking index, which indicates a slight increment in its value. △ Less

Submitted 15 June, 2018; v1 submitted 10 May, 2018; originally announced May 2018.

Comments: 19 pages, 11 figures

Journal ref: Phys. Rev. D 98, 023010 (2018)

Neutron star mergers as a probe of modifications of general relativity with finite-range scalar forces

Authors: Laura Sagunski, Jun Zhang, Matthew C. Johnson, Luis Lehner, Mairi Sakellariadou, Steven L. Liebling, Carlos Palenzuela, David Neilsen

Abstract: Observations of gravitational radiation from compact binary systems provide an unprecedented opportunity to test General Relativity in the strong field dynamical regime. In this paper, we investigate how future observations of gravitational radiation from binary neutron star mergers might provide constraints on finite-range forces from a universally coupled massive scalar field. Such scalar degree… ▽ More Observations of gravitational radiation from compact binary systems provide an unprecedented opportunity to test General Relativity in the strong field dynamical regime. In this paper, we investigate how future observations of gravitational radiation from binary neutron star mergers might provide constraints on finite-range forces from a universally coupled massive scalar field. Such scalar degrees of freedom are a characteristic feature of many extensions of General Relativity. For concreteness, we work in the context of metric $f(R)$ gravity, which is equivalent to General Relativity and a universally coupled scalar field with a non-linear potential whose form is fixed by the choice of $f(R)$. In theories where neutron stars (or other compact objects) obtain a significant scalar charge, the resulting attractive finite-range scalar force has implications for both the inspiral and merger phases of binary systems. We first present an analysis of the inspiral dynamics in Newtonian limit, and forecast the constraints on the mass of the scalar and charge of the compact objects for the Advanced LIGO gravitational wave observatory. We then perform a comparative study of binary neutron star mergers in General Relativity with those of a one-parameter model of $f(R)$ gravity using fully relativistic hydrodynamical simulations. These simulations elucidate the effects of the scalar on the merger and post-merger dynamics. We comment on the utility of the full waveform (inspiral, merger, post-merger) to probe different regions of parameter space for both the particular model of $f(R)$ gravity studied here and for finite-range scalar forces more generally. △ Less

Submitted 19 February, 2018; v1 submitted 19 September, 2017; originally announced September 2017.

Comments: 22 pages, 8 figures. Typos fixed, references added. Updated to match the version accepted for publishing in PRD

Journal ref: Phys. Rev. D 97, 064016 (2018)

Echoes of ECOs: gravitational-wave signatures of exotic compact objects and of quantum corrections at the horizon scale

Authors: Vitor Cardoso, Seth Hopper, Caio F. B. Macedo, Carlos Palenzuela, Paolo Pani

Abstract: Gravitational waves from binary coalescences provide one of the cleanest signatures of the nature of compact objects. It has been recently argued that the post-merger ringdown waveform of exotic ultracompact objects is initially identical to that of a black-hole, and that putative corrections at the horizon scale will appear as secondary pulses after the main burst of radiation. Here we extend thi… ▽ More Gravitational waves from binary coalescences provide one of the cleanest signatures of the nature of compact objects. It has been recently argued that the post-merger ringdown waveform of exotic ultracompact objects is initially identical to that of a black-hole, and that putative corrections at the horizon scale will appear as secondary pulses after the main burst of radiation. Here we extend this analysis in three important directions: (i) we show that this result applies to a large class of exotic compact objects with a photon sphere for generic orbits in the test-particle limit; (ii) we investigate the late-time ringdown in more detail, showing that it is universally characterized by a modulated and distorted train of "echoes" of the modes of vibration associated with the photon sphere; (iii) we study for the first time equal-mass, head-on collisions of two ultracompact boson stars and compare their gravitational-wave signal to that produced by a pair of black-holes. If the initial objects are compact enough as to mimic a binary black-hole collision up to the merger, the final object exceeds the maximum mass for boson stars and collapses to a black-hole. This suggests that - in some configurations - the coalescence of compact boson stars might be almost indistinguishable from that of black-holes. On the other hand, generic configurations display peculiar signatures that can be searched for in gravitational-wave data as smoking guns of exotic compact objects. △ Less

Submitted 13 October, 2016; v1 submitted 30 August, 2016; originally announced August 2016.

Comments: 13 pages, RevTex4. v2: typo in equation 7 corrected, references added, to appear in PRD

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

Electromagnetic Luminosity of the Coalescence of Charged Black Hole Binaries

Authors: Steven L. Liebling, Carlos Palenzuela

Abstract: The observation of a possible electromagnetic counterpart by the Fermi GBM group to the aLIGO detection of the merger of a black hole binary has spawned a number of ideas about its source. Furthermore, observations of fast radio bursts (FRBs) have similarly resulted in a range of new models that might endow black hole binaries with electromagnetic signatures. In this context, even the unlikely ide… ▽ More The observation of a possible electromagnetic counterpart by the Fermi GBM group to the aLIGO detection of the merger of a black hole binary has spawned a number of ideas about its source. Furthermore, observations of fast radio bursts (FRBs) have similarly resulted in a range of new models that might endow black hole binaries with electromagnetic signatures. In this context, even the unlikely idea that astrophysical black holes may have significant charge is worth exploring, and here we present results from the simulation of weakly charged black holes as they orbit and merge. Our simulations suggest that a black hole binary with mass comparable to that observed in GW150914 could produce the level of electromagnetic luminosity observed by Fermi GBM ($10^{49}$ ergs/s) with a non-dimensional charge of $q \equiv Q/M = 10^{-4}$ assuming good radiative efficiency. However even a charge such as this is difficult to imagine avoiding neutralization long enough for the binary to produce its electromagnetic counterpart, and so this value would likely serve simply as an upper bound. On the other hand, one can equivalently consider the black holes as having acquired a magnetic monopole charge that would be easy to maintain and would generate an identical electromagnetic signature as the electric charges. The observation of such a binary would have significant cosmological implications, not the least of which would be an explanation for the quantization of charge itself. We also study such a magnetically charged binary in the force-free regime and find it much more radiative, reducing even further the requirements to produce the counterpart. △ Less

Submitted 10 September, 2016; v1 submitted 7 July, 2016; originally announced July 2016.

Comments: 9 pages, 8 figures, Updated to be consistent with the published version

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

The $m=1$ instability \& gravitational wave signal in binary neutron star mergers

Authors: Luis Lehner, Steven L. Liebling, Carlos Palenzuela, Patrick Motl

Abstract: We examine the development and detectability of the $m=1$ instability in the remnant of binary neutron star mergers. The detection of the gravitational mode associated with the $m=1$ degree of freedom could potentially reveal details of the equation of state. We analyze the post-merger epoch of simulations of both equal and non-equal mass neutron star mergers using three realistic, microphysical e… ▽ More We examine the development and detectability of the $m=1$ instability in the remnant of binary neutron star mergers. The detection of the gravitational mode associated with the $m=1$ degree of freedom could potentially reveal details of the equation of state. We analyze the post-merger epoch of simulations of both equal and non-equal mass neutron star mergers using three realistic, microphysical equations of state and neutrino cooling. Our studies show such an instability develops generically and within a short dynamical time to strengths that are comparable or stronger than the $m=2$ mode which is the strongest during the early post-merger stage. We estimate the signal to noise ratio that might be obtained for the $m=1$ mode and discuss the prospects for observing this signal with available Earth-based detectors. Because the $m=1$ occurs at roughly half the frequency of the more powerful $m=2$ signal and because it can potentially be long-lived, targeted searches could be devised to observe it. We estimate that with constant amplitude direct detection of the mode could occur up to a distance of roughly $14\,\mathrm{Mpc}$ whereas a search triggered by the inspiral signal could extend this distance to roughly $100\,\mathrm{Mpc}$. △ Less

Submitted 19 August, 2016; v1 submitted 8 May, 2016; originally announced May 2016.

Comments: 9 pages, 5 figures. Update to match published version; in particular the previous version used an incorrect noise curve which affected slightly the conclusions. This is now fixed

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

Unequal mass binary neutron star mergers and multimessenger signals

Authors: Luis Lehner, Steven L. Liebling, Carlos Palenzuela, O. L. Caballero, Evan O'Connor, Matthew Anderson, David Neilsen

Abstract: We study the merger of binary neutron stars with different mass ratios adopting three different realistic, microphysical nuclear equations of state, as well as incorporating neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with neutrino cooling, in determining the properties of… ▽ More We study the merger of binary neutron stars with different mass ratios adopting three different realistic, microphysical nuclear equations of state, as well as incorporating neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with neutrino cooling, in determining the properties of the resulting hypermassive neutron star, of the neutrinos produced, and of the ejected material. The ejecta we find are consistent with other recent studies that find that small mass ratios produce more ejecta than equal mass cases (up to some limit) and this ejecta is more neutron rich. This trend indicates the importance with future kilonovae observations of measuring the individual masses of an associated binary neutron star system, presumably from concurrent gravitational wave observations, in order to be able to extract information about the nuclear equation of state △ Less

Submitted 1 March, 2016; originally announced March 2016.

Comments: 16 pages, 14 figures (some figures reduced in quality due to size constraints)

Interaction between bosonic dark matter and stars

Authors: Richard Brito, Vitor Cardoso, Caio F. B. Macedo, Hirotada Okawa, Carlos Palenzuela

Abstract: We provide a detailed analysis of how bosonic dark matter "condensates" interact with compact stars, extending significantly the results of a recent Letter. We focus on bosonic fields with mass $m_B$, such as axions, axion-like candidates and hidden photons. Self-gravitating bosonic fields generically form "breathing" configurations, where both the spacetime geometry and the field oscillate, and c… ▽ More We provide a detailed analysis of how bosonic dark matter "condensates" interact with compact stars, extending significantly the results of a recent Letter. We focus on bosonic fields with mass $m_B$, such as axions, axion-like candidates and hidden photons. Self-gravitating bosonic fields generically form "breathing" configurations, where both the spacetime geometry and the field oscillate, and can interact and cluster at the center of stars. We construct stellar configurations formed by a perfect fluid and a bosonic condensate, and which may describe the late stages of dark-matter accretion onto stars, in dark matter-rich environments. These composite stars oscillate at a frequency which is a multiple of $f=2.5\times 10^{14}\,\left(m_{B}c^2/eV\right)\,{\rm Hz}$. Using perturbative analysis and Numerical Relativity techniques, we show that these stars are generically stable, and we provide criteria for instability. Our results also indicate that the growth of the dark matter core is halted close to the Chandrasekhar limit. We thus dispel a myth concerning dark matter accretion by stars: dark mater accretion does not necessarily lead to the destruction of the star, nor to collapse to a black hole. Finally, we argue that stars with long-lived bosonic cores may also develop in other theories with effective mass couplings, such as (massless) scalar-tensor theories. △ Less

Submitted 30 January, 2016; v1 submitted 1 December, 2015; originally announced December 2015.

Comments: 26 pages, RevTex4. v3. References added, minor corrections implemented. To appear in Physical Review D as an Editor's Suggestion

Journal ref: Phys. Rev. D 93, 044045 (2016)

Effects of the microphysical Equation of State in the mergers of magnetized Neutron Stars With Neutrino Cooling

Authors: C. Palenzuela, S. L. Liebling, D. Neilsen, L. Lehner, O. L. Caballero, E. O'Connor, M. Anderson

Abstract: We study the merger of binary neutron stars using different realistic, microphysical nuclear equations of state, as well as incorporating magnetic field and neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with cooling and electromagnetic effects, in determining the properties… ▽ More We study the merger of binary neutron stars using different realistic, microphysical nuclear equations of state, as well as incorporating magnetic field and neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with cooling and electromagnetic effects, in determining the properties of the hypermassive neutron star resulting from the merger, the production of neutrinos, and the characteristics of ejecta from the system. The ejecta we find are consistent with other recent studies that find soft equations of state produce more ejecta than stiffer equations of state. Moreover, the degree of neutron richness increases for softer equations of state. In light of reported kilonova observations (associated to GRB~130603B and GRB~060614) and the discovery of relatively low abundances of heavy, radioactive elements in deep sea deposits (with respect to possible production via supernovae), we speculate that a soft EoS might be preferred---because of its significant production of sufficiently neutron rich ejecta---if such events are driven by binary neutron star mergers. We also find that realistic magnetic field strengths, obtained with a sub-grid model tuned to capture magnetic amplification via the Kelvin-Helmholtz instability at merger, are generally too weak to affect the gravitational wave signature post-merger within a time scale of $\approx 10$~ms but can have subtle effects on the post-merger dynamics. △ Less

Submitted 7 May, 2015; originally announced May 2015.

Comments: 22 pages, 24 figures

Journal ref: Phys. Rev. D 92, 044045 (2015)

Electromagnetic outflows in a class of scalar-tensor theories: binary neutron star coalescence

Authors: Marcelo Ponce, Carlos Palenzuela, Enrico Barausse, Luis Lehner

Abstract: As we showed in previous work, the dynamics and gravitational emission of binary neutron star systems in certain scalar-tensor theories can differ significantly from that expected from General Relativity in the coalescing stage. In this work we examine whether the characteristics of the electromagnetic counterparts to these binaries -- driven by magnetosphere interactions prior to the merger event… ▽ More As we showed in previous work, the dynamics and gravitational emission of binary neutron star systems in certain scalar-tensor theories can differ significantly from that expected from General Relativity in the coalescing stage. In this work we examine whether the characteristics of the electromagnetic counterparts to these binaries -- driven by magnetosphere interactions prior to the merger event -- can provide an independent way to test gravity in the most strongly dynamical stages of binary mergers. We find that the electromagnetic flux emitted by binaries in these scalar-tensor theories can show deviations from the GR prediction in particular cases. These differences are quite subtle, thus requiring delicate measurements to differentiate between GR and the type of scalar-tensor theories considered in this work using electromagnetic observations alone. However, if coupled with a gravitational-wave detection, electromagnetic measurements might provide a way to increase the confidence with which GR will be confirmed (or ruled out) by gravitational observations. △ Less

Submitted 19 April, 2015; v1 submitted 2 October, 2014; originally announced October 2014.

Comments: Improved presentation and figures; matches version accepted in PRD

Journal ref: Phys. Rev. D 91, 084038 (2015)

Interaction of misaligned magnetospheres in the coalescence of binary neutron stars

Authors: Marcelo Ponce, Carlos Palenzuela, Luis Lehner, Steven L. Liebling

Abstract: We study the dependence of the electromagnetic luminosity --produced by interactions of force-free magnetospheres-- on dipole inclinations in binary neutron star systems. We show that this interaction extracts kinetic energy from the system and powers a Poynting flux with a strong dependence on the dipole orientations. This dependence can be linked to the reconnection and redistribution of magneti… ▽ More We study the dependence of the electromagnetic luminosity --produced by interactions of force-free magnetospheres-- on dipole inclinations in binary neutron star systems. We show that this interaction extracts kinetic energy from the system and powers a Poynting flux with a strong dependence on the dipole orientations. This dependence can be linked to the reconnection and redistribution of magnetic field as the stars interact. Although the details of the Poynting luminosity are very much dependent on the orientation, all the cases considered here nevertheless radiate a large Poynting flux. This robust emission suggests that the pre-merger stage of binary neutron star systems can yield interesting electromagnetic counterparts to gravitational wave events. △ Less

Submitted 5 August, 2014; v1 submitted 2 April, 2014; originally announced April 2014.

Comments: 12 pages, 9 figures; matches version published in PRD

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

Magnetized Neutron Stars With Realistic Equations of State and Neutrino Cooling

Authors: David Neilsen, Steven L. Liebling, Matthew Anderson, Luis Lehner, Evan O'Connor, Carlos Palenzuela

Abstract: We incorporate realistic, tabulated equations of state into fully relativistic simulations of magnetized neutron stars along with a neutrino leakage scheme which accounts for cooling via neutrino emission. Both these improvements utilize open-source code (GR1D) and tables from http://stellarcollapse.org. Our implementation makes use of a novel method for the calculation of the optical depth which… ▽ More We incorporate realistic, tabulated equations of state into fully relativistic simulations of magnetized neutron stars along with a neutrino leakage scheme which accounts for cooling via neutrino emission. Both these improvements utilize open-source code (GR1D) and tables from http://stellarcollapse.org. Our implementation makes use of a novel method for the calculation of the optical depth which simplifies its use with distributed adaptive mesh refinement. We present various tests with and without magnetization and preliminary results both from single stars and from the merger of a binary system. △ Less

Submitted 9 May, 2014; v1 submitted 14 March, 2014; originally announced March 2014.

Comments: 16 pages, 10 figures; Minor modifications consistent with version accepted by PRD

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

Dynamical scalarization of neutron stars in scalar-tensor gravity theories

Authors: Carlos Palenzuela, Enrico Barausse, Marcelo Ponce, Luis Lehner

Abstract: We present a framework to study generic neutron-star binaries in scalar-tensor theories of gravity. Our formalism achieves this goal by suitably interfacing a post-Newtonian orbital evolution (described by a set of ordinary differential equations) with a set of non-linear algebraic equations, which provide a description of the scalar charge of each binary's component along the evolution in terms o… ▽ More We present a framework to study generic neutron-star binaries in scalar-tensor theories of gravity. Our formalism achieves this goal by suitably interfacing a post-Newtonian orbital evolution (described by a set of ordinary differential equations) with a set of non-linear algebraic equations, which provide a description of the scalar charge of each binary's component along the evolution in terms of isolated-star data. We validate this semi-analytical procedure by comparing its results to those of fully general-relativistic simulations, and use it to investigate the behavior of binary systems in large portions of the parameter space of scalar-tensor theories. This allows us to shed further light on the phenomenon of "dynamical scalarization", which we uncovered in [Barausse, Palenzuela, Ponce and Lehner, Phys. Rev. D 87, 081506(R) (2013)] and which takes place in tight binaries, even for stars that have exactly zero scalar charge in isolation. We also employ our formalism to study representative binary systems, obtain their gravitational-wave signals and discuss the extent to which deviations from General Relativity can be detected. The insights gained by this framework allow us to additionally show that eccentric binaries can undergo scalarization/de-scalarization phenomena. △ Less

Submitted 6 March, 2014; v1 submitted 16 October, 2013; originally announced October 2013.

Comments: 17 pages, 20 figures; matches version published in PRD

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

Linking electromagnetic and gravitational radiation in coalescing binary neutron stars

Authors: Carlos Palenzuela, Luis Lehner, Steven L. Liebling, Marcelo Ponce, Matthew Anderson, David Neilsen, Patrick Motl

Abstract: We expand on our study of the gravitational and electromagnetic emissions from the late stage of an inspiraling neutron star binary as presented in Ref. \cite{Palenzuela:2013hu}. Interactions between the stellar magnetospheres, driven by the extreme dynamics of the merger, can yield considerable outflows. We study the gravitational and electromagnetic waves produced during the inspiral and merger… ▽ More We expand on our study of the gravitational and electromagnetic emissions from the late stage of an inspiraling neutron star binary as presented in Ref. \cite{Palenzuela:2013hu}. Interactions between the stellar magnetospheres, driven by the extreme dynamics of the merger, can yield considerable outflows. We study the gravitational and electromagnetic waves produced during the inspiral and merger of a binary neutron star system using a full relativistic, resistive MHD evolution code. We show that the interaction between the stellar magnetospheres extracts kinetic energy from the system and powers radiative Poynting flux and heat dissipation. These features depend strongly on the configuration of the initial stellar magnetic moments. Our results indicate that this power can strongly outshine pulsars in binaries and have a distinctive angular and time-dependent pattern. Our discussion provides more detail than Ref. \cite{Palenzuela:2013hu}, showing clear evidence of the different effects taking place during the inspiral. Our simulations include a few milliseconds after the actual merger and study the dynamics of the magnetic fields during the formation of the hypermassive neutron star. We also briefly discuss the possibility of observing such emissions. △ Less

Submitted 19 August, 2013; v1 submitted 28 July, 2013; originally announced July 2013.

Comments: 15 pages, 15 figures; Updated with spelling/grammar corrections, updated references, and better printing figures

Journal ref: Phys. Rev. D 88, 043011 (2013)

Gravitational and electromagnetic outputs from binary neutron star mergers

Authors: Carlos Palenzuela, Luis Lehner, Marcelo Ponce, Steven L. Liebling, Matthew Anderson, David Neilsen, Patrick Motl

Abstract: The late stage of an inspiraling neutron star binary gives rise to strong gravitational wave emission due to its highly dynamic, strong gravity. Moreover, interactions between the stellar magnetospheres can produce considerable electromagnetic radiation. We study this scenario using fully general relativistic, resistive magneto-hydrodynamics simulations. We show that these interactions extract kin… ▽ More The late stage of an inspiraling neutron star binary gives rise to strong gravitational wave emission due to its highly dynamic, strong gravity. Moreover, interactions between the stellar magnetospheres can produce considerable electromagnetic radiation. We study this scenario using fully general relativistic, resistive magneto-hydrodynamics simulations. We show that these interactions extract kinetic energy from the system, dissipate heat, and power radiative Poynting flux, as well as develop current sheets. Our results indicate that this power can: (i) outshine pulsars in binaries, (ii) display a distinctive angular- and time-dependent pattern, and (iii) radiate within large opening angles. These properties suggest that some binary neutron star mergers are ideal candidates for multimessenger astronomy. △ Less

Submitted 24 July, 2013; v1 submitted 29 January, 2013; originally announced January 2013.

Comments: 6 pages, 3 figures. Updated version to appear in Phys.Rev.Letters

Journal ref: PhysRevLett 111.061105 (2013)

Neutron-star mergers in scalar-tensor theories of gravity

Authors: Enrico Barausse, Carlos Palenzuela, Marcelo Ponce, Luis Lehner

Abstract: Scalar-tensor theories of gravity are natural phenomenological alternatives to General Relativity, where the gravitational interaction is mediated by a scalar degree of freedom, besides the usual tensor gravitons. In regions of the parameter space of these theories where constraints from both solar system experiments and binary-pulsar observations are satisfied, we show that binaries of neutron st… ▽ More Scalar-tensor theories of gravity are natural phenomenological alternatives to General Relativity, where the gravitational interaction is mediated by a scalar degree of freedom, besides the usual tensor gravitons. In regions of the parameter space of these theories where constraints from both solar system experiments and binary-pulsar observations are satisfied, we show that binaries of neutron stars present marked differences from General Relativity in both the late-inspiral and merger phases. In particular, phenomena related to the spontaneous scalarization of isolated neutron stars take place in the late stages of the evolution of binary systems, with important effects in the ensuing dynamics. We comment on the relevance of our results for the upcoming Advanced LIGO/Virgo detectors. △ Less

Submitted 24 April, 2013; v1 submitted 20 December, 2012; originally announced December 2012.

Comments: 5 pages, 4 figures. One simulation highlighting a new physical effect has been added (cf. Fig 4) and the text has been modified, to match the version accepted for publication in PRD Rapid Communications

Journal ref: Phys. Rev. D 87, 081506(R) (2013)

Modeling magnetized neutron stars using resistive MHD

Authors: Carlos Palenzuela

Abstract: This work presents an implementation of the resistive MHD equations for a generic algebraic Ohm's law which includes the effects of finite resistivity within full General Relativity. The implementation naturally accounts for magnetic-field-induced anisotropies and, by adopting a phenomenological current, is able to accurately describe electromagnetic fields in the star and in its magnetosphere. We… ▽ More This work presents an implementation of the resistive MHD equations for a generic algebraic Ohm's law which includes the effects of finite resistivity within full General Relativity. The implementation naturally accounts for magnetic-field-induced anisotropies and, by adopting a phenomenological current, is able to accurately describe electromagnetic fields in the star and in its magnetosphere. We illustrate the application of this approach in interesting systems with astrophysical implications; the aligned rotator solution and the collapse of a magnetized rotating neutron star to a black hole. △ Less

Submitted 1 December, 2012; originally announced December 2012.

Comments: 14 pages, 5 figures

Journal ref: Mon. Not. R. Aston. Soc., Volume 431, Issue 2, pag.1853-1865 (2013)

Dynamical evolution of fermion-boson stars

Authors: Susana Valdez-Alvarado, Carlos Palenzuela, Daniela Alic, L. Arturo Ureña-López

Abstract: Compact objects, like neutron stars and white dwarfs, may accrete dark matter, and then be sensitive probes of its presence. These compact stars with a dark matter component can be modeled by a perfect fluid minimally coupled to a complex scalar field (representing a bosonic dark matter component), resulting in objects known as fermion-boson stars. We have performed the dynamical evolution of thes… ▽ More Compact objects, like neutron stars and white dwarfs, may accrete dark matter, and then be sensitive probes of its presence. These compact stars with a dark matter component can be modeled by a perfect fluid minimally coupled to a complex scalar field (representing a bosonic dark matter component), resulting in objects known as fermion-boson stars. We have performed the dynamical evolution of these stars in order to analyze their stability, and to study their spectrum of normal modes, which may reveal the amount of dark matter in the system. Their stability analysis shows a structure similar to that of an isolated (fermion or boson) star, with equilibrium configurations either laying on the stable or on the unstable branch. The analysis of the spectrum of normal modes indicates the presence of new oscillation modes in the fermionic part of the star, which result from the coupling to the bosonic component through the gravity. △ Less

Submitted 1 November, 2015; v1 submitted 8 October, 2012; originally announced October 2012.

Comments: 13 pages, 9 figures. Replaced with the journal final version

Journal ref: Phys. Rev. D 87, 084040 (2013)

Dynamical Boson Stars

Authors: Steven L. Liebling, Carlos Palenzuela

Abstract: The idea of stable, localized bundles of energy has strong appeal as a model for particles. In the 1950s John Wheeler envisioned such bundles as smooth configurations of electromagnetic energy that he called {\em geons}, but none were found. Instead, particle-like solutions were found in the late 1960s with the addition of a scalar field, and these were given the name {\em boson stars}. Since then… ▽ More The idea of stable, localized bundles of energy has strong appeal as a model for particles. In the 1950s John Wheeler envisioned such bundles as smooth configurations of electromagnetic energy that he called {\em geons}, but none were found. Instead, particle-like solutions were found in the late 1960s with the addition of a scalar field, and these were given the name {\em boson stars}. Since then, boson stars find use in a wide variety of models as sources of dark matter, as black hole mimickers, in simple models of binary systems, and as a tool in finding black holes in higher dimensions with only a single killing vector. We discuss important varieties of boson stars, their dynamic properties, and some of their uses, concentrating on recent efforts. △ Less

Submitted 11 February, 2023; v1 submitted 26 February, 2012; originally announced February 2012.

Comments: 106 pages, 35 figures, invited review for Living Reviews in Relativity; major revision in 2022. Comments welcome

Journal ref: Living Rev. Relativity 26, (2023), 1

Intense Electromagnetic Outbursts from Collapsing Hypermassive Neutron Stars

Authors: Luis Lehner, Carlos Palenzuela, Steven L. Liebling, Christopher Thompson, Chad Hanna

Abstract: We study the gravitational collapse of a magnetized neutron star using a novel numerical approach able to capture both the dynamics of the star and the behavior of the surrounding plasma. In this approach, a fully general relativistic magnetohydrodynamics implementation models the collapse of the star and provides appropriate boundary conditions to a force-free model which describes the stellar ex… ▽ More We study the gravitational collapse of a magnetized neutron star using a novel numerical approach able to capture both the dynamics of the star and the behavior of the surrounding plasma. In this approach, a fully general relativistic magnetohydrodynamics implementation models the collapse of the star and provides appropriate boundary conditions to a force-free model which describes the stellar exterior. We validate this strategy by comparing with known results for the rotating monopole and aligned rotator solutions and then apply it to study both rotating and non-rotating stellar collapse scenarios, and contrast the behavior with what is obtained when employing the electrovacuum approximation outside the star. The non-rotating electrovacuum collapse is shown to agree qualitatively with a Newtonian model of the electromagnetic field outside a collapsing star. We illustrate and discuss a fundamental difference between the force-free and electrovacuum solutions, involving the appearance of large zones of electric-dominated field in the vacuum case. This provides a clear demonstration of how dissipative singularities appear generically in the non-linear time-evolution of force-free fluids. In both the rotating and non-rotating cases, our simulations indicate that the collapse induces a strong electromagnetic transient. In the case of sub-millisecond rotation, the magnetic field experiences strong winding and the transient carries much more energy. This result has important implications for models of gamma-ray bursts. △ Less

Submitted 5 September, 2012; v1 submitted 12 December, 2011; originally announced December 2011.

Comments: 28 pages, 20 figures (quality lowered to reduce sizes). Improved initial data and matching condition results in a lower, but still important, energy emission. Added appendix with a discussion on effects of transition layer

Journal ref: Phys. Rev. D86,104035 (2012)

On the detectability of dual jets from binary black holes

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

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

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

Comments: 4 pages, 3 figures

Journal ref: ApJ, 749, L32, 2012

Robustness of the Blandford-Znajek mechanism

Authors: Carlos Palenzuela, Carles Bona, Luis Lehner, Oscar Reula

Abstract: The Blandford-Znajek mechanism has long been regarded as a key ingredient in models attempting to explain powerful jets in AGNs, quasars, blazzars etc. In such mechanism, energy is extracted from a rotating black hole and dissipated at a load at far distances. In the current work we examine the behaviour of the BZ mechanism with respect to different boundary conditions, revealing the mechanism rob… ▽ More The Blandford-Znajek mechanism has long been regarded as a key ingredient in models attempting to explain powerful jets in AGNs, quasars, blazzars etc. In such mechanism, energy is extracted from a rotating black hole and dissipated at a load at far distances. In the current work we examine the behaviour of the BZ mechanism with respect to different boundary conditions, revealing the mechanism robustness upon variation of these conditions. Consequently, this work closes a gap in our understanding of this important scenario. △ Less

Submitted 17 February, 2011; originally announced February 2011.

Comments: 7 pages, accepted in CQG

Journal ref: Class. Quantum Grav. 28 (2011) 134007

Boosting jet power in black hole spacetimes

Authors: David Neilsen, Luis Lehner, Carlos Palenzuela, Eric W. Hirschmann, Steven L. Liebling, Patrick M. Motl, T. Garrett

Abstract: The extraction of rotational energy from a spinning black hole via the Blandford-Znajek mechanism has long been understood as an important component in models to explain energetic jets from compact astrophysical sources. Here we show more generally that the kinetic energy of the black hole, both rotational and translational, can be tapped, thereby producing even more luminous jets powered by the i… ▽ More The extraction of rotational energy from a spinning black hole via the Blandford-Znajek mechanism has long been understood as an important component in models to explain energetic jets from compact astrophysical sources. Here we show more generally that the kinetic energy of the black hole, both rotational and translational, can be tapped, thereby producing even more luminous jets powered by the interaction of the black hole with its surrounding plasma. We study the resulting Poynting jet that arises from single boosted black holes and binary black hole systems. In the latter case, we find that increasing the orbital angular momenta of the system and/or the spins of the individual black holes results in an enhanced Poynting flux. △ Less

Submitted 11 January, 2011; v1 submitted 27 December, 2010; originally announced December 2010.

Comments: 7 pages, 5 figures

Journal ref: Proceedings of the National Academy of Sciences, Vol. 108, no. 31, pp. 12641-12646 (August 2011)

Magnetospheres of Black Hole Systems in Force-Free Plasma

Authors: Carlos Palenzuela, Travis Garrett, Luis Lehner, Steven L. Liebling

Abstract: The interaction of black holes with ambient magnetic fields is important for a variety of highly energetic astrophysical phenomena. We study this interaction within the force-free approximation in which a tenuous plasma is assumed to have zero inertia. Blandford and Znajek (BZ) used this approach to demonstrate the conversion of some of the black hole's energy into electromagnetic Poynting flux in… ▽ More The interaction of black holes with ambient magnetic fields is important for a variety of highly energetic astrophysical phenomena. We study this interaction within the force-free approximation in which a tenuous plasma is assumed to have zero inertia. Blandford and Znajek (BZ) used this approach to demonstrate the conversion of some of the black hole's energy into electromagnetic Poynting flux in stationary and axisymmetric single black hole systems. We adopt this approach and extend it to examine asymmetric and, most importantly, dynamical systems by implementing the fully nonlinear field equations of general relativity coupled to Maxwell's equations. For single black holes, we study in particular the dependence of the Poynting flux and show that, even for misalignments between the black hole spin and the direction of the asymptotic magnetic field, a Poynting flux is generated with a luminosity dependent on such misalignment. For binary black hole systems, we show both in the head-on and orbiting cases that the moving black holes generate a Poynting flux. △ Less

Submitted 7 July, 2010; originally announced July 2010.

Comments: 13 pages, 11 figures

Journal ref: Phys.Rev.D82:044045,2010

Dual Jets from Binary Black Holes

Authors: Carlos Palenzuela, Luis Lehner, Steven L. Liebling

Abstract: Supermassive black holes are found at the centers of most galaxies and their inspiral is a natural outcome when galaxies merge. The inspiral of these systems is of utmost astrophysical importance as prodigious producers of gravitational waves and in their possible role in energetic electromagnetic events. We study such binary black hole coalescence under the influence of an external magnetic field… ▽ More Supermassive black holes are found at the centers of most galaxies and their inspiral is a natural outcome when galaxies merge. The inspiral of these systems is of utmost astrophysical importance as prodigious producers of gravitational waves and in their possible role in energetic electromagnetic events. We study such binary black hole coalescence under the influence of an external magnetic field produced by the expected circumbinary disk surrounding them. Solving the Einstein equations to describe the spacetime and using the force-free approach for the electromagnetic fields and the tenuous plasma, we present numerical evidence for possible jets driven by these systems. Extending the process described by Blandford and Znajek for a single spinning black hole, the picture that emerges suggests the electromagnetic field extracts energy from the orbiting black holes, which ultimately merge and settle into the standard Blandford-Znajek scenario. Emissions along dual and single jets would be expected that could be observable to large distances. △ Less

Submitted 6 May, 2010; originally announced May 2010.

Comments: 6 pages, 4 figures

Report number: pi-other-170

Journal ref: Science V329, 927 (2010)

EM counterparts of recoiling black holes: general relativistic simulations of non-Keplerian discs

Authors: Olindo Zanotti, Luciano Rezzolla, Luca Del Zanna, Carlos Palenzuela

Abstract: We investigate the dynamics of a circumbinary disc that responds to the loss of mass and to the recoil velocity of the black hole produced by the merger of a binary system of supermassive black holes. We perform the first two-dimensional general relativistic hydrodynamics simulations of \textit{extended} non-Keplerian discs and employ a new technique to construct a "shock detector", thus determini… ▽ More We investigate the dynamics of a circumbinary disc that responds to the loss of mass and to the recoil velocity of the black hole produced by the merger of a binary system of supermassive black holes. We perform the first two-dimensional general relativistic hydrodynamics simulations of \textit{extended} non-Keplerian discs and employ a new technique to construct a "shock detector", thus determining the precise location of the shocks produced in the accreting disc by the recoiling black hole. In this way we can study how the properties of the system, such as the spin, mass and recoil velocity of the black hole, affect the mass accretion rate and are imprinted on the electromagnetic emission from these sources. We argue that the estimates of the bremsstrahlung luminosity computed without properly taking into account the radiation transfer yield cooling times that are unrealistically short. At the same time we show, through an approximation based on the relativistic isothermal evolution, that the luminosity produced can reach a peak value above $L \simeq 10^{43} \ {\rm erg/s} $ at about $\sim 30\,{\rm d}$ after the merger of a binary with total mass $M\simeq 10^6 M_\odot$ and persist for several days at values which are a factor of a few smaller. If confirmed by more sophisticated calculations such a signal could indeed lead to an electromagnetic counterpart of the merger of binary black-hole system. △ Less

Submitted 12 September, 2010; v1 submitted 22 February, 2010; originally announced February 2010.

Comments: 17 pages, 11 figures, accepted by A&A, movies available at http://numrel.aei.mpg.de/Visualisations/Archive/BinaryBlackHoles/EMCounterparts/EMCounterparts.html

Journal ref: A&A 523, A8 (2010)

Evolutions of Magnetized and Rotating Neutron Stars

Authors: Steven L. Liebling, Luis Lehner, David Neilsen, Carlos Palenzuela

Abstract: We study the evolution of magnetized and rigidly rotating neutron stars within a fully general relativistic implementation of ideal magnetohydrodynamics with no assumed symmetries in three spatial dimensions. The stars are modeled as rotating, magnetized polytropic stars and we examine diverse scenarios to study their dynamics and stability properties. In particular we concentrate on the stability… ▽ More We study the evolution of magnetized and rigidly rotating neutron stars within a fully general relativistic implementation of ideal magnetohydrodynamics with no assumed symmetries in three spatial dimensions. The stars are modeled as rotating, magnetized polytropic stars and we examine diverse scenarios to study their dynamics and stability properties. In particular we concentrate on the stability of the stars and possible critical behavior. In addition to their intrinsic physical significance, we use these evolutions as further tests of our implementation which incorporates new developments to handle magnetized systems. △ Less

Submitted 24 May, 2010; v1 submitted 4 January, 2010; originally announced January 2010.

Comments: 12 pages, 8 figures

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

Vacuum Electromagnetic Counterparts of Binary Black-Hole Mergers

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

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

Submitted 11 December, 2009; originally announced December 2009.

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

Understanding possible electromagnetic counterparts to loud gravitational wave events: Binary black hole effects on electromagnetic fields

Authors: Carlos Palenzuela, Luis Lehner, Shin Yoshida

Abstract: In addition to producing loud gravitational waves (GW), the dynamics of a binary black hole system could induce emission of electromagnetic (EM) radiation by affecting the behavior of plasmas and electromagnetic fields in their vicinity. We here study how the electromagnetic fields are affected by a pair of orbiting black holes through the merger. In particular, we show how the binary's dynamics… ▽ More In addition to producing loud gravitational waves (GW), the dynamics of a binary black hole system could induce emission of electromagnetic (EM) radiation by affecting the behavior of plasmas and electromagnetic fields in their vicinity. We here study how the electromagnetic fields are affected by a pair of orbiting black holes through the merger. In particular, we show how the binary's dynamics induce a variability in possible electromagnetically induced emissions as well as an enhancement of electromagnetic fields during the late-merge and merger epochs. These time dependent features will likely leave their imprint in processes generating detectable emissions and can be exploited in the detection of electromagnetic counterparts of gravitational waves. △ Less

Submitted 19 November, 2009; originally announced November 2009.

Comments: 12 pages

Report number: pi-other-166

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