MHD simulations of the collapsar model for GRBs

We present results from axisymmetric, time-dependent magnetohydrodynamic (MHD) simulations of the collapsar model for gamma-ray bursts. Our main conclusion is that, within the collapsar model, MHD effects alone are able to launch, accelerate and sustain a strong polar outflow. We also find that the outflow is Poynting flux-dominated, and note that this provides favorable initial conditions for the subsequent production of a baryon-poor fireball.Comment: 4 pages, to appear in proceedings of "2003 GRB Conference" (Santa Fe, NM, September 8-12, 2003), needs aipprocs LaTeX class, movies are available at http://rocinante.colorado.edu/~proga

The effect of alongcoast advection on pacific northwest shelf and slope water properties in relation to upwelling variability

The Northern California Current System experiences highly variable seasonal upwelling in addition to larger basin-scale variability, both of which can significantly affect its water chemistry. Salinity and temperature fields from a 7 year ROMS hindcast model of this region (43°N-50°N), along with extensive particle tracking, were used to study interannual variability in water properties over both the upper slope and the midshelf bottom. Variation in slope water properties was an order of magnitude smaller than on the shelf. Furthermore, the primary relationship between temperature and salinity anomalies in midshelf bottom water consisted of variation in density (cold/salty versus warm/fresh), nearly orthogonal to the anomalies along density levels (cold/fresh versus warm/salty) observed on the upper slope. These midshelf anomalies were well-explained (R2=0.6) by the combination of interannual variability in local and remote alongshore wind stress, and depth of the California Undercurrent (CUC) core. Lagrangian analysis of upper slope and midshelf bottom water shows that both are affected simultaneously by large-scale alongcoast advection of water through the northern and southern boundaries. The amplitude of anomalies in bottom oxygen and dissolved inorganic carbon (DIC) on the shelf associated with upwelling variability are larger than those associated with typical variation in alongcoast advection, and are comparable to observed anomalies in this region. However, a large northern intrusion event in 2004 illustrates that particular, large-scale alongcoast advection anomalies can be just as effective as upwelling variability in changing shelf water properties on the interannual scale

Stellar Explosions by Magnetic Towers

We propose a magnetic mechanism for the collimated explosion of a massive star relevant for GRBs, XRFs and asymmetric supernovae. We apply Lynden-Bell's magnetic tower scenario to the interior of a massive rotating star after the core has collapsed to form a black hole with an accretion disk or a millisecond magnetar acting as a central engine. We solve the force-free Grad-Shafranov equation to calculate the magnetic structure and growth of a tower embedded in a stellar environment. The pressure of the toroidal magnetic field, continuously generated by differential rotation of the central engine, drives a rapid expansion which becomes vertically collimated after lateral force balance with the surrounding gas pressure is reached. The collimation naturally occurs because hoop stress concentrates magnetic field toward the rotation axis and inhibits lateral expansion. This leads to the growth of a self-collimated magnetic tower. When embedded in a massive star, the supersonic expansion of the tower drives a strong bow shock behind which an over-pressured cocoon forms. The cocoon confines the tower by supplying collimating pressure and provides stabilization against disruption due to MHD instabilities. Because the tower consists of closed field lines starting and ending on the central engine, mixing of baryons from the cocoon into the tower is suppressed. The channel cleared by the growing tower is thus plausibly free of baryons and allows the escape of magnetic energy from the central engine through the star. While propagating down the stellar density gradient, the tower accelerates and becomes relativistic. During the expansion, fast collisionless reconnection becomes possible resulting in dissipation of magnetic energy which may be responsible for GRB prompt emission.Comment: 19 pages, 8 figures, accepted to ApJ, updated references and additional discussion adde

HI Observations of the Supermassive Binary Black Hole System in 0402+379

We have recently discovered a supermassive binary black hole system with a projected separation between the two black holes of 7.3 parsecs in the radio galaxy 0402+379. This is the most compact supermassive binary black hole pair yet imaged by more than two orders of magnitude. We present Global VLBI observations at 1.3464 GHz of this radio galaxy, taken to improve the quality of the HI data. Two absorption lines are found toward the southern jet of the source, one redshifted by 370 +/- 10 km/s and the other blueshifted by 700 +/- 10 km/s with respect to the systemic velocity of the source, which, along with the results obtained for the opacity distribution over the source, suggests the presence of two mass clumps rotating around the central region of the source. We propose a model consisting of a geometrically thick disk, of which we only see a couple of clumps, that reproduces the velocities measured from the HI absorption profiles. These clumps rotate in circular Keplerian orbits around an axis that crosses one of the supermassive black holes of the binary system in 0402+379. We find an upper limit for the inclination angle of the twin jets of the source to the line of sight of 66 degrees, which, according to the proposed model, implies a lower limit on the central mass of ~7 x 10^8 Msun and a lower limit for the scale height of the thick disk of ~12 pc .Comment: 20 pages, 7 figures. Accepted on the Astrophysical Journa

A Decline in the X-ray through Radio Emission from GW170817 Continues to Support an Off-Axis Structured Jet

We present new observations of the binary neutron star merger GW170817 at $\Delta t\approx 220-290$ days post-merger, at radio (Karl G. Jansky Very Large Array; VLA), X-ray (Chandra X-ray Observatory) and optical (Hubble Space Telescope; HST) wavelengths. These observations provide the first evidence for a turnover in the X-ray light curve, mirroring a decline in the radio emission at $\gtrsim5\sigma$ significance. The radio-to-X-ray spectral energy distribution exhibits no evolution into the declining phase. Our full multi-wavelength dataset is consistent with the predicted behavior of our previously published models of a successful structured jet expanding into a low-density circumbinary medium, but pure cocoon models with a choked jet cannot be ruled out. If future observations continue to track our predictions, we expect that the radio and X-ray emission will remain detectable until $\sim 1000$ days post-merger.Comment: Accepted to ApJL. Updated version includes new VLA observations extending through 2018 June

The Binary Neutron Star event LIGO/VIRGO GW170817 a hundred and sixty days after merger: synchrotron emission across the electromagnetic spectrum

We report deep Chandra, HST and VLA observations of the binary neutron star event GW170817 at $t<160$ d after merger. These observations show that GW170817 has been steadily brightening with time and might have now reached its peak, and constrain the emission process as non-thermal synchrotron emission where the cooling frequency $\nu_c$ is above the X-ray band and the synchrotron frequency $\nu_m$ is below the radio band. The very simple power-law spectrum extending for eight orders of magnitude in frequency enables the most precise measurement of the index $p$ of the distribution of non-thermal relativistic electrons $N(\gamma)\propto \gamma^{-p}$ accelerated by a shock launched by a NS-NS merger to date. We find $p=2.17\pm0.01$, which indicates that radiation from ejecta with $\Gamma\sim3-10$ dominates the observed emission. While constraining the nature of the emission process, these observations do \emph{not} constrain the nature of the relativistic ejecta. We employ simulations of explosive outflows launched in NS ejecta clouds to show that the spectral and temporal evolution of the non-thermal emission from GW170817 is consistent with both emission from radially stratified quasi-spherical ejecta traveling at mildly relativistic speeds, \emph{and} emission from off-axis collimated ejecta characterized by a narrow cone of ultra-relativistic material with slower wings extending to larger angles. In the latter scenario, GW170817 harbored a normal SGRB directed away from our line of sight. Observations at $t\le 200$ days are unlikely to settle the debate as in both scenarios the observed emission is effectively dominated by radiation from mildly relativistic material.Comment: Updated with the latest VLA and Chandra dat

Axisymmetric general relativistic hydrodynamics: Long-term evolution of neutron stars and stellar collapse to neutron stars and black holes

We report a new implementation for axisymmetric simulation in full general relativity. In this implementation, the Einstein equations are solved using the Nakamura-Shibata formulation with the so-called cartoon method to impose an axisymmetric boundary condition, and the general relativistic hydrodynamic equations are solved using a high-resolution shock-capturing scheme based on an approximate Riemann solver. As tests, we performed the following simulations: (i) long-term evolution of non-rotating and rapidly rotating neutron stars, (ii) long-term evolution of neutron stars of a high-amplitude damping oscillation accompanied with shock formation, (iii) collapse of unstable neutron stars to black holes, and (iv) stellar collapses to neutron stars. The tests (i)--(iii) were carried out with the $\Gamma$-law equation of state, and the test (iv) with a more realistic parametric equation of state for high-density matter. We found that this new implementation works very well: It is possible to perform the simulations for stable neutron stars for more than 10 dynamical time scales, to capture strong shocks formed at stellar core collapses, and to accurately compute the mass of black holes formed after the collapse and subsequent accretion. In conclusion, this implementation is robust enough to apply to astrophysical problems such as stellar core collapse of massive stars to a neutron star and black hole, phase transition of a neutron star to a high-density star, and accretion-induced collapse of a neutron star to a black hole. The result for the first simulation of stellar core collapse to a neutron star started from a realistic initial condition is also presented.Comment: 28 pages, to appear in PRD 67, 0440XX (2003

Can a Large Neutron Excess Help Solve the Baryon Loading Problem in Gamma-Ray Burst Fireballs?

We point out that the baryon-loading problem in Gamma-Ray Burst (GRB) models can be amelioriated if a significant fraction of the baryons which inertially confine the fireball are converted to neutrons. A high neutron fraction in some circumstances can result in a reduced transfer of energy from relativistic light particles in the fireball to baryons. The energy needed to produce the required relativistic flow in the GRB is consequently reduced, in some cases by orders of magnitude. This could be relevant to GRB models because a high neutron-to-proton ratio has been calculated in neutron star-merger fireball environments. Significant neutron excess also could occur near compact objects with high neutrino fluxes.Comment: 5 pages, 2 figures, to appear in Phys. Rev. Let