The fast variability observed in the X-ray emission from black-hole binaries has a very complex phenomenology, but offers the possibility to investigate directly the properties of the inner accretion flow. In particular, type-B oscillations in the 2-8 Hz range, observed in the Soft-Intermediate state, have been associated to the emission from a relativistic jet. We present the results of the timing and spectral analysis of a set of observations of the bright transient MAXI J1348-630 made with the NICER telescope. The observations are in the brightest part of the outburst and all feature a strong type-B QPO at ~4.5 Hz. We compute the energy dependence of the fractional rms and the phase lags at the QPO frequency, obtaining high signal-to-noise data and sampling for the first time at energies below 2 keV. The fractional rms decreases from more than 10% at 9 keV to 0.6% at 1.5 keV, and is constant below that energy. Taking the 2-3 keV band as reference, photons at all energies show a hard lag, increasing with the distance from the reference band. The behaviour below 2 keV has never been observed before, due to the higher energy bandpass of previous timing instruments. The energy spectrum can be fitted with a standard model for this state, consisting of a thin disc component and a harder power law, plus an emission line between 6 and 7 keV. We discuss the results, concentrating on the phase lags, and show that they can be interpreted within a Comptonization model.

Recently, we reported an observational correlation between a) the time-lag of the hard (9 - 15 keV) with respect to the soft (2 - 5 keV) X-ray photons in black-hole X-ray binaries (BHXRBs) and b) the power-law photon index $\Gamma$ of the X-ray spectrum. This was physically explained with a simple jet model, i.e., a model where the Comptonization (the Compton upscattering of soft photons) happens in the jet. Here, we report the inclination dependence of this correlation, which we also explain with our jet model. Photons that emerge at different polar angles from the jet axis have different spectra and different time-lags. Because of this, we can explain quantitatively the type-B QPOs of GX 339-4 as resulting from a precessing jet.

Galactic black-hole X-ray binaries (BHBs) emit a compact, optically thick, mildly relativistic radio jet when they are in the hard and hard-intermediate states. In these states, BHBs exhibit a correlation between the time lag of hard photons with respect to softer ones and the photon index of the power law component that characterizes the X-ray spectral continuum above $\sim$ 10 keV. The correlation, however, shows large scatter. Our objective is to investigate the role that the inclination of the system plays on the correlation between the time lag and the photon index. We find that the correlation between the time lag and the photon index is tight in low-inclination systems and becomes weaker in high-inclination systems. The amplitude of the lags is also larger at low and intermediate inclination angles than at high inclination. Our jet model that reproduces the process of Comptonization in an extended jet can account for the observations remarkably well.

We present the design and performance of RoboPol, a four-channel optical polarimeter operating at the Skinakas Observatory in Crete, Greece. RoboPol is capable of measuring both relative linear Stokes parameters $q$ and $u$ (and the total intensity $I$) in one sky exposure. Though primarily used to measure the polarization of point sources in the R-band, the instrument features additional filters (B, V and I), enabling multi-wavelength imaging polarimetry over a large field of view (13.6' $\times$ 13.6'). We demonstrate the accuracy and stability of the instrument throughout its five years of operation. Best performance is achieved within the central region of the field of view and in the R band. For such measurements the systematic uncertainty is below 0.1% in fractional linear polarization, $p$ (0.05% maximum likelihood). Throughout all observing seasons the instrumental polarization varies within 0.1% in $p$ and within 1$^\circ$ in polarization angle.

Black-hole transients exhibit a correlation between the time lag of hard photons with respect to softer ones and the photon index of the hard X-ray power law. The correlation is not very tight and therefore it is necessary to examine it source by source. The objective of the present work is to investigate in detail the time-lag -- photon-index correlation in GX 339-4. We have obtained RXTE energy spectra and light curves and have computed the photon index and the time lag of the $9 - 15$ keV photons with respect to the $2 - 6$ keV ones. The observations cover the first stages of the hard state, the pure hard state, and the hard-intermediate state. At low $\Gamma$, the correlation is positive and it becomes negative at large $\Gamma$. By assuming that the hard X-ray power law index $\Gamma$ is produced by inverse Compton scattering of soft disk photons in the jet, we have reproduced the entire correlation by varying two parameters in the jet: the radius of the jet at its base $R_0$ and the Thomson optical depth along the jet $\tau_\parallel$. We have found that, as the luminosity of the source increases, $R_0$ initially increases and then decreases. This behavior is expected in the context of the Cosmic Battery. As a further test of our model, we predict the break frequency in the radio spectrum as a function of the photon index during the rising part of an outburst.

We use results of our 3 year polarimetric monitoring program to investigate the previously suggested connection between rotations of the polarization plane in the optical emission of blazars and their gamma-ray flares in the GeV band. The homogeneous set of 40 rotation events in 24 sources detected by {\em RoboPol} is analysed together with the gamma-ray data provided by {\em Fermi}-LAT. We confirm that polarization plane rotations are indeed related to the closest gamma-ray flares in blazars and the time lags between these events are consistent with zero. Amplitudes of the rotations are anticorrelated with amplitudes of the gamma-ray flares. This is presumably caused by higher relativistic boosting (higher Doppler factors) in blazars that exhibit smaller amplitude polarization plane rotations. Moreover, the time scales of rotations and flares are marginally correlated.

We have performed a timing and spectral analysis of a set of black-hole binaries to study the correlation between the photon index and the time lag of the hard photons with respect to the soft ones. We provide further evidence that the timing and spectral properties in black-hole X-ray binaries are coupled. In particular, we find that the average time lag increases as the X-ray emission becomes softer. Although a correlation between the hardness of the X-ray spectrum and the time (or phase) lag has been reported in the past for a handful of systems, our study confirms that this correlated behaviour is a global property of black-hole X-ray binaries. We also demonstrate that the photon-index - time-lag correlation can be explained as a result of inverse Comptonization in a jet.

Black-hole and neutron-star X-ray binaries exhibit compact radio jets, when they are in the so called quiescent, hard, or hard intermediate states. The radio spectrum in these states is flat to slightly inverted. It is widely accepted that the energy distribution of the electrons, in the rest frame of the jet, is a power law with index p. A power-law energy distribution of the electrons in the jet is sufficient to explain the flat to slightly inverted spectrum emitted by the jet from radio to near infrared wavelengths, but is it necessary? Contrary to what our thinking was decades ago, now we know that the jets originate in the hot inner flow around black holes and neutron stars. Thus, we have investigated the spectrum that is emitted by a thermal jet with kT in the range 100-250 keV. We have computed the emitted spectrum from radio to near infrared using either a thermal distribution of electrons or a power-law one. We have found that parabolic jets with a thermal distribution of electrons give inverted spectra with alpha in the range 0-0.4, while jets with a power-law distribution of electrons give inverted spectra with alpha in the range 0-0.2. The rest of the parameters are kept the same in the two cases. The break frequency, which marks the transition from optically thick to optically thin synchrotron emission, is comparable for the two electron energy distributions. Our conclusion is that, contrary to common belief, it is not necessary to invoke a power-law energy distribution of the electrons in a jet to explain its flat to slightly inverted radio spectrum. A relativistic Maxwellian produces similar spectra. Thus, the widely invoked corona around black holes in X-ray binaries may actually be the jet.

We present average R-band optopolarimetric data, as well as variability parameters, from the first and second RoboPol observing season. We investigate whether gamma- ray--loud and gamma-ray--quiet blazars exhibit systematic differences in their optical polarization properties. We find that gamma-ray--loud blazars have a systematically higher polarization fraction (0.092) than gamma-ray--quiet blazars (0.031), with the hypothesis of the two samples being drawn from the same distribution of polarization fractions being rejected at the 3{\sigma} level. We have not found any evidence that this discrepancy is related to differences in the redshift distribution, rest-frame R-band lu- minosity density, or the source classification. The median polarization fraction versus synchrotron-peak-frequency plot shows an envelope implying that high synchrotron- peaked sources have a smaller range of median polarization fractions concentrated around lower values. Our gamma-ray--quiet sources show similar median polarization fractions although they are all low synchrotron-peaked. We also find that the random- ness of the polarization angle depends on the synchrotron peak frequency. For high synchrotron-peaked sources it tends to concentrate around preferred directions while for low synchrotron-peaked sources it is more variable and less likely to have a pre- ferred direction. We propose a scenario which mediates efficient particle acceleration in shocks and increases the helical B-field component immediately downstream of the shock.

We present a new set of optical polarization plane rotations in blazars, observed during the third year of operation of RoboPol. The entire set of rotation events discovered during three years of observations is analysed with the aim of determining whether these events are inherent in all blazars. It is found that the frequency of the polarization plane rotations varies widely among blazars. This variation cannot be explained either by a difference in the relativistic boosting or by selection effects caused by a difference in the average fractional polarization. We conclude that the rotations are characteristic of a subset of blazars and that they occur as a consequence of their intrinsic properties.

The stages before the formation of stars in molecular clouds are poorly understood. Insights can be gained by studying the properties of quiescent clouds, such as their magnetic field structure. The plane-of-the-sky orientation of the field can be traced by polarized starlight. We present the first extended, wide-field ($\sim$10 $\rm deg^2$) map of the Polaris Flare cloud in dust-absorption induced optical polarization of background stars, using the RoboPol polarimeter at the Skinakas Observatory. This is the first application of the wide-field imaging capabilities of RoboPol. The data were taken in the R-band and analysed with the automated reduction pipeline of the instrument. We present in detail optimizations in the reduction pipeline specific to wide-field observations. Our analysis resulted in reliable measurements of 641 stars with median fractional linear polarization 1.3%. The projected magnetic field shows a large scale ordered pattern. At high longitudes it appears to align with faint striations seen in the Herschel-SPIRE map of dust emission (250 $\mu m$), while in the central 4-5 deg$^2$ it shows an eddy-like feature. The overall polarization pattern we obtain is in good agreement with large scale measurements by Planck of the dust emission polarization in the same area of the sky.

Neutron star X-ray binaries emit a compact, optically thick, relativistic radio jet during low-luminosity, usually hard states, as Galactic black-hole X-ray binaries do. When radio emission is bright, a hard power-law tail without evidence for an exponential cutoff is observed in most systems. We have developed a jet model that explains many spectral and timing properties of black-hole binaries in the states where a jet is present. Our goal is to investigate whether our jet model can reproduce the hard tail, with the correct range of photon index and the absence of a high-energy cutoff, in neutron-star X-ray binaries. We have performed Monte Carlo simulations of the Compton upscattering of soft, accretion-disk or boundary layer photons, in the jet and computed the emergent energy spectra, as well as the time lag of hard photons with respect to softer ones as a function of Fourier frequency. We demonstrate that our jet model explains the observed power-law distribution with photon index in the range 1.8-3. With an appropriate choice of the parameters, the cutoff expected from Comptonization is shifted to energies above ~300 keV, producing a pure power law without any evidence for a rollover, in agreement with the observations. Our results reinforce the idea that the link between the outflow (jet) and inflow (disk) in X-ray binaries does not depend on the nature of the compact object, but on the process of accretion. Furthermore, we address the differences of jets in black-hole and neutron-star X-ray binaries and predict that the break frequency in the spectral energy distribution of neutron-star X-ray binaries, as a class, will be lower than that of black-hole binaries.

We present first results on polarization swings in optical emission of blazars obtained by RoboPol, a monitoring program of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events. A possible connection of polarization swing events with periods of high activity in gamma rays is investigated using the dataset obtained during the first season of operation. It was found that the brightest gamma-ray flares tend to be located closer in time to rotation events, which may be an indication of two separate mechanisms responsible for the rotations. Blazars with detected rotations have significantly larger amplitude and faster variations of polarization angle in optical than blazars without rotations. Our simulations show that the full set of observed rotations is not a likely outcome (probability $\le 1.5 \times 10^{-2}$) of a random walk of the polarization vector simulated by a multicell model. Furthermore, it is highly unlikely ($\sim 5 \times 10^{-5}$) that none of our rotations is physically connected with an increase in gamma-ray activity.

We present measurements of rotations of the optical polarization of blazars during the second year of operation of RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events, and we analyse the large set of rotation events discovered in two years of observation. We investigate patterns of variability in the polarization parameters and total flux density during the rotation events and compare them to the behaviour in a non-rotating state. We have searched for possible correlations between average parameters of the polarization-plane rotations and average parameters of polarization, with the following results: (1) there is no statistical association of the rotations with contemporaneous optical flares; (2) the average fractional polarization during the rotations tends to be lower than that in a non-rotating state; (3) the average fractional polarization during rotations is correlated with the rotation rate of the polarization plane in the jet rest frame; (4) it is likely that distributions of amplitudes and durations of the rotations have physical upper bounds, so arbitrarily long rotations are not realised in nature.

Galactic black-hole X-ray binaries emit a compact, optically thick, mildy relativistic radio jet when they are in the hard and hard-intermediate states. In a series of papers, we have developed a jet model and have shown, through Monte Carlo simulations, that our model can explain many observational results. In this work, we investigate one more constraining relationship between the cutoff energy and the phase lag during the early stages of an X-ray outburst of the black-hole X-ray binary GX 339-4: the cutoff energy decreases while the phase lag increases during the brightening of the hard state. We demonstrate that our jet model naturally explains the above correlation, with a minor modification consisting of introducing an acceleration zone at the base of the jet. The observed correlation between the cutoff energy and the phase lag suggests that the lags are produced by the hard component. Here we show that this correlation arises naturally if Comptonization in the jet produces these two quantities.

We present the results from the analysis of the broad-band X-ray spectra of 5 Anomalous X-ray Pulsars (AXPs) and Soft $\gamma$-ray Repeaters (SGRs). We fit their Suzaku and INTEGRAL spectra with models appropriate for the X-ray emission from the accretion flow onto a pulsar. We find that their X-ray spectra can be well described with this model. In particular we find that: (a) the radius of the accretion column is $\sim150-350$ m resulting in a transverse optical depth of $\sim 1$; (b) the vertical Thompson optical depth is $\approx 50-400$, and (c) their luminosity translates in accretion rates $\approx10^{15}\rm{g\, s^{-1}}$. These results are in good agreement with the predictions from the fall-back disk model, providing further support in the interpretation of AXPs and SGRs as accreting pulsars.

Aims: We summarize the current observational picture of the outbursts of black-hole X-ray transients (BHTs), based on the evolution traced in a hardness-luminosity diagram (HLD), and we offer a physical interpretation. Methods: The basic ingredient in our interpretation is the Poynting-Robertson Cosmic Battery (PRCB, Contopoulos & Kazanas 1998), which provides locally the poloidal magnetic field needed for the ejection of the jet. In addition, we make two assumptions, easily justifiable. The first is that the mass-accretion rate to the black hole in a BHT outburst has a generic bell-shaped form. This is guaranteed by the observational fact that all BHTs start their outburst and end it at the quiescent state. The second assumption is that at low accretion rates the accretion flow is geometrically thick, ADAF-like, while at high accretion rates it is geometrically thin. Results: Both, at the beginning and the end of an outburst, the PRCB establishes a strong poloidal magnetic field in the ADAF-like part of the accretion flow, and this explains naturally why a jet is always present in the right part of the HLD. In the left part of the HLD, the accretion flow is in the form of a thin disk, and such a disk cannot sustain a strong poloidal magnetic filed. Thus, no jet is expected in this part of the HLD. The counterclockwise traversal of the HLD is explained as follows: the poloidal magnetic field in the ADAF forces the flow to remain ADAF and the source to move upwards in the HLD rather than to turn left. Thus, the history of the system determines the counterclockwise traversal of the HLD. As a result, no BHT is expected to ever traverse the entire HLD curve in the clockwise direction. Conclusions: We offer a physical interpretation of accretion and ejection in BHTs with only one parameter, the mass transfer rate.

Polarization in classical Be stars results from Thomson scattering of the unpolarized light from the Be star in the circumstellar disc. Theory and observations agree that the maximum degree of polarization from isolated Be stars is < 4%. We report on the first optical polarimetric observations of the Be/X-ray binary EXO\,2030+375. We find that the optical (R band) light is strongly linearly polarized with a degee of polarization of 19%, the highest ever measured either in a classical or Be/X-ray binary. We argue that the interstellar medium cannot account for this high polarization degree and that a substantial amount must be intrinsic to the source. We propose that it may result from the alignment of non-spherical ferromagnetic grains in the Be star disc due to the strong neutron star magnetic field.

We report the polarized optical light curve of a gamma-ray burst afterglow obtained using the RoboPol instrument. Observations began 655 seconds after the initial burst of gamma-rays from GRB131030A, and continued uninterrupted for 2 hours. The afterglow displayed a low, constant fractional linear polarization of $p = (2.1 \pm 1.6)\,\%$ throughout, which is similar to the interstellar polarization measured on nearby stars. The optical brightness decay is consistent with a forward-shock propagating in a medium of constant density, and the low polarization fraction indicates a disordered magnetic field in the shock front. This supports the idea that the magnetic field is amplified by plasma instabilities on the shock front. These plasma instabilities produce strong magnetic fields with random directions on scales much smaller than the total observable region of the shock, and the resulting randomly-oriented polarization vectors sum to produce a low net polarization over the total observable region of the shock.

We present first results from RoboPol, a novel-design optical polarimeter operating at the Skinakas Observatory in Crete. The data, taken during the May - June 2013 commissioning of the instrument, constitute a single-epoch linear polarization survey of a sample of gamma-ray - loud blazars, defined according to unbiased and objective selection criteria, easily reproducible in simulations, as well as a comparison sample of, otherwise similar, gamma-ray - quiet blazars. As such, the results of this survey are appropriate for both phenomenological population studies and for tests of theoretical population models. We have measured polarization fractions as low as $0.015$ down to $R$ magnitude of 17 and as low as $0.035$ down to 18 magnitude. The hypothesis that the polarization fractions of gamma-ray - loud and gamma-ray - quiet blazars are drawn from the same distribution is rejected at the $10^{-3}$ level. We therefore conclude that gamma-ray - loud and gamma-ray - quiet sources have different optical polarization properties. This is the first time this statistical difference is demonstrated in optical wavelengths. The polarization fraction distributions of both samples are well-described by exponential distributions with averages of $\langle p \rangle =6.4 ^{+0.9}_{-0.8}\times 10^{-2}$ for gamma-ray--loud blazars, and $\langle p \rangle =3.2 ^{+2.0}_{-1.1}\times 10^{-2}$ for gamma-ray--quiet blazars. The most probable value for the difference of the means is $3.4^{+1.5}_{-2.0}\times 10^{-2}$. The distribution of polarization angles is statistically consistent with being uniform.

We describe the data reduction pipeline and control system for the RoboPol project. The RoboPol project is monitoring the optical $R$-band magnitude and linear polarization of a large sample of active galactic nuclei that is dominated by blazars. The pipeline calibrates and reduces each exposure frame, producing a measurement of the magnitude and linear polarization of every source in the $13'\times 13'$ field of view. The control system combines a dynamic scheduler, real-time data reduction, and telescope automation to allow high-efficiency unassisted observations.

In the fallback disk model for the persistent emission of Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs), the hard X-ray emission arises from bulk- and thermal Comptonization of bremsstrahlung photons, which are generated in the accretion column. The relatively low X-ray luminosity of these sources implies a moderate transverse optical depth to electron scattering, with photons executing a small number of shock crossings before escaping sideways. We explore the range of spectral shapes that can be obtained with this model and characterize the most important parameter dependencies. We use a Monte Carlo code to study the crisscrossing of photons in a radiative shock in an accretion column and compute the resulting spectrum. As expected, high-energy power-law X-ray spectra are produced in radiative shocks with photon-number spectral index larger than or about 0.5. We find that the required transverse optical depth is between 1 and 7. Such spectra are observed in low-luminosity X-ray pulsars. We demonstrate here with a simple model that Compton upscattering in the radiative shock in the accretion column can produce hard X-ray spectra similar to those seen in the persistent and transient emission of AXPs and SGRs. In particular, one can obtain a high-energy power-law spectrum, with photon-number spectral index ~ 1 and a cutoff at 100 - 200 keV, with a transverse Thomson optical depth of ~ 5, which is shown to be typical in AXPs/SGRs.

We investigate the X-ray enhancement and the long-term evolution of the recently discovered, second "low-B magnetar" Swift J1822.3-1606 in the frame of the fallback disk model. During a soft gamma burst episode, the inner disk matter is pushed back to larger radii forming a density gradient at the inner disk. Subsequent relaxation of the inner disk could account for the observed X-ray enhancement light curve of Swift J1822.3-1606. We obtain model fits to the X-ray data with basic disk parameters similar to those employed to explain the X-ray outburst light curves of other AXPs and SGRs. The long period (8.4 s) of the neutron star can be reached by the effect of the disk torques in the long-term accretion phase (1-3 x 10^5 yrs). The currently ongoing X-ray enhancement could be due to a transient accretion epoch, or the source could still be in the accretion phase in quiescence. Considering these different possibilities, we determine the model curves that could represent the long-term rotational and the X-ray luminosity evolution of Swift J1822.3-1606, which constrain the strength of the magnetic dipole field to the range of 1-2 x 10^12 G on the surface of the neutron star.

Recent measurements showed that the period derivative of the 'high-B' radio pulsar PSR J1734-3333 is increasing with time. For neutron stars evolving with fallback disks, this rotational behavior is expected in certain phases of the long-term evolution. Using the same model as employed earlier to explain the evolution of anomalous X-ray pulsars and soft gamma-ray repeaters, we show that the period, the first and second period derivatives and the X-ray luminosity of this source can simultaneously acquire the observed values for a neutron star evolving with a fallback disk. We find that the required strength of the dipole field that can produce the source properties is in the range of 10^{12} - 10^{13} G on the pole of the neutron star. When the model source reaches the current state properties of PSR J1734-3333, accretion onto the star has not started yet, allowing the source to operate as a regular radio pulsar. Our results imply that PSR J1734-3333 is at an age of ~ 3 x 10^4 - 2 x 10^5 years. Such sources will have properties like the X-ray dim isolated neutron stars or transient AXPs at a later epoch of weak accretion from the diminished fallback disk.

We propose that the quiescent emission of AXPs/SGRs is powered by accretion from a fallback disk, requiring magnetic dipole fields in the range 10^{12}-10^{13} G, and that the luminous hard tails of their X-ray spectra are produced by bulk-motion Comptonization in the radiative shock near the bottom of the accretion column. This radiation escapes as a fan beam, which is partly absorbed by the polar cap photosphere, heating it up to relatively high temperatures. The scattered component and the thermal emission from the polar cap form a polar beam. We test our model on the well-studied AXP 4U 0142+61, whose energy-dependent pulse profiles show double peaks, which we ascribe to the fan and polar beams. The temperature of the photosphere (kT~0.4 keV) is explained by the heating effect. The scattered part forms a hard component in the polar beam. We suggest that the observed high temperatures of the polar caps of AXPs/SGRs, compared with other young neutron stars, are due to the heating by the fan beam. Using beaming functions for the fan beam and the polar beam and taking gravitational bending into account, we fit the energy-dependent pulse profiles and obtain the inclination angle and the angle between the spin axis and the magnetic dipole axis, as well as the height of the radiative shock above the stellar surface. We do not explain the high luminosity bursts, which may be produced by the classical magnetar mechanism operating in super-strong multipole fields.

Neutron-star and black-hole X-ray binaries (XRBs) exhibit radio jets, whose properties depend on the X-ray spectral state and history of the source. In particular, black-hole XRBs emit compact, steady radio jets when they are in the so-called hard state, the jets become eruptive as the sources move toward the soft state, disappear in the soft state, and re-appear when the sources return to the hard state. On the other hand, jets from neutron-star X-ray binaries are typically weaker radio emitters than the black-hole ones at the same X-ray luminosity and in some cases radio emission is detected in the soft state. Significant phenomenology has been accumulated so far regarding the spectral states of neutron-star and black-hole XRBs, and there is general agreement about the type of the accretion disk around the compact object in the various spectral states. Our aim is to investigate whether the phenomenology regarding the X-ray emission on one hand and the jet appearance and disappearance on the other can be put together in a consistent physical picture. It has been shown that the so-called Poynting-Robertson Cosmic Battery (PRCB) explains in a natural way the formation of magnetic fields in the disks of AGN and the ejection of jets. We investigate whether the PRCB can also explain the formation, destruction, and variability of jets in XRBs. We find excellent agreement between the conditions under which the PRCB is efficient (i.e., the type of the accretion disk) and the emission or destruction of the radio jet. The disk-jet connection in XRBs is explained in a natural way using the PRCB.

We propose to perform a continuously scanning all-sky survey from 200 keV to 80 MeV achieving a sensitivity which is better by a factor of 40 or more compared to the previous missions in this energy range. The Gamma-Ray Imaging, Polarimetry and Spectroscopy (GRIPS) mission addresses fundamental questions in ESA's Cosmic Vision plan. Among the major themes of the strategic plan, GRIPS has its focus on the evolving, violent Universe, exploring a unique energy window. We propose to investigate $\gamma$-ray bursts and blazars, the mechanisms behind supernova explosions, nucleosynthesis and spallation, the enigmatic origin of positrons in our Galaxy, and the nature of radiation processes and particle acceleration in extreme cosmic sources including pulsars and magnetars. The natural energy scale for these non-thermal processes is of the order of MeV. Although they can be partially and indirectly studied using other methods, only the proposed GRIPS measurements will provide direct access to their primary photons. GRIPS will be a driver for the study of transient sources in the era of neutrino and gravitational wave observatories such as IceCUBE and LISA, establishing a new type of diagnostics in relativistic and nuclear astrophysics. This will support extrapolations to investigate star formation, galaxy evolution, and black hole formation at high redshifts.

We present a self-consistent model of the spectral energy distributions (SEDs) of spiral galaxies from the ultraviolet (UV) to the mid-infrared (MIR)/far-infrared (FIR)/submillimeter (submm) based on a full radiative transfer calculation of the propagation of starlight in galaxy disks. This model predicts not only the total integrated energy absorbed in the UV/optical and re-emitted in the infrared/submm, but also the colours of the dust emission based on an explicit calculation of the strength and colour of the UV/optical radiation fields heating the dust, and incorporating a full calculation of the stochastic heating of small dust grains and PAH molecules. The geometry of the translucent components of the model is empirically constrained using the results from the radiation transfer analysis of Xilouris et al. on spirals in the middle range of the Hubble sequence, while the geometry of the optically thick components is constrained from physical considerations with a posteriori checks of the model predictions with observational data. These geometrical constraints enable the dust emission to be predicted in terms of a minimum set of free parameters: the central face-on dust opacity in the B-band tau^f_B, a clumpiness factor F for the star-forming regions, the star-formation rate SFR, the normalised luminosity of the old stellar population old and the bulge-to-disk ratio B/D. We show that these parameters are almost orthogonal in their predicted effect on the colours of the dust/PAH emission. The results of the calculations are made available in the form of a large library of simulated dust emission SEDs spanning the whole parameter space of our model, together with the corresponding library of dust attenuation calculated using the same model. (see full abstract in the paper)

AXPs and SGRs constitute a special population of young neutron stars, which are thought to be magnetars, i.e., neutron stars with super-strong magnetic fields (10^14 - 10^15 G). Assuming that AXPs and SGRs accrete matter from a fallback disk, we attempt to explain the energy-dependent pulse profiles of AXP 4U 0142+61, as well as its phase-dependent energy spectra. In the fallback disk model, the Thomson optical depth along the accretion funnel is significant and bulk-motion Comptonization operates efficiently. This is enhanced by resonant cyclotron scattering. The thus scattered photons escape mainly sideways and produce a fan beam, which is detected as a main pulse up to energies of ~160 keV. The approximately isotropic emission from the stellar surface (soft thermal photons and reflected hard X-ray ones) is detected as a secondary pulse. This secondary pulse shows a bump at an energy of ~60 keV, which may be interpreted as resonant cyclotron scattering of fan-beam photons at the neutron-star surface, implying a dipole magnetic field strength B ~7 x 10^12 (1+z) G, where z is the gravitational redshift. Our model explains the soft and hard X-ray spectra of 4U 0142+61 and its energy dependent pulse profiles of the quiescent emission, while the short bursts are due to magnetar-type processes taking place in superstrong multiple fields.

Assuming that AXPs and SGRs accrete matter from a fallback disk, we attempt to explain both the soft and the hard X-ray emission as the result of the accretion process. We also attempt to explain their radio emission or the lack of it. We test the hypothesis that the power-law, hard X-ray spectra are produced in the accretion flow mainly by bulk-motion Comptonization of soft photons emitted at the neutron star surface. Fallback disk models invoke surface dipole magnetic fields of $10^{12} - 10^{13}$ G, which is what we assume here. Unlike normal X-ray pulsars, for which the accretion rate is highly super-Eddington, the accretion rate is approximately Eddington in AXPs and SGRs and thus the bulk-motion Comptonization operates efficiently. As an illustrative example we reproduce both the hard and the soft X-ray spectra of AXP 4U 0142+61 well using the XSPEC package compTB. Our model seems to explain both the hard and the soft X-ray spectra of AXPs and SGRs, as well as their radio emission or the lack of it, in a natural way. It might also explain the short bursts observed in these sources. On the other hand, it cannot explain the giant X-ray outbursts observed in SGRs, which may result from the conversion of magnetic energy in local multipole fields.

Some recent observational results impose significant constraints on all the models that have been proposed to explain the Galactic black-hole X-ray sources in the hard state. In particular, it has been found that during the hard state of Cyg X-1 the power-law photon number spectral index is correlated with the average time lag between hard and soft X-rays. Furthermore, the peak frequencies of the four Lorentzians that fit the observed power spectra are correlated with both the photon index and the time lag. We performed Monte Carlo simulations of Compton upscattering of soft, accretion-disk photons in the jet and computed the time lag between hard and soft photons and the power-law index of the resulting photon number spectra. We demonstrate that our jet model naturally explains the above correlations, with no additional requirements and no additional parameters.

A jet model for Galactic black-hole X-ray binaries will be presented that appears to explain several observational characteristics. In particular, it explains the energy spectrum from radio to hard X-rays, the time-lags as a function of Fourier frequency, the increase of the variability amplitude (QPO and high frequency) with increasing photon energy, and the narrowing of the autocorrelation function with increasing photon energy. On the other hand, there are additional observational constraints that no model has tried to explain yet. It is important that we all try to address these constraints if we are to make any progress in understanding black-hole X-ray sources

We present detailed modeling of the spectral energy distribution (SED) of the spiral galaxies NGC 891, NGC 4013, and NGC 5907 in the far-infrared and submm wavelengths. The model takes into account the emission produced by the diffuse dust and the star forming HII complexes. The dust mass is constrained by radiative transfer simulations in the optical (Xilouris et al. 1999). We find that the submm emission predicted by our model cannot account for the observed fluxes. Two scenarios may account for the "missing" submm flux. In the first scenario (Popescu et al. 2000), additional dust (to that derived from the optical, and associated with young stars) is embedded in the galaxy in the form of a thin disk and gives rise to additional submm emission. The other scenario investigates whether the average submm emissivity of the dust grains is higher than the values widely used in Galactic environments. In this case, the dust mass is equal to that derived from the optical observations, and the submm emissivity is treated as a free parameter calculated by fitting our model to the observed SED. We find the submm emissivity value to be ~3 times that often used for our Galaxy. While both scenarios reproduce the observed 850 micron surface brightness, the extra embedded dust model is not supported by the near infrared observations. We, thus, find that the enhanced dust submm emissivity scenario is the most plausible. [abridged]

The relationship between star formation and infrared emission in galaxies will be investigated. If galaxies were simple objects and young stars were completely covered with dust, then all the absorbed light of the young stars would be re-emitted in the infrared and from the infrared emission of galaxies we would infer the star formation rate (SFR) in them accurately. To show the complexities involved in real galaxies, we will use as a case study the late-type spiral galaxies. We will show that the heating of the dust is done mainly by the UV radiation of the young stars and therefore the infrared emission reveals the SFR in them. With a realistic model and its application to a number of galaxies, tight correlations are derived between SFR and total far infrared luminosity on one hand, and dust mass and 850 micron flux on the other. Other diagnostics of the SFR are examined and it is shown that there is consistency among them. Thus, the SFR for galaxies of all Hubble types has been determined as well as for interacting starburst galaxies. Combining different methods, the star-formation history of the universe has been determined and will be shown. Finally, some early results from the Spitzer Space Telescope will be presented.

The dust in galaxies makes radiation transport calculations in them absolutely necessary. It is not only common practice in Astrophysics, but also wisdom, to try to make as simple models as possible to simulate physical systems. For spiral galaxies, however, this turned out to be catastrophic. For years, the major question of the opacity of spiral galaxies kept the community divided, because the models were too simple. A spiral galaxy appears, to first order, to have exponential distributions of stars and dust, which cannot be approximated with uniform distributions. We will review the radiative transfer methods used in galaxies and we will comment on their pluses and minuses. We will also present some of the main results of the application of one of the methods to the observations.

We present detailed modelling of the spectral energy distribution (SED) of the spiral galaxies NGC 891, NGC 4013, and NGC 5907 in the far-infrared (FIR) and sub-millimeter (submm) wavelengths. The model takes into account the emission of the diffuse dust component, which is heated by the UV and optical radiation fields produced by the stars, as well as the emission produced locally in star forming HII complexes. The radiative transfer simulations of Xilouris et al. (1999) in the optical bands are used to constrain the stellar and dust geometrical parameters, as well as the total amount of dust. We find that the submm emission predicted by our model can not account for the observed fluxes at these wavelengths. We examine two cases, one having more dust embedded in a second thin disk and another allowing for an enhanced submillimeter emissivity of the dust grains. We argue that both cases can equally well reproduce the observed SED. The case of having more dust embedded in a second disk though, is not supported by the near-infrared observations and thus more realistic distributions of the dust (i.e., in spiral arms and clumps) have to be examined in order to better fit the surface brightness of each galaxy.

we propose a jet model for the low/hard state of black-hole X-ray sources which explains a) the X-ray spectra, b) the timelag spectra, c) the increase in the amplitude (QPO and high frequency) with increasing photon energy, and d) the narrowing of the autocorrelation function with increasing photon energy. The model (in its simplest form) assumes that i) there is a uniform magnetic field along the axis of the jet, ii) the electron density in the jet is inversely proportional to distance and iii) the jet is "hotter" near its center than at its periphery. We have performed Monte Carlo simulations of Compton upscattering of soft photons from the accretion disk and have found power-law high-energy spectra with photon number index in the range 1.5-2, power-law timelags versus Fourier frequency with index ~0.8, and an increase of the rms amplitude of the variability and a narrowing of the autocorrelation function with photon energy as they have been observed in Cygnus X-1.

We present new calculations of the attenuation of stellar light from spiral galaxies using geometries for stars and dust which can reproduce the entire spectral energy distribution from the UV to the FIR/submm and can also account for the surface brightness distribution in both the optical/NIR and FIR/submm. The calculations are based on the model of Popescu et al. (2000), which incorporates a dustless stellar bulge, a disk of old stars with associated diffuse dust, a thin disk of young stars with associated diffuse dust, and a clumpy dust component associated with star-forming regions in the thin disk. The attenuations, which incorporate the effects of multiple anisotropic scattering, are derived separately for each stellar component, and presented in the form of easily accessible polynomial fits as a function of inclination, for a grid in optical depth and wavelength. The wavelength range considered is between 912 AA and 2.2 micron, sampled such that attenuation can be conveniently calculated both for the standard optical bands and for the bands covered by GALEX. The attenuation characteristics of the individual stellar components show marked differences between each other. A general formula is given for the calculation of composite attenuation, valid for any combination of the bulge-to-disk ratio and amount of clumpiness. As an example, we show how the optical depth derived from the variation of attenuation with inclination depends on the bulge-to-disk ratio. Finally, a recipe is given for a self-consistent determination of the optical depth from the Halpha/Hbeta line ratio.

We address the problem of modeling the far-infrared (FIR) spectrum and deriving the star-formation rate (SFR) and the dust mass of spiral galaxies. We use the realistic physical model of Popescu et al. (2000) to describe the overall ultra-violet (UV), optical and FIR spectral energy distribution (SED) of a spiral galaxy. The model has only three free parameters: SFR, dust mass, and the fraction of the UV radiation which is absorbed locally by dense dust in the HII regions. Using this model, we are able to fit well the FIR spectra of 62 bright, IRAS galaxies from the "SCUBA Local Universe Galaxy Survey" of Dunne et al.(2000). As a result, we are able to determine, among others, their SFR and dust mass. We find that, on average, the SFR (in absolute units), the star-formation efficiency, the SFR surface density, and the ratio of FIR luminosity over the total intrinsic luminosity, are larger than the respective values of typical spiral galaxies of the same morphological type. We also find that the mean gas-to-dust mass ratio is close to the Galactic value, while the average central face-on optical depth of these galaxies in the V band is 2.3. Finally, we find a strong correlation between SFR or dust mass and observed FIR quantities like total FIR luminosity or FIR luminosity at 100 and 850 $\mathrm{\mu m}$. These correlations yield well defined relations, which can be used to determine a spiral galaxy's SFR and dust-mass content from FIR observations.

We present deep maps of dust emission from the edge-on spiral galaxy NGC 891, obtained with the ISOPHOT instrument on board the Infrared Space Observatory in broad band filters with reference wavelengths centered on 170 and 200 micron. Using new processing methods to remove the effects of detector transients from the data, we detect cold dust at high dynamic range. The observed surface brightness distribution and colour profile of the far-infrared (FIR) emission are found to be in good agreement with predictions for their counterparts derived from the model of Popescu et al. (2000a). Thus, NGC 891 is the first galaxy for which an intrinsic distribution of dust and stars could be found which simultaneously accounts for both the optical/near-IR and FIR morphologies.

Most, probably all, accreting binaries that are believed to contain a black-hole emit radio waves when they are in the low/hard state. Whenever this radio emission has been resolved, a jet-like structure has become apparent. We propose that Compton upscattering of low-energy photons in the jet can explain both the energy spectra and the time lags versus Fourier frequency observed in the low/hard state of black-hole systems. The soft photons originate in the inner part of the accretion disk. We have performed Monte Carlo simulations of Compton upscattering in a jet and have found that for a rather wide range of values of the parameters we can obtain power-law high-energy X-ray spectra with photon-number index in the range 1.5 - 2 and power-law time lags versus Fourier frequency with index ~ 0.7. The black-hole source Cyg X-1 in the low/hard state is well described by our model.

We have obtained discrete Fourier power spectra of a sample of persistent low-mass neutron-star X-ray binaries using long-term light curves from the All Sky Monitor on board the Rossi X-ray Timing Explorer. Our aim is to investigate their aperiodic variability at frequencies in the range 1 x 10^{-7}-5 x 10^{-6} Hz and compare their properties with those of the black-hole source Cyg X-1. We find that the classification scheme that divides LMXBs into Z and atoll sources blurs at very low frequencies. Based on the long-term (~ years) pattern of variability and the results of power-law fits (P ~ v^{-a}) to the 1 x 10^{-7}-5 x 10^{-6} Hz power density spectra, low-mass neutron-star binaries fall into three categories. Type I includes all Z sources, except Cyg X-2, and the atoll sources GX9+1 and GX13+1. They show relatively flat power spectra (a < 0.9) and low variability (rms < 20%). Type II systems comprise 4U 1636-53, 4U 1735-44 and GX3+1. They are more variable (20% < rms < 30%) and display steeper power spectra (0.9 < a < 1.2) than Type I sources. Type III systems are the most variable (rms > 30%) and exhibit the steepest power spectra (a > 1.2). The sources 4U 1705-44, GX354-0 and 4U 1820-30 belong to this group. GX9+9 and Cyg X-2 appear as intermediate systems in between Type I and II and Type II and III sources, respectively. We speculate that the differences in these systems may be caused by the presence of different types of mass-donor companions. Other factors, like the size of the accretion disc and/or the presence of weak magnetic fields, are also expected to affect their low-frequency X-ray aperiodic varibility.

We have carried out a timing analysis of Cyg X-1 and GRS 1915+105 using more than 5 years worth of data from the All Sky Monitor (ASM) on board RXTE. We have obtained for the first time power density spectra, colour-colour and colour-intensity plots, cross-correlation functions and phase-lag diagrams to investigate the variability of the sources at frequencies < 10^{-5} Hz. We find that the PDS are not flat but consistent with a power-law of index ~-1. This is the same slope as that found in the power spectra of Cyg X-1 during the soft state. In fact, the power spectrum of Cyg X-1 in the frequency range 10^{-7}-10^{-5} Hz appears as a continuation of the 10^{-3}-10 Hz power spectrum during the soft state. Significant variability, expressed as the fractional rms in the frequency range 2 x 10^{-7}-1 x 10^{-5} Hz, is detected at a level of 27% and 21% for GRS 1915+105 and Cyg X-1, respectively. Our results confirm previous suggestions that the innermost regions of the accretion disc are affected by variations occurring at larger radii from the black hole, presumably due to changes in the mass accretion rate. We also observe strong spectral changes: in Cyg X--1 the spectrum softens as the flux increases when it is in the low/hard state, and hardens with increasing flux in the high/soft state. GRS 1915+105 follows the same trend as Cyg X-1 in the soft state, but the hardness ratios show a larger amplitude of variation. The cross-correlation shows that the variations occur simultaneously in all energy bands. However, the cross-correlation functions are asymmetric toward negative lags. The phase spectrum also reveals negative lags at periods larger than a few days, suggesting that the long term variations in the soft energy band are delayed with respect to the harder bands.

We have performed Monte Carlo simulations of Compton upscattering of low-energy photons in an accretion disk around a Schwarzschild black hole. The photons gain energy from the rotational motion of the electrons in the disk. The upscattering occurs near the black hole horizon, where the flow velocity of the electrons approaches the speed of light. We show that this type of bulk-flow Comptonization can produce power-law X-ray spectra similar to the ones observed in black-hole X-ray transients in the high/soft state, i.e., a soft bump dominating the spectrum below ~ 10 keV and a power-law tail with photon index in the range 2-3. In order to reproduce the observed hard to soft flux ratio the disk has to have vertical optical depth above ~ 3 at the last stable orbit. We conclude that the power-law component of the high/soft state of black-hole transients may be due to an intrinsically cool disk extending all the way to the hole, without a separate hot plasma component.

Using tools previously described and applied to the prototype galaxy NGC 891, we model the optical to far-infrared spectral energy distributions (SED) of four additional edge-on spiral galaxies, namely NGC 5907, NGC 4013, UGC 1082 and UGC 2048. Comparing the model predictions with IRAS and, where available, sub-millimeter and millimeter observations, we determine the respective roles of the old and young stellar populations in grain heating. In all cases, the young population dominates, with the contribution of the old stellar population being at most 40%, as previously found for NGC 891. After normalization to the disk area, the massive star-formation rate (SFR) derived using our SED modeling technique, which is primarily sensitive to the non-ionizing ultraviolet output from the young stellar population, lies in the range 7e-4 - 2e-2 M_sun * yr^-1 * kpc^-2. This is consistent with normalized SFRs derived for face-on galaxies of comparable surface gas densities from H_alpha observations. Though the most active star-forming galaxy of the five in absolute terms, NGC 891 is not an exceptional system in terms of its surface density in SFR.

We describe a new tool for the analysis of the UV to the sub-millimeter (sub-mm) spectral energy distribution (SED) of spiral galaxies. We use a consistent treatment of grain heating and emission, solve the radiation transfer problem for a finite disk and bulge, and self-consistently calculate the stochastic heating of grains placed in the resulting radiation field. We use this tool to analyse the well-studied nearby edge-on spiral galaxy NGC891. First we investigate whether the old stellar population in NGC891, along with a reasonable assumption about the young stellar population, can account for the heating of the dust and the observed far-infrared (FIR) and sub-mm emission. The dust distribution is taken from the model of Xilouris et al. (1999), who used only optical and near-IR observations to determine it. We have found that such a simple model cannot reproduce the SED of NGC891, especially in the sub-mm range. It underestimates by a factor of 2-4 the observed sub-mm flux. A number of possible explanations exist for the missing sub-mm flux. We investigate a few of them and demonstrate that one can reproduce the observed SED in the FIR and the sub-mm quite well, as well as the observed radial profile at 850 microns. For the models calculated we give the relative proportion of the dust radiation powered by the old and young stellar populations as a function of FIR/sub-mm wavelength. In all models we find that the dust is predominantly heated by the young stellar population.

We have used Rossi X-ray Timing Explorer data to measure the lags between soft (2-5 keV) and hard (5-13 keV) photons and to study the aperiodic variability of the superluminal black hole candidate GRS 1915+105 during low-flux states. The power density spectra exhibit quasi-periodic oscillations (QPO) whose frequency increases with increasing count rate and varies in the frequency range 0.6-8 Hz. A correlation between the QPO frequency and the phase lag spectra is reported for the first time. This correlation is found for both the phase lag continuum and the phase lag at the QPO frequency. We find that as the QPO frequency moves to higher values the phase lags reverse sign from positive to negative. The absolute value of the lag always increases with photon energy. The negative (soft) lags are associated with a softer energy spectrum, whereas the positive (hard) lags are seen when the source is harder. We describe a possible scenario that may account for the change in the sign of the lags.

We compare 850um SCUBA images of NGC 891 with the corresponding V-band optical depth predicted from radiation transfer simulations. These two tracers of dust show a very similar distribution along the minor axis and a reasonable agreement along the major axis. Assuming that the grains responsible for optical extinction are also the source of 850um emission we derive a submillimeter emissivity (emission efficiency) for dust in the NGC 891 disk. This quantity is found to be a factor of 2-3 higher than the generally-accepted (but highly uncertain) values adopted for the Milky Way. It should be stated, however, that if a substantial fraction of dust in NGC 891 is clumped, the emissivity in the two galaxies may be quite similar. We use our newly-acquired emissivity to convert our 850um images into detailed maps of dust mass and, utilizing 21cm and CO-emission data for NGC 891, derive the gas-to-dust ratio along the disk. We compute an average ratio of 260 -- a value consistent with the Milky Way and external spirals within the uncertainties in deriving both the dust mass and the quantity of molecular gas. The bulk of dust in NGC 891 appears to be closely associated with the molecular gas phase although it may start to follow the distribution of atomic hydrogen at radii >9 kpc (i.e. >0.5 R_25). Using the optical depth of the NGC 891 disk, we quantify how light emitted at high redshift is attenuated by dust residing in foreground spirals. For B-band observations of galaxies typically found in the Hubble Deep Field, the amount of light lost is expected to be small (~ 5%). This value depends critically on the maximum radial extent of cold dust in spiral disks (which is poorly known). It may also represent a lower limit if galaxies expel dust over time into the intergalactic medium.

We demonstrate that the observed characteristics of the 5 cm OH masers in star-forming regions can be explained with the same model and the same parameters as the 18 cm and the 6 cm OH masers. In our already published study of the 18 cm and the 6 cm OH masers in star-forming regions we had examined the pumping of the 5 cm masers, but did not report the results we had found because of some missing collision rate coefficients, which in principle could be important. The recently published observations on the 5 cm masers of OH encourage us to report our old calculations along with some new ones that we have performed. These calculations, in agreement with the observations, reveal the main lines at 5 cm as strong masers, the 6049 MHz satellite line as a weak maser, and the 6017 MHz satellite line as never inverted for reasonable values of the parameters

We investigate how significant the spiral structure is on calculations concerning radiative transfer in dusty spiral galaxies seen edge-on. The widely adopted exponential disk model (i.e. both the stars and the dust are distributed exponentially in the radial direction and also perpendicular to the plane of the disk) is now subject to a detailed comparison with a realistic model that includes spiral structure for the stars and the dust in the disk. In particular, model images of galaxies with logarithmic spiral arms are constructed, such that the azimuthally averaged disk is exponential in radius and in height, as the observations suggest. Then, pure exponential disk models (i.e. with no spiral structure) are used to fit the edge-on appearance of the model images. As a result, the parameters derived after the fit are compared to the real values used to create the spiral-structured images. It turns out that the plain exponential disk model is able to give a very good description of the galactic disk with its parameters varying only by a few percent from their true values.

The opacity of spiral galaxies is examined by modelling the dust and stellar content of individual galaxies. The model is applied to five late-type spiral galaxies (NGC 4013, IC 2531, UGC 1082, NGC 5529 and NGC 5907). Having analyzed a total of seven galaxies thus far, the five galaxies mentioned above plus UGC 2048 and NGC 891 presented in (Xilouris et al. 1997, 1998), we are able to draw some general conclusions, the most significant of which are: 1) The face-on central optical depth is less than one in all optical bands indicating that typical spiral galaxies like the ones that we have modelled would be completely transparent if they were to be seen face-on. 2) The dust scaleheight is about half that of the stars, which means that the dust is more concentrated near the plane of the disk. 3) The dust scalelength is about 1.4 times larger than that of the stars and the dust is more radially extended than the stars. 4) The dust mass is found to be about an order of a magnitude more than previously measured using the IRAS fluxes, indicating the existence of a cold dust component. The gas-to-dust mass ratio calculated is close to the value derived for our Galaxy. 5) The derived extinction law matches quite well the Galactic extinction law, indicating a universal dust behaviour.

We study Compton upscattering of low-frequency photons in a converging flow of thermal plasma. The photons escape diffusively and electron scattering is the dominant source of opacity. We solve numerically and approximately analytically the equation of radiative transfer in the case of spherical, steady state accretion into black holes. Unlike previous work on this subject, we consider the inner boundary at a finite radius and this has a significant effect on the emergent spectrum. It is shown that the bulk motion of the converging flow is more efficient in upscattering photons than thermal Comptonization, provided that the electron temperature in the flow is of order a few keV or less. In this case, the spectrum observed at infinity consists of a soft component coming from those input photons which escaped after a few scatterings without any significant energy change and of a power law which extends to high energies and is made of those photons which underwent significant upscattering. The luminosity of the power law is relatively small compared to that of the soft component. The more reflective the inner boundary is, the flatter the power-law spectrum becomes. The spectral energy power-law index for black-hole accretion is always higher than 1 and it is approximately 1.5 for high accretion rates. This result tempts us to say that bulk motion Comptonization might be the mechanism behind the power-law spectra seen in black-hole X-ray sources.

Spectral formation in steady state, spherical accretion onto neutron stars and black holes is examined by solving numerically and analytically the equation of radiative transfer. The photons escape diffusively and their energy gains come from their scattering off thermal electrons in the converging flow of the accreting gas. We show that the bulk motion of the flow is more efficient in upscattering photons than thermal Comptonization in the range of non-relativistic electron temperatures. The spectrum observed at infinity is a power law with an exponential turnover at energies of order the electron rest mass. Especially in the case of accretion into a black hole, the spectral energy power-law index is distributed around 1.5. Because bulk motion near the horizon (1-5 Schwarzschild radii) is most likely a necessary characteristic of accretion into a black hole, we claim that observations of an extended power law up to about the electron rest mass, formed as a result of bulk motion Comptonization, is a real observational evidence for the existence of an underlying black hole.