COMPUTATION OF SYNTHETIC SPECTRA FROM SIMULATIONS OF RELATIVISTIC SHOCKS

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

RADIATIVE SIGNATURES OF RELATIVISTIC SHOCKS

Particle-in-cell simulations of relativistic, weakly magnetized collisionless shocks show that particles can gain energy by repeatedly crossing the shock front. This requires scattering off self-generated small length-scale magnetic fluctuations. The radiative signature of this first-order Fermi acceleration mechanism is important for models of both the prompt and afterglow emission in gamma-ray bursts and depends on the strength parameter a = lambda e/delta B/mc(2) of the fluctuations (lambda is the length scale and vertical bar delta B vertical bar is the magnitude of the fluctuations). For electrons (and positrons), acceleration saturates when the radiative losses produced by the scattering cannot be compensated by the energy gained on crossing the shock. We show that this sets an upper limit on both the electron Lorentz factor gamma <10(6) (n/1 cm(-3))(-1/6)(-1/6) and on the energy of the photons radiated during the scattering process h omega(max) <40Max(a, 1)(n/1 cm(-3))(1/6)(-1/6) eV, where n is the number density of the plasma and (gamma) over bar is the thermal Lorentz factor of the downstream plasma, provided a <a(crit) similar to 10(6). This rules out "jitter" radiation on self-excited fluctuations with a <I as a source of gamma rays, although high-energy photons might still be produced when the jitter photons are upscattered in an analog of the synchrotron self-Compton process. In fluctuations with a > 1, radiation is generated by the standard synchrotron mechanism, and the maximum photon energy rises linearly with a, until saturating at 70 MeV, when a = a(crit).

TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF THE NONRESONANT, COSMIC-RAY-DRIVEN INSTABILITY IN SUPERNOVA REMNANT SHOCKS

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

A current-driven instability in parallel, relativistic shocks

Recently, Bell ( 2004 Mon. Not. R. Astron. Soc. 353 550) 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 by Achterberg ( 1983 Astron. Astrophys. 119 274) and McKenzie and Volk ( 1982 Astron. Astrophys. 116 191), 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 generalized 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 observed 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.

On the maximum energy of shock-accelerated cosmic rays at ultra-relativistic shocks

The maximum energy to which cosmic rays can be accelerated at weakly magnetised ultra-relativistic shocks is investigated. We demonstrate that for such shocks, in which the scattering of energetic particles is mediated exclusively by ion skin-depth scale structures, as might be expected for a Weibel-mediated shock, there is an intrinsic limit on the maximum energy to which particles can be accelerated. This maximum energy is determined from the requirement that particles must be isotropized in the downstream plasma frame before the mean field transports them far downstream, and falls considerably short of what is required to produce ultra-high-energy cosmic rays. To circumvent this limit, a highly disorganized field is required on larger scales. The growth of cosmic ray-induced instabilities on wavelengths much longer than the ion-plasma skin depth, both upstream and downstream of the shock, is considered. While these instabilities may play an important role in magnetic field amplification at relativistic shocks, on scales comparable to the gyroradius of the most energetic particles, the calculated growth rates have insufficient time to modify the scattering. Since strong modification is a necessary condition for particles in the downstream region to re-cross the shock, in the absence of an alternative scattering mechanism, these results imply that acceleration to higher energies is ruled out. If weakly magnetized ultra-relativistic shocks are disfavoured as high-energy particle accelerators in general, the search for potential sources of ultra-high-energy cosmic rays can be narrowed.

Dust-acoustic shocks in strongly coupled dusty plasmas

Electrostatic dust-acoustic shock waves are investigated in a viscous, complex plasma consisting of dust particles, electrons, and ions. The system is modelled using the generalized hydrodynamic equations, with strong coupling between the dust particles being accounted for by employing the effective electrostatic temperature approach. Using a reductive perturbation method, it is demonstrated that this model predicts the existence of weakly nonlinear dust-acoustic shock waves, arising as solutions to Burgers's equation, in which the nonlinear forces are balanced by dissipative forces, in this case, associated with viscosity. The evolution and stability of dust-acoustic shocks is investigated via a series of numerical simulations, which confirms our analytical predictions on the shock characteristics.

Shocks in unmagnetized plasma with a shear flow: Stability and magnetic field generation

A pair of curved shocks in a collisionless plasma is examined with a two-dimensional particle-in-cell simulation. The shocks are created by the collision of two electron-ion clouds at a speed that exceeds everywhere the threshold speed for shock formation. A variation of the collision speed along the initially planar collision boundary, which is comparable to the ion acoustic speed, yields a curvature of the shock that increases with time. The spatially varying Mach number of the shocks results in a variation of the downstream density in the direction along the shock boundary. This variation is eventually equilibrated by the thermal diffusion of ions. The pair of shocks is stable for tens of inverse ion plasma frequencies. The angle between the mean flow velocity vector of the inflowing upstream plasma and the shock's electrostatic field increases steadily during this time. The disalignment of both vectors gives rise to a rotational electron flow, which yields the growth of magnetic field patches that are coherent over tens of electron skin depths.

Self-regulated Shocks in Massive Star Binary Systems

In an early-type, massive star binary system, X-ray bright shocks result from the powerful collision of stellar winds driven by radiation pressure on spectral line transitions. We examine the influence of the X-rays from the wind-wind collision shocks on the radiative driving of the stellar winds using steady-state models that include a parameterized line force with X-ray ionization dependence. Our primary result is that X-ray radiation from the shocks inhibits wind acceleration and can lead to a lower pre-shock velocity, and a correspondingly lower shocked plasma temperature, yet the intrinsic X-ray luminosity of the shocks, L X, remains largely unaltered, with the exception of a modest increase at small binary separations. Due to the feedback loop between the ionizing X-rays from the shocks and the wind driving, we term this scenario as self-regulated shocks. This effect is found to greatly increase the range of binary separations at which a wind-photosphere collision is likely to occur in systems where the momenta of the two winds are significantly different. Furthermore, the excessive levels of X-ray ionization close to the shocks completely suppress the line force, and we suggest that this may render radiative braking less effective. Comparisons of model results against observations reveal reasonable agreement in terms of log (L X/L bol). The inclusion of self-regulated shocks improves the match for kT values in roughly equal wind momenta systems, but there is a systematic offset for systems with unequal wind momenta (if considered to be a wind-photosphere collision).

The Impact of Uncertainty Shocks under Measurement Error: A Proxy SVAR approach

A growing literature considers the impact of uncertainty using SVAR models that include proxies for uncertainty shocks as endogenous variables. In this paper we consider the impact of measurement error in these proxies on the estimated impulse responses. We show via a Monte-Carlo experiment that measurement error can result in attenuation bias in impulse responses. In contrast, the proxy SVAR that uses the uncertainty shock proxy as an instrument does not su¤er from this bias. Applying this latter method to the Bloom (2009) data-set results in impulse responses to uncertainty shocks that are larger in magnitude and more persistent than those obtained from a recursive SVAR.

The role of investment specific technological shocks in explaining macroeconomics fluctuations.

Rapport de recherche

Habit Formation and the Persistence of Monetary Shocks

This paper studies the persistent effects of monetary shocks on output. Previous empirical literature documents this persistence, but standard general equilibrium models with sticky prices fail to generate output responses beyond the duration of nominal contracts. This paper constructs and estimates a general equilibrium model with price rigidities, habit formation, and costly capital adjustment. The model is estimated via Maximum Likelihood using US data on output, the real money stock, and the nominal interest rate. Econometric results suggest that habit formation and adjustment costs to capital play an important role in explaining the output effects of monetary policy. In particular, impulse response analysis indicates that the model generates persistent, hump-shaped output responses to monetary shocks.

A Systematic Framework for Analyzing the Dynamic Effects of Permanent and Transitory Shocks.

This paper proposes a systematic framework for analyzing the dynamic effects of permanent and transitory shocks on a system of \"n\" economic variables.

A Systematic Framework for Analyzing the Dynamic Effects of Permanent and Transitory Shocks.

This paper proposes a systematic framework for analyzing the dynamic effects of permanent and transitory shocks on a system of n economic variables.

Cross-sectoral variation in the volatility of plant-level idiosyncratic shocks

We estimate the volatility of plant–level idiosyncratic shocks in the U.S. manufacturing sector. Our measure of volatility is the variation in Revenue Total Factor Productivity which is not explained by either industry– or economy–wide factors, or by establishments’ characteristics. Consistent with previous studies, we find that idiosyncratic shocks are much larger than aggregate random disturbances, accounting for about 80% of the overall uncertainty faced by plants. The extent of cross–sectoral variation in the volatility of shocks is remarkable. Plants in the most volatile sector are subject to about six times as much idiosyncratic uncertainty as plants in the least volatile. We provide evidence suggesting that idiosyncratic risk is higher in industries where the extent of creative destruction is likely to be greater.

Short Term Health Shocks and School Attendance: The Case of a Dengue Fever Outbreak in Colombia.

This paper makes use of a short, sharp, unexpected health shock in the form of the 2010 Colombian Dengue outbreak to examine the direct and indirect impact of negative health shocks on behaviour of households in affected areas. Our analysis combines data from several sources in order to obtain a comprehensive picture of the influence of the outbreak, and furthermore to understand the underlying mechanisms driving the effects. Our initial analysis indicates that the outbreak had a substantial negative effect on the health status of adults and adversely affected their ability to function as usual in their daily lives. In our aggregated school data, in areas with high levels of haemorrhagic Dengue we observe a reduction innational exam attendance (last year of secondary school) and on enrolment rates in primary education. Further analysis aims to exploit detailed individual level data to gain a more in depth understanding of the precise channels through which this disease influenced the behaviour and outcomes of the poor in Colombia.