Source author record

Miguel Zilhao

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

14 published item(s)

preprint2022arXiv

Gauge structure of the Einstein field equations in Bondi-like coordinates

The characteristic initial (boundary) value problem has numerous applications in general relativity (GR) involving numerical studies, and is often formulated using Bondi-like coordinates. Recently it was shown that several prototype formulations of this type are only weakly hyperbolic. Presently we examine the root cause of this result. In a linear analysis we identify the gauge, constraint and physical blocks in the principal part of the Einstein field equations in such a gauge, and show that the subsystem related to the gauge variables is only weakly hyperbolic. Weak hyperbolicity of the full system follows as a consequence in many cases. We demonstrate this explicitly in specific examples, and thus argue that Bondi-like gauges result in weakly hyperbolic free evolution systems under quite general conditions. Consequently the characteristic initial (boundary) value problem of GR in these gauges is rendered ill-posed in the simplest norms one would like to employ. The possibility of finding good alternative norms, in which well-posedness is achieved, is discussed. So motivated, we present numerical convergence tests with an implementation of full GR which demonstrate the effect of weak hyperbolicity in practice.

preprint2022arXiv

New Horizons for Fundamental Physics with LISA

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas.

preprint2021arXiv

Black hole binaries and light fields: Gravitational molecules

We show that light scalars can form quasibound states around binaries. In the nonrelativistic regime, these states are formally described by the quantum-mechanical Schrodinger equation for a one-electron heteronuclear diatomic molecule. We performed extensive numerical simulations of scalar fields around black hole binaries showing that a scalar structure condenses around the binary -- we dub these states "gravitational molecules". We further show that these are well described by the perturbative, nonrelativistic description.

preprint2020arXiv

Collective scalarization or tachyonization: when averaging fails

Certain scalar-tensor theories of gravity provide negative-energy, tachyonic modes to a fundamental scalar inside matter, giving rise to non-perturbative phenomena around compact stars. Studies of this and other tachyonic instabilities always average over local matter properties. We use elementary, flat space models to understand possible collective effects and the accuracy of the averaging procedure. In particular, we consider bodies made of elementary constituents which do not, in isolation, scalarize because their compactness ${\cal C}$ is too small, ${\cal C}\lesssim {\cal C}_{\rm crit}$. We show that when the individual constituents have compactness smaller but close to the threshold, one is able to scalarize composite bodies through collective effects, and the compactness of the composite body can be made arbitrarily small. On the other hand, our results suggest that when the fundamental building blocks have very low compactness, then scalarization of the composite body requires a global compactness ${\cal C}_{\rm global}\gtrsim {\cal C}_{\rm crit}$. Thus, our results rule out scalarization of dilute bodies via collective effects.

preprint2020arXiv

Crossing a large-$N$ phase transition at finite volume

The existence of phase-separated states is an essential feature of infinite-volume systems with a thermal, first-order phase transition. At energies between those at which the phase transition takes place, equilibrium homogeneous states are either metastable or suffer from a spinodal instability. In this range the stable states are inhomogeneous, phase-separated states. We use holography to investigate how this picture is modified at finite volume in a strongly coupled, four-dimensional gauge theory. We work in the planar limit, $N \to \infty$, which ensures that we remain in the thermodynamic limit. We uncover a rich set of inhomogeneous states dual to lumpy black branes on the gravity side, as well as first- and second-order phase transitions between them. We establish their local (in)stability properties and show that fully non-linear time evolution in the bulk takes unstable states to stable ones.

preprint2020arXiv

Dynamics of Phase Separation from Holography

We use holography to develop a physical picture of the real-time evolution of the spinodal instability of a four-dimensional, strongly-coupled gauge theory with a first-order, thermal phase transition. We numerically solve Einstein's equations to follow the evolution, in which we identify four generic stages: A first, linear stage in which the instability grows exponentially; a second, non-linear stage in which peaks and/or phase domains are formed; a third stage in which these structures merge; and a fourth stage in which the system finally relaxes to a static, phase-separated configuration. On the gravity side the latter is described by a static, stable, inhomogeneous horizon. We conjecture and provide evidence that all static, non-phase separated configurations in large enough boxes are dynamically unstable. We show that all four stages are well described by the constitutive relations of second-order hydrodynamics that include all second-order gradients that are purely spatial in the local rest frame. In contrast, a Müller-Israel-Stewart-type formulation of hydrodynamics fails to provide a good description for two reasons. First, it misses some large, purely-spatial gradient corrections. Second, several second-order transport coefficients in this formulation, including the relaxation times $τ_π$ and $τ_Π$, diverge at the points where the speed of sound vanishes.

preprint2016arXiv

Inspiraling black-hole binary spacetimes: Challenges in transitioning from analytical to numerical techniques

We explore how a recently developed analytical black-hole binary spacetime can be extended using numerical simulations to go beyond the slow-inspiral phase. The analytic spacetime solves the Einstein field equations approximately, with the approximation error becoming progressively smaller the more separated the binary. To continue the spacetime beyond the slow-inspiral phase, we need to transition. Such a transition was previously explored at smaller separations. Here, we perform this transition at a separation of D=20M (large enough that the analytical metric is expected to be accurate), and evolve for six orbits. We find that small constraint violations can have large dynamical effects, but these can be removed by using a constraint-damping system like the conformal covariant formulation of the Z4 system. We find agreement between the subsequent numerical spacetime and the predictions of post-Newtonian theory for the waveform and inspiral rate that is within the post-Newtonian truncation error.

preprint2015arXiv

Testing General Relativity with Present and Future Astrophysical Observations

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

preprint2014arXiv

Higher dimensional Numerical Relativity: code comparison

The nonlinear behavior of higher dimensional black hole spacetimes is of interest in several contexts, ranging from an understanding of cosmic censorship to black hole production in high-energy collisions. However, nonlinear numerical evolutions of higher dimensional black hole spacetimes are tremendously complex, involving different diagnostic tools and "dimensional reduction methods". In this work we compare two different successful codes to evolve Einstein's equations in higher dimensions, and show that the results of such different procedures agree to numerical precision, when applied to the collision from rest of two equal-mass black holes. We calculate the total radiated energy to be E/M=9x10^{-4} in five dimensions and E/M=8.1x10^{-4} in six dimensions.

preprint2012arXiv

Dynamics of black holes in de Sitter spacetimes

Nonlinear dynamics in cosmological backgrounds has the potential to teach us immensely about our universe, and also to serve as prototype for nonlinear processes in generic curved spacetimes. Here we report on dynamical evolutions of black holes in asymptotically de Sitter spacetimes. We focus on the head-on collision of equal mass binaries and for the first time compare analytical and perturbative methods with full blown nonlinear simulations. Our results include an accurate determination of the merger/scatter transition (consequence of an expanding background) for small mass binaries and a test of the Cosmic Censorship conjecture, for large mass binaries. We observe that, even starting from small separations, black holes in large mass binaries eventually lose causal contact, in agreement with the conjecture.

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.

preprint2010arXiv

Black holes in a box: towards the numerical evolution of black holes in AdS

The evolution of black holes in "confining boxes" is interesting for a number of reasons, particularly because it mimics the global structure of Anti-de Sitter geometries. These are non-globally hyperbolic space-times and the Cauchy problem may only be well defined if the initial data is supplemented by boundary conditions at the time-like conformal boundary. Here, we explore the active role that boundary conditions play in the evolution of a bulk black hole system, by imprisoning a black hole binary in a box with mirror-like boundary conditions. We are able to follow the post-merger dynamics for up to two reflections off the boundary of the gravitational radiation produced in the merger. We estimate that about 15% of the radiation energy is absorbed by the black hole per interaction, whereas transfer of angular momentum from the radiation to the black hole is only observed in the first interaction. We discuss the possible role of superradiant scattering for this result. Unlike the studies with outgoing boundary conditions, both the Newman-Penrose scalars Ψ_4 and Ψ_0 are non-trivial in our setup, and we show that the numerical data verifies the expected relations between them.

preprint2010arXiv

Numerical relativity for D dimensional axially symmetric space-times: formalism and code tests

The numerical evolution of Einstein's field equations in a generic background has the potential to answer a variety of important questions in physics: from applications to the gauge-gravity duality, to modelling black hole production in TeV gravity scenarios, analysis of the stability of exact solutions and tests of Cosmic Censorship. In order to investigate these questions, we extend numerical relativity to more general space-times than those investigated hitherto, by developing a framework to study the numerical evolution of D dimensional vacuum space-times with an SO(D-2) isometry group for D\ge 5, or SO(D-3) for D\ge 6. Performing a dimensional reduction on a (D-4)-sphere, the D dimensional vacuum Einstein equations are rewritten as a 3+1 dimensional system with source terms, and presented in the Baumgarte, Shapiro, Shibata and Nakamura (BSSN) formulation. This allows the use of existing 3+1 dimensional numerical codes with small adaptations. Brill-Lindquist initial data are constructed in D dimensions and a procedure to match them to our 3+1 dimensional evolution equations is given. We have implemented our framework by adapting the LEAN code and perform a variety of simulations of non-spinning black hole space-times. Specifically, we present a modified moving puncture gauge which facilitates long term stable simulations in D=5. We further demonstrate the internal consistency of the code by studying convergence and comparing numerical versus analytic results in the case of geodesic slicing for D=5,6.

preprint2010arXiv

Numerical relativity for D dimensional space-times: head-on collisions of black holes and gravitational wave extraction

Black objects in higher dimensional space-times have a remarkably richer structure than their four dimensional counterparts. They appear in a variety of configurations (e.g. black holes, black branes, black rings, black Saturns), and display complex stability phase diagrams. They might also play a key role in high energy physics: for energies above the fundamental Planck scale, gravity is the dominant interaction which, together with the hoop-conjecture, implies that the trans-Planckian scattering of point particles should be well described by black hole scattering. Higher dimensional scenarios with a fundamental Planck scale of the order of TeV predict, therefore, black hole production at the LHC, as well as in future colliders with yet higher energies. In this setting, accurate predictions for the production cross-section and energy loss (through gravitational radiation) in the formation of black holes in parton-parton collisions is crucial for accurate phenomenological modelling in Monte Carlo event generators. In this paper, we use the formalism and numerical code reported in arXiv:1001.2302 to study the head-on collision of two black holes. For this purpose we provide a detailed treatment of gravitational wave extraction in generic D-dimensional space-times, which uses the Kodama-Ishibashi formalism. For the first time, we present the results of numerical simulations of the head-on collision in five space-time dimensions, together with the relevant physical quantities. We show that the total radiated energy, when two black holes collide from rest at infinity, is approximately (0.089\pm 0.006)% of the centre of mass energy, slightly larger than the 0.055% obtained in the four dimensional case, and that the ringdown signal at late time is in very good agreement with perturbative calculations.