BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrew MacFadyen

Andrew MacFadyen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.COphysics.comp-phastro-ph.GAphysics.flu-dynastro-ph.SRastro-ph.EPgr-qcphysics.plasm-ph

Catalog footprint

What is connected

22works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

22 published item(s)

preprint2022arXiv

Ellipsars: Ring-like explosions from flattened stars

The stellar cataclysms producing astronomical transients have long been modeled as either a point-like explosion or jet-like engine ignited at the center of a spherically symmetric star. However, many stars are observed, or are expected on theoretical grounds, not to be precisely spherically symmetric, but rather to have a slightly flattened geometry similar to that of an oblate spheroid. Here we present axisymmetric two-dimensional hydrodynamical simulations of the dynamics of point-like explosions initiated at the center of an aspherical massive star with a range of oblateness. We refer to these exploding aspherical stars as "ellipsars" in reference to the elliptical shape of the iso-density contours of their progenitors in the two-dimensional axisymmetric case. We find that ellipsars are capable of accelerating expanding rings of relativistic ejecta which may lead to the production of astronomical transients including low-luminosity GRBs, relativistic supernovae, and Fast Blue Optical Transients (FBOTs.)

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Equilibrium eccentricity of accreting binaries

Using high-resolution hydrodynamics simulations, we show that equal-mass binaries accreting from a circumbinary disk evolve toward an orbital eccentricity of $e \simeq 0.45$, unless they are initialized on a nearly circular orbit with $e \lesssim 0.08$, in which case they further circularize. The implied bi-modal eccentricity distribution resembles that seen in post-AGB stellar binaries. Large accretion spikes around periapse impart a tell-tale, quasi-periodic, bursty signature on the light curves of eccentric binaries. We predict that intermediate-mass and massive black hole binaries at $z \lesssim 10$ entering the \emph{LISA} band will have measurable eccentricities in the range $e \simeq 10^{-3} - 10^{-2}$, if they have experienced a gas-driven phase. On the other hand, GW190521 would have entered the LIGO/Virgo band with undetectable eccentricity $\sim 10^{-6}$ if it had been driven into the gravitational wave regime by a gas disk.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Gas-driven inspiral of binaries in thin accretion disks

Numerical studies of gas accretion onto supermassive black hole binaries (SMBHBs) have generally been limited to conditions where the circumbinary disk (CBD) is 10-100 times thicker than expected for disks in active galactic nuclei (AGN). This discrepancy arises from technical limitations, and also from publication bias toward replicating fiducial numerical models. Here we present the first systematic study of how the binary's orbital evolution varies with disk scale height. We report three key results: (1) Binary orbital evolution switches from outspiralling for warm disks (aspect ratio ~0.1), to inspiralling for more realistic cooler, thinner disks at a critical aspect ratio ~0.04, corresponding to orbital Mach number ~25. (2) The net torque on the binary arises from a competition between positive torque from gas orbiting close to the black holes, and negative torque from the inner edge of the CBD, which is denser for thinner disks. This leads to increasingly negative net torques on the binary for increasingly thin disks. (3) The accretion rate is modestly suppressed with increasing Mach number. We discuss how our results may influence modeling of the nano-Hz gravitational wave background, as well as estimates of the LISA merger event rate.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Missing-Link in the Supernova-GRB Connection: The Case of SN 2012ap

Gamma Ray Bursts (GRBs) are characterized by ultra-relativistic outflows, while supernovae are generally characterized by non-relativistic ejecta. GRB afterglows decelerate rapidly usually within days, because their low-mass ejecta rapidly sweep up a comparatively larger mass of circumstellar material. However supernovae, with heavy ejecta, can be in nearly free expansion for centuries. Supernovae were thought to have non-relativistic outflows except for few relativistic ones accompanied by GRBs. This clear division was blurred by SN 2009bb, the first supernova with a relativistic outflow without an observed GRB. Yet the ejecta from SN 2009bb was baryon loaded, and in nearly-free expansion for a year, unlike GRBs. We report the first supernova discovered without a GRB, but with rapidly decelerating mildly relativistic ejecta, SN 2012ap. We discovered a bright and rapidly evolving radio counterpart driven by the circumstellar interaction of the relativistic ejecta. However, we did not find any coincident GRB with an isotropic fluence of more than a sixth of the fluence from GRB 980425. This shows for the first time that central engines in type Ic supernovae, even without an observed GRB, can produce both relativistic and rapidly decelerating outflows like GRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Analysis of Chandra Deep Follow-up GRBs: Implications for Off-Axis Jets

We present a sample of 27 GRBs with detailed Swift light curves supplemented by late time Chandra observations. To answer the missing jet-break problem in general, we develop a numerical simulation based model which can be directly fit to the data using Monte Carlo methods. Our numerical model takes into account all the factors that can shape a jet-break: (i) lateral expansion (ii) edge effects and (iii) off-axis effects. Our results provide improved fits to the light curves and constraints on physical parameters. More importantly, our results suggest that off-axis effects are important and must be included in interpretations of GRB jet breaks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gamma Ray Bursts Are Observed Off-Axis

We use the ScaleFit package to perform Markov chain Monte Carlo light curve fitting on a large sample of Swift-XRT gamma-ray burst afterglows. The ScaleFit model uses scaling relations in the hydrodynamic and radiation equations to compute synthetic light curves directly from a set of high resolution two-dimensional relativistic blast wave simulations. The data sample consists of all Swift-XRT afterglows from 2005 to 2012 with sufficient coverage and a known redshift, 188 bursts in total. We find the jet half opening angle varies widely but is commonly less than 0.1 radians. The distribution of the electron spectral index is also broad, with a median at 2.23. This approach allows, for the first time, for the off-axis observer angle to be inferred directly from the light curve. We find the observer angle to have a median value of 0.57 of the jet opening angle over our sample, which has profound consequences for the predicted rate of observed jet breaks and effects the beaming corrected total energies of gamma-ray bursts.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Accretion into the Central Cavity of a Circumbinary Disk

A near-equal mass binary black hole can clear a central cavity in a circumbinary accretion disk; however, previous works have revealed accretion streams entering this cavity. Here we use 2D hydrodynamical simulations to study the accretion streams and their periodic behavior. In particular, we perform a suite of simulations, covering different binary mass ratios $q=M_2/M_1$ in the range $0.01 \leq q \leq 1$. In each case, we follow the system for several thousand binary orbits, until it relaxes to a stable accretion pattern. We find the following results: (i) while the binary is efficient in maintaining a low-density cavity, the time-averaged mass accretion rate into the cavity, through narrow coherent accretion streams, is suppressed by at most a factor of $\sim 5$, compared to a disk with a single BH with the same mass; (ii) the largest suppression occurs for $q\approx 0.05$; binaries whose mass ratios are either lower or higher both suppress accretion less significantly; (iii) for $q \gtrsim 0.05$, the accretion rate is strongly modulated by the binary, and depending on the precise value of $q$, the power spectrum of the accretion rate shows either one, two, or three distinct periods; and (v) for $q \lesssim 0.05$, the accretion rate becomes steady, with no time-variations. Most binaries produced in galactic mergers are expected to have $q\gtrsim 0.05$. If the luminosity of these binaries tracks their accretion rate, then a periodogram of their light-curve could help in their identification, and to constrain their mass ratio and disk properties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fitting Afterglows With Multi-Dimensional Simulations

We present preliminary data fit results of synthetic light curves computed from numerical afterglow blast wave simulations. Our technique uses Markov chain Monte Carlo (MCMC) in a new data analysis tool, ScaleFit. Scaling relations in both the hydrodynamics and radiation equations allow light curves to be parameterized by a small set of scale-invariant characteristic quantities. These quantities have been calculated and tabulated from high resolution two-dimensional hydrodynamic simulations. Producing a light curve from the characteristics takes only a millisecond, allowing for the use of MCMC data fitting techniques which can require millions of iterations. ScaleFit is a portable, lightweight, python package which performs this analysis on afterglow light curves. Using the set of Swift-XRT light curves from 2011 & 2012 with known redshifts, we find ScaleFit can measure the jet opening angle, observer angle, and spectral index of most afterglows. Globally we find gamma-ray burst afterglows tend to be observed off axis, at a significant fraction of the jet opening angle.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gamma-ray burst afterglow light curves from a Lorentz-boosted simulation frame and the shape of the jet break

GRB afterglow jets have been notoriously difficult to resolve numerically using 2D hydrodynamical simulations due to high outflow Lorentz factors. By performing simulations in a boosted frame, it is possible to calculate light curves from numerically computed flows in sufficient detail to accurately quantify the shape of the jet break and the post-break steepening of the light curve. We study jet breaks for jets with opening angles of theta_0 = 0.05, 0.1 and 0.2 radians decelerating in a constant density medium, observed at angles (theta_obs) ranging from on-axis to the jet edge. We present a single set of scale-invariant functions describing the time evolution of synchrotron spectral break frequencies and peak flux, depending only on theta_0 and theta_obs, sufficient to reconstruct light curves for arbitrary explosion energy, circumburst density and synchrotron slope p. We compare our light curves to fit functions using sharp breaks and smooth transitions and confirm that the measured jet break time can be delayed due to theta_obs. We find that the difference in temporal indices across the jet break is larger than theoretically anticipated leading to post-break slopes of roughly 0.25 - 1.3 p, sufficiently steep to not easily reconcile with post-break slopes measured for the Swift sample, suggesting that most Swift GRBs do not explode in homogeneous media or that jet breaks are hidden by additional physics such as prolonged energy injection or viewing angle effects. Light curves and spectral break and peak flux evolution functions will be made publicly available on-line at http://cosmo.nyu.edu/afterglowlibrary .

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

GRB Afterglow Blast Wave Encountering Sudden Circumburst Density Change Produces No Flares

Afterglows of gamma-ray bursts are observed to produce light curves with the flux following power law evolution in time. However, recent observations reveal bright flares at times on the order of minutes to days. One proposed explanation for these flares is the interaction of a relativistic blast wave with a circumburst density transition. In this paper, we model this type of interaction computationally in one and two dimensions, using a relativistic hydrodynamics code with adaptive mesh refinement called ram, and analytically in one dimension. We simulate a blast wave traveling in a stellar wind environment that encounters a sudden change in density, followed by a homogeneous medium, and compute the observed radiation using a synchrotron model. We show that flares are not observable for an encounter with a sudden density increase, such as a wind termination shock, nor for an encounter with a sudden density decrease. Furthermore, by extending our analysis to two dimensions, we are able to resolve the spreading, collimation, and edge effects of the blast wave as it encounters the change in circumburst medium. In all cases considered in this paper, we find that a flare will not be observed for any of the density changes studied.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

No flares from Gamma-Ray Burst afterglow blast waves encountering sudden circumburst density change

Afterglows of gamma-ray bursts are observed to produce light curves with the flux following power law evolution in time. However, recent observations reveal bright flares at times on the order of minutes to days. One proposed explanation for these flares is the interaction of a relativistic blast wave with a circumburst density transition. In this paper, we model this type of interaction computationally in one and two dimensions, using a relativistic hydrodynamics code with adaptive mesh refinement called RAM, and analytically in one dimension. We simulate a blast wave traveling in a stellar wind environment that encounters a sudden change in density, followed by a homogeneous medium, and compute the observed radiation using a synchrotron model. We show that flares are not observable for an encounter with a sudden density increase, such as a wind termination shock, nor for an encounter with a sudden density decrease. Furthermore, by extending our analysis to two dimensions, we are able to resolve the spreading, collimation, and edge effects of the blast wave as it encounters the change in circumburst medium. In all cases considered in this paper, we find that a flare will not be observed for any of the density changes studied.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spectral and Intermittency Properties of Relativistic Turbulence

High-resolution numerical simulations are utilized to examine isotropic turbulence in a compressible fluid when long wavelength velocity fluctuations approach light speed. Spectral analysis reveals an inertial sub-range of relativistic motions with a broadly 5/3 index. The use of generalized Lorentz-covariant structure functions based on the four-velocity is proposed. These structure functions extend the She-Leveque model for intermittency into the relativistic regime.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Transient Gravitational-Wave Sky

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical Simulations of Driven Relativistic MHD Turbulence

A wide variety of astrophysical phenomena involve the flow of turbulent magnetized gas with relativistic velocity or energy density. Examples include gamma-ray bursts, active galactic nuclei, pulsars, magnetars, micro-quasars, merging neutron stars, X-ray binaries, some supernovae, and the early universe. In order to elucidate the basic properties of the relativistic magnetohydrodynamical (RMHD) turbulence present in these systems, we present results from numerical simulations of fully developed driven turbulence in a relativistically warm, weakly magnetized and mildly compressible ideal fluid. We have evolved the RMHD equations for many dynamical times on a uniform grid with 1024^3 zones using a high order Godunov code. We observe the growth of magnetic energy from a seed field through saturation at about 1% of the total fluid energy. We compute the power spectrum of velocity and density-weighted velocity and conclude that the inertial scaling is consistent with a slope of -5/3. We compute the longitudinal and transverse velocity structure functions of order p up to 11, and discuss their possible deviation from the expected scaling for non-relativistic media. We also compute the scale-dependent distortion of coherent velocity structures with respect to the local magnetic field, finding a weaker scale dependence than is expected for incompressible non-relativistic flows with a strong mean field.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical Simulations of Driven Supersonic Relativistic MHD Turbulence

Models for GRB outflows invoke turbulence in relativistically hot magnetized fluids. In order to investigate these conditions we have performed high-resolution three-dimensional numerical simulations of relativistic magneto-hydrodynamical (RMHD) turbulence. We find that magnetic energy is amplified to several percent of the total energy density by turbulent twisting and folding of magnetic field lines. Values of epsilon_B near 1% are thus naturally expected. We study the dependence of saturated magnetic field energy fraction as a function of Mach number and relativistic temperature. We then present power spectra of the turbulent kinetic and magnetic energies. We also present solenoidal (curl-like) and dilatational (divergence-like) power spectra of kinetic energy. We propose that relativistic effects introduce novel couplings between these spectral components. The case we explore in most detail is for equal amounts of thermal and rest mass energy, corresponding to conditions after collisions of shells with relative Lorentz factors of several. These conditions are relevant in models for internal shocks, for the late afterglow phase, for cocoon material along the edge of a relativistic jet as it propagates through a star, as well neutron stars merging with each other and with black hole companions. We find that relativistic turbulence decays extremely quickly, on a sound crossing time of an eddy. Models invoking sustained relativistic turbulence to explain variability in GRB prompt emission are thus strongly disfavored unless a persistant driving of the turbulence is maintained for the duration of the prompt emission.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Slightly Two or Three Dimensional Self-Similar Solutions

Self similarity allows for analytic or semi-analytic solutions to many hydrodynamics problems. Most of these solutions are one dimensional. Using linear perturbation theory, expanded around such a one-dimensional solution, we find self-similar hydrodynamic solutions that are two- or three-dimensional. Since the deviation from a one-dimensional solution is small, we call these slightly two-dimensional and slightly three-dimensional self-similar solutions, respectively. As an example, we treat strong spherical explosions of the second type. A strong explosion propagates into an ideal gas with negligible temperature and density profile of the form rho(r,theta,phi)=r^{-omega}[1+sigma*F(theta,phi)], where omega>3 and sigma << 1. Analytical solutions are obtained by expanding the arbitrary function F(theta,phi) in spherical harmonics. We compare our results with two dimensional numerical simulations, and find good agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Flare-less long Gamma-ray Bursts and the properties of their massive star progenitors

While there is mounting evidence that long Gamma-Ray Bursts (GRBs) are associated with the collapse of massive stars, the detailed structure of their pre-supernova stage is still debatable. Particularly uncertain is the degree of mixing among shells of different composition, and hence the role of magnetic torques and convection in transporting angular momentum. Here we show that early-time afterglow observations with the Swift satellite place constraints on the allowed GRB pre-supernova models. In particular, they argue against pre-supernova models in which different elemental shells are unmixed. These types of models would produce energy injections into the GRB engine on timescales between several hundreds of seconds to a few hours. Flaring activity has {\em not} been observed in a large fraction of well-monitored long GRBs. Therefore, if the progenitors of long GRBs have common properties, then the lack of flares indicates that the massive stars which produce GRBs are mostly well mixed, as expected in low-metallicity, rapidly rotating massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Off-Axis Gamma-Ray Burst Afterglow Modeling Based On A Two-Dimensional Axisymmetric Hydrodynamics Simulation

Starting as highly relativistic collimated jets, gamma-ray burst outflows gradually decelerate and become non-relativistic spherical blast waves. Although detailed analytical solutions describing the afterglow emission received by an on-axis observer during both the early and late phases of the outflow evolution exist, a calculation of the received flux during the intermediate phase and for an off-axis observer requires either a more simplified analytical model or direct numerical simulations of the outflow dynamics. In this paper we present light curves for off-axis observers covering the long-term evolution of the blast wave calculated from a high resolution two-dimensional relativistic hydrodynamics simulation using a synchrotron radiation model. We compare our results to earlier analytical work and calculate the consequence of the observer angle with respect to the jet axis both for the detection of orphan afterglows and for jet break fits to the observational data. We find that observable jet breaks can be delayed for up to several weeks for off-axis observers, potentially leading to overestimation of the beaming corrected total energy. When using our off-axis light curves to create synthetic Swift X-ray data, we find that jet breaks are likely to remain hidden in the data. We also confirm earlier results in the literature finding that only a very small number of local Type Ibc supernovae can harbor an orphan afterglow.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hydrodynamical Response of a Circumbinary Gas Disk to Black Hole Recoil and Mass Loss

Finding electromagnetic (EM) counterparts of future gravitational wave (GW) sources would bring rich scientific benefits. A promising possibility, in the case of the coalescence of a super-massive black hole binary (SMBHB), is that prompt emission from merger-induced disturbances in a supersonic circumbinary disk may be detectable. We follow the post-merger evolution of a thin, zero-viscosity circumbinary gas disk with two-dimensional simulations, using the hydrodynamic code FLASH. We analyze perturbations arising from the 530 km/s recoil of a 10^6 M_sun binary, oriented in the plane of the disk, assuming either an adiabatic or a pseudo-isothermal equation of state for the gas. We find that a single-armed spiral shock wave forms and propagates outward, sweeping up about 20% of the mass of the disk. The morphology and evolution of the perturbations agrees well with those of caustics predicted to occur in a collisionless disk. Assuming that the disk radiates nearly instantaneously to maintain a constant temperature, we estimate the amount of dissipation and corresponding post-merger light-curve. The luminosity rises steadily on the time-scale of months, and reaches few times 10^{43} erg/s, corresponding to about 10% of the Eddington luminosity of the central SMBHB. We also analyze the case in which gravitational wave emission results in a 5% mass loss in the merger remnant. The mass-loss reduces the shock overdensities and the overall luminosity of the disk by 15-20%, without any other major effects on the spiral shock pattern.

0 citations0 reviewsNo code linkedOpen access