The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

The tomographic Alcock-Paczynski (AP) method can result in tight cosmological constraints by using small and intermediate clustering scales of the large scale structure (LSS) of the galaxy distribution. By focusing on the redshift dependence, the AP distortion can be distinguished from the distortions produced by the redshift space distortions (RSD). In this work, we combine the tomographic AP method with other recent observational datasets of SNIa+BAO+CMB+$H_0$ to reconstruct the dark energy equation-of-state $w$ in a non-parametric form. The result favors a dynamical DE at $z\lesssim1$, and shows a mild deviation ($\lesssim2\sigma$) from $w=-1$ at $z=0.5-0.7$. We find the addition of the AP method improves the low redshift ($z\lesssim0.7$) constraint by $\sim50\%$.

The tomographic Alcock-Paczynski (AP) test is a robust large-scale structure (LSS) measurement that receives little contamination from the redshift space distortion (RSD). It has placed tight cosmological constraints by using small and intermediate clustering scales of the LSS data. However, previous works have neglected the cross-correlation among different redshift bins, which could cause the statistical uncertainty being underestimated by $\sim$20\%. In this work, we further improve this method by including this multi-redshifts full correlation. We apply it to the SDSS DR12 galaxies sample and find out that, for $\Lambda$CDM, the combination of AP with the Planck+BAO dataset slightly reduces (within 1-$\sigma$) $\Omega_m$ to $0.304\pm0.007$ (68.3\% CL). This then leads to a larger $H_0$ and also mildly affects $\Omega_b h^2$, $n_s$ and the derived parameters $z_*$, $r_*$, $z_{re}$ but not $\tau$, $A_s$ and $\sigma_8$. For the flat $w$CDM model, our measurement gives $\Omega_m=0.301\pm 0.010$ and $w=-1.090\pm 0.047$, where the additional AP measurement reduces the error budget by $\sim 25\%$. When including more parameters into the analysis, the AP method also improves the constraints on $\Omega_k$, $\sum m_\mu$, $N_{\rm eff}$ by $20-30\%$. Early universe parameters such as $dn_s/d{\rm ln}k$ and $r$, however, are unaffected. Assuming the dark energy equation of state $w=w_0+w_a \frac{z}{1+z}$, the Planck+BAO+SNIa+$H_0$+AP datasets prefer a dynamical dark energy at $\approx1.5 \sigma$ CL. Finally, we forecast the cosmological constraints expected from the DESI galaxy survey and find that combining AP with CMB+BAO method would improve the $w_0$-$w_a$ constraint by a factor of $\sim 10$.

DArk Matter Particle Explorer (DAMPE) is a space-borne apparatus for detecting the high-energy cosmic-rays like electrons, $\gamma$-rays, protons and heavy-ions. Plastic Scintillator Detector (PSD) is the top-most sub-detector of the DAMPE. The PSD is designed to measure the charge of incident high-energy particles and it also serves as a veto detector for discriminating $\gamma$-rays from charged particles. In this paper, PSD on-orbit calibration procedure is described, which includes five steps of pedestal, dynode correlation, response to minimum-ionizing particles (MIPs), light attenuation function and energy reconstruction. A method for reconstructing the charge of incident high energy cosmic-ray particles is introduced. The detection efficiency of each PSD strip is verified to be above 99.5%, the total efficiency of the PSD for charged particles is above 99.99%.

The $\beta$-skeleton is a mathematical method to construct graphs from a set of points that has been widely applied in the areas of image analysis, machine learning, visual perception, and pattern recognition. In this work, we apply the $\beta$-skeleton to study the cosmic web. We use this tool on observed and simulated data to identify the filamentary structures and characterize the statistical properties of the skeleton. In particular, we compare the $\beta$-skeletons built from SDSS-III galaxies to those obtained from MD-PATCHY mocks, and also to mocks directly built from the Big MultiDark $N$-body simulation. We find that the $\beta$-skeleton is able to reveal the underlying structures in observed and simulated samples without any parameter fine-tuning. A different degree of sparseness can be obtained by adjusting the value of $\beta$; in addition, the statistical properties of the length and direction of the skeleton connections show a clear dependence on redshift space distortions (RSDs), cosmological effects and galaxy bias. We also find that the $N$-body simulation accurately reproduces the RSD effect in the data, while the MD-PATCHY mocks appear to underestimate its magnitude. Our proof-of-concept study shows that the statistical properties of the $\beta$-skeleton can be used to probe cosmological parameters and galaxy evolution.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

he DArk Matter Particle Explorer (DAMPE) is a general purposed satellite-borne high energy $\gamma-$ray and cosmic ray detector, and among the scientific objectives of DAMPE are the searches for the origin of cosmic rays and an understanding of Dark Matter particles. As one of the four detectors in DAMPE, the Plastic Scintillator Detector (PSD) plays an important role in the particle charge measurement and the photons/electrons separation. The PSD has 82 modules, each consists of a long organic plastic scintillator bar and two PMTs at both ends for readout, in two layers and covers an overall active area larger than 82 cm $\times$ 82 cm. It can identify the charge states for relativistic ions from H to Fe, and the detector efficiency for Z=1 particles can reach 0.9999. The PSD has been successfully launched with DAMPE on Dec. 17, 2015. In this paper, the design, the assembly, the qualification tests of the PSD and some of the performance measured on the ground have been described in detail.

We present the first observation, analysis and modeling of solar coronal twin jets, which occurred after a preceding jet. Detailed analysis on the kinetics of the preceding jet reveals its blowout-jet nature, which resembles the one studied in Liu et al. 2014. However the erupting process and kinetics of the twin jets appear to be different from the preceding one. In lack of the detailed information on the magnetic fields in the twin jet region, we instead use a numerical simulation using a three-dimensional (3D) MHD model as described in Fang et al. 2014, and find that in the simulation a pair of twin jets form due to reconnection between the ambient open fields and a highly twisted sigmoidal magnetic flux which is the outcome of the further evolution of the magnetic fields following the preceding blowout jet. Based on the similarity between the synthesized and observed emission we propose this mechanism as a possible explanation for the observed twin jets. Combining our observation and simulation, we suggest that with continuous energy transport from the subsurface convection zone into the corona, solar coronal twin jets could be generated in the same fashion addressed above.

We carry out a magneto-hydrodynamic (MHD) simulation of convective dynamo in the rotating solar convective envelope driven by the solar radiative diffusive heat flux. The simulation is similar to that reported in Fan & Fang (2014) but with further reduced viscosity and magnetic diffusion. The resulting convective dynamo produces a large scale mean field that exhibits similar irregular cyclic behavior and polarity reversals, and self-consistently maintains a solar-like differential rotation. The main driver for the solar-like differential rotation (with faster rotating equator) is a net outward transport of angular momentum away from the rotation axis by the Reynolds stress, and we found that this transport is enhanced with reduced viscosity and magnetic diffusion.

$\delta$-sunspots, with highly complex magnetic structures, are very productive in energetic eruptive events, such as X-class flares and homologous eruptions. We here study the formation of such complex magnetic structures by numerical simulations of magnetic flux emergence from the convection zone into the corona in an active-region-scale domain. In our simulation, two pairs of bipolar sunspots form on the surface, originating from two buoyant segments of a single subsurface twisted flux rope, following the approach of Toriumi et al. (2014). Expansion and rotation of the emerging fields in the two bipoles drive the two opposite polarities into each other with apparent rotating motion, producing a compact $\delta$-sunspot with a sharp polarity inversion line. The formation of the $\delta$-sunspot in such a realistic-scale domain produces emerging patterns similar to those formed in observations, e.g. the inverted polarity against Hale's law, the curvilinear motion of the spot, strong transverse field with highly sheared magnetic and velocity fields at the PIL. Strong current builds up at the PIL, giving rise to reconnection, which produces a complex coronal magnetic connectivity with non-potential fields in the Delta-spot overlaid by more relaxed fields connecting the two polarities at the two ends.

We study the formation of coronal jets through numerical simulation of the emergence of a twisted magnetic flux rope into a pre-existing open magnetic field. Reconnection inside the emerging flux rope in addition to that between the emerging and pre-existing fields give rise to the violent eruption studied. The simulated event closely resembles the coronal jets ubiquitously observed by Hinode/XRT and demonstrates that heated plasma is driven into the extended atmosphere above. Thermal conduction implemented in the model allows us to qualitatively compare simulated and observed emission from such events. We find that untwisting field lines after the reconnection drive spinning outflows of plasma in the jet column. The Poynting flux in the simulated jet is dominated by the untwisting motions of the magnetic fields loaded with high-density plasma. The simulated jet is comprised of spires of untwisting field that are loaded with a mixture of cold and hot plasma and exhibit rotational motion of order 20 km/s and match contemporary observations.

We report the results of a magneto-hydrodynamic (MHD) simulation of a convective dynamo in a model solar convective envelope driven by the solar radiative diffusive heat flux. The convective dynamo produces a large-scale mean magnetic field that exhibits irregular cyclic behavior with oscillation time scales ranging from about 5 to 15 years and undergoes irregular polarity reversals. The mean axisymmetric toroidal magnetic field is of opposite signs in the two hemispheres and is concentrated at the bottom of the convection zone. The presence of the magnetic fields is found to play an important role in the self-consistent maintenance of a solar-like differential rotation in the convective dynamo model. Without the magnetic fields, the convective flows drive a differential rotation with a faster rotating polar region. In the midst of magneto-convection, we found emergence of strong super-equipartition flux bundles at the surface, exhibiting properties that are similar to emerging solar active regions.

We describe a 3D finite-difference spherical anelastic MHD (FSAM) code for modeling the subsonic dynamic processes in the solar convective envelope. A comparison of this code with the widely used global spectral anlastic MHD code, ASH (Anelastic Spherical Harmonics), shows that FSAM produces convective flows with statistical properties and mean flows similar to the ASH results. Using FSAM, we first simulate the rotating solar convection in a partial spherical shell domain and obtain a statistically steady, giant-cell convective flow with a solar-like differential rotation. We then insert buoyant toroidal flux tubes near the bottom of the convecting envelope and simulate the rise of the flux tubes in the presence of the giant cell convection. We find that for buoyant flux tubes with an initial field strength of 100 kG, the magnetic buoyancy largely determines the rise of the tubes although strong down flows produce significant undulation and distortion to the shape of the emerging $\Omega$-shaped loops. The convective flows significantly reduce the rise time it takes for the apex of the flux tube to reach the top. For the weakly twisted and untwisted cases we simulated, the apex portion is found to rise nearly radially to the top in about a month, and produce an emerging region (at a depth of about 30 Mm below the photosphere) with an overall tilt angle consistent with the active region tilts, although the emergence pattern is more complex compared to the case without convection. Near the top boundary at a depth of about 30 Mm, the emerging flux shows a retrograde zonal flow of about 345 m/s relative to the mean flow at that latitude.

We examine a simulation of flux emergence and cancellation, which shows a complex sequence of processes that accumulate free magnetic energy in the solar corona essential for the eruptive events such as coronal mass ejections (CMEs), filament eruptions and flares. The flow velocity at the surface and in the corona shows a consistent shearing pattern along the polarity inversion line (PIL), which together with the rotation of the magnetic polarities, builds up the magnetic shear. Tether-cutting reconnection above the PIL then produces longer sheared magnetic field lines that extend higher into the corona, where a sigmoidal structure forms. Most significantly, reconnection and upward energy-flux transfer are found to occur even as magnetic flux is submerging and appears to cancel at the photosphere. A comparison of the simulated coronal field with the corresponding coronal potential field graphically shows the development of nonpotential fields during the emergence of the magnetic flux and formation of sunspots.

We simulate the buoyant rise of a magnetic flux rope from the solar convection zone into the corona to better understand the energetic coupling of the solar interior to the corona. The magnetohydrodynamic model addresses the physics of radiative cooling, coronal heating and ionization, which allow us to produce a more realistic model of the solar atmosphere. The simulation illustrates the process by which magnetic flux emerges at the photosphere and coalesces to form two large concentrations of opposite polarities. We find that the large-scale convective motion in the convection zone is critical to form and maintain sunspots, while the horizontal converging flows in the near surface layer prevent the concentrated polarities from separating. The foot points of the sunspots in the convection zone exhibit a coherent rotation motion, resulting in the increasing helicity of the coronal field. Here, the local configuration of the convection causes the convergence of opposite polarities of magnetic flux with a shearing flow along the polarity inversion line. During the rising of the flux rope, the magnetic energy is first injected through the photosphere by the emergence, followed by energy transport by horizontal flows, after which the energy is subducted back to the convection zone by the submerging flows.

We present the broad-band Spectral Energy Distributions (SEDs) of the largest available highly (72%) complete spectroscopic sample of mid-infrared (MIR) selected galaxies and AGN at intermediate redshift. The sample contains 203 extragalactic sources from the 15-micron survey in the ELAIS-SWIRE field S1, all with measured spectroscopic redshift. Most of these sources have full multi-wavelength coverage from the far-UV to the far-infrared and lie in the redshift range 0.1<z<1.3. Due to its size, this sample allows us for the first time to characterise the spectral properties of the sources responsible for the strong evolution observed in the MIR. Based on SED-fitting technique we have classified the MIR sources, identifying AGN signatures in about 50% of them. This fraction is significantly higher than that derived from optical spectroscopy (~29%) and is due in particular to the identification of AGN activity in objects spectroscopically classified as galaxies. It is likely that in most of our objects, the AGN is either obscured or of low-luminosity, and thus it does not dominate the energetic output at any wavelength, except in the MIR, showing up just in the range where the host galaxy SED has a minimum. The fraction of AGN strongly depends on the flux density, with that derived through the SED-fitting being about 20% at S(15)~0.5-1 mJy and gradually increasing up to 100% at S(15)>10 mJy, while that obtained from optical spectroscopy never being >30%, even at the higher flux densities. The results of this work will be very useful for updating all the models aimed at interpreting the deep infrared survey data and, in particular, for constraining the nature and the role of dust-obscured systems in the intermediate/high-redshift Universe.

We present the optical-to-submm spectral energy distributions for 33 radio & mid-IR identified submillimetre galaxies discovered via the SHADES 850 micron SCUBA imaging in the Subaru-XMM Deep Field (SXDF). Optical data for the sources comes from the Subaru-XMM Deep Field (SXDF) and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6 and 4.5 micron fluxes. We then fit spectral energy distribution (SED) templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The infrared template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submm sources and on the SHADES SXDF field. Our template fitting finds that AGN, present in about 10% of sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fit by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submm galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.

This paper presents galaxy source counts at 24 microns in the six Spitzer Wide-field InfraRed Extragalactic (SWIRE) fields. The source counts are compared to counts in other fields, and to model predictions that have been updated since the launch of Spitzer. This analysis confirms a very steep rise in the Euclidean-normalized differential number counts between 2 mJy and 0.3 mJy. Variations in the counts between fields show the effects of sample variance in the flux range 0.5-10 mJy, up to 100% larger than Poisson errors. Nonetheless, a "shoulder" in the normalized counts persists at around 3 mJy. The peak of the normalized counts at 0.3 mJy is higher and narrower than most models predict. In the ELAIS N1 field, the 24 micron data are combined with Spitzer-IRAC data and five-band optical imaging, and these bandmerged data are fit with photometric redshift templates. Above 1 mJy the counts are dominated by galaxies at z less than 0.3. By 300 microJy, about 25% are between z ~ 0.3-0.8, and a significant fraction are at z ~ 1.3-2. At low redshifts the counts are dominated by spirals, and starbursts rise in number density to outnumber the spirals' contribution to the counts below 1 mJy.

We investigate the clustering of galaxies selected in the 3.6 micron band of the Spitzer Wide-area Infrared Extragalactic (SWIRE) legacy survey. The angular two-point correlation function is calculated for eleven samples with flux limits of S_3.6 > 4-400 mujy, over an 8 square degree field. The angular clustering strength is measured at >5-sigma significance at all flux limits, with amplitudes of A=(0.49-29)\times10^{-3} at one degree, for a power-law model, A\theta^{-0.8}. We estimate the redshift distributions of the samples using phenomological models, simulations and photometric redshifts, and so derive the spatial correlation lengths. We compare our results with the GalICS (Galaxies In Cosmological Simulations) models of galaxy evolution and with parameterized models of clustering evolution. The GalICS simulations are consistent with our angular correlation functions, but fail to match the spatial clustering inferred from the phenomological models or the photometric redshifts. We find that the uncertainties in the redshift distributions of our samples dominate the statistical errors in our estimates of the spatial clustering. At low redshifts (median z<0.5) the comoving correlation length is approximately constant, r_0=6.1\pm0.5h^{-1} Mpc, and then decreases with increasing redshift to a value of 2.9\pm0.3h^{-1} Mpc for the faintest sample, for which the median redshift is z=1. We suggest that this trend can be attributed to a decrease in the average galaxy and halo mass in the fainter flux-limited samples, corresponding to changes in the relative numbers of early- and late-type galaxies. However, we cannot rule out strong evolution of the correlation length over 0.5<z<1.

We use a simple optical/infrared (IR) photometric selection of high-redshift QSOs that identifies a Lyman Break in the optical photometry and requires a red IR color to distinguish QSOs from common interlopers. The search yields 100 z~3 (U-dropout) QSO candidates with 19<r'<22 over 11.7 deg^2 in the ELAIS-N1 (EN1) and ELAIS-N2 (EN2) fields of the Spitzer Wide-area Infrared Extragalactic (SWIRE) Legacy Survey. The z~3 selection is reliable, with spectroscopic follow-up of 10 candidates confirming they are all QSOs at 2.83<z<3.44. We find that our z~4$ (g'-dropout) sample suffers from both unreliability and incompleteness but present 7 previously unidentified QSOs at 3.50<z<3.89. Detailed simulations show our z~3 completeness to be ~80-90% from 3.0<z<3.5, significantly better than the ~30-80% completeness of the SDSS at these redshifts. The resulting luminosity function extends two magnitudes fainter than SDSS and has a faint end slope of beta=-1.42 +- 0.15, consistent with values measured at lower redshift. Therefore, we see no evidence for evolution of the faint end slope of the QSO luminosity function. Including the SDSS QSO sample, we have now directly measured the space density of QSOs responsible for ~70% of the QSO UV luminosity density at z~3. We derive a maximum rate of HI photoionization from QSOs at z~3.2, Gamma = 4.8x10^-13 s^-1, about half of the total rate inferred through studies of the Ly-alpha forest. Therefore, star-forming galaxies and QSOs must contribute comparably to the photoionization of HI in the intergalactic medium at z~3.

By combining the VVDS with the SWIRE data, we have built the currently largest spectroscopic sample of galaxies selected in the rest-frame near-infrared. These allow us to investigate, for the first time using spectroscopic redshifts, the clustering evolution of galaxies selected from their rest-frame near-infrared luminosity in the redshift range 0.2<z<2.1. Therefore we use the projected two-point correlation function wp(rp) to study the three dimensional clustering properties of galaxies detected at 3.6 and 4.5 microns with IRAC and for which we have spectroscopic redshifts from first epoch VVDS. We find that in the flux limited samples at 3.6 and 4.5 microns, the clustering length does not change from redshift ~2 to the present. In addition, we find that the clustering lengths are systematically higher when galaxy samples are selected from increasingly redder wavelengths. These values are larger than the typical values found for I-band selected galaxies in the same redshift range, but we find that the difference in clustering length between I-band and 3.6/4.5 microns selected samples is decreasing with increasing redshift to become comparable at z~1.5. We interpret this as evidence that galaxies with older stellar populations and galaxies actively forming stars reside in comparably overdense environments at epochs earlier than z~1.5. The increasing difference in clustering length observed between rest-frame UV-optical and infrared selected samples could then be an indication that star formation is gradually shifting to lower density regions as cosmic time increases.

We study the ultraviolet to far-infrared (hereafter UV-to-IR) SEDs of a sample of intermediate redshift (0.2 < z < 0.7) UV-selected galaxies from the ELAIS-N1 and ELAIS-N2 fields by fitting a multi-wavelength dataset to a library of GRASIL templates. Star formation related properties of the galaxies are derived from the library of models by using the Bayesian statistics. We find a decreasing presence of galaxies with low attenuation and low total luminosity as redshift decreases, which does not hold for high total luminosity galaxies. In addition the dust attenuation of low mass galaxies increases as redshift decreases, and this trend seems to disappear for galaxies with M* > 10^11 M_sun. This result is consistent with a mass dependent evolution of the dust to gas ratio, which could be driven by a mass dependent efficiency of star formation in star forming galaxies. The specific star formation rates (SSFR) decrease with increasing stellar mass at all redshifts, and for a given stellar mass the SSFR decreases with decreasing redshift. The differences in the slope of the M*--SSFR relation found between this work and others at similar redshift could be explained by the adopted selection criteria of the samples which, for a UV selected sample, favours blue, star forming galaxies.

We present measurements of the spatial clustering of ultraluminous infrared galaxies in two redshift intervals, 1.5<z<2.0 and 2<z<3. Both samples cluster strongly, with r_0=14.40+/-1.99 h^-1 Mpc for the 2<z<3 sample, and r_0=9.40+/-2.24 h^-1 Mpc for the 1.5<z<2.0 sample, making them among the most biased galaxies at these epochs. These clustering amplitudes are consistent with both populations residing in dark matter haloes with masses of ~6x10^13 solar masses. We infer that a minimum dark matter halo mass is an important factor for all forms of luminous, obscured activity in galaxies at z>1. Adopting plausible models for the growth of DM haloes with redshift, then the haloes hosting the 2<z<3 sample will likely host the richest clusters of galaxies at z=0, whereas the haloes hosting the 1.5<z<2.0 sample will likely host poor to rich clusters at z=0.

Our goal is to probe the populations of obscured and unobscured AGN investigating their optical-IR and X-ray properties as a function of X-ray flux, luminosity and redshift within a hard X-ray selected sample of 136 X-ray sources in the XMM Medium Deep Survey (XMDS) with wide multiwavelength coverage. The XMDS area is covered with optical photometry from the VVDS and CFHTLS surveys and infrared Spitzer data. Based on the X-ray luminosity and X-ray to optical ratio, 132 sources are likely AGN, of which 122 have unambiguous optical - IR identification. The observed optical and IR spectral energy distributions of sources are fitted with AGN/galaxy templates in order to classify them and compute photometric redshifts. 70% of the AGN are fitted by a type 2 AGN or a star forming galaxy template and are grouped together in a single class of ``optically obscured'' AGN. They have ``red'' optical colors and generally show significant X-ray absorption from X-ray spectra or hardness ratios (N$_H > 10^{22}$ cm$^{-2}$). Sources with SEDs typical of type 1 AGN have ``blue'' optical colors and exhibit X-ray absorption in about 30% of cases. We performed a stacking analysis for obscured and type 1 AGN. The stacked X-ray spectrum of obscured AGN is flatter than that of type 1 AGN and has an average spectral slope of Gamma = 1.6. The subsample of objects fitted by a galaxy template has an even harder stacked spectrum, with Gamma = 1.2 - 1.3. The obscured fraction is larger at lower fluxes, lower redshifts and lower luminosities. X-ray absorption is less common than ``optical'' obscuration and its incidence is nearly constant with redshift and luminosity. This implies that X-ray absorption is not necessarily related to optical obscuration.

We present the spectral energy distributions (SEDs) of a hard X-ray selected sample. The sample contains 136 sources with F(2-10 keV)>10^-14 erg/cm^2/s and 132 are AGNs. The sources are detected in a 1 square degree area of the XMM-Newton-Medium Deep Survey where optical data from the VVDS, CFHTLS surveys, and infrared data from the SWIRE survey are available. Based on a SED fitting technique we derive photometric redshifts with sigma(1+z)=0.11 and 6% of outliers and identify AGN signatures in 83% of the objects. This fraction is higher than derived when a spectroscopic classification is available. The remaining 17+9-6% of AGNs shows star-forming galaxy SEDs (SF class). The sources with AGN signatures are divided in two classes, AGN1 (33+6-1%) and AGN2 (50+6-11). The AGN1 and AGN2 classes include sources whose SEDs are fitted by type 1 and type 2 AGN templates, respectively. On average, AGN1s show soft X-ray spectra, consistent with being unabsorbed, while AGN2s and SFs show hard X-ray spectra, consistent with being absorbed. The analysis of the average SEDs as a function of X-ray luminosity shows a reddening of the IR SEDs, consistent with a decreasing contribution from the host galaxy at higher luminosities. The AGNs in the SF classes are likely obscured in the mid-infrared, as suggested by their low L(3-20micron)/Lcorr(0.5-10 keV) ratios. We confirm the previously found correlation for AGNs between the radio luminosity and the X-ray and the mid-infrared luminosities. The X-ray-radio correlation can be used to identify heavily absorbed AGNs. However, the estimated radio fluxes for the missing AGN population responsible for the bulk of the background at E>10 keV are too faint to be detected even in the deepest current radio surveys.

We study dust attenuation and stellar mass of $\rm z\sim 0.6$ star-forming galaxies using new SWIRE observations in IR and GALEX observations in UV. Two samples are selected from the SWIRE and GALEX source catalogs in the SWIRE/GALEX field ELAIS-N1-00 ($\Omega = 0.8$ deg$^2$). The UV selected sample has 600 galaxies with photometric redshift (hereafter photo-z) $0.5 \leq z \leq 0.7$ and NUV$\leq 23.5$ (corresponding to $\rm L_{FUV} \geq 10^{9.6} L_\sun$). The IR selected sample contains 430 galaxies with $f_{24\mu m} \geq 0.2$ mJy ($\rm L_{dust} \geq 10^{10.8} L_\sun$) in the same photo-z range. It is found that the mean $\rm L_{dust}/L_{FUV}$ ratios of the z=0.6 UV galaxies are consistent with that of their z=0 counterparts of the same $\rm L_{FUV}$. For IR galaxies, the mean $\rm L_{dust}/L_{FUV}$ ratios of the z=0.6 LIRGs ($\rm L_{dust} \sim 10^{11} L_\sun$) are about a factor of 2 lower than local LIRGs, whereas z=0.6 ULIRGs ($\rm L_{dust} \sim 10^{12} L_\sun$) have the same mean $\rm L_{dust}/L_{FUV}$ ratios as their local counterparts. This is consistent with the hypothesis that the dominant component of LIRG population has changed from large, gas rich spirals at z$>0.5$ to major-mergers at z=0. The stellar mass of z=0.6 UV galaxies of $\rm L_{FUV} \leq 10^{10.2} L_\sun$ is about a factor 2 less than their local counterparts of the same luminosity, indicating growth of these galaxies. The mass of z=0.6 UV lunmous galaxies (UVLGs: $\rm L_{FUV} > 10^{10.2} L_\sun$) and IR selected galaxies, which are nearly exclusively LIRGs and ULIRGs, is the same as their local counterparts.

We present details of the discovery of XLSSJ022303.0-043622, a z=1.2 cluster of galaxies. This cluster was identified from its X-ray properties and selected as a z>1 candidate from its optical/near-IR characteristics in the XMM Large-Scale Structure Survey (XMM-LSS). It is the most distant system discovered in the survey to date. We present ground-based optical and near IR observations of the system carried out as part of the XMM-LSS survey. The cluster has a bolometric X-ray luminosity of 1.1 +/- 0.7 x 10^44 erg/s, fainter than most other known z>1 X-ray selected clusters. In the optical it has a remarkably compact core, with at least a dozen galaxies inside a 125 kpc radius circle centred on the X-ray position. Most of the galaxies within the core, and those spectroscopically confirmed to be cluster members, have stellar masses similar to those of massive cluster galaxies at low redshift. They have colours comparable to those of galaxies in other z>1 clusters, consistent with showing little sign of strong ongoing star formation. The bulk of the star formation within the galaxies appears to have ceased at least 1.5 Gyr before the observed epoch. Our results are consistent with massive cluster galaxies forming at z>1 and passively evolving thereafter. We also show that the system is straightforwardly identified in Spitzer/IRAC 3.6 and 4.5 micron data obtained by the SWIRE survey emphasising the power and utility of joint XMM and Spitzer searches for the most distant clusters.

We present the rest-frame optical and infrared colours of a complete sample of 1114 z<0.3 galaxies from the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) and the Sloan Digital Sky Survey (SDSS). We discuss the optical and infrared colours of our sample and analyse in detail the contribution of dusty star-forming galaxies and AGN to optically selected red sequence galaxies. We propose that the optical (g-r) colour and infrared log(L_{24}/L_{3.6}) colour of galaxies in our sample are determined primarily by a bulge-to-disk ratio. The (g-r) colour is found to be sensitive to the bulge-to-disk ratio for disk-dominated galaxies, whereas the log(L_{24}/L_{3.6}) colour is more sensitive for bulge-dominated systems. We identify ~18% (195 sources) of our sample as having red optical colours and infrared excess. Typically, the infrared luminosities of these galaxies are found to be at the high end of star-forming galaxies with blue optical colours. Using emission line diagnostic diagrams, 78 are found to have an AGN contribution, and 117 are identified as star-forming systems. The red (g-r) colour of the star-forming galaxies could be explained by extinction. However, their high optical luminosities cannot. We conclude that they have a significant bulge component. The number densities of optically red star-forming galaxies are found to correspond to ~13% of the total number density of our sample. In addition, these systems contribute ~13% of the total optical luminosity density, and 28% of the total infrared luminosity density of our SWIRE/SDSS sample. These objects may reduce the need for "dry-mergers".

We construct rest-frame luminosity functions at 3.6, 4.5, 5.8, 8 and 24 microns over the redshift range 0<z<2 for galaxies and 0<z<4 for optical QSOs, using optical and infrared data from the Spitzer Wide-area InfraRed Extragalactic survey. The 3.6 and 4.5 micron galaxy LFs show evidence for moderate positive luminosity evolution up to z~1.5, consistent with the passive ageing of evolved stellar populations. Their comoving luminosity density was found to evolve passively, gradually increasing out to z~0.5-1 but flattening, or even declining, at higher redshift. Conversely, the 24 micron galaxy LF, which is more sensitive to obscured star formation and/or AGN activity, undergoes strong positive evolution, with the derived IR energy density and SFR density proportional to (1+z)^gamma with gamma=4.5^{+0.7}_{-0.6} and the majority of this evolution occurring since z~1. Optical QSOs, however, show positive luminosity evolution in all bands, out to the highest redshifts (3<z<4). Modelling as L* proportional to (1+z)^gamma gave gamma=1.3^{+0.1}_{-0.1} at 3.6 micron, gamma=1.0^{+0.1}_{-0.1} at 4.5 micron and stronger evolution at the longer wavelengths (5.8, 8 and 24 micron), of gamma~3. Comparison of the galaxy LFs to predictions from a semi-analytic model based on CDM indicate that an IMF skewed towards higher mass star formation in bursts compared to locally is preferred. As a result the currently inferred massive star formation rates in distant sub-mm sources may require substantial downwards revision.

We fit a parametric model comprising a mixture of multi-dimensional Gaussian functions to the 3.6 to 8um colour and optical photo-z distribution of galaxy populations in the ELAIS-N1 and Lockman Fields of SWIRE. For 16,698 sources in ELAIS-N1 we find our data are best modelled (in the sense of the Bayesian Information Criterion) by the sum of four Gaussian distributions or modes (C_a, C_b, C_c and C_d). We compare the fit of our empirical model with predictions from existing semi-analytic and phenomological models. We infer that our empirical model provides a better description of the mid-infrared colour distribution of the SWIRE survey than these existing models. This colour distribution test is thus a powerful model discriminator and complementary to comparisons of number counts. We use our model to provide a galaxy classification scheme and explore the nature of the galaxies in the different modes of the model. C_a consists of dusty star-forming systems such as ULIRG's. Low redshift late-type spirals are found in C_b, where PAH emission dominates at 8um. C_c consists of dusty starburst systems at intermediate redshifts. Low redshift early-type spirals and ellipticals dominate C_d. We thus find a greater variety of galaxy types than one can with optical photometry alone. Finally we develop a new technique to identify unusual objects, and find a selection of outliers with very red IRAC colours. These objects are not detected in the optical, but have very strong detections in the mid-infrared. These sources are modelled as dust-enshrouded, strongly obscured AGN, where the high mid-infrared emission may either be attributed to dust heated by the AGN or substantial star-formation. These sources have z_ph ~ 2-4, making them incredibly infrared luminous, with a L_IR ~ 10^(12.6-14.1) L_sun.

We present measurements of the spatial clustering of galaxies with stellar masses >10^11Msun, infrared luminosities >10^12 Lsun, and star formation rates >200Msun per year in two redshift intervals; 1.5<z<2.0 and 2<z<3. Both samples cluster very strongly, with spatial correlation lengths of r_0=14.40+/-1.99 h^-1Mpc for the 2<z<3 sample, and r_0=9.40+/-2.24 h^-1Mpc for the 1.5<z<2.0 sample. These clustering amplitudes are consistent with both populations residing in dark matter haloes with masses of ~6x10^13Msun, making them among the most biased galaxies at these epochs. We infer, from this and previous results, that a minimum dark matter halo mass is an important factor for all forms of luminous, obscured activity in galaxies at z>1, both starbursts and AGN. Adopting plausible models for the growth of DM haloes with redshift, then the haloes hosting the 2<z<3 sample will likely host the richest clusters of galaxies at z=0, whereas the haloes hosting the 1.5<z<2.0 sample will likely host poor to rich clusters at z=0. We conclude that ULIRGs at z>1 signpost stellar buildup in galaxies that will reside in clusters at z=0, with ULIRGs at increasing redshifts signposting the buildup of stars in galaxies that will reside in increasingly rich clusters.

We use a simple optical/infrared photometric selection of high redshift QSOs which identifies a Lyman Break in the optical and requires a red IRAC color to distinguish QSOs from common interlopers. We find 100 U-dropout (z=3) QSO candidates with r'<22 within 11.2 deg^2 in the ELAIS-N1 & ELAIS-N2 fields in the Spitzer Wide-Area Infrared Extragalactic (SWIRE) Legacy Survey. Spectroscopy of 10 candidates shows that they are all QSOs with 2.83<z<3.44. We use detailed simulations which incorporate variations in QSO SEDs, IGM transmission and imaging depth to derive a completeness of 85-90% between 3.0<z<3.4. The resulting luminosity function is two magnitudes fainter than SDSS and, when combined with those data, gives a faint end slop \beta = 1.62 \pm 0.18, consistent with measurements at z<2 and steeper than initial measurements at the same redshift.

We discuss spectral energy distributions, photometric redshifts, redshift distributions, luminosity functions, source-counts and the far infrared to optical luminosity ratio for sources in the SWIRE Legacy Survey. The spectral energy distributions of selected SWIRE sources are modelled in terms of a simple set of galaxy and quasar templates in the optical and near infrared, and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst, and AGN dust torus) in the mid infrared. The optical data, together with the IRAC 3.6 and 4.5 mu data, have been used to determine photometric redshifts. For galaxies with known spectroscopic redshifts there is a notable improvement in the photometric redshift when the IRAC data are used, with a reduction in the rms scatter from 10% in (1+z) to 5%. While further spectroscopic data are needed to confirm this result, the prospect of determining good photometric redshifts for the 2 million extragalactic objects in SWIRE is excellent. The distribution of the different infrared sed types in the L{ir}/L{opt} versus L{ir} plane, where L{ir} and L{opt} are the infrared and optical bolometric luminosities, is discussed. Source-counts at 24, 70 and 160 mu are discussed, and luminosity functions at 3.6 and 24 mu are presented.

We present the reduction of the 24 micron data obtained during the first cosmological survey performed by the Spitzer Space Telescope (First Look Survey, FLS). The survey consists of a shallow observation of 2.5x2 sq deg centered at 17h18m +59d30m (main survey) and a deeper observation of 1x0.5 sq deg centered at 17h17m +59d45m(verification survey). Issues with the reduction of the 24 micron MIPS data are discussed and solutions to attenuate instrumental effects are proposed and applied to the data. Approximately 17000 sources are extracted with a SNR greater than five. The photometry of the point sources is evaluated through PSF fitting using an empirical PSF derived from the data. Aperture corrections and the absolute calibration have been checked using stars in the field. Astrometric and photometric errors depend on the SNR of the source varying between 0.35-1 arcsec and 5-15%, respectively, for sources detected at 20-5 sigma. The flux of the 123 extended sources have been estimated through aperture photometry. The extended sources cover less than 0.3% of the total area of the survey. Based on simulations, the main and verification surveys are 50% complete at 0.3 and 0.15 mJy, respectively. Counterparts have been searched for in optical and radio catalogs. More than 80% of the 24 micron sources have a reliable optical counterpart down to R=25.5. 16% of the sources have a 20 cm counterpart down to 0.1 mJy and ~ 80% of the radio-infrared associations have a reliable optical counterpart. A residual map is obtained by subtracting point sources detected at the 3-sigma level and interpolating the regions occupied by extended sources. Several galactic clouds with low and intermediate velocities are identified by comparison with neutral Hydrogen data from this field.

This work focuses on the properties of dusty tori in active galactic nuclei (AGN) derived from the comparison of SDSS type 1 quasars with mid-Infrared (MIR) counterparts and a new, detailed torus model. The infrared data were taken by the Spitzer Wide-area InfraRed Extragalactic (SWIRE) Survey. Basic model parameters are constraint, such as the density law of the graphite and silicate grains, the torus size and its opening angle. A whole variety of optical depths is supported. The favoured models are those with decreasing density with distance from the centre, while there is no clear tendency as to the covering factor, ie small, medium and large covering factors are almost equally distributed. Based on the models that better describe the observed SEDs, properties such as the accretion luminosity, the mass of dust, the inner to outer radius ratio and the hydrogen column density are computed.

Using the large multi-wavelength data set in the chandra/SWIRE Survey (0.6 square degrees in the Lockman Hole), we show evidence for the existence of highly obscured (Compton-thick) AGN, estimate a lower limit to their surface density and characterize their multi-wavelength properties. Two independent selection methods based on the X-ray and infrared spectral properties are presented. The two selected samples contain 1) 5 X-ray sources with hard X-ray spectra and column densities > 10^24 cm-2, and 2) 120 infrared sources with red and AGN-dominated infrared spectral energy distributions (SEDs). We estimate a surface density of at least 25 Compton-thick AGN per square degree detected in the infrared in the chandra/SWIRE field of which ~40% show distinct AGN signatures in their optical/near-infrared SEDs, the remainings being dominated by the host-galaxy emission. Only ~33% of all Compton-thick AGN are detected in the X-rays at our depth (F(0.3-8 keV)>10^-15 erg/cm2/s. We report the discovery of two sources in our sample of Compton-thick AGN, SWIRE_J104409.95+585224.8 (z=2.54) and SWIRE_J104406.30+583954.1 (z=2.43), which are the most luminous Compton-thick AGN at high-z currently known. The properties of these two sources are discussed in detail with an analysis of their spectra, SEDs, luminosities and black-hole masses.

The ESO-Spitzer extragalactic Imaging Survey (ESIS) is the optical follow up of the Spitzer Wide-Area InfraRed Extragalactic (SWIRE) survey in the ELAIS-S1 area. This paper presents B, V, R Wide Field Imager observations of the first 1.5 square degree of the ESIS survey. Data reduction is described including astrometric calibration, illumination and color corrections, completeness and photometric accuracy estimates. Number counts and color distributions are compared to literature observational and theoretical data, including non-evolutionary, PLE, evolutionary and semi-analytic Lambda-CDM galaxy models, as well as Milky Way stellar predictions. ESIS data are in good agreement with previous works and are best reproduced by evolutionary and hierarchical Lambda-CDM scenarios. The ELAIS-S1 area benefits from extensive follow-up from X-ray to radio frequencies: some potential uses of the multi-wavelength observations are illustrated. Optical-Spitzer color-color plots promise to be very powerful tools to disentangle different classes of sources (e.g. AGNs, starbursts, quiescent galaxies). Ultraviolet GALEX data are matched to optical and Spitzer samples, leading to a discussion of galaxy properties in the UV-to-24 microns color space. The spectral energy distribution of a few objects, from the X-rays to the far-IR are presented as examples of the multi-wavelength study of galaxy emission components in different spectral domains.

We introduce an information-theoretic measure, the Renyi information, to describe the galaxy distribution in space. We discuss properties of the information measure, and demonstrate its relationship with the probability distribution function and multifractal descriptions. Using the First Look Survey galaxy samples observed by the Infrared Array Camera onboard Spitzer Space Telescope, we present measurements of the Renyi information, as well as the counts-in-cells distribution and multifractal properties of galaxies in mid-infrared wavelengths. Guided by multiplicative cascade simulation based on a binomial model, we verify our measurements, and discuss the spatial selection effects on measuring information of the spatial structures. We derive structure scan functions at scales where selection effects are small for the Spitzer samples. We discuss the results, and the potential of applying the Renyi information to measuring other spatial structures.

We address the question of how to deal with confusion limited surveys in the mid-infrared domain by using informations from higher frequency observations over the same sky regions. Such informations, once applied to apparently extended mid-infrared sources, which are indeed ``blends'' of two or more different sources, allow us to disentangle the single counterparts and to split the measured flux density into different components. We present the application of this method to the 24 micron Spitzer archival data in the GOODS-EN1 test field, where apparently extended, ``blended'' sources constitute about 20% of a reliable sample of 983 sources detected above the 5-sigma threshold down to 23 microJy. As higher frequency data-set we have considered the public IRAC images and catalogues on the same field. We show that the 24 micron sample is almost unbiased down to ~40 microJy and the careful application of the deblending procedure does not require any statistical completeness correction (at least at the flux level considered). This is probed by direct comparison of our results with those of Chary et al. (2004), who analysed the same data-set through extensive Monte Carlo simulations. The deblending procedure reduces of about 30% the confusion limit of the MIPS 24 micron survey, allowing one to obtain reliable source counts down to ~40 microJy. The extrapolation of the source counts down to fainter fluxes suggests that our 24 micron sample is able to resolve ~62% of the cosmic background down to a flux level of 38 microJy.

We present ground-based optical and Spitzer infrared imaging observations of the interacting galaxy UGC 10214, the "Tadpole Galaxy" (z = 0.0310), focusing on the star formation activity in the nuclear, disk, spiral arms and tidal tail regions. The major findings of this study are that the Tadpole is actively forming stars in the main disk outside of the nucleus and in the tidal plume, with an estimated mean star formation rate of ~2 to 4 M_sun/yr. The most prominent sites of mid-infrared emission define a "ring" morphology that, combined with the overall morphology of the system, suggest the interaction may belong to the rare class of off-center collisional ring systems that form both shock-induced rings of star formation and tidal plumes. The nuclear emission is solely powered by older stars, with little evidence for ongoing star formation at the center of the Tadpole. Extra-nuclear star formation accounts for >50% of the total star formation in the disk and spiral arms, featuring infrared-bright 'hot spots' that exhibit strong PAH emission, whose band strength is comparable to that of late-type star-forming disk galaxies. The tidal tail, which extends 2 arcmin (~75 kpc) into the intergalactic medium, is populated by super massive star clusters likely triggered by the galaxy-galaxy interaction that has distorted UGC 10214 into its current "tadpole" shape.

We present Spitzer 70um and 160um observations of the Spitzer extragalactic First Look Survey (xFLS). The data reduction techniques and the methods for producing co-added mosaics and source catalogs are discussed. Currently, 26% of the 70um sample and 49% of the 160um-selected sources have redshifts. The majority of sources with redshifts are star-forming galaxies at z<0.5, while about 5% have infrared colors consistent with AGN. The observed infrared colors agree with the spectral energy distribution (SEDs) of local galaxies previously determined from IRAS and ISO data. The average 160um/70um color temperature for the dust is Td~= 30+/-5 K, and the average 70um/24um spectral index is alpha~= 2.4+/-0.4. The observed infrared to radio correlation varies with redshift as expected out to z~1 based on the SEDs of local galaxies. The xFLS number counts at 70um and 160um are consistent within uncertainties with the models of galaxy evolution, but there are indications that the current models may require slight modifications. Deeper 70um observations are needed to constrain the models, and redshifts for the faint sources are required to measure the evolution of the infrared luminosity function.

We present results of a morphological analysis of a small subset of the Spitzer Wide-area InfraRed Extragalactic survey (SWIRE) galaxy population. The analysis is based on public ACS data taken inside the SWIRE N1 field, which are the deepest optical high-resolution imaging available within the SWIRE fields as of today. Our reference sample includes 156 galaxies detected by both ACS and SWIRE. Among the various galaxy morphologies, we disentangle two main classes, spheroids (or bulge-dominated galaxies) and disk-dominated ones, for which we compute the number counts as a function of flux. We then limit our sample to objects with IRAC fluxes brighter than 10 microJy, estimated ~90% completeness limit of the SWIRE catalogues, and compare the observed counts to model predictions. We find that the observed counts of the spheroidal population agree with the expectations of a hierarchical model while a monolithic scenario predicts steeper counts. Both scenaria, however, under-predict the number of late-type galaxies. These observations show that the large majority (close to 80 per cent) of the 3.6 and 4.5 micron galaxy population, even at these moderately faint fluxes, is dominated by spiral and irregular galaxies or mergers.

We present Infrared Array Camera (IRAC) data and source catalogs from the Spitzer Space Telescope Extragalactic First Look Survey. The data were taken in four broad bands centered at nominal wavelengths of 3.6, 4.5, 5.8 and 8.0 microns. A set of mosaics and catalogs have been produced which are ~80% complete and ~99% reliable to their chosen flux density limits. The main field survey covers 3.8 deg^2, and has flux density limits of 20muJy, 25muJy, 100muJy and 100muJy at wavelengths of 3.6, 4.5, 5.8 and 8.0 microns,respectively. The deeper ``verification'' survey covers 0.25 deg^2 with limits of 10muJy, 10muJy, 30muJy and 30muJy, respectively. We also include deep data in the ELAIS-N1 field which covers 0.041deg^2 with limits of 4muJy, 3muJy, 10muJy and 10muJy, respectively, but with only two wavelength coverage at a given sky position. The final bandmerged catalogs contain 103193 objects in the main field, 12224 in the verification field and 5239 in ELAIS-N1. Flux densities of high signal-to-noise objects are accurate to about 10%, and the residual systematic error in the absolute flux density scale is ~2-3%. We have successfully extracted sources at source densities as high as 100000 deg^-2 in our deepest 3.6 and 4.5 micron data. The mosaics and source catalogs will be made available through the Spitzer Science Center archive and the Infrared Science Archive.

We present the results of using three methods, the probability distribution function, the Renyi information, and the multi-fractals, to measure the high moments of galaxy angular distribution in the Infrared Array Camera data of the Spitzer First Look Survey. This is the first time that Renyi information is explicitly used to describe the large-scale galaxy spatial distribution, and our galaxy samples are also the first of the kind in their wavelengths and sensitivities.

We discuss optical associations, spectral energy distributions and photometric redshifts for SWIRE sources in the ELAIS-N1 area and the Lockman Validation Field. The band-merged IRAC (3.6, 4.5, 5.8 and 8.0 mu) and MIPS (24, 70, 160 mu) data have been associated with optical UgriZ data from the INT Wide Field Survey in ELAIS-N1, and with our own optical Ugri data in Lockman-VF. The spectral energy distributions of selected ELAIS sources in N1 detected by SWIRE, most with spectroscopic redshifts, are modelled in terms of a simple set of galaxy and quasar templates in the optical and near infrared, and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst, and AGN dust torus) in the mid infrared. The optical data, together with the IRAC 3.6 and 4.5 mu data, have been used to determine photometric redshifts. For galaxies with known spectroscopic redshifts there is a notable improvement in the photometric redshift when the IRAC data are used, with a reduction in the rms scatter from 10% in (1+z) to 7%. The photometric redshifts are used to derive the 3.6 and 24 mu redshift distribution and to compare this with the predictions of models. For those sources with a clear mid infrared excess, relative to the galaxy starlight model used for the optical and near infrared, the mid and far infrared data are modelled in terms of the same dust emission templates. The proportions found of each template type are: cirrus 31%, M82 29%, Arp 220 10%, AGN dust tori 29%. The distribution of the different infrared sed types in the L_{ir}/L_{opt} versus L_{ir} plane, where L_{ir} and L_{opt} are the infrared and optical bolometric luminosities, is discussed.

The 24 micron array on board the Spitzer Space Telescope is one of three arrays in the Multi-band Imaging Photometer for Spitzer (MIPS) instrument. It provides 5.3 x 5.3 arcmin images at a scale of ~2.5 arcsec per pixel corresponding to sampling of the point spread function which is slightly better than critical (~0.4\lambda/D). A scan-mirror allows dithering of images on the array without the overhead of moving and stabilizing the spacecraft. It also enables efficient mapping of large areas of sky without significant compromise in sensitivity. We present an overview of the pipeline flow and reduction steps involved in the processing of image data acquired with the 24 micron array. Residual instrumental signatures not yet removed in automated processing and strategies for hands-on mitigation thereof are also given.

We present a mid-infrared analysis of 35 quasars with spectroscopic redshifts selected from the Spitzer Wide-area InfraRed Extragalactic Survey (SWIRE). We discuss their optical and mid-infrared (MIR) colors, and show that these quasars occupy well defined regions in MIR color-color space. We examine the issue of type-I AGN candidate selection in detail and propose new selection methods based on mid-IR colors. The available multi-band data allows us to construct two new, well-sampled quasar templates, covering wavelengths from the ultraviolet to the MIR.

We present the first results of the angular auto-correlation function of the galaxies detected by the Infrared Array Camera (IRAC) instrument in the First Look Survey (FLS) of the Spitzer Space Telescope. We detect significant signals of galaxy clustering within the survey area. The angular auto-correlation function of the galaxies detected in each of the four IRAC instrument channels is consistent with a power-law form $w(\theta)=A\theta^{1-\gamma}$ out to $\theta = 0.2\arcdeg$, with the slope ranging from $\gamma = 1.5$ to 1.8. We estimate the correlation amplitudes $A$ to be $2.95 \times 10^{-3}$, $2.03 \times 10^{-3}$, $4.53 \times 10^{-3}$, and $2.34 \times 10^{-3}$ at $\theta=1\arcdeg$ for galaxies detected in the IRAC 3.6$\mu$m, 4.5$\mu$m, 5.8$\mu$m, and 8.0$\mu$m instrument channels, respectively. We compare our measurements at 3.6$\mu$m with the previous K-band measurements, and discuss the implications of these results.

We report on submillimeter and infrared observations of 28 radio-selected galaxies in the Spitzer First Look Survey Verification field (FLSV). All of the radio-selected galaxies that show evidence for emission at 850um with SCUBA have Spitzer counterparts at 24um, while only half of the radio-selected galaxies without 850um emission have detectable counterparts at 24um. The data show a wide range of infrared colors (S70/S24 < 5--30, S8/S3.6 < 0.3--4), indicative of a mixture of infrared-warm AGN and cooler starburst dominated sources. The galaxies showing 850um emission have Spitzer flux densities and flux density ratios consistent with the range of values expected for high-redshift (z=1--4) ultraluminous infrared galaxies.

We characterize the SWIRE galaxy populations in the SWIRE validation field within the Lockman Hole, based on the 3.6-24$\mu$ Spitzer data and deep U,g',r',r' optical imaging within an area ~1/3 sq. deg for ~16,000 Spitzer-SWIRE sources. The entire SWIRE survey will discover over 2.3 million galaxies at 3.6$\mu$m and almost 350,000 at 24$\mu$m; ~70,000 of these will be 5-band 3.6-24$\mu$ detections. The colors cover a broad range, generally well represented by redshifted spectral energy distributions of known galaxy populations, however significant samples of unusually blue objects in the [3.6-4.5]$\mu$m color are found, as well as many objects very red in the 3.6-24$\mu$m mid-IR. Nine of these are investigated and are interpreted as star-forming systems, starbursts and AGN from z=0.37 to 2.8, with luminosities from L$_{IR}$=10$^{10.3}$ to 10$^{13.7}$ L$_{\odot}$

Using the SPITZER SPACE TELESCOPE and the VLA, we present the first DIRECT evidence that the well-known far-infrared/radio correlation is valid to cosmologically significant redshift. We also confirm, with improved statistics compared with previous surveys, a similar result for the Mid-IR/radio correlation. We explore the dependence of monochromatic q_{24} and q_{70} on z. The results were obtained by matching Spitzer sources at 24 and 70 microns with VLA 1.4 GHz micro-Jy radio sources obtained for the SPITZER FLS. Spectroscopic redshifts have been obtained for over 500 matched IR/radio sources using observations at WIYN, Keck and archival SDSS data extending out to z > 2. We find that q_{24} shows significantly more dispersion than q_{70}. By comparing the observed fluxes at 70, 24 and 4.5 microns with a library of SED templates, we find that the larger dispersion in q_{24} is predictable in terms of systematic variations in SED shape throughout the population. Although the models are not able to encompass the full range of observed behavior (both the presence of either extremely flat or extremely steep IR SEDs), the fitting parameters were used to `k-correct' the higher-z galaxies which resulted in a reduced scatter in q. For comparison, we also corrected these data using the SED for M82. The results for 24 and 70 microns provide strong consistent evidence for the universality of the mid-IR/radio and far-IR/radio correlations out to redshifts of at least z = 1.

We present the Spitzer MIPS 24 micron source counts in the Extragalactic First Look Survey main, verification and ELAIS-N1 fields. Spitzer's increased sensitivity and efficiency in large areal coverage over previous infrared telescopes, coupled with the enhanced sensitivity of the 24 micron band to sources at intermediate redshift, dramatically improve the quality and statistics of number counts in the mid-infrared. The First Look Survey observations cover areas of, respectively, 4.4, 0.26 and 0.015 sq.deg. and reach 3-sigma depths of 0.11, 0.08 and 0.03 mJy. The extragalactic counts derived for each survey agree remarkably well. The counts can be fitted by a super-Euclidean power law of index alpha=-2.9 from 0.2 to 0.9 mJy, with a flattening of the counts at fluxes fainter than 0.2 mJy. Comparison with infrared galaxy evolution models reveals a peak's displacement in the 24 micron counts. This is probably due to the detection of a new population of galaxies with redshift between 1 and 2, previously unseen in the 15 micron deep counts.

Selection of active galactic nuclei (AGN) in the infrared allows the discovery of AGN whose optical emission is extinguished by dust. In this paper, we use the Spitzer Space Telescope First Look Survey (FLS) to assess what fraction of AGN with mid-infrared luminosities comparable to quasars are missed in optical quasar surveys due to dust obscuration. We begin by using the Sloan Digital Sky Survey (SDSS) database to identify 54 quasars within the 4 deg^2 extragalactic FLS. These quasars occupy a distinct region in mid-infrared color space by virtue of their strong, red, continua. This has allowed us to define a mid-infrared color criterion for selecting AGN candidates. About 2000 FLS objects have colors consistent with them being AGN, but most are much fainter in the mid-infrared than the SDSS quasars, which typically have 8 micron flux densities, S(8.0), ~1 mJy. We have investigated the properties of the 43 objects with S(8.0) >= 1 mJy satisfying our AGN color selection. This sample should contain both unobscured quasars, and AGN which are absent from the SDSS survey due to extinction in the optical. After removing 16 known quasars, three probable normal quasars, and eight spurious or confused objects from the initial sample of 43, we are left with 16 objects which are likely to be obscured quasars or luminous Seyfert-2 galaxies. This suggests the numbers of obscured and unobscured AGN are similar in samples selected in the mid-infrared at S(8.0)~1 mJy.

We present the first analysis of large-scale clustering from the Spitzer Wide-area InfraRed Extragalactic legacy survey: SWIRE. We compute the angular correlation function of galaxies selected to have 3.6 micron fluxes brighter than 32 muJy in three fields totaling two square degrees in area. In each field we detect clustering with a high level of significance. The amplitude and slope of the correlation function is consistent between the three fields and is modeled as $w(theta)=A\theta^{1-\gamma}$ with $A=(0.6\pm0.3)\times10^{-3}, \gamma=2.03\pm0.10$. With a fixed slope of $\gamma=1.8$, we obtain an amplitude of $A=(1.7\pm0.1)\times10^{-3}$. Assuming an equivalent depth of K~18.7 mag we find our errors are smaller but our results are consistent with existing clustering measurements in K-band surveys and with stable clustering models. We estimate our median redshift $z\simeq 0.75$ and this allows us to obtain an estimate of the three-dimensional correlation function $\xi(r)$, for which we find $r_0=4.4\pm0.1 h^{-1}$ Mpc.

The SIRTF Wide-area InfraRed Extragalactic survey (SWIRE), the largest SIRTF Legacy program, is a wide-area, imaging survey to trace the evolution of dusty, star-forming galaxies, evolved stellar populations, and AGN as a function of environment, from redshifts z~3 to the current epoch. SWIRE will survey 7 high-latitude fields, totaling 60 - 65 sq. deg. in all 7 SIRTF bands: IRAC 3.6, 4.5, 5.6, 8 microns and MIPS 24, 70, 160 microns. The Legacy Extragalactic Catalog may contain in excess of 2 million IR-selected galaxies, dominated by (1) ~150,000 luminous infrared galaxies (LIRGs: L{FIR}>10^11 L_sun), ~7000 of these with z>2; (2) 1 million early-type galaxies, ~10,000 with z>2; and (3) \~20,000 classical AGN, plus significantly more dust-obscured QSO/AGN among the LIRGs. SWIRE will provide an unprecedented view of the evolution of galaxies, structure, and AGN. The key scientific goals of SWIRE are: (1) to determine the evolution of actively star-forming and passively evolving galaxies in order to understand the history of galaxy formation in the context of cosmic structure formation; (2) to determine the evolution of the spatial distribution and clustering of evolved galaxies, starbursts and AGN in the key redshift range, 0.5<z<3, over which much of cosmic evolution has occurred; (3) to determine the evolutionary relationship between ``normal galaxies'' and AGN, and the contribution of AGN accretion energy vs stellar nucleosynthesis to the cosmic backgrounds. SWIRE's large area is important to create statistically significant population samples over enough volume cells that we can resolve the star formation history as a function of epoch and environment The large volume is also optimised for finding rare objects.

The Very Large Array (VLA) has been used in C-configuration to map an area \~0.3 square degrees at 1.4 GHz with 5sigma sensitivities of 0.305, 0.325, 0.380 and 0.450 mJy/beam over four equal subareas. Radio properties are presented for 62 detected sources. Deep optical imaging to Gunn r~25 mag using the Hale 5-m telescope covering ~0.21 square degrees is reported for a subset of 43 sources. This optical follow-up is much deeper than that of existing larger area radio surveys of similar radio sensitivity. Archival J, H and K-band photometry from the Two-Micron All Sky Survey (2MASS) is also presented. Using a robust likelihood-ratio technique, we optically identified 26 radio sources with probabilities >80%. Nine are uncertain/ambiguous detections and eight are empty fields. Comparisons with a stellar synthesis model that includes radio emission and dust reddening suggests that the near-infrared--to--optical emission in a small, bright sub-sample is reddened by `optically-thin' dust with absorption A_V ~ 2-2.5 mag, regardless of morphological type. This is consistent with other, more direct determinations of absorption. The radio--to--optical(--near-infrared) flux ratios of early-type galaxies require significant contamination in the radio by an active galactic nucleus (AGN), consistent with the current paradigm. Using our simple modeling approach, we also discuss a potential diagnostic for selecting Ultraluminous Infrared Galaxies (ULIGS) to z~1.6 from microJansky radio surveys.

In this work we incorporate the newest ISO results on the mid-infrared spectral-energy-distributions (MIR SEDs) of galaxies into models for the number counts and redshift distributions of MIR surveys. A three-component model, with empirically determined MIR SED templates of (1) a cirrus/PDR component (2) a starburst component and (3) an AGN component, is developed for infrared (3--120\micron) SEDs of galaxies. The model includes a complete IRAS 25\micron selected sample of 1406 local galaxies ($z \leq 0.1$; Shupe et al. 1998a). Results based on these 1406 spectra show that the MIR emission features cause significant effects on the redshift dependence of the K-corrections for fluxes in the WIRE 25\micron band and ISOCAM 15\micron band. This in turn will affect deep counts and redshift distributions in these two bands, as shown by the predictions of two evolution models (a luminosity evolution model with $L\propto (1+z)^3$ and a density evolution model with $\rho\propto (1+z)^4$). The dips-and-bumps on curves of MIR number counts, caused by the emission features, should be useful indicators of evolution mode. The strong emission features at $\sim 6$--8\micron will help the detections of relatively high redshift ($z\sim 2$) galaxies in MIR surveys. On the other hand, determinations of the evolutionary rate based on the slope of source counts, and studies on the large scale structures using the redshift distribution of MIR sources, will have to treat the effects of the MIR emission features carefully. We have also estimated a 15\micron local luminosity function from the predicted 15\micron fluxes of the 1406 galaxies using the bivariate (15\micron vs. 25\micron luminosities) method. This luminosity function will improve our understanding of the ISOCAM 15\micron surveys.

We have established a model to systematically estimate the contribution of the mid-infrared emission features between 3 and 11.6 micron to the IRAS in-band fluxes, using the results of ISO PHT-S observation of 16 galaxies by Lu et al. (1997). The model is used to estimate more properly the k-corrections for calculating the restframe 12 and 25 micron fluxes and luminosities of IRAS galaxies. We have studied the 12-25 micron color-luminosity relation for a sample of galaxies selected at 25 microns. The color is found to correlate well with the 25 micron luminosity, the far-infrared and the blue luminosities. The relations with the mid-infrared luminosities are more sensitive to different populations of galaxies, while a single relation of the 12-25 micron color vs. the ratio of the far-infrared and the blue luminosities applies equally well to these different populations. The luminous and ultraluminous infrared galaxies have redder 12-25 micron colors than those of the quasars. These relations provide powerful tools to differentiate different populations of galaxies. We have also selected a sample of 668 galaxies from the IRAS Faint Source Survey flux-density limited at 200 mJy at 12 microns. A 12 micron local luminosity function is derived and, for the first time in the literature, effects of density variation in the local universe are considered and corrected in the calculation of the 12 micron luminosity function. It is also found that the 12 micron-selected sample are dominated by quasars and active galaxies, which therefore strongly affect the 12 micron luminosity function at high luminosities. The ultraluminous infrared galaxies are relatively rare at 12 micron comparing with a 25 micron sample.

The local luminosity function at 25 $\mu$m provides the basis for interpreting the results of deep mid-infrared surveys planned or in progress with space astrophysics missions including ISO, WIRE and SIRTF. We have selected a sample of 1458 galaxies from the IRAS Faint Source Survey with a flux density limit of 250 mJy at 25 $\mu$m. The local luminosity function is derived using both parametric and non-parametric maximum-likelihood techniques, and the classical $1/V_{max}$ estimator. Comparison of these results shows that the $1/V_{max}$ estimate of the luminosity function is significantly affected by the Local Supercluster. A maximum-likelihood fit to the radial density shows no systematic increase that would be caused by density evolution of the galaxy population. The density fit is used to correct the $1/V_{max}$ estimate. We also demonstrate the high quality and completeness of our sample by a variety of methods. The luminosity function derived from this sample is compared to previously published estimates, showing the prior estimates to have been strongly affected by the Local Supercluster. Our new luminosity function leads to lower estimates of mid-infrared backgrounds and number counts.