Source author record

Paul Duffell

Paul Duffell 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

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

10 published item(s)

preprint2016arXiv

A transition in circumbinary accretion discs at a binary mass ratio of 1:25

We study circumbinary accretion discs in the framework of the restricted three-body problem (R3Bp) and via numerically solving the height-integrated equations of viscous hydrodynamics. Varying the mass ratio of the binary, we find a pronounced change in the behaviour of the disc near mass ratio $q \equiv M_s/M_p \sim 0.04$. For mass ratios above $q=0.04$, solutions for the hydrodynamic flow transition from steady, to strongly-fluctuating; a narrow annular gap in the surface density around the secondary's orbit changes to a hollow central cavity; and a spatial symmetry is lost, resulting in a lopsided disc. This phase transition is coincident with the mass ratio above which stable orbits do not exist around the L4 and L5 equilibrium points of the R3B problem. Using the DISCO code, we find that for thin discs, for which a gap or cavity can remain open, the mass ratio of the transition is relatively insensitive to disc viscosity and pressure. The $q=0.04$ transition has relevance for the evolution of massive black hole binary+disc systems at the centers of galactic nuclei, as well as for young stellar binaries and possibly planets around brown dwarfs.

preprint2016arXiv

Synchrotron Magnetic Fields from Rayleigh-Taylor Instability in Supernovae

Synchrotron emission from a supernova necessitates a magnetic field, but it is unknown how strong the relevant magnetic fields are, and what mechanism generates them. In this study, we perform high-resolution numerical gas dynamics calculations to determine the growth of turbulence due to Rayleigh-Taylor instability, and the resulting kinetic energy in turbulent fluctuations, to infer the strength of magnetic fields amplified by this turbulence. We find that Rayleigh-Taylor instability can produce turbulent fluctuations strong enough to amplify magnetic fields to a few percent of equipartition with the thermal energy. This turbulence stays concentrated near the reverse shock, but averaging this magnetic energy throughout the shocked region (weighting by emissivity) sets the magnetic fields at a minimum of 0.3 percent of equipartition. This suggests a minimum effective magnetic field strength ($ε_B > 0.003$) which should be present in all interacting supernovae.

preprint2015arXiv

A reduced orbital period for the supermassive black hole binary candidate in the quasar PG 1302-102?

Graham et al. (2015) have detected a 5.2 year periodic optical variability of the quasar PG 1302-102 at redshift $z=0.3$, which they interpret as the redshifted orbital period $(1+z)t_{\rm bin}$ of a putative supermassive black hole binary (SMBHB). Here we consider the implications of a $3-8$ times shorter orbital period, suggested by hydrodynamical simulations of circumbinary discs (CBDs) with nearly equal--mass SMBHBs ($q\equiv M_2/M_1\gtrsim 0.3$). With the corresponding $2-4$ times tighter binary separation, PG 1302 would be undergoing gravitational wave dominated inspiral, and serve as a proof that the BHs can be fueled and produce bright emission even in this late stage of the merger. The expected fraction of binaries with the shorter $t_{\rm bin}$, among bright quasars, would be reduced by 1-2 orders of magnitude, compared to the 5.2 year period, in better agreement with the rarity of candidates reported by Graham et al. (2015). Finally, shorter periods would imply higher binary speeds, possibly imprinting periodicity on the light curves from relativistic beaming, as well as measurable relativistic effects on the Fe K$α$ line. The CBD model predicts additional periodic variability on timescales of $t_{\rm bin}$ and $\approx 0.5 t_{\rm bin}$, as well as periodic variation of broad line widths and offsets relative to the narrow lines, which are consistent with the observations. Future observations will be able to test these predictions and hence the binary+CBD hypothesis for PG 1302.

preprint2015arXiv

Producing Magnetar Magnetic Fields in the Merger of Binary Neutron Stars

The merger of binary neutron stars (BNSs) can lead to large amplifications of the magnetic field due to the development of turbulence and instabilities in the fluid, such as the Kelvin-Helmholtz shear instability, which drive small-scale dynamo activity. In order to properly resolve such instabilities and obtain the correct magnetic field amplification, one would need to employ resolutions that are currently unfeasible in global general relativistic magnetohydrodynamic (GRMHD) simulations of BNS mergers. Here, we present a subgrid model that allows global simulations to take into account the small-scale amplification of the magnetic field which is caused by the development of turbulence during BNS mergers. Assuming dynamo saturation, we show that magnetar-level fields ($\sim 10^{16}\,{\rm G}$) can be easily reached, and should therefore be expected from the merger of magnetized BNSs. The total magnetic energy can reach values up to $\sim 10^{51}\,{\rm erg}$ and the post-merger remnant can therefore emit strong electromagnetic signals and possibly produce short gamma-ray bursts.

preprint2014arXiv

Balancing the load: A Voronoi based scheme for parallel computations

The use of numerical simulations in science is ever increasing and with it the computational size. In many cases single processors are no longer adequate and simulations are run on multiple core machines or supercomputers. One of the key issues when running a simulation on multiple CPUs is maintaining a proper load balance throughout the run and minimizing communications between CPUs. We propose a novel method of utilizing a Voronoi diagram to achieve a nearly perfect load balance without the need of any global redistributions of data. As a show case, we implement our method in RICH, a 2D moving mesh hydrodynamical code, but it can be extended trivially to other codes in 2D or 3D. Our tests show that this method is indeed efficient and can be used in a large variety of existing hydrodynamical codes as well as other applications.

preprint2014arXiv

Binary Black Hole Accretion During Inspiral and Merger

We present the results of 2D, moving mesh, viscous hydrodynamical simulations of accretion onto merging supermassive black hole (SMBH) binaries. We include viscous heating, shock heating, and radiative cooling, and simulate the transition from the "pre-decoupling" epoch, where the inspiral timescale is longer than the viscous timescale, to the "post-decoupling" epoch, where the inspiral timescale is shorter than the viscous timescale. We find that there is no abrupt halt to the accretion at decoupling, but rather the accretion shows a slow decay, with significant accretion well after the expected decoupling. Moreover, we find that the luminosity in X-rays is significantly higher prior to the merger, as orbital energy from the SMBH binary is converted to heat via strong shocks inside the cavity, and radiated away. Following the merger, the cavity refills viscously and the accretion rate relaxes to the Shakura-Sunyaev value, while the X-ray luminosity drops as the shocks quickly dissipate.

preprint2014arXiv

Characteristic Signatures in the Thermal Emission from Accreting Binary Black Holes

We present the results of a calculation of the thermal spectrum from a 2D, moving mesh, high-accuracy, viscous hydrodynamical simulation of an accreting supermassive black hole binary. We include viscous heating, shock heating, and radiative cooling, evolving for longer than a viscous time so that we reach a quasi-steady accretion state. In agreement with previous work, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We also find that emission from the shock-heated mini-disks and accretion streams prevents any deficit in high-energy emission that may be expected inside the circumbinary cavity, and instead leads to a characteristic brightening of the spectrum beginning in soft X-rays.

preprint2014arXiv

Shock Corrugation by Rayleigh-Taylor Instability in GRB Afterglow Jets

Afterglow jets are Rayleigh-Taylor unstable and therefore turbulent during the early part of their deceleration. There are also several processes which actively cool the jet. In this letter, we demonstrate that if cooling significantly increases the compressibility of the flow, the turbulence collides with the forward shock, destabilizing and corrugating it. In this case, the forward shock is turbulent enough to produce the magnetic fields responsible for synchrotron emission via small scale turbulent dynamo. We calculate light curves assuming the magnetic field is in energy equipartition with the turbulent kinetic energy and discover that dynamic magnetic fields are well-approximated by a constant magnetic-to- thermal energy ratio of 1%, though there is a sizeable delay in the time of peak flux as the magnetic field turns on only after the turbulence has activated. The reverse shock is found to be significantly more magnetized than the forward shock, with a magnetic-to-thermal energy ratio of order 10%. This work motivates future Rayleigh-Taylor calculations using more physical cooling models.

preprint2013arXiv

Binary Black Hole Accretion From a Circumbinary Disk: Gas Dynamics Inside the Central Cavity

We present the results of 2D hydrodynamical simulations of circumbinary disk accretion using the finite-volume code DISCO. This code solves the 2D viscous Navier-Stokes equations on a high-resolution moving mesh which shears with the fluid flow, greatly reducing advection errors in comparison with a fixed grid. We perform a series of simulations for binary mass ratios in the range 0.026 < q < 1.0, each lasting longer than a viscous time so that we reach a quasi-steady accretion state. In each case, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We find that for q > 0.1, the binary excites eccentricity in the inner region of the circumbinary disk, creating an overdense lump which gives rise to enhanced periodicity in the accretion rate. The dependence of the periodicity on mass ratio may provide a method for observationally inferring mass ratios from measurements of the accretion rate. We also find that for all mass ratios studied, the magnitude of the accretion onto the secondary is sufficient to drive the binary toward larger mass ratio. This suggests a mechanism for biasing mass ratio distributions toward equal mass.

preprint2013arXiv

The fate of fallback matter around newly born compact objects

The presence of fallback disks around young neutron stars has been invoked over the years to explain a large variety of phenomena. Here we perform a numerical investigation of the formation of such disks during a supernova explosion, considering both neutron star (NS) and black hole (BH) remnants. Using the public code MESA, we compute the angular momentum distribution of the pre-supernova material, for stars with initial masses M in the range 13 - 40 Msun, initial surface rotational velocities vsurf between 25% and 75% of the critical velocity, and for metallicities Z of 1%, 10% and 100% of the solar value. These pre SN models are exploded with energies E varying between 10^{50} - 3x10^{52} ergs, and the amount of fallback material is computed. We find that, if magnetic torques play an important role in angular momentum transport, then fallback disks around NSs, even for low-metallicity main sequence stars, are not an outcome of SN explosions. Formation of such disks around young NSs can only happen under the condition of negligible magnetic torques and a fine-tuned explosion energy. For those stars which leave behind BH remnants, disk formation is ubiquitous if magnetic fields do not play a strong role; however, unlike the NS case, even with strong magnetic coupling in the interior, a disk can form in a large region of the {Z,M,vsurf,E} parameter space. Together with the compact, hyperaccreting fallback disks widely discussed in the literature, we identify regions in the above parameter space which lead to extended, long-lived disks around BHs. We find that the physical conditions in these disks may be conducive to planet formation, hence leading to the possible existence of planets orbiting black holes.