ALMA products are stored in the Science Archive in the form of FITS images. It is a common idea that the FITS image headers should collect in their keywords all the information that an archive User might want to search for in order to quickly select, compare, or discard datasets. With this perspective in mind, we first present a short description of the current status of the ALMA FITS archive and images. We realized that at the moment most of the parameters that could be useful for a general User are still missing in the archived data. We then provide a CASA task generating the image header keywords that we suggest to be relevant for the scientific exploitation of the ALMA archival data. The proposed tool could be also applied to several types of interferometer data and part of it is implemented in a web interface. An example of the scientific application of the keywords is also discussed.

Since its inception in the early 2000, the Virtual Observatory (VO), developed as a collaboration of many national and international projects, has become a major factor in the discovery and dissemination of astronomical information worldwide. The IVOA has been coordinating all these efforts worldwide to ensure a common VO framework that enables transparent access to and interoperability of astronomy resources (data and software) around the world. The VO is not a magic solution to all astronomy data management challenges but it does bring useful solutions in many areas borne out by the fact that VO interfaces are broadly found in astronomy major data centres and projects worldwide. Astronomy data centres have been building VO services on top of their existing data services to increase interoperability with other VO-compliant data resources to take advantage of the continuous and increasing development of VO applications. VO applications have made multi-instrument and multi-wavelength science, a difficult and fruitful part of astronomy, somewhat easier. More recently, several major new astronomy projects have been directly adopting VO standards to build their data management infrastructure, giving birth to VO built-in archives. Embracing the VO framework from the beginning brings the double gain of not needing to reinvent the wheel and ensuring from the start interoperability with other astronomy VO resources. Some of the IVOA standards are also starting to be used by neighbour disciplines like planetary sciences. There is still quite a lot to be done on the VO, in particular tackling the upcoming big data challenge and how to find interoperable solutions to the new data analysis paradigm of bringing and running the software close to the data.

This document presents HiPS, a hierarchical scheme for the description, storage and access of sky survey data. The system is based on hierarchical tiling of sky regions at finer and finer spatial resolution which facilitates a progressive view of a survey, and supports multi-resolution zooming and panning. HiPS uses the HEALPix tessellation of the sky as the basis for the scheme and is implemented as a simple file structure with a direct indexing scheme that leads to practical implementations.

The Advanced Telescope for High ENergy Astrophysics (Athena) is the X-ray observatory mission selected by ESA within its Cosmic Vision 2015-2025 programme to address the Hot and Energetic Universe scientific theme. The ESO-Athena Synergy Team (EAST) has been tasked to single out the potential scientific synergies between Athena and optical/near-infrared (NIR) and sub/mm ground based facilities, in particular those of ESO (i.e., the VLT and ELT, ALMA and APEX), by producing a White Paper to identify and develop the: 1. needs to access ESO ground-based facilities to achieve the formulated Athena science objectives; 2. needs to access Athena to achieve the formulated science objectives of ESO facilities contemporary to Athena; 3. science areas where the synergetic use of Athena and ESO facilities in the late 2020s will result in scientific added value. Community input to the process happened primarily via a dedicated ESO - Athena Synergy Workshop that took place on Sept. 14 - 16, 2016 at ESO, Garching. This White Paper presents the results of the EAST's work, sorted by synergy area, and deals with the following topics: 1. the Hot Universe: Early groups and clusters and their evolution, Physics of the Intracluster medium, Missing baryons in cosmic filaments; 2. the Energetic Universe: Supermassive black hole (SMBH) history, SMBH accretion disks, Active Galactic Nuclei feedback - Molecular outflows, Ultra-fast outflows, Accretion Physics, Transient Science; 3. Observatory Science: Star Formation, Stars. It then discusses the optical-NIR-sub-mm perspective by providing details on VLT/MOONS, the E-ELT instruments, in particular the MOS, VISTA/4MOST, the ESO and ALMA archives, future ALMA and ESO developments, and finally the (likely) ESO - Athena astronomical scene in the 2020s. (abridged)

With high over-subscription rates and significant operational costs, observatories must ensure that their operations are efficient and effective. A number of key performance indicators are generally used to evaluate the observatory's performance among which are the numbers of publications and citations of refereed journal articles to measure the overall scientific impact. Those measures, however, are broad and can not assess whether the observatory was successful on a project-by-project basis to deliver data to the PIs enabling them to carry out their science and to publish their results. In particular the reasons that prevented PIs from publishing remain hidden. Understanding and acting upon those reasons, however, have the potential to substantially improve the observatory's operational model. Of course not every approved project even should lead to a publication. Indeed, the risk of not finding the expected (or any unexpected) science in the data the PI receives is an inherent and indispensable part of the scientific process. But even here, measuring the fraction of such projects can lead to valuable insights which might then be used to instruct future proposal review committees. To fully close the loop on the end-to-end data-flow, ALMA has started in March 2015 to send survey questions to PIs where two years after the end of the proprietary period no publication making use of the delivered data could be identified. We describe our method as well as the type of conclusions we hope to be able to draw once a statistically relevant sample of answers has been received.

The increasing volumes of astronomical data require practical methods for data exploration, access and visualisation. The Hierarchical Progressive Survey (HiPS) is a HEALPix based scheme that enables a multi-resolution approach to astronomy data from the individual pixels up to the whole sky. We highlight the decisions and approaches that have been taken to make this scheme a practical solution for managing large volumes of heterogeneous data. Early implementors of this system have formed a network of HiPS nodes, with some 250 diverse data sets currently available, with multiple mirror implementations for important data sets. This hierarchical approach can be adapted to expose Big Data in different ways. We describe how the ease of implementation, and local customisation of the Aladin Lite embeddable HiPS visualiser have been keys for promoting collaboration on HiPS.

The scientific impact of a facility is the most important measure of its success. Monitoring and analysing the scientific return can help to modify and optimise operations and adapt to the changing needs of scientific research. The methodology that we have developed to monitor the scientific productivity of the ALMA Observatory, as well as the first results, are described. We focus on the outcome of the first cycle (Cycle 0) of ALMA Early Science operations. Despite the fact that only two years have passed since the completion of Cycle 0 and operations have already changed substantially, this analysisconfirms the effectiveness of the underlying concepts. We find that ALMA is fulfilling its promise as a transformational facility for the observation of the Universe in the submillimetre.

The Atacama Large Millimeter/submillimeter Array (ALMA), an international partnership of Europe, North America and East Asia in cooperation with the Republic of Chile, is the largest astronomical project in existence. While ALMA's capabilities are ramping up, Early Science observations have started. The ALMA Archive is at the center of the operations of the telescope array and is designed to manage the 200 TB of data that will be taken each year, once the observatory is in full operations. We briefly describe design principles. The second part of this paper focuses on how astronomy is likely to evolve as the amount and complexity of data taken grows. We argue that in the future observatories will compete for astronomers to work with their data, that observatories will have to reorient themselves to from providing good data only to providing an excellent end-to-end user-experience with all its implications, that science-grade data-reduction pipelines will become an integral part of the design of a new observatory or instrument and that all this evolution will have a deep impact on how astronomers will do science. We show how ALMA's design principles are in line with this paradigm.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

Telescope bibliographies are integral parts of observing facilities. They are used to associate the published literature with archived observational data, to measure an observatory's scientific output through publication and citation statistics, and to define guidelines for future observing strategies. The ESO and NRAO librarians as well as NAOJ jointly maintain the ALMA (Atacama Large Millimeter/submillimeter Array) bibliography, a database of refereed papers that use ALMA data. In this paper, we illustrate how relevant articles are identified, which procedures are used to tag entries in the database and link them to the correct observations, and how results are communicated to ALMA stakeholders and the wider community. Efforts made to streamline the process will be explained and evaluated, and a first analysis of ALMA papers published after two years of observations will be given.

A public release of slitless spectra, obtained with ACS/WFC and the G800L grism, is presented. Spectra were automatically extracted in a uniform way from 153 archival fields (or "associations") distributed across the two Galactic caps, covering all observations to 2008. The ACS G800L grism provides a wavelength range of 0.55-1.00 \mu$m, with a dispersion of $40 \ \AA / pixel$ and a resolution of $\sim 80\ \AA$ for point-like sources. The ACS G800L images and matched direct images were reduced with an automatic pipeline that handles all steps from archive retrieval, alignment and astrometric calibration, direct image combination, catalogue generation, spectral extraction and collection of metadata. The large number of extracted spectra (73,581) demanded automatic methods for quality control and an automated classification algorithm was trained on the visual inspection of several thousand spectra. The final sample of quality controlled spectra includes 47,919 datasets (65% of the total number of extracted spectra) for $32,149$ unique objects, with a median $i_{\rm AB}$-band magnitude of 23.7, reaching 26.5 AB for the faintest objects. Each released dataset contains science-ready 1D and 2D spectra, as well as multi-band image cutouts of corresponding sources and a useful preview page summarising the direct and slitless data, astrometric and photometric parameters. In order to characterize the slitless spectra, emission-line flux and equivalent width sensitivity of the ACS data were compared with public ground-based spectra in the GOODS-South field. An example list of emission line galaxies with two or more identified lines is also included, covering the redshift range $0.2-4.6$.

We use a hybrid approach that combines high-resolution simulations of the formation of a Milky Way-like halo with a semi-analytic model of galaxy formation to study the mass content of dwarf galaxies in the concordance $\Lambda$CDM cosmology. We find that the mass within 600 pc of dark matter haloes hosting luminous satellites has a median value of $\sim 3.2\times 10^7$ Msun with very little object-to-object scatter. In contrast, the present day total luminosities of the model satellites span nearly five orders of magnitude. These findings are in very good agreement with the results recently reported in the literature for the dwarf spheroidal galaxies of the Milky Way. In our model, dwarf irregular galaxies like the Small Magellanic Cloud, are predicted to have similar or slightly larger dark matter mass within 600 pc.

The Hubble Legacy Archive (HLA) aims to create calibrated science data from the Hubble Space Telescope archive and make them accessible via user-friendly and Virtual Observatory (VO) compatible interfaces. It is a collaboration between the Space Telescope Science Institute (STScI), the Canadian Astronomy Data Centre (CADC) and the Space Telescope - European Coordinating Facility (ST-ECF). Data produced by the Hubble Space Telescope (HST) instruments with slitless spectroscopy modes are among the most difficult to extract and exploit. As part of the HLA project, the ST-ECF aims to provide calibrated spectra for objects observed with these HST slitless modes. In this paper, we present the HLA NICMOS G141 grism spectra. We describe in detail the calibration, data reduction and spectrum extraction methods used to produce the extracted spectra. The quality of the extracted spectra and associated direct images is demonstrated through comparison with near-IR imaging catalogues and existing near-IR spectroscopy. The output data products and their associated metadata are publicly available through a web form at http://hla.stecf.org and via VO interfaces. In total, 2470 spectra of 1923 unique targets are included in the current release.

We present the transverse flux correlation function of the Lyman-alpha forest in quasar absorption spectra at z~2.1 from VLT-FORS and VLT-UVES observations of a total of 32 pairs of quasars; 26 pairs with separations in the range 0.6<theta<4 arcmin and 6 pairs with 4<theta<10 arcmin. Correlation is detected at the 3sigma level up to separations of the order of ~4 arcmin (or ~4.4 h^{-1} Mpc comoving at z=2.1 for Omega_{m}=0.3 and Omega_{Lambda}=0.7). We have, furthermore, measured the longitudinal correlation function at a somewhat higher mean redshift (z=2.39) from 20 lines of sight observed with high spectral resolution and high signal-to-noise ratio with VLT-UVES. We compare the observed transverse and longitudinal correlation functions to that obtained from numerical simulations and illustrate the effect of spectral resolution, thermal broadening and peculiar motions. The shape and correlation length of the correlation functions are in good agreement with those expected from absorption by the filamentary and sheet-like structures in the photoionized warm intergalactic medium predicted in CDM-like models for structures formation. Using a sample of 139 C_{iv} systems detected along the lines of sight toward the pairs of quasars we also investigate the transverse correlation of metals on the same scales. The observed transverse correlation function of intervening C_{iv} absorption systems is consistent with that of a randomly distributed population of absorbers. This is likely due to the small number of pairs with separation less than 2 arcmin. We detect, however, a significant overdensity of systems in the sightlines towards the quartet Q0103-294A&B, Q0102-2931 and Q0102-293 which extends over the redshift range 1.5<z<2.2 and an angular scale larger than 10 arcmin.

The kinematical properties of elliptical galaxies formed during the mergers of equal mass, stars+gas+dark matter spiral galaxies are compared to the observed low velocity dispersions found for planetary nebulae on the outskirts of ellipticals, which have been interpreted as pointing to a lack of dark matter in ellipticals (which poses a problem for the standard model of galaxy formation). We find that the velocity dispersion profiles of the stars in the simulated ellipticals match well the observed ones. The low outer stellar velocity dispersions are mainly caused by the radial orbits of the outermost stars, which, for a given binding energy must have low angular momentum to reach their large radial distances, usually driven out along tidal tails.

Elements of kinematical and dynamical modeling of elliptical galaxies (Es) are presented. In projection, NFW models resemble Sersic models, but with a very narrow range of shapes (m=3+/-1). The total density profile of Es cannot be NFW-like because the predicted local M/L and aperture velocity dispersion within an effective radius (Re) are much lower than observed. Stars must then dominate Es out to a few Re. Fitting an NFW model to the total density profile of Sersic+NFW (stars+dark matter [DM]) Es results in very high concentration parameters, as found by X-ray observers. Kinematical modeling of Es assuming an isotropic NFW DM model underestimates M/L at the virial radius by a factor of 1.6 to 2.4, because dissipationless LCDM halos have slightly different density profiles and slightly radial velocity anisotropy. In N-body+gas simulations of Es as merger remnants of spirals embedded in DM halos, the slope of the DM density profile is steeper when the initial spiral galaxies are gas-rich. The Hansen & Moore (2006) relation between anisotropy and the slope of the density profile breaks down for gas and DM, but the stars follow an analogous relation with slightly less radial anisotropies for a given density slope. Using kurtosis (h_4) to infer anisotropy in Es is dangerous, as h_4 is also sensitive to small levels of rotation. The stationary Jeans equation provides accurate masses out to 8 Re. The discrepancy between the modeling of Romanowsky et al. (2003), indicating a dearth of DM in Es, and the simulations analyzed by Dekel et al. (2005), which match the spectroscopic observations of Es, is partly due to radial anisotropy and to observing oblate Es face-on. However, one of the 15 solutions to the orbit modeling of Romanowsky et al. is found to have an amount and concentration of DM consistent with LCDM predictions.

We use a high-resolution simulation of a galaxy-sized dark matter halo, published simulated data as well as four cluster-sized haloes from Fukushige, Kawai & Makino to study the inner halo structure in a Lambda cold dark matter cosmology. We find that the circular velocity curves are substantially better described by Stoehr et al. (SWTS) profiles than by Navarro, Frenk & White (NFW) or Moore et al. profiles. Our findings confirm that no asymptotic slope is reached and that the profiles are nearly universal, but not perfectly. The velocity profiles curve at a_constant_ rate in log(r) over the full converged range in radii and the corresponding extrapolated density profiles reach a finite maximum density. We find, that the claim of a strong discrepancy between the observed inner slopes of the density profiles of low surface brightness galaxies (LSB) and their simulated counterparts on the grounds of currently available observations and simulations is unfounded. In addition, if the SWTS profile turns out to be a good description of the halo profile for the regions that cannot be probed with simulations of today, then even in these regions the agreeement between simulations and observations is very reasonable.

One of the most debated issues in the theoretical modeling of cosmic reionization is the impact of small-mass gravitationally-bound structures. We carry out the first numerical investigation of the role of such sterile `minihaloes', which serve as self-shielding screens of ionizing photons. Minihaloes are too small to be properly resolved in current large-scale cosmological simulations, and thus we estimate their effects using a sub-grid model, considering two cases that bracket their effect within this framework. In the `extreme suppression' case in which minihalo formation ceases once a region is partially ionized, their effect on cosmic reionization is modest, reducing the volume-averaged ionization fraction by an overall factor of less than 15%. In the other extreme, in which minihalo formation is never suppressed, they delay complete reionization as much as Delta z~2, in rough agreement with the results from a previous semi-analytical study by the authors. Thus, depending on the details of the minihalo formation process, their effect on the overall progress of reionization can range from modest to significant, but the minihalo photon consumption is by itself insufficient to force an extended reionization epoch.

We derive the transverse flux correlation function in the Lyman-alpha forest at z~2.1 from VLT-FORS observations of a total of 32 pairs of quasars. The shape and correlation length of the transverse correlation function are in good agreement with the paradigm of intergalactic medium predicted in CDM-like models for structures formation. Using a sample of 139 Civ systems detected along the lines of sight toward the pairs of quasars we investigate the transverse correlation of metals on the same scales. We find that the correlation function is consistent with that of a randomly distributed population of Civ systems. However, we detect an important overdensity of systems in front of a quartet.

While kinematical modelling of the low PN velocity dispersions observed in the outer regions of elliptical galaxies suggest a lack of dark matter around these galaxies, we report on an analysis of a suite of $N$-body simulations (with gas) of major mergers of spiral galaxies embedded in dark matter halos, and find that the outer velocity dispersions are as low as observed for the PNe. The inconsistency between our dynamical modelling and previous kinematical modelling is caused by very radial stellar orbits and projection effects when viewing face-on oblate ellipticals. Our simulations (weakly) suggest the youth of PNe around ellipticals, and we propose that the universality of the PN luminosity function may be explained if the bright PNe in ellipticals are formed after the regular accretion of very low mass gas-rich galaxies.

The cosmic star formation histories are evaluated for different minimum masses of the initial halo structures, with allowance for realistic gas outflows. With a minimum halo mass of 10^{7} - 10^{8} M_odot and a moderate outflow efficiency, we reproduce both the current baryon fraction and the early chemical enrichment of the IGM. The intensity of the formation rate of ``normal'' stars is also well constrained by the observations: it has to be dominated by star formation in elliptical galaxies, except perhaps at very low redshift. The fraction of baryons in stars is predicted as are also the type Ia and II supernova event rates. Comparison with SN observations in the redshift range z=0-2 allows us to set strong constraints on the time delay of type Ia supernovae (a total delay of \sim 4 Gyr is required to fit the data), the lower end of the mass range of the progenitors (2 - 8 M_odot) and the fraction of white dwarfs that reproduce the type Ia supernova (about 1 per cent). The intensity in the initial starburst of zero metallicity stars below 270 M_\odot must be limited in order to avoid premature overenrichment of the IGM. Only about 10 - 20 % of the metals present in the IGM at z = 0 have been produced by population III stars at very high z. The remaining 80 - 90 % are ejected later by galaxies forming normal stars, with a maximum outflow efficiency occurring at a redshift of about 5. We conclude that 10^{-3} of the mass in baryons must lie in first massive stars in order to produce enough ionizing photons to allow early reionization of the IGM by z \sim 15.

There is strong evidence that the mass in the Universe is dominated by dark matter, which exerts gravitational attraction but whose exact nature is unknown. In particular, all galaxies are believed to be embedded in massive haloes of dark matter. This view has recently been challenged by surprisingly low random stellar velocities in the outskirts of ordinary elliptical galaxies, which were interpreted as indicating a lack of dark matter (Mendez et al. 2001; Romanowsky et al. 2003). Here we show that the low velocities are in fact compatible with galaxy formation in dark-matter haloes. Using numerical simulations of disc-galaxy mergers, we find that the stellar orbits in the outer regions of the resulting ellipticals are very elongated. These stars were torn by tidal forces from their original galaxies during the first close passage and put on outgoing trajectories. The elongated orbits, combined with the steeply falling density profile of the observed tracers, explain the observed low velocities even in the presence of large amounts of dark matter. Projection effects when viewing a triaxial elliptical can lead to even lower observed velocities along certain lines of sight.

We present a semi-analytic treatment of galactic winds within high resolution, large scale cosmological N-body simulations of a LCDM Universe. The evolution of winds is investigated by following the expansion of supernova driven superbubbles around the several hundred thousand galaxies that form in an approximately spherical region of space with diameter 52 Mpc/h and mean density close to the mean density of the Universe. We focus our attention on the impact of winds on the diffuse intergalactic medium. Initial conditions for mass loss at the base of winds are taken from Shu, Mo and Mao (2003). Results are presented for the volume filling factor and the mass fraction of the IGM affected by winds and their dependence on the model parameters is carefully investigated. The mass loading efficiency of bubbles is a key factor to determine the evolution of winds and their global impact on the IGM: the higher the mass loading, the later the IGM is enriched with metals. Galaxies with 10^9 < M_* < 10^10 M_sun are responsible for most of the metals ejected into the IGM at z=3, while galaxies with M_* < 10^9 M_sun give a non negligible contribution only at higher redshifts, when larger galaxies have not yet assembled. We find a higher mean IGM metallicity than Lyalpha forest observations suggest and we argue that the discrepancy may be explained by the high temperatures of a large fraction of the metals in winds, which may not leave detectable imprints in absorption in the Lyalpha forest.

We use a large cosmological simulation and a variety of high resolution resimulations of individual cluster and galaxy haloes to study the systematics of subhalo populations over ranges of 1000 in halo mass and 1000 in the ratio of subhalo to parent halo mass. The subhalo populations of different haloes are not scaled copies of each other, but vary systematically with halo properties. On average, the amount of substructure increases with halo mass. At fixed mass, it decreases with halo concentration and with halo formation redshift. These trends are comparable in size to the scatter in subhalo abundance between similar haloes. Averaged over all haloes of given mass, the abundance of low mass subhaloes per unit parent halo mass is independendent of parent mass.It is very similar to the abundance per unit mass of low mass haloes in the universe as a whole, once differing boundary definitions for subhaloes and haloes are accounted for. The radial distribution of subhaloes within their parent haloes is substantially less centrally concentrated than that of the dark matter. It varies at most weakly with the mass(or concentration) of the parent halo and not at all with subhalo mass. It does depend on the criteria used to define the subhalo population considered. About 90% of present-day subhaloes were accreted after z=1 and about 70% after z=0.5. Only about 8% the total mass of all haloes accreted at z=1 survives as bound subhaloes at z=0. For haloes accreted at z=2, the survival mass fraction is just 2%. Subhaloes seen near the centre of their parent typically were accreted earlier and retain less of their original mass than those seen near the edge. These strong systematics mean that comparison with galaxies in real clusters is only possible if the formation of the luminous component is modelled appropriately.

We study the effect of galactic outflows on the Lya forest. The winds are modelled as fully ionised spherical bubbles centered around galactic haloes. The observed flux probability distribution and flux power spectrum limit the volume filling factor of bubbles to be less than 10%. We have compared the mean flux as a function of distance from haloes with the Adelberger etal (ASSP) measurement. A model with bubbles of constant radius of 1.5Mpc/h surrounding the most massive haloes matches the high transmissivity at separations <0.5Mpc/h, but the increase of the transmissivity at small scales is more gradual than observed. At separations <0.5Mpc/h, our estimate of the cosmic variance error is $\Delta F\sim 0.3$, 30% higher than that of ASSP. The difficulty in matching the rise in the transmissivity at small separations is caused by residual absorption of neutral hydrogen lying physically outside the bubbles but having a redshift position similar to the haloes. The flux level is thus sensitive to the amplitude of the coherent velocity shear near haloes which increases with halo mass in the simulation. A model where LBGs are starbursts in small mass haloes matches the observations with smaller bubble size than one where more massive haloes host the LBGs. A starburst model with a bubble radius of 1Mpc/h and a volume filling factor of 2% is consistent with the ASSP measurements at the 1-1.5$\sigma$ level. If this model is correct the sharp rise of the transmissivity at small separation in the ASSP sample is due to cosmic variance and is expected to become more moderate for a larger sample.

If the dark matter in the Universe is made of weakly self-interacting particles, they may self-annihilate and emit gamma-rays. We use high resolution numerical simulations to estimate directly the annihilation flux from the central regions of the Milky Way and from dark matter substructures in its halo. Although such estimates remain uncertain because of their strong dependence on the structure of the densest regions, our numerical experiments suggest that less direct calculations have overestimated the emission both from the centre and from halo substructure. We estimate a maximal enhancement of at most a factor of a few with respect to a smooth spherical halo of standard Navarro-Frenk-White (NFW) structure. We discuss detection strategies for the next generation of gamma-ray detectors and find that the annihilation flux may be detectable, regardless of uncertainties about the densest regions, for the annihilation cross-sections predicted by currently popular elementary particle models for the dark matter.

We analyse the properties of substructures within dark matter halos (subhalos) using a set of high-resolution numerical simulations of the formation of structure in a Lambda-CDM Universe. Our simulation set includes 11 high-resolution simulations of massive clusters as well as a region of mean density, allowing us to study the spatial and mass distribution of substructures down to a mass resolution limit of 10^9 h^(-1)Mo. We also investigate how the properties of substructures vary as a function of the mass of the `parent' halo in which they are located. We find that the substructure mass function depends at most weakly on the mass of the parent halo and is well described by a power-law. The radial number density profiles of substructures are steeper in low mass halos than in high mass halos. More massive substructures tend to avoid the centres of halos and are preferentially located in the external regions of their parent halos. We also study the mass accretion and merging histories of substructures, which we find to be largely independent of environment. We find that a significant fraction of the substructures residing in clusters at the present day were accreted at redshifts z < 1. This implies that a significant fraction of present-day `passive' cluster galaxies should have been still outside the cluster progenitor and more active at z~1.

This paper illustrates how mock observational samples of high-redshift galaxies with sophisticated selection criteria can be extracted from the predictions of GALICS, a hybrid model of hierarchical galaxy formation that couples the outputs of large cosmological simulations and semi-analytic recipes to describe dark matter collapse and the physics of baryons respectively. As an example of this method, we focus on the properties of Lyman Break Galaxies at redshift 3. With the MOMAF software package described in a companion paper, we generate a mock observational sample with selection criteria as similar as possible to those implied in the actual observations of z = 3 LBGs by Steidel et al.(1995). Our model predictions are in good agreement with the observed number density and 2D correlation function. We investigate the optical/IR luminosity budget as well as several other physical properties of LBGs and find them to be in general agreement with observed values. Looking into the future of these LBGs we predict that 75% of them end up as massive ellipticals today, even though only 35% of all our local ellipticals are predicted to have a LBG progenitor. In spite of some shortcomings, this new 'mock observation' method clearly represents a necessary first step toward a more accurate comparison between hierarchical models of galaxy formation and real observational surveys.

We have studied the IGM reionization process in its full cosmological context including structure evolution and a realistic galaxy population. We have used a combination of high-resolution N-body simulations (to describe the dark matter and diffuse gas component), a semi-analytic model of galaxy formation (to track the evolution of the sources of ionizing radiation) and the Monte Carlo radiative transfer code CRASH (to follow the propagation of ionizing photons into the IGM). The process has been followed in the largest volume ever used for this kind of study, a field region of the universe with a comoving length of L~20/h Mpc, embedded in a much larger cosmological simulation. To assess the effect of environment on the reionization process, the same radiative transfer simulations have been performed on a 10/h Mpc comoving box, centered on a clustered region. We find that, to account for the all ionizing radiation, objects with total masses of M~10^9 Msun must be resolved. In this case, the simulated stellar population produces a volume averaged ionization fraction x_v=0.999 by z~8, consistent with observations without requiring any additional sources of ionization. We also find that environment substantially affects the reionization process. In fact, although the simulated proto-cluster occupies a smaller volume and produces a higher number of ionizing photons, it gets totally ionized later. This is because high density regions, which are more common in the proto-cluster, are difficult to ionize because of their high recombination rates.

We compare the structure and kinematics of the 11 known satellites of the Milky Way with high resolution simulations of the formation of its dark halo in a LambdaCDM universe. In contrast to earlier work, we find excellent agreement. The observed kinematics are exactly those predicted for stellar populations with the observed spatial structure orbiting within the most massive ``satellite'' substructures in our simulations. Less massive substructures have weaker potential wells than those hosting the observed satellites. If there is a halo substructure ``problem'', it consists in understanding why halo substructures have been so inefficient in making stars. Suggested modifications of dark matter properties (for example, self-interacting or warm dark matter) may well spoil the good agreement found for standard Cold Dark Matter.

We compare two techniques for following the cooling of gas and its condensation into galaxies within high resolution simulations of cosmologically representative regions. Both techniques treat the dark matter using N-body methods. One follows the gas using smoothed particle hydrodynamics (SPH) while the other uses simplified recipes from semi-analytic (SA) models. We compare the masses and locations predicted for dense knots of cold gas (the `galaxies') when the two techniques are applied to evolution from the same initial conditions and when the additional complications of star formation and feedback are ignored. We find that above the effective resolution limit of the two techniques, they give very similar results both for global quantities such as the total amount of cooled gas and for the properties of individual `galaxies'. The SA technique has systematic uncertainties arising from the simplified cooling model adopted, while details of the SPH implementation can produce substantial systematic variations in the galaxy masses it predicts. Nevertheless, for the best current SPH methods and the standard assumptions of the SA model, systematic differences between the two techniques are remarkably small. The SA technique gives adequate predictions for the condensation of gas into `galaxies' at less than one percent of the computational cost of obtaining similar results at comparable resolution using SPH.