BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Thomas W. Baumgarte

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

ResearcherUnclaimed source record
gr-qcastro-phastro-ph.SRastro-ph.HEphysics.comp-ph

Catalog footprint

What is connected

28works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2023arXiv

Bona-Massó slices of Reissner-Nordström spacetimes

Motivated by recent numerical relativity simulations of charged black holes and their interactions, we explore the properties of common slicing conditions in Reissner-Nordström spacetimes. Specifically, we consider different choices for the so-called Bona-Massó function, and construct static and spherically symmetric slices of the Reissner-Nordström spacetime satisfying the corresponding slicing conditions. For some of these functions the construction is entirely analytical, while for others we use numerical root-finding to solve quartic equations. Our solutions are parameterized by the charge-to-mass ratio $λ= Q/M$ and approach a unique slice, independent of the Bona-Massó functions considered here, in the extremal limit $λ\to 1$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Bona-Masso slicing conditions and the lapse close to black-hole punctures

We consider several families of functions $f(α)$ that appear in the Bona-Masso slicing condition for the lapse function $α$. Focusing on spherically symmetric and time-independent slices we apply these conditions to the Schwarzschild spacetime in order to construct analytical expressions for the lapse $α$ in terms of the areal radius $R$. We then transform to isotropic coordinates and determine the dependence of $α$ on the isotropic radius $r$ in the vicinity of the black-hole puncture. We propose generalizations of previously considered functions $f(α)$ for which, to leading order, the lapse is proportional to $r$ rather than a non-integer power of $r$. We also perform dynamical simulations in spherical symmetry and demonstrate advantages of the above choices in numerical simulations employing spectral methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Shock-avoiding slicing conditions: tests and calibrations

While the 1+log slicing condition has been extremely successful in numerous numerical relativity simulations, it is also known to develop "gauge-shocks" in some examples. Alternative "shock-avoiding" slicing conditions suggested by Alcubierre prevent these pathologies in those examples, but have not yet been explored and tested very broadly. In this paper we compare the performance of shock-avoiding slicing conditions with those of 1+log slicing for a number of "text-book" problems, including black holes and relativistic stars. While, in some simulations, the shock-avoiding slicing conditions feature some unusual properties and lead to more "gauge-dynamics" than the 1+log slicing condition, we find that they perform quite similarly in terms of stability and accuracy, and hence provide a very viable alternative to 1+log slicing.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Accretion onto black holes inside neutron stars with piecewise-polytropic equations of state: analytic and numerical treatments

We consider spherically symmetric accretion onto a small, possibly primordial, black hole residing at the center of a neutron star governed by a cold nuclear equation of state (EOS). We generalize the relativistic Bondi solution for such EOSs, approximated by piecewise polytropes, and thereby obtain analytical expressions for the steady-state matter profiles and accretion rates. We compare these rates with those found by time-dependent, general relativistic hydrodynamical simulations upon relaxation and find excellent agreement. We consider several different candidate EOSs, neutron star masses and central densities and find that the accretion rates vary only little, resulting in an accretion rate that depends primarily on the black hole mass, and only weakly on the properties of the neutron star.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical relativity in spherical coordinates: A new dynamical spacetime and general relativistic MHD evolution framework for the Einstein Toolkit

We present SphericalNR, a new framework for the publicly available Einstein Toolkit that numerically solves the Einstein field equations coupled to the equations of general relativistic magnetohydrodynamics (GRMHD) in a 3+1 split of spacetime in spherical coordinates without symmetry assumptions. The spacetime evolution is performed using reference-metric versions of either the Baumgarte-Shapiro-Shibata-Nakamura equations or the fully covariant and conformal Z4 system with constraint damping. We have developed a reference-metric version of the Valencia formulation of GRMHD with a vector potential method, guaranteeing the absence of magnetic monopoles during the evolution. In our framework, every dynamical field (both spacetime and matter) is evolved using its components in an orthonormal basis with respect to the spherical reference-metric. Furthermore, all geometric information about the spherical coordinate system is encoded in source terms appearing in the evolution equations. This allows for the straightforward extension of Cartesian high-resolution shock-capturing finite volume codes to use spherical coordinates with our framework. To this end, we have adapted GRHydro, a Cartesian finite volume GRMHD code already available in the Einstein Toolkit, to use spherical coordinates. We present the full evolution equations of the framework, as well as details of its implementation in the Einstein Toolkit. We validate SphericalNR by demonstrating it passes a variety of challenging code tests in static and dynamical spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Critical phenomena in the aspherical gravitational collapse of radiation fluids

We study critical phenomena in the gravitational collapse of a radiation fluid. We perform numerical simulations in both spherical symmetry and axisymmetry, and observe critical scaling in both supercritical evolutions, which lead to the formation of a black hole, and subcritical evolutions, in which case the fluid disperses to infinity and leaves behind flat space. We identify the critical solution in spherically symmetric collapse, find evidence for its universality, and study the approach to this critical solution in the absence of spherical symmetry. For the cases that we consider, aspherical deviations from the spherically symmetric critical solution decay in damped oscillations in a manner that is consistent with the behavior found by Gundlach in perturbative calculations. Our simulations are performed with an unconstrained evolution code, implemented in spherical polar coordinates, and adopting "moving-puncture" coordinates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical Relativity in Spherical Polar Coordinates: Off-center Simulations

We have recently presented a new approach for numerical relativity simulations in spherical polar coordinates, both for vacuum and for relativistic hydrodynamics. Our approach is based on a reference-metric formulation of the BSSN equations, a factoring of all tensor components, as well as a partially implicit Runge-Kutta method, and does not rely on a regularization of the equations, nor does it make any assumptions about the symmetry across the origin. In order to demonstrate this feature we present here several off-centered simulations, including simulations of single black holes and neutron stars whose center is placed away from the origin of the coordinate system, as well as the asymmetric head-on collision of two black holes. We also revisit our implementation of relativistic hydrodynamics and demonstrate that a reference-metric formulation of hydrodynamics together with a factoring of all tensor components avoids problems related to the coordinate singularities at the origin and on the axes. As a particularly demanding test we present results for a shock wave propagating through the origin of the spherical polar coordinate system.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Trumpet Slices in Kerr Spacetimes

We introduce a new time-independent family of analytical coordinate systems for the Kerr spacetime representing rotating black holes. We also propose a (2+1)+1 formalism for the characterization of trumpet geometries. Applying this formalism to our new family of coordinate systems we identify, for the first time, analytical and stationary trumpet slices for general rotating black holes, even for charged black holes in the presence of a cosmological constant. We present results for metric functions in this slicing and analyze the geometry of the rotating trumpet surface.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Simple Family of Analytical Trumpet Slices of the Schwarzschild Spacetime

We describe a simple family of analytical coordinate systems for the Schwarzschild spacetime. The coordinates penetrate the horizon smoothly and are spatially isotropic. Spatial slices of constant coordinate time $t$ feature a trumpet geometry with an asymptotically cylindrical end inside the horizon at a prescribed areal radius $R_0$ (with $0<R_{0}\leq M$) that serves as the free parameter for the family. The slices also have an asymptotically flat end at spatial infinity. In the limit $R_{0}=0$ the spatial slices lose their trumpet geometry and become flat -- in this limit, our coordinates reduce to Painlevé-Gullstrand coordinates.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fully covariant and conformal formulation of the Z4 system in a reference-metric approach: comparison with the BSSN formulation in spherical symmetry

We adopt a reference-metric approach to generalize a covariant and conformal version of the Z4 system of the Einstein equations. We refer to the resulting system as ``fully covariant and conformal", or fCCZ4 for short, since it is well suited for curvilinear as well as Cartesian coordinates. We implement this fCCZ4 formalism in spherical polar coordinates under the assumption of spherical symmetry using a partially-implicit Runge-Kutta (PIRK) method and show that our code can evolve both vacuum and non-vacuum spacetimes without encountering instabilities. Our method does not require regularization of the equations to handle coordinate singularities, nor does it depend on constraint-preserving outer boundary conditions, nor does it need any modifications of the equations for evolutions of black holes. We perform several tests and compare the performance of the fCCZ4 system, for different choices of certain free parameters, with that of BSSN. Confirming earlier results we find that, for an optimal choice of these parameters, and for neutron-star spacetimes, the violations of the Hamiltonian constraint can be between 1 and 3 orders of magnitude smaller in the fCCZ4 system than in the BSSN formulation. For black-hole spacetimes, on the other hand, any advantages of fCCZ4 over BSSN are less evident.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

General relativistic hydrodynamics in curvilinear coordinates

In this paper we report on what we believe is the first successful implementation of relativistic hydrodynamics, coupled to dynamical spacetimes, in spherical polar coordinates without symmetry assumptions. We employ a high-resolution shock-capturing scheme, which requires that the equations be cast in flux-conservative form. One example of such a form is the :Valencia" formulation, which has been adopted in numerous applications, in particular in Cartesian coordinates. Here we generalize this formulation to allow for a reference-metric approach, which provides a natural framework for calculations in curvilinear coordinates. In spherical polar coordinates, for example, it allows for an analytical treatment of the singular r and sin(θ) terms that appear in the equations. We experiment with different versions of our generalized Valencia formulation in numerical implementations of relativistic hydrodynamics for both fixed and dynamical spacetimes. We consider a number of different tests -- non-rotating and rotating relativistic stars, as well as gravitational collapse to a black hole -- to demonstrate that our formulation provides a promising approach to performing fully relativistic astrophysics simulations in spherical polar coordinates.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Collapse of Nonlinear Gravitational Waves in Moving-Puncture Coordinates

We study numerical evolutions of nonlinear gravitational waves in moving-puncture coordinates. We adopt two different types of initial data -- Brill and Teukolsky waves -- and evolve them with two independent codes producing consistent results. We find that Brill data fail to produce long-term evolutions for common choices of coordinates and parameters, unless the initial amplitude is small, while Teukolsky wave initial data lead to stable evolutions, at least for amplitudes sufficiently far from criticality. The critical amplitude separates initial data whose evolutions leave behind flat space from those that lead to a black hole. For the latter we follow the interaction of the wave, the formation of a horizon, and the settling down into a time-independent trumpet geometry. We explore the differences between Brill and Teukolsky data and show that for less common choices of the parameters -- in particular negative amplitudes -- Brill data can be evolved with moving-puncture coordinates, and behave similarly to Teukolsky waves.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An alternative approach to solving the Hamiltonian constraint

Solving Einstein's constraint equations for the construction of black hole initial data requires handling the black hole singularity. Typically, this is done either with the excision method, in which the black hole interior is excised from the numerical grid, or with the puncture method, in which the singular part of the conformal factor is expressed in terms of an analytical background solution, and the Hamiltonian constraint is then solved for a correction to the background solution that, usually, is assumed to be regular everywhere. We discuss an alternative approach in which the Hamiltonian constraint is solved for an inverse power of the conformal factor. This new function remains finite everywhere, so that this approach requires neither excision nor a split into background and correction. In particular, this method can be used without modification even when the correction to the conformal factor is singular itself. We demonstrate this feature for rotating black holes in the trumpet topology.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Analytical Tendex and Vortex Fields for Perturbative Black Hole Initial Data

Tendex and vortex fields, defined by the eigenvectors and eigenvalues of the electric and magnetic parts of the Weyl curvature tensor, form the basis of a recently developed approach to visualizing spacetime curvature. In particular, this method has been proposed as a tool for interpreting results from numerical binary black hole simulations, providing a deeper insight into the physical processes governing the merger of black holes and the emission of gravitational radiation. Here we apply this approach to approximate but analytical initial data for both single boosted and binary black holes. These perturbative data become exact in the limit of small boost or large binary separation. We hope that these calculations will provide additional insight into the properties of tendex and vortex fields, and will form a useful test for future numerical calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravity darkening and brightening in binaries

We apply a von Zeipel gravity darkening model to corotating binaries to obtain a simple, analytical expression for the emergent radiative flux from a tidally distorted primary orbiting a point-mass secondary. We adopt a simple Roche model to determine the envelope structure of the primary, assumed massive and centrally condensed, and use the results to calculate the flux. As for single rotating stars, gravity darkening reduces the flux along the stellar equator of the primary, but, unlike for rotating stars, we find that gravity brightening enhances the flux in a region around the stellar poles. We identify a critical limiting separation beyond which hydrostatic equilibrium no longer is possible, whereby the flux vanishes at the point on the stellar equator of the primary facing the companion. For equal-mass binaries, the total luminosity is reduced by about 13 % when this limiting separation is reached.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Invariants for Tendex and Vortex Fields

Tendex and vortex fields, defined by the eigenvectors and eigenvalues of the electric and magnetic parts of the Weyl curvature tensor, form the basis of a recently developed approach to visualizing spacetime curvature. In analogy to electric and magnetic fields, these fields are coordinate-dependent. However, in a further analogy, we can form invariants from the tendex and vortex fields that are invariant under coordinate transformations, just as certain combinations of the electric and magnetic fields are invariant under coordinate transformations. We derive these invariants, and provide a simple, analytical demonstration for non-spherically symmetric slices of a Schwarzschild spacetime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical Relativity in Spherical Polar Coordinates: Evolution Calculations with the BSSN Formulation

In the absence of symmetry assumptions most numerical relativity simulations adopt Cartesian coordinates. While Cartesian coordinates have some desirable properties, spherical polar coordinates appear better suited for certain applications, including gravitational collapse and supernova simulations. Development of numerical relativity codes in spherical polar coordinates has been hampered by the need to handle the coordinate singularities at the origin and on the axis, for example by careful regularization of the appropriate variables. Assuming spherical symmetry and adopting a covariant version of the BSSN equations, Montero and Cordero-Carrión recently demonstrated that such a regularization is not necessary when a partially implicit Runge-Kutta (PIRK) method is used for the time evolution of the gravitational fields. Here we report on an implementation of the BSSN equations in spherical polar coordinates without any symmetry assumptions. Using a PIRK method we obtain stable simulations in three spatial dimensions without the need to regularize the origin or the axis. We perform and discuss a number of tests to assess the stability, accuracy and convergence of the code, namely weak gravitational waves, "hydro-without-hydro" evolutions of spherical and rotating relativistic stars in equilibrium, and single black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Oppenheimer-Snyder Collapse in Moving-Puncture Coordinates

Moving-puncture coordinates are commonly used in numerical simulations of black holes. Their properties for vacuum Schwarzschild black holes have been analyzed in a number of studies. The behavior of moving-puncture coordinates in spacetimes containing matter, however, is less well understood. In this paper we explore the behavior of these coordinates for Oppenheimer-Snyder collapse, i.e., the collapse of a uniform density, pressureless sphere of dust initially at rest to a black hole. Oppenheimer-Snyder collapse provides a stringent test of the singularity-avoiding properties of moving-puncture coordinates, since the singularity can form more quickly than it would for matter with pressure. Our results include analytical expressions for the matter density, lapse function, and mean curvature at early times, as well as interesting limits for later times. We also carry out numerical simulations to obtain the full solution and these show that, even in the absence of pressure, moving-puncture coordinates are able to avoid the singularity. At late times the geometry settles down to a trumpet slice of a vacuum black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Puncture black hole initial data in the conformal thin-sandwich formalism

We revisit the construction of puncture black hole initial data in the conformal thin-sandwich decomposition of Einstein's constraint equations. It has been shown previously that this approach cannot yield quasiequilibrium wormhole data, which connect two asymptotically flat spatial infinities. This argument does not apply to trumpet data, which connect the spatial infinity in one universe with the future timelike infinity of another. As a numerical demonstration we present results for a single boosted trumpet-puncture black holes, constructed in the original version of the conformal thin-sandwich formalism.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Trumpet slices of the Schwarzschild-Tangherlini spacetime

We study families of time-independent maximal and 1+log foliations of the Schwarzschild-Tangherlini spacetime, the spherically-symmetric vacuum black hole solution in D spacetime dimensions, for D >= 4. We identify special members of these families for which the spatial slices display a trumpet geometry. Using a generalization of the 1+log slicing condition that is parametrized by a constant n we recover the results of Nakao, Abe, Yoshino and Shibata in the limit of maximal slicing. We also construct a numerical code that evolves the BSSN equations for D=5 in spherical symmetry using moving-puncture coordinates, and demonstrate that these simulations settle down to the trumpet solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Merger of white dwarf-neutron star binaries: Prelude to hydrodynamic simulations in general relativity

White dwarf-neutron star binaries generate detectable gravitational radiation. We construct Newtonian equilibrium models of corotational white dwarf-neutron star (WDNS) binaries in circular orbit and find that these models terminate at the Roche limit. At this point the binary will undergo either stable mass transfer (SMT) and evolve on a secular time scale, or unstable mass transfer (UMT), which results in the tidal disruption of the WD. The path a given binary will follow depends primarily on its mass ratio. We analyze the fate of known WDNS binaries and use population synthesis results to estimate the number of LISA-resolved galactic binaries that will undergo either SMT or UMT. We model the quasistationary SMT epoch by solving a set of simple ordinary differential equations and compute the corresponding gravitational waveforms. Finally, we discuss in general terms the possible fate of binaries that undergo UMT and construct approximate Newtonian equilibrium configurations of merged WDNS remnants. We use these configurations to assess plausible outcomes of our future, fully relativistic simulations of these systems. If sufficient WD debris lands on the NS, the remnant may collapse, whereby the gravitational waves from the inspiral, merger, and collapse phases will sweep from LISA through LIGO frequency bands. If the debris forms a disk about the NS, it may fragment and form planets.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Einstein constraints: uniqueness and non-uniqueness in the conformal thin sandwich approach

We study the appearance of multiple solutions to certain decompositions of Einstein's constraint equations. Pfeiffer and York recently reported the existence of two branches of solutions for identical background data in the extended conformal thin-sandwich decomposition. We show that the Hamiltonian constraint alone, when expressed in a certain way, admits two branches of solutions with properties very similar to those found by Pfeiffer and York. We construct these two branches analytically for a constant-density star in spherical symmetry, but argue that this behavior is more general. In the case of the Hamiltonian constraint this non-uniqueness is well known to be related to the sign of one particular term, and we argue that the extended conformal thin-sandwich equations contain a similar term that causes the breakdown of uniqueness.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Quasi-equilibrium Binary Black Hole Initial Data for Dynamical Evolutions

We present a formalism for constructing quasi-equilibrium binary black hole initial data suitable for numerical evolution. We construct quasi-equilibrium models by imposing an approximate helical Killing symmetry appropriate for quasi-circular orbits. We use the sum of two Kerr-Schild metrics as our background metric, thereby improving on conformally flat backgrounds that do not accommodate rotating black holes and providing a horizon-penetrating lapse, convenient for implementing black hole excision. We set inner boundary conditions at an excision radius well inside the apparent horizon and construct these boundary conditions to incorporate the quasi-equilibrium condition and recover the solution for isolated black holes in the limit of large separation. We use our formalism both to generate initial data for binary black hole evolutions and to construct a crude quasi-equilibrium, inspiral sequence for binary black holes of fixed irreducible mass.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

One-Armed Spiral Instability in Differentially Rotating Stars

We investigate the dynamical instability of the one-armed spiral m=1 mode in differentially rotating stars by means of 3+1 hydrodynamical simulations in Newtonian gravitation. We find that both a soft equation of state and a high degree of differential rotation in the equilibrium star are necessary to excite a dynamical m=1 mode as the dominant instability at small values of the ratio of rotational kinetic to potential energy, T/|W|. We find that this spiral mode propagates outward from its point of origin near the maximum density at the center to the surface over several central orbital periods. An unstable m=1 mode triggers a secondary m=2 bar mode of smaller amplitude, and the bar mode can excite gravitational waves. As the spiral mode propagates to the surface it weakens, simultaneously damping the emitted gravitational wave signal. This behavior is in contrast to waves triggered by a dynamical m=2 bar instability, which persist for many rotation periods and decay only after a radiation-reaction damping timescale.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Effects of Differential Rotation on the Maximum Mass of Neutron Stars

The merger of binary neutron stars is likely to lead to differentially rotating remnants. In this paper we numerically construct models of differentially rotating neutron stars in general relativity and determine their maximum allowed mass. We model the stars adopting a polytropic equation of state and tabulate maximum allowed masses as a function of differential rotation and stiffness of the equation of state. We also provide a crude argument that yields a qualitative estimate of the effect of stiffness and differential rotation on the maximum allowed mass.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

On the Maximum Mass of Differentially Rotating Neutron Stars

We construct relativistic equilibrium models of differentially rotating neutron stars and show that they can support significantly more mass than their nonrotating or uniformly rotating counterparts. We dynamically evolve such ``hypermassive'' models in full general relativity and show that there do exist configurations which are dynamically stable against radial collapse and bar formation. Our results suggest that the remnant of binary neutron star coalescence may be temporarily stabilized by differential rotation, leading to delayed collapse and a delayed gravitational wave burst.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Numerical Evolution of Black Holes with a Hyperbolic Formulation of General Relativity

We describe a numerical code that solves Einstein's equations for a Schwarzschild black hole in spherical symmetry, using a hyperbolic formulation introduced by Choquet-Bruhat and York. This is the first time this formulation has been used to evolve a numerical spacetime containing a black hole. We excise the hole from the computational grid in order to avoid the central singularity. We describe in detail a causal differencing method that should allow one to stably evolve a hyperbolic system of equations in three spatial dimensions with an arbitrary shift vector, to second-order accuracy in both space and time. We demonstrate the success of this method in the spherically symmetric case.

0 citations0 reviewsNo code linkedOpen access