BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Olivier Sarbach

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

ResearcherUnclaimed source record
gr-qchep-thmath-phmath.MPcond-mat.stat-mechastro-ph.COastro-ph.SRastro-ph.HEphysics.comp-phphysics.data-anphysics.hist-phphysics.ins-det

Catalog footprint

What is connected

38works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

38 published item(s)

preprint2022arXiv

An introduction to the relativistic kinetic theory on curved spacetimes

This article provides a self-contained pedagogical introduction to the relativistic kinetic theory of a dilute gas propagating on a curved spacetime manifold (M,g) of arbitrary dimension. Special emphasis is made on geometric aspects of the theory in order to achieve a formulation which is manifestly covariant on the relativistic phase space. Whereas most previous work has focussed on the tangent bundle formulation, here we work on the cotangent bundle associated with (M,g) which is more naturally adapted to the Hamiltonian framework of the theory. In the first part of this work we discuss the relevant geometric structures of the cotangent bundle T*M, starting with the natural symplectic form on T*M, the one-particle Hamiltonian and the Liouville vector field, defined as the corresponding Hamiltonian vector field. Next, we discuss the Sasaki metric on T*M and its most important properties, including the role it plays for the physical interpretation of the one-particle distribution function. In the second part of this work we describe the general relativistic theory of a collisionless gas, starting with the derivation of the collisionless Boltzmann equation for a neutral simple gas. Subsequently, the description is generalized to a charged gas consisting of several species of particles and the general relativistic Vlasov-Maxwell equations are derived for this system. The last part of this work is devoted to a transparent derivation of the collision term, leading to the general relativistic Boltzmann equation on (M,g). The meaning of global and local equilibrium and the stringent restrictions for the existence of the former on a curved spacetime are discussed. We close this article with an application of our formalism to the expansion of a homogeneous and isotropic universe filled with a collisional simple gas and its behavior in the early and late epochs. [abbreviated version]

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extreme $\ell$-boson stars

A new class of complex scalar field objects, which generalize the well known boson stars, was recently found as solutions to the Einstein-Klein-Gordon system. The generalization consists in incorporating some of the effects of angular momentum, while still maintaining the spacetime's spherical symmetry. These new solutions depend on an (integer) angular parameter $\ell$, and hence were named $\ell$-boson stars. Like the standard $\ell=0$ boson stars these configurations admit a stable branch in the solution space; however, contrary to them they have a morphology that presents a shell-like structure with a "hole" in the internal region. In this article we perform a thorough exploration of the parameter space, concentrating particularly on the extreme cases with large values of $\ell$. We show that the shells grow in size with the angular parameter, doing so linearly for large values, with the size growing faster than the thickness. Their mass also increases with $\ell$, but in such a way that their compactness, while also growing monotonically, converges to a finite value corresponding to about one half of the Buchdahl limit for stable configurations. Furthermore, we show that $\ell$-boson stars can be highly anisotropic, with the radial pressure diminishing relative to the tangential pressure for large $\ell$, reducing asymptotically to zero, and with the maximum density also approaching zero. We show that these properties can be understood by analyzing the asymptotic limit $\ell\rightarrow\infty$ of the field equations and their solutions. We also analyze the existence and characteristics of both timelike and null circular orbits, especially for very compact solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Kinetic gas disks surrounding Schwarzschild black holes

We describe stationary and axisymmetric gas configurations surrounding black holes. They consist of a collisionless relativistic kinetic gas of identical massive particles following bound orbits in a Schwarzschild exterior spacetime and are modeled by a one-particle distribution function which is the product of a function of the energy and a function of the orbital inclination associated with the particle's trajectory. The morphology of the resulting configuration is analyzed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Choked accretion onto a Kerr black hole

The choked accretion model consists of a purely hydrodynamical mechanism in which, by setting an equatorial to polar density contrast, a spherically symmetric accretion flow transitions to an inflow-outflow configuration. This scenario has been studied in the case of a (non-rotating) Schwarzschild black hole as central accretor, as well as in the non-relativistic limit. In this article, we generalize these previous works by studying the accretion of a perfect fluid onto a (rotating) Kerr black hole. We first describe the mechanism by using a steady-state, irrotational analytic solution of an ultrarelativistic perfect fluid, obeying a stiff equation of state. We then use hydrodynamical numerical simulations in order to explore a more general equation of state. Analyzing the effects of the black hole's rotation on the flow, we find in particular that the choked accretion inflow-outflow morphology prevails for all possible values of the black hole's spin parameter, showing the robustness of the model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gauge-invariant spherical linear perturbations of wormholes in Einstein gravity minimally coupled to a self-interacting phantom scalar field

Recently, there has been quite a lot of interest in static, spherical wormhole spacetimes and the question of their stability with respect to time-dependent perturbations. The consideration of linearized perturbations usually leads to a master wave equation with effective potential which can then be analyzed using standard tools from quantum mechanics. However, in the wormhole case, particular care must be taken with the gauge conditions when formulating the master equation. A poor coordinate choice, based for example on fixing the areal radial coordinate, may lead to singularities at the throat which complicate the stability analysis or might even lead to erroneous conclusions regarding the stability of the underlying wormhole configuration. In this work, we present a general method for deriving a gauge-invariant wave system of linearized perturbation equations in the spherically symmetric case, assuming that the matter supporting the wormhole is a phantom scalar field, that is, a self-interacting scalar field whose kinetic energy has the reversed sign. We show how this system can be decoupled and reduced to a single master wave equation with a regular potential, with no intermediate steps involving singularities at the throat. Two applications of our formalism are given. First, we rederive the master equation for the linearly perturbed Ellis-Bronnikov wormhole using our new, singularity-free method. Second, we derive the master equation describing the linear perturbations of a certain Anti de Sitter wormhole, provide a detailed analysis of the spectral properties of the underlying operator and prove that this wormhole is linearly unstable. In the final part of the paper, we consider a wormhole with de Sitter-type ends, whose spacetime presents horizons and admits a nonstatic extension beyond them; for this system we derive partial results of linear instability.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Phase space mixing in an external gravitational central potential

This article is devoted to the study of the dynamical behavior of a collisionless kinetic gas in d=1,2,3 space dimensions which is trapped in a rotationally symmetric potential well. Although at the microscopic level the trajectories of individual gas particles are quasi-periodic and characterized by their d fundamental frequencies, at the macroscopic level the gas relaxes in time to a stationary state, provided the potential satisfies a certain non-degeneracy condition. In this article, we provide a mathematically precise formulation for this relaxation process which is due to phase space mixing. In particular, we prove that a physically relevant class of macroscopic observables computed from the one-particle distribution function, such as particle and energy densities, pressure and stress tensors, converge in time to the corresponding observables associated with an averaged distribution function. The latter can be determined from the initial datum and depends only on integrals of motion. Thus, the final state of the gas is described by an effective distribution function depending only on integrals of motion, which considerably reduces the degrees of freedom of the gas configuration. We discuss some applications to gravitational physics, including the propagation of a collisionless gas in typical potentials arising in stellar dynamics and the modeling of dark matter halos, and we also generalize our results to a relativistic gas whose individual particles follow bound timelike trajectories in the exterior region of a static, spherically symmetric black hole spacetime.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Dynamical evolutions of $\ell$-boson stars in spherical symmetry

In previous work, we have found new static, spherically symmetric boson star solutions which generalize the standard boson stars by allowing a particular superposition of scalar fields in which each of the fields is characterized by a fixed value of its non-vanishing angular momentum number $\ell$. We call such solutions "$\ell$-boson stars". Here, we perform a series of fully non-linear dynamical simulations of perturbed $\ell$-boson stars in order to study their stability, and the final fate of unstable configurations. We show that for each value of $\ell$, the configuration of maximum mass separates the parameter space into stable and unstable regions. Stable configurations, when perturbed, oscillate around the unperturbed solution and very slowly return to a stationary configuration. Unstable configurations, in contrast, can have three different final states: collapse to a black hole, migration to the stable branch, or explosion (dissipation) to infinity. Just as it happens with $\ell=0$ boson stars, migration to the stable branch or dissipation to infinity depends on the sign of the total binding energy of the star: bound unstable stars collapse to black holes or migrate to the stable branch, whereas unbound unstable stars either collapse to a black hole or explode to infinity. Thus, the parameter $\ell$ allows us to construct a new set of stable configurations. All our simulations are performed in spherical symmetry, leaving a more detailed stability analysis including non-spherical perturbations for future work.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Black holes in nonlinear electrodynamics: quasi-normal spectra and parity splitting

We discuss the quasi-normal oscillations of black holes which are sourced by a nonlinear electrodynamic field. While previous studies have focused on the computation of quasi-normal frequencies for the wave or higher spin equation on a fixed background geometry described by such black holes, here we compute for the first time the quasi-normal frequencies for the coupled electromagnetic-gravitational linear perturbations. To this purpose, we consider a parametrized family of Lagrangians for the electromagnetic field which contains the Maxwell Lagrangian as a special case. In the Maxwell case, the unique spherically symmetric black hole solutions are described by the Reissner-Nordström family and in this case it is well-known that the quasi-normal spectra in the even- and odd-parity sectors are identical to each other. However, when moving away from the Maxwell case, we obtain deformed Reissner-Nordström black holes, and we show that in this case there is a parity splitting in the quasi-normal mode spectra. A partial explanation for this phenomena is provided by considering the eikonal (high-frequency) limit.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Emergence of Gravitational Wave Science: 100 Years of Development of Mathematical Theory, Detectors, Numerical Algorithms, and Data Analysis Tools

On September 14, 2015, the newly upgraded Laser Interferometer Gravitational-wave Observatory (LIGO) recorded a loud gravitational-wave (GW) signal, emitted a billion light-years away by a coalescing binary of two stellar-mass black holes. The detection was announced in February 2016, in time for the hundredth anniversary of Einstein's prediction of GWs within the theory of general relativity (GR). The signal represents the first direct detection of GWs, the first observation of a black-hole binary, and the first test of GR in its strong-field, high-velocity, nonlinear regime. In the remainder of its first observing run, LIGO observed two more signals from black-hole binaries, one moderately loud, another at the boundary of statistical significance. The detections mark the end of a decades-long quest, and the beginning of GW astronomy: finally, we are able to probe the unseen, electromagnetically dark Universe by listening to it. In this article, we present a short historical overview of GW science: this young discipline combines GR, arguably the crowning achievement of classical physics, with record-setting, ultra-low-noise laser interferometry, and with some of the most powerful developments in the theory of differential geometry, partial differential equations, high-performance computation, numerical analysis, signal processing, statistical inference, and data science. Our emphasis is on the synergy between these disciplines, and how mathematics, broadly understood, has historically played, and continues to play, a crucial role in the development of GW science. We focus on black holes, which are very pure mathematical solutions of Einstein's gravitational-field equations that are nevertheless realized in Nature, and that provided the first observed signals.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Michel accretion of a polytropic fluid with adiabatic index gamma > 5/3: Global flows versus homoclinic orbits

We analyze the properties of a polytropic fluid which is radially accreted into a Schwarzschild black hole. The case where the adiabatic index gamma lies in the range 1 < gamma <= 5/3 has been treated in previous work. In this article we analyze the complementary range 5/3 < gamma <= 2. To this purpose, the problem is cast into an appropriate Hamiltonian dynamical system whose phase flow is analyzed. While for 1 < gamma <= 5/3 the solutions are always characterized by the presence of a unique critical saddle point, we show that when 5/3 < gamma <= 2, an additional critical point might appear which is a center point. For the parametrization used in this paper we prove that whenever this additional critical point appears, there is a homoclinic orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Observational distinction between black holes and naked singularities: the role of the redshift function

We suggest that the redshift of photons traveling from past to future null infinity through a collapsing object could provide an observational signature capable of differentiating between the formation of a globally naked singularity and the formation of an event horizon. Supporting evidence for this idea is drawn from the analysis of photons with zero angular momentum through the center of a collapsing spherical dust cloud. We show that the frequency shift as a function of proper time with respect to stationary observers has distinct features depending on whether the object collapses to a black hole or a naked singularity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quasi-normal acoustic oscillations in the Michel flow

We study spherical and nonspherical linear acoustic perturbations of the Michel flow, which describes the steady radial accretion of a perfect fluid into a nonrotating black hole. The dynamics of such perturbations are governed by a scalar wave equation on an effective curved background geometry determined by the acoustic metric, which is constructed from the spacetime metric and the particle density and four-velocity of the fluid. For the problem under consideration in this article the acoustic metric has the same qualitative features as an asymptotically flat, static and spherically symmetric black hole, and thus it represents a natural astrophysical analogue black hole. As for the case of a scalar field propagating on a Schwarzschild background, we show that acoustic perturbations of the Michel flow exhibit quasi-normal oscillations. Based on a new numerical method for determining the solutions of the radial mode equation, we compute the associated frequencies and analyze their dependency on the radii of the event and sonic horizons and the angular momentum number. Our results for the fundamental frequencies are compared to those obtained from an independent numerical Cauchy evolution, finding good agreement between the two approaches. When the radius of the sonic horizon is large compared to the event horizon radius, we find that the quasi-normal frequencies scale approximately like the surface gravity associated with the sonic horizon.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quasi-normal acoustic oscillations in the transonic Bondi flow

In recent work, we analyzed the dynamics of spherical and nonspherical acoustic perturbations of the Michel flow, describing the steady radial accretion of a relativistic perfect fluid into a nonrotating black hole. We showed that such perturbations undergo quasi-normal oscillations and computed the corresponding complex frequencies as a function of the black hole mass M and the radius r_c of the sonic horizon. It was found that when r_c is much larger than the Schwarzschild radius r_H = 2GM/c^2 of the black hole, these frequencies scale like the surface gravity of the analogue black hole associated with the acoustic metric. In this work, we analyze the Newtonian limit of the Michel solution and its acoustic perturbations. In this limit, the flow outside the sonic horizon reduces to the transonic Bondi flow, and the acoustic metric reduces to the one introduced by Unruh in the context of experimental black hole evaporation. We show that for the transonic Bondi flow, Unruh's acoustic metric describes an analogue black hole and compute the associated quasi-normal frequencies. We prove that they do indeed scale like the surface gravity of the acoustic black hole, thus providing an explanation for our previous results in the relativistic setting.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Radial accretion flows on static spherically symmetric black holes

We analyze the steady radial accretion of matter into a nonrotating black hole. Neglecting the self-gravity of the accreting matter, we consider a rather general class of static, spherically symmetric and asymptotically flat background spacetimes with a regular horizon. In addition to the Schwarzschild metric, this class contains certain deformation of it which could arise in alternative gravity theories or from solutions of the classical Einstein equations in the presence of external matter fields. Modeling the ambient matter surrounding the black hole by a relativistic perfect fluid, we reformulate the accretion problem as a dynamical system, and under rather general assumptions on the fluid equation of state, we determine the local and global qualitative behavior of its phase flow. Based on our analysis and generalizing previous work by Michel, we prove that for any given positive particle density number at infinity, there exists a unique radial, steady-state accretion flow which is regular at the horizon. We determine the physical parameters of the flow, including its accretion and compression rates, and discuss their dependency on the background metric.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shadow of a naked singularity

We analyze the redshift suffered by photons originating from an external source, traversing a collapsing dust cloud and finally being received by an asymptotic observer. In addition, we study the shadow that the collapsing cloud casts on the sky of the asymptotic observer. We find that the resulting redshift and properties of the shadow depend crucially on whether the final outcome of the complete gravitational collapse is a black hole or a naked singularity. In the black hole case, the shadow is due to the high redshift acquired by the photons as they approach the event horizon, implying that their energy is gradually redshifted toward zero within a few crossing times associated with the event horizon radius. In contrast to this, a naked singularity not only absorbs photons originating from the source, but it also emits infinitely redshifted photons with and without angular momenta. This emission introduces an abrupt cutoff in the frequency shift of the photons detected in directions close to the radial one, and it is responsible for the shadow masking the source in the naked singularity case. Furthermore, even though the shadow forms and begins to grow immediately after the observer crosses the Cauchy horizon, it takes many more crossing times than in the black hole case for the source to be occulted from the observer's eyes. We discuss possible implications of our results for testing the weak cosmic censorship hypothesis. Even though at late times the image of the source perceived by the observer looks the same in both cases, the dynamical formation of the shadow and the redshift images has distinct features and time scales in the black hole versus the naked singularity case. For stellar collapse, these time scales seem to be too short to be resolved with existing technology. However, our results may be relevant for the collapse of seeds leading to supermassive black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Schwarzschild scalar wigs: spectral analysis and late time behavior

Using the Green's function representation technique, the late time behavior of localized scalar field distributions on Schwarzschild spacetimes is studied. Assuming arbitrary initial data we perform a spectral analysis, computing the amplitude of each excited quasi-bound mode without the necessity of performing dynamical evolutions. The resulting superposition of modes is compared with a traditional numerical evolution with excellent agreement; therefore, we have an efficient way to determine final black hole wigs. The astrophysical relevance of the quasi-bound modes is discussed in the context of scalar field dark matter models and the axiverse.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The evolution of a spatially homogeneous and isotropic universe filled with a collisionless gas

We review the evolution of a spatially homogeneous and isotropic universe described by a Friedmann-Robertson-Walker spacetime filled with a collisionless, neutral, simple, massive gas. The gas is described by a one-particle distribution function which satisfies the Liouville equation and is assumed to be homogeneous and isotropic. Making use of the isometries of the spacetime, we define precisely the homogeneity and isotropicity property of the distribution function, and based on this definition we give a concise derivation of the most general family of such distribution functions. For this family, we construct the particle current density and the stress-energy tensor and consider the coupled Einstein-Liouville system of equations. We find that as long as the distribution function is collisionless, homogenous and isotropic, the evolution of a Friedmann-Robertson-Walker universe exhibits a singular origin. Its future development depends upon the curvature of the spatial sections: spatially flat or hyperboloid universes expand forever and this expansion dilutes the energy density and pressure of the gas, while a universe with compact spherical sections reaches a maximal expansion, after which it reverses its motion and recollapses to a final crunch singularity where the energy density and isotropic pressure diverge. Finally, we analyze the evolution of the universe filled with the collisionless gas once a cosmological constant is included.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The geometry of the tangent bundle and the relativistic kinetic theory of gases

This article discusses the relativistic kinetic theory for a simple collisionless gas from a geometric perspective. We start by reviewing the rich geometrical structure of the tangent bundle TM of a given spacetime manifold, including the splitting of the tangent spaces of TM into horizontal and vertical subspaces and the natural metric and symplectic structure it induces on TM. Based on these structures we introduce the Liouville vector field L and a suitable Hamiltonian function H on TM. The Liouville vector field turns out to be the Hamiltonian vector field associated to H. On the other hand, H also defines the mass shells as Lorentzian submanifolds of the tangent bundle. A simple collisionless gas is described by a distribution function on a particular mass shell, satisfying the Liouville equation. Together with the Liouville vector field the distribution function can be thought of as defining a fictitious incompressible fluid on the mass shells, with associated conserved current density. Flux integrals of this current density provide the averaged properties of the gas, while suitable fibre integrals of the distribution function define divergence-free tensor fields on the spacetime manifold such as the current density and stress-energy tensor. Finally, we discuss the relationship between symmetries of the spacetime manifold and symmetries of the distribution function. As a first application of our formalism we derive the most general spherically symmetric distribution function on any spherically symmetric spacetime and write the Einstein-Liouville equations as effective field equations on the two-dimensional radial manifold. As a second application we derive the most general collisionless distribution function on a Kerr black hole spacetime background.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The gravitational redshift of photons traversing a collapsing dust cloud and observable consequences

We analyze the frequency shift of photons propagating on an asymptotically flat spacetime describing a collapsing, spherical dust cloud. We focus on the case where the interaction of the photons with the matter can be neglected. Under fairly general assumptions on the initial data characterizing the collapse, we show that photons with zero angular momentum which travel from past to future null infinity, crossing the collapsing cloud through its center, are always redshifted with respect to stationary observers. We compute this redshift as a function of proper time of a distant stationary observer and discuss its dependency on the mass distribution of the cloud. Possible implications of this redshift effect for weak cosmic censorship and light propagation in cosmological spacetimes are also briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cauchy horizon stability in a collapsing spherical dust cloud I: geometric optic approximation and spherically symmetric test fields

A spherical dust cloud which is initially at rest and which has a monotonously decaying density profile collapses and forms a shell-focussing singularity. Provided the density profile is not too flat, meaning that its second radial derivative is negative at the center, this singularity is visible to local, and sometimes even to global observers. According to the strong cosmic censorship conjecture, such naked singularities should be unstable under generic, nonspherical perturbations of the initial data or when more realistic matter models are considered. In an attempt to gain some understanding about this stability issue, in this work we initiate the analysis of a simpler but related problem. We discuss the stability of test fields propagating in the vicinity of the Cauchy horizon associated to the naked central singularity. We first study the high-frequency limit and show that the fields undergo a blueshift as they approach the Cauchy horizon. However, in contrast to what occurs at inner horizons of black holes, we show that the blueshift is uniformly bounded along incoming and outgoing null rays. Motivated by this boundedness result, we take a step beyond the geometric optic approximation and consider the Cauchy evolution of spherically symmetric test scalar fields. We prove that under reasonable conditions on the initial data a suitable rescaled field can be continuously extended to the Cauchy horizon. In particular, this result implies that the physical field is everywhere finite on the Cauchy horizon away from the central singularity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Linear perturbations of self-gravitating spherically symmetric configurations

We present a new covariant, gauge-invariant formalism describing linear metric perturbation fields on any spherically symmetric background in general relativity. The advantage of this formalism relies in the fact that it does not require a decomposition of the perturbations into spherical tensor harmonics. Furthermore, it does not assume the background to be vacuum, nor does it require its staticity. In the particular case of vacuum perturbations, we derive two master equations describing the propagation of arbitrary linear gravitational waves on a Schwarzschild black hole. When decomposed into spherical harmonics, they reduce to covariant generalizations of the well-known Regge-Wheeler and Zerilli equations. Next, we discuss the general case where the metric perturbations are coupled to matter fields and derive a new constrained wave system describing the propagation of three gauge-invariant scalars from which the complete metric perturbations can be reconstructed. We apply our formalism to the Einstein-Euler system, dividing the fluid perturbations into two parts. The first part, which decouples from the metric perturbations, obeys simple advection equations along the background flow and describes the propagation of the entropy and the vorticity. The second part describes a perturbed potential flow, and together with the metric perturbations it forms a closed wave system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relativistic Kinetic Theory: An Introduction

We present a brief introduction to the relativistic kinetic theory of gases with emphasis on the underlying geometric and Hamiltonian structure of the theory. Our formalism starts with a discussion on the tangent bundle of a Lorentzian manifold of arbitrary dimension. Next, we introduce the Poincare one-form on this bundle, from which the symplectic form and a volume form are constructed. Then, we define an appropriate Hamiltonian on the bundle which, together with the symplectic form yields the Liouville vector field. The corresponding flow, when projected onto the base manifold, generates geodesic motion. Whenever the flow is restricted to energy surfaces corresponding to a negative value of the Hamiltonian, its projection describes a family of future-directed timelike geodesics. A collisionless gas is described by a distribution function on such an energy surface, satisfying the Liouville equation. Fibre integrals of the distribution function determine the particle current density and the stress-energy tensor. We show that the stress-energy tensor satisfies the familiar energy conditions and that both the current and stress-energy tensor are divergence-free. Our discussion also includes the generalization to charged gases, a summary of the Einstein-Maxwell-Vlasov system in any dimensions, as well as a brief introduction to the general relativistic Boltzmann equation for a simple gas.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Tangent bundle formulation of a charged gas

We discuss the relativistic kinetic theory for a simple, collisionless, charged gas propagating on an arbitrary curved spacetime geometry. Our general relativistic treatment is formulated on the tangent bundle of the spacetime manifold and takes advantage of its rich geometric structure. In particular, we point out the existence of a natural metric on the tangent bundle and illustrate its role for the development of the relativistic kinetic theory. This metric, combined with the electromagnetic field of the spacetime, yields an appropriate symplectic form on the tangent bundle. The Liouville vector field arises as the Hamiltonian vector field of a natural Hamiltonian. The latter also defines natural energy surfaces, called mass shells, which turn out to be smooth Lorentzian submanifolds. A simple, collisionless, charged gas is described by a distribution function which is defined on the mass shell and satisfies the Liouville equation. Suitable fibre integrals of the distribution function define observable fields on the spacetime manifold, such as the current density and stress-energy tensor. Finally, the geometric setting of this work allows us to discuss the relationship between the symmetries of the electromagnetic field, those of the spacetime metric, and the symmetries of the distribution function. Taking advantage of these symmetries, we construct the most general solution of the Liouville equation an a Kerr-Newman black hole background.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Conformal diagrams for the gravitational collapse of a spherically symmetric dust cloud

We present an algorithm for the construction of conformal coordinates in the interior of a spherically symmetric, collapsing matter cloud in general relativity. This algorithm is based on the numerical integration of radial null geodesics. As an application we generate conformal diagrams for collapsing spherical dust clouds and analyze the causal structure of the resulting spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Continuum and Discrete Initial-Boundary-Value Problems and Einstein's Field Equations

Many evolution problems in physics are described by partial differential equations on an infinite domain; therefore, one is interested in the solutions to such problems for a given initial dataset. A prominent example is the binary black hole problem within Einstein's theory of gravitation, in which one computes the gravitational radiation emitted from the inspiral of the two black holes, merger and ringdown. Powerful mathematical tools can be used to establish qualitative statements about the solutions, such as their existence, uniqueness, continuous dependence on the initial data, or their asymptotic behavior over large time scales. However, one is often interested in computing the solution itself, and unless the partial differential equation is very simple, or the initial data possesses a high degree of symmetry, this computation requires approximation by numerical discretization. When solving such discrete problems on a machine, one is faced with a finite limit to computational resources, which leads to the replacement of the infinite continuum domain with a finite computer grid. This, in turn, leads to a discrete initial-boundary value problem. The hope is to recover, with high accuracy, the exact solution in the limit where the grid spacing converges to zero with the boundary being pushed to infinity. The goal of this article is to review some of the theory necessary to understand the continuum and discrete initial-boundary value problems arising from hyperbolic partial differential equations and to discuss its applications to numerical relativity; in particular, we present well-posed initial and initial-boundary value formulations of Einstein's equations, and we discuss multi-domain high-order finite difference and spectral methods to solve them.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical simulations with a first order BSSN formulation of Einstein's field equations

We present a new fully first order strongly hyperbolic representation of the BSSN formulation of Einstein's equations with optional constraint damping terms. We describe the characteristic fields of the system, discuss its hyperbolicity properties, and present two numerical implementations and simulations: one using finite differences, adaptive mesh refinement and in particular binary black holes, and another one using the discontinuous Galerkin method in spherical symmetry. The results of this paper constitute a first step in an effort to combine the robustness of BSSN evolutions with very high accuracy numerical techniques, such as spectral collocation multi-domain or discontinuous Galerkin methods.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Schwarzschild black holes can wear scalar wigs

We study the evolution of a massive scalar field surrounding a Schwarzschild black hole and find configurations that can survive for arbitrarily long times, provided the black hole or the scalar field mass is small enough. In particular, both ultra-light scalar field dark matter around supermassive black holes and axion-like scalar fields around primordial black holes can survive for cosmological times. Moreover, these results are quite generic, in the sense that fairly arbitrary initial data evolves, at late times, as a combination of those long-lived configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The propagation of particles and fields in wormhole geometries

We discuss several properties of static, spherically symmetric wormholes with particular emphasis on the behavior of causal geodesics and the propagation of linear fields. We show there always exist null geodesics which are trapped in a region close to the throat. Depending upon the detailed structure of the wormhole geometry, these trapped geodesics can be stable, unlike the case of the Schwarzschild black hole. We also show that test scalar fields propagating on such wormholes are stable. However, when a mixture of ghost and Klein-Gordon scalar fields is used as a source of the Einstein equations we prove that the resulting static, spherically symmetric wormhole configurations are linearly unstable.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Are black holes a serious threat to scalar field dark matter models?

Classical scalar fields have been proposed as possible candidates for the dark matter component of the universe. Given the fact that super-massive black holes seem to exist at the center of most galaxies, in order to be a viable candidate for the dark matter halo a scalar field configuration should be stable in the presence of a central black hole, or at least be able to survive for cosmological time-scales. In the present work we consider a scalar field as a test field on a Schwarzschild background, and study under which conditions one can obtain long-lived configurations. We present a detailed study of the Klein-Gordon equation in the Schwarzschild spacetime, both from an analytical and numerical point of view, and show that indeed there exist quasi-stationary solutions that can remain surrounding a black hole for large time-scales.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Conformal diagrams for the gravitational collapse of a spherical dust cloud

We present an algorithm for the construction of conformal coordinates in the interior of a spherically symmetric, collapsing matter cloud in general relativity. This algorithm is based on the numerical integration of the radial null geodesics and a local analysis of their behavior close to the singularity. As an application, we consider a collapsing spherical dust cloud, generate the corresponding conformal diagram and analyze the structure of the resulting singularity. A new bound on the initial data which guarantees that the singularity is visible from future null infinity is also obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Tetrad formalism for numerical relativity on conformally compactified constant mean curvature hypersurfaces

We present a new evolution system for Einstein's field equations which is based on tetrad fields and conformally compactified hyperboloidal spatial hypersurfaces which reach future null infinity. The boost freedom in the choice of the tetrad is fixed by requiring that its timelike leg be orthogonal to the foliation, which consists of constant mean curvature slices. The rotational freedom in the tetrad is fixed by the 3D Nester gauge. With these conditions, the field equations reduce naturally to a first-order constrained symmetric hyperbolic evolution system which is coupled to elliptic equations for the gauge variables. The conformally rescaled equations are given explicitly, and their regularity at future null infinity is discussed. Our formulation is potentially useful for high accuracy numerical modeling of gravitational radiation emitted by inspiraling and merging black hole binaries and other highly relativistic isolated systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Boundary conditions for the Baumgarte-Shapiro-Shibata-Nakamura formulation of Einstein's field equations

We discuss the initial-boundary value problem for the Baumgarte-Shapiro-Shibata-Nakamura evolution system of Einstein's field equations which has been used extensively in numerical simulations of binary black holes and neutron stars. We specify nine boundary conditions for this system with the following properties: (i) they impose the momentum constraint at the boundary, which is shown to preserve all the constraints throughout evolution, (ii) they approximately control the incoming gravitational degrees of freedom by specifying the Weyl scalar Psi_0 at the boundary, (iii) they control the gauge freedom by requiring a Neumann boundary condition for the lapse, by setting the normal component of the shift to zero, and by imposing a Sommerfeld-like condition on the tangential components of the shift, (iv) they are shown to yield a well-posed problem in the limit of weak gravity. Possible numerical applications of our results are also discussed briefly.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nonlinear instability of wormholes supported by exotic dust and a magnetic field

Recently, spherically symmetric, static wormholes supported by exotic dust and a radial magnetic field have been derived and argued to be stable with respect to linear radial fluctuations. In this report we point out that these wormholes are unstable due to the formation of shell-crossing singularities when the nonlinearities of the theory are taken into account.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Turduckening black holes: an analytical and computational study

We provide a detailed analysis of several aspects of the turduckening technique for evolving black holes. At the analytical level we study the constraint propagation for a general family of BSSN-type formulation of Einstein's field equations and identify under what conditions the turducken procedure is rigorously justified and under what conditions constraint violations will propagate to the outside of the black holes. We present high-resolution spherically symmetric studies which verify our analytical predictions. Then we present three dimensional simulations of single distorted black holes using different variations of the turduckening method and also the puncture method. We study the effect that these different methods have on the coordinate conditions, constraint violations, and extracted gravitational waves. We find that the waves agree up to small but non-vanishing differences, caused by escaping superluminal gauge modes. These differences become smaller with increasing detector location.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

A model problem for the initial-boundary value formulation of Einstein's field equations

In many numerical implementations of the Cauchy formulation of Einstein's field equations one encounters artificial boundaries which raises the issue of specifying boundary conditions. Such conditions have to be chosen carefully. In particular, they should be compatible with the constraints, yield a well posed initial-boundary value formulation and incorporate some physically desirable properties like, for instance, minimizing reflections of gravitational radiation. Motivated by the problem in General Relativity, we analyze a model problem, consisting of a formulation of Maxwell's equations on a spatially compact region of spacetime with timelike boundaries. The form in which the equations are written is such that their structure is very similar to the Einstein-Christoffel symmetric hyperbolic formulations of Einstein's field equations. For this model problem, we specify a family of Sommerfeld-type constraint-preserving boundary conditions and show that the resulting initial-boundary value formulations are well posed. We expect that these results can be generalized to the Einstein-Christoffel formulations of General Relativity, at least in the case of linearizations about a stationary background.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Recent analytical and numerical techniques applied to the Einstein equations

Combining deeper insight of Einstein's equations with sophisticated numerical techniques promises the ability to construct accurate numerical implementations of these equations. We illustrate this in two examples, the numerical evolution of ``bubble'' and single black hole spacetimes. The former is chosen to demonstrate how accurate numerical solutions can answer open questions and even reveal unexpected phenomena. The latter illustrates some of the difficulties encountered in three-dimensional black hole simulations, and presents some possible remedies.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Convergence and stability in numerical relativity

It is often the case in numerical relativity that schemes that are known to be convergent for well posed systems are used in evolutions of weakly hyperbolic (WH) formulations of Einstein's equations. Here we explicitly show that with several of the discretizations that have been used through out the years, this procedure leads to non-convergent schemes. That is, arbitrarily small initial errors are amplified without bound when resolution is increased, independently of the amount of numerical dissipation introduced. The lack of convergence introduced by this instability can be particularly subtle, in the sense that it can be missed by several convergence tests, especially in 3+1 dimensional codes. We propose tests and methods to analyze convergence that may help detect these situations.

0 citations0 reviewsNo code linkedOpen access