With the goal of providing more accurate and realistic estimates of the secular behavior of the mass accretion and drag rates in the "common-envelope" scenario encountered when a black hole or a neutron star moves in the stellar envelope of a red supergiant star, we have carried out the first general relativistic simulations of the accretion flow onto a nonrotating black hole moving supersonically in a medium with regular but different density gradients. The simulations reveal that the supersonic motion always rapidly reaches a stationary state and it produces a shock cone in the downstream part of the flow. In the absence of density gradients we recover the phenomenology already observed in the well-known Bondi--Hoyle--Lyttleton accretion problem, with super-Eddington mass accretion rate and a shock cone whose axis is stably aligned with the direction of motion. However, as the density gradient is made stronger, the accretion rate also increases and the shock cone is progressively and stably dragged toward the direction of motion. With sufficiently large gradients, the shock-cone axis can become orthogonal to the direction, or even move in the upstream region of the flow in the case of the largest density gradient. Together with the phenomenological aspects of the accretion flow, we have also quantified the rates of accretion of mass and momentum onto the black hole. Simple analytic expressions have been found for the rates of accretion of mass, momentum, drag force, and bremsstrahlung luminosity, all of which have been employed in the astrophysical modelling of the secular evolution of a binary system experiencing a common-envelope evolution...

It is well known that the equation of state (EoS) of compact objects like neutron and quark stars is not determined despite there are several sophisticated models to describe it. From the electromagnetic observations, summarized in \cite{Lattimer01}, and the recent observation of gravitational waves from binary neutron star inspiral GW170817 \cite{Abbott2017_etal} and GW190425 \cite{Abbott2019}, it is possible to make an estimation of the range of masses and so constraint the mass of the neutron and quark stars, determining not only the best approximation for the EoS, but which kind of stars we would be observing. In this paper we explore several configurations of neutron stars assuming a simple polytropic equation of state, using a single layer model without crust. In particular, when the EoS depends on the mass rest density, $p=K \rho_{0}^{\Gamma}$, and when it depends on the energy density $p=K \rho^{\Gamma}$, considerable differences in the mass-radius relationships are found. On the other hand, we also explore quark stars models using the MIT bag EoS for different values of the vacuum energy density $B$.

We build equilibrium solutions of magnetised thick discs around a highly spinning Kerr black hole and evolve these initial data up to a final time of about 100 orbital periods. The numerical simulations reported in this paper solve the general relativistic magnetohydrodynamics equations using the BHAC code and are performed in axisymmetry. Our study assumes non-self-gravitating, polytropic, constant angular momentum discs endowed with a purely toroidal magnetic field. In order to build the initial data we consider three approaches, two of which incorporate the magnetic field in a self-consistent way and a third approach in which the magnetic field is included as a perturbation on to an otherwise purely hydrodynamical solution. To test the dependence of the evolution on the initial data, we explore four representative values of the magnetisation parameter spanning from almost hydrodynamical discs to very strongly magnetised tori. The initial data are perturbed to allow for mass and angular momentum accretion on to the black hole. Notable differences are found in the long-term evolutions of the initial data. In particular, our study reveals that highly magnetised discs are unstable, and hence prone to be fully accreted and expelled, unless the magnetic field is incorporated into the initial data in a self-consistent way.

We study the relativistic Bondi-Hoyle-Lyttleton accretion onto a Schwarzschild black hole (BH), which is surrounded by rigid and small, randomly distributed, bodies. These bodies are idealized representations of substructure like stars passing close to the BH, bubbles created by stellar winds or cold clumps.We explore cases where the filling factor of these bodies is small. The flow is assumed to be adiabatic and move supersonically towards the black hole. The interaction with these rigid obstacles transforms ram pressure of the flow into thermal pressure through bow shocks, slowing down the flow and making the accreting gas turbulent. As a consequence, although the flow reaches a statistically-steady state, the accretion rate presents some variability. For a flow Mach number at infinity of 4, a few of these objects (5 - 10) are enough to increase the accretion rate about 50% over the accretion rate without bodies, even though the gas is adiabatic and the filling factor of the obstacles is small.

We model, for the first time, the Bondi-Hoyle accretion of a fluid with velocity gradients onto a Kerr black hole, by numerically solving the fully relativistic hydrodynamics equations. Specifically, we consider a supersonic ideal gas, which has velocity gradients perpendicular to the relative motion. We measure the mass and specific angular accretion rates to illustrate whether the fluid presents unstable patterns or not. The initial parameters, we consider in this work, are the velocity gradient $\epsilon_{v}$, the black hole spin $a$, the asymptotic Mach number ${\cal M}_{\infty}$ and adiabatic index $\Gamma$. We show that the flow accretion reaches a fairly stationary regime, unlike in the Newtonian case, where significant fluctuations of the mass and angular momentum accretion rates are found. On the other hand, we consider a special case where both density and velocity gradients of the fluid are taken into account. The spin of the black hole and the asymptotic Newtonian Mach number, for this case, are $a=0.98$ and ${\cal M}_{\infty}=1$, respectively. A kind of flip-flop behavior is found at the early times; nevertheless, the system also reaches a steady state.

We present a new code designed to solve the equations of classical ideal magneto-hydrodynamics (MHD) in three dimensions, submitted to a constant gravitational field. The purpose of the code centers on the analysis of solar phenomena within the photosphere-corona region. We present 1D and 2D standard tests to demonstrate the quality of the numerical results obtained with our code. As solar tests we present the transverse oscillations of Alfvenic pulses in coronal loops using a 2.5D model, and as 3D tests we present the propagation of impulsively generated MHD-gravity waves and vortices in the solar atmosphere. The code is based on high-resolution shock-capturing methods, uses the HLLE flux formula combined with Minmod, MC and WENO5 reconstructors. The divergence free magnetic field constraint is controlled using the Flux Constrained Transport method.

In this work, we present, for the first time, a numerical study of the Bondi-Hoyle accretion with density gradients in the fully relativistic regime. In this context, we consider accretion onto a Kerr Black Hole (BH) of a supersonic ideal gas, which has density gradients perpendicular to the relative motion. The set of parameters of interest in this study are the Mach number, M, the spin of the BH, a, and the density-gradient parameter of the gas, {\rho} . We show that, unlike in the Newtonian case, all the studied cases, especially those with density gradient, approach a stationary flow pattern. To illustrate that the system reaches steady state we calculate the mass and angular momentum accretion rates on a spherical surface located almost at the event horizon. In the particular case of M = 1, {\rho} = 0.5 and BH spin a = 0.5, we observe a disk-like configuration surrounding the BH. Finally, we present the gas morphology and some of its properties.

We introduce CAFE, a new independent code designed to solve the equations of Relativistic ideal Magnetohydrodynamics (RMHD) in 3D. We present the standard tests for a RMHD code and for the Relativistic Hydrodynamics (RHD) regime since we have not reported them before. The tests include the 1D Riemann problems related to blast waves, head-on collision of streams and states with transverse velocities, with and without magnetic field, which is aligned or transverse, constant or discontinuous across the initial discontinuity. Among the 2D and 3D tests, without magnetic field we include the 2D Riemann problem, a one dimensional shock tube along a diagonal, the high speed Emery wind tunnel, the Kelvin-Helmholtz instability, a set of jets and a 3D spherical blast wave, whereas in the presence of a magnetic field we show the magnetic rotor, the cylindrical explosion, a case of Kelvin-Helmholtz instability and a 3D magnetic field advection loop. The code uses High Resolution Shock Capturing methods and we present the error analysis for a combination that uses the HLLE flux formula combined with linear, PPM and fifth order WENO reconstructors. We use the flux-CT and the divergence cleaning methods to control the divergence free magnetic field constraint.

We model the radial accretion of radiation on Primordial Black Holes (PBH) by numerically solving Einstein's equations coupled to an ultrarelativistic ideal gas with equation of state $p=\rho/3$. We calculate the final mass of a black hole by the integration of the accreted radiation energy density during the leptonic era between $t\sim10^{-4}s$ to $t\sim 10^2s$ after the Big Bang. Our results indicate that small PBHs with initial masses between $10^{-4}$ to $1M_{\odot}$ may grow up to hundreds of solar masses, and thus can be SMBH seeds. On the other hand, PBHs formed at $t\sim 1s$ with initial mass between 900 and $\sim 980M_{\odot}$, by the time $t\sim 100s$ show masses of $10^4$ to $10^6M_{\odot}$ which are masses of seeds or already formed SMBHs. The fact that we consider only radial flow implies that our results work well as limiting cases, and it is expected that under more general scenarios the accretion rates may change significantly. Nevertheless we show that it is possible that SMBHs can be PBHs that grew due to the accretion of radiation.

In this work, we track the evolution of an ultrarelativistic fluid onto a Kerr black hole, on the equatorial plane. In this treatment, we consider the limit where the rest mass density is neglected, that is, the approximation is valid in the regime where the internal energy dominates over the rest mass density. We particularly concentrate in the case of a gas with $\Gamma$ = 4/3, which corresponds to a radiation fluid. We show, as in several cases, that a shock cone appears when the asymptotic velocity of the fluid is larger than the asymptotic relativistic sound speed of the gas. On the other hand, in order to show the system approaches to steady state, we calculate the accreted total energy rate on a spherical surface. Finally, we also show the gas distribution and various of its properties.

We study numerically the relativistic Bondi-Hoyle accretion of an ideal gas onto a Kerr fixed background space-time on the equatorial plane with s-lab symmetry. We use both Kerr-Schild (KS) and Boyer-Lindquist (BL) coordinates. We particularly focus on the study of the flip-flop motion of the shock cone formed when the gas is injected at supersonic speed. The development of the flip-flop instability of the shock cone in the relativistic regime was reported recently for the first time. We reproduce the flip-flop behaviour found in the past when BL coordinates are used, and perform similar numerical experiments using horizon penetrating KS coordinates. We find that when using KS coordinates the shock cone oscillates, however such oscillations are not of the flip-flop type and their amplitude decrease with resolution.

We present the numerical evolution of a massive test scalar fields around a Schwarzschild space-time. We proceed by using hyperboloidal slices that approach future null infinity, which is the boundary of scalar fields, and also demand the slices to penetrate the event horizon of the black hole. This approach allows the scalar field to be accreted by the black hole and to escape toward future null infinity. We track the evolution of the energy density of the scalar field, which determines the rate at which the scalar field is being diluted. We find polynomial decay of the energy density of the scalar field, and use it to estimate the rate of dilution of the field in time. Our findings imply that the energy density of the scalar field decreases even five orders of magnitude in time scales smaller than a year. This implies that if a supermassive black hole is the Schwarzschild solution, then scalar field dark matter would be diluted extremely fast

We numerically solve the mass-less test scalar field equation on the space-time background of boson stars and black holes. In order to do so, we use a numerical domain that contains future null infinity. We achieve this construction using a scri-fixing conformal compactification technique based on hyperboloidal constant mean curvature foliations of the space-time and solve the conformally invariant wave equation. We present two results: the scalar field shows oscillations of the quasi- normal-mode type found for black holes only for boson star configurations that are compact, and no signs of tail decay is found in the parameter space we explored. Even though our results do not correspond to the master equation of perturbations of boson star solutions, they indicate that the parameter space of boson stars as black hole mimickers is restricted to compact configurations.