Source author record

F. S. Guzman

F. S. Guzman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2016arXiv

Accretion of Supersonic Winds on Boson Stars

We present the evolution of a supersonic wind interacting with a Boson Star (BS) and compare the resulting wind density profile with that of the shock cone formed when the wind is accreted by a non-rotating Black Hole (BH) of the same mass. The physical differences between these accretors are that a BS, unlike a BH has no horizon, it does not have a mechanical surface either and thus the wind is expected to trespass the BS. Despite these conditions, on the BS space-time the gas achieves a stationary flux with the gas accumulating in a high density elongated structure comparable to the shock cone formed behind a BH. The highest density resides in the center of the BS whereas in the case of the BH it is found on the downstream part of the BH near the event horizon. The maximum density of the gas is smaller in the BS than in the BH case. Our results indicate that the highest density of the wind is more similar on the BS to that on the BH when the BS has high self-interaction, when it is more compact and when the wind velocity is higher. We expect this and similar analyses help to know whether BSs can still be Black Hole mimickers or can be ruled out.

preprint2016arXiv

Behavior of luminous matter in the head-on encounter of two ultralight BEC dark matter halos

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.

preprint2016arXiv

Modeling long GRBs using a single shock with relativistic radiation hydrodynamics

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.

preprint2015arXiv

Accretion of supersonic winds onto black holes in 3D: stability of the shock cone

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.

preprint2015arXiv

CAFE: A New Relativistic MHD Code

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.

preprint2015arXiv

Estimating the contribution of Alfvén waves to the process of heating the quiet solar corona

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én 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én 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én 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én waves, even in the optimistic scenario of having two consecutive Alfvén wave pulses, cannot transport enough energy as to heat the quiet corona.

preprint2015arXiv

Newtonian CAFE: a new ideal MHD code to study the solar atmosphere

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.

preprint2015arXiv

Rotation curves of ultralight BEC dark matter halos with rotation

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.

preprint2014arXiv

Horizon growth of supermassive black hole seeds fed with collisional dark matter

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.

preprint2013arXiv

Exact solution of the 1D Riemann problem in Newtonian and relativistic hydrodynamics

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.

preprint2013arXiv

PBH mass growth through radial accretion during the radiation dominated era

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=ρ/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.

preprint2013arXiv

Rotation curves of rotating galactic BEC dark matter halos

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.

preprint2013arXiv

Stability of BEC galactic dark matter halos

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.

preprint2013arXiv

The shocks during the accretion of an ultrarelativistic supersonic gas onto a rotating black hole

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

preprint2012arXiv

Axisymmetric Bondi-Hoyle accretion onto a Schwarzschild Black Hole: shock cone vibrations

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.

preprint2012arXiv

Behavior of Phantom Scalar Fields near Black Holes

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 $ω= p/ρ>-1$, as opposed to a phantom fluid at cosmic scales that allows $ω< -1$.

preprint2012arXiv

Is the flip-flop behaviour of accretion shock cones on to black holes an effect of coordinates?

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.

preprint2012arXiv

Revisiting spherically symmetric relativistic hydrodynamics

In this paper we revise two classical examples of Relativistic Hydrodynamics in order to illustrate in detail the numerical methods commonly used in fluid dynamics, specifically those designed to deal with shocks, which are based on a finite volume approximation. The two cases we consider are the relativistic blast wave problem and the evolution of a Tolman-Oppenheimer-Volkoff star model, in spherical symmetry. In the first case we illustrate the implementation of relativistic Euler's equations on a fixed background space-time, whereas in the second case we also show how to couple the evolution of the fluid to the evolution of the space-time.

preprint2012arXiv

Spherical non-linear absorption of cosmological scalar fields onto a black hole

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.

preprint2011arXiv

Collisional dark matter density profiles around supermassive black holes

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^κ$. 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 $κ<1.5$ and in the limit of the pressure-less case $κ\rightarrow 1.5$; on the other hand, in region two, the value of $κ<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.

preprint2011arXiv

Exploring the effects of pressure on the radial accretion of dark matter by a Schwarzschild supermassive black hole

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.

preprint2011arXiv

Interference pattern in the collision of structures in the BEC dark matter model: comparison with fluids

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ödinger-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.

preprint2011arXiv

Scalar Field Dark Matter: behavior around black holes

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

preprint2010arXiv

Numerical solution of the wave equation on particular space-times using CMC slices and scri-fixing conformal compactification

In this paper we present in detail the numerical solution of the conformally invariant wave equation on top of a fixed background space-time corresponding to two different cases: i) 1+1 Minkowski space-time in Cartesian coordinates and ii) Schwarzschild space-time. In both cases we use hyperboloidal constant mean curvature slices and scri-fixing conformal compactification, and solve the wave equation on the conformal space-time. In the case of the Schwarzschild space-time we study the quasinormal mode oscillations and the late-time polynomial tail decay exponents corresponding to a mass-less scalar field. We also present general formulas to construct hyperboloidal constant mean curvature slicings of spherically symmetric, static, space-times in spherical coordinates.

preprint2010arXiv

Spherical Boson Stars as Black Hole mimickers

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.

preprint2009arXiv

Accretion of phantom scalar field into a black hole

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.

preprint2009arXiv

On the instability of charged wormholes supported by a ghost scalar field

In previous work, we analyzed the linear and nonlinear stability of static, spherically symmetric wormhole solutions to Einstein's field equations coupled to a massless ghost scalar field. Our analysis revealed that all these solutions are unstable with respect to linear and nonlinear spherically symmetric perturbations and showed that the perturbation causes the wormholes to either decay to a Schwarzschild black hole or undergo a rapid expansion. Here, we consider charged generalization of the previous models by adding to the gravitational and ghost scalar field an electromagnetic one. We first derive the most general static, spherically symmetric wormholes in this theory and show that they give rise to a four-parameter family of solutions. This family can be naturally divided into subcritical, critical and supercritical solutions depending on the sign of the sum of the asymptotic masses. Then, we analyze the linear stability of these solutions. We prove that all subcritical and all critical solutions possess one exponentially in time growing mode. It follows that all subcritical and critical wormholes are linearly unstable. In the supercritical case we provide numerical evidence for the existence of a similar unstable mode.

preprint2009arXiv

On Wormholes supported by phantom energy

By a combination of analytical and numerical techniques, we demonstrate the existence of spherical, asymptotically flat traversable wormholes supported by exotic matter whose stress tensor relative to the orthonormal frame of Killing observers takes the form of a perfect fluid possessing anisotropic pressures and subject to linear equations of state: $τ=λρc^{2}$, $P=μρc^{2}$. We show that there exists a four parameter family of asymptotically flat spherical wormholes parametrized by the area of the throat A(0), the gradient $Λ(0)$ of the red shift factor evaluated at the throat as well as the values of $(λ, μ)$. The latter are subject to restrictions: $λ>1$ and $2μ>λ$ or $λ<0$ and $2μ<-|λ|$. For particular values of $(λ, μ)$, the stress tensor may be interpreted as representing a phantom configuration, while for other values represents exotic matter. All solutions have the property that the two asymptotically flat ends posses finite ADM mass.

preprint2001arXiv

Scalar Field Dark Matter

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 $Φ$ with the scalar potential $V(Φ)=V_{0}(cosh {(λ\sqrt{κ_{0}}Φ)}-1)$ and the dark energy is modeled by a scalar field $Ψ$, endowed with the scalar potential $\tilde{V}(Ψ)= \tilde{V_{0}}(\sinh{(α\sqrt{κ_{0}}Ψ)})^β$, 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.