We present a system of X-ray photometry for the Chandra satellite. X-ray photometry can be a powerful tool to obtain flux estimates, hardness ratios, and colors unbiased by assumptions about spectral shape and independent of temporal and spatial changes in instrument characteristics. The system we have developed relies on our knowledge of effective area and the energy-to-channel conversion to construct filters similar to photometric filters in the optical bandpass. We show that the filters are well behaved functions of energy and that this X-ray photometric system is able to reconstruct fluxes to within about 20%, without color corrections, for non-pathological spectra. Even in the worst cases it is better than 50%. Our method also treats errors in a consistent manner, both statistical as well as systematic.

We discuss the dynamics of microquasar jets in the interstellar medium, with specific focus on the effects of the X-ray binaries' space velocity with respect to the local Galactic standard of rest. We argue that, during late stages in the evolution of large scale radio nebulae around microquasars, the ram pressure of the interstellar medium due to the microquasar's space velocity becomes important and that microquasars with high velocities form the Galactic equivalent of extragalactic head-tail sources, i.e., that they leave behind trails of stripped radio plasma. Because of their higher space velocities, low-mass X-ray binaries are more likely to leave trails than high-mass X-ray binaries. We show that the volume of radio plasma released by microquasars over the history of the Galaxy is comparable to the disk volume and argue that a fraction of a few percent of the radio plasma left behind by the X-ray binary is likely mixed with the neutral phases of the ISM before the plasma is removed from the disk by buoyancy. Because the formation of microquasars is an unavoidable by-product of star formation, and because they can travel far from their birth places, their activity likely has important consequences for the evolution of magnetic fields in forming galaxies. We show that radio emission from the plasma inside the trail should be detectable at low frequencies. We suggest that LMXBs with high detected proper motions like XTE J1118+480 will be the best candidates for such a search.

In this paper we investigate the X-ray spectra and X-ray spectral variability of compact X-ray sources for 3 Chandra observations of the Local Group galaxy M33. The observations are centered on the nucleus and the star forming region NGC 604. In the observations 261 sources have been detected. For a total of 43 sources the number of net counts is above 100, sufficient for a more detailed spectral fitting. Of these sources, 25 have been observed in more than one observation, allowing the study of spectral variability on ~months timescales. A quarter of the sources are found to be variable between observations. However, except for two foreground sources, no source is variable within any observation above the 99% confidence level. Only six sources show significant spectral variability between observations. A comparison of N_H values with HI observations shows that X-ray absorption values are consistent with Galactic X-ray binaries and most sources in M33 are intrinsically absorbed. The pattern of variability and the spectral parameters of these sources are consistent with the M33 X-ray source population being dominated by X-ray binaries: Two thirds of the 43 bright sources have spectral and timing properties consistent with X-ray binaries; we also find two candidates for super-soft sources and two candidates for quasi-soft sources.

By combining the recently derived X-ray luminosity function for Galactic X-ray binaries (XRBs) by Grimm et al. (2002) and the radio-X-ray-mass relation of accreting black holes found by Merloni et al. (2003), we derive predictions for the radio luminosity function and radio flux distribution (logN/logS) for XRBs. Based on the interpretation that the radio-X-ray-mass relation is an expression of an underlying relation between jet power and nuclear radio luminosity, we derive the kinetic luminosity function for Galactic black hole jets, up to a normalization constant in jet power. We present estimates for this constant on the basis of known ratios of jet power to core flux for AGN jets and available limits for individual XRBs. We find that, if XRB jets do indeed fall on the same radio flux--kinetic power relation as AGN jets, the estimated mean kinetic luminosity of typical low/hard state jets is of the order of <W_{XRB}> ~ 2x10^37 ergs/s, with a total integrated power output of W \~ 5.5x10^38 ergs/s. We find that the power carried in transient jets should be of comparable magnitude to that carried in low/hard state jets. Including neutron star systems increases this estimate to W ~ 9x10^38 ergs/s. We estimate the total kinetic energy output from low/hard state jets over the history of the Galaxy to be E_{XRB} ~ 7x10^56 ergs.

In this paper we present the source list for 3 Chandra observations of the Local Group galaxy M33. The observations are centered on the nucleus and on the star forming region NGC 604. We detect a total of 261 sources in an area of \~0.2 square degree down to a flux limit of 3*10^{-16} erg/s/cm^2 which corresponds to a luminosity of ~2*10^34 erg/s at a distance of 840 kpc. From the source list we construct the luminosity functions of sources observed in M33. Taking into account background contamination the luminosity functions are consistent with those of other star forming galaxies. In addition, the combination of X-ray color analysis and the existence of ``blue'' optical counterparts strongly indicate that the X-ray point source population in M33 consists of young objects. Above 3*10^35 erg/s there are few X-ray sources in the locus of the X-ray hardness ratio diagram that is generally populated by LMXBs.

We study the statistical properties of the combined emission of a population of discrete sources (e.g. X-ray emission of a galaxy due to its X-ray binaries population). Namely, we consider the dependence of their total luminosity L_tot=SUM(L_k) and of fractional rms_tot of their variability on the number of sources N or, equivalently, on the normalization of the luminosity function. We show that due to small number statistics a regime exists, in which L_tot grows non-linearly with N, in an apparent contradiction with the seemingly obvious prediction <L_tot>=integral(dN/dL*L*dL) ~ N. In this non-linear regime, the rms_tot decreases with N significantly more slowly than expected from the rms ~ 1/sqrt(N) averaging law. For example, for a power law luminosity function with a slope of a=3/2, in the non-linear regime, L_tot ~ N^2 and the rms_tot does not depend at all on the number of sources N. Only in the limit of N>>1 do these quantities behave as intuitively expected, L_tot ~ N and rms_tot ~ 1/sqrt(N). We give exact solutions and derive convenient analytical approximations for L_tot and rms_tot. Using the total X-ray luminosity of a galaxy due to its X-ray binary population as an example, we show that the Lx-SFR and Lx-M* relations predicted from the respective ``universal'' luminosity functions of high and low mass X-ray binaries are in a good agreement with observations. Although caused by small number statistics the non-linear regime in these examples extends as far as SFR<4-5 Msun/yr and log(M*/Msun)<10.0-10.5, respectively.

We compare the results of Grimm et al. (2003) and Ranalli et al. (2003) on the Lx-SFR relation in normal galaxies. Based on the Lx-stellar mass dependence for LMXBs, we show, that low SFR (SFR<1 Msun/year) galaxies in the Ranalli et al. sample are contaminated by the X-ray emission from low mass X-ray binaries, unrelated to the current star formation activity. The most important conclusion from our comparison is, however, that after the data are corrected for the ``LMXB contamination'', the two datasets become consistent with each other, despite of their different content, variability effects, difference in the adopted source distances, X-ray flux and star formation rate determination and in the cosmological parameters used in interpreting the HDF-N data. They also agree well, both in the low and high SFR regimes, with the predicted Lx-SFR dependence derived from the parameters of the ``universal'' HMXB luminosity function. This encouraging result emphasizes the potential of the X-ray luminosity as an independent star formation rate indicator for normal galaxies.

Based on recent X-ray observations of the Milky Way, Magellanic Clouds and nearby starburst galaxies we study population of high mass X-ray binaries, their connection with ultra-luminous X-ray sources and relation to the star formation. Although more subtle SFR dependent effects are likely to exist, the data in the lg(Lx)~36-40.5 luminosity range are broadly consistent with existence of a ``universal'' luminosity function of HMXBs, which can be roughly described as a power law with differential slope of ~1.6 and a cutoff at lg(Lx)~40.5. The ULX sources found in many starburst galaxies occupy the high luminosity end of this single slope power law distribution, whereas its low luminosity part is composed of ``ordinary'' high mass X-ray binaries, observed, e.g. in the Milky Way and Magellanic Clouds. As the normalization of the ``universal'' luminosity function is proportional to the star formation rate, the number and/or the collective X-ray luminosity of HMXBs can be used to measure the current value of SFR in the host galaxy. Distant (unresolved) starburst galaxies observed by Chandra at redshifts of z~0.2-1.3 obey the same Lx-SFR relation as local galaxies, indicating that the ULXs at these redshifts were not significantly more luminous than those found in nearby galaxies.

Based on CHANDRA observations of nearby starburst galaxies and RXTE/ASM, ASCA and MIR-KVANT/TTM studies of high mass X-ray binary (HMXB) populations in the Milky Way and Magellanic Clouds, we propose that the number and/or the collective X-ray luminosity of HMXBs can be used to measure the star formation rate (SFR) of a galaxy. We show that, within the accuracy of the presently available data, a linear relation between HMXB number and the star formation rate exists. The relation between SFR and collective luminosity of HMXBs is non-linear in the low SFR regime, $L_X\propto \SFR^{\approx 1.7}$, and becomes linear only for sufficiently high star formation rate, when the total number of HMXB sources becomes sufficiently large. Such behaviour is caused by the fact, that we measure collective luminosity of a population of the discrete sources. Although more subtle SFR dependent effects are likely to exist, the data are broadly consistent with the existence of a universal luminosity function of HMXBs which can be roughly described as a power law with a differential slope of $\sim 1.6$, a cutoff at $L_X \sim few \times 10^{40}$ erg/sec and a normalisation proportional to the star formation rate. We apply our results to (spatially unresolved) starburst galaxies observed by CHANDRA at redshifts up to $z\sim 1$ in the Hubble Deep Field North and show that the calibration of the collective luminosity of HMXBs as a SFR indicator based on the local sample agrees well with the SFR estimates obtained for these distant galaxies with conventional methods.

We study the Log(N)-Log(S) and X-ray luminosity function in the 2-10 keV energy band, and the spatial (3-D) distribution of bright, log(L_X) > 34-35 erg/s, X-ray binaries in the Milky Way. In agreement with theoretical expectations and earlier results we found significant differences between the spatial distributions of low (LMXB) and high (HMXB) mass X-ray binaries. The volume density of LMXB sources peaks strongly at the Galactic Bulge. HMXBs tend to avoid the inner 3-4 kpc of the Galaxy, HMXBs are more concentrated towards the Galactic Plane and show clear signatures of the spiral structure in their spatial distribution. LMXB sources have a flatter Log(N)-Log(S) distribution and luminosity function than HMXBs. The integrated 2-10 keV luminosities of X-ray binaries, averaged over 1996--2000, are 2-3 * 10^39 (LMXB) and 2-3 * 10^38 (HMXB) erg/s. Normalised to the stellar mass and the star formation rate, respectively, these correspond to 5 * 10^28 erg/s/M_sol for LMXBs and 5 * 10^37 erg/s/(M_sol/yr) for HMXBs. Due to the shallow slopes of the luminosity functions the integrated emission of X-ray binaries is dominated by the 5-10 most luminous sources which determine the appearance of the Milky Way in the standard X-ray band for an outside observer. In particular variability of individual sources or an outburst of a bright transient source can increase the integrated luminosity of the Milky Way by as much as a factor of ~2. Although the average LMXB luminosity function shows a break near the Eddington luminosity for a 1.4 M_sol neutron star, at least 11 sources showed episodes of super-Eddington luminosity during ASM observations. We provide the maps of distribution of X-ray binaries in the Milky Way in various projections, which can be compared to images of nearby galaxies taken by CHANDRA and XMM-Newton.