We perform numerical evolutions of the fully non-linear Einstein-(complex, massive)Klein-Gordon and Einstein-(complex)Proca systems, to assess the formation and stability of spinning bosonic stars. In the scalar/vector case these are known as boson/Proca stars. Firstly, we consider the formation scenario. Starting with constraint-obeying initial data, describing a dilute, axisymmetric cloud of spinning scalar/Proca field, gravitational collapse towards a spinning star occurs, via gravitational cooling. In the scalar case the formation is transient, even for a non-perturbed initial cloud; a non-axisymmetric instability always develops ejecting all the angular momentum from the scalar star. In the Proca case, by contrast, no instability is observed and the evolutions are compatible with the formation of a spinning Proca star. Secondly, we address the stability of an existing star, a stationary solution of the field equations. In the scalar case, a non-axisymmetric perturbation develops collapsing the star to a spinning black hole. No such instability is found in the Proca case, where the star survives large amplitude perturbations; moreover, some excited Proca stars decay to, and remain as, fundamental states. Our analysis suggests bosonic stars have different stability properties in the scalar/vector case, which we tentatively relate to their toroidal/spheroidal morphology. A parallelism with instabilities of spinning fluid stars is briefly discussed.

The evolution of planets, stars and even galaxies is driven, to a large extent, by dynamical friction of gravitational origin. There is now a good understanding of the friction produced by extended media, either collisionless of fluid-like. However, the physics of accretion or protoplanetary disks, for instance, is described by slab-like geometries instead, compact in one spatial direction. Here, we find, for the first time, the gravitational wake due to a massive perturber moving through a slab-like medium, describing e.g. accretion disks with sharp transitions. We show that dynamical friction in such environments can be substantially reduced relatively to spatially extended profiles. Finally, we provide simple and accurate expressions for the gravitational drag force felt by the perturber, in both the subsonic and supersonic regime.

Spherically symmetric bosonic stars are one of the few examples of gravitating solitons that are known to form dynamically, via a classical process of (incomplete) gravitational collapse. As stationary solutions of the Einstein--Klein-Gordon or the Einstein--Proca theory, bosonic stars may also become sufficiently compact to develop light rings and hence mimic, in principle, gravitational-wave observational signatures of black holes (BHs). In this paper, we discuss how these horizonless ultra-compact objects (UCOs) are actually distinct from BHs, both phenomenologically and dynamically. In the electromagnetic channel, the light ring associated phenomenology reveals remarkable lensing patterns, quite distinct from a standard BH shadow, with an infinite number of Einstein rings accumulating in the vicinity of the light ring, both inside and outside the latter. The strong lensing region, moreover, can be considerably smaller than the shadow of a BH with a comparable mass. Dynamically, we investigate the fate of such UCOs under perturbations, via fully non-linear numerical simulations and observe that, in all cases, they decay into a Schwarzschild BH within a time scale of $\mathcal{O}(M)$, where $M$ is the mass of the bosonic star. Both these studies reinforce how difficult it is for horizonless UCOs to mimic BH phenomenology and dynamics, in all its aspects.

In this paper we evolve magnetized and unmagnetized circumbinary accretion disks around supermassive black hole binaries in the relativistic regime. We use a post-Newtonian expansion to construct an analytical spacetime and determine how the order of the post-Newtonian (PN) expansion affects the dynamics of the gas. We find very small differences in the late-time bulk dynamics of non-magnetized hydrodynamic evolutions between the two spacetimes down to separations of approximately $40GM/c^2$ where $M$ is the total mass of the binary. For smaller separations, the differences due to PN-order become comparable to differences caused by using initial data further from equilibrium. For magnetized gas, MHD stresses, which drives the accretion dynamics, tends to mask all higher order PN effects even at separations of $20GM/c^2$, leading to essentially the same observed electromagnetic luminosity. This implies that our calculations of the EM signal may be robust down to small binary separations. Our investigation is the first to demonstrate how the level of PN accuracy affects a circumbinary disk's evolution and informs us of the range in separation within which to trust the PN approximation for this kind of study. We also address the influence the initial conditions and binary separation have on simulation predictions.

We present a new warped gridding scheme adapted to simulating gas dynamics in binary black hole spacetimes. The grid concentrates grid points in the vicinity of each black hole to resolve the smaller scale structures there, and rarefies grid points away from each black hole to keep the overall problem size at a practical level. In this respect, our system can be thought of as a "double" version of the fisheye coordinate system, used before in numerical relativity codes for evolving binary black holes. The gridding scheme is constructed as a mapping between a uniform coordinate system--in which the equations of motion are solved--to the distorted system representing the spatial locations of our grid points. Since we are motivated to eventually use this system for circumbinary disc calculations, we demonstrate how the distorted system can be constructed to asymptote to the typical spherical polar coordinate system, amenable to efficiently simulating orbiting gas flows about central objects with little numerical diffusion. We discuss its implementation in the Harm3d code, tailored to evolve the magnetohydrodynamics equations in curved spacetimes. We evaluate the performance of the system's implementation in Harm3d with a series of tests, such as the advected magnetic field loop test, magnetized Bondi accretion, and evolutions of hydrodynamic discs about a single black hole and about a binary black hole. Like we have done with Harm3d, this gridding scheme can be implemented in other unigrid codes as a (possibly) simpler alternative to adaptive mesh refinement.

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.