David Neilsen

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.

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.

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.

Mergers of magnetized neutron stars with spinning black holes: disruption, accretion, and fallback.

We investigate the merger of a neutron star in orbit about a spinning black hole in full general relativity with a mass ratio of 5:1, allowing the star to have an initial magnetization of 10(12)  G. We present the resulting gravitational waveform and analyze the fallback accretion as the star is disrupted. We see no significant dynamical effects in the simulations or changes in the gravitational waveform resulting 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.

The transient gravitational-wave sky

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

Effects of the microphysical Equation of State in the mergers of magnetized Neutron Stars With Neutrino Cooling

We study the merger of binary neutron stars using different realistic, microphysical nuclear equations of state, as well as incorporating magnetic field and 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 cooling and electromagnetic effects, in determining the properties of the hypermassive neutron star resulting from the merger, the production of neutrinos, and the characteristics of ejecta from the system. The ejecta we find are consistent with other recent studies that find soft equations of state produce more ejecta than stiffer equations of state. Moreover, the degree of neutron richness increases for softer equations of state. In light of reported kilonova observations (associated to GRB~130603B and GRB~060614) and the discovery of relatively low abundances of heavy, radioactive elements in deep sea deposits (with respect to possible production via supernovae), we speculate that a soft EoS might be preferred---because of its significant production of sufficiently neutron rich ejecta---if such events are driven by binary neutron star mergers. We also find that realistic magnetic field strengths, obtained with a sub-grid model tuned to capture magnetic amplification via the Kelvin-Helmholtz instability at merger, are generally too weak to affect the gravitational wave signature post-merger within a time scale of

Electromagnetic and gravitational outputs from binary-neutron-star coalescence.

\approx 10

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

~ms but can have subtle effects on the post-merger dynamics.

Boosting jet power in black hole spacetimes

Carlos Palenzuela, Luis Lehner, Marcelo Ponce, Steven L. Liebling, Matthew Anderson, David Neilsen, and Patrick Motl Canadian Institute for Theoretical Astrophysics, Toronto, Ontario M5S 3H8, Canada, Perimeter Institute for Theoretical Physics,Waterloo, Ontario N2L 2Y5, Canada Department of Physics, University of Guelph, Guelph, Ontario N1G 2W1, Canada, Department of Physics, Long Island University, New York 11548, USA 5 Pervasive Technology Institute, Indiana University, Bloomington, IN 47405, USA Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA, Department of Science, Mathematics and Informatics, Indiana University Kokomo, Kokomo, IN 46904, USA,