Eric W. Hirschmann

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.

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.

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.

Simulating binary neutron stars: Dynamics and gravitational waves

We model two mergers of orbiting binary neutron stars, the first forming a black hole and the second a differentially rotating neutron star. We extract gravitational waveforms in the wave zone. Comparisons to a post-Newtonian analysis allow us to compute the orbital kinematics, including trajectories and orbital eccentricities. We verify our code by evolving single stars and extracting radial perturbative modes, which compare very well to results from perturbation theory. The Einstein equations are solved in a first-order reduction of the generalized harmonic formulation, and the fluid equations are solved using a modified convex essentially non-oscillatory method. All calculations are done in three spatial dimensions without symmetry assumptions. We use the had computational infrastructure for distributed adaptive mesh refinement.

An axisymmetric gravitational collapse code

We present a new numerical code designed to solve the Einstein field equations for axisymmetric spacetimes. The long-term goal of this project is to construct a code that will be capable of studying many problems of interest in axisymmetry, including gravitational collapse, critical phenomena, investigations of cosmic censorship and head-on black-hole collisions. Our objective here is to detail the (2+1)+1 formalism we use to arrive at the corresponding system of equations and the numerical methods we use to solve them. We are able to obtain stable evolution, despite the singular nature of the coordinate system on the axis, by enforcing appropriate regularity conditions on all variables and by adding numerical dissipation to hyperbolic equations.

Magnetic solutions to 2+1 gravity.

We report on a new solution to the Einstein-Maxwell equations in 2+1 dimensions with a negative cosmological constant. The solution is static, rotationally symmetric and has a non-zero magnetic field. The solution can be interpreted as a monopole with an everywhere finite energy density.

Perturbed disks get shocked: Binary black hole merger effects on accretion disks

The merger process of a binary black hole system can have a strong impact on a circumbinary disk. In the present work we study the effect of both central mass reduction (due to the energy loss through gravitational waves) and a possible black hole recoil (due to asymmetric emission of gravitational radiation). For the mass reduction case and recoil directed along the disks angular momentum, oscillations are induced in the disk which then modulate the internal energy and bremsstrahlung luminosities. On the other hand, when the recoil direction has a component orthogonal to the disks angular momentum, the disks dynamics are strongly impacted, giving rise to relativistic shocks. The shock heating leaves its signature in our proxies for radiation, the total internal energy and bremsstrahlung luminosity. Interestingly, for cases where the kick velocity is below the smallest orbital velocity in the disk (a likely scenario in real active galactic nuclei), we observe a common, characteristic pattern in the internal energy of the disk. Variations in kick velocity simply provide a phase offset in the characteristic pattern implying that observations of such a signature could yield a measure of the kick velocity through electromagnetic signals alone.

Collapse of a scalar field in 2 + 1 gravity

We consider the problem of critical gravitational collapse of a scalar field in 2 + 1 dimensions with spherical (circular) symmetry. After surveying all the analytic, continuously self-similar solutions and considering their global structure, we examine their perturbations with the intent of understanding which are the critical solutions with a single unstable mode. The critical solution which we find is the one which agrees most closely with that found in numerical evolutions. However, the critical exponent which we find does not seem to agree with the numerical result.

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. @S0556-2821~99!03024-6#

Criticality and bifurcation in the gravitational collapse of a self-coupled scalar field

We examine the gravitational collapse of a non-linear sigma model in spherical symmetry. There exists a family of continuously self-similar solutions parameterized by the coupling constant of the theory. These solutions are calculated together with the critical exponents for black hole formation of these collapse models. We also find that the sequence of solutions exhibits a Hopf-type bifurcation as the continuously self-similar solutions become unstable to perturbations away from self-similarity.