Search

The radio-quiet quasar SDSS J1240+1455 lies at a redshift of z=3.11, is surrounded by a Ly-alpha blob (LAB), and is absorbed by a proximate damped Ly-alpha system. In order to better define the morphology of the blob and determine its emission mechanism, we gathered deep narrow-band images isolating the Ly-alpha line of this object in linearly polarized light. We provide a deep intensity image of the blob, showing a filamentary structure extending up to 16'' (or ~122 physical kpc) in diameter. No significant polarization signal could be extracted from the data, but 95% probability upper limits were defined through simulations. They vary between ~3% in the central 0.75'' disk (after subtraction of the unpolarized quasar continuum) and ~10% in the 3.8-5.5'' annulus. The low polarization suggests that the Ly-alpha photons are emitted mostly in situ, by recombination and de-excitation in a gas largely ionized by the quasar ultraviolet light, rather than by a central source and scattered subsequently by neutral hydrogen gas. This blob shows no detectable polarization signal, contrary to LAB1, a brighter and more extended blob that is not related to the nearby active galactic nucleus (AGN) in any obvious way, and where a significant polarization signal of about 18% was detected., Comment: 15 pages, 8 figures; accepted for publication in A&A

This document summarizes the efforts of the EMMI Rapid Reaction Task Force on "Suppression and (re)generation of quarkonium in heavy-ion collisions at the LHC", centered around their 2019 and 2022 meetings. It provides a review of existing experimental results and theoretical approaches, including lattice QCD calculations and semiclassical and quantum approaches for the dynamical evolution of quarkonia in the quark-gluon plasma as probed in high-energy heavy-ion collisions. The key ingredients of the transport models are itemized to facilitate comparisons of calculated quantities such as reaction rates, binding energies, and nuclear modification factors. A diagnostic assessment of the various results is attempted and coupled with an outlook for the future., Comment: 54 pages, 19 figures

We present the Lensed Lyman-Alpha MUSE Arcs Sample (LLAMAS) selected from MUSE and HST observations of 17 lensing clusters. The sample consists of 603 continuum-faint (-23

The neutral intergalactic medium above redshift 6 is opaque to ionizing radiation, therefore one needs indirect measurements of the escape fraction of ionizing photons from galaxies of this epoch. Low-ionization state absorption lines are a common feature in the spectrum of galaxies, showing a diversity of strengths and shapes. Since these lines indicate the presence of neutral gas in front of the stars, they have been proposed to carry information on the escape of ionizing radiation from galaxies. We study which processes are responsible for the shape of the absorption lines, to better understand their origin. We then explore whether the absorption lines can be used to predict the escape fractions. Using a radiation-hydrodynamical zoom-in simulation and the radiative transfer code RASCAS, we generate mock CII 1334 and LyB lines of a virtual galaxy at redshift 3 as seen from many directions of observation. We also compute the escape fraction of ionizing photons in those directions and look for correlations between the lines and the escape fractions. We find that the resulting mock absorption lines are comparable to observations and that the lines and the escape fractions vary strongly depending on the direction of observation. Gas velocity and dust always affect the absorption profile significantly. We find no strong correlations between observable LyB or CII 1334 and the escape fraction. After correcting the continuum for attenuation by dust to recover the intrinsic continuum, the residual flux of CII 1334 correlates well with the escape fraction for directions with a dust corrected residual flux larger than 30%. For other directions, the relations have a strong dispersion, and the residual flux overestimates the escape fraction for most cases. Concerning LyB, the residual flux after dust correction does not correlate with the escape fraction but can be used as a lower limit. (abridged), Comment: 25 pages, 22 figures, accepted for publication in Astronomy and Astrophysics. New title, shorter abstract, corrected typos and url for Krome website

We present a study of the intersection between the populations of star forming galaxies (SFGs) selected as either Lyman break galaxies (LBGs) or Lyman-alpha emitters (LAEs) in the redshift range 2.9 - 6.7, within the same volume of universe sampled by the Multi-Unit Spectroscopic Explorer (MUSE) behind the Hubble Frontier Fields lensing cluster A2744. We define three samples of star-forming galaxies: LBG galaxies with an LAE counterpart (92 galaxies), LBG galaxies without LAE counterpart (408 galaxies) and LAE galaxies without an LBG counterpart (46 galaxies). All these galaxies are intrinsically faint due to the lensing nature of the sample (Muv $\ge$ -20.5). The fraction of LAEs among all selected SFGs increases with redshift up to z $\sim$ 6 and decreases for higher redshifts. The evolution of LAE/LBG populations with UV magnitude and Lya luminosity shows that the LAE selection is able to identify intrinsically UV faint galaxies with Muv $\ge$ -15 that are typically missed in the deepest lensing photometric surveys. The LBG population seems to fairly represent the total population of SFGs down to Muv$\sim$-15. Galaxies with Muv$<-17$ tend to have SFRLya$<$SFRuv, whereas the opposite trend is observed within our sample for faint galaxies with Muv$>-17$, including galaxies only detected by their Lya emission, with a large scatter. These trends, previously observed in other samples of SFGs at high-$z$, are seen here for very faint Muv$\sim -15$ galaxies, much fainter than in previous studies. There is no clear evidence, based on the present results, for an intrinsic difference on the properties of the two populations selected as LBG and/or LAE. The observed trends could be explained by a combination of several facts, like the existence of different star-formation regimes, the dust content, the relative distribution and morphology of dust and stars, or the stellar populations

We simulate the flux emitted from galaxy halos in order to quantify the brightness of the circumgalactic medium (CGM). We use dedicated zoom-in cosmological simulations with the hydrodynamical Adaptive Mesh Refinement code RAMSES, which are evolved down to z=0 and reach a maximum spatial resolution of 380 $h^{-1}$pc and a gas mass resolution up to 1.8$\times 10^{5} h^{-1} \rm{M}_{\odot}$ in the densest regions. We compute the expected emission from the gas in the CGM using CLOUDY emissivity models for different lines (e.g. Ly$\alpha$, CIV, OVI, CVI, OVIII) considering UV background fluorescence, gravitational cooling and continuum emission. In the case of Ly$\alpha$ we additionally consider the scattering of continuum photons. We compare our predictions to current observations and find them to be in good agreement at any redshift after adjusting the Ly$\alpha$ escape fraction. We combine our mock observations with instrument models for FIREBall-2 (UV balloon spectrograph) and HARMONI (visible and NIR IFU on the ELT) to predict CGM observations with either instrument and optimise target selections and observing strategies. Our results show that Ly$\alpha$ emission from the CGM at a redshift of 0.7 will be observable with FIREBall-2 for bright galaxies (NUV$\sim$18 mag), while metal lines like OVI and CIV will remain challenging to detect. HARMONI is found to be well suited to study the CGM at different redshifts with various tracers., Comment: accepted for publication in MNRAS

We present the deepest study to date of the Lya luminosity function (LF) in a blank field using blind integral field spectroscopy from MUSE. We constructed a sample of 604 Lya emitters (LAEs) across the redshift range 2.91 < z < 6.64 using automatic detection software in the Hubble Ultra Deep Field. We calculate accurate total Lya fluxes capturing low surface brightness extended Lya emission now known to be a generic property of high-redshift star-forming galaxies. We simulated realistic extended LAEs to characterise the selection function of our samples, and performed flux-recovery experiments to test and correct for bias in our determination of total Lya fluxes. We find an accurate completeness correction accounting for extended emission reveals a very steep faint-end slope of the LF, alpha, down to luminosities of log10 L erg s^-1< 41.5, applying both the 1/Vmax and maximum likelihood estimators. Splitting the sample into three broad redshift bins, we see the faint-end slope increasing from -2.03+1.42-inf at z ~ 3.44 to -2.86+0.76-inf at z ~ 5.48, however no strong evolution is seen between the 68% confidence regions in L*-alpha parameter space. Using the Lya line flux as a proxy for star formation activity, and integrating the observed LFs, we find that LAEs' contribution to the cosmic SFRD rises with redshift until it is comparable to that from continuum-selected samples by z ~ 6. This implies that LAEs may contribute more to the star-formation activity of the early Universe than previously thought - any additional interglactic medium correction would act to further boost the Lya luminosities. Finally, assuming fiducial values for the escape of Lya and LyC radiation, and the clumpiness of the IGM, we integrated the maximum likelihood LF at 5.00 < z < 6.64 and find we require only a small extrapolation beyond the data (< 1 dex in L) for LAEs alone to maintain an ionised IGM at z ~ 6., Comment: 16 pages, 10 figures. Accepted for publication in A&A special issue: "The MUSE Hubble Ultra Deep Field Survey" (MUSE UDF Series Paper VI)

GalICS 2.0 is a new semianalytic code to model the formation and evolution of galaxies in a cosmological context. N-body simulations based on a Planck cosmology are used to construct halo merger trees, track subhaloes, compute spins and measure concentrations. The accretion of gas onto galaxies and the morphological evolution of galaxies are modelled with prescriptions derived from hydrodynamic simulations. Star formation and stellar feedback are described with phenomenological models (as in other semianalytic codes). GalICS 2.0 computes rotation speeds from the gravitational potential of the dark matter, the disc and the central bulge. As the rotation speed depends not only on the virial velocity but also on the ratio of baryons to dark matter within a galaxy, our calculation predicts a different Tully-Fisher relation from models in which the rotation speed is proportional to the virial velocity. This is why GalICS 2.0 is able to reproduce the galaxy stellar mass function and the Tully-Fisher relation simultaneously. Our results are also in agreement with halo masses from weak lensing and satellite kinematics, gas fractions, the relation between star formation rate (SFR) and stellar mass, the evolution of the cosmic SFR density, bulge-to-disc ratios, disc sizes and the Faber-Jackson relation., Comment: 31 pages, 19 figures, 1 table, submitted to MNRAS

[Abridged] We use the final data of the VIMOS Public Extragalactic Redshift Survey (VIPERS) to investigate the effect of environment on the evolution of galaxies between $z=0.5$ and $z=0.9$. We characterise local environment in terms of the density contrast smoothed over a cylindrical kernel, the scale of which is defined by the distance to the $5^{th}$ nearest neighbour. We find that more massive galaxies tend to reside in higher-density environments over the full redshift range explored. Defining star-forming and passive galaxies through their (NUV$-r$) vs ($r-K$) colours, we then quantify the fraction of star-forming over passive galaxies, $f_{\rm ap}$, as a function of environment at fixed stellar mass. $f_{\rm ap}$ is higher in low-density regions for galaxies with masses ranging from $\log(\mathcal{M}/\mathcal{M}_\odot)=10.38$ (the lowest value explored) to at least $\log(\mathcal{M}/\mathcal{M}_\odot)\sim11.3$, although with decreasing significance going from lower to higher masses. This is the first time that environmental effects on high-mass galaxies are clearly detected at redshifts as high as $z\sim0.9$. We compared these results to VIPERS-like galaxy mock catalogues based on the galaxy formation model of De Lucia & Blaizot. The model correctly reproduces $f_{\rm ap}$ in low-density environments, but underpredicts it at high densities. The discrepancy is particularly strong for the lowest-mass bins. We find that this discrepancy is driven by an excess of low-mass passive satellite galaxies in the model. Looking at the accretion history of these model galaxies, i.e. the times when they become satellites, a better (yet not perfect) agreement with observations can be obtained in high density regions by assuming either that a not-negligible fraction of satellites is destroyed, or that their quenching time-scale is longer than $\sim 2$ Gyr., Comment: 20 pages, 13 figures, accepted to be published in A&A

Lyman-$\alpha$ nebulae are usually found in massive environments at high redshift ($z > 2$). The origin of their Lyman-$\alpha$ (Lya) emission remains debated. Recent polarimetric observations showed that at least some Lya sources are polarized. This is often interpreted as a proof that the photons are centrally produced, and opposed to the scenario in which the Lya emission is the cooling radiation emitted by gas heated during the accretion onto the halo. We suggest that this scenario is not incompatible with the polarimetric observations. In order to test this idea, we post-process a radiative hydrodynamics simulation of a blob with the MCLya Monte Carlo transfer code. We compute radial profiles for the surface brightness and the degree of polarization and compare them to existing observations. We find that both are consistent with a significant contribution of the extragalactic gas to the Lya emission. Most of the photons are centrally emitted and scattered inside the filament afterwards, producing the observed high level of polarization. We argue that the contribution of the extragalactic gas to the Lya emission does not prevent polarization to arise. On the contrary, we find that pure galactic emission causes the polarization profile to be too steep to be consistent with observations., Comment: 14 pages, 9 figures (+4 in appendices). Accepted for publication in Astronomy & Astrophysics

Small galaxies are thought to be the main contributors to the ionising budget of the Universe before reionisation was complete. There have been a number of numerical studies trying to quantify their ionising efficiency through the escape fraction $f_{esc}$. While there is a clear trend that $f_{esc}$ is higher for smaller haloes, there is a large scatter in the distribution of $f_{esc}$ for a single halo mass. We propose that this is due to the intrinsic burstiness of star formation in low mass galaxies. We performed high resolution radiative hydrodynamics simulations with Ramses-RT to model the evolution of three galaxies and their ionising efficiency. We found that the variability of $f_{esc}$ follows that of the star formation rate. We then discuss the consequences of this variability on the observability of such galaxies by JWST., Comment: 4 pages, 3 figures, to appear in the proceedings of SF2A 2015

We report the detection of extended Ly alpha emission around individual star-forming galaxies at redshifts z = 3-6 in an ultradeep exposure of the Hubble Deep Field South obtained with MUSE on the ESO-VLT. The data reach a limiting surface brightness (1sigma) of ~1 x 10^-19 erg s^-1 cm^-2 arcsec^-2 in azimuthally averaged radial profiles, an order of magnitude improvement over previous narrowband imaging. Our sample consists of 26 spectroscopically confirmed Ly alpha-emitting, but mostly continuum-faint (m_AB >~ 27) galaxies. In most objects the Ly alpha emission is considerably more extended than the UV continuum light. While 5 of the faintest galaxies in the sample show no significantly detected Ly alpha haloes, the derived upper limits suggest that this is just due to insufficient S/N. Ly alpha haloes therefore appear to be (nearly) ubiquitous even for low-mass (~10^8-10^9 M_sun) star-forming galaxies at z>3. We decompose the Ly alpha emission of each object into a compact `continuum-like' and an extended halo component, and infer sizes and luminosities of the haloes. The extended Ly alpha emission approximately follows an exponential surface brightness distribution with a scale length of a few kpc. While these haloes are thus quite modest in terms of their absolute sizes, they are larger by a factor of 5-15 than the corresponding rest-frame UV continuum sources as seen by HST. They are also much more extended, by a factor ~5, than Ly alpha haloes around low-redshift star-forming galaxies. Between ~40% and >90% of the observed Ly alpha flux comes from the extended halo component, with no obvious correlation of this fraction with either the absolute or the relative size of the Ly alpha halo. Our observations provide direct insights into the spatial distribution of at least partly neutral gas residing in the circumgalactic medium of low to intermediate mass galaxies at z > 3., Comment: accepted by Astronomy & Astrophysics

We explore the properties of high-redshift Lyman-alpha emitters (LAE), and their link with the Lyman-Break galaxy population (LBG), using a semi-analytic model of galaxy formation that takes into account resonant scattering of Lya photons in gas outflows. We can reasonably reproduce the abundances of LAEs and LBGs from redshift 3 to 7, as well as most UV LFs of LAEs. The stronger dust attenuation for (resonant) Lya photons compared to UV continuum photons in bright LBGs provides a natural interpretation to the increase of the LAE fraction in LBG samples, X_LAE, towards fainter magnitudes. The redshift evolution of X_LAE seems however very sensitive to UV magnitudes limits and EW cuts. In spite of the apparent good match between the statistical properties predicted by the model and the observations, we find that the tail of the Lya equivalent width distribution (EW > 100 A) cannot be explained by our model, and we need to invoke additional mechanisms. We find that LAEs and LBGs span a very similar dynamical range, but bright LAEs are about 4 times rarer than LBGs in massive halos. Moreover, massive halos mainly contain weak LAEs in our model, which might introduce a bias towards low-mass halos in surveys which select sources with high EW cuts. Overall, our results are consistent with the idea that LAEs and LBGs make a very similar galaxy population. Their apparent differences seem mainly due to EW selections, UV detection limits, and a decreasing Lya-to-UV escape fraction ratio in high SFR galaxies., Comment: accepted in MNRAS

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field., Comment: 21 pages, 20 figures. Note that the catalogues, data cube and associated spectra will be released upon acceptance of the paper

We develop a Principal Component Analysis aimed at classifying a sub-set of 27,350 spectra of galaxies in the range 0.4 < z < 1.0 collected by the VIMOS Public Extragalactic Redshift Survey (VIPERS). We apply an iterative algorithm to simultaneously repair parts of spectra affected by noise and/or sky residuals, and reconstruct gaps due to rest-frame transformation, and obtain a set of orthogonal spectral templates that span the diversity of galaxy types. By taking the three most significant components, we find that we can describe the whole sample without contamination from noise. We produce a catalogue of eigen-coefficients and template spectra that will be part of future VIPERS data releases. Our templates effectively condense the spectral information into two coefficients that can be related to the age and star formation rate of the galaxies. We examine the spectrophotometric types in this space and identify early, intermediate, late and starburst galaxies., Comment: 15 pages, 20 images, accepted for publication in MNRAS: MN-12-1739-MJ.R1

Simulations of galaxy evolution aim to capture our current understanding as well as to make predictions for testing by future experiments. Simulations and observations are often compared in an indirect fashion: physical quantities are estimated from the data and compared to models. However, many applications can benefit from a more direct approach, where the observing process is also simulated and the models are seen fully from the observer's perspective. To facilitate this, we have developed the Millennium Run Observatory (MRObs), a theoretical virtual observatory which uses virtual telescopes to `observe' semi-analytic galaxy formation models based on the suite of Millennium Run dark matter simulations. The MRObs produces data that can be processed and analyzed using the standard software packages developed for real observations. At present, we produce images in forty filters from the rest-frame UV to IR for two stellar population synthesis models, three different models of IGM absorption, and two cosmologies (WMAP1/7). Galaxy distributions for a large number of mock lightcones can be `observed' using models of major ground- and space-based telescopes. The data include lightcone catalogues linked to structural properties of galaxies, pre-observation model images, mock telescope images, and Source Extractor products that can all be traced back to the higher level dark matter, semi-analytic galaxy, and lightcone catalogues available in the Millennium database. Here, we describe our methods and announce a first public release of simulated surveys (e.g., SDSS, CFHT-LS, GOODS, GOODS/ERS, CANDELS, and HUDF). The MRObs browser, an online tool, further facilitates exploration of the simulated data. We demonstrate the benefits of a direct approach through a number of example applications (galaxy number counts in CANDELS, clusters, morphologies, and dropout selections)., Comment: MNRAS, in press. Millennium Run Observatory data products, online tools, and more available through http://galformod.mpa-garching.mpg.de/mrobs/

[Abridged] We perform on galaxy mock catalogues the same colour-density analysis made by Cucciati et al. (2006) on a 5 Mpc/h scale using the VVDS-Deep survey, and compare the results from mocks with observed data. We use mocks with the same flux limits (I=24) as the VVDS (CMOCKS), built using the semi- analytic model by De Lucia & Blaizot (2007) applied to the Millennium Simulation. From CMOCKS, we extracted samples of galaxies mimicking the VVDS observational strategy (OMOCKS). We computed the B-band Luminosity Function LF and the colour-density relation (CDR) in the mocks. We find that the LF in mocks roughly agrees with the observed LF, but at z<0.8 the faint-end slope of the model LF is steeper than the VVDS one. Computing the LF for early and late type galaxies, we show that mocks have an excess of faint early-type and of bright late-type galaxies with respect to data. We find that the CDR in OMOCKS is in excellent agreement with the one in CMOCKS. At z~0.7, the CDR in mocks agrees with the VVDS one (red galaxies reside mainly in high densities). Yet, the strength of the CDR in mocks does not vary within 0.2

{Abridged} We investigate the observability of cold accretion streams at redshift 3 via Lyman-alpha (Lya) emission and the feasibility of cold accretion as the main driver of Lya blobs (LABs). We run cosmological zoom simulations focusing on 3 halos spanning two orders of magnitude in mass, roughly from 10^11 to 10^13 solar masses. We use a version of the Ramses code that includes radiative transfer of UV photons, and we employ a refinement strategy that allows us to resolve accretion streams in their natural environment to an unprecedented level. For the first time, we self-consistently model self-shielding in the cold streams from the cosmological UV background, which enables us to predict their temperatures, ionization states and Lya luminosities with improved accuracy. We find the efficiency of gravitational heating in cold streams in a ~10^11 solar mass halo is around 10-20% throughout most of the halo but reaching much higher values close to the center. As a result most of the Lya luminosity comes from gas which is concentrated at the central 20% of the halo radius, leading to Lya emission which is not extended. In more massive halos, of >10^12 solar masses, cold accretion is complex and disrupted, and gravitational heating does not happen as a steady process. Ignoring the factors of Lya scattering, local UV enhancement, and SNe feedback, cold accretion alone in these massive halos can produce LABs that largely agree with observations in terms of morphology, extent, and luminosity. Our simulations slightly and systematically over-predict LAB abundances, perhaps hinting that the interplay of these ignored factors may have a negative net effect on extent and luminosity. We predict that a factor of a few increase in sensitivity from current observational limits should unambiguously reveal continuum-free accretion streams around massive galaxies at z=3., Comment: 25 pages, 21 figures, 3 tables. Accepted for publication in MNRAS

We present for the first time predictions for UV line emission of intergalactic and circumgalactic gas from Adaptive Mesh Resolution (AMR) Large Scale Structure (LSS) simulations at redshifts 0.3

We present in detail the implementation of the Blaizot-M\'endez-Wschebor (BMW) approximation scheme of the nonperturbative renormalization group, which allows for the computation of the full momentum dependence of correlation functions. We discuss its signification and its relation with other schemes, in particular the derivative expansion. Quantitative results are presented for the testground of scalar O(N) theories. Besides critical exponents which are zero-momentum quantities, we compute in three dimensions in the whole momentum range the two-point function at criticality and, in the high temperature phase, the universal structure factor. In all cases, we find very good agreement with the best existing results., Comment: 21 pages, 7 figures. Added some minor corrections

We have developed a method for detecting clusters in large imaging surveys, based on the detection of structures in galaxy density maps made in slices of photometric redshifts. This method was first applied to the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Deep 1 field by Mazure et al. (2007), then to all the Deep and Wide CFHTLS fields available in the T0004 data release by Adami et al. (2010). The validity of the cluster detection rate was estimated by applying the same procedure to galaxies from the Millennium simulation. Here we analyse with the same method the full CFHTLS Wide survey, based on the T0006 data release. In a total area of 154 deg2, we have detected 4061 candidate clusters at 3sigma or above (6802 at 2sigma and above), in the redshift range 0.1<=z<=1.15, with estimated mean masses between 1.3 10^14 and 12.6 10^14 M_solar. This catalogue of candidate clusters will be available online via VizieR. We compare our detections with those made in various CFHTLS analyses with other methods. By stacking a subsample of clusters, we show that this subsample has typical cluster characteristics (colour-magnitude relation, galaxy luminosity function). We also confirm that the cluster-cluster correlation function is comparable to that obtained for other cluster surveys and analyze large scale filamentary galaxy distributions. We have increased the number of known optical high redshift cluster candidates by a large factor, an important step towards obtaining reliable cluster counts to measure cosmological parameters. The clusters that we detect behave as expected for a sample of clusters fed by filaments at the intersection of which they are located., Comment: Accepted for publication in Astronomy \& Astrophysics

This paper presents a detailed comparison between high-redshift observations from the VIMOS-VLT Deep Survey (VVDS) and predictions from the Munich semi-analytical model of galaxy formation. In particular, we focus this analysis on the magnitude, redshift, and colour distributions of galaxies, as well as their clustering properties. We constructed 100 quasi-independent mock catalogues, using the output of the semi-analytical model presented in De Lucia & Blaizot (2007).We then applied the same observational selection function of the VVDS-Deep survey, so as to carry out a fair comparison between models and observations. We find that the semi-analytical model reproduces well the magnitude counts in the optical bands. It tends, however, to overpredict the abundance of faint red galaxies, in particular in the i' and z' bands. Model galaxies exhibit a colour bimodality that is only in qualitative agreement with the data. In particular, we find that the model tends to overpredict the number of red galaxies at low redshift and of blue galaxies at all redshifts probed by VVDS-Deep observations, although a large fraction of the bluest observed galaxies is absent from the model. In addition, the model overpredicts by about 14 per cent the number of galaxies observed at 0.2

We propose a model for the propagation of a heavy-quark in a hot plasma, to be viewed as a first step towards a full description of the dynamics of heavy quark systems in a quark-gluon plasma, including bound state formation. The heavy quark is treated as a non relativistic particle interacting with a fluctuating field, whose correlator is determined by a hard thermal loop approximation. This approximation, which concerns only the medium in which the heavy quark propagates, is the only one that is made, and it can be improved. The dynamics of the heavy quark is given exactly by a quantum mechanical path integral that is calculated in this paper in the Euclidean space-time using numerical Monte Carlo techniques. The spectral function of the heavy quark in the medium is then reconstructed using a Maximum Entropy Method. The path integral is also evaluated exactly in the case where the mass of the heavy quark is infinite; one then recovers known results concerning the complex optical potential that controls the long time behavior of the heavy quark. The heavy quark correlator and its spectral function is also calculated semi-analytically at the one-loop order, which allows for a detailed description of the coupling between the heavy quark and the plasma collective modes.

We study the gluon distribution in nucleus-nucleus collisions in the framework of the Color-Glass-Condensate. Approximate analytical solutions are compared to numerical solutions of the non-linear Yang-Mills equations. We find that the full numerical solution can be well approximated by taking the full initial condition of the fields in Coulomb gauge and using a linearized solution for the time evolution. We also compare kt-factorized approximations to the full solution., Comment: 13 pages, 12 figures, RevTeX

[Abridged] We present a homogeneous and complete catalogue of optical groups identified in the purely flux limited (17.5<=I<=24.0) VIMOS-VLT Deep Survey (VVDS). We use mock catalogues extracted from the MILLENNIUM simulation, to correct for potential systematics that might affect the overall distribution as well as the individual properties of the identified systems. Simulated samples allow us to forecast the number and properties of groups that can be potentially found in a survey with VVDS-like selection functions. We use them to correct for the expected incompleteness and also to asses how well galaxy redshifts trace the line-of-sight velocity dispersion of the underlying mass overdensity. In particular, we train on these mock catalogues the adopted group-finding technique (the Voronoi-Delaunay Method, VDM). The goal is to fine-tune its free parameters, recover in a robust and unbiased way the redshift and velocity dispersion distributions of groups and maximize the level of completeness (C) and purity (P) of the group catalogue. We identify 318 VVDS groups with at least 2 members within 0.2<=z<=1.0, among which 144 (/30) with at least 3 (/5) members. The sample has globally C=60% and P=50%. Nearly 45% of the groups with at least 3 members are still recovered if we run the algorithm with a parameter set which maximizes P (75%). We exploit the group sample to study the redshift evolution of the fraction f_b of blue galaxies (U-B<=1) within 0.2<=z<=1. We find that f_b is significantly lower in groups than in the whole ensemble of galaxies irrespectively of their environment. These quantities increase with redshift, with f_b in groups showing a marginally significant steeper increase. We also confirm that, at any explored redshift, f_b decreases for increasing group richness, and we extend towards fainter luminosities the magnitude range over which this result holds., Comment: Submitted to A&A, revised version after referee comments, Table 5 added

In order to enlarge publicly available optical cluster catalogs, in particular at high redshift, we have performed a systematic search for clusters of galaxies in the CFHTLS. We used the Le Phare photometric redshifts for the galaxies detected with magnitude limits of i'=25 and 23 for the Deep and Wide fields respectively. We then constructed galaxy density maps in photometric redshift bins of 0.1 based on an adaptive kernel technique and detected structures with SExtractor. In order to assess the validity of our cluster detection rates, we applied a similar procedure to galaxies in Millennium simulations. We measured the correlation function of our cluster candidates. We analyzed large scale properties and substructures by applying a minimal spanning tree algorithm both to our data and to the Millennium simulations. We have detected 1200 candidate clusters with various masses (minimal masses between 1.0 10$^{13}$ and 5.5 10$^{13}$ and mean masses between 1.3 10$^{14}$ and 12.6 10$^{14}$ M$_\odot$), thus notably increasing the number of known high redshift cluster candidates. We found a correlation function for these objects comparable to that obtained for high redshift cluster surveys. We also show that the CFHTLS deep survey is able to trace the large scale structure of the universe up to z$\geq$1. Our detections are fully consistent with those made in various CFHTLS analyses with other methods. We now need accurate mass determinations of these structures to constrain cosmological parameters., Comment: Accepted in A&A

The heavy-quark spectral function in a hot plasma is reconstructed from the corresponding euclidean propagator. The latter is evaluated through a path-integral simulation. A weak-coupling calculation is also performed, allowing to interpret the qualitative behavior of the spectral function in terms of quite general physical processes., Comment: 4 pages, 3 figures - To appear in the conference proceedings for Quark Matter 2009, March 30 - April 4, Knoxville, Tennessee

Although a fair amount of work has been devoted to growing Monte-Carlo merger trees which resemble those built from an N-body simulation, comparatively little effort has been invested in quantifying the caveats one necessarily encounters when one extracts trees directly from such a simulation. To somewhat revert the tide, this paper seeks to provide its reader with a comprehensive study of the problems one faces when following this route. The first step to building merger histories of dark matter haloes and their subhaloes is to identify these structures in each of the time outputs (snapshots) produced by the simulation. Even though we discuss a particular implementation of such an algorithm (called AdaptaHOP) in this paper, we believe that our results do not depend on the exact details of the implementation but extend to most if not all (sub)structure finders. We then highlight different ways to build merger histories from AdaptaHOP haloes and subhaloes, contrasting their various advantages and drawbacks. We find that the best approach to (sub)halo merging histories is through an analysis that goes back and forth between identification and tree building rather than one which conducts a straightforward sequential treatment of these two steps. This is rooted in the complexity of the merging trees which have to depict an inherently dynamical process from the partial temporal information contained in the collection of instantaneous snapshots available from the N-body simulation., Comment: 19 pages, 28 figures

We demonstrate the power of a recently-proposed approximation scheme for the non-perturbative renormalization group that gives access to correlation functions over their full momentum range. We solve numerically the leading-order flow equations obtained within this scheme, and compute the two-point functions of the O(N) theories at criticality, in two and three dimensions. Excellent results are obtained for both universal and non-universal quantities at modest numerical cost., Comment: 4 pages, 1 figure

We have investigated the dependence of galaxy clustering on their stellar mass at z~1, using the data from the VIMOS-VLT Deep Survey (VVDS). We have measured the projected two-point correlation function of galaxies, wp(rp) for a set of stellar mass selected samples at an effective redshift =0.85. We have control and quantify all effects on galaxy clustering due to the incompleteness of our low mass samples. We find that more massive galaxies are more clustered. When compared to similar results at z~0.1 in the SDSS, we observed no evolution of the projected correlation function for massive galaxies. These objects present a stronger linear bias at z~1 with respect to low mass galaxies. As expected, massive objects at high redshift are found in the highest pics of the dark matter density field., Comment: 4 pages, 2 figures, 43rd Rencontres de Moriond - March 15-22, 2008 - La Thuile (Val d'Aosta, Italy)

The VVDS-Wide survey has been designed with the general aim of tracing the large-scale distribution of galaxies at z~1 on comoving scales reaching ~100Mpc/h, while providing a good control of cosmic variance over areas as large as a few square degrees. This is achieved by measuring redshifts with VIMOS at the ESO VLT to a limiting magnitude I_AB=22.5, targeting four independent fields with size up to 4 sq.deg. each. The whole survey covers 8.6 sq.deg., here we present the general properties of the current redshift sample. This includes 32734 spectra in the four regions (19977 galaxies, 304 type I AGNs, and 9913 stars), covering a total area of 6.1 sq.deg, with a sampling rate of 22 to 24%. The redshift success rate is above 90% independently of magnitude. It is the currently largest area coverage among redshift surveys reaching z~1. We give the mean N(z) distribution averaged over 6.1 sq.deg. Comparing galaxy densities from the four fields shows that in a redshift bin Deltaz=0.1 at z~1 one still has factor-of-two variations over areas as large as ~0.25 sq.deg. This level of cosmic variance agrees with that obtained by integrating the galaxy two-point correlation function estimated from the F22 field alone, and is also in fairly good statistical agreement with that predicted by the Millennium mocks. The variance estimated over the survey fields shows explicitly how clustering results from deep surveys of even ~1 sq.deg. size should be interpreted with caution. This paper accompanies the public release of the first 18143 redshifts of the VVDS-Wide survey from the 4 sq.deg. contiguous area of the F22 field at RA=22h, publicly available at http://cencosw.oamp.fr, Comment: Accepted for publication on Astronomy & Astrophysics

Observations of distant supernovae indicate that the Universe is now in a phase of accelerated expansion the physical cause of which is a mystery. Formally, this requires the inclusion of a term acting as a negative pressure in the equations of cosmic expansion, accounting for about 75 per cent of the total energy density in the Universe. The simplest option for this "dark energy" corresponds to a cosmological constant, perhaps related to the quantum vacuum energy. Physically viable alternatives invoke either the presence of a scalar field with an evolving equation of state, or extensions of general relativity involving higher-order curvature terms or extra dimensions. Although they produce similar expansion rates, different models predict measurable differences in the growth rate of large-scale structure with cosmic time. A fingerprint of this growth is provided by coherent galaxy motions, which introduce a radial anisotropy in the clustering pattern reconstructed by galaxy redshift surveys. Here we report a measurement of this effect at a redshift of 0.8. Using a new survey of more than 10,000 faint galaxies, we measure the anisotropy parameter b = 0.70 +/- 0.26, which corresponds to a growth rate of structure at that time of f = 0.91 +/- 0.36. This is consistent with the standard cosmological-constant model with low matter density and flat geometry, although the error bars are still too large to distinguish among alternative origins for the accelerated expansion. This could be achieved with a further factor-of-ten increase in the sampled volume at similar redshift., Comment: One PDF file including both main paper and Supplementary Information (28 pages, 3+2 figures). Published version available at http://www.nature.com/nature/journal/v451/n7178/abs/nature06555.html

We investigate the behavior of a pair of heavy fermions, denoted by $Q$ and $\bar{Q}$, in a hot/dense medium. Although we have in mind the situation where $Q$ and $\bar{Q}$ denote heavy quarks, our treatment will be limited to simplified models, which bear only some general similarities with QCD. We study in particular the limiting case where the mass of the heavy fermions is infinite. Then a number of results can be derived exactly: a Schr\"odinger equation can be established for the correlator of the heavy quarks; the interaction effects exponentiate, leading to a simple instantaneous effective potential for this Schr\"odinger equation. We consider simple models for the medium in which the $Q\bar Q$ pair propagates. In the case where the medium is a plasma of photons and light charged fermions, an imaginary part develops in this effective potential. We discuss the physical interpretation of this imaginary part in terms of the collisions between the heavy particles and the light fermions of the medium; the same collisions also determine the damping rate of the heavy fermions. Finally we study the connection between the real-time propagator of the heavy fermion pair and its Euclidean counterpart, and show that the real part of the potential entering the Schr\"odinger equation for the real-time propagator is the free energy calculated in the imaginary-time formalism., Comment: 32 pages, 8 figures

Aims: We use the VVDS-Deep first-epoch data to measure the dependence of galaxy clustering on galaxy stellar mass, at z~0.85. Methods: We measure the projected correlation function wp(rp) for sub-samples with 0.510^9 Msun to r0 ~ 4.28 h^-1 Mpc for galaxies more massive than 10^10.5 Msun. At the same time, the slope increases from ~ 1.67 to ~ 2.28. A comparison of the observed wp(rp) to local measurements by the SDSS shows that the evolution is faster for objects less massive than ~10^10.5 Msun. This is interpreted as a higher dependence on redshift of the linear bias b_L for the more massive objects. While for the most massive galaxies b_L decreases from 1.5+/-0.2 at z~0.85 to 1.33+/-0.03 at z~0.15, the less massive population maintains a virtually constant value b_L~1.3. This result is in agreement with a scenario in which more massive galaxies formed at high redshift in the highest peaks of the density field, while less massive objects form at later epochs from the more general population of dark-matter halos., Comment: 13 pages, 10 figures, accepted in A&A

We apply the previously proposed scheme of approximately self-consistent hard-thermal-loop resummations in the entropy of high-temperature QCD to N=4 supersymmetric Yang-Mills (SYM) theories and compare with a (uniquely determined) R[4,4] Pad\'e approximant that interpolates accurately between the known perturbative result and the next-to-leading order strong-coupling result obtained from AdS/CFT correspondence. We find good agreement up to couplings where the entropy has dropped to about 85% of the Stefan-Boltzmann value. This is precisely the regime which in purely gluonic QCD corresponds to temperatures above 2.5 times the deconfinement temperature and for which this method of hard-thermal-loop resummation has given similar good agreement with lattice QCD results. This suggests that in this regime the entropy of both QCD and N=4 SYM is dominated by effectively weakly coupled hard-thermal-loop quasiparticle degrees of freedom. In N=4 SYM, strong-coupling contributions to the thermodynamic potential take over when the entropy drops below 85% of the Stefan-Boltzmann value., Comment: 14 pages, 2 figures, JHEP3. v2: revised and expanded, with unchanged HTL results but corrected NLO strong-coupling result from AdS/CFT (which is incorrectly reproduced in almost all previous papers comparing weak and strong coupling results of N=4 SYM) and novel (unique) Pade approximant interpolating between weak and strong coupling results

We investigate a zero-dimensional toy model originally introduced by Mueller and Salam which mimics high-energy scattering in QCD in the presence of both gluon saturation and gluon number fluctuations, and hence of Pomeron loops. Unlike other toy models of the reaction-diffusion type, the model studied in this paper is consistent with boost invariance and, related to that, it exhibits a mechanism for particle saturation close to that of the JIMWLK equation in QCD, namely the saturation of the emission rate due to high-density effects. Within this model, we establish the dominant high-energy behaviour of the S-matrix element for the scattering between a target obtained by evolving one particle and a projectile made with exactly n particles. Remarkably, we find that all such matrix elements approach the black disk limit S=0 at high rapidity Y, with the same exponential law: ~ exp(-Y) for all values of n. This is so because the S-matrix is dominated by rare target configurations which involve only few particles. We also find that the bulk distribution for a saturated system is of the Poisson type., Comment: 34 pages, 9 figures. Some explanations added on the frame-dependence of the relevant configurations (new section 3.3)

We show how the Landau-Pomeranchuk-Migdal effect on photon production rates in a quark-gluon plasma can be derived via the usual Boltzmann equation. To do this, we first derive the electromagnetic polarization tensor using linear response theory, and then formulate the Boltzmann equation including the collisions mediated by soft gluon exchanges. We then identify the resulting expression for the production rate with that obtained by the field-theoretic formalism recently proposed by Arnold, Moore and Yaffe. To illustrate the LPM effect we solve the Boltzmann equation in the diffusion approximation., Comment: 17 pages, 4 figures Minor rewriting. In particular, the specific small angle behavior of the collision term involved in the LPM effect is better explained. Version accepted for publication in Physics Letters B

[Abridged] We use mock catalogues based on the GALICS model (Hatton et al. 03) to explore the nature of galaxy clustering observed in the SDSS. We measure low and high order angular clustering statistic from these mock catalogues, after selecting galaxies the same way as for observations, and compare them directly to estimates from SDSS data. Note that we also present measurements of S3-S5 on the SDSS DR1. We find that our model is in general good agreement with observations in the scale/luminosity range where we can trust the predictions. This range is found to be limited (i) by the size of the dark matter simulation used -- which introduces finite volume effects at large scales -- and by the mass resolution of this simulation -- which introduces incompleteness at apparent magnitudes fainter than $r\sim 20$. We then focus on the small scale clustering properties of galaxies and investigate the behaviour of three different prescriptions for positioning galaxies within haloes of dark matter. We show that galaxies are poor tracers both of DM particles or DM sub-structures, within groups and clusters. Instead, SDSS data tells us that the distribution of galaxies lies somewhat in between these two populations. This confirms the general theoretical expectation from numerical simulations and semi-analytic modelling., Comment: MNRAS, in press

We reproduce the blue and red sequences in the observed joint distribution of colour and magnitude for galaxies at low and high redshifts using hybrid N-body/semi-analytic simulations of galaxy formation. The match of model and data is achieved by mimicking the effects of cold flows versus shock heating coupled to feedback from active galactic nuclei (AGNs), as predicted by Dekel & Birnboim (2006). After a critical epoch z=3, only haloes below a critical shock-heating mass 10^12MSun enjoy gas supply by cold flows and form stars, while cooling and star formation are shut down abruptly above this mass. The shock-heated gas is kept hot because being dilute it is vulnerable to feedback from energetic sources such as AGNs in their self-regulated mode. The shutdown explains in detail the bright-end truncation of the blue sequence at ~L*, the appearance of luminous red-and-dead galaxies on the red sequence starting already at z~2, the colour bimodality, its strong dependence on environment density and its correlations with morphology and other galaxy properties. Before z~2-3, even haloes above the shock-heating mass form stars by cold streams penetrating through the hot gas. This explains the bright star-forming galaxies at z~3-4, the early appearance of massive galaxies on the red sequence, the high cosmological star-formation rate at high redshifts and the subsequent low rate at low redshifts., Comment: 17 pages, 12 figures, submitted to MNRAS

We present an approximation scheme to solve the Non Perturbative Renormalization Group equations and obtain the full momentum dependence of the $n$-point functions. It is based on an iterative procedure where, in a first step, an initial ansatz for the $n$-point functions is constructed by solving approximate flow equations derived from well motivated approximations. These approximations exploit the derivative expansion and the decoupling of high momentum modes. The method is applied to the O($N$) model. In leading order, the self energy is already accurate both in the perturbative and the scaling regimes. A stringent test is provided by the calculation of the shift $\Delta T_c$ in the transition temperature of the weakly repulsive Bose gas, a quantity which is particularly sensitive to all momentum scales. The leading order result is in agreement with lattice calculations, albeit with a theoretical uncertainty of about 25%., Comment: 48 pages, 15 figures A few minor corrections. A reference added

The Active Galactic Nuclei In Cosmological Simulations (AGNICS) project includes AGN into the GalICS hybrid (N-body + semi-analytic) model to explore the quasar-galaxy link in a cosmological perspective. The key problems are the quasar fuelling mechanism, the origin of the BH to bulge mass relation, the causal and chronological link between BH growth and galaxy formation, the properties of quasar hosts and the role of AGN feedback in galaxy formation. This first paper has two goals: to describe the general structure and the assumptions of the method, and to apply AGNICS to studying the joint formation of BHs and spheroids in galaxy mergers. We investigate under what conditions this scenario can reproduce the local BH distribution and the evolution of the quasar population. AGNICS contains two star formation modes: a quiescent one in discs and a starburst one in protospheroids, the latter triggered by mergers and disc instabilities. We assume that BH growth is linked to the starburst mode. The simplest version of this scenario, in which the BH accretion rate and the star formation rate in the starburst component are simply related by a constant of proportionality, does not to reproduce the cosmic evolution of the quasar population. A model in which this relation is modulated by a scaling with the square root of the starburst gas density can explain the evolution of the quasar luminosity function in B-band and X-rays (accounting for the obscuration inferred from X-ray studies). The scatter and the tilt that this model introduces in the BH-to-bulge mass relation are within the observational constraints. The quasar contribution grows with bolometric luminosity and for a given bulge mass the most massive BHs are in the bulges with the oldest stars., Comment: 20 pages, 11 figures, submitted to MNRAS on 1 Dec 2004, accepted

We argue that RHIC data, in particular those on the anisotropic flow coefficients v_2 and v_4, suggest that the matter produced in the early stages of nucleus-nucleus collisions is incompletely thermalized. We interpret the parameter (1/S)(dN/dy), where S is the transverse area of the collision zone and dN/dy the multiplicity density, as an indicator of the number of collisions per particle at the time when elliptic flow is established, and hence as a measure of the degree of equilibration. This number serves as a control parameter which can be varied experimentally by changing the system size, the centrality of the collision, or the beam energy. We provide predictions for Cu-Cu collisions at RHIC as well as for Pb-Pb collisions at the LHC., Comment: 5 pages, 2 figures, minor corrections

We present the first measurements of the Probability Distribution Function (PDF) of galaxy fluctuations in the VIMOS-VLT Deep Survey (VVDS) cone, covering 0.4x0.4 deg between 0.40.8; ii) the formation of bright galaxies is inhibited below a characteristic mass-overdensity threshold whose amplitude increases with redshift and luminosity; iii) the biasing function is non linear in all the redshift bins investigated with non-linear effects of the order of a few to 10% on scales >5Mpc., Comment: 30 pages, 17 figs, Accepted by A&A

We propose a method to solve the Non Perturbative Renormalization Group equations for the $n$-point functions. In leading order, it consists in solving the equations obtained by closing the infinite hierarchy of equations for the $n$-point functions. This is achieved: i) by exploiting the decoupling of modes and the analyticity of the $n$-point functions at small momenta: this allows us to neglect some momentum dependence of the vertices entering the flow equations; ii) by relating vertices at zero momenta to derivatives of lower order vertices with respect to a constant background field. Although the approximation is not controlled by a small parameter, its accuracy can be systematically improved. When it is applied to the O(N) model, its leading order is exact in the large $N$ limit; in this case, one recovers known results in a simple and direct way, i.e., without introducing an auxiliary field., Comment: Minor changes. Version to be published

We propose an effective theory which governs Pomeron dynamics in QCD at high energy, in the leading logarithmic approximation, and in the limit where N_c, the number of colors, is large. In spite of its remarkably simple structure, this effective theory generates precisely the evolution equations for scattering amplitudes that have been recently deduced from a more complete microscopic analysis. It accounts for the BFKL evolution of the Pomerons together with their interactions: dissociation (one Pomeron splitting into two) and recombination (two Pomerons merging into one). It is constructed by exploiting a duality principle relating the evolutions in the target and the projectile, more precisely, splitting and merging processes, or fluctuations in the dilute regime and saturation effects in the dense regime. The simplest Pomeron loop calculated with the effective theory is free of both ultraviolet or infrared singularities., Comment: 13 pages, 1 figure

We discuss proton-nucleus collisions in the framework of the color glass condensate. By assuming that the proton can be described as a low density color source, we solve exactly the Yang-Mills equations corresponding to this type of collision, and then use this solution in order to calculate inclusive gluon production or quark-antiquark production. Our result shows that k_T-factorization, while valid for gluon production, is violated for quark pair production in proton-nucleus collisions., Comment: Talk given at SEWM2004, Helsinki, June 2004

This paper presents the evolution of the clustering of the main population of galaxies from z=2.1 to z=0.2, from the first epoch VIMOS VLT Deep Survey (VVDS), a magnitude limited sample with 17.5<=I_{AB}<=24. We have computed the correlation functions \xi(r_p,\pi) and w_p(r_p), and the correlation length r_0(z), for the VVDS-02h and VVDS-CDFS fields, for a total of 7155 galaxies in a 0.61deg^2 area. We find that the correlation length in this sample stays roughly constant from z=0.5 to z=1.1, with r_0(z)=2.5-2.8 h^{-1} Mpc (comoving), for galaxies comparable in luminosity to the local 2dFGRS and SDSS samples, indicating that the amplitude of the correlation function was ~2.5x lower at z~1 than observed locally. The correlation length in our lowest redshift bin z=[0.2,0.5] is r_0=2.4 h^{-1} Mpc, lower than for any other population at the same redshift, indicating the low clustering of very low luminosity galaxies, 1.5 magnitudes fainter than in the 2dFGRS or SDSS. The correlation length is increasing to r_0~3.0 h^{-1} Mpc at higher redshifts z=[1.3,2.1], as we are observing increasingly brighter galaxies, comparable to galaxies with MB_AB=-20.5 locally. We compare our measurement to the DEEP2 measurements in the range z=[0.7,1.35] \citep{coil} on the population selected applying the same magnitude and color selection criteria as in their survey, and find comparable results. The slowly varying clustering of VVDS galaxies as redshift increases is markedly different from the predicted evolution of the clustering of dark matter, indicating that bright galaxies are already tracing the large scale structures emerging from the dark matter distribution 9-10 billion years ago, a supporting evidence for a strong evolution of the galaxy vs. dark matter bias., Comment: 11 pages, accepted 11-Apr-05 in A&A

We are presenting in this paper a detailed account of the methods used to compute the three-dimensional two-point galaxy correlation function in the VIMOS-VLT deep survey (VVDS). We investigate how instrumental selection effects and observational biases affect the measurements and identify the methods to correct them. We quantify the accuracy of our correction method using an ensemble of fifty mock galaxy surveys generated with the GalICS semi-analytic model of galaxy formation which incorporate the same selection biases and tiling strategy as the real data does. We demonstrate that we are able to recover the real-space two-point correlation function xi(s) to an accuracy better than 10% on scales larger than 1 h^{-1} Mpc, and of about 30% on scales below 1 h^{-1} Mpc, with the sampling strategy used for the first epoch VVDS data. The projected correlation function w_p(r_p) is recovered with an accuracy better than 10% on all scales 0.1 <= r <= 10 h^{-1} Mpc. There is a tendency for a small but systematic under-estimate of the correlation length derived from w_p(r_p) of 6% on average, remaining after our correction process. The large number of simulated surveys allows us to provide a reliable estimate of the cosmic error on the measurements of the correlation length r_0, of about 15-20% for the first epoch VVDS observation (Le Fevre et al.2004, astro-ph/0409133). The error estimation and measurement techniques outlined in this paper are being used in several studies which investigate in detail the clustering properties of galaxies in the VVDS data., Comment: 16 pages, accepted 10-Mar-05 in A&A