We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Gaia Data Release 2 (Gaia DR2) contains results for 1693 million sources in the magnitude range 3 to 21 based on observations collected by the European Space Agency Gaia satellite during the first 22 months of its operational phase. We describe the input data, models, and processing used for the astrometric content of Gaia DR2, and the validation of these results performed within the astrometry task. Some 320 billion centroid positions from the pre-processed astrometric CCD observations were used to estimate the five astrometric parameters (positions, parallaxes, and proper motions) for 1332 million sources, and approximate positions at the reference epoch J2015.5 for an additional 361 million mostly faint sources. Special validation solutions were used to characterise the random and systematic errors in parallax and proper motion. For the sources with five-parameter astrometric solutions, the median uncertainty in parallax and position at the reference epoch J2015.5 is about 0.04 mas for bright (G<14 mag) sources, 0.1 mas at G=17 mag, and 0.7 mas at G=20 mag. In the proper motion components the corresponding uncertainties are 0.05, 0.2, and 1.2 mas/yr, respectively. The optical reference frame defined by Gaia DR2 is aligned with ICRS and is non-rotating with respect to the quasars to within 0.15 mas/yr. From the quasars and validation solutions we estimate that systematics in the parallaxes depending on position, magnitude, and colour are generally below 0.1 mas, but the parallaxes are on the whole too small by about 0.03 mas. Significant spatial correlations of up to 0.04 mas in parallax and 0.07 mas/yr in proper motion are seen on small (<1 deg) and intermediate (20 deg) angular scales. Important statistics and information for the users of the Gaia DR2 astrometry are given in the appendices.

Accurate tracking of probes is one of the key points of space exploration. Range and Doppler techniques are the most commonly used. In this paper we analyze the impact of the transponder delay, $i.e.$ the processing time between reception and re-emission of a two-way tracking link at the satellite, on tracking observables and on spacecraft orbits. We show that this term, only partially accounted for in the standard formulation of computed space observables, can actually be relevant for future missions with high nominal tracking accuracies or for the re-processing of old missions. We present several applications of our formulation to Earth flybys, the NASA GRAIL and the ESA BepiColombo missions.

A key objective of the ESA Gaia satellite is the realization of a quasi-inertial reference frame at visual wavelengths by means of global astrometric techniques. This requires an accurate mathematical and numerical modeling of relativistic light propagation, as well as double-blind-like procedures for the internal validation of the results, before they are released to the scientific community at large. Aim of this work is to specialize the Time Transfer Functions (TTF) formalism to the case of the Gaia observer and prove its applicability to the task of Global Sphere Reconstruction (GSR), in anticipation of its inclusion in the GSR system, already featuring the suite of RAMOD models, as an additional semi-external validation of the forthcoming Gaia baseline astrometric solutions. We extend the current GSR framework and software infrastructure (GSR2) to include TTF relativistic observation equations compatible with Gaia's operations. We use simulated data generated by the Gaia Data Reduction and Analysis Consortium (DPAC) to obtain different least-squares estimations of the full stellar spheres and gauge results. These are compared to analogous solutions obtained with the current RAMOD model in GSR2 and to the catalog generated with GREM, the model baselined for Gaia and used to generate the DPAC synthetic data. Linearized least-squares TTF solutions are based on spheres of about 132,000 primary stars uniformly distributed on the sky and simulated observations spanning the entire 5-yr range of Gaia's nominal operational lifetime. The statistical properties of the results compare well with those of GREM. Finally, comparisons to RAMOD@GSR2 solutions confirmed the known lower accuracy of that model and allowed us to establish firm limits on the quality of the linearization point outside of which an iteration for non-linearity is required for its proper convergence.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

The Einstein Equivalence Principle (EEP) is one of the foundations of the theory of General Relativity and several alternative theories of gravitation predict violations of the EEP. Experimental constraints on this fundamental principle of nature are therefore of paramount importance. The EEP can be split in three sub-principles: the Universality of Free Fall (UFF), the Local Lorentz Invariance (LLI) and the Local Position Invariance (LPI). In this paper we propose to use stable clocks in eccentric orbits to perform a test of the gravitational redshift, a consequence of the LPI. The best test to date was performed with the Gravity Probe A (GP-A) experiment in 1976 with an uncertainty of $1.4\times10^{-4}$. Our proposal considers the opportunity of using Galileo satellites 5 and 6 to improve on the GP-A test uncertainty. We show that considering realistic noise and systematic effects, and thanks to a highly eccentric orbit, it is possible to improve on the GP-A limit to an uncertainty around $(3-4)\times 10^{-5}$ after one year of integration of Galileo 5 and 6 data.

Determining range, Doppler and astrometric observables is of crucial interest for modelling and analyzing space observations. We recall how these observables can be computed when the travel time of a light ray is known as a function of the positions of the emitter and the receiver for a given instant of reception (or emission). For a long time, such a function--called a reception (or emission) time transfer function--has been almost exclusively calculated by integrating the null geodesic equations describing the light rays. However, other methods avoiding such an integration have been considerably developped in the last twelve years. We give a survey of the analytical results obtained with these new methods up to the third order in the gravitational constant $G$ for a mass monopole. We briefly discuss the case of quasi-conjunctions, where higher-order enhanced terms must be taken into account for correctly calculating the effects. We summarize the results obtained at the first order in $G$ when the multipole structure and the motion of an axisymmetric body is taken into account. We present some applications to on-going or future missions like Gaia and Juno. We give a short review of the recent works devoted to the numerical estimates of the time transfer functions and their derivatives.

Given the extreme accuracy of modern space science, a precise relativistic modeling of observations is required. We use the Time Transfer Functions formalism to study light propagation in the field of uniformly moving axisymmetric bodies, which extends the field of application of previous works. We first present a space-time metric adapted to describe the geometry of an ensemble of uniformly moving bodies. Then, we show that the expression of the Time Transfer Functions in the field of a uniformly moving body can be easily derived from its well-known expression in a stationary field by using a change of variables. We also give a general expression of the Time Transfer Function in the case of an ensemble of arbitrarily moving point masses. This result is given in the form of an integral easily computable numerically. We also provide the derivatives of the Time Transfer Function in this case, which are mandatory to compute Doppler and astrometric observables. We particularize our results in the case of moving axisymmetric bodies. Finally, we apply our results to study the different relativistic contributions to the range and Doppler tracking for the JUNO mission in the Jovian system.

During the last decade, the precision in the tracking of spacecraft has constantly improved. The discovery of few astrometric anomalies, such as the Pioneer and Earth flyby anomalies, stimulated further analysis of the operative modeling currently adopted in Deep Space Navigation (DSN). Our study shows that some traditional approximations lead to neglect tiny terms that could have consequences in the orbit determination of a probe in specific configurations such as during an Earth flyby. Therefore, we suggest here a way to improve the light time calculation used for probe tracking.

Given the extreme accuracy of modern space science, a precise relativistic modeling of observations is required. In particular, it is important to describe properly light propagation through the Solar System. For two decades, several modeling efforts based on the solution of the null geodesic equations have been proposed but they are mainly valid only for the first order Post-Newtonian approximation. However, with the increasing precision of ongoing space missions as Gaia, GAME, BepiColombo, JUNO or JUICE, we know that some corrections up to the second order have to be taken into account for future experiments. We present a procedure to compute the relativistic coordinate time delay, Doppler and astrometric observables avoiding the integration of the null geodesic equation. This is possible using the Time Transfer Function formalism, a powerful tool providing key quantities such as the time of flight of a light signal between two point-events and the tangent vector to its null-geodesic. Indeed we show how to compute the Time Transfer Functions and their derivatives (and thus range, Doppler and astrometric observables) up to the second post-Minkowskian order. We express these quantities as quadratures of some functions that depend only on the metric and its derivatives evaluated along a Minkowskian straight line. This method is particularly well adapted for numerical estimations. As an illustration, we provide explicit expressions in static and spherically symmetric space-time up to second post-Minkowskian order. Then we give the order of magnitude of these corrections for the range/Doppler on the BepiColombo mission and for astrometry in a GAME-like observation.

Given the extreme accuracy of modern space astrometry, a precise relativistic modeling of observations is required. Moreover, the availability of several models formulated in different and independent ways is a security against the presence of systematic errors in the analysis of future experimental results, like in the case of the Gaia mission. In this work, we simulate a series of observations using the two models to be used for the data analysis of Gaia, the Gaia RElativistic Model (GREM) and the Relativistic Astrometric MODel (RAMOD), and we compare them with the results of our astrometric model based on the Time Transfer Functions.

We review the mathematical models available for relativistic astrometry, discussing the different approaches and their accuracies in the context of the modern experiments from space like Gaia and GAME, and we show how these models can be applied to the real world, and their consequences from the mathematical and numerical point of view, with specific reference to the case of Gaia, whose launch is due before the end of the year.

Given the extreme accuracy reached by future global space astrometry, one needs a global relativistic modeling of observations. A relativistic definition of astrometric observables is then essential to find uniquely coordinates, parallax and proper motion of a stellar object. The standard procedure is to explicitly solve the null geodesic equations to describe the trajectory of a photon emitted by a celestial object and received by a moving observing satellite. However, this task can be avoided if one builds an astrometric set up using the recent formalism of time transfer functions.

In this work we investigate the cooling channels of diffuse gas (i.e. n_H<0.1 cm^-3) in cosmology. We aim to identify the wavelengths where most of the energy is radiated in the form of emission lines or continuum radiation, and the main elements and ions responsible for the emission. We use a subset of cosmological, hydrodynamical runs from the OWLS project to calculate the emission of diffuse gas and its evolution with time. We find that at z=0 (z=2) about 70 (80) per cent of the energy emitted by diffuse gas is carried by emission lines, with the continuum radiation contributing the remainder. Hydrogen lines in the Lyman series are the primary contributors to the line emission, with a share of 16 (20) per cent. Oxygen lines are the main metal contributors at high redshift, while silicon, carbon and iron lines are strongest at low redshift, when the contributions of AGB stars and supernova Ia explosions to the metal budget become important and when there is more hot gas. The ionic species carrying the most energy are OIII, CII, CIII, SiII, SiIII, FeII and SIII. The great majority of energy is emitted in the UV band (lambda=100-4000 A), both as continuum radiation and line emission. With almost no exception, all the strongest lines fall in this band. At high energies, continuum radiation is dominant (e.g., 80 per cent in the X-ray band), while lines contribute progressively more at lower energies. While the results do depend on the details of the numerical implementation of the physical processes modelled in the simulations, the comparison of results from different simulations demonstrates that the variations are overall small, and that the conclusions are fairly robust. Given the overwhelming importance of UV emission for the cooling of diffuse gas, it is desirable to build instruments dedicated to the detection and characterisation of diffuse UV emission.

A lot of fundamental tests of gravitational theories rely on highly precise measurements of the travel time and/or the frequency shift of electromagnetic signals propagating through the gravitational field of the Solar System. In practically all of the previous studies, the explicit expressions of such travel times and frequency shifts as predicted by various metric theories of gravity are derived from an integration of the null geodesic differential equations. However, the solution of the geodesic equations requires heavy calculations when one has to take into account the presence of mass multipoles in the gravitational field or the tidal effects due to the planetary motions, and the calculations become quite complicated in the post-post-Minkowskian approximation. This difficult task can be avoided using the time transfer function's formalism. We present here our last advances in the formulation of the one-way frequency shift using this formalism up to the post-post-Minkowskian approximation.

Rest-frame UV emission lines offer the possibility to directly image the gas around high-redshift galaxies with upcoming optical instruments. We use a suite of large, hydrodynamical simulations to predict the nature and detectability of emission lines from the intergalactic medium at 2<z<5. The brightest emission comes from HI Ly-alpha and the strongest metal line, CIII, is about an order of magnitude fainter, although HI Ly-alpha may be fainter if the gas is self-shielded to the UV background or if dust is important. The highest surface brightness regions for CIV, SiIII, SiIV and OVI are fainter than CIII by factors of a few. The NV and NeVIII lines, as well as HeII H-alpha, are substantially weaker but their maximum surface brightnesses still exceed 100 photon/cm^2/s/sr at z=2 (for 2" pixels). Lower ionisation lines arise in denser and colder gas that produces clumpier emission. The brightest HI Ly-alpha emission arises in highly overdense gas, but the highest surface brightness emission from high-ionisation metal lines traces a wider range of overdensities. Bright metal-line emission traces gas with temperatures close to the peak of the corresponding emissivity curve. While HI Ly-alpha, HeII H-alpha, CIII, SiIII, and SiIV are excellent probes of cold accretion flows and the colder parts of outflows, CIV, NV, OVI, and NeVIII are powerful tracers of the diffuse WHIM and galactic winds. A comparison of results from simulations with varying physical prescriptions demonstrates that the predictions for the brighter metal-line emission are robust to within factors of a few. Several emission lines from the high-redshift IGM will become detectable in the near future, possibly starting with the Cosmic Web Imager on Palomar. MUSE and the Keck Cosmic Web Imager have the potential to revolutionise studies of the interactions between high-redshift galaxies and their environment. (Abridged)

Approximately half the baryons in the local Universe are thought to reside in the warm-hot intergalactic medium (WHIM). Emission lines from metals in the UV band are excellent tracers of the cooler fraction of this gas. We present predictions for the surface brightness of a sample of UV lines that could potentially be observed by the next generation of UV telescopes at z<1. We use a subset of simulations from the OWLS project to create emission maps and to investigate the effect of varying the physical prescriptions for star formation, supernova and AGN feedback, chemodynamics and radiative cooling. Most models produce results in agreement within a factor of a few, indicating that the predictions are robust. Of the lines we consider, C III is the strongest line, but it typically traces gas colder than 10^5 K. The same is true for Si IV. The second strongest line, C IV, traces circum-galactic gas with T~10^5 K. O VI and Ne VIII probe the warmer (T~10^5.5 K and T~10^6 K, respectively) and more diffuse gas that may be a better tracer of the large scale structure. N V emission is intermediate between C IV and O VI. The intensity of all emission lines increases strongly with gas density and metallicity, and for the bright emission it is tightly correlated with the temperature for which the line emissivity is highest. In particular, the C III, C IV, Si IV and O VI emission that is sufficiently bright to be potentially detectable in the near future (>10^3 photon/s/cm^2/sr), comes from relatively dense (rho>10^2 rho_mean) and metal rich (Z>0.1 Z_sun) gas. As such, emission lines are highly biased tracers of the missing baryons and are not an optimal tool to close the baryon budget. However, they do provide a powerful means to detect the gas cooling onto or flowing out of galaxies and groups. (Abridged)

Emission lines from metals offer one of the most promising ways to detect the elusive warm-hot intergalactic medium (WHIM; 10^5 K<T<10^7 K), which is thought to contain a substantial fraction of the baryons in the low-redshift Universe. We present predictions for the soft X-ray line emission from the WHIM using a subset of cosmological simulations from the OverWhelmingly Large Simulations (OWLS) project. We use the OWLS models to test the dependence of the predicted emission on a range of physical prescriptions, such as cosmology, gas cooling and feedback from star formation and accreting black holes. Provided that metal-line cooling is taken into account, the models give surprisingly similar results, indicating that the predictions are robust. Soft X-ray lines trace the hotter part of the WHIM (T>10^6 K). We find that the OVIII 18.97A is the strongest emission line, with a predicted maximum surface brightness of ~10^2 photon/s/cm^2/sr, but a number of other lines are only slightly weaker. All lines show a strong correlation between the intensity of the observed flux and the density and metallicity of the gas responsible for the emission. On the other hand, the potentially detectable emission consistently corresponds to the temperature at which the emissivity of the electronic transition peaks. The emission traces neither the baryonic nor the metal mass. In particular, the emission that is potentially detectable with proposed missions, traces overdense (rho>10^2rho_mean) and metal-rich (Z>0.1Z_sun) gas in and around galaxies and groups. While soft X-ray line emission is therefore not a promising route to close the baryon budget, it does offer the exciting possibility to image the gas accreting onto and flowing out of galaxies.

We investigate the physics driving the cosmic star formation (SF) history using the more than fifty large, cosmological, hydrodynamical simulations that together comprise the OverWhelmingly Large Simulations (OWLS) project. We systematically vary the parameters of the model to determine which physical processes are dominant and which aspects of the model are robust. Generically, we find that SF is limited by the build-up of dark matter haloes at high redshift, reaches a broad maximum at intermediate redshift, then decreases as it is quenched by lower cooling rates in hotter and lower density gas, gas exhaustion, and self-regulated feedback from stars and black holes. The higher redshift SF is therefore mostly determined by the cosmological parameters and to a lesser extent by photo-heating from reionization. The location and height of the peak in the SF history, and the steepness of the decline towards the present, depend on the physics and implementation of stellar and black hole feedback. Mass loss from intermediate-mass stars and metal-line cooling both boost the SF rate at late times. Galaxies form stars in a self-regulated fashion at a rate controlled by the balance between, on the one hand, feedback from massive stars and black holes and, on the other hand, gas cooling and accretion. Paradoxically, the SF rate is highly insensitive to the assumed SF law. This can be understood in terms of self-regulation: if the SF efficiency is changed, then galaxies adjust their gas fractions so as to achieve the same rate of production of massive stars. Self-regulated feedback from accreting black holes is required to match the steep decline in the observed SF rate below redshift two, although more extreme feedback from SF, for example in the form of a top-heavy IMF at high gas pressures, can help.

We present a detailed analysis of predicted galaxy-galaxy merger fractions and rates in the Millennium simulation and compare these with the most up to date observations of the same quantities up to z~3. We carry out our analysis by considering the predicted merger history in the Millennium simulation within a given time interval, as a function of stellar mass. This method, as opposed to pair fraction counts, considers mergers that have already taken place, and allows a more direct comparison with the observed rates and fractions measured with the concentration-asymmetry-clumpiness (CAS) method. We examine the evolution of the predicted merger fraction and rate in the Millennium simulation for galaxies with stellar masses M_* ~ 10^9 - 10^12 M_sun. We find that the predicted merger rates and fractions match the observations well for galaxies with M_* > 10^11 M_sun at z<2, while significant discrepancies occur at lower stellar masses, and at z>2 for M_* > 10^11 M_sun systems. At z>2 the simulations underpredict the observed merger fractions by a factor of 4-10. The shape of the predicted merger fraction and rate evolutions are similar to the observations up to z~2, and peak at 1<z<2 in almost all mass bins. The exception is the merger rate of galaxies with M_* > 10^11 M_sun. We discuss possible reasons for these discrepancies, and compare different realisations of the Millennium simulation to understand the effect of varying the physical implementation of feedback. We conclude that the comparison is potentially affected by a number of issues, including uncertainties in interpreting the observations and simulations in terms of the assumed merger mass ratios and merger time-scales. (abridged)

In this paper we review the current predictions of numerical simulations for the origin and observability of the warm hot intergalactic medium (WHIM), the diffuse gas that contains up to 50 per cent of the baryons at z~0. During structure formation, gravitational accretion shocks emerging from collapsing regions gradually heat the intergalactic medium (IGM) to temperatures in the range T~10^5-10^7 K. The WHIM is predicted to radiate most of its energy in the ultraviolet (UV) and X-ray bands and to contribute a significant fraction of the soft X-ray background emission. While O VI and C IV absorption systems arising in the cooler fraction of the WHIM with T~10^5-10^5.5 K are seen in FUSE and HST observations, models agree that current X-ray telescopes such as Chandra and XMM-Newton do not have enough sensitivity to detect the hotter WHIM. However, future missions such as Constellation-X and XEUS might be able to detect both emission lines and absorption systems from highly ionised atoms such as O VII, O VIII and Fe XVII.

We present results for a galaxy formation model that includes a simple treatment for the disruption of dwarf galaxies by gravitational forces and galaxy encounters within galaxy clusters. This is implemented a posteriori in a semi-analytic model by considering the stability of cluster dark matter sub-haloes at z=0. We assume that a galaxy whose dark matter substructure has been disrupted will itself disperse, while its stars become part of the population of intracluster stars responsible for the observed intracluster light. Despite the simplicity of this assumption, our results show a substantial improvement over previous models and indicate that the inclusion of galaxy disruption is indeed a necessary ingredient of galaxy formation models. We find that galaxy disruption suppresses the number density of dwarf galaxies by about a factor of two. This makes the slope of the faint end of the galaxy luminosity function shallower, in agreement with observations. In particular, the abundance of faint, red galaxies is strongly suppressed. As a result, the luminosity function of red galaxies and the distinction between the red and the blue galaxy populations in colour-magnitude relationships are correctly predicted. Finally, we estimate a fraction of intracluster light comparable to that found in clusters of galaxies.

We present results of a new feedback scheme implemented in the Munich galaxy formation model. The new scheme includes a dynamical treatment of galactic winds powered by supernovae explosions and stellar winds in a cosmological context. We find that such a scheme is a good alternative to empirically-motivated recipes for feedback in galaxy formation. Model results are in good agreement with the observed luminosity functions and stellar mass function for galaxies in the local Universe. In particular, the new scheme predicts a number density of dwarfs that is lower than in previous models. This is a consequence of a new feature of the model, that allows an estimate of the amount of mass and metals that haloes can permanently deposit into the IGM. This loss of material leads to the suppression of star formation in small haloes and therefore to the decrease in the number density of dwarf galaxies. The model is able to reproduce the observed mass-stellar metallicity and luminosity-gas metallicity relationships. This demonstrates that our scheme provides a significant improvement in the treatment of the feedback in dwarf galaxies. Despite these successes, our model does not reproduce the observed bimodality in galaxy colours and predicts a larger number of bright galaxies than observed. Finally, we investigate the efficiency of metal injection in winds and in the intergalactic medium. We find that galaxies that reside in haloes with M_vir < 10^12 M_sun/h may deposit most of their metal mass into the intergalactic medium, while groups and clusters at z=0 have lost at most a few percent of their metals before the bulk of the halo mass was accreted.

The origin of intergalactic magnetic fields is still a mystery and several scenarios have been proposed so far: among them, primordial phase transitions, structure formation shocks and galactic outflows. In this work we investigate how efficiently galactic winds can provide an intense and widespread "seed" magnetisation. This may be used to explain the magnetic fields observed today in clusters of galaxies and in the intergalactic medium (IGM). We use semi-analytic simulations of magnetised galactic winds coupled to high resolution N-body simulations of structure formation to estimate lower and upper limits for the fraction of the IGM which can be magnetised up to a specified level. We find that galactic winds are able to seed a substantial fraction of the cosmic volume with magnetic fields. Most regions affected by winds have magnetic fields in the range -12 < Log B < -8 G, while higher seed fields can be obtained only rarely and in close proximity to wind-blowing galaxies. These seed fields are sufficiently intense for a moderately efficient turbulent dynamo to amplify them to the observed values. The volume filling factor of the magnetised regions strongly depends on the efficiency of winds to load mass from the ambient medium. However, winds never completely fill the whole Universe and pristine gas can be found in cosmic voids and regions unaffected by feedback even at z=0. This means that, in principle, there might be the possibility to probe the existence of primordial magnetic fields in such regions.

We investigate the effect of galactic winds on the Lyalpha forest in cosmological simulations of structure and galaxy formation. We combine high resolution N-body simulations of the evolution of the dark matter with a semi-analytic model for the formation and evolution of galaxies which includes detailed prescriptions for the long-term evolution of galactic winds. This model is the first to describe the evolution of outflows as a two-phase process (an adiabatic bubble followed by a momentum--driven shell) and to include metal--dependent cooling of the outflowing material. We find that the main statistical properties of the Lyalpha forest, namely the flux power spectrum P(k) and the flux probability distribution function (PDF), are not significantly affected by winds and so do not significantly constrain wind models. Winds around galaxies do, however, produce detectable signatures in the forest, in particular, increased flux transmissivity inside hot bubbles, and narrow, saturated absorption lines caused by dense cooled shells. We find that the Lyalpha flux transmissivity is highly enhanced near strongly wind-blowing galaxies, almost half of all high-redshift galaxies in our sample, in agreement with the results of Adelberger et al. (2005). Finally, we propose a new method to identify absorption lines potentially due to wind shells in the Lyalpha forest: we calculate the abundance of saturated regions in spectra as a function of region width and we find that the number with widths smaller than about 1 Angstrom at z=3 and 0.6 Angstrom at z=2 may be more than doubled. This should be detectable in real spectra.

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.