Search

In this paper we study the linear stability of selfinteracting boson stars in the nonrelativistic limit of the Einstein-Klein-Gordon theory. For this purpose, based on a combination of analytic and numerical methods, we determine the behavior of general linear perturbations around the stationary and spherically symmetric solutions of the Gross-Pitaevskii-Poisson system. In particular, we conclude that ground state configurations are linearly stable if the selfinteraction is repulsive, whereas there exist a state of maximum mass that divides the stable and the unstable branches in case the selfinteraction is attractive. Regarding the excited states, they are in general unstable under generic perturbations, although we identify a stability band in the first excited states of the repulsive theory. This result is independent of the mass of the scalar field and the details of the selfinteraction potential, and it is in contrast to the situation of vanishing selfinteraction, in which excited states are always unstable., Comment: 15+1 pages and 10 figures. Minor changes were made. It was accepted in journal Physical Review D

In previous work we analyzed the linear stability of non-relativistic $\ell$-boson stars with respect to radial modes and showed that ground state configurations are stable with respect to these modes, whereas excited states are unstable. In this work we extend the analysis to non-spherical linear mode perturbations. To this purpose, we expand the wave function in terms of tensor spherical harmonics which allows us to decouple the perturbation equations into a family of radial problems. By using a combination of analytic and numerical methods, we show that ground state configurations with $\ell > 1$ possess exponentially in time growing non-radial modes, whereas only oscillating modes are found for $\ell=0$ and $\ell=1$. This leads us to conjecture that nonrelativistic $\ell$-boson stars in their ground state are stable for $\ell=1$ as well as $\ell=0$, while ground state and excited configurations with $\ell > 1$ are unstable., Comment: 21 pages, 5 figures, 2 tables. Minor changes were made. It was accepted in journal Physical Review D

We study the dynamics of a collisionless kinetic gas whose particles follow future-directed timelike and spatially bound geodesics in the exterior of a sub-extremal Kerr black hole spacetime. Based on the use of generalized action-angle variables, we analyze the large time asymptotic behavior of macroscopic observables associated with the gas. We show that, as long as the fundamental frequencies of the system satisfy a suitable non-degeneracy condition, these macroscopic observables converge in time to the corresponding observables determined from an averaged distribution function. In particular, this implies that the final state is characterized by a distribution function which is invariant with respect to the full symmetry group of the system, that is, it is stationary, axisymmetric and Poisson-commutes with the integral of motion associated with the Carter constant. As a corollary of our result, we demonstrate the validity of the strong Jeans theorem in our setting, stating that the distribution function belonging to a stationary state must be a function which is independent of the generalized angle variables. An analogous theorem in which the assumption of stationarity is replaced with the requirement of invariance with respect to the Carter flow is also proven. Finally, we prove that the aforementioned non-degeneracy condition holds. This is achieved by providing suitable asymptotic expansions for the energy and Carter constant in terms of action variables for orbits having sufficiently large radii, and by exploiting the analytic dependency of the fundamental frequencies on the integrals of motion., Comment: 52 pages, 6 figures

We study the linear stability of nonrelativistic $\ell$-boson stars, describing static, spherically symmetric configurations of the Schr\"odinger-Poisson system with multiple wave functions having the same value of the angular momentum $\ell$. In this work we restrict our analysis to time-dependent perturbations of the radial profiles of the $2\ell+1$ wave functions, keeping their angular dependency fixed. Based on a combination of analytic and numerical methods, we find that for each $\ell$, the ground state is linearly stable, whereas the $n$'th excited states possess $2n$ unstable (exponentially in time growing) modes. Our results also indicate that all excited states correspond to saddle points of the conserved energy functional of the theory., Comment: 21+1 pages, 7 figures, 5 tables

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], Comment: 79 pages, 12 figures. Typos corrected, references added

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., Comment: 33 pages, 9 figures. Additional references, conclusions expanded, and typos corrected. To appear in Classical and Quantum Gravity

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., Comment: 30 pages, 1 figure, new section discussing a wormhole with de Sitter-type ends. To appear in Physical Review D

We analyze the stability of the Cauchy horizon associated with a globally naked, shell-focussing singularity arising from the complete gravitational collapse of a spherical dust cloud. In previous work, we have studied the dynamics of spherical test scalar fields on such a background. In particular, we proved that such fields cannot develop any divergences which propagate along the Cauchy horizon. In the present work, we extend our analysis to the more general case of test fields without symmetries and to linearized gravitational perturbations with odd parity. To this purpose, we first consider test fields possessing a divergence-free stress-energy tensor satisfying the dominant energy condition, and we prove that a suitable energy norm is uniformly bounded in the domain of dependence of the initial slice. In particular, this result implies that free-falling observers co-moving with the dust particles measure a finite energy of the field, even as they cross the Cauchy horizon at points lying arbitrarily close to the central singularity. Next, for the case of Klein-Gordon fields, we derive point-wise bounds from our energy estimates which imply that the scalar field cannot diverge at the Cauchy horizon, except possibly at the central singular point. Finally, we analyze the behaviour of odd-parity, linear gravitational and dust perturbations of the collapsing spacetime. Similarly to the scalar field case, we prove that the relevant gauge-invariant combinations of the metric perturbations stay bounded away from the central singularity, implying that no divergences can propagate in the vacuum region. Our results are in accordance with previous numerical studies and analytic work in the self-similar case., Comment: 27 pages, 3 figures. Conclusion improved. References added. Matches published version

We provide a systematic study for the accretion of a collisionless, relativistic kinetic gas into a nonrotating black hole. To this end, we first solve the relativistic Liouville equation on a Schwarzschild background spacetime. The most general solution for the distribution function is given in terms of appropriate symplectic coordinates on the cotangent bundle, and the associated observables, including the particle current density and stress energy-momentum tensor, are determined. Next, we explore the case where the flow is steady-state and spherically symmetric. Assuming that in the asymptotic region the gas is described by an equilibrium distribution function, we determine the relevant parameters of the accretion flow as a function of the particle density and the temperature of the gas at infinity. In particular, we find that in the low temperature limit the tangential pressure at the horizon is about an order of magnitude larger than the radial one, showing explicitly that a collisionless gas, despite exerting kinetic pressure, behaves very differently than an isotropic perfect fluid, and providing a partial explanation for the known fact that the accretion rate is much lower than in the hydrodynamic case of Bondi-Michel accretion. Finally, we establish the asymptotic stability of the steady-state spherical flows by proving pointwise convergence results which show that a large class of (possibly nonstationary and nonspherical) initial conditions for the distribution function lead to solutions of the Liouville equation which relax in time to a steady-state, spherically symmetric configuration., Comment: 38 pages, 6 figures. Updated version with typos corrected. To appear in Classical and Quantum Gravity

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., Comment: To appear in Living Reviews in Relativity, 191 pages

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., Comment: Corrected factor of 2 errors in Eqs. (A8) and (A9) and a few typos; final version