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$.

The origin of the magnetic field in some systems like AGNs or GRBs is still an open question in astrophysics. A possible mechanism to explain the magnetic field generation is the Kelvin-Helmholtz instability, since it is able to transform the kinetic energy, in a shear flow, into magnetic energy. Through the present work, we investigate the linear and non linear effects produced by the magnetic susceptibility in the development of the Kelvin-Helmholtz instability in a relativistic plasma. The magnetic field in the system is parallel to the flows and the susceptibility is assumed to be homogeneous, constant in time, and equal in both fluids. In particular, we analyze the instability in three different cases, when the fluids are diamagnetic, paramagnetic, and when the susceptibility is zero. We compute the dispersion relation in the linear regime and found that the interface between diamagnetic fluids is more stable than the original case without magnetic susceptibility. On the contrary, the paramagnetism makes the interface more unstable. This behavior is increasingly important the higher the magnetization parameter. Additionally, the interval of relativistic Mach numbers for which the interface is unstable grows with the magnetic susceptibility, i.e., it is larger for paramagnetic fluids. We check these analytical results with numerical simulations, and explore the effect of the magnetic polarization in the non linear regime. We found that the magnetic field is more amplified in paramagnetic fluids than in diamagnetic ones. Surprisingly, the effect of the susceptibility in the amplification is stronger when the magnetization parameter is smaller. These results make the KH instability a more efficient and effective amplification mechanism of seed magnetic fields when considering the susceptibility of matter.

We numerically simulate emergence of a magnetic plasma blob into a solar coronal hole. This blob may be associated with granulation and therefore it has a weak magnetic field. Two-dimensional simulations are performed using the MAGNUS code which solves magnetohydrodynamic equations, taking into account magnetic resistivity and thermal conduction. As a result of the interaction of the emerging blob with the ambient plasma, the magnetic lines experience reconnection with the blob getting flattened and deformed with time. Additionally, this process launches a vertical outflow of hot plasma and the chromosphere in its response increases its temperature. We perform parametric studies by varying the magnitude of the magnetic field of the blob and observing the net heating of the chromosphere. These studies are inspired by realistic simulations of granulation made with the use of two-fluid JOANNA code. In these simulations a number of magnetic blobs are detected in the convection zone and in the photosphere. From the numerical results, we conclude that as a result of granulation operating in a solar quiet region the emerging blob may trigger very complex dynamics in the upper regions of the solar atmosphere, and the associated outflows may be a source of heating of the chromosphere and possibly the solar corona.

We present the first family of magnetically polarized equilibrium tori around a Kerr black hole. The models were obtained in the test fluid approximation by assuming that the tori is a linear media, making it is possible to characterize the magnetic polarization of the fluid through the magnetic susceptibility $\chi_{m}$. The magnetohydrodynamic (MHD) structure of the models was solved by following the Komissarov approach, but with the aim of including the magnetic polarization of the fluid, the integrability condition for the magnetic counterpart was modified. We build two kinds of magnetized tori depending on whether the magnetic susceptibility is constant in space or not. In the models with constant $\chi_{m}$, we find that the paramagnetic tori ($\chi_{m}>0$) are more dense and less magnetized than the diamagnetic ones ($\chi_{m}<0$) in the region between the inner edge, $r_{in}$, and the center of the disk, $r_{c}$; however, we find the opposite behavior for $r>r_{c}$. Now, in the models with non-constant $\chi_{m}$, the tori become more magnetized than the Komissarov solution in the region where $\partial\chi_{m}/\partial r<0$, and less magnetized when $\partial\chi_{m}/\partial r>0$. Nevertheless, it is worth mentioning that in all solutions presented in this paper the magnetic pressure is greater than the hydrodynamic pressure. These new equilibrium tori can be useful for studying the accretion of a magnetic media onto a rotating black hole.

The magnetically polarized matter in astrophysical systems may be relevant in some magnetically dominated regions. For instance, the funnel that is generated in some highly magnetized disks configurations whereby relativistic jets are thought to spread, or in pulsars where the fluids are subject to very intense magnetic fields. With the aim of dealing with magnetic media in the astrophysical context, we present for the first time the conservative form of the GRMHD equations with a non-zero magnetic polarization vector $m^{\mu}$. Then, we follow the Anile method to compute the eigenvalue structure in the case where the magnetic polarization is parallel to the magnetic field, and it is parametrized by the magnetic susceptibility $\chi_m$. This approximation allows us to describe diamagnetic fluids, for which $\chi_m<0$, and paramagnetic fluids where $\chi_m>0$. The theoretical results were implemented in the CAFE code to study the role of the magnetic polarization in some 1D Riemann problems. We found that independently of the initial condition, the first waves that appear in the numerical solutions are faster in diamagnetic materials than in paramagnetic ones. Moreover, the constant states between the waves change notably for different magnetic susceptibilities. All these effects are more appreciable if the magnetic pressure is much bigger than the fluid pressure. Additionally, with the aim of analysing a magnetic media in a strong gravitational field, we carry out for the first time the magnetized Michel accretion of a magnetically polarized fluid. With this test, we found that the numerical solution is effectively maintained over time ($t>4000$), and that the global convergence of the code is $\gtrsim$ 2 for $\chi_{m}\lesssim 0.005$, for all the magnetic field strength $\beta$ we considered.

The origin of magnetic fields in the universe still remains unknown and constitutes one of the most intriguing questions in astronomy and astrophysics. Their significance is enormous since they have a strong influence on many astrophysical phenomena. In regards of this motivation, theoretical models of galactic disks with sources of magnetic field may contribute to understand the physics behind them. Inspired by this, we present a new family of analytical models for thin disks composed by magnetized material. The solutions are axially symmetric, conformastatic and are obtained by solving the Einstein-Maxwell Field Equations for continuum media without the test field approximation, and assuming that the sources are razor-thin disk of magnetically polarized matter. We find analytical expressions for the surface energy density, the pressure, the polarization vector, the electromagnetic fields, the mass and the rotational velocity for circular orbits, for two particular solutions. In each case, the energy-momentum tensor agrees with the energy conditions and also the convergence of the mass for all the solutions is proved. Since the solutions are well-behaved, they may be used to model astrophysical thin disks, and also may contribute as initial data in numerical simulations. In addition, the process to obtain the solutions is described in detail, which may be used as a guide to find solutions with magnetized material in General Relativity.

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 present a new magnetohydrodynamic (MHD) code for the simulation of wave propagation in the solar atmosphere, under the effects of electrical resistivity, but not dominant, and heat transference in a uniform 3D grid. The code is based on the finite volume method combined with the HLLE and HLLC approximate Riemann solvers, which use different slope limiters like MINMOD, MC, and WENO5. In order to control the growth of the divergence of the magnetic field, due to numerical errors, we apply the Flux Constrained Transport method, which is described in detail to understand how the resistive terms are included in the algorithm. In our results, it is verified that this method preserves the divergence of the magnetic fields within the machine round-off error. For the validation of the accuracy and efficiency of the schemes implemented in the code, we present some numerical tests in 1D and 2D for the ideal MHD. Later, we show one test for the resistivity in a magnetic reconnection process and one for the thermal conduction, where the temperature is advected by the magnetic field lines. Moreover, we display two numerical problems associated with the MHD wave propagation. The first one corresponds to a 3D evolution of a vertical velocity pulse at the photosphere-transition-corona region, while the second one consists in a 2D simulation of a transverse velocity pulse in a coronal loop.

An infinite family of exact solutions of the electrovacuum Einstein-Maxwell equations is presented. The family is static, axially symmetric and describe thin disks composed by electrically polarized material in a conformastatic spacetime. The form of the conformastatic metric allows us to write down the metric functions and the electromagnetic potentials in terms of a solution of the Laplace equation. We find a general expression for the surface energy density of the disk, the pressure, the polarization vector, the electromagnetic fields and the velocity rotation for circular orbits. As an example, we present the first model of the family and show the behavior of the different physical variables.

In the present paper, using the first-order approximation of the $n$-body Lagrangian (derived on the basis of the post-Newtonian gravitational theory of Einstein, Infeld, and Hoffman), we explicitly write down the equations of motion for the planar circular restricted three-body problem. Additionally, with some simplified assumptions, we obtain two formulas for estimating the values of the mass/distance and velocity/speed of light ratios appropriate for a given post-Newtonian approximation. We show that the formulas derived in the present study, lead to a good numerical conservation of the Jacobi constant and allow for an approximate equivalence between the Lagrangian and Hamiltonian approaches at the same post-Newtonian order. Accordingly, the dynamics of the system is analyzed in terms of the Poincar\'e sections method and Lyapunov exponents, finding that for specific values of the Jacobi constant the dynamics can be either chaotic or regular. Our results suggest that the chaoticity of the post-Newtonian system is slightly in- creased in comparison with its Newtonian counterpart.

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.

Removing outbursts from multiwavelength light curves of the blazar Mrk~421, we construct outburstless time series for this system. A model-independent power spectrum light curve analysis in the optical, hard X-ray and gamma-rays for this outburstless state and also the full light-curves, show clear evidence for a periodicity of ~ 310 days across all wavelengths studied. A subsequent full maximum likelihood analysis fitting an eclipse model confirms a periodicity of 310 \pm 1 days. The power spectrum of the signal in the outburstless state of the source does not follow a flicker noise behaviour and so, the system producing it is not self-organised. This and the fact that the periodicity is better defined in the outburstless state, strongly suggests that it is not produced by any internal physical processes associated to the central engine. The simplest physical mechanism to which this periodicity could be ascribed is a dynamical effect produced by an orbiting supermassive black hole companion eclipsing the central engine. Interestingly, the optimal eclipse model infers a brightness enhancement of (136.4 \pm 20 )%, suggesting an eclipse resulting in a gravitational lensing brightening. Consisting with this interpretation, the eclipse occurs for only ( 9.7 \pm 0.2)% of the orbital period.

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 study the rotation curves of ultralight BEC dark matter halos. These halos are long lived solutions of initially rotating BEC fluctuations. In order to study the implications of the rotation characterizing these long-lived configurations we consider the particular case of a boson mass $m=10^{-23}\mathrm{eV/c}^2$ and no self-interaction. We find that these halos successfully fit samples of rotation curves (RCs) of LSB galaxies.

We present the accretion of collisional dark matter on a supermassive black hole seed. The analysis is based on the numerical solution of the fully coupled system of Einstein-Euler equations for spherically symmetric flow, where the dark matter is modeled as a perfect fluid that obeys an ideal gas equation of state. As the black hole actually grows, the accretion rate of dark matter corresponds to the black hole apparent horizon growth rate. We analyse cases with infall velocity as high as $0.5c$ and an environment density of $100M_{\odot}/\mathrm{pc}^3$, which are rather extreme conditions. Being the radial flux the maximum accretion case, our results show that the accretion of an ideal gas, eventually collisional dark matter, does not contribute significantly to SMBH masses. This result favors models predicting SMBHs were formed already with supermasses. We show that despite the fact that we are solving the full general relativistic system, for the parameter space studied our results are surprisingly similar to those obtained using the Bondi formula, which somehow certifies its use as a good approximation of a fully evolving space-time with spherical symmetry at short scales at least for dark matter densities. Additionally, we study the density profile of the gas and find that the presence of SMBHs redistributes the gas near the event horizon with a cuspy profile, whereas beyond a small fraction of a parsec it is not-cuspy anymore.

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.

We present the dynamics of rotating Bose Condensate galactic dark matter halos, made of an ultralight spinless boson. We restrict to the case of adding axisymmetric rigid rotation to initially spherically symmetric structures and show there are three regimes: i) small angular momentum, that basically retains the drawbacks of spherically symmetric halos related to compactness and failure at explaining galactic RCs, ii) an intermediate range of values of angular momentum that allow the existence of long-lived structures with acceptable RC profiles, and iii) high angular momentum, in which the structure is dispersed away by rotation. We also present in detail the new code used to solve the Gross-Pitaevskii Poisson system of equations in three dimensions.

In this paper we show that spherically symmetric BEC dark matter halos, with the $\sin r/r$ density profile, that accurately fit galactic rotation curves and represent a potential solution to the cusp-core problem are unstable. We do this by introducing back the density profiles into the fully time-dependent Gross-Pitaevskii-Poisson system of equations. Using numerical methods to track the evolution of the system, we found that these galactic halos lose mass at an approximate rate of half of its mass in a time scale of dozens of Myr. We consider this time scale is enough as to consider these halos are unstable and unlikely to be formed. We provide some arguments to show that this behavior is general and discuss some other drawbacks of the model that restrict its viability.

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.

Some of the most interesting scenarios that can be studied in astrophysics, contain fluids and plasma moving under the influence of strong gravitational fields. To study these problems it is required to implement numerical algorithms robust enough to deal with the equations describing such scenarios, which usually involve hydrodynamical shocks. It is traditional that the first problem a student willing to develop research in this area is to numerically solve the one dimensional Riemann problem, both Newtonian and relativistic. Even a more basic requirement is the construction of the exact solution to this problem in order to verify that the numerical implementations are correct. We describe in this paper the construction of the exact solution and a detailed procedure of its implementation.

We study numerically the axisymmetric relativistic Bondi-Hoyle accretion of a supersonic ideal gas onto a fixed Schwarzschild background space-time described with horizon penetrating coordinates. We verify that a nearly stationary shock cone forms and that the properties of the shock cone are consistent with previous results in Newtonian gravity and former relativistic studies. The fact that the evolution of the gas is tracked on a spatial domain that contains a portion of the inner part of the black hole avoids the need to impose boundary conditions on a time-like boundary as done in the past. Thus, our approach contributes to the solution to the Bondi-Hoyle accretion problem at the length scale of the accretor in the sense that the gas is genuinely entering the black hole. As an astrophysical application, we study for a set of particular physical parameters, the spectrum of the shock cone vibrations and their potential association with QPO sources.

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 accretion of a phantom scalar field into a black hole for various scalar field potentials in the full non-linear regime. Our results are based on the use of numerical methods and show that for all the cases studied the black hole's apparent horizon mass decreases. We explore a particular subset of the parameter space and from our results we conclude that this is a very efficient black hole shrinking process because the time scales of the area reduction of the horizon are short. We show that the radial equation of state of the scalar field depends strongly on the space and time, with the condition $\omega = p/\rho>-1$, as opposed to a phantom fluid at cosmic scales that allows $\omega < -1$.

In this paper we track the non-linear spherical evolution of a massless scalar field onto a Schwarzschild black hole space-time as a first approximation to the accretion of cosmologically motivated classical scalar fields. We perform an analysis related to wave packets described by wave number and width. We study various values of the wave number k, and found that for k = 0 and width packets bigger than the Schwarzschild radius, the absorption is not total. In the cases we studied for k > 0, the black hole absorbs the total amount of energy density of the scalar field moving toward the horizon. Our results indicate that assuming spherical symmetry, in the non-linear regime, there are cases for which scalar fields are allowed to survive outside black holes and may eventually have life-times consistent with cosmological time scales.

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 solve the spherically symmetric time dependent relativistic Euler equations on a Schwarzschild background space-time for a perfect fluid, where the perfect fluid models the dark matter and the space-time background is that of a non-rotating supermassive black hole. We consider the fluid obeys an ideal gas equation of state as a simple model of dark matter with pressure. Assuming out of equilibrium initial conditions we search for late-time attractor type of solutions, which we found to show a constant accretion rate for the non-zero pressure case, that is, the pressure itself suffices to produce stationary accretion regimes. We then analyze the resulting density profile of such late-time solutions with the function $A/r^{\kappa}$. For different values of the adiabatic index we find different slopes of the density profile, and we study such profile in two regions: a region one near the black hole, located from the horizon up to 50$M$ and a region two from $\sim 800M$ up to $\sim 1500M$, which for a black hole of $10^{9}M_{\odot}$ corresponds to $\sim 0.1$pc. The profile depends on the adiabatic index or equivalently on the pressure of the fluid and our findings are as follows: in the near region the density profile shows values of $\kappa <1.5$ and in the limit of the pressure-less case $\kappa \rightarrow 1.5$; on the other hand, in region two, the value of $\kappa<0.3$ in all the cases we studied. If these results are to be applied to the dark matter problem, the conclusion is that, in the limit of pressure-less gas the density profile is cuspy only near the black hole and approaches a non-cuspy profile at bigger scales within 1pc. These results show on the one hand that pressure suffices to provide flat density profiles of dark matter and on the other hand show that the presence of a central black hole does not distort the density profile of dark matter at scales of 0.1pc.

Based on the numerical solution of the time-dependent relativistic Euler equations onto a fixed Schwarzschild background space-time, we estimate the accretion rate of radial flow toward the horizon of a test perfect fluid obeying an ideal gas equation of state. We explore the accretion rate in terms of the initial density of the fluid for various values of the inflow velocity in order to investigate whether or not sufficiently arbitrary initial conditions allow a steady state accretion process depending on the values of the pressure. We extrapolate our results to the case where the fluid corresponds to dark matter and the black hole is a supermassive black hole seed. Then we estimate the equation of state parameters that provide a steady state accretion process. We found that when the pressure of the dark matter is zero, the black hole's mass grows up to values that are orders of magnitude above $10^{9}M_{\odot}$ during a lapse of 10Gyr, whereas in the case of the accretion of the ideal gas dark matter with non zero pressure the accreted mass can be of the order of $\sim 1M_{\odot}/10Gyr$ for black holes of $10^{6}M_{\odot}$. This would imply that if dark matter near a supermassive black hole acquires an equation of state with non trivial pressure, the contribution of accreted dark matter to the supermassive black hole growth could be small, even though only radial accretion is considered.

We present the first analytical superposition of a charged black hole with an annular disk of extremal dust. In order to obtain the solutions, we first solve the Einstein-Maxwell field equations for sources that represent disk-like configurations of matter in confomastatic spacetimes by assuming a functional dependence among the metric function, the electric potential and an auxiliary function,which is taken as a solution of the Laplace equation. We then employ the Lord Kelvin Inversion Method applied to models of finite extension in order to obtain annular disks. The structures obtained extend to infinity, but their total masses are finite and all the energy conditions are satisfied. Finally, we observe that the extremal Reissner-Nordstr\"{o}m black hole can be embedded into the center of the disks by adding a boundary term in the inversion.

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.