Long Term Evolution of Magnetic Turbulence in Relativistic Collisionless Shocks

We study the long term evolution of magnetic fields generated by an initially unmagnetized collisionless relativistic $e^+e^-$ shock. Our 2D particle-in-cell numerical simulations show that downstream of such a Weibel-mediated shock, particle distributions are approximately isotropic, relativistic Maxwellians, and the magnetic turbulence is highly intermittent spatially, nonpropagating, and decaying. Using linear kinetic theory, we find a simple analytic form for these damping rates. Our theory predicts that overall magnetic energy decays like $(\omega_p t)^{-q}$ with $q \sim 1$, which compares favorably with simulations, but predicts overly rapid damping of short wavelength modes. Magnetic trapping of particles within the magnetic structures may be the origin of this discrepancy. We conclude that initially unmagnetized relativistic shocks in electron-positron plasmas are unable to form persistent downstream magnetic fields. These results put interesting constraints on synchrotron models for the prompt and afterglow emission from GRBs.Comment: 4 pages, 3 figures, contributed talk at the workshop: High Energy Phenomena in Relativistic Outflows (HEPRO), Dublin, 24-28 September 2007; Downsampled version for arXiv. Full resolution version available at http://astro.berkeley.edu/~pchang/proceedings.pd

Time Variability in the X-ray Nebula Powered by Pulsar B1509-58

We use new and archival Chandra and ROSAT data to study the time variability of the X-ray emission from the pulsar wind nebula (PWN) powered by PSR B1509-58 on timescales of one week to twelve years. There is variability in the size, number, and brightness of compact knots appearing within 20" of the pulsar, with at least one knot showing a possible outflow velocity of ~0.6c (assuming a distance to the source of 5.2 kpc). The transient nature of these knots may indicate that they are produced by turbulence in the flows surrounding the pulsar. A previously identified prominent jet extending 12 pc to the southeast of the pulsar increased in brightness by 30% over 9 years; apparent outflow of material along this jet is observed with a velocity of ~0.5c. However, outflow alone cannot account for the changes in the jet on such short timescales. Magnetohydrodynamic sausage or kink instabilities are feasible explanations for the jet variability with timescale of ~1.3-2 years. An arc structure, located 30"-45" north of the pulsar, shows transverse structural variations and appears to have moved inward with a velocity of ~0.03c over three years. The overall structure and brightness of the diffuse PWN exterior to this arc and excluding the jet has remained the same over the twelve year span. The photon indices of the diffuse PWN and possibly the jet steepen with increasing radius, likely indicating synchrotron cooling at X-ray energies.Comment: accepted to ApJ, 14 pages, 8 figure

Experimentally Observed Instability of a Laminar Ekman Flow in a Rotating Basin

In studying the axi-symmetric flow induced by source-sink distributions in a rotating cylindrical basin in the absence of radial barriers, a highly organized pattern of instability has been observed in the laminar Ekman layer along the bottom of the basin. The instability manifests itself in the form of almost perfectly concentric cylindrical sheets or curtains of water which rise as sharply defined vertical jets from the Ekman layer and penetrate the entire depth of fluid. A less sharply defined downward motion between the curtains completes the circulation celis thus developed. At some maximum critical radius, the curtains usually disappear, and the flow at larger radii is a stable, laminar Ekman flow. Quantitative observations of ring spacing and critical radius are reported for experiments in which angular velocity, flow rate, viscosity and total depth of water were varied over experimentally available ranges

Magnetar Driven Bubbles and the Origin of Collimated Outflows from GRBs

We model the interaction between the wind from a newly formed rapidly rotating magnetar and the surrounding progenitor. In the first few seconds after core collapse the magnetar inflates a bubble of plasma and magnetic fields behind the supernova shock, which expands asymmetrically because of the pinching effect of the toroidal magnetic field, as in PWNe, even if the host star is spherically symmetric. The degree of asymmetry depends on the ratio of the magnetic energy to the total energy in the bubble. We assume that the wind by newly formed magnetars inflating these bubbles is more magnetized than for PWNe. We show that for a magnetic to total power supplied by the central magnetar $\sim 0.1$ the bubble expands relatively spherically while for values greater than 0.3, most of the pressure in the bubble is exerted close to the rotation axis, driving a collimated outflow out through the host star. This can account for the collimation inferred from observations of long-duration gamma-ray bursts (GRBs). Given that the wind magnetization increases in time, we thus suggest that the magnetar-driven bubble initially expands relatively spherically (enhancing the energy of the associated supernova) while at late times it becomes progressivelymore collimated (producing the GRB). Similar processes may operate in more modestly rotating neutron stars to produce asymmetric supernovae and lower energy transients such as X-ray flashes.Comment: Proceeding of the conference "Astrophysics of Compact Objects", 1-7 July, Huangshan, Chin

Magnetic Confinement, MHD Waves, and Smooth Line Profiles in AGN

In this paper, we show that if the broad line region clouds are in approximate energy equipartition between the magnetic field and gravity, as hypothesized by Rees, there will be a significant effect on the shape and smoothness of broad emission line profiles in active galactic nuclei. Line widths of contributing clouds or flow elements are much wider than their thermal widths, due to the presence of non-dissipative MHD waves, and their collective contribution produce emission line profiles broader and smoother than would be expected if a magnetic field were not present. As an illustration, a simple model of isotropically emitting clouds, normally distributed in velocity, is used to show that smoothness can be achieved for less than 80,000 clouds and may even be as low as a few hundred. We conclude that magnetic confinement has far reaching consequences for observing and modeling active galactic nuclei.Comment: to appear in MNRA

Stokes tomography of radio pulsar magnetospheres. II. Millisecond pulsars

The radio polarization characteristics of millisecond pulsars (MSPs) differ significantly from those of non-recycled pulsars. In particular, the position angle (PA) swings of many MSPs deviate from the S-shape predicted by the rotating vector model, even after relativistic aberration is accounted for, indicating that they have non-dipolar magnetic geometries, likely due to a history of accretion. Stokes tomography uses phase portraits of the Stokes parameters as a diagnostic tool to infer a pulsar's magnetic geometry and orientation. This paper applies Stokes tomography to MSPs, generalizing the technique to handle interpulse emission. We present an atlas of look-up tables for the Stokes phase portraits and PA swings of MSPs with current-modified dipole fields, filled core and hollow cone beams, and two empirical linear polarization models. We compare our look-up tables to data from 15 MSPs and find that the Stokes phase portraits for a current-modified dipole approximately match several MSPs whose PA swings are flat or irregular and cannot be reconciled with the standard axisymmetric rotating vector model. PSR J1939+2134 and PSR J0437$-$4715 are modelled in detail. The data from PSR J1939+2134 at 0.61\,GHz can be fitted well with a current-modified dipole at $(\alpha, i) = (22 \pm 2^\circ, 80 \pm 1^\circ)$ and emission altitude 0.4 $r_\text{LC}$. The fit is less accurate for PSR J1939+2134 at 1.414\,GHz, and for PSR J0437$-$4715 at 1.44\,GHz, indicating that these objects may have a more complicated magnetic field geometry, such as a localized surface anomaly or a polar magnetic mountain.Comment: 38 pages, 33 figures, accepted for publication by MNRA

High Energy Neutrinos and Photons from Curvature Pions in Magnetars

We discuss the relevance of the curvature radiation of pions in strongly magnetized pulsars or magnetars, and their implications for the production of TeV energy neutrinos detectable by cubic kilometer scale detectors, as well as high energy photons.Comment: 19 pages, 4 figures, to appear in JCA

Thermal emission from low-field neutron stars

We present a new grid of LTE model atmospheres for weakly magnetic (B<=10e10G) neutron stars, using opacity and equation of state data from the OPAL project and employing a fully frequency- and angle-dependent radiation transfer. We discuss the differences from earlier models, including a comparison with a detailed NLTE calculation. We suggest heating of the outer layers of the neutron star atmosphere as an explanation for the featureless X-ray spectra of RX J1856.5-3754 and RX J0720.4-3125 recently observed with Chandra and XMM.Comment: 8 pages A&A(5)-Latex, 6 Figures, A&A in press. The model spectra presented here are available as XSPEC tables at http://www.astro.soton.ac.uk/~btg/outgoing/nsspec

Study of the accreting pulsar 4U 0115+634 with a bulk and thermal Comptonization model

Highly magnetized pulsars accreting matter in a binary system are bright sources in the X-ray band (0.1-100 keV). Despite the early comprehension of the basic emission mechanism, their spectral energy distribution is generally described by phenomenological or simplified models. We propose a study of the spectral emission from the high mass X-ray binary pulsar 4U 0115+634 by means of thermal and bulk Comptonization models based on the physical properties of such objects. For this purpose, we analyze the BeppoSAX data in the energy range 0.7-100 keV of the 1999 giant outburst, 12 days after the maximum. We model the spectral energy distribution of the system using a two-component continuum. At higher energy, above ~7 keV, the emission is due to thermal and bulk Comptonization of the seed photons produced by cyclotron cooling of the accretion column, and at lower energy, the emission is due to thermal Comptonization of a blackbody source in a diffuse halo close to the stellar surface. From the best fit parameters, we argue that the cyclotron emission is produced ~1.7 km above the stellar surface, and escapes from the column near its base, where the absorption features are generated by the interaction with the magnetic field in a surrounding halo. We find that in 4U 0115+634, the observed spectrum is dominated by reprocessed cyclotron radiation, whereas in other bright sources with stronger magnetic fields such as Her X-1, the spectrum is dominated by reprocessed bremsstrahlung.Comment: 12 pages, accepted for publication on Astronomy and Astrophysic

Using Flow Specifications of Parameterized Cache Coherence Protocols for Verifying Deadlock Freedom

We consider the problem of verifying deadlock freedom for symmetric cache coherence protocols. In particular, we focus on a specific form of deadlock which is useful for the cache coherence protocol domain and consistent with the internal definition of deadlock in the Murphi model checker: we refer to this deadlock as a system- wide deadlock (s-deadlock). In s-deadlock, the entire system gets blocked and is unable to make any transition. Cache coherence protocols consist of N symmetric cache agents, where N is an unbounded parameter; thus the verification of s-deadlock freedom is naturally a parameterized verification problem. Parametrized verification techniques work by using sound abstractions to reduce the unbounded model to a bounded model. Efficient abstractions which work well for industrial scale protocols typically bound the model by replacing the state of most of the agents by an abstract environment, while keeping just one or two agents as is. However, leveraging such efficient abstractions becomes a challenge for s-deadlock: a violation of s-deadlock is a state in which the transitions of all of the unbounded number of agents cannot occur and so a simple abstraction like the one above will not preserve this violation. In this work we address this challenge by presenting a technique which leverages high-level information about the protocols, in the form of message sequence dia- grams referred to as flows, for constructing invariants that are collectively stronger than s-deadlock. Efficient abstractions can be constructed to verify these invariants. We successfully verify the German and Flash protocols using our technique