Source author record

J. G. Kirk

J. G. Kirk appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

23 published item(s)

preprint2016arXiv

Radiative trapping in intense laser beams

The dynamics of electrons in counter-propagating, circularly polarized laser beams are shown to exhibit attractors whose ability to trap particles depends on the ratio of the beam intensities and a single parameter describing radiation reaction. Analytical expressions are found for the underlying limit cycles and the parameter range in which they are stable. In high-intensity optical pulses, where radiation reaction strongly modifies the trajectories, the production of collimated gamma-rays and the initiation of non-linear cascades of electron-positron pairs can be optimized by a suitable choice of the intensity ratio.

preprint2015arXiv

Quantum radiation reaction in laser-electron beam collisions

It is possible using current high intensity laser facilities to reach the quantum radiation reaction regime for energetic electrons. An experiment using a wakefield accelerator to drive GeV electrons into a counterpropagating laser pulse would demonstrate the increase in the yield of high energy photons caused by the stochastic nature of quantum synchrotron emission: we show that a beam of $10^9$ 1 GeV electrons colliding with a 30 fs laser pulse of intensity $10^{22}~\text{Wcm}^{-2}$ will emit 6300 photons with energy greater than 700 MeV, $60\times$ the number predicted by classical theory.

preprint2013arXiv

Dense electron-positron plasmas and bursts of gamma-rays from laser-generated QED plasmas

In simulations of a 12.5PW laser (focused intensity I = 4x10^23W/cm^2) striking a solid aluminium target 10% of the laser energy is converted to gamma-rays. A dense electron-positron plasma is generated with a maximum density of 10^26/m^3; seven orders of magnitude denser than pure e-e+ plasmas generated with 1PW lasers. When the laser power is increased to 320PW (I = 10^25W/cm^2) 40% of the laser energy is converted to gamma-ray photons and 10% to electron-positron pairs. In both cases there is strong feedback between the QED emission processes and the plasma physics; the defining feature of the new `QED-plasma' regime reached in these interactions.

preprint2013arXiv

Modelling Gamma Ray Emission and Pair Production in High-Intensity Laser-Matter Interactions

In high-intensity (> 10^21W/cm^2) laser-matter interactions gamma-ray photon emission by the electrons can strongly affect the electron's dynamics and copious numbers of electron-positron pairs can be produced by the emitted photons. We show how these processes can be included in simulations by coupling a Monte-Carlo algorithm describing the emission to a particle-in-cell code. The Monte-Carlo algorithm includes quantum corrections to the photon emission, which we show must be included if the pair production rate is to be correctly determined. The accuracy, convergence and energy conservation properties of the Monte-Carlo algorithm are analysed in simple test problems.

preprint2013arXiv

Pair plasma cushions in the hole-boring scenario

Pulses from a 10 PW laser are predicted to produce large numbers of gamma-rays and electron-positron pairs on hitting a solid target. However, a pair plasma, if it accumulates in front of the target, may partially shield it from the pulse. Using stationary, one-dimensional solutions of the two-fluid (electron-positron) and Maxwell equations, including a classical radiation reaction term, we examine this effect in the hole-boring scenario. We find the collective effects of a pair plasma "cushion" substantially reduce the reflectivity, converting the absorbed flux into high-energy gamma-rays. There is also a modest increase in the laser intensity needed to achieve threshold for a non-linear pair cascade.

preprint2012arXiv

Dense Electron-Positron Plasmas and Ultra-Intense Bursts of Gamma-Rays from Laser-Irradiated Solids

In simulations of a 10PW laser striking a solid we demonstrate the possibility of producing a pure electron-positron plasma by the same processes as those thought to operate in high-energy astrophysical environments. A maximum positron density of 10^26/m^3 is achieved, seven orders of magnitude greater than achieved in previous experiments. Additionally, 35% of the laser energy is converted to a burst of gamma-rays of intensity 10^22W/cm^2, potentially the most intense gamma-ray source available in the laboratory. This absorption results in a strong feedback between both pair and gamma-ray production and classical plasma physics in the new `QED-plasma' regime.

preprint2010arXiv

Waves in pulsar winds

The radio, optical, X-ray and gamma-ray nebulae that surround many pulsars are thought to arise from synchrotron and inverse Compton emission. The energy powering this emission, as well as the magnetic fields and relativistic particles, are supplied by a "wind" driven by the central object. The inner parts of the wind can be described using the equations of MHD, but these break down in the outer parts, when the density of charge carriers drops below a critical value. This paper reviews the wave properties of the inner part (striped wind), and uses a relativistic two-fluid model (cold electrons and positrons) to re-examine the nonlinear electromagnetic modes that propagate in the outer parts. It is shown that in a radial wind, two solutions exist for circularly polarised electromagnetic modes. At large distances one of them turns into a freely expanding flow containing a vacuum wave, whereas the other decelerates, corresponding to a confined flow.

preprint2008arXiv

Steady-State Solutions in Nonlinear Diffusive Shock Acceleration

Stationary solutions to the equations of non-linear diffusive shock acceleration play a fundamental role in the theory of cosmic-ray acceleration. Their existence usually requires that a fraction of the accelerated particles be allowed to escape from the system. Because the scattering mean-free-path is thought to be an increasing function of energy, this condition is conventionally implemented as an upper cut-off in energy space -- particles are then permitted to escape from any part of the system, once their energy exceeds this limit. However, because accelerated particles are responsible for substantial amplification of the ambient magnetic field in a region upstream of the shock front, we examine an alternative approach in which particles escape over a spatial boundary. We use a simple iterative scheme that constructs stationary numerical solutions to the coupled kinetic and hydrodynamic equations. For parameters appropriate for supernova remnants, we find stationary solutions with efficient acceleration when the escape boundary is placed at the point where growth and advection of strongly driven non-resonant waves are in balance. We also present the energy dependence of the distribution function close to the energy where it cuts off - a diagnostic that is in principle accessible to observation.

preprint2007arXiv

Automatic quenching of high energy gamma-ray sources by synchrotron photons

We investigate a magnetized plasma in which injected high energy gamma-rays annihilate on a soft photon field, that is provided by the synchrotron radiation of the created pairs. For a very wide range of magnetic fields, this process involves gamma-rays between 0.3 GeV and 30 TeV. We derive a simple dynamical system for this process, analyze its stability to runaway production of soft photons and paris, and find conditions for it to automatically quench by reaching a steady state with an optical depth to photon-photon annihilation larger than unity. We discuss applications to broad-band gamma-ray emitters, in particular supermassive black holes. Automatic quenching limits the gamma-ray luminosity of these objects and predicts substantial pair loading of the jets of less active sources.

preprint2007arXiv

The theory of pulsar winds and nebulae

We review current theoretical ideas on pulsar winds and their surrounding nebulae. Relativistic MHD models of the wind of the aligned rotator, and of the striped wind, together with models of magnetic dissipation are discussed. It is shown that the observational signature of this dissipation is likely to be point-like, rather than extended, and that pulsed emission may be produced. The possible pulse shapes and polarisation properties are described. Particle acceleration at the termination shock of the wind is discussed, and it is argued that two distinct mechanisms must be operating, with the first-order Fermi mechanism producing the high-energy electrons (above 1 TeV) and either magnetic annihilation or resonant absorption of ion cyclotron waves responsible for the 100 MeV to 1 TeV electrons. Finally, MHD models of the morphology of the nebula are discussed and compared with observation.

preprint2005arXiv

Gamma-rays from PSR B1259 -63

The high-energy gamma-ray emission discovered using the H.E.S.S. telescopes from the binary system PSR B1259 -63, is modelled using an extension of the approach that successfully predicted it. We find that the simultaneous INTEGRAL and H.E.S.S. data permit both a model with dominant radiative losses, high pulsar wind Lorentz factor and modest efficiency as well as one with dominant adiabatic losses, a slower wind and higher efficiency. Additional, simultaneous, X-ray and TeV data sets are needed to lift this degeneracy.

preprint2005arXiv

The high-energy gamma-ray light curve of PSR B1259 -63

The high-energy gamma-ray light curve of the binary system PSR B1259 -63, is computed using the approach that successfully predicted the spectrum at periastron. The simultaneous INTEGRAL and H.E.S.S. spectra taken 16 days after periastron currently permit both a model with dominant radiative losses, high pulsar wind Lorentz factor and modest efficiency as well as one with dominant adiabatic losses, a slower wind and higher efficiency. In this paper we shown how the long-term light curve may help to lift this degeneracy.

preprint2001arXiv

The beaming pattern and spectrum of inverse Compton radiation in blazars

By including Klein-Nishina effects, we generalize previous calculations of the beaming pattern of photons produced by inverse Compton scattering. For an isotropic distribution of soft photons upscattered by nonthermal electrons with a power-law density distribution embedded in a plasma moving with relativistic bulk speed we show that the observed radiation intensity is proportional to D^(3+p), where D is the Doppler boosting factor, and p the electron distribution index. This agrees with previous computations performed in the Thomson limit, where the observed spectral index is a=(p-1)/2 and the beaming pattern D^(4+2a). Independent of D, Klein-Nishina effects limit the location of the peak energy of the observed spectral energy distribution such that e_{peak} less than 1/e_0, where e_0 is the energy of the seed photons in units of mc^2. Assuming that the seed photons originate in the broad line region, we demonstrate that the GeV emission of blazars is significantly modified by Klein-Nishina effects, the spectrum being softer than that calculated in the Thomson limit. We further show that the change in spectral index of the inverse Compton emission across the peak energy can exceed the value of 0.5 predicted by computations performed in the Thomson limit. The model spectra agree with OSSE and COMPTEL limits on this break without invoking the effects of differential absorption at the edge of a gamma-ray photosphere.

preprint1999arXiv

Particle acceleration at oblique shocks and discontinuities of the density profile

In the theory of diffusive acceleration at oblique shock fronts the question of the existence of a discontinuity of energetic particle density is contentious. The resolution of this problem is interesting from a theoretical point of view, and potentially for the interpretation of observations of particle densities at heliospheric shocks and of high-resolution radio observations of the rims of supernova remnants. It can be shown analytically that an isotropic particle distribution at a shock front implies continuity of the particle density -- whether or not the shock is oblique. However, if the obliquity of the shock induces an anisotropy, a jump is permitted. Both semi-analytic computations and Monte-Carlo simulations are used to show that, for interesting parameter ranges, a jump is indeed produced, with accelerated particles concentrated in a precursor ahead of the shock front.

preprint1998arXiv

Inverse Compton Emission of TeV Gamma Rays from PSR B1259-63

We derive light curves for the hard gamma-ray emission, at energies up to several TeV, expected from the unique pulsar/Be-star binary system PSR B1259-63. This is the only known system in our galaxy in which a radio pulsar is orbiting a main sequence star. We show that inverse Compton emission from the electrons and positrons in the shocked pulsar wind, scattering target photons from the Be star, produces a flux of hard gamma-rays that should be above the sensitivity threshold of present day atmospheric Cerenkov detectors. Furthermore, we predict that the flux of hard gamma-rays produced via this mechanism has a characteristic variation with orbital phase that should be observable, and which is not expected from any other mechanism.

preprint1998arXiv

Particle acceleration and synchrotron emission in blazar jets

We model the acceleration of electrons at a shock front in a relativistic blazar jet and compute the radiation they emit in a post-shock region which contains a homogeneous magnetic field. The full space, time and momentum dependence of the electron distribution is used in this calculation. It is shown that the "homogeneous" synchrotron model is recovered, provided the downstream speed of the plasma away from the shock front is nonrelativistic, and provided that the light travel times across the face of the shock front is unimportant. By varying the rate at which particles are picked up by the acceleration process, we calculate the time-dependence of the spectra. Since the magnetic field strength is assumed constant within the emission region, each frequency band can be identified with electrons of a particular energy. We find that for a band in which the electrons are accelerated rapidly compared to the rate at which they cool, the spectra typically harden during phases of rising flux, and soften during phases of falling flux, as has been observed in the objects PKS 2155-304 and Mkn 421. However, in a frequency band in which the timescales are comparable, the reverse behaviour is to be expected. We discuss the extent to which observations of both the stationary spectrum and the spectral variability of the synchrotron component of blazar emission can be used to constrain the model.

preprint1997arXiv

Anomalous transport of cosmic ray electrons

Anomalous transport processes in which the variance of the distance travelled does not necessarily increase linearly with time are modelled using the formalism of continuous time random walks. We compute particle propagators which have the required dependence on space and time and use these to find the spatial dependence of the synchrotron radiation emitted by a population of continuously injected electrons. As the electrons are transported away from the source they cool, and the synchrotron spectrum softens. Sub-diffusive transport -- corresponding to stochastic trapping, or restriction of the transport across the average direction of a stochastic magnetic field -- produces a much slower rate of change of spectral index than does supra-diffusion -- which occurs when particles move almost without scattering, in a field containing large ordered regions. Application to the diffuse emission of the outer parts of the Coma cluster favours an interpretation in terms of supra-diffusion.

preprint1996arXiv

On the production of very high energy beamed gamma-rays in blazars

The variable flux of TeV gamma-rays detected from Mkn 421 and Mkn 501 requires the presence of high energy electrons, which could in principle produce large numbers of electron/positron pairs, leading to an electromagnetic cascade. We point out that this scenario can be avoided if electrons are accelerated to high energy rectilinearly, rather than being injected isotropically into a blob, as in most of the models of the GeV gamma-ray emission. By balancing linear acceleration by an electric field against inverse Compton losses in the radiation field of the accretion disk we calculate the emitted spectra and find the conditions which must be fulfilled in order to exclude the development of electromagnetic cascades during acceleration. Assuming these to be fulfilled, we show that the maximum possible photon energy is approximately $10M_8^{2/5}\,$TeV, where $M_8$ is the mass of the central black hole in units of $10^8\,\msolar$. In addition we compute the optical depth to absorption of TeV photons on a possible isotropic scattered component and on the observed nonthermal radiation (in the case of Mkn 421) and find that TeV photons can escape provided the nonthermal X-rays originate in a jet moving with a Lorentz factor of around 8.

preprint1996arXiv

Variability in the synchrotron self-Compton model of blazar emission

We present a model of the spectra of gamma-ray emitting blazars in which a single homogeneous emission region both emits synchrotron photons directly and scatters them to high (gamma-ray) energy before emission (a ``synchrotron self-Compton'' or SSC model). In contrast to previous work, we follow the full time dependent evolution of the electron and photon spectra, assuming a power-law form of the electron injection and examine the predictions of the model with regard to variability of the source. We apply these computations to the object Mkn 421, which displayed rapid variability in its X-ray and TeV emission during a multiwavelength campaign in 1994. This observation strongly implies that the same population of electrons produces the radiation in both energy bands. By fitting first the observed quiescent spectrum over all 18 orders of magnitude in frequency, we show that the time dependence of the keV/TeV flare could have been the result of a sudden increase in the maximum energy of the injected electrons. We show also that different types of flare may occur in this object and others, and that the energy band most sensitive to the properties of the acceleration mechanism is the X-ray band.

preprint1995arXiv

A MODEL FOR THE RADIO EMISSION FROM SNR 1987A

The observations of radio emission from SNR~1987A can be accounted for on the basis of diffusive shock acceleration of electrons by the supernova blast wave. However, with this interpretation the observed spectral index implies that the compression ratio of the gas subshock is roughly $2.7$ rather than the value of $4$ expected of a strong shock front. We propose that in SNR~1987A, ions also undergo diffusive acceleration at the shock, a process that is likely to be rapid. Unlike the electron population, the accelerated ions can have an important effect on the gas dynamics. We calculate this coupled gas and energetic particle dynamics on the basis of the two-fluid model, in which the accelerated ions provide an additional component to the total pressure acting on the fluid. By accelerating and possibly heating the upstream plasma, the initially strong shock is modified and a weaker subshock with an upstream precursor results. The electrons behave as test particles. They are accelerated at the evolving subshock, escape downstream, and emit synchrotron radiation in the swept up magnetic field. Two models are considered for the surroundings of the progenitor: that of a freely expanding wind of number density $n\propto r^{-2}$, and that of a wind confined by a shell of denser material, creating a stagnation zone of roughly constant density beyond the standing shock which terminates the free wind. We model the observed radio light curves and the relatively steep spectrum of SNR~1987A using similar values for the ion acceleration parameters to those used in models of cosmic ray acceleration in older SNRs which can also contain high Mach number shocks, and find a good fit for the case in which the termination shock is located at about $2\times 10^{15}\,$m from the progenitor.

preprint1995arXiv

Anomalous Transport and Particle Acceleration at Shocks

The theory of first order Fermi acceleration at shocks assumes that particles diffuse due to scattering off slow-moving magnetic irregularities. However, cosmic rays are closely tied to magnetic field lines, and the transport process, particularly across the direction of the field, is likely to be more complicated. To describe cross field transport we employ recent extensions of the Rechester-Rosenbluth theory in localised stochastic regions of magnetic field. During acceleration at a shock and when the motion along the field is diffusive, there is a transition at a critical energy from \lq\lq sub-diffusive\rq\rq\ motion, where the mean square displacement of a particle increases with time as $t^{1/2}$, to compound diffusion, a combined process involving diffusion along a magnetic field which is itself wandering. Requiring this critical energy to be less than the system cut-off in a SNR of radius $R$ places an upper limit on the coherence length, $λ$, of the magnetic field for which diffusive shock acceleration can occur at quasi-perpendicular shocks when field line wandering dominates the effective transport of particles: $λ<R(U/v)^{1/2}$ where $U$ is the shock speed and $v$ the particle speed. The acceleration rate for particles in the sub-diffusive regime is derived and its implication for the origin of high energy cosmic rays is discussed.

preprint1995arXiv

The Acceleration of Electrons In Radio Supernova SN1986J

We propose a model for radio supernovae (RSN) based on the synchrotron emission from relativistic electrons which are diffusively accelerated at the expanding supernova shock. This model was originally developed for application to the optically thin emission observed from SN1987A. Here we generalise it by including the effects of free-free absorption from both an external screen and from material internal to the source, and by relaxing the restriction to an azimuthal B-field. We find a good fit to the entire set of radio data for the best observed highly-luminous RSN -- SN1986J -- with a reduced χ^2 of 3.85. Applying the new model to SN1988Z, another intrinsically bright RSN, also yields a good fit (χ^2_{\rm red}\approx 2) but this is less significant, because of the limited data on this distant (z=0.02) source. These fits suggest that the shock expands at constant speed, that the magnetic field within the source decreases with time according to t^{-2}, and that the compression ratio of the shock front is close to the value expected of a strong shock in an ideal gas of adiabatic index 5/3 -- indicating a relatively low value of the cosmic ray pressure compared with SN1987A. In the case of SN1986J we derive an explosion date in August/September 1982, a magnetic field at the position of the shock 1000 days after explosion of 4nT and a spatial diffusion coefficient of the electrons of 4 10^{19} cm^2/s, four orders of magnitude greater than the Bohm value. In addition, we obtain the optical depths to external and internal absorption, and derive an estimate of the mass-loss rate.

preprint1994arXiv

Radio Supernovae as Tev Gamma-Ray Sources

When applied to the blast wave formed by the explosion of a massive star as a supernova (SN), the theory of diffusive particle acceleration at shock fronts predicts a very high energy density in cosmic rays. Almost immediately after particles begin to be injected into the process, the cosmic ray pressure rises until comparable to the ram-pressure encountered by the shock front. Those supernovae which are observed in the radio band i.e., radio supernovae (RSNe), provide direct evidence of particle acceleration in the form of synchrotron emitting electrons. Furthermore, these objects are particularly interesting, since they are usually surrounded by a relatively dense confining medium. The acceleration of cosmic rays can then lead to the production of very high energy (VHE) gamma-rays which arise from collisions between energetic particles and target nuclei. We estimate the cosmic ray energy density assuming a fraction $ϕ\lesim1$ of the energy available at the shock front is converted into cosmic rays. Combining this with the parameters describing the environment of the SN progenitor, as deduced from observations, and from more detailed modelling, we compute the flux at Earth $\tevflux$ of photons of energy above $1\,$TeV. For the relatively weak but nearby supernova SN1987A we predict $\tevflux=2\times10^{-13}\,\fluxunits$ before the shock front encounters the ring of dense matter seen by the Hubble Space Telescope. Subsequently, the flux is expected to rise further. The medium around SN1993J in M81 is thought to have a density