In 2014 the NGC 5548 Space Telescope and Optical Reverberation Mapping campaign discovered a two-month anomaly when variations in the absorption and emission lines decorrelated from continuum variations. During this time the soft X-ray part of the intrinsic spectrum had been strongly absorbed by a line-of-sight (LOS) obscurer, which was interpreted as the upper part of a disk wind. Our first paper showed that changes in the LOS obscurer produce the decorrelation between the absorption lines and the continuum. A second study showed that the base of the wind shields the BLR, leading to the emission-line decorrelation. In that study, we proposed the wind is normally transparent with no effect on the spectrum. Changes in the wind properties alter its shielding and affect the SED striking the BLR, producing the observed decorrelations. In this work, we investigate the impact of a translucent wind on the emission lines. We simulate the obscuration using XMM-Newton, NuSTAR, and HST observations to determine the physical characteristics of the wind. We find that a translucent wind can contribute a part of the He II and Fe K? emission. It has a modest optical depth to electron scattering, which explains the fainter far-side emission in the observed velocity delay maps. The wind produces the very broad base seen in the UV emission lines and may also be present in the Fe K? line. Our results highlight the importance of accounting for the effects of such winds in the analysis of the physics of the central engine.

The 180-day Space Telescope and Optical Reverberation Mapping campaign on NGC 5548 discovered an anomalous period, the broad-line region (BLR) holiday, in which the emission lines decorrelated from the continuum variations. This is important since the correlation between the continuum-flux variations and the emission-line response is the basic assumption for black hole (BH) mass determinations through reverberation mapping. During the BLR holiday, the high-ionization intrinsic absorption lines also decorrelated from the continuum as a result of variable covering factor of the line of sight (LOS) obscurer. The emission lines are not confined to the LOS, so this does not explain the BLR holiday. If the LOS obscurer is a disk wind, its streamlines must extend down to the plane of the disk and the base of the wind would lie between the BH and the BLR, forming an equatorial obscurer. This obscurer can be transparent to ionizing radiation, or can be translucent, blocking only parts of the SED, depending on its density. An emission-line holiday is produced if the wind density increases only slightly above its transparent state. Both obscurers are parts of the same wind, so they can have associated behavior in a way that explains both holidays. A very dense wind would block nearly all ionizing radiation, producing a Seyfert 2 and possibly providing a contributor to the changing-look AGN phenomenon. Disk winds are very common and we propose that the equatorial obscurers are too, but mostly in a transparent state.

In this manuscript the three types of late-time behavior of spherically symmetric Boson Stars are presented, these are: stable configurations, unstable bounded that collapse to form black holes and unstable unbounded that explode. The results are found by solving the spherically symmetric Einstein-Klein-Gordon system of equations for a complex scalar field with initial conditions corresponding to a Boson Star. The solution is based on a constrained evolution that uses the method of lines for the scalar field and solves the constraint equations for the geometry on the fly.

We study binary mergers of ultralight bosonic dark matter cores by solving the Gross-Pitaevskii-Poisson system of equations. The analysis centers on the dynamics of the relaxation process and the behavior of the configuration resulting from the merger, including the Gravitational Cooling with its corresponding emission of mass and angular momentum. The oscillations of density and size of the final configuration are characterized, indicating that for the equal mass case the dependency of the amplitude and frequency of these oscillations on the impact parameter of the pre-merger configuration is linear. The amplitude of these oscillations changes by a factor of two or more indicating the final configuration does not approach a clear stationary state even though it oscillates around a virialized state. For the unequal mass case, global quantities also indicate the final configuration oscillates around a virialized state, although the density does not show a dominant oscillation mode. The evolution of the angular momentum prior and post merger is analyzed in all cases.

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches.

The flux variations in the emission lines in active galactic nuclei (AGNs) are driven by variations in the ionizing continuum flux --which are usually reflected in the observable UV-optical continuum. The "Reverberation mapping" technique measures the delay between line and continuum variations to determine the size of the line emitting region, this is the basis for measurements of the central black hole mass in AGNs. The Space Telescope and Optical Reverberation Mapping Project (AGN STORM) on NGC 5548 in 2014 is the most intensive multi-wavelength AGN monitoring campaign ever. For most of the campaign, the emission-line variations followed changes in the continuum with a time lag, as expected. However, the lines varied independently of the observed UV-optical continuum during a 60 -- 70 day "holiday." To understand this remarkable phenomenon, we study the intrinsic absorption lines present in NGC 5548. We identify a novel cycle that reproduces the absorption line variability and thus identify the physics that allows the holiday to occur. In our model, variations in this obscurer's line-of-sight covering factor modify the soft X-ray continuum. This leads to changes in the ionization of helium gas in the broad-line region. Ionizing radiation produced by recombining helium then affects the ionization of other species as observed during the AGN STORM holiday. It is likely that any other model which selectively changes the soft X-ray part of the continuum during the holiday can also explain the anomalous emission line behavior observed.

Structure formation simulations of ultralight bosonic dark matter, ruled by the Gross-Pitaevskii-Poisson (GPP) system of equations, show that dark matter clumps to form structures with density profiles consisting of a core surrounded by a power-law distribution. The core has a density profile similar to that of spherically symmetric equilibrium configurations of the GPP system. These configurations have been shown to be stable under a variety of perturbations and to have attractor properties. It is interesting to know the dominant frequencies of oscillation of these configurations. One reason is that in galaxies the effects of perturbations can trigger observable effects for specific frequency modes. Another reason is that during the process of structure formation, the oscillation modes can help to characterize a core. Based on these motivations, in this manuscript we present a systematic numerical analysis of the reaction of equilibrium configurations to various axi-symmetric perturbations with modes $l=0,1,2,3,4$. We then calculate the first few oscillation frequencies of equilibrium configurations for each mode.

We present a method for estimating the velocity of a wandering black hole and the equation of state for the gas around, based on a catalog of numerical simulations. The method uses machine learning methods based on convolutional neural networks applied to the classification of images resulting from numerical simulations. Specifically we focus on the supersonic velocity regime and choose the direction of the black hole to be parallel to its spin. We build a catalog of 900 simulations by numerically solving Euler's equations onto the fixed space-time background of a black hole, for two parameters: the adiabatic index $\Gamma$ with values in the range [1.1, 5/3], and the asymptotic relative velocity of the black hole with respect to the surroundings $v_{\infty}$, with values within $[0.2, 0.8]c$. For each simulation we produce a 2D image of the gas density once the process of accretion has approached a stationary regime. The results obtained show that the implemented Convolutional Neural Networks are capable to classify correctly the adiabatic index $87.78\%$ of the time within an uncertainty of $\pm 0.0284$ while the prediction of the velocity is correct $96.67\%$ of the times within an uncertainty of $\pm 0.03c$. We expect that this combination of a massive number of numerical simulations and machine learning methods will help analyze more complicated scenarios related to future high resolution observations of black holes, like those from the Event Horizon Telescope.

Within the context of ultralight BEC dark matter, we analyze the head-on encounters of two structures. These structures are made of a BEC component, which is a ground state equilibrium solution of the Gross-Pitaevskii-Poisson system, together with a component of luminous matter. The evolution of the Condensate dark matter is carried out by solving the time dependent GPP equations, whereas the luminous matter is modeled with particles interacting gravitationally on top of the BEC dark matter halos. We track the evolution of frontal encounters for various values of the collision velocity and analyze the regime of high velocity regime showing solitonic behavior of the BEC halos and that of slow velocities producing a single final structure. We measure the relative velocity of the dark matter with respect to the luminous matter after the encounters in the solitonic case and track the evolution of luminous matter in the case of merger.

We explore the possibility that a single relativistic shock, where the gas dynamics is coupled with radiation, can fit the light curves of long GRBs. For this we numerically solve the one dimensional relativistic radiation hydrodynamics equations with a single initial shock. We calculate light curves due to the evolution of this shock in terms of the velocity of the shock, the opacity of the gas, mass density and density of radiated energy. We explore how the variation of each of these parameters provides different features in the light curves. As examples we include the fitting of two long GRBs.

Using numerical simulations we present the accretion of supersonic winds onto a rotating black hole in three dimensions. We study five representative directions of the wind with respect to the axis of rotation of the black hole and focus on the evolution and stability of the high density shock cone that is formed during the process. We explore both, the regime in which the shock cone is expected to be stable in order to confirm previous results obtained with two dimensional simulations, and the regime in which the shock cone is expected to show a flip-flop type of instability. The methods used to attempt triggering the instability were first the accumulation of numerical errors and second the explicit application of a perturbation on the velocity field after the shock-cone was formed. The result is negative, that is, we did not find the flip-flop instability within the parameter space we explored, which includes cases that are expected to be unstable.

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.

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 solve numerically the ideal MHD equations with an external gravitational field in 2D in order to study the effects of impulsively generated linear and non-linear Alfv\'en waves into isolated solar arcades and coronal funnels. We analyze the region containing the interface between the photosphere and the corona. The main interest is to study the possibility that Alfv\'en waves triggers the energy flux transfer toward the quiet solar corona and heat it, including the case that two consecutive waves can occur. We find that in the case of arcades, short or large, the transferred fluxes by Alfv\'en waves are sufficient to heat the quiet corona only during a small lapse of time and in a certain region. In the case of funnels the threshold is achieved only when the wave is faster than 10 km/s, which is extremely high. We conclude from our analysis, that Alfv\'en waves, even in the optimistic scenario of having two consecutive Alfv\'en wave pulses, cannot transport enough energy as to heat the quiet corona.

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.

The status of boson stars as black hole mimickers is presented. We focus on the analysis of the emission spectrum of a simple accretion disk model. We describe the free parameters that allow a boson star to become a black hole mimicker and present an example of a particular astrophysical case.

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.

In order to explore nonlinear effects on the distribution of matter during collisions within the Bose-Einstein condensate (BEC) dark matter model driven by the Schr\"odinger-Poisson system of equations, we study the head-on collision of structures and focus on the interference pattern formation in the density of matter during the collision process. We explore the possibility that the collision of two structures of fluid matter modeled with an ideal gas equation of state also forms interference patterns and found a negative result. Given that a fluid is the most common flavor of dark matter models, we conclude that one fingerprint of the BEC dark matter model is the pattern formation in the density during a collision of structures.

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 spherically symmetric boson stars as black hole mimickers based on the power spectrum of a simple accretion disk model. The free parameters of the boson star are the mass of the boson and the fourth order self-interaction coefficient in the scalar field potential. We show that even if the mass of the boson is the only free parameter it is possible to find a configuration that mimics the power spectrum of the disk due to a black hole of the same mass. We also show that for each value of the self-interaction a single boson star configuration can mimic a black hole at very different astrophysical scales in terms of the mass of the object and the accretion rate. In order to show that it is possible to distinguish one of our mimickers from a black hole we also study the deflection of light.

Using numerical methods we present the first full nonlinear study of phantom scalar field accreted into a black hole. We study different initial configurations and find that the accretion of the field into the black hole can reduce its area down to 50 percent within time scales of the order of few masses of the initial horizon. The analysis includes the cases where the total energy of the space-time is positive or negative. The confirmation of this effect in full nonlinear general relativity implies that the accretion of exotic matter could be considered an evaporation process. We speculate that if this sort of exotic matter has some cosmological significance, this black hole area reduction process might have played a crucial role in black hole formation and population.

In this manuscript we track the evolution of a system consisting of two self-gravitating virialized objects made of a scalar field in the newtonian limit. The Schr\"odinger-Poisson system contains a potential with self-interaction of the Gross-Pitaevskii type for Bose Condensates. Our results indicate that solitonic behavior is allowed in the scalar field dark matter model when the total energy of the system is positive, that is, the two blobs pass through each other as should happen for solitons; on the other hand, there is a true collision of the two blobs when the total energy is negative.

We show the evolution of non-spherically symmetric balls of a self-gravitating scalar field in the Newtonian regime or equivalently an ideal self-gravitating condensed Bose gas. In order to do so, we use a finite differencing approximation of the Shcr\"odinger-Poisson (SP) system of equations with axial symmetry in cylindrical coordinates. Our results indicate: 1) that spherically symmetric ground state equilibrium configurations are stable against non-spherical perturbations and 2) that such configurations of the SP system are late-time attractors for non-spherically symmetric initial profiles of the scalar field, which is a generalization of such behavior for spherically symmetric initial profiles. Our system and the boundary conditions used, work as a model of scalar field dark matter collapse after the turnaround point. In such case, we have found that the scalar field overdensities tolerate non-spherical contributions to the profile of the initial fluctuation.

Equilibrium configurations for a self-gravitating scalar field with self-interaction are constructed. The corresponding Schr\"odinger-Poisson (SP) system is solved using finite differences assuming spherical symmetry. It is shown that equilibrium configurations of the SP system are late-time attractor solutions for initially quite arbitrary density profiles, which relax and virialize through the emission of scalar field bursts; a process dubbed gravitational cooling. Among other potential applications, these results indicate that scalar field dark matter models (in its different flavors) tolerate the introduction of a self-interaction term in the SP equations. This study can be useful in exploring models in which dark matter in galaxies is not point-like.

The emission spectrum from a simple accretion disk model around a compact object is compared for the cases of a black hole (BH) and a boson star (BS) playing the role of the central object. It was found in the past that such a spectrum presents a hardening at high frequencies; however, here it is shown that the self-interaction and compactness of the BS have the effect of softening the spectrum, the less compact the star is, the softer the emission spectrum at high frequencies. Because the mass of the boson fixes the mass of the star and the self-interaction the compactness of the star, we find that, for certain values of the BS parameters, it is possible to produce similar spectra to those generated when the central object is a BH. This result presents two important implications: (i) using this simple accretion model, a BS can supplant a BH in the role of compact object accreting matter, and (ii) within the assumptions of the present accretion disk model we do not find a prediction that could help distinguish a BH from a BS with appropriate parameters of mass and self-interaction.

Using numerical techniques, we study the collapse of a scalar field configuration in the Newtonian limit of the spherically symmetric Einstein--Klein--Gordon (EKG) system, which results in the so called Schr\"odinger--Newton (SN) set of equations. We present the numerical code developed to evolve the SN system and topics related, like equilibrium configurations and boundary conditions. Also, we analyze the evolution of different initial configurations and the physical quantities associated to them. In particular, we readdress the issue of the gravitational cooling mechanism for Newtonian systems and find that all systems settle down onto a 0--node equilibrium configuration.

The present effort addresses the question about the existence of a well-defined thermodynamic limit for the astrophysical systems with the following power law form: to tend the number of particles, N, the total energy, E, and the characteristic linear dimension of the system, L, to infinity, keeping constant E/N^{\Lambda_{E}} and L/N^{\Lambda_{L}}, being \Lambda_{E} and \Lambda_{L} certain scaling exponent constant. This study is carried out for a system constituted by a non-rotating fluid under the influence of its own Newtonian gravitational interaction. The analysis yields that a thermodynamic limit of the above form will only appear when the local pressure depends on the energy density of fluid as $p=\gamma\epsilon$, being $\gamma$ certain constant. Therefore, a thermodynamic limit with a power law form can be only satisfied by a reduced set of models, such as the selfgravitating gas of fermions and the Antonov isothermal model.

Astrophysical systems will never be in a real Thermodynamic equilibrium: they undergo an evaporation process due to the fact that the gravity is not able to confine the particles. Ordinarily, this difficulty is overcome by enclosing the system in a rigid container which avoids the evaporation. We proposed an energetic prescription which is able to confine the particles, leading in this way to an alternative version of the Antonov isothermal model which unifies the well-known isothermal and polytropic profiles. Besides of the main features of the isothermal sphere model: the existence of the gravitational collapse and the energetic region with a negative specific heat, this alternative model has the advantage that the system size naturally appears as a consequence of the particles evaporation.

We reconsider some general aspects about the mean field thermodynamical description of the astrophysical systems based on the microcanonical ensemble. Starting from these basis, we devote a special attention to the analysis of the scaling laws of the thermodynamical variables and potentials in the thermodynamic limit. Geometrical considerations motivate a way by means of which could be carried out a well-defined generalized canonical-like description for this kind of systems, even being nonextensive. This interesting possibility allows us to extend the applicability of the Standard Thermodynamic methods, even in the cases in which the system exhibits a negative specific heat. As example of application, we reconsider the classical Antonov problem of the isothermal spheres.

According to self-similarity hypothesis, the thermodynamic limit could be defined from the scaling laws for the system self-similarity by using the microcanonical ensemble. This analysis for selfgravitating systems yields the following thermodynamic limit: send N to infinity, keeping constant E/N^{(7/3)} and LN^{(1/3)}, in which is ensured the extensivity of the Boltzmann entropy S_{B}=lnW(E,N). It is shown how the consideration of this thermodynamic limit allows us to explain the origin of dynamical anomalies in numerical simulations of selfgravitating systems.

In this letter, we develop a Newtonian approach to the collapse of galaxy fluctuations of scalar field dark matter under initial conditions inferred from simple assumptions. The full relativistic system, the so called Einstein-Klein-Gordon, is reduced to the Schr\"odinger-Newton one in the weak field limit. The scaling symmetries of the SN equations are exploited to track the non-linear collapse of single scalar matter fluctuations. The results can be applied to both real and complex scalar fields.

In an arbitrary axisymmetric stationary spacetime, we determine the expression for the tangential velocity of test objects following a circular stable geodesic motion in the equatorial plane, as function of the metric coefficients. Next, we impose the condition, observed in large samples of disks galaxies, that the magnitude of such tangential velocity be radii independent, obtaining a constraint equation among the metric coefficients, and thus arriving to an iff condition: The tangential velocity of test particles is radii independent iff the metric coefficients satisfied the mentioned constraint equation. Furthermore, for the static case, the constraint equation can be easily integrated, leaving the spacetime at the equatorial plane essentially with only one independent metric coefficient. With the geometry thus fixed, we compute the Einstein tensor and equate it to and arbitrary stress energy tensor, in order to determined the type of energy-matter which could produce such a geometry. Within an approximation, we deduce a constraint equation among the components of the stress energy tensor. We test in that constraint equation several well known types of matter, which have been proposed as dark matter candidates and are able to point for possible right ones. Finally, we also present the spherically symmetric static case and apply the mentioned procedure to perfect fluid stress energy tensor, recovering the Newtonian result as well as the one obtained in the axisymmetric case. We also present arguments on the need to use GR to study non dust types of matter.

We present a general description of the scalar field dark matter (SFDM) hypothesis in the cosmological context. The scenario of structure formation under such a hypothesis is based on Jeans instabilities of fluctuations of the scalar field. It is shown that it is possible to form stable long lived objects consisting of a wide range of typical galactic masses around $10^{12}M_{\odot}$ once the parameters of the effective theory are fixed with the cosmological constraints. The energy density at the origin of such an object is smooth as it should.

We present a scenario for galaxy formation based on the hypothesis of scalar field dark matter. We interpret galaxy formation through the collapse of a scalar field fluctuation. We find that a cosh potential for the self-interaction of the scalar field provides a reasonable scenario for galactic formation, which is in agreement with cosmological observations and phenomenological studies in galaxies.

We present an exact solution of the averaged Einstein's field equations in the presence of two real scalar fields and a component of dust with spherical symmetry. We suggest that the space-time found provides the characteristics required by a galactic model that could explain the supermassive central object and the dark matter halo at once, since one of the fields constitutes a central oscillaton surrounded by the dust and the other scalar field distributes far from the coordinate center and can be interpreted as a halo. We show the behavior of the rotation curves all along the background. Thus, the solution could be a first approximation of a ``long exposition photograph'' of a galaxy.

This work is a review of the last results of research on the Scalar Field Dark Matter model of the Universe at cosmological and at galactic level. We present the complete solution to the scalar field cosmological scenario in which the dark matter is modeled by a scalar field $\Phi$ with the scalar potential $V(\Phi)=V_{0}(cosh {(\lambda \sqrt{\kappa_{0}}\Phi)}-1)$ and the dark energy is modeled by a scalar field $\Psi$, endowed with the scalar potential $\tilde{V}(\Psi)= \tilde{V_{0}}(\sinh{(\alpha \sqrt{\kappa_{0}}\Psi)})^{\beta}$, which together compose the 95% of the total matter energy in the Universe. The model presents successfully deals with the up to date cosmological observations, and is a good candidate to treat the dark matter problem at the galactic level.

We study a spherically symmetric fluctuation of scalar dark matter in the cosmos and show that it could be the dark matter in galaxies, provided that the scalar field has an exponential potential whose overall sign is negative and whose exponent is constrained observationally by the rotation velocities of galaxies. The local space-time of the fluctuation contains a three dimensional space-like hypersurface with surplus of angle.

An exact, axially symmetric solution to the Einstein-Klein-Gordon field equations is employed to model the dark matter in spiral galaxies. The extended rotation curves from a previous analysis are used to fit the model and a very good agreement is found. It is argued that, although our model possesses three parameters to be fitted, it is better than the non-relativistic alternatives in the sense that it is not of a phenomenological nature, since the dark matter would consist entirely of a scalar field.

Through the geodesic analysis of a static and axially symmetric space time, we present conditions on the state equation of an isotropic perfect fluid $p=\omega d$, when it is considered as dark matter in spiral galaxies. The main conclusion is that it can be an exotic fluid ($-1<\omega <-1/3$) as it is found for Quintessence at cosmological scale.

We investigate the hypothesis that the scalar field is the dark matter and the dark energy in the Cosmos, wich comprises about 95% of the matter of the Universe. We show that this hypothesis explains quite well the recent observations on type Ia supernovae.

Recently it has been proposed that the main contributor to the dark energy of the Universe is a dynamical, slow evolving, spatially inhomogeneous scalar field called quintessence. We investigate the behavior of this scalar field at galactic level by assuming that it is the dark matter compossing the halos of galaxies. Using an exact solution of the Einstein's equations we find an excellent concordance between our results and observations.

We present a model for the dark matter in spiral galaxies, which is a result of a static and axial symmetric exact solution of the Einstein-Dilaton theory. We suposse that dark matter is a scalar field endowed with a scalar potential. We obtain that a) the effective energy density goes like $1/(r^2+r_{c}^{2})$ and b) the resulting circular velocity profile of tests particles is in good agreement with the observed one.