Observational evidence suggests a link between long duration gamma ray bursts (LGRBs) and Type Ic supernovae. Here, we propose a potential mechanism for Type Ic supernovae in LGRB progenitors powered solely by accretion energy. We present spherically-symmetric hydrodynamic simulations of the long-term accretion of a rotating gamma-ray burst progenitor star, a "collapsar," onto the central compact object, which we take to be a black hole. The simulations were carried out with the adaptive mesh refinement code FLASH in one spatial dimension and with rotation, an explicit shear viscosity, and convection in the mixing length theory approximation. Once the accretion flow becomes rotationally supported outside of the black hole, an accretion shock forms and traverses the stellar envelope. Energy is carried from the central geometrically thick accretion disk to the stellar envelope by convection. Energy losses through neutrino emission and nuclear photodisintegration are calculated but do not seem important following the rapid early drop of the accretion rate following circularization. We find that the shock velocity, energy, and unbound mass are sensitive to convective efficiency, effective viscosity, and initial stellar angular momentum. Our simulations show that given the appropriate combinations of stellar and physical parameters, explosions with energies ~5x10^50 ergs, velocities 3000 km/s, and unbound material masses >6 solar masses are possible in a rapidly rotating 16 solar mass main sequence progenitor star. Further work is needed to constrain the values of these parameters, to identify the likely outcomes in more plausible and massive LRGB progenitors, and to explore nucleosynthetic implications.

The association of long-duration gamma-ray bursts (LGRBs) with Type Ic supernovae presents a challenge to supernova explosion models. In the collapsar model for LGRBs, gamma rays are produced in an ultrarelativistic jet launching from the magnetosphere of the black hole that forms in the aftermath of the collapse of a rotating progenitor star. The jet is collimated along the star's rotation axis, but the concomitant luminous supernova should be relatively--though certainly not entirely--spherical, and should synthesize a substantial mass of 56Ni. Our goal is to provide a qualitative assessment of the possibility that accretion of the progenitor envelope onto the black hole, which powers the LGRB, could also deposit sufficient energy and nickel mass in the envelope to produce a luminous supernova. For this, the energy dissipated near the black hole during accretion must be transported outward, where it can drive a supernova-like shockwave. Here we suggest that the energy is transported by convection and develop an analytical toy model, relying on global mass and energy conservation, for the dynamics of stellar collapse. The model suggests that a ~10,000 km/s shock can be driven into the envelope and that ~10^51 erg explosions are possible. The efficiency with which the accretion energy is being transferred to the envelope is governed by the competition of advection and convection at distances ~100-1,000 km from the black hole and is sensitive to the values of the convective mixing length, the magnitude of the effective viscous stress, and the specific angular momentum of the infalling envelope. Substantial masses of 56Ni may be synthesized in the convective accretion flow over the course of tens of seconds from the initial circularization of the infalling envelope around the black hole. The synthesized nickel is convectively mixed with a much larger mass of unburned ejecta.

We present axisymmetric hydrodynamical simulations of the long-term accretion of a rotating GRB progenitor star, a "collapsar," onto the central compact object. The simulations were carried out with the adaptive mesh refinement code FLASH in two spatial dimensions and with an explicit shear viscosity. The evolution of the central accretion rate exhibits phases reminiscent of the long GRB gamma-ray and X-ray light curve, which lends support to the proposal that the luminosity is modulated by the central accretion rate. After a few tens of seconds, an accretion shock sweeps outward through the star. The formation and outward expansion of the accretion shock is accompanied with a sudden and rapid power-law decline in the central accretion rate Mdot ~ t^{-2.8}, which resembles the L_X ~ t^{-3} decline observed in the X-ray light curves. The collapsed, shock-heated stellar envelope settles into a thick, low-mass equatorial disk embedded within a massive, pressure-supported atmosphere. After a few hundred seconds, the inflow of low-angular-momentum material in the axial funnel reverses into an outflow from the surface of the thick disk. Meanwhile, the rapid decline of the accretion rate slows down, or even settles a in steady state with Mdot ~ 5x10^{-5} Msun/s, which resembles the "plateau" phase in the X-ray light curve. While the duration of the "prompt" phase depends on the resolution in our simulations, we provide an analytical model taking into account neutrino losses that estimates the duration to be ~20 s. The model suggests that the steep decline in GRB X-ray light curves is triggered by the circularization of the infalling stellar envelope at radii where the virial temperature is below ~10^{10} K, such that neutrino cooling shuts off and an outward expansion of the accretion shock becomes imminent.

In recent work we presented the first results of global general relativistic magnetohydrodynamic (GRMHD) simulations of tilted (or misaligned) accretion disks around rotating black holes. The simulated tilted disks showed dramatic differences from comparable untilted disks, such as asymmetrical accretion onto the hole through opposing "plunging streams" and global precession of the disk powered by a torque provided by the black hole. However, those simulations used a traditional spherical-polar grid that was purposefully underresolved along the pole, which prevented us from assessing the behavior of any jets that may have been associated with the tilted disks. To address this shortcoming we have added a block-structured "cubed-sphere" grid option to the Cosmos++ GRMHD code, which will allow us to simultaneously resolve the disk and polar regions. Here we present our implementation of this grid and the results of a small suite of validation tests intended to demonstrate that the new grid performs as expected. The most important test in this work is a comparison of identical tilted disks, one evolved using our spherical-polar grid and the other with the cubed-sphere grid. We also demonstrate an interesting dependence of the early-time evolution of our disks on their orientation with respect to the grid alignment. This dependence arises from the differing treatment of current sheets within the disks, especially whether they are aligned with symmetry planes of the grid or not.

We present two-dimensional spectroscopy covering the rest-frame wavelengths of strong optical emission lines in six luminous submillimetre galaxies at z=1.3--2.5. Together with HST ACS and NICMOS imaging we map the dynamics and morphologies of these systems on scales from 4--11kpc. Four of the systems show multiple components in their spatially-resolved spectra with average velocity offsets of ~180km/s across 8kpc in projection. From the ensemble properties from our survey and the literature, we estimate the typical dynamical masses of bright sub-mm galaxies as (5+/-3) x10^{11}Mo. This is similar to recent estimates of their stellar masses -- suggesting that the dynamics of the central regions of these galaxies are baryon dominated, with a substantial fraction of those baryons in stars by the epoch of observation. Combining our dynamical mass estimates with stellar luminosities for this population we investigate whether submillimetre galaxies can evolve onto the Faber-Jackson relation for local ellipticals. Adopting a typical lifetime of tau_{burst}~300Myr for the submillimetre-luminous phase we find that the stellar populations of sub-mm galaxies should fade to place them on the Faber-Jackson relation, at M_K ~ -25.1. Furthermore, using the same starburst lifetime we correct the observed space density of submillimetre galaxies for the duty cycle to derive a volume density of the progenitors of ~1x10^{-4} Mpc^{-3}. This is consistent with the space density of local luminous early-type galaxies with M_K ~ -25.1, indicating that submillimetre galaxies can evolve onto the scaling relations observed for local early-type galaxies, and the observed population at z~2 is then sufficient to account for the formation of the whole population of >~3 L^{*} ellipticals seen at z~0.