The effect of dwarf galaxy disruption in semi-analytic models

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

Galactic mergers and gravitationally unbound populations

Motivated by the observations on the intracluster light and intergalactic stellar populations, N -body simulations are used to model the galactic merging events as a goal to investigate the production and distribution of gravitational unbound populations (GUPs). Both the parabolic and hyperbolic mergers are considered, and each category includes six models with different relative orientations between two galaxies. Our results show that there are more (about a factor of 2) GUPs after a hyperbolic merging event than after a parabolic one. In general, depending on the relative orientation and also on the relative velocity of the two galaxies in a merging pair, a head-on collision of a galaxy pair would only make a tiny fraction (less than 1 per cent) of the initial stellar mass luminous GUP, but a considerable fraction (8–14 per cent) of the dark matter becomes dark GUP.     

Galaxy cluster masses without non-baryonic dark matter

We apply the modified acceleration law obtained from Einstein gravity coupled to a massive skew symmetric field,   F _μνλ  , to the problem of explaining X-ray galaxy cluster masses without exotic dark matter. Utilizing X-ray observations to fit the gas mass profile and temperature profile of the hot intracluster medium (ICM) with King 'β-models', we show that the dynamical masses of the galaxy clusters resulting from our modified acceleration law fit the cluster gas masses for our sample of 106 clusters without the need of introducing a non-baryonic dark matter component. We are further able to show for our sample of 106 clusters that the distribution of gas in the ICM as a function of radial distance is well fitted by the dynamical mass distribution arising from our modified acceleration law without any additional dark matter component. In a previous work, we applied this theory to galaxy rotation curves and demonstrated good fits to our sample of 101 low surface brightness, high surface brightness and dwarf galaxies including 58 galaxies that were fitted photometrically with the single-parameter mass-to-light ratio ( M / L )_stars. The results obtained there were qualitatively similar to those obtained using Milgrom's phenomenological Modified Newtonian Dynamics (MOND) model, although the determined galaxy masses were quantitatively different, and MOND does not show a return to Keplerian behaviour at extragalactic distances. The results obtained here are compared to those obtained using Milgrom's phenomenological MOND model which does not fit the X-ray galaxy cluster masses unless an auxiliary dark matter component is included.     

Regularized orbit models unveiling the stellar structure and dark matter halo of the Coma elliptical NGC 4807

This is the second in a series of papers dedicated to unveiling the mass structure and orbital content of a sample of flattened early-type galaxies in the Coma cluster. The ability of our orbit libraries to reconstruct internal stellar motions and the mass composition of a typical elliptical in the sample is investigated by means of Monte Carlo simulations of isotropic rotator models. The simulations allow a determination of the optimal amount of regularization needed in the orbit superpositions. It is shown that under realistic observational conditions and with the appropriate regularization, internal velocity moments can be reconstructed to an accuracy of ≈15 per cent; the same accuracy can be achieved for the circular velocity and dark matter fraction. In contrast, the flattening of the halo remains unconstrained. Regularized orbit superpositions are applied to a first galaxy in our sample, NGC 4807, for which stellar kinematical observations extend to  3  r _eff  . The galaxy seems dark-matter dominated outside   r > 2  r _eff  . Logarithmic dark matter potentials are consistent with the data, as well as NFW profiles, mimicking logarithmic potentials over the observationally sampled radial range. In both cases, the derived stellar mass-to-light ratio ϒ agrees well with independently obtained mass-to-light ratios from stellar population analysis. The achieved accuracy is  Δϒ≈ 0.5  . Kinematically, NGC 4807 is characterized by mild radial anisotropy outside   r > 0.5  r _eff  , becoming isotropic towards the centre. Our orbit models hint at either a distinct stellar component or weak triaxiality in the outer parts of the galaxy.     

Mass distribution in nearby Abell clusters

We study the mass distribution in six nearby  ( z < 0.06)  relaxed Abell clusters of galaxies A0262, A0496, A1060, A2199, A3158 and A3558. Given the dominance of dark matter in galaxy clusters, we approximate their total density distribution by the Navarro, Frenk & White (NFW) formula characterized by virial mass and concentration. We also assume that the anisotropy of galactic orbits is reasonably well described by a constant and that galaxy distribution traces that of the total density. Using the velocity and position data for 120–420 galaxies per cluster we calculate, after removal of interlopers, the profiles of the lowest order even velocity moments, dispersion and kurtosis. We then reproduce the velocity moments by jointly fitting the moments to the solutions of the Jeans equations. Including the kurtosis in the analysis allows us to break the degeneracy between the mass distribution and anisotropy and constrain the anisotropy as well as the virial mass and concentration. The method is tested in detail on mock data extracted from the N -body simulations of dark matter haloes. We find that the best-fitting Galactic orbits are remarkably close to isotropic in most clusters. Using the fitted pairs of mass and concentration parameters for the six clusters, we conclude that the trend of decreasing concentration for higher masses found in the cosmological N -body simulations is consistent with the data. By scaling the individual cluster data by mass, we combine them to create a composite cluster with 1465 galaxies and perform a similar analysis on such sample. The estimated concentration parameter then lies in the range  1.5 < c < 14  and the anisotropy parameter in the range  −1.1 < β < 0.5  at the 95 per cent confidence level.     

Dynamical friction and galaxy merging time-scales

The time-scale for galaxies within merging dark matter haloes to merge with each other is an important ingredient in galaxy formation models. Accurate estimates of merging time-scales are required for predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar mass in central galaxies and the statistical properties of satellite galaxies within dark matter haloes. In this paper, we study the merging time-scales of extended dark matter haloes using N -body simulations. We compare these results to standard estimates based on the Chandrasekhar theory of dynamical friction. We find that these standard predictions for merging time-scales, which are often used in semi-analytic galaxy formation models, are systematically shorter than those found in simulations. The discrepancy is approximately a factor of 1.7 for M _sat/ M _host≈ 0.1 and becomes larger for more disparate satellite-to-host mass ratios, reaching a factor of ∼3.3 for M _sat/ M _host≈ 0.01. Based on our simulations, we propose a new, easily implementable fitting formula that accurately predicts the time-scale for an extended satellite to sink from the virial radius of a host halo down to the halo's centre for a wide range of M _sat/ M _host and orbits. Including a central bulge in each galaxy changes the merging time-scale by ≲10 per cent. To highlight one concrete application of our results, we show that merging time-scales often used in the literature overestimate the growth of stellar mass by satellite accretion by ≈40 per cent, with the extra mass gained in low mass ratio mergers.     

Populating a cluster of galaxies – I. Results at

We simulate the assembly of a massive rich cluster and the formation of its constituent galaxies in a flat, low-density universe. Our most accurate model follows the collapse, the star formation history and the orbital motion of all galaxies more luminous than the Fornax dwarf spheroidal, while dark halo structure is tracked consistently throughout the cluster for all galaxies more luminous than the SMC. Within its virial radius this model contains about     dark matter particles and almost 5000 distinct dynamically resolved galaxies. Simulations of this same cluster at a variety of resolutions allow us to check explicitly for numerical convergence both of the dark matter structures produced by our new parallel N -body and substructure identification codes, and of the galaxy populations produced by the phenomenological models we use to follow cooling, star formation, feedback and stellar aging. This baryonic modelling is tuned so that our simulations reproduce the observed properties of isolated spirals outside clusters. Without further parameter adjustment our simulations then produce a luminosity function, a mass-to-light ratio, luminosity, number and velocity dispersion profiles, and a morphology–radius relation which are similar to those observed in real clusters. In particular, since our simulations follow galaxy merging explicitly, we can demonstrate that it accounts quantitatively for the observed cluster population of bulges and elliptical galaxies.     

Polarization of resonance X-ray lines from clusters of galaxies

It is known that resonant scattering can distort the surface brightness profiles of clusters of galaxies in X-ray lines. We demonstrate that the scattered line emission should be polarized and possibly detectable with future X-ray polarimeters. Spectrally resolved mapping of a galaxy cluster in polarized X-rays could provide valuable independent information on the physical conditions, in particular element abundances and the characteristic velocity of small-scale turbulent motions, in the intracluster gas. The expected degree of polarization is of the order of 10 per cent for the richest regular clusters (e.g. Coma) and clusters whose X-ray emission is dominated by a central cooling flow (such as Perseus and M87/Virgo).     

Modelling gaseous and stellar kinematics in the disc galaxies NGC 772, 3898 and 7782

We present V -band surface photometry and major-axis kinematics of stars and ionized gas of three early-type spiral galaxies, namely NGC 772, 3898 and 7782. For each galaxy we present a self-consistent Jeans model for the stellar kinematics, adopting the light distribution of bulge and disc derived by means of a two-dimensional parametric photometric decomposition. This allows us to investigate the presence of non-circular gas motions, and derive the mass distribution of luminous and dark matter in these objects.
NGC 772 and 7782 have apparently normal kinematics with the ionized gas tracing the gravitational equilibrium circular speed. This is not true in the innermost region (| r |≲8 arcsec) of NGC 3898, where the ionized gas is rotating more slowly than the circular velocity predicted by dynamical modelling. This phenomenon is common in the bulge-dominated galaxies for which dynamical modelling enables us to make the direct comparison between the gas velocity and the circular speed, and it poses questions about the reliability of galaxy mass distributions derived by the direct decomposition of the observed ionized-gas rotation curve into the contributions of luminous and dark matter.     

The effects of photoionization on galaxy formation – I. Model and results at z=0

We develop a coupled model for the evolution of the global properties of the intergalactic medium (IGM) and the formation of galaxies, in the presence of a photoionizing background due to stars and quasars. We use this model to predict the thermodynamic history of the IGM when photoionized by galaxies forming in a cold dark matter (CDM) universe. The evolution of the galaxies is calculated using a semi-analytical model, including a detailed treatment of the effects of tidal stripping and dynamical friction on satellite galaxies orbiting inside larger dark matter haloes. We include in the model the negative feedback on galaxy formation from the photoionizing background. Photoionization inhibits galaxy formation in low-mass dark matter haloes in two ways: (i) heating of the IGM and inhibition of the collapse of gas into dark haloes by the IGM pressure, and (ii) reduction in the rate of radiative cooling of gas within haloes. The result of our method is a self-consistent model of galaxy formation and the IGM. The IGM is reheated twice (during reionization of H  i and He  ii ), and we find that the star formation rate per unit volume is slightly suppressed after each episode of reheating. We find that galaxies brighter than L _★ are mostly unaffected by reionization, while the abundance of faint galaxies is significantly reduced, leading to present-day galaxy luminosity functions with shallow faint-end slopes, in good agreement with recent observational data. Reionization also affects other properties of these faint galaxies, in a readily understandable way.     

Numerical simulations of hot halo gas in galaxy mergers

Galaxy merger simulations have explored the behaviour of gas within the galactic disc, yet the dynamics of hot gas within the galaxy halo have been neglected. We report on the results of high-resolution hydrodynamic simulations of colliding galaxies with metal-free hot halo gas. To isolate the effect of the halo gas, we simulate only the dark matter halo and the hot halo gas over a range of mass ratios, gas fractions and orbital configurations to constrain the shocks and gas dynamics within the progenitor haloes. We find that (i) a strong shock is produced in the galaxy haloes before the first passage, increasing the temperature of the gas by almost an order of magnitude to   T ∼ 10^6.3 K  . (ii) The X-ray luminosity of the shock is strongly dependent on the gas fraction; it is  ≳10³⁹ erg s⁻¹  for halo gas fractions larger than 10 per cent. (iii) The hot diffuse gas in the simulation produces X-ray luminosities as large as  10⁴² erg s⁻¹  . This contributes to the total X-ray background in the Universe. (iv) We find an analytic fit to the maximum X-ray luminosity of the shock as a function of merger parameters. This fit can be used in semi-analytic recipes of galaxy formation to estimate the total X-ray emission from shocks in merging galaxies. (v) ∼10–20 per cent of the initial gas mass is unbound from the galaxies for equal-mass mergers, while 3–5 per cent of the gas mass is released for the 3:1 and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible mechanism to enrich the intergalactic medium with metals.     

Ram pressure stripping the hot gaseous haloes of galaxies in groups and clusters

We use a large suite of carefully controlled full hydrodynamic simulations to study the ram pressure stripping of the hot gaseous haloes of galaxies as they fall into massive groups and clusters. The sensitivity of the results to the orbit, total galaxy mass, and galaxy structural properties is explored. For typical structural and orbital parameters, we find that ∼30 per cent of the initial hot galactic halo gas can remain in place after 10 Gyr. We propose a physically simple analytic model that describes the stripping seen in the simulations remarkably well. The model is analogous to the original formulation of Gunn & Gott, except that it is appropriate for the case of a spherical (hot) gas distribution (as opposed to a face-on cold disc) and takes into account that stripping is not instantaneous but occurs on a characteristic time-scale. The model reproduces the results of the simulations to within ≈10 per cent at almost all times for all the orbits, mass ratios, and galaxy structural properties we have explored. The one exception involves unlikely systems where the orbit of the galaxy is highly non-radial and its mass exceeds about 10 per cent of the group or cluster into which it is falling (in which case the model underpredicts the stripping following pericentric passage). The proposed model has several interesting applications, including modelling the ram pressure stripping of both observed and cosmologically simulated galaxies and as a way to improve present semi-analytic models of galaxy formation. One immediate consequence is that the colours and morphologies of satellite galaxies in groups and clusters will differ significantly from those predicted with the standard assumption of complete stripping of the hot coronae.     

Ram pressure stripping of spiral galaxies in clusters

We use three-dimensional SPH/ N -body simulations to study ram pressure stripping of gas from spiral galaxies orbiting in clusters. We find that the analytic expectation of Gunn & Gott, relating the gravitational restoring force provided by the disc to the ram pressure force, provides a good approximation to the radius at which gas will be stripped from a galaxy. However, at small radii it is also important to consider the potential provided by the bulge component. A spiral galaxy passing through the core of a rich cluster, such as Coma, will have its gaseous disc truncated to ∼4 kpc, thus losing ∼80 per cent of its diffuse gas mass. The time-scale for this to occur is a fraction of a crossing time ∼10⁷ yr. Galaxies orbiting within poorer clusters, or inclined to the direction of motion through the intracluster medium, will lose significantly less gas. We conclude that ram pressure alone is insufficient to account for the rapid and widespread truncation of star formation observed in cluster galaxies, or the morphological transformation of Sabs to S0s that is necessary to explain the Butcher–Oemler effect.     

Galaxies in clusters: the observational characteristics of bow shocks, wakes and tails

The dynamical signatures of the interaction between galaxies in clusters and the intracluster medium (ICM) can potentially yield significant information about the structure and dynamical history of clusters. To develop our understanding of this phenomenon we present results from numerical modelling of the galaxy–ICM interaction, as the galaxy moves through the cluster. The simulations have been performed for a broad range of ICM temperatures ( kT _cl=1, 4 and 8 keV), representative of poor clusters or groups through to rich clusters.
There are several dynamical features that can be identified in these simulations. For supersonic galaxy motion, a leading bow shock is present, and also a weak gravitationally focused wake or tail behind the galaxy (analogous to Bondi–Hoyle accretion). For galaxies with higher mass replenishment rates and a denser interstellar medium (ISM), the dominant feature is a dense ram-pressure stripped tail. In line with other simulations, we find that the ICM/galaxy–ISM interaction can result in complex time-dependent dynamics, with ram-pressure stripping occurring in an episodic manner.
In order to facilitate this comparison between the observational consequences of dynamical studies and X-ray observations we have calculated synthetic X-ray flux and hardness maps from these simulations. These calculations predict that the ram-pressure stripped tail will usually be the most visible feature, though in nearby galaxies the bow shock preceding the galaxy should also be apparent in deeper X-ray observations. We briefly discuss these results and compare them with X-ray observations of galaxies where there is evidence of such interactions.     

Galaxy–dark matter correlations applied to galaxy–galaxy lensing: predictions from the semi-analytic galaxy formation models

We use semi-analytic models of galaxy formation combined with high-resolution N -body simulations to make predictions for galaxy–dark matter correlations and apply them to galaxy–galaxy lensing. We analyse cross-power spectra between the dark matter and different galaxy samples selected by luminosity, colour or star formation rate. We compare the predictions with the recent detection by the Sloan Digital Sky Survey (SDSS). We show that the correlation amplitude and the mean tangential shear depend strongly on the luminosity of the sample on scales below 1  h ⁻¹ Mpc, reflecting the correlation between the galaxy luminosity and the halo mass. The cross-correlation cannot, however, be used to infer the halo profile directly because different halo masses dominate on different scales and because not all galaxies are at the centres of the corresponding haloes. We compute the redshift evolution of the cross-correlation amplitude and compare it with those of galaxies and dark matter. We also compute the galaxy–dark matter correlation coefficient and show that it is close to unity on scales above 1  h ⁻¹ Mpc for all considered galaxy types. This would allow one to extract the bias and the dark matter power spectrum on large scales from the galaxy and galaxy–dark matter correlations.     

Ram-pressure histories of cluster galaxies

Ram-pressure stripping can remove significant amounts of gas from galaxies that orbit in clusters and massive groups, and thus has a large impact on the evolution of cluster galaxies. In this paper, we reconstruct the present-day distribution of ram pressure and the ram-pressure histories of cluster galaxies. To this aim, we combine the Millennium Simulation and an associated semi-analytic model of galaxy evolution with analytic models for the gas distribution in clusters. We find that about one quarter of galaxies in massive clusters are subject to strong ram pressures that are likely to cause an expedient loss of all gas. Strong ram pressures occur predominantly in the inner core of the cluster, where both the gas density and the galaxy velocity are higher. Since their accretion on to a massive system, more than 64 per cent of galaxies that reside in a cluster today have experienced strong ram pressures of  >10⁻¹¹ dyn cm⁻²  which most likely led to a substantial loss of the gas.     

Adiabatic scaling relations of galaxy clusters

The aim of this work is to show that, contrary to popular belief, galaxy clusters are not expected to be self-similar, even when the only energy sources available are gravity and shock-wave heating. In particular, we investigate the scaling relations between mass, luminosity and temperature of galaxy groups and clusters in the absence of radiative processes. Theoretical expectations are derived from a polytropic model of the intracluster medium and compared with the results of high-resolution adiabatic gasdynamical simulations. It is shown that, in addition to the well-known relation between the mass and concentration of the dark matter halo, the effective polytropic index of the gas also varies systematically with cluster mass, and therefore neither the dark matter nor the gas profiles are exactly self-similar. It is remarkable, though, that the effects of concentration and polytropic index tend to cancel each other, leading to scaling relations whose logarithmic slopes roughly match the predictions of the most-basic self-similar models. We provide a phenomenological fit to the relation between polytropic index and concentration, as well as a self-consistent scheme to derive the non-linear scaling relations expected for any cosmology and the best-fitting normalizations of the M – T , L – T and F – T relations appropriate for a Λ cold dark matter universe. The predicted scaling relations reproduce observational data reasonably well for massive clusters, where the effects of cooling and star formation are expected to play a minor role.     

Bulges versus discs: the evolution of angular momentum in cosmological simulations of galaxy formation

We investigate the evolution of angular momentum in simulations of galaxy formation in a cold dark matter universe. We analyse two model galaxies generated in the N -body/hydrodynamic simulations of Okamoto et al. Starting from identical initial conditions, but using different assumptions for the baryonic physics, one of the simulations produced a bulge-dominated galaxy and the other one a disc-dominated galaxy. The main difference is the treatment of star formation and feedback, both of which were designed to be more efficient in the disc-dominated object. We find that the specific angular momentum of the disc-dominated galaxy tracks the evolution of the angular momentum of the dark matter halo very closely: the angular momentum grows as predicted by linear theory until the epoch of maximum expansion and remains constant thereafter. By contrast, the evolution of the angular momentum of the bulge-dominated galaxy resembles that of the central, most bound halo material: it also grows at first according to linear theory, but 90 per cent of it is rapidly lost as pre-galactic fragments, into which gas had cooled efficiently, merge, transferring their orbital angular momentum to the outer halo by tidal effects. The disc-dominated galaxy avoids this fate because the strong feedback reheats the gas, which accumulates in an extended hot reservoir and only begins to cool once the merging activity has subsided. Our analysis lends strong support to the classical theory of disc formation whereby tidally torqued gas is accreted into the centre of the halo conserving its angular momentum.     

Effects of cluster galaxies on arc statistics

We present the results of a set of numerical simulations evaluating the effect of cluster galaxies on arc statistics.
We perform a first set of gravitational lensing simulations using three independent projections for each of nine different galaxy clusters obtained from N -body simulations. The simulated clusters consist of dark matter only. We add a population of galaxies to each cluster, mimicking the observed luminosity function and the spatial galaxy distribution, and repeat the lensing simulations including the effects of cluster galaxies, which themselves act as individual lenses. Each galaxy is represented by a spherical Navarro, Frenk & White density profile.
We consider the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of lengths, widths, curvature radii and length-to-width ratios of long arcs. We find some changes to the distribution of short-arc properties in the presence of cluster galaxies. The differences appear in the distribution of curvature radii for arc lengths smaller than 12 arcsec, while the distributions of lengths, widths and length-to-width ratios are significantly changed only for arcs shorter than 4 arcsec.