Source author record

Matthew W. Choptuik

Matthew W. Choptuik 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

14works
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

14 published item(s)

preprint2023arXiv

Stability of Non-Minimally Coupled Topological-Defect Boson Stars

As shown by Marunovic and Murkovic, non-minimal d-stars, composite structures consisting of a boson star and a global monopole non-minimally coupled to the general relativistic field, can have extremely high gravitational compactness. In a previous paper we demonstrated that these ground-state stationary solutions are sometimes additionally characterized by shells of bosonic matter located far from the center of symmetry. In order to investigate the question of stability posed by Marunovic and Murkovic, we investigate the stability of several families of d-stars using both numerical simulations and linear perturbation theory. For all families investigated, we find that the most highly compact solutions, along with those solutions exhibiting shells of bosonic matter, are unstable to radial perturbations and are therefore poor candidates for astrophysically-relevant black hole mimickers or other highly compact stable objects.

preprint2016arXiv

Black hole formation in Randall-Sundrum II braneworlds

We present the first numerical study of the full dynamics of a braneworld scenario, working within the framework of the single brane model of Randall and Sundrum (RSII). In particular, we study the process of gravitational collapse driven by a massless scalar field which is confined to the brane. Imposing spherical symmetry on the brane, we show that the evolutions of sufficiently strong initial configurations of the scalar field result in black holes that have finite extension into the bulk. Furthermore, we find preliminary evidence that the black holes generated form a unique sequence, irrespective of the details of the initial data. The black hole solutions we obtain from dynamical evolutions are consistent with those previously computed from a static vacuum ansatz.

preprint2015arXiv

Black hole critical behaviour with the generalized BSSN formulation

The development of hyperbolic formulations of Einstein's equations has revolutionized our ability to perform long-time, stable, accurate numerical simulations of strong field gravitational phenomena. However, hyperbolic methods have seen relatively little application in one area of interest, type II critical collapse, where the challenges for a numerical code are particularly severe. Using the critical collapse of a massless scalar field in spherical symmetry as a test case, we study a generalization of the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation due to Brown that is suited for use with curvilinear coordinates. We adopt standard dynamical gauge choices, including 1+log slicing and a shift that is either zero or evolved by a Gamma-driver condition. With both choices of shift we are able to evolve sufficiently close to the black hole threshold to (1) unambiguously identify the discrete self-similarity of the critical solution (2) determine an echoing exponent consistent with previous calculations and (3) measure a mass scaling exponent, also in accord with prior computations. Our results can be viewed as an encouraging first step towards the use of hyperbolic formulations in more generic type II scenarios, including the as yet unresolved problem of critical collapse of axisymmetric gravitational waves.

preprint2015arXiv

Driven neutron star collapse: Type~I critical phenomena and the initial black hole mass distribution

We study the general relativistic collapse of neutron star (NS) models in spherical symmetry. Our initially stable models are driven to collapse by the addition of one of two things: an initially in-going velocity profile, or a shell of minimally coupled, massless scalar field that falls onto the star. Tolman-Oppenheimer-Volkoff (TOV) solutions with an initially isentropic, gamma-law equation of state serve as our NS models. The initial values of the velocity profile's amplitude and the star's central density span a parameter space which we have surveyed extensively and which we find provides a rich picture of the possible end states of NS collapse. This parameter space survey elucidates the boundary between Type I and Type II critical behavior in perfect fluids which coincides, on the subcritical side, with the boundary between dispersed and bound end states. For our particular model, initial velocity amplitudes greater than 0.3c are needed to probe the regime where arbitrarily small black holes can form. In addition, we investigate Type I behavior in our system by varying the initial amplitude of the initially imploding scalar field. In this case we find that the Type I critical solutions resemble TOV solutions on the 1-mode unstable branch of equilibrium solutions, and that the critical solutions' frequencies agree well with the fundamental mode frequencies of the unstable equilibria. Additionally, the critical solution's scaling exponent is shown to be well approximated by a linear function of the initial star's central density.

preprint2014arXiv

Critical collapse in the spherically-symmetric Einstein-Vlasov model

We solve the coupled Einstein-Vlasov system in spherical symmetry using direct numerical integration of the Vlasov equation in phase space. Focusing on the case of massless particles we study critical phenomena in the model, finding strong evidence for generic type I behaviour at the black hole threshold that parallels what has previously been observed in the massive sector. For differing families of initial data we find distinct critical solutions, so there is no universality of the critical configuration itself. However we find indications of at least a weak universality in the lifetime scaling exponent, which is yet to be understood. Additionally, we clarify the role that angular momentum plays in the critical behaviour in the massless case.

preprint2010arXiv

Generalized harmonic formulation in spherical symmetry

In this pedagogically structured article, we describe a generalized harmonic formulation of the Einstein equations in spherical symmetry which is regular at the origin. The generalized harmonic approach has attracted significant attention in numerical relativity over the past few years, especially as applied to the problem of binary inspiral and merger. A key issue when using the technique is the choice of the gauge source functions, and recent work has provided several prescriptions for gauge drivers designed to evolve these functions in a controlled way. We numerically investigate the parameter spaces of some of these drivers in the context of fully non-linear collapse of a real, massless scalar field, and determine nearly optimal parameter settings for specific situations. Surprisingly, we find that many of the drivers that perform well in 3+1 calculations that use Cartesian coordinates, are considerably less effective in spherical symmetry, where some of them are, in fact, unstable.

preprint2010arXiv

Ultra Relativistic Particle Collisions

We present results from numerical solution of the Einstein field equations describing the head-on collision of two solitons boosted to ultra relativistic energies. We show, for the first time, that at sufficiently high energies the collision leads to black hole formation, consistent with hoop conjecture arguments. This implies that the non-linear gravitational interaction between the kinetic energy of the solitons causes gravitational collapse, and that arguments for black hole formation in super-Planck scale particle collisions are robust.

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.

preprint1999arXiv

Critical phenomena in perfect fluids

We investigate the gravitational collapse of a spherically symmetric, perfect fluid with equation of state P = (Gamma -1)rho. We restrict attention to the ultrarelativistic (``kinetic-energy-dominated'', ``scale-free'') limit where black hole formation is anticipated to turn on at infinitesimal black hole mass (Type II behavior). Critical solutions (those which sit at the threshold of black hole formation in parametrized families of collapse) are found by solving the system of ODEs which result from a self-similar ansatz, and by solving the full Einstein/fluid PDEs in spherical symmetry. These latter PDE solutions (``simulations'') extend the pioneering work of Evans and Coleman (Gamma = 4/3) and verify that the continuously self-similar solutions previously found by Maison and Hara et al for $1.05 < Gamma < 1.89 are (locally) unique critical solutions. In addition, we find strong evidence that globally regular critical solutions do exist for 1.89 < Gamma <= 2$, that the sonic point for Gamma_dn = 1.8896244 is a degenerate node, and that the sonic points for Gamma >Gamma_dn are nodal points rather than focal points as previously reported. We also find a critical solution for Gamma = 2, and present evidence that it is continuously self-similar and Type II. Mass-scaling exponents for all of the critical solutions are calculated by evolving near-critical initial data, with results which confirm and extend previous calculations based on linear perturbation theory. Finally, we comment on critical solutions generated with an ideal-gas equation of state.

preprint1999arXiv

New Critical Behavior in Einstein-Yang-Mills Collapse

We extend the investigation of the gravitational collapse of a spherically symmetric Yang-Mills field in Einstein gravity and show that, within the black hole regime, a new kind of critical behavior arises which separates black holes formed via Type I collapse from black holes formed through Type II collapse. Further, we provide evidence that these new attracting critical solutions are in fact the previously discovered colored black holes with a single unstable mode.

preprint1999arXiv

Ultrarelativistic fluid dynamics

This is the first of two papers examining the critical collapse of spherically symmetric perfect fluids with the equation of state P = (Gamma -1)rho. Here we present the equations of motion and describe a computer code capable of simulating the extremely relativistic flows encountered in critical solutions for Gamma <= 2. The fluid equations are solved using a high-resolution shock-capturing scheme based on a linearized Riemann solver.

preprint1997arXiv

Instability of an "Approximate Black Hole"

We investigate the stability of a family of spherically symmetric static solutions in vacuum Brans-Dicke theory (with $ω=0$) recently described by van Putten. Using linear perturbation theory, we find one exponentially growing mode for every member of the family of solutions, and thus conclude that the solutions are not stable. Using a previously constructed code for spherically symmetric Brans-Dicke, additional evidence for instability is provided by directly evolving the static solutions with perturbations. The full non-linear evolutions also suggest that the solutions are black-hole-threshold critical solutions.

preprint1996arXiv

Black Hole--Scalar Field Interactions in Spherical Symmetry

We examine the interactions of a black hole with a massless scalar field using a coordinate system which extends ingoing Eddington-Finkelstein coordinates to dynamic spherically symmetric-spacetimes. We avoid problems with the singularity by excising the region of the black hole interior to the apparent horizon. We use a second-order finite difference scheme to solve the equations. The resulting program is stable and convergent and will run forever without problems. We are able to observe quasi-normal ringing and power-law tails as well an interesting nonlinear feature.