We performed 2D, axisymmetric, MHD simulations with Cosmos++ in order to examine the growth of the magnetorotational instability (MRI) in core--collapse supernovae. We have initialized a non--rotating 15 solar mass progenitor, infused with differential rotation and poloidal magnetic fields. The collapse of the iron core is simulated with the Shen EOS, and the parametric Ye and entropy evolution. The wavelength of the unstable mode in the post--collapse environment is expected to be only ~ 200 m. In order to achieve the fine spatial resolution requirement, we employed remapping technique after the iron core has collapsed and bounced. The MRI unstable region appears near the equator and angular momentum and entropy are transported outward. Higher resolution remap run display more vigorous overturns and stronger transport of angular momentum and entropy. Our results are in agreement with the earlier work by Akiyama et al. (2003) and Obergaulinger et al. (2009).

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.

This paper presents a continuation of our efforts to numerically study accretion disks that are misaligned (tilted) with respect to the rotation axis of a Kerr black hole. Here we present results of a global numerical simulation which fully incorporates the effects of the black hole spacetime as well as magnetorotational turbulence that is the primary source of angular momentum transport in the flow. This simulation shows dramatic differences from comparable simulations of untilted disks. Accretion onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is due to the aspherical nature of the gravitational spacetime around the rotating black hole. These plunging streams start from a larger radius than would be expected for an untilted disk. In this regard the tilted black hole effectively acts like an untilted black hole of lesser spin. Throughout the duration of the simulation, the main body of the disk remains tilted with respect to the symmetry plane of the black hole; thus there is no indication of a Bardeen-Petterson effect in the disk at large. The torque of the black hole instead principally causes a global precession of the main disk body. In this simulation the precession has a frequency of $3 (M_\odot/M)$ Hz, a value consistent with many observed low-frequency quasi-periodic oscillations. However, this value is strongly dependent on the size of the disk, so this frequency may be expected to vary over a large range.

We posit that accreting compact objects, including stellar mass black holes and neutron stars as well as supermassive black holes, may undergo extended periods of accretion during which the angular momentum of the disk at large scales is misaligned with that of the compact object. In such a scenario, Lense-Thirring precession caused by the rotating compact object can dramatically affect the disk. In this presentation we describe results from a three-dimensional relativistic magnetohydrodynamic simulation of an MRI turbulent disk accreting onto a tilted rapidly rotating black hole. For this case, the disk does not achieve the commonly described Bardeen-Petterson configuration; rather, it remains nearly planar, undergoing a slow global precession. Accretion from the disk onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is a consequence of the non-sphericity of the gravitational spacetime of the black hole.

We use the relativistic hydrodynamics code Cosmos++ to model the evolution of the radio nebula triggered by the Dec. 27, 2004 giant flare event of soft gamma repeater 1806-20. We primarily focus on the rebrightening and centroid motion occurring subsequent to day 20 following the flare event. We model this period as a mildly relativistic (gamma ~ 1.07 - 1.67) jetted outflow expanding into the interstellar medium (ISM). We demonstrate that a jet with total energy ~ 10^46 ergs confined to a half opening angle ~ 20 degrees fits the key observables of this event, e.g. the flux lightcurve, emission map centroid position, and aspect ratio. In particular, we find excellent agreement with observations if the rebrightening is due to the jet, moving at 0.5c and inclined ~ 0 - 40 degrees toward the observer, colliding with a density discontinuity in the ISM at a radius of several 10^16 cm. We also find that a jet with a higher velocity, >~ 0.7c, and larger inclination, >~ 70 degrees, moving into a uniform ISM can fit the observations in general, but tends to miss the details of rebrightening. The latter, uniform ISM model predicts an ISM density more than 100 times lower than that of the former model, and thus suggests an independent test which might discriminate between the two. One of the strongest constraints of both models is that the data seems to require a non-uniform jet in order to be well fit.

A new code and methodology are introduced for solving the general relativistic magnetohydrodynamic (GRMHD) equations in fixed background spacetimes using time-explicit, finite-volume discretization. The code has options for solving the GRMHD equations using traditional artificial-viscosity (AV) or non-oscillatory central difference (NOCD) methods, or a new extended AV (eAV) scheme using artificial-viscosity together with a dual energy-flux-conserving formulation. The dual energy approach allows for accurate modeling of highly relativistic flows at boost factors well beyond what has been achieved to date by standard artificial viscosity methods. It provides the benefit of Godunov methods in capturing high Lorentz boosted flows but without complicated Riemann solvers, and the advantages of traditional artificial viscosity methods in their speed and flexibility. Additionally, the GRMHD equations are solved on an unstructured grid that supports local adaptive mesh refinement using a fully threaded oct-tree (in three dimensions) network to traverse the grid hierarchy across levels and immediate neighbors. A number of tests are presented to demonstrate robustness of the numerical algorithms and adaptive mesh framework over a wide spectrum of problems, boosts, and astrophysical applications, including relativistic shock tubes, shock collisions, magnetosonic shocks, Alfven wave propagation, blast waves, magnetized Bondi flow, and the magneto-rotational instability in Kerr black hole spacetimes.

We present a numerical investigation of the contribution of the presupernova ejecta of Wolf-Rayet stars to the environment surrounding gamma-ray bursts (GRBs), and describe how this external matter can affect the observable afterglow characteristics. An implicit hydrodynamic calculation for massive stellar evolution is used here to provide the inner boundary conditions for an explicit hydrodynamical code to model the circumstellar gas dynamics. The resulting properties of the circumstellar medium are then used to calculate the deceleration of a relativistic, gas-dynamic jet and the corresponding afterglow light curve produced as the shock wave propagates through the shocked-wind medium. We find that variations in the stellar wind drive instabilities that may produce radial filaments in the shocked-wind region. These comet-like tails of clumps could give rise to strong temporal variations in the early afterglow lightcurve. Afterglows may be expected to differ widely among themselves, depending on the angular anisotropy of the jet and the properties of the stellar progenitor; a wide diversity of behaviors may be the rule, rather than the exception.

Funded by $1.2M in grants and donations, we are now building PROMPT at CTIO. When completed in late 2005, PROMPT will consist of six 0.41-meter diameter Ritchey-Chretien telescopes on rapidly slewing mounts that respond to GRB alerts within seconds, when the afterglow is potentially extremely bright. Each mirror and camera coating is being optimized for a different wavelength range and function, including a NIR imager, two red-optimized imagers, a blue-optimized imager, an UV-optimized imager, and an optical polarimeter. PROMPT will be able to identify high-redshift events by dropout and distinguish these events from the similar signatures of extinction. In this way, PROMPT will act as a distance-finder scope for spectroscopic follow up on the larger 4.1-meter diameter SOAR telescope, which is also located at CTIO. When not chasing GRBs, PROMPT serves broader educational objectives across the state of North Carolina. Enclosure construction and the first two telescopes are now complete and functioning: PROMPT observed Swift's first GRB in December 2004. We upgrade from two to four telescope in February 2005 and from four to six telescopes in mid-2005.

We numerically compute light and polarisation curves of gamma-ray burst afterglows for various configurations of the jet luminosity structure and for different dynamical evolutions. We especially consider the standard homogeneous ``top hat'' jet and the ``universal structured jet'' with power-law wings. We also investigate a possible more physical variation of the ``top hat'' model: the ``Gaussian jet''. The polarisation curves for the last two jet types are shown here for the first time together with the computation of X-ray and radio polarised fluxes. We show that the lightcurves of the total flux from these configurations are very similar to each other, and therefore only very high quality data could allow us to pin down the underlying jet structure. We demonstrate instead that polarisation curves are a powerful means to solve the jet structure, since the predicted behaviour of polarisation and its position angle at times around the jet break are very different if not opposite. We conclude that the afterglow polarisation measurements provide clear footprints of any outflow energy distribution (unlike the lightcurves of the total flux) and the joint analysis of the total and polarised flux should reveal GRBs jet structure.

The polarization curve of GRB 020813 is discussed and compared to different models for the structure, evolution and magnetisation properties of the jet and the interstellar medium onto which the fireball impacts. GRB 020813 is best suited for this kind of analysis for the smoothness of its afterglow light curve, ensuring the applicability of current models. The polarization dataset allows us to rule out the standard GRB jet, in which the energy and Lorentz factor have a well defined value inside the jet opening angle and the magnetic field is generated at the shock front. We explore alternative models finding that a structured jet or a jet with a toroidal component of the magnetic field can fit equally well the polarization curve. Stronger conclusions cannot be drawn due to the incomplete sampling of the polarization curve. A more dense sampling, especially at early times, is required to pin down the structure of the jet and the geometry of its magnetic field.

Herein I present numerical calculations of lightcurves of homogeneous and structured afterglows with various lateral expansion rates as seen from any vantage point. Such calculations allow for direct simulation of observable quantities for complex afterglows with arbitrary energy distributions and lateral expansion paradigms. A simple, causal model is suggested for lateral expansion of the jet as it evolves; namely, that the lateral expansion kinetic energy derives from the forward kinetic energy. As such the homogeneous jet model shows that lateral expansion is important at all times in the afterglow evolution and that analytical scaling laws do a poor job at describing the afterglow decay before and after the break. In particular, I find that lateral expansion does not cause a break in the lightcurve as had been predicted. A primary purpose of this paper is to study structured afterglows, which do a good job of reproducing global relationships and correlations in the data and thus suggest the possibility of a universal afterglow model. Simulations of structured jets show a general trend in which jet breaks become more pronounced with increasing viewing angle with respect to the jet axis. In fact, under certain conditions a bump can occur in the lightcurve at the jet break time. I derive scaling relations for this bump and suggest that it may be a source of some bumps in observed lightcurves such as that of GRB 000301C. A couple of lateral expansion models are tested over a range of efficiencies and viewing angles and it is found that lateral expansion can, in some cases, substantially sharpen the jet break. I show flux surface contour maps and simulated images of the afterglows which give insight into how they evolve and determine their lightcurves.

It has been recently proposed that current burst and afterglow observations can be accounted for by a top hat jet structure as well as by a more universal structured outflow, in which the energy carried by the flow is maximum at the jet spine, decreasing outwards as a power-law with the off-axis angle. To pin down the GRBs jet structure is important in order to quantify the total energy released by the unknown engine and put constrains on its nature. In paper we will review the most popular GRBs jet structures, focusing on the predicted observed properties that could discriminate between them. In particular the behaviour of linearly polarised emission during the afterglow phase seems very promising, since the degree and temporal dependence are predicted to be very different from model to model. In addition, polarization measurements and detections are accumulating fast enouqgh to lead soon to a comprehension of magnetic field and GRBs jet configurations.

In this paper we present a model for the short (< 1 second) population of gamma-ray bursts. In this model heated neutron stars in a close binary system near its last stable orbit emit a large amount of neutrinos (~10^53 ergs). A fraction of these neutrinos will annihilate to form an e+e- pair plasma wind which will, in turn, expand and recombine to photons which make the gamma-ray burst. We study neutrino annihilation and show that a substantial fraction (~1/2) of energy deposited into e+e- pairs comes from inter-star neutrinos, where each member of the neutrino pair originates from each neutron star. Thus, in addition to the annihilation of neutrinos blowing off of a single star, there is a new source of baryon-free plasma that is deposited between the stars. To model the e+e- pair plasma wind between stars, we do three-dimensional relativistic numerical hydrodynamic calculations. We find that the time scale for these bursts, deriving from the baryon-free plasma, is less than one second and they will have a hot spectrum ~5 MeV. The energy in bursts is the order of 10^52 ergs.

A strong correlation is reported between gamma-ray burst (GRB) pulse lags and afterglow jet-break times for the set of bursts (seven) with known redshifts, luminosities, pulse lags, and jet-break times. This may be a valuable clue toward understanding the connection between the burst and afterglow phases of these events. The relation is roughly linear (i.e. doubling the pulse lag in turn doubles the jet break time) and thus implies a simple relationship between these quantities. We suggest that this correlation is due to variation among bursts of emitter Doppler factor. Specifically, an increased speed or decreased angle of velocity, with respect to the observed line-of-site, of burst ejecta will result in shorter perceived pulse lags in GRBs as well as quicker evolution of the external shock of the afterglow to the time when the jet becomes obvious, i.e. the jet-break time. Thus this observed variation among GRBs may result from a perspective effect due to different observer angles of a morphologically homogeneous populations of GRBs. Also, a conjecture is made that peak luminosities not only vary inversely with burst timescale, but also are directly proportional to the spectral break energy. If true, this could provide important information for explaining the source of this break.

In this paper we present a model for the short (< second) population of gamma-ray bursts (GRBs). In this model heated neutron stars in a close binary system near their last stable orbit emit neutrinos at large luminosities (~ 10^53 ergs/sec). A fraction of these neutrinos will annihilate to form an electron-positron pair plasma wind which will, in turn, expand and recombine to photons which make the gamma-ray burst. We study neutrino annihilation and show that a substantial fraction (~ 50%) of energy deposited comes from inter-star neutrinos, where each member of the neutrino pair originates from each neutron star. Thus, in addition to the annihilation of neutrinos blowing off of a single star, we have a new source of baryon free energy that is deposited between the stars. To model the pair plasma wind between stars, we do three-dimensional relativistic numerical hydrodynamic calculations. Preliminary results are also presented of new, fully general relativistic calculations of gravitationally attracting stars falling from infinity with no angular momentum. These simulations exhibit a compression effect.

A correlation is presented between the pulse lag and the jet-break time for seven BATSE gamma-ray bursts with known redshifts. This is, to our best knowledge, the first known direct tight correlation between a property of the gamma-ray burst phase (the pulse lag) and the afterglow phase (the jet-break time). As pulse lag and luminosity have been found to be correlated this also represents a correlation between peak luminosity and jet-break time. Observed timescales (variability or spectral lags) as well as peak luminosity naturally have a strong dependence on the Lorentz factor of the outflow and so we propose that much of the variety among GRBs has a purely kinematic origin (the speed or direction of the outflow). We explore a model in which the variation among GRBs is due to a variation in jet-opening angles, and find that the narrowest jets have the fastest outflows. We also explore models in which the jets have similar morphology and size, and the variation among bursts is caused by variation in viewing angle and/or due to a velocity profile. The relations between luminosity, variability, spectral lag and jet-break time can be qualitatively understood from models in which the Lorentz factor decreases as a function of angle from the jet axis. One expects to see high luminosities, short pulse lags and high variability as well as an early jet-break time for bursts viewed on axis, while higher viewing inclinations will yield lower luminosities, longer pulse lags, smoother bursts and later jet-break times.

We calculate the neutrino pair annihilation rate into electron pairs between two neutron stars in a binary system. We present a closed formula for the energy deposition rate at any point between the stars, where each neutrino of a pair derives from each star, and compare this result to that where all neutrinos derive from a single neutron star. An approximate generalization of this formula is given to include the relativistic effects of gravity. We find that this inter-star neutrino annihilation is a significant contributor to the energy deposition between heated neutron star binaries. In particular, for two neutron stars near their last stable orbit, inter-star neutrino annihilation energy deposition is almost equal to that of single star energy deposition.

A model is proposed for gamma-ray bursts based upon a neutrino burst of about 10^52 ergs lasting a few seconds above a heated collapsing neutron star. This type of thermal neutrino burst is suggested by relativistic hydrodynamic studies of the compression, heating, and collapse of close binary neutron stars as they approach their last stable orbit, but may arise from other sources as well. We present a hydrodynamic simulation of the formation and evolution of the pair plasma associated with such a neutrino burst. This pair plasma leads to the production of ~10^51 - 10^52 ergs in gamma-rays with spectral and temporal properties consistent with observed gamma-ray bursts.

In this paper I put forward a model in which GRB980425 is both associated with SN1998bw and is also a standard canonical (long; ~seconds) gamma-ray burst. Herein it is argued that if gamma-ray bursts are relativistic jets with the fastest moving material at the core, then the range of observed jet inclinations to the line-of-sight produces a range in the observed properties of GRBs, i.e. the lag-luminosity relationship. In particular, if the jet inclination is high enough, the observed emitter will move slowly enough to render relativistic beaming ineffective, thus distinguishing the jet from apparent isotropic emission. Thus we expect a break in the lag-luminosity relationship. I propose that GRB980425 defines that break. The position of this break gives important physical parameters such as the Lorentz factor (\gamma_max ~ 1000), the jet opening angle (~1 degree), and thus the beaming fraction (~10^-4). Estimates of burst rates are consistent with observation. If correct, this model is evidence in favor of the collapsar mode as the progenitor of cosmological, long gamma-ray bursts.

This paper presents an interpretation based on gamma-ray burst source kinematics for the relationship found by Norris et al. between peak luminosity and energy-dependent pulse lag. I argue that the correlation should instead be between "number" luminosity and pulse lag. This interpretation improves the least-squares fit of this correlation for the known bursts by 25 percent or more. It also suggests a distance estimation scheme. I propose that this relationship is due to the variation in line-of-sight velocity among bursts. This interpretation allows one to speculate on the range of gamma-ray burst expansion velocities or the size of their jet opening angles.

The interaction of an expanding Pair-Electromagnetic pulse (PEM pulse) with a shell of baryonic matter surrounding a Black Hole with electromagnetic structure (EMBH) is analyzed for selected values of the baryonic mass at selected distances well outside the dyadosphere of an EMBH. The dyadosphere, the region in which a super critical field exists for the creation of electron-positron pairs, is here considered in the special case of a Reissner-Nordstrom geometry. The interaction of the PEM pulse with the baryonic matter is described using a simplified model of a slab of constant thickness in the laboratory frame (constant-thickness approximation) as well as performing the integration of the general relativistic hydrodynamical equations. The validation of the constant-thickness approximation, already presented in a previous paper Ruffini, et al.(1999) for a PEM pulse in vacuum, is here generalized to the presence of baryonic matter. It is found that for a baryonic shell of mass-energy less than 1% of the total energy of the dyadosphere, the constant-thickness approximation is in excellent agreement with full general relativistic computations. The approximation breaks down for larger values of the baryonic shell mass, however such cases are of less interest for observed Gamma Ray Bursts (GRBs). On the basis of numerical computations of the slab model for PEM pulses, we describe (i) the properties of relativistic evolution of a PEM pulse colliding with a baryonic shell; (ii) the details of the expected emission energy and observed temperature of the associated GRBs for a given value of the EMBH mass; 10^3 solar masses, and for baryonic mass-energies in the range 10^{-8} to 10^{-2} the total energy of the dyadosphere.

In this paper we model the emission from a relativistically expanding electron-positron pair plasma fireball originating near the surface of a heated neutron star. This pair fireball is deposited via the annihilation of neutrino pairs emanating from the surface of the hot neutron star. The heating of neutron stars may occur in close neutron star binary systems near their last stable orbit. We model the relativistic expansion and subsequent emission of the plasma and find 10^51 to 10^52 ergs in gamma-rays are produced with spectral and temporal properties consistent with observed gamma-ray bursts.

General relativistic calculations are made of neutrino-antineutrino annihilation into electron-positron pairs near the surface of a neutron star. It is found that the efficiency of this process is enhanced over the Newtonian values up to a factor of more than 4 in the regime applicable to Type II supernovae and by up to a factor of 30 for collapsing neutron stars.

We study the relativistically expanding electron-positron pair plasma formed by the process of vacuum polarization around an electromagnetic black hole (EMBH). Such processes can occur for EMBH's with mass all the way up to $6\cdot 10^5M_\odot$. Beginning with a idealized model of a Reissner-Nordstrom EMBH with charge to mass ratio $\xi=0.1$, numerical hydrodynamic calculations are made to model the expansion of the pair-electromagnetic pulse (PEM pulse) to the point that the system is transparent to photons. Three idealized special relativistic models have been compared and contrasted with the results of the numerically integrated general relativistic hydrodynamic equations. One of the three models has been validated: a PEM pulse of constant thickness in the laboratory frame is shown to be in excellent agreement with results of the general relativistic hydrodynamic code. It is remarkable that this precise model, starting from the fundamental parameters of the EMBH, leads uniquely to the explicit evaluation of the parameters of the PEM pulse, including the energy spectrum and the astrophysically unprecedented large Lorentz factors (up to $6\cdot 10^3$ for a $10^3 M_{\odot}$ EMBH). The observed photon energy at the peak of the photon spectrum at the moment of photon decoupling is shown to range from 0.1 MeV to 4 MeV as a function of the EMBH mass. Correspondingly the total energy in photons is in the range of $10^{52}$ to $10^{54}$ ergs, consistent with observed gamma-ray bursts. In these computations we neglect the presence of baryonic matter which will be the subject of forthcoming publications.

Using hydrodynamic computer codes, we study the possible patterns of relativistic expansion of an enormous pair-electromagnetic-pulse (P.E.M. pulse); a hot, high density plasma composed of photons, electron-positron pairs and baryons deposited near a charged black hole (EMBH). On the bases of baryon-loading and energy conservation, we study the bulk Lorentz factor of expansion of the P.E.M. pulse by both numerical and analytical methods.

In this paper we present the preliminary results of a model for the production of gamma-ray bursts (GRBs) through the compressional heating of binary neutron stars near their last stable orbit prior to merger.

We report on general relativistic hydrodynamic studies which indicate several new physical processes which may contribute to powering gamma-ray bursts in neutron star binaries. Relativistically driven compression, heating, and collapse of the individual stars can occur many seconds before inspiral and merger. This compression may produce a neutrino burst of $\sim 10^{53}$ ergs lasting several seconds. The associated thermal neutrino emission produces an $e^+-e^-$ pair plasma by $\nu \bar \nu$ annihilation. We show first results of a simulated burst which produces $\sim 10^{51}$ erg in $\gamma$-rays. We also discuss a preliminary study of the evolution of the magnetic field lines attached to the fluid as the stars orbit. We show that the relativistically driven fluid motion might lead to the formation of extremely strong magnetic fields ($\sim 10^{17}$ gauss) in and around the stars which could affect to the formation and evolution of a gamma-ray burst.