We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys.

The COSMOS-Legacy survey is a 4.6 Ms Chandra program that has imaged 2.2 deg$^2$ of the COSMOS field with an effective exposure of $\simeq$160 ks over the central 1.5 deg$^2$ and of $\simeq$80 ks in the remaining area. The survey is the combination of 56 new observations, obtained as an X-ray Visionary Project, with the previous C-COSMOS survey. We describe the reduction and analysis of the new observations and the properties of 2273 point sources detected above a spurious probability of 2$\times 10^{-5}$. We also present the updated properties of the C-COSMOS sources detected in the new data. The whole survey includes 4016 point sources (3814, 2920 and 2440 in the full, soft and hard band). The limiting depths are 2.2 $\times$ 10$^{-16}$, 1.5 $\times$ 10$^{-15}$ and 8.9$\times$ 10$^{-16}$ ${\rm erg~cm}^{-2}~{\rm s}^{-1}$ in the 0.5-2, 2-10 and 0.5-10 keV bands, respectively. The observed fraction of obscured AGN with column density $> 10^{22}$ cm$^{-2}$ from the hardness ratio (HR) is $\sim$50$^{+17}_{-16}$%. Given the large sample, we compute source number counts in the hard and soft bands, significantly reducing the uncertainties of 5-10%. For the first time, we compute number counts for obscured (HR$>$-0.2) and unobscured (HR$<$-0.2) sources and find significant differences between the two populations in the soft band. Due to the un-precedent large exposure, COSMOS-Legacy area is 3 times larger than surveys at similar depth and its depth is 3 times fainter than surveys covering similar area. The area-flux region occupied by COSMOS-Legacy is likely to remain unsurpassed for years to come.

Limited studies have been performed on the radio-loud fraction in X-ray selected type 1 AGN samples. The consistency between various radio-loudness definitions also needs to be checked. We measure the radio-loudness of the 407 type 1 AGNs in the XMM-COSMOS quasar sample using nine criteria from the literature (six defined in the rest-frame and three defined in the observed frame): $R_L=\log(L_{5GHz}/L_B)$, $q_{24}=\log(L_{24\mu m}/L_{1.4GHz})$, $R_{uv}=\log(L_{5GHz}/L_{2500\AA})$, $R_{i}=\log(L_{1.4GHz}/L_i)$, $R_X=\log(\nu L_{\nu}(5GHz)/L_X)$, $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$, $R_{L,obs}=\log(f_{1.4GHz}/f_B)$ (observed frame), $R_{i,obs}=\log(f_{1.4GHz}/f_i)$ (observed frame), and $q_{24, obs}=\log(f_{24\mu m}/f_{1.4GHz})$ (observed frame). Using any single criterion defined in the rest-frame, we find a low radio-loud fraction of $\lesssim 5\%$ in the XMM-COSMOS type 1 AGN sample, except for $R_{uv}$. Requiring that any two criteria agree reduces the radio-loud fraction to $\lesssim 2\%$ for about 3/4 of the cases. The low radio-loud fraction cannot be simply explained by the contribution of the host galaxy luminosity and reddening. The $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$ gives the smallest radio-loud fraction. Two of the three radio-loud fractions from the criteria defined in the observed frame without k-correction ($R_{L,obs}$ and $R_{i,obs}$) are much larger than the radio-loud fractions from other criteria.

The mid-infrared to ultraviolet (0.1 -- 10 $\mu m$) spectral energy distribution (SED) shapes of 407 X-ray-selected radio-quiet type 1 AGN in the wide-field ``Cosmic Evolution Survey" (COSMOS) have been studied for signs of evolution. For a sub-sample of 200 radio-quiet quasars with black hole mass estimates and host galaxy corrections, we studied their mean SEDs as a function of a broad range of redshift, bolometric luminosity, black hole mass and Eddington ratio, and compared them with the Elvis et al. (1994, E94) type 1 AGN mean SED. We found that the mean SEDs in each bin are closely similar to each other, showing no statistical significant evidence of dependence on any of the analyzed parameters. We also measured the SED dispersion as a function of these four parameters, and found no significant dependencies. The dispersion of the XMM-COSMOS SEDs is generally larger than E94 SED dispersion in the ultraviolet, which might be due to the broader ``window function'' for COSMOS quasars, and their X-ray based selection.

We define a quasar-galaxy mixing diagram using the slopes of their spectral energy distributions (SEDs) from 1\mu m to 3000\AA\ and from 1\mu m to 3\mu m in the rest frame. The mixing diagram can easily distinguish among quasar-dominated, galaxy-dominated and reddening-dominated SED shapes. By studying the position of the 413 XMM selected Type 1 AGN in the wide-field "Cosmic Evolution Survey" (COSMOS) in the mixing diagram, we find that a combination of the Elvis et al. (1994, hereafter E94) quasar SED with various contributions from galaxy emission and some dust reddening is remarkably effective in describing the SED shape from 0.3-3\mu m for large ranges of redshift, luminosity, black hole mass and Eddington ratio of type 1 AGN. In particular, the location in the mixing diagram of the highest luminosity AGN is very close (within 1\sigma) to that of the E94 SED. The mixing diagram can also be used to estimate the host galaxy fraction and reddening in quasar. We also show examples of some outliers which might be AGN in different evolutionary stages compared to the majority of AGN in the quasar-host galaxy co-evolution cycle.

We review our knowledge of the most basic properties of the AGN obscuring region - its location, scale, symmetry, and mean covering factor - and discuss new evidence on the distribution of covering factors in a sample of ~9000 quasars with WISE, UKIDSS, and SDSS photometry. The obscuring regions of AGN may be in some ways more complex than we thought - multi-scale, not symmetric, chaotic - and in some ways simpler - with no dependence on luminosity, and a covering factor distribution that may be determined by the simplest of considerations - e.g. random misalignments.

We investigate the range of covering factors (determined from the ratio of IR to UV/optical luminosity) seen in luminous type 1 quasars using a combination of data from the WISE, UKIDSS and SDSS surveys. Accretion disk (UV/optical) and obscuring dust (IR) luminosities are measured via the use of a simple three component SED model. We use these estimates to investigate the distribution of covering factors and its relationship to both accretion luminosity and IR SED shape. The distribution of covering factors (f_C) is observed to be log-normal, with a bias-corrected mean of <log10 f_C>=-0.41 and standard deviation of 0.2. The fraction of IR luminosity emitted in the near-IR (1--5 micron) is found to be high (~40 per cent), and strongly dependant on covering factor.

The "Cosmic Evolution Survey" (COSMOS) enables the study of the Spectral Energy Distributions (SEDs) of Active Galactic Nuclei (AGN) because of the deep coverage and rich sampling of frequencies from X-ray to radio. Here we present a SED catalog of 413 X-ray (\xmm) selected type 1 (emission line FWHM$>2000$ km s$^{-1}$) AGN with Magellan, SDSS or VLT spectrum. The SEDs are corrected for the Galactic extinction, for broad emission line contributions, constrained variability, and for host galaxy contribution. We present the mean SED and the dispersion SEDs after the above corrections in the rest frame 1.4 GHz to 40 keV, and show examples of the variety of SEDs encountered. In the near-infrared to optical (rest frame $\sim 8\mu m$-- 4000\AA), the photometry is complete for the whole sample and the mean SED is derived from detections only. Reddening and host galaxy contamination could account for a large fraction of the observed SED variety. The SEDs are all available on-line.

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 0.5-2 keV band, 7.3 10^-16 erg cm^-2 s^-1 in the 2-10 keV band, and 5.7 10^-16 erg cm^-2 s-1 in the 0.5-10 keV band. In this paper we report the i, K and 3.6micron 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. Making use of the large number of X-ray sources, we update the "classic locus" of AGN and define a new locus containing 90% of the AGN in the survey with full band luminosity >10^42 erg/s. 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.

Close encounters between galaxies are expected to be a viable mechanism, as predicted by numerical simulations, by which accretion onto supermassive black holes can be initiated. To test this scenario, we construct a sample of 562 galaxies (M*>2.5x10^10 M_sun) in kinematic pairs over the redshift range 0.25 < z < 1.05 that are more likely to be interacting than a well-matched control sample of 2726 galaxies not identified as being in a pair, both from the zCOSMOS 20k spectroscopic catalog. Galaxies that harbor an active galactic nucleus (AGN) are identified on the basis of their X-ray emission (L_x>2x10^42 erg s^-1) detected by Chandra. We find a higher fraction of AGN in galaxies in pairs relative to isolated galaxies of similar stellar mass. Our result is primarily due to an enhancement of AGN activity, by a factor of 1.9 (observed) and 2.6 (intrinsic), for galaxies in pairs of projected separation less than 75 kpc and line-of-sight velocity offset less than 500 km s^-1. This study demonstrates that close kinematic pairs are conducive environments for black hole growth either indicating a causal physical connection or an inherent relation, such as, to enhanced star formation. In the Appendix, we describe a method to estimate the intrinsic fractions of galaxies (either in pairs or the field) hosting an AGN with confidence intervals, and an excess fraction in pairs. We estimate that 17.8_{-7.4}^{+8.4}% of all moderate-luminosity AGN activity takes place within galaxies undergoing early stages of interaction that leaves open the question as to what physical processes are responsible for fueling the remaining ~80% that may include late-stage mergers.

We present the high redshift (3<z<5.3) 0.5-2 keV number counts and the 2-10 keV (rest frame) space density of X-ray selected AGNs detected in the Chandra COSMOS survey. The sample comprises 81 X-ray detected sources with available spectroscopic (31) and photometric (50) redshifts plus 20 sources with a formal z_phot<3 but with a broad photometric redshift probability distribution, such that z_phot+1sigma>3. 81 sources are selected in the 0.5-2 keV band, 14 are selected in the 2-10 keV and 6 in the 0.5-10 keV bands. We sample the high luminosity (log L_(2-10 keV)>44.15 erg/s) space density up to z~5 and a fainter luminosity range (43.5<log L(2-10 keV)<44.15 erg/s) than previous studies, up to z=3.5. We weighted the contribution to the number counts and the space density of the sources with photometric redshift by using their probability of being at z>3. We find that the space density of high-luminosity AGNs declines exponentially at all the redshifts, confirming the trend observed for optically selected quasars. At lower luminosity, the measured space density is not conclusive, and a larger sample of faint sources is needed. Comparisons with optical luminosity functions and BH formation models are presented together with prospects for future surveys.

We show that the Hot-Dust-Poor (HDP) quasars, originally found in the X-ray selected XMM-COSMOS type 1 AGN sample, are just as common in two samples selected at optical/infrared wavelengths: the Richards et al. Spitzer/SDSS sample ($8.7%\pm2.2%$), and the PG-quasar dominated sample of Elvis et al. ($9.5%\pm 5.0%$). The properties of the HDP quasars in these two samples are consistent with the XMM-COSMOS sample, except that, at the $99% (\sim 2.5\sigma)$ significance, a larger proportion of the HDP quasars in the Spitzer/SDSS sample have weak host galaxy contributions, probably due to the selection criteria used. Either the host-dust is destroyed (dynamically or by radiation), or is offset from the central black hole due to recoiling. Alternatively, the universality of HDP quasars in samples with different selection methods and the continuous distribution of dust covering factor in type 1 AGNs, suggest that the range of SEDs could be related to the range of tilts in warped fueling disks, as in the model of Lawrence and Elvis (2010), with HDP quasars having relatively small warps.

We show how accretion rate governs the physical properties of a sample of unobscured broad-line, narrow-line, and lineless active galactic nuclei (AGNs). We avoid the systematic errors plaguing previous studies of AGN accretion rate by using accurate accretion luminosities (L_int) from well-sampled multiwavelength SEDs from the Cosmic Evolution Survey (COSMOS), and accurate black hole masses derived from virial scaling relations (for broad-line AGNs) or host-AGN relations (for narrow-line and lineless AGNs). In general, broad emission lines are present only at the highest accretion rates (L_int/L_Edd > 0.01), and these rapidly accreting AGNs are observed as broad-line AGNs or possibly as obscured narrow-line AGNs. Narrow-line and lineless AGNs at lower specific accretion rates (L_int/L_Edd < 0.01) are unobscured and yet lack a broad line region. The disappearance of the broad emission lines is caused by an expanding radiatively inefficient accretion flow (RIAF) at the inner radius of the accretion disk. The presence of the RIAF also drives L_int/L_Edd < 10^-2 narrow-line and lineless AGNs to 10 times higher ratios of radio to optical/UV emission than L_int/L_Edd > 0.01 broad-line AGNs, since the unbound nature of the RIAF means it is easier to form a radio outflow. The IR torus signature also tends to become weaker or disappear from L_int/L_Edd < 0.01 AGNs, although there may be additional mid-IR synchrotron emission associated with the RIAF. Together these results suggest that specific accretion rate is an important physical "axis" of AGN unification, described by a simple model.

We report a sizable class of type 1 active galactic nuclei (AGNs) with unusually weak near-infrared (1-3{\mu}m) emission in the XMM-COSMOS type 1 AGN sample. The fraction of these "hot-dust-poor" AGNs increases with redshift from 6% at lowredshift (z < 2) to 20% at moderate high redshift (2 < z < 3.5). There is no clear trend of the fraction with other parameters: bolometric luminosity, Eddington ratio, black hole mass, and X-ray luminosity. The 3{\mu}m emission relative to the 1{\mu}m emission is a factor of 2-4 smaller than the typical Elvis et al. AGN spectral energy distribution (SED), which indicates a "torus" covering factor of 2%-29%, a factor of 3-40 smaller than required by unified models. The weak hot dust emission seems to expose an extension of the accretion disk continuum in some of the source SEDs. We estimate the outer edge of their accretion disks to lie at $(0.3-2.0) /times 10^4$ Schwarzschild radii, ~10-23 times the gravitational stability radii. Formation scenarios for these sources are discussed.

We report the final optical identifications of the medium-depth (~60 ksec), contiguous (2 deg^2) XMM-Newton survey of the COSMOS field. XMM-Newton has detected ~800 X-ray sources down to limiting fluxes of ~5x10^{-16}, ~3x10^{-15}, and ~7x10^{-15} erg/cm2/s in the 0.5-2 keV, 2-10 keV and 5-10 keV bands, respectively. The work is complemented by an extensive collection of multi-wavelength data from 24 micron 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 24micron 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 logL_X>44.5 AGN is at z~2. Spectroscopically-identified obscured and unobscured AGN, as well as normal and starforming 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 distribution 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.

We present a detailed study of a peculiar source in the COSMOS survey at z=0.359. Source CXOCJ100043.1+020637 (CID-42) presents two compact optical sources embedded in the same galaxy. The distance between the 2, measured in the HST/ACS image, is 0.495" that, at the redshift of the source, corresponds to a projected separation of 2.46 kpc. A large (~1200 km/s) velocity offset between the narrow and broad components of Hbeta has been measured in three different optical spectra from the VLT/VIMOS and Magellan/IMACS instruments. CID-42 is also the only X-ray source having in its X-ray spectra a strong redshifted broad absorption iron line, and an iron emission line, drawing an inverted P-Cygni profile. The Chandra and XMM data show that the absorption line is variable in energy by 500 eV over 4 years and that the absorber has to be highly ionized, in order not to leave a signature in the soft X-ray spectrum. That these features occur in the same source is unlikely to be a coincidence. We envisage two possible explanations: (1) a gravitational wave recoiling black hole (BH), caught 1-10 Myr after merging, (2) a Type 1/ Type 2 system in the same galaxy where the Type 1 is recoiling due to slingshot effect produced by a triple BH system. The first possibility gives us a candidate gravitational waves recoiling BH with both spectroscopic and imaging signatures. In the second case, the X-ray absorption line can be explained as a BAL-like outflow from the foreground nucleus (a Type 2 AGN) at the rearer one (a Type 1 AGN), which illuminates the otherwise undetectable wind, giving us the first opportunity to show that fast winds are present in obscured AGN.

(Abriged) We report on the measurement of the rest frame K-band luminosity and total stellar mass of the hosts of 89 broad line Active Galactic Nuclei detected in the zCOSMOS survey in the redshift range 1<z<2.2. The unprecedented multiwavelength coverage of the survey field allows us to disentangle the emission of the host galaxy from that of the nuclear black hole in their Spectral Energy Distributions. We derive an estimate of black hole masses through the analysis of the broad Mg II emission lines observed in the medium-resolution spectra taken with VIMOS/VLT as part of the zCOSMOS project. We found that, as compared to the local value, the average black hole to host galaxy mass ratio appears to evolve positively with redshift, with a best fit evolution of the form (1+z)^{0.68 \pm0.12 +0.6 -0.3}, where the large asymmetric systematic errors stem from the uncertainties in the choice of IMF, in the calibration of the virial relation used to estimate BH masses and in the mean QSO SED adopted. A thorough analysis of observational biases induced by intrinsic scatter in the scaling relations reinforces the conclusion that an evolution of the MBH-M* relation must ensue for actively growing black holes at early times: either its overall normalization, or its intrinsic scatter (or both) appear to increase with redshift. This can be interpreted as signature of either a more rapid growth of supermassive black holes at high redshift, a change of structural properties of AGN hosts at earlier times, or a significant mismatch between the typical growth times of nuclear black holes and host galaxies.

We present black hole masses and accretion rates for 182 Type 1 AGN in COSMOS. We estimate masses using the scaling relations for the broad Hb, MgII, and CIV emission lines in the redshift ranges 0.16<z<0.88, 1<z<2.4, and 2.7<z<4.9. We estimate the accretion rate using an Eddington ratio L_I/L_Edd estimated from optical and X-ray data. We find that very few Type 1 AGN accrete below L_I/L_Edd ~ 0.01, despite simulations of synthetic spectra which show that the survey is sensitive to such Type 1 AGN. At lower accretion rates the BLR may become obscured, diluted or nonexistent. We find evidence that Type 1 AGN at higher accretion rates have higher optical luminosities, as more of their emission comes from the cool (optical) accretion disk with respect to shorter wavelengths. We measure a larger range in accretion rate than previous works, suggesting that COSMOS is more efficient at finding low accretion rate Type 1 AGN. However the measured range in accretion rate is still comparable to the intrinsic scatter from the scaling relations, suggesting that Type 1 AGN accrete at a narrow range of Eddington ratio, with L_I/L_Edd ~ 0.1.

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra} program that has imaged the central 0.5 sq.deg of the COSMOS field (centered at 10h, +02deg) with an effective exposure of ~160ksec, and an outer 0.4sq.deg. area with an effective exposure of ~80ksec. The limiting source detection depths are 1.9e-16 erg cm(-2) s(-1) in the Soft (0.5-2 keV) band, 7.3e(-16) erg cm^-2 s^-1 in the Hard (2-10 keV) band, and 5.7e(-16) erg cm(-2) s(-1) in the Full (0.5-10 keV) band. Here we describe the strategy, design and execution of the C-COSMOS survey, and present the catalog of 1761 point sources detected at a probability of being spurious of <2e(-5) (1655 in the Full, 1340 in the Soft, and 1017 in the Hard bands). By using a grid of 36 heavily (~50%) overlapping pointing positions with the ACIS-I imager, a remarkably uniform (to 12%) exposure across the inner 0.5 sq.deg field was obtained, leading to a sharply defined lower flux limit. The widely different PSFs obtained in each exposure at each point in the field required a novel source detection method, because of the overlapping tiling strategy, which is described in a companion paper. (Puccetti et al. Paper II). This method produced reliable sources down to a 7-12 counts, as verified by the resulting logN-logS curve, with sub-arcsecond positions, enabling optical and infrared identifications of virtually all sources, as reported in a second companion paper (Civano et al. Paper III). The full catalog is described here in detail, and is available on-line.

We present optical spectroscopy for an X-ray and optical flux-limited sample of 677 XMM-Newton selected targets covering the 2 deg^2 COSMOS field, with a yield of 485 high-confidence redshifts. The majority of the spectra were obtained over three seasons (2005-2007) with the IMACS instrument on the Magellan (Baade) telescope. We also include in the sample previously published Sloan Digital Sky Survey spectra and supplemental observations with MMT/Hectospec. We detail the observations and classification analyses. The survey is 90% complete to flux limits of f_{0.5-10 keV}>8 x 10^-16 erg cm^-2 s^-1 and i_AB+<22, where over 90% of targets have high-confidence redshifts. Making simple corrections for incompleteness due to redshift and spectral type allows for a description of the complete population to $i_AB+<23. The corrected sample includes 57% broad emission line (Type 1, unobscured) AGN at 0.13<z<4.26, 25% narrow emission line (Type 2, obscured) AGN at 0.07<z<1.29, and 18% absorption line (host-dominated, obscured) AGN at 0<z<1.22 (excluding the stars that made up 4% of the X-ray targets). We show that the survey's limits in X-ray and optical flux include nearly all X-ray AGN (defined by L_{0.5-10 keV}>3 x 10^42 erg s^-1) to z<1, of both optically obscured and unobscured types. We find statistically significant evidence that the obscured to unobscured AGN ratio at z<1 increases with redshift and decreases with luminosity.

We report on the discovery of strongly variable FeII and MgII absorption lines seen at z=1.48 in the spectra of the z=4.05 GRB 060206 obtained between 4.13 to 7.63 hours (observer frame) after the burst. In particular, the FeII line equivalent width (EW) decayed rapidly from 1.72+-0.25 AA to 0.28+-0.21 AA, only to increase to 0.96+-0.21 AA in a later date spectrum. The MgII doublet shows even more complicated evolution: the weaker line of the doublet drops from 2.05+-0.25 AA to 0.92+-0.32 AA, but then more than doubles to 2.47+-0.41 AA in later data. The ratio of the EWs for the MgII doublet is also variable, being closer to 1:1 (saturated regime) when the lines are stronger and becoming closer to 2:1 (unsaturated regime) when the lines are weaker, consistent with expectations based on atomic physics. We have investigated and rejected the possibility of any instrumental or atmospheric effects causing the observed strong variations. Our discovery of clearly variable intervening FeII and MgII lines lends very strong support to their scenario, in which the characteristic size of intervening patches of MgII ``clouds'' is comparable to the GRB beam size, i.e, about 10^16 cm. We discuss various implications of this discovery, including the nature of the MgII absorbers, the physics of GRBs, and measurements of chemical abundances from GRB and quasar absorption lines.

We report on two recent z~4 gamma-ray bursts (GRBs), GRB 060206 and GRB 060210, for which we have obtained well-sampled optical light curves. Our data, combined with early optical data reported in the literature, shows unusual behavior for both afterglows. In R-band GRB 060206 (z=4.045) experienced a slow early decay, followed by a rapid increase in brightness by factor ~2.5 about 1 hour after the burst. Its afterglow then faded in a broken power-law fashion, with a smooth break at t_b=0.6 days, but with additional, less dramatic (~10%) ``bumps and wiggles'', well detected in the densely sampled light curve. The R-band afterglow of GRB 060210 (z=3.91) is also unusual: the light curves was more or less flat between 60 and 300 sec after the burst, followed by ~70% increase at ~600 sec after the burst, after which the light curve declined as a t^{-1.3} power-law. Despite earlier reports to the contrary, we find that for GRB 060206 X-rays also more or less follow the optical decay, but with significant variations on short timescales. However, the X-ray afterglow is contaminated by a nearby, variable source, which especially at late times obscures the behavior of the X-ray afterglow. We argue that ``anomalous'' optical afterglows are likely to be the norm, and that the rapid variations often seen in Swift-XRT data would also be seen in the optical light curves, given good enough sampling. As a result, some of the often employed procedures, such as deriving the jet opening angle using a smooth broken power-law fit to the optical light curves, in many cases might have a poor statistical significance. We argue that the early increase in brighness for both bursts might be due to the turn-on of the external shock. Existence of such features could provide valuable additional information about the burst. (Abridged)

We present early photometric and spectroscopic data on the afterglow of GRB 060218 and report the evolution of the underlying supernova 2006aj. Our data span a time-range of 4 days to 22 days after the GRB and clearly establish that SN 2006aj is a fast-evolving broad-lined Type Ic SN with an extremely short rise-time (~ 10 days) and a large optical luminosity (M_V = -18.7 mag). The SN properties are deduced well since the GRB afterglow does not contribute a significant amount to the total light output. The spectra show broad lines indicative of large expansion velocities, but are better matched by those of SN 2002ap and SN 1997ef (that are not associated with a GRB) than those of the proto-typical GRB-related SN 1998bw. We refine the redshift estimate to z = 0.0335 +/- 0.00007. The host-galaxy is a low-metallicity dwarf galaxy (with M_V \~ -16.0 mag), similar to host-galaxies of other GRB-associated SNe.

Numerical simulation is an important tool to help us understand the process of structure formation in the universe. However many simulation results of cold dark matter (CDM) halos on small scale are inconsistent with observations: the central density profile is too cuspy and there are too many substructures. Here we point out that these two problems may be connected with a hitherto unrecognized bias in the simulation halos. Although CDM halos in nature and in simulation are both virialized systems of collisionless CDM particles, gravitational encounter cannot be neglected in the simulation halos because they contain much less particles. We demonstrate this by two numerical experiments, showing that there is a difference on the microcosmic scale between the natural and simulation halos. The simulation halo is more akin to globular clusters where gravitational encounter is known to lead to such drastic phenomena as core collapse. And such artificial core collapse process appears to link the two problems together in the bottom-up scenario of structure formation in the $\Lambda$CDM universe. The discovery of this bias also has implications on the applicability of the Jeans Theorem in Galactic Dynamics.

Numerical simulation plays an important role in the study of structure formation of the universe. However, the mass resolution in current simulations is still poor. Due to technical difficulties, it is necessary to use both greatly reduced number density of particles and greatly raised unit particle mass. Consequently, the particle masses used in cosmological simulations are about $10^{70}$ times larger than the $Gev$ candidates in particle physics. This is a huge physical bias that cannot be neglected when interpreting the results of the simulations. Here we discuss how such a bias affects Cold Dark Matter (CDM) cosmological simulations. We find that the small scale properties of the CDM particle system are changed in two aspects. 1) An upper limit is imposed on the spatial resolution of the simulation results. 2) Most importantly, an unexpected short mean free path is produced, and the corresponding two body scattering cross section is close to the value expected in the Self-Interaction Dark Matter(SIDM) model. Since the mean free path of real CDM particle systems is much longer than that in the simulations, our results imply 1) that there is probably no 'cusp problem' in real CDM halos, and 2) that a much longer time is needed to form new virialized halos in real CDM particle systems than in the simulations. This last result can help us understand the 'substructure problem'. Our discussion can also explain why the massive halos in the simulations may have smaller concentration coefficients.