The thermal and chemical properties of the hot diffuse intragroup medium (IGrM) provide important constraints on the feedback processes associated with massive galaxy formation and evolution. Here we explore these constraints via a detailed analysis of the global properties of simulated z<3 galaxy groups from a cosmological simulation that includes a well-constrained prescription for stellar/supernovae-powered galactic outflows but no AGN feedback. Our aims are to (a) establish a baseline against which we will compare future models; (b) identify model successes due to stellar/supernovae-powered outflows; and (c) pinpoint features that signal the need for, and constrain the nature of, AGN feedback. Our simulation successfully reproduces key observed z=0 group IGrM properties, including the various X-ray Lx - Tx - entropy scaling relations, for all but the most massive groups. The z<1 redshift evolution of these also agree with observations. Contrary to expectations, the simulated groups' IGrM does not suffer catastrophic cooling. Yet, the z=0 group stellar mass is ~ 2X too large. This is due to the build-up of cold gas in the massive galaxies before they are incorporated inside groups. This in turn indicates that other feedback mechanisms must activate in real galaxies once their stellar masses grow to a few X 10^{10} M_sun. We show that these must be powerful enough to expel a significant fraction of the gas from the galactic halos. Gentle maintenance-mode (quenching) AGN feedback, as seen in galaxy clusters, will not do. Just as importantly, we find that the stellar/supernovae-powered winds are essential for understanding the IGrM metal abundances. Our simulation is able to reproduce the observed relationship between the global IGrM iron and silicon abundance and the group X-ray temperature, and these results ought to be relatively insensitive to the addition of AGN feedback.

The evolution of the baryon distribution in different phases, derived from cosmological simulations, are here reported. These computations indicate that presently most of baryons are in a warm-hot intergalactic (WHIM) medium (about 43%) while at z = 2.5 most of baryons constitute the diffuse medium (about 74%). Stars and the cold gas in galaxies represent only 14% of the baryons at z = 0. For z < 4 about a half of the metals are locked into stars while the fraction present in the WHIM and in the diffuse medium increases with a decreasing redshift. In the redshift range 0 < z < 2.5, the amount of metals in the WHIM increases from 4% to 22% while in the diffuse medium it increases from 0.6% to 4%. This enrichment process is due essentially to a turbulent diffusion mechanism associated to mass motions driven by supernova explosions. At z = 0, simulated blue (late type) galaxies show a correlation of the oxygen abundance present in the cold gas with the luminosity of the considered galaxy that agrees quite well with data derived from HII regions.

We perform simulations of feedback from supernovae and black holes with smoothed particle hydrodynamics (SPH). Such strong perturbations are inaccurately handled with standard time integration schemes, leading to poor energy conservation, a problem that is commonly overlooked. We show for the first time that, in the absence of radiative cooling, concordance of thermal and kinetic feedback are achieved when using an accurate time integration. In order to preserve the concordance of feedback methods when using a more efficient time integration scheme - as for instance the hierarchical time-step scheme - we implement a modified version of the time-step limiter proposed by Saitoh & Makino (2009). We apply the limiter to general test cases, and first show that this scheme violates energy conservation up to almost four orders of magnitude when energy is injected at random times. To tackle this issue, we find necessary, not only to ensure a fast information propagation, but also to enforce a prompt response of the system to the energy perturbation. The method proposed here to handle strong feedback events enables us to achieve energy conservation at percent level in all tests, even if all the available energy is injected into only one particle. We argue that concordance of feedback methods can be achieved in numerical simulations only if the time integration scheme preserve a high energy conservation level.

The coalescence history of massive black holes has been derived from cosmological simulations, in which the evolution of those objects and that of the host galaxies are followed in a consistent way. The present study indicates that supermassive black holes having masses greater than $\sim 10^{9} M_{\odot}$ underwent up to 500 merger events along their history. The derived coalescence rate per comoving volume and per mass interval permitted to obtain an estimate of the expected detection rate distribution of gravitational wave signals ("ring-down") along frequencies accessible by the planned interferometers either in space (LISA) or in the ground (Einstein). For LISA, in its original configuration, a total detection rate of about $15 yr^{-1}$ is predicted for events having a signal-to-noise ratio equal to 10, expected to occur mainly in the frequency range $4-9 mHz$. For the Einstein gravitational wave telescope, one event each 14 months down to one event each 4 years is expected with a signal-to-noise ratio of 5, occurring mainly in the frequency interval $10-20 Hz$. The detection of these gravitational signals and their distribution in frequency would be in the future an important tool able to discriminate among different scenarios explaining the origin of supermassive black holes.

One possible scenario for the formation of massive black holes (BHs) in the early Universe is from the direct collapse of primordial gas in atomic-cooling dark matter haloes in which the gas is unable to cool efficiently via molecular transitions. We study the formation of such BHs, as well as the accretion of gas onto these objects and the high energy radiation emitted in the accretion process, by carrying out cosmological radiation hydrodynamics simulations. In the absence of radiative feedback, we find an upper limit to the accretion rate onto the central object which forms from the initial collapse of hot (~ 10^4 K) gas of the order of 0.1 MSun per year. This is high enough for the formation of a supermassive star, the immediate precursor of a BH, with a mass of the order of 10^5 MSun. Assuming that a fraction of this mass goes into a BH, we track the subsequent accretion of gas onto the BH self-consistently with the high energy radiation emitted from the accretion disk. Using a ray-tracing algorithm to follow the propagation of ionizing radiation, we model in detail the evolution of the photoionized region which forms around the accreting BH. We find that BHs with masses of the order of 10^4 MSun initially accrete at close to the Eddington limit, but that the accretion rate drops to of order 1 percent of the Eddington limit after ~ 10^6 yr, due to the expansion of the gas near the BH in response to strong photoheating and radiation pressure. One signature of the accretion of gas onto BHs formed by direct collapse, as opposed to massive Pop III star formation, is an extremely high ratio of the luminosity emitted in He II 1640 to that emitted in H_alpha (or Ly_alpha); this could be detected by the James Webb Space Telescope. Finally, we briefly discuss implications for the coevolution of BHs and their host galaxies.

The correlations between the mass of supermassive black holes and properties of their host galaxies are investigated through cosmological simulations. Black holes grow from seeds of 100 solar masses inserted into density peaks present in the redshift range 12-15. Seeds grow essentially by accreting matter from a nuclear disk and also by coalescences resulting from merger episodes. At z=0, our simulations reproduce the black hole mass function and the correlations of the black hole mass both with stellar velocity dispersion and host dark halo mass. Moreover, the evolution of the black hole mass density derived from the present simulations agrees with that derived from the bolometric luminosity function of quasars, indicating that the average accretion history of seeds is adequately reproduced . However, our simulations are unable to form black holes with masses above $10^9 M_{\odot}$ at $z\sim 6$, whose existence is inferred from the bright quasars detected by the Sloan survey in this redshift range.

Cosmological N-body simulations were performed to study the evolution of the phase-space density Q = rho/sigma^3 of dark matter halos. No significant differences in the scale relations Q ~ sigma^(-2.1) or Q ~ M^(-0.82) are seen for "cold" or "warm" dark matter models. The follow up of individual halos from z = 10 up to the present time indicate the existence of two main evolutionary phases: an early and fast one (10 > z > 6.5), in which Q decreases on the average by a factor of 40 as a consequence of the randomization of bulk motions and a late and long one (6.5 > z > 0), in which Q decreases by a factor of 20 because of mixing induced by merger events. The study of these halos has also evidenced that rapid and positive variations of the velocity dispersion, induced by merger episods, are related to a fast decrease of the phase density Q.