Showing 1–39 of 39 results for author: Duffell, P C

Sculpting the Morphology of Supernova Remnant Pa 30 via Efficient Ejecta Cooling

Authors: Paul C. Duffell, Abigail Polin, Soham Mandal

Abstract: We demonstrate in a proof-of-concept numerical hydrodynamics calculation that the narrow radial filamentary structures seen in Pa 30 could be generated through highly efficient cooling (e.g. via line emission) in the ejecta. Efficient cooling in the ejecta causes a drop of pressure support in Rayleigh-Taylor fingers, leading them to be compressed, and suppressing the growth of Kelvin-Helmholtz ins… ▽ More We demonstrate in a proof-of-concept numerical hydrodynamics calculation that the narrow radial filamentary structures seen in Pa 30 could be generated through highly efficient cooling (e.g. via line emission) in the ejecta. Efficient cooling in the ejecta causes a drop of pressure support in Rayleigh-Taylor fingers, leading them to be compressed, and suppressing the growth of Kelvin-Helmholtz instability. Following this result, we make three predictions that could determine whether this is the mechanism responsible for shaping Pa 30: First, we predict very strong emission lines, strong enough to cool a significant fraction of the shock energy in an expansion time. Secondly, we predict that the forward shock should be highly corrugated on small scales, with the shock front closely following the structure of the filaments. Third, we predict that these filaments should be nearly ballistic, with velocities around 90% of the free-expansion velocity ($v \approx 0.9 ~r/t$). These predictions should be falsifiable in follow-up observations of this remnant. △ Less

Submitted 13 May, 2024; v1 submitted 20 March, 2024; originally announced March 2024.

Comments: Rejected from both ApJ Letters and MNRAS letters. I'm giving up on this one; I'll just leave the updated version here while we work on the follow-up. Please let me know if you have any comments

Measurement of anisotropies in Supernova Remnant observations and their interpretation using numerical models

Authors: Soham Mandal, Paul C. Duffell, Abigail Polin, Dan Milisavljevic

Abstract: Supernova remnants (SNRs) exhibit varying degrees of anisotropy, which have been extensively modeled using numerical methods. We implement a technique to measure anisotropies in SNRs by calculating power spectra from their high-resolution images. To test this technique, we develop 3D hydrodynamical models of supernova remnants and generate synthetic x-ray images from them. Power spectra extracted… ▽ More Supernova remnants (SNRs) exhibit varying degrees of anisotropy, which have been extensively modeled using numerical methods. We implement a technique to measure anisotropies in SNRs by calculating power spectra from their high-resolution images. To test this technique, we develop 3D hydrodynamical models of supernova remnants and generate synthetic x-ray images from them. Power spectra extracted from both the 3D models and the synthetic images exhibit the same dominant angular scale, which separates large scale features from small scale features due to hydrodynamic instabilities. The angular power spectrum at small length scales during relatively early times is too steep to be consistent with Kolmogorov turbulence, but it transitions to Kolmogorov turbulence at late times. As an example of how this technique can be applied to observations, we extract a power spectrum from a \textit{Chandra} observation of Tycho's SNR and compare with our models. Our predicted power spectrum picks out the angular scale of Tycho's fleece-like structures and also agrees with the small-scale power seen in Tycho. We use this to extract an estimate for the density of the circumstellar gas ($n \sim 0.28/\mathrm{cm^3}$), consistent with previous measurements of this density by other means. The power spectrum also provides an estimate of the density profile of the outermost ejecta. Moreover, we observe additional power at large scales which may provide important clues about the explosion mechanism itself. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 13 pages, 8 figures, submitted to ApJ. Comments are most welcome

Fast Methods for Computing Photometric Variability of Eccentric Binaries: Boosting, Lensing, and Variable Accretion

Authors: Daniel J. D'Orazio, Paul C. Duffell, Christopher Tiede

Abstract: We analyze accretion-rate time series for equal-mass binaries in co-planar gaseous disks spanning a continuous range of orbital eccentricities up to 0.8, for both prograde and retrograde systems. The dominant variability timescales match that of previous investigations; the binary orbital period is dominant for prograde binaries with $e \gtrsim 0.1$, with a 5 times longer "lump" period taking over… ▽ More We analyze accretion-rate time series for equal-mass binaries in co-planar gaseous disks spanning a continuous range of orbital eccentricities up to 0.8, for both prograde and retrograde systems. The dominant variability timescales match that of previous investigations; the binary orbital period is dominant for prograde binaries with $e \gtrsim 0.1$, with a 5 times longer "lump" period taking over for $e\lesssim 0.1$. This lump period fades and drops from 5 times to 4.5 times the binary period as $e$ approaches 0.1, where it vanishes. For retrograde orbits, the binary orbital period dominates at $e \lesssim 0.55$ and is accompanied by a 2 times longer-timescale periodicity at higher eccentricities. The shape of the accretion-rate time series varies with binary eccentricity. For prograde systems, the orientation of an eccentric disk causes periodic trading of accretion between the binary components in a ratio that we report as a function of binary eccentricity. We present a publicly available tool, binlite, that can rapidly ($\lesssim 0.01$~sec) generate templates for the accretion-rate time series, onto either binary component, for choice of binary eccentricity below 0.8. As an example use-case, we build lightcurve models where the accretion rate through the circumbinary disk and onto each binary component sets contributions to the emitted specific flux. We combine these rest-frame, accretion-variability lightcurves with observer-dependent Doppler boosting and binary self-lensing. This allows a flexible approach to generating lightcurves over a wide range of binary and observer parameter space. We envision binlite as the access point to a living database that will be updated with state-of-the-art hydrodynamical calculations as they advance. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Comments: Comments welcome

jetsimpy: A Highly Efficient Hydrodynamic Code for Gamma-ray Burst Afterglow

Authors: Hao Wang, Ranadeep G. Dastidar, Dimitrios Giannios, Paul C. Duffell

Abstract: Gamma-ray burst (GRB) afterglows are emissions from ultrarelativistic blast waves produced by a narrow jet interacting with surrounding matter. Since the first multimessenger observation of a neutron star merger, hydrodynamic modeling of GRB afterglows for structured jets with smoothly varying angular energy distributions has gained increased interest. While the evolution of a jet is well describe… ▽ More Gamma-ray burst (GRB) afterglows are emissions from ultrarelativistic blast waves produced by a narrow jet interacting with surrounding matter. Since the first multimessenger observation of a neutron star merger, hydrodynamic modeling of GRB afterglows for structured jets with smoothly varying angular energy distributions has gained increased interest. While the evolution of a jet is well described by self-similar solutions in both ultrarelativistic and Newtonian limits, modeling the transitional phase remains challenging. This is due to the nonlinear spreading of a narrow jet to a spherical configuration and the breakdown of self-similar solutions. Analytical models are limited in capturing these nonlinear effects, while relativistic hydrodynamic simulations are computationally expensive, which restricts the exploration of various initial conditions. In this work, we introduce a reduced hydrodynamic model that approximates the blast wave as an infinitely thin two-dimensional surface. Further assuming axial symmetry, this model simplifies the simulation to one dimension and drastically reduces the computational costs. We have compared our modeling to relativistic hydrodynamic simulations and semianalytic methods, and applied it to fit the light curve and flux centroid motion of GRB 170817A. These comparisons demonstrate good agreement and validate our approach. We have developed this method into a numerical tool, \texttt{jetsimpy}, which models the synchrotron GRB afterglow emission from a blast wave with arbitrary angular energy and Lorentz factor distribution. Although the code is built with GRB afterglow in mind, it applies to any relativistic jet. This tool is particularly useful in Markov Chain Monte Carlo studies and is provided to the community. △ Less

Submitted 11 July, 2024; v1 submitted 29 February, 2024; originally announced February 2024.

Comments: 24 pages, 13 figures. Published in ApJS

The Santa Barbara Binary-Disk Code Comparison

Authors: Paul C. Duffell, Alexander J. Dittmann, Daniel J. D'Orazio, Alessia Franchini, Kaitlin M. Kratter, Anna B. T. Penzlin, Enrico Ragusa, Magdalena Siwek, Christopher Tiede, Haiyang Wang, Jonathan Zrake, Adam M. Dempsey, Zoltan Haiman, Alessandro Lupi, Michal Pirog, Geoffrey Ryan

Abstract: We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed… ▽ More We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, $0.2 (a^2 Ω_B /ν)$ binary orbits, where $ν$ is the kinematic viscosity. We did not investigate dependence on binary mass ratio, eccentricity, disk temperature, or disk viscosity; therefore, these benchmarks may act as guides towards expanding converged solutions to the wider parameter space but might need to be updated in a future study that investigates dependence on system parameters. We find the most major discrepancies between codes resulted from the dimensionality of the setup (3D vs 2D disks). Beyond this, we find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D. △ Less

Submitted 20 June, 2024; v1 submitted 20 February, 2024; originally announced February 2024.

Comments: ApJ Accepted

Disk induced binary precession: Implications for dynamics and multi-messenger observations of black hole binaries

Authors: Christopher Tiede, Daniel J. D'Orazio, Lorenz Zwick, Paul C. Duffell

Abstract: Many studies have recently documented the orbital response of eccentric binaries accreting from thin circumbinary disks, characterizing the change in binary semi-major axis and eccentricity. We extend these calculations to include the precession of the binary's longitude of periapse induced by the circumbinary disk, and we characterize this precession continuously with binary eccentricity $e_b$ fo… ▽ More Many studies have recently documented the orbital response of eccentric binaries accreting from thin circumbinary disks, characterizing the change in binary semi-major axis and eccentricity. We extend these calculations to include the precession of the binary's longitude of periapse induced by the circumbinary disk, and we characterize this precession continuously with binary eccentricity $e_b$ for equal mass components. This disk-induced apsidal precession is prograde with a weak dependence on binary eccentricity when $e_b \lesssim 0.4$ and decreases approximately linearly for $e_b \gtrsim 0.4$; yet at all $e_b$ binary precession is faster than the rates of change to the semi-major axis and eccentricity by an order of magnitude. We estimate that such precession effects are likely most important for sub-parsec separated binaries with masses $\lesssim 10^7 M_\odot$, like LISA precursors. We find that accreting, equal-mass LISA binaries with $M < 10^6 M_\odot$ (and the most massive $M \sim 10^7 M_\odot$ binaries out to $z \sim 3$) may acquire a detectable phase offset due to the disk-induced precession. Moreover, disk-induced precession can compete with General Relativistic precession in vacuum, making it important for observer-dependent electromagnetic searches for accreting massive binaries -- like Doppler boost and binary self-lensing models -- after potentially only a few orbital periods. △ Less

Submitted 9 April, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 16 pages and 10 figures

Journal ref: The Astrophysical Journal, Volume 964, pg 46, March 20 2024

The effect of relativistic precession on light curves of tidal disruption events

Authors: Diego Calderón, Ondřej Pejcha, Brian D. Metzger, Paul C. Duffell

Abstract: The disruption of a star by the tidal forces of a spinning black hole causes the stellar stream to precess affecting the conditions for triggering the tidal disruption event (TDE). In this work, we study the effect that precession imprints on TDE light curves due to the interaction of the TDE wind and luminosity with the stream wrapped around the black hole. We perform two-dimensional radiation-hy… ▽ More The disruption of a star by the tidal forces of a spinning black hole causes the stellar stream to precess affecting the conditions for triggering the tidal disruption event (TDE). In this work, we study the effect that precession imprints on TDE light curves due to the interaction of the TDE wind and luminosity with the stream wrapped around the black hole. We perform two-dimensional radiation-hydrodynamic simulations using the moving-mesh hydrodynamic code JET with its radiation treatment module. We study the impact of black hole mass, accretion efficiency, and inclination between the orbital and spin planes. From our results, we identified two behaviours: $i)$ models with low-mass black holes ($M_\text{h}\sim10^6~\text{M}_{\odot}$), low inclination ($i\sim0$), and low accretion efficiency ($η\sim0.01$) show light curves with a short early peak caused by the interaction of the wind with the inner edge of the stream. The line of sight has little effect on the light curve, since the stream covers a small fraction of the solid angle due to the precession occurring in the orbital plane; $ii)$ models with high-mass black holes ($M_\text{h}\gtrsim10^7~\text{M}_{\odot}$), high inclination ($i\sim90^{\circ}$), and high accretion efficiency ($η\sim0.1$) produce light curves with luminosity peaks that can be delayed by up to 50-100 d depending on the line of sight due to presence of the precessed stream blocking the radiation in the early phase of the event. Our results show that black hole spin and misalignment do not imprint recognisable features on the light curves but rather can add complications to their analysis. △ Less

Submitted 18 January, 2024; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: 19 pages (+1 in Appendix), 13 figures (+1 in Appendix), 1 table. Accepted for publication in MNRAS

SPROUT: A moving mesh hydro code using a uniformly expanding Cartesian grid

Authors: Soham Mandal, Paul C. Duffell

Abstract: We present the publicly available moving-mesh hydrodynamics code Sprout. Sprout solves the equations of ideal hydrodynamics on an expanding Cartesian mesh. The expanding mesh can follow fluid outflows for several orders of magnitude with very little numerical diffusion, thereby capturing shocks and fine structures accurately. Following the bulk flow accurately also allows for longer timesteps in g… ▽ More We present the publicly available moving-mesh hydrodynamics code Sprout. Sprout solves the equations of ideal hydrodynamics on an expanding Cartesian mesh. The expanding mesh can follow fluid outflows for several orders of magnitude with very little numerical diffusion, thereby capturing shocks and fine structures accurately. Following the bulk flow accurately also allows for longer timesteps in general. This makes Sprout particularly suitable for studying expanding outflows such as supernova remnants and active galactic nuclei. Relative to other moving mesh codes, the simple mesh structure in Sprout is also convenient for implementing additional physics or algorithms. Many code tests are performed to test the accuracy and performance of the numerical scheme. △ Less

Submitted 6 October, 2023; v1 submitted 25 July, 2023; originally announced July 2023.

Comments: 13 pages, 12 figures; accepted for publication in ApJS

A 3D Numerical Study of Anisotropies in Supernova Remnants

Authors: Soham Mandal, Paul C. Duffell, Abigail Polin, Dan Milisavljevic

Abstract: We develop a suite of 3D hydrodynamic models of supernova remnants (SNRs) expanding against the circumstellar medium (CSM). We study the Rayleigh-Taylor Instability (RTI) forming at the expansion interface by calculating an angular power spectrum for each of these models. The power spectra of young SNRs is seen to exhibit a dominant angular mode, which is a diagnostic of their ejecta density profi… ▽ More We develop a suite of 3D hydrodynamic models of supernova remnants (SNRs) expanding against the circumstellar medium (CSM). We study the Rayleigh-Taylor Instability (RTI) forming at the expansion interface by calculating an angular power spectrum for each of these models. The power spectra of young SNRs is seen to exhibit a dominant angular mode, which is a diagnostic of their ejecta density profile as found by previous studies. The steep scaling of power at smaller modes and the time evolution of the spectra is indicative of absence of a turbulent cascade. Instead, as the time evolution of the spectra suggests, they may be governed by an angular mode dependent net growth rate. We also study the impact of anisotropies in the ejecta as well as in the CSM on the power spectra of velocity and density. We confirm that perturbations in the density field (whether imposed on the ejecta or the CSM) do not influence the anisotropy of the remnant significantly unless they have a very large amplitude and form large-scale coherent structures. In any case, these clumps can only affect structures on large angular scales. The power spectra on small angular scales is completely independent of the initial clumpiness and only governed by the growth and saturation of the Rayleigh-Taylor instability. △ Less

Submitted 6 October, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 13 pages, 9 figures; accepted for publication in ApJ

Hydrodynamic mixing of accretion disk outflows in collapsars: implications for r-process signatures

Authors: Jennifer Barnes, Paul C. Duffell

Abstract: The astrophysical environments capable of triggering heavy-element synthesis via rapid neutron capture (the r-process) remain uncertain. While binary neutron star mergers (NSMs) are known to forge r-process elements, certain rare supernovae (SNe) have been theorized to supplement, or even dominate, r-production by NSMs. However, the most direct evidence for such SNe, unusual reddening of the emiss… ▽ More The astrophysical environments capable of triggering heavy-element synthesis via rapid neutron capture (the r-process) remain uncertain. While binary neutron star mergers (NSMs) are known to forge r-process elements, certain rare supernovae (SNe) have been theorized to supplement, or even dominate, r-production by NSMs. However, the most direct evidence for such SNe, unusual reddening of the emission caused by the high opacities of r-process elements, has not been observed. Recent work identified the distribution of r-process material within the SN ejecta as a key predictor of the ease with which signals associated with r-process enrichment could be discerned. Though this distribution results from hydrodynamic processes at play during the SN explosion, thus far it has been treated only in a parameterized way. We use hydrodynamic simulations to model how disk winds, the alleged locus of r-production in rare SNe, mix with initially r-process-free ejecta. We study mixing as a function of the wind mass and duration and of the initial SN explosion energy, and find that it increases with the first two of these and decreases with the third. This suggests that SNe accompanying the longest long-duration gamma-ray bursts are promising places to search for signs of r-process enrichment. We use semianalytic radiation transport to connect hydrodynamics to electromagnetic observables, allowing us to assess the mixing level at which the presence of r-process material can be diagnosed from SN light curves. Analytic arguments constructed atop this foundation imply that a wind-driven r-process-enriched SN model is unlikely to explain standard energetic SNe. △ Less

Submitted 7 June, 2023; v1 submitted 28 April, 2023; originally announced May 2023.

Comments: V2 accepted to ApJ. Absolute banger. The authors encourage the widest possible adoption of the term "collapsaR-process.''

Scary Barbie: An Extremely Energetic, Long-Duration Tidal Disruption Event Candidate Without a Detected Host Galaxy at z = 0.995

Authors: Bhagya M. Subrayan, Dan Milisavljevic, Ryan Chornock, Raffaella Margutti, Kate D. Alexander, Vandana Ramakrishnan, Paul C. Duffell, Danielle A. Dickinson, Kyoung-Soo Lee, Dimitrios Giannios, Geoffery Lentner, Mark Linvill, Braden Garretson, Matthew J. Graham, Daniel Stern, Daniel Brethauer, Tien Duong, Wynn Jacobson-Gal�n, Natalie LeBaron, David Matthews, Huei Sears, Padma Venkatraman

Abstract: We report multi-wavelength observations and characterization of the ultraluminous transient AT 2021lwx (ZTF20abrbeie; aka ``Barbie'') identified in the alert stream of the Zwicky Transient Facility (ZTF) using a Recommender Engine For Intelligent Transient Tracking (REFITT) filter on the ANTARES alert broker. From a spectroscopically measured redshift of 0.995, we estimate a peak observed pseudo-b… ▽ More We report multi-wavelength observations and characterization of the ultraluminous transient AT 2021lwx (ZTF20abrbeie; aka ``Barbie'') identified in the alert stream of the Zwicky Transient Facility (ZTF) using a Recommender Engine For Intelligent Transient Tracking (REFITT) filter on the ANTARES alert broker. From a spectroscopically measured redshift of 0.995, we estimate a peak observed pseudo-bolometric luminosity of log (L$_{\text{max}} / [\text{erg}/\text{s}]$) = 45.7 from slowly fading ztf-$\it{g}$ and ztf-$r$ light curves spanning over 1000 observer-frame days. The host galaxy is not detected in archival Pan-STARRS observations ($g > 23.3$ mag), implying a lower limit to the outburst amplitude of more than 5 mag relative to the quiescent host galaxy. Optical spectra from Lick and Keck Observatories exhibit strong emission lines with narrow cores from the H Balmer series and ultraviolet semi-forbidden lines of Si III] $λ$1892, C III] $λ$1909, and C II] $λ$2325. Typical nebular lines in AGN spectra from ions such as [O II] and [O III] are not detected. These spectral features, along with the smooth light curve that is unlike most AGN flaring activity, and the luminosity that exceeds any observed or theorized supernova, lead us to conclude that AT 2021lwx is most likely an extreme tidal disruption event (TDE). Modeling of ZTF photometry with MOSFiT suggests that the TDE was between a $\approx 14 M_{\odot}$ star and a supermassive black hole of mass $M_{\text{BH}} \sim$ $10^{8} M_{\odot}$. Continued monitoring of the still-evolving light curve along with deep imaging of the field once AT 2021lwx has faded can test this hypothesis and potentially detect the host galaxy. △ Less

Submitted 8 June, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: 15 pages, 4 figures, 1 Table; Version as published in The Astrophysical Journal Letters. Observations of AT 2021lwx published in the paper can be found at https://bsubraya.github.io/research/

Numerical Investigation of Dynamical and Morphological Trends in Relativistic Jets

Authors: Soham Mandal, Paul C. Duffell, Yuan Li

Abstract: Active galactic nuclei (AGN) show a range of morphologies and dynamical properties, which are determined not only by parameters intrinsic to the central engine but also their interaction with the surrounding environment. We investigate the connection of kiloparsec scale AGN jet properties to their intrinsic parameters and surroundings. This is done using a suite of 40 relativistic hydrodynamic sim… ▽ More Active galactic nuclei (AGN) show a range of morphologies and dynamical properties, which are determined not only by parameters intrinsic to the central engine but also their interaction with the surrounding environment. We investigate the connection of kiloparsec scale AGN jet properties to their intrinsic parameters and surroundings. This is done using a suite of 40 relativistic hydrodynamic simulations spanning a wide range of engine luminosities and opening angles. We explore AGN jet propagation with different ambient density profiles, including $r^{-2}$ (self-similar solution) and $r^{-1}$, which is more relevant for AGN host environments. The Fanaroff-Riley (FR) morphological dichotomy arises naturally in our models. Jets with low energy density compared to the ambient medium produce a center-brightened emissivity distribution, while emissivity from relatively higher energy density jets is dominated by a terminal bright spot. We observe recollimation shocks in our simulations that can generate bright spots along the spine of the jet, providing a possible explanation for "knots" observed in AGN jets. We additionally find a scaling relation between the number of knots and the jet-head-to-surroundings energy density ratio. This scaling relation is generally consistent with the observations of the jets in M87 and Cygnus A. Our model also correctly predicts M87 as FR I and Cygnus A as FR II. Our model can be used to relate jet dynamical parameters such as jet head velocity, jet opening angle, and external pressure to jet power and ambient density estimates. △ Less

Submitted 14 August, 2022; v1 submitted 26 March, 2022; originally announced March 2022.

Comments: 12 pages, 7 figures and 2 tables. Published in ApJ

Journal ref: ApJ Vol. 935, No. 1, Jan 2022, pp 42

Moving-mesh radiation-hydrodynamic simulations of wind-reprocessed transients

Authors: Diego Calderón, Ondřej Pejcha, Paul C. Duffell

Abstract: Motivated by recent theoretical work on tidal disruption events and other peculiar transients, we present moving-mesh radiation-hydrodynamic simulations of radiative luminosity emitted by a central source being reprocessed by a wind-like outflow. We couple the moving-mesh hydrodynamic code JET with our newly-developed radiation module based on mixed-frame grey flux-limited diffusion with implicit… ▽ More Motivated by recent theoretical work on tidal disruption events and other peculiar transients, we present moving-mesh radiation-hydrodynamic simulations of radiative luminosity emitted by a central source being reprocessed by a wind-like outflow. We couple the moving-mesh hydrodynamic code JET with our newly-developed radiation module based on mixed-frame grey flux-limited diffusion with implicit timestep update. This allows us to study the self-consistent multi-dimensional radiation-hydrodynamic evolution over more than ten orders of magnitude in both space and time in a single run. We simulate an optically-thick spherical wind with constant or evolving mass-loss rate, which is irradiated by a central isotropic or angularly-dependent radiation source. Our spherically-symmetric simulations confirm previous analytic results by identifying different stages of radiation reprocessing: radiation trapped in the wind, diffusing out through the wind, and reaching constant maximum attenuation. We find that confining the central radiation source in a cone with moderate opening angles decrease up to one order of magnitude the early flux along sightlines oriented away from the direction of radiation injection but that the reprocessed radiation becomes isotropic roughly after one lateral diffusion time through the ejecta. We discuss further applications and guidelines for the use of our novel radiation-hydrodynamics tool in the context of transient modelling. △ Less

Submitted 29 July, 2021; v1 submitted 18 May, 2021; originally announced May 2021.

Comments: 10 pages (+6 in appendix), 6 figures (+3 in appendix), 1 table. Accepted for publication in MNRAS

Proto-magnetar jets as central engines for broad-lined type Ic supernovae

Authors: Swapnil Shankar, Philipp Mösta, Jennifer Barnes, Paul C. Duffell, Daniel Kasen

Abstract: A subset of type Ic supernovae (SNe Ic), broad-lined SNe Ic (SNe Ic-bl), show unusually high kinetic energies ($\sim 10^{52}$ erg) which cannot be explained by the energy supplied by neutrinos alone. Many SNe Ic-bl have been observed in coincidence with long gamma-ray bursts (GRBs) which suggests a connection between SNe and GRBs. A small fraction of core-collapse supernovae (CCSNe) form a rapidly… ▽ More A subset of type Ic supernovae (SNe Ic), broad-lined SNe Ic (SNe Ic-bl), show unusually high kinetic energies ($\sim 10^{52}$ erg) which cannot be explained by the energy supplied by neutrinos alone. Many SNe Ic-bl have been observed in coincidence with long gamma-ray bursts (GRBs) which suggests a connection between SNe and GRBs. A small fraction of core-collapse supernovae (CCSNe) form a rapidly-rotating and strongly-magnetized protoneutron star (PNS), a proto-magnetar. Jets from such magnetars can provide the high kinetic energies observed in SNe Ic-bl and also provide the connection to GRBs. In this work we use the jetted outflow produced in a 3D CCSN simulation from a consistently formed proto-magnetar as the central engine for full-star explosion simulations. We extract a range of central engine parameters and find that the extracted engine energy is in the range of $6.231 \times 10^{51}-1.725 \times 10^{52}$ erg, the engine time-scale in the range of $0.479-1.159$ s and the engine half-opening angle in the range of $\sim 9-19^{\circ}$. Using these as central engines, we perform 2D special-relativistic (SR) hydrodynamic (HD) and radiation transfer simulations to calculate the corresponding light curves and spectra. We find that these central engine parameters successfully produce SNe Ic-bl which demonstrates that jets from proto-magnetars can be viable engines for SNe Ic-bl. We also find that only the central engines with smaller opening angles ($\sim 10^{\circ}$) form a GRB implying that GRB formation is likely associated with narrower jet outflows and Ic-bl's without GRBs may be associated with wider outflows. △ Less

Submitted 17 May, 2021; originally announced May 2021.

Comments: 13 pages, 12 figures

Orbital Evolution of Equal-mass Eccentric Binaries due to a Gas Disk: Eccentric Inspirals and Circular Outspirals

Authors: Daniel J. D'Orazio, Paul C. Duffell

Abstract: We solve the equations of two-dimensional hydrodynamics describing a circumbinary disk accreting onto an eccentric, equal-mass binary. We compute the time rate of change of the binary semimajor axis $a$ and eccentricity $e$ over a continuous range of eccentricities spanning $e=0$ to $e=0.9$. We find that binaries with initial eccentricities $e_0 \lesssim 0.1$ tend to $e=0$, where the binary semima… ▽ More We solve the equations of two-dimensional hydrodynamics describing a circumbinary disk accreting onto an eccentric, equal-mass binary. We compute the time rate of change of the binary semimajor axis $a$ and eccentricity $e$ over a continuous range of eccentricities spanning $e=0$ to $e=0.9$. We find that binaries with initial eccentricities $e_0 \lesssim 0.1$ tend to $e=0$, where the binary semimajor axis expands. All others are attracted to $e \approx 0.4$, where the binary semimajor axis decays. The $e \approx 0.4$ attractor is caused by a rapid change in the disk response from a nearly origin-symmetric state to a precessing asymmetric state. The state change causes the time rates of change $\dot{a}$ and $\dot{e}$ to steeply change sign at the same critical eccentricity resulting in an attracting solution where $\dot{a} = \dot{e} = 0$. This does not, however, result in a stalled, eccentric binary. The finite-transition time between disk states causes the binary eccentricity to evolve beyond the attracting eccentricity in both directions resulting in oscillating orbital parameters and a drift of the semimajor axis. For the chosen disk parameters, binaries with $e_0 \gtrsim 0.1$ evolve toward and then oscillate around $e \approx 0.4$ where they shrink in semimajor axis. Because unequal mass binaries grow toward equal mass through preferential accretion, our results are applicable to a wide range of initial binary mass ratios. Hence, these findings merit further investigations of this disk transition; understanding its dependence on disk parameters is vital for determining the fate of binaries undergoing orbital evolution with a circumbinary disk. △ Less

Submitted 15 June, 2021; v1 submitted 16 March, 2021; originally announced March 2021.

Comments: Published in The Astrophysical Journal Letters

Journal ref: ApJL 914 L21 (2021)

The Effect of Jet-Ejecta Interaction on the Viewing Angle Dependence of Kilonova Light Curves

Authors: Hannah Klion, Paul C. Duffell, Daniel Kasen, Eliot Quataert

Abstract: The merger of two neutron stars produces an outflow of radioactive heavy nuclei. Within a second of merger, the central remnant is expected to also launch a relativistic jet, which shock-heats and disrupts a portion of the radioactive ejecta. Within a few hours, emission from the radioactive material gives rise to an ultraviolet, optical, and infrared transient (a kilonova). We use the endstates o… ▽ More The merger of two neutron stars produces an outflow of radioactive heavy nuclei. Within a second of merger, the central remnant is expected to also launch a relativistic jet, which shock-heats and disrupts a portion of the radioactive ejecta. Within a few hours, emission from the radioactive material gives rise to an ultraviolet, optical, and infrared transient (a kilonova). We use the endstates of a suite of 2D relativistic hydrodynamic simulations of jet-ejecta interaction as initial conditions for multi-dimensional Monte Carlo radiation transport simulations of the resulting viewing angle-dependent light curves and spectra starting at 1.5h after merger. We find that on this timescale, jet shock heating does not affect the kilonova emission. However, the jet disruption to the density structure of the ejecta does change the light curves. The jet carves a channel into the otherwise spheroidal ejecta, revealing the hot, inner regions. As seen from near ($\lesssim 30 �$) the jet axis, the kilonova is brighter by a factor of a few and bluer. The strength of this effect depends on the jet parameters, since the light curves of more heavily disrupted ejecta are more strongly affected. The light curves and spectra are also more heavily modified in the ultraviolet than in the optical. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: 11 pages, 10 figures. Submitted to MNRAS

Circumbinary Disks: Accretion and Torque as a Function of Mass Ratio and Disk Viscosity

Authors: Paul C. Duffell, Daniel D'Orazio, Andrea Derdzinski, Zoltan Haiman, Andrew MacFadyen, Anna L. Rosen, Jonathan Zrake

Abstract: Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with previous studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary towards equal mass. Accretion variability depends sen… ▽ More Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with previous studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary towards equal mass. Accretion variability depends sensitively on the numerical sink prescription, but the torque and relative accretion onto each component do not depend on the sink timescale. Positive torque and highly variable accretion occurs only for mass ratios greater than around $0.05$. This means that for mass ratios below $0.05$, the binary would migrate inward until the secondary accreted sufficient mass, after which it would execute a U-turn and migrate outward. We explore a range of viscosities, from $α= 0.03$ to $α= 0.15$, and find that this outward torque is proportional to the viscous torque, simply proportional to viscosity in this range. Dependence of accretion and torque on mass ratio is explored in detail, densely sampling mass ratios between $0.01$ and unity. For mass ratio $q > 0.6$, accretion variability is found to exhibit a distinct sawtooth pattern, typically with a five-orbit cycle that provides a "smoking gun" prediction for variable quasars observed over long periods, as a potential means to confirm the presence of a binary. △ Less

Submitted 31 July, 2020; v1 submitted 13 November, 2019; originally announced November 2019.

Comments: ApJ Accepted

How Dense a CSM is Sufficient to Choke a Jet?

Authors: Paul C. Duffell, Anna Y. Q. Ho

Abstract: The progenitor stars of stripped-envelope high-velocity supernovae (Ic-BL SNe) can explode inside a dense circumstellar medium (CSM) that extends out to many times the progenitor radius. This complicates the question of whether all Ic-BL SNe harbor a jet, which can tunnel through the star and be viewed on-axis as a long-duration gamma-ray burst (GRB). More specifically, a sufficiently dense CSM mi… ▽ More The progenitor stars of stripped-envelope high-velocity supernovae (Ic-BL SNe) can explode inside a dense circumstellar medium (CSM) that extends out to many times the progenitor radius. This complicates the question of whether all Ic-BL SNe harbor a jet, which can tunnel through the star and be viewed on-axis as a long-duration gamma-ray burst (GRB). More specifically, a sufficiently dense CSM might "choke" the jet, redistributing its energy quasi-spherically. In this study, we numerically calculate the CSM density necessary for jet-choking. For typical GRBs, we determine the jet is not choked in the CSM unless $ρr^2 > 4 \times 10^{19}$ g/cm; this requires several solar masses of CSM to be situated within $10^{13}$ cm of the progenitor, a much higher density than any CSM observed. We conclude that typical GRB jets are not choked in the CSM. However, in many cases the CSM has sufficient mass to decelerate the jet to a modest Lorentz factor ($Γ\sim 10$), which should lead to a long coasting phase for the jet, observable as a long plateau (potentially up to a few days) in the afterglow light curve. For extreme cases of low-energy GRBs in a high-mass CSM, the jet will decelerate to nonrelativistic velocities, causing it to spread modestly to a larger opening angle ($θ_j \approx 20$ degrees) before breaking out of the CSM. Even in these extreme examples, the jet does not have time to redistribute its energy quasi-spherically in the CSM before breakout. △ Less

Submitted 8 July, 2019; originally announced July 2019.

An empirically-derived formula for the shape of planet-induced gaps in protoplanetary disks

Authors: Paul C. Duffell

Abstract: This study uses numerical hydrodynamics calculations and a novel method for densely sampling parameter space to measure the precise shape of a gap opened by a planet in a gaseous disk, as a function of planet-to-star mass ratio, disk Mach number, and disk viscosity. Formulas for gap depth and width are determined, which are combined to form a complete formula for surface density as a function of r… ▽ More This study uses numerical hydrodynamics calculations and a novel method for densely sampling parameter space to measure the precise shape of a gap opened by a planet in a gaseous disk, as a function of planet-to-star mass ratio, disk Mach number, and disk viscosity. Formulas for gap depth and width are determined, which are combined to form a complete formula for surface density as a function of radius in the disk. This new analytical formula is compared with numerically-derived gaps opened by planets ranging from very low masses up to a few times Jupiter's mass, and excellent agreement is found over a wide range of parameter space relevant to planet-disk interactions. A simple-to-use code is presented to rapidly generate synthetic disk profiles. △ Less

Submitted 26 June, 2019; originally announced June 2019.

Gas Flows Within Cavities of Circumbinary Discs in Eccentric Binary Protostellar Systems

Authors: Philipp M�sta, Ronald E. Taam, Paul C. Duffell

Abstract: The structure and evolution of gas flows within the cavity of a circumbinary disk (CBD) surrounding the stellar components in eccentric binaries are examined via two-dimensional hydrodynamical simulations. The degree to which gas fills the cavity between the circumstellar disks (CSDs) and the CBD is found to be greater for highly eccentric systems, in comparison to low-eccentricity systems, reflec… ▽ More The structure and evolution of gas flows within the cavity of a circumbinary disk (CBD) surrounding the stellar components in eccentric binaries are examined via two-dimensional hydrodynamical simulations. The degree to which gas fills the cavity between the circumstellar disks (CSDs) and the CBD is found to be greater for highly eccentric systems, in comparison to low-eccentricity systems, reflecting the spatial extent over which mass enters into the cavity throughout the orbit. The pattern of the gas flow in the cavity differs for eccentric binaries from that of binaries in a circular orbit. In particular, the former reveals tightly wound gas streams and figure-eight-like structures for systems characterized by eccentricities, $e \ge 0.4$, whereas the latter only reveal relatively loosely bent streams from the CBD to the CSDs. Hence, the description of the stream structures can be a probe of sufficient non-circularity of the binary orbital motion. Given that the inner edge of the CBD is not very well defined for highly eccentric systems due to the complex gas structures, it is suggested that the area of the cavity for high-sensitivity imaging observations may prove to be a more useful diagnostic for probing the effectiveness of CBD clearing in the future. △ Less

Submitted 22 April, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 8 pages, 5 figures

Journal ref: ApJL, 875, L21 (2019)

Jet Dynamics in Compact Object Mergers: GW170817 Likely had a Successful Jet

Authors: Paul C. Duffell, Eliot Quataert, Daniel Kasen, Hannah Klion

Abstract: We use relativistic hydrodynamic numerical calculations to study the interaction between a jet and a homologous outflow produced dynamically during binary neutron star mergers. We quantify how the thermal energy supplied by the jet to the ejecta and the ability of a jet to escape the homologous ejecta depend on the parameters of the jet engine and the ejecta. For collimated jets initiated at early… ▽ More We use relativistic hydrodynamic numerical calculations to study the interaction between a jet and a homologous outflow produced dynamically during binary neutron star mergers. We quantify how the thermal energy supplied by the jet to the ejecta and the ability of a jet to escape the homologous ejecta depend on the parameters of the jet engine and the ejecta. For collimated jets initiated at early times compared to the engine duration, we show that successful breakout of the forward cocoon shock necessitates a jet that successfully escapes the ejecta. This is because the ejecta is expanding and absorbing thermal energy, so that the forward shock from a failed jet stalls before it reaches the edge of the ejecta. This conclusion can be circumvented only for very energetic wide angle jets, with parameters that are uncomfortable given short-duration GRB observations. For successful jets, we find two regimes of jet breakout from the ejecta, early breakout on timescales shorter than the engine duration, and late breakout well after the engine shuts off. A late breakout can explain the observed delay between gravitational waves and gamma rays in GW 170817. We show that for the entire parameter space of jet parameters surveyed here (covering energies $\sim 10^{48}-10^{51}$ ergs and opening angles $θ_j \sim 0.07-0.4$) the thermal energy deposited into the ejecta by the jet propagation is less than that produced by r-process heating on second timescales by a factor of $\gtrsim 10$. Shock heating is thus energetically subdominant in setting the luminosity of thermally powered transients coincident with neutron star mergers (kilonovae). For typical short GRB jet parameters, our conclusion is stronger: there is little thermal energy in the cocoon, much less than what is needed to explain the early blue component of the kilonova in GW 170817. △ Less

Submitted 27 June, 2018; originally announced June 2018.

Interaction of a Supernova with a Circumstellar Disk

Authors: Austin T. McDowell, Paul C. Duffell, Daniel Kasen

Abstract: Interaction between supernova (SN) ejecta and a dense circumstellar medium (CSM) can power a luminous light curve and create narrow emission lines in the spectra. While theoretical studies of interaction often assume a spherically symmetric CSM, there are observational indications that the gas surrounding some SN has a disk-like geometry. Here, we use moving-mesh hydrodynamics simulations to study… ▽ More Interaction between supernova (SN) ejecta and a dense circumstellar medium (CSM) can power a luminous light curve and create narrow emission lines in the spectra. While theoretical studies of interaction often assume a spherically symmetric CSM, there are observational indications that the gas surrounding some SN has a disk-like geometry. Here, we use moving-mesh hydrodynamics simulations to study the interaction of a SN with a disk and determine how the dynamics and observable signatures may depend on the disk mass, thickness, and radial extent. We find that simple modifications to standard spherically-symmetric scaling laws can be used to describe the propagation and heating rate of the interaction shock. We use the resulting shock heating rates to derive approximate bolometric light curves, and provide analytic formulas that can be used to generate simple synthetic light curves for general supernova-disk interactions. For certain disk parameters and explosion energies, we are able to produce luminosities akin to those seen in super-luminous SN. Because the SN ejecta can flow around and engulf the CSM disk, the interaction region may become embedded and from certain viewing angles the narrow emission lines indicative of interaction may be hidden. △ Less

Submitted 14 February, 2018; originally announced February 2018.

On the Deceleration and Spreading of Relativistic Jets I: Jet Dynamics

Authors: Paul C. Duffell, Tanmoy Laskar

Abstract: Jet breaks in gamma ray burst (GRB) afterglows provide a direct probe of their collimation angle. Modeling a jet break requires an understanding of the "jet spreading" process, whereby the jet transitions from a collimated outflow into the spherical Sedov-Taylor solution at late times. Currently, direct numerical calculations are the most accurate way to capture the deceleration and spreading proc… ▽ More Jet breaks in gamma ray burst (GRB) afterglows provide a direct probe of their collimation angle. Modeling a jet break requires an understanding of the "jet spreading" process, whereby the jet transitions from a collimated outflow into the spherical Sedov-Taylor solution at late times. Currently, direct numerical calculations are the most accurate way to capture the deceleration and spreading process, as analytical models have previously given inaccurate descriptions of the dynamics. Here (in paper I) we present a new, semi-analytical model built empirically by performing relativistic numerical jet calculations and measuring the relationship between Lorentz factor and opening angle. We then calculate the Lorentz factor and jet opening angle as a function of shock radius and compare to the numerical solutions. Our analytic model provides an efficient means of computing synthetic GRB afterglow light curves and spectra, which is the focus of paper II. △ Less

Submitted 19 October, 2017; originally announced October 2017.

Comments: ApJ Submitted

A GRB and Broad-lined Type Ic Supernova from a Single Central Engine

Authors: Jennifer Barnes, Paul C. Duffell, Yuqian Liu, Maryam Modjaz, Federica B. Bianco, Daniel Kasen, Andrew I. MacFadyen

Abstract: Unusually high velocities (< ~0.1c) and correspondingly high kinetic energies have been observed in a subset of Type Ic supernovae (so-called "broad-lined Ic" supernovae; SNe Ic-BL), prompting a search for a central engine model capable of generating such energetic explosions. A clue to the explosion mechanism may lie in the fact that all supernovae that accompany long-duration gamma-ray bursts be… ▽ More Unusually high velocities (< ~0.1c) and correspondingly high kinetic energies have been observed in a subset of Type Ic supernovae (so-called "broad-lined Ic" supernovae; SNe Ic-BL), prompting a search for a central engine model capable of generating such energetic explosions. A clue to the explosion mechanism may lie in the fact that all supernovae that accompany long-duration gamma-ray bursts belong to the SN Ic-BL class. Using a combination of two-dimensional relativistic hydrodynamics and radiation transport calculations, we demonstrate that the central engine responsible for long gamma-ray bursts can also trigger a SN Ic-BL. We find that a reasonable gamma-ray burst engine injected into a stripped Wolf-Rayet progenitor produces a relativistic jet with energy ~10^51 ergs, as well as a SN whose synthetic light curves and spectra are fully consistent with observed SNe Ic-BL during the photospheric phase. As a result of the jet's asymmetric energy injection, the SN spectra and light curves depend on viewing angle. The impact of viewing angle on the spectrum is particularly pronounced at early times, while the viewing angle dependence for the light curves (~10% variation in bolometric luminosity) persists throughout the photospheric phase. △ Less

Submitted 8 August, 2017; originally announced August 2017.

Comments: Submitted to ApJ; comments welcome

DISCO: A 3D Moving-Mesh Magnetohydrodynamics Code Designed for the Study of Astrophysical Disks

Authors: Paul C. Duffell

Abstract: This work presents the publicly available moving-mesh magnetohydrodynamics code DISCO. DISCO is efficient and accurate at evolving orbital fluid motion in two and three dimensions, especially at high Mach number. DISCO employs a moving-mesh approach utilizing a dynamic cylindrical mesh that can shear azimuthally to follow the orbital motion of the gas. The moving mesh removes diffusive advection e… ▽ More This work presents the publicly available moving-mesh magnetohydrodynamics code DISCO. DISCO is efficient and accurate at evolving orbital fluid motion in two and three dimensions, especially at high Mach number. DISCO employs a moving-mesh approach utilizing a dynamic cylindrical mesh that can shear azimuthally to follow the orbital motion of the gas. The moving mesh removes diffusive advection errors and allows for longer timesteps than a static grid. Magnetohydrodynamics is implemented in DISCO using an HLLD Riemann solver and a novel constrained transport scheme which is compatible with the mesh motion. DISCO is tested against a wide variety of problems, which are designed to test its stability, accuracy and scalability. In addition, several magnetohydrodynamics tests are performed which demonstrate the accuracy and stability of the new constrained transport approach, including two tests of the magneto-rotational instability (MRI); one testing the linear growth rate and the other following the instability into the fully turbulent regime. △ Less

Submitted 23 June, 2016; v1 submitted 11 May, 2016; originally announced May 2016.

Comments: Accepted to ApJS

Eccentric Jupiters via Disk-Planet Interactions

Authors: Paul C. Duffell, Eugene Chiang

Abstract: Numerical hydrodynamics calculations are performed to determine conditions under which giant planet eccentricities can be excited by parent gas disks. Unlike in other studies, Jupiter-mass planets are found to have their eccentricities amplified --- provided their orbits start eccentric. We disentangle the web of co-rotation, co-orbital, and external resonances to show that this finite-amplitude i… ▽ More Numerical hydrodynamics calculations are performed to determine conditions under which giant planet eccentricities can be excited by parent gas disks. Unlike in other studies, Jupiter-mass planets are found to have their eccentricities amplified --- provided their orbits start eccentric. We disentangle the web of co-rotation, co-orbital, and external resonances to show that this finite-amplitude instability is consistent with that predicted analytically. Ellipticities can grow until they reach of order the disk's aspect ratio, beyond which the external Lindblad resonances that excite eccentricity are weakened by the planet's increasingly supersonic epicyclic motion. Forcing the planet to still larger eccentricities causes catastrophic eccentricity damping as the planet collides into gap walls. For standard parameters, the range of eccentricities for instability is modest; the threshold eccentricity for growth ($\sim$$0.04$) is not much smaller than the final eccentricity to which orbits grow ($\sim$$0.07$). If this threshold eccentricity can be lowered (perhaps by non-barotropic effects), and if the eccentricity driving documented here survives in 3D, it may robustly explain the low-to-moderate eccentricities $\lesssim 0.1$ exhibited by many giant planets (including Jupiter and Saturn), especially those without planetary or stellar companions. △ Less

Submitted 3 September, 2015; v1 submitted 30 July, 2015; originally announced July 2015.

Comments: Accepted to ApJ with added references and minor revisions

A Narrow Short-Duration GRB Jet from a Wide Central Engine

Authors: Paul C. Duffell, Eliot Quataert, Andrew I. MacFadyen

Abstract: We use two-dimensional relativistic hydrodynamic numerical calculations to show that highly collimated relativistic jets can be produced in neutron star merger models of short-duration gamma ray bursts without the need for a highly directed engine or a large net magnetic flux. Even a hydrodynamic engine generating a very wide sustained outflow on small scales can in principle produce a highly coll… ▽ More We use two-dimensional relativistic hydrodynamic numerical calculations to show that highly collimated relativistic jets can be produced in neutron star merger models of short-duration gamma ray bursts without the need for a highly directed engine or a large net magnetic flux. Even a hydrodynamic engine generating a very wide sustained outflow on small scales can in principle produce a highly collimated relativistic jet, facilitated by a dense surrounding medium which provides a cocoon surrounding the jet core. An oblate geometry to the surrounding gas significantly enhances the collimation process. Previous numerical simulations have shown that the merger of two neutron stars produces an oblate, expanding cloud of dynamical ejecta. We show that this gas can efficiently collimate the central engine power much as the surrounding star does in long-duration GRB models. For typical short-duration GRB central engine parameters, we find jets with opening angles of order 10 degrees in which a large fraction of the total outflow power of the central engine resides in highly relativistic material. These results predict large differences in the opening angles of outflows from binary neutron star mergers versus neutron star-black hole mergers. △ Less

Submitted 20 May, 2015; originally announced May 2015.

A Simple Analytical Model for Gaps in Protoplanetary Disks

Authors: Paul C. Duffell

Abstract: An analytical model is presented for calculating the surface density as a function of radius $Σ(r)$ in protoplanetary disks in which a planet has opened a gap. This model is also applicable to circumbinary disks with extreme binary mass ratios. The gap profile can be solved for algebraically, without performing any numerical integrals. In contrast with previous one-dimensional gap models, this mod… ▽ More An analytical model is presented for calculating the surface density as a function of radius $Σ(r)$ in protoplanetary disks in which a planet has opened a gap. This model is also applicable to circumbinary disks with extreme binary mass ratios. The gap profile can be solved for algebraically, without performing any numerical integrals. In contrast with previous one-dimensional gap models, this model correctly predicts that low-mass (sub-Jupiter) planets can open gaps in sufficiently low-viscosity disks, and it correctly recovers the power-law dependence of gap depth on planet-to-star mass ratio $q$, disk aspect ratio $h/r$, and dimensionless viscosity $α$ found in previous numerical studies. Analytical gap profiles are compared with numerical calculations over a range of parameter space in $q$, $h/r$, and $α$, demonstrating accurate reproduction of the "partial gap" regime, and general agreement over a wide range of parameter space. △ Less

Submitted 13 May, 2015; originally announced May 2015.

A One-Dimensional Model for Rayleigh-Taylor Instability in Supernova Remnants

Authors: Paul C. Duffell

Abstract: This study presents a method for approximating the multidimensional effects of Rayleigh-Taylor instability as a modification of the one-dimensional hydro equations. This modification is similar to the Shakura-Sunyaev α prescription for modeling the coarse-grained effects of turbulence in astrophysical disks. The model introduces several dimensionless tunable parameters that are calibrated by compa… ▽ More This study presents a method for approximating the multidimensional effects of Rayleigh-Taylor instability as a modification of the one-dimensional hydro equations. This modification is similar to the Shakura-Sunyaev α prescription for modeling the coarse-grained effects of turbulence in astrophysical disks. The model introduces several dimensionless tunable parameters that are calibrated by comparing with high-resolution two-dimensional axisymmetric numerical calculations of Rayleigh-Taylor unstable flows. A complete description of the model is presented, along with a handful of test problems that demonstrate the extent to which the one-dimensional model is able to reproduce multidimensional effects. △ Less

Submitted 2 March, 2016; v1 submitted 9 April, 2015; originally announced April 2015.

Comments: ApJ Accepted

Halting Migration: Numerical Calculations of Corotation Torques in the Weakly Nonlinear Regime

Authors: Paul C. Duffell

Abstract: Planets in their formative years can migrate due to the influence of gravitational torques in the protoplanetary disk they inhabit. For low-mass planets in an isothermal disk, it is known that there is a strong negative torque on the planet due to its linear perturbation to the disk, causing fast inward migration. The current investigation demonstrates that in these same isothermal disks, for inte… ▽ More Planets in their formative years can migrate due to the influence of gravitational torques in the protoplanetary disk they inhabit. For low-mass planets in an isothermal disk, it is known that there is a strong negative torque on the planet due to its linear perturbation to the disk, causing fast inward migration. The current investigation demonstrates that in these same isothermal disks, for intermediate-mass planets, there is a strong positive nonlinear corotation torque due to the effects of gas being pulled through a gap on horseshoe orbits. For intermediate-mass planets, this positive torque can partially or completely cancel the linear (Type I) torque, leading to slower or outward migration, even in an isothermal disk. The effect is most significant for Super-Earth and Sub-Jovian planets, during the transition from a low-mass linear perturber to a non-linear gap-opening planet, when the planet has opened a so-called 'partial gap'. In this study, numerical calculations of planet-disk interactions calculate these torques explicitly, and scalings are empirically constructed for migration rates in this weakly nonlinear regime. These results find outward migration is possible for planets with masses in the 20 - 100 Earth Mass range. △ Less

Submitted 4 May, 2015; v1 submitted 27 December, 2014; originally announced December 2014.

Comments: ApJ Accepted

Shallow Cavities in Multiple-Planet Systems

Authors: Paul C. Duffell, Ruobing Dong

Abstract: Large cavities are often observed in protoplanetary disks, which might suggest the presence of planets opening gaps in the disk. Multiple planets are necessary to produce a wide cavity in the gas. However, multiple planets may also be a burden to the carving out of very deep gaps. When additional planets are added to the system, the time-dependent perturbations from these additional satellites can… ▽ More Large cavities are often observed in protoplanetary disks, which might suggest the presence of planets opening gaps in the disk. Multiple planets are necessary to produce a wide cavity in the gas. However, multiple planets may also be a burden to the carving out of very deep gaps. When additional planets are added to the system, the time-dependent perturbations from these additional satellites can stir up gas in the gap, suppressing cavity opening. In this study, we perform two-dimensional numerical hydro calculations of gap opening for single and multiple planets, showing the effect that additional planets have on the gap depths. We show that multiple planets produce much shallower cavities than single planets, so that more massive planets are needed in the multiple-planet case to produce an equivalent gap depth as in the single-planet case. To deplete a gap by a factor of 100 for the parameters chosen in this study, one only requires $M_p \approx 3.5M_J$ in the single-planet case, but much more massive planets, $M_p \approx 7M_J$ are required in the multiple-planet case. This requirement of high-mass planets implies that such planets may be detectable in the next generation of direct imaging projects, in gaps whose depths are constrained to be sufficiently deep by ALMA. △ Less

Submitted 9 April, 2015; v1 submitted 11 December, 2014; originally announced December 2014.

Comments: ApJ Accepted

From Engine to Afterglow: Collapsars Naturally Produce Top-Heavy Jets and Early-Time Plateaus in Gamma Ray Burst Afterglows

Authors: Paul C. Duffell, Andrew I. MacFadyen

Abstract: We demonstrate that the steep decay and long plateau in the early phases of gamma ray burst (GRB) X-ray afterglows are naturally produced in the collapsar model, by a means ultimately related to the dynamics of relativistic jet propagation through a massive star. We present two-dimensional axisymmetric hydrodynamical simulations which start from a collapsar engine and evolve all the way through th… ▽ More We demonstrate that the steep decay and long plateau in the early phases of gamma ray burst (GRB) X-ray afterglows are naturally produced in the collapsar model, by a means ultimately related to the dynamics of relativistic jet propagation through a massive star. We present two-dimensional axisymmetric hydrodynamical simulations which start from a collapsar engine and evolve all the way through the late afterglow phase. The resultant outflow includes a jet core which is highly relativistic after breaking out of the star, but becomes baryon-loaded after colliding with a massive outer shell, corresponding to mass from the stellar atmosphere of the progenitor star which became trapped in front of the jet core at breakout. The prompt emission produced before or during this collision would then have the signature of a high Lorentz factor jet, but the afterglow is produced by the amalgamated post-collision ejecta which has more inertia than the original highly relativistic jet core and thus has a delayed deceleration. This naturally explains the early light curve behavior discovered by Swift, including a steep decay and a long plateau, without invoking late-time energy injection from the central engine. The numerical simulation is performed continuously from engine to afterglow, covering a dynamic range of over ten orders of magnitude in radius. Light curves calculated from the numerical output demonstrate that this mechanism reproduces basic features seen in early afterglow data. Initial steep decays are produced by internal shocks, and the plateau corresponds to the coasting phase of the outflow. △ Less

Submitted 26 April, 2015; v1 submitted 30 July, 2014; originally announced July 2014.

Comments: ApJ Accepted

High-Frequency Voronoi Noise Reduced by Smoothed Mesh Motion

Authors: Paul C. Duffell, Andrew I. MacFadyen

Abstract: We describe a technique for improving the performance of hydrodynamics codes which employ a moving Voronoi mesh. Currently, such codes are susceptible to high-frequency noise produced by rapid adjustments in the grid topology on the smallest scales. The treatment for this grid noise is simple; instead of moving the mesh-generating marker points with the local fluid velocity, this velocity field is… ▽ More We describe a technique for improving the performance of hydrodynamics codes which employ a moving Voronoi mesh. Currently, such codes are susceptible to high-frequency noise produced by rapid adjustments in the grid topology on the smallest scales. The treatment for this grid noise is simple; instead of moving the mesh-generating marker points with the local fluid velocity, this velocity field is smoothed on small scales, so that neighboring marker points generally have similar velocities. We demonstrate significant improvement gained by this adjustment in several code tests relevant to the physics which moving-mesh codes are designed to capture. △ Less

Submitted 9 April, 2015; v1 submitted 27 July, 2014; originally announced July 2014.

Comments: MNRAS Accepted

The Migration of Gap-Opening Planets is not Locked to Viscous Disk Evolution

Authors: Paul C. Duffell, Zoltan Haiman, Andrew I. MacFadyen, Daniel J. D'Orazio, Brian D. Farris

Abstract: Most standard descriptions of Type II migration state that massive, gap-opening planets must migrate at the viscous drift rate. This is based on the idea that the disk is separated into an inner and outer region and gas is considered unable to cross the gap. In fact, gas easily crosses the gap on horseshoe orbits, nullifying this necessary premise which would set the migration rate. In this work,… ▽ More Most standard descriptions of Type II migration state that massive, gap-opening planets must migrate at the viscous drift rate. This is based on the idea that the disk is separated into an inner and outer region and gas is considered unable to cross the gap. In fact, gas easily crosses the gap on horseshoe orbits, nullifying this necessary premise which would set the migration rate. In this work, it is demonstrated using highly accurate numerical calculations that the actual migration rate is dependent on disk and planet parameters, and can be significantly larger or smaller than the viscous drift rate. In the limiting case of a disk much more massive than the secondary, the migration rate saturates to a constant which is sensitive to disk parameters and is not necessarily of order viscous rate. In the opposite limit of a low-mass disk, the migration rate decreases linearly with disk mass. Steady-state solutions in the low disk mass limit show no pile-up outside the secondary's orbit, and no corresponding drainage of the inner disk. △ Less

Submitted 1 August, 2014; v1 submitted 14 May, 2014; originally announced May 2014.

Comments: ApJ Accepted

A "Boosted Fireball" Model for Structured Relativistic Jets

Authors: Paul C. Duffell, Andrew I. MacFadyen

Abstract: We present a model for relativistic jets which generates a particular angular distribution of Lorentz factor and energy per solid angle. We consider a fireball with specific internal energy E/M launched with bulk Lorentz factor γ_B. This "boosted fireball" model is motivated by the phenomenology of collapsar jets, but is applicable to a wide variety of relativistic flows. In its center-of-momentum… ▽ More We present a model for relativistic jets which generates a particular angular distribution of Lorentz factor and energy per solid angle. We consider a fireball with specific internal energy E/M launched with bulk Lorentz factor γ_B. This "boosted fireball" model is motivated by the phenomenology of collapsar jets, but is applicable to a wide variety of relativistic flows. In its center-of-momentum frame the fireball expands isotropically, converting its internal energy into radially expanding flow with asymptotic Lorentz factor η_0 ~ E/M. In the lab frame the flow is beamed, expanding with Lorentz factor Γ= 2 η_0 γ_B in the direction of its initial bulk motion and with characteristic opening angle θ_0 ~ 1/γ_B. The flow is jet-like with Γθ_0 ~ 2 η_0 such that jets with Γ> 1/θ_0 are naturally produced. The choice η_0 ~ γ_B ~ 10 yields a jet with Γ~ 200 on-axis and angular structure characterized by opening angle θ_0 ~ 0.1 of relevance for cosmological GRBs, while γ_B >~ 1 may be relevant for low-luminosity GRBs. The model produces a family of outflows, of relevance for different relativistic phenomena with structures completely determined by η_0 and γ_B. We calculate the energy per unit solid angle for the model and use it to compute light curves for comparison with the widely used top-hat model. The jet break in the boosted fireball light curve is greatly subdued when compared to the top-hat model because the edge of the jet is smoother than for a top-hat. This may explain missing jet breaks in afterglow light curves. △ Less

Submitted 2 October, 2013; v1 submitted 7 August, 2013; originally announced August 2013.

Comments: ApJ Accepted

Rayleigh-Taylor Instability in a Relativistic Fireball on a Moving Computational Grid

Authors: Paul C. Duffell, Andrew I. MacFadyen

Abstract: We numerically calculate the growth and saturation of the Rayleigh-Taylor instability caused by the deceleration of relativistic outflows with Lorentz factor Γ = 10, 30, and 100. The instability generates turbulence whose scale exhibits strong dependence on Lorentz factor, as only modes with angular size smaller than 1/Γ can grow. We develop a simple diagnostic to measure the kinetic energy in tur… ▽ More We numerically calculate the growth and saturation of the Rayleigh-Taylor instability caused by the deceleration of relativistic outflows with Lorentz factor Γ = 10, 30, and 100. The instability generates turbulence whose scale exhibits strong dependence on Lorentz factor, as only modes with angular size smaller than 1/Γ can grow. We develop a simple diagnostic to measure the kinetic energy in turbulent fluctuations, and calculate a ratio of turbulent kinetic energy to thermal energy of .03 in the region affected by the instability. Although our numerical calculation does not include magnetic fields, we argue that small scale turbulent dynamo amplifies magnetic fields to nearly this same fraction, giving a ratio of magnetic to thermal energy of ~ .01, to within a factor of two. The instability completely disrupts the contact discontinuity between the ejecta and the swept up circumburst medium. The reverse shock is stable, but is impacted by the Rayleigh-Taylor instability, which strengthens the reverse shock and pushes it away from the forward shock. The forward shock front is unaffected by the instability, but Rayleigh-Taylor fingers can penetrate of order 10% of the way into the energetic region behind the shock during the two-shock phase of the explosion. We calculate afterglow emission from the explosion and find the reverse shock emission peaks at a later time due to its reduced Lorentz factor and modified density and pressure at the shock front. These calculations are performed using a novel numerical technique that includes a moving computational grid. The moving grid is essential as it maintains contact discontinuities to high precision and can easily evolve flows with extremely large Lorentz factors. △ Less

Submitted 2 October, 2013; v1 submitted 28 February, 2013; originally announced February 2013.

Comments: ApJ Accepted

Gap Opening by Extremely Low Mass Planets in a Viscous Disk

Authors: Paul C. Duffell, Andrew I MacFadyen

Abstract: By numerically integrating the compressible Navier-Stokes equations in two dimensions, we calculate the criterion for gap formation by a very low mass (q ~10^{-4}) protoplanet on a fixed orbit in a thin viscous disk. In contrast with some previously proposed gap-opening criteria, we find that a planet can open a gap even if the Hill radius is smaller than the disk scale height. Moreover, in the lo… ▽ More By numerically integrating the compressible Navier-Stokes equations in two dimensions, we calculate the criterion for gap formation by a very low mass (q ~10^{-4}) protoplanet on a fixed orbit in a thin viscous disk. In contrast with some previously proposed gap-opening criteria, we find that a planet can open a gap even if the Hill radius is smaller than the disk scale height. Moreover, in the low-viscosity limit, we find no minimum mass necessary to open a gap for a planet held on a fixed orbit. In particular, a Neptune-mass planet will open a gap in a minimum mass solar nebula with suitably low viscosity (α<10^{-4}). We find that the mass threshold scales as the square root of viscosity in the low mass regime. This is because the gap width for critical planet masses in this regime is a fixed multiple of the scale height, not of the Hill radius of the planet. △ Less

Submitted 5 April, 2013; v1 submitted 7 February, 2013; originally announced February 2013.

Comments: ApJ accepted

Global Calculations of Density Waves and Gap Formation in Protoplanetary Disks using a Moving Mesh

Authors: Paul C. Duffell, Andrew I. MacFadyen

Abstract: We calculate the global quasi-steady state of a thin disk perturbed by a low-mass protoplanet orbiting at a fixed radius using extremely high-resolution numerical integrations of Euler's equations in two dimensions. The calculations are carried out using a moving computational domain, which greatly reduces advection errors and allows for much longer time-steps than a fixed grid. We calculate the a… ▽ More We calculate the global quasi-steady state of a thin disk perturbed by a low-mass protoplanet orbiting at a fixed radius using extremely high-resolution numerical integrations of Euler's equations in two dimensions. The calculations are carried out using a moving computational domain, which greatly reduces advection errors and allows for much longer time-steps than a fixed grid. We calculate the angular momentum flux and the torque density as a function of radius and compare them with analytical predictions. We discuss the quasi-steady state after 100 orbits and the prospects for gap formation by low mass planets. △ Less

Submitted 4 June, 2012; v1 submitted 25 February, 2012; originally announced February 2012.

Comments: Accepted to ApJ

TESS: A Relativistic Hydrodynamics Code on a Moving Voronoi Mesh

Authors: Paul C. Duffell, Andrew I. MacFadyen

Abstract: We have generalized a method for the numerical solution of hyperbolic systems of equations using a dynamic Voronoi tessellation of the computational domain. The Voronoi tessellation is used to generate moving computational meshes for the solution of multi-dimensional systems of conservation laws in finite-volume form. The mesh generating points are free to move with arbitrary velocity, with the ch… ▽ More We have generalized a method for the numerical solution of hyperbolic systems of equations using a dynamic Voronoi tessellation of the computational domain. The Voronoi tessellation is used to generate moving computational meshes for the solution of multi-dimensional systems of conservation laws in finite-volume form. The mesh generating points are free to move with arbitrary velocity, with the choice of zero velocity resulting in an Eulerian formulation. Moving the points at the local fluid velocity makes the formulation effectively Lagrangian. We have written the TESS code to solve the equations of compressible hydrodynamics and magnetohydrodynamics for both relativistic and non-relativistic fluids on a dynamic Voronoi mesh. When run in Lagrangian mode, TESS is significantly less diffusive than fixed mesh codes and thus preserves contact discontinuities to high precision while also accurately capturing strong shock waves. TESS is written for Cartesian, spherical and cylindrical coordinates and is modular so that auxilliary physics solvers are readily integrated into the TESS framework and so that the TESS framework can be readily adapted to solve general systems of equations. We present results from a series of test problems to demonstrate the performance of TESS and to highlight some of the advantages of the dynamic tessellation method for solving challenging problems in astrophysical fluid dynamics. △ Less

Submitted 8 November, 2011; v1 submitted 18 April, 2011; originally announced April 2011.

Comments: ApJS, 197, 15