F. R. Bouchet

galics– I. A hybrid N-body/semi-analytic model of hierarchical galaxy formation

This is the first paper of a series that describes the methods and basic results of the galics model (Galaxies In Cosmological Simulations). galics is a hybrid model for hierarchical galaxy formation studies, combining the outputs of large cosmological N-body simulations with simple, semi-analytic recipes to describe the fate of the baryons within dark matter haloes. The simulations produce a detailed merging tree for the dark matter haloes, including complete knowledge of the statistical properties arising from the gravitational forces. We intend to predict the overall statistical properties of galaxies, with special emphasis on the panchromatic spectral energy distribution emitted by galaxies in the ultraviolet/optical and infrared/submillimetre wavelength ranges. In this paper, we outline the physically motivated assumptions and key free parameters that go into the model, comparing and contrasting with other parallel efforts. We specifically illustrate the success of the model in comparison with several data sets, showing how it is able to predict the galaxy disc sizes, colours, luminosity functions from the ultraviolet to far infrared, the Tully–Fisher and Faber–Jackson relations, and the fundamental plane in the local Universe. We also identify certain areas where the model fails, or where the assumptions needed to succeed are at odds with observations, and pay special attention to understanding the effects of the finite resolution of the simulations on the predictions made. Other papers in this series will take advantage of different data sets available in the literature to extend the study of the limitations and predictive power of galics, with particular emphasis put on high-redshift galaxies.

The structure and dynamical evolution of dark matter haloes

We use N-body simulations to investigate the structure and dynamical evolution of dark matter halos in clusters of galaxies. Our sample consists of nine massive halos from an Einstein-De Sitter universe with scale free power spectrum and spectral index n = 1. Halos are resolved by 20000 particles each, on average, and have a dynamical resolution of 20-25 kpc, as shown by extensive tests. Large scale tidal fields are included up to a scale L = 150 Mpc using background particles. We find that the halo formation process can be characterized by the alternation of two dynamical configurations: a merging phase and a relaxation phase, defined by their signature on the evolution of the total mass and root mean square (rms) velocity. Halos spend on average one third of their evolution in the merging phase and two thirds in the relaxation phase. Using this definition, we study the density profiles and show how they change during the halo dynamical history. In particular, we find that the average density profiles of our halos are fitted by the Navarro, Frenk & White (1995) analytical model with an rms residual of 17% between the virial radius Rv and 0.01Rv. The Hernquist (1990) analytical density profiles fits the same halos with an rms residual of 26%. The trend with mass of the scale radius of these fits is marginally consistent with that found by Cole & Lacey (1996): compared to their results our halos are more centrally concentrated, and the relation between scale radius and halo mass is slightly steeper. We find a moderately large scatter in this relation, due both to dynamical evolution within halos and to fluctuations in the halo population. We analyze the dynamical equilibrium of our halos using the Jeans’ equation, and find that on average they are approximately in equilibrium within their virial radius. Finally, we find that the projected mass profiles of our simulated halos are in very good agreement with the profiles of three rich galaxy clusters derived from strong and weak gravitational lensing observations.

Semi-analytic modelling of galaxy evolution in the IR/submm range

This paper proposes a new semi-analytic modelling of galaxy properties in the IR/submm wavelength range, which is explicitly set in a cosmological framework. We start from a description of the non-dissipative and dissipative collapses of primordial perturbations, and add star formation, stellar evolution and feedback, as well as the absorption of starlight by dust and its re-emission in the IR and submm. This type of approach has had some success in reproducing the optical properties of galaxies. We hereafter propose a simple extension to the IR/submm range. The growth of structures is followed according to the standard cold dark matter model. We assume that star formation proceeds either in a ‘quiescent’ mode, e.g., as in discs, or in a ‘burst’ mode with 10 times shorter time-scales. In order to reproduce the current data on the evolution of the comoving cosmic star formation rate and gas densities, we need to introduce a mass fraction involved in the ‘burst’ mode strongly increasing with redshift, probably reflecting the increase of interaction and merging activity. We estimate the IR/submm luminosities of these ‘mild starburst’ and ‘luminous UV/IR galaxies’, and we explore how much star formation could be hidden in heavily extinguished, ‘ultraluminous IR galaxies’ by designing a family of evolutionary scenarios which are consistent with the current status of the ‘cosmic constraints’, as well as with the IRAS 60-μm luminosity function and faint counts, but with different high-z IR luminosity densities. However, these scenarios generate a cosmic infrared background whose spectrum falls within the ±1σ range of the isotropic IR component detected by Puget et al. and revisited by Guiderdoni et al. We give predictions for the faint galaxy counts and redshift distributions at IR and submm wavelengths. The submm range is very sensitive to the details of the evolutionary scenarios. As a result, the ongoing and forthcoming observations with ISO and SCUBA (and later with SIRTF, SOFIA, FIRST and PLANCK) will put strong constraints on the evolution of galaxies at z∼1 and beyond.

Dancing in the dark: galactic properties trace spin swings along the cosmic web

A large-scale hydrodynamical cosmological simulation, Horizon-AGN , is used to investigate the alignment between the spin of galaxies and the large-scale cosmic filaments above redshift one. The analysis of more than 150 000 galaxies with morphological diversity in a 100h −1 Mpc comoving box size shows that the spin of low-mass, rotationdominated, blue, star-forming galaxies is preferentially aligned with their neighbouring filaments. High-mass, dispersion-dominated, red, quiescent galaxies tend to have a spin perpendicular to nearby filaments. The reorientation of the spin of massive galaxies is provided by galaxy mergers which are significant in the mass build up of high-mass galaxies. We find that the stellar mass transition from alignment to misalignment happens around 3×10 10 M⊙. This is consistent with earlier findings of a dark matter mass transition for the orientation of the spin of halos (5 × 10 11 M⊙ at the same redshift from Codis et al. 2012). With these numerical evidence, we advocate a scenario in which galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the physical and morphological properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as Active Galactic Nuclei is mandatory to quench in situ star formation in massive galaxies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping galaxy spins misaligned with cosmic filaments. It also promotes diversity amongst galaxy properties.

Planck pre-launch status: The HFI instrument, from specification to actual performance

Context. The High Frequency Instrument (HFI) is one of the two focal instruments of the Planck mission. It will observe the whole sky in six bands in the 100 GHz-1 THz range. Aims: The HFI instrument is designed to measure the cosmic microwave background (CMB) with a sensitivity limited only by fundamental sources: the photon noise of the CMB itself and the residuals left after the removal of foregrounds. The two high frequency bands will provide full maps of the submillimetre sky, featuring mainly extended and point source foregrounds. Systematic effects must be kept at negligible levels or accurately monitored so that the signal can be corrected. This paper describes the HFI design and its characteristics deduced from ground tests and calibration. Methods: The HFI instrumental concept and architecture are feasible only by pushing new techniques to their extreme capabilities, mainly: (i) bolometers working at 100 mK and absorbing the radiation in grids; (ii) a dilution cooler providing 100 mK in microgravity conditions; (iii) a new type of AC biased readout electronics and (iv) optical channels using devices inspired from radio and infrared techniques. Results: The Planck-HFI instrument performance exceeds requirements for sensitivity and control of systematic effects. During ground-based calibration and tests, it was measured at instrument and system levels to be close to or better than the goal specification.

Cosmological constraints from Archeops

We analyze the cosmological constraints that Archeops places on adiabatic cold dark matter models with passive power-law initial fluctuations. Because its angular power spectrum has small bins in l and large l coverage down to COBE scales, Archeops provides a precise determination of the first acoustic peak in terms of position at multipole l_peak=220 +- 6, height and width. An analysis of Archeops data in combination with other CMB datasets constrains the baryon content of the Universe, Omega(b)h^2 = 0.022 (+0.003,-0.004), compatible with Big-Bang nucleosynthesis and with a similar accuracy. Using cosmological priors obtainedfrom recent non-CMB data leads to yet tighter constraints on the total density, e.g. Omega(tot)=1.00 (+0.03,-0.02) using the HST determination of the Hubble constant. An excellent absolute calibration consistency is found between Archeops and other CMB experiments, as well as with the previously quoted best fit model.The spectral index n is measured to be 1.04 (+0.10,-0.12) when the optical depth to reionization, tau, is allowed to vary as a free parameter, and 0.96 (+0.03,-0.04) when tau is fixed to zero, both in good agreement with inflation.

The cosmic microwave background anisotropy power spectrum measured by archeops

We present a determination by the Archeops experiment of the angular power spectrum of the cosmic microwave background anisotropy in 16 bins over the multipole range l=15-350. Archeops was conceived as a precursor of the Planck HFI instrument by using the same optical design and the same technology for the detectors and their cooling. Archeops is a balloon-borne instrument consisting of a 1.5 m aperture diameter telescope and an array of 21 photometers maintained at ~100 mK that are operating in 4 frequency bands centered at 143, 217, 353 and 545 GHz. The data were taken during the Arctic night of February 7, 2002 after the instrument was launched by CNES from Esrange base (Sweden). The entire data cover ~ 30% of the sky.This first analysis was obtained with a small subset of the dataset using the most sensitive photometer in each CMB band (143 and 217 GHz) and 12.6% of the sky at galactic latitudes above 30 degrees where the foreground contamination is measured to be negligible. The large sky coverage and medium resolution (better than 15 arcminutes) provide for the first time a high signal-to-noise ratio determination of the power spectrum over angular scales that include both the first acoustic peak and scales probed by COBE/DMR. With a binning of Delta(l)=7 to 25 the error bars are dominated by sample variance for l below 200. A companion paper details the cosmological implications.

Skewness induced by gravity

Large-scale structures, observed today, are generally believed to have grown from random, small-amplitude inhomogeneities, present in the early universe. We investigate how gravitational instability drives the distribution of these fluctuations away from the initial state, assumed to be Gaussian. Using second-order perturbation theory, we calculate the skewness factor, S 3 ≡ / 2 . Here the angle brackets, , denote an ensemble average, and δ is the density contrast field, smoothed with a low-pass spatial filter. We show that S 3 decreases with the slope of the fluctuation power spectrum; it depends only weakly on Ω, the cosmological density parameter

The pre-launch Planck Sky Model: a model of sky emission at submillimetre to centimetre wavelengths

We present the Planck Sky Model (PSM), a parametric model for generating all-sky, few arcminute resolution maps of sky emission at submillimetre to centimetre wavelengths, in both intensity and polarisation. Several options are implemented to model the cosmic microwave background, Galactic diffuse emission (synchrotron, free-free, thermal and spinning dust, CO lines), Galactic HII regions, extragalactic radio sources, dusty galaxies, and thermal and kinetic Sunyaev-Zeldovich signals from clusters of galaxies. Each component is simulated by means of educated interpolations/extrapolations of data sets available at the time of the launch of the Planck mission, complemented by state-of-the-art models of the emission. Distinctive features of the simulations are spatially varying spectral properties of synchrotron and dust; different spectral parameters for each point source; modelling of the clustering properties of extragalactic sources and of the power spectrum of fluctuations in the cosmic infrared background. The PSM enables the production of random realisations of the sky emission, constrained to match observational data within their uncertainties. It is implemented in a software package that is regularly updated with incoming information from observations. The model is expected to serve as a useful tool for optimising planned microwave and sub-millimetre surveys and testing data processing and analysis pipelines. It is, in particular, used to develop and validate data analysis pipelines within the Planck collaboration. A version of the software that can be used for simulating the observations for a variety of experiments is made available on a dedicated website.

Weakly nonlinear gravitational instability for arbitrary Omega

The weakly nonlinear evolution of particle trajectories in Friedmann-Lemaitre models with zero cosmological constant is investigated. The matter is assumed to be a nonrelativistic pressureless fluid. A perturbative expansion in Lagrangian coordinates is used and analytic expressions for the second-order solutions with arbitary density parameter Ω are derived. This perturbative expansion is valid provided the gradients of the displacement field are small, a much weaker condition than the usual Eulerian requirement of smallness of the density perturbations