Source author record

Carsten Gundlach

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

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

15 published item(s)

preprint2022arXiv

Critical collapse of an axisymmetric ultrarelativistic fluid in $2+1$ dimensions

We carry out numerical simulations of the gravitational collapse of a rotating perfect fluid with the ultrarelativistic equation of state $P=κρ$, in axisymmetry in $2+1$ spacetime dimensions with $Λ<0$. We show that for $κ\lesssim 0.42$, the critical phenomena are type I and the critical solution is stationary. The picture for $κ\gtrsim 0.43$ is more delicate: for small angular momenta, we find type II phenomena and the critical solution is quasistationary, contracting adiabatically. The spin-to-mass ratio of the critical solution increases as it contracts, and hence so does that of the black hole created at the end as we fine-tune to the black-hole threshold. Forming extremal black holes is avoided because the contraction of the critical solution smoothly ends as extremality is approached.

preprint2016arXiv

The Einstein-Vlasov system in spherical symmetry: reduction of the equations of motion and classification of single-shell static solutions, in the limit of massless particles

We express the Einstein-Vlasov system in spherical symmetry in terms of a dimensionless momentum variable $z$ (radial over angular momentum). This regularises the limit of massless particles, and in that limit allows us to obtain a reduced system in independent variables $(t,r,z)$ only. Similarly, in this limit the Vlasov density function $f$ for static solutions depends on a single variable $Q$ (energy over angular momentum). This reduction allows us to show that any given static metric which has vanishing Ricci scalar, is vacuum at the centre and for $r>3M$ and obeys certain energy conditions uniquely determines a consistent $f=\bar k(Q)$ (in closed form). Vice versa, any $\bar k(Q)$ within a certain class uniquely determines a static metric (as the solution of a system of two first-order quasilinear ODEs). Hence the space of static spherically symmetric solutions of Einstein-Vlasov is locally a space of functions of one variable. For a simple 2-parameter family of functions $\bar k(Q)$, we construct the corresponding static spherically symmetric solutions, finding that their compactness is in the interval $0.7\lesssim \rm max_r(2M/r)\le 8/9$. This class of static solutions includes one that agrees with the approximately universal type-I critical solution recently found by Akbarian and Choptuik (AC) in numerical time evolutions. We speculate on what singles it out as the critical solution found by fine-tuning generic data to the collapse threshold, given that AC also found that {\em all} static solutions are one-parameter unstable and sit on the threshold of collapse.

preprint2015arXiv

Scalar field critical collapse in 2+1 dimensions

We carry out numerical experiments in the critical collapse of a spherically symmetric massless scalar field in 2+1 spacetime dimensions in the presence of a negative cosmological constant and compare them against a new theoretical model. We approximate the true critical solution as the $n=4$ Garfinkle solution, matched at the lightcone to a Vaidya-like solution, and corrected to leading order for the effect of $Λ<0$. This approximation is only $C^3$ at the lightcone and has three growing modes. We {\em conjecture} that pointwise it is a good approximation to a yet unknown true critical solution that is analytic with only one growing mode (itself approximated by the top mode of our amended Garfinkle solution). With this conjecture, we predict a Ricci-scaling exponent of $γ=8/7$ and a mass-scaling exponent of $δ=16/23$, compatible with our numerical experiments.

preprint2012arXiv

Critical phenomena at the threshold of immediate merger in binary black hole systems: the extreme mass ratio case

In numerical simulations of black hole binaries, Pretorius and Khurana [Class. Quant. Grav. {\bf 24}, S83 (2007)] have observed critical behaviour at the threshold between scattering and immediate merger. The number of orbits scales as $n\simeq -γ\ln|p-p_*|$ along any one-parameter family of initial data such that the threshold is at $p=p_*$. Hence they conjecture that in ultrarelavistic collisions almost all the kinetic energy can be converted into gravitational waves if the impact parameter is fine-tuned to the threshold. As a toy model for the binary, they consider the geodesic motion of a test particle in a Kerr black hole spacetime, where the unstable circular geodesics play the role of critical solutions, and calculate the critical exponent $γ$. Here, we incorporate radiation reaction into this model using the self-force approximation. The critical solution now evolves adiabatically along a sequence of unstable circular geodesic orbits under the effect of the self-force. We confirm that almost all the initial energy and angular momentum are radiated on the critical solution. Our calculation suggests that, even for infinite initial energy, this happens over a finite number of orbits given by $n_\infty\simeq 0.41/η$, where $η$ is the (small) mass ratio. We derive expressions for the time spent on the critical solution, number of orbits and radiated energy as functions of the initial energy and impact parameter.

preprint2012arXiv

Shock formation in stellar perturbations and tidal shock waves in binaries

We investigate whether tidal forcing can result in sound waves steepening into shocks at the surface of a star. To model the sound waves and shocks, we consider adiabatic non-spherical perturbations of a Newtonian perfect fluid star. Because tidal forcing of sounds waves is naturally treated with linear theory, but the formation of shocks is necessarily nonlinear, we consider the perturbations in two regimes. In most of the interior, where tidal forcing dominates, we treat the perturbations as linear, while in a thin layer near the surface we treat them in full nonlinearity but in the approximation of plane symmetry, fixed gravitational field and a barotropic equation of state. Using a hodograph transformation, this nonlinear regime is also described by a linear equation. We show that the two regimes can be matched to give rise to a single mode equation which is linear but models nonlinearity in the outer layers. This can then be used to obtain an estimate for the critical mode amplitude at which a shock forms near the surface. As an application, we consider the tidal waves raised by the companion in an irrotational binary system in circular orbit. We find that shocks form at the same orbital separation where Roche lobe overflow occurs, and so shock formation is unlikely to occur.

preprint2011arXiv

A conservation law formulation of nonlinear elasticity in general relativity

We present a practical framework for ideal hyperelasticity in numerical relativity. For this purpose, we recast the formalism of Carter and Quintana as a set of Eulerian conservation laws in an arbitrary 3+1 split of spacetime. The resulting equations are presented as an extension of the standard Valencia formalism for a perfect fluid, with additional terms in the stress-energy tensor, plus a set of kinematic conservation laws that evolve a configuration gradient. We prove that the equations can be made symmetric hyperbolic by suitable constraint additions, at least in a neighbourhood of the unsheared state. We discuss the Newtonian limit of our formalism and its relation to a second formalism also used in Newtonian elasticity. We validate our framework by numerically solving a set of Riemann problems in Minkowski spacetime, as well as Newtonian ones from the literature.

preprint2011arXiv

Universality in the run-up of shock waves to the surface of a star

We investigate the run-up of a shock wave from inside to the surface of a perfect fluid star in equilibrium and bounded by vacuum. Near the surface we approximate the fluid motion as plane-symmetric and the gravitational field as constant. We consider the "hot" equation of state $P=(Γ-1)ρe$ and its "cold" (fixed entropy, barotropic) form $P=K_0ρ^Γ$ (the latter does not allow for shock heating). We find numerically that the evolution of generic initial data approaches universal similarity solutions sufficiently near the surface, and we construct these similarity solutions explicitly. The two equations of state show very different behaviour, because shock heating becomes the dominant effect when it is allowed. In the barotropic case, the fluid velocity behind the shock approaches a constant value, while the density behind the shock approaches a power law in space, as the shock approaches the surface. In the hot case with shock heating,the density jumps by a constant factor through the shock, while the sound speed and fluid velocity behind the shock diverge in a whiplash effect. We tabulate the similarity exponents as a function of the equation of state parameter $Γ$and the stratification index $n_*$.

preprint2009arXiv

Generic behaviour of nonlinear sound waves near the surface of a star: smooth solutions

We are interested in the generic behaviour of nonlinear sound waves as they approach the surface of a star, here assumed to have the polytropic equation of state $P=Kρ^Γ$. Restricting to spherical symmetry, and considering only the region near the surface, we generalise the methods of Carrier and Greenspan (1958) for the shallow water equations on a sloping beach to this problem. We give a semi-quantitative criterion for a shock to form near the surface during the evolution of generic initial data with support away from the surface. We show that in smooth solutions the velocity and the square of the sound speed remain regular functions of Eulerian radius at the surface.

preprint2007arXiv

Critical phenomena in gravitational collapse

As first discovered by Choptuik, the black hole threshold in the space of initial data for general relativity shows both surprising structure and surprising simplicity. Universality, power-law scaling of the black hole mass, and scale echoing have given rise to the term "critical phenomena". They are explained by the existence of exact solutions which are attractors within the black hole threshold, that is, attractors of codimension one in phase space, and which are typically self-similar. Critical phenomena give a natural route from smooth initial data to arbitrarily large curvatures visible from infinity, and are therefore likely to be relevant for cosmic censorship, quantum gravity, astrophysics, and our general understanding of the dynamics of general relativity.

preprint2005arXiv

Constraint damping in the Z4 formulation and harmonic gauge

We show that by adding suitable lower-order terms to the Z4 formulation of the Einstein equations, all constraint violations except constant modes are damped. This makes the Z4 formulation a particularly simple example of a lambda-system as suggested by Brodbeck et al. We also show that the Einstein equations in harmonic coordinates can be obtained from the Z4 formulation by a change of variables that leaves the implied constraint evolution system unchanged. Therefore the same method can be used to damp all constraints in the Einstein equations in harmonic gauge.

preprint2002arXiv

Black-hole threshold solutions in stiff fluid collapse

Numerical studies of the gravitational collapse of a stiff (P=rho) fluid have found the now familiar critical phenomena, namely scaling of the black hole mass with a critical exponent and continuous self-similarity at the threshold of black hole formation. Using the equivalence of an irrotational stiff fluid to a massless scalar field, we construct the critical solution as a scalar field solution by making a self-similarity ansatz. We find evidence that this solution has exactly one growing perturbation mode; both the mode and the critical exponent, gamma ~ 0.94, derived from its eigenvalue agree with those measured in perfect fluid collapse simulations. We explain why this solution is seen as a critical solution in stiff fluid collapse but not in scalar field collapse, and conversely why the scalar field critical solution is not seen in stiff fluid collapse, even though the two systems are locally equivalent.

preprint2002arXiv

Perturbations of an exact solution for 2+1 dimensional critical collapse

We find the perturbation spectrum of a family of spherically symmetric and continuously self-similar (CSS) exact solutions that appear to be relevant for the critical collapse of scalar field matter in 2+1 spacetime dimensions. The rate of exponential growth of the unstable perturbation yields the critical exponent. Our results are compared to the numerical simulations of Pretorius and Choptuik and are inconclusive: We find a CSS solution with exactly one unstable mode, which suggests that it may be the critical solution, but another CSS solution which has three unstable modes fits the numerically found critical solution better.

preprint2000arXiv

Gauge-invariant and coordinate-independent perturbations of stellar collapse II: matching to the exterior

In Paper I in this series we constructed evolution equations for the complete gauge-invariant linear perturbations of a time-dependent spherically symmetric perfect fluid spacetime. A key application of this formalism is the interior of a collapsing star. Here we derive boundary conditions at the surface of the star, matching the interior perturbations to the well-known perturbations of the vacuum Schwarzschild spacetime outside the star.

preprint1993arXiv

The Isaacson expansion in quantum cosmology

This paper is an application of the ideas of the Born-Oppenheimer (or slow/fast) approximation in molecular physics and of the Isaacson (or short-wave) approximation in classical gravity to the canonical quantization of a perturbed minisuperspace model of the kind examined by Halliwell and Hawking. Its aim is the clarification of the role of the semiclassical approximation and the backreaction in such a model. Approximate solutions of the quantum model are constructed which are not semiclassical, and semiclassical solutions in which the quantum perturbations are highly excited.