Source author record

F. D. Lora-Clavijo

F. D. Lora-Clavijo 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

28works
10topics
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

28 published item(s)

preprint2022arXiv

OSIRIS: A New Code for Ray Tracing Around Compact Objects

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.

preprint2020arXiv

Neutron and quark stars: constraining the parameters for simple EoS using the GW170817

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 ρ_{0}^Γ$, and when it depends on the energy density $p=K ρ^Γ$, 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$.

preprint2020arXiv

sPad of Super Accuracy and Geometry in Old Babylon. -- sPad de Súper Exactitud y Geometría en la Antigua Babilonia

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 $ε_m$. Additionally, we present the physical and astrophysical constants most used in science and engineering in sexagesimal base. ----- Reciéntemente se ha descubierto que en una tableta de piedra de más de 3800 años de antigüedad, la Tabla Plimpton 322, están ta\-lla\-das las relaciones geométricas que existen entre los lados de 15 triángulos rectángulos escogidos de manera muy especial. Debido a su propiedad como herramienta de cálculo súper preciso, en este trabajo la hemos llamado sPad por \emph{stone pad}, y hemos calculado su exactitud de máquina $ε_{m}$. Adicionalmente presentamos las constantes físicas y astrofísicas más usadas en ciencia e ingeniería en base sexagesimal.

preprint2016arXiv

An eclipsing binary black hole candidate system in the blazar Mrk 421

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.

preprint2016arXiv

Non Axisymmetric Relativistic Wind Accretion with Velocity Gradients onto a Rotating Black Hole

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 $ε_{v}$, the black hole spin $a$, the asymptotic Mach number ${\cal M}_{\infty}$ and adiabatic index $Γ$. 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.

preprint2016arXiv

Pseudo-Newtonian planar circular restricted 3-body problem

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és 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é section method and the Lyapunov exponents, where the introduction of a new parameter $ε$, 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.

preprint2016arXiv

The Energy-Momentum Tensor for a Dissipative Fluid in General Relativity

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.

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

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

Relativistic Bondi-Hoyle-Lyttleton Accretion onto a Rotating Black-Hole: Density Gradients

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

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

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

Charged Annular Disks and Reissner-Nordström Type Black Holes from Extremal Dust

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öm black hole can be embedded into the center of the disks by adding a boundary term in the inversion.

preprint2010arXiv

Evolution of a mass-less test scalar field on Boson Stars space-times

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.

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.