Intergalactic shock acceleration and the cosmic gamma-ray background

We investigate numerically the contribution to the cosmic gamma-ray background from cosmic-ray ions and electrons accelerated at intergalactic shocks associated with cosmological structure formation. We show that the kinetic energy of accretion flows in the low-redshift intergalactic medium is thermalized primarily through moderately strong shocks, which allow for an efficient conversion of shock ram pressure into cosmic-ray pressure. Cosmic rays accelerated at these shocks produce a diffuse gamma-ray flux which is dominated by inverse Compton emission from electrons scattering off cosmic microwave background photons. Decay of neutral π mesons generated in p–p inelastic collisions of the ionic cosmic-ray component with the thermal gas contribute about 30 per cent of the computed emission. Based on experimental upper limits on the photon flux above 100 MeV from nearby clusters we constrain the efficiency of conversion of shock ram pressure into relativistic CR electrons to  ≲1 per cent  . Thus, we find that cosmic rays of cosmological origin can generate an overall significant fraction of order 20 per cent and no more than 30 per cent of the measured gamma-ray background.     

On the role of injection in kinetic approaches to non-linear particle acceleration at non-relativistic shock waves

The dynamical reaction of the particles accelerated at a shock front by the first-order Fermi process can be determined within kinetic models that account for both the hydrodynamics of the shocked fluid and the transport of the accelerated particles. These models predict the appearance of multiple solutions, all physically allowed. We discuss here the role of injection in selecting the real solution, in the framework of a simple phenomenological recipe, which is a variation of what is sometimes referred to as thermal leakage. In this context we show that multiple solutions basically disappear and when they are present they are limited to rather peculiar values of the parameters. We also provide a quantitative calculation of the efficiency of particle acceleration at cosmic ray modified shocks and we identify the fraction of energy which is advected downstream and that of particles escaping the system from upstream infinity at the maximum momentum. The consequences of efficient particle acceleration for shock heating are also discussed.     

Electron surfing acceleration in oblique magnetic fields

Initially, inhomogeneous plasma jets, ejected by active galactic nuclei and associated with gamma-ray bursts, are thermalized by the formation of internal shocks. Jet subpopulations can hereby collide at Lorentz factors of a few. As the resulting relativistic shock expands into the upstream plasma, a significant fraction of the upstream ions is reflected. These ions, together with downstream ions that leak through the shock, form relativistic beams of ions that outrun the shock. The thermalization of these beams via the two-stream instability is thought to contribute significantly to plasma heating and particle acceleration by the shock. Here, the capability of a two-stream instability to generate relativistic field-aligned and cross-field electron flow, is examined for a magnetized plasma by means of a particle-in-cell (PIC) simulation. The electrons interact with the developing quasi-electrostatic waves and oblique magnetic fields. The simulation results bring forward evidence that such waves, by their non-linear interactions with the plasma, produce a highly relativistic field-aligned electron flow and electron energies, which could contribute to the radio synchrotron emissions from astrophysical jets, to ultrarelativistic leptonic subpopulations propagating with the jet and to the halo particles surrounding the accretion disc of the black hole.     

Non-linear particle acceleration at non-relativistic shock waves in the presence of self-generated turbulence

Particle acceleration at astrophysical shocks may be very efficient if magnetic scattering is self-generated by the same particles. This non-linear process adds to the non-linear modification of the shock due to the dynamical reaction of the accelerated particles on the shock. Building on a previous general solution of the problem of particle acceleration with arbitrary diffusion coefficients, we present here the first semi-analytical calculation of particle acceleration with both effects taken into account at the same time; charged particles are accelerated in the background of Alfvén waves that they generate due to the streaming instability, and modify the dynamics of the plasma in the shock vicinity.     

A multiple shooting approach for the numerical treatment of stellar structure and evolution

We present a new numerical method for solving the system of partial differential equations describing the structure and evolution of a spherically symmetric star. As usual, we employ the transversal method of lines in order to split the equations into a coupled spatial and temporal part. The novel features of the algorithm are the following: (a) Instead of using the Lagrangian picture we formulate the system of partial differential equations in the Eulerian picture. (b) We reformulate the equations of stellar structure as a multipoint boundary-value problem. By means of this reformulation the rather clumsy iterative matching procedure of stellar atmosphere and interior is avoided. (c) The multipoint boundary-value problem is solved by the multiple shooting method. This approach not only ensures a high accuracy of the stellar models calculated at each time step but also allows the free boundaries inside the star due to different energy transport mechanisms to be located exactly. (d) The time derivatives involved in the stellar-structure equations are discretized implicitly to second order accuracy. Moreover, at each time step, the chemical abundances are determined by using a sophisticated update procedure. In this way, a high accuracy is achieved with respect to the integration in time. The algorithm has turned out to be exceedingly reliable and numerically accurate. This is shown by the evolution of a 1 M_⊙ star up to the hydrogen-shell burning phase. In this example, the virial theorem, the law of mass conservation, and the law of energy conservation is fulfilled to a hitherto unattainable degree of accuracy. Since the multiple shooting method, which is at the heart of our approach, is a perfect example of a parallel algorithm, the computational speed of the algorithm might be substantially improved provided easy-to-program, high-performance parallel computers with sufficiently many processors become available in the near future.     

A fitting formula for the non-linear evolution of the bispectrum

We present a fitting formula for the non-linear evolution of the bispectrum in cold dark matter (CDM) models, obtained from measurements in high-resolution numerical simulations. The formula interpolates between the perturbative and highly non-linear regimes, and generalizes previous results obtained for scale-free initial conditions.     

The shock-reprocessing model of electron acceleration in impulsive solar flares

We propose a new two-stage model for acceleration of electrons in solar flares. In the first stage, electrons are accelerated stochastically in a post-reconnection turbulent downflow. The second stage is the reprocessing of a subset of these electrons as they pass through a weakly compressive fast shock above the apex of the closed flare loop on their way to the chromosphere. We call this the 'shock-reprocessing' model. The model reproduces the sign and magnitude of the energy-dependent arrival time delays for both the pulsed and smooth component of impulsive solar flare X-rays, but requires either enhanced cooling or the presence of a loop-top trap to explain the concavity of the observed time delay energy relation for the smooth component. The model also predicts an emission site above the loop-top, as seen in the Masuda flare. The loop-top source distinguishes the shock-reprocessing model from previous models. The model makes testable predictions for the energy dependence of footpoint pulse strengths and the location and spectrum of the loop-top emission, and can account for the observed soft-hard-soft trend in the spectral evolution of footpoint emission. The model also highlights the concept that magnetic reconnection provides an environment which permits multiple acceleration processes. Which combination of processes operates within a particular flare may depend on the initial conditions that determine, for example, whether the reconnection downflow is turbulent or laminar. The shock-reprocessing model comprises one such combination.     

Studying magneto‐convection by numerical simulation

Following a brief overview of the two main approaches to investigate the interaction between magnetic fields and convective flows near the solar surface layers by numerical simulation, namely idealized model problems and ‘realistic’ large‐eddy simulations, we present first results obtained with a newly developed MHD code. The first example concerns the realistic simulation of the magnetic field dynamics in a solar plage region while the second example demonstrates small‐scale dynamo action in idealized compressible convection.     

Cosmic ray acceleration to very high energy through the non-linear amplification by cosmic rays of the seed magnetic field

The maximum energy for cosmic ray acceleration at supernova shock fronts is usually thought to be limited to around 10¹⁴–10¹⁵ eV by the size of the shock and the time for which it propagates at high velocity. We show that the magnetic field can be amplified non-linearly by the cosmic rays to many times the pre-shock value, thus increasing the acceleration rate and facilitating acceleration to energies well above 10¹⁵ eV. A supernova remnant expanding into a uniform circumstellar medium may accelerate protons to 10¹⁷ eV and heavy ions, with charge Ze , to Z ×10¹⁷ eV. Expansion into a pre-existing stellar wind may increase the maximum cosmic ray energy by a further factor of 10.     

初始谱指数为-1的非线性重谱(Bispectrum)

重谱即傅立叶空间的三点相关函数，是研究由引力不稳定性产生的大尺度结构的最低非高斯阶统计性质．在傅立叶空间中，发现二维和三维空间的统计性质是等价的，相对三维快速傅立叶，二维快速傅立叶变换可以节省很多计算机资源．结合二维和三维的快速傅立叶变换，利用高精度的初始谱指数为-1的N体数值模拟，测量了从弱线性区域到强非线性区域的重谱．在测量过程中，考虑并严格校正了作用于功率谱和重谱上的数值效应，如离散效应（particle discreteness）、软化因子（force softening）、模拟盒子的大小、混淆效应（alias effect）．发现约化重谱即使在强非线性区域仍然依赖于三角形的形状和尺度，并随着波数的增加而增加，这和前人的假设完全不同．在测量结果的基础上，拟合了重谱的一个适用于初始谱指数为-1的宇宙学模型的经验公式．     

High-energy particle acceleration and production of ultra-high-energy cosmic rays in the giant lobes of Centaurus A

The nearby radio galaxy Centaurus A is poorly studied at high frequencies with conventional radio telescopes because of its very large angular size, but is one of a very few extragalactic objects to be detected and resolved by the Wilkinson Microwave Anisotropy Probe ( WMAP ). We have used the five-year WMAP data for Cen A to constrain the high-frequency radio spectra of the 10° giant lobes and to search for spectral changes as a function of position along the lobes. We show that the high-frequency radio spectra of the northern and southern giant lobes are significantly different: the spectrum of the southern lobe steepens monotonically (and is steeper further from the active nucleus) whereas the spectrum of the northern lobe remains consistent with a power law. The inferred differences in the northern and southern giant lobes may be the result of real differences in their high-energy particle acceleration histories, perhaps due to the influence of the northern middle lobe, an intermediate-scale feature which has no detectable southern counterpart. In light of these results, we discuss the prospects for Fermi Gamma-ray Space Telescope detections of inverse-Compton emission from the giant lobes and the lobes' possible role in the production of the ultra-high-energy cosmic rays (UHECR) detected by the Pierre Auger Observatory. We show that the possibility of a Fermi detection depends sensitively on the physical conditions in the giant lobes, with the northern lobe more likely to be detected, and that any emission observed by Fermi is likely to be dominated by photons at the soft end of the Fermi energy band. On the other hand, we argue that the estimated conditions in the giant lobes imply that UHECRs can be accelerated there, with a potentially detectable γ-ray signature at TeV energies.