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.Comment: 7 pages, 5 figure

The Surface Brightness Fluctuations and Globular Cluster Populations of M87 and its Companions

Using the surface brightness fluctuations in HST WFPC-2 images, we determine that M87, NGC 4486B, and NGC 4478 are all at a distance of ~16 Mpc, while NGC 4476 lies in the background at ~21 Mpc. We also examine the globular clusters of M87 using archived HST fields. We detect the bimodal color distribution, and find that the amplitude of the red peak relative to the blue peak is greatest near the center. This feature is in good agreement with the merger model of elliptical galaxy formation, where some of the clusters originated in progenitor galaxies while other formed during mergers.Comment: 5 pages, 2 figure

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 $\nabla\cdot {\bf 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.Comment: 24 pages, 14 figures, 3 table

Critical Phenomena in Neutron Stars I: Linearly Unstable Nonrotating Models

We consider the evolution in full general relativity of a family of linearly unstable isolated spherical neutron stars under the effects of very small, perturbations as induced by the truncation error. Using a simple ideal-fluid equation of state we find that this system exhibits a type-I critical behaviour, thus confirming the conclusions reached by Liebling et al. [1] for rotating magnetized stars. Exploiting the relative simplicity of our system, we are able carry out a more in-depth study providing solid evidences of the criticality of this phenomenon and also to give a simple interpretation of the putative critical solution as a spherical solution with the unstable mode being the fundamental F-mode. Hence for any choice of the polytropic constant, the critical solution will distinguish the set of subcritical models migrating to the stable branch of the models of equilibrium from the set of subcritical models collapsing to a black hole. Finally, we study how the dynamics changes when the numerically perturbation is replaced by a finite-size, resolution independent velocity perturbation and show that in such cases a nearly-critical solution can be changed into either a sub or supercritical. The work reported here also lays the basis for the analysis carried in a companion paper, where the critical behaviour in the the head-on collision of two neutron stars is instead considered [2].Comment: 15 pages, 9 figure

A joint Chandra and Swift view of the 2015 X-Ray dust-scattering echo of V404 Cygni

We present a combined analysis of the Chandra and Swift observations of the 2015 X-ray echo of V404 Cygni. Using a stacking analysis, we identify eight separate rings in the echo. We reconstruct the soft X-ray light curve of the 2015 June outburst using the high-resolution Chandra images and cross-correlations of the radial intensity profiles, indicating that about 70% of the outburst fluence occurred during the bright flare at the end of the outburst on MJD 57199.8. By deconvolving the intensity profiles with the reconstructed outburst light curve, we show that the rings correspond to eight separate dust concentrations with precise distance determinations. We further show that the column density of the clouds varies significantly across the field of view, with the centroid of most of the clouds shifted toward the Galactic plane, relative to the position of V404 Cyg, invalidating the assumption of uniform cloud column typically made in attempts to constrain dust properties from light echoes. We present a new XSPEC spectral dust-scattering model that calculates the differential dust-scattering cross section for a range of commonly used dust distributions and compositions and use it to jointly fit the entire set of Swift echo data. We find that a standard Mathis-Rumpl-Nordsieck model provides an adequate fit to the ensemble of echo data. The fit is improved by allowing steeper dust distributions, and models with simple silicate and graphite grains are preferred over models with more complex composition. © 2016. The American Astronomical Society. All rights reserved

Coherent Ro-vibrational Revivals in a Thermal Molecular Ensemble

We report an experimental and theoretical study of the evolution of vibrational coherence in a thermal ensemble of nitrogen molecules. Rotational dephasing and rephasing of the vibrational coherence is detected by coherent anti-Stokes Raman scattering. The existence of ro-vibrational coupling and the discrete energy spectrum of the rotational bath lead to a whole new class of full and fractional ro-vibrational revivals. Following the rich ro-vibrational dynamics on a nanosecond time scale with sub-picosecond time resolution enables us to determine the second-order ro-vibrational constant $gamma_e$ and assess new possibilities of controlling decoherence.Comment: submitted at Physical Review

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 disk's 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 disk's angular momentum, the disk's 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 AGN), 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.Comment: 10 pages, 13 figures. v2: Minor changes, version to be published in PR

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.Comment: 14 pages, 16 figures. Added one figure from previous version; corrected typo

WhiskyMHD: a new numerical code for general relativistic magnetohydrodynamics

The accurate modelling of astrophysical scenarios involving compact objects and magnetic fields, such as the collapse of rotating magnetized stars to black holes or the phenomenology of gamma-ray bursts, requires the solution of the Einstein equations together with those of general-relativistic magnetohydrodynamics. We present a new numerical code developed to solve the full set of general-relativistic magnetohydrodynamics equations in a dynamical and arbitrary spacetime with high-resolution shock-capturing techniques on domains with adaptive mesh refinements. After a discussion of the equations solved and of the techniques employed, we present a series of testbeds carried out to validate the code and assess its accuracy. Such tests range from the solution of relativistic Riemann problems in flat spacetime, over to the stationary accretion onto a Schwarzschild black hole and up to the evolution of oscillating magnetized stars in equilibrium and constructed as consistent solutions of the coupled Einstein-Maxwell equations.Comment: minor changes to match the published versio

Bondian frames to couple matter with radiation

A study is presented for the non linear evolution of a self gravitating distribution of matter coupled to a massless scalar field. The characteristic formulation for numerical relativity is used to follow the evolution by a sequence of light cones open to the future. Bondian frames are used to endow physical meaning to the matter variables and to the massless scalar field. Asymptotic approaches to the origin and to infinity are achieved; at the boundary surface interior and exterior solutions are matched guaranteeing the Darmois--Lichnerowicz conditions. To show how the scheme works some numerical models are discussed. We exemplify evolving scalar waves on the following fixed backgrounds: A) an atmosphere between the boundary surface of an incompressible mixtured fluid and infinity; B) a polytropic distribution matched to a Schwarzschild exterior; C) a Schwarzschild- Schwarzschild spacetime. The conservation of energy, the Newman--Penrose constant preservation and other expected features are observed.Comment: 20 pages, 6 figures; to appear in General Relativity and Gravitatio