Hot, tenuous advection-dominated accretion flows around black holes are ideal sites for the Fermi acceleration of relativistic particles at standing shock waves in the accretion disc. Previous work has demonstrated that the shock-acceleration process can be efficient enough to power the observed, strong outflows in radio-loud active galaxies such as M87. However, the dynamical effect (back-reaction) on the flow, exerted by the pressure of the relativistic particles, has not been previously considered, and this effect can have a significant influence on the disc structure. We reexamine the problem by developing a new, two-fluid model for the structure of the accretion disc that includes the dynamical effect of the relativistic particle pressure, combined with the pressure of the background (thermal) gas. The new model is analogous to the two-fluid model of cosmic ray acceleration in supernova-driven shock waves. As part of the model, we also develop a new set of shock jump conditions, which are solved along with the hydrodynamic conservation equations to determine the structure of the accretion disc. The solutions include the formation of a mildly relativistic outflow (jet) at the shock radius, driven by the relativistic particles accelerated in the disc. One of our main conclusions is that in the context of the new two-fluid accretion model, global smooth (shock-free) solutions do not exist, and the disc must always contain a standing shock wave, at least in the inviscid case considered here.

The multi-wavelength spectrum observed from M87 extends from radio wavelengths up to TeV gamma-ray energies. The radio through GeV components have been interpreted successfully using SSC models based on misaligned blazar jets, but the origin of the intense TeV emission detected during flares in 2004, 2005, and 2010 remains puzzling. It has been previously suggested that the TeV flares are produced when a relativistic proton jet originating in the core of M87 collides with a molecular cloud (or stellar atmosphere) located less than one parsec from the central black hole. We explore this scenario in detail here using a self-consistent model for the acceleration of relativistic protons in a shocked, two-fluid ADAF accretion disc. The relativistic protons accelerated in the disc escape to power the observed jet outflows. The distribution function for the jet protons is used to compute the TeV emission produced when the jet collides with a cloud or stellar atmosphere. The simulated broadband radiation spectrum includes radio, X-ray, and GeV components generated via synchrotron, as well as TeV emission generated via the production and decay of muons, positrons, and electrons. The self-consistency of the model is verified by computing the relativistic particle pressure using the distribution function, and comparing it with the relativistic particle pressure obtained from the hydrodynamical model. We demonstrate that the model is able to reproduce the multi-wavelength spectrum from M87 observed by VERITAS and HESS during the high-energy flares in 2004, 2005, and 2010.

Multiwavelength monitoring of the blazar 3C 279 observed a very bright, 12-hour, orphan gamma-ray flare on 20 Dec 2013 with a uniquely hard Fermi-LAT spectrum and high Compton dominance. We work with a one-zone, leptonic model with both first- and second-order Fermi acceleration, which now reproduces the unique flaring behavior. We present a simplified analytic electron energy distribution to provide intuition about how particle acceleration shapes multi-wavelength blazar jet emission spectra. The contributions of individual processes in relativistic jets is fundamental to understanding the particle energy budget in the formation and propagation of astrophysical jets. We show that first- and second-order Fermi acceleration are sufficient to explain the flare, and that magnetic reconnection is not needed. Our analysis suggests that the flare is initiated by an increase in the particle energies due to shock acceleration, which also increases the stochastic acceleration. The higher energy particle preferentially occupy the outer jet, along the sheath, which decreases the apparent magnetic field and synchrotron radiation, while increasing electron exposure to the broad line region photon fields, driving up the external Compton emission.

We present an analysis of X-ray spectra of the high mass X-ray binary 4U 0115+634 as observed with Suzaku and RXTE in 2011 July, during the fading phase of a giant X-ray outburst. We used a continuum model consisting of an absorbed cutoff power-law and an ad-hoc Gaussian emission feature centered around 8.5 keV, which we discuss to be due to cyclotron emission. Our results are consistent with a fundamental cyclotron absorption line centered at ${\sim}10.2$ keV for all observed flux ranges. At the same time we rule out significant influence of the 8.5 kev Gaussian on the CRSF parameters, which are not consistent with the cyclotron line energies and depths of previously reported flux-dependent descriptions. We also show that some continuum models can lead to artificial line-like residuals in the analyzed spectra, which are then misinterpreted as unphysically strong cyclotron lines. Specifically, our results do not support the existence of a previously claimed additional cyclotron feature at ${\sim}15$ keV. Apart from these features, we find for the first time evidence for a He-like Fe XXV emission line at ${\sim}6.7$ keV and weak H-like Fe XXVI emission close to ${\sim}7.0$ keV.

The broadband spectrum from the 2013 December 20 $\gamma$-ray flare from 3C~279 is analyzed with our previously-developed one-zone blazar jet model. We are able to reproduce two SEDs, a quiescent and flaring state, the latter of which had an unusual SED, with hard $\gamma$-ray spectrum, high Compton dominance, and short duration. Our model suggests that there is insufficient energy for a comparable X-ray flare to have occurred simultaneously, which is an important constraint given the lack of X-ray data. We show that first- and second-order Fermi acceleration are sufficient to explain the flare, and that magnetic reconnection is not needed. The model includes particle acceleration, escape, and adiabatic and radiative energy losses, including the full Compton cross-section, and emission from the synchrotron, synchrotron self-Compton, and external Compton processes. We provide a simple analytic approximation to the electron distribution solution to the transport equation that may be useful for simplified modeling in the future.

Studying the physical processes occurring in the region just above the magnetic poles of strongly magnetized, accreting binary neutron stars is essential to our understanding of stellar and binary system evolution. Perhaps more importantly, it provides us with a natural laboratory for studying the physics of high temperature and high density plasmas exposed to extreme radiation, gravitational, and magnetic fields. Observations over the past decade have shed new light on the manner in which plasma falling at velocities near the speed of light onto a neutron star surface is halted. Recent advances in modeling these processes have resulted in direct measurement of the magnetic fields and plasma properties. On the other hand, numerous physical processes have been identified that challenge our current picture of how the accretion process onto neutron stars works. Observation and theory are our essential tools in this regime because the extreme conditions cannot be duplicated on Earth. This white paper gives an overview of the current theory, the outstanding theoretical and observational challenges, and the importance of addressing them in contemporary astrophysics research.

In 2007 September, the Crab Nebula exhibited a bright gamma-ray flare in the GeV energy range that was detected by AGILE. The observed emission at >160 MeV indicates that the radiating electrons had energies above the classical synchrotron radiation-reaction limit, thus presenting a serious challenge to classical models for electron acceleration in astrophysical environments. In this paper, we apply our recently developed time-dependent self-similar analytical model describing electrostatic acceleration in the explosive reconnection region around the pulsar wind termination shock to the 2007 September flare. This event was unique in that it displayed both long-duration "wave" and short-duration "sub-flare" features. The unusual temporal variation makes this flare an especially interesting test for our model. We demonstrate that our model can reproduce the time-dependent gamma-ray spectrum for this event, as well as the associated gamma-ray light curve, obtained by integrating the spectrum for photon energies >100 MeV. This establishes that our time-dependent electrostatic acceleration model can explain both wave and sub-flare transients, which lends further support to the theoretical framework we have developed. We also further examine the validity of the self-similar electric and magnetic field evolution implied by our model. We conclude that strong electrostatic acceleration driven by shock-induced magnetic reconnection is able to power the Crab Nebula gamma-ray flares by energizing the electrons on sub-Larmor timescales.

The $\gamma$-ray flares from the Crab nebula observed by {\it AGILE} and {\it Fermi}-LAT between 2007-2013 reached GeV photon energies and lasted several days. The strongest emission, observed during the 2011 April "super-flare," exceeded the quiescent level by more than an order of magnitude. These observations challenge the standard models for particle acceleration in pulsar wind nebulae, because the radiating electrons have energies exceeding the classical radiation-reaction limit for synchrotron. Particle-in-cell simulations have suggested that the classical synchrotron limit can be exceeded if the electrons also experience electrostatic acceleration due to shock-driven magnetic reconnection. In this paper, we revisit the problem using an analytic approach based on solving a fully time-dependent electron transport equation describing the electrostatic acceleration, synchrotron losses, and escape experienced by electrons in a magnetically confined plasma "blob" as it encounters and passes through the pulsar-wind termination shock. We show that our model can reproduce the $\gamma$-ray spectra observed during the rising and decaying phases of each of the two sub-flare components of the 2011 April super-flare. We integrate the spectrum for photon energies $\ge 100\,$MeV to obtain the light curve for the event, which agrees with the observations. We find that strong electrostatic acceleration occurs on both sides of the termination shock, driven by magnetic reconnection. We also find that the dominant mode of particle escape changes from diffusive escape to advective escape as the blob passes through the shock.

We introduce a new theoretical model to describe the emitting region in a blazar jet. We assume a one-zone leptonic picture, and construct the particle transport equation for a plasma blob experiencing low-energy, monoenergetic particle injection, energy dependent particle escape, shock acceleration, adiabatic expansion, stochastic acceleration, synchrotron radiation, and external Compton radiation from the dust torus and broad line region. We demonstrate that a one-zone leptonic model is able to explain the IR though {\gamma}-ray spectrum for 3C 279 in 2008-2009. We determine that the broad-line region seed photons cannot be adequately described by a single average distribution, but rather we find that a stratified broad line region provides an improvement in the estimation of the distance of the emitting region from the black hole. We calculate that the jet is not always in equipartition between the particles and magnetic field, and find that stochastic acceleration provides more energy to the particles than does shock acceleration, where the latter is also overshadowed by adiabatic losses. We further introduce a novel technique to implement numerical boundary conditions and determine the global normalization for the electron distribution, based on analysis of stiff ordinary differential equations. Our astrophysical results are compared with those obtained by previous authors.

Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons off quantized plasma electrons in the strong magnetic field (of the order 10^12 G) close to the surface of an accreting X-ray pulsar. The line profiles of CRSFs cannot be described by an analytic expression. Numerical methods such as Monte Carlo (MC) simulations of the scattering processes are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time consuming simulations possible on convenient time scales. Furthermore, it can generate responses to mono-energetic photon injections, producing Green's functions, which can be used later to generate spectra for arbitrary continua. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green's function tables are available online at http://www.sternwarte.uni-erlangen.de/research/cyclo .

The availability of the unprecedented spectral resolution provided by modern X-ray observatories is opening up new areas for study involving the coupled formation of the continuum emission and the cyclotron absorption features in accretion-powered X-ray pulsar spectra. Previous research focusing on the dynamics and the associated formation of the observed spectra has largely been confined to the single-fluid model, in which the super-Eddington luminosity inside the column decelerates the flow to rest at the stellar surface, while the dynamical effect of gas pressure is ignored. In a companion paper, we have presented a detailed analysis of the hydrodynamic and thermodynamic structure of the accretion column obtained using a new self-consistent model that includes the effects of both gas and radiation pressure. In this paper, we explore the formation of the associated X-ray spectra using a rigorous photon transport equation that is consistent with the hydrodynamic and thermodynamic structure of the column. We use the new model to obtain phase-averaged spectra and partially-occulted spectra for Her X-1, Cen X-3, and LMC X-4. We also use the new model to constrain the emission geometry, and compare the resulting parameters with those obtained using previously published models. Our model sheds new light on the structure of the column, the relationship between the ionized gas and the photons, the competition between diffusive and advective transport, and the magnitude of the energy-averaged cyclotron scattering cross section.

Previous research centered on the hydrodynamics in X-ray pulsar accretion columns has largely focused on the single-fluid model, in which the super-Eddington luminosity inside the column decelerates the flow to rest at the stellar surface. This type of model has been relatively successful in describing the overall properties of the accretion flows, but it does not account for the possible dynamical effect of the gas pressure. On the other hand, the most successful radiative transport models for pulsars generally do not include a rigorous treatment of the dynamical structure of the column, instead assuming an ad hoc velocity profile. In this paper, we explore the structure of X-ray pulsar accretion columns using a new, self-consistent, "two-fluid" model, which incorporates the dynamical effect of the gas and radiation pressure, the dipole variation of the magnetic field, the thermodynamic effect of all of the relevant coupling and cooling processes, and a rigorous set of physical boundary conditions. The model has six free parameters, which we vary in order to approximately fit the phase-averaged spectra in Her X-1, Cen X-3, and LMC X-4. In this paper, we focus on the dynamical results, which shed new light on the surface magnetic field strength, the inclination of the magnetic field axis relative to the rotation axis, the relative importance of gas and radiation pressure, and the radial variation of the ion, electron, and inverse-Compton temperatures. The results obtained for the X-ray spectra are presented in a separate paper.

The {\gamma}-ray flares from the Crab nebula observed by AGILE and Fermi-LAT reaching GeV energies and lasting several days challenge the standard models for particle acceleration in pulsar wind nebulae, because the radiating electrons have energies exceeding the classical radiation-reaction limit for synchrotron. Previous modeling has suggested that the synchrotron limit can be exceeded if the electrons experience electrostatic acceleration, but the resulting spectra do not agree very well with the data. As a result, there are still some important unanswered questions about the detailed particle acceleration and emission processes occurring during the flares. We revisit the problem using a new analytical approach based on an electron transport equation that includes terms describing electrostatic acceleration, stochastic wave-particle acceleration, shock acceleration, synchrotron losses, and particle escape. An exact solution is obtained for the electron distribution, which is used to compute the associated {\gamma}-ray synchrotron spectrum. We find that in our model the {\gamma}-ray flares are mainly powered by electrostatic acceleration, but the contributions from stochastic and shock acceleration play an important role in producing the observed spectral shapes. Our model can reproduce the spectra of all the Fermi-LAT and AGILE flares from the Crab nebula, using magnetic field strengths in agreement with the multi-wavelength observational constraints. We also compute the spectrum and duration of the synchrotron afterglow created by the accelerated electrons, after they escape into the region on the downstream side of the pulsar wind termination shock. The afterglow is expected to fade over a maximum period of about three weeks after the {\gamma}-ray flare.

Electron cyclotron resonant scattering features (CRSFs) are observed as absorption-like lines in the spectra of X-ray pulsars. A significant fraction of the computing time for Monte Carlo simulations of these quantum mechanical features is spent on the calculation of the mean free path for each individual photon before scattering, since it involves a complex numerical integration over the scattering cross section and the (thermal) velocity distribution of the scattering electrons. We aim to numerically calculate interpolation tables which can be used in CRSF simulations to sample the mean free path of the scattering photon and the momentum of the scattering electron. The tables also contain all the information required for sampling the scattering electron's final spin. The tables were calculated using an adaptive Simpson integration scheme. The energy and angle grids were refined until a prescribed accuracy is reached. The tables are used by our simulation code to produce artificial CRSF spectra. The electron momenta sampled during these simulations were analyzed and justified using theoretically determined boundaries. We present a complete set of tables suited for mean free path calculations of Monte Carlo simulations of the cyclotron scattering process for conditions expected in typical X-ray pulsar accretion columns (0.01<B/B_{crit}<=0.12, where B_{crit}=4.413x10^{13} G and 3keV<=kT<15keV). The sampling of the tables is chosen such that the results have an estimated relative error of at most 1/15 for all points in the grid. The tables are available online at http://www.sternwarte.uni-erlangen.de/research/cyclo.

We report new spectral modeling of the accreting X-ray pulsar Hercules X- 1. Our radiation-dominated radiative shock model is an implementation of the analytic work of Becker & Wolff on Comptonized accretion flows onto magnetic neutron stars. We obtain a good fit to the spin-phase averaged 4 to 78 keV X-ray spectrum observed by the Nuclear Spectroscopic Telescope Array during a main- on phase of the Her X-1 35-day accretion disk precession period. This model allows us to estimate the accretion rate, the Comptonizing temperature of the radiating plasma, the radius of the magnetic polar cap, and the average scattering opacity parameters in the accretion column. This is in contrast to previous phenomenological models that characterized the shape of the X-ray spectrum but could not determine the physical parameters of the accretion flow. We describe the spectral fitting details and discuss the interpretation of the accretion flow physical parameters.

Depending on the values of the energy and angular momentum per unit mass in the gas supplied at large radii, inviscid advection-dominated accretion flows can display velocity profiles with either preshock deceleration or preshock acceleration. Nakayama has shown that these two types of flow configurations are expected to have different stability properties. By employing the Chevalier \& Imamura linearization method and the Nakayama instability boundary conditions, we discover that there are regions of parameter space where disks/shocks with outflows can be stable or unstable. In regions of instability, we find that preshock deceleration is always unstable to the zeroth mode with zero frequency of oscillation, but is always stable to the fundamental mode and overtones. Furthermore, we also find that preshock acceleration is always unstable to the zeroth mode and that the fundamental mode and overtones become increasingly less stable as the shock location moves away from the horizon when the disk half-height expands above $\sim 12$ gravitational radii at the shock radius. In regions of stability, we demonstrate the zeroth mode to be stable for the velocity profiles that exhibit preshock acceleration and deceleration. Moreover, for models that are linearly unstable, our model suggests the possible existence of quasi-periodic oscillations (QPOs) with ratios 2:3 and 3:5. These ratios are believed to occur in stellar and supermassive black hole candidates, for example, in GRS 1915+105 and Sgr A*, respectively. We expect that similar QPO ratios also exist in regions of stable shocks.

Electromagnetic radiation from blazar jets often displays strong variability, extending from radio to $\gamma$-ray frequencies. In a few cases, this variability has been characterized using Fourier time lags, such as those detected in the X-rays from Mrk~421 using BeppoSAX. The lack of a theoretical framework to interpret the data has motivated us to develop a new model for the formation of the X-ray spectrum and the time lags in blazar jets based on a transport equation including terms describing stochastic Fermi acceleration, synchrotron losses, shock acceleration, adiabatic expansion, and spatial diffusion. We derive the exact solution for the Fourier transform of the electron distribution, and use it to compute the Fourier transform of the synchrotron radiation spectrum and the associated X-ray time lags. The same theoretical framework is also used to compute the peak flare X-ray spectrum, assuming that a steady-state electron distribution is achieved during the peak of the flare. The model parameters are constrained by comparing the theoretical predictions with the observational data for Mrk~421. The resulting integrated model yields, for the first time, a complete first-principles physical explanation for both the formation of the observed time lags and the shape of the peak flare X-ray spectrum. It also yields direct estimates of the strength of the shock and the stochastic MHD wave acceleration components in the Mrk~421 jet.

Many accreting black holes manifest time lags during outbursts, in which the hard Fourier component typically lags behind the soft component. Despite decades of observations of this phenomenon, the underlying physical explanation for the time lags has remained elusive, although there are suggestions that Compton reverberation plays an important role. However, the lack of analytical solutions has hindered the interpretation of the available data. In this paper, we investigate the generation of X-ray time lags in Compton scattering coronae using a new mathematical approach based on analysis of the Fourier-transformed transport equation. By solving this equation, we obtain the Fourier transform of the radiation Green's function, which allows us to calculate the exact dependence of the time lags on the Fourier frequency, for both homogeneous and inhomogeneous coronal clouds. We use the new formalism to explore a variety of injection scenarios, including both monochromatic and broadband (bremsstrahlung) seed photon injection. We show that our model can successfully reproduce both the observed time lags and the time-averaged (quiescent) X-ray spectra for Cyg~X-1 and GX~339-04, using a single set of coronal parameters for each source. The time lags are the result of impulsive bremsstrahlung injection occurring near the outer edge of the corona, while the time-averaged spectra are the result of continual distributed injection of soft photons throughout the cloud.

Accreting highly magnetized pulsars in binary systems are among the brightest X-ray emitters in our Galaxy. Although a number of high statistical quality broad-band (0.1-100 keV) X-ray observations are available, the spectral energy distribution of these sources is usually investigated by adopting pure phenomenological models, rather than models linked to the physics of accretion. In this paper, a detailed spectral study of the X-ray emission recorded from the high-mass X-ray binary pulsars Cen X-3, 4U 0115+63, and Her X-1 is carried out by using BeppoSAX and joined Suzaku+NuStar data, together with an advanced version of the compmag model. The latter provides a physical description of the high energy emission from accreting pulsars, including the thermal and bulk Comptonization of cyclotron and bremsstrahlung seed photons along the neutron star accretion column. The compmag model is based on an iterative method for solving second-order partial differential equations, whose convergence algorithm has been improved and consolidated during the preparation of this paper. Our analysis shows that the broad-band X-ray continuum of all considered sources can be self-consistently described by the compmag model. The cyclotron absorption features, not included in the model, can be accounted for by using Gaussian components. From the fits of the compmag model to the data we inferred the physical properties of the accretion columns in all sources, finding values reasonably close to those theoretically expected according to our current understanding of accretion in highly magnetized neutron stars. The updated version of the compmag model has been tailored to the physical processes that are known to occur in the columns of highly magnetized accreting neutron stars and it can thus provide a better understanding of the high energy radiation from these sources.

We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 d orbit with 2-3 outbursts per orbit during phases of activity. We improve the orbital solution using data from multiple instruments. The X-ray spectrum can be described by an absorbed Fermi-Dirac cutoff power law model along with a narrow Fe K line at 6.4 keV and a weak Cyclotron Resonance Scattering Feature (CRSF) at ~35 keV. The Suzaku data are consistent with the previously observed continuum flux versus iron line flux correlation expected from fluorescence emission along the line of sight. However, the observed iron line flux is slightly higher, indicating the possibility of a higher iron abundance or the presence of non-uniform material. We argue that the source most likely has only been observed in the subcritical (non-radiation dominated) state since its pulse profile is stable over all observed luminosities and the energy of the CRSF is approximately the same at the highest (~$5 \times 10^{37} $erg s$^{-1}$) and lowest (~$5 \times 10^{36} $erg s$^{-1}$) observed 3-60 keV luminosities.

The strong variability of blazars can be characterized by power spectral densities (PSDs) and Fourier frequency-dependent time lags. In previous work, we created a new theoretical formalism for describing the PSDs and time lags produced via a combination of stochastic particle injection and emission via the synchrotron, synchrotron self-Compton, and external Compton (EC) processes. This formalism used the Thomson cross section and simple $\delta$-function approximations to model the synchrotron and Compton emissivities. Here we expand upon this work, using the full Compton cross section and detailed and accurate emissivities. Our results indicate good agreement between the PSDs computed using the $\delta$-function approximations and those computed using the accurate expressions, provided the observed photons are produced primarily by electrons with energies exceeding the lower limit of the injected particle population. Breaks are found in the PSDs at frequencies corresponding to the cooling timescales of the electrons primarily responsible for the observed emission, and the associated time lags are related to the difference in electron cooling timescales between the two energy channels, as expected. If the electron cooling timescales can be determined from the observed time lags and/or the observed EC PSDs, then one could in principle use the method developed here to determine the energy of the external seed photon source for EC, which is an important unsolved problem in blazar physics.

This is an overview of pulsar accretion modeling. The physics of pulsar accretion, i.e., the process of plasma flow onto the neutron star surface, can be constrained from the spectral properties of the X-ray source. We discuss a new implementation of the physical continuum model developed by Becker and Wolff (2007, ApJ 654, 435). The model incorporates Comptonized blackbody, bremsstrahlung, and cyclotron emission. We discuss preliminary results of applying the new tool to the test cases of Suzaku data of Cen X-3 and XTE J1946+274. Cen X-3 is a persistent accreting pulsar with an O-star companion observed during a bright period. XTE J1946+274 is a transient accreting pulsar with a Be companion observed during a dim period. Both sources show spectra that are well described with an empirical Fermi Dirac cutoff power law model. We extend the spectral analysis by making the first steps towards a physical description of Cen X-3 and XTE J1946+274.

Blazars display strong variability on multiple timescales and in multiple radiation bands. Their variability is often characterized by power spectral densities (PSDs) and time lags plotted as functions of the Fourier frequency. We develop a new theoretical model based on the analysis of the electron transport (continuity) equation, carried out in the Fourier domain. The continuity equation includes electron cooling and escape, and a derivation of the emission properties includes light travel time effects associated with a radiating blob in a relativistic jet. The model successfully reproduces the general shapes of the observed PSDs and predicts specific PSD and time lag behaviors associated with variability in the synchrotron, synchrotron self-Compton (SSC), and external Compton (EC) emission components, from sub-mm to gamma-rays. We discuss applications to BL Lacertae objects and to flat-spectrum radio quasars (FSRQs), where there are hints that some of the predicted features have already been observed. We also find that FSRQs should have steeper PSD power-law indices than BL Lac objects at Fourier frequencies < 10^{-4} Hz, in qualitative agreement with previously reported observations by the Fermi Large Area Telescope.

Accreting black hole sources show a wide variety of rapid time variability, including the manifestation of time lags during X-ray transients, in which a delay (phase shift) is observed between the Fourier components of the hard and soft spectra. Despite a large body of observational evidence for time lags, no fundamental physical explanation for the origin of this phenomenon has been presented. We develop a new theoretical model for the production of X-ray time lags based on an exact analytical solution for the Fourier transform describing the diffusion and Comptonization of seed photons propagating through a spherical corona. The resulting Green's function can be convolved with any source distribution to compute the associated Fourier transform and time lags, hence allowing us to explore a wide variety of injection scenarios. We show that thermal Comptonization is able to self-consistently explain both the X-ray time lags and the steady-state (quiescent) X-ray spectrum observed in the low-hard state of Cyg X-1. The reprocessing of bremsstrahlung seed photons produces X-ray time lags that diminish with increasing Fourier frequency, in agreement with the observations for a wide range of sources.

Context: Accretion-powered X-ray pulsars show highly energy-dependent and complex pulse-profile morphologies. Significant deviations from the average pulse profile can appear, in particular close to the cyclotron line energies. These deviations can be described as energy-dependent phase lags, that is, as energy-dependent shifts of main features in the pulse profile. Aims: Using a numerical study we explore the effect of cyclotron resonant scattering on observable, energy-resolved pulse profiles. Methods: We generated the observable emission as a function of spin phase, using Monte Carlo simulations for cyclotron resonant scattering and a numerical ray-tracing routine accounting for general relativistic light-bending effects on the intrinsic emission from the accretion columns. Results: We find strong changes in the pulse profile coincident with the cyclotron line energies. Features in the pulse profile vary strongly with respect to the average pulse profile with the observing geometry and shift and smear out in energy additionally when assuming a non-static plasma. Conclusions: We demonstrate how phase lags at the cyclotron energies arise as a consequence of the effects of angular redistribution of X-rays by cyclotron resonance scattering in a strong magnetic field combined with relativistic effects. We also show that phase lags are strongly dependent on the accretion geometry. These intrinsic effects will in principle allow us to constrain a system's accretion geometry.

We report on an outburst of the high mass X-ray binary 4U 0115+63 with a pulse period of 3.6s in 2008 March/April as observed with RXTE and INTEGRAL. During the outburst the neutron star's luminosity varied by a factor of 10 in the 3--50\,keV band. In agreement with earlier work we find evidence for five cyclotron resonance scattering features at ~10.7, 21.8, 35.5, 46.7, and 59.7keV. Previous work had found an anti-correlation between the fundamental cyclotron line energy and the X-ray flux. We show that this apparent anti-correlation is probably due to the unphysical interplay of parameters of the cyclotron line with the continuum models used previously, e.g., the negative and positive exponent power law (NPEX). For this model, we show that cyclotron line modeling erroneously leads to describing part of the exponential cutoff and the continuum variability, and not the cyclotron lines. When the X-ray continuum is modeled with a simple exponentially cutoff power law modified by a Gaussian emission feature around 10keV, the correlation between the line energy and the flux vanishes and the line parameters remain virtually constant over the outburst. We therefore conclude that the previously reported anti-correlation is an artifact of the assumptions adopted in the modeling of the continuum.

Accreting X-ray pulsars are among the best observed objects of X-ray astronomy with a rich data set of observational phenomena in the spectral and timing domain. While the general picture for these sources is well established, the detailed physics behind the observed phenomena are often subject of debate. We present recent observational, theoretical and modeling results for the structure and dynamics of the accretion column in these systems. Our results indicate the presence of different accretion regimes and possible explanations for observed variations of spectral features with luminosity.

Accretion-powered X-ray pulsars exhibit significant variability of the Cyclotron Resonance Scattering Feature (CRSF) centroid energy on pulse-to-pulse timescales, and also on much longer timescales. Two types of spectral variability are observed. For sources in group 1, the CRSF energy is negatively correlated with the variable source luminosity, and for sources in group 2, the opposite behavior is observed. The physical basis for this bimodal behavior is currently not understood. We explore the hypothesis that the accretion dynamics in the group 1 sources is dominated by radiation pressure near the stellar surface, and that Coulomb interactions decelerate the gas to rest in the group 2 sources. We derive a new expression for the critical luminosity such that radiation pressure decelerates the matter to rest in the supercritical sources. The formula for the critical luminosity is evaluated for 5 sources, using the maximum value of the CRSF centroid energy to estimate the surface magnetic field strength. The results confirm that the group 1 sources are supercritical and the group 2 sources are subcritical, although the situation is less clear for those highly variable sources that cross over the critical line. We also explain the variation of the CRSF energy with luminosity as a consequence of the variation of the characteristic emission height. The sign of the height variation is opposite in the supercritical and subcritical cases, hence creating the observed bimodal behavior.

Recently observed minute timescale variability of blazar emission at TeV energies has imposed severe constraints on jet models and TeV emission mechanisms. We focus on a robust jet instability to explain this variability. As a consequence of the bulk outflow of the jet plasma, the pressure is likely to be anisotropic, with the parallel pressure $P_{||}$ in the forward jet direction exceeding the perpendicular pressure $P_{\perp}$. Under these circumstances, the jet is susceptible to the firehose instability, which can cause disruptions in the large scale jet structure and result in variability of the observed radiation. For a realistic range of parameters, we find that the growth timescale of the firehose instability is $\approx$ a few minutes, in good agreement with the observed TeV variability timescales for Mrk 501 (Albert et al. 2007) and PKS 2155-304 (Aharonian et al. 2007).

High mass X-ray binaries are among the brightest objects of our Galaxy in the high energy domain (0.1-100 keV). Despite our relatively good knowledge of their basic emission mechanisms, the complex problem of understanding their time and energy dependent X-ray emission is not completely solved. <P> In this paper, we study the energy dependent pulse profiles of the high mass X-ray binary pulsar 4U 0115+63 to investigate how they are affected by cyclotron resonant scattering. <P> We analyze archival BeppoSAX and RXTE observations performed during the giant outburst of the source which occurred in 1999. We exploit a cross correlation technique to compare the pulse profiles in different energy ranges and develop a relativistic ray-tracing model to interpret our findings. We also study the phase dependency of the cyclotron absorption features by performing phase resolved spectroscopy. <P> The pulse profiles of 4U 0115+63 displayed clear "phase-lags" at energies close to those of the cyclotron absorption features that characterize the X-ray emission of the source. We reproduce qualitatively this phenomenon by assuming an energy dependent beaming of the emission from the column surface and verify that our model is also compatible with the results of phase resolved spectral analysis. <P> We showed that cyclotron resonant scattering affects the pulse profile formation mechanisms in a complex way, which necessitates both improvements in the modeling and the study of other sources to be better understood.

We develop and discuss global accretion solutions for viscous ADAF disks containing centrifugally supported isothermal shock waves. The fact that such shocks can exist at all in ADAF disks is a new result. Interestingly, we find that isothermal shocks can form even when the level of viscous dissipation is relatively high. In order to better understand this phenomenon, we explore all possible combinations of the fundamental flow parameters, such as specific energy, specific angular momentum, and viscosity, to obtain the complete family of global solutions. This procedure allows us to identify the region of the parameter space where isothermal shocks can exist in viscous ADAF disks. The allowed region is maximized in the inviscid case, and it shrinks as the level of viscous dissipation increases. Adopting the canonical value gamma=1.5 for the ratio of specific heats, we find that the shock region disappears completely when the Shakura-Sunyaev viscosity parameter alpha exceeds the critical value ~0.27. This establishes for the first time that steady ADAF disks containing shocks can exist even for relatively high levels of viscous dissipation. If an isothermal shock is present in the disk, it would have important implications for the acceleration of energetic particles that can escape to power the relativistic jets commonly observed around underfed, radio-loud black holes. In two specific applications, we confirm that the kinetic luminosity lost from the disk at the isothermal shock location is sufficient to power the observed relativistic outflows in M87 and Sgr A*.

In this Letter, we present a new self-consistent theory for the production of the relativistic outflows observed from radio-loud black hole candidates and active galaxies as a result of particle acceleration in hot, viscous accretion disks containing standing, centrifugally-supported isothermal shocks. This is the first work to obtain the structure of such disks for a relatively large value of the Shakura-Sunyaev viscosity parameter ($\alpha=0.1$), and to consider the implications of the shock for the acceleration of relativistic particles in viscous disks. In our approach, the hydrodynamics and the particle acceleration are coupled and the solutions are obtained self-consistently based on a rigorous mathematical method. We find that particle acceleration in the vicinity of the shock can provide enough energy to power the observed relativistic jet in M87.

The distribution function describing the acceleration of relativistic particles in an advection-dominated accretion disk is analyzed using a transport formalism that includes first-order Fermi acceleration, advection, spatial diffusion, and the escape of particles through the upper and lower surfaces of the disk. When a centrifugally-supported shock is present in the disk, the concentrated particle acceleration occurring in the vicinity of the shock channels a significant fraction of the binding energy of the accreting gas into a population of relativistic particles. These high-energy particles diffuse vertically through the disk and escape, carrying away both energy and entropy and allowing the remaining gas to accrete. The dynamical structure of the disk/shock system is computed self-consistently using a model previously developed by the authors that successfully accounts for the production of the observed relativistic outflows (jets) in M87 and \SgrA. This ensures that the rate at which energy is carried away from the disk by the escaping relativistic particles is equal to the drop in the radial energy flux at the shock location, as required for energy conservation. We investigate the influence of advection, diffusion, and acceleration on the particle distribution by computing the nonthermal Green's function, which displays a relatively flat power-law tail at high energies. We also obtain the energy distribution for the particles escaping from the disk, and we conclude by discussing the spectrum of the observable secondary radiation produced by the escaping particles.

Highly magnetized pulsars accreting matter in a binary system are bright sources in the X-ray band (0.1-100 keV). Despite the early comprehension of the basic emission mechanism, their spectral energy distribution is generally described by phenomenological or simplified models. We propose a study of the spectral emission from the high mass X-ray binary pulsar 4U 0115+634 by means of thermal and bulk Comptonization models based on the physical properties of such objects. For this purpose, we analyze the BeppoSAX data in the energy range 0.7-100 keV of the 1999 giant outburst, 12 days after the maximum. We model the spectral energy distribution of the system using a two-component continuum. At higher energy, above ~7 keV, the emission is due to thermal and bulk Comptonization of the seed photons produced by cyclotron cooling of the accretion column, and at lower energy, the emission is due to thermal Comptonization of a blackbody source in a diffuse halo close to the stellar surface. From the best fit parameters, we argue that the cyclotron emission is produced ~1.7 km above the stellar surface, and escapes from the column near its base, where the absorption features are generated by the interaction with the magnetic field in a surrounding halo. We find that in 4U 0115+634, the observed spectrum is dominated by reprocessed cyclotron radiation, whereas in other bright sources with stronger magnetic fields such as Her X-1, the spectrum is dominated by reprocessed bremsstrahlung.

We aim to illustrate the role of hot protons in enhancing the magnetorotational instability (MRI) via the ``hybrid'' viscosity, which is due to the redirection of protons interacting with static magnetic field perturbations, and to establish that it is the only relevant mechanism in this situation. It has recently been shown by Balbus \cite{PBM1} and Islam & Balbus \cite{PBM11} using a fluid approach that viscous momentum transport is key to the development of the MRI in accretion disks for a wide range of parameters. However, their results do not apply in hot, advection-dominated disks, which are collisionless. We develop a fluid picture using the hybrid viscosity mechanism, that applies in the collisionless limit. We demonstrate that viscous effects arising from this mechanism can significantly enhance the growth of the MRI as long as the plasma $\beta \gapprox 80$. Our results facilitate for the first time a direct comparison between the MHD and quasi-kinetic treatments of the magnetoviscous instability in hot, collisionless disks.

Accretion-powered X-ray pulsars are among the most luminous X-ray sources in the Galaxy. However, despite decades of theoretical and observational work since their discovery, no satisfactory model for the formation of the observed X-ray spectra has emerged. In this paper, we report on a self-consistent calculation of the spectrum emerging from a pulsar accretion column that includes an explicit treatment of the bulk and thermal Comptonization occurring in the radiation-dominated shocks that form in the accretion flows. Using a rigorous eigenfunction expansion method, we obtain a closed-form expression for the Green's function describing the upscattering of monochromatic radiation injected into the column. The Green's function is convolved with bremsstrahlung, cyclotron, and blackbody source terms to calculate the emergent photon spectrum. We show that energization of photons in the shock naturally produces an X-ray spectrum with a relatively flat continuum and a high-energy exponential cutoff. Finally, we demonstrate that our model yields good agreement with the spectra of the bright pulsar Her X-1 and the low luminosity pulsar X Per.

We develop a new theoretical model for the spectral formation process in accretion-powered X-ray pulsars based on a detailed treatment of the bulk and thermal Comptonization occurring in the accreting, shocked gas. A rigorous eigenfunction expansion method is employed to obtain the analytical solution for the Green's function describing the scattering of radiation injected into the column from a monochromatic source located at an arbitrary height above the stellar surface. The emergent spectrum is calculated by convolving the Green's function with source terms corresponding to bremsstrahlung, cyclotron, and blackbody emission. The energization of the photons in the shock, combined with cyclotron absorption, naturally produces an X-ray spectrum with a relatively flat continuum shape and a high-energy quasi-exponential cutoff. We demonstrate that the new theory successfully reproduces the phase-averaged spectra of the bright pulsars Her X-1, LMC X-4, and Cen X-3. In these luminous sources, it is shown that the emergent spectra are dominated by Comptonized bremsstrahlung emission.

Stochastic acceleration of charged particles due to interactions with magnetohydrodynamic (MHD) plasma waves is the dominant process leading to the formation of the high-energy electron and ion distributions in a variety of astrophysical systems. Collisions with the waves influence both the energization and the spatial transport of the particles, and therefore it is important to treat these two aspects of the problem in a self-consistent manner. We solve the representative Fokker-Planck equation to obtain a new, closed-form solution for the time-dependent Green's function describing the acceleration and escape of relativistic ions interacting with Alfven or fast-mode waves characterized by momentum diffusion coefficient $D(p)\propto p^q$ and mean particle escape timescale $t_esc(p) \propto p^{q-2}$, where $p$ is the particle momentum and $q$ is the power-law index of the MHD wave spectrum. In particular, we obtain solutions for the momentum distribution of the ions in the plasma and also for the momentum distribution of the escaping particles, which may form an energetic outflow. The general features of the solutions are illustrated via examples based on either a Kolmogorov or Kraichnan wave spectrum. The new expressions complement the results obtained by Park and Petrosian, who presented exact solutions for the hard-sphere scattering case ($q=2$) in addition to other scenarios in which the escape timescale has a power-law dependence on the momentum. Our results have direct relevance for models of high-energy radiation and cosmic-ray production in astrophysical environments such as $\gamma$-ray bursts, active galaxies, and magnetized coronae around black holes.

We reexamine the energetics of nonthermal electron acceleration in solar noise storms. A new result is obtained for the minimum nonthermal electron number density required to produce a Langmuir wave population of sufficient intensity to power the noise storm emission. We combine this constraint with the stochastic electron acceleration formalism developed by Subramanian & Becker (2005) to derive a rigorous estimate for the efficiency of the overall noise storm emission process, beginning with nonthermal electron acceleration and culminating in the observed radiation. We also calculate separate efficiencies for the electron acceleration -- Langmuir wave generation stage and the Langmuir wave -- noise storm production stage. In addition, we obtain a new theoretical estimate for the energy density of the Langmuir waves in noise storm continuum sources.

The nature of the viscosity operative in hot, two-temperature accretion disks around AGN has been a long-standing, unsolved problem. It has been previously suggested that protons, in conjunction with the turbulent magnetic field that is likely to exist in the accretion disk, might be crucial in providing this viscosity. Several authors have recently determined diffusion coefficients for charged particles (cosmic rays) propagating in turbulent magnetic fields by means of extensive Monte Carlo simulations. We use the diffusion coefficients for protons determined by these simulations to find the effective mean free path for protons in hot accretion disks. This in turn yields good estimates of the viscosity due to energetic protons embedded in the turbulent magnetic field of a hot, two-temperature accretion disk. We combine this with a simple two-temperature accretion disk model to determine the Shakura-Sunyaev $\alpha$ viscosity parameter arising out of this mechanism. We find that protons diffusing in the turbulent magentic field embedded in a hot accretion d isk provide a physically reasonable source of viscosity in hot accretion disks around AGN.

Relativistic outflows (jets) of matter are commonly observed from systems containing black holes. The strongest outflows occur in the radio-loud systems, in which the accretion disk is likely to have an advection-dominated structure. In these systems, it is clear that the binding energy of the accreting gas is emitted primarily in the form of particles rather than radiation. However, no comprehensive model for the disk structure and the associated outflows has yet been produced. In particular, none of the existing models establishes a direct physical connection between the presence of the outflows and the action of a microphysical acceleration mechanism operating in the disk. In this paper we explore the possibility that the relativistic protons powering the jet are accelerated at a standing, centrifugally-supported shock in the underlying accretion disk via the first-order Fermi mechanism. The theoretical analysis employed here parallels the early studies of cosmic-ray acceleration in supernova shock waves, and the particle acceleration and disk structure are treated in a coupled, self-consistent manner based on a rigorous mathematical approach. We find that first-order Fermi acceleration at standing shocks in advection-dominated disks proves to be a very efficient means for accelerating the jet particles. Using physical parameters appropriate for M87 and SgrA*, we verify that the jet kinetic luminosities computed using our model agree with estimates based on observations of the sources.

We have observed an unusually strong X-ray burst as a part of our regular eclipse timing observations of the low mass binary system EXO0748-676. The burst peak flux was 5.2x10^-8 ergs cm^-2 s^-1, approximately five times the normal peak X-ray burst flux observed from this source by RXTE. Spectral fits to the data strongly suggest that photospheric radius expansion occurred during the burst. In this Letter we examine the properties of this X-ray burst, which is the first example of a radius expansion burst from EXO0748-676 observed by RXTE. We find no evidence for coherent burst oscillations. Assuming that the peak burst luminosity is the Eddington luminosity for a 1.4 solar mass neutron star we derive a distance to EXO0748-676 of 7.7 kpc for a helium-dominated burst photosphere and 5.9 kpc for a hydrogen-dominated burst photosphere.

Accretion-powered X-ray pulsars are among the most luminous X-ray sources in the Galaxy. However, despite decades of theoretical and observational work since their discovery, no satisfactory model for the formation of the observed X-ray spectra has emerged. In particular, the previously available theories are unable to reproduce the power-law variation observed at high energies in many sources. In this paper, we present the first self-consistent calculation of the spectrum emerging from a pulsar accretion column that includes an explicit treatment of the energization occurring in the shock. Using a rigorous eigenfunction expansion method based on the exact dynamical solution for the velocity profile in the column, we obtain a closed-form expression for the Green's function describing the upscattering of radiation injected into the column from a monochromatic source located at the top of the thermal mound, near the base of the flow. The Green's function is convolved with a Planck distribution to calculate the radiation spectrum resulting from the reprocessing of blackbody photons emitted by the thermal mound. We demonstrate that the energization of the photons in the shock naturally produces an X-ray spectrum with a power-law shape at high energies and a blackbody shape at low energies, in agreement with many observations of accreting X-ray pulsars.

We present the first self-consistent model for the dynamics and the radiative transfer occurring in bright X-ray pulsar accretion columns, with a special focus on the role of the shock in energizing the emerging X-rays. The pressure inside the accretion column of a luminous X-ray pulsar is dominated by the photons, and consequently the equations describing the coupled radiative-dynamical structure must be solved simultaneously. Spectral formation in these sources is therefore a complex, nonlinear phenomenon. We obtain the analytical solution for the Green's function describing the upscattering of monochromatic radiation injected into the column from the thermal mound located near the base of the flow. The Green's function is convolved with a Planck distribution to model the X-ray spectrum resulting from the reprocessing of blackbody photons produced in the thermal mound. These photons diffuse through the infalling gas and eventually escape out the walls of the column, forming the observed X-ray spectrum. We show that the resulting column-integrated, phase-averaged spectrum has a power-law shape at high energies and a blackbody shape at low energies, in agreement with the observational data for many X-ray pulsars.

It has recently been demonstrated that the Shakura-Sunyaev prescription for the kinematic viscosity in an advection-dominated accretion disk yields physically reasonable solutions for the structure of the inflow close to the event horizon. In particular, no violations of relativistic causality occur at the horizon. This is somewhat surprising considering the diffusive nature of the angular momentum transport in the Shakura-Sunyaev scenario, and it is therefore natural to ask whether one can also obtain acceptable solutions for the disk structure based on the various alternative models for the viscosity that have been proposed, including the "deterministic" forms. In this paper we perform a rigorous asymptotic analysis of the structure of an advection-dominated accretion disk close to the event horizon of a nonrotating black hole based on three of the alternative prescriptions for the viscosity that have been suggested in the literature. We constrain the physical disk model by stipulating that the stress must vanish at the horizon, which is the fundamental inner boundary condition imposed by general relativity. Surprisingly, we find that none of the three alternative viscosity prescriptions yield physically acceptable disk structures close to the horizon when the zero-torque condition is applied, whether the flow is in vertical hydrostatic equilibrium or free-fall. Hence we conclude that the original Shakura-Sunyaev prescription is the only one proposed so far that is physically consistent close to the event horizon.

In this Letter, we suggest that the relativistic protons powering the outflows emanating from radio-loud systems containing black holes are accelerated at standing, centrifugally-supported shocks in hot, advection-dominated accretion disks. Such disks are ideal sites for first-order Fermi acceleration at shocks because the gas is tenuous, and consequently the mean free path for particle-particle collisions generally exceeds the thickness of the disk. The accelerated particles are therefore able to avoid thermalization, and as a result a small fraction of them achieve very high energies and escape from the disk. In our approach the hydrodynamics and the particle acceleration are coupled and the solutions are obtained self-consistently based on a rigorous mathematical treatment. The theoretical analysis of the particle transport parallels the early studies of cosmic-ray acceleration in supernova shock waves. We find that particle acceleration in the vicinity of the shock can extract enough energy to power a relativistic jet. Using physical parameters appropriate for M87 and Sgr A*, we confirm that the jet kinetic luminosities predicted by the theory agree with the observational estimates.

We reexamine arguments advanced by Hayashi & Matsuda (2001), who claim that several simple, physically motivated derivations based on mean free path theory for calculating the viscous torque in a quasi-Keplerian accretion disk yield results that are inconsistent with the generally accepted model. If correct, the ideas proposed by Hayashi & Matsuda would radically alter our understanding of the nature of the angular momentum transport in the disk, which is a central feature of accretion disk theory. However, in this paper we point out several fallacies in their arguments and show that there indeed exists a simple derivation based on mean free path theory that yields an expression for the viscous torque that is proportional to the radial derivative of the angular velocity in the accretion disk, as expected. The derivation is based on the analysis of the epicyclic motion of gas parcels in adjacent eddies in the disk.

We use a generic stochastic acceleration formalism to examine the power $L_{\rm in}$ (${\rm erg s^{-1}}$) input to nonthermal electrons that cause noise storm continuum emission. The analytical approach includes the derivation of the Green's function for a general second-order Fermi process, and its application to obtain the particular solution for the nonthermal electron distribution resulting from the acceleration of a Maxwellian source in the corona. We compare $L_{\rm in}$ with the power $L_{\rm out}$ observed in noise storm radiation. Using typical values for the various parameters, we find that $L_{\rm in} \sim 10^{23-26}$ ${\rm erg s^{-1}}$, yielding an efficiency estimate $\eta \equiv L_{\rm out}/L_{\rm in}$ in the range $10^{-10} \lsim \eta \lsim 10^{-6}$ for this nonthermal acceleration/radiation process. These results reflect the efficiency of the overall process, starting from electron acceleration and culminating in the observed noise storm emission.

We formulate and solve the diffusion problem of line photon propagation in a bulk outflow from a compact object (black hole or neutron star) using a generic assumption regarding the distribution of line photons within the outflow. Thomson scattering of the line photons within the expanding flow leads to a decrease of their energy which is of first order in v/c, where v is the outflow velocity and c is the speed of light. We demonstrate that the emergent line profile is closely related to the time distribution of photons diffusing through the flow (the light curve) and consists of a broad redshifted feature. We analyzed the line profiles for the general case of outflow density distribution. We emphasize that the redshifted lines are intrinsic properties of the powerful outflow that are supposed to be in many compact objects.

We obtain an exact, closed-form expression for the time-dependent Green's function solution to the Kompaneets equation. The result, which is expressed as the integral of a product of two Whittaker functions, describes the evolution in energy space of a photon distribution that is initially monoenergetic. Effects of spatial transport within a homogeneous scattering cloud are also included within the formalism. The Kompaneets equation that we solve includes both the recoil and energy diffusion terms, and therefore our solution for the Green's function approaches the Wien spectrum at large times. We show that the Green's function can be used to generate all of the previously known steady-state and time-dependent solutions to the Kompaneets equation. The new solution allows the direct determination of the spectrum, without the need to numerically solve the partial differential equation. Based upon the Green's function, we obtain a new time-dependent solution for the photon distribution resulting from the reprocessing of an optically thin bremsstrahlung initial spectrum with a low-energy cutoff. The new bremsstrahlung solution possesses a finite photon number density, and therefore it displays proper equilibration to a Wien spectrum at large times. The relevance of our results for the interpretation of emission from variable X-ray sources is discussed, with particular attention to the production of hard X-ray time lags, and the Compton broadening of narrow features such as iron lines.

The structure of the inner region of an advection-dominated accretion disk around a nonrotating black hole is explored by applying asymptotic analysis in the region just outside the event horizon. We assume that the viscous transport is described by the standard Shakura-Sunyaev prescription throughout the disk, including the inner region close to the horizon. One of our goals is to explore the self-consistency of this assumption by analyzing the causality of the viscous transport near the black hole. The effects of general relativity are incorporated in an approximate manner by utilizing a pseudo-Newtonian gravitational potential. Analysis of the conservation equations yields unique asymptotic forms for the behaviors of the radial inflow velocity, the density, the sound speed, and the angular velocity. The specific behaviors are determined by three quantities; namely, the accreted specific energy, the accreted specific angular momentum, and the accreted specific entropy. The additional requirement of passage through a sonic point further constrains the problem, leaving only two free parameters. Our detailed results confirm that the Shakura-Sunyaev viscosity yields a well-behaved flow structure in the inner region that satisfies the causality constraint. We also show that the velocity distribution predicted by our pseudo-Newtonian model agrees with general relativity in the vicinity of the horizon. The asymptotic expressions we derive provide convenient boundary conditions for the development of global models via numerical integration of the conservation equations.

Using a large set of optically detected hotspots in powerful extragalactic double radio sources, we examine the basic question of whether the detection of optical synchrotron radiation requires in-situ acceleration of relativistic electrons within the hotspots/lobes. For this, we take into account the jet's bulk relativistic motion, as well as its likely misalignment from the plane of the sky. Together, both these factors can drastically reduce the apparent range of the ultra-relativistic electrons ejected from the nucleus in the form of a jet. The conventionally adopted parameter space for the fundamental variables, namely, the hotspot magnetic field, radio source orientation angle relative to the line-of-sight and the bulk speed of the jet plasma, is considered. We find that the observed optical/near-IR synchrotron emission of the hotspots can be explained even if the radiating relativistic electrons were accelerated exclusively within the nuclear region, provided the energy losses incurred by the electrons during their transport down the jet are dominated by inverse Compton upscatterings of the cosmic microwave background photons. Under this circumstance, in-situ acceleration of relativistic electrons inside the hotspots or lobes is not found to be mandated by their reported optical/near-infrared detections.

The acceleration of relativistic particles due to repeated scattering across a shock wave remains the most attractive model for the production of energetic cosmic rays. This process has been analyzed extensively during the past two decades using the ``two-fluid'' model of diffusive shock acceleration. It is well known that 1, 2, or 3 distinct solutions for the flow structure can be found depending on the upstream parameters. The precise nature of the critical conditions delineating the number and character of shock transitions has remained unclear, mainly due to the inappropriate choice of parameters used in the determination of the upstream boundary conditions. We derive the exact critical conditions by reformulating the upstream boundary conditions in terms of two individual Mach numbers defined with respect to the cosmic-ray and gas sound speeds, respectively. The gas and cosmic-ray adiabatic indices are assumed to remain constant throughout the flow, although they may have arbitrary, independent values. Our results provide for the first time a complete, analytical classification of the parameter space of shock transitions in the two-fluid model. When multiple solutions are possible, we propose using the associated entropy distributions as a means for indentifying the most stable configuration.

The Newtonian ADAF model has been generalized recently to obtain the ADIOS model that includes outflows of energy and angular momentum. The utilization of a Newtonian gravitational potential limits the ability of the ADIOS model to describe the inner region of the disk, where any relativistic outflows are likely to originate. In this paper we modify the ADIOS scenario to incorporate a pseudo-Newtonian potential, which approximates the effects of general relativity. The analysis yields a unique, self-similar solution for the structure of the coupled disk/wind system. Interesting features of the new solution include the relativistic character of the outflow in the vicinity of the radius of marginal stability, which may help to explain the formation of relativistic jets in active galaxies. The unique density variation obtained agrees with that implied by the dependence of the X-ray hard time lags on the Fourier frequency for a number of accreting galactic black hole candidates.

We examine the possibility that the relativistic jets observed in many active galactic nuclei may be powered by the Fermi acceleration of protons in a tenuous corona above a two-temperature accretion disk. In this picture the acceleration arises as a consequence of the shearing motion of the magnetic field in the corona, which is anchored in the underlying Keplerian disk. The protons in the corona have a power-law distribution because the density there is too low for proton-proton collisions to thermalize the energy supplied via Fermi acceleration. The same shear acceleration mechanism also operates in the disk itself, however, there the density is high enough for thermalization to occur and consequently the disk protons have a Maxwellian distribution. Particle acceleration in the corona leads to the development of a pressure-driven wind that passes through a critical point and subsequently transforms into a relativistic jet at large distances from the black hole. We combine the critical conditions for the wind with the structure equations for the disk and the corona to obtain a coupled disk/corona/wind model. Using the coupled model we compute the asymptotic Lorentz factor $\Gamma_\infty$ of the jet as a function of the cylindrical starting radius at the base of the outflow, in the corona. Our results suggest that $\Gamma_\infty \lapprox 10$, which is consistent with observations of superluminal motion in blazars. We show that collisions between the jet and broad-line emission clouds can produce high-energy radiation with a luminosity sufficient to power the $\gamma$-rays observed from blazars. Subject headings: radiation mechanisms: non-thermal, accretion, accretion disks, acceleration of particles, gamma rays: theory

We develop a general expression for the gamma-gamma absorption coefficient for gamma-rays propagating in an arbitrary direction at an arbitrary point in space above an X-ray emitting accretion disk. The X-ray intensity is assumed to vary as a power law in energy and radius between the outer disk radius and the inner radius, which is the radius of marginal stability for a Schwarzschild black hole. We use our result for the absorption coefficient to calculate the gamma-gamma optical depth for gamma-rays created at an arbitrary height and propagating at an arbitrary angle relative to the disk axis. As an application, we use our formalism to compute the minimum distance between the central black hole and the site of production of the gamma-rays detected by EGRET during the June 1991 flare of 3C 279. Our results indicate that the ``focusing'' of the gamma-rays along the disk axis due to pair production is strong enough to explain the observed degree of alignment in blazar sources. If the gamma-rays are produced isotropically in gamma-ray blazars, then these objects should appear as bright MeV sources when viewed along off-axis lines of sight.

We derive expressions for a ``hybrid'' viscosity arising from Coulomb collisions between ions in the presence of tangled magnetic fields. Such magnetic fields are expected to be embedded in accretion flows onto black holes and other compact objects. Ion viscosity in the absence of magnetic fields has been shown to appreciable in hot accretion disks, but it was thought that the presence of the smallest of magnetic fields would effectively quench it. We show that the hybrid viscosity in hot (ion temperatures around 10^11-10^12 K) accretion disks results in values of the Shakura Sunyaev \alpha parameter around 0.01 to 0.1-0.5. We do not model the generation of the tangled magnetic fields; we merely assume its presence and assume it to be arbitrarily strong.