Researcher profile

John G. Kirk

John G. Kirk contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
18works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2015arXiv

Rapid cosmic-ray acceleration at perpendicular shocks in supernova remnants

Perpendicular shocks are shown to be rapid particle accelerators that perform optimally when the ratio $u_{\rm s}$ of the shock speed to the particle speed roughly equals the ratio $1/η$ of the scattering rate to the gyro frequency. We use analytical methods and Monte-Carlo simulations to solve the kinetic equation that governs the anisotropy generated at these shocks, and find, for $ηu_{\rm s}\approx1$, that the spectral index softens by unity and the acceleration time increases by a factor of two compared to the standard result of diffusive shock acceleration theory. These results provide a theoretical basis for the thirty-year-old conjecture that a supernova exploding into the wind of a Wolf-Rayet star may accelerate protons to an energy exceeding $10^{15}\,$eV.

preprint2013arXiv

High-energy emission from pulsar binaries

Unpulsed, high-energy emission from pulsar binaries can be attributed to the interaction of a pulsar wind with that of a companion star. At the shock between the outflows, particles carried away from the pulsar magnetosphere are accelerated and radiate both in synchrotron and inverse Compton processes. This emission constitutes a significant fraction of the pulsar spin-down luminosity. It is not clear however, how the highly magnetized pulsar wind could convert its mainly electromagnetic energy into the particles with such high efficiency. Here we investigate a scenario in which a pulsar striped wind converts into a strong electromagnetic wave before reaching the shock. This mode can be thought of as a shock precursor that is able to accelerate particles to ultrarelativistic energies at the expense of the electromagnetic energy it carries. Radiation of the particles leads to damping of the wave. The efficiency of this process depends on the physical conditions imposed by the external medium. Two regimes can be distinguished: a high density one, where the EM wave cannot be launched at all and the current sheets in the striped wind are first compressed by an MHD shock and subsequently dissipate by reconnection, and a low density one, where the wind can first convert into an electromagnetic wave in the shock precursor, which then damps and merges into the surroundings. Shocks in binary systems can transit from one regime to another according to binary phase. We discuss possible observational implications for these objects.

preprint2013arXiv

Propagation and stability of superluminal waves in pulsar winds

Nonlinear electromagnetic waves with superluminal phase velocity can propagate in the winds around isolated pulsars, and around some pulsars in binary systems. Using a short-wavelength approximation, we find and analyze an integrable system of equations that govern their evolution in spherical geometry. A confined mode is identified that stagnates to finite pressure at large radius and can form a precursor to the termination shock. Using a simplified criterion, we find this mode is stable for most isolated pulsars, but may be unstable if the external pressure is high, such as in the pulsar wind nebulae in starburst galaxies and in W44. Pulsar winds in eccentric binary systems, such as PSR 1259-63, may go through phases with stable and unstable electromagnetic precursors, as well as phases in which the density is too high for these modes to propagate.

preprint2013arXiv

Radiative damping and emission signatures of strong superluminal waves in pulsar winds

We analyse the damping by radiation reaction and by Compton drag of strong, superluminal electromagnetic waves in the context of pulsar winds. The associated radiation signature is found by estimating the efficiency and the characteristic radiation frequencies. Applying these estimates to the gamma-ray binary containing PSR B1259-63, we show that the GeV flare observed by Fermi-LAT can be understood as inverse Compton emission by particles scattering photons from the companion star, if the pulsar wind termination shock acquires a precursor of superluminal waves roughly 30 days after periastron. This constrains the mass-loading factor of the wind $μ=L/\dot{N}mc^2$ (where $L$ is the luminosity and $\dot{N}$ the rate of loss of electrons and positrons) to be roughly $6\times 10^4$.

preprint2013arXiv

Superluminal Waves and the Structure of Pulsar Wind Termination Shocks

The termination shock of a pulsar wind is located roughly where the ram pressure matches that of the surrounding medium. Downstream of the shock, MHD models of the diffuse nebular emission suggest the plasma is weakly magnetized. However, the transition from a Poynting-dominated MHD wind to a particle-dominated flow is not well understood. We discuss a solution of this "sigma problem" in which a striped wind converts into a strong, superluminal electromagnetic wave. This mode slows down as it propagates radially, and its ram pressure tends to a constant value at large radius, a property we use to match the solution to the surrounding nebula. The wave thus forms a pre-cursor to the termination shock, which occurs at the point where the wave dissipates. Possible damping and dissipation mechanisms are discussed qualitatively.

preprint2013arXiv

The Role of Superluminal Electromagnetic Waves in Pulsar Wind Termination Shocks

The dynamics of a standing shock front in a Poynting-flux dominated relativistic flow is investigated by using a one-dimensional, relativistic, two-fluid simulation. An upstream flow containing a circularly polarized, sinusoidal magnetic shear wave is considered, mimicking a wave driven by an obliquely rotating pulsar. It is demonstrated that this wave is converted into large amplitude electromagnetic waves with superluminal phase speeds by interacting with the shock when the shock-frame frequency of the wave exceeds the proper plasma frequency. The superluminal waves propagate in the upstream, modify the shock structure substantially, and form a well-developed precursor region ahead of a subshock. Dissipation of Poynting flux occurs in the precursor as well as in the downstream region through a parametric instability driven by the superluminal waves. The Poynting flux remaining in the downstream region is carried entirely by the superluminal waves. The downstream plasma is therefore an essentially unmagnetized, relativistically hot plasma with a non-relativistic flow speed, as suggested by observations of pulsar wind nebulae.

preprint2013arXiv

The stability of strong waves and its implications for pulsar wind shocks

Strong waves can mediate a shock transition between a pulsar wind and its surroundings, playing the role of an extended precursor, in which the energy is effectively transferred from fields to non-thermal particles. The damping of such precursors results in an essentially unmagnetized shock near the equator. In this context, we discuss the stability of strong waves and its implications for the properties of shocks. Those with stable precursors can exist in the winds of most of isolated pulsars, but the precursors may be unstable if the external pressure in the nebula is high, as in Vela-like pulsars. Pulsar wind shocks in eccentric binary systems, such as B1259-63, can acquire precursors only at certain orbital phases, and this process should be accompanied by enhanced synchro-Compton and inverse Compton emission from the precursor. The same scenario may be at work in the binary HESS J0632+057.

preprint2012arXiv

Linearly polarized superluminal waves in pulsar winds

Pulsar winds are the ideal environment for the study of non-linear electromagnetic waves. It is generally thought that a pulsar launches a striped wind, a magnetohydrodynamic entropy wave, where plasma sheets carried along with the flow separate regions of alternating magnetic field. But when the density drops below a critical value, or equivalently for distances from the pulsar greater than a critical radius, a strong superluminal wave can also propagate. In this contribution we discuss the conversion of the equatorial striped wind into a linearly polarized superluminal wave, and we argue that this mode is important for the conversion of Poynting flux to kinetic energy flux before the outflow reaches the termination shock.

preprint2011arXiv

Charge-starved, relativistic jets and blazar variability

High energy emission from blazars is thought to arise in a relativistic jet launched by a supermassive black hole. The emission site must be far from the hole and the jet relativistic, in order to avoid absorption of the photons. In extreme cases, rapid variability of the emission suggests that structures of length-scale smaller than the gravitational radius of the central black hole are imprinted on the jet as it is launched, and modulate the radiation released after it has been accelerated to high Lorentz factor. We propose a mechanism which can account for the acceleration of the jet, and for the rapid variability of the radiation, based on the propagation characteristics of large-amplitude waves in charge-starved, polar jets. Using a two-fluid (electron-positron) description, we find the outflows exhibit a delayed acceleration phase, that starts at roughly 1pc, where the inertia associated with the wave currents becomes important. The time-structure imprinted on the jet at launch modulates photons produced by the accelerating jet provided the pair multiplicity in the black-hole magnetosphere is sufficiently small, suggesting that very rapid variability is confined to sources in which the electromagnetic cascade in the black-hole magnetosphere is not prolific.

preprint2011arXiv

Superluminal waves in pulsar winds

The energy lost by a rotation-powered pulsar is carried by a relativistic flow containing a mixture of electromagnetic fields and particles. In the inner regions, this is thought to be a magnetically dominated, cold, electron-positron wind that is well described by the MHD equations. However, beyond a critical radius r_{cr}, the same particle, energy and momentum fluxes can be transported by a strong, transverse electromagnetic wave with superluminal phase speed. We analyze the nonlinear dispersion relation of these waves for linear and circular polarization, and find the dependence of r_{cr} on the mass-loading, magnetization and luminosity of the flow, as well as on the net magnetic flux. We show that, for most isolated pulsars, the wind lies well outside r_{cr}, and speculate that superluminal modes play an important role in the dissipation of electromagnetic energy into nonthermal particles at the termination shock.

preprint2010arXiv

Computation of synthetic spectra from simulations of relativistic shocks

Particle-in-cell (PIC) simulations of relativistic shocks are in principle capable of predicting the spectra of photons that are radiated incoherently by the accelerated particles. The most direct method evaluates the spectrum using the fields given by the Lienard-Wiechart potentials. However, for relativistic particles this procedure is computationally expensive. Here we present an alternative method, that uses the concept of the photon formation length. The algorithm is suitable for evaluating spectra both from particles moving in a specific realization of a turbulent electromagnetic field, or from trajectories given as a finite, discrete time series by a PIC simulation. The main advantage of the method is that it identifies the intrinsic spectral features, and filters out those that are artifacts of the limited time resolution and finite duration of input trajectories.

preprint2010arXiv

Linear acceleration emission in pulsar magnetospheres

(abbrev.) Linear acceleration emission occurs when a charged particle is accelerated parallel to its velocity. We evaluate the spectral and angular distribution of this radiation for several special cases, including constant acceleration (hyperbolic motion) of finite duration. Based on these results, we find the following general properties of the emission from an electron in a linear accelerator that can be characterized by an electric field E acting over a distance L: (i) the spectrum extends to a cut-off photon energy ~ LE^2 MeV, where E is in units of the Schwinger critical field and L in units of the Compton wavelength of the electron. (ii) the total energy emitted by a particle traversing the accelerator is in agreement with the standard Larmor formula (iii) the low frequency spectrum is flat for hyperbolic trajectories, but in general depends on the details of the accelerator. We also show that linear acceleration emission complements curvature radiation in the strongly magnetized pair formation regions in pulsar magnetospheres. It dominates when the length L of the accelerator is less than the formation length of curvature photons, which is given by the ratio of the radius of curvature of the magnetic field lines to the Lorentz factor of the particle. In standard static models of pair creating regions linear acceleration emission is negligible, but it is important in more realistic dynamical models in which the accelerating field fluctuates on a short length-scale.

preprint2010arXiv

Radiative Signatures of Relativistic Shocks

(Abbreviated) Particle-in-cell simulations of relativistic, weakly magnetized collisionless shocks show that particles can gain energy by repeatedly crossing the shock front. This requires scattering off self-generated small length-scale magnetic fluctuations. The radiative signature of this first-order Fermi acceleration mechanism is important for models of both the prompt and afterglow emission in gamma-ray bursts and depends on the strength parameter &#34;a&#34; of the fluctuations. For electrons (and positrons), acceleration saturates when the radiative losses produced by the scattering cannot be compensated by the energy gained on crossing the shock. We show that this sets an upper limit on both the electron Lorentz factor and on the energy of the photons radiated during the scattering process. This rules out &#34;jitter&#34; radiation on self-excited fluctuations with a < 1 as a source of gamma-rays, although high-energy photons might still be produced when the jitter photons are upscattered in an analog of the synchrotron self-Compton process. In fluctuations with a > 1, radiation is generated by the standard synchrotron mechanism, and the maximum photon energy rises linearly with a, until saturating at approximately 70 MeV.

preprint2009arXiv

Two-Dimensional particle-in-cell simulations of the nonresonant, cosmic-ray driven instability in SNR shocks

In supernova remnants, the nonlinear amplification of magnetic fields upstream of collisionless shocks is essential for the acceleration of cosmic rays to the energy of the &#34;knee&#34; at 10^{15.5}eV. A nonresonant instability driven by the cosmic ray current is thought to be responsible for this effect. We perform two-dimensional, particle-in-cell simulations of this instability. We observe an initial growth of circularly polarized non-propagating magnetic waves as predicted in linear theory. It is demonstrated that in some cases the magnetic energy density in the growing waves, can grow to at least 10 times its initial value. We find no evidence of competing modes, nor of significant modification by thermal effects. At late times we observe saturation of the instability in the simulation, but the mechanism responsible is an artefact of the periodic boundary conditions and has no counterpart in the supernova-shock scenario.

preprint2008arXiv

Gamma-ray polarization constraints on Planck scale violations of special relativity

Using recent polarimetric observations of the Crab Nebula in the hard X-ray band by INTEGRAL, we show that the absence of vacuum birefringence effects constrains O(E/M) Lorentz violation in QED to the level |ξ| < 9x10^{-10} at three sigma CL, tightening by more than three orders of magnitude previous constraints. We show that planned X-ray polarimeters have the potential the potential to probe |ξ|~ 10^{-16} by detecting polarization in active galaxies at red-shift ~1.

preprint2001arXiv

Particle acceleration by ultra-relativistic shocks: theory and simulations

We consider the acceleration of charged particles near ultra-relativistic shocks, with Lorentz factor Gamma_s >> 1. We present simulations of the acceleration process and compare these with results from semi-analytical calculations. We show that the spectrum that results from acceleration near ultra-relativistic shocks is a power law, N(E) \propto E^{-s}, with a nearly universal value s \approx 2.2 - 2.3 for the slope of this power law. We confirm that the ultra-relativistic equivalent of Fermi acceleration at a shock differs from its non-relativistic counterpart by the occurence of large anisotropies in the distribution of the accelerated particles near the shock. In the rest frame of the upstream fluid, particles can only outrun the shock when their direction of motion lies within a small loss cone of opening angle theta_c \approx 1/Gamma_s around the shock normal. We also show that all physically plausible deflection or scattering mechanisms can change the upstream flight direction of relativistic particles originating from downstream by only a small amount: Delta theta ~ 1/Gamma_s. This limits the energy change per shock crossing cycle to Delta E ~ E, except for the first cycle where particles originate upstream. In that case the upstream energy is boosted by a factor ~ Gamma_s^2 for those particles that are scattered back across the shock into the upstream region.

preprint1999arXiv

Evidence for Intergalactic Absorption in the TeV Gamma-Ray Spectrum of Mkn 501

The recent HEGRA observations of the blazar Mkn 501 show strong curvature in the very high energy gamma-ray spectrum. Applying the gamma-ray opacity derived from an empirically based model of the intergalactic infrared background radiation field (IIRF), to these observations, we find that the intrinsic spectrum of this source is consistent with a power-law: dN/dE~ E^-alpha with alpha=2.00 +/- 0.03 over the range 500 GeV - 20 TeV. Within current synchrotron self-Compton scenarios, the fact that the TeV spectral energy distribution of Mkn 501 does not vary with luminosity, combined with the correlated, spectrally variable emission in X-rays, as observed by the BeppoSAX and RXTE instruments, also independently implies that the intrinsic spectrum must be close to alpha=2. Thus, the observed curvature in the spectrum is most easily understood as resulting from intergalactic absorption.