The magnetic structures of electron phase-space holes formed in the electron two-stream instability

It is well known that the parallel cuts of the parallel and perpendicular electric field in electron phase-space holes (electron holes) have bipolar and unipolar structures, respectively. Recently, electron holes in the Earth’s plasma sheet have been observed by THEMIS satellites to have detectable fluctuating magnetic field with regular structures. Du et al. (2011) investigated the evolution of a one-dimensional (1D) electron hole with two-dimensional (2D) electromagnetic particle-in-cell (PIC) simulations in weakly magnetized plasma (Ω _e <ω _pe , where Ω _e and ω _pe are the electron gyrofrequency and electron plasma frequency, respectively), which initially exists in the simulation domain. The electron hole is unstable to the transverse instability and broken into several 2D electron holes. They successfully explained the observations by THEMIS satellites based on the generated magnetic structures associated with these 2D electron holes. In this paper, 2D electromagnetic particle-in-cell (PIC) simulations are performed in the x–y plane to investigate the nonlinear evolution of the electron two-stream instability in weakly magnetized plasma, where the background magnetic field (B₀ = B₀[(e)\vec] _x)(\mathbf{B}_{0} =B_{0}\vec{\mathbf{e}} _{x}) is along the x direction. Several 2D electron holes are formed during the nonlinear evolution, where the parallel cuts of E _x and E _y have bipolar and unipolar structures, respectively. Consistent with the results of Du et al. (2011), we found that the current along the z direction is generated by the electric field drift motion of the trapped electrons in the electron holes due to the existence of E _y , which produces the fluctuating magnetic field δB _x and δB _y in the electron holes. The parallel cuts of δB _x and δB _y in the electron holes have unipolar and bipolar structures, respectively.     

Simulation studies of electron acceleration by ion ring distributions in solar flares

Using a 2 1/2-D fully relativistic electromagnetic particle-in-cell code (PIC) we have investigated a potential electron acceleration mechanism in solar flares. The free energy is provided by ions which have a ring velocity distribution about the magnetic field direction. Ion rings may be produced by perpendicular shocks, which could in turn be generated by the super-Alfvénic motion of magnetic flux tubes emerging from the photosphere or by coronal mass ejections (CMEs). Such ion distributions are known to be unstable to the generation of lower hybrid waves, which have phase velocities in excess of the electron thermal speed parallel to the field and can, therefore, resonantly accelerate electrons in that direction. The simulations show the transfer of perpendicular ion energy to energetic electrons via lower hybrid wave turbulence. With plausible ion ring velocities, the process can account for the observationally inferred fluxes and energies of non-thermal electrons during the impulsive phase of flares. Our results also show electrostatic wave generation close to the plasma frequency: we suggest that this is due to a bump-in-tail instability of the electron distribution.     

Cross-scale coupling at a perpendicular collisionless shock

A full particle simulation study is carried out on a perpendicular collisionless shock with a relatively low Alfven Mach number (M_A = 5). Recent self-consistent hybrid and full particle simulations have demonstrated ion kinetics are essential for the non-stationarity of perpendicular collisionless shocks, which means that physical processes due to ion kinetics modify the shock jump condition for fluid plasmas. This is a cross-scale coupling between fluid dynamics and ion kinetics. On the other hand, it is not easy to study cross-scale coupling of electron kinetics with ion kinetics or fluid dynamics, because it is a heavy task to conduct large-scale full particle simulations of collisionless shocks. In the present study, we have performed a two-dimensional (2D) electromagnetic full particle simulation with a “shock-rest-frame model”. The simulation domain is taken to be larger than the ion inertial length in order to include full kinetics of both electrons and ions. The present simulation result has confirmed the transition of shock structures from the cyclic self-reformation to the quasi-stationary shock front. During the transition, electrons and ions are thermalized in the direction parallel to the shock magnetic field. Ions are thermalized by low-frequency electromagnetic waves (or rippled structures) excited by strong ion temperature anisotropy at the shock foot, while electrons are thermalized by high-frequency electromagnetic waves (or whistler mode waves) excited by electron temperature anisotropy at the shock overshoot. Ion acoustic waves are also excited at the shock overshoot where the electron parallel temperature becomes higher than the ion parallel temperature. We expect that ion acoustic waves are responsible for parallel diffusion of both electrons and ions, and that a cross-scale coupling between an ion-scale mesoscopic instability and an electron-scale microscopic instability is important for structures and dynamics of a collisionless perpendicular shock.     

3-D Rpic Simulations of Relativistic Jets: Particle Acceleration, Magnetic Field Generation, and Emission

We have applied numerical simulations and modeling to the particle acceleration, magnetic field generation, and emission from relativistic shocks. We investigate the nonlinear stage of theWeibel instability and compare our simulations with the observed gamma-ray burst emission. In collisionless shocks, plasma waves and their associated instabilities (e.g., the Weibel, Buneman and other two-stream instabilities) are responsible for particle (electron, positron, and ion) acceleration and magnetic field generation. 3-D relativistic electromagnetic particle (REMP) simulations with three different electron-positron jet velocity distributions and also with an electron-ion plasma have been performed and show shock processes including spatial and temporal evolution of shocks in unmagnetized ambient plasmas. The growth time and nonlinear saturation levels depend on the initial jet parallel velocity distributions. Simulations show that the Weibel instability created in the collisionless shocks accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The nonlinear fluctuation amplitude of densities, currents, electric, and magnetic fields in the electron-positron shocks are larger for smaller jet Lorentz factor. This comes from the fact that the growth time of the Weibel instability is proportional to the square of the jet Lorentz factor. We have performed simulations with broad Lorentz factor distribution of jet electrons and positrons, which is assumed to be created by photon annihilation. Simulation results with this broad distribution show that the Weibel instability is excited continuously by the wide-range of jet Lorentz factor from lower to higher values. In all simulations the Weibel instability is responsible for generating and amplifying magnetic fields perpendicular to the jet propagation direction, and contributes to the electron’s (positron’s) transverse deflection behind the jet head. This small scale magnetic field structure contributes to the generation of “jitter” radiation from deflected electrons (positrons), which is different from synchrotron radiation in uniform magnetic fields. The jitter radiation resulting from small scale magnetic field structures may be important for understanding the complex time structure and spectral evolution observed in gamma-ray bursts or other astrophysical sources containing relativistic jets and relativistic collisionless shocks. The detailed studies of shock microscopic process evolution may provide some insights into early and later GRB afterglows.     

The electron–cyclotron maser for astrophysical application

The electron–cyclotron maser is a process that generates coherent radiation from plasma. In the last two decades, it has gained increasing attention as a dominant mechanism of producing high-power radiation in natural high-temperature magnetized plasmas. Originally proposed as a somewhat exotic idea and subsequently applied to include non-relativistic plasmas, the electron–cyclotron maser was considered as an alternative to turbulent though coherent wave–wave interaction which results in radio emission. However, when it was recognized that weak relativistic corrections had to be taken into account in the radiation process, the importance of the electron–cyclotron maser rose to the recognition it deserves. Here we review the theory and application of the electron–cyclotron maser to the directly accessible plasmas in our immediate terrestrial and planetary environments. In situ access to the radiating plasmas has turned out to be crucial in identifying the conditions under which the electron–cyclotron maser mechanism is working. Under extreme astrophysical conditions, radiation from plasmas may provide a major energy loss; however, for generating the powerful radiation in which the electron–cyclotron maser mechanism is capable, the plasma must be in a state where release of susceptible amounts of energy in the form of radiation is favorable. Such conditions are realized when the plasma is unable to digest the available free energy that is imposed from outside and stored in its particle distribution. The lack of dissipative processes is a common property of collisionless plasmas. When, in addition, the plasma density becomes so low that the amount of free energy per particle is large, direct emission becomes favorable. This can be expressed as negative absorption of the plasma which, like in conventional masers, leads to coherent emission even though no quantum correlations are involved. The physical basis of this formal analogy between a quantum maser and the electron–cyclotron maser is that in the electron–cyclotron maser the free-space radiation modes can be amplified directly. Several models have been proposed for such a process. The most famous one is the so-called loss-cone maser. However, as argued in this review, the loss-cone maser is rather inefficient. Available in situ measurements indicate that the loss-cone maser plays only a minor role. Instead, the main source for any strong electron–cyclotron maser is found in the presence of a magnetic-field-aligned electric potential drop which has several effects: (1) it dilutes the local plasma to such an extent that the plasma enters the regime in which the electron–cyclotron maser becomes effective; (2) it generates energetic relativistic electron beams and field-aligned currents; (3) it deforms, together with the magnetic mirror force, the electron distribution function, thereby mimicking a high energy level sufficiently far above the Maxwellian ground state of an equilibrium plasma; (4) it favors emission in the free-space RX mode in a direction roughly perpendicular to the ambient magnetic field; (5) this emission is the most intense, since it implies the coherent resonant contribution of a maximum number of electrons in the distribution function to the radiation (i.e., to the generation of negative absorption); (6) it generates a large number of electron holes via the two-stream instability, and ion holes via the current-driven ion-acoustic instability which manifest themselves as subtle fine structures moving across the radiation spectrum and being typical for the electron–cyclotron maser emission process. These fine structures can thus be taken as the ultimate identifier of the electron–cyclotron maser. The auroral kilometric radiation of Earth is taken here as the paradigm for other manifestations of intense radio emissions such as the radiation from other planets in the solar system, from exoplanets, the Sun and other astrophysical objects.     

Numerical simulation of electromagnetic emissions in the solar wind plasma with non-thermal electron distribution

Dynamics of fundamental and second harmonic electromagnetic emissions are simulated in the solar wind plasma in the presence of non-thermal electron distribution function in which primary Langmuir waves are driven by an electron beam. The electron velocity distribution function is separated into two distributions representing the distribution of the ambient electrons (Maxwellian) and the suprathermal electrons (non-thermal electrons). The effects of the non-thermal electrons on the generation of primary Langmuir waves, emission rates of the fundamental (F) and harmonic waves (H) and their distributions are investigated. The both of the F and H emissions are sensitive to the characterizes of the non-thermal electrons. It is found that in the presence of non-thermal electrons the production of the Langmuir waves decreases and consequently the levels of fundamental and second harmonic waves are reduced. The emission rate of the fundamental transverse waves decreases and its peak moves slightly toward smaller wave-numbers.     

Observational Evidence of Dissipative Small Scale Processes: Geotail Spacecraft Observation and Simulation of Electrostatic Solitary Waves

In recent spacecraft observations, coherent microscale structures such as electrostatic solitary waves are observed in various regions of the magnetosphere. The Geotail spacecraft observation has shown that these solitary waves are associated with high energy non-thermal electrons flowing along the magnetic field. The solitary structures are generated as a result of a long time evolution of coherent nonlinear trapping of electrons as found in bump-on-tail, bi-stream and Buneman instabilities. It is noted that these solitary waves can be generated at distant regions far away from the spacecraft locations, because these trapped electrons, or electron holes, are drifting much faster than the local thermal plasmas. Some of the solitary waves are accompanied by perpendicular electric fields indicating that two-or three-dimensional potential structures are passing by the spacecraft. Depending on the local plasma parameters, these multi-dimensional solitary structures couple with perpendicular modes such as electrostatic whistler modes and lower-hybrid modes. In a long time evolution, these perpendicular modes are dissipated via self-organization of small solitary potentials, leading to formation of one-dimensional potential troughs as observed in the deep magnetotail. The above dissipative small-scale processes are reproduced in particle simulations, and they can be used for diagnostics of electron dynamics from spacecraft observation of multi-dimensional solitary waves in various regions of the magnetosphere. This revised version was published online in July 2006 with corrections to the Cover Date.     

Spherical electron acoustic solitary waves in plasma with suprathermal electrons

A reductive perturbation technique is employed to solve the fluid-Poisson equations in spherical geometry describing a weakly nonlinear electron–acoustic (EA) waves in unmagnetized plasma consisting of stationary ions, cold electrons and kappa distributed hot electrons. It is shown that a variable coefficient Kadomtsev–Petviashvili (KP) equation governs the evolution of scalar potential describing propagation of EA waves. The influence of suprathermality and geometry effects on propagation of EA solitary waves is investigated. We found that when electrons evolve toward their thermodynamic equilibrium, EA solitons are generated with large amplitudes. Also it is shown that EA solitary structures can be significantly modified by transverse perturbations.     

太阳大气中动力学阿尔文波的产生与加热研究进展

动力学阿尔文波是垂直波长接近离子回旋半径或者电子惯性长度的色散阿尔文波.由于波的尺度接近粒子的动力学尺度,动力学阿尔文波在太阳和空间等离子体加热、加速等能化现象中起重要作用.因此,动力学阿尔文波通常被认为是日冕加热的候选者.本研究工作深入、系统地调研了太阳大气中动力学阿尔文波的激发和耗散机制.基于日冕等离子体环境,介绍了几种常见的动力学阿尔文波激发机制:温度各向异性不稳定性、场向电流不稳定性、电子束流不稳定性、密度非均匀不稳定性以及共振模式转换.还介绍了太阳大气中动力学阿尔文波的耗散机制,并讨论了这些耗散机制对黑子加热、冕环加热以及冕羽加热的影响.不仅为认识太阳大气中动力学阿尔文波的驱动机制、动力学演化特征以及波粒相互作用提供合理的理论依据,而且有助于揭示日冕等离子体中能量储存和释放、粒子加热等能化现象的微观物理机制.     

Ion-temperature-gradient driven modes in dust-contaminated plasma with nonthermal electron distribution and dust charge fluctuations

The effect of nonthermal electrons on ion-temperature-gradient (ITG) driven modes is investigated in the presence of variable dust charge and ion shear flow. The dust charge fluctuating expression is obtained in the presence of kappa distributed electrons. A dispersion relation is derived and analyzed numerically by choosing space plasma parameters of Jupiter/Saturn magnetospheres. It is found that the presence of nonthermal electrons population reduces the growth rate of ITG mode driven instability. The effects of ion temperature, electron density and magnetic field variation on the growth rate of ITG instability are presented numerically. It is also pointed out that the present results will be useful to understand the ITG driven modes with variable dust charge and kappa distributed electrons, present in most of the space plasma environments.     

The runaway instability of thick discs around black holes I. The constant angular momentum case

We present results from a numerical study of the runaway instability of thick discs around black holes. This instability is an important issue for most models of cosmic gamma-ray bursts, where the central engine responsible for the initial energy release is such a system consisting of a thick disc surrounding a black hole. We have carried out a comprehensive number of time-dependent simulations aimed at exploring the appearance of the instability. Our study has been performed using a fully relativistic hydrodynamics code. The general relativistic hydrodynamic equations are formulated as a hyperbolic flux-conservative system and solved using a suitable Godunov-type scheme. We build a series of constant angular momentum discs around a Schwarzschild black hole. Furthermore, the self-gravity of the disc is neglected and the evolution of the central black hole is assumed to be that of a sequence of exact Schwarzschild black holes of varying mass. The black hole mass increase is thus determined by the mass accretion rate across the event horizon. In agreement with previous studies based on stationary models, we find that by allowing the mass of the black hole to grow the disc becomes unstable. Our hydrodynamical simulations show that for all disc-to-hole mass ratios considered (between 1 and 0.05), the runaway instability appears very fast on a dynamical time-scale of a few orbital periods, typically a few 10 ms and never exceeding 1 s for our particular choice of the mass of the black hole (2.5 M) and a large range of mass fluxes  ( m 10^-3 M s^-1)  . The implications of our results in the context of gamma-ray bursts are briefly discussed.     

Pair emission from bare magnetized strange stars

The dominant emission from bare strange stars is thought to be electron–positron pairs, produced through spontaneous pair creation (SPC) in a surface layer of electrons tied to the star by a superstrong electric field. The positrons escape freely, but the electrons are directed towards the star and quickly fill all available states, such that their degeneracy suppresses further SPC. An electron must be reflected and gain energy in order to escape, along with the positron. Each escaping electron leaves a hole that is immediately filled by another electron through SPC. We discuss the collisional processes that produce escaping electrons. When the Landau quantization of the motion perpendicular to the magnetic field is taken into account, electron–electron collisions can lead to an escaping electron only through a multistage process involving higher Landau levels. Although the available estimates of the collision rate are deficient in several ways, it appears that the rate is too low for electron–electron collisions to be effective. A simple kinetic model for electron–quark collisions leads to an estimate of the rate of pair production that is analogous to thermionic emission, but the work function is poorly determined.     

Enhanced energetic electron intensities at 100 km altitude and a whistler propagating through the plasmasphere

During the flight of a Petrel rocket, instrumented by the SRC Radio and Space Research Station with Geiger counters and launched westwards from South Uist, Outer Hebrides, Scotland (L=3.38), a transient increase was observed in the intensity of energetic electrons having pitch angles between 60 and 120°. The increase, by a factor of 20 above the quasi-steady intensity observed throughout the remainder of the flight, occurred in 0.8 sec and was simultaneous for both >45 keV and >110 keV electrons. Recorded ～0.5 sec later, on the ground, was a two-hop whistler. During the enhanced electron intensity event, the entire duration of which was ～6 sec, the four-, six- and eight-hop whistlers were also received. From an analysis of the whistlers' spectrogram, it is concluded that the whistlers were ducted through the magnetosphere along the L=3.3 ±0.1 field line; the electron density in the equatorial plane is found to be 330 ±10 cm^?3, a value characteristic of conditions within the plasmapause. It is suggested that these temporally and/or spatially associated phenomena, rather than arising by a chance coincidence, were the result of a gyroresonant interaction between energetic electrons and whistler mode waves moving in opposite directions. For gyroresonance on this field line at the equator, the parallel component of energy of the electrons is 25 keV at 3 kHz in the whistler band, or 100 keV at 1 kHz below it. It is suggested that a magnetospheric event occurred, causing both sudden enhanced electron precipitation and favourable conditions for the propagation and/or amplification of whistlers. A possible explanation is that energetic electrons, having a sufficiently anisotropic distribution function and associated with those injected during an earlier auroral substorm, become unstable via the transverse resonance instability when they drift into the plasmasphere, a region of high density thermal plasma.     

Collective plasma effects and the electron number problem in solar hard X-ray bursts

Due to the relatively high stream densities involved, collective interactions with the ambient plasma are likely to be important for the electrons producing solar hard X-ray bursts. In thick- and thin-target bremsstrahlung models the most relevant process is limitation of the invoked electron beams by ion sound wave generation in the neutralizing reverse current established in the atmosphere. For the thick target model it is shown that typical electron fluxes are near the maximum permitted by stability of the reverse current so that ion-sound wave generation may be the process which limits the electron injection rate. On the other hand the chromospheric reverse current is sufficient to supply the large total number of electrons which have to be accelerated in the corona. For the thin target the low density of the corona severely limits the possible reverse current so that the maximum upward flux of fast electrons is probably much too small to explain X-ray bursts but compatible with observations of interplanetary electrons.A distinct class of model postulates a small number of electrons confined by resonant scattering in a dense coronal slab surrounding a current sheet with continuous stochastic acceleration offsetting collisional losses. The energetic aspects of such a situation described by Hoyng (1975) are developed here by addition of equations describing the slab geometry in terms of electron diffusion by whistler scattering and of the collisional damping of the accelerating Langmuir waves. Solution of these equations results in values for the fieldB(70–350 G), densityn ₀(2–5 × 10¹² cm ^–3), slab dimensions (10¹⁸ km² × 0.3–3 km) and relative Langmuir energy density (10^–3 – 10^–2) required to produce the observed range of bursts. It is pointed out, however, that there may be no real gain in electron number requirements since the fast electrons in the emitting slab would be constantly swept out along with the frozen-in plasma as dissipation proceeds so that a large total number of electrons is still required. It could in fact be that just such a coronal region is the injection mechanism for the thick-target model.On leave from Department of Astronomy, University of Glasgow, Scotland.     

Hydromagnetic wave propagation and the Kelvin-Helmholtz instability in an arbitrary-β collisional inhomogeneous plasma

The instability of an inhomogeneous arbitrary- plasma occurring due to the transverse velocity shear, has been studied to analyse the effects of collisional thermal transfer. The dissipation of hydromagnetic waves in such a plasma has also been discussed. It has been found that the thermal forces modify the instability criteria in several limiting cases. Numerical solutions have also been obtained to investigate the effects of various physical parameters for a non-isothermal plasma with different adiabticity of two species, viz., electrons and ions.     

Some ideas concerning the problem of anomalous polar E region heating due to Farley-Buneman turbulence

The energetics of the excitation of the Farley-Buneman instability is considered, which is recently observed in the auroral and equatorial E regions of the Earth's ionosphere at altitudes between 100 km and 120 km. In the magnetic field of the Earth the Farley-Buneman instability is excited under the condition of a strong enough external electric field in the case of ion-neutral collisions with frequencies much larger than the ion gyrofrequency and electron-neutral collisions with frequencies much below the electron gyrofrequency. It is shown that the linear increase of the wave amplitudes is caused by a small disbalance between the processes of nonlinear energy pumping into the wave from an external electric field and the energy loss because of the collisions of the electrons and ions with the neutral particles. During the nonlinear energy pumping energy of the external electric field is transferred into a nonlinear current of second order, which is connected with the oscillating motion of the electrons in the wave. The oscillating electron motion takes place perpendicular to wave propagation. From the estimations follows that the energy pumped into a Farley-Buneman wave during one period of pulsation is much larger than the wave energy itself. A new and simply to understand derivation of the anomalous diffusion coefficient is presented, related to the study of the behaviour of a test wave with frequency much above the frequencies of the Farley-Buneman turbulence in developed stage can cause an additional macroscopic nonlinear Pedersen current directed along the external electric field. It is found that the nonlinear Pedersen current can reach the order of the usual Pedersen current and should contribute to the effective heating of the ionospheric plasma.     

On the propagation of the electron streams generating type III bursts

The suggestion is explored that the two-stream instability has little effect on the propagation of the electron streams which generate type III bursts because the time required (ti) for development of the instability is comparable with or greater than the time available (Δt) for growth of the waves. Inferred parameters for streams in the corona and measured parameters for streams at the orbit of the Earth are compatible with this suggestion. Quasi-linear relaxation, which should occur as the stream forms, ensures that equality ti = Δt is set up initially, and restricts the number of escaping electrons to N _s ≈ 10³¹. The minimum density requirement on the stream for the two-stream instability to occur is found to be much less restrictive than the requirement that there should be many streaming electrons per Debye sphere.     

10–100 keV electron acceleration and emission from solar flares

We present an analysis of spacecraft observations of non-thermal X-rays and escaping electrons for 5 selected small solar flares in 1967. OSO-3 multi-channel energetic X-ray measurements during the non-thermal component of the solar flare X-ray bursts are used to derive the parent electron spectrum and emission measure. IMP-4 and Explorer-35 observations of > 22 keV and > 45 keV electrons in the interplanetary medium after the flares provide a measure of the total number and spectrum of the escaping particles. The ratio of electron energy loss due to collisions with the ambient solar flare gas to the energy loss due to bremsstrahlung is derived. The total energy loss due to collisions is then computed from the integrated bremsstrahlung energy loss during the non-thermal X-ray burst. For > 22 keV flare electrons the total energy loss due to collisions is found to be 10⁴ times greater than the bremsstrahlung energy loss and 10² times greater than the energy loss due to escaping electrons. Therefore the escape of electrons into the interplanetary medium is a negligible energetic electron loss mechanism and cannot be a substantial factor in the observed decay of the non-thermal X-ray burst for these solar flares.We present a picture of electron acceleration, energy loss and escape consistent with previous observations of an inverse relationship between rise and decay times of the non-thermal X-ray burst and X-ray energy. In this picture the acceleration of electrons occurs throughout the 10–100 sec duration of the non-thermal X-ray burst and determines the time profile of the burst. The average energy of the accelerated electrons first rises and then falls through the burst. Collisions with the ambient gas provide the dominant energetic electron loss mechanism with a loss time of 1 sec. This picture is consistent with the ratio of the total number of energetic electrons accelerated in the flare to the maximum instantaneous number of electrons in the flare region. Typical values for the parameters derived from the X-ray and electron observations are: total energy in > 22 keV electrons total energy lost by collisions = 10^28–29 erg, total number of electrons accelerated above 22 keV = 10³⁶, total energy lost by non-thermal bremsstrahlung = 10²⁴erg, total energy lost in escaping > 22 keV electrons = 10²⁶erg, total number of > 22 keV electrons escaping = 10^33–34.The total energy in electrons accelerated above 22 keV is comparable to the energy in the optical or quasi-thermal flare, implying a flare mechanism with particle acceleration as one of the dominant modes of energy dissipation.The overall efficiency for electron escape into the interplanetary medium is 0.1–1% for these flares, and the spectrum of escaping electrons is found to be substantially harder than the X-ray producing electrons.Currently at Tokyo Astronomical Observatory, Mitaka, Tokyo, Japan.     

Effects of non-thermal tails in Maxwellian electron distributions on synchrotron and Compton processes

We investigate how the presence of a non-thermal tail beyond a Maxwellian electron distribution affects the synchrotron process as well as Comptonization in plasmas with parameters typical for accretion flows on to black holes. We find that the presence of the tail can significantly increase the net (after accounting for self-absorption) cyclo-synchrotron emission of the plasma, which then provides seed photons for Compton upscattering. Thus, the luminosity in the thermally Comptonized spectrum is enhanced as well. The importance of these effects increases with both increasing Eddington ratio and black hole mass. The enhancement of the Comptonized synchrotron luminosity can be as large as ∼10³ and ∼10⁵ for stellar and supermassive black holes, respectively, when the energy content in the non-thermal tail is 1 per cent.
The presence of the tail only weakly hardens the thermal Comptonization spectrum but it leads to the formation of a high-energy tail beyond the thermal cut-off, which two effects are independent of the nature of the seed photons. Since observations of high-energy tails in Comptonization spectra can constrain the non-thermal tails in the electron distribution and thus the Comptonized synchrotron luminosity, they provide upper limits on the strength of magnetic fields in accretion flows. In particular, the measurement of an MeV tail in the hard state of Cyg X-1 by McConnell et al. implies the magnetic field strength in this source to be at most an order of magnitude below equipartition.     

The runaway instability of self-gravitating tori with non-constant specific angular momentum around black holes

We discuss the runaway instability of axisymmetric tori with non-constant specific angular momentum around black holes, taking into account self-gravity of the tori. The distribution of specific angular momentum of the tori is assumed to be a positive power law with respect to the distance from the rotational axis. By employing the pseudo-Newtonian potential for the gravity of the spherical black hole, we have found that self-gravity of the tori causes a runaway instability if the amount of the mass which is transferred from the torus to the black hole exceeds a critical value, i.e. 3 per cent of the mass of the torus. This has been shown by two different approaches: (1) by using equilibrium models and (2) by dynamical simulations. In particular, dynamical simulations using an SPH code have been carried out for both self-gravitating and non-self-gravitating tori. For non-self-gravitating models, all tori are runaway stable. Therefore we come to the conclusion that self-gravity of the tori has a stronger destabilizing effect than the stabilizing effect of the positive power-law distribution of the angular momentum.