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

Luis Lehner

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

ResearcherUnclaimed source record
gr-qchep-thastro-ph.HEastro-ph.COastro-ph.SRnlin.CDastro-phcond-mat.stat-mechhep-phastro-ph.GAcond-mat.otherMachine Learningnucl-thphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

68works
15topics
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

68 published item(s)

preprint2022arXiv

Gamma-radiation sky maps from compact binaries

We study sky maps and light curves of gamma-ray emission from neutron stars in compact binaries, and in isolation. We briefly review some gamma-ray emission models, and reproduce sky maps from a standard isolated pulsar in the Separatrix Layer model. We consider isolated pulsars with several variations of a dipole magnetic field, including superpositions, and predict their gamma-ray emission. Our results provide new heuristics on what can and cannot be inferred about the magnetic field configuration of pulsars from high-energy observations. We find that typical double-peak light curves can be produced by pulsars with significant multipole structure beyond a single dipole. For binary systems, we also present a simple approximation that is useful for rapid explorations of binary magnetic field structure. Finally, we predict the gamma-ray emission pattern from a compact black hole-neutron star binary moments before merger by applying the Separatrix Layer model to data simulated in full general relativity; we find that face-on observers receive little emission, equatorial observers see one broad peak, and more generic observers typically see two peaks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Using machine learning to parametrize postmerger signals from binary neutron stars

There is growing interest in the detection and characterization of gravitational waves from postmerger oscillations of binary neutron stars. These signals contain information about the nature of the remnant and the high-density and out-of-equilibrium physics of the postmerger processes, which would complement any electromagnetic signal. However, the construction of binary neutron star postmerger waveforms is much more complicated than for binary black holes: (i) there are theoretical uncertainties in the neutron-star equation of state and other aspects of the high-density physics, (ii) numerical simulations are expensive and available ones only cover a small fraction of the parameter space with limited numerical accuracy, and (iii) it is unclear how to parametrize the theoretical uncertainties and interpolate across parameter space. In this work, we describe the use of a machine-learning method called a conditional variational autoencoder (CVAE) to construct postmerger models for hyper/massive neutron star remnant signals based on numerical-relativity simulations. The CVAE provides a probabilistic model, which encodes uncertainties in the training data within a set of latent parameters. We estimate that training such a model will ultimately require $\sim 10^4$ waveforms. However, using synthetic training waveforms as a proof-of-principle, we show that the CVAE can be used as an accurate generative model and that it encodes the equation of state in a useful latent representation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nonlinear effects in the black hole ringdown: absorption-induced mode excitation

Gravitational-wave observations of black hole ringdowns are commonly used to characterize binary merger remnants and to test general relativity. These analyses assume linear black hole perturbation theory, in particular that the ringdown can be described in terms of quasinormal modes even for times approaching the merger. Here we investigate a nonlinear effect during the ringdown, namely how a mode excited at early times can excite additional modes as it is absorbed by the black hole. This is a third-order secular effect: the change in the black-hole mass causes a shift in the mode spectrum, so that the original mode is projected onto the new ones. Using nonlinear simulations, we study the ringdown of a spherically-symmetric scalar field around an asymptotically anti-de Sitter black hole, and we find that this "absorption-induced mode excitation" (AIME) is the dominant nonlinear effect. We show that this effect takes place well within the nonadiabatic regime, so we can analytically estimate it using a sudden mass-change approximation. Adapting our estimation technique to asymptotically-flat Schwarzschild black holes, we expect AIME to play a role in the analysis and interpretation of current and future gravitational wave observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Critical collapse of a scalar field in semiclassical loop quantum gravity

We study the collapse in spherical symmetry of a massless scalar field minimally coupled to gravity using the semiclassical equations that are expected from loop quantum gravity. We find critical behavior of the mass as a function of the parameters of the initial data similar to that found by Choptuik in classical general relativity for a large set of initial data and values of the polymerization parameter. Contrary to wide expectations for quantum gravity, our semiclassical field equations have an exact scale invariance, as do the classical field equations. As one would then expect, we numerically find that the phase transition is second order, again as in the classical case.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Excited hairy black holes: dynamical construction and level transitions

We study the dynamics of unstable Reissner-Nordström anti-de Sitter black holes under charged scalar field perturbations in spherical symmetry. We unravel their general behavior and approach to the final equlibrium state. In the first part of this work, we present a numerical analysis of massive charged scalar field quasinormal modes. We identify the known mode families - superradiant modes, zero-damped modes, AdS modes, and the near-horizon mode - and we track their migration under variation of the black hole and field parameters. We show that the zero-damped modes become superradiantly unstable for large RNAdS with large gauge coupling; the leading unstable mode is identified with the near-horizon condensation instability. In the second part, we present results of numerical simulations of perturbed large RNAdS, showing the nonlinear development of these unstable modes. For generic initial conditions, charge and mass are ransferred from the black hole to the scalar field, until an equilibrium solution with a scalar condensate is reached. We use results from the linear analysis, however, to select special initial data corresponding to an unstable overtone mode. We find that these data evolve to produce a new equilibrium state - an excited hairy black hole with the scalar condensate in an overtone configuration. This state is, however, unstable, and the black hole eventually decays to the generic end state. Nevertheless, this demonstrates the potential relevance of overtone modes as transients in black hole dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

LRP2020: The cosmic origin and evolution of the elements

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Large-Misalignment Mechanism for the Formation of Compact Axion Structures: Signatures from the QCD Axion to Fuzzy Dark Matter

Axions are some of the best motivated particles beyond the Standard Model. We show how the attractive self-interactions of dark matter (DM) axions over a broad range of masses, from $10^{-22}$ eV to $10^7$ GeV, can lead to nongravitational growth of density fluctuations and the formation of bound objects. This structure formation enhancement is driven by parametric resonance when the initial field misalignment is large, and it affects axion density perturbations on length scales of order the Hubble horizon when the axion field starts oscillating, deep inside the radiation-dominated era. This effect can turn an otherwise nearly scale-invariant spectrum of adiabatic perturbations into one that has a spike at the aforementioned scales, producing objects ranging from dense DM halos to scalar-field configurations such as solitons and oscillons. We call this class of cosmological scenarios for axion DM production "the large-misalignment mechanism." We explore observational consequences of this mechanism for axions with masses up to $10$ eV. For axions heavier than $10^{-5}$ eV, the compact axion halos are numerous enough to significantly impact Earth-bound direct detection experiments, yielding intermittent but coherent signals with repetition rates exceeding one per decade and crossing times less than a day. These episodic increases in the axion density and kinematic coherence suggest new approaches for axion DM searches, including for the QCD axion. Dense structures made up of axions from $10^{-22}$ eV to $10^{-5}$ eV are detectable through gravitational lensing searches, and their gravitational interactions can also perturb baryonic structures and alter star formation. At very high misalignment amplitudes, the axion field can undergo self-interaction-induced implosions long before matter-radiation equality, producing potentially-detectable low-frequency stochastic gravitational waves.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gravitational action with null boundaries

We present a complete discussion of the boundary term in the action functional of general relativity when the boundary includes null segments in addition to the more usual timelike and spacelike segments. We confirm that ambiguities appear in the contribution from a null segment, because it depends on an arbitrary choice of parametrization for the generators. We also show that similar ambiguities appear in the contribution from a codimension-two surface at which a null segment is joined to another (spacelike, timelike, or null) segment. The parametrization ambiguity can be tamed by insisting that the null generators be affinely parametrized; this forces each null contribution to the boundary action to vanish, but leaves intact the fredom to rescale the affine parameter by a constant factor on each generator. Once a choice of parametrization is made, the ambiguity in the joint contributions can be eliminated by formulating well-motivated rules that ensure the additivity of the gravitational action. Enforcing these rules, we calculate the time rate of change of the action when it is evaluated for a so-called "Wheeler-deWitt patch" of a black hole in asymptotically-anti de Sitter space. We recover a number of results cited in the literature, obtained with a less complete analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlinear Evolution and Final Fate of Charged Anti-de Sitter Black Hole Superradiant Instability

We describe the full nonlinear development of the superradiant instability for a charged massless scalar field, coupled to general relativity and electromagnetism, in the vicinity of a Reissner-Nordstrom-AdS black hole. The presence of the negative cosmological constant provides a natural context for considering perfectly reflecting boundary conditions and studying the dynamics as the scalar field interacts repeatedly with the black hole. At early times, small superradiant perturbations grow as expected from linearized studies. Backreaction then causes the black hole to lose charge and mass until the perturbation becomes nonsuperradiant, with the final state described by a stable hairy black hole. For large gauge coupling, the instability extracts a large amount of charge per unit mass, resulting in greater entropy increase. We discuss the implications of the observed behavior for the general problem of superradiance in black hole spacetimes.

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

Towards an understanding of the force-free magnetosphere of rapidly spinning black holes

The ability of a plasma surrounding spinning black holes to extract rotational energy and power energetic emissions has been recognized as a key astrophysical phenomenon. Important insights into the nature of this process are obtained through the analysis of the interplay between a force-free magnetosphere and the black hole. This task involves solving a complicated system of equations, often requiring complex numerical simulations. Recent analytical attempts at tackling this problem have exploited the fact that the near horizon region of extreme Kerr (NHEK) is endowed with an enhanced symmetry group. We continue in this direction and show that for some conformally self-similar solutions, the NHEK force-free equations reduce to a single non-linear ordinary differential equation which is difficult to solve with straightforward integration. We here introduce a new approach specifically tailored to this problem and describe how one can obtain physically meaningful solutions.

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

Conserved quantities and dual turbulent cascades in Anti-de Sitter spacetime

We consider the dynamics of a spherically symmetric massless scalar field coupled to general relativity in Anti--de Sitter spacetime in the small-amplitude limit. Within the context of our previously developed two time framework (TTF) to study the leading self-gravitating effects, we demonstrate the existence of two new conserved quantities in addition to the known total energy $E$ of the modes: The particle number $N$ and Hamiltonian $H$ of our TTF system. Simultaneous conservation of $E$ and $N$ implies that weak turbulent processes undergo dual cascades (direct cascade of $E$ and inverse cascade of $N$ or vice versa). This partially explains the observed dynamics of 2-mode initial data. In addition, conservation of $E$ and $N$ limits the region of phase space that can be explored within the TTF approximation and in particular rules out equipartion of energy among the modes for general initial data. Finally, we discuss possible effects of conservation of $N$ and $E$ on late time dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coupled Oscillator Model for Nonlinear Gravitational Perturbations

Motivated by the gravity/fluid correspondence, we introduce a new method for characterizing nonlinear gravitational interactions. Namely we map the nonlinear perturbative form of the Einstein equation to the equations of motion of a collection of nonlinearly-coupled harmonic oscillators. These oscillators correspond to the quasinormal or normal modes of the background spacetime. We demonstrate the mechanics and the utility of this formalism within the context of perturbed asymptotically anti-de Sitter black brane spacetimes. We confirm in this case that the boundary fluid dynamics are equivalent to those of the hydrodynamic quasinormal modes of the bulk spacetime. We expect this formalism to remain valid in more general spacetimes, including those without a fluid dual. In other words, although borne out of the gravity/fluid correspondence, the formalism is fully independent and it has a much wider range of applicability. In particular, as this formalism inspires an especially transparent physical intuition, we expect its introduction to simplify the often highly technical analytical exploration of nonlinear gravitational dynamics.

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
preprint2015arXiv

Islands of stability and recurrence times in AdS

We study the stability of anti-de Sitter (AdS) spacetime to spherically symmetric perturbations of a real scalar field in general relativity. Further, we work within the context of the "two time framework" (TTF) approximation, which describes the leading nonlinear effects for small amplitude perturbations, and is therefore suitable for studying the weakly turbulent instability of AdS---including both collapsing and non-collapsing solutions. We have previously identified a class of quasi-periodic (QP) solutions to the TTF equations, and in this work we analyze their stability. We show that there exist several families of QP solutions that are stable to linear order, and we argue that these solutions represent islands of stability in TTF. We extract the eigenmodes of small oscillations about QP solutions, and we use them to predict approximate recurrence times for generic non-collapsing initial data in the full (non-TTF) system. Alternatively, when sufficient energy is driven to high-frequency modes, as occurs for initial data far from a QP solution, the TTF description breaks down as an approximation to the full system. Depending on the higher order dynamics of the full system, this often signals an imminent collapse to a black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetosphere of a Kerr black hole immersed in magnetized plasma and its perturbative mode structure

This work studies jet-like electromagnetic configurations surrounding a slowly-spinning black-hole immersed in a uniformly magnetized force-free plasma. In the first part of this work, we present a family of stationary solutions that are jet-capable. While these solutions all satisfy the force-free equations and the appropriate boundary conditions, our numerical experiments show a unique relaxed state starting from different initial data, and so one member of the family is likely preferred over the others. In the second part of this work, we analyze the perturbations of this family of jet-like solutions, and show that the perturbative modes exhibit a similar split into the trapped and traveling categories previously found for perturbed Blandford-Znajek solutions. In the eikonal limit, the trapped modes can be identified with the fast magnetosonic waves in the force-free plasma and the traveling waves are essentially the Alfven waves. Moreover, within the scope of our analysis, we have not seen signs of unstable modes at the light-crossing timescale of the system, within which the numerical relaxation process occurs. This observation disfavors mode instability as the selection mechanism for picking out a preferred solution. Consequently, our analytical study is unable to definitively select a particular solution out of the family to serve as the aforementioned preferred final state. This remains an interesting open problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Scaling Relations in Two-Dimensional Relativistic Hydrodynamic Turbulence

We derive exact scaling relations for two-dimensional relativistic hydrodynamic turbulence in the inertial range of scales. We consider both the energy cascade towards large scales and the enstrophy cascade towards small scales. We illustrate these relations by numerical simulations of turbulent weakly compressible flows. Intriguingly, the fluid-gravity correspondence implies that the gravitational field in spacetimes with anti-de Sitter asymptotics should exhibit similar scaling relations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Small black holes in $AdS_5\times S^5$

We consider small black holes in $AdS_5\times S^5$, smeared on $S^5$. We compute the spectrum of $\ell \in [1,10]$ $S^5$-quasinormal modes corresponding to fluctuations leading to localization of these black holes on $S^5$. We recover the zero mode found by Hubeny and Rangamani (HR) previously \cite{Hubeny:2002xn}, and explicitly demonstrate that a Gregory-Laflamme type instability is at play in this system. The instability is associated with the expectation value of a dimension-5 operator.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Universality of non-equilibrium dynamics of CFTs from holography

Motivated by a low-energy effective description of gauge theory/string theory duality, we conjecture that the dynamics of $SO(4)$-invariant states in a large class of four-dimensional conformal gauge theories on $S^3$ with non-equal central charges $c\ne a$ are universal on time scales $t_{\rm universal}\propto ({\cal E}-{\cal E}_{\rm vacuum})^{-1}$, in the limit where the energy ${\cal E}\to {\cal E}_{\rm vacuum}$. We show that low-energy excitations in $c\ne a$ CFTs do not thermalize in this limit. The holographic universality conjecture then implies that within the Einstein-scalar field system (dual to theories with $c=a$), $AdS_5$ is stable to spherically symmetric perturbations against formation of trapped surfaces within time scales $t_{\rm universal}$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Collisions of oppositely charged black holes

The first fully non-linear numerical simulations of colliding charged black holes in D=4 Einstein-Maxwell theory were recently reported arXiv:1205.1063. These collisions were performed for black holes with equal charge-to-mass ratio, for which initial data can be found in closed analytic form. Here we generalize the study of collisions of charged black holes to the case of unequal charge-to-mass ratios. We focus on oppositely charged black holes, as to maximize acceleration-dependent effects. As |Q|/M increases from 0 to 0.99, we observe that the gravitational radiation emitted increases by a factor of ~ 2.7; the electromagnetic radiation emission becomes dominant for |Q|/M >~ 0.37 and at |Q|/M=0.99 is larger, by a factor of ~ 5.8, than its gravitational counterpart. We observe that these numerical results exhibit a precise and simple scaling with the charge. Furthermore, we show that the results from the numerical simulations are qualitatively captured by a simple analytic model that computes the electromagnetic dipolar radiation and the gravitational quadrupolar radiation of two non-relativistic interacting particles in Minkowski spacetime.

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

Holographic Thermalization, stability of AdS, and the Fermi-Pasta-Ulam-Tsingou paradox

For a real massless scalar field in general relativity with a negative cosmological constant, we uncover a large class of spherically symmetric initial conditions that are close to AdS, but whose numerical evolution does not result in black hole formation. According to the AdS/CFT dictionary, these bulk solutions are dual to states of a strongly interacting boundary CFT that fail to thermalize at late times. Furthermore, as these states are not stationary, they define dynamical CFT configurations that do not equilibrate. We develop a two-timescale perturbative formalism that captures both direct and inverse cascades of energy and agrees with our fully nonlinear evolutions in the appropriate regime. We also show that this formalism admits a large class of quasi-periodic solutions. Finally, we demonstrate a striking parallel between the dynamics of AdS and the classic Fermi-Pasta-Ulam-Tsingou problem.

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

Numerical Relativity and Astrophysics

Throughout the Universe many powerful events are driven by strong gravitational effects that require general relativity to fully describe them. These include compact binary mergers, black hole accretion and stellar collapse, where velocities can approach the speed of light, and extreme gravitational fields --$Φ_{\rm Newt}/c^2 \simeq 1$-- mediate the interactions. Many of these processes trigger emission across a broad range of the electromagnetic spectrum. Compact binaries further source strong gravitational wave emission that could directly be detected in the near future. This feat will open up a gravitational wave window into our Universe and revolutionize its understanding. Describing these phenomena requires general relativity, and --where dynamical effects strongly modify gravitational fields-- the full Einstein equations coupled to matter sources. Numerical relativity is a field within general relativity concerned with studying such scenarios that cannot be accurately modeled via perturbative or analytical calculations. In this review, we examine results obtained within this discipline, with a focus on its impact in astrophysics.

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
preprint2014arXiv

Quasinormal Modes Beyond Kerr

The quasinormal modes (QNMs) of a black hole spacetime are the free, decaying oscillations of the spacetime, and are well understood in the case of Kerr black holes. We discuss a method for computing the QNMs of spacetimes which are slightly deformed from Kerr. We mention two example applications: the parametric, turbulent instability of scalar fields on a background which includes a gravitational QNM, and the shifts to the QNM frequencies of Kerr when the black hole is weakly charged. This method may be of use in studies of black holes which are deformed by external fields or are solutions to alternative theories of gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating the universe(s): from cosmic bubble collisions to cosmological observables with numerical relativity

The theory of eternal inflation in an inflaton potential with multiple vacua predicts that our universe is one of many bubble universes nucleating and growing inside an ever-expanding false vacuum. The collision of our bubble with another could provide an important observational signature to test this scenario. We develop and implement an algorithm for accurately computing the cosmological observables arising from bubble collisions directly from the Lagrangian of a single scalar field. We first simulate the collision spacetime by solving Einstein's equations, starting from nucleation and ending at reheating. Taking advantage of the collision's hyperbolic symmetry, simulations are performed with a 1+1-dimensional fully relativistic code that uses adaptive mesh refinement. We then calculate the comoving curvature perturbation in an open Friedmann-Robertson-Walker universe, which is used to determine the temperature anisotropies of the cosmic microwave background radiation. For a fiducial Lagrangian, the anisotropies are well described by a power law in the cosine of the angular distance from the center of the collision signature. For a given form of the Lagrangian, the resulting observational predictions are inherently statistical due to stochastic elements of the bubble nucleation process. Further uncertainties arise due to our imperfect knowledge about inflationary and pre-recombination physics. We characterize observational predictions by computing the probability distributions over four phenomenological parameters which capture these intrinsic and model uncertainties. This represents the first fully-relativistic set of predictions from an ensemble of scalar field models giving rise to eternal inflation, yielding significant differences from previous non-relativistic approximations. Thus, our results provide a basis for a rigorous confrontation of these theories with cosmological data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Testing the nonlinear stability of Kerr-Newman black holes

The nonlinear stability of Kerr-Newman black holes (KNBHs) is investigated by performing numerical simulations within the full Einstein-Maxwell theory. We take as initial data a KNBH with mass $M$, angular momentum to mass ratio $a$ and charge $Q$. Evolutions are performed to scan this parameter space within the intervals $0\le a/M\le 0.994$ and $0\le Q/M\le 0.996$, corresponding to an extremality parameter $a/a_{\rm max}$ ($a_{\rm max} \equiv \sqrt{M^2-Q^2}$) ranging from $0$ to $0.995$. These KNBHs are evolved, together with a small bar-mode perturbation, up to a time of order $120M$. Our results suggest that for small $Q/a$, the quadrupolar oscillation modes depend solely on $a/a_{\rm max}$, a universality also apparent in previous perturbative studies in the regime of small rotation. Using as a stability criterion the absence of significant relative variations in the horizon areal radius and BH spin, we find no evidence for any developing instability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Turbulent Black Holes

We show that rapidly-spinning black holes can display turbulent gravitational behavior which is mediated by a new type of parametric instability. This instability transfers energy from higher temporal and azimuthal spatial frequencies to lower frequencies--- a phenomenon reminiscent of the inverse energy cascade displayed by 2+1-dimensional turbulent fluids. Our finding reveals a path towards gravitational turbulence for perturbations of rapidly-spinning black holes, and provides the first evidence for gravitational turbulence in an asymptotically flat spacetime. Interestingly, this finding predicts observable gravitational wave signatures from such phenomena in black hole binaries with high spins and gives a gravitational description of turbulence relevant to the fluid-gravity duality.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Holographic Path to the Turbulent Side of Gravity

We study the dynamics of a 2+1 dimensional relativistic viscous conformal fluid in Minkowski spacetime. Such fluid solutions arise as duals, under the "gravity/fluid correspondence", to 3+1 dimensional asymptotically anti-de Sitter (AAdS) black brane solutions to the Einstein equation. We examine stability properties of shear flows, which correspond to hydrodynamic quasinormal modes of the black brane. We find that, for sufficiently high Reynolds number, the solution undergoes an inverse turbulent cascade to long wavelength modes. We then map this fluid solution, via the gravity/fluid duality, into a bulk metric. This suggests a new and interesting feature of the behavior of perturbed AAdS black holes and black branes, which is not readily captured by a standard quasinormal mode analysis. Namely, for sufficiently large perturbed black objects (with long-lived quasinormal modes), nonlinear effects transfer energy from short to long wavelength modes via a turbulent cascade within the metric perturbation. As long wavelength modes have slower decay, this lengthens the overall lifetime of the perturbation. We also discuss various implications of this behavior, including expectations for higher dimensions, and the possibility of predicting turbulence in more general gravitational scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Post-Newtonian approach to black hole-fluid systems

This work devises a formalism to obtain the equations of motion for a black hole-fluid configuration. Our approach is based on a Post-Newtonian expansion and adapted to scenarios where obtaining the relevant dynamics requires long time-scale evolutions. These systems are typically studied with Newtonian approaches, which have the advantage that larger time-steps can be employed than in full general-relativistic simulations, but have the downside that important physical effects are not accounted for. The formalism presented here provides a relatively straightforward way to incorporate those effects in existing implementations, up to 2.5PN order, with lower computational costs than fully relativistic simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Boson Stars in AdS

We construct boson stars in global Anti de Sitter (AdS) space and study their stability. Linear perturbation results suggest that the ground state along with the first three excited state boson stars are stable. We evolve some of these solutions and study their nonlinear stability in light of recent work \cite{Bizon:2011gg} arguing that a weakly turbulent instability drives scalar perturbations of AdS to black hole formation. However evolutions suggest that boson stars are nonlinearly stable and immune to the instability for sufficiently small perturbation. Furthermore, these studies find other families of initial data which similarly avoid the instability for sufficiently weak parameters. Heuristically, we argue that initial data families with widely distributed mass-energy distort the spacetime sufficiently to oppose the coherent amplification favored by the instability. From the dual CFT perspective our findings suggest that there exist families of rather generic initial conditions in strongly coupled CFT (with large number of degrees of freedom) that do not thermalize in the infinite future.

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

Holographic thermalization, quasinormal modes and superradiance in Kerr-AdS

Black holes in anti-de Sitter (AdS) backgrounds play a pivotal role in the gauge/gravity duality where they determine, among other things, the approach to equilibrium of the dual field theory. We undertake a detailed analysis of perturbed Kerr-AdS black holes in four- and five-dimensional spacetimes, including the computation of its quasinormal modes, hydrodynamic modes and superradiantly unstable modes. Our results shed light on the possibility of new black hole phases with a single Killing field, possible new holographic phenomena and phases in the presence of a rotating chemical potential, and close a crucial gap in our understanding of linearized perturbations of black holes in anti-de Sitter scenarios.

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
preprint2013arXiv

Quantum quenches of holographic plasmas

We employ holographic techniques to study quantum quenches at finite temperature, where the quenches involve varying the coupling of the boundary theory to a relevant operator with an arbitrary conformal dimension $2\leq\D\leq4$. The evolution of the system is studied by evaluating the expectation value of the quenched operator and the stress tensor throughout the process. The time dependence of the new coupling is characterized by a fixed timescale and the response of the observables depends on the ratio of the this timescale to the initial temperature. The observables exhibit universal scaling behaviours when the transitions are either fast or slow, i.e. when this ratio is very small or very large. The scaling exponents are smooth functions of the operator dimension. We find that in fast quenches, the relaxation time is set by the thermal timescale regardless of the operator dimension or the precise quenching rate.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Transient Gravitational-Wave Sky

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Collisions of charged black holes

We perform fully non-linear numerical simulations of charged-black-hole collisions, described by the Einstein-Maxwell equations, and contrast the results against analytic expectations. We focus on head-on collisions of non-spinning black holes, starting from rest and with the same charge to mass ratio, Q/M. The addition of charge to black holes introduces a new interesting channel of radiation and dynamics, most of which seem to be captured by Newtonian dynamics and flat-space intuition. The waveforms can be qualitatively described in terms of three stages; (i) an infall phase prior to the formation of a common apparent horizon; (ii) a nonlinear merger phase which corresponds to a peak in gravitational and electromagnetic energy; (iii) the ringdown marked by an oscillatory pattern with exponentially decaying amplitude and characteristic frequencies that are in good agreement with perturbative predictions. We observe that the amount of gravitational-wave energy generated throughout the collision decreases by about three orders of magnitude as the charge-to-mass ratio Q/M is increased from 0 to 0.98. We interpret this decrease as a consequence of the smaller accelerations present for larger values of the charge. In contrast, the ratio of energy carried by electromagnetic to gravitational radiation increases, reaching about 22% for the maximum Q/M ratio explored, which is in good agreement with analytic predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Determining the outcome of cosmic bubble collisions in full General Relativity

Cosmic bubble collisions provide an important possible observational window on the dynamics of eternal inflation. In eternal inflation, our observable universe is contained in one of many bubbles formed from an inflating metastable vacuum. The collision between bubbles can leave a detectable imprint on the cosmic microwave background radiation. Although phenomenological models of the observational signature have been proposed, to make the theory fully predictive one must determine the bubble collision spacetime, and thus the cosmological observables, from a scalar field theory giving rise to eternal inflation. Because of the intrinsically non-linear nature of the bubbles and their collision, this requires a numerical treatment incorporating General Relativity. In this paper, we present results from numerical simulations of bubble collisions in full General Relativity. These simulations allow us to accurately determine the outcome of bubble collisions, and examine their effect on the cosmology inside a bubble universe. We confirm the validity of a number of approximations used in previous analytic work, and identify qualitatively new features of bubble collision spacetimes. Both vacuum bubbles and bubbles containing a realistic inflationary cosmology are studied. We identify the constraints on the scalar field potential that must be satisfied in order to obtain collisions that are consistent with our observed cosmology, yet leave detectable signatures.

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

NR/HEP: roadmap for the future

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Scalar Collapse in AdS

Recently, studies of the gravitational collapse of a scalar field within spherically symmetric AdS spacetimes was presented in \cite{Bizon:2011gg,Jalmuzna:2011qw} which showed an instability of pure AdS to black hole formation. In particular, the work showed that arbitrarily small initial configurations of scalar field evolved through some number of reflections off the AdS boundary until a black hole forms. We consider this same system, extended to include a complex scalar field, and reproduce this phenomena. We present tests of our numerical code that demonstrate convergence and consistency. We study the properties of the evolution as the scalar pulse becomes more compact examining the asymptotic behavior of the scalar field, an observable in the corresponding boundary CFT. We demonstrate that such BH formation occurs even when one places a reflecting boundary at finite radius indicating that the sharpening is a property of gravity in a bounded domain, not of AdS itself. We examine how the initial energy is transferred to higher frequencies --which leads to black hole formation-- and uncover interesting features of this transfer.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermal quenches in N=2* plasmas

We exploit gauge/gravity duality to study `thermal quenches' in a plasma of the strongly coupled N=2* gauge theory. Specifically, we consider the response of an initial thermal equilibrium state of the theory under variations of the bosonic or fermionic mass, to leading order in m/T<<1. When the masses are made to vary in time, novel new counterterms must be introduced to renormalize the boundary theory. We consider transitions the conformal super-Yang-Mills theory to the mass deformed gauge theory and also the reverse transitions. By construction, these transitions are controlled by a characteristic time scale \calt and we show how the response of the system depends on the ratio of this time scale to the thermal time scale 1/T. The response shows interesting scaling behaviour both in the limit of fast quenches with T\calt<<1 and slow quenches with T\calt>>1. In the limit that T\calt\to\infty, we observe the expected adiabatic response. For fast quenches, the relaxation to the final equilibrium is controlled by the lowest quasinormal mode of the bulk scalar dual to the quenched operator. For slow quenches, the system relaxes with a (nearly) adiabatic response that is governed entirely by the late time profile of the mass. We describe new renormalization scheme ambiguities in defining gauge invariant observables for the theory with time dependant couplings.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Turbulent flows for relativistic conformal fluids in 2+1 dimensions

We demonstrate that relativistic conformal hydrodynamics in 2+1 dimensions displays a turbulent behaviour which cascades energy to longer wavelengths on both flat and spherical manifolds. Our motivation for this study is to understand the implications for gravitational solutions through the AdS/CFT correspondence. The observed behaviour implies gravitational perturbations of the corresponding black brane/black hole spacetimes (for sufficiently large scales/temperatures) will display a similar cascade towards longer wavelengths.

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

Final State of Gregory-Laflamme Instability

We describe the behavior of a perturbed 5-dimensional black string subject to the Gregory-Laflamme instability. We show that the horizon evolves in a self-similar manner, where at any moment in the late-time development of the instability the horizon can be described as a sequence of 3-dimensional spherical black holes of varying size, joined by black string segments of similar radius. As with the initial black string, each local string segment is itself unstable, and this fuels the self-similar cascade to (classically) arbitrarily small scales; in the process the horizon develops a fractal structure. In finite asymptotic time, the remaining string segments shrink to zero-size, yielding a naked singularity. Since no fine-tuning is required to excite the instability, this constitutes a generic violation of cosmic censorship. We further discuss how this behavior is related to satellite formation in low-viscosity fluid streams subject to the Rayleigh-Plateau instability, and estimate the fractal dimension of the horizon prior to formation of the naked singularity.

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

Black Strings, Low Viscosity Fluids, and Violation of Cosmic Censorship

We describe the behavior of 5-dimensional black strings, subject to the Gregory-Laflamme instability. Beyond the linear level, the evolving strings exhibit a rich dynamics, where at intermediate stages the horizon can be described as a sequence of 3-dimensional spherical black holes joined by black string segments. These segments are themselves subject to a Gregory-Laflamme instability, resulting in a self-similar cascade, where ever-smaller satellite black holes form connected by ever-thinner string segments. This behavior is akin to satellite formation in low-viscosity fluid streams subject to the Rayleigh-Plateau instability. The simulation results imply that the string segments will reach zero radius in finite asymptotic time, whence the classical space-time terminates in a naked singularity. Since no fine-tuning is required to excite the instability, this constitutes a generic violation of cosmic censorship.

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

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

Mergers of Magnetized Neutron Stars with Spinning Black Holes: Disruption, Accretion and Fallback

We investigate the merger of a neutron star (of compaction ratio $0.1$) in orbit about a spinning black hole in full general relativity with a mass ratio of $5:1$, allowing for the star to have an initial magnetization of $10^{12} {\rm Gauss}$. We present the resulting gravitational waveform and analyze the fallback accretion as the star is disrupted. The evolutions suggest no significant effects from the initial magnetization. We find that only a negligible amount of matter becomes unbound; $99%$ of the neutron star material has a fallback time of 10 seconds or shorter to reach the region of the central engine and that $99.99%$ of the star will interact with the central disk and black hole within 3 hours.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Post-merger electromagnetic emissions from disks perturbed by binary black holes

We simulate the possible emission from a disk perturbed by a recoiling super-massive black hole. To this end, we study radiation transfer from the system incorporating bremsstrahlung emission from a Maxwellian plasma and absorption given by Kramer's opacity law modified to incorporate blackbody effects. We employ this model in the radiation transfer integration to compute the luminosity at several frequencies, and compare with previous bremsstrahlung luminosity estimations from a transparent limit (in which the emissivity is integrated over the computational domain and over all frequencies) and with a simple thermal emission model. We find close agreement between the radiation transfer results and the estimated bremsstrahlung luminosity from previous work for electromagnetic signals above $10^{14}$ Hz. For lower frequencies, we find a self-eclipsing behavior in the disk, resulting in a strong intensity variability connected to the orbital period of the disk.

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
preprint2005arXiv

Geometrically motivated hyperbolic coordinate conditions for numerical relativity: Analysis, issues and implementations

We study the implications of adopting hyperbolic driver coordinate conditions motivated by geometrical considerations. In particular, conditions that minimize the rate of change of the metric variables. We analyze the properties of the resulting system of equations and their effect when implementing excision techniques. We find that commonly used coordinate conditions lead to a characteristic structure at the excision surface where some modes are not of outflow-type with respect to any excision boundary chosen inside the horizon. Thus, boundary conditions are required for these modes. Unfortunately, the specification of these conditions is a delicate issue as the outflow modes involve both gauge and main variables. As an alternative to these driver equations, we examine conditions derived from extremizing a scalar constructed from Killing's equation and present specific numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Recent analytical and numerical techniques applied to the Einstein equations

Combining deeper insight of Einstein's equations with sophisticated numerical techniques promises the ability to construct accurate numerical implementations of these equations. We illustrate this in two examples, the numerical evolution of ``bubble'' and single black hole spacetimes. The former is chosen to demonstrate how accurate numerical solutions can answer open questions and even reveal unexpected phenomena. The latter illustrates some of the difficulties encountered in three-dimensional black hole simulations, and presents some possible remedies.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Numerical Relativity: A review

Computer simulations are enabling researchers to investigate systems which are extremely difficult to handle analytically. In the particular case of General Relativity, numerical models have proved extremely valuable for investigations of strong field scenarios and been crucial to reveal unexpected phenomena. Considerable efforts are being spent to simulate astrophysically relevant simulations, understand different aspects of the theory and even provide insights in the search for a quantum theory of gravity. In the present article I review the present status of the field of Numerical Relativity, describe the techniques most commonly used and discuss open problems and (some) future prospects.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Exact Solutions for the Intrinsic Geometry of Black Hole Coalescence

We describe the null geometry of a multiple black hole event horizon in terms of a conformal rescaling of a flat space null hypersurface. For the prolate spheroidal case, we show that the method reproduces the pair-of-pants shaped horizon found in the numerical simulation of the head-on-collision of black holes. For the oblate case, it reproduces the initially toroidal event horizon found in the numerical simulation of collapse of a rotating cluster. The analytic nature of the approach makes further conclusions possible, such as a bearing on the hoop conjecture. From a time reversed point of view, the approach yields a description of the past event horizon of a fissioning white hole, which can be used as null data for the characteristic evolution of the exterior space-time.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Matching characteristic codes: exploiting two directions

Combining incoming and outgoing characteristic formulations can provide numerical relativists with a natural implementation of Einstein's equations that better exploits the causal properties of the spacetime and gives access to both null infinity and the interior region simultaneously (assuming the foliation is free of caustics and crossovers). We discuss how this combination can be performed and illustrate its behavior in the Einstein-Klein-Gordon field in 1D.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

High-powered Gravitational News

We describe the computation of the Bondi news for gravitational radiation. We have implemented a computer code for this problem. We discuss the theory behind it as well as the results of validation tests. Our approach uses the compactified null cone formalism, with the computational domain extending to future null infinity and with a worldtube as inner boundary. We calculate the appropriate full Einstein equations in computational eth form in (a) the interior of the computational domain and (b) on the inner boundary. At future null infinity, we transform the computed data into standard Bondi coordinates and so are able to express the news in terms of its standard $N_{+}$ and $N_{\times}$ polarization components. The resulting code is stable and second-order convergent. It runs successfully even in the highly nonlinear case, and has been tested with the news as high as 400, which represents a gravitational radiation power of about $10^{13}M_{\odot}/sec$.

0 citations0 reviewsNo code linkedOpen access