The Astropy Project supports and fosters the development of open-source and openly developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package, as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of interoperable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy Project.

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as flexible image transport system (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

We report the final optical identifications of the medium-depth (∼60 ks), contiguous (2 deg 2 ) XMM-Newton survey of the COSMOS field. XMM-Newton has detected ∼1800 X-ray sources down to limiting fluxes of ∼5 × 10 −16 , ∼3 × 10 −15 , and ∼7 × 10 −15 erg cm −2 s −1 in the 0.5-2 keV, 2-10 keV, and 5-10 keV bands, respectively (∼1 × 10 −15 , ∼6 × 10 −15 , and ∼1 × 10 −14 erg cm −2 s −1 , in the three bands, respectively, over 50% of the area). The work is complemented by an extensive collection of multiwavelength data from 24 μm to UV, available from the COSMOS survey, for each of the X-ray sources, including spectroscopic redshifts for 50% of the sample, and high-quality photometric redshifts for the rest. The XMM and multiwavelength flux limits are well matched: 1760 (98%) of the X-ray sources have optical counterparts, 1711 (∼95%) have IRAC counterparts, and 1394 (∼78%) have MIPS 24 μm detections. Thanks to the redshift completeness (almost 100%) we were able to constrain the high-luminosity tail of the X-ray luminosity function confirming that the peak of the number density of log L X > 44.5 active galactic nuclei (AGNs) is at z ∼ 2. Spectroscopically identified obscured and unobscured AGNs, as well as normal and star-forming galaxies, present well-defined optical and infrared properties. We devised a robust method to identify a sample of ∼150 high-redshift (z > 1), obscured AGN candidates for which optical spectroscopy is not available. We were able to determine that the fraction of the obscured AGN population at the highest (L X > 10 44 erg s −1 ) X-ray luminosity is ∼15%-30% when selection effects are taken into account, providing an important observational constraint for X-ray background synthesis. We studied in detail the optical spectrum and the overall spectral energy 348 No. 1, 2010 THE XMM-NEWTON WIDE-FIELD SURVEY IN THE COSMOS FIELD 349distribution of a prototypical Type 2 QSO, caught in a stage transitioning from being starburst dominated to AGN dominated, which was possible to isolate only thanks to the combination of X-ray and infrared observations.

Context. The COSMOS survey is a multiwavelength survey aimed to study the evolution of galaxies, AGN and large scale structures. Within this survey XMM-COSMOS a powerful tool to detect AGN and galaxy clusters. The XMM-COSMOS is a deep X-ray survey over the full 2 deg 2 of the COSMOS area. It consists of 55 XMM-Newton pointings for a total exposure of ∼1.5 Ms with an average vignetting-corrected depth of 40 ks across the field of view and a sky coverage of 2.13 deg 2 . Aims. We present the catalogue of point-like X-ray sources detected with the EPIC CCD cameras, the log N − log S relations and the X-ray colour-colour diagrams. Methods. The analysis was performed using the XMM-SAS data analysis package in the 0.5-2 keV, 2-10 keV and 5-10 keV energy bands. Source detection has been performed using a maximum likelihood technique especially designed for raster scan surveys. The completeness of the catalogue as well as log N − log S and source density maps have been calibrated using Monte Carlo simulations. Results. The catalogs contains a total of 1887 unique sources detected in at least one band with likelihood parameter det_ml > 10. The survey, which shows unprecedented homogeneity, has a flux limit of ∼1.7×10 −15 erg cm −2 s −1 , ∼9.3×10 −15 erg cm −2 s −1 and ∼1.3×10 −14 erg cm −2 s −1 over 90% of the area (1.92 deg 2 ) in the 0.5-2 keV, 2-10 keV and 5-10 keV energy band, respectively. Thanks to the rather homogeneous exposure over a large area, the derived log N − log S relations are very well determined over the flux range sampled by XMM-COSMOS. These relations have been compared with XRB synthesis models, which reproduce the observations with an agreement of ∼10% in the 5-10 keV and 2-10 keV band, while in the 0.5-2 keV band the agreement is of the order of ∼20%. The hard X-ray colors confirmed that the majority of the extragalactic sources in a bright subsample are actually type I or type II AGN. About 20% of the sources have a X-ray luminosity typical of AGN (L X > 10 42 erg/s) although they do not show any clear signature of nuclear activity in the optical spectrum.

A large population of heavily obscured, Compton-thick AGN is predicted by AGN synthesis models for the cosmic X-ray background and by the "relic" super-massive black-hole mass function measured from local bulges. However, even the deepest X-ray surveys are inefficient to search for these elusive AGN. Alternative selection criteria, combining mid-infrared with near-infrared and optical photometry, have instead been successful to pin-point a large population of Compton thick AGN. We take advantage of the deep Chandra and Spitzer coverage of a large area (more than 10 times the area covered by the Chandra deep fields, CDFs) in the COSMOS field, to extend the search of highly obscured, Compton-thick active nuclei to higher luminosity. These sources have low surface density and large samples can be provided only through large area surveys, like the COSMOS survey. We analyze the X-ray properties of COSMOS MIPS sources with 24µm fluxes higher than 550µJy. For the MIPS sources not directly detected in the Chandra images we produce stacked images in soft and hard X-rays bands. To estimate the fraction of Compton-thick AGN in the MIPS source population we compare the observed stacked count rates and hardness ratios to those predicted by detailed Monte Carlo simulations including both obscured AGN and star-forming galaxies. The volume density of Compton thick QSOs (logL(2-10keV)=44-45 ergs s −1 , or logλL λ (5.8µm)=44.79-46.18 ergs s −1 for a typical infrared to X-ray luminosity ratio) evaluated in this way is (4.8 ± 1.1) × 10 −6 Mpc −3 in the redshift bin 1.2-2.2. This density is ∼ 44% of that of all X-ray selected QSOs in the same redshift and luminosity bin, and it is consistent with the expectation of most up-to-date AGN synthesis models for the Cosmic X-ray background ). The density of lower luminosity Compton-thick AGN (logL(2-10keV)=43.5-44) at z=0.7-1.2 is (3.7 ± 1.1) × 10 −5 Mpc −3 , corresponding to ∼ 67% of that of X-ray selected AGN. The comparison between the fraction of infrared selected, Compton thick AGN to the X-ray selected, unobscured and moderately obscured AGN at high and low luminosity suggests that Compton-thick AGN follow a luminosity dependence similar to that discovered for Compton-thin AGN, becoming relatively rarer at high luminosities. We estimate that the fraction of AGN (unobscured, moderately obscured and Compton thick) to the total MIPS source population is 49 ± 10%, a value significantly higher than that previously estimated at similar 24µm fluxes. We discuss how our findings can constrain AGN feedback models.

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra program that has imaged the central 0.9 deg 2 of the COSMOS field down to limiting depths of 1.9×10 −16 erg cm −2 s −1 in the soft (0.5-2 keV) band, 7.3×10 −16 erg cm −2 s −1 in the hard (2-10 keV) band, and 5.7×10 −16 erg cm −2 s −1 in the full (0.5-10 keV) band. In this paper we report the i, K and 3.6µm identifications of the 1761 X-ray point sources. We use the likelihood ratio technique to derive the association of optical/infrared counterparts for 97% of the X-ray sources. For most of the remaining 3% , the presence of multiple counterparts or the faintness of the possible counterpart prevented a unique association. For only 10 X-ray sources we were not able to associate a counterpart, mostly due to the presence of a very bright field source close by. Only 2 sources are truly empty fields. The full catalog, including spectroscopic and photometric redshifts and classification described here in detail, is available on-line. Making use of the large number of X-ray sources, we update the "classic locus" of AGN defined 20 years ago in soft X-ray surveys and define a new locus containing 90% of the AGN in the survey with full band luminosity >10 42 erg s −1 . We present the linear fit between the total i band magnitude and the X-ray flux in the soft and hard band, drawn over 2 orders of magnitude in X-ray flux, obtained using the combined C-COSMOS and XMM-COSMOS samples. We focus on the X-ray to optical flux ratio (X/O) and we test its known correlation with redshift and luminosity, and a recently introduced anti-correlation with the concentration index (C). We find a strong anti-correlation (though the dispersion is of the order of 0.5 dex) between C and X/O, computed in the hard band, and that 90% of the obscured AGN in the sample with morphological information live in galaxies with regular morphology (bulgy and disky/spiral), suggesting that secular processes govern a significant fraction of the BH growth at X-ray luminosities of 10 43 − 10 44.5 erg s −1 . We also investigate the degree of obscuration of the sample, using the hardness ratio and we compare the X-ray color with the near-infrared-to-optical one.

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra program that covers the central contiguous ∼0.92 deg 2 of the COSMOS field. C-COSMOS is the result of a complex tiling, with every position being observed in up to six overlapping pointings (four overlapping pointings in most of the central ∼0.45 deg 2 area with the best exposure, and two overlapping pointings in most of the surrounding area, covering an additional ∼0.47 deg 2 ). Therefore, the full exploitation of the C-COSMOS data requires a dedicated and accurate analysis focused on three main issues: (1) maximizing the sensitivity when the point-spread function (PSF) changes strongly among different observations of the same source (from ∼1 arcsec up to ∼10 arcsec half-power radius); (2) resolving close pairs; and (3) obtaining the best source localization and count rate. We present here our treatment of four key analysis items: source detection, localization, photometry, and survey sensitivity. Our final procedure consists of a two step procedure: (1) a wavelet detection algorithm to find source candidates and (2) a maximum likelihood PSF fitting algorithm to evaluate the source count rates and the probability that each source candidate is a fluctuation of the background. We discuss the main characteristics of this procedure, which was the result of detailed comparisons between different detection algorithms and photometry tools, calibrated with extensive and dedicated simulations.

We study the spectral energy distributions and evolution of a large sample of optically selected quasars from the Sloan Digital Sky Survey that were observed in 323 Chandra images analyzed by the Chandra Multiwavelength Project. Our highest-confidence matched sample includes 1135 X-ray detected quasars in the redshift range 0.2 < z < 5.4, representing some 36 Msec of effective exposure. We provide catalogs of QSO properties, and describe our novel method of calculating X-ray flux upper limits and effective sky coverage. Spectroscopic redshifts are available for about 1/3 of the detected sample; elsewhere, redshifts are estimated photometrically. We detect 56 QSOs with redshift z > 3, substantially expanding the known sample. We find no evidence for evolution out to z ∼ 5 for either the X-ray photon index Γ or for the ratio of optical/UV to X-ray flux α ox . About 10% of detected QSOs show best-fit intrinsic absorbing columns greater than 10 22 cm −2 , but the fraction might reach ∼1/3 if most nondetections are absorbed. We confirm a significant correlation between α ox and optical luminosity, but it flattens or disappears for fainter (M B −23) active galactic nucleus (AGN) alone. We report significant hardening of Γ both toward higher X-ray luminosity, and for relatively X-ray loud quasars. These trends may represent a relative increase in nonthermal X-ray emission, and our findings thereby strengthen analogies between Galactic black hole binaries and AGN. For uniformly selected subsamples of narrow-line Seyfert 1s and narrow absorption line QSOs, we find no evidence for unusual distributions of either α ox or Γ.