X-ray bright knots of kpc-scale jets of several radio loud quasars have been an actively discussed issue. Among various models to explain observations, synchrotron radiation from the electron population different from radio to IR emitting electrons is promising. However, the origin of this electron population has been debated. Recently, we proposed that this electron population is produced by proton-photon collisions (mainly, Bethe-Heitler process), and we applied this model to PKS 0637-752. We found that this model works if the proton power is by an order of magnitude larger than the Eddington power. In this paper we apply this model to the X-ray emission in the knots of 3C 273 and PKS 1136-135. The target photons for electron-positron pair production are supplied by synchrotron radiation at radio-IR by primary electrons and by the active galactic nucleus (AGN) core as well as cosmic microwave background (CMB) radiation. The effects of the AGN photons are included for the first time in the hadronic model. Though the observed X-ray flux is obtained with the contribution of the AGN photons, the required proton power turns out to be highly super-Eddington. However, we find that our model works for a nearly Eddington proton power, if the photon density of the AGN is enhanced. This can occur if the AGN photons are more beamed toward the X-ray knots than toward our line of sight and the AGN photon frequency is shifted by the Doppler effect.

The Blandford-Znajek process, the steady electromagnetic energy extraction from a rotating black hole, is widely believed to work for driving relativistic jets, although it is still under debate where the electric current is driven. We address this issue analytically by investigating the time-dependent state in the Boyer-Lindquist and Kerr-Schild coordinate systems. This analysis suggests that a non-ideal magnetohydrodynamic region is required in the time-dependent state, while not in the steady state.

Strong X-ray emission from large scale jets of radio loud quasars still remains an open problem. Models based on inverse Compton scattering off CMB photons by relativistically beamed jets have recently been ruled out, since Fermi LAT observations for 3C 273 and PKS 0637-752 give the upper limit far below the model prediction. Synchrotron emission from a separate electron population with multi-hundred TeV energies remains a possibility although its origin is not well known. We examine a photo-hadronic origin of such high energy electrons/positrons, assuming that protons are accelerated up to $10^{19}$ eV and produce electrons/positrons through Bethe-Heitler process and photo-pion production. These secondary electrons/positrons are injected at sufficiently high energies and produce X-rays and $\gamma$-rays by synchrotron radiation without conflicting with the Fermi LAT upper limits. We find that the resultant spectrum well reproduces the X-ray observations from PKS 0637-752, if the proton power is at least $10^{49} \text{erg} \, \text{s}^{-1}$, which is highly super-Eddington. It is noted that the X-ray emission originates primarily from leptons through Bethe-Heitler process, while leptons from photo-pion origin lose energy directly through synchrotron emission of multi-TeV photons rather than cascading. To avoid the overproduction of the optical flux, optical emission is primarily due to synchrotron emission of secondary leptons rather than primary electrons, or a mild degree of beaming of the jet is needed, if it is owing to the primary electrons. Proton synchrotron luminosity is a few orders of magnitude smaller.

Blandford-Znajek process, the steady electromagnetic energy extraction from a rotating black hole (BH), is widely believed to work for driving relativistic jets in active galactic nuclei, gamma-ray bursts and Galactic microquasars, although it is still under debate how the Poynting flux is causally produced and how the rotational energy of the BH is reduced. We generically discuss the Kerr BH magnetosphere filled with a collisionless plasma screening the electric field along the magnetic field, extending the arguments of Komissarov and our previous paper, and propose a new picture for resolving the issues. For the magnetic field lines threading the equatorial plane in the ergosphere, we find that the inflow of particles with negative energy as measured in the coordinate basis is generated near that plane as a feedback from the Poynting flux production, which appears to be a similar process to the mechanical Penrose process. For the field lines threading the event horizon, we first show that the concept of the steady inflow of negative electromagnetic energy is not physically essential, partly because the sign of the electromagnetic energy density depends on the coordinates. Then we build an analytical toy model of a time-dependent process both in the Boyer-Lindquist and Kerr-Schild coordinate systems in which the force-free plasma injected continuously is filling a vacuum, and suggest that the structure of the steady outward Poynting flux is causally constructed by the displacement current and the cross-field current at the in-going boundary between the plasma and the vacuum. In the steady state, the Poynting flux is maintained without any electromagnetic source.

We examine the plasma composition of relativistic jets in four FRII radio galaxies by analyzing the total cocoon pressure in terms of partial pressures of thermal and non-thermal electrons/positrons and protons. The total cocoon pressure is determined by cocoon dynamics via comparison of theoretical model with the observed cocoon shape. By inserting the observed number density of non-thermal electrons/positrons and the upper limit of thermal electron/positron number density into the equation of state, the number density of protons is constrained. We apply this method to four FRII radio galaxies (Cygnus A, 3C 219, 3C 223 and 3C 284), for which the total cocoon pressures have been already evaluated. We find that the positron-free plasma comprising of protons and electrons is ruled out, when we consider plausible particle distribution functions. In other words, the mixture of positrons is required for all four FRII radio galaxies; the number density ratio of electrons/positrons to protons is larger than two. Thus, we find that the plasma composition is independent of the jet power and the size of cocoons. We also investigate the additional contribution of thermal electrons/positrons and protons on the cocoon dynamics. When thermal electrons/positrons are absent, the cocoon is supported by the electron/ proton plasma pressure, while both electron/positron pressure supported and electron/proton plasma pressure supported cocoons are allowed if the number density of thermal electrons/positrons is about 10 times larger than that of non-thermal ones.

We explore the degree of magnetization at the jet base of M87 by using the observational data of the event horizon telescope (EHT) at 230~GHz obtained by Doeleman et al. By utilizing the method in Kino et al., we derive the energy densities of magnetic fields ($U_{B}$) and electrons and positrons ($U_{\pm}$) in the compact region detected by EHT (EHT-region) with its full-width-half-maximum size $40~{\rm \mu as}$. First, we assume that an optically-thick region for synchrotron self absorption (SSA) exists in the EHT-region. Then, we find that the SSA-thick region should not be too large not to overproduce the Poynting power at the EHT-region. The allowed ranges of the angular size and the magnetic field strength of the SSA-thick region are $21~{\rm \mu as} \le \theta_{\rm thick}\le 26.3~{\rm \mu as}$ and $50~{\rm G}\le B_{\rm tot}\le 124~{\rm G}$, respectively. Correspondingly $U_{B}\gg U_{\pm} $ is realized in this case. We further examine the composition of plasma and energy density of protons by utilizing the Faraday rotation measurement ($RM$) at 230~GHz obtained by Kuo et al. Then, we find that $U_{B}\gg U_{\pm}+U_{p} $ still holds in the SSA-thick region. Second, we examine the case when EHT-region is fully SSA-thin. Then we find that $U_{B}\gg U_{\pm}$ still holds unless protons are relativistic. Thus, we conclude that magnetically driven jet scenario in M87 is viable in terms of energetics close to ISCO scale unless the EHT-region is fully SSA-thin and relativistic protons dominated.

We examine the applicability of the stochastic electron acceleration to two high synchrotron peaked blazars, Mrk 421 and Mrk 501, assuming synchrotron self-Compton emission of gamma-rays. Our model considers an emitting region moving at relativistic speed, where non-thermal electrons are accelerated and attain a steady-state energy spectrum together with the photons they emit. The kinetic equations of the electrons and photons are solved numerically, given a stationary wave number spectrum of the magnetohydrodynamic (MHD) disturbances, which are responsible for the electron acceleration and escape. Our simple formulation appears to reproduce the two well-sampled, long-term averaged photon spectra. In order to fit the model to the emission component from the radio to the X-ray bands, we need both a steeper wave spectral index than the Kolmogorov spectrum and efficient particle escape. Although the model provides a natural explanation for the high-energy cutoff of the electron energy distribution, the derived physical parameters raise a problem with an energy budget if the MHD waves with the Alfv{\'e}n velocity are assumed to be the acceleration agent.

We study effects of high-energy particles on the accretion flows onto a supermassive black hole and luminosities of escaping particles such as protons, neutrons, gamma-rays, and neutrinos. We formulate a one-dimensional model of the two-component accretion flow consisting of thermal particles and high-energy particles, supposing that some fraction of the released energy is converted to the acceleration of the high-energy particles. The thermal component is governed by fluid dynamics while the high-energy particles obey the moment equations of the diffusion-convection equation. By solving the time evolution of these equations, we obtain advection dominated flows as the steady state solutions. Effects of the high-energy particles on the flow structures turn out to be small even if the pressure of the high-energy particles dominates over the thermal pressure. For a model in which the escaping protons take away almost all the released energy, the high-energy particles have large influence enough to make the flow have the Keplerian angular velocity at the inner region. We calculate the luminosities of the escaping particles for these steady solutions. The escaping particles can extract the energy from about $10^{-4}\dot M c^2$ to $10^{-2}\dot M c^2$, where $\dot M$ is the mass accretion rates. The luminosities of the escaping particles depend on the parameters such as the injection Lorentz factors, the mass accretion rates, and the diffusion coefficients. We also discuss some implications on the relativistic jet production by the escaping particles.

One of the mechanisms widely considered for driving relativistic jets in active galactic nuclei, galactic microquasars, and gamma-ray bursts is the electromagnetic extraction of the rotational energy of a central black hole, i.e., the Blandford-Znajek process, although the origin of the electromotive force in this process is still under debate. We study this process as the steady unipolar induction in the Kerr black hole magnetosphere filled with a collisionless plasma screening the electric field (the D field) along the magnetic field (the B field), i.e., D dot B = 0. We extend the formulations and arguments made by Komissarov, and generally show that the origin of the electromotive force is ascribed to the ergosphere. It is explicitly shown that open magnetic field lines penetrating the ergosphere have a region where the D field is stronger than the B field in the ergosphere, and it keeps driving the poloidal currents and generating the electromotive force and the outward Poynting flux. The range of the possible value of the so-called angular velocity of the magnetic field line Omega_F is deduced for the field lines threading the equatorial plane in the ergosphere. We briefly discuss the relation between our conclusion and the ideal magnetohydrodynamic condition.

We explore energy densities of magnetic field and relativistic electrons in the M87 jet. Since the radio core at the jet base is identical to the optically thick surface against synchrotron self absorption (SSA), the observing frequency is identical to the SSA turnover frequency. As a first step, we assume the radio core as a simple uniform sphere geometry. Using the observed angular size of the radio core measured by the Very Long Baseline Array at 43 GHz, we estimate the energy densities of magnetic field ($U_{B}$) and relativistic electrons ($U_{e}$) based on the standard SSA formula. Imposing the condition that the Poynting power and relativistic electron one should be smaller than the total power of the jet, we find that (i) the allowed range of the magnetic field strength ($B_{tot}$) is from 1 G to 15 G, and that (ii) $1 times 10^{-5} < U_{e}/U_{B} < 6 times 10^{2}$ holds. The uncertainty of $U_{e}/U_{B}$ comes from the strong dependence on the angular size of the radio core and the minimum Lorentz factor of non-thermal electrons ($gamma_{e,min}$) in the core. It is still open that the resultant energetics is consistent with either the magnetohydrodynamic jet or with kinetic power dominated jet even on ~10 Schwarzschild radii scale.

We examine the radiation spectra from relativistic electrons moving in a Langmuir turbulence expected to exist in high energy astrophysical objects by using numerical method. The spectral shape is characterized by the spatial scale {\lambda}, field strength {\sigma}, and frequency of the Langmuir waves, and in term of frequency they are represented by {\omega}_0 = 2{\pi}c/{\lambda}, {\omega}_st = e{\sigma}/mc, and {\omega}_p, respectively. We normalize {\omega}_st and {\omega}_p by {\omega}_0 as \a \equiv {\omega}_st/{\omega}_0 and \b \equiv{\omega}_p/{\omega}_0, and examine the spectral shape in the a-b plane. An earlier study based on Diffusive Radiation in Langmuir turbulence (DRL) theory by Fleishman and Toptygin showed that the typical frequency is {\gamma}^2{\omega}_p and that the low frequency spectrum behaves as F_{\omega} pronto {\omega}^1 for b > 1 irrespective of a. Here, we adopt the first principle numerical approach to obtain the radiation spectra in more detail. We generate Langmuir turbulence by superposing Fourier modes,inject monoenergetic electrons, solve the equation of motion, and calculate the radiation spectra using Lienard-Wiechert potential. We find different features from the DRL theory for a>b>1. The peak frequency turns out to be {\gamma}^2{\omega}_st which is higher than {\gamma}^2{\omega}_p predicted in the DRL theory, and the spectral index of low frequency region is not 1 but 1/3. It is because the typical deflection angle of electrons is larger than the angle of the beaming cone sim 1/{\gamma}. We call the radiation for this case "Wiggler Radiation in Langmuir turbulence" (WRL).

The second-order Fermi acceleration (Fermi-II) driven by turbulence may be responsible for the electron acceleration in blazar jets. We test this model with time-dependent simulations. The hard electron spectrum predicted by the Fermi-II process agrees with the hard photon spectrum of 1ES 1101-232. For other blazars that show softer spectra, the Fermi-II model requires radial evolution of the electron injection rate and/or diffusion coefficient in the outflow. Such evolutions can yield a curved electron spectrum, which can reproduce the synchrotron spectrum of Mrk 421 from the radio to the X-ray regime.The photon spectrum in the GeV energy range of Mrk 421 is hard to fit with a synchrotron self-Compton model. However, if we introduce an external radio photon field with a luminosity of $4.9 \times 10^{38}~\mbox{erg}~\mbox{s}^{-1}$, GeV photons are successfully produced via inverse Compton scattering. The temporal variability of the diffusion coefficient or injection rate causes flare emission. The observed synchronicity of X-ray and TeV flares implies a decrease of the magnetic field in the flaring source region.

Pulsar winds have longstanding problems in energy conversion and pair cascade processes which determine the magnetization $\sigma$, the pair multiplicity $\kappa$ and the bulk Lorentz factor $\gamma$ of the wind. We study induced Compton scattering by a relativistically moving cold plasma to constrain wind properties by imposing that radio pulses from the pulsar itself are not scattered by the wind as was first studied by Wilson & Rees. We find that relativistic effects cause a significant increase or decrease of the scattering coefficient depending on scattering geometry. Applying to the Crab, we consider uncertainties of an inclination angle of the wind velocity with respect to the radio beam $\theta_{\rm pl}$ and the emission region size $r_{\rm e}$ which determines an opening angle of the radio beam. We obtain the lower limit $\gamma\gtrsim10^{1.7}r^{1/2}_{\rm e,3}\theta^{-1}_{\rm pl}(1+\sigma)^{-1/4}$ ($r_{\rm e}=10^3r_{\rm e,3}$ cm) at the light cylinder $r_{\rm LC}$ for an inclined wind $\theta_{\rm pl}>10^{-2.7}$. For an aligned wind $\theta_{\rm pl}<10^{-2.7}$, we require $\gamma>10^{2.7}$ at $r_{\rm LC}$ and an additional constraint $\gamma>10^{3.4}r^{1/5}_{\rm e,3}(1+\sigma)^{-1/10}$ at the characteristic scattering radius $r_{\rm c}=10^{9.6}r^{2/5}_{\rm e,3}$ cm within which the `lack of time' effect prevents scattering. Considering the lower limit $\kappa\gtrsim10^{6.6}$ suggested by recent studies of the Crab Nebula, for $r_{\rm e}=10^3$ cm, we obtain the most optimistic constraint $10^{1.7}\lesssim\gamma\lesssim10^{3.9}$ and $10^{6.6}\lesssim\kappa\lesssim10^{8.8}$ which are independent of $r$ when $\theta_{\rm pl}\sim1$ and $1+\sigma\sim1$ at $r_{\rm LC}$.

Relativistic jets in active galactic nuclei, galactic microquasars, and gamma-ray bursts are widely considered to be magnetohydrodynamically driven by black hole accretion systems, although conversion mechanism from Poynting into particle kinetic energy flux is still open. Recent detailed numerical and analytical studies of global structures of steady, axisymmetric magnetohydrodynamic (MHD) flows with specific boundary conditions have not reproduced as rapid an energy conversion as required by observations. In order to find more suitable boundary conditions, we focus on the flow along a poloidal magnetic field line just inside the external boundary, without treating transfield force balance in detail. We find some examples of the poloidal field structure and corresponding external pressure profile for an efficient and rapid energy conversion as required by observations, and that the rapid acceleration requires a rapid decrease of the external pressure above the accretion disk. We also clarify the differences between the fast magnetosonic point of the MHD flow and the sonic point of de Laval nozzle.

Among dozens young pulsar wind nebulae, some have been detected in TeV \gamma-rays (TeV PWNe), while others have not (non-TeV PWNe). The TeV emission detectability is not correlated either with the spin-down power or with the characteristic age of their central pulsars, and it is an open problem what determines the detectability. To study this problem, we investigate spectral evolution of five young non-TeV PWNe, 3C58, G310.6-1.6, G292.0+1.8, G11.2-0.3 and SNR B0540-69.3. We use a spectral evolution model which has been developed to be applied to young TeV PWNe in our previous works. TeV \gamma-ray flux upper limits of non-TeV PWNe give upper or lower limits on parameters, such as the age of the PWN and the fraction of the spin-down power going to the magnetic energy injection (the fraction parameter). Combined with other independent observational and theoretical studies, we can guess a plausible value of the parameters for each object. For 3C58, we prefer the parameters with an age of 2.5 kyr old and the fraction parameter of 3.0x10^{-3}, although the spectral modeling alone does not rule out a shorter age and a higher fraction parameter. The fraction parameter of 3.0x10^{-3} is also consistent for other non-TeV PWNe and then the value is regarded as common to young PWNe including TeV PWNe. Moreover, we find that the intrinsic properties of the central pulsars are similar, 10^{48-50}erg for the initial rotational energy and 10^{42-44}erg for the magnetic energy (2x10^{12} - 3x10^{13}G for the dipole magnetic field strength at their surfaces). The TeV detectability is correlated with the total injected energy and the energy density of the interstellar radiation field around PWNe. Except for G292.0+1.8, a broken power-law injection of the particles well reproduces the broadband emission from non-TeV PWNe.

The gamma ray flares of the Crab nebula detected by Fermi and AGILE satellites challenge our understanding of physics of pulsars and their nebulae. The central problem is that the peak energy of the flares exceeds the maximum energy E_{\mathrm{c}} determined by synchrotron radiation loss. However, when there exist turbulent magnetic fields with scales \lambda_{\mathrm{B}} smaller than 2\pi mc^2/eB, jitter radiation can emit photons with energy higher than E_{\mathrm{c}}. The scale required for the Crab flares is about two orders of magnitude less than the wavelength of the striped wind. We discuss the model in which the flares are triggered by plunging of the high density blobs into the termination shock. The observed hard spectral shape may be explained by jitter mechanism. We make three observational predictions: firstly the polarization degree will become lower in flares, secondly, no counterpart will be seen in TeV-PeV range, and thirdly the flare spectrum will not be harder than \nu F_\nu \propto \nu^1.

We propose a model for the gamma-ray binary LS 5039 in which the X-ray emission is due to the inverse Compton (IC) process instead of the synchrotron radiation. Although the synchrotron model has been discussed in previous studies, it requires a strong magnetic field which leads to a severe suppression of the TeV gamma-ray flux in conflict with H.E.S.S. observations. In this paper, we calculate the IC emission by low energy electrons (\gamma_e \lesssim 10^3) in the Thomson regime. We find that IC emission of the low energy electrons can explain the X-ray flux and spectrum observed with Suzaku if the minimum Lorentz factor of injected electrons \gamma_min is around 10^3. In addition, we show that the Suzaku light curve is well reproduced if \gamma_min varies in proportion to the Fermi flux when the distribution function of injected electrons at higher energies is fixed. We conclude that the emission from LS 5039 is well explained by the model with the IC emission from electrons whose injection properties are dependent on the orbital phase. Since the X-ray flux is primarily determined by the total number of cooling electrons, this conclusion is rather robust, although some mismatches between the model and observations at the GeV band remain in the present formulation.

Plasmas of geometrically thick, black hole (BH) accretion flows in active galactic nuclei (AGNs) are generally collisionless for protons, and involve magnetic field turbulence. Under such conditions a fraction of protons can be accelerated stochastically and create relativistic neutrons via nuclear collisions. These neutrons can freely escape from the accretion flow and decay into protons in dilute polar region above the rotating BH to form relativistic jets. We calculate geometric efficiencies of the neutron energy and mass injections into the polar region, and show that this process can deposit luminosity as high as L_j ~ 2e-3 dot{M} c^2 and mass loading dot{M}_j ~ 6e-4 dot{M} for the case of the BH mass M ~ 1e8 M_sun, where dot{M} is mass accretion rate. The terminal Lorentz factors of the jets are Gamma ~ 3, and they may explain the AGN jets having low luminosities. For higher luminosity jets, which can be produced by additional energy inputs such as Poynting flux, the neutron decay still can be a dominant mass loading process, leading to e.g., Gamma ~ 50 for L_{j,tot} ~ 3e-2 dot{M}c^2.

We examine plasma composition of jets in active galactic nuclei through the comparison of the total pressure ($P$) with partial pressures of electrons and protons in a cocoon. The total pressure is estimated from the analysis of an expanding cocoon dynamics. We determine the average kinetic energy per particle for several representative cases of particle energy distribution such as one- and two-temperature thermal plasmas and non-thermal electrons by evaluating the dissipation of total kinetic energy of the jet into the internal energy of cocoon plasma. The number density of the total electrons/positrons ($n_{\pm}$) in the cocoon is constrained by using the particle supply from hot spots and the absence of thermal bremsstrahlung emission from radio lobes. By inserting $P$, $n_{\pm}$ and the particle energy of each population into the equation of state, the number density ($n_{p}$) and pressure ($P_{p}$) of protons in the cocoon can be constrained. Applying this method to Cygnus A, we find that (i) electron/positron ($e^{\pm}$) pairs always dominate in terms of number density, but that (ii) either an "$e^{\pm}$-supported cocoon (i.e., $P_{\pm} >P_{p}$)" or "proton-supported one (i.e, $P_{\pm} <P_{p}$)" is possible.

Recent observations of Sgr A$^*$ by Fermi and HESS have detected steady gamma-ray emission in the GeV and TeV bands. We present a new model to explain the GeV gamma-ray emission by inverse Compton scattering by nonthermal electrons supplied by the NIR/X-ray flares of Sgr A$^*$. The escaping electrons from the flare regions accumulate in a region with a size of $\sim 10^{18}$ cm and magnetic fields of $\lesssim 10^{-4}$ G. Those electrons produce gamma-rays by inverse Compton scattering off soft photons emitted by stars and dust around the central black hole. By fitting the GeV spectrum, we find constraints on the magnetic field and the energy density of optical-UV radiation in the central 1 pc region around the supermassive black hole. While the GeV spectrum is well fitted by our model, the TeV $\gamma$-rays, whose spectral index is different from that of the GeV emission, may be from different sources such as pulsar wind nebulae.

We study four young Pulsar Wind Nebulae (PWNe) detected in TeV gamma-rays, G21.5-0.9, G54.1+0.3, Kes 75, and G0.9+0.1, using the spectral evolution model developed and applied to the Crab Nebula in our previous work. We model the evolution of magnetic field and particle distribution function inside a uniformly expanding PWN considering a time-dependent injection from the pulsar and radiative and adiabatic losses. Considering uncertainties in the interstellar radiation field (ISRF) and their distance, we study two cases for each PWN. Because TeV PWNe have a large TeV gamma-rays to X-rays flux ratio, the magnetic energy of the PWNe accounts for only a small fraction of the total energy injected (typically a few x 10^{-3}). The gamma-ray emission is dominated by inverse Compton scattering off the infrared photons of the ISRF. A broken power-law distribution function for the injected particles reproduces the observed spectrum well, except for G0.9+0.1. For G0.9+0.1, we do not need a low energy counterpart because adiabatic losses alone are enough to reproduce the radio observations. High energy power-law indices at injection are similar (2.5 -- 2.6), while low energy power-law indices range from 1.0 to 1.6. The lower limit of the particle injection rate indicates that the pair multiplicity is larger than 10^4. The corresponding upper limit of the bulk Lorentz factor of the pulsar winds is close to the break energy of the broken power-law injection, except for Kes 75.

We calculate numerically the radiation spectrum from relativistic electrons moving in small scale turbulent magnetic fields expected in high energy astrophysical sources. Such radiation spectrum is characterized by the strength parameter $a = \lambda_{B}} e|B|/mc^2$, where $\lambda_{B}}$ is the length scale of the turbulent field. When $a$ is much larger than the Lorentz factor of a radiating electron $\gamma$, synchrotron radiation is realized, while $a\ll 1$ corresponds to the so-called jitter radiation regime. Because for $1<a<\gamma$ we cannot use either approximations, we should have recourse to the Lienard-Wiechert potential to evaluate the radiation spectrum, which is performed in this paper. We generate random magnetic fields assuming Kolmogorov turbulence, inject monoenergetic electrons, solve the equation of motion, and calculate the radiation spectrum. We perform numerical calculations for several values of $a$ with $\gamma = 10$. We obtain various types of spectra ranging between jitter radiation and synchrotron radiation. For $a\sim 7$, the spectrum turns out to take a novel shape which has not been noticed up to now. It is like a synchrotron spectrum in the middle energy region, but in the low frequency region it is a broken power law and in the high frequency region an extra power law component appears beyond the synchrotron cutoff. We give a physical explanation of these features.

We study the escape of cosmic-ray (CR) protons accelerated at a supernova remnant (SNR) by numerically solving a diffusion-convection equation from the vicinity of the shock front to the region far away from the front. We consider the amplifications of Alfven waves generated by the escaping CR particles and their effects on CR escape into interstellar medium (ISM). We find that the amplification of the waves significantly delays the escape of the particles even far away from the shock front (on a scale of the SNR). This means that the energy spectrum of CR particles measured through gamma-ray observations at molecular clouds around SNRs is seriously affected by the particle scattering by the waves.

Sagittarius A$^*$ in the Galactic center harbors a supermassive black hole and exhibits various active phenomena. Besides quiescent emission in radio and submillimeter radiation, flares in the near infrared (NIR) and X-ray bands are observed to occur frequently. We study a time-dependent model of the flares, assuming that the emission is from a blob ejected from the central object. Electrons obeying a power law with the exponential cutoff are assumed to be injected in the blob for a limited time interval. The flare data of 2007 April 4 were used to determine the values of model parameters. The spectral energy distribution of flare emission is explained by nonthermal synchrotron radiation in the NIR and X-ray bands. The model light curves suggest that electron acceleration is still underway during the rising phase of the flares. GeV gamma-rays are also emitted by synchrotron self-Compton scattering, although its luminosity is not strictly constrained by the current model. If the GeV emission is faint, the plasma blob is dominated by the magnetic energy density over the electron kinetic energy density. Observations in the GeV band will clarify the origin of the blob.

H-alpha emission from supernova remnants (SNRs) implies the existence of neutral hydrogens in the ambient medium. In the precursor of an SNR shock modified by cosmic rays (CRs), upstream plasmas are pushed by the CR pressure, but neutral particles are not, so that the relative velocity appears and some neutral particles become pickup ions by the charge exchange process in the precursor. We investigate how the pickup ions generated in the precursor affect the shock structure and the particle acceleration. If the CR pressure is larger than 20% of the shock ram pressure, the compression of the subshock becomes smaller than that without pickup ions because of the pressure of the pickup ions. Moreover, even if the shock is modified by CRs, the total compression ratio can be smaller than 4. In addition, the pickup ions play an important role for the injection into the particle acceleration. If the shock is a quasi-perpendicular shock and if the the multiply reflected ion acceleration occurs, the CR spectrum can be harder than that of the test particle diffusive shock acceleration below GeV.

The emission mechanism of the gamma-ray binary LS5039 in energy bands of TeV, GeV, and X-ray is investigated. Observed light curves in LS5039 show that TeV and GeV fluxes anticorrelate and TeV and X-ray fluxes correlate. However, such correlated variations have not been explained yet reasonably at this stage. Assuming that relativistic electrons are injected constantly at the location of the compact object as a point source, and that they lose energy only by the inverse Compton (IC) process, we calculate gamma-ray spectra and light curves by the Monte Carlo method, including the full electromagnetic cascade process. Moreover, we calculated X-ray spectra and light curves by using the resultant electron distribution. As a result, we are able to reproduce qualitatively spectra and light curves observed by HESS, Fermi, and Suzaku for the inclination angle i = 30 dig and the index of injected electron distribution p = 2.5. We conclude that TeV-GeV anticorrelation is due to anisotropic IC scattering and anisotropic gamma-gamma absorption, and that TeV-X correlation is due to the dependence of IC cooling time on orbital phases. In addition, the constraint on the inclination angle implies that the compact object in LS5039 is a black hole.

We study the spectral evolution of PWNe taking into account the energy injected when they are young. We model the evolution of the magnetic field inside a uniformly expanding PWN. Considering time dependent injection from the pulsar and coolings by radiative and adiabatic losses, we solve the evolution of the particle distribution function. The model is calibrated by fitting the calculated spectrum to the observations of the Crab Nebula at an age of a thousand years. The spectral evolution of the Crab Nebula in our model shows that the flux ratio of TeV gamma-rays to X-rays increases with time, which implies that old PWNe are faint in X-rays, but not in TeV gamma-rays. The increase of this ratio is because the magnetic field decreases with time and is not because the X-ray emitting particles are cooled more rapidly than the TeV gamma-ray emitting particles. Our spectral evolution model matches the observed rate of the radio flux decrease of the Crab Nebula. This result implies that our magnetic field evolution model is close to the reality. Finally, from the viewpoint of the spectral evolution, only a small fraction of the injected energy from the Crab Pulsar needs to go to the magnetic field, which is consistent with previous studies.

We study the escape of cosmic-ray protons accelerated at a supernova remnant (SNR). We are interested in their propagation in interstellar medium (ISM) after they leave the shock neighborhood where they are accelerated, but when they are still near the SNR with their energy density higher than that in the average ISM. Using Monte-Carlo simulations, we found that the cosmic-rays with energies of <~TeV excite Alfven waves around the SNR on a scale of the SNR itself if the ISM is highly ionized. Thus, even if the cosmic-rays can leave the shock, scattering by the waves prevents them from moving further away from the SNR. The cosmic-rays form a slowly expanding cosmic-ray bubble, and they spend a long time around the SNR. This means that the cosmic-rays cannot actually escape from the SNR until a fairly late stage of the SNR evolution. This is consistent with some results of Fermi and H.E.S.S. observations.

Recent observations suggest that many old pulsar wind nebulae (PWNe) are bright TeV gamma-ray sources without a strong X-ray counterpart. In this paper, we study the spectral evolution of PWNe taking into account the energy which was injected when they were young for old PWNe. We model the evolution of the magnetic field and solve for the particle distribution inside a uniformly expanding PWN. The model is calibrated by fitting the calculated spectrum to the observations of the Crab Nebula at an age of a thousand years. We find that only a small fraction of the injected energy from the Crab Pulsar goes to the magnetic field, consistent with previous studies. The spectral evolution model of the Crab Nebula shows that the flux ratio of TeV gamma-rays to X-rays increases with time, which implies that old PWNe are faint at X-rays, but not at TeV gamma-rays. The increase of this ratio is primarily because the magnetic field decreases with time and is not because the X-ray emitting particles are cooled more rapidly than the TeV gamma-ray emitting particles. Our spectral evolution model matches the observed rate of the radio flux decrease of the Crab Nebula.

We study cosmic-ray acceleration in a supernova remnant (SNR) and the escape from it. We model nonthermal particle and photon spectra for the hidden SNR in the open cluster Westerlund 2, and the old-age mixed-morphology SNR W 28. We assume that the SNR shock propagates in a low-density cavity, which is created and heated through the activities of the progenitor stars and/or previous supernova explosions. We indicate that the diffusion coefficient for cosmic-rays around the SNRs is less than ~1% of that away from them. We compare our predictions with the gamma-ray spectra of molecular clouds illuminated by the cosmic-rays (Fermi and H.E.S.S.). We found that the spectral indices of the particles are ~2.3. This may be because the particles were accelerated at the end of the Sedov phase, and because energy dependent escape and propagation of particles did not much affect the spectrum.

H-alpha emission from supernova remnants (SNRs) implies the existence of neutral hydrogen in the circumstellar medium. Some of the neutral particles penetrating the shock are ionized by the charge exchange process and make a cold ion beam in the shock downstream region. We perform linear analyses of collisionless plasma instabilities between the cold beam and the hot downstream plasma. We find that, under typical SNR conditions, either the resonant instability or the Weibel instability is the most unstable. This mechanism may amplify the magnetic field to more than 100 micro G and changes the shock structure. As a result, the radio spectrum and the large magnetic field can be explained, apart from the widely discussed Bell's mechanism.

We present the analysis of Suzaku observations of the young open cluster Westerlund 2, which is filled with diffuse X-ray emission. We found that the emission consists of three thermal components or two thermal and one non-thermal components. The upper limit of the energy flux of the non-thermal component is smaller than that in the TeV band observed with H.E.S.S. This may indicate that active particle acceleration has stopped in this cluster, and that the accelerated electrons have already cooled. The gamma-ray emission observed with H.E.S.S. is likely to come from high-energy protons, which hardly cool in contrast with electrons. Metal abundances of the diffuse X-ray gas may indicate the explosion of a massive star in the past.

We investigate apparent internal structure of kiloparsec scale jets of AGNs arising from the energy dependent cooling of accelerated electrons and light travel time effect for relativistically moving sources. Using a simple cylindrical shell model, we find that the offsets between the peaks of X-ray, optical and radio brightness distributions observed for many cases are basically explained with this model. Assuming that electrons in the moving shell are instantaneously accelerated, X-rays are emitted for a shorter time scale and observed at the side nearer to the nucleus, while radio emission continues to the far side of the nucleus because of longer cooling time. The resultant offset turns out to be order of the jet diameter owing to the light travel time effect.

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 non-propagating 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 artefact of the periodic boundary conditions and has no counterpart in the supernova-shock scenario.

Through detailed numerical simulations, we demonstrate that relativistic outflows (Lorentz factor $\Gamma \sim 7$) of electron-positron pairs can be produced by radiative acceleration even when the flow starts from a nearly pair equilibrium state at subrelativistic temperatures. Contrary to the expectation that pairs annihilate during an expansion stage for such low temperatures, we find that most pairs can survive for the situations obtained in our previous work. This is because in the outflow-generating region the dynamical timescale is short enough even though the fireball is optically thick to scattering. Several problems that should be solved to apply to actual active galactic nucleus jets are discussed.

Electrostatic behavior of a collisionless plasma in the foot region of high Mach number perpendicular shocks is investigated through the two-dimensional linear analysis and electrostatic particle-in-cell (PIC) simulation. The simulations are double periodic and taken as a proxy for the situation in the foot. The linear analysis for relatively cold unmagnetized plasmas with a reflected proton beam shows that obliquely propagating Buneman instability is strongly excited. We also found that when the electron temperature is much higher than the proton temperature, the most unstable mode is the highly obliquely propagating ion two-stream instability excited through the resonance between ion plasma oscillations of the background protons and of the beam protons, rather than the ion acoustic instability that is dominant for parallel propagation. To investigate nonlinear behavior of the ion two-stream instability, we have made PIC simulations for the shock foot region in which the initial state satisfies the Buneman instability condition. In the first phase, electrostatic waves grow two-dimensionally by the Buneman instability to heat electrons. In the second phase, highly oblique ion two-stream instability grows to heat mainly ions. This result is in contrast to previous studies based on one-dimensional simulations, for which ion acoustic instability further heats electrons. The present result implies that overheating problem of electrons for shocks in supernova remnants is resolved by considering ion two-stream instability propagating highly obliquely to the shock normal and that multi-dimensional analysis is crucial to understand the particle heating and acceleration processes in shocks.

H.E.S.S. observed TeV blazar PKS 2155--304 in a strong flare state in 2006 July. The TeV flux varied on timescale as short as a few minutes, which sets strong constraints on the properties of the emission region. By use of the synchrotron self-Compton model, we found that models with the bulk Lorentz factor $\sim 100$, the size of the emission region $\sim 10^{15}$ cm, and magnetic field $\sim 0.1$ G explain the observed spectral energy distribution and the flare timescale $\sim$ a few minutes. This model with a large value of $\Gamma$ accounts for the emission spectrum not only in the TeV band but also in the X-ray band. The major cooling process of electrons/positrons in the jet is inverse Compton scattering off synchrotron photons. The energy content of the jet is highly dominated by particle kinetic energy over magnetic energy.

Invisible plasma content in blazar jets such as protons and/or thermal electron-positron ($e^{\pm}$) pairs is explored through combined arguments of dynamical and radiative processes. By comparing physical quantities required by the internal shock model with those obtained through the observed radio-to-gamma-ray spectra for Mrk 421, we find the existence of a copious amount of invisible plasma in the jet. We speculate that the blazar sequence could arise from variations of total amount and/or blending ratio of $e^{\pm}$ pair and electron-proton plasma.

Invisible plasma content in blazar jets such as protons and/or thermal electron-positron ($e^{\pm}$) pairs is explored through combined arguments of dynamical and radiative processes. By comparing physical quantities required by the internal shock model with those obtained through the observed broadband spectra for Mrk 421, we obtain that the ratio of the Lorentz factors of a pair of cold shells resides in about $2\sim 20$, which implies that the shocks are at most mildly relativistic. Using the obtained Lorentz factors, the total mass density $\rho$ in the shocked shells is investigated. The upper limit of $\rho$ is obtained from the condition that thermal bremsstrahlung emission should not exceed the observed $\gamma$-ray luminosity, whilst the lower limit is constrained from the condition that the energy density of non-thermal electrons is smaller than that of the total plasma. Then we find $\rho$ is $10^2$-$10^3$ times heavier than that of non-thermal electrons for pure $e^{\pm}$ pairs, while $10^2$-$10^6$ times heavier for pure electron-proton ($e/p$) content, implying the existence of a large amount of invisible plasma. The origin of the continuous blazar sequence is shortly discussed and we speculate that the total mass density and/or the blending ratio of $e^{\pm}$ pairs and $e/p$ plasma could be new key quantities for the origin of the sequence.

Electron acceleration in high Mach number perpendicular shocks is investigated through two-dimensional electrostatic particle-in-cell (PIC) simulation. We simulate the shock foot region by modeling particles that consist of three components such as incident protons and electrons and reflected protons in the initial state which satisfies the Buneman instability condition. In contrast to previous one-dimensional simulations in which strong surfing acceleration is realized, we find that surfing acceleration does not occur in two-dimensional simulation. This is because excited electrostatic potentials have a two-dimensional structure that makes electron trapping impossible. Thus, the surfing acceleration does not work either in itself or as an injection mechanism for the diffusive shock acceleration. We briefly discuss implications of the present results on the electron heating and acceleration by shocks in supernova remnants.

We numerically calculate the energy and momentum transfer rates due to Coulomb scattering between two fluids moving with a relative velocity. The results are fitted by simple functions. The fitting formulae are useful to simulate outflows from active galactic nuclei and compact high energy sources.

Motivated by relativistic jets observed in active galactic nuclei (AGN), we simulate outflows of electron-positron pairs strongly coupled with photons from normal electron-proton plasmas. Using multi-fluid approximation and a Monte Carlo method of radiative transfer, we obtain spherically symmetric, steady solutions of radiation and pair outflows for the luminosity $L \leq 10^{47}$ erg ${\rm s^{-1}}$. For microphysics, Coulomb scattering, Compton scattering, bremsstrahlung, electron-positron pair annihilation and creation are taken into account. Although a significant amount of pairs outflow by powerful radiative force with a mildly relativistic velocity, the temperature is not high enough to avoid pair annihilation before the fireball becomes optically thin to scattering. Several caveats in the simulations are also discussed.

We study the emission from an old supernova remnant (SNR) with an age of around 10^5 yrs and that from a giant molecular cloud (GMC) encountered by the SNR. When the SNR age is around 10^5 yrs, proton acceleration is efficient enough to emit TeV gamma-rays both at the shock of the SNR and that in the GMC. The maximum energy of primarily accelerated electrons is so small that TeV gamma-rays and X-rays are dominated by hadronic processes, pi^0-decay and synchrotron radiation from secondary electrons, respectively. However, if the SNR is older than several 10^5 yrs, there are few high-energy particles emitting TeV gamma-rays because of the energy loss effect and/or the wave damping effect occurring at low-velocity isothermal shocks. For old SNRs or SNR-GMC interacting systems capable of generating TeV gamma-ray emitting particles, we calculated the ratio of TeV gamma-ray (1-10 TeV) to X-ray (2-10 keV) energy flux and found that it can be more than ~10^2. Such a source showing large flux ratio may be a possible origin of recently discovered unidentified TeV sources.

The emission spectra of TeV blazars extend up to tens of TeV and the emission mechanism of the TeV $\gamma$-rays is explained by synchrotron self-Compton scattering in leptonic models. In these models the time variabilities of X-rays and TeV $\gamma$-rays are correlated. However, recent observations of 1ES 1959+650 and Mrk 421 have found the ``orphan'' TeV $\gamma$-ray flares, i.e., TeV $\gamma$-ray flares without simultaneous X-ray flares. In this paper we propose a model for the ``orphan'' TeV $\gamma$-ray flares, employing an inhomogeneous leptonic jet model. After a primary flare that accompanies flare-up both in X-rays and TeV $\gamma$-rays, radiation propagates in various directions in the comoving frame of the jet. When a dense region in the jet receives the radiation, X-rays are scattered by relativistic electrons/positrons to become TeV $\gamma$-rays. These $\gamma$-ray photons are observed as an ``orphan'' TeV $\gamma$-ray flare. The observed delay time between the primary and ``orphan'' flares is about two weeks and this is accounted for in our model for parameters such as $\Gamma = 20$, $d = 4 \times 10^{17}$cm, $\alpha = 3$, and $\eta = 1$, where $\Gamma$ is the bulk Lorentz factor of the jet, $d$ is the distance between the central black hole and the primary flare site, $\alpha/\Gamma$ is the angle between the jet axis and the direction of the motion of the dense region that scatters incoming X-rays produced by the primary flare, and $\eta/\Gamma$ is the angle between the jet axis and the line of sight.

We report detailed studies of the X-ray emission from the kpc scale jet in the nearest active galaxy, Cen A. 41 compact sources were found within the jet, 13 of which were newly identified. We construct the luminosity function for the detected jet-knots and argue that the remaining emission is most likely to be truly diffuse, rather than resulting from the pile-up of unresolved faint knots. The transverse jet profile reveals that the extended emission has the intensity peak at the jet boundaries. We note that limb-brightened jet morphologies have been observed previously at radio frequencies in some jet sources, but never so clearly at higher photon energies. Our result therefore supports a stratified jet model, consisting of a relativistic outflow including a boundary layer with a velocity shear. In addition, we found that the X-ray spectrum of the diffuse component is almost uniform across and along the jet. We discuss this spectral behavior within a framework of shock and stochastic particle acceleration processes. We note some evidence for a possible spectral hardening at the outer sheath of the jet. Due to the limited photon statistics of the present data, further deep observations of Cen A are required to determine the reality of this finding, however we note that the existence of the hard X-ray features at outer jet boundaries would provide an important challenge to theories for the evolution of ultra-relativistic particles within the jets.

H1426+428 is one of blazars that are observed by $\gamma$-rays in the TeV region. Because TeV $\gamma$-rays from distant sources are subject to attenuation by the extragalactic background light (EBL) via electron-positron pair production, the intrinsic spectrum of the TeV $\gamma$-rays should be inferred by using the models of radiation processes and EBL spectrum. We set constraints on the physical condition of H1426+428 with the synchrotron-self-Compton model applying several EBL models. We find that the emission region of H1426+428 is moving toward us with the bulk Lorentz factor of $\sim 20$ and that its magnetic field strength is $\sim 0.1$ G. These properties are similar to other TeV blazars such as Mrk 421 and Mrk 501. However, the ratio of the energy density of nonthermal electrons to that of the magnetic fields is about 190 and fairly larger than those of Mrk 421 and Mrk 501, which are about 5 -- 20. It is also found that the intensity of EBL in the middle and near infrared wavelengths should be low, i.e., the intensity at 10 $\mu$m is about 1 nW m$^{-2}$ sr$^{-1}$ to account for the observed TeV $\gamma$-ray flux. Because the spectral data of H1426+428 in X-rays and $\gamma$-rays used in our analysis were not obtained simultaneously, further observations of TeV blazars are necessary to make the constraints on EBL more stringent.

In the proton counterflow model of a pulsar magnetosphere that we have recently proposed, non-relativistic protons are supplied from the magnetosphere to flow toward the pulsar surface and screen an electric field above the polar cap region. In this Letter, we show that the proton counterflow is also suitable for the bunching of pair plasma. The two-stream instability is easily excited and can produce bunches of pairs with a relevant length scale to emit coherent curvature radiation.

Japanese magnetospheric explorer GEOTAIL recorded a detailed light curve during the initial 600 msec of a giant flare from SGR 1806-20 on December 27, 2004. We show that the observed light curve is well explained by an emission from relativistically expanding fireballs, like those of gamma-ray bursts (GRBs). Especially, the observed rapid fading after 500msec suggests that ejecta is collimated in a jet. We derive an upper limit on the jet opening half-angle of 0.2 radian that is as narrow as those of GRBs.

In blazars such as 3C 279, GeV gamma-rays are thought to be produced by inverse Compton scattering of soft photons injected from external sources into the jet. Because of the large bulk Lorentz factor of the jet, the energy of soft photons is Doppler shifted in the comoving frame of the jet, and the scattering is likely to occur in the Klein-Nishina regime. Although the Klein-Nishina effects are well known, the properties of the electron and emission spectra have not been studied in detail in the environment of blazars. We solve the kinetic equation of electrons with the spatial escape term of the electrons to obtain the electron energy spectrum in the jet and calculated the observed emission spectrum. In calculations of the Compton losses in the Klein-Nishina regime, we use the discrete loss formalism to take into account the significant energy loss in a single scattering. Although the scattering cross section decreases because of the Klein-Nishina effects, ample gamma rays are emitted by inverse Compton scattering. When the injection spectrum of electrons obeys a power law, the electron spectrum does not follow a broken power law, as a result of the Klein-Nishina effects, and a large number of high-energy electrons remain in the emitting region.

We propose a new mechanism to screen the electric field in the pulsar polar cap. Previous studies have shown that if an electron beam from the stellar surface is accelerated to energies high enough to create electron-positron pairs, the required electric field parallel to the magnetic field lines is too strong to be screened out by the produced pairs. We argue here that if non-relativistic protons are supplied from the magnetosphere to flow towards the stellar surface, they can provide an anode to screen out such a strong electric field. Injected electron-positron pairs yield an asymmetry of the electrostatic potential around the screening point. The required pair creation rate in this model is consistent with the conventional models.

Motivated by the particle acceleration problem in pulsars, we numerically investigate electrostatic instability of electron-positron pairs injected in an external electric field. The electric field is expected to be so strong that we cannot neglect effects of spatial variation in the 0-th order distribution functions on the scale of the plasma oscillation. We assume that pairs are injected mono-energetically with 4-velocity $u_0>0$ in a constant external electric field by which electrons (positrons) are accelerated (decelerated). By solving linear perturbations of the field and distribution functions of pairs, we find a new type of electrostatic instability. The properties of the instability are characterized by $u_0$ and the ratio $R$ of the braking time-scale (determined by the external electric field) to the time-scale of the plasma oscillation. The growth rate is as large as a few times the plasma frequency. We discuss the possibility that the excited waves prevent positrons from returning to the stellar surface.

We explore the energetics and plasma composition in FR II sources using a new simple method of combining shock dynamics and radiation spectrum. The hot spots are identified with the reverse shocked region of jets. With the one-dimensional shock jump conditions taking account of the finite pressure of hot ICM, we estimate the rest mass and energy densities of the sum of thermal and non-thermal particles in hot spots. Independently, based on the Synchrotron Self-Compton (SSC) model, we estimate the number and energy densities of {\it non-thermal} electrons using the multi-frequency radiation spectrum of hot spots. We impose the condition that the obtained rest mass, internal energy, and number densities of non-thermal electrons should be lower than those of the total particles determined by shock dynamics. We apply this method to Cygnus A. We examine three extreme cases of pure electron-positron pair plasma (Case I), pure electron-proton plasma with separate thermalization (Case II), and pure electron-proton plasma in thermal-equilibrium (Case III). By detailed SSC analysis for Cygnus A and 3C123, we find that the energy density of non-thermal electrons is about 10 times larger than that of magnetic field. We find that the Case III is not acceptable because predicted photon spectra do not give a good fit to the observed one. We find that Case II can also be ruled out since the number density of non-thermal electrons exceeds that of the total number density. Hence, we find that only pure $e^{\pm}$ plasma (Case I) is acceptable among the three cases. Total kinetic power of jet and electron acceleration efficiency are also constrained by internal energy densities of non-thermal and total particles.

We study steady and spherically symmetric outflows of pure electron-positron pair plasma as a possible acceleration mechanism of relativistic jets up to the bulk Lorentz factor of greater than 10. These outflows are initiated by the ``Wien fireball'', which is optically thick to Compton scattering but thin to absorption and in a Wien equilibrium state between pairs and photons at a relativistic temperature.

We report the Chandra detection of an X-ray jet in 3C 15. The peak of the X-ray emission in the jet is 4.1'' (a projected distance of 5.1 kpc) from the nucleus, and coincident with a component previously identified in the radio and optical jets. We examine four models for the X-ray jet emission: (I) weak synchrotron cooling in equip., (II) moderate synchrotron cooling in equip., (III) weak synchrotron plus SSC cooling, and (IV) moderate synchrotron plus SSC cooling. We argue that case (II) can most reasonably explain the overall emission from knot C. Case (III) is also possible, but requires a large departure from equipartition and for the jet power to be comparable to that of the brightest quasars. Diffuse X-ray emission has also been detected, distributed widely over the full extent (63kpc x 25kpc) of the radio lobes. We compare the total energy contained in the lobes with the jet power estimated from knot C, and discuss the energetic link between the jet and the lobes. We argue that the fueling time (t_fuel) and the source age (t_src) are comparable for case (II), whereas t_fuel << t_src is likely for case (III). The latter may imply that the jet has a very small filling factor, ~10^{-3}. We consider the pressure balance between the thermal galaxy halo and non-thermal relativistic electrons in the radio lobes. Finally, we show that the X-ray emission from the nucleus is not adequately fitted by a simple absorbed power-law model, but needs an additional power-law with heavy absorption intrinsic to the source. Such a high column density is consistent with the presence of a dense, dusty torus which obscures the quasar nucleus.

The emission mechanisms of the blazar 3C 279 are studied by solving the kinetic equations of electrons and photons in a relativistically moving blob. The gamma-ray spectral energy distribution (SED) is fitted by inverse Compton scattering of external photons. The bulk Lorentz factor of the emitting blob is found to be 25, and the magnetic field is found to be 0.3 G. GeV gamma-rays are well explained by inefficiently cooled electrons because of the Klein-Nishina effects. The electron spectrum is not a broken power law with a steeper spectrum above a break energy, which is often used to fit the observed SED. The kinetic energy density of the nonthermal electrons dominates the magnetic energy density; this result is qualitatively the same as that for TeV blazars such as Mrk 421 and Mrk 501. The gamma-ray luminosity of 3C 279 is often observed to increase rapidly. We show that one of the better sampled gamma-ray flares can be explain by the internal shock model.

Gamma-ray bursts are generally considered to be the result of internal shocks generated in an inhomogeneous relativistic outflow that arises from a fireball. In such shocks, the Fermi acceleration of protons is naturally expected to be at least as efficient as that of electrons. We investigate the consequences of proton acceleration in the residual thermal photon field of a fireball, especially those on the emission spectra of photons. In contrast to most other studies, we do not invoke a direct electron acceleration in shock waves. We show that the residual radiation field of the fireball can ignite the photopion production and subsequent cascade, and that the photons emitted in this process further enhance the photon-initiated cascade. We find that Fermi accelerated protons with $\gtrsim 10^{13}$ eV efficiently bring their energy into pion production and subsequent photon emission. The particles cascading from the pions emit photons over a wide energy range. The photons emitted from electron--positron pairs distribute continuously from the GeV range down to the X-ray range, while muons can emit gamma-rays by synchrotron radiation with a break at around 1--10 MeV. We also discuss several radiation processes which may possibly produce a break feature in the MeV range, as observed in addition to muon synchrotron radiation.

We give a new coherent description of the first-order Fermi acceleration of particles in shock waves from the point of view of stochastic process of the individual particles, under the test particle approximation. The time development of the particle distribution function can be dealt with by this description, especially for relativistic shocks. We formulate the acceleration process of a particle as a two-dimensional Markov process in a logarithmic momentum-time space, and relate the solution of the Markov process with the particle distribution function at the shock front, for both steady and time-dependent case. For the case where the probability density function of the energy gain and cycle-time at each shock crossing of the particles obeys a scaling law in momentum, which is usually assumed in the literature, it is confirmed in more general form that the energy distribution of particles has the power-law feature in steady state. The equation to determine the exact power-law index which is applicable for any shock speed is derived and it is shown that the power-law index, in general, depends on the shape of the probability density function of the energy gain at each shock crossing; in particular for relativistic shocks, the dispersion of the energy gain can influence the power-law index. It is also shown that the time-dependent solution has a self-similarity for the same case.

We report the detection at X-rays of the radio/optical hotspot and knots of 3C 303 from a short (15 ksec) Chandra exposure in 2001 March. The X-ray morphology is similar to that of the radio/optical emission with peaks in the X-ray emission found at 5.5''(knot B), 9'' (knot C) and 17'' (hotspot) from the core of 3C 303. Despite the limited signal-to-noise ratio of the short Chandra exposure, the X-ray photon spectrum was measured for the hotspot. We construct the spectral energy distribution (SED) and find that the X-ray flux is well below the extrapolation of the radio-to-optical continuum, which we interpret as resulting from the production of X-rays via inverse Compton scattering of both synchrotron photons (SSC) and cosmic microwave background photons (EC/CMB). The magnetic field strength, region size, and the maximum energy of electrons are self-consistently determined for the hotspot to be B ~ 4.3 \mu G, R ~ 6.5x10^{21} cm, and G_max ~ 1.4x10^7. This implies a magnetic field strength a factor of ~30 below the equipartition value; B_eq ~ 150 \mu G. The origin of this large departure from equipartition is still uncertain, but the discrepancy is reduced if the plasma in the hotspot is moving with mildly relativistic speeds. Our observation of 3C 303, as well as recent Chandra detections of large scale jets and hotspots in a number of radio galaxies, confirm that particles are accelerated very efficiently in radio galaxies.

We investigate the relation between the properties of Brightest Cluster Galaxies (BCGs) and those of their host clusters. To quantify the properties of cluster hot gas, we employ the parameter $Z$ of the fundamental plane of X-ray clusters. It is found that the offset of the BCG from the peak of cluster X-ray emission is larger for smaller $Z$ clusters. The parameter $Z$ (not the redshift {\it z}), which mainly depends on virial density $\rho_{\rm {vir}}$, is considered to represent the formation epoch of a cluster. We thus consider that the offset of the BCG is correlated with the dynamical equilibrium state of its host cluster. On the contrary, no significant correlation is found between the absolute optical magnitude of the BCG and the parameter $Z$. If the extreme brightness of the BCG is mainly acquired in the course of cluster evolution by environmental effect, BCGs are expected to be brighter in large $Z$ clusters. Our result is not consistent with this simplified view. On the contrary, it is possible that the extreme brightness of the BCG is likely to be determined in the early history of cluster collapse.

We study self-consistency for the pulsar polar cap model. In general, there will remain an unscreened electric field at the pair production front (PPF). We have derived the condition of electric field screening by pairs in the presence of returning particles. A previous belief that pairs with a density higher than the Goldreich-Julian density immediately screen out the electric field beyond PPF is unjustified. Pairs have little contribution on screening the electric field, and the geometrical screening is the only known way of getting a finite potential drop.

We study on the self-consistency of the pulsar polar cap model, i.e., the problem of whether the field-aligned electric field is screened by electron-positron pairs that are injected beyond the pair production front. We solve the one-dimensional Poisson equation along a magnetic field line, both analytically and numerically, for a given current density incorporating effects of returning positrons, and we obtain the conditions for the electric-field screening. The formula which we obtained gives the screening distance and the return flux for given primary current density, field geometry and pair creation rate at the pair production front. If the geometrical screening is not possible, for instance, on field lines with a super-Goldreich-Julian current, then the electric field at the pair production front is constrained to be fairly small in comparison with values expected typically by the conventional polar cap models. This is because (1) positive space charge by pair polarization is limited to a small value, and (2) returning of positrons leave pair electrons behind. A previous belief that pair creation with a pair density higher than the Goldreich-Julian density immediately screens out the electric field is unjustified at least for for a super Goldreich-Julian current density. We suggest some possibilities to resolve this difficulty.

We present the results of a long-look monitoring of 3C273 with RXTE between 1996 and 2000. A total of 230 observations amounts to a net exposure of 845ksec, with this spectral and variability analysis of 3C273 covering the longest observation period available at hard X-ray energies. Flux variations by a factor of 4 have been detected over 4years, whereas less than 30% flux variations have been observed for individual flares on time-scales of 3 days. Two temporal methods, the power spectrum density (PSD) and the structure function (SF), have been used to study the variability characteristics of 3C273. The hard X-ray photon spectra generally show a power-law shape with a differential photon index of G = 1.6+-0.1. In 10 of 261 data segments, exceptions to power-law behaviour have been found: (i) an additional soft excess below 4keV, and (ii) a broad Fe fluorescent line feature with EW = 100-200 eV. Our new observations of these previously reported X-ray features may imply that 3C273 is a unique object whose hard X-ray emission occasionally contains a component which is not related to a beamed emission (Seyfert like), but most hard X-rays are likely to originate in inverse Compton radiation from the relativistic jet (blazar like). Multi-frequency spectra from radio to gamma-ray are presented in addition to our RXTE results. The X-ray time variability and spectral evolution are discussed in the framework of beamed, synchrotron self-Compton picture. We consider the ``power balance'' (both radiative and kinetic) between the accretion disk, sub-pc-scale jet, and the 10kpc-scale jet.

We investigate a novel mechanism of bulk acceleration for relativistic jets of pure electron-positron pair plasma.We assumed that relativistically hot plasma is optically thick for Compton scattering and then pair plasma coupled dynamically to photons may establish quasi-Wien equilibria. We name such an electron-positron pair plasma the "Wien-fireball". Assuming there is a Wien-fireball, it is easy to form highly relativistic jets not only because coupled photons do not cause radiation drag effect, but also because coupled photons prevent catastrophic annihilation of pair plasma. It is found that the Wien-fireball provides a very effective mechanism to form high power jets. And the radiation of the Wien-fireball in blazars should be observed as MeV peaked emission at infinity with an order of kinetic power of jets.

Using multi-frequency spectra from TeV blazars in quiescent states, we obtain the physical parameters of the emission region of blazars within the framework of the one-zone synchrotron self-Compton (SSC) model. We numerically calculate the steady-state energy spectra of electrons by self-consistently taking into account the effects of radiative cooling with a proper account of the Klein-Nishina effects. Here electrons are assumed to be injected with a power-law spectrum and to escape on a finite time scale, which naturally leads to the existence of a break energy scale. Although we do not use time variabilities but utilize a model of electron escape to constrain the size of the emission region, the resultant size turns out to be similar to that obtained based on time variabilities. Through detailed comparison of the predicted emission spectra with observations, we find that for Mrk 421, Mrk 501, and PKS 2155--304, the energy density of relativistic electrons is about an order of magnitude larger than that of magnetic fields with an uncertainty within a factor of a few.

We have studied the rapid X-ray variability of three extragalactic TeV gamma-ray sources; Mrk 421, Mrk 501 and PKS 2155-304. Analyzing the X-ray light curves obtained from ASCA and/or RXTE observations between 1993 and 1998, we have investigated the variability in the time domain from 10^3 to 10^8 sec. For all three sources, both the power spectrum density (PSD) and the structure function (SF) show a roll-over with a time-scale of the order of 1 day or longer, which may be interpreted as the typical time-scale of successive flare events. Although the exact shape of turn-over is not well constrained and the low-frequency (long time-scale) behavior is still unclear, the high-frequency (short time-scale) behavior is clearly resolved. We found that, on time-scales shorter than 1 day, there is only small power in the variability, as indicated by a steep power spectrum density of f^{-2~-3}. This is very different from other types of mass-accreting black-hole systems for which the short time-scale variability is well characterized by a fractal, flickering-noise PSD (f^{-1~-2}). The steep PSD index and the characteristic time-scale of flares imply that the X-ray emitting site in the jet is of limited spatial extent; D > 10^{17} cm distant from the base of the jet, which corresponds to > 100 Schwarzschild radii for 10^{7-10} M_solar black-hole systems.

We study the temporal/spectral variability of two extragalactic TeV sources, Mrk 421 and Mrk 501, based on 5-year observations with the ASCA and RXTE satellites. We found that the peak of the synchrotron emission exists just in the X-ray band and its position shifted from lower to higher energy when the source became brighter. Relationship between the peak energy and peak luminosity showed quite different behavior in the two sources; Mrk 421 showed very little change in the peak position (0.5-2 keV), while Mrk 501 revealed the largest shift ever observed in blazars (1-100 keV). We analyze these X-ray data with the flux changes in TeV band, which are obtained from 35 truly simultateneous observations. Very different spectral evolution of both objects indicates some differences in the electron acceleration mechanism at work during the flares. We argue that the flux variability of Mrk 421 is associated with an increase in the number of electrons, while the flare of Mrk 501 is mostly due to the large changes in maximum energy of electrons. We also discuss the characteristic temporal variability of TeV sources, and implications for the X-ray emitting site in the relativistic jet.

We investigate statistical properties of galaxy clusters in the context of hierarchical clustering scenario, taking account of their formation epoch distribution, motivated by the recent finding by Fujita and Takahara that X-ray clusters form a fundamental plane, where the mass and formation epoch are regarded as two independent parameters. Using the formalism which discriminates between major merger and accretion, the epoch of a cluster formation is identified with that of the last major merger. Since tiny mass accretion after the formation does not much affect the core structure of clusters, the properties of X-ray emission from clusters are determined by the total mass and density at their formation time. Under these assumptions, we calculate X-ray luminosity and temperature functions of galaxy clusters. We find that the behavior of luminosity function is different from the model which does not take account of formation epoch distribution, while the behavior of temperature function is not much changed. In our model, luminosity function is shifted to a higher luminosity and shows no significant evolution up to $z \sim 1$, independent of cosmological models. The clusters are populated on the temperature-luminosity plane with a finite dispersion. Since the simple scaling model in which the gas temperature is equal to the virial temperature fails to reproduce the observed luminosity-temperature relation, we also consider a model which takes the effects of preheating into account. The preheating model reproduces the observations much better.

We formulate the first order Fermi acceleration in parallel shock waves in terms of the random walk theory. The formulation is applicable to any value of the shock speed and the particle speed, in particular to the acceleration in relativistic shocks, as long as large angle scattering is suitable for the scattering process of particles. We first show that the trajectory of a particle suffering from large angle scattering can be treated as a random walk in a moving medium with an absorbing boundary (e.g., the shock front). We derive an integral equation to determine the density of scattering points of the random walk, and by solving it approximately we obtain approximate solutions of the probability density of pitch angle at and the return probability after the shock crossing in analytical form. These approximate solutions include corrections of several non-diffusive effects to the conventional diffusion approximation and we show that they agree well with the Monte Carlo results for isotropic scattering model for any shock speed and particle speed. Finally, we give an analytical expression of the spectral index of accelerated particles in parallel shocks valid for arbitrary shock speed.

We conducted a multi-frequency campaign for the TeV blazar Mrk~421 in 1998 April. The campaign started from a pronounced high amplitude flare recorded by SAX and Whipple; ASCA observation started three days later. In the X-ray data, we detected multiple flares, occuring on time scales of about one day. ASCA data clearly reveal spectral variability. The comparison of the data from ASCA, EUVE and RXTE indicates that the variability amplitudes in the low energy synchrotron component are larger at higher photon energies. In TeV Gamma-rays, large intra-day variations -- which were correlated with the X-ray flux -- were observed when results from three Cherenkov telescopes are combined. The RMS variability of TeV Gamma--rays was similar to that observed in hard X-rays, above 10 keV. The X-ray light curve reveals flares which are almost symmetric for most of cases, implying the dominant time scale is the light crossing time through the emitting region. The structure function analysis based on the continuous X-ray light curve of seven days indicates that the characteristic time scale is ~0.5 day. The analysis of ASCA light curves in various energy bands appears to show both soft (positive) and hard (negative) lags. These may not be real, as systematic effects could also produce these lags, which are all much smaller than an orbit. If the lags of both signs are real, these imply that the particle acceleration and X-ray cooling time scales are similar.

Blazars are known to emit a broad band emission from radio to gamma-rays with rapid time variations, particularly, in X- and gamma-rays. Synchrotron radiation and inverse Compton scattering are thought to play an important role in emission and the time variations are likely related to the acceleration of nonthermal electrons. As simultaneous multiwavelength observations with continuous time spans are recently available, some characteristics of electron acceleration are possibly inferred from the spectral changes of high energy emission. In order to make such inferences, we solve the time-dependent kinetic equations of electrons and photons simultaneously using a simple model for electron acceleration. We then show how the time variations of emission are dependent on electron acceleration. We also present a simple model for a flare in X-rays and TeV gamma-rays by temporarily changing the acceleration timescale. Our model will be used, in future, to analyze observed data in detail to obtain information on electron acceleration in blazars.

We investigate the origin of the variation of the gas mass fraction in the core of galaxy clusters, which was indicated by our work on the X-ray fundamental plane. The adopted model supposes that the gas distribution characterized by the slope parameter is related to the preheated temperature. Comparison with observations of relatively hot (~> 3 keV) and low redshift clusters suggests that the preheated temperature is about 0.5-2 keV, which is higher than expected from the conventional galactic wind model and possibly suggests the need for additional heating such as quasars or gravitational heating on the largest scales at high redshift. The dispersion of the preheated temperature may be attributed to the gravitational heating in subclusters. We calculate the central gas fraction of a cluster from the gas distribution, assuming that the global gas mass fraction is constant within a virial radius at the time of the cluster collapse. We find that the central gas density thus calculated is in good agreement with the observed one, which suggests that the variation of gas mass fraction in cluster cores appears to be explained by breaking the self-similarity in clusters due to preheated gas. We also find that this model does not change major conclusions on the fundamental plane and its cosmological implications obtained in previous papers, which strongly suggests that not only for the dark halo but also for the intracluster gas the core structure preserves information about the cluster formation.

Using X-ray data, we investigate interrelations between gas density rho_0, virial density rho_vir, core radius R, temperature T, entropy S_gas, and metal abundance Z in the core region of clusters of galaxies. First, we confirm that fundamental relations among rho_0, R, and T found by Fujita and Takahara are reproduced by another data catalogue. Second, we find that, when clusters have two components in their surface brightness distribution, the inner components also satisfy the same fundamental relations on the assumption that the average temperature of the inner component is the same as that of the outer component. These results strengthen our interpretation that clusters form a two parameter family in terms of mass and rho_vir; larger rho_vir corresponds to earlier formation epoch. We argue that the inner components represent distinct dark matter components which collapsed ahead of the outer components. Third, we also find a tight relation between S_gas and rho_vir both for the outer and inner components; S_gas is smaller for larger rho_vir but is larger than that produced through gravitational collapse alone for larger rho_vir. Although radiative cooling affects the thermal evolution, the tight relationship discovered suggests the existence of stable heating sources or stable energy transmission mechanisms. Finally, we find that the iron abundance at the centers of clusters is correlated with rho_vir for the inner components. This implies that iron produced by Type Ia supernovae has more accumulated for clusters formed earlier. We briefly discuss the implications of these findings.

Based on the two-parameter family nature of X-ray clusters of galaxies obtained in a separate paper, we discuss the formation history of clusters and cosmological parameters of the universe. Utilizing the spherical collapse model of cluster formation, and assuming that the cluster X-ray core radius is proportional to the virial radius at the time of the cluster collapse, the observed relations among the density, radius, and temperature of clusters imply that cluster formation occurs in a wide range of redshift. The observed relations favor the low-density universe. Moreover, we find that the model of $n\sim -1$ is preferable.

We analyze the relations among central gas density, core radius, and temperature of X-ray clusters by plotting the observational data in the three-dimensional ($\log \rho_0$, $\log R$, and $\log T$) space and find that the data lie on a 'fundamental plane'. Its existence implies that the clusters form a two-parameter family. The data on the plane still has a correlation and form a band on the plane. The observed relation $L_{\rm X} \propto T^3$ turns out to be the cross section of the band perpendicular to the major axis, while the major axis is found to describe the virial density. We discuss implications of this two-parameter family nature of X-ray clusters.

We present a new method for recovering the kinetic temperature of the intervening diffuse gas to an accuracy of 10%. The method is based on the comparison of unsaturated absorption-line profiles of two species with different atomic weights. The species are assumed to have the same temperature and bulk motion within the absorbing region. The computational technique involves the Fourier transform of the absorption profiles and the consequent Entropy-Regularized chi^2-Minimization [ERM] to estimate the model parameters. The procedure is tested using synthetic spectra of CII, SiII and FeII ions. The comparison with the standard Voigt fitting analysis is performed and it is shown that the Voigt deconvolution of the complex absorption-line profiles may result in estimated temperatures which are not physical. We also successfully analyze Keck telescope spectra of CII1334 and SiII1260 lines observed at the redshift z = 3.572 toward the quasar Q1937--1009 by Tytler {\it et al.}.

We present data for 18 blazars observed with ASCA, half of which were also observed contemporaneously with EGRET as parts of multi-wavelength campaigns. The observations show a clear difference in the spectra between three subclasses of blazars, namely the High-energy peaked BL Lac objects (HBLs), Low-energy peaked BL Lac objects (LBLs), and quasar-hosted blazars (QHBs). We find that the radiation process responsible for the HE peak for HBLs can be explained solely by Synchrotron-Self-Compton (SSC) emission, with the Doppler factor consistent with the VLBI and/or gamma-ray variability data. For many QHBs, on the other hand, the gamma-rays cannot be solely due to the SSC mechanism unless the Doppler factor is significantly in excess of that inferred from VLBI data. We consider an alternative scenario consistent with the measured values of the Doppler factor, where the SSC component is still present in QHBs and it dominates in the X-ray band, but it is below the observed gamma-ray spectrum. With an assumption that the peak of the SSC emission is on the extrapolation of the X-ray spectrum, and adopting the Doppler factor of 10, we infer the magnetic field to be 0.1 - 1 Gauss, and Lorentz factors of electrons radiating at the peak of the nu F(nu) spectrum of 10^3 for QHBs; this is much lower than 10^5 for HBLs. This difference is most likely due to the large photon density expected in QHBs.

Deuterium abundance re-measurements by Burles and Tytler (1998; hereafter BT) yielded D/H = (3.3 +/- 0.3) 10^{-5} and the robust upper limit D/H < 3.9 10^{-5} from the z_a = 3.572 system toward Q1937-1009. In this new analysis BT adopted multicomponent microturbulent models together with the possibility to vary freely the local continuum level around each HI line to improve the fit. The procedure failed, however, to fit adequately D Ly-beta without recourse to an additional H Ly-alpha contamination at the position of D Ly-beta. We show that this obstacle may be successfully overcome within the framework of the mesoturbulent model accounting (in contrast to the microturbulent approximation) for a correlated structure of the large scale velocity field. Using the same observational data and the original continuum as determined by Tytler et al. (1996), we obtained good fits. The one-component mesoturbulent models provide D/H in the range (3.2 - 4.8) 10^{-5} and the total hydrogen column density N(HI) = (5.6 - 7.0) 10^{17} cm^{-2}. This result is consistent with that found by us from the z_a = 2.504 and z_a = 0.701 systems toward Q1009+2956 and Q1718+4807, respectively. The range for D/H common to all three analyses is D/H = (4.1 - 4.6) 10^{-5}. This value is consistent with standard big bang nucleosynthesis [SBBN] if the baryon-to-photon ratio, \eta, is in the range 4.2 10^{-10} <= \eta <= 4.6 10^{-10}, implying 0.0155 <= \Omega_b h^2_{100} <= 0.0167.

The apparent discrepancy between low and high D abundances derived from QSO spectra may be caused by spatial correlations in the stochastic velocity field. If one accounts for such correlations, one finds good agreement between different observations and the theoretical predictions for standard big bang nucleosynthesis (SBBN). In particular, we show that the H+D Ly-alpha profile observed at z = 0.7 toward Q1718+4807 is compatible with 4.1*10^{-5} <= D/H <= 4.7*10^{-5}. This result is consistent with our previous D/H determination for the z = 2.504 system toward Q1009+2956 and, thus, supports SBBN.

Deuterium abundances measured recently from QSO absorption-line systems lie in the range from 3 10^{-5} to 3 10^{-4}, which shed some questions on standard big bang theory. We show that this discordance may simply be an artifact caused by inadequate analysis ignoring spatial correlations in the velocity field in turbulent media. The generalized procedure (accounting for such correlations) is suggested to reconcile the D/H measurements. An example is presented based on two high-resolution observations of Q1009+2956 (low D/H) [1,2] and Q1718+4807 (high D/H) [8,9]. We show that both observations are compatible with D/H = 4.1 - 4.6 10^{-5}, and thus support SBBN. The estimated mean value <D/H> = 4.4 10^{-5} corresponds to the baryon-to-photon ratio during SBBN eta = 4.4 10^{-10} which yields the present-day baryon density Omega_b h^2 = 0.015.

The current status of extragalactic deuterium abundance is discussed using two examples of `low' and `high' D/H measurements. We show that the discordance of these two types of D abundances may be a consequence of the spatial correlations in the stochastic velocity field. Within the framework of the generalized procedure (accounting for such effects) one finds good agreement between different observations and the theoretical predictions for standard big bang nucleosynthesis (SBBN). In particular, we show that the deuterium absorption seen at z = 2.504 toward Q1009+2956 and the H+D Ly-alpha profile observed at z = 0.701 toward Q1718+4807 are compatible with D/H $\sim 4.1 - 4.6\times10^{-5}$. This result supports SBBN and, thus, no inhomogeneity is needed. The problem of precise D/H measurements is discussed.

A new method based on a Reverse Monte Carlo [RMC] technique and aimed at the inverse problem in the analysis of interstellar (intergalactic) absorption lines is presented. The line formation process in chaotic media with a finite correlation length $(l > 0)$ of the stochastic velocity field (mesoturbulence) is considered. This generalizes the standard assumption of completely uncorrelated bulk motions $(l \equiv 0)$ in the microturbulent approximation which is used for the data analysis up-to-now. It is shown that the RMC method allows to estimate from an observed spectrum the proper physical parameters of the absorbing gas and simultaneously an appropriate structure of the velocity field parallel to the line-of-sight. The application to the analysis of the H+D Ly$\alpha$ profile is demonstrated using Burles & Tytler [B&T] data for QSO 1009+2956 where the DI Ly$\alpha$ line is seen at $z_a = 2.504$. The results obtained favor a low D/H ratio in this absorption system, although our upper limit for the hydrogen isotopic ratio of about $4.5\times10^{-5}$ is slightly higher than that of B&T (D/H = $3.0^{+0.6}_{-0.5} \times 10^{-5}$). We also show that the D/H and N(HI) values are, in general, correlated, i.e. the derived D-abundance may be badly dependent on the assumed hydrogen column density. The corresponding confidence regions for an arbitrary and a fixed stochastic velocity field distribution are calculated.

A new method, based on a Reverse Monte Carlo technique and aimed at the inverse problem in the analysis of interstellar (intergalactic) absorption lines is presented. We consider the process of line formation in media with a stochastic velocity field accounting for the effects of a finite correlation length (mesoturbulence). This approach generalizes the standard microturbulent approximation, which is commonly used to model the formation of absorption spectra in turbulent media. The method allows to estimate from an observed spectrum both, the physical parameters of the absorbing gas, and an appropriate structure of the velocity distribution parallel to the line of sight. The procedure is applied to a template H+D Ly-alpha profile which reproduces the Burles and Tytler [BT] spectrum of the quasar Q 1009+2956 with the DI Ly-alpha line seen at the redshift z = 2.504. The results obtained favor a low D/H ratio in this particular absorption system, although the inferred upper limit for the hydrogen isotopic ratio of about 4.6 10^{-5} is slightly higher than that of BT [log (D/H) = -4.60 +/- 0.08 +/- 0.06,i.e. D/H <= 3.47 10^{-5}]. This revision of the BT data in the framework of the mesoturbulent model leads to limits on D/H consistent with standard big bang nucleosynthesis [SBBN] predictions and observational constraints on both extra-galactic 4He mass fraction Y_p (Izotov et al. 1997) and 7Li abundance in the atmospheres of population II (halo) stars (Bonifacio and Molaro 1997). SBBN restricts the common interval of the baryon to photon ratio \eta for which there is concordance between the various abundances to 3.9 10^{-10} <= \eta <= 5.3 10^{-10}. It implies that 0.014 <= \Omega_b h^2 <= 0.020, where h is the Hubble constant in units 100 km/s/Mpc.

A large-scale gradient in the metal abundance has been detected with ASCA from an X-ray bright cluster of galaxies AWM7. The metal abundance shows a peak of 0.5 solar at the center and smoothly declines to <~ 0.2 solar at a radius of 500 kpc. The gas temperature is found to be constant at 3.8 keV. The radial distribution of iron can be fit with a beta-model with beta ~ 0.8 assuming the same core radius (115 kpc) as that of the intracluster medium. The metal distribution in AWM7 suggests that the gas injected from galaxies is not efficiently mixed in the cluster space and traces the distribution of galaxies.

By using a Reverse Monte Carlo technique aimed at the inverse problem in the analysis of the H+D Ly-a absorption from QSO spectra we have estimated the physical parameters of the gas and the large-scale line of sight velocity structure towards Q1009+2956 at z = 2.504 (observations by Burles and Tytler, astro-ph/9603070). We show that for a precise evaluation of the column densities, the detailed velocity field structure is to be determined simultaneously. Our study of the z = 2.504 system yields D/H = (3.75 +/- 0.85)*10^{-5} (2sigma) which gives good agreement between the standard BBN predictions and observational constraints on extra-galactic 4He mass fraction and 7Li abundances in the atmospheres of population II (halo) starts. We conclude that the discordance of D/H with the light element abundances discussed in the literature is a consequence of the assumption of microturbulence (i.e. completely uncorrelated bulk motions) in the analysis of line profiles. The generalized model accounting for a finite correlation length leads to the common eta-interval of (3.9-5.3)*10^{-10}, where eta is the baryon to photon ratio.

Recently, Jakobsen et al. reported that the HeII Gunn-Peterson optical depth is greater than 1.7 at $ z=3.3$. When we consider this observation combined with the HI Gunn-Peterson optical depths and the intensity of the UV background at high redshifts, we can put constraints on the spectrum of the UV background and the mass fraction of intergalactic medium. We find that the spectrum must be steep ($\alpha \gsim 3; J_{\nu} \propto \nu^{-\alpha} $) to be consistent with observations, and that the mass fraction of the intergalactic medium must be greater than 0.3 (0.1) in the case $ \alpha=3 (4)$ for baryon density consistent with the primordial nucleosynthesis. We also find that HeII fraction is as low as $10^{-3} \sim 10^{-1} $ even for this steep spectrum.