Search

A line of first-order phase transitions is conjectured in the phase diagram of Quantum Chromodynamics at non-zero baryon density. If this is the case, numerical simulations of neutron star mergers suggest that various regions of the stars may cross this line multiple times. This results in the nucleation of bubbles of the preferred phase, which subsequently expand and collide. The resulting gravitational wave spectrum is highly sensitively to the velocity of the bubble walls. We use holography to perform the first microscopic simulation of bubble dynamics in a theory that qualitatively mirrors the expected phase diagram of Quantum Chromodynamics. We determine the wall velocity in the metastable regions and we compare it to theoretical estimates. We discuss implications for gravitational wave production., Comment: 31 pages, 21 figures

Recently, the piercing of a mini boson star by a black hole was studied, with tidal capture and the discovery of a "gravitational atom" being reported ( arXiv:2206.00021 [gr-qc] ). Building on this research, we extend the study by including a hexic solitonic potential and explore the piercing of a solitonic boson star by a black hole. Notably, the solitonic boson star can reach higher compactness, which one might expect could alter the dynamics in this context. Our findings suggest that even when the black hole's size approaches the test particle limit, the solitonic boson star is easily captured by the black hole due to an extreme tidal capture process. Regardless of the black hole initial mass and velocity, our results indicate that over 85% of the boson star material is accreted. Thus, the self-interaction does not alter the qualitative behavior of the system., Comment: 11 pages, 11 figures, revised to match the published version

Fuzzy dark matter is an exciting alternative to the standard cold dark matter paradigm, reproducing its large scale predictions, while solving most of the existing tension with small scale observations. These models postulate that dark matter is constituted by light bosons and predict the condensation of a solitonic core -- also known as boson star, supported by wave pressure -- at the center of halos. However, solitons which host a \emph{parasitic} supermassive black hole are doomed to be swallowed by their guest. It is thus crucial to understand in detail the accretion process. In this work, we use numerical relativity to self-consistently solve the problem of accretion of a boson star by a central black hole, in spherical symmetry. We identify three stages in the process, a {\it boson-quake}, a {\it catastrophic stage} and a linear phase, as well as a general accurate expression for the lifetime of a boson star with an endoparasitic black hole. Lifetimes of these objects can be large enough to allow them to survive until the present time., Comment: v2. Overall improvement, added Supplemental Material. To appear in Physical Review D Letters. Movies are also available on https://centra.tecnico.ulisboa.pt/network/grit/files/movies/

We explore the dynamics of boson stars in the presence of axionic couplings through nonlinear evolutions of Einstein's field equations. We show that, for large axionic couplings, isolated boson stars become unstable, and decay via a large burst of electromagnetic radiation, becoming less massive and more dilute. Our full nonlinear results are in good agreement with flat-space estimates for the critical couplings. We then consider head-on collisions of sub-critical boson stars and study the electromagnetic and gravitational signal. Boson stars cluster around the critical point via interactions, and we argue that mergers will generically be also sources of electromagnetic radiation, in addition to gravitational waves, which can be used to place constraints on the axionic coupling if such multimessenger signals are detected., Comment: 12 pages, 10 figures

New light fundamental fields are natural candidates for all or a fraction of dark matter. Self-gravitating structures of such fields might be common objects in the universe, and could comprise even galactic halos. These structures would interact gravitationally with black holes, a process of the utmost importance since it dictates their lifetime, the black hole motion, and possible gravitational radiation emission. Here, we study the dynamics of a black hole piercing through a much larger fully relativistic boson star, made of a complex minimally coupled massive scalar without self-interactions. As the black hole pierces through the bosonic structure, it is slowed down by accretion and dynamical friction, giving rise to gravitational-wave emission. Since we are interested in studying the interaction with large and heavy scalar structures, we consider mass ratios up to $q\sim 10$ and length ratios ${\cal L} \sim 62$. Somewhat surprisingly, for all our simulations, the black hole accretes more than 95% of the boson star material, even if an initially small black hole collides with large velocity. This is a consequence of an extreme "tidal capture" process, which binds the black hole and the boson star together, for these mass ratios. We find evidence of a "gravitational atom" left behind as a product of the process., Comment: 13 pages, 11 figures, 2 movies in ancillary files. revised to match the published version. Movies are also available on https://centra.tecnico.ulisboa.pt/network/grit/files/movies/

Cosmological phase transitions can proceed via the nucleation of bubbles that subsequently expand and collide. The resulting gravitational wave spectrum depends crucially on the properties of these bubbles. We extend our previous holographic work on planar bubbles to circular bubbles in a strongly-coupled, non-Abelian, four-dimensional gauge theory. This extension brings about two new physical properties. First, the existence of a critical bubble, which we determine. Second, the bubble profile at late times exhibits a richer self-similar structure, which we verify. These results require a new 3+1 evolution code called Jecco that solves the Einstein equations in the characteristic formulation in asymptotically AdS spaces. Jecco is written in the Julia programming language and is freely available. We present an outline of the code and the tests performed to assess its robustness and performance., Comment: 63 pages, 23 figures. minor typo fixed

We uncover a new gravitational-wave production mechanism in cosmological, first-order, thermal phase transitions. These are usually assumed to proceed via the nucleation of bubbles of the stable phase inside the metastable phase. However, if the nucleation rate is sufficiently suppressed, then the Universe may supercool all the way down the metastable branch and enter the spinodal region. In this case the transition proceeds via the exponential growth of unstable modes and the subsequent formation, merging and relaxation of phase domains. We use holography to follow the real-time evolution of this process in a strongly coupled, four-dimensional gauge theory. The resulting gravitational wave spectrum differs qualitatively from that in transitions mediated by bubble nucleation. We discuss the possibility that the spinodal dynamics may be preceded by a period of thermal inflation., Comment: 45 pages, 19 figures

We use holography to study collisions of phase domains formed in a four-dimensional, strongly-coupled gauge theory with a first-order, thermal phase transition. We find three qualitatively different dynamical regimes depending on the collision velocity. For low velocities the domains slow down before the collision and subsequently merge and relax to equilibrium. For intermediate velocities no slow down is present before the merger. For high enough velocities the domains can collide and break apart several times before they finally merge. These features leave an imprint on the time evolution of the entropy of the system, which we compute from the area of the dual horizon on the gravity side., Comment: 24 pages, 23 figures

Cosmological phase transitions proceed via the nucleation of bubbles that subsequently expand and collide. The resulting gravitational wave spectrum depends crucially on the bubble wall velocity. Microscopic calculations of this velocity are challenging even in weakly coupled theories. We use holography to compute the wall velocity from first principles in a strongly coupled, non-Abelian, four-dimensional gauge theory. The wall velocity is determined dynamically in terms of the nucleation temperature. We find an approximately linear relation between the velocity and the ratio $\Delta \mathcal{P}/\mathcal{E}$, with $\Delta \mathcal{P}$ the pressure difference between the inside and the outside of the bubble and $\mathcal{E}$ the energy density outside the bubble. Up to a rescaling, the wall profile is well approximated by that of an equilibrium, phase-separated configuration at the critical temperature. We verify that ideal hydrodynamics provides a good description of the system everywhere except near the wall., Comment: 6 pages, 9 figures. Matches published version in PRD

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., Comment: 16 pages, 12 figures. Matches published version in PRD

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., Comment: 6 pages, version to appear in Physical Review Letters

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\"uller-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 $\tau_\pi$ and $\tau_\Pi$, diverge at the points where the speed of sound vanishes., Comment: 45 pages, 35 figures; published version

We use holography to analyse relativistic collisions in a one-parameter family of strongly coupled gauge theories with thermal phase transitions. For a critical value of the parameter the transition is second order, for subcritical values it is first order, and for super-critical values it is a smooth crossover. We extract the gauge theory stress tensor from collisions of gravitational shock waves on the dual geometries. Regardless of the nature of the transition, for values of the parameter close to the critical value almost all the energy of the projectiles is deposited into a long-lived, quasi-static blob of energy at mid-rapidity. This configuration is 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\"uller-Israel-Stewart-type formulation of hydrodynamics fails to provide a good description. We discuss possible implications for searches of the QCD critical point., Comment: 6 pages, 3 figures

We extend our previous analysis of holographic heavy ion collisions in non-conformal theories. We provide a detailed description of our numerical code. We study collisions at different energies in gauge theories with different degrees of non-conformality. We compare four relaxation times: the hydrodynamization time (when hydrodynamics becomes applicable), the EoSization time (when the average pressure approaches its equilibrium value), the isotropization time (when the longitudinal and transverse pressures approach each other) and the condensate relaxation time (when the expectation value of a scalar operator approaches its equilibrium value). We find that these processes can occur in several different orderings. In particular, the condensate can remain far from equilibrium even long after the plasma has hydrodynamized and EoSized. We also explore the rapidity distribution of the energy density at hydrodynamization. This is far from boost-invariant and its width decreases as the non-conformality increases. Nevertheless, the velocity field at hydrodynamization is almost exactly boost-invariant regardless of the non-conformality. This result may be used to constrain the initialization of hydrodynamic fields in heavy ion collisions., Comment: 45 pages, 17 figures. Expanded some sections. Matches published version in JHEP

We use holography to study the spinodal instability of a four-dimensional, strongly-coupled gauge theory with a first-order thermal phase transition. We place the theory on a cylinder in a set of homogeneous, unstable initial states. The dual gravity configurations are black branes afflicted by a Gregory-Laflamme instability. We numerically evolve Einstein's equations to follow the instability until the system settles down to a stationary, inhomogeneous black brane. The dual gauge theory states have constant temperature but non-constant energy density. We show that the time evolution of the instability and the final states are accurately described by second-order hydrodynamics. In the static limit, the latter reduces to a single, second-order, non-linear differential equation from which the inhomogeneous final states can be derived., Comment: 11 pages, 6 figures. Matches published version in JHEP

We numerically simulate gravitational shock wave collisions in a holographic model dual to a non-conformal four-dimensional gauge theory. We find two novel effects associated to the non-zero bulk viscosity of the resulting plasma. First, the hydrodynamization time increases. Second, if the bulk viscosity is large enough then the plasma becomes well described by hydrodynamics before the energy density and the average pressure begin to obey the equilibrium equation of state. We discuss implications for the quark-gluon plasma created in heavy ion collision experiments., Comment: 10 pages, 3 figures; published version

We study the equilibrium and near-equilibrium properties of a holographic five-dimensional model consisting of Einstein gravity coupled to a scalar field with a non-trivial potential. The dual four-dimensional gauge theory is not conformal and, at zero temperature, exhibits a renormalisation group flow between two different fixed points. We quantify the deviations from conformality both in terms of thermodynamic observables and in terms of the bulk viscosity of the theory. The ratio of bulk over shear viscosity violates Buchel's bound. We study relaxation of small-amplitude, homogeneous perturbations by computing the quasi-normal modes of the system at zero spatial momentum. In this approximation we identify two different relaxation channels. At high temperatures, the different pressures first become approximately equal to one another, and subsequently this average pressure evolves towards the equilibrium value dictated by the equation of state. At low temperatures, the average pressure first evolves towards the equilibrium pressure, and only later the different pressures become approximately equal to one another., Comment: 22 pages, 29 figures; published version

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., Comment: 188 pages, 46 figures, 6 tables, 903 references. Matches version published in Classical and Quantum Gravity. Supplementary data files available at http://www.phy.olemiss.edu/~berti/research/ and http://centra.tecnico.ulisboa.pt/network/grit/files/

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., Comment: 7 pages, RevTex4

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., Comment: 68 pages, 4 Figures

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., Comment: 31 pages, 12 figures, RevTex4. v2: Published version. Further details can be found at http://blackholes.ist.utl.pt/

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 \Psi_4 and \Psi_0 are non-trivial in our setup, and we show that the numerical data verifies the expected relations between them., Comment: REvTex4, 17 pages, 12 Figs. v2: Minor improvements. Published version. Animation of a black hole binary in a box can be found at http://blackholes.ist.utl.pt/

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., Comment: 31 pages, 6 figures; v2 Minor changes and added two references. Matches the published version in PRD.