BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Carlos Palenzuela

Carlos Palenzuela appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.COhep-thastro-ph.SRhep-phastro-phastro-ph.GA

Catalog footprint

What is connected

41works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2022arXiv

No evidence of kinetic screening in simulations of merging binary neutron stars beyond general relativity

We have conducted fully relativistic simulations in a class of scalar-tensor theories with derivative self-interactions and screening of local scales. By using high-resolution shock-capturing methods and a non-vanishing shift vector, we have managed to avoid issues plaguing similar attempts in the past. We have first confirmed recent results by ourselves in spherical symmetry, obtained with an approximate approach and pointing at a partial breakdown of the screening in black-hole collapse. Then, we considered the late inspiral and merger of binary neutron stars. We found that screening tends to suppress the (subdominant) dipole scalar emission, but not the (dominant) quadrupole scalar mode. Our results point at quadrupole scalar signals as large as (or even larger than) in Fierz-Jordan-Brans-Dicke theories with the same conformal coupling, for strong-coupling scales in the MeV range that we can simulate.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dynamics of Screening in Modified Gravity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GRMHD large eddy simulations with gradient subgrid-scale model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Introduction to Numerical Relativity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electromagnetic Luminosity of the Coalescence of Charged Black Hole Binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Interaction between bosonic dark matter and stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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}$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unequal mass binary neutron star mergers and multimessenger signals

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraining scalar-tensor theories of gravity from the most massive neutron stars

Scalar-tensor~(ST) theories of gravity are natural phenomenological extensions to general relativity. Although these theories are severely constrained both by solar system experiments and by binary pulsar observations, a large set of ST families remain consistent with these observations. Recent work has suggested probing the unconstrained region of the parameter space of ST theories based on the stability properties of highly compact neutron stars. Here, the dynamical evolution of very compact stars in a fully nonlinear code demonstrates that the stars do become unstable and that the instability, in some cases, drives the stars to collapse. We discuss the implications of these results in light of recent observations of the most massive neutron star yet observed. In particular, such observations suggest that such a star would be subject to the instability for a certain regime; its existence therefore supports a bound on the ST parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamical evolution of fermion-boson stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamical scalarization of neutron stars in scalar-tensor gravity theories

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interaction of misaligned magnetospheres in the coalescence of binary neutron stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetized Neutron Stars With Realistic Equations of State and Neutrino Cooling

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Projected Constraints on Scalarization with Gravitational Waves from Neutron Star Binaries

Certain scalar-tensor theories have the property of endowing stars with scalar hair, sourced either by the star's own compactness (spontaneous scalarization) or, for binary systems, by the companion's scalar hair (induced scalarization) or by the orbital binding energy (dynamical scalarization). Scalarized stars in binaries present different conservative dynamics than in General Relativity, and can also excite a scalar mode in the metric perturbation that carries away dipolar radiation. As a result, the binary orbit shrinks faster than predicted in General Relativity, modifying the rate of decay of the orbital period. In spite of this, scalar-tensor theories can pass existing binary pulsar tests, because observed pulsars may not be compact enough or sufficiently orbitally bound to activate scalarization. Gravitational waves emitted during the last stages of compact binary inspirals are thus ideal probes of scalarization effects. For the standard projected sensitivity of advanced LIGO, we here show that, if neutron stars are sufficiently compact to enter the detector's sensitivity band already scalarized, then gravitational waves could place constraints at least comparable to binary pulsars. If the stars dynamically scalarize while inspiraling in band, then constraints are still possible provided the scalarization occurs sufficiently early in the inspiral, roughly below an orbital frequency of 50Hz. In performing these studies, we derive an easy-to-calculate data analysis measure, an integrated phase difference between a General Relativistic and a modified signal, that maps directly to the Bayes factor so as to determine whether a modified gravity effect is detectable. Finally, we find that custom-made templates are equally effective as model-independent, parameterized post-Einsteinian waveforms at detecting such modified gravity effects at realistic signal-to-noise ratios.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Almost-Killing conserved currents: a general mass function

A new class of conserved currents, describing non-gravitational energy-momentum density, is presented. The proposed currents do not require the existence of a (timelike) Killing vector, and are not restricted to spherically symmetric spacetimes (or similar ones, in which the Kodama vector can be defined). They are based instead on almost-Killing vectors, which could in principle be defined on generic spacetimes. We provide local arguments, based on energy density profiles in highly simplified (stationary, rigidly-rotating) star models, which confirm the physical interest of these 'almost-Killing currents'. A mass function is defined in this way for the spherical case, qualitatively different from the Hernández-Misner mass function. An elliptic equation determining the new mass function is derived for the Tolman-Bondi spherically symmetric dust metrics, including a simple solution for the Oppenheimer-Schneider collapse. The equations for the non-symmetric case are shown to be of a mixed elliptic-hyperbolic nature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

General-relativistic resistive magnetohydrodynamics in three dimensions: Formulation and tests

We present a new numerical implementation of the general-relativistic resistive magnetohydrodynamics (MHD) equations within the Whisky code. The numerical method adopted exploits the properties of implicit-explicit Runge-Kutta numerical schemes to treat the stiff terms that appear in the equations for large electrical conductivities. Using tests in one, two, and three dimensions, we show that our implementation is robust and recovers the ideal-MHD limit in regimes of very high conductivity. Moreover, the results illustrate that the code is capable of describing scenarios in a very wide range of conductivities. In addition to tests in flat spacetime, we report simulations of magnetized nonrotating relativistic stars, both in the Cowling approximation and in dynamical spacetimes. Finally, because of its astrophysical relevance and because it provides a severe testbed for general-relativistic codes with dynamical electromagnetic fields, we study the collapse of a nonrotating star to a black hole. We show that also in this case our results on the quasinormal mode frequencies of the excited electromagnetic fields in the Schwarzschild background agree with the perturbative studies within 0.7% and 5.6% for the real and the imaginary part of the l=1 mode eigenfrequency, respectively. Finally we provide an estimate of the electromagnetic efficiency of this process.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gravitational and electromagnetic outputs from binary neutron star mergers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Linking electromagnetic and gravitational radiation in coalescing binary neutron stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Neutron-star mergers in scalar-tensor theories of gravity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Conformal and covariant formulation of the Z4 system with constraint-violation damping

We present a new formulation of the Einstein equations based on a conformal and traceless decomposition of the covariant form of the Z4 system. This formulation combines the advantages of a conformal decomposition, such as the one used in the BSSNOK formulation (i.e. well-tested hyperbolic gauges, no need for excision, robustness to imperfect boundary conditions) with the advantages of a constraint-damped formulation, such as the generalized harmonic one (i.e. exponential decay of constraint violations when these are produced). We validate the new set of equations through standard tests and by evolving binary black hole systems. Overall, the new conformal formulation leads to a better behavior of the constraint equations and a rapid suppression of the violations when they occur. The changes necessary to implement the new conformal formulation in standard BSSNOK codes are very small as are the additional computational costs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Intense Electromagnetic Outbursts from Collapsing Hypermassive Neutron Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modeling magnetized neutron stars using resistive MHD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the detectability of dual jets from binary black holes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The role of the ergosphere in the Blandford-Znajek process

The Blandford-Znajek process, one of the most promising model for powering the relativistic jets from black holes, was initially introduced as a mechanism in which the magnetic fields extract energy from a rotating black hole. We study the evolution of force-free electromagnetic fields on regular spacetimes with an ergosphere, which are generated by rapidly rotating stars. Our conclusive results confirm previous works, claiming that the Blandford-Znajek mechanism is not directly related to the horizon of the black hole. We also show that the radiated energy depends exponentially on the compactness of the star.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Boosting jet power in black hole spacetimes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Robustness of the Blandford-Znajek mechanism

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dual Jets from Binary Black Holes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evolutions of Magnetized and Rotating Neutron Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Magnetospheres of Black Hole Systems in Force-Free Plasma

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multi-state Boson Stars

Motivated by the increasing interest in models which consider scalar fields as viable dark matter candidates, we have constructed a generalization of relativistic Boson Stars (BS) composed of two coexisting states of the scalar field, the ground state and the first excited state. We have studied the dynamical evolution of these Multi-state Boson Stars (MSBS) under radial perturbations, using numerical techniques. We show that stable MSBS can be constructed, when the number of particles in the first excited state, N2, is smaller than the number of particles in the ground state, N1. On the other hand, when N2 > N1, the configurations are initially unstable. However, they evolve and settle down into stable configurations. In the stabilization process, the initially ground state is excited and ends in a first excited state, whereas the initially first excited state ends in a ground state. During this process, both states emit scalar field radiation, decreasing their number of particles. This behavior shows that even though BS in the first excited state are intrinsically unstable under finite perturbations, the configuration resulting from the combination of this state with the ground state produces stable objects. Finally we show in a qualitative way, that stable MSBS could be realistic models of dark matter galactic halos, as they produce rotation curves that are flatter at large radii than the rotation curves produced by BS with only one state.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

One dimensional description of the gravitational perturbation in a Kerr background

The perturbation equation in a Kerr background is written as a coupled system of one dimensional equations for the different modes in the time domain. Numerical simulations show that the dominant mode in the gravitational response is the one corresponding to the mode of the initial perturbation, allowing us to conjecture that the coupling among the modes has a weak influence in our system of equations. We conclude that by neglecting the coupling terms it can be obtained a one dimensional harmonic equation which indeed describes with good accuracy the gravitational response from the Kerr black hole with low spin, while only few couplings are necessary to describe a high spin one. This result may help to understand the structure of test fields in a Kerr background and even to generate accurate waveforms for various cases in an efficient manner.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Signatures of the sources in the gravitational waves of a perturbed Schwarzschild black hole

The explicit form of perturbation equation for the $Ψ_4$ Weyl scalar, containing the matter source terms, is derived for general type D spacetimes. It is described in detail the particular case of the Schwarzschild spacetime using in-going penetrating coordinates. As a practical application, we focused on the emission of gravitational waves when a black hole is perturbed by a surrounding dust-like fluid matter. The symmetries of the spacetime and the simplicity of the matter source allow, by means of a spherical harmonic decomposition, to study the problem by means of a one dimensional numerical code.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Vacuum Electromagnetic Counterparts of Binary Black-Hole Mergers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Beyond ideal MHD: towards a more realistic modeling of relativistic astrophysical plasmas

Many astrophysical processes involving magnetic fields and quasi-stationary processes are well described when assuming the fluid as a perfect conductor. For these systems, the ideal-magnetohydrodynamics (MHD) description captures the dynamics effectively and a number of well-tested techniques exist for its numerical solution. Yet, there are several astrophysical processes involving magnetic fields which are highly dynamical and for which resistive effects can play an important role. The numerical modeling of such non-ideal MHD flows is significantly more challenging as the resistivity is expected to change of several orders of magnitude across the flow and the equations are then either of hyperbolic-parabolic nature or hyperbolic with stiff terms. We here present a novel approach for the solution of these relativistic resistive MHD equations exploiting the properties of implicit-explicit (IMEX) Runge Kutta methods. By examining a number of tests we illustrate the accuracy of our approach under a variety of conditions and highlight its robustness when compared with alternative methods, such as the Strang-splitting. Most importantly, we show that our approach allows one to treat, within a unified framework, both those regions of the flow which are fluid-pressure dominated (such as in the interior of compact objects) and those which are instead magnetic-pressure dominated (such as in their magnetospheres)

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Method to estimate ISCO and ring-down frequencies in binary systems and consequences for gravitational wave data analysis

Recent advances in the description of compact binary systems have produced gravitational waveforms that include inspiral, merger and ring-down phases. Comparing results from numerical simulations with those of post-Newtonian (PN), and related, expansions has provided motivation for employing PN waveforms in near merger epochs when searching for gravitational waves and has encouraged the development of analytic fits to full numerical waveforms. The models and simulations do not yet cover the full binary coalescence parameter space. For these yet un-simulated regions, data analysts can still conduct separate inspiral, merger and ring-down searches. Improved knowledge about the end of the inspiral phase, the beginning of the merger, and the ring-down frequencies could increase the efficiency of both coherent inspiral-merger-ring-down (IMR) searches and searches over each phase separately. Insight can be gained for all three cases through a recently presented theoretical calculation, which, corroborated by the numerical results, provides an implicit formula for the final spin of the merged black holes, accurate to within 10% over a large parameter space. Knowledge of the final spin allows one to predict the end of the inspiral phase and the quasinormal mode ring-down frequencies, and in turn provides information about the bandwidth and duration of the merger. In this work we will discuss a few of the implications of this calculation for data analysis.

0 citations0 reviewsNo code linkedOpen access