Source author record

Steven L. Liebling

Steven L. Liebling 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

31works
10topics
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

31 published item(s)

preprint2022arXiv

Conformally curved initial data for charged, spinning black hole binaries on arbitrary orbits

We present a method to construct conformally curved initial data for charged black hole binaries with spin on arbitrary orbits. We generalize the superposed Kerr-Schild, extended conformal thin sandwich construction from [Lovelace et al., Phys. Rev. D {78}, 084017 (2008)] to use Kerr-Newman metrics for the superposed black holes and to solve the electromagnetic constraint equations. We implement the construction in the pseudospectral code SGRID. The code thus provides a complementary and completely independent excision-based construction, compared to the existing charged black hole initial data constructed using the puncture method [Bozzola and Paschalidis, Phys. Rev. D {99}, 104044 (2019)]. It also provides an independent implementation (with some small changes) of the Lovelace et al. vacuum construction. We construct initial data for different configurations of orbiting binaries, e.g., with black holes that are highly charged or rapidly spinning (90 and 80 percent of the extremal values, respectively, for this initial test, though the code should be able to produce data with even higher values of these parameters using higher resolutions), as well as for generic spinning, charged black holes. We carry out exploratory evolutions with the finite difference, moving punctures codes BAM (in the vacuum case) and HAD (for head-on collisions including charge), filling inside the excision surfaces. In the charged case, evolutions of these initial data provide a proxy for binary black hole waveforms in modified theories of gravity. Moreover, the generalization of the construction to Einstein-Maxwell-dilaton theory should be straightforward.

preprint2022arXiv

Gamma-radiation sky maps from compact binaries

We study sky maps and light curves of gamma-ray emission from neutron stars in compact binaries, and in isolation. We briefly review some gamma-ray emission models, and reproduce sky maps from a standard isolated pulsar in the Separatrix Layer model. We consider isolated pulsars with several variations of a dipole magnetic field, including superpositions, and predict their gamma-ray emission. Our results provide new heuristics on what can and cannot be inferred about the magnetic field configuration of pulsars from high-energy observations. We find that typical double-peak light curves can be produced by pulsars with significant multipole structure beyond a single dipole. For binary systems, we also present a simple approximation that is useful for rapid explorations of binary magnetic field structure. Finally, we predict the gamma-ray emission pattern from a compact black hole-neutron star binary moments before merger by applying the Separatrix Layer model to data simulated in full general relativity; we find that face-on observers receive little emission, equatorial observers see one broad peak, and more generic observers typically see two peaks.

preprint2020arXiv

Accessing the axion via compact object binaries

Black holes in binaries with other compact objects can provide natural venues for indirect detection of axions or other ultralight fields. The superradiant instability associated with a rapidly spinning black hole leads to the creation of an axion cloud which carries energy and angular momentum from the black hole. This cloud will then decay via gravitational wave emission. We show that the energy lost as a result of this process tends toward an outspiraling of the binary orbit. A given binary system is sensitive to a narrow range of axion masses, determined by the mass of the black hole. Pulsar-black hole binaries, once detected in the electromagnetic band, will allow high-precision measurements of black hole mass loss via timing measurements of the companion pulsar. This avenue of investigation is particularly promising in light of the recent preliminary announcements of two candidate black hole-neutron star mergers by LIGO/VIRGO (#S190814bv and #S190426c). We demonstrate that for such a binary system with a typical millisecond pulsar and a 3-solar-mass black hole, axions with masses between $2.7 \times 10^{-12}$ eV and $3.2 \times 10^{-12}$ eV are detectable. Recent gravitational wave observations by LIGO/VIRGO of binary black hole mergers imply that, for these binaries, gravitational radiation from the rotating quadrupole moment is a dominant effect, causing an inspiraling orbit. With some reasonable assumptions about the period of the binary when it formed and the spins of the black holes, these observations rule out possible axion masses between $3 \times 10^{-13}$ eV and $6 \times 10^{-13}$ eV. Future binary black hole observations, for example by LISA, are expected to provide more robust bounds. In some circumstances, neutron stars may also undergo superradiant instabilities, and binary pulsars could be used to exclude axions with certain masses and matter couplings.

preprint2018arXiv

Probing Near-Horizon Fluctuations with Black Hole Binary Mergers

The strong version of the nonviolent nonlocality proposal of Giddings predicts "strong but soft" quantum metric fluctuations near black hole horizons in an attempt to resolve the information paradox. To study observable signatures of this proposal, we numerically solve Einstein's equations modified by these fluctuations and analyze the gravitational wave signal from the inspiral and merger of two black holes. In a model of evolution for such fluctuations, we show that they lead to significant deviations in the observed waveform, even when the black holes are still well separated, and could potentially be observed by aLIGO.

preprint2016arXiv

Electromagnetic Luminosity of the Coalescence of Charged Black Hole Binaries

The observation of a possible electromagnetic counterpart by the Fermi GBM group to the aLIGO detection of the merger of a black hole binary has spawned a number of ideas about its source. Furthermore, observations of fast radio bursts (FRBs) have similarly resulted in a range of new models that might endow black hole binaries with electromagnetic signatures. In this context, even the unlikely idea that astrophysical black holes may have significant charge is worth exploring, and here we present results from the simulation of weakly charged black holes as they orbit and merge. Our simulations suggest that a black hole binary with mass comparable to that observed in GW150914 could produce the level of electromagnetic luminosity observed by Fermi GBM ($10^{49}$ ergs/s) with a non-dimensional charge of $q \equiv Q/M = 10^{-4}$ assuming good radiative efficiency. However even a charge such as this is difficult to imagine avoiding neutralization long enough for the binary to produce its electromagnetic counterpart, and so this value would likely serve simply as an upper bound. On the other hand, one can equivalently consider the black holes as having acquired a magnetic monopole charge that would be easy to maintain and would generate an identical electromagnetic signature as the electric charges. The observation of such a binary would have significant cosmological implications, not the least of which would be an explanation for the quantization of charge itself. We also study such a magnetically charged binary in the force-free regime and find it much more radiative, reducing even further the requirements to produce the counterpart.

preprint2016arXiv

The $m=1$ instability \& gravitational wave signal in binary neutron star mergers

We examine the development and detectability of the $m=1$ instability in the remnant of binary neutron star mergers. The detection of the gravitational mode associated with the $m=1$ degree of freedom could potentially reveal details of the equation of state. We analyze the post-merger epoch of simulations of both equal and non-equal mass neutron star mergers using three realistic, microphysical equations of state and neutrino cooling. Our studies show such an instability develops generically and within a short dynamical time to strengths that are comparable or stronger than the $m=2$ mode which is the strongest during the early post-merger stage. We estimate the signal to noise ratio that might be obtained for the $m=1$ mode and discuss the prospects for observing this signal with available Earth-based detectors. Because the $m=1$ occurs at roughly half the frequency of the more powerful $m=2$ signal and because it can potentially be long-lived, targeted searches could be devised to observe it. We estimate that with constant amplitude direct detection of the mode could occur up to a distance of roughly $14\,\mathrm{Mpc}$ whereas a search triggered by the inspiral signal could extend this distance to roughly $100\,\mathrm{Mpc}$.

preprint2016arXiv

Unequal mass binary neutron star mergers and multimessenger signals

We study the merger of binary neutron stars with different mass ratios adopting three different realistic, microphysical nuclear equations of state, as well as incorporating neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with neutrino cooling, in determining the properties of the resulting hypermassive neutron star, of the neutrinos produced, and of the ejected material. The ejecta we find are consistent with other recent studies that find that small mass ratios produce more ejecta than equal mass cases (up to some limit) and this ejecta is more neutron rich. This trend indicates the importance with future kilonovae observations of measuring the individual masses of an associated binary neutron star system, presumably from concurrent gravitational wave observations, in order to be able to extract information about the nuclear equation of state

preprint2015arXiv

Conserved quantities and dual turbulent cascades in Anti-de Sitter spacetime

We consider the dynamics of a spherically symmetric massless scalar field coupled to general relativity in Anti--de Sitter spacetime in the small-amplitude limit. Within the context of our previously developed two time framework (TTF) to study the leading self-gravitating effects, we demonstrate the existence of two new conserved quantities in addition to the known total energy $E$ of the modes: The particle number $N$ and Hamiltonian $H$ of our TTF system. Simultaneous conservation of $E$ and $N$ implies that weak turbulent processes undergo dual cascades (direct cascade of $E$ and inverse cascade of $N$ or vice versa). This partially explains the observed dynamics of 2-mode initial data. In addition, conservation of $E$ and $N$ limits the region of phase space that can be explored within the TTF approximation and in particular rules out equipartion of energy among the modes for general initial data. Finally, we discuss possible effects of conservation of $N$ and $E$ on late time dynamics.

preprint2015arXiv

Islands of stability and recurrence times in AdS

We study the stability of anti-de Sitter (AdS) spacetime to spherically symmetric perturbations of a real scalar field in general relativity. Further, we work within the context of the "two time framework" (TTF) approximation, which describes the leading nonlinear effects for small amplitude perturbations, and is therefore suitable for studying the weakly turbulent instability of AdS---including both collapsing and non-collapsing solutions. We have previously identified a class of quasi-periodic (QP) solutions to the TTF equations, and in this work we analyze their stability. We show that there exist several families of QP solutions that are stable to linear order, and we argue that these solutions represent islands of stability in TTF. We extract the eigenmodes of small oscillations about QP solutions, and we use them to predict approximate recurrence times for generic non-collapsing initial data in the full (non-TTF) system. Alternatively, when sufficient energy is driven to high-frequency modes, as occurs for initial data far from a QP solution, the TTF description breaks down as an approximation to the full system. Depending on the higher order dynamics of the full system, this often signals an imminent collapse to a black hole.

preprint2015arXiv

Universality of non-equilibrium dynamics of CFTs from holography

Motivated by a low-energy effective description of gauge theory/string theory duality, we conjecture that the dynamics of $SO(4)$-invariant states in a large class of four-dimensional conformal gauge theories on $S^3$ with non-equal central charges $c\ne a$ are universal on time scales $t_{\rm universal}\propto ({\cal E}-{\cal E}_{\rm vacuum})^{-1}$, in the limit where the energy ${\cal E}\to {\cal E}_{\rm vacuum}$. We show that low-energy excitations in $c\ne a$ CFTs do not thermalize in this limit. The holographic universality conjecture then implies that within the Einstein-scalar field system (dual to theories with $c=a$), $AdS_5$ is stable to spherically symmetric perturbations against formation of trapped surfaces within time scales $t_{\rm universal}$.

preprint2014arXiv

Holographic Thermalization, stability of AdS, and the Fermi-Pasta-Ulam-Tsingou paradox

For a real massless scalar field in general relativity with a negative cosmological constant, we uncover a large class of spherically symmetric initial conditions that are close to AdS, but whose numerical evolution does not result in black hole formation. According to the AdS/CFT dictionary, these bulk solutions are dual to states of a strongly interacting boundary CFT that fail to thermalize at late times. Furthermore, as these states are not stationary, they define dynamical CFT configurations that do not equilibrate. We develop a two-timescale perturbative formalism that captures both direct and inverse cascades of energy and agrees with our fully nonlinear evolutions in the appropriate regime. We also show that this formalism admits a large class of quasi-periodic solutions. Finally, we demonstrate a striking parallel between the dynamics of AdS and the classic Fermi-Pasta-Ulam-Tsingou problem.

preprint2014arXiv

Interaction of misaligned magnetospheres in the coalescence of binary neutron stars

We study the dependence of the electromagnetic luminosity --produced by interactions of force-free magnetospheres-- on dipole inclinations in binary neutron star systems. We show that this interaction extracts kinetic energy from the system and powers a Poynting flux with a strong dependence on the dipole orientations. This dependence can be linked to the reconnection and redistribution of magnetic field as the stars interact. Although the details of the Poynting luminosity are very much dependent on the orientation, all the cases considered here nevertheless radiate a large Poynting flux. This robust emission suggests that the pre-merger stage of binary neutron star systems can yield interesting electromagnetic counterparts to gravitational wave events.

preprint2014arXiv

Magnetized Neutron Stars With Realistic Equations of State and Neutrino Cooling

We incorporate realistic, tabulated equations of state into fully relativistic simulations of magnetized neutron stars along with a neutrino leakage scheme which accounts for cooling via neutrino emission. Both these improvements utilize open-source code (GR1D) and tables from http://stellarcollapse.org. Our implementation makes use of a novel method for the calculation of the optical depth which simplifies its use with distributed adaptive mesh refinement. We present various tests with and without magnetization and preliminary results both from single stars and from the merger of a binary system.

preprint2014arXiv

Simulating the universe(s): from cosmic bubble collisions to cosmological observables with numerical relativity

The theory of eternal inflation in an inflaton potential with multiple vacua predicts that our universe is one of many bubble universes nucleating and growing inside an ever-expanding false vacuum. The collision of our bubble with another could provide an important observational signature to test this scenario. We develop and implement an algorithm for accurately computing the cosmological observables arising from bubble collisions directly from the Lagrangian of a single scalar field. We first simulate the collision spacetime by solving Einstein's equations, starting from nucleation and ending at reheating. Taking advantage of the collision's hyperbolic symmetry, simulations are performed with a 1+1-dimensional fully relativistic code that uses adaptive mesh refinement. We then calculate the comoving curvature perturbation in an open Friedmann-Robertson-Walker universe, which is used to determine the temperature anisotropies of the cosmic microwave background radiation. For a fiducial Lagrangian, the anisotropies are well described by a power law in the cosine of the angular distance from the center of the collision signature. For a given form of the Lagrangian, the resulting observational predictions are inherently statistical due to stochastic elements of the bubble nucleation process. Further uncertainties arise due to our imperfect knowledge about inflationary and pre-recombination physics. We characterize observational predictions by computing the probability distributions over four phenomenological parameters which capture these intrinsic and model uncertainties. This represents the first fully-relativistic set of predictions from an ensemble of scalar field models giving rise to eternal inflation, yielding significant differences from previous non-relativistic approximations. Thus, our results provide a basis for a rigorous confrontation of these theories with cosmological data.

preprint2013arXiv

Boson Stars in AdS

We construct boson stars in global Anti de Sitter (AdS) space and study their stability. Linear perturbation results suggest that the ground state along with the first three excited state boson stars are stable. We evolve some of these solutions and study their nonlinear stability in light of recent work \cite{Bizon:2011gg} arguing that a weakly turbulent instability drives scalar perturbations of AdS to black hole formation. However evolutions suggest that boson stars are nonlinearly stable and immune to the instability for sufficiently small perturbation. Furthermore, these studies find other families of initial data which similarly avoid the instability for sufficiently weak parameters. Heuristically, we argue that initial data families with widely distributed mass-energy distort the spacetime sufficiently to oppose the coherent amplification favored by the instability. From the dual CFT perspective our findings suggest that there exist families of rather generic initial conditions in strongly coupled CFT (with large number of degrees of freedom) that do not thermalize in the infinite future.

preprint2013arXiv

Gravitational and electromagnetic outputs from binary neutron star mergers

The late stage of an inspiraling neutron star binary gives rise to strong gravitational wave emission due to its highly dynamic, strong gravity. Moreover, interactions between the stellar magnetospheres can produce considerable electromagnetic radiation. We study this scenario using fully general relativistic, resistive magneto-hydrodynamics simulations. We show that these interactions extract kinetic energy from the system, dissipate heat, and power radiative Poynting flux, as well as develop current sheets. Our results indicate that this power can: (i) outshine pulsars in binaries, (ii) display a distinctive angular- and time-dependent pattern, and (iii) radiate within large opening angles. These properties suggest that some binary neutron star mergers are ideal candidates for multimessenger astronomy.

preprint2013arXiv

Linking electromagnetic and gravitational radiation in coalescing binary neutron stars

We expand on our study of the gravitational and electromagnetic emissions from the late stage of an inspiraling neutron star binary as presented in Ref. \cite{Palenzuela:2013hu}. Interactions between the stellar magnetospheres, driven by the extreme dynamics of the merger, can yield considerable outflows. We study the gravitational and electromagnetic waves produced during the inspiral and merger of a binary neutron star system using a full relativistic, resistive MHD evolution code. We show that the interaction between the stellar magnetospheres extracts kinetic energy from the system and powers radiative Poynting flux and heat dissipation. These features depend strongly on the configuration of the initial stellar magnetic moments. Our results indicate that this power can strongly outshine pulsars in binaries and have a distinctive angular and time-dependent pattern. Our discussion provides more detail than Ref. \cite{Palenzuela:2013hu}, showing clear evidence of the different effects taking place during the inspiral. Our simulations include a few milliseconds after the actual merger and study the dynamics of the magnetic fields during the formation of the hypermassive neutron star. We also briefly discuss the possibility of observing such emissions.

preprint2013arXiv

Nonlinear Collapse in the Semilinear Wave Equation in AdS

Previous studies of the semilinear wave equation in Minkowski space have shown a type of critical behavior in which large initial data collapse to singularity formation due to nonlinearities while small initial data does not. Numerical solutions in spherically symmetric Anti-de Sitter (AdS) are presented here which suggest that, in contrast, even small initial data collapse eventually. Such behavior appears analogous to the recent result of Ref. [1] that found that even weak, scalar initial data collapse gravitationally to black hole formation via a weakly turbulent instability. Furthermore, the imposition of a reflecting boundary condition in the bulk introduces a cut-off, below which initial data fails to collapse. This threshold appears to arise because of the dispersion introduced by the boundary condition.

preprint2012arXiv

Intense Electromagnetic Outbursts from Collapsing Hypermassive Neutron Stars

We study the gravitational collapse of a magnetized neutron star using a novel numerical approach able to capture both the dynamics of the star and the behavior of the surrounding plasma. In this approach, a fully general relativistic magnetohydrodynamics implementation models the collapse of the star and provides appropriate boundary conditions to a force-free model which describes the stellar exterior. We validate this strategy by comparing with known results for the rotating monopole and aligned rotator solutions and then apply it to study both rotating and non-rotating stellar collapse scenarios, and contrast the behavior with what is obtained when employing the electrovacuum approximation outside the star. The non-rotating electrovacuum collapse is shown to agree qualitatively with a Newtonian model of the electromagnetic field outside a collapsing star. We illustrate and discuss a fundamental difference between the force-free and electrovacuum solutions, involving the appearance of large zones of electric-dominated field in the vacuum case. This provides a clear demonstration of how dissipative singularities appear generically in the non-linear time-evolution of force-free fluids. In both the rotating and non-rotating cases, our simulations indicate that the collapse induces a strong electromagnetic transient. In the case of sub-millisecond rotation, the magnetic field experiences strong winding and the transient carries much more energy. This result has important implications for models of gamma-ray bursts.

preprint2012arXiv

Scalar Collapse in AdS

Recently, studies of the gravitational collapse of a scalar field within spherically symmetric AdS spacetimes was presented in \cite{Bizon:2011gg,Jalmuzna:2011qw} which showed an instability of pure AdS to black hole formation. In particular, the work showed that arbitrarily small initial configurations of scalar field evolved through some number of reflections off the AdS boundary until a black hole forms. We consider this same system, extended to include a complex scalar field, and reproduce this phenomena. We present tests of our numerical code that demonstrate convergence and consistency. We study the properties of the evolution as the scalar pulse becomes more compact examining the asymptotic behavior of the scalar field, an observable in the corresponding boundary CFT. We demonstrate that such BH formation occurs even when one places a reflecting boundary at finite radius indicating that the sharpening is a property of gravity in a bounded domain, not of AdS itself. We examine how the initial energy is transferred to higher frequencies --which leads to black hole formation-- and uncover interesting features of this transfer.

preprint2011arXiv

Boosting jet power in black hole spacetimes

The extraction of rotational energy from a spinning black hole via the Blandford-Znajek mechanism has long been understood as an important component in models to explain energetic jets from compact astrophysical sources. Here we show more generally that the kinetic energy of the black hole, both rotational and translational, can be tapped, thereby producing even more luminous jets powered by the interaction of the black hole with its surrounding plasma. We study the resulting Poynting jet that arises from single boosted black holes and binary black hole systems. In the latter case, we find that increasing the orbital angular momenta of the system and/or the spins of the individual black holes results in an enhanced Poynting flux.

preprint2010arXiv

Dual Jets from Binary Black Holes

Supermassive black holes are found at the centers of most galaxies and their inspiral is a natural outcome when galaxies merge. The inspiral of these systems is of utmost astrophysical importance as prodigious producers of gravitational waves and in their possible role in energetic electromagnetic events. We study such binary black hole coalescence under the influence of an external magnetic field produced by the expected circumbinary disk surrounding them. Solving the Einstein equations to describe the spacetime and using the force-free approach for the electromagnetic fields and the tenuous plasma, we present numerical evidence for possible jets driven by these systems. Extending the process described by Blandford and Znajek for a single spinning black hole, the picture that emerges suggests the electromagnetic field extracts energy from the orbiting black holes, which ultimately merge and settle into the standard Blandford-Znajek scenario. Emissions along dual and single jets would be expected that could be observable to large distances.

preprint2010arXiv

Dynamics of Rotating, Magnetized Neutron Stars

Using a fully general relativistic implementation of ideal magnetohydrodynamics with no assumed symmetries in three spatial dimensions, the dynamics of magnetized, rigidly rotating neutron stars are studied. Beginning with fully consistent initial data constructed with Magstar, part of the Lorene project, we study the dynamics and stability of rotating, magnetized polytropic stars as models of neutron stars. Evolutions suggest that some of these rotating, magnetized stars may be minimally unstable occurring at the threshold of black hole formation.

preprint2010arXiv

Evolutions of Magnetized and Rotating Neutron Stars

We study the evolution of magnetized and rigidly rotating neutron stars within a fully general relativistic implementation of ideal magnetohydrodynamics with no assumed symmetries in three spatial dimensions. The stars are modeled as rotating, magnetized polytropic stars and we examine diverse scenarios to study their dynamics and stability properties. In particular we concentrate on the stability of the stars and possible critical behavior. In addition to their intrinsic physical significance, we use these evolutions as further tests of our implementation which incorporates new developments to handle magnetized systems.

preprint2010arXiv

Magnetospheres of Black Hole Systems in Force-Free Plasma

The interaction of black holes with ambient magnetic fields is important for a variety of highly energetic astrophysical phenomena. We study this interaction within the force-free approximation in which a tenuous plasma is assumed to have zero inertia. Blandford and Znajek (BZ) used this approach to demonstrate the conversion of some of the black hole's energy into electromagnetic Poynting flux in stationary and axisymmetric single black hole systems. We adopt this approach and extend it to examine asymmetric and, most importantly, dynamical systems by implementing the fully nonlinear field equations of general relativity coupled to Maxwell's equations. For single black holes, we study in particular the dependence of the Poynting flux and show that, even for misalignments between the black hole spin and the direction of the asymptotic magnetic field, a Poynting flux is generated with a luminosity dependent on such misalignment. For binary black hole systems, we show both in the head-on and orbiting cases that the moving black holes generate a Poynting flux.

preprint2010arXiv

Mergers of Magnetized Neutron Stars with Spinning Black Holes: Disruption, Accretion and Fallback

We investigate the merger of a neutron star (of compaction ratio $0.1$) in orbit about a spinning black hole in full general relativity with a mass ratio of $5:1$, allowing for the star to have an initial magnetization of $10^{12} {\rm Gauss}$. We present the resulting gravitational waveform and analyze the fallback accretion as the star is disrupted. The evolutions suggest no significant effects from the initial magnetization. We find that only a negligible amount of matter becomes unbound; $99%$ of the neutron star material has a fallback time of 10 seconds or shorter to reach the region of the central engine and that $99.99%$ of the star will interact with the central disk and black hole within 3 hours.

preprint2004arXiv

The Nonlinear Sigma Model With Distributed Adaptive Mesh Refinement

An adaptive mesh refinement (AMR) scheme is implemented in a distributed environment using Message Passing Interface (MPI) to find solutions to the nonlinear sigma model. Previous work studied behavior similar to black hole critical phenomena at the threshold for singularity formation in this flat space model. This work is a follow-up describing extensions to distribute the grid hierarchy and presenting tests showing the correctness of the model.

preprint2000arXiv

Critical Phenomena in Nonlinear Sigma Models

We consider solutions to the nonlinear sigma model (wave maps) with target space S^3 and base space 3+1 Minkowski space, and we find critical behavior separating singular solutions from nonsingular solutions. For families of solutions with localized spatial support a self-similar solution is found at the boundary. For other families, we find that a static solution appears to sit at the boundary. This behavior is compared to the black hole critical phenomena found by Choptuik.

preprint1998arXiv

Multiply Unstable Black Hole Critical Solutions

The gravitational collapse of a complex scalar field in the harmonic map is modeled in spherical symmetry. Previous work has shown that a change of stability of the attracting critical solution occurs in parameter space from the discretely self-similarity critical (DSS) solution originally found by Choptuik to the continuously self-similar (CSS) solution found by Hirschmann and Eardley. In the region of parameter space in which the DSS is the attractor, a family of initial data is found which finds the CSS as its critical solution despite the fact that it has more than one unstable mode. An explanation of this is proposed in analogy to families that find the DSS in the region where the CSS is the attractor.