Researcher profile

Thomas W. Baumgarte

Thomas W. Baumgarte contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 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$.

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.

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.

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.

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.

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.

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.

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.

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.