Search

In this work, we study the impact of anisotropy on slowly rotating neutron stars by extending the Hartle-Thorne formalism in general relativity to include anisotropy in pressure up to second order in the angular velocity. We assess the presence of anisotropy within the star by employing a quasi-local relationship. Our results show that the ratio between the gravitational mass of the fastest anisotropic rotating configurations and the corresponding non-rotating ones ranges from $1.12$ to $1.25$, consistent with recent findings. We develop universal relations for the moment of inertia, binding energy, and quadrupole moment of the rotating stars. These relations are tested against various equations of state, which were modeled by a piecewise polytropic function with continuous sound speed., Comment: 12 pages and 6 figures, accepted in Phys. Rev. D

In this paper, we study the impact of anisotropy on neutron stars with different equations of state, which have been modeled by a piecewise polytropic function with continuous sound speed. Anisotropic pressure in neutron stars is often attributed to interior magnetic fields, rotation, and the presence of exotic matter or condensates. We quantify the presence of anisotropy within the star by assuming a quasi-local relationship. We find that the radial and tangential sound velocities constrain the range of anisotropy allowed within the star. As expected, the anisotropy affects the macroscopic properties of stars, and it can be introduced to reconcile them with astrophysical observations. For instance, the maximum mass of anisotropic neutron stars can be increased by up to 15\% compared to the maximum mass of the corresponding isotropic configuration. This allows neutron stars to reach masses greater than $2.5M_\odot$, which may explain the secondary compact object of the GW190814 event. Additionally, we propose a universal relation for the binding energy of an anisotropic neutron star as a function of the star's compactness and the degree of anisotropy., Comment: 12 pages, 7 figures

In this work, we investigate whether the compact object at the center of the Milky Way is a naked singularity described by the $\textit{q}$-metric spacetime. Our fitting of the astrometric and spectroscopic data for the S2 star implies that similarly to the Schwarzschild black hole, the $\textit{q}$-metric naked singularity offers a satisfactory fit to the observed measurements. Additionally, it is shown that the shadow produced by the naked singularity is consistent with the shadow observed by the Event Horizon Telescope collaboration for Sgr-A*. It is worth mentioning that the spatial distribution of the S-stars favors the notion that the compact object at the center of our Galaxy can be described by an almost static spacetime. Based on these findings, the $\textit{q}$-metric naked singularity turns up as a compelling candidate for further investigation., Comment: Accepted for publication in Classical and Quantum Gravity

From a theoretical perspective, matter accretion processes around compact objects are highly relevant as they serve as a natural laboratory to test general relativity in the strong field regime. This enables us to validate fundamental concepts such as the no-hair theorem, the cosmic censorship hypothesis, and the existence of alternative solutions to Einstein's equations that mimic the effects of black holes. In this study, we analyze the emission spectra of geometrically thick accretion disks, referred to as Polish doughnuts, around naked singularities described by the $q$-metric. To begin, we revisit the construction of equilibrium configurations of magnetized tori in this spacetime and evaluate the role of the deformation parameter over these configurations. Once we have systematically studied the disks in this spacetime, we use the \texttt{OSIRIS} code to perform a backward ray-tracing method, resulting in the first simulations of the intensity map and emission profiles of magnetized tori within this metric. Furthermore, we validate the effect of both the quadrupole moment and the angular momentum on observable quantities such as flux and intensity for optically thin and thick disks, since for values of $ q < 0$, which correspond to objects with prolate deformation, and which in turn, are constructed with higher values of angular momentum, the emission spectrum exhibits higher intensity than that obtained for Schwarzschild's spacetime. Hence, we find a first differential feature that distinguishes tori formed around naked singularities from those around static black holes., Comment: 27 pages, 9 figures. Accepted for publication in Classical and Quantum Gravity

The radiation observed in quasars and active galactic nuclei is mainly produced by a relativistic plasma orbiting close to the black hole event horizon, where strong gravitational effects are relevant. The observational data of such systems can be compared with theoretical models to infer the black hole and plasma properties. In the comparison process, ray tracing algorithms are essential to computing the trajectories followed by the photons from the source to our telescopes. In this paper, we present OSIRIS: a new stable FORTRAN code capable of efficiently computing null geodesics around compact objects, including general relativistic effects such as gravitational lensing, redshift, and relativistic boosting. The algorithm is based on the Hamiltonian formulation and uses different integration schemes to evolve null geodesics while tracking the error in the Hamiltonian constrain to ensure physical results. We found from an error analysis that the integration schemes are all stable, and the best one maintains an error below $10^{-11}$. Particularly, to test the robustness and ability of the code to evolve geodesics in curved spacetime, we compute the shadow and Einstein rings of a Kerr black hole with different rotation parameters and obtain the image of a thin Keplerian accretion disk around a Schwarzschild black hole. Although OSIRIS is parallelized neither with MPI nor with CUDA, the computation times are of the same order as those reported by other codes with these types of parallel computing platforms., Comment: 13 pages, 13 figures, accepted for publication in The European Physical Journal C

The magneto-rotational instability (MRI) is the most likely mechanism for transportation of angular momentum and dissipation of energy within hot, ionized accretion discs. This instability is produced through the interactions of a differentially rotating plasma with an embedded magnetic field. Like all substances in nature, the plasma in an accretion disc has the potential to become magnetically polarized when it interacts with the magnetic field. In this paper, we study the effect of this magnetic susceptibility, parameterized by $\chi_m$, on the MRI, specifically within the context of black hole accretion. We find from a linear analysis within the Newtonian limit that the minimum wavelength of the first unstable mode and the wavelength of the fastest growing mode are shorter in paramagnetic ($\chi_m>0$) than in diamagnetic ($\chi_m<0$) discs, all other parameters being equal. Furthermore, the magnetization parameter (ratio of gas to magnetic pressure) in the saturated state should be smaller when the magnetic susceptibility is positive than when it is negative. We confirm this latter prediction through a set of numerical simulations of magnetically polarized black hole accretion discs. We additionally find that the vertically integrated stress and mass accretion rate are somewhat larger when the disc is paramagnetic than when it is diamagnetic. If astrophysical discs are able to become magnetically polarized to any significant degree, then our results would be relevant to properly interpreting observations., Comment: 12 pages, 8 figures, accepted for publication in MNRAS

One crucial problem in relativistic astrophysics is that of the nature of black hole candidates. It is usually assumed that astrophysical black holes are described by the Schwarzschild or Kerr space-times; however, there is no direct evidence to assert this. Moreover, there are various solutions in general relativity that can be alternatives to black holes, usually called black hole mimickers. In this work, we study the shadow produced by a compact object described by the q-metric, which is the simplest static and axially symmetric solution of Einstein equations with a non-vanishing quadrupole moment. This particular spacetime has the property of containing an independent parameter $q$, which is related to the compact object deformation. The solution corresponds to naked singularities for some specific values of this parameter. Additionally, we analyze the eigenvalues of the Riemann tensor using the $SO(3,C)$ representation, which allows us to find, in an invariant way, regions where there may be repulsive effects. Furthermore, we numerically solve the motion equations to show the shadow, the Einstein ring, and the gravitational lensing to establish a possible signature of such repulsive effects. We found that as $q$ is smaller, the Einstein ring decreases, but the shape is the same as the Schwarzschild black hole case. However, for values of $q$ lower or equal than $-0.5$, repulsive gravitational effects appear in the gravitational lensing close to the compact object, where a strong dependence of the system to the initial conditions seems to take place., Comment: Pages 20, 32 figures. Accepted for publication in Classical and Quantum Gravity (CQG)

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$., Comment: 7 pages, 4 figures and 1 table

Recently it has been discovered that on a stone tablet over 3800 years old, the Plimpton-322 table, are carved the geometric relations that exist between the sides of 15 right triangles chosen in a very special way. Due to its property as a super accuracy calculation tool, in this work we have called it sPad by stone pad, and we have calculated its machine accuracy $\epsilon_m$. Additionally, we present the physical and astrophysical constants most used in science and engineering in sexagesimal base. ----- Reci\'entemente se ha descubierto que en una tableta de piedra de m\'as de 3800 a\~nos de antig\"uedad, la Tabla Plimpton 322, est\'an ta\-lla\-das las relaciones geom\'etricas que existen entre los lados de 15 tri\'angulos rect\'angulos escogidos de manera muy especial. Debido a su propiedad como herramienta de c\'alculo s\'uper preciso, en este trabajo la hemos llamado sPad por \emph{stone pad}, y hemos calculado su exactitud de m\'aquina $\epsilon_{m}$. Adicionalmente presentamos las constantes f\'isicas y astrof\'isicas m\'as usadas en ciencia e ingenier\'ia en base sexagesimal., Comment: in Spanish

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., Comment: 12 pages, 30 figures, 1 table. Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

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., Comment: 6 pages, 14 png figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

A deep exploration of the parameter space that relates the interacting equation of state with the bag constant B, and the interaction parameter a, is fundamental for the construction of diverse models of quark stars. In particular, the anisotropy of quark stars with a well-motivated quantum chromodynamics (QCD) equation of state is presented here. The contribution of the fourth order corrections parameter ($\mathrm{a}$) of the QCD perturbation on the radial and tangential pressure generate significant effects on the mass-radius relation and the stability of the quark star. An adequate set of solutions for several values of the bag factor and the interaction parameter are used in order to calculate the relation between the mass, radius, density, compactness, and consequently the maximum masses and the stability. Therefore, while the more interactive quark solution leads to higher masses, the weak interaction among quarks give solutions similar to the widely known MIT bag model., Comment: All comments and suggestions are welcome

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., Comment: 7 pages, 15 .pdf figures, 1 table, Accepted for publication in Astronomy & Astrophysics (A&A)

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., Comment: 18 pages, 11 .pdf figures, Accepted for publication in General Relativity and Gravitation

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., Comment: 20 pages, 53 .pdf figures, 1 table, Accepted for publication in The Astrophysical Journal

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., Comment: 8 pages, 13 figures, 1 table. Accepted for publication in Monthly Notices of the Royal Astronomical Society

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., Comment: 15 Pages, 33 eps figures, 4 tables, Accepted for publication in Astrophysical Journal Supplement Series

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., Comment: 7 pages, 4 figures, 70 and 70 Gravitation Fest, 28 September 2016, Cartagena, Colombia

Nowadays, the magnetic and radiation fields are very important to understand the matter accretion into compact objects, the dynamics of binary systems, the equilibrium configurations of neutron stars, the photon diffusion, etc. The energy and the momentum associated to these fields, along with the matter one, need to satisfy some conditions that guarantee an appropriate physical behavior of the source and its gravitational field. Based on this fact, we present the energy conditions for a perfect fluid with magnetic and radiation field, in which the radiation part of the energy-momentum tensor is assumed to be approximately isotropic, in accordance with the optically thick regime. In order to find these conditions, the stress tensor of the system is written in an orthonormal basis in which it becomes diagonal, and the energy conditions are computed through contractions of the energy-momentum tensor with the four velocity vector of an arbitrary observer. Finally, the conditions for a magnetized fluid are presented as a particular case in which the radiation contribution is zero., Comment: 12 pages

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., Comment: 13 pages, 6 figures. This version has been substantially revised. Some new results were added

We study the dynamics of the planar circular restricted three-body problem in the context of a pseudo-Newtonian approximation. By using the Fodor-Hoenselaers-Perj\'es procedure, we perform an expansion in the mass potential of a static massive spherical source up to the first non-Newtonian term, giving place to a gravitational potential that includes first-order general relativistic effects. With this result, we model a system composed by two pseudo-Newtonian primaries describing circular orbits around their common center of mass, and a test particle orbiting the system in the equatorial plane. The dynamics of the new system of equations is studied in terms of the Poincar\'e section method and the Lyapunov exponents, where the introduction of a new parameter $\epsilon$, allows us to observe the transition from the Newtonian to the pseudo-Newtonian regime. We show that when the Jacobian constant is fixed, a chaotic orbit in the Newtonian regime can be either chaotic or regular in the pseudo-Newtonian approach. As a general result, we find that most of the pseudo-Newtonian configurations are less stable than their Newtonian equivalent., Comment: 11 pages, 2 figures. Accepted for publication in Physics Letters A, In Press

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., Comment: 11 pages, 20 figures, 1 table. Accepted for publication in MNRAS

Considering the growing interest of the astrophysicist community in the study of dissipative fluids with the aim of getting a more realistic description of the universe, we present in this paper a physical analysis of the energy-momentum tensor of a viscous fluid with heat flux. We introduce the general form of this tensor and, using the approximation of small velocity gradients, we relate the stresses of the fluid with the viscosity coefficients, the shear tensor and the expansion factor. Exploiting these relations, we can write the stresses in terms of the extrinsic curvature of the normal surface to the 4-velocity vector of the fluid, and we can also establish a connection between the perfect fluid and the symmetries of the spacetime. On the other hand, we calculate the energy conditions for a dissipative fluid through contractions of the energy-momentum tensor with the 4-velocity vector of an arbitrary observer. This method is interesting because it allows us to compute the conditions in a reasonable easy way and without considering any approximation or restriction on the energy-momentum tensor., Comment: Accepted for publication in General Relativity and Gravitation, 8 pages, Minor changes

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., Comment: 5 pages, 2 figures

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., Comment: 17 pages, 48 png fiugres,3 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

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., Comment: 14 pages, 60 png figures, 1 table. Accepted for publication in the Astrophysical Journal Supplement

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., Comment: 7 pages, 10 eps figures, 1 tables. Accepted for publication in General Relativity and Gravitation

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., Comment: 30 pages, 110 png figures, 4 tables. Accepted for publication in the Astrophysical Journal Supplement. More numerical details, tests and additional references

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., Comment: 11 pages, 6 eps figures, 3 tables. Accepted for publication in MNRAS

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., Comment: 15 pages, 6 eps figures. Accepted for publication in JCAP

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., Comment: 10 revtex pages, 30 eps figures

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., Comment: 11 pages, 8 eps figures. Accepted for publication in JCAP

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., Comment: 8 pages, 5 eps figures, Prepared for the conference proceedings of the IX Mexican School on Gravitation and Mathematical Physics: Cosmology for the XXI Century

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., Comment: Prepared with educative purposes, 15 pages, 34 eps figures. Accepted for publication in the educative section of the Mexican Journal of Physics

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., Comment: 12 pages, 34 eps figures, accepted for publication in MNRAS

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., Comment: Prepared with educative purposes, 15 pages, 34 eps figures

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., Comment: 8 pages, 9 eps figures, accepted for publication in MNRAS

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$., Comment: 8 pages, 16 eps figures, Proceedings of the VIII Mexican School on Gravitation and Mathematical Physics. Gravitational physics: testing gravity from submillimeter to cosmic

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., Comment: 7 revtex pages, accepted for publication in Phys. Rev. D

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., Comment: 7 pages, 8 eps figures, accepted for publication in MNRAS

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., Comment: 9 pages, 24 eps figures, 2 tables. Accepted for publication in MNRAS

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., Comment: 17 revtex pages, 8 eps figures

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, Comment: 15 pages, 21 eps figures. Appendix added, accepted for publication in JCAP

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., Comment: 16 revtex4 pages, 32 eps figures, may work as an educative paper

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., Comment: 9 pages, 15 eps figures, revtex4