-
Equilibriums of extremely magnetized compact stars with force-free magnetotunnels
Authors:
Koji Uryu,
Shijun Yoshida,
Eric Gourgoulhon,
Charalampos Markakis,
Kotaro Fujisawa,
Antonios Tsokaros,
Keisuke Taniguchi,
Mina Zamani
Abstract:
We present numerical solutions for stationary and axisymmetric equilibriums of compact stars associated with extremely strong magnetic fields. The interior of the compact stars is assumed to satisfy ideal magnetohydrodynamic (MHD) conditions, while in the region of negligible mass density the force-free conditions or electromagnetic vacuum are assumed. Solving all components of Einstein's equation…
▽ More
We present numerical solutions for stationary and axisymmetric equilibriums of compact stars associated with extremely strong magnetic fields. The interior of the compact stars is assumed to satisfy ideal magnetohydrodynamic (MHD) conditions, while in the region of negligible mass density the force-free conditions or electromagnetic vacuum are assumed. Solving all components of Einstein's equations, Maxwell's equations, ideal MHD equations, and force-free conditions, equilibriums of rotating compact stars associated with mixed poloidal and toroidal magnetic fields are obtained. It is found that in the extreme cases the strong mixed magnetic fields concentrating in a toroidal region near the equatorial surface expel the matter and form a force-free toroidal magnetotunnel. We also introduce a new differential rotation law for computing solutions associated with force-free magnetosphere, and present other extreme models without the magnetotunnel.
△ Less
Submitted 31 March, 2023;
originally announced March 2023.
-
Peeling at extreme black hole horizons
Authors:
Jack Borthwick,
Eric Gourgoulhon,
Jean-Philippe Nicolas
Abstract:
The starting point of this work was an intriguing similarity between the behaviour of fields near a degenerate horizon and near the infinity of an asymptotically flat spacetime, as revealed by the scattering theory for Dirac fields in the ``exterior'' region of the extreme Kerr - de Sitter black hole, developed by one of the authors (JB). However, in that situation, the comparison was somewhat clo…
▽ More
The starting point of this work was an intriguing similarity between the behaviour of fields near a degenerate horizon and near the infinity of an asymptotically flat spacetime, as revealed by the scattering theory for Dirac fields in the ``exterior'' region of the extreme Kerr - de Sitter black hole, developed by one of the authors (JB). However, in that situation, the comparison was somewhat clouded by some of the analytical techniques used in intermediate steps of the proof. The aim of the present work is to clarify the comparison further by studying instead the peeling behaviour of solutions to the wave equation at an extremal horizon. We focus first on the extreme Reissner-Nordström black hole, for which the Couch-Torrence inversion (a global conformal isometry that exchanges the horizon and infinity) makes the analogy explicit. Then, we explore more general spherically symmetric situations using the Couch-Torrence inversion outside of its natural context.
△ Less
Submitted 25 March, 2023;
originally announced March 2023.
-
High precision numerical sequences of rotating hairy black holes
Authors:
Gustavo García,
Eric Gourgoulhon,
Philippe Grandclément,
Marcelo Salgado
Abstract:
We analyze numerically the existence of regular stationary rotating hairy black holes within the framework of general relativity, which are the result of solving the Einstein-Klein-Gordon system for a complex-valued scalar field under suitable boundary (regularity and asymptotically flat) conditions. To that aim we solve the corresponding system of elliptic partial differential equations using spe…
▽ More
We analyze numerically the existence of regular stationary rotating hairy black holes within the framework of general relativity, which are the result of solving the Einstein-Klein-Gordon system for a complex-valued scalar field under suitable boundary (regularity and asymptotically flat) conditions. To that aim we solve the corresponding system of elliptic partial differential equations using spectral methods which are specially suited for such a numerical task. In order to obtain such system of equations we employ a parametrization for the metric that corresponds to quasi-isotropic coordinates (QIC) that have been used in the past for analyzing different kinds of stationary rotating relativistic systems. Our findings are in agreement with those reported originally by Herdeiro \& Radu ite{Herdeiro2014,Herdeiro2015}. The method is submitted to several analytic and numerical tests, which include the recovery of the Kerr solution in QIC and the cloud solutions in the Kerr background. We report different global quantities that allow us to determine the contribution of the boson hair to the spacetime, as well as relevant quantities at the horizon, like the surface gravity. The latter indicates to what extent the hairy solutions approach the extremal limit, noting that for this kind of solutions the ratio of the angular momentum per squared mass $J_\infty/M^2_{\rm ADM}$ can be larger than unity due to the contribution of the scalar hair, a situation which differs from the Kerr metric where this parameter is bounded according to $0\leq |J/M^2| \leq 1$ with the upper bound corresponding to the extremal case.
△ Less
Submitted 13 February, 2023;
originally announced February 2023.
-
Equatorial orbits and imaging of hairy cubic Galileon black holes
Authors:
Karim Van Aelst,
Eric Gourgoulhon,
Frederic H. Vincent
Abstract:
Null and timelike equatorial orbits are investigated in a family of hairy black holes in the cubic Galileon theory. These include rotating generalizations of static black hole metrics supporting a time-dependent scalar field. Depending on the coupling and rotation, the properties of the geodesics expectedly deviate from general relativity. In particular, it is found that stable circular geodesics…
▽ More
Null and timelike equatorial orbits are investigated in a family of hairy black holes in the cubic Galileon theory. These include rotating generalizations of static black hole metrics supporting a time-dependent scalar field. Depending on the coupling and rotation, the properties of the geodesics expectedly deviate from general relativity. In particular, it is found that stable circular geodesics only exist below a critical coupling, which is related to the existence of an outermost stable circular orbit. Focusing on the strong-field region, images of an accretion disk are also produced to highlight tendencies that would constrain the model given further accurate observations of supermassive black holes.
△ Less
Submitted 26 January, 2022; v1 submitted 2 March, 2021;
originally announced March 2021.
-
Geometric modeling of M87* as a Kerr black hole or a non-Kerr compact object
Authors:
F. H. Vincent,
M. Wielgus,
M. A. Abramowicz,
E. Gourgoulhon,
J. -P. Lasota,
T. Paumard,
G. Perrin
Abstract:
Context. The Event Horizon Telescope (EHT) collaboration recently obtained first images of the surroundings of the supermassive compact object M87* at the center of the galaxy M87. Aims. We want to develop a simple analytic disk model for the accretion flow of M87*. Compared to general-relativistic magnetohydrodynamic (GRMHD) models, it has the advantage of being independent of the turbulent chara…
▽ More
Context. The Event Horizon Telescope (EHT) collaboration recently obtained first images of the surroundings of the supermassive compact object M87* at the center of the galaxy M87. Aims. We want to develop a simple analytic disk model for the accretion flow of M87*. Compared to general-relativistic magnetohydrodynamic (GRMHD) models, it has the advantage of being independent of the turbulent character of the flow, and controlled by only few easy-to-interpret, physically meaningful parameters. We want to use this model to predict the image of M87* assuming that it is either a Kerr black hole, or an alternative compact object. Methods. We compute the synchrotron emission from the disk model and propagate the resulting light rays to the far-away observer by means of relativistic ray tracing. Such computations are performed assuming different spacetimes (Kerr, Minkowski, non-rotating ultracompact star, rotating boson star or Lamy spinning wormhole). We perform numerical fits of these models to the EHT data. Results. We discuss the highly-lensed features of Kerr images and show that they are intrinsically linked to the accretion-flow properties, and not only to gravitation. This fact is illustrated by the notion of secondary ring that we introduce. Our model of spinning Kerr black hole predicts mass and orientation consistent with the EHT interpretation. The non-Kerr images result in similar quality of the numerical fits and may appear very similar to Kerr images, once blurred to the EHT resolution. This implies that a strong test of the Kerr spacetime may be out of reach with the current data. We notice that future developments of the EHT could alter this situation. Conclusions. Our results show the importance of studying alternatives to the Kerr spacetime in order to be able to test the Kerr paradigm unambiguously.
△ Less
Submitted 27 October, 2020; v1 submitted 21 February, 2020;
originally announced February 2020.
-
Prospects for Fundamental Physics with LISA
Authors:
Enrico Barausse,
Emanuele Berti,
Thomas Hertog,
Scott A. Hughes,
Philippe Jetzer,
Paolo Pani,
Thomas P. Sotiriou,
Nicola Tamanini,
Helvi Witek,
Kent Yagi,
Nicolas Yunes,
T. Abdelsalhin,
A. Achucarro,
K. V. Aelst,
N. Afshordi,
S. Akcay,
L. Annulli,
K. G. Arun,
I. Ayuso,
V. Baibhav,
T. Baker,
H. Bantilan,
T. Barreiro,
C. Barrera-Hinojosa,
N. Bartolo
, et al. (296 additional authors not shown)
Abstract:
In this paper, which is of programmatic rather than quantitative nature, we aim to further delineate and sharpen the future potential of the LISA mission in the area of fundamental physics. Given the very broad range of topics that might be relevant to LISA, we present here a sample of what we view as particularly promising directions, based in part on the current research interests of the LISA sc…
▽ More
In this paper, which is of programmatic rather than quantitative nature, we aim to further delineate and sharpen the future potential of the LISA mission in the area of fundamental physics. Given the very broad range of topics that might be relevant to LISA, we present here a sample of what we view as particularly promising directions, based in part on the current research interests of the LISA scientific community in the area of fundamental physics. We organize these directions through a "science-first" approach that allows us to classify how LISA data can inform theoretical physics in a variety of areas. For each of these theoretical physics classes, we identify the sources that are currently expected to provide the principal contribution to our knowledge, and the areas that need further development. The classification presented here should not be thought of as cast in stone, but rather as a fluid framework that is amenable to change with the flow of new insights in theoretical physics.
△ Less
Submitted 27 April, 2020; v1 submitted 27 January, 2020;
originally announced January 2020.
-
Hairy rotating black holes in cubic Galileon theory
Authors:
Karim Van Aelst,
Eric Gourgoulhon,
Philippe Grandclément,
Christos Charmousis
Abstract:
Numerical solutions for asymptotically flat rotating black holes in the cubic Galileon theory are presented. These black holes are endowed with a nontrivial scalar field and exhibit a non-Schwarzschild behaviour: faster than $1/r$ convergence to Minkowski spacetime at spatial infinity and hence vanishing of the Komar mass. The metrics are compared with the Kerr metric for various couplings and ang…
▽ More
Numerical solutions for asymptotically flat rotating black holes in the cubic Galileon theory are presented. These black holes are endowed with a nontrivial scalar field and exhibit a non-Schwarzschild behaviour: faster than $1/r$ convergence to Minkowski spacetime at spatial infinity and hence vanishing of the Komar mass. The metrics are compared with the Kerr metric for various couplings and angular velocities. Their physical properties are extracted and show significant deviations from the Kerr case.
△ Less
Submitted 23 January, 2020; v1 submitted 18 October, 2019;
originally announced October 2019.
-
New code for equilibriums and quasiequilibrium initial data of compact objects. IV. Rotating relativistic stars with mixed poloidal and toroidal magnetic fields
Authors:
Koji Uryu,
Shijun Yoshida,
Eric Gourgoulhon,
Charalampos Markakis,
Kotaro Fujisawa,
Antonios Tsokaros,
Keisuke Taniguchi,
Yoshiharu Eriguchi
Abstract:
A new code for computing fully general relativistic solutions of strongly magnetized rapidly rotating compact stars is developed as a part of the COCAL (Compact Object CALculator) code. The full set of Einstein's equations, Maxwell's equations and magnetohydrodynamic equations are consistently solved assuming perfect conductivity, stationarity, and axisymmetry, and strongly magnetized solutions as…
▽ More
A new code for computing fully general relativistic solutions of strongly magnetized rapidly rotating compact stars is developed as a part of the COCAL (Compact Object CALculator) code. The full set of Einstein's equations, Maxwell's equations and magnetohydrodynamic equations are consistently solved assuming perfect conductivity, stationarity, and axisymmetry, and strongly magnetized solutions associated with mixed poloidal and toroidal components of magnetic fields are successfully obtained in generic (non-circular) spacetimes. We introduce the formulation of the problem and the numerical method in detail, then present examples of extremely magnetized compact star solutions and their convergence tests. It is found that, in extremely magnetized stars, the stellar matter can be expelled from the region of strongest toroidal fields. Hence we conjecture that a toroidal electro-vacuum region may appear inside of the extremely magnetized compact stars, which may seem like the neutron star becoming the strongest toroidal solenoid coil in the universe.
△ Less
Submitted 25 June, 2019;
originally announced June 2019.
-
Gravitational waves from bodies orbiting the Galactic Center black hole and their detectability by LISA
Authors:
Eric Gourgoulhon,
Alexandre Le Tiec,
Frederic H. Vincent,
Niels Warburton
Abstract:
We present the first fully relativistic study of gravitational radiation from bodies in circular equatorial orbits around the massive black hole at the Galactic Center, Sgr A* and we assess the detectability of various kinds of objects by the gravitational wave detector LISA. Our computations are based on the theory of perturbations of the Kerr spacetime and take into account the Roche limit induc…
▽ More
We present the first fully relativistic study of gravitational radiation from bodies in circular equatorial orbits around the massive black hole at the Galactic Center, Sgr A* and we assess the detectability of various kinds of objects by the gravitational wave detector LISA. Our computations are based on the theory of perturbations of the Kerr spacetime and take into account the Roche limit induced by tidal forces in the Kerr metric. The signal-to-noise ratio in the LISA detector, as well as the time spent in LISA band, are evaluated. We have implemented all the computational tools in an open-source SageMath package, within the Black Hole Perturbation Toolkit framework. We find that white dwarfs, neutrons stars, stellar black holes, primordial black holes of mass larger than $10^{-4} M_\odot$, main-sequence stars of mass lower than $\sim 2.5\, M_\odot$ and brown dwarfs orbiting Sgr A* are all detectable in one year of LISA data with a signal-to-noise ratio above 10 during at least $10^5$ years in the slow inspiral towards either the innermost stable circular orbit (compact objects) or the Roche limit (main-sequence stars and brown dwarfs). The longest times in-band, of the order of $10^6$ years, are achieved for primordial black holes of mass $\sim 10^{-3} M_\odot$ down to $10^{-5} M_\odot$, depending on the spin of Sgr A*, as well as for brown dwarfs, just followed by white dwarfs and low mass main-sequence stars. The long time in-band of these objects makes Sgr A* a valuable target for LISA. We also consider bodies on close circular orbits around the massive black hole in the nucleus of the nearby galaxy M32 and find that, among them, compact objects and brown dwarfs stay for $10^3$ to $10^4$ years in LISA band with a 1-year signal-to-noise ratio above 10.
△ Less
Submitted 9 June, 2019; v1 submitted 5 March, 2019;
originally announced March 2019.
-
Symbolic tensor calculus on manifolds: a SageMath implementation
Authors:
Eric Gourgoulhon,
Marco Mancini
Abstract:
These lecture notes present a method for symbolic tensor calculus that (i) runs on fully specified smooth manifolds (described by an atlas), (ii) is not limited to a single coordinate chart or vector frame, (iii) runs even on non-parallelizable manifolds and (iv) is independent of the symbolic backend used to perform calculus at the level of coordinate expressions. In addition to the main ideas, w…
▽ More
These lecture notes present a method for symbolic tensor calculus that (i) runs on fully specified smooth manifolds (described by an atlas), (ii) is not limited to a single coordinate chart or vector frame, (iii) runs even on non-parallelizable manifolds and (iv) is independent of the symbolic backend used to perform calculus at the level of coordinate expressions. In addition to the main ideas, we discuss some details of the implementation in the open-source mathematics software system SageMath, which has been performed via the SageManifolds project.
△ Less
Submitted 12 November, 2018; v1 submitted 19 April, 2018;
originally announced April 2018.
-
Imaging a non-singular rotating black hole at the center of the Galaxy
Authors:
Frédéric Lamy,
Eric Gourgoulhon,
Thibaut Paumard,
Frédéric H. Vincent
Abstract:
We show that the rotating generalization of Hayward's non-singular black hole previously studied in the literature is geodesically incomplete, and that its straightforward extension leads to a singular spacetime. We present another extension, which is devoid of any curvature singularity. The obtained metric depends on three parameters and, depending on their values, yields an event horizon or not.…
▽ More
We show that the rotating generalization of Hayward's non-singular black hole previously studied in the literature is geodesically incomplete, and that its straightforward extension leads to a singular spacetime. We present another extension, which is devoid of any curvature singularity. The obtained metric depends on three parameters and, depending on their values, yields an event horizon or not. These two regimes, named respectively regular rotating Hayward black hole and naked rotating wormhole, are studied both numerically and analytically. In preparation for the upcoming results of the Event Horizon Telescope, the images of an accretion torus around Sgr A*, the supermassive object at the center of the Galaxy, are computed. These images contain, even in the absence of a horizon, a central faint region which bears a resemblance to the shadow of Kerr black holes and emphasizes the difficulty of claiming the existence of an event horizon from the analysis of strong-field images. The frequencies of the co- and contra-rotating orbits at the innermost stable circular orbit (ISCO) in this geometry are also computed, in the hope that quasi-periodic oscillations may permit to compare this model with Kerr's black hole on observational grounds.
△ Less
Submitted 12 May, 2018; v1 submitted 5 February, 2018;
originally announced February 2018.
-
Accurate ray tracing of realistic neutron star atmospheres for constraining their parameters
Authors:
F. H. Vincent,
M. Bejger,
A. Rozanska,
O. Straub,
T. Paumard,
M. Fortin,
J. Madej,
A. Majczyna,
E. Gourgoulhon,
P. Haensel,
L. Zdunik,
B. Beldycki
Abstract:
Thermal dominated X-ray spectra of neutron stars in quiescent transient X-ray binaries and neutron stars that undergo thermonuclear bursts are sensitive to mass and radius. The mass-radius relation of neutron stars depends on the equation of state that governs their interior. Constraining this relation accurately is thus of fundamental importance to understand the nature of dense matter. In this c…
▽ More
Thermal dominated X-ray spectra of neutron stars in quiescent transient X-ray binaries and neutron stars that undergo thermonuclear bursts are sensitive to mass and radius. The mass-radius relation of neutron stars depends on the equation of state that governs their interior. Constraining this relation accurately is thus of fundamental importance to understand the nature of dense matter. In this context we introduce a pipeline to calculate realistic model spectra of rotating neutron stars with hydrogen and helium atmospheres. An arbitrarily fast rotating neutron star with a given equation of state generates the spacetime in which the atmosphere emits radiation. We use the Lorene/nrotstar code to compute the spacetime numerically and the ATM24 code to solve the radiative transfer equations self-consistently. Emerging specific intensity spectra are then ray-traced through the neutron star's spacetime from the atmosphere to a distant observer with the Gyoto code. Here, we present and test our fully relativistic numerical pipeline. To discuss and illustrate the importance of realistic atmosphere models we compare our model spectra to simpler models like the commonly used isotropic color-corrected blackbody emission. We highlight the importance of considering realistic model-atmosphere spectra together with relativistic ray tracing to obtain accurate predictions. We also insist on the crucial impact of the star's rotation on the observables. Finally, we close a controversy that has been appearing in the literature in the recent years regarding the validity of the ATM24 code.
△ Less
Submitted 21 February, 2018; v1 submitted 7 November, 2017;
originally announced November 2017.
-
Comparing timelike geodesics around a Kerr black hole and a boson star
Authors:
M. Grould,
Z. Meliani,
F. H. Vincent,
P. Grandclément,
E. Gourgoulhon
Abstract:
The second-generation beam combiner at the Very Large Telescope (VLT), GRAVITY, observes the stars orbiting the compact object located at the center of our galaxy, with an unprecedented astrometric accuracy of 10 $μ$as. The nature of this compact source is still unknown since black holes are not the only candidates explaining the four million solar masses at the Galactic center. Boson stars are su…
▽ More
The second-generation beam combiner at the Very Large Telescope (VLT), GRAVITY, observes the stars orbiting the compact object located at the center of our galaxy, with an unprecedented astrometric accuracy of 10 $μ$as. The nature of this compact source is still unknown since black holes are not the only candidates explaining the four million solar masses at the Galactic center. Boson stars are such an alternative model to black holes. This paper focuses on the study of trajectories of stars orbiting a boson star and a Kerr black hole. We put in light strong differences between orbits obtained in both metrics when considering stars with sufficiently close pericenters to the compact object, typically $\lesssim 30~M$. Discovery of closer stars to the Galactic center than the S2 star by the GRAVITY instrument would thus be a powerful tool to possibly constrain the nature of the central source.
△ Less
Submitted 18 September, 2017;
originally announced September 2017.
-
Conservation laws and evolution schemes in geodesic, hydrodynamic and magnetohydrodynamic flows
Authors:
Charalampos Markakis,
Kōji Uryū,
Eric Gourgoulhon,
Jean-Philippe Nicolas,
Nils Andersson,
Athina Pouri,
Vojtech Witzany
Abstract:
Carter and Lichnerowicz have established that barotropic fluid flows are conformally geodesic and obey Hamilton's principle. This variational approach can accommodate neutral, or charged and poorly conducting, fluids. We show that, unlike what has been previously thought, this approach can also accommodate perfectly conducting magnetofluids, via the Bekenstein-Oron description of ideal magnetohydr…
▽ More
Carter and Lichnerowicz have established that barotropic fluid flows are conformally geodesic and obey Hamilton's principle. This variational approach can accommodate neutral, or charged and poorly conducting, fluids. We show that, unlike what has been previously thought, this approach can also accommodate perfectly conducting magnetofluids, via the Bekenstein-Oron description of ideal magnetohydrodynamics. When Noether symmetries associated with Killing vectors or tensors are present in geodesic flows, they lead to constants of motion polynomial in the momenta. We generalize these concepts to hydrodynamic flows. Moreover, the Hamiltonian descriptions of ideal magnetohydrodynamics allow one to cast the evolution equations into a hyperbolic form useful for evolving rotating or binary compact objects with magnetic fields in numerical general relativity. Conserved circulation laws, such as those of Kelvin, Alfvén and Bekenstein-Oron, emerge simply as special cases of the Poincaré-Cartan integral invariant of Hamiltonian systems. We use this approach to obtain an extension of Kelvin's theorem to baroclinic (non-isentropic) fluids, based on a temperature-dependent time parameter. We further extend this result to perfectly or poorly conducting baroclinic magnetoflows. Finally, in the barotropic case, such magnetoflows are shown to also be geodesic, albeit in a Finsler (rather than Riemann) space.
△ Less
Submitted 29 December, 2016;
originally announced December 2016.
-
Astrophysical imaging of Kerr black holes with scalar hair
Authors:
F. H. Vincent,
E. Gourgoulhon,
C. Herdeiro,
E. Radu
Abstract:
We address the astrophysical imaging of a family of deformed Kerr black holes (BHs). These are stationary, asymptotically flat black hole (BH) spacetimes, that are solutions of General Relativity minimally coupled to a massive, complex scalar field: Kerr BHs with scalar hair (KBHsSH). Such BHs bifurcate from the vacuum Kerr solution and can be regarded as a horizon within a rotating boson star. In…
▽ More
We address the astrophysical imaging of a family of deformed Kerr black holes (BHs). These are stationary, asymptotically flat black hole (BH) spacetimes, that are solutions of General Relativity minimally coupled to a massive, complex scalar field: Kerr BHs with scalar hair (KBHsSH). Such BHs bifurcate from the vacuum Kerr solution and can be regarded as a horizon within a rotating boson star. In a recent letter, it was shown that KBHsSH can exhibit very distinct shadows from the ones of their vacuum counterparts. The setup therein, however, considered the light source to be a celestial sphere sufficiently far away from the BH. Here, we analyse KBHsSH surrounded by an emitting torus of matter, simulating a more realistic astrophysical environment, and study the corresponding lensing of light as seen by a very far away observer, to appropriately model ground-based observations of Sgr A*. We find that the differences in imaging between KBHsSH and comparable vacuum Kerr BHs remain, albeit less dramatic than those observed for the corresponding shadows in the previous setup. In particular, we highlight two observables that might allow differentiating KBHsSH and Kerr BHs. The first is the angular size of the photon ring (in a Kerr spacetime) or lensing ring (in a KBHSH spacetime), the latter being significantly smaller for sufficiently non-Kerr-like spacetimes. The second is the existence of an edge in the intensity distribution (the photon ring in Kerr spacetime). This edge can disappear for very non-Kerr-like KBHsSH. It is plausible, therefore, that sufficiently precise Very Long Baseline Interferometric observations of BH candidates can constrain this model.
△ Less
Submitted 19 October, 2016; v1 submitted 14 June, 2016;
originally announced June 2016.
-
Thermalization of holographic Wilson loops in spacetimes with spatial anisotropy
Authors:
Dmitry S. Ageev,
Irina Ya. Aref'eva,
Anastasia A. Golubtsova,
Eric Gourgoulhon
Abstract:
In this paper, we study behaviour of Wilson loops in the boost-invariant nonequilibrium anisotropic quark-gluon plasma produced in heavy-ion collisions within the holographic approach. We describe the thermalization studying the evolution of the Vaidya metric in the boost-invariant and spatially anisotropic background. To probe the system during this process we calculate rectangular Wilson loops o…
▽ More
In this paper, we study behaviour of Wilson loops in the boost-invariant nonequilibrium anisotropic quark-gluon plasma produced in heavy-ion collisions within the holographic approach. We describe the thermalization studying the evolution of the Vaidya metric in the boost-invariant and spatially anisotropic background. To probe the system during this process we calculate rectangular Wilson loops oriented in different spatial directions. We find that anisotropic effects are more visible for the Wilson loops lying in the transversal plane unlike the Wilson loops with partially longitudinal orientation. In particular, we observe that the Wilson loops can thermalizes first unlike to the order of the isotropic model. We see that Wilson loops on transversal contours have the shortest thermalization time. We also calculate the string tension and the pseudopotential at different temperatures for the static quark-gluon plasma. We show that the pseudopotential related to the configuration on the transversal plane has the screened Cornell form. We also show that the jet-quenching parameter related with the average of the light-like Wilson loop exhibits the dependence on orientations.
△ Less
Submitted 22 May, 2018; v1 submitted 13 June, 2016;
originally announced June 2016.
-
Analytic black branes in Lifshitz-like backgrounds and thermalization
Authors:
Irina Ya. Aref'eva,
Anastasia A. Golubtsova,
Eric Gourgoulhon
Abstract:
Using black brane solutions in 5d Lifshitz-like backgrounds with arbitrary dynamical exponent $ν$, we construct the Vaidya geometry, asymptoting to the Lifshitz-like spacetime, which represents a thin shell infalling at the speed of light. We apply the new Lifshitz-Vaidya background to study the thermalization process of the quark-gluon plasma via the thin shell approach previously successfully us…
▽ More
Using black brane solutions in 5d Lifshitz-like backgrounds with arbitrary dynamical exponent $ν$, we construct the Vaidya geometry, asymptoting to the Lifshitz-like spacetime, which represents a thin shell infalling at the speed of light. We apply the new Lifshitz-Vaidya background to study the thermalization process of the quark-gluon plasma via the thin shell approach previously successfully used in several backgrounds. We find that the thermalization depends on the chosen direction because of the spatial anisotropy. The plasma thermalizes thus faster in the transversal direction than in the longitudinal one. To probe the system described by the Lifshitz-like backgrounds, we also calculate the holographic entanglement entropy for the subsystems delineated along both transversal and longitudinal directions. We show that the entropy has some universality in the behavior for both subsystems. At the same time, we find that certain characteristics strongly depend on the critical exponent $ν$.
△ Less
Submitted 22 October, 2016; v1 submitted 22 January, 2016;
originally announced January 2016.
-
Circular geodesics and thick tori around rotating boson stars
Authors:
Z. Meliani,
F. H. Vincent,
P. Grandclément,
E. Gourgoulhon,
R. Monceau-Baroux,
O. Straub
Abstract:
Accretion disks play an important role in the evolution of their relativistic inner compact objects. The emergence of a new generation of interferometers will allow to resolve these accretion disks and provide more information about the properties of the central gravitating object. Due to this instrumental leap forward it is crucial to investigate the accretion disk physics near various types of i…
▽ More
Accretion disks play an important role in the evolution of their relativistic inner compact objects. The emergence of a new generation of interferometers will allow to resolve these accretion disks and provide more information about the properties of the central gravitating object. Due to this instrumental leap forward it is crucial to investigate the accretion disk physics near various types of inner compact objects now to deduce later constraints on the central objects from observations. A possible candidate for the inner object is the boson star. Here, we will try to analyze the differences between accretion structures surrounding boson stars and black holes. We aim at analysing the physics of circular geodesics around boson stars and study simple thick accretion tori (so-called Polish doughnuts) in the vicinity of these stars. We realize a detailed study of the properties of circular geodesics around boson stars. We then perform a parameter study of thick tori with constant angular momentum surrounding boson stars. This is done using the boson star models computed by a code constructed with the spectral solver library KADATH. We demonstrate that all the circular stable orbits are bound. In the case of a constant angular momentum torus, a cusp in the torus surface exists only for boson stars with a strong gravitational scalar field. Moreover, for each inner radius of the disk, the allowed specific angular momentum values lie within a constrained range which depends on the boson star considered. We show that the accretion tori around boson stars have different characteristics than in the vicinity of a black hole. With future instruments it could be possible to use these differences to constrain the nature of compact objects.
△ Less
Submitted 14 October, 2015;
originally announced October 2015.
-
Imaging a boson star at the Galactic center
Authors:
F. H. Vincent,
Z. Meliani,
P. Grandclement,
E. Gourgoulhon,
O. Straub
Abstract:
Millimeter very long baseline interferometry will soon produce accurate images of the closest surroundings of the supermassive compact object at the center of the Galaxy, Sgr A*. These images may reveal the existence of a central faint region, the so-called shadow, which is often interpreted as the observable consequence of the event horizon of a black hole. In this paper, we compute images of an…
▽ More
Millimeter very long baseline interferometry will soon produce accurate images of the closest surroundings of the supermassive compact object at the center of the Galaxy, Sgr A*. These images may reveal the existence of a central faint region, the so-called shadow, which is often interpreted as the observable consequence of the event horizon of a black hole. In this paper, we compute images of an accretion torus around Sgr A* assuming this compact object is a boson star, i.e. an alternative to black holes within general relativity, with no event horizon and no hard surface. We show that very relativistic rotating boson stars produce images extremely similar to Kerr black holes, showing in particular shadow-like and photon-ring-like structures. This result highlights the extreme difficulty of unambiguously telling the existence of an event horizon from strong-field images.
△ Less
Submitted 16 March, 2016; v1 submitted 14 October, 2015;
originally announced October 2015.
-
A characterization of 3+1 spacetimes via the Simon-Mars tensor
Authors:
C. Somé,
P. Grandclément,
E. Gourgoulhon
Abstract:
We present the 3+1 decomposition of the Simon-Mars tensor, which has the property of being identically zero for a vacuum and asymptotically flat spacetime if and only if the latter is locally isometric to the Kerr spacetime. Using this decomposition we form two dimensionless scalar fields. Computing these scalars provides a simple way of comparing locally a generic (even non vacuum and non analyti…
▽ More
We present the 3+1 decomposition of the Simon-Mars tensor, which has the property of being identically zero for a vacuum and asymptotically flat spacetime if and only if the latter is locally isometric to the Kerr spacetime. Using this decomposition we form two dimensionless scalar fields. Computing these scalars provides a simple way of comparing locally a generic (even non vacuum and non analytic) stationary spacetime to Kerr. As an illustration, we evaluate the Simon-Mars scalars for numerical solutions of the Einstein equations generated by boson stars and neutron stars, for analytic solutions of the Einstein equations such as Curzon-Chazy spacetime and $δ=2$ Tomimatsu-Sato spacetime, and for an approximate solution of the Einstein equations : the modified Kerr metric, which is an example of a parametric deviation from Kerr spacetime.
△ Less
Submitted 19 December, 2014;
originally announced December 2014.
-
Tensor calculus with open-source software: the SageManifolds project
Authors:
Eric Gourgoulhon,
Michal Bejger,
Marco Mancini
Abstract:
The SageManifolds project aims at extending the mathematics software system Sage towards differential geometry and tensor calculus. Like Sage, SageManifolds is free, open-source and is based on the Python programming language. We discuss here some details of the implementation, which relies on Sage's parent/element framework, and present a concrete example of use.
The SageManifolds project aims at extending the mathematics software system Sage towards differential geometry and tensor calculus. Like Sage, SageManifolds is free, open-source and is based on the Python programming language. We discuss here some details of the implementation, which relies on Sage's parent/element framework, and present a concrete example of use.
△ Less
Submitted 21 December, 2014; v1 submitted 15 December, 2014;
originally announced December 2014.
-
Equilibrium solutions of relativistic rotating stars with mixed poloidal and toroidal magnetic fields
Authors:
Koji Uryu,
Eric Gourgoulhon,
Charalampos Markakis,
Kotaro Fujisawa,
Antonios Tsokaros,
Yoshiharu Eriguchi
Abstract:
Stationary and axisymmetric solutions of relativistic rotating stars with strong mixed poloidal and toroidal magnetic fields are obtained numerically. Because of the mixed components of the magnetic field, the underlying stationary and axisymmetric spacetimes are no longer circular. These configurations are computed from the full set of the Einstein-Maxwell equations, Maxwell's equations and from…
▽ More
Stationary and axisymmetric solutions of relativistic rotating stars with strong mixed poloidal and toroidal magnetic fields are obtained numerically. Because of the mixed components of the magnetic field, the underlying stationary and axisymmetric spacetimes are no longer circular. These configurations are computed from the full set of the Einstein-Maxwell equations, Maxwell's equations and from first integrals and integrability conditions of the magnetohydrodynamic equilibrium equations. After a brief introduction of the formulation of the problem, we present the first results for highly deformed magnetized rotating compact stars.
△ Less
Submitted 14 October, 2014;
originally announced October 2014.
-
Models of rotating boson stars and geodesics around them: new type of orbits
Authors:
Philippe Grandclement,
Claire Somé,
Eric Gourgoulhon
Abstract:
We have developed a highly accurate numerical code capable of solving the coupled Einstein-Klein-Gordon system, in order to construct rotating boson stars in general relativity. Free fields and self-interacting fields, with quartic and sextic potentials, are considered. In particular, we present the first numerical solutions of rotating boson stars with rotational quantum number $k=3$ and $k=4$, a…
▽ More
We have developed a highly accurate numerical code capable of solving the coupled Einstein-Klein-Gordon system, in order to construct rotating boson stars in general relativity. Free fields and self-interacting fields, with quartic and sextic potentials, are considered. In particular, we present the first numerical solutions of rotating boson stars with rotational quantum number $k=3$ and $k=4$, as well as the first determination of the maximum mass of free-field boson stars with $k=2$. We have also investigated timelike geodesics in the spacetime generated by a rotating boson star for $k=1$, $2$ and $3$. A numerical integration of the geodesic equation has enabled us to identify a peculiar type of orbits: the zero-angular-momentum ones. These orbits pass very close to the center and are qualitatively different from orbits around a Kerr black hole. Should such orbits be observed, they would put stringent constraints on astrophysical compact objects like the Galactic center.
△ Less
Submitted 19 July, 2014; v1 submitted 19 May, 2014;
originally announced May 2014.
-
Extracting black-hole rotational energy: The generalized Penrose process
Authors:
J. -P. Lasota,
E. Gourgoulhon,
M. Abramowicz,
A. Tchekhovskoy,
R. Narayan
Abstract:
In the case involving particles the necessary and sufficient condition for the Penrose process to extract energy from a rotating black hole is absorption of particles with negative energies and angular momenta. No torque at the black-hole horizon occurs. In this article we consider the case of arbitrary fields or matter described by an unspecified, general energy-momentum tensor $T_{μν}$ and show…
▽ More
In the case involving particles the necessary and sufficient condition for the Penrose process to extract energy from a rotating black hole is absorption of particles with negative energies and angular momenta. No torque at the black-hole horizon occurs. In this article we consider the case of arbitrary fields or matter described by an unspecified, general energy-momentum tensor $T_{μν}$ and show that the necessary and sufficient condition for extraction of a black hole's rotational energy is analogous to that in the mechanical Penrose process: absorption of negative energy and negative angular momentum. We also show that a necessary condition for the Penrose process to occur is for the Noether current (the conserved energy-momentum density vector) to be spacelike or past directed (timelike or null) on some part of the horizon. In the particle case, our general criterion for the occurrence of a Penrose process reproduces the standard result. In the case of relativistic jet-producing "magnetically arrested disks" we show that the negative energy and angular-momentum absorption condition is obeyed when the Blandford-Znajek mechanism is at work, and hence the high energy extraction efficiency up to $\sim 300\%$ found in recent numerical simulations of such accretion flows results from tapping the black hole's rotational energy through the Penrose process. We show how black-hole rotational energy extraction works in this case by describing the Penrose process in terms of the Noether current.
△ Less
Submitted 2 February, 2014; v1 submitted 28 October, 2013;
originally announced October 2013.
-
Constraining the nature of the Galactic center black hole Sgr A* with present and future observations
Authors:
F. H. Vincent,
E. Gourgoulhon,
O. Straub,
M. Abramowicz,
J. Novak,
T. Paumard,
G. Perrin
Abstract:
The Galactic center is an ideal laboratory to study strong-field general relativistic phenomena, as the supermassive black hole Sgr A* has the biggest angular Schwarzschild radius among all black holes. This article presents three different ways of using the immediate surroundings of Sgr A* as a probe in order to either constrain its spin, or even test the very nature of this compact object.
The Galactic center is an ideal laboratory to study strong-field general relativistic phenomena, as the supermassive black hole Sgr A* has the biggest angular Schwarzschild radius among all black holes. This article presents three different ways of using the immediate surroundings of Sgr A* as a probe in order to either constrain its spin, or even test the very nature of this compact object.
△ Less
Submitted 3 September, 2012;
originally announced September 2012.
-
3+1 geodesic equation and images in numerical spacetimes
Authors:
Frederic H. Vincent,
Eric Gourgoulhon,
Jérôme Novak
Abstract:
The equations governing null and timelike geodesics are derived within the 3+1 formalism of general relativity. In addition to the particle's position, they encompass an evolution equation for the particle's energy leading to a 3+1 expression of the redshift factor for photons. An important application is the computation of images and spectra in spacetimes arising from numerical relativity, via th…
▽ More
The equations governing null and timelike geodesics are derived within the 3+1 formalism of general relativity. In addition to the particle's position, they encompass an evolution equation for the particle's energy leading to a 3+1 expression of the redshift factor for photons. An important application is the computation of images and spectra in spacetimes arising from numerical relativity, via the ray-tracing technique. This is illustrated here by images of numerically computed stationary neutron stars and dynamical neutron stars collapsing to a black hole.
△ Less
Submitted 24 October, 2012; v1 submitted 20 August, 2012;
originally announced August 2012.
-
Modelling the black hole silhouette in Sgr A* with ion tori
Authors:
Odele Straub,
Frederic H. Vincent,
Marek A. Abramowicz,
Eric Gourgoulhon,
Thibaut Paumard
Abstract:
We calculate the "observed at infinity" image and spectrum of the accretion structure in Sgr A*, by modelling it as an optically thin, constant angular momentum ion torus in hydrodynamic equilibrium. The physics we consider includes a two-temperature plasma, a toroidal magnetic field, as well as radiative cooling by bremsstrahlung, synchrotron and inverse Compton processes. Our relativistic model…
▽ More
We calculate the "observed at infinity" image and spectrum of the accretion structure in Sgr A*, by modelling it as an optically thin, constant angular momentum ion torus in hydrodynamic equilibrium. The physics we consider includes a two-temperature plasma, a toroidal magnetic field, as well as radiative cooling by bremsstrahlung, synchrotron and inverse Compton processes. Our relativistic model has the virtue of being fully analytic and very simple, depending only on eight tunable parameters: the black hole spin and the inclination of the spin axis to our line of sight, the torus angular momentum, the polytropic index, the magnetic to total pressure ratio, the central values of density and electron temperature and the ratio of electron to ion temperatures. The observed image and spectrum are calculated numerically using the ray-tracing code GYOTO. Our results demonstrate that the ion torus model is able to account for the main features of the accretion structure surrounding Sgr A*.
△ Less
Submitted 12 March, 2012;
originally announced March 2012.
-
Quasi-equilibrium models of magnetized compact objects
Authors:
Charalampos Markakis,
Koji Uryu,
Eric Gourgoulhon
Abstract:
We report work towards a relativistic formulation for modeling strongly magnetized neutron stars, rotating or in a close circular orbit around another neutron star or black hole, under the approximations of helical symmetry and ideal MHD. The quasi-stationary evolution is governed by the first law of thermodynamics for helically symmetric systems, which is generalized to include magnetic fields. T…
▽ More
We report work towards a relativistic formulation for modeling strongly magnetized neutron stars, rotating or in a close circular orbit around another neutron star or black hole, under the approximations of helical symmetry and ideal MHD. The quasi-stationary evolution is governed by the first law of thermodynamics for helically symmetric systems, which is generalized to include magnetic fields. The formulation involves an iterative scheme for solving the Einstein-Maxwell and relativistic MHD-Euler equations numerically. The resulting configurations for binary systems could be used as self-consistent initial data for studying their inspiral and merger.
△ Less
Submitted 6 November, 2011;
originally announced November 2011.
-
GYOTO: a new general relativistic ray-tracing code
Authors:
Frederic H. Vincent,
Thibaut Paumard,
Eric Gourgoulhon,
Guy Perrin
Abstract:
GYOTO, a general relativistic ray-tracing code, is presented. It aims at computing images of astronomical bodies in the vicinity of compact objects, as well as trajectories of massive bodies in relativistic environments. This code is capable of integrating the null and timelike geodesic equations not only in the Kerr metric, but also in any metric computed numerically within the 3+1 formalism of g…
▽ More
GYOTO, a general relativistic ray-tracing code, is presented. It aims at computing images of astronomical bodies in the vicinity of compact objects, as well as trajectories of massive bodies in relativistic environments. This code is capable of integrating the null and timelike geodesic equations not only in the Kerr metric, but also in any metric computed numerically within the 3+1 formalism of general relativity. Simulated images and spectra have been computed for a variety of astronomical targets, such as a moving star or a toroidal accretion structure. The underlying code is open source and freely available. It is user-friendly, quickly handled and very modular so that extensions are easy to integrate. Custom analytical metrics and astronomical targets can be implemented in C++ plug-in extensions independent from the main code.
△ Less
Submitted 22 September, 2011;
originally announced September 2011.
-
Magnetohydrodynamics in stationary and axisymmetric spacetimes: a fully covariant approach
Authors:
Eric Gourgoulhon,
Charalampos Markakis,
Koji Uryu,
Yoshiharu Eriguchi
Abstract:
A fully geometrical treatment of general relativistic magnetohydrodynamics (GRMHD) is developed under the hypotheses of perfect conductivity, stationarity and axisymmetry. The spacetime is not assumed to be circular, which allows for greater generality than the Kerr-type spacetimes usually considered in GRMHD. Expressing the electromagnetic field tensor solely in terms of three scalar fields relat…
▽ More
A fully geometrical treatment of general relativistic magnetohydrodynamics (GRMHD) is developed under the hypotheses of perfect conductivity, stationarity and axisymmetry. The spacetime is not assumed to be circular, which allows for greater generality than the Kerr-type spacetimes usually considered in GRMHD. Expressing the electromagnetic field tensor solely in terms of three scalar fields related to the spacetime symmetries, we generalize previously obtained results in various directions. In particular, we present the first relativistic version of the Soloviev transfield equation, subcases of which lead to fully covariant versions of the Grad-Shafranov equation and of the Stokes equation in the hydrodynamical limit. We have also derived, as another subcase of the relativistic Soloviev equation, the equation governing magnetohydrodynamical equilibria with purely toroidal magnetic fields in stationary and axisymmetric spacetimes.
△ Less
Submitted 4 May, 2011; v1 submitted 18 January, 2011;
originally announced January 2011.
-
Thermodynamics of magnetized binary compact objects
Authors:
Koji Uryu,
Eric Gourgoulhon,
Charalampos Markakis
Abstract:
Binary systems of compact objects with electromagnetic field are modeled by helically symmetric Einstein-Maxwell spacetimes with charged and magnetized perfect fluids. Previously derived thermodynamic laws for helically-symmetric perfect-fluid spacetimes are extended to include the electromagnetic fields, and electric currents and charges; the first law is written as a relation between the change…
▽ More
Binary systems of compact objects with electromagnetic field are modeled by helically symmetric Einstein-Maxwell spacetimes with charged and magnetized perfect fluids. Previously derived thermodynamic laws for helically-symmetric perfect-fluid spacetimes are extended to include the electromagnetic fields, and electric currents and charges; the first law is written as a relation between the change in the asymptotic Noether charge $\dl Q$ and the changes in the area and electric charge of black holes, and in the vorticity, baryon rest mass, entropy, charge and magnetic flux of the magnetized fluid. Using the conservation laws of the circulation of magnetized flow found by Bekenstein and Oron for the ideal magnetohydrodynamic (MHD) fluid, and also for the flow with zero conducting current, we show that, for nearby equilibria that conserve the quantities mentioned above, the relation $\dl Q=0$ is satisfied. We also discuss a formulation for computing numerical solutions of magnetized binary compact objects in equilibrium with emphasis on a first integral of the ideal MHD-Euler equation.
△ Less
Submitted 21 October, 2010;
originally announced October 2010.
-
An introduction to the theory of rotating relativistic stars
Authors:
Eric Gourgoulhon
Abstract:
These lecture notes are intended to introduce the theory of rotating stars in general relativity. The focus is put on the theoretical foundations, with a detailed discussion of the spacetime symmetries, the choice of coordinates and the derivation of the equations of structure from the Einstein equation. The global properties of rotating stars (mass, angular momentum, redshifts, orbits, etc.) are…
▽ More
These lecture notes are intended to introduce the theory of rotating stars in general relativity. The focus is put on the theoretical foundations, with a detailed discussion of the spacetime symmetries, the choice of coordinates and the derivation of the equations of structure from the Einstein equation. The global properties of rotating stars (mass, angular momentum, redshifts, orbits, etc.) are also introduced.
△ Less
Submitted 20 September, 2011; v1 submitted 25 March, 2010;
originally announced March 2010.
-
Mass and Angular Momentum in General Relativity
Authors:
J. L. Jaramillo,
E. Gourgoulhon
Abstract:
We present an introduction to mass and angular momentum in General Relativity. After briefly reviewing energy-momentum for matter fields, first in the flat Minkowski case (Special Relativity) and then in curved spacetimes with or without symmetries, we focus on the discussion of energy-momentum for the gravitational field. We illustrate the difficulties rooted in the Equivalence Principle for de…
▽ More
We present an introduction to mass and angular momentum in General Relativity. After briefly reviewing energy-momentum for matter fields, first in the flat Minkowski case (Special Relativity) and then in curved spacetimes with or without symmetries, we focus on the discussion of energy-momentum for the gravitational field. We illustrate the difficulties rooted in the Equivalence Principle for defining a local energy-momentum density for the gravitational field. This leads to the understanding of gravitational energy-momentum and angular momentum as non-local observables that make sense, at best, for extended domains of spacetime. After introducing Komar quantities associated with spacetime symmetries, it is shown how total energy-momentum can be unambiguously defined for isolated systems, providing fundamental tests for the internal consistency of General Relativity as well as setting the conceptual basis for the understanding of energy loss by gravitational radiation. Finally, several attempts to formulate quasi-local notions of mass and angular momentum associated with extended but finite spacetime domains are presented, together with some illustrations of the relations between total and quasi-local quantities in the particular context of black hole spacetimes. This article is not intended to be a rigorous and exhaustive review of the subject, but rather an invitation to the topic for non-experts. In this sense we follow essentially the expositions in Szabados 2004, Gourgoulhon 2007, Poisson 2004 and Wald 84, and refer the reader interested in further developments to the existing literature, in particular to the excellent and comprehensive review by Szabados (2004).
△ Less
Submitted 30 April, 2010; v1 submitted 29 January, 2010;
originally announced January 2010.
-
Non-conformally flat initial data for binary compact objects
Authors:
Koji Uryu,
Francois Limousin,
John L. Friedman,
Eric Gourgoulhon,
Masaru Shibata
Abstract:
A new method is described for constructing initial data for a binary neutron-star (BNS) system in quasi-equilibrium circular orbit. Two formulations for non-conformally flat data, waveless (WL) and near-zone helically symmetric (NHS), are introduced; in each formulation, the Einstein-Euler system, written in 3+1 form on an asymptotically flat spacelike hypersurface, is exactly solved for all met…
▽ More
A new method is described for constructing initial data for a binary neutron-star (BNS) system in quasi-equilibrium circular orbit. Two formulations for non-conformally flat data, waveless (WL) and near-zone helically symmetric (NHS), are introduced; in each formulation, the Einstein-Euler system, written in 3+1 form on an asymptotically flat spacelike hypersurface, is exactly solved for all metric components, including the spatially non-conformally flat potentials, and for irrotational flow. A numerical method applicable to both formulations is explained with an emphasis on the imposition of a spatial gauge condition. Results are shown for solution sequences of irrotational BNS with matter approximated by parametrized equations of state that use a few segments of polytropic equations of state. The binding energy and total angular momentum of solution sequences computed within the conformally flat -- Isenberg-Wilson-Mathews (IWM) -- formulation are closer to those of the third post-Newtonian (3PN) two point particles up to the closest orbits, for the more compact stars, whereas sequences resulting from the WL/NHS formulations deviate from the 3PN curve even more for the sequences with larger compactness. We think it likely that this correction reflects an overestimation in the IWM formulation as well as in the 3PN formula, by $\sim 1$ cycle in the gravitational wave phase during the last several orbits. The work suggests that imposing spatial conformal flatness results in an underestimate of the quadrupole deformation of the components of binary neutron-star systems in the last few orbits prior to merger.
△ Less
Submitted 8 October, 2009; v1 submitted 5 August, 2009;
originally announced August 2009.
-
Improved constrained scheme for the Einstein equations: An approach to the uniqueness issue
Authors:
Isabel Cordero-Carrión,
Pablo Cerdá-Durán,
Harald Dimmelmeier,
José Luis Jaramillo,
Jérôme Novak,
Eric Gourgoulhon
Abstract:
Uniqueness problems in the elliptic sector of constrained formulations of Einstein equations have a dramatic effect on the physical validity of some numerical solutions, for instance when calculating the spacetime of very compact stars or nascent black holes. The fully constrained formulation (FCF) proposed by Bonazzola, Gourgoulhon, Grandclément, and Novak is one of these formulations. It conta…
▽ More
Uniqueness problems in the elliptic sector of constrained formulations of Einstein equations have a dramatic effect on the physical validity of some numerical solutions, for instance when calculating the spacetime of very compact stars or nascent black holes. The fully constrained formulation (FCF) proposed by Bonazzola, Gourgoulhon, Grandclément, and Novak is one of these formulations. It contains, as a particular case, the approximation of the conformal flatness condition (CFC) which, in the last ten years, has been used in many astrophysical applications. The elliptic part of the FCF basically shares the same differential operators as the elliptic equations in CFC scheme. We present here a reformulation of the elliptic sector of CFC that has the fundamental property of overcoming the local uniqueness problems. The correct behavior of our new formulation is confirmed by means of a battery of numerical simulations. Finally, we extend these ideas to FCF, complementing the mathematical analysis carried out in previous studies.
△ Less
Submitted 23 January, 2009; v1 submitted 13 September, 2008;
originally announced September 2008.
-
New theoretical approaches to black holes
Authors:
Eric Gourgoulhon,
Jose Luis Jaramillo
Abstract:
Quite recently, some new mathematical approaches to black holes have appeared in the literature. They do not rely on the classical concept of event horizon -- which is very global, but on the local concept of hypersurfaces foliated by trapped surfaces. After a brief introduction to these new horizons, we focus on a viscous fluid analogy that can be developed to describe their dynamics, in a fash…
▽ More
Quite recently, some new mathematical approaches to black holes have appeared in the literature. They do not rely on the classical concept of event horizon -- which is very global, but on the local concept of hypersurfaces foliated by trapped surfaces. After a brief introduction to these new horizons, we focus on a viscous fluid analogy that can be developed to describe their dynamics, in a fashion similar to the membrane paradigm introduced for event horizons in the seventies, but with a significant change of sign of the bulk viscosity.
△ Less
Submitted 20 March, 2008; v1 submitted 20 March, 2008;
originally announced March 2008.
-
Mathematical Issues in a Fully-Constrained Formulation of Einstein Equations
Authors:
I. Cordero-Carrión,
J. M. Ibáñez,
E. Gourgoulhon,
J. L. Jaramillo,
J. Novak
Abstract:
Bonazzola, Gourgoulhon, Grandclément, and Novak [Phys. Rev. D {\bf 70}, 104007 (2004)] proposed a new formulation for 3+1 numerical relativity. Einstein equations result, according to that formalism, in a coupled elliptic-hyperbolic system. We have carried out a preliminary analysis of the mathematical structure of that system, in particular focusing on the equations governing the evolution for…
▽ More
Bonazzola, Gourgoulhon, Grandclément, and Novak [Phys. Rev. D {\bf 70}, 104007 (2004)] proposed a new formulation for 3+1 numerical relativity. Einstein equations result, according to that formalism, in a coupled elliptic-hyperbolic system. We have carried out a preliminary analysis of the mathematical structure of that system, in particular focusing on the equations governing the evolution for the deviation of a conformal metric from a flat fiducial one. The choice of a Dirac's gauge for the spatial coordinates guarantees the mathematical characterization of that system as a (strongly) hyperbolic system of conservation laws. In the presence of boundaries, this characterization also depends on the boundary conditions for the shift vector in the elliptic subsystem. This interplay between the hyperbolic and elliptic parts of the complete evolution system is used to assess the prescription of inner boundary conditions for the hyperbolic part when using an excision approach to black hole spacetime evolutions.
△ Less
Submitted 8 April, 2008; v1 submitted 21 February, 2008;
originally announced February 2008.
-
From Geometry to Numerics: interdisciplinary aspects in mathematical and numerical relativity
Authors:
J. L. Jaramillo,
J. A. Valiente Kroon,
E. Gourgoulhon
Abstract:
This article reviews some aspects in the current relationship between mathematical and numerical General Relativity. Focus is placed on the description of isolated systems, with a particular emphasis on recent developments in the study of black holes. Ideas concerning asymptotic flatness, the initial value problem, the constraint equations, evolution formalisms, geometric inequalities and quasi-…
▽ More
This article reviews some aspects in the current relationship between mathematical and numerical General Relativity. Focus is placed on the description of isolated systems, with a particular emphasis on recent developments in the study of black holes. Ideas concerning asymptotic flatness, the initial value problem, the constraint equations, evolution formalisms, geometric inequalities and quasi-local black hole horizons are discussed on the light of the interaction between numerical and mathematical relativists.
△ Less
Submitted 30 January, 2008; v1 submitted 14 December, 2007;
originally announced December 2007.
-
Shells around black holes: the effect of freely specifiable quantities in Einstein's constraint equations
Authors:
Keith Matera,
Thomas W. Baumgarte,
Eric Gourgoulhon
Abstract:
We solve Einstein's constraint equations in the conformal thin-sandwich decomposition to model thin shells of non-interacting particles in circular orbit about a non-rotating black hole. We use these simple models to explore the effects of some of the freely specifiable quantities in this decomposition on the physical content of the solutions. Specifically, we adopt either maximal slicing or Ker…
▽ More
We solve Einstein's constraint equations in the conformal thin-sandwich decomposition to model thin shells of non-interacting particles in circular orbit about a non-rotating black hole. We use these simple models to explore the effects of some of the freely specifiable quantities in this decomposition on the physical content of the solutions. Specifically, we adopt either maximal slicing or Kerr-Schild slicing, and make different choices for the value of the lapse on the black hole horizon. For one particular choice of these quantities the resulting equations can be solved analytically; for all others we construct numerical solutions. We find that these different choices have no effect on our solutions when they are expressed in terms of gauge-invariant quantities.
△ Less
Submitted 13 November, 2007;
originally announced November 2007.
-
Trapping Horizons as inner boundary conditions for black hole spacetimes
Authors:
J. L. Jaramillo,
E. Gourgoulhon,
I. Cordero-Carrion,
J. M. Ibanez
Abstract:
We present a set of inner boundary conditions for the numerical construction of dynamical black hole space-times, when employing a 3+1 constrained evolution scheme and an excision technique. These inner boundary conditions are heuristically motivated by the dynamical trapping horizon framework and are enforced in an elliptic subsystem of the full Einstein equation. In the stationary limit they r…
▽ More
We present a set of inner boundary conditions for the numerical construction of dynamical black hole space-times, when employing a 3+1 constrained evolution scheme and an excision technique. These inner boundary conditions are heuristically motivated by the dynamical trapping horizon framework and are enforced in an elliptic subsystem of the full Einstein equation. In the stationary limit they reduce to existing isolated horizon boundary conditions. A characteristic analysis completes the discussion of inner boundary conditions for the radiative modes.
△ Less
Submitted 18 September, 2007;
originally announced September 2007.
-
Construction of initial data for 3+1 numerical relativity
Authors:
Eric Gourgoulhon
Abstract:
This lecture is devoted to the problem of computing initial data for the Cauchy problem of 3+1 general relativity. The main task is to solve the constraint equations. The conformal technique, introduced by Lichnerowicz and enhanced by York, is presented. Two standard methods, the conformal transverse-traceless one and the conformal thin sandwich, are discussed and illustrated by some simple exam…
▽ More
This lecture is devoted to the problem of computing initial data for the Cauchy problem of 3+1 general relativity. The main task is to solve the constraint equations. The conformal technique, introduced by Lichnerowicz and enhanced by York, is presented. Two standard methods, the conformal transverse-traceless one and the conformal thin sandwich, are discussed and illustrated by some simple examples. Finally a short review regarding initial data for binary systems (black holes and neutron stars) is given.
△ Less
Submitted 7 November, 2007; v1 submitted 2 April, 2007;
originally announced April 2007.
-
3+1 Formalism and Bases of Numerical Relativity
Authors:
Eric Gourgoulhon
Abstract:
These lecture notes provide some introduction to the 3+1 formalism of general relativity, which is the foundation of most modern numerical relativity. The text is rather self-contained, with detailed calculations and numerous examples. Contents: 1. Introduction, 2. Geometry of hypersurfaces, 3. Geometry of foliations, 4. 3+1 decomposition of Einstein equation, 5. 3+1 equations for matter and ele…
▽ More
These lecture notes provide some introduction to the 3+1 formalism of general relativity, which is the foundation of most modern numerical relativity. The text is rather self-contained, with detailed calculations and numerous examples. Contents: 1. Introduction, 2. Geometry of hypersurfaces, 3. Geometry of foliations, 4. 3+1 decomposition of Einstein equation, 5. 3+1 equations for matter and electromagnetic field, 6. Conformal decomposition, 7. Asymptotic flatness and global quantities, 8. The initial data problem, 9. Choice of foliation and spatial coordinates, 10. Evolution schemes.
△ Less
Submitted 6 March, 2007;
originally announced March 2007.
-
Area evolution, bulk viscosity and entropy principles for dynamical horizons
Authors:
Eric Gourgoulhon,
Jose Luis Jaramillo
Abstract:
We derive from Einstein equation an evolution law for the area of a trapping or dynamical horizon. The solutions to this differential equation show a causal behavior. Moreover, in a viscous fluid analogy, the equation can be interpreted as an energy balance law, yielding to a positive bulk viscosity. These two features contrast with the event horizon case, where the non-causal evolution of the a…
▽ More
We derive from Einstein equation an evolution law for the area of a trapping or dynamical horizon. The solutions to this differential equation show a causal behavior. Moreover, in a viscous fluid analogy, the equation can be interpreted as an energy balance law, yielding to a positive bulk viscosity. These two features contrast with the event horizon case, where the non-causal evolution of the area and the negative bulk viscosity require teleological boundary conditions. This reflects the local character of trapping horizons as opposed to event horizons. Interpreting the area as the entropy, we propose to use an area/entropy evolution principle to select a unique dynamical horizon and time slicing in the Cauchy evolution of an initial marginally trapped surface.
△ Less
Submitted 12 September, 2006; v1 submitted 12 July, 2006;
originally announced July 2006.
-
Viscosity driven instability in rotating relativistic stars
Authors:
Motoyuki Saijo,
Eric Gourgoulhon
Abstract:
We investigate the viscosity driven instability in rotating relativistic stars by means of an iterative approach. We focus on polytropic rotating equilibrium stars and impose an m=2 perturbationin the lapse. We vary both the stiffness of the equation of state and the compactness of the star to study those effects on the value of the threshold. For a uniformly rotating star, the criterion T/W, wh…
▽ More
We investigate the viscosity driven instability in rotating relativistic stars by means of an iterative approach. We focus on polytropic rotating equilibrium stars and impose an m=2 perturbationin the lapse. We vary both the stiffness of the equation of state and the compactness of the star to study those effects on the value of the threshold. For a uniformly rotating star, the criterion T/W, where T is the rotational kinetic energy and W is the gravitational binding energy, mainly depends on the compactness of the star and takes values around 0.13 ~ 0.16, which differ slightly from that of Newtonian incompressible stars (~ 0.14). For differentially rotating stars, the critical value of T/W is found to span the range 0.17 - 0.25. This is significantly larger than the uniformly rotating case with the same compactness of the star. Finally we discuss a possibility of detecting gravitational waves from viscosity driven instability with ground-based interferometers.
△ Less
Submitted 16 October, 2006; v1 submitted 22 June, 2006;
originally announced June 2006.
-
An introduction to relativistic hydrodynamics
Authors:
Eric Gourgoulhon
Abstract:
This lecture provides some introduction to perfect fluid dynamics within the framework of general relativity. The presentation is based on the Carter-Lichnerowicz approach. It has the advantage over the more traditional approach of leading very straightforwardly to important conservation laws, such as the relativistic generalizations of Bernoulli's theorem or Kelvin's circulation theorem. It als…
▽ More
This lecture provides some introduction to perfect fluid dynamics within the framework of general relativity. The presentation is based on the Carter-Lichnerowicz approach. It has the advantage over the more traditional approach of leading very straightforwardly to important conservation laws, such as the relativistic generalizations of Bernoulli's theorem or Kelvin's circulation theorem. It also permits to get easily first integrals of motion which are particularly useful for computing equilibrium configurations of relativistic stars in rotation or in binary systems. The presentation is relatively self-contained and does not require any a priori knowledge of general relativity. In particular, the three types of derivatives involved in relativistic hydrodynamics are introduced in detail: this concerns the Lie, exterior and covariant derivatives.
△ Less
Submitted 5 March, 2006;
originally announced March 2006.
-
Binary neutron stars: Equilibrium models beyond spatial conformal flatness
Authors:
Koji Uryu,
Francois Limousin,
John L. Friedman,
Eric Gourgoulhon,
Masaru Shibata
Abstract:
Equilibria of binary neutron stars in close circular orbits are computed numerically in a waveless formulation: The full Einstein-relativistic-Euler system is solved on an initial hypersurface to obtain an asymptotically flat form of the 4-metric and an extrinsic curvature whose time derivative vanishes in a comoving frame. Two independent numerical codes are developed, and solution sequences th…
▽ More
Equilibria of binary neutron stars in close circular orbits are computed numerically in a waveless formulation: The full Einstein-relativistic-Euler system is solved on an initial hypersurface to obtain an asymptotically flat form of the 4-metric and an extrinsic curvature whose time derivative vanishes in a comoving frame. Two independent numerical codes are developed, and solution sequences that model inspiraling binary neutron stars during the final several orbits are successfully computed. The binding energy of the system near its final orbit deviates from earlier results of third post-Newtonian and of spatially conformally flat calculations. The new solutions may serve as initial data for merger simulations and as members of quasiequilibrium sequences to generate gravitational wave templates, and may improve estimates of the gravitational-wave cutoff frequency set by the last inspiral orbit.
△ Less
Submitted 21 October, 2006; v1 submitted 25 November, 2005;
originally announced November 2005.
-
A generalized Damour-Navier-Stokes equation applied to trapping horizons
Authors:
Eric Gourgoulhon
Abstract:
An identity is derived from Einstein equation for any hypersurface H which can be foliated by spacelike two-dimensional surfaces. In the case where the hypersurface is null, this identity coincides with the two-dimensional Navier-Stokes-like equation obtained by Damour in the membrane approach to a black hole event horizon. In the case where H is spacelike or null and the 2-surfaces are marginal…
▽ More
An identity is derived from Einstein equation for any hypersurface H which can be foliated by spacelike two-dimensional surfaces. In the case where the hypersurface is null, this identity coincides with the two-dimensional Navier-Stokes-like equation obtained by Damour in the membrane approach to a black hole event horizon. In the case where H is spacelike or null and the 2-surfaces are marginally trapped, this identity applies to Hayward's trapping horizons and to the related dynamical horizons recently introduced by Ashtekar and Krishnan. The identity involves a normal fundamental form (normal connection 1-form) of the 2-surface, which can be viewed as a generalization to non-null hypersurfaces of the Hajicek 1-form used by Damour. This 1-form is also used to define the angular momentum of the horizon. The generalized Damour-Navier-Stokes equation leads then to a simple evolution equation for the angular momentum.
△ Less
Submitted 19 October, 2005; v1 submitted 31 July, 2005;
originally announced August 2005.
-
A 3+1 perspective on null hypersurfaces and isolated horizons
Authors:
E. Gourgoulhon,
J. L. Jaramillo
Abstract:
The isolated horizon formalism recently introduced by Ashtekar et al. aims at providing a quasi-local concept of a black hole in equilibrium in an otherwise possibly dynamical spacetime. In this formalism, a hierarchy of geometrical structures is constructed on a null hypersurface. On the other side, the 3+1 formulation of general relativity provides a powerful setting for studying the spacetime…
▽ More
The isolated horizon formalism recently introduced by Ashtekar et al. aims at providing a quasi-local concept of a black hole in equilibrium in an otherwise possibly dynamical spacetime. In this formalism, a hierarchy of geometrical structures is constructed on a null hypersurface. On the other side, the 3+1 formulation of general relativity provides a powerful setting for studying the spacetime dynamics, in particular gravitational radiation from black hole systems. Here we revisit the kinematics and dynamics of null hypersurfaces by making use of some 3+1 slicing of spacetime. In particular, the additional structures induced on null hypersurfaces by the 3+1 slicing permit a natural extension to the full spacetime of geometrical quantities defined on the null hypersurface. This 4-dimensional point of view facilitates the link between the null and spatial geometries. We proceed by reformulating the isolated horizon structure in this framework. We also reformulate previous works, such as Damour's black hole mechanics, and make the link with a previous 3+1 approach of black hole horizon, namely the membrane paradigm. We explicit all geometrical objects in terms of 3+1 quantities, putting a special emphasis on the conformal 3+1 formulation. This is in particular relevant for the initial data problem of black hole spacetimes for numerical relativity. Illustrative examples are provided by considering various slicings of Schwarzschild and Kerr spacetimes.
△ Less
Submitted 26 October, 2005; v1 submitted 29 March, 2005;
originally announced March 2005.
-
The final phase of inspiral of neutron stars: realistic equations of state
Authors:
Dorota Gondek-Rosinska,
Michal Bejger,
Tomek Bulik,
Eric Gourgoulhon,
Pawel Haensel,
Francois Limousin,
Keisuke Taniguchi,
Leszek Zdunik
Abstract:
Coalescing compact star binaries are expected to be among the strongest sources of gravitational radiation to be seen by laser interferometers. We present calculations of the final phase of inspiral of equal mass irrotational neutron star binaries and strange quark star binaries. Six types of equations of state at zero temperature are used - three realistic nuclear equations of state of various…
▽ More
Coalescing compact star binaries are expected to be among the strongest sources of gravitational radiation to be seen by laser interferometers. We present calculations of the final phase of inspiral of equal mass irrotational neutron star binaries and strange quark star binaries. Six types of equations of state at zero temperature are used - three realistic nuclear equations of state of various softness and three different MIT bag models of strange quark matter. We study the precoalescing stage within the Isenberg-Wilson-Mathews approximation of general relativity using a multidomain spectral method. The gravitational-radiation driven evolution of the binary system is approximated by a sequence of quasi-equilibrium configurations at fixed baryon number and decreasing separation. We find that the innermost stable circular orbit (ISCO) is given by an orbital instability for binary strange quark stars and by the mass-shedding limit for neutron star binaries. The gravitational wave frequency at the ISCO, which marks the end of the inspiral phase, is found to be around 1100-1460 Hz for two 1.35 solar masses irrotational strange stars described by the MIT bag model and between 800 Hz and 1230 Hz for neutron stars.
△ Less
Submitted 10 January, 2007; v1 submitted 2 December, 2004;
originally announced December 2004.
-
Quasiequilibrium sequences of binary strange quark stars in general relativity
Authors:
Francois Limousin,
Dorota Gondek-Rosinska,
Eric Gourgoulhon
Abstract:
Inspiraling compact binaries are expected to be the strongest sources of gravitational waves for VIRGO, LIGO and other laser interferometers. We present the first computations of quasi-equilibrium sequences of compact binaries containing two strange quark stars (which are currently considered as a possible alternative to neutron stars). We study a precoalescing stage in the conformal flatness ap…
▽ More
Inspiraling compact binaries are expected to be the strongest sources of gravitational waves for VIRGO, LIGO and other laser interferometers. We present the first computations of quasi-equilibrium sequences of compact binaries containing two strange quark stars (which are currently considered as a possible alternative to neutron stars). We study a precoalescing stage in the conformal flatness approximation of general relativity using a multidomain spectral method. A hydrodynamical treatment is performed under the assumption that the flow is irrotational.
△ Less
Submitted 30 November, 2004;
originally announced November 2004.