The very high energy (VHE) gamma-ray emission reported from a number of pulsar wind nebulae (PWNe) is naturally explained by the inverse Compton (IC) scattering of multi-TeV electrons. On the other hand, the physical dimensions of some gamma-ray emitting PWNe significantly exceed the scales anticipated by the standard hydrodynamical paradigm of PWN formation. The most "disturbing" case in this regard is HESS J1825-137 which extends to distances $r\approx35$ pc from the central pulsar. If the gamma-ray emission is indeed produced inside the plerion, but not by electrons which escaped the nebula and diffuse in the interstellar medium (ISM), the formation of such an anomalously extended plerion could be realized, in an extremely diluted environment with the hydrogen number density \(n_\textsc{ism}\le10^{-2}\rm\,cm^{-3}\). In this paper, we explore an alternative scenario assuming that the pulsar responsible for the formation of the nebula initially had a very short rotation period. In this case, the sizes of both the PWN and the surrounding SNR depend on the initial pulsar period, the braking index, and the ISM density. To check the feasibility of this scenario, we study the parameter space that would reproduce the size of HESS J1825-137. We show that this demand can be achieved if the braking index is small, \(n\leq2\) and the pulsar birth period is short, \(P_{\rm b}\simeq1\rm\,ms\). Remarkably, this scenario can reproduce the size of HESS J1825-137 only in a dense environment with \(n_\textsc{ism}\geq\rm1\,cm^{-3}\). The requirement of the dense surrounding gas is supported by the presence of molecular clouds found in the source vicinity.

High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment.

Stars and their winds can contribute to the non-thermal (NT) emission in extragalactic jets. Given the complexity of jet-star interactions, the properties of the resulting emission are strongly linked to those of the emitting flows. We simulate the interaction between a stellar wind and a relativistic extragalactic jet and use the hydrodynamic results to compute the NT emission under different conditions. We perform relativistic axisymmetric hydrodynamical simulations of a relativistic jet interacting with a supersonic, non-relativistic stellar wind. We compute the corresponding streamlines out of the simulation results, and calculate the injection, evolution, and emission of NT particles accelerated in the jet shock, focusing on electrons or $e^\pm$-pairs. Several cases are explored, considering different jet-star interaction locations, magnetic fields and observer lines of sight. The jet luminosity and star properties are fixed, but the results are easily scalable under changes of these parameters. Individual jet-star interactions produce synchrotron and inverse Compton (IC) emission that peaks from X-rays to MeV energies (depending on the magnetic field), and at $\sim 100-1000$ GeV (depending on the stellar type), respectively. The radiation spectrum is hard in the scenarios explored here due to non-radiative cooling dominance, as low-energy electrons are efficiently advected even under relatively high magnetic fields. Interactions of jets with cold stars lead to an even harder IC spectrum because of the Klein-Nishina effect in the cross-section. Doppler boosting has a strong impact on the observer luminosity. The emission levels for individual interactions found here are in the line of previous, more approximate, estimates, strengthening the hypothesis that collective jet-star interactions could significantly contribute at high energies under efficient particle acceleration.

Binaries hosting a massive star and a non-accreting pulsar are powerful non-thermal emitters due to the interaction of the pulsar and the stellar wind. The winds of massive stars are thought to be inhomogeneous, which could have an impact on the non-thermal emission. We study numerically the impact of the presence of inhomogeneities or clumps in the stellar wind on the high-energy non-thermal radiation of high-mass binaries hosting a non-accreting pulsar. We compute the trajectories and physical properties of the streamlines in the shocked pulsar wind without clumps, with a small clump, and with a large one. This information is used to compute the synchrotron and inverse Compton emission from the non-thermal populations, accounting also for the effect of gamma-ray absorption through pair creation. A specific study is done for PSR B1259-63/LS2883. When stellar wind clumps perturb the two-wind interaction region, the associated non-thermal radiation in the X-ray band,of synchrotron origin, and in the GeV-TeV band, of inverse Compton origin, is affected by several effects: (i) strong changes in the the plasma velocity direction that result in Doppler boosting factor variations, (ii) strengthening of the magnetic field that mainly enhances the synchrotron radiation, (iii) strengthening of the pulsar wind kinetic energy dissipation at the shock, potentially available for particle acceleration, and (iv) changes in the rate of adiabatic losses that affect the lower energy part of the non-thermal particle population. The radiation above 100 GeV detected, presumably, during the post-periastron crossing of the Be star disc in PSR B1259-63/LS2883, can be roughly reproduced assuming that the crossing of the disc is modeled as the encounter with a large inhomogeneity.

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions developed by the Institute of Space and Astronautical Science (ISAS), with a planned launch in 2015. The ASTRO-H mission is equipped with a suite of sensitive instruments with the highest energy resolution ever achieved at E > 3 keV and a wide energy range spanning four decades in energy from soft X-rays to gamma-rays. The simultaneous broad band pass, coupled with the high spectral resolution of Delta E < 7 eV of the micro-calorimeter, will enable a wide variety of important science themes to be pursued. ASTRO-H is expected to provide breakthrough results in scientific areas as diverse as the large-scale structure of the Universe and its evolution, the behavior of matter in the gravitational strong field regime, the physical conditions in sites of cosmic-ray acceleration, and the distribution of dark matter in galaxy clusters at different redshifts.

The inverse Compton (IC) scattering of relativistic electrons is one of the major gamma-ray production mechanisms in different environments. Often the target photons for the IC scattering are dominated by black (or grey) body radiation. In this case, the precise treatment of the characteristics of IC radiation requires numerical integrations over the Planckian distribution. Formally, analytical integrations are also possible but they result in series of several special functions; this limits the efficiency of usage of these expressions. The aim of this work is the derivation of approximate analytical presentations which would provide adequate accuracy for the calculations of the energy spectra of up-scattered radiation, the rate of electron energy losses, and the mean energy of emitted photons. Such formulae have been obtained by merging the analytical asymptotic limits. The coefficients in these expressions are calculated via the least square fitting of the results of numerical integrations. The simple analytical presentations, obtained for both the isotropic and anisotropic target radiation fields, provide adequate (as good as $1\%$) accuracy for broad astrophysical applications.

We develop a Monte Carlo Comptonization model for the X-ray spectrum of accretion-powered pulsars. Simple, spherical, thermal Comptonization models give harder spectra for higher optical depth, while the observational data from Vela X-1 show that the spectra are harder at higher luminosity. This suggests a physical interpretation where the optical depth of the accreting plasma increases with mass accretion rate. We develop a detailed Monte-Carlo model of the accretion flow, including the effects of the strong magnetic field ($\sim 10^{12}$ G) both in geometrically constraining the flow into an accretion column, and in reducing the cross section. We treat bulk-motion Comptonization of the infalling material as well as thermal Comptonization. These model spectra can match the observed broad-band {\it Suzaku} data from Vela X-1 over a wide range of mass accretion rates. The model can also explain the so-called "low state", in which the uminosity decreases by an order of magnitude. Here, thermal Comptonization should be negligible, so the spectrum instead is dominated by bulk-motion Comptonization.

In a small scale turbulent medium, when the nonrelativistic Larmor radius $R_{\rm L}=mc^2/eB$ exceeds the correlation length $\lambda$ of the magnetic field, the magnetic bremsstrahlung of charged relativistic particles unavoidably proceeds in the so-called jitter radiation regime. The cooling timescale of parent particles is identical to the synchrotron cooling time, thus this radiation regime can be produced with very high efficiency in different astrophysical sources characterized by high turbulence. The jitter radiation has distinct spectral features shifted, compared to synchrotron radiation, towards high energies. This makes the jitter mechanism an attractive broad-band gamma-ray production channel which in highly magnetized and turbulent environments can compete or even dominate over other high energy radiation mechanisms. In this paper we present a novel study on spectral properties of the jitter radiation performed within the framework of perturbation theory. The derived general expression for the spectral power of radiation is presented in a compact and convenient for numerical calculations form.

We have analyzed the time variability of the wide-band X-ray spectrum of Vela X-1, the brightest wind-fed accreting neutron star, on a short timescale of 2 ks by using {\it Suzaku} observations with an exposure of 100 ks. During the observation, the object showed strong variability including several flares and so-called "low states", in which the X-ray luminosity decreases by an order of magnitude. Although the spectral hardness increases with the X-ray luminosity, the majority of the recorded flares do not show any significant changes of circumstellar absorption. However, a sign of heavy absorption was registered immediately before one short flare that showed a significant spectral hardening. In the low states, the flux level is modulated with the pulsar spin period, indicating that even at this state the accretion flow reaches the close proximity of the neutron star. Phenomenologically, the broad-band X-ray spectra, which are integrated over the entire spin phase, are well represented by the "NPEX" function (a combination of negative and positive power laws with an exponential cutoff by a common folding energy) with a cyclotron resonance scattering feature at 50 keV. Fitting of the data allowed us to infer a correlation between the photon index and X-ray luminosity. Finally, the circumstellar absorption shows a gradual increase in the orbital phase interval 0.25--0.3, which can be interpreted as an impact of a bow shock imposed by the motion of the compact object in the supersonic stellar wind.

The high and very-high energy spectrum of gamma-ray binaries has become a challenge for all theoretical explanations since the detection of powerful, persistent GeV emission from LS 5039 and LS I +61 303 by Fermi/LAT. The spectral cutoff at a few GeV indicates that the GeV component and the fainter, hard TeV emission above 100 GeV are not directly related. We explore the possible origins of these two emission components in the framework of a young, non-accreting pulsar orbiting the massive star, and initiating the non-thermal emission through the interaction of the stellar and pulsar winds. The pulsar/stellar wind interaction in a compact orbit binary gives rise to two potential locations for particle acceleration: the shocks at the head-on collision of the winds and the termination shock caused by Coriolis forces on scales larger than the binary separation. We explore the suitability of these two locations to host the GeV and TeV emitters, respectively, through the study of their non-thermal emission along the orbit. We focus on the application of this model to LS 5039 given its well determined stellar wind with respect to other gamma-ray binaries. The application of the proposed model to LS 5039 indicates that these two potential emitter locations provide the necessary conditions for reproduction of the two-component high-energy gamma-ray spectrum of LS 5039. In addition, the ambient postshock conditions required at each of the locations are consistent with recent hydrodynamical simulations. The scenario based on the interaction of the stellar and pulsar winds is compatible with the GeV and TeV emission observed from gamma-ray binaries with unknown compact objects, such as LS 5039 and LS I +61 303.

Strong iron fluorescence at 6.4 keV and hard-X-ray emissions from giant molecular clouds in the Galactic center region have been interpreted as reflections of a past outburst of the Sgr A* supermassive black hole. Careful treatment of multiple interactions of photons in a complicated geometry is essential to modeling the reprocessed emissions from the dense clouds. We develop a new calculation framework of X-ray reflection from molecular clouds based on Monte Carlo simulations for accurate interpretation of high-quality observational data. By utilizing this simulation framework, we present the first calculations of morphologies and spectra of the reflected X-ray emission for several realistic models of Sgr B2, which is the most massive molecular cloud in our Galaxy. The morphology of scattered hard X-rays above 20 keV is significantly different from that of iron fluorescence due to their large penetrating power into dense regions of the cloud, probing the structure of the cloud. High-resolution spectra provide quantitative evaluation of the iron line including its Compton shoulder to constrain the mass and the chemical composition of the cloud as well as the luminosity of the illuminating source. These predictions can be checked in the near future with future X-ray missions such as NuStar (hard X-rays) and ASTRO-H (both iron lines and hard X-rays).

We study the applicability of the idea of internal absorption of gamma-rays produced through synchrotron radiation of ultrarelativistic protons in highly magnetized blobs to 1ES 0229+200 and 3C 66A, the two TeV blazars which show unusually hard intrinsic gamma-ray spectra after being corrected for the intergalactic absorption. We show that for certain combinations of reasonable model parameters, even with quite modest energy requirements, the scenario allows a self-consistent explanation of the non-thermal emission of these objects in the keV, GeV, and TeV energy bands.

Several binary systems that contain a massive star have been detected in both the radio band and at very high energies. In the dense stellar photon field of these sources, gamma-ray absorption and pair creation are expected to occur, and the radiation from these pairs may contribute significantly to the observed radio emission. We aim at going deeper in the study of the properties, and in particular the morphology, of the pair radio emission in gamma-ray binaries. We apply a Monte-Carlo code that computes the creation location, the spatial trajectory and the energy evolution of the pairs produced in the binary system and its surroundings. The radio emission produced by these pairs, with its spectral, variability and spatial characteristics, is calculated as it would be seen from a certain direction. A generic case is studied first, and then the specific case of LS 5039 is also considered. We find that, confirming previous results, the secondary radio emission should appear as an extended radio structure of a few milliarcseconds size. This radiation would be relatively hard, with fluxes up to ~ 10 mJy. Modulation is expected depending on the gamma-ray production luminosity, system eccentricity, and wind ionization fraction, and to a lesser extent on the magnetic field structure. In gamma-ray binaries in general, the pairs created due to photon-photon interactions can contribute significantly to the core, and generate an extended structure. In the case of LS 5039, the secondary radio emission is likely to be a significant fraction of the detected core flux, with a marginal extension.

Binary pulsar systems emit potentially detectable components of gamma ray emission due to Comptonization of the optical radiation of the companion star by relativistic electrons of the pulsar wind, both before and after termination of the wind. The recent optical observations of binary pulsar system PSR B1259-63/LS 2883 revealed radiation properties of the companion star which differ significantly from previous measurements. In this paper we study the implications of these observations for the interaction rate of the unshocked pulsar wind with the stellar photons and the related consequences for fluxes of high energy (HE) and very high energy (VHE) gamma rays. We show that the signal should be strong enough to be detected with Fermi close to the periastron passage, unless the pulsar wind is strongly anisotropic or the Lorentz factor of the wind is smaller than $10^3$ or larger that $10^5$. The higher luminosity of the optical star also has two important implications: (i) attenuation of gamma rays due to photon-photon pair production, and (ii) Compton drag of the unshocked wind. While the first effect has an impact on the lightcurve of VHE gamma rays, the second effect may significantly decrease the energy available for particle acceleration after termination of the wind.

Only a few binary systems with compact objects display TeV emission. The physical properties of the companion stars represent basic input to understand the physical mechanisms behind the particle acceleration, emission, and absorption processes in these so-called gamma-ray binaries. Here we present high-resolution and high signal-to-noise optical spectra of LS2883, the Be star forming a gamma-ray binary with the young non-accreting pulsar PSR B1259-63, showing it to rotate faster and be significantly earlier and more luminous than previously thought. Analysis of the interstellar lines suggest that the system is located at the same distance as (and thus is likely a member of) CenOB1. Taking the distance to the association, d=2.3kpc, and a color excess of E(B-V)=0.85 for LS2883, results in Mv=-4.4. Because of fast rotation, LS2883 is oblate (R_eq=9.7R_sun and R_pole=8.1R_sun) and presents a temperature gradient (T_eq=27500K, log g_eq=3.7; T_pole=34000K, log g_pole=4.1). If the star did not rotate, it would have parameters corresponding to a late O-type star. We estimate its luminosity at log(L*/L_sun)=4.79, and its mass at M=30M_sun. The mass function then implies an inclination of the binary system i_orb=23deg,slightly smaller than previous estimates. We discuss the implications of these new astrophysical parameters of LS2883 for the production of high energy and very high energy gamma rays in the PSR B1259-63/LS2883 gamma-ray binary system. In particular, the stellar properties are very important for prediction of the line-like bulk Comptonization component from the unshocked ultra-relativistic pulsar wind.

Microquasars are X-ray binaries with relativistic jets. These jets are powerful energy carriers, thought to be fed by accretion, which produce non-thermal emission at different energy bands. To date, several Galactic sources showing extended radio emission, among them at least one confirmed microquasar, Cygnus X-1, have been detected in the TeV range. All of them show complex patterns of spectral and temporal behavior. In this work, we discuss the physics behind the very high-energy emission in microquasars. In concrete, we focus on the microquasar Cygnus X-1, and also in the other two TeV binaries with detected extended outflows, LS 5039 and LS I +61 303, pointing out relevant aspects of the complex phenomena occurring in them. We conclude that the TeV emission is likely of leptonic origin, although hadrons cannot be discarded. In addition, efficient electromagnetic cascades can hardly develop since even relatively low magnetic fields suppress them. Also, the modeling of the radiation from some of the detected sources points to them as either extremely efficient accelerators, and/or having the TeV emitter at a distance from the compact object of about 10e12 cm. Finally, we point out that the role of a massive and hot stellar companion, due to its strong photon field and wind, cannot be neglected when trying to understand the behavior of microquasars at high and very high energies. The complexity of microquasars precludes straightforward generalizations to a whole population, and are better studied presently in a source by source base.

LS 5039 is an X-ray binary detected at very high energies. Along the orbit, there is a significant detection even during the superior conjunction of the compact object, when very large gamma-ray opacities are expected. Electromagnetic cascades, which may make the system more transparent to gamma-rays, are hardly efficient for reasonable magnetic fields in the massive star surroundings. A jet-like flow could transport energy to regions where the photon-photon absorption is much lower and the TeV radiation is not so severely absorbed. Otherwise, in the standard pulsar scenario for LS 5039, the emitter would be located between the star and the compact object, which would imply the violation of the observational constraints at X-rays.

The recent detections of TeV gamma-rays from compact binary systems show that relativistic outflows (jets or winds) are sites of effective acceleration of particles up to multi-TeV energies. In this paper, we discuss the conditions of acceleration and radiation of ultra-relativistic electrons in LS 5039, the gamma-ray emitting binary system for which the highest quality TeV data are available. Assuming that the gamma-ray emitter is a jet-like structure, we performed detailed numerical calculations of the energy spectrum and lightcurves accounting for the acceleration efficiency, the location of the accelerator, the speed of the emitting flow, the inclination angle of the system, as well as specific features related to anisotropic inverse Compton scattering and pair production. We conclude that the accelerator should not be deep inside the binary system unless we assume a very efficient acceleration rate. We show that within the IC scenario both the gamma-ray spectrum and flux are strongly orbital phase dependent. Formally, our model can reproduce, for specific sets of parameter values, the energy spectrum of gamma-rays reported by HESS for wide orbital phase intervals. However, the physical properties of the source can be constrained only by observations capable of providing detailed energy spectra for narrow orbital phase intervals ($\Delta\phi\ll 0.1$).

The inverse Compton scattering of ultrarelativistic electrons accelerated at the pulsar wind termination shock is believed to be responsible for TeV gamma-ray signal recently reported from the binary system PSR B1259-63/SS2883. While this process can explain the energy spectrum of the observed TeV emission, the detected gamma-ray fluxes do not agree with the published theoretical predictions of the TeV lightcurve. The main objective of this paper is to show that the HESS results can be explained, under certain reasonable assumptions, by inverse Compton scenarios of gamma-ray production in the system. In this paper we study evolution of the energy spectra of relativistic electrons under different assumptions about the acceleration and energy-loss rates of electrons, and the impact of these processes on the lightcurve of IC gamma-rays. We demonstrate that the observed TeV lightcurve can be explained (i) by adiabatic losses which dominate over the entire trajectory of the pulsar with a significant increase towards the periastron, or (ii) by the sub-TeV cutoffs in the energy spectra of electrons due to the enhanced rate of Compton losses close to the periastron. The Compton deceleration of the pulsar wind contributes to the decrease of the nonthermal power released in the accelerated electrons, and thus to the reduction of the IC and synchrotron components of radiation close to the periastron. Although this effect alone cannot explain the observed lightcurves, the Comptonization of the pulsar wind leads to the formation of gamma-radiation with a line-type energy spectrum. While the HESS data already constrain the Lorentz factor of the wind, $\Gamma \le 10^6$, future observations of this object with GLAST should allow a deep probe of the wind Lorentz factor in the range between $10^4$ and $10^6$.

The recent detection of TeV gamma-rays from the microquasar LS 5039 by HESS is one of the most exciting discoveries of observational gamma-ray astronomy in the very high energy regime. This result clearly demonstrates that X-ray binaries with relativistic jets (microquasars) are sites of effective acceleration of particles (electrons and/or protons) to multi-TeV energies. Whether the gamma-rays are of hadronic or leptonic origin is a key issue related to the origin of Galactic Cosmic Rays. We discuss different possible scenarios for the production of gamma-rays, and argue in favor of hadronic origin of TeV photons, especially if they are produced within the binary system. If so, the detected gamma-rays should be accompanied by a flux of high energy neutrinos emerging from the decays of \pi^\pm mesons produced at pp and/or p \gamma interactions. The flux of TeV neutrinos, which can be estimated on the basis of the detected TeV gamma-ray flux, taking into account the internal \gamma \gamma \to e^+e^- absorption, depends significantly on the location of gamma-ray production region(s). The minimum neutrino flux above 1 TeV is expected to be at the level of 10^{-12} cm^{-2} s^{-1}; however, it could be up to a factor of 100 larger. The detectability of the signal of multi-TeV neutrinos significantly depends on the high energy cutoff in the spectrum of parent protons; if the spectrum of accelerated protons continues to 1 PeV and beyond, the predicted neutrino fluxes can be probed by the planned km^3-scale neutrino detector.