With a total integration time of 168 hours and a narrowband (NB) filter tuned to Lyman alpha at z = 8.8, the UltraVISTA survey has set out to find some of the most distant galaxies, on the verge of the Epoch of Reionization. Previous calculations of the expected number of detected Lya-emitting galaxies (LAEs) at this redshift did not explicitly take into account the radiative transfer (RT) of Lya. In this work we combine a theoretical model for the halo mass function with numerical results from high-res cosmological hydrosimulations with LyC+Lya RT, assessing the visibility of LAEs residing in these halos. Uncertainties such as cosmic variance and the anisotropic escape of Lya are taken into account, and it is predicted that once the survey has finished, the probabilities of detecting none, one, or more than one are ~90%, ~10%, and ~1%; a significantly smaller success rate compared to earlier predictions, due to the combined effect of a highly neutral IGM scattering Lya to such large distances from the galaxy that they fall outside the observational aperture, and to the actual depth of the survey being less than predicted. Because the IGM affects NB and broadband (BB) magnitudes differently, we argue for a relaxed color selection criterion of NB - BB ~ +0.85. But since the flux is continuum-dominated, even if a galaxy is detectable in the NB its probability of being selected as a NB excess object is <~35%. Various properties of galaxies at this redshift are predicted, e.g. UV and Lya LFs, M*-Mh relation, spectral shape, optimal aperture, and the anisotropic escape of Lya through both a dusty ISM and a partly neutral IGM. Finally, we describe and publish a fast numerical code for adding numbers with asymmetric uncertainties ("x_{-sigma_1}^{+sigma_2}") proving to be significantly better than the standard, but wrong, way of adding upper and lower uncertainties in quadrature separately.

We study the evolution of the Red Sequence (RS) scatter in galaxy clusters and groups simultaneously using predictions from our simulations (cosmological hydrodynamic + semi-analytical) as well as observational data from the HAWK-I Cluster Survey (HCS), a sample of galaxy clusters at redshifts $0.8 < z < 1.5$. We analyze the intrinsic scatter of the RS to investigate whether the stellar age can be its main effective driver, at the same time assessing the role of metallicity variations in shaping the RS building at around epoch $z\sim$1 and beyond. To this purpose we rely on various methods to derive the average age and age spread from the RS colour scatter of the HCS sample, with the aid of population synthesis models. The RS scatter predicted by the models at $z< 0.7$ is found to not depend on the star formation history adopted, whilst at $z>the correlation found between age and rest-frame colour scatters is quite robust, although all age scatter estimations ultimately depend on the definition of RS as well as on the completeness limits adopted. We find that the age spread of RS galaxies predicted by both hydrodynamical simulations and SAM increases with cosmic epoch, while the ratio between the age spread and the average age remains approximately constant. Both trends are in agreement with observational results from both the HCS and other literature samples.

The quasar Q0918+1636 (z=3.07) has an intervening high-metallicity Damped Lyman-alpha Absorber (DLA) along the line of sight, at a redshift of z=2.58. The DLA is located at a large impact parameter of 16.2 kpc, and has an almost solar metallicity. It is shown, that a novel type of cosmological galaxy formation models, invoking a new SNII feedback prescription, the Haardt & Madau (2012) UVB field and explicit treatment of UVB self-shielding, can reproduce the observed characteristics of the DLA. UV radiation from young stellar populations in the galaxy, in particular in the photon energy range 10.36-13.61 eV (relating to Sulfur II abundance), are also considered in the analysis. It is found that a) for L~L* galaxies (at z=2.58), about 10% of the sight-lines through the galaxies at impact parameter 16.2 kpc will display a Sulfur II column density N(SII)$>$ 10$^{15.82}$ cm$^{-2}$ (the observed value for the DLA), and b) considering only cases where a near-solar metallicity will be detected at 16.2 kpc impact parameter, the probability distribution of galaxy SFR peaks near the value observed for the DLA galaxy counterpart of ~27 Msun/yr. It is argued, that the bulk of the alpha-elements, like Sulfur, traced by the high metal column density, b=16.2 kpc absorption lines, were produced by evolving young stars in the inner galaxy, and later transported outward by galactic winds.

We present SIGAME (SImulator of GAlaxy Millimetre/submillimetre Emission), a new numerical code designed to simulate the 12CO rotational line emission spectrum of galaxies. Using sub-grid physics recipes to post-process the outputs of smoothed particle hydrodynamics (SPH) simulations, a molecular gas phase is condensed out of the hot and partly ionized SPH gas. The gas is subjected to far-UV radiation fields and cosmic ray ionization rates which are set to scale with the local star formation rate volume density. Level populations and radiative transport of the CO lines are solved with the 3-D radiative transfer code LIME. We have applied SIGAME to cosmological SPH simulations of three disc galaxies at z=2 with stellar masses in the range ~(0.5-2)x10^11 Msun and star formation rates ~40-140 Msun/yr. Global CO luminosities and line ratios are in agreement with observations of disc galaxies at z~2 up to and including J=3-2 but falling short of the few existing J=5-4 observations. The central 5 kpc regions of our galaxies have CO 3-2/1-0 and 7-6/1-0 brightness temperature ratios of ~0.55-0.65 and ~0.02-0.08, respectively, while further out in the disc the ratios drop to more quiescent values of ~0.5 and <0.01. Global CO-to-H2 conversion (alpha_CO) factors are ~=1.5 Msun*pc^2/(K km s/1), i.e. ~2-3 times below typically adopted values for disc galaxies, and alpha_CO increases with radius, in agreement with observations of nearby galaxies. Adopting a top-heavy Giant Molecular Cloud (GMC) mass spectrum does not significantly change the results. Steepening the GMC density profile leads to higher global line ratios for J_up>=3 and CO-to-H2 conversion factors [~=3.6 Msun*pc^2/(K km/s)].

By means of our own cosmological-hydrodynamical simulation and semi-analytical model we studied galaxy population properties in clusters and groups, spanning over 10 different bands from UV to NIR, and their evolution since redshift z=2. We compare our results in terms of galaxy red/blue fractions and luminous-to-faint ratio (LFR) on the Red Sequence (RS) with recent observational data reaching beyond z=1.5. Different selection criteria were tested in order to retrieve galaxies belonging to the RS: either by their quiescence degree measured from their specific SFR ("Dead Sequence"), or by their position in a colour-colour plane which is also a function of sSFR. In both cases, the colour cut and the limiting magnitude threshold were let evolving with redshift, in order to follow the natural shift of the characteristic luminosity in the LF. We find that the Butcher-Oemler effect is wavelength-dependent, with the fraction of blue galaxies increasing steeper in optical colours than in NIR. Besides, only when applying a lower limit in terms of fixed absolute magnitude, a steep BO effect can be reproduced, while the blue fraction results less evolving when selecting samples by stellar mass or an evolving magnitude limit. We then find that also the RS-LFR behaviour, highly debated in the literature, is strongly dependent on the galaxy selection function: in particular its very mild evolution recovered when measured in terms of stellar mass, is in agreement with values reported for some of the highest redshift confirmed (proto)clusters. As to differences through environments, we find that normal groups and (to a lesser extent) cluster outskirts present the highest values of both star forming fraction and LFR at low z, while fossil groups and cluster cores the lowest: this separation among groups begins after z~0.5, while earlier all group star forming properties are undistinguishable.

We present results from SPH-cosmological simulations, including self-consistent modelling of SN feedback and chemical evolution, of galaxies belonging to two clusters and twelve groups. We reproduce the mass-metallicity (ZM) relation of galaxies classified in two samples according to their star-forming activity, as parametrized by their sSFR, across a redshift range up to z=2. Its slope shows irrelevant evolution in the passive sample, being steeper in groups than in clusters. However, the sub-sample of high-mass passive galaxies only is characterized by a steep increase of the slope with redshift, from which it can be inferred that the bulk of the slope evolution of the ZM relation is driven by the more massive passive objects. (...ABRIDGED...) The ZM relation for the star-forming sample reveals an increasing scatter with redshift, indicating that it is still being built at early epochs. The star-forming galaxies make up a tight sequence in the SFR-M_* plane at high redshift, whose scatter increases with time alongside with the consolidation of the passive sequence. We also confirm the anti-correlation between sSFR and stellar mass, pointing at a key role of the former in determining the galaxy downsizing, as the most significant means of diagnostics of the star formation efficiency. Likewise, an anti-correlation between sSFR and metallicity can be established for the star-forming galaxies, while on the contrary more active galaxies in terms of simple SFR are also metal-richer. We discuss these results in terms of the mechanisms driving the evolution within the high- and low-mass regimes at different epochs: mergers, feedback-driven outflows and the intrinsic variation of the star formation efficiency.

We perform hydrodynamical simulations of a young galactic disc embedded in a hot gaseous halo using parameters typical for Lyman break galaxies (LBGs). We take into account the (static) gravitational potentials due to a dark matter halo, a stellar bulge and a disc of stars and gas. Star formation is treated by a local Kennicutt-Schmidt law. We simplify the structure of the interstellar medium by restricting the computational domain to a 25th of the full azimuthal angle, effectively assuming large-scale axisymmetry and neglecting any effects of spiral structure, and focus on the large-scale ISM drivers, the superbubbles. Supernovae are triggered randomly and have preset event sizes of several tens to hundreds. We further investigate different halo gas pressures and energy injection methods. Many of our simulated galaxies, but not all, develop bipolar outflows. We characterise the strength of the outflow by mass and energy outflow rates, and investigate the effect of changes to the details of the model. We find that supernovae are more effective if comprised into larger superbubbles. The weight and the pressure of the halo gas is able to quench galactic outflows. The wind emerges from a series of superbubbles in regions where a critical star formation density is exceeded. The superbubbles expand into the gaseous halo at slightly supersonic speed, producing radiative shock waves with similar characteristics as the absorptions systems observed around LBGs.

We use numerical simulations to study the effects of the patchiness of a partly reionized intergalactic medium (IGM) on the observability of Ly-alpha emitters (LAEs) at high redshifts (z ~ 6). We present a new model that divides the Ly-alpha radiative transfer into a (circum-)galactic and an extragalactic (IGM) part, and investigate how the choice of intrinsic line model affects the IGM transmission results. We use our model to study the impact of neutral hydrogen on statistical observables such as the Ly-alpha restframe equivalent width (REW) distribution, the LAE luminosity function and the two-point correlation function. We find that if the observed changes in LAE luminosity functions and equivalent width distributions between z ~ 6 and z ~ 7 are to be explained by an increased IGM neutral fraction alone, we require an extremely late and rapid reionization scenario, where the Universe was ~ 40 % ionized at z = 7, ~ 50 % ionized at z = 6.5 and ~ 100 % ionized at z = 6. This is in conflict with other observations, suggesting that intrinsic LAE evolution at z > 6 cannot be completely neglected. We show how the two-point correlation function can provide more robust constraints once future observations obtain larger LAE samples, and provide predictions for the sample sizes needed to tell different reionization scenarios apart.

We study the intergalactic transmission of radiation in the vicinity of the Ly{\alpha} wavelength. Simulating sightlines through the intergalactic medium (IGM) in detailed cosmological hydrosimulations, the impact of the IGM on the shape of the line profile from Ly{\alpha} emitting galaxies at redshifts 2.5 to 6.5 is investigated. In particular we show that taking into account the correlation of the density and velocity fields of the IGM with the galaxies, the blue part of the spectrum may be appreciably reduced, even at relatively low redshifts. This may in some cases provide an alternative to the often-invoked outflow scenario, although it is concluded that this model is still a plausible explanation of the many asymmetric Ly{\alpha} profiles observed. Applying the calculated wavelength dependent transmission to simulated spectra from Ly{\alpha} emitting galaxies, we derive the fraction of photons that are lost in the IGM, in addition to what is absorbed internally in the galaxies due to dust. Moreover, by comparing the calculated transmission of radiation blueward of the Ly{\alpha} line, the total optical depth to Thomson scattering of cosmic microwave background, with corresponding observations, we are able to constrain the epoch when the Universe was reionized to z <~ 8.5.

The structure of a sample of high-redshift (z=2), rotating galaxies with high star formation rates and turbulent gas velocities of sigma=40-80 km/s is investigated. Fitting the observed disk rotational velocities and radii with a Mo, Mao, White (1998) (MMW) model requires unusually large disk spin parameters lambda_d>0.1 and disk-to-dark halo mass fraction m_d=0.2, close to the cosmic baryon fraction. The galaxies segregate into dispersion-dominated systems with 1<vmax/sigma<3, maximum rotational velocities vmax<200 km/s and disk half-light radii rd=1-3 kpc and rotation-dominated systems with vmax>200 km/s, vmax/sigma>3 and rd=4-8 kpc. For the dispersion-dominated sample, radial pressure gradients partly compensate the gravitational force, reducing the rotational velocities. Including this pressure effect in the MMW model, dispersion-dominated galaxies can be fitted well with spin parameters lf lambda_d=0.03-0.05 for high disk mass fractions of m_d=0.2 and with lambda_d=0.01-0.03 for m_d=0.05. These values are in good agreement with cosmological expectations. For the rotation-dominated sample however pressure effects are small and better agreement with theoretically expected disk spin parameters can only be achieved if the dark halo mass contribution in the visible disk regime (2-3*rd) is smaller than predicted by the MMW model. We argue that these galaxies can still be embedded in standard cold dark matter halos if the halos did not contract adiabatically in response to disk formation. It is shown that the observed high turbulent gas motions of the galaxies are consistent with a Toomre instability parameter Q=1 which is equal to the critical value, expected for gravitational disk instability to be the major driver of turbulence. The dominant energy source of turbulence is then the potential energy of the gas in the disk.

The age and metallicity gradients for a sample of group and cluster galaxies from N-body+hydrodynamical simulation are analyzed in terms of galaxy stellar mass. Dwarf galaxies show null age gradient with a tail of high and positive values for systems in groups and cluster outskirts. Massive systems have generally zero age gradients which turn to positive for the most massive ones. Metallicity gradients are distributed around zero in dwarf galaxies and become more negative with mass; massive galaxies have steeper negative metallicity gradients, but the trend flatten with mass. In particular, fossil groups are characterized by a tighter distribution of both age and metallicity gradients. We find a good agreement with both local observations and independent simulations. The results are also discussed in terms of the central age and metallicity, as well as the total colour, specific star formation and velocity dispersion.

Numerical models of disc galaxy formation predict the existence of extended, hot ( T~10^6 K) gas haloes around present day spirals. The X-ray luminosity of these haloes is predicted to increase strongly with galaxy mass. However, searches for their X-ray emission have not been successful so far. We calculate the all sky O VII column density distributions for the haloes of three Milky Way like disc galaxies, resulting from cosmological high-resolution, N-body/gasdynamical simulations. We perform calculations both including the disc gas and without it, so the disc contribution to the column density is quantified. It is found that the column densities estimated for Milky Way-like galaxies are just below the observational upper limit, making a test of the hot halo paradigm likely within observational reach.

We compute the escape of ionizing radiation from galaxies in the redshift interval z=4-10, i.e., during and after the epoch of reionization, using a high-resolution set of galaxies, formed in fully cosmological simulations. The simulations invoke early, energetic feedback, and the galaxies evolve into a realistic population at z=0. Our galaxies cover nearly four orders of magnitude in masses (10^{7.8}-10^{11.5}\msun) and more than five orders in star formation rates (10^{-3.5}-10^{1.7}\msun\yr^{-1}), and we include an approximate treatment of dust absorption. We show that the source-averaged Lyman-limit escape fraction at z=10.4 is close to 80% declining monotonically with time as more massive objects build up at lower redshifts. Although the amount of dust absorption is uncertain to 1-1.5 dex, it is tightly correlated with metallicity; we find that dust is unlikely to significantly impact the observed UV output. These results support reionization by stellar radiation from low-luminosity dwarf galaxies and are also compatible with Lyman continuum observations and theoretical predictions at z\sim3-4.

In order to investigate the structure and dynamics of the recently discovered massive (M_* > 10^11 M_sun) compact z~2 galaxies, cosmological hydrodynamical/N-body simulations of a proto-cluster region have been undertaken. At z=2, the highest resolution simulation contains ~5800 resolved galaxies, of which 509, 27 and 5 have M_* > 10^10 M_sun, > 10^11 M_sun and > 4x10^11 M_sun, respectively. Effective radii and characteristic stellar densities have been determined for all galaxies. At z=2, for the definitely well resolved mass range of M_* > 10^11 Msun, the mass-size relation is consistent with observational findings for the most compact z~2 galaxies. The very high velocity dispersion recently measured for a compact z~2 galaxy (~510 km/s; van Dokkum et al 2009) can be matched at about the 1-sigma level, although a somewhat larger mass than the estimated M_* ~ 2 x 10^11 M_sun is indicated. For the above mass range, the galaxies have an average axial ratio <b/a> = 0.64 +/- 0.02 with a dispersion of 0.1, an average rotation to 1D velocity dispersion ratio <v/sigma> = 0.46 +/- 0.06 with a dispersion of 0.3, and a maximum value of v/sigma ~ 1.1. Rotation and velocity anisotropy both contribute in flattening the compact galaxies. Some of the observed compact galaxies appear flatter than any of the simulated galaxies. Finally, it is found that the massive compact galaxies are strongly baryon dominated in their inner parts, with typical dark matter mass fractions of order only 20% inside of r=2R_eff.

We present results obtained with a 3D, Ly alpha radiative transfer code, applied to a fully cosmological galaxy formation simulation. The developed Monte Carlo code is capable of treating an arbitrary distribution of source Ly alpha emission, neutral hydrogen density, temperature, and peculiar velocity of the interstellar medium. We investigate the influence of resonant scattering on the appearance and properties of young galaxies by applying the code to a simulated "Lyman Break Galaxy" at redshift z = 3.6, and of star formation rate 22 M_sun/yr and total Ly alpha luminosity 2.0 X 10^43 erg/s. It is found that resonant scattering of Ly alpha radiation can explain that young galaxies frequently are observed to be more extended on the sky in Ly alpha than in the optical. Moreover, it is shown that, for the system investigated, due to the anisotropic escape of the photons, the appearent maximum surface brightness can differ by a factor of ~15, and the total derived luminosity by a factor of ~4, depending on the orientation of the system relative to the observer.

A numerical code for solving various Lyman alpha (Lya) radiative transfer (RT) problems is presented. The code is suitable for an arbitrary, three-dimensional distribution of Lya emissivity, gas temperature, density, and velocity field. Capable of handling Lya RT in an adaptively refined grid-based structure, it enables detailed investigation of the effects of clumpiness of the interstellar (or intergalactic) medium. The code is tested against various geometrically and physically idealized configurations for which analytical solutions exist, and subsequently applied to three "Lyman-break galaxies", extracted from high-resolution cosmological simulations at redshift z = 3.6. Proper treatment of the Lya scattering reveals a diversity of surface brightness (SB) and line profiles. Specifically, for a given galaxy the maximum observed SB can vary by an order of magnitude, and the total flux by a factor of 3 - 6, depending on the viewing angle. This may provide an explanation for differences in observed properties of high-redshift galaxies, and in particular a possible physical link between Lyman-break galaxies and regular Lya emitters.

The Lyman alpha emission line is an essential diagnostic tool for probing galaxy formation and evolution. Not only is it commonly the strongest observable line from high-redshift galaxies but from its shape detailed information about its host galaxy can be revealed. However, due to the scattering nature of Lya photons increasing their path length in a non-trivial way, if dust is present in the galaxy the line may be severely suppressed and its shape altered. In order to interpret observations correctly, it is thus of crucial significance to know how much of the emitted light actually escapes the galaxy. In the present work, using a combination of high-resolution cosmological hydro-simulations and an adaptively refinable Monte Carlo Lya radiative transfer code including an advanced model of dust, the escape fractions f_esc of Lya radiation from high-redshift (z = 3.6) galaxies are calculated. In addition to the average escape fraction, the variation of f_esc in different directions and from different parts of the galaxies is investigated, as well as the effect on the emergent spectrum. Escape fractions from a sample of simulated galaxies of representative physical properties are found to decrease for increasing galaxy virial mass M_vir, from f_esc approaching unity for M_vir ~ 10^9 M_sun to f_esc less than 10% for M_vir ~ 10^12 M_sun. In spite of the dust being almost grey, it is found that the emergent spectrum is affected non-uniformly, with the escape fraction of photons close to the line center being much higher than of those in the wings, thus effectively narrowing the Lya line.

To assess the effect of baryonic ``pinching'' of galaxy cluster dark matter (DM) haloes, cosmological (LCDM) TreeSPH simulations of the formation and evolution of two galaxy clusters have been performed, with and without baryons included. The simulations with baryons invoke star formation, chemical evolution with non-instantaneous recycling, metallicity dependent radiative cooling, strong star-burst, driven galactic super-winds and the effects of a meta-galactic UV field, including simplified radiative transfer. The two clusters have T_X~3 and 6 keV, respectively, and, at z~0, both host a prominent, central cD galaxy. Comparing the simulations without and with baryons, it is found for the latter that the inner DM density profiles, r<50-100 kpc, steepen considerably: Delta(alpha)~0.5-0.6, where -alpha is the logarithmic DM density gradient. This is mainly due to the central stellar cDs becoming very massive, as a consequence of the onset of late time cooling flows and related star formation. Once these spurious cooling flows have been corrected for, and the cluster gravitational potentials dynamically adjusted, much smaller pinching effects are found: Delta(alpha)~0.1. Including the effects of baryonic pinching, central slopes of alpha~1.0 and 1.2 are found for the DM in the two clusters, interestingly close to recent observational findings. For the simulations with baryons, the inner density profile of DM+ICM gas combined is found to be only very marginally steeper than that of the DM, Delta(alpha)<0.05. However, the total matter inner density profiles are found to be Delta(alpha)~0.5 steeper than the inner profiles in the dark matter only simulations.

Models of disk galaxy formation commonly predict the existence of an extended reservoir of accreted hot gas surrounding massive spirals at low redshift. As a test of these models, we use X-ray and H-alpha data of the two massive, quiescent edge-on spirals NGC 5746 and NGC 5170 to investigate the amount and origin of any hot gas in their halos. Contrary to our earlier claim, the Chandra analysis of NGC 5746, employing more recent calibration data, does not reveal any significant evidence for diffuse X-ray emission outside the optical disk, with a 3-sigma upper limit to the halo X-ray luminosity of 4e39 erg/s. An identical study of the less massive NGC 5170 also fails to detect any extraplanar X-ray emission. By extracting hot halo properties of disk galaxies formed in cosmological hydrodynamical simulations, we compare these results to expectations for cosmological accretion of hot gas by spirals. For Milky Way-sized galaxies, these high-resolution simulations predict hot halo X-ray luminosities which are lower by a factor of ~2 compared to our earlier results reported by Toft et al. (2002). We find the new simulation predictions to be consistent with our observational constraints for both NGC 5746 and NGC 5170, while also confirming that the hot gas detected so far around more actively star-forming spirals is in general probably associated with stellar activity in the disk. Observational results on quiescent disk galaxies at the high-mass end are nevertheless providing powerful constraints on theoretical predictions, and hence on the assumed input physics in numerical studies of disk galaxy formation and evolution.

We analyze high resolution, N-body hydrodynamical simulations of fiducial galaxy clusters to probe tidal stripping of the dark matter subhalos. These simulations include a prescription for star formation allowing us to track the fate of the stellar component as well. We investigate the effect of tidal stripping on cluster galaxies hosted in these dark matter subhalos as a function of cluster-centric radius. To quantify the extent of the dark matter halos of cluster galaxies, we introduce the half mass radius r_half as a diagnostic, and study its evolution with projected cluster-centric distance R as a function of redshift. We find a well defined trend for (r_half,R): the closer the galaxies are to the center of the cluster, the smaller the half mass radius. Interestingly, this trend is inferred in all redshift frames examined in this work ranging from z=0 to z=0.7. At z=0, galaxy halos in the central regions of clusters are found to be highly truncated, with the most compact half mass radius of 10 kpc. We also find that r_half depends on luminosity and we present scaling relations of r_half with galaxy luminosity. The corresponding total mass of the cluster galaxies is also found to increase with projected cluster-centric distance and luminosity, but with more scatter than the (r_half,R) trend. Comparing the distribution of stellar mass to total mass for cluster galaxies, we find that the dark matter component is preferentially stripped, whereas the stellar component remains protected by the halo and is much less affected by tidal forces. We compare these results with galaxy-galaxy lensing probes of r_half and find qualitative agreement. (Abridged)

We present the modeling of SINFONI integral field dynamics of 18 star forming galaxies at z ~ 2 from Halpha line emission. The galaxies are selected from the larger sample of the SINS survey, based on the prominence of ordered rotational motions with respect to more complex merger induced dynamics. The quality of the data allows us to carefully select systems with kinematics dominated by rotation, and to model the gas dynamics across the whole galaxy using suitable exponential disk models. We obtain a good correlation between the dynamical mass and the stellar mass, finding that large gas fractions Mgas~M*) are required to explain the difference between the two quantities. We use the derived stellar mass and maximum rotational velocity Vmax from the modeling to construct for the first time the stellar mass Tully-Fisher relation at z ~ 2.2. The relation obtained shows a slope similar to what is observed at lower redshift, but we detect an evolution of the zero point. We find that at z ~ 2.2 there is an offset in log(M*) for a given rotational velocity of 0.41+-0.11 with respect to the local Universe. This result is consistent with the predictions of the latest N-body/hydrodynamical simulations of disk formation and evolution, which invoke gas accretion onto the forming disk in filaments and cooling flows. This scenario is in agreement with other dynamical evidence from SINS, where gas accretion from the halo is required to reproduce the observed properties of a large fraction of the z ~ 2 galaxies.

In this article, we apply strong lensing techniques in Abell 1703, a massive X-ray luminous galaxy cluster at z=0.28. Our analysis is based on imaging data both from space and ground in 8 bands, complemented with a spectroscopic survey. Abell 1703 looks rather circular from the general shape of its multiply imaged systems and present a dominant giant elliptical cD galaxy in its centre. This cluster exhibits a remarkable bright 'central ring' formed by 4 bright images at z_{spec}=0.888 located very close to the cD galaxy, providing observational constraints that are potentially very interesting to probe the central mass distribution. The stellar contribution from the cD galaxy (~1.25 10^{12} M_{sun} within 7") is accounted for in our parametric mass modelling, and the underlying smooth dark matter component distribution is described using a generalized NFW profile parametrized with a central logarithmic slope \alpha. We find that within the range where observational constraints are present (from ~5" to ~50"), the slope of the dark matter distribution in Abell1703 is equal to 1.09^{+0.05}_{-0.11} (3\sigma confidence level). The concentration parameter is equal to c_{200} ~ 3.5, and the scale radius is constrained to be larger than the region where observational constraints are available. Within this radius, the 2D mass is equal to M(50")=2.4 10^{14} M_{\sun}. We cannot draw any conclusions on cosmological models at this point since we lack results from realistic numerical simulations containing baryons to make a proper comparison. We advocate the need for a sample of observed and simulated unimodal relaxed galaxy clusters in order to make reliable comparisons, and potentially provide a test of cosmological models.

N-body + hydrodynamical simulations of the formation and evolution of galaxy groups and clusters in a LambdaCDM cosmology are used in order to follow the building-up of the colour-magnitude relation in two clusters and in 12 groups. We have found that galaxies, starting from the more massive, move to the Red Sequence (RS) as they get aged over times and eventually set upon a ``dead sequence'' (DS) once they have stopped their bulk star formation activity. Fainter galaxies keep having significant star formation out to very recent epochs and lie broader around the RS. Environment plays a role as galaxies in groups and cluster outskirts hold star formation activity longer than the central cluster regions. However galaxies experiencing infall from the outskirts to the central parts keep star formation on until they settle on to the DS of the core galaxies. Merging contributes to mass assembly until z~1, after which major events only involve the brightest cluster galaxies. The emerging scenario is that the evolution of the colour-magnitude properties of galaxies within the hierarchical framework is mainly driven by star formation activity during dark matter halos assembly. Galaxies progressively quenching their star formation settle to a very sharp ``red and dead'' sequence, which turns out to be universal, its slope and scatter being almost independent of the redshift (since at least z~1.5) and environment. Differently from the DS, the operatively defined RS evolves more evidently with z, the epoch when it changes its slope being closely corresponding to that at which the passive galaxies population takes over the star forming one: this goes from z~1 in clusters down to 0.4 in normal groups.

We test the hypothesis that the host galaxies of long-duration gamma-ray bursts (GRBs) as well as quasar-selected damped Lyman-alpha (DLA) systems are drawn from the population of UV-selected star-forming, high-z galaxies (generally referred to as Lyman-break galaxies). Specifically, we compare the metallicity distributions of the GRB and DLA populations to simple models where these galaxies are drawn randomly from the distribution of star-forming galaxies according to their star-formation rate and HI cross-section respectively. We find that it is possible to match both observational distributions assuming very simple and constrained relations between luminosity, metallicity and HI sizes. The simple model can be tested by observing the luminosity distribution of GRB host galaxies and by measuring the luminosity and impact parameters of DLA selected galaxies as a function of metallicity. Our results support the expectation that GRB and DLA samples, in contrast to magnitude limited surveys, provide an almost complete census of z=3 star-forming galaxies that are not heavily dust-obscured.

The median observed velocity width v_90 of low-ionization species in damped Ly-alpha systems is close to 90 km/s, with approximately 10% of all systems showing v_90 > 210 km/s at z=3. We show that a relative shortage of such high-velocity neutral gas absorbers in state-of-the-art galaxy formation models is a fundamental problem, present both in grid-based and particle-based numerical simulations. Using a series of numerical simulations of varying resolution and box size to cover a wide range of halo masses, we demonstrate that energy from gravitational infall alone is insufficient to produce the velocity dispersion observed in damped Ly-alpha systems, nor does this dispersion arise from an implementation of star formation and feedback in our highest resolution (~ 45 pc) models, if we do not put any galactic winds into our models by hand. We argue that these numerical experiments highlight the need to separate dynamics of different components of the multiphase interstellar medium at z=3.

[Abridged] The detection of the rotational lines of CO in proto-galaxies in the early Universe provides one of the most promising ways of probing the fundamental physical properties of a galaxy, such as its size, dynamical mass, gas density, and temperature. While such observations are currently limited to the most luminous galaxies, the advent of ALMA will change the situation dramatically, resulting in the detection of numerous normal galaxies at high redshifts. Maps and spectra of rotational CO line emission were calculated from a cosmological N-body/hydrodynamical TreeSPH simulation of a z ~ 3 "Lyman break galaxy'' of UV luminosity about one order of magnitude below L*. To simulate a typical observation of our system with ALMA, we imposed characteristic noise, angular, and spectral resolution constraints. The CO line properties predicted by our simulation are in good agreement with the two Lyman break systems detected in CO to date. We find that while supernovae explosions from the ongoing star formation carve out large cavities in the molecular ISM, they do not generate large enough gas outflows to make a substantial imprint on the CO line profile. This implies that for most proto-galaxies - except possibly the most extreme cases - stellar feedback effects do not affect CO as a dynamical mass tracer. Detecting CO in sub-L* galaxies at z ~3 will push ALMA to the limits of its cababilities, and whether a source is detected or not may depend critically on its inclination angle. Both these effects (sensitivity and inclination) will severely impair the ability of ALMA to infer the gas kinematics and dynamical masses using line observations.

The Sunyaev-Zeldovich (SZ) effect gives a measure of the thermal energy and electron pressure in groups and clusters of galaxies. In the near future SZ surveys will map hundreds of systems, shedding light on the pressure distribution in the systems. The thermal energy is related to the total mass of a system of galaxies, but it is only a projection that is observed through the SZ effect. A model for the 3D distribution of pressure is needed to link the SZ signal to the total mass of the system. In this work we construct an empirical model for the 2D and 3D SZ profile, and compare it to a set of realistic high resolution SPH simulations of galaxy clusters and groups, and to a stacked SZ profile for massive clusters derived from WMAP data. Furthermore, we combine observed temperature profiles with dark matter potentials to yield an additional constraint, under the assumption of hydrostatic equilibrium. We find a very tight correlation between the characteristic scale in the model, the integrated SZ signal, and the total mass in the systems with a scatter of only 4%. The model only contains two free parameters, making it readily applicable even to low resolution SZ observations of galaxy clusters. A fitting routine for the model that can be applied to observed or simulated data can be found at http://www.phys.au.dk/~haugboel/software.shtml

The temperature distribution of the cosmic gas-phase oxygen at z~3 is determined by combining high resolution cosmological simulations with the observed, (rest-frame) V-band galaxy LF. The simulations invoke three different stellar IMFs, a Kroupa (K98), a Salpeter (S) and an Arimoto-Yoshii (AY), spanning a range of a factor of five in chemical yield and specific SNII energy feedback. Oxygen is found in all T phases (cold: logT<~4, warm: logT~4.5, hot: logT>~5), in particular for the (top-heavy) AY IMF. For the K98 and S IMFs, the most important phase is the cold one. Moreover, the cold phase alone contains 1-3 times the mass of oxygen in galactic stars for the three IMFs. The implications of this in relation to DLA studies are discussed. In relation to ``missing metals'' it is found that a significant fraction of the oxygen is located in a warm/hot phase that may be very difficult to detect. Moreover, it is found that less than about 20-25% of the cosmic oxygen is associated with galaxies brighter than M_V~-22, the faintest LBGs probed by current metallicity determinations. From the LBG based, lambda~1500 AA luminosity density history at z>3, the mean oxygen density at z=3 is inferred. Comparing this to the models, it is found that the (solar neighbourhood type) K98 IMF is strongly excluded (too small yield), the S is marginally excluded, and the AY matches the constraint well. The optimal IMF has a yield between the S and AY. A similar conclusion is reached on the basis of the C IV abundance of the moderate density IGM.

Ongoing accretion onto galactic disks has been recently theorized to progress via the unstable cooling of the baryonic halo into condensed clouds. These clouds have been identified as analogous to the High-Velocity Clouds (HVCs) observed in HI in our Galaxy. Here we compare the distribution of HVCs observed around our own Galaxy and extra-planar gas around the Andromeda galaxy to these possible HVC analogs in a simulation of galaxy formation that naturally generates these condensed clouds. We find a very good correspondence between these observations and the simulation, in terms of number, angular size, velocity distribution, overall flux and flux distribution of the clouds. We show that condensed cloud accretion only accounts for ~ 0.2 M_solar / year of the current overall Galactic accretion in the simulations. We also find that the simulated halo clouds accelerate and become more massive as they fall toward the disk. The parameter space of the simulated clouds is consistent with all of the observed HVC complexes that have distance constraints, except the Magellanic Stream which is known to have a different origin. We also find that nearly half of these simulated halo clouds would be indistinguishable from lower-velocity gas and that this effect is strongest further from the disk of the galaxy, thus indicating a possible missing population of HVCs. These results indicate that the majority of HVCs are consistent with being infalling, condensed clouds that are a remnant of Galaxy formation.

Based on cosmological simulations, we model Lyman continuum emission from a sample of 11 high-redshift star forming galaxies spanning a mass range of a factor 20. Each of the 11 galaxies has been simulated both with a Salpeter and a Kroupa initial mass function (IMF). We find that the Lyman continuum (LyC) luminosity of an average star forming galaxy in our sample declines from z=3.6 to 2.4 due to the steady gas infall and higher gas clumping at lower redshifts, increasingly hampering the escape of ionizing radiation. The galaxy-to-galaxy variation of apparent LyC emission at a fixed redshift is caused in approximately equal parts by the intrinsic variations in the LyC emission and by orientation effects. The combined scatter of an order of magnitude can explain the variance in the far-UV spectra of high-redshift galaxies detected by Shapley et al. (2006). Our results imply that the cosmic galactic ionizing UV luminosity would be monotonically decreasing from z=3.6 to 2.4, curiously anti-correlated with the star formation rate in the smaller galaxies, which on average rises during this redshift interval.

We report coincident damped Lyman alpha (DLA) and sub-DLA absorption at z = 2.66 and z = 2.94 towards the z ~ 3 13.8 arcsecond separation binary quasar SDSS 1116+4118 AB. At the redshifts of the absorbers, this angular separation corresponds to a proper transverse separation of ~ 110 kpc. A third absorber, a sub-DLA at z = 2.47, is detected towards SDSS 1116+4118 B, but no corresponding high column density absorber is present towards SDSS 1116+4118 A. We use high resolution galaxy simulations and a clustering analysis to interpret the coincident absorption and its implications for galaxy structure at z ~ 3. We conclude that the common absorption in the two lines of sight is unlikely to arise from a single galaxy, or a galaxy plus satellite system, and is more feasibly explained by a group of two or more galaxies with separations ~ 100 kpc. The impact of these findings on single line of sight observations is also discussed; we show that abundances of DLAs may be affected by up to a few tenths of a dex by line of sight DLA blending. From a Keck ESI spectrum of the two quasars, we measure metal column densities for all five absorbers and determine abundances for the three absorbers with log N(HI) > 20. For the two highest N(HI) absorbers, we determine high levels of metal enrichment, corresponding to 1/3 and 1/5 solar. These metallicities are amongst the highest measured for DLAs at any redshift and are consistent with values measured in Lyman break galaxies at 2 < z < 3. For the DLA at z = 2.94 we also infer an approximately solar ratio of alpha-to-Fe peak elements from [S/Zn] = +0.05, and measure an upper limit for the molecular fraction in this particular line of sight of log f(H_2)< -5.5.

High resolution cosmological N-body simulations of four galaxy-scale dark matter halos are compared to corresponding N-body/hydrodynamical simulations containing dark matter, stars and gas. The simulations without baryons share features with others described in the literature in that the dark matter density slope continuously decreases towards the center, with a density \rho~r^{-1.3+/-0.2}, at about 1% of the virial radius for our Milky Way sized galaxies. The central cusps in the simulations which also contain baryons steepen significantly, to \rho~r^{-1.9+/-0.2}, with an indication of the inner logarithmic slope converging. Models of adiabatic contraction of dark matter halos due to the central build-up of stellar/gaseous galaxies are examined. The simplest and most commonly used model, by Blumenthal et al., is shown to overestimate the central dark matter density considerably. A modified model proposed by Gnedin et al. is shown to be a considerable improvement, but not perfect. Moreover it is found that the contraction parameters not only depend on the orbital structure of the dark-matter-only halos but also on the stellar feedback prescription which is most relevant for the baryonic distribution. Implications for dark matter annihilation at the galactic center are discussed and it is found that although our simulations show a considerable reduced halo contraction as compared to the Blumenthal et al. model, the fluxes from dark matter annihilation is still expected to be enhanced by at least a factor of a hundred as compared to dark-matter-only halos. Finally, it is shown that while dark-matter-only halos are typically prolate, the dark matter halos containing baryons are mildly oblate with minor-to-major axis ratios of c/a=0.73+/-0.11, with their flattening aligned with the central baryonic disks.

N--body + hydrodynamical simulations of formation and evolution of galaxy groups in a LambdaCDM cosmology have been performed. The properties of the galaxy populations in 12 groups are here discussed, with focus on the colour-magnitude relation in both normal and fossil groups.

We present predictions on the evolution of the Tully-Fisher (TF) relation with redshift, based on cosmological N-body/hydrodynamical simulations of disc galaxy formation and evolution. The simulations invoke star formation and stellar feedback, chemical evolution with non-instantaneous recycling, metallicity dependent radiative cooling and effects of a meta-galactic UV field, including simplified radiative transfer. At z=0, the simulated and empirical TF relations are offset by about 0.4 magnitudes (1 sigma) in the B and I bands. The origin of these offsets is somewhat unclear, but it may not necessarily be a problem of the simulations only. As to evolution, we find a brightening of the TF relation between z=0 and z=1 of about 0.85 mag in rest-frame B band, with a non-evolving slope. The brightening we predict is intermediate between the (still quite discrepant) observational estimates. This evolution is primarily a luminosity effect, while the stellar mass TF relation shows negligible evolution. The individual galaxies do gain stellar mass between z=1 and z=0, by a 50-100%; but they also correspondingly increase their characteristic circular speed. As a consequence, individually they mainly evolve ALONG the stellar mass TF relation, while the relation as such does not show any significant evolution.

Models of disk galaxy formation commonly predict the existence of an extended reservoir of hot gas surrounding massive spirals at low redshift. As a test of these models, we have obtained X-ray and optical data of the two massive edge-on spirals NGC 5746 and NGC 5170, in order to investigate the amount and origin of hot gas in their disks and halos. Chandra observations of NGC 5746 reveal evidence for diffuse X-ray emission with a total luminosity of ~7 x 10^39 erg/s surrounding this galaxy out to at least ~20 kpc from the disk, whereas an identical study of the less massive NGC 5170 fails to detect any extraplanar X-ray emission. Unlike the case for other disk galaxies with detected X-ray halos, the halo emission around NGC 5746 is not accompanied by extraplanar H-alpha or radio emission, and there is no evidence for significant nuclear or starburst activity in the disk. In contrast to these other cases, the emission around NGC 5746 therefore appears to arise from the cooling of externally accreted material rather than from disk outflows. To verify this idea, we present results of cosmological simulations of galaxy formation and evolution, showing our observations to be in good agreement with expectations for cosmological accretion, while also confirming that the X-ray halos of other spirals do not fit well into an accretion scenario. We find that the estimated cooling rate of hot halo gas around NGC 5746 would provide sufficient material for star formation in the disk to proceed at its present rate. This lends support to the idea that a supply of hot ambient gas is potentially available as fuel for star formation in massive, nearby spirals, and suggests that accretion of hot gas could be important for maintaining the stellar disks of such galaxies.

We report broadband UBV and/or BVRI CCD photometry for a total of 1857 stars in the thick-disk and halo populations of the Galaxy. The majority of our targets were selected as candidate field horizontal-branch or other A-type stars (FHB/A, N = 576), or candidate low-metallicity stars (N = 1221), from the HK and Hamburg/ESO objective-prism surveys. Similar data for a small number of additional stars from other samples are also reported. These data are being used for several purposes. In the case of the FHB/A candidates they are used to accurately separate the lower-gravity FHB stars from various higher-gravity A-type stars, a subsample that includes the so-called Blue Metal Poor stars, halo and thick-disk blue stragglers, main-sequence A-type dwarfs, and Am and Ap stars. These data are also being used to derive photometric distance estimates to high-velocity hydrogen clouds in the Galaxy and for improved measurements of the mass of the Galaxy. Photometric data for the metal-poor candidates are being used to refine estimates of stellar metallicity for objects with available medium-resolution spectroscopy, to obtain distance estimates for kinematic analyses, and to establish initial estimates of effective temperature for analysis of high-resolution spectroscopy of the stars for which this information now exists.

Using results from high-resolution galaxy formation simulations in a standard Lambda-CDM cosmology and a fully conservative multi-resolution radiative transfer code around point sources, we compute the energy-dependent escape fraction of ionizing photons from a large number of star forming regions in two galaxies at five different redshifts from z=3.8 to 2.39. All escape fractions show a monotonic decline with time, from (at the Lyman-limit) ~6-10% at z=3.6 to ~1-2% at z=2.39, due to higher gas clumping at lower redshifts. It appears that increased feedback can lead to higher f_esc at z>3.4 via evacuation of gas from the vicinity of star forming regions and to lower f_esc at z<2.39 through accumulation of swept-up shells in denser environments. Our results agree well with the observational findings of \citet{inoue..06} on redshift evolution of f_esc in the redshift interval z=2-3.6.

We have performed a series of N-body/hydrodynamical (TreeSPH) simulations of clusters and groups of galaxies, selected from cosmological N-body simulations within a $\Lambda$CDM framework: these objects have been re-simulated at higher resolution to $z$=0, in order to follow also the dynamical, thermal and chemical input on to the ICM from stellar populations within galaxies. The simulations include metal dependent radiative cooling, star formation according to different IMFs, energy feedback as strong starburst-driven galactic super-winds, chemical evolution with non-instantaneous recycling of gas and heavy elements, effects of a meta-galactic UV field and thermal conduction in the ICM. In this Paper I of a series of three, we derive results, mainly at $z=0$, on the temperature and entropy profiles of the ICM, its X-ray luminosity, the cluster cold components (cold fraction as well as mass--to--light ratio) and the metal distribution between ICM and stars. In general, models with efficient super-winds, along with a top-heavy stellar IMF, are able to reproduce fairly well the observed $L_X-T$ relation, the entropy profiles and the cold fraction. Observed radial ICM temperature profiles can be matched, except for the gradual decline in temperature inside of $r\sim$~0.1$R_{\rm{vir}}$. Metal enrichment of the ICM gives rise to somewhat steep inner iron gradients; yet, the global level of enrichment compares well to observational estimates after correcting for the stars formed at late times at the base of the cooling flows; also the metal partition between stars and ICM gets into good agreement with observations.

Extended gas haloes around galaxies are a ubiquitous prediction of galaxy formation scenarios. However, the density profiles of this hot halo gas is virtually unknown, although various profiles have been suggested on theoretical grounds. In order to quantitatively address the gas profile, we compare galaxies from direct cosmological simulations with analytical solutions of the underlying gas equations. We find remarkable agreement between simulations and theoretical predictions. We present an expression for this gas profile with a non-trivial dependence on the total mass profile. This expression is useful when setting up equilibrium galaxy models for numerical experiments.

We report on the discovery of a z = 3.16 Lyman-alpha emitting blob in the GOODS South field. The blob has a total Ly-alpha luminosity of ~ 10^(43) erg s^(-1) and a diameter larger than 60 kpc. The available multi-wavelength data in the GOODS field consists of 13 bands from X-rays (Chandra) to infrared (Spitzer). Unlike other discovered Ly-alpha blobs, this blob shows no obvious continuum counter-part in any of the broad-bands. In particular, no optical counter-parts are found in the deep HST/ACS imaging available. For previously published blobs, AGN (Active Galactic Nuclei) or 'superwind' models have been found to provide the best match with the data. We here argue that the most probable origin of the extended Ly-alpha emission from the blob in the GOODS South field is cold accretion onto a dark matter halo.

Cosmological LCDM simulations of 12 M_vir~10^14 Msun galaxy groups have been performed, invoking star formation, chemical evolution with non-instantaneous recycling, metallicity dependent radiative cooling, strong star-burst driven galactic super-winds and effects of a meta-galactic UV field. At z=0, intra-group light (IGL) fractions are found to be 12-45%. Low values refer to groups with only a small difference between the R-band magnitudes of the first and second ranked group galaxy, large are typical of "fossil" groups (FGs). The IG stars in the 4 FGs are 0.3-0.5 Gyr older than in the 8 nonFGs. For the IGL, B-R=~1.4, in good agreement with observations. For FGs/nonFGs the iron abundance of the IG stars is slightly sub-solar in the central parts (r~100 kpc) decreasing to about 40% solar at about 0.5 r_vir The IG stars are alpha-element enhanced with [O/Fe] increasing with r, and an overall [O/Fe]~0.45, indicating predominant SNII enrichment. The velocity distributions of the IG stars and group galaxies are, at r>~30 kpc, significantly more radially anisotropic for FGs than for nonFGs. So a characteristic of FG formation, apart from formation time (D'Onghia et al.), may be the "initial" velocity distribution of the group galaxies. For FGs one can dynamically infer the (dark matter dominated) mass distribution of the groups all the way to r_vir, from the kinematics of the IG stars or group galaxies. For the nonFGs this method overestimates the group mass at r>~200 kpc, by up to a factor of two at r_vir. This is interpreted as FGs being, in general, more relaxed than nonFGs. Finally, FGs of the above M_vir should host ~500 planetary nebulae at projected distances between 100 and 1000 kpc from the first ranked galaxy. All results appear consistent with the FG formation scenario of D'Onghia et al.

Cosmological (LCDM) TreeSPH simulations of galaxy clusters have been performed, invoking star-formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star-burst and/or AGN driven super-winds, effects of a meta-galactic UV field and thermal conduction. Results for two clusters, one Virgo-like (T~3 keV) and one (sub) Coma-like (T~6 keV), are presented. At z=0, the stellar contents of both clusters consist of a central dominant (cD) galaxy surrounded by galaxies and intracluster (IC) stars. The IC stars are found to contribute 20-40% of the total cluster B-band light and to form at a mean redshift z_f~3, on average about 0.5 Gyr prior to the stars in galaxies. UBVRIJHK surface brightness profiles of the IC star populations are presented; the profile of the larger cluster matches the observed V-band profile of the cD in Abell 1413 (T~8 keV). The IC stellar populations have B-R=1.4-1.5, comparable to larger E galaxies. The mean [Fe/H]~0.0 in the central part of the cluster (r~100 kpc) decreasing to about -0.3 at r_vir. The IC stars are alpha-element enhanced, with a weak trend of [O/Fe] increasing with r, and an overall [O/Fe]~0.4 dex. The IC stars are kinematically colder than the galaxies: The velocity dispersions of the IC stars are at r~100-500 kpc only about half of those of the galaxies increasing to about 70% at r=1-2 Mpc. The projected velocity dispersion in the Virgo-like cluster at R>~50 kpc is 300-600 km/s depending on orientation and R. Rotation is dynamically insignificant for the IC stars. The velocity distributions of IC stars and galaxies are in one cluster highly radially anisotropic, in the other close to isotropic.

We have performed TreeSPH simulations of disk galaxy formation in various warm dark matter (WDM) cosmologies. Our results indicate that for a range of WDM free-streaming masses, the disk galaxy formation angular momentum problem can be completely resolved by going to the WDM structure formation scenario, without having to invoke stellar feedback processes at all. We also confirm our previous suspicion, that part of the angular momentum problem is due to numerical effects, most likely related to the shock capturing, artificial viscosity used in SPH. Furthermore we find that we can match the observed I-band Tully-Fisher (TF) relation, provided that the I-band mass-to-light ratio of disk galaxies is about 0.8. We argue that this is quite a reasonable value in comparison with various dynamical and spectrophotometric estimates, including one given in this paper. We speculate that our success in matching the TF relation may be due to WDM halos being less centrally concentrated than CDM halos and suggest to check this exciting possibility with high resolution simulations, in particular in low Omega_M, WDM cosmologies. Finally, we discuss possible physical candidates for WDM particles extensively. We find that the most promising are neutrinos with weaker or stronger interactions than normal, majorons (light pseudogoldstone bosons) or mirror or shadow world neutrinos.

Based on 19 high-resolution N-body/gas-dynamical galaxy formation simulations in the LCDM cosmology it is shown, that for a galaxy like the Milky Way, in addition to the baryonic mass of the galaxy itself, about 70% extra baryonic mass should reside around the galaxy (inside of the virial radius), chiefly in the form of hot gas. Averaging over the entire field galaxy population, this ``external'' component amounts to 64-85% of the baryonic mass of the population itself. These results are supported by the recent detection of very extended, soft X-ray emission from the halo of the quiescent, massive disk galaxy NGC 5746. Some of the hot gas may, by thermal instability, have condensed into mainly pressure supported, warm clouds, similar to the Galactic High Velocity Clouds (HVCs). Based on an ultra-high resolution cosmological test simulation of a Milky Way like galaxy (with a gas particle mass and gravity softening length of only 7600 h^-1 Msun and 83 h^-1 pc,respectively), it is argued, that the hot gas phase dominates over the warm gas phase, in the halo. Finally, an origin of HVCs as ``leftovers'' from filamentary, ``cold'' accretion events, mainly occurring early in the history of galaxies, is proposed.

Hot gaseous haloes surrounding galaxies and extending well beyond the distribution of stars are a ubiquitous prediction of galaxy formation scenarios. The haloes are believed to consist of gravitationally trapped gas with a temperature of millions of Kelvin. The existence of such hot haloes around massive elliptical galaxies has been established through their X-ray emission. While gas out-flowing from starburst spiral galaxies has been detected, searches for hot haloes around normal, quiescent spiral galaxies have so far failed, casting doubts on the fundamental physics in galaxy formation models. Here we present the first detection of a hot, large-scale gaseous halo surrounding a normal, quiescent spiral galaxy, NGC 5746, alleviating a long-standing problem for galaxy formation models. In contrast to starburst galaxies, where the X-ray halo can be powered by the supernova energy, there is no such power source in NGC 5746. The only compelling explanation is that we are here witnessing a galaxy forming from gradually in-flowing hot and dilute halo gas.

Extended, soft X-ray emission from the halo of a very large disk galaxy has been detected out to r~20-40 kpc. The luminosity and surface brightness distribution is in very good agreement with predictions by recent, cosmological galaxy formation models. Predictions of Ly-alpha emission, associated with ``cold'' accretion of filamentary gas onto galaxies, are presented. In particular these predictions are compared to the properties of a recently detected Ly-alpha ``blob'', which shows no obvious continuum counter-part in any of 12 broad-bands, from X-rays (Chandra) to infrared (Spitzer). Finally, the predicted evolution of the Tully-Fisher relation, going from z=0 to 1, is discussed in relation to recent observations.

We use a set of twelve high-resolution N-body/hydrodynamical simulations in the $\Lambda$CDM cosmology to investigate the origin and formation rate of fossil groups (FGs), which are X-ray bright galaxy groups dominated by a large elliptical galaxy, with the second brightest galaxy being at least two magnitudes fainter. The simulations invoke star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst driven galactic super winds, effects of a meta-galactic UV field and full stellar population synthesis. We find an interesting correlation between the magnitude gap between the first and second brightest galaxy and the formation time of the group. It is found that FGs have assembled half of their final dark matter mass already at $z\ga1$, and subsequently typically grow by minor merging only, wheras non-FGs on average form later. The early assembly of FGs leaves sufficient time for galaxies of $L \sim L_*$ to merge into the central one by dynamical friction, resulting in the large magnitude gap at $z=0$. A fraction of 33$\pm$16% of the groups simulated are found to be fossil, whereas the observational estimate is $\sim$10-20%. The FGs are found to be X-ray over-luminous relative to non-FGs of the same optical luminosity, in qualitative agreement with observations. Finally, from a dynamical friction analysis is found that only because infall of $L \sim L_*$ galaxies happens along filaments with small impact parameters do FGs exist at all.

Cosmological (LambdaCDM) TreeSPH simulations of the formation and evolution of galaxy groups and clusters have been performed. The simulations invoke star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst and (optionally) AGN driven galactic super winds, effects of a meta-galactic UV field and thermal conduction. The properties of the galaxy populations in two clusters, one Virgo-like (T~3 keV) and one (sub) Coma-like (T~6 keV), are discussed. The global star formation rates of the cluster galaxies are found to decrease very significantly with time from redshift z=2 to 0, in agreement with observations. The total K-band luminosity of the cluster galaxies correlates tightly with total cluster mass, and for models without additional AGN feedback, the zero point of the relation matches the observed one fairly well. The match to observed galaxy luminosity functions is reasonable, except for a deficiency of bright galaxies (M_B < -20), which becomes increasingly significant with super-wind strength. Results of a high resolution test indicate that this deficiency is not due to ``over--merging''. The redshift evolution of the luminosity functions from z=1 to 0 is mainly driven by luminosity evolution, but also by merging of bright galaxies with the cD. The colour--magnitude relation of the cluster galaxies matches the observed "red sequence" very well and, on average, galaxy metallicity increases with luminosity. As the brighter galaxies are essentially coeval, the colour--magnitude relation results from metallicity rather than age effects, as observed.

Cosmological LambdaCDM TreeSPH simulations of the formation and evolution of galaxy groups and clusters have been performed. The simulations include: star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst and (optionally) AGN driven galactic super winds, effects of a meta-galactic UV field and thermal conduction. We report results on the temperature and entropy profiles of the ICM, the X-ray luminosity, cold fraction, M/L and IMLR ratios in gas and stars, metal abundances and gradients. Besides, the properties of the galaxy populations in the two richest clusters are discussed: global star formation rates of the cluster galaxies, the total K-band luminosity, the galaxy luminosity functions at z=0 and their redshift evolution, the colour-magnitude relation (``red sequence'') as resulting from metallicity effects, and the role of the IMF in reproducing colours and abundances of the stellar populations. Finally the contribute (20-40 %) to cluster light from the intra-cluster stars and the cD galaxy has been investigated: surface brightness profiles, mean colours and chemical abundances, kinematics (velocity dispersion and distributions) of the IC stellar populations.

The X-ray luminosities of the hot halo gas around simulated, Milky Way like disk galaxies have been determined, as a function of redshift. The X-ray luminosity increases significantly with redshift, in some cases as much as a factor 30 going from z=0 to 2. Consequently, the optimal detection redshift can be >~0.5. Results of fully cosmological simulations of galaxy groups and clusters, incorporating star formation, chemical evolution with non-instantaneous recycling, metal-dependent radiative cooling, galactic super-winds and thermal conduction are presented. X-ray luminosities at z=0 are somewhat high and central entropies somewhat low compared to observations. This is likely a combined effect of chemical evolution and metal-dependent, radiative cooling. Central ICM abundance profiles are somewhat steep, and the observed level of ICM enrichment can only be reproduced with IMFs more top-heavy than the Salpeter IMF. In agreement with observations it is found that the iron in the ICM is in place already at z~1. The [Si/Fe] of the ICM decreases with time, and increases slightly with radius at a given time. The cluster galaxies match the observed "red sequence" very well, and the metallicity of cluster galaxies increases with galaxy mass, as observed.

Cosmological simulations of disc galaxy formation, when compared to the observed Tully-Fisher relation, suggest a low Mass-to-Light (M/L) ratio for the stellar component in spirals. We show that a number of "bottom-light" Initial Mass Functions (IMFs) suggested independently in literature, do imply M/L ratios as low as required, at least for late type spirals (Sbc/Sc). However the typical M/L ratio, and correspondingly the zero-point of the Tully-Fisher relation, is expected to vary considerably with Hubble type. Bottom-light IMFs tend to have a metal production in excess of what is tipically estimated for spiral galaxies. Suitable tuning of the IMF slope and mass limits, post-supernova fallback of metals onto black holes or metal outflows must then be invoked, to reproduce the observed chemical properties of disc galaxies.

X-ray luminosities and surface brightness profiles of the hot gas haloes of simulated disc galaxies at redshifts z=0-2 are presented. The galaxies are extracted from fully cosmological simulations and correspond in mass to the Milky Way. We find that the bolometric X-ray luminosities of the haloes decrease by a factor 4-10 from z~1 to z~0, reflecting the decrease in the rate at which hot halo gas cools out on to the disc. At all redshifts, most of the emission is found to originate within 10-15 kpc of the disc. When combined with models in which the evolution of disc X-ray luminosity is dominated by X-ray binaries, the predicted halo luminosities at z~1 show good agreement with constraints from spiral galaxies in Chandra Deep Field data. There is an indication that haloes with a metal abundance of 0.3 solar overpredict observed X-ray luminosities at z~1, suggesting that halo metallicities are lower than this value. Prospects for direct detection of the haloes of Milky Way-sized galaxies with current and future X-ray instrumentation are discussed. It is found that XEUS should be able to single out the halo emission of highly inclined Milky Way-sized disc galaxies out to approximately z=0.3. For such galaxies in this redshift interval, we estimate a lower limit to the surface density of detectable haloes on the sky of ~10 deg^(-2). More generally, owing to their luminosity evolution, the optimum redshifts at which to observe such haloes could be 0.5<z<1, depending on their assembly history.

It is frequently debated in literature whether a "standard" Initial Mass Function (IMF) - meaning an IMF of the kind usually adopted to explain the chemical evolution in the local Solar Neighbourhood - can account for the observed metal enrichment and Iron Mass-to-Light Ratio in clusters of galaxies. We address this problem by means of straightforward estimates that should hold independently of the details of chemical evolution models. It is crucial to compute self-consistently the amount of mass and metals locked-up in stars, by accounting for the stellar mass-to-light ratio predicted by a given IMF. It becomes then clear that a "standard" Solar Neighbourhood IMF cannot provide enough metals to account for the observed chemical properties in clusters: clusters of galaxies and the local environment must be characterized by different IMFs. Alternatively, if we require the IMF to be universal, in order to explain clusters such an IMF must be much more efficient in metal production than usually estimated for the Solar Vicinity. In this case, substantial loss of metals is required from the Solar Neighbourhood and from disc galaxies in general. This "non-standard" scenario for the local chemical evolution would challenge our present understanding of the Milky Way and of disc galaxy formation.

TreeSPH simulations of galaxy formation in a LCDM cosmology, with star formation, energetic stellar feedback and a meta-galactic UV field have been performed, resulting in realistic disk, S0 and E galaxies at z=0. The disk galaxies are deficient in angular momentum by only about a factor of two compared to observations for runs where fairly strong star-bursts in early, proto-galactic clouds lead to "blow-away" of the remaining gas. The surface density profiles of the stellar disks are approximately exponential and those of the bulges range from exponential to r^(1/4). B/D ratios and integrated B-V colours are consistent with observations. The observed I-band TF relation can be matched with M/L_I ~ 0.8, in fair agreement with recent determinations. The (E/S0)s have approximately r^(1/4) profiles, non disk-like kinematics and are flattened due to non-isotropic stellar velocity distributions. We predict hot halo gas to cool out and accrete onto the Galactic disk at 0.5-1 Msun/yr at z=0, consistent with upper limits from FUSE observations of O VI. We analyzed two Milky Way-sized galaxies and find accretion rates, and X-ray halo luminosities, 6-7 times larger at z=1 than at z=0. The gas infall declines nearly exponentially with time, supporting the approximation often used in chemical evolution models. The infall time-scales are comparable to what is used to solve the "G-dwarf problem". One disk forms "inside-out", the other "outside-in". For both, the mean stellar ages in the outskirts agree with findings for the disk of M31. The amount of hot gas in disk galaxy haloes is consistent with observational upper limits, as are dispersion measures to pulsars in the globular cluster M53 and the LMC, which were "inserted" in the disk galaxy haloes.

A low mass-to-light ratio for the stellar component of spiral galaxies (M/L<1 in the I-band) is advocated by various dynamical arguments and by recent cosmological simulations of the formation of these systems. We discuss this possibility by means of chemo-photometric models for galactic discs, adopting different Initial Mass Functions. We show that a number of "bottom--light" IMFs (namely, with less mass locked in low-mass stars than the standard Salpeter IMF), suggested independently in recent literature, do imply M/L ratios as low as mentioned above, at least for late type spirals (Sbc/Sc). This conclusion still holds when the bulge contribution to mass and light is included. We also predict the typical stellar M/L ratio, and correspondingly the zero-point of the Tully-Fisher relation, to vary considerably with Hubble type (about 0.5-0.7 mag in the red bands, from Sa to Sc type). For some of the bottom-light IMFs considered, the efficiency of metal production tends to exceed what is typically estimated for spiral galaxies. Suitable tuning of the IMF mass limits, post-supernova fallback of metals onto black holes or metal outflows must then be invoked, to reproduce the observed chemical properties of disc galaxies. In the appendix we provide M/L-colour relations to estimate the stellar M/L ratio of a galaxy on the base of its colours, for several IMFs.

For the Warm Dark Matter (WDM) cosmological model the implications of strongly inhomogenous, primordial baryon distribution on sub-galactic scales for Big Bang Nucleosynthesis, Cosmic Microwave Background anisotropies and Galaxy Formation (including fully non-linear evolution to z=0) are discussed, and the inflationary theory leading to such distributions is briefly reviewed. It is found that Big Bang Nucleosynthesis is essentially unaffected relative to SBBN and that the change in recombination history at z~1500-700 relative to ``standard'' theory leads to differences in the anisotropy and polarization power spectra, which should be detectable by the Planck satellite provided systematic effects can be accounted for. Moreover, it is shown by fully cosmological, hydro/gravity simulations that the formation of galactic discs is only weakly affected by going from smooth to highly non-homogenous, initial baryon distributions. In particular, the final disc angular momenta at z=0 are as large as for the standard case and the ``disc angular momentum problem'' is solved to within a factor of two or better without invoking (hypothetical) energetic feedback events. A very desirable difference relative to the the standard WDM model, however, is that the on-set of star (and AGN) formation happens earlier. For the ``optimal'' free-streaming mass scale of M_f~1.5x10^11 h^-1 Msun the redshift of formation of the first stars increases from z_*=4-5 to >6.5, in much better agreement with observational data on high-redshift galaxies and QSOs. It will, however, not be possible to push z_* above ~10, because at higher redshifts the gas velocity field is nowhere compressive. Probing the ``dark ages'' will hence enable a direct test of this theory.

We measure the average gravitational shear profile of 6 massive clusters (M_vir ~ 10^15 M_sun) at z=0.3 out to a radius ~2h^-1 Mpc. The measurements are fitted to a generalized NFW-like halo model \rho(r) with an arbitrary r -> 0 slope \alpha. The data are well fitted by such a model with a central cusp with \alpha ~ 0.9 - 1.6 (68% confidence interval). For the standard-NFW case \alpha = 1.0, we find a concentration parameter c_vir that is consistent with recent predictions from high-resolution CDM N-body simulations. Our data are also well fitted by an isothermal sphere model with a softened core. For this model, our 1\sigma upper limit for the core radius corresponds to a limit \sigma_star \leq 0.1 cm^2 g^-1 on the elastic collision cross-section in a self-interacting dark matter model.

We present warm dark matter (WDM) as a possible solution to the missing satellites and angular momentum problem in galaxy formation and introduce improved initial conditions for numerical simulations of WDM models, which avoid the formation of unphysical haloes found in earlier simulations. There is a hint, that the mass function of satellite haloes has been overestimated so far, pointing to higher values for the WDM particle mass.

The X-ray properties of the haloes of disk galaxies formed in fully cosmological, hydro/gravity simulations are discussed. The results are found to be consistent with observational X-ray detections and upper limits. Disk galaxy haloes are predicted to be about an order of magnitude brighter in soft X-rays at z~1 than at z=0. The Ly-alpha and H-alpha surface brightness of an edge-on, Milky Way like model galaxy has been determined. The emission is found to be quite extended, with a scale height of about 600 pc, neglecting extinction corrections.

A low mass-to-light ratio for the baryonic component of spiral galaxies is advocated by a number of dynamical studies and by cosmological simulations of galaxy formation. We discuss the possibility of obtaining low mass-to-light ratios for the stellar component in discs, by changing the Initial Mass Function and the Star Formation History.

We have performed TreeSPH simulations of galaxy formation in a standard LCDM cosmology, including effects of star formation, energetic stellar feedback processes and a meta-galactic UV field, and obtain a mix of disk, lenticular and elliptical galaxies. The disk galaxies are deficient in angular momentum by only about a factor of two compared to observed disk galaxies. The stellar disks have approximately exponential surface density profiles, and those of the bulges range from exponential to r^{1/4}, as observed. The bulge-to-disk ratios of the disk galaxies are consistent with observations and likewise are their integrated B-V colours, which have been calculated using stellar population synthesis techniques. Furthermore, we can match the observed I-band Tully-Fisher (TF) relation, provided that the mass-to-light ratio of disk galaxies, (M/L_I), is about 0.8. The ellipticals and lenticulars have approximately r^{1/4} stellar surface density profiles, are dominated by non-disklike kinematics and flattened due to non-isotropic stellar velocity distributions, again consistent with observations.

Bolometric and 0.2-2 keV X-ray luminosities of the hot gas haloes of simulated disc galaxies have been calculated at redshift z=0. The TreeSPH simulations are fully cosmological and the sample of 44 disc galaxies span a range in characteristic circular speeds of V_c = 130-325 km/s. The galaxies have been obtained in simulations with a considerable range of physical parameters, varying the baryonic fraction, the gas metallicity, the meta-galactic UV field, the cosmology, the dark matter type, and also the numerical resolution. The models are found to be in agreement with the (few) relevant X-ray observations available at present. The amount of hot gas in the haloes is also consistent with constraints from pulsar dispersion measures in the Milky Way. Forthcoming XMM and Chandra observations should enable much more stringent tests and provide constraints on the physical parameters. We find that simple cooling flow models over-predict X-ray luminosities by up to two orders of magnitude for high (but still realistic) cooling efficiencies relative to the models presented here. Our results display a clear trend that increasing cooling efficiency leads to decreasing X-ray luminosities at z=0. The reason is found to be that increased cooling efficiency leads to a decreased fraction of hot gas relative to total baryonic mass inside of the virial radius at present. At gas metal abundances of a third solar this hot gas fraction becomes as low as just a few percent. We also find that most of the X-ray emission comes from the inner parts (inner about 20 kpc) of the hot galactic haloes. Finally, we find for realistic choices of the physical parameters that disc galaxy haloes possibly were more than one order of magnitude brighter in soft X-ray emission at z=1, than at present.

Two ways of possibly solving the angular momentum problem plaguing cold dark matter (CDM) ab initio simulations of disk galaxy formation are discussed: 1) Stellar feedback processes and 2) Warm dark matter (WDM) rather than CDM. In relation to the chemical evolution of disk galaxies our simulations indicate that in case 1) the first generation of disk stars formed in disk galaxies like the Milky Way should have an abundance about two dex below solar, in fairly good agreement with the lowest observed abundance of the metal-weak tail of the Galactic thick disk. For the second case no such statements can be made without further assumptions about the star-formation history of the galaxies. We find that the I-band Tully-Fisher relation can be matched by WDM disk galaxy formation simulations provided (M/L_I) is about 0.8 for disk galaxies, which Sommer-Larsen & Dolgov (1999) argue is a reasonable value. Finally it is discussed how the magnetic field strengths observed in galactic disks can be obtained through disk galaxy formation, as an alternative to the conventional dynamo hypothesis.

A few, simple and qualitative examples of the potential of galaxy clusters as diagnostics of cosmology are presented. In relation to these we discuss briefly three ongoing or forthcoming cluster surveys in the optical/NIR, X-rays and the cosmic microwave radiation background.

We point out that polarization measurements as planned for the upcoming PLANCK mission can significantly enhance the accuracy of cosmic parameter estimation compared to the temperature anisotropy spectrum alone. In order to illustrate this, we consider a standard cosmological model and several modifications that adjust one parameter each to fit the recently published Maxima-1 data. While all models produce acceptable fits as far as the power spectrum is concerned, their corresponding polarization spectra differ widely. The strongest differences are expected for a model with delayed recombination, reflecting the fact that polarization measurements are most sensitive to the processes governing the epoch of recombination.

We use the recent CMB power spectrum measurement by the Maxima experiment (Hanany et al. 2000) to constrain the redshift of reionization z_re. This becomes possible by combining the CMB data with cosmological parameters from various independent measurements, including Big Bang Nucleosynthesis (BBN) and X-ray cluster data. Most notably, our results provide a robust lower bound on z_re. We find that z_re>15 (8) at the 68% (95%) confidence level, unless the Hubble constant is larger than 75 km/s/Mpc.

A numerical double burst model of the chemical evolution of gas-rich dwarf galaxies has been developed. The model is fitted to a sample of N/O, O/H, Y and gas fraction observations, where N/O and O/H are the relative abundances by number of nitrogen to oxygen and oxygen to hydrogen, respectively. Y is the abundance by mass of helium. Closed models as well as models including enriched outflow, ordinary outflow and ordinary outflow combined with inflow are considered.The bursts are assumed to be instantaneous but ordered in pairs to explain the scatter in N/O-O/H. The method of gas fraction fitting is revised, and it is found that it is very important to specify whether dwarf irregulars (dIrrs) or blue compact galaxies (BCGs) are considered. Effective enriched winds fail when fitting N/O, whereas closed models, models with ordinary winds or a combination of ordinary winds and inflow are all viable.

Ibata et al (1999) have recently discovered faint, moving objects in the Hubble Deep Field. The quantity, magnitudes and proper motions of these objects are consistent with old white dwarfs making up the Galactic dark halo. We review a number of ground-based proper motion surveys in which nearby dark halo white dwarfs might be present, if they have the colours and absolute magnitudes proposed. No such objects have been found, whereas we argue here that several times more would be expected than in the Hubble Deep Field. We conclude it is unlikely that hydrogen atmosphere white dwarfs make up a significant fraction of the halo dark matter. No limits can be placed yet on helium atmosphere dwarfs from optical searches.

The hypothesis that the Damped Ly-alpha systems (DLAs) are large, galactic disks (Milky Way sized) is tested by confronting predictions of models of the formation and evolution of (large) disk galaxies with observations, in particular the Zinc abundance distribution with neutral hydrogen column density found for DLAs. A pronounced mismatch is found strongly hinting that the majority of DLAs may not be large, galactic disks.

When only cooling processes are included, smoothed-particle hydrodynamical simulations of galaxy formation in a Cold Dark Matter hierarchical clustering scenario consistently produce collapsed objects that are deficient in angular momentum by a factor of about 25 relative to the disks of observed spiral galaxies. It is widely hoped that proper allowance for star formation feedback effects will resolve this discrepancy. We explore and compare the effects of including two different types of feedback event: uniform reheating of the entire Universe to 5e5 K at a redshift z~6, and gas blow-out from pregalactic gas clouds comparable to present-day dwarf galaxies at a somewhat lower redshift. We find that blow-out is far more successful than early reheating, and that it may even be sufficient to solve the angular momentum problem. We see indications that the remaining angular momentum deficit (a factor of 5 in the blow-out models) is due to current limitations in the numerical method (SPH) used in our and other authors' work. Our most successful models distinguish themselves by the fact that a large fraction of the gas is accreted very gradually, in a cooling flow from a surrounding hot phase, rather than by mergers of massive cold clumps.

The structure, kinematics and dynamics of the Galactic stellar halo are reviewed including evidence of substructure in the spatial distribution and kinematics of halo stars. Implications for galaxy formation theory are subsequently discussed; in particular it is argued that the observed kinematics of stars in the outer Galactic halo can be used as an important constraint on viable galaxy formation scenarios.

We investigate the claim by Burkert and Silk (1997) that the observed rotation curve of the dwarf irregular galaxy DDO154 cannot be reconciled with the universal CDM halo profile of Navarro, Frenk & White (1996,1997) even when allowance is made for the effect of violent gas outflow events on the structure of the galaxy. By means of N-body simulations we show that under certain conditions it is possible to obtain a reasonable fit to the observed rotation curve without invoking Burkert & Silk's proposed spheroidal MACHO component. We are able to best reproduce the observed decline in the rotation curve by postulating additional hidden disc mass, in an amount that is compatible with disc stability requirements. In the process we improve upon the results of Navarro, Eke & Frenk (1996) on the formation of halo cores by mass loss by using actual haloes from Cold Dark Matter simulations instead of Hernquist (1990) distributions.

Effects of gravitational softening on the global structure of self-gravitating disks in centrifugal equillibrium are examined in relation to hydrodynamical/gravitational simulations. The one-parameter spline softening proposed by Hernquist & Katz is used. It is found that if the characteristic size of a disk, r, is comparable or less than the gravitational softening length, epsilon, then the cross section of the simulated disk is significantly larger than that of a no-softening (Newtonian) disk with the same mass and angular momentum. We furthermore demonstrate that if r is less than or about epsilon/2 then the scaling relation r proportional to epsilon^3/4 holds for a given mass and specific angular momentum distribution with mass. Finally we compare some of the theoretical results obtained in this and a previous paper with the results of numerical Tree-SPH simulations and find qualitative agreement.

New X-ray observations of the sparse NGC 3258 group of galaxies made by the ASCA satellite with good spectral and spatial resolution has revealed that this group has a gravitational potential deep enough to prohibit significant mass removal from the system. The baryonic fraction within 240 kpc is found to be 0.065 +0.051/-0.020 for h_{50}=1, where h_{50}=H_0/50 km/s/Mpc, in good agreement with the universal value of 0.05 +/-0.01 predicted by standard Big Bang nucleosynthesis for a Universe with Omega_0=1 and h_{50}=1. Since the deep potential of the NGC 3258 group ensures that all pristine intragroup gas has been retained, the baryonic fraction of the NGC 3258 group is indicative of the universal value. Consequently it seems premature to rule out a critical Universe.

We study the accuracy of inference of the massive halo objects' (MHO or `Macho') mass function from microlensing events observed toward LMC. Assuming the spatial distribution and kinematics of the objects are known, the slope and the range of the MHO mass function (modeled here by a simple power law) will be possible to determine from 100-1000 detected events if the slope \alpha is in the range [-2.5, -0.5], with the statistical errors reaching their minima at \alpha =-1.5. Outside this range the errors grow rapidly making the inference difficult even at very large numbers of events (N\approx 10000). On the other hand, the average mass of the MHOs will be determined to better than about 30% accuracy from N\approx 100 events for any slope. We find that the accuracy of inference at fixed N will not be strongly affected by the presently available detection efficiencies if the typical MHO masses are in the range indicated by the events detected so far. We also estimate the effects of uncertainty of the MHOs' spatial distribution and kinematics on the determination of their mass function. The halo models considered are all spherical but we allow for various density profiles and a radius-dependent, anisotropic velocity dispersion. We find that, while the mass function slope and range are weakly affected for -2 < \alpha < 0, the error in the average mass due to the halo uncertainties could be reduced to less than 50% only through the detection of about 1000 or more events. Reliable estimates of the halo structure itself [density profile and velocity dispersion profile] can start only at very large numbers of detections (N ~ 10000).

We re-analyse the kinematics of the system of blue horizontal branch field (BHBF) stars in the Galactic halo (in particular the outer halo), fitting the kinematics with the model of radial and tangential velocity dispersions in the halo as a function of galactocentric distance r proposed by Sommer-Larsen, Flynn & Christensen (1994), using a much larger sample (almost 700) of BHBF stars. The basic result is that the character of the stellar halo velocity ellipsoid changes markedly from radial anisotropy at the sun to tangential anisotropy in the outer parts of the Galactic halo (r greater than approx 20 kpc). Specifically, the radial component of the stellar halo's velocity ellipsoid decreases fairly rapidly beyond the solar circle, from approx 140 +/- 10 km/s at the sun, to an asymptotic value of 89 +/- 19 km/s at large r. The rapid decrease in the radial velocity dispersion is matched by an increase in the tangential velocity dispersion, with increasing r. Our results may indicate that the Galaxy formed hierarchically (partly or fully) through merging of smaller subsystems - the 'bottom-up' galaxy formation scenario, which for quite a while has been favoured by most theorists and recently also has been given some observational credibility by HST observations of a potential group of small galaxies, at high redshift, possibly in the process of merging to a larger galaxy (Pascarelle et al 1996).

A theory for the structure of isothermal, self-gravitating gas spheres in pressure equilibrium in a softened gravitational field is developed. The one parameter spline softening proposed by Hernquist & Katz (1989) is used. We show that the addition of this extra scale parameter implies that the set of equilibrium solutions constitute a one-parameter family, rather than the one and only one isothermal sphere solution for Newtonian gravity. We demonstrate the perhaps somewhat surprising result that for any finite choice of softening length and temperature, it is possible to deposit an arbitrarily large mass of gas in pressure equilibrium and with a non-singular density distribution inside of r_0 for any r_0 > 0. The theoretical predictions of our models are compared with the properties of the small, massive, quasi-isothermal gas clumps which typically form in numerical Tree-SPH simulations of 'passive' galaxy formation of Milky Way sized galaxies. We find reasonable agreement despite the neglect of rotational support in the models. We comment on whether the hydrodynamical resolution in our numerical simulation of galaxy formation is sufficient, and finally we conclude that one should be cautious, when comparing results of numerical simulations involving gravitational softening and hydrodynamical smoothing, with reality.

We have tested whether the simple model for the kinematics of the Galactic stellar halo (in particular the outer halo) proposed by Sommer-Larsen, Flynn and Christensen (SLFC) is physically realizable, by directly integrating particles in a 3-D model of the Galactic potential. We are able to show that the SLFC solution can be realized in terms of a distribution of particles with stationary statistical properties in phase-space. Hence, the SLFC model, which shows a notable change in the anisotropy from markedly radial at the sun to markedly tangential beyond about Galactocentric radius r=20 kpc, seems a tenable description of outer halo kinematics.