Steven L. Liebling

Dynamical boson stars

The idea of stable, localized bundles of energy has strong appeal as a model for particles. In the 1950s, John Wheeler envisioned such bundles as smooth configurations of electromagnetic energy that he called geons, but none were found. Instead, particle-like solutions were found in the late 1960s with the addition of a scalar field, and these were given the name boson stars. Since then, boson stars find use in a wide variety of models as sources of dark matter, as black hole mimickers, in simple models of binary systems, and as a tool in finding black holes in higher dimensions with only a single Killing vector. We discuss important varieties of boson stars, their dynamic properties, and some of their uses, concentrating on recent efforts.

Boson stars in AdS spacetime

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.

Holographic Thermalization, Stability of Anti–de Sitter Space, and the Fermi-Pasta-Ulam 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 anti-de Sitter space (AdS) but whose numerical evolution does not result in black hole formation. According to the AdS/conformal field theory (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-time-scale 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 quasiperiodic solutions. Finally, we demonstrate a striking parallel between the dynamics of AdS and the classic Fermi-Pasta-Ulam-Tsingou problem.

Magnetized Neutron-Star Mergers and Gravitational-Wave Signals

We investigate the influence of magnetic fields upon the dynamics of, and resulting gravitational waves from, a binary neutron-star merger in full general relativity coupled to ideal magnetohydrodynamics. We consider two merger scenarios: one where the stars have aligned poloidal magnetic fields and one without. Both mergers result in a strongly differentially rotating object. In comparison to the nonmagnetized scenario, the aligned magnetic fields delay the full merger of the stars. During and after merger we observe phenomena driven by the magnetic field, including Kelvin-Helmholtz instabilities in shear layers, winding of the field lines, and transition from poloidal to toroidal magnetic fields. These effects not only mediate the production of electromagnetic radiation, but also can have a strong influence on the gravitational waves. Thus, there are promising prospects for studying such systems with both types of waves.

Dual Jets from Binary Black Holes

When Black Holes Collide When galaxies merge, their central black holes interact, initially coming together in a binary system and eventually coalescing into one single black hole. The dynamics of orbiting binary black holes in vacuum is well understood; however, when black holes merge, their accretion disks are expected to combine into a circumbinary disk anchoring a magnetic field. Numerical simulations from Palenzuela et al. (p. 927; see the Perspective by Yunes) now show that when this surrounding environment is taken into account, the black holes effectively stir the plasma that surrounds them, generating collimated beams of electromagnetic radiation that then transition to a single jet as the black holes merge due to the emission of gravitational waves. Simulations show that merging black holes may produce a detectable electromagnetic signature. The coalescence of supermassive black holes—a natural outcome when galaxies merge—should produce gravitational waves and would likely be associated with energetic electromagnetic events. We have studied the coalescence of such binary black holes within an external magnetic field produced by the expected circumbinary disk surrounding them. Solving the Einstein equations to describe black holes interacting with surrounding 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 that the electromagnetic field extracts energy from the orbiting black holes, which ultimately merge and settle into the standard Blandford-Znajek scenario. Emissions along these jets could potentially be observable at large distances.

Scalar collapse in AdS spacetimes

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.

Critical collapse of the massless scalar field in axisymmetry

We present the results from a numerical study of critical gravitational collapse of axisymmetric distributions of massless scalar field energy. We find threshold behavior that can be described by the spherically symmetric critical solution with axisymmetric perturbations. However, we see indications of a growing, nonspherical mode about the spherically symmetric critical solution. The effect of this instability is that the small asymmetry present in what would otherwise be a spherically symmetric self-similar solution grows. This growth continues until a bifurcation occurs and two distinct regions form on the axis, each resembling the spherically symmetric self-similar solution. The existence of a nonspherical unstable mode is in conflict with previous perturbative results, and we therefore discuss whether such a mode exists in the continuum limit, or whether we are instead seeing a marginally stable mode that is rendered unstable by numerical approximation.

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 used this approach to demonstrate the conversion of some of the black holes 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 Maxwells 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.

Relativistic MHD with adaptive mesh refinement

This paper presents a new computer code to solve the general relativistic magnetohydrodynamics (GRMHD) equations using distributed parallel adaptive mesh refinement (AMR). The fluid equations are solved using a finite difference convex ENO method (CENO) in 3 + 1 dimensions, and the AMR is Berger–Oliger. Hyperbolic divergence cleaning is used to control the ∇ ⋅ B = 0 constraint. We present results from three flat space tests, and examine the accretion of a fluid onto a Schwarzschild black hole, reproducing the Michel solution. The AMR simulations substantially improve performance while reproducing the resolution equivalent unigrid simulation results. Finally, we discuss strong scaling results for parallel unigrid and AMR runs.

Binary black holes' effects on electromagnetic fields.

In addition to producing gravitational waves, the dynamics of a binary black hole system could induce emission of electromagnetic radiation by affecting the behavior of plasmas and electromagnetic fields in their vicinity. We here study how the electromagnetic fields are affected by a pair of orbiting black holes through the merger. In particular, we show how the binarys dynamics induce a variability in possible electromagnetically induced emissions as well as a possible enhancement of electromagnetic fields during the late-merge and merger epochs. These time dependent features will likely leave their imprint in processes generating detectable emissions and can be exploited in the detection of electromagnetic counterparts of gravitational waves.