We present a new model for computing the spectral evolution of stellar populations at ages between 100,000 yr and 20 Gyr at a resolution of 3 A across the whole wavelength range from 3200 to 9500 A for a wide range of metallicities. These predictions are based on a newly available library of observed stellar spectra. We also compute the spectral evolution across a larger wavelength range, from 91 A to 160 micron, at lower resolution. The model incorporates recent progress in stellar evolution theory and an observationally motivated prescription for thermally-pulsing stars on the asymptotic giant branch. The latter is supported by observations of surface brightness fluctuations in nearby stellar populations. We show that this model reproduces well the observed optical and near-infrared colour-magnitude diagrams of Galactic star clusters of various ages and metallicities. Stochastic fluctuations in the numbers of stars in different evolutionary phases can account for the full range of observed integrated colours of star clusters in the Magellanic Clouds. The model reproduces in detail typical galaxy spectra from the Early Data Release (EDR) of the Sloan Digital Sky Survey (SDSS). We exemplify how this type of spectral fit can constrain physical parameters such as the star formation history, metallicity and dust content of galaxies. Our model is the first to enable accurate studies of absorption-line strengths in galaxies containing stars over the full range of ages. Using the highest-quality spectra of the SDSS EDR, we show that this model can reproduce simultaneously the observed strengths of those Lick indices that do not depend strongly on element abundance ratios [abridged].

https://hal.archives-ouvertes.fr/hal-00012948/document

Stellar population synthesis at the resolution of 2003

Statistics for Spatial Data.

We discuss the cosmological simulation code GADGET-2, a new massively parallel TreeSPH code, capable of following a collisionless fluid with the N-body method, and an ideal gas by means of smoothed particle hydrodynamics (SPH). Our implementation of SPH manifestly conserves energy and entropy in regions free of dissipation, while allowing for fully adaptive smoothing lengths. Gravitational forces are computed with a hierarchical multipole expansion, which can optionally be applied in the form of a TreePM algorithm, where only short-range forces are computed with the ‘tree’ method while long-range forces are determined with Fourier techniques. Time integration is based on a quasi-symplectic scheme where long-range and short-range forces can be integrated with different time-steps. Individual and adaptive short-range time-steps may also be employed. The domain decomposition used in the parallelization algorithm is based on a space-filling curve, resulting in high flexibility and tree force errors that do not depend on the way the domains are cut. The code is efficient in terms of memory consumption and required communication bandwidth. It has been used to compute the first cosmological N-body simulation with more than 10 10 dark matter particles, reaching a homogeneous spatial dynamic range of 10 5 per dimension in a three-dimensional box. It has also been used to carry out very large cosmological SPH simulations that account for radiative cooling and star formation, reaching total particle numbers of more than 250 million. We present the algorithms used by the code and discuss their accuracy and performance using a number of test problems. GADGET-2 is publicly released to the research community. Ke yw ords: methods: numerical ‐ galaxies: interactions ‐ dark matter.

/pdf/the-cosmological-simulation-code-gadget-2-48ba4eyb8k.pdf

The Cosmological simulation code GADGET-2

This paper describes the Seventh Data Release of the Sloan Digital Sky Survey (SDSS), marking the completion of the original goals of the SDSS and the end of the phase known as SDSS-II. It includes 11,663 deg^2 of imaging data, with most of the ~2000 deg^2 increment over the previous data release lying in regions of low Galactic latitude. The catalog contains five-band photometry for 357 million distinct objects. The survey also includes repeat photometry on a 120° long, 2°.5 wide stripe along the celestial equator in the Southern Galactic Cap, with some regions covered by as many as 90 individual imaging runs. We include a co-addition of the best of these data, going roughly 2 mag fainter than the main survey over 250 deg^2. The survey has completed spectroscopy over 9380 deg^2; the spectroscopy is now complete over a large contiguous area of the Northern Galactic Cap, closing the gap that was present in previous data releases. There are over 1.6 million spectra in total, including 930,000 galaxies, 120,000 quasars, and 460,000 stars. The data release includes improved stellar photometry at low Galactic latitude. The astrometry has all been recalibrated with the second version of the USNO CCD Astrograph Catalog, reducing the rms statistical errors at the bright end to 45 milliarcseconds per coordinate. We further quantify a systematic error in bright galaxy photometry due to poor sky determination; this problem is less severe than previously reported for the majority of galaxies. Finally, we describe a series of improvements to the spectroscopic reductions, including better flat fielding and improved wavelength calibration at the blue end, better processing of objects with extremely strong narrow emission lines, and an improved determination of stellar metallicities.

/pdf/the-seventh-data-release-of-the-sloan-digital-sky-survey-5cm9fqkek5.pdf

The Seventh Data Release of the Sloan Digital Sky Survey

The cold dark matter model has become the leading theoretical picture for the formation of structure in the Universe. This model, together with the theory of cosmic inflation, makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a simulation of the growth of dark matter structure using 2,1603 particles, following them from redshift z = 127 to the present in a cube-shaped region 2.230 billion lightyears on a side. In postprocessing, we also follow the formation and evolution of the galaxies and quasars. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with future generations of observational surveys of galaxies.

/pdf/simulations-of-the-formation-evolution-and-clustering-of-1yahvmv5qk.pdf

Simulations of the formation, evolution and clustering of galaxies and quasars

Using high resolution N-body simulations we address the problem of emptiness of giant 20 Mpc/h diameter voids found in the distribution of bright galaxies. Are the voids filled by dwarf galaxies? Do cosmological models predict too many small dark matter haloes inside the voids? Can the problems of cosmological models on small scales be addressed by studying the abundance of dwarf galaxies inside voids? We find that voids in the distribution of 10^12 Msun/h haloes (expected galactic magnitudes ~ M_*) are almost the same as the voids in 10^11 Msun/h haloes. Yet, much smaller haloes with masses 10^9 Msun/h and circular velocities v_circ about 20 km/s readily fill the voids: there should be almost 1000 of these haloes in a 20 Mpc/h void. A typical void of diameter 20 Mpc/h contains about 50 haloes with v_circ > 50 km/s. The haloes are arranged in a pattern, which looks like a miniature Universe: it has the same structural elements as the large-scale structure of the galactic distribution of the Universe. There are filaments and voids; larger haloes are at the intersections of filaments. The only difference is that all masses are four orders of magnitude smaller. There is severe (anti)bias in the distribution of haloes, which depends on halo mass and on the distance from the centre of the void. Large haloes are more antibiased and have a tendency to form close to void boundaries. The mass function of haloes in voids is different from the ``normal'' mass function. It is much steeper for high masses resulting in very few M33-type galaxies (v_circ about 100 km/s). We present an analytical approximation for the mass function of haloes in voids.

The structure of voids

Using high resolution N-body simulations of stellar disks embedded in cosmologically motivated dark matter halos, we study the evolution of bars and the transfer of angular momentum between halos and bars We find that dynamical friction results in some transfer of angular momentum to the halo, but the effect is much smaller than previously found in low resolution simulations and is incompatible with early analytical estimates After 5 Gyrs of evolution the stellar component loses only 5% – 7% of its initial angular momentum Mass and force resolutions are crucial for the modeling of bar dynamics In low resolution (300 – 500 pc) simulations we find that the bar slows down and angular momentum is lost relatively fast In simulations with millions of particles reaching a resolution of 20-40 pc, the pattern speed may not change over billions of years Our high resolution models produce bars which are fast rotators, where the ratio of the corotation radius to the bar major semi-axis lies in the range R = 12 17, marginally compatible with observational results In contrast to many previous simulations, we find that bars are relatively short As in many observed cases, the bar major semi-axis is close to the exponential length of the disk The transfer of angular momentum between inner and outer parts of the disk plays a very important role in the secular evolution of the disk and the bar The bar formation increases the exponential length of the disk by a factor of 12 -15 The transfer substantially increases the stellar mass in the center of the galaxy and decreases the dark matter-to-baryons ratio As the result, the central 2 kpc region is always strongly dominated by the baryonic component At intermediate (3 – 10 kpc) scales the disk is sub-dominant These models demonstrate that the efficiency of angular momentum transfer to the dark matter has been greatly overestimated More realistic models produce bar structure in striking agreement with observational results

Secular bar formation in galaxies with a significant amount of dark matter

We report results from high-resolution particle-mesh (PM) N-body simulations of structure formation in an $\Omega=1$ cosmological model with a mixture of Cold plus Hot Dark Matter (C+HDM) having $\Omega_{\rm cold}=0.6$, $\Omega_\nu=0.3$, and $\Omega_{\rm baryon}=0.1$. We present analytic fits to the C+HDM power spectra for both cold and hot ($\nu$) components, which provide initial conditions for our nonlinear simulations. In order to sample the neutrino velocities adequately, these simulations included six times as many neutrino particles as cold particles. Our simulation boxes were 14, 50, and 200~Mpc cubes (with $H_0=50$ km s$^{-1}$ Mpc$^{-1}$); we also did comparison simulations for Cold Dark Matter (CDM) in a 50~Mpc box. C+HDM with linear bias factor $b=1.5$ is consistent both with the COBE data and with the galaxy correlations we calculate. We find the number of halos as a function of mass and redshift in our simulations; our results for both CDM and C+HDM are well fit by a Press-Schechter model. The number density of galaxy-mass halos is smaller than for CDM, especially at redshift $z>2$, but the numbers of cluster-mass halos are comparable. We also find that on galaxy scales the neutrino velocities and flatter power spectrum in C+HDM result in galaxy pairwise velocities that are in good agreement with the data, and about 30\% smaller than in CDM with the same biasing factor. On scales of several tens of Mpc, the C+HDM streaming velocities are considerably larger than CDM. Thus C+HDM looks promising as a model of structure formation.

Structure Formation with Cold + Hot Dark Matter

We present the online MultiDark Database – a Virtual Observatory-oriented, relational database for hosting various cos-mological simulations. The data is accessible via an SQL (Structured Query Language) query interface, which also allows users to directly pose scientific questions, as shown in a number of examples in this paper. Further examples for the usage of the database are given in its extensive online documentation. The database is based on the same technology as the Millennium Database, a fact that will greatly facilitate the usage of both suites of cosmological simulations. The first release of the MultiDark Database hosts two 8.6 billion particle cosmological N-body simulations: the Bolshoi (250 h–1 Mpc simulation box, 1 h–1 kpc resolution) and MultiDark Run1 simulation (MDR1, or BigBolshoi, 1000 h–1 Mpc simulation box, 7 h–1 kpc resolution). The extraction methods for halos/subhalos from the raw simulation data, and how this data is structured in the database are explained in this paper. With the first data release, users get full access to halo/subhalo catalogs, various profiles of the halos at redshifts z = 0–15, and raw dark matter data for one time-step of the Bolshoi and four time-steps of the MultiDark simulation. Later releases will also include galaxy mock catalogs and additional merger trees for both simulations as well as new large volume simulations with high resolution. This project is further proof of the viability to store and present complex data using relational database technology. We encourage other simulators to publish their results in a similar manner. (© 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

The MultiDark Database: release of the Bolshoi and MultiDark cosmological simulations

We report results from high-resolution particle-mesh (PM) N-body simulations of structure formation in an Ω=1 cosmological model with a mixture of cold plus hot dark matter (C+HDM) having Ω cold =0.6, Ω ν =0.3, and Ω baryon =0.1. We present analytic fits to the C+HDM power spectra for both cold and hot (neutrino) components, which provide initial conditions for our nonlinear simulations. In order to sample the neutrino velocities adequately, these simulations included 6 times as many neutrino particles as cold particles. Our simulations boxes were 14, 50, and 200 Mpc cubes (with H 0 =50 km s −1 Mpc −1 ); we also did comparison simulations for cold dark matter (CDM) in a 50 Mpc box

Anatoly Klypin

Papers

The structure of voids

Secular bar formation in galaxies with a significant amount of dark matter

Structure Formation with Cold + Hot Dark Matter

The MultiDark Database: release of the Bolshoi and MultiDark cosmological simulations

Structure formation with cold plus hot dark matter