-
Thermodynamic limits on oxygenic photosynthesis around M-dwarf stars: Generalized models and strategies for optimization
Authors:
Samir Chitnavis,
Thomas J. Haworth,
Edward Gillen,
Conrad W. Mullineaux,
Christopher D. P. Duffy
Abstract:
We explore the feasibility and potential characteristics of photosynthetic light-harvesting on exo-planets orbiting in the habitable zone of low mass stars ($< 1$ M$_{\odot}$). As stellar temperature, $T_{s}$, decreases, the irradiance maximum red-shifts out of the $400 \textrm{nm} \leq λ< 750$ nm range of wavelengths that can be utilized by \emph{oxygenic} photosynthesis on Earth. However, limite…
▽ More
We explore the feasibility and potential characteristics of photosynthetic light-harvesting on exo-planets orbiting in the habitable zone of low mass stars ($< 1$ M$_{\odot}$). As stellar temperature, $T_{s}$, decreases, the irradiance maximum red-shifts out of the $400 \textrm{nm} \leq λ< 750$ nm range of wavelengths that can be utilized by \emph{oxygenic} photosynthesis on Earth. However, limited irradiance in this region does not preclude oxygenic photosynthesis and Earth's plants, algae and cyanobacteria all possess very efficient \emph{light-harvesting antennae} that facilitate photosynthesis in very low light. Here we construct general models of photosynthetic light-harvesting structures to determine how an oxygenic photosystem would perform in different irradiant spectral fluxes. We illustrate that the process of light-harvesting, capturing energy over a large antenna and concentrating it into a small \emph{reaction centre}, must overcome a fundamental \emph{entropic barrier}. We show that a plant-like antenna cannot be adapted to the light from stars of $T_{s}<3400$ K, as increasing antenna size offers diminishing returns on light-harvesting. This can be overcome if one introduces a slight \emph{enthalpic gradient}, to the antenna. Interestingly, this strategy appears to have been adopted by Earth's oxygenic cyanobacteria, and we conclude that \emph{bacterial} oxygenic photosynthesis is feasible around even the lowest mass M-dwarf stars.
△ Less
Submitted 22 September, 2023;
originally announced September 2023.
-
Photosynthesis Under a Red Sun: Predicting the absorption characteristics of an extraterrestrial light-harvesting antenna
Authors:
Christopher D. P. Duffy,
Gregoire Canchon,
Thomas J. Haworth,
Edward Gillen,
Samir Chitnavis,
Conrad W. Mullineaux
Abstract:
Here we discuss the feasibility of photosynthesis on Earth-like rocky planets in close orbit around ultra-cool red dwarf stars. Stars of this type have very limited emission in the \textit{photosynthetically active} region of the spectrum ($400 - 700$ nm), suggesting that they may not be able to support oxygenic photosynthesis. However, photoautotrophs on Earth frequently exploit very dim environm…
▽ More
Here we discuss the feasibility of photosynthesis on Earth-like rocky planets in close orbit around ultra-cool red dwarf stars. Stars of this type have very limited emission in the \textit{photosynthetically active} region of the spectrum ($400 - 700$ nm), suggesting that they may not be able to support oxygenic photosynthesis. However, photoautotrophs on Earth frequently exploit very dim environments with the aid of highly structured and extremely efficient antenna systems. Moreover, the anoxygenic photosynthetic bacteria, which do not need to oxidize water to source electrons, can exploit far red and near infrared light. Here we apply a simple model of a photosynthetic antenna to a range of model stellar spectra, ranging from ultra-cool (2300 K) to Sun-like (5800 K). We assume that a photosynthetic organism will evolve an antenna that maximizes the rate of energy input while also minimizing fluctuations. The latter is the 'noise cancelling' principle recently reported by Arp et al. 2020. Applied to the Solar spectrum this predicts optimal antenna configurations in agreement with the chlorophyll Soret absorption bands. Applied to cooler stars, the optimal antenna peaks become redder with decreasing stellar temperature, crossing to the typical wavelength ranges associated with anoxygenic photoautotrophs at $\sim 3300$ K. Lastly, we compare the relative input power delivered by antennae of equivalent size around different stars and find that the predicted variation is within the same order of magnitude. We conclude that low-mass stars do not automatically present light-limiting conditions for photosynthesis but they may select for anoxygenic organisms.
△ Less
Submitted 3 May, 2023;
originally announced May 2023.
-
Particle Acceleration in Relativistic Shearing Flows: Energy Spectrum
Authors:
Frank M. Rieger,
Peter Duffy
Abstract:
We consider the acceleration of charged particles in relativistic shearing flows, with Lorentz factor up to $Γ_0 \sim 20$. We present numerical solutions to the particle transport equation and compare these with results from analytical calculations. We show that in the highly relativistic limit the particle energy spectrum that results from acceleration approaches a power law,…
▽ More
We consider the acceleration of charged particles in relativistic shearing flows, with Lorentz factor up to $Γ_0 \sim 20$. We present numerical solutions to the particle transport equation and compare these with results from analytical calculations. We show that in the highly relativistic limit the particle energy spectrum that results from acceleration approaches a power law, $N(E)\propto E^{-\tilde{q}}$, with a universal value $\tilde{q}=(1+α)$ for the slope of this power law, where $α$ parameterizes the power-law momentum dependence of the particle mean free path. At mildly relativistic flow speeds, the energy spectrum becomes softer and sensitive to the underlying flow profile. We explore different flow examples, including Gaussian and power-law-type velocity profiles, showing that the latter yield comparatively harder spectra, producing $\tilde{q}\simeq 2$ for $Γ_0 \simeq 3$ and Kolmogorov turbulence. We provide a comparison with a simplified leaky-box approach and derive an approximate relation for estimating the spectral index as a function of the maximum shear flow speed. These results are of relevance for jetted, high-energy astrophysical sources such as active galactic nuclei, since shear acceleration is a promising mechanism for the acceleration of charged particles to relativistic energies and is likely to contribute to the high-energy radiation observed.
△ Less
Submitted 27 June, 2022;
originally announced June 2022.
-
Progenitor, environment, and modelling of the interacting transient, AT 2016jbu (Gaia16cfr)
Authors:
S. J. Brennan,
M. Fraser,
J. Johansson,
A. Pastorello,
R. Kotak,
H. F. Stevance,
T. -W. Chen,
J. J. Eldridge,
S. Bose,
P. J. Brown,
E. Callis,
R. Cartier,
M. Dennefeld,
Subo Dong,
P. Duffy,
N. Elias-Rosa,
G. Hosseinzadeh,
E. Hsiao,
H. Kuncarayakti,
A. Martin-Carrillo,
B. Monard,
G. Pignata,
D. Sand,
B. J. Shappee,
S. J. Smartt
, et al. (45 additional authors not shown)
Abstract:
We present the bolometric lightcurve, identification and analysis of the progenitor candidate, and preliminary modelling of AT2016jbu (Gaia16cfr). We find a progenitor consistent with a $\sim$22--25~$M_{\odot}$ yellow hypergiant surrounded by a dusty circumstellar shell, in agreement with what has been previously reported. We see evidence for significant photometric variability in the progenitor,…
▽ More
We present the bolometric lightcurve, identification and analysis of the progenitor candidate, and preliminary modelling of AT2016jbu (Gaia16cfr). We find a progenitor consistent with a $\sim$22--25~$M_{\odot}$ yellow hypergiant surrounded by a dusty circumstellar shell, in agreement with what has been previously reported. We see evidence for significant photometric variability in the progenitor, as well as strong H$α$ emission consistent with pre-existing circumstellar material. The age of the environment as well as the resolved stellar population surrounding AT2016jbu, support a progenitor age of $>$10 Myr, consistent with a progenitor mass of $\sim$22~$M_{\odot}$. A joint analysis of the velocity evolution of AT2016jbu, and the photospheric radius inferred from the bolometric lightcurve shows the transient is consistent with two successive outbursts/explosions. The first outburst ejected material with velocity $\sim$650$kms^{-1}$, while the second, more energetic event, ejected material at $\sim$4500$kms^{-1}$. Whether the latter is the core-collapse of the progenitor remains uncertain. We place a limit on the ejected $^{56}$Ni mass of $<$0.016$M_{\odot}$. Using the BPASS code, we explore a wide range of possible progenitor systems, and find that the majority of these are in binaries, some of which are undergoing mass transfer or common envelope evolution immediately prior to explosion. Finally, we use the SNEC code to demonstrate that the low-energy explosion within some of these binary systems, together with sufficient CSM, can reproduce the overall morphology of the lightcurve of AT2016jbu.
△ Less
Submitted 27 April, 2022; v1 submitted 18 February, 2021;
originally announced February 2021.
-
Photometric and spectroscopic evolution of the interacting transient AT 2016jbu (Gaia16cfr)
Authors:
S. J. Brennan,
M. Fraser,
J. Johansson,
A. Pastorello,
R. Kotak,
H. F. Stevance,
T. -W. Chen,
J. J. Eldridge,
S. Bose,
P. J. Brown,
E. Callis,
R. Cartier,
M. Dennefeld,
Subo Dong,
P. Duffy,
N. Elias-Rosa,
G. Hosseinzadeh,
E. Hsiao,
H. Kuncarayakti,
A. Martin-Carrillo,
B. Monard,
A. Nyholm,
G. Pignata,
D. Sand,
B. J. Shappee
, et al. (46 additional authors not shown)
Abstract:
We present the results from a high cadence, multi-wavelength observation campaign of AT 2016jbu (aka Gaia16cfr), an interacting transient. This dataset complements the current literature by adding higher cadence as well as extended coverage of the lightcurve evolution and late-time spectroscopic evolution. Photometric coverage reveals that AT 2016jbu underwent significant photometric variability f…
▽ More
We present the results from a high cadence, multi-wavelength observation campaign of AT 2016jbu (aka Gaia16cfr), an interacting transient. This dataset complements the current literature by adding higher cadence as well as extended coverage of the lightcurve evolution and late-time spectroscopic evolution. Photometric coverage reveals that AT 2016jbu underwent significant photometric variability followed by two luminous events, the latter of which reached an absolute magnitude of M$_V\sim$-18.5 mag. This is similar to the transient SN 2009ip whose nature is still debated. Spectra are dominated by narrow emission lines and show a blue continuum during the peak of the second event. AT 2016jbu shows signatures of a complex, non-homogeneous circumstellar material (CSM). We see slowly evolving asymmetric hydrogen line profiles, with velocities of 500km$s^{-1}$ seen in narrow emission features from a slow moving CSM, and up to 10,000km$s^{-1}$ seen in broad absorption from some high velocity material. Late-time spectra ($\sim$+1 year) show a lack of forbidden emission lines expected from a core-collapse supernova and are dominated by strong emission from H, He i and Ca ii. Strong asymmetric emission features, a bumpy lightcurve, and continually evolving spectra suggest an inhibit nebular phase. We compare the evolution of H$α$ among SN 2009ip-like transients and find possible evidence for orientation angle effects. The light-curve evolution of AT 2016jbu suggests similar, but not identical, circumstellar environments to other SN 2009ip-like transients.
△ Less
Submitted 27 April, 2022; v1 submitted 18 February, 2021;
originally announced February 2021.
-
Turbulence and Particle Acceleration in Shearing Flows
Authors:
Frank M. Rieger,
Peter Duffy
Abstract:
We explore constraints imposed by shear-driven instabilities on the acceleration of energetic particles in relativistic shearing flows. We show that shearing layers in large-scale AGN jets are likely to encompass a sizeable fraction ($\geq 0.1$) of the jet radius, requiring seed injection of GeV electrons for efficient acceleration. While the diffusion process may depend on pre-developed turbulenc…
▽ More
We explore constraints imposed by shear-driven instabilities on the acceleration of energetic particles in relativistic shearing flows. We show that shearing layers in large-scale AGN jets are likely to encompass a sizeable fraction ($\geq 0.1$) of the jet radius, requiring seed injection of GeV electrons for efficient acceleration. While the diffusion process may depend on pre-developed turbulence if injection occurs at higher energies, electron acceleration to PeV and proton acceleration to EeV energies appears possible within the constraints imposed by jet stability.
△ Less
Submitted 15 January, 2021;
originally announced January 2021.
-
Particle Acceleration in Shearing Flows: Efficiencies and Limits
Authors:
Frank M. Rieger,
Peter Duffy
Abstract:
We examine limits to the efficiency for particles acceleration in shearing flows, showing that relativistic flow speeds are required for efficient gradual shear acceleration. We estimate maximum achievable particle energies for parameters applicable to relativistic AGN jets. The implications of our estimates is that if large-scale jets are relativistic, then efficient electron acceleration up to s…
▽ More
We examine limits to the efficiency for particles acceleration in shearing flows, showing that relativistic flow speeds are required for efficient gradual shear acceleration. We estimate maximum achievable particle energies for parameters applicable to relativistic AGN jets. The implications of our estimates is that if large-scale jets are relativistic, then efficient electron acceleration up to several PeV, and proton acceleration up to several EeV energies appears feasible. This suggests that shear particle acceleration could lead to a continued energization of synchrotron X-ray emitting electrons, and be of relevance for the production of ultra-high-energy cosmic-ray particles.
△ Less
Submitted 13 November, 2019;
originally announced November 2019.
-
Thermal emission from bow shocks I: 2D Hydrodynamic Models of the Bubble Nebula
Authors:
Samuel Green,
Jonathan Mackey,
Thomas J. Haworth,
Vasilii V. Gvaramadze,
Peter Duffy
Abstract:
The Bubble Nebula (or NGC 7635) is a parsec-scale seemingly spherical wind-blown bubble around the relatively unevolved O star BD+60$^\circ$2522. The small dynamical age of the nebula and significant space velocity of the star suggest that the Bubble Nebula might be a bow shock. We have run 2D hydrodynamic simulations to model the interaction of the central star's wind with the interstellar medium…
▽ More
The Bubble Nebula (or NGC 7635) is a parsec-scale seemingly spherical wind-blown bubble around the relatively unevolved O star BD+60$^\circ$2522. The small dynamical age of the nebula and significant space velocity of the star suggest that the Bubble Nebula might be a bow shock. We have run 2D hydrodynamic simulations to model the interaction of the central star's wind with the interstellar medium (ISM). The models cover a range of possible ISM number densities of $n=50-200 {\rm cm}^{-3}$ and stellar velocities of $v_{\star}=20-40$ km s$^{-1}$. Synthetic H$α$ and 24 $μ$m emission maps predict the same apparent spherical bubble shape with quantitative properties similar to observations. The synthetic maps also predict a maximum brightness similar to that from the observations and agree that the maximum brightness is at the apex of the bow shock. The best-matching simulation had $v_{\star}\approx20$ km s$^{-1}$ into an ISM with $n\sim100 {\rm cm}^{-3}$, at an angle of 60$^\circ$ with respect to the line of sight. Synthetic maps of soft ($0.3-2$ keV) and hard ($2-10$ keV) X-ray emission show that the brightest region is in the wake behind the star and not at the bow shock itself. The unabsorbed soft X-rays have luminosity $\sim10^{32}-10^{33}$ erg s$^{-1}$. The hard X-rays are fainter, luminosity $\sim 10^{30} - 10^{31}$ erg s$^{-1}$, and may be too faint for current X-ray instruments to successfully observe. Our results imply that the O star creates a bow shock as it moves through the ISM and in turn creates an asymmetric bubble visible at optical and infrared wavelengths, and predicted to be visible in X-rays. The Bubble Nebula does not appear to be unique, it could be just a favourably oriented very dense bow shock. The dense ISM surrounding BD+60$^\circ$2522 and its strong wind suggest that it could be a good candidate for detecting non-thermal emission.
△ Less
Submitted 13 March, 2019;
originally announced March 2019.
-
Shear Acceleration in Expanding Flows
Authors:
F. M. Rieger,
P. Duffy
Abstract:
Shear flows are naturally expected to occur in astrophysical environments and potential sites of continuous non-thermal Fermi-type particle acceleration. Here we investigate the efficiency of expanding relativistic outflows to facilitate the acceleration of energetic charged particles to higher energies. To this end, the gradual shear acceleration coefficient is derived based on an analytical trea…
▽ More
Shear flows are naturally expected to occur in astrophysical environments and potential sites of continuous non-thermal Fermi-type particle acceleration. Here we investigate the efficiency of expanding relativistic outflows to facilitate the acceleration of energetic charged particles to higher energies. To this end, the gradual shear acceleration coefficient is derived based on an analytical treatment. The results are applied to the context of the relativistic jets of active galactic nuclei. The inferred acceleration timescale is investigated for a variety of conical flow profiles (i.e., power law, Gaussian, Fermi-Dirac) and compared to the relevant radiative and non-radiative loss timescales. The results exemplify that relativistic shear flows are capable of boosting cosmic-rays to extreme energies. Efficient electron acceleration, on the other hand, requires weak magnetic fields and may thus be accompanied by a delayed onset of particle energization and affect the overall jet appearance (e.g., core, ridge line and limb-brightening).
△ Less
Submitted 14 November, 2016;
originally announced November 2016.
-
Charge Exchange Losses and Stochastic Acceleration in the Solar Wind
Authors:
Ciaran Kenny,
Peter Duffy
Abstract:
Stochastic acceleration of particles under a pressure balance condition can accommodate the universal $p^{-5}$ spectra observed under many different conditions in the inner heliosphere. In this model, in order to avoid an infinite build up of particle pressure, a relationship between the momentum diffusion of particles and the adiabatic deceleration in the solar wind must exist. This constrains bo…
▽ More
Stochastic acceleration of particles under a pressure balance condition can accommodate the universal $p^{-5}$ spectra observed under many different conditions in the inner heliosphere. In this model, in order to avoid an infinite build up of particle pressure, a relationship between the momentum diffusion of particles and the adiabatic deceleration in the solar wind must exist. This constrains both the spatial and momentum diffusion coefficients and results in the $p^{-5}$ spectrum in the presence of adiabatic losses in the solar wind. However, this theory cannot explain the presence of such spectra beyond the termination shock, where adiabatic deceleration is negligible. To explain this apparent discrepancy, we include the effect of charge exchange losses, resulting in new forms of both the spatial and momentum diffusion coefficients that have not previously been considered. Assuming that the turbulence is of a large-scale compressible nature, we find that a balance between momentum diffusion and losses can still readily lead to the creation of $p^{-5}$ suprathermal tails, including those found in the outer heliosphere.
△ Less
Submitted 5 September, 2016;
originally announced September 2016.
-
A New Time-Dependent Finite Difference Method for Relativistic Shock Acceleration
Authors:
Sean Delaney,
Paul Dempsey,
Peter Duffy,
Turlough P. Downes
Abstract:
We present a new approach to calculate the particle distribution function about relativistic shocks including synchrotron losses using the method of lines with an explicit finite difference scheme. A steady, continuous, one dimensional plasma flow is considered to model thick (modified) shocks, leading to a calculation in three dimensions plus time, the former three being momentum, pitch angle and…
▽ More
We present a new approach to calculate the particle distribution function about relativistic shocks including synchrotron losses using the method of lines with an explicit finite difference scheme. A steady, continuous, one dimensional plasma flow is considered to model thick (modified) shocks, leading to a calculation in three dimensions plus time, the former three being momentum, pitch angle and position. The method accurately reproduces the expected power law behaviour in momentum at the shock for upstream flow speeds ranging from 0.1c to 0.995c (1 < Γ< 10). It also reproduces approximate analytical results for the synchrotron cutoff shape for a non-relativistic shock, demonstrating that the loss process is accurately represented. The algorithm has been implemented as a hybrid OpenMP--MPI parallel algorithm to make efficient use of SMP cluster architectures and scales well up to many hundreds of CPUs.
△ Less
Submitted 24 November, 2011;
originally announced November 2011.
-
The non-thermal emission of extended radio galaxy lobes with curved electron spectra
Authors:
Peter Duffy,
Katherine M. Blundell
Abstract:
The existing theoretical framework for the energies stored in the synchrotron-emitting lobes of radio galaxies and quasars doesn't properly account for the curved spectral shape that many of them exhibit. We characterise these spectra using parameters that are straightforwardly observable in the era of high-resolution, low-frequency radio astronomy: the spectral curvature and the turnover in the f…
▽ More
The existing theoretical framework for the energies stored in the synchrotron-emitting lobes of radio galaxies and quasars doesn't properly account for the curved spectral shape that many of them exhibit. We characterise these spectra using parameters that are straightforwardly observable in the era of high-resolution, low-frequency radio astronomy: the spectral curvature and the turnover in the frequency spectrum. This characterisation gives the Lorentz factor at the turnover in the energy distribution (we point out that this is distinctly different from the Lorentz factor corresponding to the turnover frequency in a way that depends on the amount of curvature in the spectrum) and readily gives the equipartition magnetic field strength and the total energy of the radiating plasma obviating the need for any assumed values of the cutoff frequencies to calculate these important physical quantities. This framework readily yields the form of the X-ray emission due to inverse-Compton (IC) scattering of Cosmic Microwave Background (CMB) photons by the electrons in the plasma having Lorentz factors of $\sim$1000. We also present the contribution to CMB anisotropies due to relativistic plasmas such as giant radio galaxy lobes, expressed in terms of the extent to which the lobes have their magnetic field and particle energies are in equipartition with one another.
△ Less
Submitted 21 November, 2011;
originally announced November 2011.
-
Steady-State Solutions in Nonlinear Diffusive Shock Acceleration
Authors:
B. Reville,
J. G. Kirk,
P. Duffy
Abstract:
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 a…
▽ More
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.
△ Less
Submitted 20 December, 2008;
originally announced December 2008.
-
Particle-acceleration timescales in TeV blazar flares
Authors:
Joni Tammi,
Peter Duffy
Abstract:
Observations of minute-scale flares in TeV Blazars place constraints on particle acceleration mechanisms in those objects. The implications for a variety of radiation mechanisms have been addressed in the literature; in this paper we compare four different acceleration mechanisms: diffusive shock acceleration, second-order Fermi, shear acceleration and the converter mechanism. When the accelerat…
▽ More
Observations of minute-scale flares in TeV Blazars place constraints on particle acceleration mechanisms in those objects. The implications for a variety of radiation mechanisms have been addressed in the literature; in this paper we compare four different acceleration mechanisms: diffusive shock acceleration, second-order Fermi, shear acceleration and the converter mechanism. When the acceleration timescales and radiative losses are taken into account, we can exclude shear acceleration and the neutron-based converted mechanism as possible acceleration processes in these systems. The first-order Fermi process and the converter mechanism working via SSC photons are still practically instantaneous, however, provided sufficient turbulence is generated on the timescale of seconds. We propose stochastic acceleration as a promising candidate for the energy-dependent time delays in recent gamma-ray flares of Markarian 501.
△ Less
Submitted 1 December, 2008; v1 submitted 21 November, 2008;
originally announced November 2008.
-
Properties of vertically self-gravitating accretion discs with a dissipative corona
Authors:
Fazeleh Khajenabi,
Peter Duffy
Abstract:
The steady-state structure of a disc with a corona is analyzed when the vertical component of the gravitational force due to the self-gravity of the disc is considered. For the energy exchange between the disc and the corona, we assume a fraction f of the dissipated energy inside the accretion disc is transported to the corona via the magnetic tubes. Analytical solutions corresponding to a presc…
▽ More
The steady-state structure of a disc with a corona is analyzed when the vertical component of the gravitational force due to the self-gravity of the disc is considered. For the energy exchange between the disc and the corona, we assume a fraction f of the dissipated energy inside the accretion disc is transported to the corona via the magnetic tubes. Analytical solutions corresponding to a prescription for f (in which this parameter directly depends on the ratio of the gas pressure to the total pressure) or free f are presented and their physical properties are studied in detail. We show that the existence of the corona not only decreases the temperature of the disc, but also increases the surface density.The vertical component of the gravitational force due to the self-gravity of the disc decreases the self-gravitating radius and the mass of the fragments at this radius. However, as more energy is transported from the disc to the corona, the effect of the vertical component of the gravitational force due to the self-gravity of the disc on the self-gravitating radius becomes weaker, though the mass of the fragments is reduced irrespective of the amount of the energy exchange from the disc to the corona.
△ Less
Submitted 30 May, 2008;
originally announced May 2008.
-
Multiwavelength study of Cygnus A II. X-ray inverse-Compton emission from a relic counterjet and implications for jet duty-cycles
Authors:
K. C. Steenbrugge,
K. M. Blundell,
P. Duffy
Abstract:
The duty-cycle of powerful radio galaxies and quasars such as the prototype Cygnus A is poorly understood. X-ray observations of inverse-Compton scattered Cosmic Microwave Background (ICCMB) photons probe lower Lorentz-factor particles than radio observations of synchrotron emission. Comparative studies of the nearer and further lobes, separated by many 10s of kpc and thus 10s of thousands of ye…
▽ More
The duty-cycle of powerful radio galaxies and quasars such as the prototype Cygnus A is poorly understood. X-ray observations of inverse-Compton scattered Cosmic Microwave Background (ICCMB) photons probe lower Lorentz-factor particles than radio observations of synchrotron emission. Comparative studies of the nearer and further lobes, separated by many 10s of kpc and thus 10s of thousands of years in light-travel time, yield additional temporal resolution in studies of the lifecycles. We have co-added all archival Chandra ACIS-I data and present a deep 200 ks image of Cygnus A. This deep image reveals the presence of X-ray emission from a counterjet i.e. a jet receding from Earth and related to a previous episode of jet activity. The non-thermal X-ray emission, we interpret as ICCMB radiation. There is an absence of any discernible X-ray emission associated with a jet flowing towards Earth. We conclude that: (1) The emission from a relic jet, indicates a previous episode of jet activity, that took place earlier than the current jet activity appearing as synchrotron radio emission. (2) The presence of X-ray emission from a relic counterjet of Cygnus A and the absence of X-ray emission associated with any relic approaching jet constrains the timescale between successive episodes of jet activity to ~10^6 years. (3) Transverse expansion of the jet causes expansion losses which shifts the energy distribution to lower energies. (4) Assuming the electrons cooled due to adiabatic expansion, the required magnetic field strength is substantially smaller than the equipartition magnetic field strength. (5) A high minimum Lorentz factor for the distribution of relativistic particles in the current jet, of a few 10^3, is ejected from the central nucleus of this active galaxy. Abridged.
△ Less
Submitted 16 May, 2008; v1 submitted 14 May, 2008;
originally announced May 2008.
-
Environmental limits on the non-resonant cosmic-ray current-driven instability
Authors:
B. Reville,
J. G. Kirk,
P. Duffy,
S. O'Sullivan
Abstract:
We investigate the so-called non-resonant cosmic-ray streaming instability, first discussed by Bell (2004). The extent to which thermal damping and ion-neutral collisions reduce the growth of this instability is calculated. Limits on the growth of the non-resonant mode in SN1006 and RX J1713.7-3946 are presented.
We investigate the so-called non-resonant cosmic-ray streaming instability, first discussed by Bell (2004). The extent to which thermal damping and ion-neutral collisions reduce the growth of this instability is calculated. Limits on the growth of the non-resonant mode in SN1006 and RX J1713.7-3946 are presented.
△ Less
Submitted 22 February, 2008;
originally announced February 2008.
-
The transport of cosmic rays in self-excited magnetic turbulence
Authors:
B. Reville,
S. O'Sullivan,
P. Duffy,
J. G. Kirk
Abstract:
The process of diffusive shock acceleration relies on the efficacy with which hydromagnetic waves can scatter charged particles in the precursor of a shock. The growth of self-generated waves is driven by both resonant and non-resonant processes. We perform high-resolution magnetohydrodynamic simulations of the non-resonant cosmic-ray driven instability, in which the unstable waves are excited b…
▽ More
The process of diffusive shock acceleration relies on the efficacy with which hydromagnetic waves can scatter charged particles in the precursor of a shock. The growth of self-generated waves is driven by both resonant and non-resonant processes. We perform high-resolution magnetohydrodynamic simulations of the non-resonant cosmic-ray driven instability, in which the unstable waves are excited beyond the linear regime. In a snapshot of the resultant field, particle transport simulations are carried out. The use of a static snapshot of the field is reasonable given that the Larmor period for particles is typically very short relative to the instability growth time. The diffusion rate is found to be close to, or below, the Bohm limit for a range of energies. This provides the first explicit demonstration that self-excited turbulence reduces the diffusion coefficient and has important implications for cosmic ray transport and acceleration in supernova remnants.
△ Less
Submitted 1 February, 2008;
originally announced February 2008.
-
Detection of a relic X-ray jet in Cygnus A
Authors:
K. C. Steenbrugge,
K. M. Blundell,
P. Duffy
Abstract:
We present a 200 ks Chandra ACIS-I image of Cygnus A, and discuss a long linear feature seen in its counterlobe. This feature has a non-thermal spectrum and lies on the line connecting the brighter hotspot on the approaching side and the nucleus. We therefore conclude that this feature is (or was) a jet. However, the outer part of this X-ray jet does not trace the current counterjet observed in…
▽ More
We present a 200 ks Chandra ACIS-I image of Cygnus A, and discuss a long linear feature seen in its counterlobe. This feature has a non-thermal spectrum and lies on the line connecting the brighter hotspot on the approaching side and the nucleus. We therefore conclude that this feature is (or was) a jet. However, the outer part of this X-ray jet does not trace the current counterjet observed in radio. No X-ray counterpart is observed on the jet side. Using light-travel time effects we conclude that this X-ray 50 kpc linear feature is a relic jet that contains enough low-energy plasma (gamma ~ 10^3) to inverse-Compton scatter cosmic microwave background photons, producing emission in the X-rays.
△ Less
Submitted 18 December, 2007;
originally announced December 2007.
-
A cosmic ray current driven instability in partially ionised media
Authors:
B. Reville,
J. G. Kirk,
P. Duffy,
S. O. Sullivan
Abstract:
We investigate the growth of hydromagnetic waves driven by streaming cosmic rays in the precursor environment of a supernova remnant shock. It is known that transverse waves propagating parallel to the mean magnetic field are unstable to anisotropies in the cosmic ray distribution, and may provide a mechanism to substantially amplify the ambient magnetic field. We quantify the extent to which te…
▽ More
We investigate the growth of hydromagnetic waves driven by streaming cosmic rays in the precursor environment of a supernova remnant shock. It is known that transverse waves propagating parallel to the mean magnetic field are unstable to anisotropies in the cosmic ray distribution, and may provide a mechanism to substantially amplify the ambient magnetic field. We quantify the extent to which temperature and ionisation fractions modify this picture. Using a kinetic description of the plasma we derive the dispersion relation for a collisionless thermal plasma with a streaming cosmic ray current. Fluid equations are then used to discuss the effects of neutral-ion collisions. We calculate the extent to which the environment into which the cosmic rays propagate influences the growth of the magnetic field, and determines the range of possible growth rates. If the cosmic ray acceleration is efficient, we find that very large neutral fractions are required to stabilise the growth of the non-resonant mode. For typical supernova parameters in our galaxy, thermal effects do not significantly alter the growth rates. For weakly driven modes, ion-neutral damping can dominate over the instability at more modest ionisation fractions. In the case of a supernova shock interacting with a molecular clouds, such as in RX J1713.7-3946, with high density and low ionisation, the modes can be rapidly damped.
△ Less
Submitted 20 September, 2007; v1 submitted 25 July, 2007;
originally announced July 2007.
-
Radiative losses and cut-offs of energetic particles at relativistic shocks
Authors:
Paul Dempsey,
Peter Duffy
Abstract:
We investigate the acceleration and simultaneous radiative losses of electrons in the vicinity of relativistic shocks. Particles undergo pitch angle diffusion, gaining energy as they cross the shock by the Fermi mechanism and also emitting synchrotron radiation in the ambient magnetic field. A semi-analytic approach is developed which allows us to consider the behaviour of the shape of the spect…
▽ More
We investigate the acceleration and simultaneous radiative losses of electrons in the vicinity of relativistic shocks. Particles undergo pitch angle diffusion, gaining energy as they cross the shock by the Fermi mechanism and also emitting synchrotron radiation in the ambient magnetic field. A semi-analytic approach is developed which allows us to consider the behaviour of the shape of the spectral cut-off and the variation of that cut-off with the particle pitch angle. The implications for the synchrotron emission of relativistic jets, such as those in gamma ray burst sources and blazars, are discussed.
△ Less
Submitted 2 April, 2007;
originally announced April 2007.
-
Particle Acceleration at Multiple Internal Relativistic Shocks
Authors:
Paul Dempsey,
Peter Duffy
Abstract:
Relativistic shocks provide an efficient method for high-energy particle acceleration in many astrophysical sources. Multiple shock systems are even more effective and of importance, for example, in the internal shock model of gamma-ray bursts. We investigate the reacceleration of pre-existing energetic particles at such relativistic internal shocks by the first order Fermi process of pitch angl…
▽ More
Relativistic shocks provide an efficient method for high-energy particle acceleration in many astrophysical sources. Multiple shock systems are even more effective and of importance, for example, in the internal shock model of gamma-ray bursts. We investigate the reacceleration of pre-existing energetic particles at such relativistic internal shocks by the first order Fermi process of pitch angle scattering. We use a well established eigenfunction method to calculate the resulting spectra for infinitely thin shocks. Implications for GRBs and relativistic jets are discussed.
△ Less
Submitted 3 April, 2007; v1 submitted 9 October, 2006;
originally announced October 2006.
-
A microscopic analysis of shear acceleration
Authors:
Frank M. Rieger,
Peter Duffy
Abstract:
A microscopic analysis of the viscous energy gain of energetic particles in (gradual) non-relativistic shear flows is presented. We extend previous work and derive the Fokker-Planck coefficients for the average rate of momentum change and dispersion in the general case of a momentum-dependent scattering time $τ(p) \propto p^α$ with $α\geq 0$. We show that in contrast to diffusive shock accelerat…
▽ More
A microscopic analysis of the viscous energy gain of energetic particles in (gradual) non-relativistic shear flows is presented. We extend previous work and derive the Fokker-Planck coefficients for the average rate of momentum change and dispersion in the general case of a momentum-dependent scattering time $τ(p) \propto p^α$ with $α\geq 0$. We show that in contrast to diffusive shock acceleration the characteristic shear acceleration timescale depends inversely on the particle mean free path which makes the mechanism particularly attractive for high energy seed particles. Based on an analysis of the associated Fokker-Planck equation we show that above the injection momentum $p_0$ power-law differential particle number density spectra $n(p) \propto p^{-(1+ α)}$ are generated for $α>0$ if radiative energy losses are negligible. We discuss the modifications introduced by synchrotron losses and determine the contribution of the accelerated particles to the viscosity of the background flow. Possible implications for the plasma composition in mildly relativistic extragalactic jet sources (WATs) are addressed.
△ Less
Submitted 6 October, 2006;
originally announced October 2006.
-
Fermi acceleration in astrophysical jets
Authors:
Frank M. Rieger,
Valenti Bosch-Ramon,
Peter Duffy
Abstract:
We consider the acceleration of energetic particles by Fermi processes (i.e., diffusive shock acceleration, second order Fermi acceleration, and gradual shear acceleration) in relativistic astrophysical jets, with particular attention given to recent progress in the field of viscous shear acceleration. We analyze the associated acceleration timescales and the resulting particle distributions, an…
▽ More
We consider the acceleration of energetic particles by Fermi processes (i.e., diffusive shock acceleration, second order Fermi acceleration, and gradual shear acceleration) in relativistic astrophysical jets, with particular attention given to recent progress in the field of viscous shear acceleration. We analyze the associated acceleration timescales and the resulting particle distributions, and discuss the relevance of these processes for the acceleration of charged particles in the jets of AGNs, GRBs and microquasars, showing that multi-component powerlaw-type particle distributions are likely to occur.
△ Less
Submitted 5 October, 2006;
originally announced October 2006.
-
A current driven instability in parallel, relativistic shocks
Authors:
B. Reville,
J. G. Kirk,
P. Duffy
Abstract:
Recently, Bell has reanalysed the problem of wave excitation by cosmic rays propagating in the pre-cursor region of a supernova remnant shock front. He pointed out a strong, non-resonant, current-driven instability that had been overlooked in the kinetic treatments, and suggested that it is responsible for substantial amplification of the ambient magnetic field. Magnetic field amplification is a…
▽ More
Recently, Bell has reanalysed the problem of wave excitation by cosmic rays propagating in the pre-cursor region of a supernova remnant shock front. He pointed out a strong, non-resonant, current-driven instability that had been overlooked in the kinetic treatments, and suggested that it is responsible for substantial amplification of the ambient magnetic field. Magnetic field amplification is also an important issue in the problem of the formation and structure of relativistic shock fronts, particularly in relation to models of gamma-ray bursts. We have therefore generalised the linear analysis to apply to this case, assuming a relativistic background plasma and a monoenergetic, unidirectional incoming proton beam. We find essentially the same non-resonant instability noticed by Bell, and show that also under GRB conditions, it grows much faster than the resonant waves. We quantify the extent to which thermal effects in the background plasma limit the maximum growth rate.
△ Less
Submitted 2 October, 2006; v1 submitted 22 August, 2006;
originally announced August 2006.
-
Particle Acceleration in Gamma-Ray Burst Jets
Authors:
Frank M. Rieger,
Peter Duffy
Abstract:
Gradual shear acceleration of energetic particles in gamma-ray burst (GRB) jets is considered. Special emphasis is given to the analysis of universal structured jets, and characteristic acceleration timescales are determined for a power-law and a Gaussian evolution of the bulk flow Lorentz factor $γ_b$ with angle $φ$ from the jet axis. The results suggest that local power-law particle distributi…
▽ More
Gradual shear acceleration of energetic particles in gamma-ray burst (GRB) jets is considered. Special emphasis is given to the analysis of universal structured jets, and characteristic acceleration timescales are determined for a power-law and a Gaussian evolution of the bulk flow Lorentz factor $γ_b$ with angle $φ$ from the jet axis. The results suggest that local power-law particle distributions may be generated and that higher energy particles are generally concentrated closer to the jet axis. Taking several constraints into account we show that efficient electron acceleration in gradual shear flows, with maximum particle energy successively decreasing with time, may be possible on scales larger than $r \sim 10^{15}$ cm, provided the jet magnetic field becomes sufficiently weak and/or decreases rapidly enough with distance, while efficient acceleration of protons to ultra-high energies $> 10^{20}$ eV may be possible under a wide range of conditions.
△ Less
Submitted 2 November, 2005; v1 submitted 16 September, 2005;
originally announced September 2005.
-
Particle acceleration timescales in relativistic shear flows
Authors:
Frank M. Rieger,
Peter Duffy
Abstract:
We review the acceleration of energetic particles in relativistic astrophysical jets characterized by a significant velocity shear. The possible formation of power-law momentum spectra is discussed and typical acceleration timescales are determined for a variety of different conditions such as parallel and azimuthal shear flows. Special attendance is given to the analysis of parallel shear flows…
▽ More
We review the acceleration of energetic particles in relativistic astrophysical jets characterized by a significant velocity shear. The possible formation of power-law momentum spectra is discussed and typical acceleration timescales are determined for a variety of different conditions such as parallel and azimuthal shear flows. Special attendance is given to the analysis of parallel shear flows with either a linear decreasing or a Gaussian-type velocity profile. It is shown that in the presence of a gradual shear flow and a particle mean free path scaling with the gyroradius, synchrotron radiation losses may no longer be able to stop the acceleration once it has started to work efficiently. Finally, the relevance of shear acceleration in small- and large-scale relativistic jets is addressed.
△ Less
Submitted 29 June, 2005;
originally announced June 2005.
-
Particle acceleration efficiencies in astrophysical shear flows
Authors:
F. M. Rieger,
P. Duffy
Abstract:
The acceleration of energetic particles in astrophysical shear flows is analyzed. We show that in the presence of a non-relativistic gradual velocity shear, power law particle momentum distributions $f(p) \propto p^{-(3+α)}$ may be generated, assuming a momentum-dependent scattering time $τ\propto p^α$, with $α> 0$. We consider possible acceleration sites in astrophysical jets and study the cond…
▽ More
The acceleration of energetic particles in astrophysical shear flows is analyzed. We show that in the presence of a non-relativistic gradual velocity shear, power law particle momentum distributions $f(p) \propto p^{-(3+α)}$ may be generated, assuming a momentum-dependent scattering time $τ\propto p^α$, with $α> 0$. We consider possible acceleration sites in astrophysical jets and study the conditions for efficient acceleration. It is shown, for example, that in the presence of a gradual shear flow and a gyro-dependent particle mean free path, synchrotron radiation losses no longer stop the acceleration once it has started to work efficiently. This suggests that shear acceleration may naturally account for a second, non-thermal population of energetic particles in addition to a shock-accelerated one. The possible relevance of shear acceleration is briefly discussed with reference to the relativistic jet in the quasar 3C 273.
△ Less
Submitted 4 February, 2005;
originally announced February 2005.
-
Particle acceleration in astrophysical shear flows
Authors:
Frank M. Rieger,
Peter Duffy
Abstract:
We consider the acceleration of particles due to a velocity shear in relativistic astrophysical flows. The basic physical picture and the formation of power law momentum spectra is discussed for a non-relativistic velocity field using a microscopic approach. We identify possible sites for shear acceleration in relativistic astrophysical jets and analyze their associated acceleration timescales.…
▽ More
We consider the acceleration of particles due to a velocity shear in relativistic astrophysical flows. The basic physical picture and the formation of power law momentum spectra is discussed for a non-relativistic velocity field using a microscopic approach. We identify possible sites for shear acceleration in relativistic astrophysical jets and analyze their associated acceleration timescales. It is shown in particular that for a mean scattering time $τ$ scaling with the gyro-radius, the acceleration timescale for gradual shear scales in the same manner as the synchrotron cooling timescale, so that losses may no longer be able to stop the acceleration once it has started to work efficiently. Finally, the possible role of shear acceleration is discussed with reference to the relativistic jet in the quasar 3C~273.
△ Less
Submitted 10 January, 2005;
originally announced January 2005.
-
Shear acceleration in relativistic astrophysical jets
Authors:
Frank M. Rieger,
Peter Duffy
Abstract:
We consider the acceleration of energetic particles by a velocity shear in the relativistic background flow containing scattering centers. Three possible acceleration sites for astrophysical jets are identified: (1) gradual velocity shear parallel to the jet axis, such as a velocity profile decreasing linearly outward with radial coordinate, (2) gradual velocity shear perpendicular to the jet ax…
▽ More
We consider the acceleration of energetic particles by a velocity shear in the relativistic background flow containing scattering centers. Three possible acceleration sites for astrophysical jets are identified: (1) gradual velocity shear parallel to the jet axis, such as a velocity profile decreasing linearly outward with radial coordinate, (2) gradual velocity shear perpendicular to the jet axis, such as intrinsic jet rotation, and (3) non-gradual/discontinuous, longitudinal velocity shear at the jet-side boundary. We determine the characteristic acceleration timescales, specify the conditions for efficient acceleration and discuss observational features with respect to each process. In particular, it is shown that in the case of (2) the higher energy emission is expected to be concentrated closer to the jet axis, while in the case of (1) and (3) the higher energy particles are likely to be located near the edges of the jet, thus possibly leading to some form of limb-brightening.
△ Less
Submitted 10 January, 2005; v1 submitted 11 October, 2004;
originally announced October 2004.
-
Jet evolution, flux ratios and light-travel time effects
Authors:
James C. A. Miller-Jones,
Katherine M. Blundell,
Peter Duffy
Abstract:
Studies of the knotty jets in both quasars and microquasars frequently make use of the ratio of the intensities of corresponding knots on opposite sides of the nucleus in order to infer the product of the intrinsic jet speed (beta) and the cosine of the inclination angle of the jet-axis (cos{theta}), via the formalism
I_{a}/I_{r} = ((1+beta cos{theta})/(1-beta cos{theta}))^{3+alpha}, where alph…
▽ More
Studies of the knotty jets in both quasars and microquasars frequently make use of the ratio of the intensities of corresponding knots on opposite sides of the nucleus in order to infer the product of the intrinsic jet speed (beta) and the cosine of the inclination angle of the jet-axis (cos{theta}), via the formalism
I_{a}/I_{r} = ((1+beta cos{theta})/(1-beta cos{theta}))^{3+alpha}, where alpha relates the intensity I_{nu} as a function of frequency nu as I_{nu} propto nu^{-alpha}. Where beta cos{theta} is determined independently, the intensity ratio of a given pair of jet to counter-jet knots is over-predicted by the above formalism compared with the intensity ratio actually measured from radio images. As an example in the case of Cygnus X-3 the original formalism predicts an intensity ratio of about 185, whereas the observed intensity ratio at one single epoch is about 3. Mirabel and Rodriguez (1999) have refined the original formalism, and suggested measuring the intensity ratio of knots when they are at equal angular separations from the nucleus. This method is only applicable where there is sufficient time-sampling with sufficient physical resolution to interpolate the intensities of the knots at equal distances from the nucleus, and can therefore be difficult to apply to microquasars and is impossible to apply to quasars. Accounting for both the light-travel time between the knots and the simple evolution of the knots themselves reconciles this over-prediction and renders the original formalism obsolete.
△ Less
Submitted 7 January, 2004;
originally announced January 2004.
-
Time-sequenced Multi-Radio-Frequency Observations of Cygnus X-3 in Flare
Authors:
James C. A. Miller-Jones,
Katherine M. Blundell,
Michael P. Rupen,
Amy J. Mioduszewski,
Peter Duffy,
Anthony J. Beasley
Abstract:
Multifrequency observations from the VLA, VLBA and OVRO Millimeter Array of a major radio outburst of Cygnus X-3 in 2001 September are presented, measuring the evolution of the spectrum of the source over three decades in frequency, over a period of six days. Following the peak of the flare, as the intensity declines the high-frequency spectrum at frequency nu steepens from nu^{-0.4} to nu^{-0.6…
▽ More
Multifrequency observations from the VLA, VLBA and OVRO Millimeter Array of a major radio outburst of Cygnus X-3 in 2001 September are presented, measuring the evolution of the spectrum of the source over three decades in frequency, over a period of six days. Following the peak of the flare, as the intensity declines the high-frequency spectrum at frequency nu steepens from nu^{-0.4} to nu^{-0.6}, after which the spectral index remains at this latter terminal value; a trend previously observed but hitherto not satisfactorily explained. VLBA observations, for the first time, track over several days the expansion of a sequence of knots whose initial diameters are approximately 8 milliarcseconds. The light-crossing time within these plasmons is of the same order as the time-scale over which the spectrum is observed to evolve. We contend that properly accounting for light-travel time effects in and between plasmons which are initially optically thick, but which after expansion become optically thin, explains the key features of the spectral evolution, for example the observed timescale. Using the VLBA images, we have directly measured for the first time the proper motions of individual knots, analysis of which shows a two-sided jet whose axis is precessing. The best-fit jet speed is roughly beta = 0.63 and the precession period is about 5 days, significantly lower than fitted for a previous flare. Extrapolation of the positions of the knots measured by the VLBA back to zero-separation shows this to occur approximately 2.5 days after the detection of the rise in flux density of Cygnus X-3.
△ Less
Submitted 12 November, 2003;
originally announced November 2003.
-
Relativistic Blastwaves and Synchrotron Emission
Authors:
T. P. Downes,
P. Duffy,
S. Komissarov
Abstract:
Relativistic shocks accelerate particles by the first order Fermi mechanism. These particles then emit synchrotron emission in the post shock gas. We have developed a numerical code which integrates the relativistic Euler equations for fluid dynamics with a general equation of state, together with the Liouville equation for the accelerated particles. We present tests of this code and, in additio…
▽ More
Relativistic shocks accelerate particles by the first order Fermi mechanism. These particles then emit synchrotron emission in the post shock gas. We have developed a numerical code which integrates the relativistic Euler equations for fluid dynamics with a general equation of state, together with the Liouville equation for the accelerated particles. We present tests of this code and, in addition, we use it to study the gamma ray burst afterglow predicted by the fireball model, along with the hydrodynamics of a relativistic blastwave.
We find that, while, broadly speaking, the behaviour of the emission is similar to that already predicted with semi-analytic approaches, the detailed behaviour is somewhat different. The ``breaks'' in the synchrotron spectrum behave differently with time, and the spectrum above the final break is harder than previously expected. These effects are due to the incorporation of the geometry of the (spherical) blastwave, along with relativistic beaming and adiabatic cooling of the energetic particles leading to a mix, in the observed spectrum, between recently injected "uncooled" particles and the older "cooled" population in different parts of the evolving, inhomogeneous flow.
△ Less
Submitted 22 January, 2002; v1 submitted 15 February, 2001;
originally announced February 2001.
-
On "box" models of shock acceleration and electron synchrotron spectra
Authors:
Luke O'C Drury,
Peter Duffy,
David Eichler,
Apostolos Mastichiadis
Abstract:
The recent detection of high energy gamma-rays coming from supernova remnants and active galactic nuclei has revived interest in the diffusive shock acceleration of electrons. In the present paper we examine the basis of the so-called ``box'' model for particle acceleration and present a more physical version of it. Using this we determine simple criteria for the conditions under which ``pile-up…
▽ More
The recent detection of high energy gamma-rays coming from supernova remnants and active galactic nuclei has revived interest in the diffusive shock acceleration of electrons. In the present paper we examine the basis of the so-called ``box'' model for particle acceleration and present a more physical version of it. Using this we determine simple criteria for the conditions under which ``pile-ups'' can occur in shock accelerated electron spectra subject to synchrotron or inverse Compton losses (the latter in the Thompson limit). An extension to include nonlinear effects is proposed.
△ Less
Submitted 14 May, 1999;
originally announced May 1999.
-
Particle acceleration and relativistic shocks
Authors:
J. G. Kirk,
P. Duffy
Abstract:
Observations of both gamma-ray burst sources and certain classes of active galaxy indicate the presence of relativistic shock waves and require the production of high energy particles to explain their emission. In this paper we review the basic theory of shock waves in relativistic hydrodynamics and magneto-hydrodynamics, emphasising the astrophysically interesting cases. This is followed by an…
▽ More
Observations of both gamma-ray burst sources and certain classes of active galaxy indicate the presence of relativistic shock waves and require the production of high energy particles to explain their emission. In this paper we review the basic theory of shock waves in relativistic hydrodynamics and magneto-hydrodynamics, emphasising the astrophysically interesting cases. This is followed by an overview of the theory of particle acceleration at such shocks. We summarise the applications to the astrophysics of relativistic jets and fireball models of gamma-ray-bursts.
△ Less
Submitted 6 May, 1999;
originally announced May 1999.
-
Stochastic particle acceleration at shocks in the presence of braided magnetic fields
Authors:
J. G. Kirk,
P. Duffy,
Y. A. Gallant
Abstract:
The theory of diffusive acceleration of energetic particles at shock fronts assumes charged particles undergo spatial diffusion in a uniform magnetic field. If, however, the magnetic field is not uniform, but has a stochastic or braided structure, the transport of charged particles across the average direction of the field is more complicated. Assuming quasi-linear behaviour of the field lines,…
▽ More
The theory of diffusive acceleration of energetic particles at shock fronts assumes charged particles undergo spatial diffusion in a uniform magnetic field. If, however, the magnetic field is not uniform, but has a stochastic or braided structure, the transport of charged particles across the average direction of the field is more complicated. Assuming quasi-linear behaviour of the field lines, the particles undergo sub-diffusion on short time scales. We derive the propagator for such motion, which differs from the Gaussian form relevant for diffusion, and apply it to a configuration with a plane shock front whose normal is perpendicular to the average field direction. Expressions are given for the acceleration time as a function of the diffusion coefficient of the wandering magnetic field lines and the spatial diffusion coefficient of the charged particles parallel to the local field. In addition we calculate the spatial dependence of the particle density in both the upstream and downstream plasmas. In contrast to the diffusive case, the density of particles at the shock front is lower than it is far downstream. This is a consequence of the partial trapping of particles by structures in the magnetic field. As a result, the spectrum of accelerated particles is a power-law in momentum which is steeper than in the diffusive case. For a phase- space density $f\propto p^{-s}$, we find $s=\sdiff[1+1/(2\rcomp)]$, where $\rcomp$ is the compression ratio of the shock front and $\sdiff$ is the standard result of diffusive acceleration: $\sdiff=3\rcomp/(\rcomp-1)$. A strong shock in a monatomic ideal gas yields a spectrum of $s=4.5$. In the case of electrons, this corresponds to a radio synchrotron spectral index of $α=0.75$.
△ Less
Submitted 11 April, 1996;
originally announced April 1996.
-
Limits on diffusive shock acceleration in dense and incompletely ionised media
Authors:
L O'C Drury,
P Duffy,
J G Kirk
Abstract:
The limits imposed on diffusive shock acceleration by upstream ion-neutral Alfven wave damping, and by ionisation and Coulomb losses of low energy particles, are calculated. Analytic solutions are given for the steady upstream wave excitation problem with ion-neutral damping and the resulting escaping upstream flux calculated. The time dependent problem is discussed and numerical solutions prese…
▽ More
The limits imposed on diffusive shock acceleration by upstream ion-neutral Alfven wave damping, and by ionisation and Coulomb losses of low energy particles, are calculated. Analytic solutions are given for the steady upstream wave excitation problem with ion-neutral damping and the resulting escaping upstream flux calculated. The time dependent problem is discussed and numerical solutions presented. Finally the significance of these results for possible observational tests of shock acceleration in supernova remnants is discussed.
△ Less
Submitted 31 January, 1996; v1 submitted 12 October, 1995;
originally announced October 1995.
-
Anomalous Transport and Particle Acceleration at Shocks
Authors:
P. Duffy,
J. G. Kirk,
Y. A. Gallant,
R. O. Dendy
Abstract:
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-Ro…
▽ More
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.
△ Less
Submitted 11 September, 1995;
originally announced September 1995.
-
A MODEL FOR THE RADIO EMISSION FROM SNR 1987A
Authors:
P. Duffy,
Lewis Ball,
J. G. Kirk
Abstract:
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 under…
▽ More
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.
△ Less
Submitted 11 January, 1995;
originally announced January 1995.
-
Radio Supernovae as Tev Gamma-Ray Sources
Authors:
J. G. Kirk,
P. Duffy,
Lewis Ball
Abstract:
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…
▽ More
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
△ Less
Submitted 21 November, 1994;
originally announced November 1994.
-
Cosmic Ray Hydrodynamics at Shock Fronts
Authors:
P. Duffy,
L. O'C. Drury,
H. J. Voelk
Abstract:
When a magnetised, thermal plasma containing a shock wave efficiently accelerates Cosmic Rays (CR) a reaction is exerted on the fluid structure. The simplest description of this process involves three conservation equations for the background fluid, with the CR pressure gradient included as an extra force, plus a hydrodynamic equation for the CR energy density. Since the CR spectrum is not deter…
▽ More
When a magnetised, thermal plasma containing a shock wave efficiently accelerates Cosmic Rays (CR) a reaction is exerted on the fluid structure. The simplest description of this process involves three conservation equations for the background fluid, with the CR pressure gradient included as an extra force, plus a hydrodynamic equation for the CR energy density. Since the CR spectrum is not determined in this approach two quantities are introduced, the ratio of CR pressure to energy density and a mean diffusion coefficient, that must be specified to close the system of equations. Physical arguments are presented for how these quantities should evolve and then compared with time dependent numerical solutions of acceleration at a plane parallel, piston driven shock.
△ Less
Submitted 3 August, 1994;
originally announced August 1994.