We present new, fully nonlinear numerical solutions to the static, spherically symmetric Einstein-Klein-Gordon system for a collection of an arbitrary odd number $N$ of complex scalar fields with an internal $U(N)$ symmetry and no self-interactions. These solutions, which we dub $\ell$-boson stars, are parametrized by an angular momentum number $\ell=(N-1)/2$, an excitation number $n$, and a continuous parameter representing the amplitude of the fields. They are regular at every point and possess a finite total mass. For $\ell = 0$ the standard spherically symmetric boson stars are recovered. We determine their generalizations for $\ell > 0$, and show that they give rise to a large class of new static configurations which might have a much larger compactness ratio than $\ell=0$ stars.

It has long been known that no static, spherically symmetric, asymptotically flat Klein-Gordon scalar field configuration surrounding a nonrotating black hole can exist in general relativity. In a series of previous papers we proved that, at the effective level, this no-hair theorem can be circumvented by relaxing the staticity assumption: for appropriate model parameters there are quasi-bound scalar field configurations living on a fixed Schwarzschild background which, although not being strictly static, have a larger lifetime than the age of the universe. This situation arises when the mass of the scalar field distribution is much smaller than the black hole mass, and following the analogies with the hair in the literature we dubbed these long-lived field configurations wigs. Here we extend our previous work to include the gravitational backreaction produced by the scalar wigs. We derive new approximate solutions of the spherically symmetric Einstein-Klein-Gordon system which represent self-gravitating scalar wigs surrounding black holes. These configurations interpolate between boson star configurations and Schwarzschild black holes dressed with the long-lived scalar test field distributions discussed in previous papers. Nonlinear numerical evolutions of initial data sets extracted from our approximate solutions support the validity of our approach. Arbitrarily large lifetimes are still possible, although for the parameter space that we analyze in this paper they seem to decay faster than the quasi-bound states. Finally, we speculate about the possibility that these configurations could describe the innermost regions of dark matter halos.

We study the evolution of a massive scalar field surrounding a Schwarzschild black hole and find configurations that can survive for arbitrarily long times, provided the black hole or the scalar field mass is small enough. In particular, both ultra-light scalar field dark matter around supermassive black holes and axion-like scalar fields around primordial black holes can survive for cosmological times. Moreover, these results are quite generic, in the sense that fairly arbitrary initial data evolves, at late times, as a combination of those long-lived configurations.

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.

The merger process of a binary black hole system can have a strong impact on a circumbinary disk. In the present work we study the effect of both central mass reduction (due to the energy loss through gravitational waves) and a possible black hole recoil (due to asymmetric emission of gravitational radiation). For the mass reduction case and recoil directed along the disk's angular momentum, oscillations are induced in the disk which then modulate the internal energy and bremsstrahlung luminosities. On the other hand, when the recoil direction has a component orthogonal to the disk's angular momentum, the disk's dynamics are strongly impacted, giving rise to relativistic shocks. The shock heating leaves its signature in our proxies for radiation, the total internal energy and bremsstrahlung luminosity. Interestingly, for cases where the kick velocity is below the smallest orbital velocity in the disk (a likely scenario in real AGN), we observe a common, characteristic pattern in the internal energy of the disk. Variations in kick velocity simply provide a phase offset in the characteristic pattern implying that observations of such a signature could yield a measure of the kick velocity through electromagnetic signals alone.