This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

We measure cosmic weak lensing shear power spectra with the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalog covering 137deg$^2$ of the sky. Thanks to the high effective galaxy number density of $\sim$17 arcmin$^{-2}$ even after conservative cuts such as magnitude cut of $i<24.5$ and photometric redshift cut of $0.3\leq z \leq 1.5$, we obtain a high significance measurement of the cosmic shear power spectra in 4 tomographic redshift bins, achieving a total signal-to-noise ratio of 16 in the multipole range $300 \leq \ell \leq 1900$. We carefully account for various uncertainties in our analysis including the intrinsic alignment of galaxies, scatters and biases in photometric redshifts, residual uncertainties in the shear measurement, and modeling of the matter power spectrum. The accuracy of our power spectrum measurement method as well as our analytic model of the covariance matrix are tested against realistic mock shear catalogs. For a flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, we find $S_8\equiv \sigma_8(\Omega_{\rm m}/0.3)^\alpha=0.800^{+0.029}_{-0.028}$ for $\alpha=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $\alpha=0.5$) from our HSC tomographic cosmic shear analysis alone. In comparison with Planck cosmic microwave background constraints, our results prefer slightly lower values of $S_8$, although metrics such as the Bayesian evidence ratio test do not show significant evidence for discordance between these results. We study the effect of possible additional systematic errors that are unaccounted in our fiducial cosmic shear analysis, and find that they can shift the best-fit values of $S_8$ by up to $\sim 0.6\sigma$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

Photometric redshifts are a key component of many science objectives in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). In this paper, we describe and compare the codes used to compute photometric redshifts for HSC-SSP, how we calibrate them, and the typical accuracy we achieve with the HSC five-band photometry (grizy). We introduce a new point estimator based on an improved loss function and demonstrate that it works better than other commonly used estimators. We find that our photo-z's are most accurate at 0.2<~zphot<~1.5, where we can straddle the 4000A break. We achieve sigma(d_zphot/(1+zphot))~0.05 and an outlier rate of about 15% for galaxies down to i=25 within this redshift range. If we limit to a brighter sample of i<24, we achieve sigma~0.04 and ~8% outliers. Our photo-z's should thus enable many science cases for HSC-SSP. We also characterize the accuracy of our redshift probability distribution function (PDF) and discover that some codes over/under-estimate the redshift uncertainties, which have implications for N(z) reconstruction. Our photo-z products for the entire area in the Public Data Release 1 are publicly available, and both our catalog products (such as point estimates) and full PDFs can be retrieved from the data release site, https://hsc-release.mtk.nao.ac.jp/.

We present the deepest optical images of the COSMOS field based on a joint dataset taken with Hyper Suprime-Cam (HSC) by the HSC Subaru Strategic Program (SSP) team and the University of Hawaii (UH). The COSMOS field is one of the key extragalactic fields with a wealth of deep, multi-wavelength data. However, the current optical data are not sufficiently deep to match with, e.g., the UltraVista data in the near-infrared. The SSP team and UH have joined forces to produce very deep optical images of the COSMOS field by combining data from both teams. The coadd images reach depths of g=27.8, r=27.7, i=27.6, z=26.8, and y=26.2 mag at 5 sigma for point sources based on flux uncertainties quoted by the pipeline and they cover essentially the entire COSMOS 2 square degree field. The seeing is between 0.6 and 0.9 arcsec on the coadds. We perform several quality checks and confirm that the data are of science quality; ~2% photometry and 30 mas astrometry. This accuracy is identical to the Public Data Release 1 from HSC-SSP. We make the joint dataset including fully calibrated catalogs of detected objects available to the community at https://hsc-release.mtk.nao.ac.jp/.

In this paper, we describe the optical imaging data processing pipeline developed for the Subaru Telescope's Hyper Suprime-Cam (HSC) instrument. The HSC Pipeline builds on the prototype pipeline being developed by the Large Synoptic Survey Telescope's Data Management system, adding customizations for HSC, large-scale processing capabilities, and novel algorithms that have since been reincorporated into the LSST codebase. While designed primarily to reduce HSC Subaru Strategic Program (SSP) data, it is also the recommended pipeline for reducing general-observer HSC data. The HSC pipeline includes high level processing steps that generate coadded images and science-ready catalogs as well as low-level detrending and image characterizations.

We describe a weak lensing view of the downsizing of star forming galaxies based on cross correlating a weak lensing ($\kappa$) map with a predicted map constructed from a redshift survey. Moderately deep and high resolution images with Subaru/Hyper Suprime-Cam covering the 4 deg^2 DLS F2 field provide a $\kappa$ map with 1 arcmin resolution. A dense complete redshift survey of the F2 field including 12,705 galaxies with $R\leq20.6$ is the basis for construction of the predicted map. The zero-lag cross-correlation between the \kappa and predicted maps is significant at the $30\sigma$ level. The width of the cross-correlation peak is comparable with the angular scale of rich cluster at $z\sim0.3$, the median depth of the redshift survey. Slices of the predicted map in $\delta{z} = 0.05$ redshift bins enable exploration of the impact of structure as a function of redshift. The zero-lag normalised cross-correlation has significant local maxima at redshifts coinciding with known massive X-ray clusters. Even in slices where there are no known massive clusters, there is significant signal in the cross-correlation originating from lower mass groups that trace the large-scale of the universe. Spectroscopic $D_n4000$ measurements enable division of the sample into star-forming and quiescent populations. The significance of the cross-correlation with structure containing star-forming galaxies increases with redshift from $5\sigma$ at $z = 0.3$ to $7 \sigma$ at $z = 0.5$. The weak lensing results are consistent with the downsizing view of galaxy evolution established on the basis of many other independent studies.

We present measurements of the evolution of normal-galaxy X-ray emission from $z \approx$ 0-7 using local galaxies and galaxy samples in the 6 Ms Chandra Deep Field-South (CDF-S) survey. The majority of the CDF-S galaxies are observed at rest-frame energies above 2 keV, where the emission is expected to be dominated by X-ray binary (XRB) populations; however, hot gas is expected to provide small contributions to the observed- frame < 1 keV emission at $z < 1$. We show that a single scaling relation between X-ray luminosity ($L_{\rm X}$) and star-formation rate (SFR) is insufficient for characterizing the average X-ray emission at all redshifts. We establish that scaling relations involving not only SFR, but also stellar mass ($M_\star$) and redshift, provide significantly improved characterizations of the average X-ray emission from normal galaxy populations at $z \approx$ 0-7. We further provide the first empirical constraints on the redshift evolution of X-ray emission from both low-mass XRB (LMXB) and high-mass XRB (HMXB) populations and their scalings with $M_\star$ and SFR, respectively. We find $L_{\rm 2-10 keV}({\rm LMXB})/M_\star \propto (1 + z)^{2-3}$ and $L_{\rm 2-10 keV}({\rm HMXB})$/SFR $\propto (1 + z)$, and show that these relations are consistent with XRB population-synthesis model predictions, which attribute the increase in LMXB and HMXB scaling relations with redshift as being due to declining host galaxy stellar ages and metallicities, respectively. We discuss how emission from XRBs could provide an important source of heating to the intergalactic medium in the early Universe, exceeding that of active galactic nuclei.

We present a comprehensive study of the total X-ray emission from the colliding galaxy pair NGC2207/IC2163, based on Chandra, Spitzer, and GALEX data. We detect 28 ultra-luminous X-ray sources (ULXs), 7 of which were not detected previously due to X-ray variability. Twelve sources show significant long-term variability, with no correlated spectral changes. Seven sources are transient candidates. One ULX coincides with an extremely blue star cluster (B-V = -0.7). We confirm that the global relation between the number and luminosity of ULXs and the integrated star formation rate (SFR) of the host galaxy also holds on local scales. We investigate the effects of dust extinction and/or age on the X-ray binary (XRB) population on sub-galactic scales. The distributions of Nx and Lx are peaked at L(IR)/L(NUV)~1, which may be associated with an age of ~10 Myr for the underlying stellar population. We find that ~1/3 of the XRBs are located in close proximity to young star complexes. The luminosity function of the X-ray binaries is consistent with that typical for high-mass X-ray binaries, and appears unaffected by variability. We disentangle and compare the X-ray diffuse spectrum with that of the bright XRBs. The hot interstellar medium dominates the diffuse X-ray emission at E<1 keV, has a temperature kT=0.28 (+0.05/-0.04) keV and intrinsic 0.5-2 keV luminosity of 7.9e+40 erg/s, a factor of ~2.3 higher than the average thermal luminosity produced per unit SFR in local star-forming galaxies. The total X-ray output of NGC2207/IC2163 is 1.5e+41 erg/s, and the corresponding total integrated SFR is 23.7 Msol/yr.

We present results from a Chandra study of ultraluminous X-ray sources (ULXs) in a sample of 17 nearby (D_L<60 Mpc) luminous infrared galaxies (LIRGs), selected to have star formation rates (SFRs) in excess of 7 M_sun yr^-1 and low foreground Galactic column densities (N_H < 5*10^20 cm^-2). A total of 53 ULXs were detected and we confirm that this is a complete catalogue of ULXs for the LIRG sample. We examine the evolution of ULX spectra with luminosity by stacking the spectra of individual objects in three luminosity bins, finding a distinct change in spectral index at luminosity ~2 *10^39 erg s^-1. This may be a change in spectrum as 10 M_sun black holes transit from a ~Eddington to a super-Eddington accretion regime, and is supported by a plausible detection of partially-ionised absorption imprinted on the spectrum of the luminous ULX (L_X ~5*10^39 erg s^-1) CXOU J024238.9-000055 in NGC 1068, consistent with the highly ionised massive wind that we would expect to see driven by a super-Eddington accretion flow. This sample shows a large deficit in the number of ULXs detected per unit SFR (0.2 ULXs M_sun^-1 yr^-1). This deficit also manifests itself as a lower differential X-ray luminosity function normalisation for the LIRG sample than for samples of other star forming galaxies. We show that it is unlikely that this deficit is a purely observational effect. Part of this deficit might be attributable to the high metallicity of the LIRGs impeding the production efficiency of ULXs and/or a lag between the star formation starting and the production of ULXs; however, we argue that the evidence -- including very low N_ULX/L_FIR, and an even lower ULX incidence in the central regions of the LIRGs -- shows that the main culprit for this deficit is likely to be the high column of gas and dust in these galaxies, that fuels the high SFR but also acts to obscure many ULXs from our view.

The cosmological 21cm signal is a physics-rich probe of the early Universe, encoding information about both the ionization and the thermal history of the intergalactic medium (IGM). The latter is likely governed by X-rays from star-formation processes inside very high redshift (z > 15) galaxies. Due to the strong dependence of the mean free path on the photon energy, the X-ray SED can have a significant impact on the interferometric signal from the cosmic dawn. Recent Chandra observations of nearby, star-forming galaxies show that their SEDs are more complicated than is usually assumed in 21cm studies. In particular, these galaxies have ubiquitous, sub-keV thermal emission from the hot interstellar medium (ISM), which generally dominates the soft X-ray luminosity (with energies < 1 keV, sufficiently low to significantly interact with the IGM). Using illustrative soft and hard SEDs, we show that the IGM temperature fluctuations in the early Universe would be substantially increased if the X-ray spectra of the first galaxies were dominated by the hot ISM, compared with X-ray binaries with harder spectra. The associated large-scale power of the 21cm signal would be higher by roughly a factor of three. More generally, we show that the peak in the redshift evolution of the large-scale (k = 0.2 1/Mpc) 21cm power is a robust probe of the soft-band SED of the first galaxies, and importantly, is not degenerate with their bolometric luminosities. On the other hand, the redshift of the peak constrains the X-ray luminosity and halo masses which host the first galaxies.

We report significant anisotropies in the projected two-dimensional (2D) spatial distributions of Globular Clusters (GCs) of the giant Virgo elliptical galaxy NGC4649 (M60). Similar features are found in the 2D distribution of low-mass X-ray binaries (LMXBs), both associated with GCs and in the stellar field. Deviations from azimuthal symmetry suggest an arc-like excess of GCs extending north at 4-15 kpc galactocentric radii in the eastern side of major axis of NGC4649. This feature is more prominent for red GCs, but still persists in the 2D distribution of blue GCs. High and low luminosity GCs also show some segregation along this arc, with high-luminosity GCs preferentially located in the southern end and low-luminosity GCs in the northern section of the arc. GC-LMXBs follow the anisotropy of red-GCs, where most of them reside; however, a significant overdensity of (high-luminosity) field LMXBs is present to the south of the GC arc. These results suggest that NGC4649 has experienced mergers and/or multiple accretions of less massive satellite galaxies during its evolution, of which the GCs in the arc may be the fossil remnant. We speculate that the observed anisotropy in the field LMXB spatial distribution indicates that these X-ray binaries may be the remnants of a star formation event connected with the merger, or maybe be ejected from the parent red GCs, if the bulk motion of these clusters is significantly affected by dynamical friction. We also detect a luminosity enhancement in the X-ray source population of the companion spiral galaxy NGC4647. We suggest that these may be younger high mass X-ray binaries formed as a result of the tidal interaction of this galaxy with NGC4649.

We investigate the radial distribution of the low-mass X-ray binary (LMXB) population in the elliptical galaxy NGC 4649, using Chandra and Hubble data to separate the field and globular cluster (GC) populations. GCs with LMXBs have the same radial distribution as the parent red and blue GCs. The radial profile of field LMXBs follows the V-band profile within the D25 of NGC 4649. Using the spatial information provided by our data, we find that the global galaxy-wide relations between cumulative number and luminosity of LMXBs and the integrated stellar mass hold on local scales within D25. An excess of field LMXBs with respect to the V-light is observed in the galaxy's outskirt, which may be partially due to unidentified GC sources or to a rejuvenated field LMXB population caused by past merging interactions.

We investigate the relation between total X-ray emission from star-forming galaxies and their star formation activity. Using nearby late-type galaxies and ULIRGs from Paper I and star-forming galaxies from Chandra Deep Fields, we construct a sample of 66 galaxies spanning the redshift range z~0-1.3 and the star-formation rate (SFR) range ~0.1-10^3 M_sun/yr. In agreement with previous results, we find that the Lx-SFR relation is consistent with a linear law both at z=0 and for the z=0.1-1.3 CDF galaxies, within the statistical accuracy of ~0.1 in the slope of the Lx-SFR relation. For the total sample, we find a linear scaling relation Lx/SFR~(4.0\pm 0.4)x10^{39}(erg/s)/(Msun/yr), with a scatter of ~0.4 dex. About ~2/3 of the 0.5-8 keV luminosity generated per unit SFR is expected to be due to HMXBs. We find no statistically significant trends in the mean Lx/SFR ratio with the redshift or star formation rate and constrain the amplitude of its variations by <0.1-0.2 dex. These properties make X-ray observations a powerful tool to measure the star formation rate in normal star-forming galaxies that dominate the source counts at faint fluxes.

We present a new method for the determination of the two-dimensional (2D) projected spatial distribution of globular clusters (GCs) in external galaxies. This method is based on the K-Nearest Neighbor density estimator of Dressler (1980), complemented by MonteCarlo simulations to establish the statistical significance of the results. We apply this method to NGC4261, a "test galaxy" where significant 2D anisotropy in the GC distribution has been reported. We confirm that the 2D distribution of GC is not azimuthally isotropic. Moreover, we demonstrate that the 2D distribution departures from the average GC radial distribution results in highly significant spiral-like or broken shell features. Overall, the same perturbations are found in "red" and "blue" GCs, but with some differences. In particular, we observe a central feature, roughly aligned with the minor axis of NGC4261, composed of red and most luminous GCs. Blue and fainter GCs are more frequent at large radial distances and follow the spiral-like features of the overall density structure. These results suggest a complex merging history for NGC4261.

The colliding galaxy pair NGC 2207/IC 2163, at a distance of ~39 Mpc, was observed with Chandra, and an analysis reveals 28 well resolved X-ray sources, including 21 ultraluminous X-ray sources (ULXs) with Lx > 10^39 erg/s, as well as the nucleus of NGC 2207. The number of ULXs is comparable with the largest numbers of ULXs per unit mass in any galaxy yet reported. In this paper we report on these sources, and quantify how their locations correlate with the local star formation rates seen in spatially-resolved star formation rate density images that we have constructed using combinations of Galex FUV and Spitzer 24um images. We show that the numbers of ULXs are strongly correlated with the local star formation rate densities surrounding the sources, but that the luminosities of these sources are not strongly correlated with star formation rate density.

We present a Subaru weak lensing measurement of ACT-CL J0022.2-0036, one of the most luminous, high-redshift (z=0.81) Sunyaev-Zel'dovich (SZ) clusters discovered in the 268 deg^2 equatorial region survey of the Atacama Cosmology Telescope. For the weak lensing analysis using i'-band images, we use a model-fitting (Gauss-Laguerre shapelet) method to measure shapes of galaxy images, where we fit galaxy images in different exposures simultaneously to obtain best-fit ellipticities taking into account the different PSFs in each exposure. We also take into account the astrometric distortion effect on galaxy images by performing the model fitting in the world coordinate system. To select background galaxies behind the cluster at z=0.81, we use photometric redshift (photo-z) estimates for every galaxy derived from the co-added images of multi-passband Br'i'z'Y, with PSF matching/homogenization. After a photo-z cut for background galaxy selection, we detect the tangential weak lensing distortion signal with a total signal-to-noise ratio of about 3.7. By fitting a Navarro-Frenk-White model to the measured shear profile, we find the cluster mass to be M_200\bar{\rho}_m = [7.5^+3.2_-2.8(stat.)^+1.3_-0.6(sys.)] x 10^14 M_\odot/h. The weak lensing-derived mass is consistent with previous mass estimates based on the SZ observation, with assumptions of hydrostatic equilibrium and virial theorem, as well as with scaling relations between SZ signal and mass derived from weak lensing, X-ray, and velocity dispersion, within the measurement errors. We also show that the existence of ACT-CL J0022.2-0036 at z=0.81 is consistent with the cluster abundance prediction of the \Lambda-dominated cold dark matter structure formation model. We thus demonstrate the capability of Subaru-type ground-based images for studying weak lensing of high-redshift clusters.

We study the emission from the hot interstellar medium in a sample of nearby late type galaxies defined in Paper I. Our sample covers a broad range of star formation rates, from ~0.1 Msun/yr to ~17 Msun/yr and stellar masses, from ~3x10^8 Msun to ~6x10^10 Msun. We take special care of systematic effects and contamination from bright and faint compact sources. We find that in all galaxies at least one optically thin thermal emission component is present in the unresolved emission, with the average temperature of <kT>= 0.24 keV. In about ~1/3 of galaxies, a second, higher temperature component is required, with the <kT>= 0.71 keV. Although statistically significant variations in temperature between galaxies are present, we did not find any meaningful trends with the stellar mass or star formation rate of the host galaxy. The apparent luminosity of the diffuse emission in the 0.5-2 keV band linearly correlates with the star formation rate with the scale factor of Lx/SFR\approx 8.3x10^38 erg/s per Msun/yr, of which in average ~30-40% is likely produced by faint compact sources of various types. We attempt to estimate the bolometric luminosity of the gas and and obtained results differing by an order of magnitude, log(Lbol/SFR)\sim39-40, depending on whether intrinsic absorption in star-forming galaxies was allowed or not. Our theoretically most accurate, but in practice the most model dependent result for the intrinsic bolometric luminosity of ISM is Lbol/SFR\sim 1.5x10^40 erg/s per Msun/yr. Assuming that core collapse supernovae are the main source of energy, it implies that \epsilon_SN\sim5x10^-2 (E_SN/10^51)^-1 of mechanical energy of supernovae is converted into thermal energy of ISM.

In this paper we report accurate Chandra positions for two ultraluminous X-ray sources: NGC 7319-X4 at Right Ascension (RA) = 339.02917(2) deg, Declination (Dec) = 33.97476(2) deg and NGC 5474-X1 at RA = 211.24859(3) deg, Dec = 53.63584(3) deg. We perform bore-sight corrections on the Chandra X-ray Satellite observations of these sources to get to these accurate positions of the X-ray sources and match these positions with archival optical data from the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope. We do not find the optical counterparts: the limiting absolute magnitudes of the observations in the WFPC2 standard magnitude system are B = -7.9, V = -8.7 and I = -9.3 for NGC 7319-X4 and U = -6.4 for NGC 5474-X1. We report on the X-ray spectral properties and we find evidence for X-ray variability in NGC 5474-X1. Finally, we briefly discuss several options for the nature of these ULXs.

Based on a homogeneous set of X-ray, infrared and ultraviolet observations from Chandra, Spitzer, GALEX and 2MASS archives, we study populations of high-mass X-ray binaries (HMXBs) in a sample of 29 nearby star-forming galaxies and their relation with the star formation rate (SFR). In agreement with previous results, we find that HMXBs are a good tracer of the recent star formation activity in the host galaxy and their collective luminosity and number scale with the SFR, in particular, Lx~2.6 10^{39} SFR. However, the scaling relations still bear a rather large dispersion of ~0.4 dex, which we believe is of a physical origin. We present the catalog of 1057 X-ray sources detected within the $D25$ ellipse for galaxies of our sample and construct the average X-ray luminosity function (XLF) of HMXBs with substantially improved statistical accuracy and better control of systematic effects than achieved in previous studies. The XLF follows a power law with slope of 1.6 in the logLx~35-40 luminosity range with a moderately significant evidence for a break or cut-off at Lx~10^{40} erg/s. As before, we did not find any features at the Eddington limit for a neutron star or a stellar mass black hole. We discuss implications of our results for the theory of binary evolution. In particular we estimate the fraction of compact objects that once upon their lifetime experienced an X-ray active phase powered by accretion from a high mass companion and obtain a rather large number, fx~0.2 (0.1 Myr/tau_x) (tau_x is the life time of the X-ray active phase). This is ~4 orders of magnitude more frequent than in LMXBs. We also derive constrains on the mass distribution of the secondary star in HMXBs.

We study the relation between the X-ray luminosity of compact sources and the SFR of the host galaxy. Our sample includes 38 galaxies for which a uniform set of X-ray, infra-red and ultraviolet data from Chandra, Spitzer and GALEX has been collected. Our primary goals are (i) to obtain a more accurate calibration of the Lx-SFR relation and (ii) to understand the origin of the dispersion in the Lx-SFR relation observed in previous studies. Preliminary results of this project are reported below.