We perform N-body simulations on some of the most massive galaxies extracted from a cosmological simulation of hierarchical structure formation with total masses in the range $10^{12} M_{\odot} < M_{tot} < 3\times 10^{13} M_{\odot}$ from $4\geq z \geq 0$. After galactic mergers, we track the dynamical evolution of the infalling black holes (BHs) around their host's central BHs. From 11 different simulations, we find that, of the 86 infalling BHs with masses > $10^4 M_{\odot}$, 36 merge with their host's central BH, 13 are ejected from their host galaxy, and 37 are still orbiting at $z=0$. Across all galaxies, 33 BHs are kicked to a higher orbit after close interactions with the central BH binary or multiple, after which only one of them merged with their hosts. These orbiting BHs should be detectable by their anomalous (not Low Mass X-ray Binary) spectra. The X-ray luminosities of the orbiting massive BHs at z=0 are in the range $10^{28}-10^{43}$ $\mathrm{erg}~\mathrm{s}^{-1}$, with a currently undetectable median value of $10^{33}$ $\mathrm{erg}~\mathrm{s}^{-1}$. However, the most luminous $\sim$5\% should be detectable by existing X-ray facilities.

Simulating thin and extended galactic disks has long been a challenge in computational astrophysics. We introduce the NIHAO-UHD suite of cosmological hydrodynamical simulations of Milky Way mass galaxies and study stellar disk properties such as stellar mass, size and rotation velocity which agree well with observations of the Milky Way and local galaxies. In particular, the simulations reproduce the age-velocity dispersion relation and a multi-component stellar disk as observed for the Milky Way. Half of our galaxies show a double exponential vertical profile, while the others are well described by a single exponential model which we link to the disk merger history. In all cases, mono-age populations follow a single exponential whose scale height varies monotonically with stellar age and radius. The scale length decreases with stellar age while the scale height increases. The general structure of the stellar disks is already set at time of birth as a result of the inside-out and upside-down formation. Subsequent evolution modifies this structure by increasing both the scale length and height of all mono-age populations. Thus, our results put tight constraints on how much dynamical memory stellar disks can retain over cosmological timescales. Our simulations demonstrate that it is possible to form thin galactic disks in cosmological simulations provided there are no significant stellar mergers at low redshifts. Most of the stellar mass is formed in-situ with only a few percent ($\lesssim5\%$) brought in by merging satellites at early times. Redshift zero snapshots and halo catalogues are publicly available.

What is the expected infrared output of elliptical galaxies? Here we report the latest findings obtained in this high time resolution (~10 years) and high spatial resolution (2.5 parsec at center) study. We add a set of grain physics to the MACER code, including (a) dust grains made in passive stellar evolution; (b) dust grain growth due to collision and sticking; (c) grain destruction due to thermal sputtering; (d) dust cooling of hot gas via inelastic collisions; and (e) radiation pressure on dust grains. The code improvements enable us to analyze the effects of dust on metal depletion and AGN obscuration, and also to assess the infrared output of the modeled galaxies. We simulate a representative massive elliptical galaxy of a central stellar velocity dispersion ~ 260 km/s and modest rotation. We find that: (1) the circumnuclear disk (of a size <~ 1 kpc) is dusty in its outer region where most of the metals are in dust grains, while in the inner disk most of the dust grains are destroyed by the AGN irradiation; (2) the dusty disk is optically thick to both the starlight within the disk and the radiation from the central AGN. Thus the AGN will be obscured behind the disk, and the latter is of a covering factor ~ 0.2; (3) the dust infrared emission is mainly due to the AGN irradiation. The median infrared luminosity is ~ 2e44 erg/s, and it can reach >~ 1e46 erg/s during outbursts; (4) the duty cycles of the AGN activities, star formation, and the dust infrared luminosity roughly match observations, e.g., in most of its lifetime, the simulated galaxy is a stereotypical "quiescent" elliptical galaxy with L_{IR} ~ 1e11*L_{solar}, while the star formation rate can exceed 250 M_{solar}/yr during central outbursts.

A hot plasma is the dominant phase of the interstellar medium of early-type galaxies. Its origin can reside in stellar mass losses, residual gas from the formation epoch, and accretion from outside of the galaxies. Its evolution is linked to the dynamical structure of the host galaxy, to the supernova and AGN feedback, and to (late-epoch) star formation, in a way that has yet to be fully understood. Important clues about the origin and evolution of the hot gas come from the abundances of heavy metals, that have been studied with increasing detail with XMM-Newton and Chandra. We present recent high resolution hydrodynamical simulations of the hot gas evolution that include the above processes, and where several chemical species, originating in AGB stars and supernovae of type Ia and II, have also been considered. The high resolution, of few parsecs in the central galactic region, allows us to track the metal enrichment, transportation and dilution throughout the galaxy. The comparison of model results with observed abundances reveals a good agreement for the region enriched by the AGN wind, but also discrepancies for the diffuse hot gas; the latter indicate the need for a revision of standard assumptions, and/or the importance of neglected effects as those due to the dust, and/or residual uncertainties in deriving abundances from the X-ray spectra.

We perform two-dimensional hydrodynamical simulations of slowly rotating accretion flows in the region of $ 0.01-7\, \mathrm{pc} $ around a supermassive black holes with $ M_\mathrm{BH} = 10^{8} M_{\odot} $. The accretion flow is irradiated by the photons from the central active galactic nucleus (AGN). In addition to the direct radiation from the AGN, we have also included the "re-radiation", i.e., the locally produced radiation by Thomson scattering, line and bremsstrahlung radiation. Compare to our previous work, we have improved the calculation of radiation force due to the Thomson scattering of X-ray photons from the central AGN. We find that this improvement can significantly increase the mass flux and velocity of outflow. We have compared the properties of outflow --- including mass outflow rate, velocity, and kinetic luminosity of outflow --- in our simulation with the observed properties of outflow in AGNs and found that they are in good consistency. This implies that the combination of line and re-radiation forces is the possible origin of observed outflow in luminous AGNs.

Depending on the value of the accretion rate, black hole accretion is divided into cold and hot modes. The two modes have distinctly different physics and correspond to two feedback modes. Most previous feedback works either focus only on one mode, or the accretion physics is not always properly adopted. Here, by performing hydrodynamical numerical simulations of AGN feedback in an elliptical galaxy, we show that including both is important, and gives different results from including just one or the other. We specifically focus on the wind and radiation (but neglecting the jet) feedback in the hot mode, to explore whether this particular mode of feedback can play any important feedback role. For this aim, we have run two test models. In one model, we always adopt the cold mode no matter what value the accretion rate is; in another model, we turn off the AGN once it enters into the hot mode. We have calculated the AGN light curves, black hole growth, star formation, and AGN energy duty-cycle; and compared the results with the model in which the two modes are correctly included. Important differences are found. For example, if we were to adopt only the cold mode, the total mass of newly formed stars would become two orders of magnitude smaller, and the fraction of energy ejected within the high accretion regime (i.e., $L_{\rm BH} > 2\% L_{\rm Edd}$) would be too small to be consistent with observations. We have also investigated the respective roles of wind and radiation in the hot mode.

The faintness of satellite systems in galaxy groups has contributed to the widely discussed "missing satellite" and "too big to fail" issues. Using techniques based on Tremaine & Richstone (1977), we show that there is no problem with the luminosity function computed from modern codes per se, but that the gap between first and second brightest systems is too big {\it given} the luminosity function, that the same large gap is found in modern, large scale baryonic $\Lambda$CDM simulations such as EAGLE and IllustrisTNG, is even greater in dark matter only simulations, and finally, that this is most likely due to gravitationally induced merging caused by classical dynamical friction. Quantitatively the gap is larger in the computed simulations than in the randomized ones by $1.79 \pm 1.04$, $1.51 \pm 0.93$, $3.43 \pm 1.44$ and $3.33 \pm 1.35$ magnitudes in the EAGLE, IllustrisTNG, and dark matter only simulations of EAGLE and IllustrisTNG respectively. Furthermore the anomalous gaps in the simulated systems are even larger than in the real data by over half a magnitude and are still larger in the dark matter only simulations. Briefly stated, $\Lambda$CDM does not have a problem with an absence of "too big to fail" galaxies. Statistically significant large gaps between first and second brightest galaxies are to be expected.

We explore the effect of active galactic nucleus (AGN) feedback from central galaxies on their satellites by comparing two sets of cosmological zoom-in runs of 27 halos with masses ranging from $10^{12}$ to $10^{13.4}$ solar masses at z=0, with (wAGN) and without (noAGN) AGN feedback. Both simulations include stellar feedback from multiple processes, including powerful winds from supernovae, stellar winds from young massive stars, AGB stars, radiative heating within Str\"omgren spheres and photoelectric heating. Our wAGN model is identical to the noAGN model except that it also includes a model for black hole seeding and accretion, as well as AGN feedback via high-velocity broad absorption line winds and Compton/photoionization heating. We show that the inclusion of AGN feedback from the central galaxy significantly affects the star formation history and the gas content of the satellite galaxies. AGN feedback starts to affect the gas content and the star formation of the satellites as early as z=2. The mean gas rich fraction of satellites at z=0 decreases from 15% in the noAGN simulation to 5% in the wAGN simulation. The difference between the two sets extends as far out as five times the virial radius of the central galaxy at z=1. We investigate the quenching mechanism by studying the physical conditions in the surroundings of pairs of satellites matched across the wAGN and noAGN simulations and find an increase in the temperature and relative velocity of the intergalactic gas.

We study the properties of rotating accretion flows onto supermassive black holes (SMBHs) using axisymmetric two-dimensional hydrodynamical simulations with radiative cooling and BH feedback. The simulations resolve the accretion dynamics of gas outside from the BH influence radius through an inner accretion disk. For lower Bondi accretion rates in units of the Eddington rate ($\dot{M}_{\rm B}\ll 10^{-3}~\dot{M}_{\rm Edd}$), the BH feeding is suppressed due to turbulent motion by several orders of magnitudes from the Bondi rate. Thus, the radiative luminosity results in as low as $\sim 10^{-10}-10^{-7}~L_{\rm Edd}$, where $L_{\rm Edd}$ is the Eddington luminosity. For higher rates of $\dot{M}_{\rm B}> 10^{-3}~\dot{M}_{\rm Edd}$, the optically-thin accreting gas cools via free-free emission and forms a geometrically-thin disk, which feeds the BH efficiently and increases the radiative luminosity to $> 10^{-3}~L_{\rm Edd}$. The transitional behavior of accreting BHs in galactic nuclei from radiatively inefficient phases to cold disk accretion naturally explains (1) the reason for the offset between the observed luminosities and theoretical predictions for nearby quiescent SMBHs, and (2) the conditions to fuel gas into the nuclear SMBH. In addition, the cold disk formed in galactic nuclei tends to be gravitationally unstable and leads to star formation when the Bondi rate is as high as $ \dot{M}_{\rm B} > 10^{-2}~M_\odot~{\rm yr}^{-1}$. This is a plausible explanation of the correlation observed between star formation rates and BH feeding rates in Seyfert galaxies.

We add a suite of chemical abundances to the MACER (Massive AGN Controlled Ellipticals Resolved) 2D code, by solving 12 additional continuity equations for H, He, C, N, O, Ne, Mg, Si, S, Ca, Fe and Ni respectively with sources from AGB stars and supernovae of type Ia and II with metal yields based on standard stellar physics. New stars, formed in Toomre unstable circumnuclear disks (of a size $\lesssim150$ parsec), are assumed to have a top-heavy initial mass function with a power index of 1.65. The metal dilution effects due to cosmic accretion are also included. With the high resolution of few parsecs in central regions, resolved black hole accretion and AGN feedback, we can track the metal enrichment, transportation and dilution throughout the modeled massive elliptical galaxy of velocity dispersion $\sim280$ km/s. We retrieve the chemical composition of the BAL winds launched by the central AGN, synthesize the X-ray features of the hot ISM, and find that (1) the simulated metallicity in the BAL winds could be up to $\sim 8 Z_\odot$, while that of the hot ISM in the host galaxy is $\sim 2.3 Z_\odot$, matching well with SDSS observations of BLR gas; (2) the X-ray emitting hot gas is metal enriched with a typical value $\sim2.5 Z_\odot$; (3) the circumunuclear cold gas disk, where the metals are condensed, further enriched and recycled, plays a critical role in the metal enrichment; (4) the black hole accretion rate $\dot{M}_{\rm BH}$ linearly correlates with the star formation rate $\dot{M}_\star^+$ in the circumnuclear disk, i.e, $\dot{M}_\star^+ \sim 7.7 \dot{M}_{\rm BH}$, but lagged in time by roughly $10^6$ years.

Although highly successful on cosmological scales, Cold Dark Matter (CDM) models predict unobserved over-dense `cusps' in dwarf galaxies and overestimate their formation rate. We consider an ultra-light axion-like scalar boson which promises to reduce these observational discrepancies at galactic scales. The model, known as Fuzzy Dark Matter (FDM), avoids cusps, suppresses small-scale power, and delays galaxy formation via macroscopic quantum pressure. We compare the substructure and density fluctuations of galactic dark matter haloes comprised of ultra-light axions to conventional CDM results. Besides self-gravitating subhaloes, FDM includes non-virialized over-dense wavelets formed by quantum interference patterns which are an efficient source of heating to galactic discs. We find that, in the solar neighborhood, wavelet heating is sufficient to give the oldest disc stars a velocity dispersion of $\sim 30 \: \mathrm{km} \: \mathrm{s}^{-1}$ within a Hubble time if energy is not lost from the disc, the velocity dispersion increasing with stellar age as $\sigma_D \propto t^{0.4}$ in agreement with observations. Furthermore, we calculate the radius-dependent velocity dispersion and corresponding scale height caused by the heating of this dynamical substructure in both CDM and FDM with the determination that these effects will produce a flaring that terminates the Milky Way disc at $15 - 20 \: \mathrm{kpc}$. Although the source of thickened discs is not known, the heating due to perturbations caused by dark substructure cannot exceed the total disc velocity dispersion. Therefore, this work provides a lower bound on the FDM particle mass of $m_a > 0.6 \times 10^{-22} \mathrm{eV}$. Furthermore, FDM wavelets with this particle mass should be considered a viable mechanism for producing the observed disc thickening with time.

Based on our previous modeling of AGN feedback in isolated elliptical galaxies (Gan et al. 2014) using the MACER (Massive AGN Controlled Ellipticals Resolved) code, we extend and improve the model to include rotation, to limit star formation to regions of high density and low temperature, to facilitate angular momentum transfer via the Toomre instability in gaseous disks, and to improve the treatment of hot mode (low accretion rate) AGN feedback. The model galaxy now has an extended dark matter profile that matches with standard observations, but has a resolution of parsecs in the inner region and resolves the Bondi radius. We find that the results agree reasonably well with a panoply of observations: (1) both AGN activity and star formation are primarily in central cold gaseous disks, are bursty and mainly driven by the Toomre instability; (2)AGN duty cycle agrees well with the Soltan argument, i.e., the AGN spends most of its lifetime when it is in low luminosity (half of time with $L/L_{Edd}<7\times10^{-5}$), while emitting most of its energy when it is in high luminosity (half of radiant energy emitted with $L/L_{Edd}>0.06$); (3) the total star formation is $\sim$ few percents of the initial stellar mass, occurring in the bursts that would be associated with the observed E+A phenomenon. Most of the star formation occurs in the circumnuclear disk of a size <1 kpc, which is in agreement with recent observations; (4) the ISM X-ray luminosity varies within a reasonable range (median $L_{\rm X, ISM}=9.1\times10^{39}$ erg/s) in agreement with observations.

Abundance matching between galaxies and halos has proven to be an informative technique, less dependent on detailed physical approximations than some other methods. We extend the discussion to the dwarf galaxies realm and to the study of the alternative hypothesis that the dark matter is composed of ultra--light particles: Fuzzy Dark Matter. We find that, given current observations, both CDM and FDM have difficulties with dwarf systems in the local group, but that, if mechanisms are proposed to alleviate these difficulties, they make radically different predictions for FDM and CDM at high redshift. These differences should allow clear observational tests with currently planned experiments, with the number of galaxies per unit volume with stellar mass over $10^6$ M$_{\odot}$ lower in FDM by factors of ($10^{-2.7},10^{-2.0},10^{-1.7},10^{-0.7},10^{-0.1}$) at $z=(10,8,6,4,2)$.

We study low-density axisymmetric accretion flows onto black holes (BHs) with two-dimensional hydrodynamical simulations, adopting the $\alpha$-viscosity prescription. When the gas angular momentum is low enough to form a rotationally supported disk within the Bondi radius ($R_{\rm B}$), we find a global steady accretion solution. The solution consists of a rotational equilibrium distribution at $r\sim R_{\rm B}$, where the density follows $\rho \propto (1+R_{\rm B}/r)^{3/2}$, surrounding a geometrically thick and optically thin accretion disk at the centrifugal radius, where thermal energy generated by viscosity is transported via strong convection. Physical properties of the inner solution agree with those expected in convection-dominated accretion flows (CDAF; $\rho \propto r^{-1/2}$). In the inner CDAF solution, the gas inflow rate decreases towards the center due to convection ($\dot{M}\propto r$), and the net accretion rate (including both inflows and outflows) is strongly suppressed by several orders of magnitude from the Bondi accretion rate $\dot{M}_{\rm B}$ The net accretion rate depends on the viscous strength, following $\dot{M}/\dot{M}_{\rm B}\propto (\alpha/0.01)^{0.6}$. This solution holds for low accretion rates of $\dot{M}_{\rm B}/\dot{M}_{\rm Edd}< 10^{-3}$ having minimal radiation cooling, where $\dot{M}_{\rm Edd}$ is the Eddington rate. In a hot plasma at the bottom ($r<10^{-3}~R_{\rm B}$), thermal conduction would dominate the convective energy flux. Since suppression of the accretion by convection ceases, the final BH feeding rate is found to be $\dot{M}/\dot{M}_{\rm B} \sim 10^{-3}-10^{-2}$. This rate is as low as $\dot{M}/\dot{M}_{\rm Edd} \sim 10^{-7}-10^{-6}$ inferred for SgrA$^*$ and the nuclear BHs in M31 and M87, and can explain the low luminosities in these sources, without invoking any feedback mechanism.

We use cosmological hydrodynamical simulations to investigate the role of feedback from accreting black holes on the evolution of sizes, compactness, stellar core density and specific star-formation of massive galaxies with stellar masses of $ M_{star} > 10^{10.9} M_{\odot}$. We perform two sets of cosmological zoom-in simulations of 30 halos to z=0: (1) without black holes and Active Galactic Nucleus (AGN) feedback and (2) with AGN feedback arising from winds and X-ray radiation. We find that AGN feedback can alter the stellar density distribution, reduce the core density within the central 1 kpc by 0.3 dex from z=1, and enhance the size growth of massive galaxies. We also find that galaxies simulated with AGN feedback evolve along similar tracks to those characterized by observations in specific star formation versus compactness. We confirm that AGN feedback plays an important role in transforming galaxies from blue compact galaxies into red extended galaxies in two ways: (1) it effectively quenches the star formation, transforming blue compact galaxies into compact quiescent galaxies and (2) it also removes and prevents new accretion of cold gas, shutting down in-situ star formation and causing subsequent mergers to be gas-poor or mixed. Gas poor minor mergers then build up an extended stellar envelope. AGN feedback also puffs up the central region through the fast AGN driven winds as well as the slow expulsion of gas while the black hole is quiescent. Without AGN feedback, large amounts of gas accumulate in the central region, triggering star formation and leading to overly massive blue galaxies with dense stellar cores.

Understanding how galaxies maintain the inefficiency of star formation with physically self-consistent models is a central problem for galaxy evolution. Although numerous theoretical models have been proposed in recent decades, the debate still exists. By means of high-resolution two-dimensional hydrodynamical simulations, we study the three feedback effects (the stellar wind heating, SNe feedback, and AGN feedback) in suppressing star formation activities on the evolution of early-type galaxies with different stellar masses. AGN feedback models are updated based on \citet{Yuan2018}. The gas sources comes exclusively from the mass losses of dying low-mass stars for most of our models. We find that SNe feedback can keep star formation at a significantly low level for low mass elliptical galaxies for a cosmological evolution time. For the high mass galaxies, AGN feedback can efficiently offset the radiative cooling and thus regulate the star formation activities. Such a suppression of star formation is extremely efficient in the inner region of the galaxies. AGB heating cannot account for this suppression for low and high mass galaxies. The X-ray temperature $T_{\rm X}$ and luminosity $L_{\rm X}$ of hot plasma can be in agreement with the observed data with the inclusion of effective feedback processes. These results thus suggest that we can use $T_{\rm X}$ and $L_{\rm X}$ to probe the role of different feedback processes. The inclusion of additional gas sources can make the mass scale between SNe and AGN feedback dominating in suppressing star formation decrease to an observationally inferred value of a few $10^{10}~M_{\odot}$.

This is the second paper of our series of works of studying the effects of active galactic nuclei (AGN) feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. In these simulations, the inner boundary is chosen so that the Bondi radius is resolved. Physical processes like star formation, SNe Ia and II are taken into account. Compared to previous works, the main improvements is that we adopt the most updated AGN physics, which is described in detail in the first paper of this series (Yuan et al. 2018, Paper I). These improvements include the discrimination of the two accretion modes of the central AGN and the most updated descriptions of the wind and radiation in the two modes. In Paper I, we consider the case that the specific angular momentum of the gas in the galaxy is very low. In this paper, we consider the case that the specific angular momentum of the gas is high. In the galactic scale, we adopt the gravitational torques raised due to non-axisymmetric structure in the galaxy as the mechanism of the transfer of angular momentum of gas, as proposed in some recent works. Since our simulations are axisymmetric, we make use of a parameterized prescription to mimic this mechanism. Same as Paper I, we investigate the AGN light curve, typical AGN lifetime, growth of the black hole mass, AGN duty-cycle, star formation, and the X-ray surface brightness of the galaxy. Special attention is paid to the effects of specific angular momentum of the galaxy on these properties. We find that some results are qualitatively similar to those shown in Paper I, while some results such as star formation and black hole growth do show a significant difference due to the mass concentration in the galactic disk as a consequence of galactic rotation.

This paper summarizes our recent works of studying AGN feedback in an isolated elliptical galaxy by performing high-resolution hydrodynamical numerical simulations. Bondi radius is resolved and the mass accretion rate of the black hole is calculated. The most updated AGN physics, namely the discrimination of cold and hot accretion modes and the exact descriptions of the AGN radiation and wind for a given accretion rate are adopted and their interaction with the gas in the host galaxy is calculated. Physical processes such as star formation and SNe feedback are taken into account. Consistent with observation, we find the AGN spends most of the time in the low-luminosity regime. AGN feedback overall suppresses the star formation; but depending on location in the galaxy and time, it can also enhance it. The light curve of specific star formation rate is not synchronous with the AGN light curve. We find that wind usually plays a dominant role in controlling the AGN luminosity and star formation, but radiation also cannot be neglected.

A large population of X-ray binaries (XRBs) was recently discovered within the central parsec of the Galaxy by Hailey et al. While the presence of compact objects on this scale due to radial mass segregation is, in itself, unsurprising, the fraction of binaries would naively be expected to be small because of how easily primordial binaries are dissociated in the dynamically hot environment of the nuclear star cluster (NSC). We propose that the formation of XRBs in the central parsec is dominated by the tidal capture of stars by black holes (BHs) and neutron stars (NSs). We model the time-dependent radial density profiles of stars and compact objects in the NSC with a Fokker-Planck approach, using the present-day stellar population and rate of in situ massive star (and thus compact object) formation as observational constraints. Of the ~10,000-40,000 BHs that accumulate in the central parsec over the age of the Galaxy, we predict that ~60 - 200 currently exist as BH-XRBs formed from tidal capture, consistent with the population seen by Hailey et al. A somewhat lower number of tidal capture NS-XRBs is also predicted. We also use our observationally calibrated models for the NSC to predict rates of other exotic dynamical processes, such as the tidal disruption of stars by the central supermassive black hole (~0.0001 per year at z=0).

We explore the effect of momentum-driven winds representing radiation pressure driven outflows from accretion onto supermassive black holes in a set of numerical hydrodynamical simulations. We explore two matched sets of cosmological zoom-in runs of 24 halos with masses ~$10^{12.0}-10^{13.4}$ M_sun run with two different feedback models. Our `NoAGN' model includes stellar feedback via UV heating, stellar winds and supernovae, photoelectric heating and cosmic X-ray background heating from a meta-galactic background. Our fiducial `MrAGN' model is identical except that it also includes a model for black hole seeding and accretion, as well as heating and momentum injection associated with the radiation from black hole accretion. Our MrAGN model launches galactic outflows which result in both `ejective' feedback - the outflows themselves which drive gas out of galaxies - and `preventative' feedback, which suppresses the inflow of new and recycling gas. As much as 80 % of outflowing galactic gas can be expelled, and accretion can be suppressed by as much as a factor of 30 in the MrAGN runs when compared with the NoAGN runs. The histories of NoAGN galaxies are recycling-dominated, with ~70% of material that leaves the galaxy eventually returning, and the majority of outflowing gas re-accretes on 1 Gyr timescales without AGN feedback. Outflowing gas in the MrAGN runs has higher characteristic velocity (500 - 1,000 km/s versus 100-300 km/s for outflowing NoAGN gas) and travels as far as a few Mpcs. Only ~10% of ejected material is re-accreted in the MrAGN galaxies.

A recent determination of the relationships between the X-ray luminosity of the ISM (Lx) and the stellar and total mass, for a sample of nearby early-type galaxies (ETGs), is used to investigate the origin of the hot gas, via a comparison with the results of hydrodynamical simulations of the ISM evolution for a large set of isolated ETGs. After the epoch of major galaxy formation (after z~2), the ISM is replenished by stellar mass losses and SN ejecta, at the rate predicted by stellar evolution, and is depleted by star formation; it is heated by the thermalization of stellar motions, SNe explosions and the mechanical (from winds) and radiative AGN feedback. The models agree well with the observed relations, even for the largely different Lx values at the same mass, thanks to the sensitivity of the gas flow to many galaxy properties; this holds for models including AGN feedback, and those without. Therefore, the mass input from the stellar population is able to account for a major part of the observed Lx; and AGN feedback, while very important to maintain massive ETGs in a time-averaged quasi-steady state, keeping low star formation and the black hole mass, does not dramatically alter the gas content originating in stellar recycled material. These conclusions are based on theoretical predictions for the stellar population contributions in mass and energy, and on a self-consistent modeling of AGN feedback.

Using few-body simulations, we investigate the evolution of supermassive black holes (SMBHs) in galaxies ($M_{\star}=10^{10}-10^{12}{\rm M}_{\odot}$ at $z=0$) at $0<z<4$. Following galaxy merger trees from the Millennium simulation, we model BH mergers with two extreme binary decay scenarios for the `hard binary' stage: a full or an empty loss cone. These two models should bracket the true evolution, and allow us to separately explore the role of dynamical friction and that of multi-body BH interactions on BH mergers. Using the computed merger rates, we infer the stochastic gravitational wave background (GWB). Our dynamical approach is a first attempt to study the dynamical evolution of multiple SMBHs in the host galaxies undergoing mergers with various mass ratios ($10^{-4} < q_{\star} < 1$). Our main result demonstrates that SMBH binaries are able to merge in both scenarios. In the empty loss cone case, we find that BHs merge via multi-body interactions, avoiding the `final parsec' problem, and entering the PTA band with substantial orbital eccentricity. Our full loss cone treatment, albeit more approximate, suggests that the eccentricity becomes even higher when GWs become dominant, leading to rapid coalescences (binary lifetime $\lesssim1 {\rm ~Gyr}$). Despite the lower merger rates in the empty loss cone case, due to their higher mass ratios and lower redshifts, the GWB in the full/empty loss cone models are comparable ($0.70\times10^{-15}$ and $0.53\times10^{-15}$ at a frequency of $1~{\rm yr}^{-1}$, respectively). Finally, we compute the effects of high eccentricities on the GWB spectrum.

Galaxies occupy different regions of the [OIII]$\lambda5007$/H$\beta$-versus-[NII]$\lambda6584$/H$\alpha$ emission-line ratio diagram in the distant and local Universe. We investigate the origin of this intriguing result by modelling self-consistently, for the first time, nebular emission from young stars, accreting black holes (BHs) and older, post-asymptotic-giant-branch (post-AGB) stellar populations in galaxy formation simulations in a full cosmological context. In post-processing, we couple new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies to explore which galaxy physical properties drive the cosmic evolution of the optical-line ratios [OIII]$\lambda5007$/H$\beta$, [NII]$\lambda6584$/H$\alpha$, [SII]$\lambda\lambda6717,6731$/H$\alpha$ and [OI]$\lambda6300$/H$\alpha$. The line ratios of simulated galaxies agree well with observations of both star-forming and active local SDSS galaxies. Towards higher redshifts, at fixed galaxy stellar mass, the average [OIII]/H$\beta$ increases and [NII]/H$\alpha$, [SII]/H$\alpha$ and [OI]/H$\alpha$ decrease -- widely consistent with observations. At fixed stellar mass, we identify star formation history, which controls nebular emission from young stars via the ionization parameter, as the primary driver of the cosmic evolution of [OIII]/H$\beta$ and [NII]/H$\alpha$. For [SII]/H$\alpha$ and [OI]/H$\alpha$, this applies only to redshifts above $z=1.5$, the evolution at lower redshift being driven in roughly equal parts by nebular emission from AGN and post-AGB stars. Instead, changes in the hardness of ionizing radiation, ionized-gas density, the prevalence of BH accretion relative to star formation and the dust-to-metal mass ratio (whose impact on the gas-phase N/O ratio we model at fixed O/H) play at most a minor role in the cosmic evolution of simulated galaxy line ratios.

We ask how the inclusion of various physical heating processes due to the metal content of gas affect the evolution of massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, turbulent diffusion and cooling rate based on individual element abundances. The feedback models include (1) AGN feedback via high velocity BAL winds and Compton/photoionization heating, (2) explicit stellar feedback from multiple processes including powerful winds from supernova events, stellar winds from young massive stars and AGB stars as well as radiative heating within Stromgren spheres around massive stars, and (3) additional heating effects due to the presence of metals including grain photoelectric heating, metallicity dependent X-ray heating by nearby accreting black holes and from the cosmic X-ray background, which are the major improvement in our feedback model. With a suite of zoom-in simulations of 30 halos with $M_{vir} \sim 10^{12-13.4}$, we show that energy and momentum budget from all feedback effects generate realistic galaxy properties. We explore the detailed role of each feedback model with three additional sets of simulations with varying input physics. We show that the metal induced heating mechanisms reduce the fraction of accreted stellar material by mainly suppressing the growth of diffuse small stellar systems at high redshift but overall have a relatively minor effect on the final stellar and gas properties of massive galaxies. The inclusion of AGN feedback significantly improves the ability of our cosmological simulations to yield realistic gas and stellar properties of massive galaxies with reasonable fraction of the final stellar mass accreted from other galaxies.

Galactic outflows are ubiquitously observed in star-forming disk galaxies and are critical for galaxy formation. Supernovae (SNe) play the key role in driving the outflows, but there is no consensus as to how much energy, mass and metal they can launch out of the disk. We perform 3D, high-resolution hydrodynamic simulations to study SNe-driven outflows from stratified media. Assuming SN rate scales with gas surface density $\Sigma_{\rm{gas}}$ as in the Kennicutt-Schmidt (KS) relation, we find the mass loading factor, $\eta_m$, defined as the mass outflow flux divided by the star formation surface density, decreases with increasing $\Sigma_{\rm{gas}}$ as $\eta_m \propto \Sigma^{-0.61}_{\rm{gas}}$. Approximately $\Sigma_{\rm{gas}}$ $\lesssim$ 50 $M_\odot/\rm{pc}^2$ marks when $\eta_m \gtrsim$1. About 10-50\% of the energy and 40-80\% of the metals produced by SNe end up in the outflows. The tenuous hot phase ($T>3\times 10^5$ K), which fills 60-80\% of the volume at mid-plane, carries the majority of the energy and metals in outflows. We discuss how various physical processes, including vertical distribution of SNe, photoelectric heating, external gravitational field and SN rate, affect the loading efficiencies. The relative scale height of gas and SNe is a very important factor in determining the loading efficiencies.

We present three-dimensional hydrodynamical simulations showing the effect of kinetic and radiative AGN feedback on a model galaxy representing a massive quiescent low-redshift early-type galaxy of $M_* = 8.41\times 10^{10} M_\odot$, harbouring a $M_\mathrm{BH} = 4\times 10^8 M_\odot $ black hole surrounded by a cooling gaseous halo. We show that, for a total baryon fraction of $\sim 20\%$ of the cosmological value, feedback from the AGN can keep the galaxy quiescent for about 4.35 Gyr and with properties consistent with black hole mass and X-ray luminosity scaling relations. However, this can only be achieved if the AGN feedback model includes both kinetic and radiative feedback modes. The simulation with only kinetic feedback fails to keep the model galaxy fully quiescent, while one with only radiative feedback leads to excessive black-hole growth. For higher baryon fractions (e.g. 50\% of the cosmological value), the X-ray luminosities exceed observed values by at least one order of magnitude, and rapid cooling results in a star-forming galaxy. The AGN plays a major role in keeping the circumgalactic gas at observed metallicities of $Z/Z_\odot \gtrsim 0.3 $ within the central $\sim 30$ kpc by venting nuclear gas enriched with metals from residual star formation activity. As indicated by previous cosmological simulations, our results are consistent with a model for which the black hole mass and the total baryon fraction are set at higher redshifts $z > 1$ and the AGN alone can keep the model galaxy on observed scaling relations. Models without AGN feedback violate both the quiescence criterion as well as CGM metallicity constraints.

An intriguing alternative to cold dark matter (CDM) is that the dark matter is a light ( $m \sim 10^{-22}$ eV) boson having a de Broglie wavelength $\lambda \sim 1$ kpc, often called fuzzy dark matter (FDM). We describe the arguments from particle physics that motivate FDM, review previous work on its astrophysical signatures, and analyze several unexplored aspects of its behavior. In particular, (i) FDM halos smaller than about $10^7 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$ do not form. (ii) FDM halos are comprised of a core that is a stationary, minimum-energy configuration called a "soliton", surrounded by an envelope that resembles a CDM halo. (iii) The transition between soliton and envelope is determined by a relaxation process analogous to two-body relaxation in gravitating systems, which proceeds as if the halo were composed of particles with mass $\sim \rho\lambda^3$ where $\rho$ is the halo density. (iv) Relaxation may have substantial effects on the stellar disk and bulge in the inner parts of disk galaxies. (v) Relaxation can produce FDM disks but an FDM disk in the solar neighborhood must have a half-thickness of at least $300 (m/10^{-22} {\rm eV})^{-2/3}$ pc. (vi) Solitonic FDM sub-halos evaporate by tunneling through the tidal radius and this limits the minimum sub-halo mass inside 30 kpc of the Milky Way to roughly $10^8 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$. (vii) If the dark matter in the Fornax dwarf galaxy is composed of CDM, most of the globular clusters observed in that galaxy should have long ago spiraled to its center, and this problem is resolved if the dark matter is FDM.

The OVI $\lambda\lambda$1032, 1038\AA\ doublet emission traces collisionally ionized gas with $T\approx 10^{5.5}$ K, where the cooling curve peaks for metal-enriched plasma. This warm-hot phase is usually not well-resolved in numerical simulations of the multiphase interstellar medium (ISM), but can be responsible for a significant fraction of the emitted energy. Comparing simulated OVI emission to observations is therefore a valuable test of whether simulations predict reasonable cooling rates from this phase. We calculate OVI $\lambda$1032\AA\ emission, assuming collisional ionization equilibrium, for our small-box simulations of the stratified ISM regulated by supernovae. We find that the agreement is very good for our solar neighborhood model, both in terms of emission flux and mean OVI density seen in absorption. We explore runs with higher surface densities and find that, in our simulations, the OVI emission from the disk scales roughly linearly with the star formation rate. Observations of OVI emission are rare for external galaxies, but our results do not show obvious inconsistency with the existing data. Assuming the solar metallicity, OVI emission from the galaxy disk in our simulations accounts for roughly 0.5\% of supernovae heating.

Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, energy conserving Sedov-Taylor, and momentum conserving snowplow$-$ with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, $f_{\mathrm{hot}}$, as $(1 - f_{\mathrm{hot}})^{-4/5}$. We also include winds from young massive stars and AGB stars, Str\"omgren sphere gas heating by massive stars, and a gas cooling limiting mechanism driven by radiative recombination of dense HII regions. We present initial tests for isolated Milky-Way like systems simulated with the GADGET based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a two order of magnitude variation in particle mass in the range (1.3$-$130)$\times 10^4$ solar masses.

The interplay between supermassive black holes (SMBHs) and their environments is believed to command an essential role in galaxy evolution. The majority of these SMBHs are in the radiative inefficient accretion phase where this interplay remains elusive, but suggestively important, due to few observational constraints. To remedy this, we directly fit 2-D hydrodynamic simulations to Chandra observations of Sgr A* with Markov Chain Monte Carlo sampling, self-consistently modelling the 2-D inflow-outflow solution for the first time. We find the temperature and density at flow onset are consistent with the origin of the gas in the stellar winds of massive stars in the vicinity of Sgr A*. We place the first observational constraints on the angular momentum of the gas and estimate the centrifugal radius, r$_c$ $\approx$ 0.056 r$_b$ $\approx8\times10^{-3}$ pc, where r$_b$ is the Bondi radius. Less than 1\% of the inflowing gas accretes onto the SMBH, the remainder being ejected in a polar outflow. For the first time, we decouple the quiescent point-like emission from the spatially extended flow. We find this point-like emission, accounting for $\sim4$\% of the quiescent flux, is spectrally too steep to be explained by unresolved flares, nor Bremsstrahlung, but is likely a combination of a relatively steep synchrotron powerlaw and the high energy tail of inverse-Compton emission. With this self-consistent model of the accretion flow structure, we make predictions for the flow dynamics and discuss how future X-ray spectroscopic observations can further our understanding of the Sgr A* accretion flow.

Numerical simulations have become a major tool for understanding galaxy formation and evolution. Over the decades the field has made significant progress. It is now possible to simulate the formation of individual galaxies and galaxy populations from well defined initial conditions with realistic abundances and global properties. An essential component of the calculation is to correctly estimate the inflow to and outflow from forming galaxies since observations indicating low formation efficiency and strong circum-glactic presence of gas are persuasive. Energetic 'feedback' from massive stars and accreting super-massive black holes - generally unresolved in cosmological simulations - plays a major role for driving galactic outflows, which have been shown to regulate many aspects of galaxy evolution. A surprisingly large variety of plausible sub-resolution models succeeds in this exercise. They capture the essential characteristics of the problem, i.e. outflows regulating galactic gas flows, but their predictive power is limited. In this review we focus on one major challenge for galaxy formation theory: to understand the underlying physical processes that regulate the structure of the interstellar medium, star formation and the driving of galactic outflows. This requires accurate physical models and numerical simulations, which can precisely describe the multi-phase structure of the interstellar medium on the currently unresolved few hundred parsecs scales of large scale cosmological simulations. Such models ultimately require the full accounting for the dominant cooling and heating processes, the radiation and winds from massive stars and accreting black holes, an accurate treatment of supernova explosions as well as the non-thermal components of the interstellar medium like magnetic fields and cosmic rays.

We present 2D hydrodynamical simulations for the evolution of early-type galaxies containing central massive black holes (MBHs), starting at age 2 Gyr. The code contains accurate and physically consistent radiative and mechanical AGN wind feedback, with parsec-scale central resolution. Mass input comes from stellar evolution; energy input includes Type Ia and II supernova and stellar heating; star-formation is included. Realistic, axisymmetric dynamical models for the galaxies are built solving the Jeans' equations. The lowest mass models (Mstar = 8 10^{10}Msun) develop global outflows sustained by SNIa's heating, ending with a significantly lower amount of hot gas and new stars. In more massive models, nuclear outbursts last to the present epoch, with large and frequent fluctuations in nuclear emission and from the gas (Lx). Each burst last ~ 10^{7.5} yr, during which (for r < 2-3 kpc) cold, inflowing, and hot, outflowing gas phases coexist. The Lx-T relation for the gas matches that of local galaxies. AGN activity causes positive feedback for star formation. Roughly half of the total mass loss is recycled into new stars (DeltaMstar), just ~ 3% of it is accreted on the MBH, the remainder being ejected from the galaxy. The ratio between the mass of gas expelled to that in to new stars, the load factor, is ~0.6. Rounder galaxies shapes lead to larger final MBH masses, DeltaMstar, and Lx. Almost all the time is spent at very low nuclear luminosities, yet one quarter of the total energy is emitted at an Eddington ratio > 0.1. The duty-cycle of AGN activity approximates 4% (Abridged).

Nuclear star clusters (NSCs) and supermassive black holes (SMBHs) both inhabit galactic nuclei, coexisting in a range of bulge masses, but excluding each other in the largest or smallest galaxies. We propose that the transformation of NSCs into SMBHs occurs via runaway tidal captures, once NSCs exceed a certain critical central density and velocity dispersion. The bottleneck in this process, as with all collisional runaways, is growing the first e-fold in black hole mass. The growth of a stellar mass black hole past this bottleneck occurs as tidally captured stars are consumed in repeated episodes of mass transfer at pericenter. Tidal captures may turn off as a growth channel once the black hole reaches a mass ~100-1000 solar masses, but tidal disruption events will continue and appear capable of growing the seed SMBH to larger sizes. The runaway slows (becomes sub-exponential) once the seed SMBH consumes the core of its host NSC. While the bulk of the cosmic mass density in SMBHs is ultimately produced (via the Soltan-Paczynski argument) by episodic gaseous accretion in very massive galaxies, the smallest SMBHs have probably grown from strong tidal encounters with NSC stars. SMBH seeds that grow for a time $t$ entirely through this channel will follow simple power law relations with the velocity dispersion, $\sigma$, of their host galaxy. In the simplest regime it is $M_\bullet \sim \sigma^{3/2}\sqrt{M_\star t / G} \sim 10^{6}M_\odot (\sigma / 50~{\rm km~s}^{-1})^{3/2}(t/10^{10}~{\rm yr})^{1/2}$, but the exponents and prefactor can differ slightly depending on the details of loss cone refilling. Current tidal disruption event rates predicted from this mechanism are consistent with observations.

We study very-high rate spherically symmetric accretion flows onto a massive black hole (BH; 10^2 < M_BH < 10^6 Msun) embedded in a dense gas cloud with a low abundance of metals, performing one-dimensional hydrodynamical simulations which include multi-frequency radiation transfer and non-equilibrium primordial chemistry. We find that rapid gas supply from the Bondi radius at a hyper-Eddington rate can occur without being impeded by radiation feedback when (n/10^5 cm^-3) > (M_BH/10^4Msun)^{-1}(T/10^4 K)^{3/2}, where n and T are the density and temperature of ambient gas outside of the Bondi radius. The resulting accretion rate in this regime is steady, and larger than 3000 times the Eddington rate. At lower Bondi rates, the accretion is episodic due to radiative feedback and the average rate is limited below the Eddington rate. For the hyper-Eddington case, the steady solution consists of two parts: a radiation-dominated central core, where photon trapping due to electron scattering is important, and an accreting envelope which follows a Bondi profile with T~8000 K. When the emergent luminosity is limited below the Eddington luminosity because of photon trapping, radiation from the central region does not affect the gas dynamics at larger scales. We apply our result to the rapid formation of massive BHs in protogalaxies with a virial temperature of T_vir> 10^4 K. Once a seed BH forms at the center of the galaxy, it can grow up to a maximum ~10^5 (T_vir/10^4 K) Msun via gas accretion independent of the initial BH mass. Finally, we discuss possible observational signatures of rapidly accreting BHs with/without allowance for dust. We suggest that these systems could explain Lya emitters without X-rays and luminous infrared sources with hot dust emission, respectively.

We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is computed using a chemical network that follows the abundances of H+, H, H2, CO, C+, and free electrons and includes (self-)shielding of the gas and dust. We find that CRs perceptibly thicken the disk with the heights of 90% (70%) enclosed mass reaching ~1.5 kpc (~0.2 kpc). The simulations indicate that CRs alone can launch and sustain strong outflows of atomic and ionized gas with mass loading factors of order unity, even in solar neighborhood conditions and with a CR energy injection per SN of 10^50 erg, 10% of the fiducial thermal energy of an SN. The CR-driven outflows have moderate launching velocities close to the midplane (~100 km/s) and are denser (\rho~1e-24 - 1e-26 g/cm^3), smoother, and colder than the (thermal) SN-driven winds. The simulations support the importance of CRs for setting the vertical structure of the disk as well as the driving of winds.

AGN feedback from supermassive black holes (SMBHs) at the center of early type galaxies is commonly invoked as the explanation for the quenching of star formation in these systems. The situation is complicated by the significant amount of mass injected in the galaxy by the evolving stellar population over cosmological times. In absence of feedback, this mass would lead to unobserved galactic cooling flows, and to SMBHs two orders of magnitude more massive than observed. By using high-resolution 2D hydrodynamical simulations with radiative transport and star formation in state-of-the-art galaxy models, we show how the intermittent AGN feedback is highly structured on spatial and temporal scales, and how its effects are not only negative (shutting down the recurrent cooling episodes of the ISM), but also positive, inducing star formation in the inner regions of the host galaxy.

Supernovae (SN), the most energetic stellar feedback mechanism, are crucial for regulating the interstellar medium (ISM) and launching galactic winds. We explore how supernova remnants (SNRs) create a multiphase medium by performing 3D hydrodynamical simulations at various SN rates, $S$, and ISM average densities, $\bar{n}$. The evolution of a SNR in a self-consistently generated three-phase ISM is qualitatively different from that in a uniform or a two-phase warm/cold medium. By travelling faster and further in the low-density hot phase, the domain of a SNR increases by $>10^{2.5}$. Varying $\bar{n}$ and $S$, we find that a steady state can only be achieved when the hot gas volume fraction $f_{\rm{V,hot}}\lesssim 0.6 \pm 0.1 $. Above that level, overlapping SNRs render connecting topology of the hot gas, and the ISM is subjected to thermal runaway. Photoelectric heating (PEH) has a surprisingly strong impact on $f_{\rm{V,hot}}$. For $\bar{n}\gtrsim 3 \cm-3 $, a reasonable PEH rate is able to suppress the thermal runaway. Overall, we determine the critical SN rate for the onset of thermal runaway to be $S_{\rm{crit}} = 200 (\bar{n}/1\cm-3)^k (E_{\rm{SN}}/10^{51}\erg)^{-1} \kpc^{-3} \myr-1$, where $k = (1.2,2.7)$ for $\bar{n} \leq 1$ and $> 1\cm-3 $, respectively. We present a fitting formula of the ISM pressure $P(\bar{n}$, $S$), which can be used as an effective equation of state in cosmological simulations. Despite the 5 orders of magnitude span of $(\bar{n},S)$, the average Mach number varies little: $\mathcal{M} \approx \ 0.5\pm 0.2, \ 1.2\pm 0.3,\ 2.3\pm 0.9$ for the hot, warm and cold phases, respectively.

We present a potential-density pair designed to model nearly isothermal star clusters (and similar self-gravitating systems) with a central core and an outer turnover radius, beyond which density falls off as $r^{-4}$. In the intermediate zone, the profile is similar to that of an isothermal sphere (density $\rho \propto r^{-2}$), somewhat less steep than the King 62 profile, and with the advantage that many dynamical quantities can be written in a simple closed form. We derive new analytic expressions for the cluster binding energy and velocity dispersion, and apply these to create toy models for cluster core collapse and evaporation. We fit our projected surface brightness profiles to observed globular and open clusters, and find that the quality of the fit is generally at least as good as that for the surface brightness profiles of King 62. This model can be used for convenient computation of the dynamics and evolution of globular and nuclear star clusters.

We describe a new method for adding a prescribed amount of kinetic energy to simulated gas modeled on a cartesian grid by directly altering grid cells' mass and velocity in a distributed fashion. The method is explored in the context of supernova feedback in high-resolution ($\sim 10$ pc) hydrodynamic simulations of galaxy formation. Resolution-dependence is a primary consideration in our application of the method and simulations of isolated explosions (performed at different resolutions) motivate a resolution-dependent scaling for the injected fraction of kinetic energy that we apply in cosmological simulations of a $10^9$ Msun dwarf halo. We find that in high density media ($\gtrsim$ 50 cm$^{-3}$) with coarse resolution ($\gtrsim 4$ pc per cell), results are sensitive to the initial kinetic energy fraction due to early and rapid cooling. In our galaxy simulations, the deposition of small amounts of supernova energy in kinetic form (as little as 1%) has a dramatic impact on the evolution of the system, resulting in an order of magnitude suppression of stellar mass. The overall behavior of the galaxy in the two highest resolution simulations we perform appears to converge. We discuss the resulting distribution of stellar metallicities, an observable sensitive to galactic wind properties, and find that while the new method demonstrates increased agreement with observed systems, significant discrepancies remain, likely due to simplistic assumptions that neglect contributions from Type Ia supernovae and stellar winds.

Halo Abundance Matching has been used to construct a one-parameter mapping between galaxies and dark matter haloes by assuming that halo mass and galaxy luminosity (or stellar mass) are monotonically related. While this approach has been reasonably successful, it is known that galaxies must be described by at least two parameters, as can be seen from the two-parameter Fundamental Plane on which massive early-type galaxies lie. In this paper, we derive a connection between initial dark matter density perturbations in the early universe and present-day virialized dark matter haloes by assuming simple spherical collapse combined with conservation of mass and energy. We find that $z = 0$ halo concentration, or alternatively the inner slope of the halo density profile $\alpha$, is monotonically and positively correlated with the collapse redshift of the halo. This is qualitatively similar to the findings of some previous works based on numerical simulations, with which we compare our results. We then describe how the halo mass and concentration (or inner slope $\alpha$) can be used as two halo parameters in combination with two parameters of early-type galaxies to create an improved abundance matching scheme.

We employ cosmological hydrodynamical simulations to investigate the effects of AGN feedback on the formation of massive galaxies with present-day stellar masses of $M_{stel} = 8.8 \times 10^{10} - 6.0 \times 10^{11} M_{sun}$. Using smoothed particle hydrodynamics simulations with a pressure-entropy formulation that allows an improved treatment of contact discontinuities and fluid mixing, we run three sets of simulations of 20 halos with different AGN feedback models: (1) no feedback, (2) thermal feedback, and (3) mechanical and radiation feedback. We assume that seed black holes are present at early cosmic epochs at the centre of emerging dark matter halos and trace their mass growth via gas accretion and mergers with other black holes. Both feedback models successfully recover the observed M_BH - sigma relation and black hole-to-stellar mass ratio for simulated central early-type galaxies. The baryonic conversion efficiencies are reduced by a factor of two compared to models without any AGN feedback at all halo masses. However, massive galaxies simulated with thermal AGN feedback show a factor of ~10-100 higher X-ray luminosities than observed. The mechanical/radiation feedback model reproduces the observed correlation between X-ray luminosities and velocity dispersion, e.g. for galaxies with sigma = 200 km/s, the X-ray luminosity is reduced from $10^{42}$ erg/s to $10^{40}$ erg/s. It also efficiently suppresses late time star formation, reducing the specific star formation rate from $10^{-10.5}$ $yr^{-1}$ to $10^{-14}$ $yr^{-1}$ on average and resulting in quiescent galaxies since z=2, whereas the thermal feedback model shows higher late time in-situ star formation rates than observed.

We use hydrodynamical simulations in a $(256\;{\rm pc})^3$ periodic box to model the impact of supernova (SN) explosions on the multi-phase interstellar medium (ISM) for initial densities $n = 0.5-30$ cm$^{-3}$ and SN rates $1-720$ Myr$^{-1}$. We include radiative cooling, diffuse heating, and the formation of molecular gas using a chemical network. The SNe explode either at random positions, at density peaks, or both. We further present a model combining thermal energy for resolved and momentum input for unresolved SNe. Random driving at high SN rates results in hot gas ($T\gtrsim 10^6$ K) filling $> 90$% of the volume. This gas reaches high pressures ($10^4 < P/k_\mathrm{B} < 10^7$ K cm$^{-3}$) due to the combination of SN explosions in the hot, low density medium and confinement in the periodic box. These pressures move the gas from a two-phase equilibrium to the single-phase, cold branch of the cooling curve. The molecular hydrogen dominates the mass ($>50$%), residing in small, dense clumps. Such a model might resemble the dense ISM in high-redshift galaxies. Peak driving results in huge radiative losses, producing a filamentary ISM with virtually no hot gas, and a small molecular hydrogen mass fraction ($\ll 1$%). Varying the ratio of peak to random SNe yields ISM properties in between the two extremes, with a sharp transition for equal contributions. The velocity dispersion in HI remains $\lesssim 10$ km s$^{-1}$ in all cases. For peak driving the velocity dispersion in H$_\alpha$ can be as high as $70$ km s$^{-1}$ due to the contribution from young, embedded SN remnants.

We investigate the evolution of stellar population gradients from $z=2$ to $z=0$ in massive galaxies at large radii ($r > 2R_{\mathrm{eff}}$) using ten cosmological zoom simulations of halos with $6 \times 10^{12} M_{\odot} < M_{\mathrm{halo}} < 2 \times 10^{13}M_{\odot}$. The simulations follow metal cooling and enrichment from SNII, SNIa and AGB winds. We explore the differential impact of an empirical model for galactic winds that reproduces the mass-metallicity relation and its evolution with redshift. At larger radii the galaxies, for both models, become more dominated by stars accreted from satellite galaxies in major and minor mergers. In the wind model, fewer stars are accreted, but they are significantly more metal poor resulting in steep global metallicity ($\langle \nabla Z_{\mathrm{stars}} \rangle= -0.35$ dex/dex) and color (e.g. $\langle \nabla g-r \rangle = -0.13$ dex/dex) gradients in agreement with observations. In contrast, colour and metallicity gradients of the models without winds are inconsistent with observations. Age gradients are in general mildly positive at $z=0$ ($\langle \nabla Age_{\mathrm{stars}} \rangle= 0.04$ dex/dex) with significant differences between the models at higher redshift. We demonstrate that for the wind model, stellar accretion is steepening existing in-situ metallicity gradients by about 0.2 dex by the present day and helps to match observed gradients of massive early-type galaxies at large radii. Colour and metallicity gradients are significantly steeper for systems which have accreted stars in minor mergers, while galaxies with major mergers have relatively flat gradients, confirming previous results. This study highlights the importance of stellar accretion for stellar population properties of massive galaxies at large radii, which can provide important constraints for formation models.

Accretion is thought to primarily contribute to the mass accumulation history of supermassive black holes throughout cosmic time. While this may be true at high redshifts, at lower redshifts and for the most massive black holes mergers themselves might add significantly to the mass budget. We evolve SMBHs from $4 > z > 0$ using merger trees derived from hydrodynamical cosmological simulations of a cluster and void region, scaled to the observed value of the stellar mass fraction to account for overcooling. Mass gains from gas accretion proportional to bulge growth and BH-BH mergers are tracked, as are black holes that remain "orbiting" due to insufficient dynamical friction in a merger remnant, as well as those that are ejected due to gravitational recoil. We find that gas accretion remains the dominant source of mass accumulation in almost all SMBHs; mergers contribute $2.5\pm0.1\%$ for all SMBHs in the cluster and $1.0\pm0.1\%$ in the void since $z = 4$. However, mergers are significant for massive SMBHs. The fraction of mass accumulated from mergers for central BHs generally increases for larger values of the host bulge mass: in the void, the fraction is $2\%$ at $M_{*, bul} = 10^{10} M_{\odot}$, increasing to $4\%$ at $M_{*, bul} \gtrsim 10^{11} M_{\odot}$, and in the cluster it is $4\%$ at $M_{*, bul} = 10^{10} M_{\odot}$ and $23\%$ at $10^{12} M_{\odot}$. We find that $40\%$ of SMBHs and $\approx 8\%$ of the total SMBH mass is found orbiting in the cluster region at $z = 0$. The existence of orbiting and ejected SMBHs requires modification of the Soltan argument. We estimate this correction to the integrated accreted mass density of SMBHs to be in the range $6-21\%$, with a mean value of $11\pm3\%$. We also calculate the total energy output and strain from gravitational waves emitted by merging SMBHs, and obtain a signal potentially detectable by pulsar timing arrays.

Based on two-dimensional high resolution hydrodynamic numerical simulation, we study the mechanical and radiative feedback effects from the central AGN on the cosmological evolution of an isolated elliptical galaxy. Physical processes such as star formation and supernovae are considered. The inner boundary of the simulation domain is carefully chosen so that the fiducial Bondi radius is resolved and the accretion rate of the black hole is determined self-consistently. In analogy to previous works, we assume that the specific angular momentum of the galaxy is low. It is well-known that when the accretion rates are high and low, the central AGNs will be in cold and hot accretion modes, which correspond to the radiative and kinetic feedback modes, respectively. The emitted spectrum from the hot accretion flows is harder than that from the cold accretion flows, which could result in a higher Compton temperature accompanied by a more efficient radiative heating, according to previous theoretical works. Such a difference of the Compton temperature between the two feedback modes, the focus of this study, has been neglected in previous works. Significant differences in the kinetic feedback mode are found as a result of the stronger Compton heating and accretion becomes more chaotic. More importantly, if we constrain models to correctly predict black hole growth and AGN duty cycle after cosmological evolution, we find that the favored model parameters are constrained: mechanical feedback efficiency diminishes with decreasing luminosity (the maximum efficiency being $\simeq 10^{-3.5}$) and X-ray Compton temperature increases with decreasing luminosity, although models with fixed mechanical efficiency and Compton temperature can be found that are satisfactory as well. We conclude that radiative feedback in the kinetic mode is much more important than previously thought.

We study the effects of a detailed dust treatment on the properties and evolution of early-type galaxies containing central black holes, as determined by AGN feedback. We find that during cooling flow episodes, radiation pressure on the dust in and interior to infalling shells of cold gas can greatly impact the amount of gas able to be accreted and therefore the frequency of AGN bursts. However, the overall hydrodynamic evolution of all models, including mass budget, is relatively robust to the assumptions on dust. We find that IR re-emission from hot dust can dominate the bolometric luminosity of the galaxy during the early stages of an AGN burst, reaching values in excess of $10^{46}$ erg/s. The AGN-emitted UV is largely absorbed, but the optical depth in the IR does not exceed unity, so the radiation momentum input never exceeds $L_{\rm BH}/c$. We constrain the viability of our models by comparing the AGN duty cycle, broadband luminosities, dust mass, black hole mass, and other model predictions to current observations. These constraints force us to models wherein the dust to metals ratios are $\simeq 1%$ of the Galactic value, and only models with a dynamic dust to gas ratio are able to produce both quiescent galaxies consistent with observations and high obscured fractions during AGN "on" phases. During AGN outbursts, we predict that a large fraction of the FIR luminosity can be attributed to warm dust emission ($\simeq100$ K) from dense dusty gas within $\leq 1$ kpc reradiating the AGN UV emission.

We study the effect of AGN mechanical and radiation feedback on the formation of bulge dominated galaxies via mergers of disc galaxies. The merging galaxies have mass-ratios of 1:1 to 6:1 and include pre-existing hot gaseous halos to properly account for the global impact of AGN feedback. Using smoothed particle hydrodynamics simulation code (GADGET-3) we compare three models with different AGN feedback models: (1) no black hole and no AGN feedback; (2) thermal AGN feedback; and (3) mechanical and radiative AGN feedback. The last model is motivated by observations of broad line quasars which show winds with initial velocities of $v_w \ge$ 10,000 km/s and also heating associated with the central AGN X-ray radiation. The primary changes in gas properties due to mechanical AGN feedback are lower thermal X-ray luminosity from the final galaxy - in better agreement with observations - and galactic outflows with higher velocity $\sim 1000$ km/s similar to recent direct observations of nearby merger remnants. The kinetic energy of the outflowing gas is a factor of $\sim$ 20 higher than in the thermal feedback case. All merger remnants with momentum-based AGN feedback with $v_w \sim 10,000$ km/s and $\epsilon_w=2 \times 10^{-3}$, independent of their progenitor mass-ratios, reproduce the observed relations between stellar velocity dispersion and black hole mass ($M_{\rm bh} - \sigma$) as well as X-ray luminosity ($L_X - \sigma$) with $10^{37.5}$ erg/s $\lesssim L_X (0.3-8~{\rm keV}) \lesssim 10^{39.2}$ erg/s for velocity dispersions in the range of 120 km/s $\lesssim \sigma \lesssim$ 190 km/s. In addition, the mechanical feedback produces a much greater AGN variability. We also show that gas is more rapidly and impulsively stripped from the galactic centres driving a moderate increase in galaxy size and decrease in central density with the mechanical AGN feedback model.

We revisit the hypothesis that dense galactic nuclei are formed from inspiraling globular clusters. Recent advances in understanding of the continuous formation of globular clusters over cosmic time and the concurrent evolution of the galaxy stellar distribution allow us to construct a simple model that matches the observed spatial and mass distributions of clusters in the Galaxy and the giant elliptical galaxy M87. In order to compare with observations, we model the effects of dynamical friction and dynamical evolution, including stellar mass loss, tidal stripping of stars, and tidal disruption of clusters by the growing galactic nucleus. We find that inspiraling globular clusters form a dense central structure, with mass and radius comparable to the typical values in observed nuclear star clusters (NSCs) in late-type and low-mass early-type galaxies. The density contrast associated with the NSC is less pronounced in giant elliptical galaxies. Our results indicate that the NSC mass as a fraction of mass of the galaxy stellar spheroid scales as M_NSC/M_* = 0.0025 M_{*,11}^{-0.5}. Thus disrupted globular clusters could contribute most of the mass of NSCs in galaxies with stellar mass below 10^{11} Msun. The inner part of the accumulated cluster may seed the growth of a central black hole via stellar dynamical core collapse, thereby relieving the problem of how to form luminous quasars at high redshift. The seed black hole may reach about 10^5 Msun within 1 Gyr of the beginning of globular cluster formation.

Using basic physical arguments, we derive by dimensional and physical analysis the characteristic masses and sizes of important objects in the Universe in terms of just a few fundamental constants. This exercise illustrates the unifying power of physics and the profound connections between the small and the large in the Cosmos we inhabit. We focus on the minimum and maximum masses of normal stars, the corresponding quantities for neutron stars, the maximum mass of a rocky planet, the maximum mass of a white dwarf, and the mass of a typical galaxy. To zeroth order, we show that all these masses can be expressed in terms of either the Planck mass or the Chandrasekar mass, in combination with various dimensionless quantities. With these examples we expose the deep interrelationships imposed by Nature between disparate realms of the Universe and the amazing consequences of the unifying character of physical law.

We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.

We perform time-dependent, 2DHD numerical simulations to study the dynamics of a slowly rotating accretion flow from sub-pc to pc scales under the irradiation from the central AGN. Compared to previous work, we improve the calculation of the radiative force due to X-rays. More importantly, in addition to radiative pressure and radiative heating/cooling directly from the central AGN, in the momentum equation we also include the force due to the scattered and reprocessed photons. We find that the accretion flow properties change significantly due to this "re-radiation" effect. The inflow rate at the inner boundary is reduced, while the outflow rate at the outer boundary is enhanced by about one order of magnitude. This effect is more significant when the density at the outer boundary is higher. The properties of outflows such as velocity, momentum and energy fluxes, and the ratio of outflow rate and the accretion rate, are calculated. We find that the efficiency of transferring the radiation power into the kinetic power of outflow is typically $10^{-3}$, far below the value of $\sim 0.05$ which is assumed in some cosmological simulations. The effect of the temperature of the gas at the outer boundary ($T_0$) is investigated. When $T_0$ is high, the emitted luminosity of the accretion flow oscillates. This is because in this case the gas around the Bondi radius can be more easily heated to be above the virial temperature due to its high internal energy. Another question we hope to address is the so-called "sub-Eddington" puzzle. Observationally, the luminosity of almost all AGNs are sub-Eddington, while theoretically the luminosity of an accretion flow can easily be super-Eddington. We find that even when the re-radiation effect is included and outflow does become much stronger, the luminosity, while reduced, can still be super-Eddington.

Accretion onto a supermassive black hole of a rotating inflow is a particularly difficult problem to study because of the wide range of length scales involved. There have been broadly utilized analytic and numerical treatments of the global properties of accretion flows, but detailed numerical simulations are required to address certain critical aspects. We use the ZEUS code to run hydrodynamical simulations of rotating, axisymmetric accretion flows with Bremsstrahlung cooling, considering solutions for which the centrifugal balance radius significantly exceeds the Schwarzschild radius, with and without viscous angular momentum transport. Infalling gas is followed from well beyond the Bondi radius down to the vicinity of the black hole. We produce a continuum of solutions with respect to the single parameter Mdot_Bondi/Mdot_Edd, and there is a sharp transition between two general classes of solutions at an Eddington ratio of Mdot_Bondi/Mdot_Edd ~ few x 10^(-2). Our high inflow solutions are very similar to the standard Shakura & Sunyaev (1973) results. But our low inflow results are to zeroth order the stationary Papaloizou and Pringle (1984) solution, which has no accretion. To next order in the small, assumed viscosity they show circulation, with disk and conical wind outflows almost balancing inflow. These solutions are characterized by hot, vertically extended disks, and net accretion proceeds at an extremely low rate, only of order alpha times the inflow rate. Our simulations have converged with respect to spatial resolution and temporal duration, and they do not depend strongly on our choice of boundary conditions.

We present a model for merger-driven evolution of the mass function for massive galaxies and their central supermassive black holes at late times. We discuss the current observational evidence in favor of merger-driven massive galaxy evolution during this epoch, and demonstrate that the observed evolution of the mass function can be reproduced by evolving an initial mass function under the assumption of negligible star formation. We calculate the stochastic gravitational wave signal from the resulting black-hole binary mergers in the low redshift universe (z <= 1) implied by this model, and find that this population has a signal-to-noise ratio as much as ~5x larger than previous estimates for pulsar timing arrays, with an expectation value for the characteristic strain h_c (f=1 yr^{-1}) = 4.1 x 10^{-15} that may already be in tension with observational constraints, and a {2-sigma, 3-sigma} lower limit within this model of h_c (f=1 yr^{-1}) = {1.1 x 10^{-15}, 6.8 x 10^{-16}}. The strength of this signal is sufficient to make it detectable with high probability under conservative assumptions within the next several years, if the principle assumption of merger-driven galaxy evolution since z = 1 holds true. For cases where a galaxy merger fails to lead to a black hole merger, we estimate the probability for a given number of satellite unmerged black holes to remain within a massive host galaxy, and interpret the result in light of ULX observations. In particular, we find that the brightest cluster galaxies should have 1-2 such sources with luminosities above 10^{39} erg/s, which is consistent with the statistics of observed ULXs.

There is observational evidence for inside-out growth of elliptical galaxies since $z \gtrsim 2-3$, which is not driven by in-situ star formation. Many systems at high redshift have small sizes $\sim 1kpc$ and surface brightness profiles with low Sersic indices n. The most likely descendants have, on average, grown by a factor of two in mass and a factor of four in size, indicating $r \propto M^{\alpha}$ with $\alpha \gtrsim 2$. They also have surface brightness profiles with $n \gtrsim 5$. This evolution can be qualitatively explained on the basis of two assumptions: compact ellipticals predominantly grow by collisionless minor or intermediate 'dry' mergers, and they are embedded in massive dark matter halos. We draw these conclusions from idealized collisionless mergers spheroidal galaxies - with and without dark matter - with mass ratios of 1:1, 1:5, and 1:10. The sizes evolve as $r \propto M^{\alpha}$ with $\alpha < 2$ for mass-ratios of 1:1. For minor mergers of galaxies embedded in dark matter halos, the sizes grow significantly faster and the profile shapes change more rapidly. Mergers with moderate mass-ratios of 1:5 give $\alpha \sim 2.3$ and a final Sersic index of $n = 9.5$ after doubling the stellar mass. This is accompanied by a significant increase of the dark matter fraction within the stellar half-mass radius, driven by the strong size increase probing larger, dark matter dominated regions. Only a few intermediate mass-ratio mergers of galaxies embedded in massive dark matter halos can result in the observed concurrent inside-out growth and the rapid evolution in profile shapes. Apart from negative stellar metallicity gradients such a 'minor' merger scenario also predicts significantly lower dark matter fractions for $z \sim 2$ compact quiescent galaxies and their rare present day analogues (abbreviated).

Recent observations of massive galaxies indicate that they double in mass and quintuple in size between redshift z = 1 and the present, despite undergoing very little star formation, suggesting that galaxy mergers drive the evolution. Since these galaxies will contain supermassive black holes, this suggests a larger black hole merger rate, and therefore a larger gravitational-wave signal, than previously expected. We calculate the merger-driven evolution of the mass function, and find that merger rates are 10 to 30 times higher and gravitational waves are 3 to 5 times stronger than previously estimated, so that the gravitational-wave signal may already be detectable with existing data from pulsar timing arrays. We also provide an explanation for the disagreement with past estimates that were based on dark matter halo simulations.

We study the growth of black holes (BHs) in galaxies using three-dimensional smoothed particle hydrodynamic simulations with new implementations of the momentum mechanical feedback, and restriction of accreted elements to those that are gravitationally bound to the BH. We also include the feedback from the X-ray radiation emitted by the BH, which heats the surrounding gas in the host galaxies, and adds radial momentum to the fluid. We perform simulations of isolated galaxies and merging galaxies and test various feedback models with the new treatment of the Bondi radius criterion. We find that overall the BH growth is similar to what has been obtained by earlier workers using the Springel, Di Matteo, & Hernquist algorithms. However, the outflowing wind velocities and mechanical energy emitted by winds are considerably higher (v_w ~ 1000-3000 km/s) compared to the standard thermal feedback model (v_w ~ 50-100 km/s). While the thermal feedback model emits only 0.1 % of BH released energy in winds, the momentum feedback model emits more than 30 % of the total energy released by the BH in winds. In the momentum feedback model, the degree of fluctuation in both radiant and wind output is considerably larger than in the standard treatments. We check that the new model of the BH mass accretion agrees with analytic results for the standard Bondi problem.

We revisit collisionless major and minor mergers of spheroidal galaxies in the context of the size evolution of elliptical galaxies. The simulations are performed as a series of mergers with mass-ratios of 1:1 and 1:10 for models representing pure bulges as well as bulges embedded in dark matter halos. For major and minor mergers, respectively, we identify and analyze two different processes, violent relaxation and stripping, leading to size evolution and a change of the dark matter fraction within the observable effective radius. Violent relaxation - which is the dominant mixing process for major mergers but less important for minor mergers - scatters relatively more dark matter particles than bulge particles to small. Stripping in minor mergers assembles stellar satellite particles at large radii in halo dominated regions of the massive host. This strongly increases the size of the bulge into regions with higher dark matter fractions leaving the inner host structure almost unchanged. A factor of two mass increase by minor mergers increases the dark matter fraction by 20 per cent. We present analytic corrections to simple one-component virial estimates for the evolution of the gravitational radii. If such a two-component system grows by minor mergers alone its size growth, $r_{\mathrm{e}} \propto M^\alpha$, reaches values of $\alpha \approx 2.4$, significantly exceeding the simple theoretical limit of $\alpha = 2$. For major mergers the sizes grow with $\alpha \lesssim 1$. Our results indicate that minor mergers of galaxies embedded in massive dark matter halos provide a potential mechanism for explaining the rapid size growth and the build-up of massive elliptical systems predicting significant dark matter fractions and radially biased velocity dispersions at large radii (abbreviated)

In this Letter, we investigate the potential power of the Cosmic Mach Number (CMN), which is the ratio between the mean velocity and the velocity dispersion of galaxies as a function of cosmic scales, to constrain cosmologies. We first measure the CMN from 5 catalogues of galaxy peculiar velocity surveys at low redshift (0.002<z<0.03), and use them to contrast cosmological models. Overall, current data is consistent with the WMAP7 LCDM model. We find that the CMN is highly sensitive to the growth of structure on scales 0.01<k<0.1 h/Mpc in Fourier space. Therefore, modified gravity models, and models with massive neutrinos, in which the structure growth generally deviates from that in the LCDM model in a scale-dependent way, can be well differentiated from the LCDM model using future CMN data.

We present a sample of nine high resolution cosmological simulations in the mass range of M_vir=7x10^11-4x10^12 M_sun starting from LambdaCDM initial conditions. Our simulations include primordial radiative cooling, photoionization, star formation, supernova II feedback, but exclude supernova driven winds and AGN feedback. The simulated galaxies assemble in two phases, with the initial growth dominated by compact (r<r_eff) in situ star formation fueled by cold, low entropy gas streams resulting in a very similar mean assembly redshift of z_{f,ins}~2.5 for the in situ stellar component in all galaxies. The late growth is dominated by accretion of old stars formed in subunits outside the main galaxy (r>r_eff) resulting in an assembly redshift of z_{f,acc}~0.5-1.5 with much larger scatter. We find a positive correlation between the fraction of accreted stars and the final mass of our galaxies. We show that gravitational feedback strongly suppresses late star formation in massive galaxies contributing to the observed galaxy color bimodality. The accretion of stellar material is also responsible for the observed size growth of early-type galaxies. In addition, we find that the dark matter fractions within the stellar half-mass radii continuously increase towards lower redshift from about f_DM~0.05 at z~3 to f_DM~0.1-0.3 at z=0. Furthermore, the logarithmic slope of the total density profile is nearly isothermal at the present-day (gamma'~1.9-2.2). Finally, the input of gravitational heating lowers the central dark matter densities in the galaxies, with the effect being smaller compared to simulations without supernova feedback.

We examine galaxy formation in a cosmological AMR simulation, which includes two high resolution boxes, one centered on a 3 \times 10^14 M\odot cluster, and one centered on a void. We examine the evolution of 611 massive (M\ast > 10^10M\odot) galaxies. We find that the fraction of the final stellar mass which is accreted from other galaxies is between 15 and 40% and increases with stellar mass. The accreted fraction does not depend strongly on environment at a given stellar mass, but the galaxies in groups and cluster environments are older and underwent mergers earlier than galaxies in lower density environments. On average, the accreted stars are ~2.5 Gyrs older, and ~0.15 dex more metal poor than the stars formed in-situ. Accreted stellar material typically lies on the outskirts of galaxies; the average half-light radius of the accreted stars is 2.6 times larger than that of the in-situ stars. This leads to radial gradients in age and metallicity for massive galaxies, in qualitative agreement with observations. Massive galaxies grow by mergers at a rate of approximately 2.6% per Gyr. These mergers have a median (mass-weighted) mass ratio less than 0.26 \pm 0.21, with an absolute lower limit of 0.20, for galaxies with M\ast ~ 10^12 M\odot. This suggests that major mergers do not dominate in the accretion history of massive galaxies. All of these results agree qualitatively with results from SPH simulations by Oser et al. (2010, 2012).

We present a new method for incorporating nonthermal pressure from bulk motions of gas into an analytic model of the intracluster medium in clusters of galaxies, which is based on a polytropic equation of state and hydrostatic equilibrium inside gravitational potential wells drawn from cosmological dark matter simulations. The pressure is allowed to have thermal and nonthermal components with different radial distributions; the overall level of nonthermal support is based on the dynamical state of the halo, such that it is lower in more relaxed clusters. This level is normalized by comparison to pressure profiles derived from X-ray observations, and to a high resolution hydrodynamical simulation. The nonthermal pressure fraction measured at r_500 is typically in the range 10-20%, increasing with cluster mass and with redshift. The resulting model cluster properties are in accord with Sunyaev-Zel'dovich (SZ) effect observations of clusters. Inclusion of nonthermal pressure reduces the expected angular power spectrum of SZ fluctuations in the microwave sky by 24%.

To facilitate the study of black hole fueling, star formation, and feedback in galaxies, we outline a method for treating the radial forces on interstellar gas due to absorption of photons by dust grains. The method gives the correct behavior in all of the relevant limits (dominated by the central point source; dominated by the distributed isotropic source; optically thin; optically thick to UV/optical; optically thick to IR) and reasonably interpolates between the limits when necessary. The method is explicitly energy conserving so that UV/optical photons that are absorbed are not lost, but are rather redistributed to the IR where they may scatter out of the galaxy. We implement the radiative transfer algorithm in a two-dimensional hydrodynamical code designed to study feedback processes in the context of early-type galaxies. We find that the dynamics and final state of simulations are measurably but only moderately affected by radiative forces on dust, even when assumptions about the dust-to-gas ratio are varied from zero to a value appropriate for the Milky Way. In simulations with high gas densities designed to mimic ULIRGs with a star formation rate of several hundred solar masses per year, dust makes a more substantial contribution to the dynamics and outcome of the simulation. We find that, despite the large opacity of dust to UV radiation, the momentum input to the flow from radiation very rarely exceeds L/c due to two factors: the low opacity of dust to the re-radiated IR and the tendency for dust to be destroyed by sputtering in hot gas environments. We also develop a simplification of our radiative transfer algorithm that respects the essential physics but is much easier to implement and requires a fraction of the computational cost.

We extend the black hole (BH) feedback models of Ciotti, Ostriker, and Proga to two dimensions. In this paper, we focus on identifying the differences between the one-dimensional and two-dimensional hydrodynamical simulations. We examine a normal, isolated $L_*$ galaxy subject to the cooling flow instability of gas in the inner regions. Allowance is made for subsequent star formation, Type Ia and Type II supernovae, radiation pressure, and inflow to the central BH from mildly rotating galactic gas which is being replenished as a normal consequence of stellar evolution. The central BH accretes some of the infalling gas and expels a conical wind with mass, momentum, and energy flux derived from both observational and theoretical studies. The galaxy is assumed to have low specific angular momentum in analogy with the existing one-dimensional case in order to isolate the effect of dimensionality. The code then tracks the interaction of the outflowing radiation and winds with the galactic gas and their effects on regulating the accretion. After matching physical modeling to the extent possible between the one-dimensional and two-dimensional treatments, we find essentially similar results in terms of BH growth and duty cycle (fraction of the time above a given fraction of the Eddington luminosity). In the two-dimensional calculations, the cool shells forming at 0.1--1 kpc from the center are Rayleigh--Taylor unstable to fragmentation, leading to a somewhat higher accretion rate, less effective feedback, and a more irregular pattern of bursting compared to the one-dimensional case.

We analyze 40 cosmological re-simulations of individual massive galaxies with present-day stellar masses of $M_{*} > 6.3 \times 10^{10} M_{\odot}$ in order to investigate the physical origin of the observed strong increase in galaxy sizes and the decrease of the stellar velocity dispersions since redshift $z \approx 2$. At present 25 out of 40 galaxies are quiescent with structural parameters (sizes and velocity dispersions) in agreement with local early type galaxies. At z=2 all simulated galaxies with $M_* \gtrsim 10^{11}M_{\odot}$ (11 out of 40) at z=2 are compact with projected half-mass radii of $\approx$ 0.77 ($\pm$0.24) kpc and line-of-sight velocity dispersions within the projected half-mass radius of $\approx$ 262 ($\pm$28) kms$^{-1}$ (3 out of 11 are already quiescent). Similar to observed compact early-type galaxies at high redshift the simulated galaxies are clearly offset from the local mass-size and mass-velocity dispersion relations. Towards redshift zero the sizes increase by a factor of $\sim 5-6$, following $R_{1/2} \propto (1+z)^{\alpha}$ with $\alpha = -1.44$ for quiescent galaxies ($\alpha = -1.12$ for all galaxies). The velocity dispersions drop by about one-third since $z \approx 2$, following $\sigma_{1/2} \propto (1+z)^{\beta}$ with $\beta = 0.44$ for the quiescent galaxies ($\beta = 0.37$ for all galaxies). The simulated size and dispersion evolution is in good agreement with observations and results from the subsequent accretion and merging of stellar systems at $z\lesssim 2$ which is a natural consequence of the hierarchical structure formation. A significant number of the simulated massive galaxies (7 out of 40) experience no merger more massive than 1:4 (usually considered as major mergers). On average, the dominant accretion mode is stellar minor mergers with a mass-weighted mass-ratio of 1:5. (abridged)

We investigate how environmental effects by gas stripping alter the growth of a super massive black hole (SMBH) and its host galaxy evolution, by means of 1D hydrodynamical simulations that include both mechanical and radiative AGN feedback effects. By changing the truncation radius of the gas distribution (R_t), beyond which gas stripping is assumed to be effective, we simulate possible environments for satellite and central galaxies in galaxy clusters and groups. The continuous escape of gas outside the truncation radius strongly suppresses star formation, while the growth of the SMBH is less affected by gas stripping because the SMBH accretion is primarily ruled by the density of the central region. As we allow for increasing environmental effects - the truncation radius decreasing from about 410 to 50 kpc - we find that the final SMBH mass declines from about 10^9 to 8 x 10^8 Msol, but the outflowing mass is roughly constant at about 2 x 10^10 Msol. There are larger change in the mass of stars formed, which declines from about 2 x 10^10 to 2 x 10^9 Msol, and the final thermal X-ray gas, which declines from about 10^9 to 5 x 10^8 Msol, with increasing environmental stripping. Most dramatic is the decline in the total time that the objects would be seen as quasars, which declines from 52 Myr (for R_t = 377 kpc) to 7.9 Myr (for R_t = 51 kpc). The typical case might be interpreted as a red and dead galaxy having episodic cooling flows followed by AGN feedback effects resulting in temporary transitions of the overall galaxy color from red to green or to blue, with (cluster) central galaxies spending a much larger fraction of their time in the elevated state than do satellite galaxies.(Abridged)

The ISM evolution of elliptical galaxies experiencing feedback from accretion onto a central black hole was studied recently with high-resolution 1D hydrodynamical simulations including radiative heating and pressure effects, a RIAF-like radiative efficiency, mechanical input from AGN winds, and accretion-driven starbursts. Here we focus on the observational properties of the models in the X-ray band (nuclear luminosity; hot ISM luminosity and temperature; temperature and brightness profiles during quiescence and during outbursts). The nuclear bursts last for ~10^7 yr, with a duty-cycle of a few X (10^-3-10^-2); the present epoch bolometric nuclear emission is very sub-Eddington. The ISM thermal luminosity \lx oscillates in phase with the nuclear one; this helps reproduce statistically the observed large \lx variation. In quiescence the temperature profile has a negative gradient; thanks to past outbursts, the brightness profile lacks the steep shape typical of inflowing models. Outbursts produce disturbances in these profiles. Most significantly, a hot bubble from shocked hot gas is inflated at the galaxy center; the bubble would be conical in shape, and show radio emission. The ISM resumes a smooth appearance on a time-scale of ~200 Myr; the duty-cycle of perturbances in the ISM is of the order of 5-10%. From the present analysis, additional input physics is important in the ISM-black hole coevolution, to fully account for the properties of real galaxies, as a confining external medium and a jet. The jet will reduce further the mass available for accretion (and then the Eddington ratio $l$), and may help, together with an external pressure, to produce flat or positive temperature gradient profiles (observed in high density environments). Alternatively, $l$ can be reduced if the switch from high to low radiative efficiency takes place at a larger $l$ than routinely assumed.

Hydrodynamic simulations of galaxies with active galactic nuclei (AGN) have typically employed feedback that is purely local: i.e., an injection of energy to the immediate neighborhood of the black hole. We perform GADGET-2 simulations of massive elliptical galaxies with an additional feedback component: an observationally calibrated X-ray radiation field which emanates from the black hole and heats gas out to large radii from the galaxy center. We find that including the heating and radiation pressure associated with this X-ray flux in our simulations enhances the effects which are commonly reported from AGN feedback. This new feedback model is twice as effective as traditional feedback at suppressing star formation, produces 3 times less star formation in the last 6 Gyr, and modestly lowers the final BH mass (30%). It is also significantly more effective than an X-ray background in reducing the number of satellite galaxies.

The importance of feedback (radiative and mechanical) from massive black holes at the centers of elliptical galaxies is not in doubt, given the well established relation among black hole mass and galaxy optical luminosity. Here, with the aid of high-resolution hydrodynamical simulations, we discuss how this feedback affects the hot ISM of isolated elliptical galaxies of different mass. The cooling and heating functions include photoionization plus Compton heating, the radiative transport equations are solved, and the mechanical feedback due to the nuclear wind is also described on a physical basis; star formation is considered. In the medium-high mass galaxies the resulting evolution is highly unsteady. At early times major accretion episodes caused by cooling flows in the recycled gas produced by stellar evolution trigger AGN flaring: relaxation instabilities occur so that duty cycles are small enough to account for the very small fraction of massive ellipticals observed to be in the QSO-phase, when the accretion luminosity approaches the Eddington luminosity. At low redshift all models are characterized by smooth, very sub-Eddington mass accretion rates. The mass accumulated by the central black hole is limited to range observed today, even though the mass lost by the evolving stellar population is roughly two order of magnitude larger than the black hole masses observed in elliptical galaxies.

Cosmological simulations of galaxy formation appear to show a two-phase character with a rapid early phase at z>2 during which in-situ stars are formed within the galaxy from infalling cold gas followed by an extended phase since z<3 during which ex-situ stars are primarily accreted. In the latter phase massive systems grow considerably in mass and radius by accretion of smaller satellite stellar systems formed at quite early times (z>3) outside of the virial radius of the forming central galaxy. These tentative conclusions are obtained from high resolution re-simulations of 39 individual galaxies in a full cosmological context with present-day virial halo masses ranging from 7e11 M_sun h^-1 < M_vir < 2.7e13 M_sun h^-1 and central galaxy masses between 4.5e10 M_sun h^-1 < M_* < 3.6e11 M_sun h^-1. The simulations include the effects of a uniform UV background, radiative cooling, star formation and energetic feedback from SNII. The importance of stellar accretion increases with galaxy mass and towards lower redshift. In our simulations lower mass galaxies ($M_* < 9e10 M_sun h^-1) accrete about 60 per cent of their present-day stellar mass. High mass galaxy ($M_* > 1.7e11 M_sun h^-1) assembly is dominated by accretion and merging with about 80 per cent of the stars added by the present-day. In general the simulated galaxies approximately double their mass since z=1. For massive systems this mass growth is not accompanied by significant star formation. The majority of the in-situ created stars is formed at z>2, primarily out of cold gas flows. We recover the observational result of archaeological downsizing, where the most massive galaxies harbor the oldest stars. We find that this is not in contradiction with hierarchical structure formation. Most stars in the massive galaxies are formed early on in smaller structures, the galaxies themselves are assembled late.

Even though the dark-matter power spectrum in the absence of biasing predicts a number density of halos n(M) ~ M^-2 (i.e. a Schechter alpha value of -2) at the low-mass end (M < 10^10 M_solar), hydrodynamic simulations have typically produced values for stellar systems in good agreement with the observed value alpha ~ -1. We explain this with a simple physical argument and show that an efficient external gas-heating mechanism (such as the UV background included in all hydro codes) will produce a critical halo mass below which halos cannot retain their gas and form stars. We test this conclusion with GADGET-2-based simulations using various UV backgrounds, and for the first time we also investigate the effect of an X-ray background. We show that at the present epoch alpha is depends primarily on the mean gas temperature at the star-formation epoch for low-mass systems (z <~ 3): with no background we find alpha ~ -1.5, with UV only alpha ~ -1.0, and with UV and X-rays alpha ~ -0.75. We find the critical final halo mass for star formation to be ~4x10^8 M_solar with a UV background and ~7x10^8 M_solar with UV and X-rays.

The deposition of mechanical feedback from a supermassive black hole (SMBH) in an active galactic nucleus (AGN) into the surrounding galaxy occurs via broad-line winds which must carry mass and radial momentum as well as energy. The effect can be summarized by the dimensionless parameter $\eta=dot{M_outflow}/dot{M_accretion}= (2 \epsilon_w c^2)/v_w^2$ where ($\epslion_w \equiv dot{E}_w/(dot{M_accretion} c^2)$) is the efficiency by which accreted matter is turned into wind energy in the disc surrounding the central SMBH. The outflowing mass and omentum are proportional to $\eta$, and many prior treatments have essentially assumed that $\eta=0$. We perform one- and two-dimensional simulations and find that the growth of the central SMBH is very sensitive to the inclusion of the mass and momentum driving but is insensitive to the assumed mechanical efficiency. For example in representative calculations, the omission of momentum and mass feedback leads to an hundred fold increase in the mass of the SMBH to over $10^{10} \Msun$. When allowance is made for momentum driving, the final SMBH mass is much lower and the wind efficiencies which lead to the most observationally acceptable results are relatively low with $\epsilon_w \lesssim 10^{-4}$.

We present templates for the Sunyaev-Zel'dovich (SZ) angular power spectrum based on four models for the nonlinear gas distribution. The frequency-dependent SZ temperature fluctuations, with thermal (TSZ) and kinetic (KSZ) contributions, are calculated by tracing through a dark matter simulation, processed to include gas in dark matter halos and in the filamentary intergalactic medium. Different halo gas models are compared to study how star formation, energetic feedback, and nonthermal pressure support influence the angular power spectrum. The standard model has been calibrated to reproduce the stellar and gas fractions and X-ray scaling relations measured from low redshift clusters and groups. The other models illustrate the current theoretical and empirical uncertainties relating to properties of the intracluster medium. Relative to the standard model, their angular power spectra differ by approximately 50% (TSZ), 20% (KSZ), and 40% (SZ at 148 GHz) for l=3000, sigma_8=0.8, and homogeneous reionization at z=10. The angular power spectrum decreases in amplitude as gas mass and binding energy is removed through star formation, and as gas is pushed out to larger radii by energetic feedback. With nonthermal pressure support, less pressure is required to maintain hydrostatic equilibrium, thus reducing the thermal contribution to the SZ power. We also calculate the SZ templates as a function of sigma_8 and quantify this dependence. Assuming C_l is proprotional to (sigma_8/0.8)^alpha, the effective scaling index ranges from 7<alpha_tsz<9, 4.5<alpha_ksz<5.5, and 6.5<alpha_sz(148 GHz)<8 at l=3000 for 0.6<sigma_8<1. The template spectra are publicly available and can be used when fitting for the SZ contribution to the cosmic microwave background on arcminute scales.

A novel statistic is proposed to examine the hypothesis that all cluster galaxies are drawn from the same luminosity distribution (LD). In such a "statistical model" of galaxy LD, the brightest cluster galaxies (BCGs) are simply the statistical extreme of the galaxy population. Using a large sample of nearby clusters, we show that BCGs in high luminosity clusters (e.g., L_tot > 4x10^11 L_sun) are unlikely (probability <3x10^-4) to be drawn from the LD defined by all red cluster galaxies more luminous than M_r=-20. On the other hand, BCGs in less luminous clusters are consistent with being the statistical extreme. Applying our method to the second brightest galaxies, we show that they are consistent with being the statistical extreme, which implies that the BCGs are also distinct from non-BCG luminous, red, cluster galaxies. We point out some issues with the interpretation of the classical tests proposed by Tremaine & Richstone (1977) that are designed to examine the statistical nature of BCGs, investigate the robustness of both our statistical test and those of TR against difficulties in photometry of galaxies of large angular size, and discuss the implication of our findings on surveys that use the luminous red galaxies to measure the baryon acoustic oscillation features in the galaxy power spectrum.

We find, from high-resolution hydro simulations, that winds from AGN effectively heat the inner parts (~100 pc) of elliptical galaxies, reducing infall to the central SMBH; and radiative (photoionization and X-ray) heating reduces cooling flows at the kpc scale. Including both types of feedback with (peak) efficiencies of 3 10^{-4} < epsilon_mech < 10^{-3} and of epsilon_rad ~10^{-1.3} respectively, produces systems having duty-cycles, central SMBH masses, X-ray luminosities, optical light profiles, and E+A spectra in accord with the broad suite of modern observations of massive elliptical systems. Our main conclusion is that mechanical feedback (including all three of energy, momentum and mass) is necessary but the efficiency, based on several independent arguments must be a factor of 10 lower than is commonly assumed. Bursts are frequent at z>1 and decline in frequency towards the present epoch as energy and metal rich gas are expelled from the galaxies into the surrounding medium. For a representative galaxy of final stellar mass ~3 10^{11} Msun, roughly 3 10^{10} Msun of recycled gas has been added to the ISM since z~2 and, of that, roughly 63% has been expelled from the galaxy, 19% has been converted into new metal rich stars in the central few hundred parsecs, and 2% has been added to the central SMBH, with the remaining 16% in the form hot X-ray emitting ISM. The bursts occupy a total time of ~170 Myr, which is roughly 1.4% of the available time. Of this time, the central SMBH would be seen as an UV or optical source for ~45% and ~71$% of the time, respectively. Restricting to the last 8.5 Gyr, the burst occupy ~44 Myr, corresponding to a fiducial duty-cycle of ~5 10^{-3}.

We examine two extreme models for the build-up of the stellar component of luminous elliptical galaxies. In one case, we assume the build-up of stars is dissipational, with centrally accreted gas radiating away its orbital and thermal energy; the dark matter halo will undergo adiabatic contraction and the central dark matter density profile will steepen. For the second model, we assume the central galaxy is assembled by a series of dissipationless mergers of stellar clumps that have formed far from the nascent galaxy. In order to be accreted, these clumps lose their orbital energy to the dark matter halo via dynamical friction, thereby heating the central dark matter and smoothing the dark matter density cusp. The central dark matter density profiles differ drastically between these models. For the isolated elliptical galaxy, NGC 4494, the central dark matter densities follow the power-laws r^(-0.2) and r^(-1.7) for the dissipational and dissipationless models, respectively. By matching the dissipational and dissipationless models to observations of the stellar component of elliptical galaxies, we examine the relative contributions of dissipational and dissipationless mergers to the formation of elliptical galaxies and look for observational tests that will distinguish between these models. Comparisons to strong lensing brightest cluster galaxies yield median M*/L_B ratios of 2.1+/-0.8 and 5.2+/-1.7 at z=0.39 for the dissipational and dissipationless models, respectively. For NGC 4494, the best-fit dissipational and dissipationless models have M*/L_B=2.97 and 3.96. Comparisons to expected stellar mass-to-light ratios from passive evolution and population syntheses appear to rule out a purely dissipational formation mechanism for the central stellar regions of giant elliptical galaxies.

By using high-resolution 1D hydrodynamical simulations, we investigate the effects of purely mechanical feedback from super massive black holes (SMBHs) in the evolution of elliptical galaxies for a broad range of feedback efficiencies and compare the results to four major observational constraints. In particular, we focus on 1) the central black hole to stellar mass ratio of the host galaxy, 2) the lifetime of the luminous quasar phase, 3) the mass of stars formed in the host galaxy within the last Gyr, and 4) the X-ray luminosity of the hot diffuse gas. As a result, we try to pin down the most successful range of mechanical feedback efficiencies. We find that while low feedback efficiencies result in too much growth of the SMBH, high efficiencies totally blow out the hot interstellar gas, and the models are characterized by very low thermal X-ray luminosity well below the observed range. The net lifetime of the quasar phase is strongly coupled to the mass ratio between SMBH and its host galaxy, while the X-ray luminosity is generally correlated to the recent star formation within the last Gyr. When considering the popularly adopted model of the constant feedback efficiency, the feedback energy deposited into the ambient medium should be more than 0.01% of the SMBH accretion energy to be consistent with the SMBH mass to stellar mass ratio in the local universe. Yet, the X-ray luminosity of the hot gas favors about 0.005% of the accretion energy as the mechanical AGN feedback energy. We conclude that the purely mechanical feedback mode is unlikely to be simultaneously compatible with all four observable tests, even allowing a broad range of feedback efficiencies, and that including both radiative and mechanical feedback together may be a solution to comply the observational constraints. [abridged]

Both radiative and mechanical feedback from Active Galactic Nuclei have been found to be important for the evolution of elliptical galaxies. We compute how a shock may be driven from a central black hole into the gaseous envelope of an elliptical galaxy by mechanical as well as radiative feedback (in the form of nuclear winds) using high resolution 1-D hydrodynamical simulations. We calculate the synchrotron emission from the electron cosmic rays accelerated by the shocks (not the jets) in the central part of elliptical galaxies, and we also study the synchrotron spectrum's evolution using the standard diffusive shock acceleration mechanism, which is routinely applied to the scaled volume case of supernova remnants. We find good agreement quantitatively between the synchrotron radio emission produced via this mechanism with extant observations of elliptical galaxies which are undergoing outbursts. Additionally, we also find that synchrotron optical and X-ray emission can co-exist inside elliptical galaxies during a certain phase of evolution subsequent to central outbursts. In fact, our calculations predict a synchrotron luminosity of $\sim 1.3\times 10^6 L_{\odot}$ at the frequency 5 GHz (radio band), of $\sim 1.1\times 10^6 L_{\odot}$ at $4.3\times10^{14}$ Hz (R band, corresponding to the absolute magnitude -10.4), and of $\sim 1.5\times 10^{7} L_{\odot}$ at $2.4\times10^{17}$ Hz (soft X-ray, 0.5 -- 2.0 keV band).

We create realistic, full-sky, half-arcminute resolution simulations of the microwave sky matched to the most recent astrophysical observations. The primary purpose of these simulations is to test the data reduction pipeline for the Atacama Cosmology Telescope (ACT) experiment; however, we have widened the frequency coverage beyond the ACT bands to make these simulations applicable to other microwave background experiments. Some of the novel features of these simulations are that the radio and infrared galaxy populations are correlated with the galaxy cluster populations, the CMB is lensed by the dark matter structure in the simulation via a ray-tracing code, the contribution to the thermal and kinetic Sunyaev-Zel'dovich (SZ) signals from galaxy clusters, groups, and the IGM has been included, and the gas prescription to model the SZ signals matches the most recent X-ray observations. Regarding the contamination of cluster SZ flux by radio galaxies, we find for 148 GHz (90 GHz) only 3% (4%) of halos have their SZ decrements contaminated at a level of 20% or more. We find the contamination levels higher for infrared galaxies. However, at 90 GHz, less than 20% of clusters with M_{200} > 2.5 x 10^{14} Msun and z<1.2 have their SZ decrements filled in at a level of 20% or more. At 148 GHz, less than 20% of clusters with M_{200} > 2.5 x 10^{14} Msun and z<0.8 have their SZ decrements filled in at a level of 50% or larger. Our models also suggest that a population of very high flux infrared galaxies, which are likely lensed sources, contribute most to the SZ contamination of very massive clusters at 90 and 148 GHz. These simulations are publicly available and should serve as a useful tool for microwave surveys to cross-check SZ cluster detection, power spectrum, and cross-correlation analyses.

We find that the amount and nature of the assumed ionizing background can strongly affect galaxy formation and evolution. Galaxy evolution simulations typically incorporate an ultraviolet background which falls off rapidly above z=3; e.g., that of Haardt & Madau (1996). However, this decline may be too steep to fit the WMAP constraints on electron scattering optical depth or observations of intermediate redshift (z ~ 2-4) Ly-alpha forest transmission. As an alternative, we present simulations of the cosmological formation of individual galaxies with UV backgrounds that decline more slowly at high redshift: both a simple intensity rescaling and the background recently derived by Faucher-Giguere (2009), which softens the spectrum at higher redshifts. We also test an approximation of the X-ray background with a similar z-dependence. We find for the test galaxies that an increase in either the intensity or hardness of ionizing radiation generically pushes star formation towards lower redshifts: although overall star formation in the simulation boxes is reduced by 10-25%, the galaxies show a factor of ~2 increase in the fraction of stars within a 30 kpc radius that are formed after z=1. Other effects include late gas inflows enhanced up to 30 times, stellar half-mass radii decreased by up to 30%, central velocity dispersions increased up to 40%, and a strong reduction in substructure. The magnitude of the effects depends on the environmental/accretion properties of the particular galaxy.

Using a high resolution hydrodynamical cosmological simulation of the formation of a massive spheroidal galaxy we show that elliptical galaxies can be very compact and massive at high redshift in agreement with recent observations. Accretion of stripped in-falling stellar material increases the size of the system with time and the central concentration is reduced by dynamical friction of the surviving stellar cores. In a specific case of a spheroidal galaxy with a final stellar mass of $1.5 \times 10^{11} M_{\odot}$ we find that the effective radius $r_e$ increases from $0.7 \pm 0.2 \rm kpc$ at z = 3 to $r_e = 2.4 \pm 0.4 \rm kpc$ at z = 0 with a concomitant decrease in the effective density of an order of magnitude and a decrease of the central velocity dispersion by approximately 20% over this time interval. A simple argument based on the virial theorem shows that during the accretion of weakly bound material (minor mergers) the radius can increase as the square of the mass in contrast to the usual linear rate of increase for major mergers. By undergoing minor mergers compact high redshift spheroids can evolve into present-day systems with sizes and concentrations similar to observed local ellipticals. This indicates that minor mergers may be the main driver for the late evolution of sizes and densities of early-type galaxies.

The state of the hot gas in clusters of galaxies is investigated with a set of model clusters, created by assuming a polytropic equation of state (Gamma=1.2) and hydrostatic equilibrium inside gravitational potential wells drawn from a dark matter simulation. Star formation, energy input, and nonthermal pressure support are included. To match the gas fractions seen in non-radiative hydrodynamical simulations, roughly 5% of the binding energy of the dark matter must be transferred to the gas during cluster formation; the presence of nonthermal pressure support increases this value. In order to match X-ray observations, scale-free behavior must be broken. This can be due to either variation of the efficiency of star formation with cluster mass M_500, or the input of additional energy proportional to the formed stellar mass M_F. These two processes have similar effects on X-ray scalings. If 9% of the gas is converted into stars, independent of cluster mass, then feedback energy input of 1.2e-5*M_Fc^2 (or ~1.0 keV per particle) is required to match observed clusters. Alternatively, if the stellar mass fraction varies as M_500^-0.26 then a lower feedback of 4e-6*M_Fc^2 is needed, and if the stellar fraction varies as steeply as M_500^-0.49 then no additional feedback is necessary. The model clusters reproduce the observed trends of gas temperature and gas mass fraction with cluster mass, as well as observed entropy and pressure profiles; thus they provide a calibrated basis with which to interpret upcoming SZ surveys. One consequence of the increased gas energy is that the baryon fraction inside the virial radius is less than roughly 90% of the cosmic mean, even for the most massive clusters.

The centers of elliptical galaxies host supermassive black holes that significantly affect the surrounding interstellar medium through feedback resulting from the accretion process. The evolution of this gas and of the nuclear emission during the galaxies' lifetime has been studied recently with high-resolution hydrodynamical simulations. These included gas cooling and heating specific for an average AGN spectral energy distribution, a radiative efficiency declining at low mass accretion rates, and mechanical coupling between the hot gas and AGN winds. Here we present a short summary of the observational properties resulting from the simulations, focussing on 1) the nuclear luminosity; 2) the global luminosity and temperature of the hot gas; 3) its temperature profile and X-ray brightness profile. These properties are compared with those of galaxies of the local universe, pointing out the successes of the adopted feedback and the needs for new input in the simulations.

The importance of the radiative feedback from SMBHs at the centers of elliptical galaxies is not in doubt, given the well established relations among electromagnetic output, black hole mass and galaxy optical luminosity. In addition, feedback due to mechanical and thermal deposition of energy from jets and winds emitted by the accretion disk around the central SMBH is also expected to occur. In this paper we improve and extend the accretion and feedback physics explored in our previous papers to include also a physically motivated mechanical feedback. We study the evolution of an isolated elliptical galaxy with the aid of a high-resolution 1-D hydrodynamical code, where the cooling and heating functions include photoionization and Compton effects, and restricting to models which include only radiative or only mechanical feedback. We confirm that for Eddington ratios above 0.01 both the accretion and radiative output are forced by feedback effects to be in burst mode, so that strong intermittencies are expected at early times, while at low redshift the explored models are characterized by smooth, very sub-Eddington mass accretion rates punctuated by rare outbursts. However, the explored models always fail some observational tests. If we assume the high mechanical efficiency of 10^{-2.3}, we find that most of the gas is ejected from the galaxy, the resulting X-ray luminosity is far less than is typically observed and little SMBH growth occurs. But models with low enough mechanical efficiency to accomodate satisfactory SMBH growth tend to allow too strong cooling flows and leave galaxies at z=0 with E+A spectra more frequently than is observed. We conclude that both types of feedback are required. Models with combined feedback are explored in a forthcoming paper [abridged]

We study the thermal formation history of four simulated galaxies that were shown in Naab et al. (2007) to reproduce a number of observed properties of elliptical galaxies. The temperature of the gas in the galaxies is steadily increasing with decreasing redshift, although much of the gas has a cooling time shorter than the Hubble time. The gas is being heated and kept hot by gravitational heating processes through the release of potential energy from infalling stellar clumps. The energy is dissipated in supersonic collisions of infalling gas lumps with the ambient gas and through the dynamical capturing of satellite systems causing gravitational wakes that transfer energy to the surrounding gas. Furthermore dynamical friction from the infalling clumps pushes out dark matter, lowering the central dark matter density by up to a factor of two from z=3 to z=0. In galaxies in which the late formation history (z<2) is dominated by minor merging and accretion the energy released (E~5x10^{59} ergs) from gravitational feedback is sufficient to form red and dead elliptical galaxies by z~1 even in the absence of supernova and AGN feedback.

Over the last decade optical spectroscopic surveys have characterized the low redshift galaxy population and uncovered populations of star-forming galaxies back to z ~ 7. This work has shown that the primary epoch of galaxy building and black hole growth occurs at redshifts of 2 to 3. The establishment of the concordance LCDM cosmology shifted the focus of galaxy population studies from constraining cosmological parameters to characterizing the processes which regulate the formation and evolution of galaxies.In the next decade, high redshift observers will attempt to formulate a coherent evolutionary picture connecting galaxies in the high redshift Universe to galaxies today. In order to link galaxy populations at different redshifts, we must not only characterize their evolution in a systematic way, we must establish which physical processes are responsible for it. Considerable progress has already been made in understanding how galaxies evolved from z ~ 1 to the present day. Large spectroscopic surveys in the near infrared are required to push these studies back towards the main epoch of galaxy building. Only then will we understand the full story of the formation of L* galaxies like our own Milky Way. A large near-IR spectroscopic survey will also provide the calibration needed to avoid systematics in the large photometric programs proposed to study the nature of dark matter and dark energy. We provide an outline design for a multi-object 0.4 to 1.8 micron spectrograph, which could be placed on an existing telescope, and which would allow a full characterization of the galaxy population out to z ~ 2. We strongly recommend a serious further study to design a real instrument, which will be required for galaxy formation studies to advance to the next frontier.

Great strides have been made in the last two decades in determining how galaxies evolve from their initial dark matter seeds to the complex structures we observe at z=0. The role of mergers has been documented through both observations and simulations, numerous satellites that may represent these initial dark matter seeds have been discovered in the Local Group, high redshift galaxies have been revealed with monstrous star formation rates, and the gaseous cosmic web has been mapped through absorption line experiments. Despite these efforts, the dark matter simulations that include baryons are still unable to accurately reproduce galaxies. One of the major problems is our incomplete understanding of how a galaxy accretes its baryons and subsequently forms stars. Galaxy formation simulations have been unable to accurately represent the required gas physics on cosmological timescales, and observations have only just begun to detect the star formation fuel over a range of redshifts and environments. How galaxies obtain gas and subsequently form stars is a major unsolved, yet tractable problem in contemporary extragalactic astrophysics. In this paper we outline how progress can be made in this area in the next decade.

A popular formation scenario for giant elliptical galaxies proposes that they might have formed from binary mergers of disk galaxies. Difficulties with the scenario that emerged from earlier studies included providing the necessary stellar mass and metallicity, maintaining the tight color-magnitude relation and avoiding phase space limits. In this paper we revisit the issue and put constraints on the binary disc merger scenario based on the stellar populations of disc galaxies. We draw the following conclusions: Low redshift collisionless or gaseous mergers of present day Milky Way like disc galaxies do not form present day elliptical galaxies. Binary mergers of the progenitors of present day Milky Way like disc galaxies can have evolved into intermediate mass elliptical galaxies ($M < M_*$) if they have merged earlier than $\approx$ 3-4 Gyrs ago. Assuming that most present day disk galaxies formed in a similar way to the Milky Way model presented here, more massive giant ellipticals in general can not have formed from binary mergers of the progenitors of present day disc galaxies. A major reason for these conclusions is that the mass in metals of typical disk galaxy is approximately a factor of 4-8 smaller than the mass in metals of a typical early-type galaxy and this ratio grows to larger values with increasing redshift.

The hot accretion flow is usually optically thin in the radial direction, therefore the photons produced at one radius can travel for a long distance without being absorbed. These photons thus can heat or cool electrons at other radii via Compton scattering. This effect has been ignored in most previous works on hot accretion flows and is the focus of this paper. If the mass accretion rate is described by $\dot{M}=\dot{M}_0(r/r_{\rm out})^{0.3}$ and $r_{\rm out}=10^4r_{\rm s}$, we find that the Compton scattering will play a cooling and heating role at $r\la 5\times 10^3 r_{\rm s}$ and $r\ga 5\times 10^3 r_{\rm s}$, respectively. Specifically, when $\dot{M}_0>0.1L_{\rm Edd}/c^2$, the Compton cooling rate is larger than the local viscous heating rate at certain radius; therefore the cooling effect is important. When $\dot{M}_0>2L_{\rm Edd}/c^2$, the heating effect at $r_{\rm out}$ is important. We can obtain the self-consistent steady solution with the global Compton effect included only if $\dot{M}_0\la L_{\rm Edd}/c^2$ for $r_{\rm out}=50r_{\rm s}$, which corresponds to $L\la 0.02L_{\rm Edd}$. Above this rate the Compton cooling is so strong at the inner region that hot solutions can not exist. On the other hand, for $r_{\rm out}= 10^5r_{\rm s}$, we can only get the self-consistent solution when $\dot{M}_0\la L_{\rm Edd}/c^2$ and $L<0.01L_{\rm Edd}$. Above this accretion rate the equilibrium temperature of electrons at $r_{\rm out}$ is higher than the virial temperature as a result of strong Compton heating, so the accretion is suppressed. In this case the activity of the black hole will likely "oscillate" between an active and an inactive phases, with the oscillation timescale being the radiative timescale of the gas at $r_{\rm out}$.

A variety of physical heating mechanisms are combined with radiative cooling to explore, via one dimensional hydrodynamic simulations, the expected thermal properties of the intracluster medium (ICM) in the context of the cooling flow problem. Energy injection from type Ia supernovae, thermal conduction, and dynamical friction (DF) from orbiting satellite galaxies are considered. The novel feature of this work is the exploration of a wide range of efficiencies of each heating process. While the latter two can provide a substantial amount of energy, neither mechanism operating alone can produce nor maintain an ICM in thermal balance over cosmological timescales, in stark contrast with observations. For simulated clusters with initially isothermal temperature profiles, both mechanisms acting in combination result in long-term thermal balance for a range of ICM temperatures and for central electron densities less than n_e~0.02 cm^-3; at greater densities catastrophic cooling invariably occurs. Furthermore, these heating mechanisms can neither produce nor maintain clusters with a declining temperature profile in the central regions, implying that the observed "cooling-core'' clusters, which have such declining temperature profiles, cannot be maintained with these mechanisms alone. Thus, while there appears to be an abundant supply of energy capable of heating the ICM in clusters, it is extremely difficult for the energy deposition to occur in such a way that the ICM remains in thermal balance over cosmological time-scales. These results strongly suggest that a more dynamic heating process such as feedback from a central black hole is required to generate the properties of observed intracluster media. (ABRIDGED)

We present results of an attempt to include more detailed gas physics motivated from hydrodynamical simulations within semi-analytic models (SAM) of galaxy formation, focusing on the role that environmental effects play. The main difference to previous SAMs is that we include 'gravitational' heating of the intra-cluster medium (ICM) by the net surplus of gravitational potential energy released from gas that has been stripped from infalling satellites. Gravitational heating appears to be an efficient heating source able to prevent cooling in environments corresponding to dark matter halos more massive than $\sim 10^{13} $M$_{\odot}$. The energy release by gravitational heating can match that by AGN-feedback in massive galaxies and can exceed it in the most massive ones. However, there is a fundamental difference in the way the two processes operate. Gravitational heating becomes important at late times, when the peak activity of AGNs is already over, and it is very mass dependent. This mass dependency and time behaviour gives the right trend to recover down-sizing in the star-formation rate of massive galaxies. Abridged...

We study the axisymmetric, time-dependent hydrodynamics of rotating flows that are under the influence of supermassive black hole gravity and radiation from an accretion disk surrounding the black hole. This work is an extension of the earlier work presented by Proga, where nonrotating flows were studied. Here, we consider effects of rotation, a position-dependent radiation temperature, density at large radii, and uniform X-ray background radiation. As in the non-rotating case, the rotating flow settles into a configuration with two components (1) an equatorial inflow and (2) a bipolar inflow/outflow with the outflow leaving the system along the pole. However, with rotation the flow does not always reach a steady state. In addition, rotation reduces the outflow collimation and the outward flux of mass and kinetic energy. Moreover rotation increases the outward flux of the thermal energy and can lead to fragmentation and time-variability of the outflow. We also show that a position-dependent radiation temperature can significantly change the flow solution. In particular, the inflow in the equatorial region can be replaced by a thermally driven outflow. Generally, as it have been discussed and shown in the past, we find that self-consistently determined preheating/cooling from the quasar radiation can significantly reduce the rate at which the central BH is fed with matter. However, our results emphasize also a little appreciated feature. Namely, quasar radiation drives a non-spherical, multi-temperature and very dynamic flow. These effects become dominant for luminosities in excess of 0.01 of the Eddington luminosity.

We investigate the effect of non-evaporating primordial black holes (PBHs) on the ionization and thermal history of the universe. X-rays emitted by gas accretion onto PBHs modify the cosmic recombination history, producing measurable effects on the spectrum and anisotropies of the Cosmic Microwave Background (CMB). Using the third-year WMAP data and FIRAS data we improve existing upper limits on the abundance of PBHs with masses >0.1 Msun by several orders of magnitude. Fitting WMAP3 data with cosmological models that do not allow for non-standard recombination histories, as produced by PBHs or other early energy sources, may lead to an underestimate of the best-fit values of the amplitude of linear density fluctuations (sigma_8) and the scalar spectral index (n_s). Cosmological parameter estimates are affected because models with PBHs allow for larger values of the Thomson scattering optical depth, whose correlation with other parameters may not be correctly taken into account when PBHs are ignored. Values of tau_e=0.2, n_s=1 and sigma_8=0.9 are allowed at 95% CF. This result that may relieve recent tension between WMAP3 data and clusters data on the value of sigma_8. PBHs may increase the primordial molecular hydrogen abundance by up to two orders of magnitude, this promoting cooling and star formation. The suppression of galaxy formation due to X-ray heating is negligible for models consistent with the CMB data. Thus, the formation rate of the first galaxies and stars would be enhanced by a population of PBHs.

The statistics of giant arcs and large separation lensed quasars provide powerful constraints for the parameters of the underlying cosmological model. So far, most investigations have been carried out using pure dark matter simulations. Here we present a recipe for including the effects of baryon cooling (i.e. large galaxy formation) in dark matter N-body simulations that is consistent with observations of massive galaxies. Then we quantitatively compare lensing with and without applying this baryon correction to the pure dark matter case. Including the baryon correction significantly increases the frequency of giant arcs and lensed quasars, particularly on scales of 10 arcsec and smaller: the overall frequency of multiple images increases by about 25% for source redshifts between z_s = 1.5 and 7.5 and splittings larger than about 3 arcsec. The baryon rearrangement also slightly increases the fraction of quadruple images over doubles.

Using the full, three-dimensional potential of galaxy cluster halos (drawn from an N-body simulation of the current, most favored cosmology), the distribution of the X-ray emitting gas is found by assuming a polytropic equation of state and hydrostatic equilibrium, with constraints from conservation of energy and pressure balance at the cluster boundary. The resulting properties of the gas for these simulated redshift zero clusters (the temperature distribution, mass-temperature and luminosity-temperature relations, and the gas fraction) are compared with observations in the X-ray of nearby clusters. The observed properties are reproduced only under the assumption that substantial energy injection from non-gravitational sources has occurred. Our model does not specify the source, but star formation and AGN may be capable of providing this energy, which amounts to 3 to 5 x10^{-5} of the rest mass in stars (assuming ten percent of the gas initially in the cluster forms stars). With the method described here it is possible to generate realistic X-ray and Sunyaev-Zel'dovich cluster maps and catalogs from N-body simulations, with the distributions of internal halo properties (and their trends with mass, location, and time) taken into account.

The importance of the radiative feedback from massive black holes at the centers of elliptical galaxies is not in doubt, given the well established relations among electromagnetic output, black hole mass and galaxy optical luminosity. We show how this AGN radiative output affects the hot ISM of an isolated elliptical galaxy with the aid of a high-resolution hydrodynamical code, where the cooling and heating functions include photoionization plus Compton heating. We find that radiative heating is a key factor in the self-regulated coevolution of massive black holes and their host galaxies and that 1) the mass accumulated by the central black hole is limited by feedback to the range observed today, and 2) relaxation instabilities occur so that duty cycles are small enough (~0.03) to account for the very small fraction of massive ellipticals observed to be in the "on" -QSO- phase, when the accretion luminosity approaches the Eddington luminosity. The duty cycle of the hot bubbles inflated at the galaxy center during major accretion episodes is of the order of 0.1-0.4. Major accretion episodes caused by cooling flows in the recycled gas produced by normal stellar evolution trigger nuclear starbursts coincident with AGN flaring. During such episodes the central sources are often obscured; but overall, in the bursting phase (1<z<3), the duty cycle of the black hole in its "on" phase is of the order of percents and it is unobscured approximately one-third of the time. Mechanical energy output from non-relativistic gas winds integrates to 2.3 10^{59} erg, with most of it caused by broadline AGN outflows. [abridged]

Cosmological shock waves result from supersonic flow motions induced by hierarchical clustering of nonlinear structures in the universe. These shocks govern the nature of cosmic plasma through thermalization of gas and acceleration of nonthermal, cosmic-ray (CR) particles. We study the statistics and energetics of shocks formed in cosmological simulations of a concordance $\Lambda$CDM universe, with a special emphasis on the effects of non-gravitational processes such as radiative cooling, photoionization/heating, and galactic superwind feedbacks. Adopting an improved model for gas thermalization and CR acceleration efficiencies based on nonlinear diffusive shock acceleration calculations, we then estimate the gas thermal energy and the CR energy dissipated at shocks through the history of the universe. Since shocks can serve as sites for generation of vorticity, we also examine the vorticity that should have been generated mostly at curved shocks in cosmological simulations. We find that the dynamics and energetics of shocks are governed primarily by the gravity of matter, so other non-gravitational processes do not affect significantly the global energy dissipation and vorticity generation at cosmological shocks. Our results reinforce scenarios in which the intracluster medium and warm-hot intergalactic medium contain energetically significant populations of nonthermal particles and turbulent flow motions.

Convection-dominated accretion flows (CDAF) are expected to have a shallower density profile and a higher radiation efficiency as compared to advection-dominated accretion flows (ADAF). Both solutions have been developed to account for the observed properties of the low luminosity, high temperature X-ray sources believed to involve accretion onto massive black holes. Self-similar CDAFs also have steeper poloidal density gradients and temperatures close to the virial temperature. All these characteristics make CDAFs more capable of producing polar outflows driven by Compton heating as compared to other classical accretion disks. We investigate the conditions for producing such outflows in CDAFs and look for the mass accretion rate, or, equally, the luminosity of CDAFs for which such outflows will exist. When the electron temperature saturates around 10^11 K at the inner region, polar outflows are probable for 8x10^-7 <~ L/L_E <~ 4x10^-5, where L_E is the Eddington luminosity. Outflows are well collimated with small opening angles. The luminosity range for which outflow solutions exist is narrower for lower electron temperature flows and disappears completely for electron temperature <~ 6x10^9 K. When the magnetic field is present, we find that outflows are possible if the magnetic field is less than from 10% to 1% of the equipartition field. We also find that outflows are more likely to be produced when the viscosity parameter alpha is small. The tendency for jet-like collimated outflows for these solutions is presumably astrophysically relevant given the high frequency of jets from AGNs.

We present a new definition of subhalos in dissipationless dark matter N-body simulations, based on the coherent identification of their dynamically bound constituents. Whereas previous methods of determining the energetically bound components of a subhalo ignored the contribution of all the remaining particles in the halo (those not geometrically or dynamically associated with the subhalo), our method allows for all the forces, both internal and external, exerted on the subhalo. We demonstrate, using the output of a simulation at different timesteps, that our new method is more accurate at identifying the bound mass of a subhalo. We then compare our new method to previously adopted means of identifying subhalos by applying each to a sample of 1838 virialized halos extracted from a high resolution cosmological simulation. We find that the subhalo distributions are similar in each case, and that the increase in the binding energy of a subhalo from including all the particles located within it is almost entirely balanced by the losses due to the external forces; the net increase in the mass fraction of subhalos is roughly 10%, and the extra substructures tending to reside in the inner parts of the system. Finally, we compare the subhalo populations of halos to the sub-subhalo populations of subhalos, finding the two distributions to be similar. This is a new and interesting result, suggesting a self-similarity in the hierarchy substructures within cluster mass halos.

We revisit the longstanding question of whether first brightest cluster galaxies are statistically drawn from the same distribution as other cluster galaxies or are "special", using the new non-parametric, empirically based model presented in Vale&Ostriker (2006) for associating galaxy luminosity with halo/subhalo masses. We introduce scatter in galaxy luminosity at fixed halo mass into this model, building a conditional luminosity function (CLF) by considering two possible models: a simple lognormal and a model based on the distribution of concentration in haloes of a given mass. We show that this model naturally allows an identification of halo/subhalo systems with groups and clusters of galaxies, giving rise to a clear central/satellite galaxy distinction. We then use these results to build up the dependence of brightest cluster galaxy (BCG) magnitudes on cluster luminosity, focusing on two statistical indicators, the dispersion in BCG magnitude and the magnitude difference between first and second brightest galaxies. We compare our results with two simple models for BCGs: a statistical hypothesis that the BCGs are drawn from a universal distribution, and a cannibalism scenario merging two galaxies from this distribution. The statistical model is known to fail from work as far back as Tremaine&Richstone (1977). We show that neither the statistical model nor the simplest possibility of cannibalism provide a good match for observations. Our CLF models both give similar results, in good agreement with observations. Specifically, we find <m_1> between -25 and -25.5 in the K-band, sigma(m_1)~0.25 and <Delta_12> between 0.6 and 0.8, for cluster luminosities in the range of 10^12 to 10^13 h^-2 L_sun.

We describe high resolution Smoothed Particle Hydrodynamics (SPH) simulations of three approximately $M_*$ field galaxies starting from \LCDM initial conditions. The simulations are made intentionally simple, and include photoionization, cooling of the intergalactic medium, and star formation but not feedback from AGN or supernovae. All of the galaxies undergo an initial burst of star formation at $z \approx 5$, accompanied by the formation of a bubble of heated gas. Two out of three galaxies show early-type properties at present whereas only one of them experienced a major merger. Heating from shocks and -PdV work dominates over cooling so that for most of the gas the temperature is an increasing function of time. By $z \approx 1$ a significant fraction of the final stellar mass is in place and the spectral energy distribution resembles those of observed massive red galaxies. The galaxies have grown from $z=1 \to 0$ on average by 25% in mass and in size by gas poor (dry) stellar mergers. By the present day, the simulated galaxies are old ($\approx 10 {\rm Gyrs}$), kinematically hot stellar systems surrounded by hot gaseous haloes. Stars dominate the mass of the galaxies up to $\approx 4$ effective radii ($\approx 10$ kpc). Kinematic and most photometric properties are in good agreement with those of observed elliptical galaxies. The galaxy with a major merger develops a counter-rotating core. Our simulations show that realistic intermediate mass giant elliptical galaxies with plausible formation histories can be formed from \LCDM initial conditions even without requiring recent major mergers or feedback from supernovae or AGN.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

We measure the three-dimensional topology of large-scale structure in the Sloan Digital Sky Survey (SDSS). This allows the genus statistic to be measured with unprecedented statistical accuracy. The sample size is now sufficiently large to allow the topology to be an important tool for testing galaxy formation models. For comparison, we make mock SDSS samples using several state-of-the-art N-body simulations: the Millennium run of Springel et al. (2005)(10 billion particles), Kim & Park (2006) CDM models (1.1 billion particles), and Cen & Ostriker (2006) hydrodynamic code models (8.6 billion cell hydro mesh). Each of these simulations uses a different method for modeling galaxy formation. The SDSS data show a genus curve that is broadly characteristic of that produced by Gaussian random phase initial conditions. Thus the data strongly support the standard model of inflation where Gaussian random phase initial conditions are produced by random quantum fluctuations in the early universe. But on top of this general shape there are measurable differences produced by non-linear gravitational effects (cf. Matsubara 1994), and biasing connected with galaxy formation. The N-body simulations have been tuned to reproduce the power spectrum and multiplicity function but not topology, so topology is an acid test for these models. The data show a ``meatball'' shift (only partly due to the Sloan Great Wall of Galaxies; this shift also appears in a sub-sample not containing the Wall) which differs at the 2.5\sigma level from the results of the Millennium run and the Kim & Park dark halo models, even including the effects of cosmic variance.

We have identified over 2000 well resolved cluster halos, and also their associated bound subhalos, from the output of a 1024^3 particle cosmological N-body simulation (of box size 320 h^-1 Mpc and softening length 3.2 h^-1 kpc). This has allowed us to measure halo quantities in a statistically meaningful way, and for the first time analyse their distribution for a large and well resolved sample. We characterize each halo in terms of its morphology, concentration, spin, circular velocity and the fraction of their mass in substructure. We also identify those halos that have not yet reached a state of dynamical equilibrium using the virial theorem with an additional correction to account for the surface pressure at the boundary. These amount to 3.4% of our initial sample. For the virialized halos, we find a median of 5.6% of halo mass is contained within substructure, with the distribution ranging between no identified subhalos to 65%. The fraction of mass in substructure increases with halo mass with logarithmic slope of 0.44 +- 0.06. Halos tend to have a prolate morphology, becoming more so with increasing mass. Subhalos have a greater orbital angular momentum per unit mass than their host halo. Furthermore, their orbital angular momentum is typically well aligned with that of their host. Overall, we find that dimensionless properties of dark matter halos do depend on their mass, thereby demonstrating a lack of self-similarity.

We discuss the possibilities for primordial black holes (PBHs) to grow via the accretion of dark matter. In agreement with previous works, we find that accretion during the radiation-dominated era does not lead to a significant mass increase. However, during matter-domination, PBHs may grow by up to two orders of magnitude in mass through the acquisition of large dark matter halos. We discuss the possibility of PBHs being an important component in dark matter halos of galaxies as well as their potential to explain the ultra-luminous x-ray sources (ULXs) observed in nearby galactic disks. We point out that although PBHs are ruled out as the dominant component of dark matter, there is still a great deal of parameter space open to them playing a role in the modern-day universe. For example, a primordial halo population of PBHs each at $10^{2.5} M_\odot$ making up 0.1% of the dark matter grow to $10^{4.5} M_\odot$ via the accumulation of dark matter halos to account for $\sim 10%$ of the dark matter mass by a redshift of $z \approx 30$. These intermediate mass black holes may then ``light up'' when passing through molecular clouds, becoming visible as ULXs at the present day, or they may form the seeds for supermassive black holes at the centers of galaxies.

Taking three independent approaches, we investigate the simultaneous constraints set on the cosmic star formation history from various observations, including stellar mass density and extragalactic background light (EBL). We compare results based on: 1) direct observations of past light-cone, 2) a model using local fossil evidence constrained by SDSS observations at z~0 (the `Fossil' model), and 3) theoretical ab initio models from three calculations of cosmic star formation history: (a) new (1024)^3 Total Variation Diminishing (TVD) cosmological hydrodynamic simulation, (b) analytic expression of Hernquist & Springel based on cosmological Smoothed Particle Hydrodynamics (SPH) simulations, and (c) semi-analytic model of Cole et al. We find good agreement among the three independent approaches up to the order of observational errors, except that all the models predict bolometric EBL of I_tot ~= 37-52 nW m^-2 sr^-1, which is at the lower edge of the the observational estimate by Hauser & Dwek. We emphasize that the Fossil model that consists of two components -- spheroids and disks --, when normalized to the local observations, provides a surprisingly simple but accurate description of the cosmic star formation history and other observable quantities. Our analysis suggests that the consensus global parameters at z=0 are: Omega_* = 0.0023+-0.0004, I_EBL = 43+-7 nW m^-2 sr^-1 rho_SFR=(1.06+-0.22)e-2 Msun yr^-1 Mpc^-3, j_bol = (3.1+-0.2)e8 Lsun Mpc^-3.

We study gravitational lensing with a multiple lens plane approach, proposing a simple analytical model for the probability distribution function (PDF) of the dark matter convergence, kappa, for the different lens planes in a given cosmology as a function of redshift and smoothing angle, theta. The model is fixed solely by the variance of kappa, which in turn is fixed by the amplitude of the power spectrum, sigma_8. We test the PDF against a high resolution Tree-Particle-Mesh simulation and find that it is far superior to the Gaussian or the lognormal, especially for small values of theta << 1 arcmin and at large values of kappa relevant to strong lensing. With this model, we predict the probabilities of strong lensing by a single plane or by multiple planes. We find that for theta ~ 10 arcsec, a single plane accounts for almost all (~ 98%) of the strong lensing cases for source redshift unity. However, for a more typical source redshift of 4, about 12% of the strong lensing cases will result from the contribution of a secondary clump of matter along the line of sight, introducing a systematic error in the determination of the surface density of clusters, typically overestimating it by about 2-5%. We also find that matter inhomogenieties introduce a dispersion in the value of the angular diameter distance about its cosmological mean. The probable error relative to the mean increases with redshift to a value of about 8% for z ~ 6 and theta ~ 10 arcsec.

Cosmological hydrodynamic simulations have significantly improved over the past several years, and we have already shown that the observed properties of Lyman-break galaxies (LBGs) at z=3 can be explained well by the massive galaxies in the simulations. Here we extend our study to z=6 and show that we obtain good agreement for the LBGs at the bright-end of the luminosity function (LF). Our simulations also suggest that the cosmic star formation rate density has a peak at z= 5-6, and that the current LBG surveys at z=6 are missing a significant number of faint galaxies that are dimmer than the current magnitude limit. Together, our results suggest that the universe could be reionized at z=6 by the Pop II stars in ordinary galaxies. We also estimate the LF of Lyman-alpha emitters (LAEs) at z=6 by relating the star formation rate in the simulation to the Ly-alpha luminosity. We find that the simulated LAE LFs agree with the observed data provided that the net escape fraction of Ly-alpha photon is f_{Ly-alpha} <= 0.1. We investigate two possible scenarios for this effect: (1) all sources in the simulation are uniformly dimmer by a factor of 10 through attenuation, and (2) one out of ten LAEs randomly lights up at a given moment. We show that the correlation strength of the LAE spatial distribution can possibly distinguish the two scenarios.

Galactic sized gravitational lenses are simulated by combining a cosmological N-body simulation and models for the baryonic component of the galaxy. The lens caustics, critical curves, image locations and magnification ratios are calculated by ray-shooting on an adaptive grid. When the source is near a cusp in a smooth lens' caustic the sum of the magnifications of the three closest images should be close to zero. It is found that in the observed cases this sum is generally too large to be consistent with the simulations implying that there is not enough substructure in the simulations. This suggests that other factors play an important role. These may include limited numerical resolution, lensing by structure outside the halo, selection bias and the possibility that a randomly selected galaxy halo may be more irregular, for example due to recent mergers, than the isolated halo used in this study. It is also shown that, with the level of substructure computed from the N-body simulations, the image magnifications of the Einstein cross type lenses are very weak functions of source size up to $\sim 1\kpc$. This is also true for the magnification ratios of widely separated images in the fold and cusp caustic lenses. This means that selected magnification ratios for different the emission regions of a lensed quasar should agree with each other, barring microlensing by stars. The source size dependence of the magnification ratio between the closest pair of images is more sensitive to substructure.

Numerical simulations of the intergalactic medium have shown that at the present epoch a significant fraction (40-50%) of the baryonic component should be found in the (T~10^6K) Warm-Hot Intergalactic Medium (WHIM) - with several recent observational lines of evidence indicating the validity of the prediction. We here recompute the evolution of the WHIM with the following major improvements: (1) galactic superwind feedback processes from galaxy/star formation are explicitly included; (2) major metal species (O V to O IX) are computed explicitly in a non-equilibrium way; (3) mass and spatial dynamic ranges are larger by a factor of 8 and 2, respectively, than in our previous simulations. Here are the major findings: (1) galactic superwinds have dramatic effects, increasing the WHIM mass fraction by about 20%, primarily through heating up warm gas near galaxies with density 10^{1.5}-10^4 times the mean density. (2) the fraction of baryons in WHIM is increased modestly from the earlier work but is ~40-50%. (3) the gas density of the WHIM is broadly peaked at a density 10-20 times the mean density, ranging from underdense regions to regions that are overdense by 10^3-10^4. (4) the median metallicity of the WHIM is 0.18 Zsun for oxygen with 50% and 90% intervals being (0.040,0.38) and (0.0017,0.83).

We present a non-parametric, empirically based, model for associating galaxy luminosities with halo/subhalo masses, based on a self-consistent treatment of subhalo mass loss and the subhalo mass function. We find that, at high mass, the mass-luminosity relation is almost independent of the actual luminosity function considered, when luminosity is scaled by the characteristic luminosity L*. Additionally, the shape of the total halo luminosity depends on the slope of the subhalo mass function. For these high mass, cluster sized haloes, we find that total luminosity scales as L_tot ~ M^0.88, while the luminosity of the first brightest galaxy has a much weaker dependence on halo mass, L_1 ~ M^0.28, in good agreement with observations and previous results. At low mass, the resulting slope of the mass-luminosity relation depends strongly of the faint end slope of the luminosity function, and we obtain a steep relation, with approximately L ~ M^4.5 in the K-band. The average number of galaxies per halo/cluster is also in very good agreement with observations, scaling as M^0.9. In general, we obtain a good agreement with several independent sets of observational data. We find that, when comparing with observations and for a flat cosmology, the model tends to prefer lower values for Omega_m and sigma_8. Within the WMAP+SDSS concordance plane of Tegmark et al. (2004), we find best agreement around Omega_m=0.25 and sigma_8=0.8, also in very good agreement with the results of the CMB+2dF study of Sanchez et al. (2005). We also check on possible corrections for observed mass based on a comparison of the equivalent number of haloes/clusters. Additionally, we include further checks on the model results based on the mass to light ratio, the occupation number, the group luminosity function and the multiplicity function. (abridged)

A simple model for the evolution of disc galaxies is presented. We adopt three numbers from observations of the Milky Way disc, the local surface mass density, the stellar scale length (of the assumedly exponential disc) and the amplitude of the (assumedly flat) rotation curve, and physically, the (local) dynamical Kennicutt star formation prescription, standard chemical evolution equations assuming and a model for spectral evolution of stellar populations. We can determine the detailed evolution of the model with only the addition of standard cosmological scalings with time of the dimensional parameters. A surprising wealth of detailed specifications follows from this prescription including the gaseous infall rate as a function of radius and time, the distribution of stellar ages and metallicities with time and radius, surface brightness profiles at different wavelengths, colours etc. At the solar neighbourhood stars start to form $\approx 10 Gyrs$ ago at an increasing rate peaking 4 billion years ago and then slowly declining in good agreement with observations. The mean age of long lived stars at the solar neighbourhood is about $4 Gyrs$. The local surface density of the stars and gas are 35 and $15 M_{\odot}pc^{-2}$, respectively. The metallicity distribution of the stars at the solar radius is narrow with a peak at $[Z/Z_{\odot}] = -0.1$.Both a Salpeter IMF and a Chabrier IMF are consistent with observations. Comparisons with the current and local fossil evidence provides support for the model which can then be used to assess other local disc galaxies, the evolution of disc galaxies in deep optical surveys and also for theoretical investigations such as simulations of merging disc galaxies (abbreviated).

Starting with the well-known NFW dark matter halo distribution, we construct a simple polytropic model for the intracluster medium which is in good agreement with high resolution numerical hydrodynamical simulations, apply this model to a very large scale concordance dark matter simulation, and compare the resulting global properties with recent observations of X-ray clusters, including the mass-temperature and luminosity-temperature relations. We make allowances for a non-negligible surface pressure, removal of low entropy (short cooling time) gas, energy injection due to feedback, and for a relativistic (non-thermal) pressure component. A polytropic index n=5 (Gamma=1.2) provides a good approximation to the internal gas structure of massive clusters (except in the very central regions where cooling becomes important), and allows one to recover the observed M_500-T, L_x-T, and T/n_e^{2/3} ~ T^0.65 relations. Using these concepts and generalizing this method so that it can be applied to fully three-dimensional N-body simulations, one can predict the global X-ray and SZE trends for any specified cosmological model. We find a good fit to observations when assuming that twelve percent of the initial baryonic mass condenses into stars, the fraction of rest mass of this condensed component transferred back to the remaining gas (feedback) is 3.9E-5, and the fraction of total pressure from a nonthermal component is near ten percent.

We present an algorithm which is designed to allow the efficient identification and preliminary dynamical analysis of thousands of structures and substructures in large N-body simulations. First we utilise a refined density gradient system (based on DENMAX) to identify the structures, and then apply an iterative approximate method to identify unbound particles, allowing fast calculation of bound substructures. After producing a catalog of separate energetically bound substructures we check to see which of these are energetically bound to adjacent substructures. For such bound complex subhalos, we combine components and check if additional free particles are also bound to the union, repeating the process iteratively until no further changes are found. Thus our subhalos can contain more than one density maximum, but the scheme is stable: starting with a small smoothing length initially produces small structures which must be combined later, and starting with a large smoothing length produces large structures within which sub-substructure is found. We apply this algorithm to three simulations. Two which are using the TPM algorithm by Bode et al. (2000) and one on a simulated halo by Diemand (2004). For all these halos we find about 5-8% of the mass in substructures.

Recent observations have revealed a population of red massive galaxies at high redshift which are challenging to explain in hierarchical galaxy formation models. We analyze this ``massive galaxy problem'' with two different types of hydrodynamic simulations -- Eulerian TVD and SPH -- of a concordance Lambda CDM universe. We consider two separate but connected aspects of the problem posed by these extremely red objects (EROs): (1) the mass-scale, and (2) red colors. We perform spectrophotometric analyses of simulated galaxies in B, z, R, I, Js, Ks, K filters, and compare their near-IR properties with observations at redshift z=1-3. We find that the simulated galaxies brighter than the magnitude limit of K_vega=20 mag have stellar masses Mstar >= 10^{11} Msun/h and a number density of a few \times 10^{-4} h^3 Mpc^{-3} at z~2, in good agreement with the observed number density in the K20 survey. Therefore, our hydrodynamic simulations do not exhibit the ``mass-scale problem''. The answer to the ``redness problem'' is less clear because of our poor knowledge of the amount of dust extinction in EROs and the uncertain fraction of star-forming EROs. However, our simulations can account for the observed comoving number density of ~ 1e-4 Mpc^{-3} at z=1-2 if we assume a uniform extinction of E(B-V)=0.4 for the entire population of simulated galaxies. The critical question of the existence of a significant population of old, passively evolving massive EROs should continue to be tested with future hydrodynamic simulations. Upcoming observations of the thermal emission of dust in 24 micron by the Spitzer Space Telescope will help to better estimate the dust content of EROs at z=1-3, and thus to further constrain the star formation history of the Universe, and theoretical models of galaxy formation. (abridged)

We present a detailed investigation into which properties of CDM halos make them effective strong gravitational lenses. Strong lensing cross sections of 878 clusters from an N-body simulation are measured by ray tracing through 13,594 unique projections. We measure concentrations, axis ratios, orientations, and the amount of substructure of each cluster, and compare the lensing weighted distribution of each quantity to that of the cluster population as a whole. The concentrations of lensing clusters are on average 34% larger than the typical cluster in the Universe. Despite this bias, the anomalously high concentrations (c >14) recently measured by several groups, appear to be inconsistent with the concentration distribution in our simulations, which predict < 2% of lensing clusters should have concentrations this high. No correlation is found between lensing cross section and the amount of substructure. We introduce several types of simplified dark matter halos, and use them to isolate which properties of CDM clusters make them effective lenses. Projections of halo substructure onto small radii and the large scale mass distribution of clusters do not significantly influence cross sections. The abundance of giant arcs is primarily determined by the mass distribution within an average overdensity of ~ 10,000. A multiple lens plane ray tracing algorithm is used to show that projections of large scale structure increase the giant arc abundance by a modest amount <7%. We revisit the question of whether there is an excess of giant arcs behind high redshift clusters in the RCS survey and find that the number of high redshift (z > 0.6) lenses is in good agreement with LCDM, although our simulations predict more low redshift (z < 0.6) lenses than were observed. (abridged)

We present a simple method for tracing the spatial distribution and predicting the physical properties of the Warm-Hot Intergalactic Medium (WHIM), from the map of galaxy light in the local universe. Under the assumption that biasing is local and monotonic we map the ~ 2 Mpc/h smoothed density field of galaxy light into the mass density field from which we infer the spatial distribution of the WHIM in the local supercluster. Taking into account the scatter in the WHIM density-temperature and density-metallicity relation, extracted from the z=0 outputs of high-resolution and large box size hydro-dynamical cosmological simulations, we are able to quantify the probability of detecting WHIM signatures in the form of absorption features in the X-ray spectra, along arbitrary directions in the sky. To illustrate the usefulness of this semi-analytical method we focus on the WHIM properties in the Virgo Cluster region.

We present predictions of Ly-alpha forest-galaxy correlations at z=3 from Eulerian simulations that include the effects of galactic winds, driven primarily by supernova explosions. Galactic winds produce expanding bubbles of shock-heated gas within ~0.5 comoving Mpc/h of luminous galaxies in the simulation, which have space density similar to that of observed Lyman break galaxies (LBGs). However, most of the low-density intergalactic gas that determines the observed properties of the Ly-alpha forest is unaffected by winds. The impact of winds on the Ly-alpha optical depth near galaxies is less dramatic than their impact on gas temperature because winds heat only a small fraction of the gas present in the turnaround regions surrounding galaxies. Hence, Ly-alpha absorption from gas outside the wind bubbles is spread out over the same velocity range occupied by the wind-heated gas. In general, Ly-alpha absorption is expected to be stronger than average near galaxies because of the high gas density. Winds result in a modest reduction of this expected increase of Ly-alpha absorption. Our predictions can be compared to future observations to detect the wind effects and infer their strength, although with the caveat that the results are still dependent on the correspondence of simulated galaxies and observed LBGs. We find that wind effects in our simulations are not strong enough to reproduce the high Ly-alpha transmission within 0.5 Mpc/h comoving of galaxies that has been suggested by recent observations; powerful galactic explosions or ejecta with hyper-escape velocities would be required, but these are unlikely to be produced by ordinary star formation and supernovae alone

It is now established that several of the Local Group dwarf Spheroidal galaxies (dSphs) have large mass-to-light ratios. We consider the possibility that the dark matter in the halos of dSphs is composed of massive black holes with masses in the range 10^5 to 10^7 solar masses. We use direct N-body simulations to determine the long term evolution of a two-component dSph composed of a pressure-supported stellar population in a black hole halo. The black holes are initially distributed according to a Navarro, Frenk & White profile. For black hole masses between 10^5 and 10^6 solar masses, the dark matter halo evolves towards a shallower inner profile in less than a Hubble time. This suggests that black holes in this mass range might provide an explanation for the origin of the dark matter cores inferred from observations of Low Surface Brightness galaxy rotation curves. We compare the simulated evolution of the stellar population with observed data for the Draco dSph. The dependence of the heating rate on the black hole mass is determined, and an upper limit of 10^5 solar masses is placed on the individual black holes comprising the dark matter halo of Draco, if its present stellar distribution is representative of the initial one. We also present a simple scaling argument which demonstrates that the dynamical heating of an initially compact, though not self-gravitating, stellar distribution might produce a remnant distribution similar in extent to Draco after 10 Gyr, even for black hole masses somewhat in excess of 10^5 solar masses.

We use cosmological simulations to study the X-ray ionisation and heating of the intergalactic medium by an early population of accreting black holes. By considering observational constraints from the X-ray background, we find an upper limit for the optical depth to Thompson scattering tau_e~0.17. The redshifted soft X-ray background from these early sources produces:(i) fully ionised atomic hydrogen in the low density intergalactic medium before redshift z~7 (consequently stellar reionisation is characterised by an instantaneous overlap phase of HII regions),(ii) a second HeII reionisation at z~3 and (iii) heats the intergalactic medium to near 10000 K at low redshifts. The typical luminosity in the soft X-ray band of the galaxies hosting the black holes is about one order of magnitude below the sensitivity limit of the Chandra deep field. About a thousand of these sources may be present per square arcmin of the sky, producing detectable fluctuations. Few rarer objects could be luminous enough to be visible in the Chandra deep field. XEUS and Con-X satellites will be able to detect more of these sources that, if radio loud, could be used to study the 21cm forest in absorption.A signature of an early X-ray preionisation is the production of secondary CMB anisotropies on small angular scales. We find that in these models the power spectrum of temperature fluctuations increases with decreasing angular scale (dT~16 muK at 1arcsec scales), while for stellar reionisation scenarios the power decreases on smaller scales. We also show that the redshifted 21 cm radiation from neutral hydrogen can be marginally detected in emission at redshifts 7<z<12. At a redshift of about z~30 a stronger and narrower (in redshift space) signal in absorption against the CMB, that is peculiar to these models,could be detectable.[abridged]

We study the properties of galaxies at z=2 in a Lambda cold dark matter universe, using two different types of hydrodynamic simulation methods -- Eulerian TVD and smoothed particle hydrodynamics (SPH) -- and a spectrophotometric analysis in the U_n, G, R filter set. The simulated galaxies at z=2 satisfy the color-selection criteria proposed by Adelberger et al. (2004) and Steidel et al. (2004) when we assume Calzetti extinction with E(B-V)=0.15. We find that the number density of simulated galaxies brighter than R<25.5 at z=2 is about 1e-2 h^3 Mpc^-3 for E(B-V)=0.15, which is roughly twice that of the number density found by Erb et al. (2004) for the UV bright sample. This suggests that roughly half of the massive galaxies with M*>10^{10} Msun/h at z=2 are UV bright population, and the other half is bright in the infra-red wavelengths. The most massive galaxies at z=2 have stellar masses >= 10^{11-12} Msun. They typically have been continuously forming stars with a rate exceeding 30 Msun/yr over a few Gyrs from z=10 to z=2, together with significant contribution by starbursts reaching up to 1000 Msun/yr which lie on top of the continuous component. TVD simulations indicate a more sporadic star formation history than the SPH simulations. Our results do not imply that hierarchical galaxy formation fails to account for the observed massive galaxies at z>=1. The global star formation rate density in our simulations peaks at z>=5, a much higher redshift than predicted by the semianalytic models. This star formation history suggests early build-up of the stellar mass density, and predicts that 70 (50, 30)% of the total stellar mass at z=0 had already been formed by z=1 (2, 3). Upcoming observations by Spitzer and Swift might help to better constrain the star formation history at high redshift.

We discuss the importance of feedback via photoionization and Compton heating on the co-evolution of massive black holes (MBHs) at the center of spheroidal galaxies and their stellar and gaseous components. We first assess the energetics of the radiative feedback from a typical quasar on the ambient interstellar medium. We then demonstrate that the observed Mbh-sigma relation could be established following the conversion of most of the gas of an elliptical progenitor into stars, specifically when the gas-to-stars mass ratio in the central regions has dropped to a low level ~0.01 or less, so that gas cooling is no longer able to keep up with the radiative heating by the growing central MBH. A considerable amount of the remaining gas will be expelled and both MBH accretion and star formation will proceed at significantly reduced rates thereafter, in agreement with observations of present day ellipticals. We find further support for this scenario by evolving over an equivalent Hubble time a simple, physically based toy model that additionally takes into account the mass and energy return for the spheroid evolving stellar population, a physical ingredient often neglected in similar approaches.

We discuss the role of feedback via photoionization and Compton heating in the co-evolution of massive black holes at the center of spheroidal galaxies and their stellar and gaseous components. We first assess the energetics of the radiative feedback from a typical quasar on the ambient interstellar gas. We then demonstrate that the observed Mbh-sigma relation could be established at a relatively early epoch in galactic evolution when the formation of the stellar bulge was almost completed and the gas-to-stars mass ratio was reduced to a low level ~0.01 such that cooling could not keep up with radiative heating. A considerable amount of gas was expelled at that time and black hole accretion proceeded at a much lower rate thereafter.

There is abundant evidence that heating processes in the central regions of elliptical galaxies has both prevented large-scale cooling flows and assisted in the expulsion of metal rich gas. We now know that each such spheroidal system harbors in its core a massive black hole weighing approximately 0.13% of the mass in stars and also know that energy was emitted by each of these black holes with an efficiency exceeding 10% of its rest mass. Since, if only 0.5% of that radiant energy were intercepted by the ambient gas, its thermal state would be drastically altered, it is worth examining in detail the interaction between the out-flowing radiation and the equilibrium or inflowing gas. On the basis of detailed hydrodynamic computations we find that relaxation oscillations are to be expected with the radiative feedback quite capable of regulating both the growth of the central black hole and also the density and thermal state of the gas in the galaxy. Mechanical input of energy by jets may assist or dominate over these radiative effects. We propose specific observational tests to identify systems which have experienced strong bursts of radiative heating from their central black holes.

We present a summary of the results obtained with a time-dependent, one-zone toy model aimed at exploring the importance of radiative feedback on the co-evolution of massive black holes (MBHs) at the center of stellar spheroids and their stellar and gaseous components. We consider cosmological infall of gas as well as the mass and energy return for the evolving stellar population. The AGN radiative heating and cooling are described by assuming photoionization equilibrium of a plasma interacting with the average quasar SED. Our results nicely support a new scenario in which the AGN accretion phase characterized by a very short duty-cycle (and now common in the Universe) is due to radiative feedback. The establishment of this phase is recorded as a fossil in the Magorrian and Mbh-sigma relations.

Observational cosmology has indeed made very rapid progress in recent years. The ability to quantify the universe has largely improved due to observational constraints coming from structure formation. The transition to precision cosmology has been spearheaded by measurements of the anisotropy in the cosmic microwave background (CMB) over the past decade. Observations of the large scale structure in the distribution of galaxies, high red-shift supernova, have provided the required complementary information. We review the current status of cosmological parameter estimates from joint analysis of CMB anisotropy and large scale structure (LSS) data. We also sound a note of caution on overstating the successes achieved thus far.

We study the properties of galaxies at z=2 in a Lambda CDM universe, using two different types of hydrodynamic simulation methods (Eulerian TVD and SPH) and a spectrophotometric analysis in the Un, G, R filter set. The simulated galaxies at z=2 satisfy the color-selection criteria proposed by Adelberger et al. (2004) when we assume Calzetti extinction with E(B-V)=0.15. We find that the number density of simulated galaxies brighter than R<25.5 at z=2 is about 2e-2 h^3/Mpc^3, roughly one order of magnitude larger than that of Lyman break galaxies at z=3. The most massive galaxies at z=2 have stellar masses >~1e11 Msun, and their observed-frame G-R colors lie in the range 0.0<G-R<1.0. They typically have been continuously forming stars with a rate exceeding 30 Msun/yr over a few Gyrs from z=10 to z=2, although the TVD simulation indicates a more sporadic star formation history than the SPH simulations. Of order half of their stellar mass was already assembled by z~4. The reddest massive galaxies at z=2 with G-R >= 1.0 and Mstar>1e10 Msun/h finished the build-up of their stellar mass by z~3. Interestingly, our study suggests that the majority of the most massive galaxies at z=2 should be detectable at rest-frame UV wavelengths, contrary to some recent claims made on the basis of near-IR studies of galaxies at the same epoch, provided the median extinction is less than E(B-V)<0.3. However, our results also suggest that the fraction of stellar mass contained in galaxies that pass the color-selection criteria could be as low as 50% of the total stellar mass in the Universe at z=2. Our simulations suggest that the missing stellar mass is contained in fainter (R>25.5) and intrinsically redder galaxies. Our results do not suggest that hierarchical galaxy formation fails to account for the massive galaxies at z>=1. (abridged)

We explore several physical effects on the power spectrum of the Lyman-alpha forest transmitted flux. The effects we investigate here are usually not part of hydrodynamic simulations and so need to be estimated separately. The most important effect is that of high column density absorbers with damping wings, which add power on large scales. We compute their effect using the observational constraints on their abundance as a function of column density. Ignoring their effect leads to an underestimation of the slope of the linear theory power spectrum. The second effect we investigate is that of fluctuations in the ionizing radiation field. For this purpose we use a very large high resolution N-body simulation, which allows us to simulate both the fluctuations in the ionizing radiation and the small scale LyaF within the same simulation. We find an enhancement of power on large scales for quasars and a suppression for galaxies. The strength of the effect rapidly increases with increasing redshift, allowing it to be uniquely identified in cases where it is significant. We develop templates which can be used to search for this effect as a function of quasar lifetime, quasar luminosity function, and attenuation length. Finally, we explore the effects of galactic winds using hydrodynamic simulations. We find the wind effects on the LyaF power spectrum to be be degenerate with parameters related to the temperature of the gas that are already marginalized over in cosmological fits. While more work is needed to conclusively exclude all possible systematic errors, our results suggest that, in the context of data analysis procedures where parameters of the LyaF model are properly marginalized over, the flux power spectrum is a reliable tracer of cosmological information.

We combine the constraints from the recent Ly-alpha forest and bias analysis of the SDSS with previous constraints from SDSS galaxy clustering, the latest supernovae, and WMAP . Combining WMAP and the lya we find for the primordial slope n_s=0.98\pm 0.02. We see no evidence of running, dn/d\ln k=-0.003\pm 0.010, a factor of 3 improvement over previous constraints. We also find no evidence of tensors, r<0.36 (95% c.l.). A positive correlation between tensors and primordial slope disfavors chaotic inflation type models with steep slopes: V \propto \phi^4 is outside the 3-sigma contour. For the amplitude we find sigma_8=0.90\pm 0.03 from the lyaf and WMAP alone. We find no evidence of neutrino mass: for the case of 3 massive neutrino families with an inflationary prior, \sum m_{\nu}<0.42eV and the mass of lightest neutrino is m_1<0.13eV at 95% c.l. For the 3 massless + 1 massive neutrino case we find m_{\nu}<0.79eV for the massive neutrino, excluding at 95% c.l. all neutrino mass solutions compatible with the LSND results. We explore dark energy constraints in models with a fairly general time dependence of dark energy equation of state, finding Omega_lambda=0.72\pm 0.02, w(z=0.3)=-0.98^{+0.10}_{-0.12}, the latter changing to w(z=0.3)=-0.92^{+0.09}_{-0.10} if tensors are allowed. We find no evidence for variation of the equation of state with redshift, w(z=1)=-1.03^{+0.21}_{-0.28}. These results rely on the current understanding of the lyaf and other probes, which need to be explored further both observationally and theoretically, but extensive tests reveal no evidence of inconsistency among different data sets used here.

We investigate the possibility that there is a first phase of partial ionisation due to X-rays produced by black hole accretion in small-mass galaxies at redshifts 7<z<20. This is followed by complete reionisation by stellar sources at z~7. This scenario is motivated by the large optical depth to Thompson scattering, tau_e=0.17, measured by WMAP. But it is also consistent with the observed Gunn-Peterson trough in the spectra of quasars at z~5-6. We use a semianalytic code to explore models with different black hole accretion histories and cosmological parameters.We find that ``preionisation'' by X-rays can increase the intergalactic medium (IGM) optical depth from tau_e~0.06 given by stellar sources only, to 0.1<tau_e<0.2, if a fraction of baryons 5x10^-5 is accreted onto seed black holes produced in the collapse of low metallicity, high mass stars before z~15. To be effective, preionisation requires a non-negligible star formation in the first small-mass galaxies in hich seed black holes are formed.The model predicts that dwarf spheroidal galaxies may host a mass in black holes that is 5-40% of their stellar mass. The redshifted X-ray background produced by this early epoch of black hole accretion constitutes about 5-10% of the X-ray background in the 2-50 keV bands and roughtly half of the currently estimated black hole mass density was formed at early times.Moreover, in most models, the photons from the redshifted background are sufficient to fully reionise HeII at redshift z~3 without any additional contribution from quasars at lower redshifts and the temperature of the mean density intergalactic medium remains close to 10^4 K down to redshift z~1.[abridged]

We investigate possible signatures of feedback from galactic superwinds on the metallicity of the Lyman alpha forest, using a set of high resolution hydrodynamic simulations of a Lambda cold dark matter model. Simulations produce metals self-consistently, based on one single parameter, the metal yield, which in turn is constrained by metallicity in the intra-cluster gas. We follow metals as a separate density species. For the metallicity of Lyman alpha clouds with column density of N_{HI}~10^{14.5}-10^{15.5}/cm^2 at z=2-4 we find reasonable agreement between simulations, both with and without GSW, and observations (Schaye et al). A unique signature and sensitive test of GSW is, however, provided by lower density regions with gas density of \rho/<\rho> = 0.01-1.0 and a corresponding column density of 10^{12}-10^{14}/cm^2. Without GSW we predict that both the mean and median metallicity of \lya clouds in this column density range at z=2-4 should have Z<10^{3}Zsun, since these small systems support little star formation. GSW contaminate these regions, however, and also there is a significant fraction (~25%) of Lyman alpha clouds in this column density range which have a high metallicity excess of 10^{-2}Zsun, resulting in a mean metallicity of ~10^{-2}Zsun. In addition, we find that there is a minimum in the median metallicity for clouds of N_{HI}~10^{13}-10^{14}/cm^2 in the case with GSW, whereas without GSW the metallicity decreases monotonically and rapidly with decreasing column density. Finally, we predict that the ratio of secondary (e.g., N) to primary metals (e.g., O,C) is expected to be smaller by a factor of 10 in clouds of N_{HI}~10^{14.5}/cm^2 compared to that in large galaxies; this factor increases to >50 for N_{HI}< 10^{13.5}/cm^2.

In this paper we present a new, essentially empirical, model for the relation between the mass of a dark matter halo/subhalo and the luminosity of a galaxy hosted in it. To estimate this, we replace the assumption of linearity between light and mass fluctuations with the assumption of monotonicity between galaxy light and halo or subhalo mass. We are enabled to proceed with this less restrictive ansatz by the availability of new, very high resolution dark matter simulations and more detailed and comprehensive global galactic luminosity functions. We find that the relation between halo/subhalo mass and hosted galaxy luminosity is fairly well fit by a double power law, while that between halo mass and group luminosity is fairly well fit by a two branch function, with both branches double power laws. Both relations asymptote to $L\propto M^4$ at low $M$, while at high mass the former follows $L\propto M^{0.28}$ and the latter $L\propto M^{0.9}$. We also derive results for the occupation number, luminosity function of cluster galaxies, group luminosity function and multiplicity function. Then, using some additional assumptions, we further derive results for biasing between mass and light and mass and galaxy number, light distribution function and the void probability distribution. Our results for the most part seem to match well with observations and previous expectations. We feel this is a potentially powerful way of modelling this relation, since the main inputs are readily testable against dark matter simulation results and galaxy surveys. (abridged)

We study the evolution of the global stellar mass density in a Lambda cold dark matter universe using two different types of hydrodynamical simulations (Eulerian TVD and SPH) and the analytical model of Hernquist & Springel (2003). We find that the theoretical calculations all predict both a higher stellar mass density at z~3 than indicated by current observations, and that the peak of the cosmic star formation rate history should lie at z~5. Such a star formation history implies that as much as (70%, 30%) of the total stellar mass density today must already have formed by z=(1, 3). Our results suggest that current observations at z~3 are missing as much as 50% of the total stellar mass density in the Universe, perhaps owing to an inadequate allowance for dust obscuration in star-forming galaxies, limited sample sizes, or cosmic variance. We also compare our results with some of the updated semi-analytic models of galaxy formation.

To date, in almost all strong gravitational lensing analyses for modeling giant arc systems and multiple quasar images, it has been assumed that all the deflecting matter is concentrated in one lens plane at a certain distance - the thin lens approximation. However, in a few observed cases, lenses at more than one redshift have been identified as contributing to the image splitting. Here we report on a quantitative investigation of the importance and frequency of significant multiple lensing agents. We use multi-lens plane simulations to evaluate how frequently two or more lens planes combined are essential for multiple imaging, as compared with the cases where a single lens plane alone provides enough focusing to be supercritical. We find that the fraction of cases for which more than one lens plane contributes significantly to a multi-image lensing situation is a strong function of source redshift. For sources at redshift unity, 95% of lenses involve only a single mass concentration, but for a more typical scenario with, e.g., a source at a redshift of z_s = 3.8, as many as 38% of the strongly lensed quasars/arcs occur because of a significant matter contribution from one or more ADDITIONAL lens planes. In the 30% to 40% of cases when additional planes make a significant contribution, the surface mass density of the primary lens will be overestimated by about 15% to 20%, if the additional contributions are not recognized.

We compute the properties of dark matter halos with mass $10^{6.5}-10^9\msun$ at redshift $z=6-11$ in the standard cold dark matter cosmological model, utilizing a very high resolution N-body simulation. We find that dark matter halos in these mass and redshift ranges are significantly biased over matter with a bias factor in the range 2-6. The dark matter halo mass function displays a slope of $2.05\pm 0.15$ at the small mass end. We do not find a universal dark matter density profile. Instead, we find a significant dependence of the central density profile of dark matter halos on halo mass and epoch with $\alpha_0=0.4-1.0$; the high-mass ($M\ge 10^8\msun$) low-redshift ($z\sim 6$) halos occupy the high end of the range and low-mass ($M\sim 10^{7}\msun$) high-redshift ($z\sim 11$) halos occupy the low end. Additionally, for fixed mass and epoch there is a significant dispersion in $\alpha_0$ due to the stochastic assembly of halos. Our results fit a relationship of the form $\alpha_0=0.75((1+z)/7.0)^{-1.25}(M/10^7\msun)^{0.11(1+z)/7.0}$ with a dispersion about this fit of $\pm 0.5$ and no systematic dependence of variance correlated with environment. The median spin parameter of dark matter halos is $0.03-0.04$ but with a large lognormal dispersion of $\sim 0.4$. Various quantities are tabulated or fitted with empirical formulae.

We investigate the effects of a top-heavy stellar initial mass function on the reionisation history of the intergalactic medium (IGM). We use cosmological simulations that include self-consistently the feedback from ionising radiation, H_2 dissociating radiation and supernova (SN) explosions. We find that it is difficult to reionise the IGM at z_rei>10 with stellar sources even after making extreme assumptions. If star formation in 10^9 M_\odot galaxies is not suppressed by SN explosions, the optical depth to Thomson scattering is tau_e< 0.13. If we allow for the normal energy input from SNe or if pair-instability SNe are dominant, we find taue<0.09. Assuming normal yields for the first stars (popIII), the mean metallicity of the IGM is already Z/Z_\odot=2x10^-3 (10^-3<Z/Z_\odot<1 in overdense regions) when the IGM mean ionisation fraction is less than 10%. For these reasons popIII stars cannot contribute significantly to reionisation unless the mechanical energy input from SNe is greatly reduced and either the metal yield or the mixing efficiency is reduced by a factor of 1000. Both problems have a solution if popIII stars collapse to black holes. This can happen if, having masses M_*< 130 M_\odot, they are characterised by heavy element fall-back or if, having masses M_*>260 M_\odot, they collapse directly onto black holes without exploding as SNe.If metal-poor stars are initially important and collapse to black holes is the typical outcome, then the secondary emission of ionising radiation from accretion on SN induced seed black holes, might be more important than the primary emission. Surprisingly, including feedback effects, we estimate that a warm dark matter scenario (with particle mass of 1.25 keV) reduces tau_e by only approximately 10%.[abridged]

We review the evidence for substructures from the anomalous flux ratios in gravitational lenses. Using high-resolution numerical simulations, we show that at typical image positions, the fraction of surface mass density in substructures is <~ 0.5% with mass above 10^{-4} virial masses in the ``concordance'' LCDM cosmology. Substructures outside the virial radius (but projected at typical lens image positions) only increase the fraction moderately. Several effects, in particular baryonic settling and the requirement of compactness, may further decrease the predictions by a factor of few. The predicted fraction with appropriate properties thus appears to be lower than that required by lensing, although both are still uncertain. More speculative substructures such as massive black holes (M ~ 10^5-10^6 solar masses) in the halo may offer viable alternatives.

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.

In standard treatments the Eddington luminosity is calculated by assuming that the electron-photon cross section is well described by the Thomson cross section which is gray (frequency independent). Here we discuss some consequence of the introduction of free-free opacity in the Eddington luminosity computation: in particular, due to the dependence of free-free emission on the square of the gas density, it follows that the associated absorption cross section increases linearly with the gas density, so that in high density environments Eddington luminosity is correspondingly reduced. We present a summary of an ongoing exploration of the parameter space of the problem, and we conclude that Eddington luminosity in high density environments can be lowered by a factor of ten or more, making it considerably easier for black holes to accelerate and eject ambient gas.

Recent observations by Chandra and XMM-Newton demonstrate that the central gas in "cooling flow" galaxy clusters has a mass cooling rate that decreases rapidly with decreasing temperature. This contrasts the predictions of a steady state cooling flow model. On the basis of these observational results, the gas can be in a steady state only if a steady temperature dependent heating mechanism is present; alternatively the gas could be in an unsteady state, i.e., heated intermittently. Intermittent heating can be produced by accretion on the supermassive black hole residing in the central cluster galaxy, via Compton heating. This mechanism can be effective provided that the radiation temperature of the emitted spectrum is higher than the gas temperature. Here we explore whether this heating mechanism can be at the origin of the enigmatic behavior of the hot gas in the central regions of ``cooling flow'' clusters. Although several characteristics of Compton heating appear attractive in this respect, we find that the fraction of absorbed heating for realistic gas and radiation temperatures falls short by two orders of magnitude of the required heating.

We use a large-scale Lambda-CDM hydrodynamical simulation to assess the dependence of the cosmic Star Formation Rate (SFR) on overdensity of luminosity. The `morphology-density relation' is qualitatively captured when the SFR is plotted as a function of redshift, divided into the quartiles of the overdensity distribution.

We measure the large-scale real-space power spectrum P(k) using a sample of 205,443 galaxies from the Sloan Digital Sky Survey, covering 2417 square degrees with mean redshift z~0.1. We employ a matrix-based method using pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 22 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.02 h/Mpc < k < 0.3h/Mpc. We pay particular attention to modeling, quantifying and correcting for potential systematic errors, nonlinear redshift distortions and the artificial red-tilt caused by luminosity-dependent bias. Our final result is a measurement of the real-space matter power spectrum P(k) up to an unknown overall multiplicative bias factor. Our calculations suggest that this bias factor is independent of scale to better than a few percent for k<0.1h/Mpc, thereby making our results useful for precision measurements of cosmological parameters in conjunction with data from other experiments such as the WMAP satellite. As a simple characterization of the data, our measurements are well fit by a flat scale-invariant adiabatic cosmological model with h Omega_m =0.201+/- 0.017 and L* galaxy sigma_8=0.89 +/- 0.02 when fixing the baryon fraction Omega_b/Omega_m=0.17 and the Hubble parameter h=0.72; cosmological interpretation is given in a companion paper.

We present the results of a series of gas dynamical cosmological simulations of the formation of individual massive field galaxies in the standard concordance \LCDM and in a \LWDM cosmology. The galaxies form in an initial burst of star formation followed by accretion of small satellites. The simulated galaxies are old ($\approx$ 10 $Gyrs$) hot stellar systems with masses of $\approx 1.7 \times 10^{11} M_{\odot}$. Baryonic matter dominates the mass in the luminous part of the galaxies up to $\approx 5$ effective radii. The \LCDM galaxy is a slowly rotating spheroidal stellar system (E2) with predominantly disky isophotes. The line-of-sight velocity distributions deviate from Gaussian shape and $h_3$ is anticorrelated with $v_{los}$. The corresponding \LWDM galaxy is more elongated (E3 - E4) and rotates faster. The anisotropy parameter $(v/\sigma)^*$ is close to unity. Within the uncertainties of $M/L$ both galaxies follow the Fundamental Plane. We conclude that the properties of the two galaxies simulated in the \LCDM and \LWDM cosmology are in good agreement with observations of intermediate mass elliptical or S0 galaxies (abstract abbreviated).

Using information on the cosmic X-ray background and the cumulative light of active galactic nuclei at infrared wavelengths, the estimated local mass density of galactic massive black holes (MBHs) and published AGN composite spectra in the optical, UV and X-ray, we compute the characteristic angular-integrated, broad-band spectral energy distribution of the average quasar in the universe. We demonstrate that the radiation from such sources can photoionize and Compton heat the plasma surrounding them up to an equilibrium Compton temperature (Tc) of 2x10^7 K. It is shown that circumnuclear obscuration cannot significantly affect the net gas Compton heating and cooling rates, so that the above Tc value is approximately characteristic of both obscured and unobscured quasars. This temperature is above typical gas temperatures in elliptical galaxies and just above the virial temperatures of giant ellipticals. The general results of this work can be used for accurate calculations of the feedback effect of MBHs on both their immediate environs and the more distant interstellar medium of their host galaxies.

Based on observed rotation curves of galaxies and theoretical simulations of dark matter halos, there are reasons for believing that at least three different types of dark matter halos exist in the Universe classified by their masses M and the inner slope of mass density -\alpha: Population A (galaxies): 10^{10} h^{-1} M_\odot < M < 2 \times 10^{13} h^{-1} M_\odot, \alpha = 2; Population B (cluster halos): M > 2 \times 10^{13} h^{-1} M_\odot, \alpha = 1.3; and Population C (dwarf halos): M < 10^{10} h^{-1} M_\odot, \alpha = 1.3. In this paper we calculate the lensing probability produced by such a compound population of dark halos, for both image separation and time delay, assuming that the mass function of halos is given by the Press-Schechter function and the Universe is described by an LCDM, OCDM, or SCDM model. The LCDM model is normalized to the WMAP observations, OCDM and SCDM models are normalized to the abundance of rich clusters. We compare the predictions of the different cosmological models with observational data and show that, both LCDM and OCDM models are marginally consistent with the current available data, but the SCDM model is ruled out. The fit of the compound model to the observed correlation between splitting angle and time delay is excellent but the fit to the number vs splitting angle relation is only adequate using the small number of sources in the objective JVAS/CLASS survey. A larger survey of the same type would have great power in discriminating among cosmological models. Furthermore, population C in an LCDM model has a unique signature in the time domain, an additional peak at ~3 seconds potentially observable in GRBs, which makes it distinguishable from variants of CDM scenarios, such as warm dark matter, repulsive dark matter, or collisional dark matter.

Utilizing a new, high mass resolution hydrodynamic simulation we compute the metallicity evolution of damped Lyman alpha systems (DLAs) and find a reasonable agreement with observations. In particular, the observed slow evolution of the DLA metallicity occurs naturally in the simulation due to the combined effects of physical and observational selection. The slow metallicity evolution is caused by the steady transformation, with increasing time, of the highest metallicity systems to "galaxies", thus depleting this category, while all the lower metallicity systems show, individually, an increase in metallicity. Although the trend of DLA metallicity with redshift is in good agreement with observations, it appears that the average metallicity of simulated DLAs is higher than observed by 0.3-0.5 dex in the probed redshift range (z=0-5). Our study indicates that this difference may be attributed to observational selection effects due to dust obscuration. If we allow for a dust obscuration effect, our model reproduces the observed metallicity evolution in both amplitude and slope. We find that DLAs are not a simple population but probe a range of different systems and the mix changes with redshift. About 50% of all metals in the gaseous phase is in DLAs at all times from z=5 to z=1, making a rapid downturn at z\le 1 to ~20% by z=0, as metals are swept into the hotter components of the IGM as well as locked up in stars. While not the primary focus of this study, we find that the model provides good matches to observations with respect to column density distribution and evolution of neutral gas content, if the same dust obscuration is taken into account. We find Omega_{DLA,comp}=(1-3)E-3, depending on the effect of dust obscuration.

While second and higher order correlations of the light distribution have received extensive study, the lowest order probability distribution function (PDF) -- the probability that a unit volume of space will emit a given amount of light -- has received very little attention. We estimate this function with the aid of hydrodynamic simulations of the Lambda CDM model, finding it significantly different from the mass density PDF, and not simply related to it by linear bias or any of the other prescriptions commonly adopted. If the optical light PDF is, in reality, similar to what we find in the simulations, then some measures of Omega_M based on mass-to-light ratio and the cosmic virial theorem will have significantly underestimated Omega_M. Basically, the problem is one of selection bias, with galaxy forming regions being unrepresentative of the dark matter distribution in a way not described by linear bias. Knowledge of the optical PDF and the plausible assumption of a log-normal distribution for the matter PDF will allow one to correct for these selection biases. We find that this correction (which amounts to 20-30%) brings the values of Omega_M estimated by using the mass-to-light ratio and the cosmic virial theorem to the range 0.2-0.3, in better agreement with the WMAP result than the uncorrected estimates. In addition, the relation between mass and light PDFs gives us insight concerning the nature of the void phenomenon. In particular our simulation indicates that 20% of mass is distributed in voids, which occupy 85% of volume in the universe.

Dark matter, proposed decades ago as a speculative component of the universe, is now known to be the vital ingredient in the cosmos, eight times more abundant than ordinary matter, one quarter of the total energy density and the component which has controlled the growth of structure in the universe. Its nature remains a mystery, but, assuming it is comprised of weakly interacting sub-atomic particles, is consistent with large scale cosmic structure. However, recent analyses of structure on galactic and sub-galactic scales have suggested discrepancies and stimulated numerous alternative proposals. We discuss how studies of the density, demography, history and environment of smaller scale structures may distinguish among these possibilities and shed new light on the nature of dark matter.

The frequency of giant arcs - highly distorted and strongly gravitationally lensed background galaxies - is a powerful test for cosmological models. Previous comparisons of arc statistics for the currently favored concordance cosmological model (a flat lambda-dominated universe) with observations have shown an apparently large discrepancy. We present here new ray-shooting results, based on a high resolution (1024^3 particles in a 320 Mpc/h box) large-scale structure simulation normalized to the WMAP observations. We follow light rays through a pseudo-3D matter distribution approximated by up to 38 lens planes, and evaluate the occurrence of arcs for various source redshifts. We find that the frequency of strongly lensed background sources is a steep function of source redshift: the optical depth for giant arcs increases by a factor of five when the background sources are moved from redshift z_s = 1.0 to z_s = 1.5. This is a consequence of a moderate decrease of the critical surface mass density for lensing, combined with the very steep cluster mass function at the high mass end. Our results are consistent with those of Bartelmann et al. (1998) if we - as they did - restrict all sources to be exactly at z_s = 1. But if we allow for a more realistic distribution of source redshifts extending to or beyond z_s > 1.5, the apparent discrepancy vanishes: the frequency of arcs is increased by about a factor of ten as compared to previous estimates, and results in roughly one arc per 20 square degrees over the sky. This prediction for an LCDM model is then in good agreement with the observed frequency of arcs. Hence we consider the ``missing arc'' problem for a concordance LCDM cosmology to be solved.

We present constraints on cosmological and star formation parameters based on combining observations of the Wilkinson Microwave Anisotropy Probe (WMAP) and high-redshift quasars from the Sloan Digital Sky Survey (SDSS). We use a semi-analytic model for reionization (Chiu and Ostriker 2000) that takes into account a number of important physical processes both within collapsing halos and in the intergalactic medium. Assuming that the efficiency of producing UV photons per baryon is constant, we derive a constraint of the form sigma_8 Omega_0^0.5~0.33 in a flat, Lambda-dominated universe with h=0.72, n=0.99, and Omega_b h^2=0.024. However, the calculated optical depth to electron scattering of tau_es~0.06 is well below the value found by WMAP of 0.17+/-(0.04~0.07) (Spergel et al 2003). Since the WMAP constraints on tau_es are somewhat degenerate with the value of the spectral index n, we then permit the primordial spectral index n to float and fix Omega_0 h^2=0.14, while normalizing the power spectrum using WMAP. In addition, we allow the UV-efficiency to have time-dependence. Combining the WMAP constraints with the quasar transmission data, our analysis then favors a model with tau_es=0.11^{+0.02}_{-0.03}, n=0.96^{+0.02}_{-0.03}$, implying sigma_8=0.83^{+0.03}_{-0.05} (95% confidence), and an effective UV-efficiency that was at least ~10x greater at z >> 6. These results indicate that the quasar and WMAP observations are consistent. If future observations confirm an optical depth to electron scattering tau_es~0.1, then it would appear that no more "exotic" sources of UV-photons, such as mini-quasars or AGNs, are necessary; but our analysis indicates that a determination of tau_es>~0.17 would require a more radical solution.

We discuss a simple model for the growth of supermassive black holes (BHs) at the center of spheroidal stellar systems. In particular, we assess the hypotheses that (1) star formation in spheroids and BH fueling are proportional to one another, and (2) the BH accretion luminosity stays near the Eddington limit during luminous quasar phases. With the aid of this simple model, we are able to interpret many properties of the QSO luminosity function, including the puzzling steep decline of the characteristic luminosity from redshift z=2 to to z=0: indeed the residual star formation in spheroidal systems is today limited to a small number of bulges, characterized by stellar velocity dispersions a factor of 2-3 smaller those of the elliptical galaxies hosting QSOs at z > 2. A simple consequence of our hypotheses is that the redshift evolution of the QSO emissivity and of the star formation history in spheroids should be roughly parallel. We find this result to be broadly consistent with our knowledge of the evolution of both the global star formation rate, and of the evolution of the QSO emissivity, but we identify interesting discrepancies at both low and high redshifts, to which we offer tentative solutions. Finally, our hypotheses allow us to present a robust method to derive the duty cycle of QSO activity, based on the observed QSO luminosity function, and on the present-day relation between the masses of supermassive BHs and those of their spheroidal host stellar systems. The duty cycle is found to be substantially less than unity, with characteristic values in the range (3-6)x10^(-3), and we compute that the average bolometric radiative efficiency is epsilon=0.07. Finally, we find that the growth in mass of individual black holes at high redshift (z>2) can be dominated by mergers, and is therefore not necessarily limited by accretion.

The recent announcement from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite experiment combined with other recent advances in observational cosmology verifies key components of the standard cosmological model. However, we argue that there remain some significant open issues regarding the basic history and composition of the Universe and uncertainties in some of the most important parameters.

An improved implementation of an N-body code for simulating collisionless cosmological dynamics is presented. TPM (Tree-Particle-Mesh) combines the PM method on large scales with a tree code to handle particle-particle interactions at small separations. After the global PM forces are calculated, spatially distinct regions above a given density contrast are located; the tree code calculates the gravitational interactions inside these denser objects at higher spatial and temporal resolution. The new implementation includes individual particle time steps within trees, an improved treatment of tidal forces on trees, new criteria for higher force resolution and choice of time step, and parallel treatment of large trees. TPM is compared to P^3M and a tree code (GADGET) and is found to give equivalent results in significantly less time. The implementation is highly portable (requiring a Fortran compiler and MPI) and efficient on parallel machines. The source code can be found at http://astro.princeton.edu/~bode/TPM/

We propose a possible way to detect baryons at low redshifts from the analysis of X-ray absorption spectra of bright AGN pairs. A simple semi-analytical model to simulate the spectra is presented. We model the diffuse warm-hot intergalactic medium (WHIM) component, responsible for the X-ray absorption, using inputs from high-resolution hydro-dynamical simulations and analytical prescriptions. We show that the number of OVII absorbers per unit redshift with column density larger than $10^{13.5}$ cm$^{-2}$ - corresponding to an equivalent width of $\sim$ 1 km/s - which will be possibly detectable by {\it XEUS}, is $\magcir 30$ per unit redshift. {\it Constellation-X} will detect $\sim 6$ OVII absorptions per unit redshift with an equivalent width of 10 km/s. Our results show that, in a $\Lambda$CDM Universe, the characteristic size of these absorbers at $z\sim 0.1$ is $\sim 1$ $h^{-1}$ Mpc. The filamentary structure of WHIM can be probed by finding coincident absorption lines in the spectra of background AGN pairs. We estimate that at least 20 AGN pairs at separation $\mincir 20$ arcmin are needed to detect this filamentary structure at a 3$\sigma$ level. Assuming observations of distant sources using {\it XEUS} for exposure times of 500 ksec, we find that the minimum source flux to probe the filamentary structure is $\sim 2\times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$, in the 0.1-2.4 keV energy band. Thus, most pairs of these extragalactic X-ray bright sources have already been identified in the {\it ROSAT} All-Sky Survey. Re-observation of these objects by future missions could be a powerful way to search for baryons in the low redshift Universe.

The mass function of clusters of galaxies is determined from 400 deg^2 of early commissioning imaging data of the Sloan Digital Sky Survey; ~300 clusters in the redshift range z = 0.1 - 0.2 are used. Clusters are selected using two independent selection methods: a Matched Filter and a red-sequence color magnitude technique. The two methods yield consistent results. The cluster mass function is compared with large-scale cosmological simulations. We find a best-fit cluster normalization relation of sigma_8*omega_m^0.6 = 0.33 +- 0.03 (for 0.1 ~< omega_m ~< 0.4), or equivalently sigma_8 = (0.16/omega_m)^0.6. The amplitude of this relation is significantly lower than the previous canonical value, implying that either omega_m is lower than previously expected (omega_m = 0.16 if sigma_8 = 1) or sigma_8 is lower than expected (sigma_8 = 0.7 if omega_m = 0.3). The best-fit mass function parameters are omega_m = 0.19 (+0.08,-0.07) and sigma_8 = 0.9 (+0.3,-0.2). High values of omega_m (>= 0.4) and low sigma_8 (=< 0.6) are excluded at >~ 2 sigma.

We consider the consequences of SIDM for a velocity dependent cross section per unit mass. Accretion of SIDM onto seed black holes can produce supermassive black holes that are too large for certain combinations of parameters,which is used to obtain a new constraint on the dark matter interaction. Constraints due to other considerations are presented and previous ones are generalized. The black hole constraint is extremely sensitive to the slope \alpha, of the inner density profile of dark halos. For the most probable value of \alpha=1.3, there exists a narrow range in parameter space, such that all constraints are satisfied. However, the adiabatic compression of the dark halo by baryons as they cool and contract in normal galaxies yields a steeper cusp, \alpha=1.7. This gives a tighter constraint, which would exclude SIDM as a possible solution to the purported problems with CDM in the absence of other dynamical processes. Nevertheless, SIDM with parameters consistent with this stronger constraint, can explain the ubiquity of supermassive black holes in the centers of galaxies. A ``best fit'' model is presented which reproduces the supermassive black hole masses and their observed correlations with the velocity dispersion of the host bulges. Specifically, the fourth power dependence of black hole mass on velocity dispersion is a direct consequence of the power spectrum having an index of n=-2. Although the dark matter collision rates for this model are too small to directly remedy problems with CDM, mergers between dark halos harboring supermassive black holes at high redshift could ameliorate the cuspy halo problem. This scenario also explains the lack of comparable supermassive black holes in bulgeless galaxies like M33.

We present a robust method to derive the duty cycle of QSO activity based on the empirical QSO luminosity function and on the present-day linear relation between the masses of supermassive black holes and those of their spheroidal host stellar systems. It is found that the duty cycle is substantially less than unity, with characteristic values in the range $3-6\times 10^{-3}$. Finally, we tested the expectation that the QSO luminosity evolution and the star formation history should be roughly parallel, as a consequence of the above--mentioned relation between BH and galaxy masses.

During the recombination epoch, cosmic background photons couple not only to free electrons through Thompson scattering, but also to the neutral hydrogen through Rayleigh scattering. This latter is ~2% effect for photons near the peak of the photon energy distribution at z=800 and a ~0.2% effect at z=1100. Including Rayleigh scattering in the calculation reduces Silk damping at fixed redshift, alters the position of the surface of last scattering and alters the propagation of acoustic waves. We estimate the amplitude of these effects. For the Microwave Anisotropy Probe (MAP), Rayleigh scattering increases the anisotropy spectrum by 0.1% at the most. For the highest frequencies of the Planck Surveyor, the effects of Rayleigh scattering are much more dramatic (decreasing the anisotropy spectrum by 3% at \nu~550GHz and l~1000). The relative difference between the spectra of low and high frequencies is imposed by an oscillation with a function of multipole l and the oscillation amplitude is up to 0.5% between 100 and 550GHz. Rayleigh scattering also slows the decoupling between radiation and matter, but the effect is undetectably small.

The luminosity function of galaxies is derived from a cosmological hydrodynamic simulation of a Lambda cold dark matter (CDM) universe with the aid of a stellar population synthesis model. At z=0, the resulting B band luminosity function has a flat faint end slope of \alpha \approx -1.15 with the characteristic luminosity and the normalization in a fair agreement with observations, while the dark matter halo mass function is steep with a slope of \alpha \approx -2. The colour distribution of galaxies also agrees well with local observations. We also discuss the evolution of the luminosity function, and the colour distribution of galaxies from z=0 to 5. A large evolution of the characteristic mass in the stellar mass function due to number evolution is compensated by luminosity evolution; the characteristic luminosity increases only by 0.8 mag from z=0 to 2, and then declines towards higher redshift, while the B band luminosity density continues to increase from z=0 to 5 (but only slowly at z>3).

Hydrodynamic simulations of growth of cosmic structure suggest that 30-50% of the total baryons at z=0 may be in a warm-hot intergalactic medium (WHIM) with temperatures ~10^5-10^7K. The O VI \lambda \lambda 1032, 1038 absorption line doublet in the FUV portion of QSO spectra provides an important probe of this gas. Utilizing recent hydrodynamic simulations, it is found that there should be ~5 O VI absorption lines per unit redshift with equivalent widths >= 35 mA, decreasing rapidly to ~0.5 per unit redshift at >= 350 mA. About 10% of the total baryonic matter or 20-30% of the WHIM is expected to be in the O VI absorption line systems with equivalent width >= 20 mA; the remaining WHIM gas may be too hot or have too low metallicity to be detected in O VI. We find that the simulation results agree well with observations with regard to the line abundance and total mass contained in these systems. Some of the O VI systems are collisionally ionized and some are photoionized, but most of the mass is in the collisionally ionized systems. We show that the gas that produces the O VI absorption lines does not reside in virialized regions such as galaxies, groups, or clusters of galaxies, but rather has an overdensity of 10-40 times the average density. These regions form a somewhat connected network of filaments. The typical metallicity of these regions is 0.1-0.3Zsun.

The dynamical evolution of galaxies in clusters is modified if dark matter is self-interacting. Heat conduction from the hot cluster halo leads to evaporation of the relatively cooler galactic halos. The stellar distribution would adiabatically expand as it readjusts to the loss of dark matter, reducing the velocity dispersion and increasing the half-light radius. If the dark matter content within that radius was f_dm = 25-50% of the total, as indicated by current observations, the ellipticals in clusters would be offset from the fundamental plane relation beyond the observational scatter. The requirement that their halos survive for a Hubble time appears to exclude just that range of the dark matter cross-section, 0.3 < sigma/m < 10^4 cm^2/g, thought to be optimal for reducing central halo cusps, unless f_dm < 15%. If the cross-section is allowed to vary with the relative velocity of dark matter particles, sigma \propto v^{-2 delta}, a new problem of evaporation of dark matter arises in the dwarf galaxies with low velocity dispersion. The halos of large galaxies in clusters and dwarf galaxies in the Local Group can both survive only if delta < 1.1 or delta > 1.8. In either case the problem of central density cusps remains.

High-resolution N-body simulations of four popular Cold Dark Matter cosmologies (LCDM, OCDM, QCDM, and tilted SCDM), each containing 10^5 clusters of galaxies in a cubic gigaparsec volume, are used to determine the evolution of the cluster mass function from z=3 to z=0. The large volume and high resolution of these simulations allow an accurate measure of the evolution of cosmologically important (but rare) massive clusters at high redshift. The simulated mass function is presented for cluster masses within several radii typically used observationally (R=0.5, 1.0, and 1.5 Mpc/h, both comoving and physical) in order to enable direct comparison with current and future observations. The simulated evolution is compared with current observations of massive clusters at redshifts 0.3<z<0.8. The Omega_m=1 tilted SCDM model, which exhibits very rapid evolution of the cluster abundance, produces too few clusters at z>0.3 and no massive clusters at z>0.5, in stark contradiction with observations. The Omega_m=0.3 models- LCDM, OCDM, and QCDM- all exhibit considerably weaker evolution and are consistent with current data. Among these low density models, OCDM evolves the least. These trends are enhanced at high redshift and can be used to discriminate between flat and open low density models. The simulated mass functions are compared with the Press-Schechter approximation. Standard Press-Schechter predicts too many low mass clusters at z=0, and too few clusters at higher redshift. We modify the approximation by a simple parameterization of the density contrast threshold for collapse, which has a redshift dependence. This modified Press-Schechter approximation provides a good fit to the simulated mass functions.

We use the semi-analytical approach to analyze gravitational lensing of quasars by dark halos in various cold dark matter (CDM) cosmologies, in order to determine the sensitivity of the prediction probabilities of images separations to the input assumptions regarding halos and cosmologies. The mass function of dark halos is assumed to be given by the Press-Schechter function. The mass density profile of dark halos is alternatively taken to be the singular isothermal sphere (SIS), the Navarro-Frenk-White (NFW) profile, or the generalized NFW profile. The cosmologies include: the Einstein-de Sitter model (SCDM), the open model (OCDM), and the flat \Lambda-model (LCDM). As expected, we find that the lensing probability is extremely sensitive to the mass density profile of dark halos, and somewhat less so to the mean mass density in the universe, and the amplitude of primordial fluctuations. NFW halos are very much less effective in producing multiple images than SIS halos. However, none of these models can completely explain the current observations: the SIS models predict too many large splitting lenses, while the NFW models predict too few small splitting lenses. This indicates that there must be at least two populations of halos in the universe. A combination of SIS and NFW halos can reasonably reproduce the current observations if we choose the mass for the transition from SIS to NFW to be ~ 10^{13} solar masses. Additionally, there is a tendency for CDM models to have too much power on small scales, i.e. too much mass concentration; and it appears that the cures proposed for other apparent difficulties of CDM would help here as well, an example being the warm dark matter (WDM) variant which is shown to produce large splitting lenses fewer than the corresponding CDM model by one order of magnitude.

Discrepancies have emerged between the predictions of standard cold dark matter (CDM) theory and observations of clustering on sub-galactic scales. Warm dark matter (WDM) is a simple modification of CDM in which the dark matter particles have initial velocities due either to their having decoupled as thermal relics, or having been formed via non-equilibrium decay. We investigate the nonlinear gravitational clustering of WDM with a high resolution N-body code, and identify a number of distinctive observational signatures. Relative to CDM, halo concentrations and core densities are lowered, core radii are increased, and large halos emerge with far fewer low mass satellites. The number of small halos is suppressed, and those present are formed by `top down' fragmentation of caustics, as part of a `cosmic web' connecting massive halos. Few small halos form outside this web. If we identify small halos with dwarf galaxies, their number, spatial distribution, and formation epoch appear in better agreement with the observations for WDM than they are for CDM.

We study the constraints that high-redshift structure formation in the universe places on warm dark matter (WDM) dominated cosmological models. We modify the extended Press-Schechter formalism to derive the halo mass function in WDM models. We show that our predictions agree with recent numerical simulations at low redshift over the halo masses of interest. Applying our model to galaxy formation at high redshift, we find that the loss of power on small scales, together with the delayed collapse of low-mass objects, results in strong limits on the root-mean-square velocity dispersion v_rms of the WDM particles at z=0. For fermions decoupling while relativistic, these limits are equivalent to constraints on the mass m_X of the particles. The presence of a 4 billion solar mass black hole at z=5.8, believed to power the quasar SDSS 1044-1215, implies m_X > 0.5 keV (or v_rms < 0.10 km/s), assuming that the quasar is unlensed and radiating at or below the Eddington limit. Reionization by z=5.8 also implies a limit on m_X. If high-redshift galaxies produce ionizing photons with an efficiency similar to their z=3 counterparts, we find m_X > 1.2 keV (or v_rms < 0.03 km/s). However, given the uncertainties in current measurements from the proximity effect of the ionizing background at redshift 3, values of m_X as low as 0.75 keV (v_rms = 0.06 km/s) are not ruled out. The limit weakens further if, instead, the ionizing-photon production efficiency is greater at high z, but this limit will tighten considerably if reionization is shown in the future to have occurred at higher redshifts. WDM models with m_X < 1 keV (v_rms > 0.04 km/s) produce a low-luminosity cutoff in the high-redshift galaxy luminosity function which is directly detectable with the Next Generation Space Telescope (abridged).

We study star formation history and stellar metallicity distribution in galaxies in a Lambda cold dark matter universe using a hydrodynamic cosmological simulation. Our model predicts star formation rate declining in time exponentially from an early epoch to the present with the time-scale of 6 Gyr, which is consistent with the empirical Madau plot with modest dust obscuration. Star formation in L^* galaxies continues intermittently to the present also with an exponentially declining rate of a similar time-scale, whereas in small galaxies star formation ceases at an early epoch. The mean age of the extant stars decreases only slowly with increasing redshift, and exceeds 1 Gyr at z=3. Normal galaxies contain stars with a wide range of metallicity and age: stars formed at z<1 have metallicity of 0.1-1.0 Zsun, while old stars take a wide range of values from 10^{-6} Zsun to 3.0 Zsun. The mean metallicity of normal galaxies is in the range 0.1-1.0 Zsun. Dwarf galaxies that contain only old stars have a wide range of mean metallicity (10^{-4}-1.0 Zsun), but on average they are metal deficient compared with normal galaxies.

It is well known that the mass function for_halos_ in CDM cosmology is a relatively steep power law for low masses, possibly too steep to be consistent with observations. But how steep is the_galaxy_ mass function? We have analyzed the stellar and gas mass functions of the first massive luminous objects formed in a \Lambda CDM universe, as calculated in the numerical simulation described in Gnedin (2000ab). We found that while the dark matter mass function is steep, the stellar and gas mass functions are flatter for low mass objects. The stellar mass function is consistently flat at the low mass end. Moreover, while the gas mass function follows the dark matter mass function until reionization at z~7, between z=7 and z=4, the gas mass function also flattens considerably at the low mass end. At z=4, the gas and stellar mass functions are fit by a Schechter function with \alpha ~ -1.2 +/- 0.1, significantly shallower than the dark matter halo mass function and consistent with some recent observations. The baryonic mass functions are shallower because (a) the dark matter halo mass function is consistent with the Press-Schechter formulation at low masses n(M) M^-2 and (b) heating/cooling and ionization processes appear to cause baryons to collect in halos with the relationship M_b M_d^4 at low masses. Combining (a) and (b) gives n(M_b) M_b^-5/4, comparable to the simulation results. Thus, the well known observational fact that low mass galaxies are underabundant as compared to expectations from numerical dark matter simulations or Press-Schechter modeling of CDM universes emerges naturally from these results, implying that perhaps no ``new physics'' beyond the standard model is needed.

Warm Dark Matter (WDM) models have recently been resurrected to resolve apparent conflicts of Cold Dark Matter (DM) models with observations. Endowing the DM particles with non-negligible velocities causes free-streaming, which suppresses the primordial power spectrum on small scales. The choice of a root-mean-square velocity dispersion v(rms) = 0.05 km/s at redshift z=0 (corresponding to a particle mass of 1 keV if the WDM particles are fermions decoupling while relativistic) helps alleviate most, but probably not all, of the small-scale problems faced by CDM. An important side-effect of the particle velocities is the severe decrease in the number of collapsed halos at high redshift. This is caused both by the loss of small-scale power, and by the delay in the collapse of the smallest individual halos (with masses near the effective Jeans mass of the DM). The presence of early halos is required in order (1) to host either early quasars or galaxies that can reionize the universe by redshift z=5.8, and (2) to allow the growth of the supermassive black hole believed to power the recently discovered quasar SDSS 1044-1215 at this redshift. We quantify these constraints using a modified Press-Schechter formalism, and find v(rms) < 0.04 km/s (or m_X > 1 keV). If future observations uncover massive black holes at z > 10, or reveal that reionization occurred at z > 10, this could conclusively rule out WDM models as the solution to the small-scale crisis of the CDM paradigm.

The Tree-Particle-Mesh (TPM) N-body algorithm couples the tree algorithm for directly computing forces on particles in an hierarchical grouping scheme with the extremely efficient mesh based PM structured approach. The combined TPM algorithm takes advantage of the fact that gravitational forces are linear functions of the density field. Thus one can use domain decomposition to break down the density field into many separate high density regions containing a significant fraction of the mass but residing in a very small fraction of the total volume. In each of these high density regions the gravitational potential is computed via the tree algorithm supplemented by tidal forces from the external density distribution. For the bulk of the volume, forces are computed via the PM algorithm; timesteps in this PM component are large compared to individually determined timesteps in the tree regions. Since each tree region can be treated independently, the algorithm lends itself to very efficient parallelization using message passing. We have tested the new TPM algorithm (a refinement of that originated by Xu 1995) by comparison with results from Ferrell & Bertschinger's P^3M code and find that, except in small clusters, the TPM results are at least as accurate as those obtained with the well-established P^3M algorithm, while taking significantly less computing time. Production runs of 10^9 particles indicate that the new code has great scientific potential when used with distributed computing resources.

We carry out an extensive study of the cosmic Mach number (\mach) on scales of R=5, 10 and 20h^-1Mpc using an LCDM hydrodynamical simulation. We particularly put emphasis on the environmental dependence of \mach on overdensity, galaxy mass, and galaxy age. We start by discussing the difference in the resulting \mach according to different definitions of \mach and different methods of calculation. The simulated Mach numbers are slightly lower than the linear theory predictions even when a non-linear power spectrum was used in the calculation, reflecting the non-linear evolution in the simulation. We find that the observed \mach is higher than the simulated mean by more than 2-standard deviations, which suggests either that the Local Group is in a relatively low-density region or that the true value of \Omega_m is ~ 0.2, significantly lower than the simulated value of 0.37. We show from our simulation that the Mach number is a weakly decreasing function of overdensity. We also investigate the correlations between galaxy age, overdensity and \mach for two different samples of galaxies --- DWARFs and GIANTs. Older systems cluster in higher density regions with lower \mach, while younger ones tend to reside in lower density regions with larger \mach, as expected from the hierarchical structure formation scenario. However, for DWARFs, the correlation is weakened by the fact that some of the oldest DWARFs are left over in low-density regions during the structure formation history. For giant systems, one expects blue-selected samples to have higher \mach than red-selected ones. We briefly comment on the effect of the warm dark matter on the expected Mach number.

At low redshifts, measurements of the total baryon content in stars, atomic and molecular hydrogen, and cluster gas fall a factor of two to four below the baryon density derived from observed light-element ratios and nucleosynthesis arguments. A possible hiding place for a significant fraction of the missing baryons is in the warm/hot gas at temperatures T = 10^{5}-10^{7}K. We present predictions of the contribution to the soft X-ray background from warm/hot gas emission calculated using new hydrodynamical simulations and discuss the possibility of detecting the spectral signature of this gas using the Chandra X-ray Observatory.

Most elliptical galaxies contain central black holes (BHs), and most also contain significant amounts of hot gas capable of accreting on to the central BH due to cooling times short compared to the Hubble time. Why therefore do we not see AGNs at the center of most elliptical galaxies rather than in only (at most) a few per cent of them? We propose here the simple idea that feedback from accretion events heats the ambient gas retarding subsequent infall. In this context, we present a class of 1D hydrodynamical evolutionary sequences for the gas flows in elliptical galaxies with a massive central BH. The resulting evolution is characterized by strong oscillations, in which very fast and energetic bursts of the BH are followed by longer periods in which the X-ray galaxy emission comes from the coronal gas. We also allow for departures from spherical symmetry by examining scenarios in which the central engine is either an ADAF or a more conventional accretion disk that is optically thick except for a polar region. In all cases the duty cycle (fraction of the time that the system will be seen as an AGN) is quite small and in the range 10^{-4} - 10^{-3}. Thus, for any reasonable value of the efficiency, the presence of a massive BH at the center of a galaxy seems to be incompatible with the presence of a long-lived cooling flow.

At low redshifts, a census of the baryons in all known reservoirs falls a factor of two to four below the total baryon density predicted from Big Bang nucleosynthesis arguments and observed light element ratios. Recent cosmological hydrodynamic simulations suggest that a significant fraction of these missing baryons could be in the form of warm/hot gas in the filaments and halos within which most field galaxies are embedded. With the release of source count results from Chandra and recent detections of this gas in O VI quasar absorption lines, it becomes interesting to examine the predictions and limits placed on this component of the X-ray background (XRB). We have used new hydrodynamical simulations to predict the total X-ray spectrum from the gas in the 100 eV to 10 keV range. We find that, when uncertainties in the normalization of the observed XRB and the value of Omega_b are taken into account, our results are consistent with current observational limits placed on the contribution of emission from gas to the XRB. In the 0.5-2 keV range, we expect the contribution from this component to be 0.63 10^{-12} erg s^-1 cm^-2 deg^-2 or between 6% and 18% of the extragalactic surface brightness. The peak fraction occurs in the 0.5-1 keV range where the predicted line emission mirrors a spectral bump seen in the latest ASCA/ROSAT XRB data.

Approximately 30-40% of all baryons in the present day universe reside in a warm-hot intergalactic medium (WHIM), with temperatures between 10^5<T<10^7 K. This is a generic prediction from six hydrodynamic simulations of currently favored structure formation models having a wide variety of numerical methods, input physics, volumes, and spatial resolutions. Most of these warm-hot baryons reside in diffuse large-scale structures with a median overdensity around 10-30, not in virialized objects such as galaxy groups or galactic halos. The evolution of the WHIM is primarily driven by shock heating from gravitational perturbations breaking on mildly nonlinear, non-equilibrium structures such as filaments. Supernova feedback energy and radiative cooling play lesser roles in its evolution. WHIM gas is consistent with observations of the 0.25 keV X-ray background without being significantly heated by non-gravitational processes because the emitting gas is very diffuse. Our results confirm and extend previous work by Cen & Ostriker and Dave' et al.

We use large-scale cosmological simulations to estimate the mass-to-light ratio of galaxy systems as a function of scale, and compare the results with observations of galaxies, groups, clusters, and superclusters of galaxies. We find remarkably good agreement between observations and simulations. Specifically, we find that the simulated mass-to-light ratio increases with scale on small scales and flattens to a constant value on large scales, as suggested by observations. We find that while mass typically follows light on large scales, high overdensity regions --- such as rich clusters and superclusters of galaxies --- exhibit higher M/L_B values than average, while low density regions exhibit lower M/L_B values; high density regions are thus antibiased in M/L_B, with mass more strongly concentrated than blue light. This is true despite the fact that the galaxy mass density is unbiased or positively biased relative to the total mass density in these regions. The M/L_B antibias is likely due to the relatively old age of the high density regions, where light has declined significantly since their early formation time, especially in the blue band which traces recent star formation. Comparing the simulated results with observations, we place a powerful constraint on the mass density of the universe; using, for the first time, the entire observed mass-to-light function, from galaxies to superclusters, we find Omega =0.16+/-0.05.

We have examined the properties of shock waves in simulations of large scale structure formation for two cosmological scenarios (a SCDM and a LCDM with Omega =1). Large-scale shocks result from accretion onto sheets, filaments and Galaxy Clusters (GCs) on a scale of circa 5 Mpc/h in both cases. Energetic motions, both residual of past accretion history and due to current asymmetric inflow along filaments, generate additional, common shocks on a scale of about 1 Mpc/h, which penetrate deep inside GCs. Also collisions between substructures inside GCs form merger shocks. Consequently, the topology of the shocks is very complex and highly connected. During cosmic evolution the comoving shock surface density decreases, reflecting the ongoing structure merger process in both scenarios. Accretion shocks have very high Mach numbers (10-10^3), when photo-heating of the pre-shock gas is not included. The typical shock speed is of order v_{sh}(z) =H(z)lambda_{NL}(z), with lambda_{NL}(z) the wavelength scale of the nonlinear perturbation at the given epoch. However, the Mach number for shocks occuring within clusters is usually smaller (3-10), due to the fact that the intracluster gas is already hot. Statistical fits of shock speed around GCs as a function of GCs temperature give power-law's in accord with 1-D predictions. However, a very different result is obtained for fits of the shock radius, reflecting the very complex shock structures forming in 3-D simulations. The in-flowing kinetic energy across such shocks, giving the power available for cosmic-ray acceleration, is comparable to the cluster X-ray luminosity emitted from a central region of radius 0.5 Mpc/h. Considering their large size and long lifetimes, those shocks are potentially interesting sites for cosmic-ray acceleration, if modest magnetic fields exist within them.

We compute the cosmic star formation rate (SFR) and the rest-frame comoving luminosity density in various pass-bands as a function of redshift using large-scale \Lambda-CDM hydrodynamical simulations with the aim of understanding their behavior as a function of redshift. To calculate the luminosity density of galaxies, we use an updated isochrone synthesis model which takes metallicity variations into account. The computed SFR and the UV-luminosity density have a steep rise from z=0 to 1, a moderate plateau between z=1 - 3, and a gradual decrease beyond z=3. The raw calculated results are significantly above the observed luminosity density, which can be explained either by dust extinction or the possibly inappropriate input parameters of the simulation. We model the dust extinction by introducing a parameter f; the fraction of the total stellar luminosity (not galaxy population) that is heavily obscured and thus only appears in the far-infrared to sub-millimeter wavelength range. When we correct our input parameters, and apply dust extinction with f=0.65, the resulting luminosity density fits various observations reasonably well, including the present stellar mass density, the local B-band galaxy luminosity density, and the FIR-to-submm extragalactic background. Our result is consistent with the picture that \sim 2/3 of the total stellar emission is heavily obscured by dust and observed only in the FIR. The rest of the emission is only moderately obscured which can be observed in the optical to near-IR wavelength range. We also argue that the steep falloff of the SFR from z=1 to 0 is partly due to the shock-heating of the universe at late times, which produces gas which is too hot to easily condense into star-forming regions.

If the cosmological dark matter is primarily in the form of an elementary particle which has cross section and mass for self-interaction having a ratio similar to that of ordinary nuclear matter, then seed black holes (formed in stellar collapse) will grow in a Hubble time, due to accretion of the dark matter, to a mass range 10^6 - 10^9 solar masses. Furthermore, the dependence of the final black hole mass on the galaxy velocity dispersion will be approximately as observed and the growth rate will show a time dependence consistent with observations. Other astrophysical consequences of collisional dark matter and tests of the idea are noted.

We use cosmological hydrodynamic simulations to investigate formation of galactic bulges within the framework of hierarchical clustering in a representative CDM cosmological model. We show that largest objects forming at cosmological redshifts z~4 resemble observed bulges of spiral galaxies or moderate size ellipticals in their general properties like sizes, shapes, and density profiles. This is consistent with observational data indicating the existence of ``old'' bulges and ellipticals at more moderate redshifts. These bulges are gas dominated at redshift z=3, with high rates of star formation and would appear to be good candidates for small blue galaxies seen in the Hubble Deep Field.

All high temperature accretion solutions including ADAF are physically thick, so outgoing radiation interacts with the incoming flow, sharing as much or more resemblance with classical spherical accretion flows as with disk flows. We examine this interaction for the popular ADAF case. We find that without allowance for Compton preheating, a very restricted domain of ADAF solution is permitted and with Compton preheating included a new high temperature PADAF branch appears in the solution space. In the absence of preheating, high temperature flows do not exist when the mass accretion rate mdot == Mdot c^2 / L_E >~ 10^-1.5. Below this mass accretion rate, a roughly conical region around the hole cannot sustain high temperature ions and electrons for all flows having mdot >~ 10^-4, which may lead to a funnel possibly filled with a tenuous hot outgoing wind. If the flow starts at large radii with the usual equilibrium temperature ~10^4 K, the critical mass accretion rate is much lower, mdot \~10^-3.7 above which level no self-consistent ADAF (without preheating) can exist. However, above this critical mass accretion rate, the flow can be self-consistently maintained at high temperature if Compton preheating is considered. These solutions constitute a new branch of solutions as in spherical accretion flows. High temperature PADAF flows can exist above the critical mass accretion rate in addition to the usual cold thin disk solutions. We also find solutions where the flow near the equatorial plane accretes normally while the flow near the pole is overheated by Compton preheating, possibly becoming, a polar wind, solutions which we designate WADAF.

We analyze a new large-scale ($100h^{-1}$Mpc) numerical hydrodynamic simulation of the popular $\Lambda$CDM cosmological model, including in our treatment dark matter, gas and star-formation, on the basis of standard physical processes. The method, applied with a numerical resolution of $<200h^{-1}$kpc (which is still quite coarse for following individual galaxies, especially in dense regions), attempts to estimate where and when galaxies form. We then compare the smoothed galaxy distribution with the smoothed mass distribution to determine the "bias" defined as $b\equiv (\delta M/M)_{gal}/(\delta M/M)_{total}$ on scales large compared with the code numerical resolution (on the basis of resolution tests given in the appendix of this paper). We find that (holding all variables constant except the quoted one) bias increases with decreasing scale, with increasing galactic age or metallicity and with increasing redshift of observations. At the $8h^{-1}$Mpc fiducial comoving scale bias (for bright regions) is 1.35 at $z=0$ reaching to 3.6 at $z=3$, both numbers being consistent with extant observations. We also find that $(10-20)h^{-1}$Mpc voids in the distribution of luminous objects are as observed (i.e., observed voids are not an argument against CDM-like models) and finally that the younger systems should show a colder Hubble flow than do the early type galaxies (a testable proposition). Surprisingly, little evolution is found in the amplitude of the smoothed galaxy-galaxy correlation function (as a function of {\it comoving} separation). Testing this prediction vs observations will allow a comparison between this work and that of Kauffmann et al which is based on a different physical modelingmethod.

We present a comprehensive study of the observational constraints on spatially flat cosmological models containing a mixture of matter and quintessence --- a time varying, spatially inhomogeneous component of the energy density of the universe with negative pressure. Our study also includes the limiting case of a cosmological constant. Low red shift constraints include the Hubble parameter, baryon fraction, cluster abundance, age of the universe, bulk velocity and shape of the mass power spectrum; intermediate red shift constraints are due to type 1a supernovae, gravitational lensing, the Ly-a forest, and the evolution of large scale structure; high red shift constraints are based on cosmic microwave background temperature anisotropy. Mindful of systematic errors, we adopt a conservative approach in applying these constraints. We determine that quintessence models in which the matter density parameter is $0.2 \ls \Omega_m \ls 0.5$ and the effective, density-averaged equation of state is $-1 \le w \ls -0.2$, are consistent with the most reliable, current low red shift and CMB observations at the $2\sigma$ level. Factoring in the constraint due to type 1a SNe, the range for the equation of state is reduced to $-1 \le w \ls -0.4$, where this range represents models consistent with each observational constraint at the 2$\sigma$ level or better (concordance analysis). A combined maximum likelihood analysis suggests a smaller range, $-1 \le w \ls -0.6$. We find that the best-fit and best-motivated quintessence models lie near $\Omega_m \approx 0.33$, $h \approx 0.65$, and spectral index $n_s=1$, with an effective equation of state $w \approx -0.65$ for ``tracker'' quintessence and $w=-1$ for ``creeper'' quintessence. (abstract shortened)

A sample of eight quasars observed at high resolution and signal-to-noise is used to determine the probability distribution function (PDF), the power spectrum, and the correlation function of the transmitted flux in the \lya forest, in three redshift bins centered at z=2.41, 3.00, and 3.89. All the results are presented in tabular form, with full error covariance matrices to allow for comparisons with any numerical simulations and with other data sets. The observations are compared with a numerical simulation of the \lya forest of a Lambda-CDM model with Omega=0.4, known to agree with other large-scale structure observational constraints. There is excellent agreement for the PDF, if the mean transmitted flux is adjusted to match the observations. A small difference between the observed and predicted PDF is found at high fluxes and low redshift, which may be due to the uncertain effects of fitting the spectral continuum. Using the numerical simulation, we show how the flux power spectrum can be used to recover the initial power spectrum of density fluctuations. From our sample of eight quasars, we measure the amplitude of the mass power spectrum to correspond to a linear variance per unit ln(k) of $\Delta^2_\rho(k)=0.72\pm0.09$ at k=0.04(km/s)^{-1} and z=3, and the slope of the power spectrum near the same k to be $n_p=-2.55\pm0.10$ (statistical error bars). The results are statistically consistent with Croft et. al. (1999), although our value for the rms fluctuation is lower by a factor 0.75. For the Lambda-CDM model we use, the implied primordial slope is $n=0.93\pm0.10$, and the normalization is $\sigma_8=0.68+1.16(0.95-n)\pm0.04$.

During hierarchical clustering, smaller masses generally collapse earlier than larger masses and so are denser on the average. The core of a small mass halo could be dense enough to resist disruption and survive undigested, when it is incorporated into a bigger object. We explore the possibility that a nested sequence of undigested cores in the center of the halo, which have survived the hierarchical, inhomogeneous collapse to form larger and larger objects, determines the halo structure in the inner regions. For a flat universe with $P(k) \propto k^n$, scaling arguments then suggest that the core density profile is, $\rho \propto r^{-\alpha}$ with $\alpha = (9+3n)/(5+n)$. But whether such behaviour obtains depends on detailed dynamics. We first examine the dynamics using a fluid approach to the self-similar collapse solutions for the dark matter phase space density, including the effect of velocity dispersions. We highlight the importance of tangential velocity dispersions to obtain density profiles shallower than $1/r^2$ in the core regions. If tangential velocity dispersions in the core are constrained to be less than the radial dispersion, a cuspy core density profile shallower than 1/r cannot obtain, in self-similar collapse. We then briefly look at the profiles of the outer halos in low density cosmological models where the total halo mass is convergent. Finally, we analyze a suite of dark halo density and velocity dispersion profiles obtained in cosmological N-body simulations of models with n= 0, -1 and -2. We find that the core-density profiles of dark halos, show considerable scatter in their properties, but nevertheless do appear to reflect a memory of the initial power spectrum, with steeper initial spectra producing flatter core profiles. (Abridged)

The "cosmic triangle" is introduced as a way of representing the past, present, and future status of the universe. Our current location within the cosmic triangle is determined by the answers to three questions: How much matter is in the universe? Is the expansion rate slowing down or speeding up? And, is the universe flat? A review of recent observations suggests a universe that is lightweight (matter density about one-third the critical value), is accelerating, and is flat. The acceleration implies the existence of cosmic dark energy that overcomes the gravitational self-attraction of matter and causes the expansion to speed up.

We present a semi-analytic model for the thermal and ionization history of the universe at 1000 >~ z >~ 3. This model incorporates much of the essential physics included in full-scale hydrodynamical simulations, such as (1) gravitational collapse and virialization; (2) star/quasar formation and subsequent ionizing radiation; (3) heating and cooling; (4) atomic and molecular physics of hydrogen; and (5) the feedback relationships between these processes. In addition, we model the process of reheating and reionization using two separate phases, self-consistently calculating the filling factor of each phase. Thus radiative transfer is treated more accurately than simulations to date have done: we allow to lowest order for both the inhomogeneity of the sources and the sinks of radiation. After calibrating and checking the results of this model against a hydrodynamical simulation, we apply our model to a variety of Gaussian and non-Gaussian CDM-dominated cosmological models. Our major conclusions include: (1) the epoch of reheating depends most strongly on the power spectrum of density fluctuations at small scales; (2) because of the effects of gas clumping, full reionization occurs at z ~ 10 in all models; (3) the CMBR polarization and the stars and quasars to baryons ratio are strong potential discrimants between different assumed power spectra; (4) the formation of galactic spheroids may be regulated by the evolution of reheating through feedback, so that the Jeams mass tracks the non-linear mass scale; and (5) the evolution of the bias of luminous objects can potentially discriminate strongly between Gaussian and non-Gaussian PDFs.

Using large-scale hydrodynamic simulations with heuristic criteria for galaxy formation, we investigate how the galaxy field is related to the mass density and the gas temperature. We find that the bias b(R) = sigma_g(R)/sigma(R), where sigma_g(R) is the variance of galaxy counts in spheres of radius R and sigma(R) is the same for mass, varies from 2.6 at 1 h^{-1} Mpc to 1.2 at 30 h^{-1} Mpc. Including the dependence of the galaxy density on local gas temperature as well as on local mass density can fully account for this scale dependence. Galaxy density depends on temperature because gas which is too hot cannot cool to form galaxies; this causes scale dependence of b(R) because local gas temperature is related to the gravitational potential, and thus contains information about the large scale density field. We show that temperature dependence generally causes b(R) to vary on quasilinear and nonlinear scales, indicating that scale dependence of bias may be a generic effect in realistic galaxy formation scenarios. We find that the relationship between the galaxy and mass density fields is also a function of galaxy age. On large scales, the older galaxies are highly biased (b ~ 1.7) and highly correlated (r ~ 1.0) with the mass density field; younger galaxies are not biased (b ~ 0.8) and are poorly correlated (r ~ 0.5) with the mass. We argue that linear bias is inadequate to describe the relationship between galaxies and mass. (Abridged)

We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.

We present new Fokker-Planck models of the evolution of globular clusters, including gravitational tidal shocks. We extend our calculations beyond the core collapse by adopting three-body binary heating. Effects of the shocks are included by adding the tidal shock diffusion coefficients to the ordinary Fokker-Planck equation: the first order heating term, <dE>, and the second order energy dispersion term, <dE^2>. As an example, we investigate the evolution of models for the globular cluster NGC 6254. Using the Hipparcos proper motions, we are now able to construct orbits of this cluster in the Galaxy. Tidal shocks accelerate significantly both core collapse and the evaporation of the cluster and shorten the destruction time from 24 Gyr to 18 Gyr. We examine various types of adiabatic corrections and find that they are critical for accurate calculation of the evolution. Without adiabatic corrections, the destruction time of the cluster is twice as short. We examine cluster evolution for a wide range of the concentration and tidal shock parameters, and determine the region of the parameter space where tidal shocks dominate the evolution. We present fitting formulae for the core collapse time and the destruction time, covering all reasonable initial conditions. In the limit of strong shocks, the typical value of the core collapse time decreases from 10 t_{rh} to 3 t_{rh} or less, while the destruction time is just twice that number. The effects of tidal shocks are rapidly self-limiting: as clusters lose mass and become more compact, the importance of the shocks diminishes. This implies that tidal shocks were more important in the past.

Numerical simulations of standard cosmological scenarios have now reached the degree of sophistication required to provide tentative answers to the fundamental question: Where and when were the heavy elements formed? Averaging globally, these simulations give a metallicity that increases from 1% of the solar value at $z=3$ to 20% at present. This conclusion is, in fact, misleading, as it masks the very strong dependency of metallicity on local density. At every epoch higher density regions have much higher metallicity than lower density regions. Moreover, the highest density regions quickly approach near solar metallicity and then saturate, while more typical regions slowly catch up. These results are much more consistent with observational data than the simpler picture (adopted by many) of gradual, quasi-uniform increase of metallicity with time.

The clustering of galaxies relative to the mass distribution declines with time because: first, nonlinear peaks become less rare events; second, the densest regions stop forming new galaxies because gas there becomes too hot to cool and collapse; third, after galaxies form, they are gravitationally ``debiased'' because their velocity field is the same as the dark matter. To show these effects, we perform a hydrodynamic cosmological simulation and examine the density field of recently formed galaxies as a function of redshift. We find the bias b_* of recently formed galaxies (the ratio of the rms fluctuations of these galaxies and mass), evolves from 4.5 at z=3 to around 1 at z=0, on 8 h^{-1} Mpc comoving scales. The correlation coefficient r_* between recently formed galaxies and mass evolves from 0.9 at z=3 to 0.25 at z=0. As gas in the universe heats up and prevents star formation, star-forming galaxies become poorer tracers of the mass density field. After galaxies form, the linear continuity equation is a good approximation to the gravitational debiasing, even on nonlinear scales. The most interesting observational consequence of the simulations is that the linear regression of the star-formation density field on the galaxy density field evolves from about 0.9 at z=1 to 0.35 at z=0. These effects also provide a possible explanation for the Butcher-Oemler effect, the excess of blue galaxies in clusters at redshift z ~ 0.5. Finally, we examine cluster mass-to-light ratio estimates of Omega, finding that while Omega(z) increases with z, one's estimate Omega_est(z) decreases. (Abridged)

We derive expressions for the tidal field exerted by a spherically symmetric galaxy having an extended mass distribution, and use our analysis to calculate tidal perturbations and heating of stars in a globular cluster or a satellite galaxy orbiting in the external potential. We consider a Hernquist model and several power-law models for the primary. We integrate the tidal force over accurate orbits in the given potential by using the impulse approximation and by introducing adiabatic corrections that allow for the time variation of the force. We compare these results with the conventional straight path approximation. We also check the results with direct N-body simulations and find very good agreement. Heating on highly eccentric orbits dominates as the adiabatic corrections strongly damp the energy change on low eccentricity orbits. The results are illustrated for the example of globular cluster NGC 6712. For the orbital eccentricities higher than 0.9 the future lifetime of NGC 6712 is less than 10^{10} yr. Our analysis can be used to study tidal interactions of star clusters and galaxies in a semi-analytical manner.

We study the reaction of a globular star cluster to a time-varying tidal perturbation (gravitational shock) using self-consistent N-body simulations and address two questions. First, to what extent is the cluster interior protected by adiabatic invariants. Second, how much further energy change does the postshock evolution of the cluster potential produce and how much does it affect the dispersion of stellar energies. We introduce the adiabatic correction} as ratio of the energy change, <dE>, to its value in the impulse approximation. When the potential is kept fixed, the numerical results for the adiabatic correction for stars with orbital frequency \omega can be approximated as (1 + \omega^2 \tau^2)^{-\gamma}. For shocks with the characteristic duration of the order the half-mass dynamical time of the cluster, \tau < t_{dyn,h}, the exponent \gamma = 5/2. For more prolonged shocks, \tau > 4 t_{dyn,h}, the adiabatic correction is shallower, \gamma = 3/2. When we allow for self-gravity and potential oscillations which follow the shock, the energy of stars in the core changes significantly, while the total energy of the system is conserved. Paradoxically, the postshock potential fluctuations reduce the total amount of energy dispersion, <dE^2>. The effect is small but real and is due to the postshock energy change being statistically anti-correlated with the shock induced heating. These results are to be applied to Fokker-Planck models of the evolution of globular clusters.

We study the thermal structure of the widely adopted two-dimensional advection dominated accretion flow (ADAF) of Narayan & Yi (1995a). The critical radius for a given mass accretion rate, outside of which the optically thin hot solutions do not exist in the equatorial plane, agrees with one-dimensional study. However, we find that, even within the critical radius, there always exists a conical region of the flow, around the pole, which cannot maintain the assumed high electron temperature, regardless of the mass accretion rate, in the absence of radiative heating. This could lead to torus-like advection inflow shape since, in general, the ions too will cool down. We also find that Compton preheating is generally important and, if the radiative efficiency, defined as the luminosity output divided by the mass accretion rate times the velocity of light squared, is above sim 4x10^-3, the polar region of the flow is preheated above the virial temperature by Compton heating and it may result in time-dependent behaviour or outflow while accretion continues in the equatorial plane. Thus, under most relevant circumstances, ADAF solutions may be expected to be accompanied by polar outflow winds. While preheating instabilities exist in ADAF, as for spherical flows, the former are to some extent protected by their characteristically higher densities and higher cooling rates, which reduce their susceptibility to Compton driven overheating.

Spherical and axisymmetric accretion onto black holes is discussed. Physical processes in various families of solutions are explained and their characteristics are summarized. Recently discovered solutions of axisymmetric flow provide us with various radiation efficiency and spectrum, which may successfully model diverse accretion systems. Possible role of preheating is also speculated. The various families of solutions can be plotted as trajectories on the plane of (l,e) [or (l,dot m)] where l is the luminosity in units of the Eddington luminosity, dot m is the similarly defined mass accretion rate and e is the efficiency defined by e equiv L/dot M c^2 = l/dot m. We discuss the domains on these planes where solutions are known or expected to be unstable to either spherical or on-spherical perturbations. Preliminary analysis indicates that a preheating instability will occur along the polar direction of the advection dominated flow for e \gsim 10^{-2}.

Why do galactic bars rotate with high pattern speeds, when dynamical friction should rapidly couple the bar to the massive, slowly rotating dark halo? This longstanding paradox may be resolved by considering the dynamical interactions between the galactic disk and structures in the dark halo. Dynamical friction between small-scale halo structure and the disk spins up and flattens the inner halo, thereby quenching the dynamical friction exerted by the halo on the bar; at the same time the halo heats and thickens the disk, perhaps forming a rapidly rotating bulge. Two possible candidates for the required halo structures are massive black holes and tidal streamers from disrupted precursor halos. More generally, gravitational scattering from phase-wrapped inhomogeneities represents a novel relaxation process in stellar systems, intermediate between violent relaxation and two-body relaxation, which can isotropize the distribution function at radii where two-body relaxation is not effective.

We perform a variety of tests to determine the numerical resolution of the cosmological TVD eulerian code developed by Ryu et al (1993). Tests include 512^3 and 256^3 simulations of a Pk=k^{-1} spectrum to check for self-similarity and comparison of results with those from higher resolution SPH and grid-based calculations (Frenk et al 1998). We conclude that in regions where density gradients are not produced by shocks the code degrades resolution with a Gaussian smoothing (radius) length of 1.7 cells. At shock caused gradients (for which the code was designed) the smoothing length is 1.1 cells. Finally, for \beta model fit clusters, we can approximately correct numerical resolution by the transformation R^2_{core}\to R^2_{core}-(C\Delta l)^2, where \Delta l is the cell size and C=1.1-1.7. When we use these corrections on our previously published computations for the SCDM and \Lambda CDM models we find luminosity weighted, zero redshift, X-ray cluster core radii of (210\pm 86, 280\pm 67)h^{-1}kpc, respectively, which are marginally consistent with observed (Jones & Forman 1992) values of 50-200h^{-1}kpc. Using the corrected core radii, the COBE normalized SCDM model predicts the number of bright L_x>10^{43}erg/s clusters too high by a factor of \sim 20 and the \Lambda CDM model is consistent with observations.

We use high-resolution hydrodynamic simulations to investigate the density profile of hot gas in clusters of galaxies, adopting a variant of cold dark matter cosmologies and employing a cosmological N-body/smoothed particle hydrodynamics code to follow the evolution of dark matter and gas. In addition to gravitational interactions, gas pressure, and shock heating, we include bremsstrahlung cooling in the computation. Dynamical time, two-body relaxation time, and cooling time in the simulations are examined to demonstrate that the results are free from artificial relaxation effects and that the time step is short enough to accurately follow the evolution of the system. In the simulation with nominal resolution of 66h^{-1} kpc the computed cluster appears normal, but in a higher (by a factor 2) resolution run, cooling is so efficient that the final gas density profile shows a steep rise toward the cluster center that is not observed in real clusters. Also, the X-ray luminosity of 7\times10^{45} ergs s^{-1} far exceeds that for any cluster of the computed temperature. The most reasonable explanation for this discrepancy is that there are some physical processes still missing in the simulations that actually mitigate the cooling effect and play a crucial role in the thermal and dynamical evolution of the gas near the center. Among the promising candidate processes are heat conduction and heat input from supernovae. We discuss the extent to which these processes can alter the evolution of gas.

New, high resolution, large-scale, cosmological hydrodynamic galaxy formation simulations of a standard cold dark matter model (with a cosmological constant) are utilized to predict the distribution of baryons at the present and at moderate redshift. It is found that the average temperature of baryons is an increasing function of time, with most of the baryons at the present time having a temperature in the range 10^{5-7} K. Thus, not only is the universe dominated by dark matter, but more than one half of the normal matter is yet to be detected. Detection of this warm/hot gas poses an observational challenge, requiring sensitive EUV and X-ray satellites. Signatures include a soft, cosmic X-ray background, apparent warm components in hot clusters due to both intrinsic warm intra-cluster gas and warm inter-cluster gas projected onto clusters along the line of sight, absorption lines in X-ray and UV quasar spectra [e.g., O VI (1032,1038)A lines, OVII 574 eV line], strong emission lines (e.g., O VIII 653 eV line) and low redshift, broad, low column density $\lya$ absorption lines. We estimate that approximately 1/4 of the extragalactic soft X-ray background (SXRB) (at 0.7 keV) arises from the warm/hot gas, half of it coming from $z<0.65$ and three-quarters from $z<1.00$, so the source regions should be identifiable on deep optical images.

The semi-analytic theory of tidal shocks proves to be a powerful tool to study tidal interactions of star clusters and satellite galaxies with their massive hosts. New models of the globular cluster evolution employ a combination of analytic estimates, solutions of the Fokker-Planck equation and direct N-body simulations. The models predict large destruction rates for the Galactic globular clusters. Those on the highly eccentric orbits around the Galactic center are much more likely to be disrupted than the ones on nearly circular orbits. The destruction rates are largely increased near the bulge. Disruption of the low-mass clusters changes the Luminosity Function of the Globular Cluster System, shifting the peak of the Luminosity Function to the brighter end.

We propose a new method of analysis for the \lya forest, namely to measure the 1-point and 2-point joint probability distribution of the transmitted flux. The results for a sample of seven observed quasars and from two simulations of structure formation are shown and compared. Statistically significant differences in the 2-point function between the results of the numerical simulations and the observations are easily found. The analysis we suggest is very simple to apply to observed data sets, and we discuss its superiority over the traditional Voigt-profile fitting algorithms for accurate comparison to the predictions of theoretical models.

We examine the correlational property of $\lya$ clouds in detail and compare it to that of mass and galaxies. We show that the correlation strength of $\lya$ clouds is somewhat weaker than that of the underlying matter, which in turn is weaker than that of galaxies (biased galaxy formation). On the scales probed, $10-300$km/s, for $\lya$ clouds, we find that higher density, higher optical depth, higher column density regions are more strongly clustered than lower density, lower optical depth, lower column density regions. Thus, a consistent picture seems to emerge: the correlation strength for a given set of objects is positively correlated with their characteristic global density and the differences among the correlations of galaxies, $\lya$ clouds and mass reflect the differences in density that each trace. Significant positive correlations with a strength of $0.1-1.0$ are found for $\lya$ clouds in the velocity range $50-300$km/s. Among the correlational measures examined, an optical depth correlation function (Equation 5) proposed here may serve as the best correlational measure. It reasonably faithfully represents the true correlation of the underlying matter, enabling a better indication of both matter correlation and the relationship between galaxies and $\lya$ clouds. Furthermore, it appears to be an alternative to the conventional line-line correlation function with the virtue that it does not require ambiguous post-observation fitting procedures such as those commonly employed in the conventional line-finding methods. Neither does it depend sensitively on the observational resolution (e.g., FWHM), insofaras the clouds are resolved Conveniently, it can be easily measured with the current observational sensitivity without being contaminated significantly by the presence of noise.

(Abridged) By comparing the frequency of typical events with that of unusual events, one can test whether the cosmological density distribution function is consistent with the normally made assumption of Gaussianity. To this end, we compare the consistency of the tail-inferred (from clusters) and measured values (from large-scale flows) of the rms level of mass fluctuations for two distribution functions: a Gaussian, and a texture (positively-skewed) PDF. Averaging the recent large-scale flow measurements, we find that observations of the rms and the tail at the 10 h^-1 Mpc scale disfavor a texture PDF at ~1.5 sigma in all cases. However, taking only the most recent measurement of the rms, that from Willick et al. (1997b), the comparison disfavors textures for low Omega_0=0.3, and disfavors Gaussian models if Omega_0=1 (again at ~1.5 sigma). Predictions for evolution of high temperature clusters can also be made for the models considered, and strongly disfavor Omega_0=1 in Gaussian models and marginally disfavor Omega_0=1 in texture models. Only Omega_0=0.3 Gaussian models are consistent with all the data considered.

We present a new class of solutions for the gas flows in elliptical galaxies containing massive central black holes (BH). Modified King model galaxies are assumed. Two source terms operate: mass loss from evolving stars, and a secularly declining heating by SNIa. Relevant atomic physical processes are modeled in detail. Like the previous models investigated by Ciotti et al. (1991), these new models first evolve through three consecutive evolutionary stages: wind, outflow, and inflow. At this point the presence of the BH alters dramatically the subsequent evolution, because the energy emitted by the BH can heat the surrounding gas to above virial temperatures, causing the formation of a hot expanding central bubble. Short and strong nuclear bursts of radiation are followed by longer periods during which the X-ray galaxy emission comes from the coronal gas (Lx). The range and approximate distribution spanned by Lx are found to be in accordance with observations of X-ray early type galaxies. Moreover, although high accretion rates occur during bursting phases when the central BH has a luminosity characteristic of QSOs, the total mass accreted is very small when compared to that predicted by stationary cooling-flow solutions and computed masses are in accord with putative BH nuclear masses. In the bursting phases Lx is low and the surface brightness profile is very low compared to pre-burst or to cooling flow models. We propose that these new models, while solving some long-standing problems of the cooling flow scenario, can provide a unified description of QSO-like objects and X-ray emitting elliptical galaxies, these being the same objects observed at two different evolutionary phases.

The standard method of measuring the galaxy pairwise velocity dispersion on small scales is heavily weighted by the densest regions in a way that is difficult to calibrate. We propose a new statistic which measures the small-scale velocity dispersion as an explicit function of density. Computing this statistic for a volume-limited subsample of the Optical Redshift Survey, we find that the small-scale velocity dispersion rises from 220 to 760 km/s as density increases. We calculate this statistic for a series of mock catalogs drawn from a hydrodynamic simulation of an \Omega h = 0.5 Cold Dark Matter universe (standard CDM), and find that the observed velocity distribution lies ~1 \sigma below the simulations in each of eight density bins, formally ruling out this model at the 7.4 \sigma level, quantifying the well-known problem that this model produces too high a velocity dispersion. This comparison is insensitive to the normalization of the power spectrum, although it is quite sensitive to the density and velocity bias of galaxies relative to dark matter on small scales. (abridged)

We have previously found in the Hubble Deep Field a significant angular correlation of faint, high color-redshift objects on scales below one arcsecond, or several kiloparsecs in metric size. We examine the correlation and nearest neighbor statistics to conclude that 38% of these objects in the HDF have a companion within one arcsecond, three times the number expected in a random distribution with the same number of objects. We examine three dynamical scenarios for these object multiplets: 1) the objects are star-forming regions within normal galaxies, whose disks have been relatively dimmed by K-correction and surface brightness dimming; 2) they are fragments merging into large galaxies; 3) they are satellites accreting onto normal L_* galaxies. We find that hypothesis 1 is most tenable. First, large galaxies in the process of a merger formation would have accumulated too much mass in their centers (5e12 M_sun inside 2 kpc) to correspond to present day objects. Second, accretion by dynamical friction occurs with a predictable density vs. radius slope, not seen among the faint HDF objects. Since the dynamical friction time is roughly (1 Gyr), a steady-state should have been reached by redshift z < 5. Star-forming regions within galaxies clearly present no dynamical problems. Since large spirals would still appear as such in the HDF, we favor a scenario in which the faint compact sources in the HDF are giant starforming regions within small normal galaxies, such as Magellanic irregulars. Finally we checked that reduction in mass-to-light from induced star-formation cannot alone explain the luminosity overdensity.

We have measured the distribution function of the flux decrement caused by Lyman alpha forest absorption in a new sample of high resolution QSO spectra. The observations are compared to the results from two simulations of the Lya forest: an Eulerian Lambda-CDM model, and an SPH standard CDM model. Good agreement between the shapes of simulated and observed distributions is achieved by globally scaling the optical depth to match the mean flux decrements. This procedure amounts to a measurement of the parameter Omega_b^2 h^3 / Gamma (where Omega_b is the baryonic matter density and Gamma is the HI ionization rate). Estimating a lower limit Gamma > 7 10^{-13} s^{-1} from the abundance of known QSOs, we derive a lower limit to the baryon density, Omega_b h^2>0.021(0.017) for the Lambda-CDM (SCDM) model. In both cases the large values are inconsistent with some recent D/H determinations (Rugers & Hogan 1996a,b), favoring a low deuterium abundance as reported by Tytler, Fan & Burles (1996). Adopting a fixed Omega_b, we can determine the evolution of the ionizing radiation field. Our models predict the intensity to be approximately constant with redshift, consistent with the assumption that the ionizing background is produced by known quasars for z < 3. However, additional sources of ionizing photons are required at higher redshift.

We simulate a plausible cosmological model in considerable physical and numerical detail through the successive phases of reheating (at 10<z<20), formation of Pop III stars at z=15 (due to molecular hydrogen cooling), with subsequent reionization at z=7. We assume an efficiency of high mass star formation appropriate to leave the universe, after it becomes transparent, with an ionizing background J_21=0.4 (at z=4), near (and perhaps slightly below) the observed value. Since the same stars produce the ionizing radiation and the first generation of heavy elements, a mean metallicity of Z/Z_\sun=1/200 is produced in this early phase, but there is a large variation about this mean, with the high density regions having Z/Z_\sun=1/30 and low density regions essentially no metals. Reionization, when it occurs, is very rapid, which will leave a signature which may be detectable by very large area meter-wavelength radio instruments. Also, the background UV radiation field will show a sharp drop from 1Ryd to 4Ryd due to absorption edges. The simulated volume is too small to form L_* galaxies, but the smaller objects which are found in the simulation obey the Faber-Jackson relation. In order to explore theoretically this domain of "the end of the dark ages" quantitatively, numerical simulations must have a mass resolution of the order of 10^{4.5} M_\sun} in baryons, high spatial resolution (1 kpc) to resolve strong clumping, and allow for detailed and accurate treatment of both the radiation field and atomic/molecular physics.

The dynamical evolution of the globular cluster systems in galaxies is predicted, based on the standard dynamical theory normalized to the example of the Milky Way. The major processes varying with the galactocentric distance are the tidal shocks and dynamical friction. Our simple model explains, on a quantitative basis, the observed differences of the inner and outer populations of globular clusters. We can thus calculate corrections for dynamical evolution for the luminosity function of globular clusters with the assumption that the initial luminosity function is identical in all galaxies (and we can test this assumption as well, in certain cases). Then we can compute the expected distribution of absolute magnitudes and compare it with the observed distribution of apparent magnitudes to estimate the distance moduli for M31 and M87. Using this new method we find dm(M31)=24.05 +- 0.23, dm(M87)=30.83 +- 0.17, as compared to current best estimates using other methods of dm(M31)=24.30 +- 0.20, dm(M87)=31.0 +- 0.1. As a check on the method we compute, and compare with observations, the differences between the inner and outer globular clusters in all three galaxies. This new method, coupled with HST observations, promises to provide an independent method of estimating distances to galaxies with recession velocities < 10,000 km/s, or D < 100 Mpc.

Gravitational lensing directly measures mass density fluctuations along the lines of sight to very distant objects. No assumptions need to be made concerning bias, the ratio of fluctuations in galaxy density to mass density. Hence, lensing is a very useful tool to study the universe at low to moderate redshifts. We describe in detail a new method to trace light rays from redshift zero through a three dimensional mass distribution to high redshift. As an example, this method is applied here to a standard cold dark matter universe. We obtain a variety of results, some of them statistical in nature, others from rather detailed case studies of individual ``lines of sight". Among the former are the frequency of multiply imaged quasars, the distribution of separation of the multiple quasars, and the redshift distribution of lenses: all that as a function of quasar redshift. We find effects from very weak lensing to highly magnified multiple images of high redshift objects, which, for extended background sources, (i.e. galaxies), range from slight deformations of the shapes through tangentially aligned arclets up to giant luminous arcs. Different cosmological models differ, increasingly with redshift, in their predictions for the mass (thus gravitational potential) distributions. Our ultimate goal is to apply this method to a number of cosmogonic models and to eliminate some models whose gravitational lensing properties are inconsistent with those observed.

We compute the two-point angular correlation function and number-magnitude relation of Hubble Deep Field sources in order to assess their nature. We find that the correlation peaks between 0.25 arcsec and 0.4 arcsec with amplitude of 2 or greater, and much more for the smallest objects. This angular scale corresponds to physical scales of order 1 kpc for redshifts z > 1. The correlation must therefore derive from objects with subgalaxian separations. At faint magnitudes, the counts satisfy the relation Number prop. to 1/flux, expected for images which are subdivisions of larger ones. A conservative conjecture may explain these results. Since high redshift space (z > 0.5) dominates the volume of the sample, observational redshift effects are important. The K-correction and surface brightness dimming of diffuse sources enhances the prominence of compact, unresolved, UV-bright objects, such as star-forming regions within normal gas-rich galaxies. If such regions appear as individual sources, they could explain the subgalaxian correlation.

We note that current observational evidence strongly favors a conventional recombination of ionized matter subsequent to redshift z=1200, followed by reionization prior to redshift z=5 and compute how this would have occurred in a standard scenario for the growth of structure. Extending prior semi-analytic work, we show by direct, high-resolution numerical simulations (of a COBE normalized CDM+Lambda model) that reheating, will occur in the interval 15>z>7, followed by reionization and accompanied by a significant increase in the Jeans mass. However, the evolution of the Jeans mass does not significantly affect star formation in dense, self-shielded clumps of gas, which are detached from the thermal evolution of the rest of the universe. On average, the growth of the Jeans mass tracks the growth of the nonlinear mass scale, a result we suspect is due to nonlinear feedback effects. Cooling on molecular hydrogen leads to a burst of star formation prior to reheating which produces Population III stars with Omega_* reaching 10^{-5.5} and Z/Z_sun reaching 10^{-3.7} by z=14. Star formation subsequently slows down as molecular hydrogen is depleted by photo-destruction and the rise of the temperature. At later times, z<10, when the characteristic virial temperature of gas clumps reach 10,000 degrees, star formation increases again as hydrogen line cooling become efficient. Objects containing Pop III stars accrete mass with time and, as soon as they reach 10,000 K virial temperature, they engage in renewed star formation and turn into normal Pop II objects having an old Pop III metal poor component.

It is demonstrated that strong magnetic fields are produced from a zero initial magnetic field during the pregalactic era, when galaxies are first forming. Their development proceeds in three phases. In the first phase, weak magnetic fields are created by the Biermann battery mechanism, acting in shocked parts of the intergalactic medium where caustics form and intersect. In the second phase, these weak magnetic fields are amplified to strong magnetic fields by the Kolmogoroff turbulence endemic to gravitational structure formation of galaxies. During this second phase, the magnetic fields reach saturation with the turbulent power, but they are coherent only on the scale of the smallest eddy. In the third phase, the magnetic field strength increases to equipartition with the turbulent energy, and the coherence length of the magnetic fields increases to the scale of the largest turbulent eddy, comparable to the scale of the entire galaxy. The resulting magnetic field represents a galactic magnetic field of primordial origin. No further dynamo action is necessary, after the galaxy forms, to explain the origin of magnetic fields. However, the magnetic field may be altered by dynamo action once the galaxy and the galactic disk have formed. It is first shown by direct numerical simulations, that thermoelectric currentsassociated with the Biermann battery, build the field up from zero to $10^{-21}$ G in the regions about to collapse into galaxies, by $z\sim3$. For weak fields, in the absence of dissipation, the cyclotron frequency ${\bf \omega_{cyc}}=e{\bf B } /m_H c $ and $ {\bf \omega}/(1+ \chi )$, where ${\bf \omega = \nabla \times v }$ is the vorticity and $\chi$ is the degree of ionization, satisfy the same equations, and initial conditions ${\bf \omega_{cyc}=\omega}=0$, so that, ${\bf \omega_{cyc}}({\bf

Weak gravitational lensing by large-scale structure affects the determination of the cosmological deceleration parameter $q_0$. We find that the lensing induced dispersions on truly standard candles are $0.04$ and $0.02$ mag at redshift $z=1$ and $z=0.5$, respectively, in a COBE-normalized cold dark matter universe with $\Omega_0=0.40$, $\Lambda_0=0.6$, $H=65$km/s/Mpc and $\sigma_8=0.79$. It is shown that one would observe $q_0=-0.44^{+0.17}_{-0.05}$ and $q_0=-0.45^{+0.10}_{-0.03}$ (the errorbars are $2\sigma$ limits) with standard candles with zero intrinsic dispersion at redshift $z=1$ and $z=0.5$, respectively, compared to the truth of $q_0=-0.40$ in this case, i.e., a 10\% error in $q_0$ will be made. A standard COBE normalized $\Omega_0=1$ CDM model would produce three times as much variance and a mixed (hot and cold) dark matter model would lead to an intermediate result. One unique signature of this dispersion effect is its non Gaussianity. Although the lensing induced dispersion at lower redshift is still significantly smaller than the currently best observed (total) dispersion of $0.12$ mag in a sample of type Ia supernovae, selected with the multicolor light curve shape method, it becomes significant at higher redshift. We show that there is an optimal redshift, in the range $z\sim 0.5-2.0$ depending on the amplitude of the intrinsic dispersion of the standard candles, at which $q_0$ can be most accurately determined.

The nonlinear clustering of dark matter particles in an expanding universe is usually studied by N-body simulations. One can gain some insight into this complex problem if simple relations between physical quantities in the linear and nonlinear regimes can be extracted from the results of N-body simulations. Hamilton et al. (1991) and Nityananda and Padmanabhan (1994) have made an attempt in this direction by relating the mean relative pair velocities to the mean correlation function in a useful manner. We investigate this relation and other closely related issues in detail for the case of six different power spectra: power laws with spectral indexes $n=-2,-1$, cold dark matter (CDM), and hot dark matter models with density parameter $\Omega=1$; CDM including a cosmological constant ($\Lambda$) with $\Omega_{CDM}=0.4$, $\Omega_{\Lambda}=0.6$; and $n=-1$ model with $\Omega=0.1$. We find that: (i) Power law spectra lead to self-similar evolution in an $\Omega=1$ universe. (ii) Stable clustering does not hold in an $\Omega=1$ universe to the extent our simulations can ascertain. (iii) Stable clustering is a better approximation in the case of $\Omega<1$ universe in which structure formation freezes out at some low redshift. (iv) The relation between dimensionless pair velocity and the mean correlation function, $\bar\xi$, is only approximately independent of the shape of the power spectrum. At the nonlinear end, the asymptotic value of the dimensionless pair velocity decreases with increasing small scale power, because the stable clustering assumption is not universally true. (v) The relation between the evolved $\bar\xi$ and the linear regime $\bar\xi$ is also not universal but shows a weak spectrum dependence. We present simple theoretical arguments for these conclusions.

We investigate the dynamical evolution of the Galactic Globular Cluster System in considerably greater detail than has been done hitherto, finding that destruction rates are significantly larger than given by previous estimates. More than half of the present clusters (52%-58% for the OC galactic model, and 75%-86% for the BSS model) will be destroyed in the next Hubble time. For the evolution of individual clusters we use a Fokker-Planck code including two-body relaxation, tidal truncation of clusters, compressive gravitational shocks while clusters pass through the Galactic disk, and tidal shocks due to passage close to the bulge. Gravitational shocks are treated comprehensively, using a recent result by Kundic & Ostriker (1995) that the shock-induced relaxation term, driving an additional dispersion of energies, is generally more important than the usual energy shift term. We discuss some implications for a past history of the Globular Cluster System, and the initial distribution of the destruction times raising the possibility that the current population is but a very small fraction of the initial population with the remnants of the destroyed clusters constituting presently a large fraction of the spheroid (bulge+halo) stellar population.

We compute, including a current state-of-the-art treatment of hydrodynamical processes, heating and cooling, a variety of cosmological models into the extreme nonlinear phase to enable comparisons with observations. First, we note the common, model independent results. All have a mean ($z=0$) temperature of $10^{4.5}-10^{5.5}$\Kel, set essentially by photoheating processes. Most gas is in one of two components: either at the photoheating floor $10^{4.5}$\Kel and primarily in low density regions or else shock heated to $10^5-10^6$\Kel and in regions of moderate overdensity (in caustics and near groups and clusters). It presents a major observational challenge to observationally detect this second, abundant component as it is neither an efficient radiator nor absorber. About 2% to 10% of the baryons cool and collapse into galaxies forming on caustics and migrating to clusters. About 1%-2% of baryons are in the very hot X-ray emitting gas near cluster cores, in good agreement with observations. These correspondances between the simulations and the real world imply that there is some significant truth to the underlying standard scenarios for the growth of structure. The differences among model predictions may help us find the path to the correct model. For COBE normalized models the most relevant differences concern epoch of structure formation. In the open variants having $\Omega=0.3$, with or without a cosmological constant, structure formation on galactic scales is well advanced at redshift z=5, and reionization occurs early. But if observations require models for which most galaxy formation occurs more recently than $z=2$, then the flat $\Omega=1$ models are to be preferred. The velocity dispersion on the $1h^{-1}$Mpc scale also provides a strong discriminant with, as expected, the $\Omega=1$ models giving a much higher

We model the $\lya$ forest in a spatially flat, CDM model with $\Omega=0.4$, with an Eulerian hydrodynamic simulation, and find that the intergalactic, photoionized gas collapses into sheet-like and filamentary structures which produce absorption lines like those observed in the $\lya$ forest. A typical filament is $\sim 1 h^{-1} \mpc$ long with thickness $\sim 50-100 h^{-1} \kpc$ (in proper units), and baryonic mass $\sim 10^{10} h^{-1} \msun$. The gas temperature increases with time as structures with larger velocities collapse gravitationally. The predicted distributions of column densities and b-parameters of the absorption lines agree with the observed ones, and their evolution is consistent with our model if the ionizing background has an approximately constant intensity between $z=2$ and $z=4$. We also predict the distribution of transmitted flux and its correlation along a spectrum and on parallel spectra, and the \heii flux decrement as a function of redshift. We predict a correlation length of $\sim 80 h^{-1}\kpc$ perpendicular to the line of sight for features in the Lyman alpha forest. In order to reproduce the observed number of lines and average flux transmission, the baryon content of the clouds may need to be significantly higher than in previous models because of the low densities we predict. If the background intensity $\jhi$ is at least that predicted from the observed quasars, $\Omega_b$ needs to be $\sim 0.025 h^{-2}$, higher than expected by light element nucleosynthesis; the model also predicts that most of the baryons at $z>2$ are in $\lya$ clouds, and that the rate at which the baryons move to more overdense regions is slow.

The observed baryon fraction and velocity--temperature relation in clusters of galaxies are compared with hydrodynamic simulations in two cosmological models : standard (Omega = 1) and a low-density flat (Omega=0.45 and \lambda=0.55) CDM models, normalized to the COBE background fluctuations. The observed properties of clusters include the velocity dispersion versus temperature relation, the gas mass versus total mass relation, and the gas mass fraction versus velocity dispersion relation. We find that, while both cosmological models reproduce well the shape of these observed functions, only low-density CDM can reproduce the observed amplitudes. The cluster gas mass fraction reflects approximately the baryon fraction in the models, with a slight anti-bias. Therefore, due to the low baryon density given by nucleosynthesis, Omega = 1 models produce too few baryons in clusters compared with observations. Scaling our results as a function of Omega, we find that a low-density CDM model, with Omega approximately 0.3 - 0.4, best reproduces the observed mean baryon fraction in clusters. The observed beta parameter of clusters, beta = sigma^2/(k T/mu m_p) = 0.94 \pm 0.08 discriminates less well between the models; it is consistent with that produced by low-density CDM (1.10 \pm 0.22), while it is slightly larger than expected but still consistent with Omega = 1 (0.70 \pm 0.14).

The topology of large scale structure is studied as a function of galaxy type using the genus statistic. In hydrodynamical cosmological CDM simulations, galaxies form on caustic surfaces (Zeldovich pancakes) then slowly drain onto filaments and clusters. The earliest forming galaxies in the simulations (defined as ``ellipticals") are thus seen at the present epoch preferentially in clusters (tending toward a meatball topology), while the latest forming galaxies (defined as ``spirals") are seen currently in a spongelike topology. The topology is measured by the genus (= number of ``donut" holes - number of isolated regions) of the smoothed density-contour surfaces. The measured genus curve for all galaxies as a function of density obeys approximately the theoretical curve expected for random-phase initial conditions, but the early forming elliptical galaxies show a shift toward a meatball topology relative to the late forming spirals. Simulations using standard biasing schemes fail to show such an effect. Large observational samples separated by galaxy type could be used to test for this effect.

We present a new treatment of two popular models for the growth of structure, examining the X-ray emission from hot gas with allowance for spectral line emission from various atomic species, primarily ``metals". The X-ray emission from the bright cluster sources is not significantly changed from prior work and, as noted earlier, shows the CDM$+\Lambda$ model (LCDM) to be consistent but the standard, COBE normalized model (SCDM) to be inconsistent with existing observations --- after allowance for still the considerable numerical modelling uncertainties. But we find one important new result. Radiation in the softer band 0.5-1.0keV is predominantly emitted by gas far from cluster centers (hence ``background"). This background emission dominates over the cluster emission below 1keV and observations of it should show clear spectral signatures indicating its origin. In particular the ``iron blend" should be seen prominantly in this spectral bin from cosmic background hot gas at high galactic latitudes and should show shadowing against the SMC indicating its extragalactic origin. Certain OVII lines also provide a signature of this gas which emits a spectrum characteristic of $10^{6.6\pm 0.6}$K gas. Recent ASCA observations of the X-ray background tentatively indicate the presence of component with exactly the spectral features we predict here.

Observational constraints guide one forcefully to examine models in which the matter density is substantially less than critical density. Particularly noteworthy are those which are consistent with inflation. For these models, microwave background anisotropy, large-scale structure measurements, direct measurements of the Hubble constant, $H_0$, and the closure parameter, $\Omega_{\rm Matter}$, ages of stars and a host of more minor facts are all consistent with a spatially flat model having significant cosmological constant $\Omega_{\Lambda} = 0.65 \pm 0.1$, $\Omega_{\rm Matter} = 1 - \Omega_{\Lambda}$ (in the form of ``cold dark matter") and a small tilt: $0.8< n < 1.2$. (Read comments at top of file to find out how to obtain figures.)

The properties of axisymmetric accretion flows of cold adiabatic gas with zero total energy in the vicinity of a Newtonian point mass are characterized by a single dimensionless parameter, the thickness of incoming flow. In the limit of thin accretion flows with vanishing thickness, we show that the governing equations become self-similar, involving no free parameters. We study numerically thin accretion flows with finite thickness as well as those with vanishing thickness. Mass elements of the incoming flow enter the computational regime as thin rings. In the case with finite thickness, after a transient period of initial adjustment, an almost steady-state accretion shock with a small oscillation amplitude forms, confirming the previous work by Molteni, Lanzafame, \& Chakrabarti (1994). The gas in the region of vorticity between the funnel wall and the accretion shock follows closed streamlines, forming a torus. This torus, in turn, behaves as an effective barrier to the incoming flow and supports the accretion shock which reflects the incoming gas away from the equatorial plane. The postshock flow, which is further accelerated by the pressure gradient behind the shock, goes through a second shock which then reflects the flow away from the symmetry axis to form a conical outgoing wind. As the thickness of the inflowing layer decreases (or if the ratio of the half thickness to the distance to the funnel wall along the equatorial plan is smaller than $\sim0.1$), the flow becomes unstable. In the case with vanishing thickness, the accretion shock formed to stop the incoming flow behind the funnel wall oscillates quasi-periodically with an amplitude comparable to the thickness. The structure between the funnel wall and the accretion shock is destroyed as the shock moves inwards toward the central mass and re-generated

If the large-scale galaxy distribution is filamentary, as suggested by some observations and recent hydrodynamical simulations, then lengthwise views of filaments will apparently produce compact groups (CGs) that are in reality stretched out along the line of sight. This possibility has been advocated recently by Hernquist, Katz \& Weinberg (1994). Here, we propose a test for this hypothesis using X--ray emission from CGs. The observable quantity $Q \equiv L_{x} a_{p}^{3}/L_{g}^2 T_{x}^{1/2}$ should be proportional to the axis ratio of the group, $a/c$, where $a$ and $c$ are the long and short axis of a prolate distribution, $a_p$ is the radius of the group projected onto the sky, $L_{x}$ is the bolometric X-ray luminosity, $L_{g}$ is the group blue luminosity, and $T_{x}$ is the gas temperature. We find that the distribution of $Q$ is consistent with the notion that many spiral-rich CGs with unusually small values of ($a/c)$ are frauds, i.e. that the values of $Q$ are anomalously small. An alternative possibility is that CGs are gas-poor relative to rich clusters; however, this can be tested using the Sunyaev-Zeldovich effect. If the groups have a normal ratio of gas-to-total mass, but are simply stretched out along the line of sight, a Sunyaev-Zeldovich signal should be detectable.

An adaptation of the Particle-Particle/Particle-Mesh (P3M) code to the special purpose hardware GRAPE is presented. The short range force is calculated by a four chip GRAPE-3A board, while the rest of the calculation is performed on a Sun Sparc 10/51 workstation. The limited precision of the GRAPE hardware and algorithm constraints introduce stochastic errors of the order of a few percent in the gravitational forces. Tests of this new P3MG3A code show that it is a robust tool for cosmological simulations. The code currently achieves a peak efficiency of one third the speed of the vectorized P3M code on a Cray C-90 and significant improvements are planned in the near future. Special purpose computers like GRAPE are therefore an attractive alternative to supercomputers for numerical cosmology.

Gravitational lensing provides a strict test of cosmogonic models. Detailed numerical propagation of light rays through a universe having a distribution of inhomogeneities derived from the standard CDM (cold dark matter) scenario with the aid of massive, fully nonlinear computer simulations is used to test the model. It predicts that we should have seen far more widely split quasar images than have been found. These and other inconsistencies rule out the COBE normalized CDM model with $\Omega=1$ and $H_0=50$km/s/Mpc; but variants of this model might be constructed, which could pass the stringent tests provided by strong gravitational lensing.

If gravitational clustering is a hierarchical process, the present large-scale structure of the galaxy distribution implies that structures on smaller scales must have formed at high redshift. We simulate the formation of small-scale structure (average cell mass: $\Delta \bar m_b=10^{4.2}$M$_\odot$) and the evolution of photoionized gas, in the specific case of a CDM model with a cosmological constant. The photoionized gas has a natural minimal scale of collapse, the Jeans scale ($m_{b,J}\simeq 10^{9}$M$_\odot$). We find that low column density ($\nhi \le 10^{14}\cm^{-2}$) lines originate in regions resembling Zel'dovich pancakes, where gas with overdensities in the range $3 - 30$ is enclosed by two shocks but is typically re-expanding at approximately the Hubble velocity. However, higher column density ($\nhi \ge 10^{15}\cm^{-2}$) lines stem from more overdense regions where the shocked gas is cooling. We show that this model can probably account for the observed number of lines, their distribution in column density and b-parameters, as well as the cloud physical sizes as observed in gravitationally lensed quasars. We find a redshift evolution that is

Lauer \& Postman (LP) observe that all Abell clusters with redshifts less than 15,000\kms\ appear to be participating in a bulk flow of 689 km s$^{-1}$ with respect to the Cosmic Microwave Background. We find this result difficult to reconcile with all popular models for large-scale structure formation that assume Gaussian initial conditions. This conclusion is based on Monte-Carlo realizations of the LP data, drawn from large Particle-Mesh $N$-body simulations. We have taken special care to treat properly the longest-wavelength components of the power spectra. Bulk flows with amplitude as large as that reported by LP are not uncommon in the Monte-Carlo datasets. However, the $\chi^2$ of the observed bulk flow, taking into account the anisotropy of the error ellipsoid, is much more difficult to match in the simulations. The models examined are ruled out at confidence levels between 94\% and 98\%. Any model that has {\it intrinsic} flows of less than 480\kms\ on the scales probed by LP scales can be ruled out at a similar level.

We present a detailed comparison of the simulation results of various cosmological hydrodynamic codes. Starting with identical initial conditions based on the Cold Dark Matter scenario for the growth of structure, we integrate from redshift $z=20$ to $z=0$ to determine the physical state within a representative volume of size $L^3$ where $L=64 h^{-1} {\rm Mpc}$. Five independent codes are compared: three of them Eulerian mesh based and two variants of the Smooth Particle Hydrodynamics "SPH" Lagrangian approach. The Eulerian codes were run at $N^3=(32^3,~64^3,~128^3,~{\rm and},~256^3)$ cells, the SPH codes at $N^3= 32^3$ and $64^3$ particles. Results were then rebinned to a $16^3$ grid with the expectation that the rebinned data should converge, by all techniques, to a common and correct result as $N \rightarrow \infty$. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainties in even primitive quantities such as $\langle T \rangle$, $\langle \rho^2\rangle^{1/2}$ persists at the 3\%-17\% level after completion of very large simulations. The two SPH codes and the two shock capturing Eulerian codes achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high density, high temperature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by $\rho^2$) is only 4.5\%. Overall the comparison allows us to better estimate errors, it points to ways of improving this current generation of hydrodynamic codes and of suiting their use to problems which exploit their individually best features.

We exmaine the distribution of hot gas in a standard CDM model of the universe using high resolution hydrodynamic simulations. Adopting standard parameters determined from COBE and light element nucleosynthesis, $\sigma_8=1.05$, $\Omega_b=0.06$ and assuming $h=0.5$, we find the X-ray emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution and estimated sizes as well as the evolution of these quantities with redshift. This standard CDM model, normalized to COBE, produces approximately 5 times too much emission from clusters having $L_x>10^{43}$erg/s, a not unexpected result. If all other parameters were unchanged, we would expect adequate agreement for $\sigma_8=0.6$. This provides a new and independent argument for lower small scale power than standard CDM at the $8h^{-1}$Mpc scale. The background radiation field at 1keV due to clusters in this model is approximately $1/3$ of the observed background which, after correction for numerical effects, again indicates approximately 5 times too much emission and the appropriateness of $\sigma_8=0.6$. If we had used the observed ratio of gas to total mass in clusters, rather than basing the mean density on light element nucleosynthesis, then the computed luminosity of each cluster would have increased still further, by a factor of approximately ten. Examining the ratio of gas to total mass in the clusters (which we find to be anti-biased by a factor of approximately 0.6), normalized to $\Omega_b h^2=0.015$, and comparing to observations, we conclude, in agreement with S. White, that the cluster observations argue for an open universe.

A new, three-dimensional, shock capturing, TVD hydrodynamic code is utilized to determine the distribution of hot gas in a CDM$+\Lambda$ model universe. This open model succeeds in matching local observations of clusters in contrast to the standard $\Omega=1$, CDM model, which fails. It predicts an order of magnitude decline in the number density of bright ($h\nu = 2-10$keV) clusters from $z=0$ to $z=2$ in contrast to a slight increase in the number density for standard $\Omega=1$, CDM model. This COBE-normalized CDM$+\Lambda$ model produces approximately the same number of X-ray clusters having $L_x>10^{43}$erg/s as observed. The number density of bright clusters increases to $z\sim 0.2-0.5$ and declines, but the luminosity per typical cluster decreases monotonically with redshift, with the result that the number density of bright clusters shows a broad peak near $z=0.5$, and then a rapid decline as $z\rightarrow 3$. The most interesting point which we find is that the temperatures of clusters in this model freeze out at later times ($z\le 0.3$), while previously we found in the CDM model that there was a steep increase during the same interval of redshift. The trend should be detectable even with a relatively ``soft" X-ray instrument such as ROSAT, providing a powerful discriminant between $\Omega=1$ and $\Omega<1$ models. Examining the ratio of gas-to-total mass in the clusters, we find a slight antibias ($b=0.9$ or $({\Omega_{gas}\over\Omega_{tot}})_{cl} =0.083\pm 0.007$), which is consistent with observations [$({\Omega_{gas}\over\Omega_{tot}})_{obs}=0.097\pm 0.019$ for the Coma cluster for the given value of $h$, \cf, White 1991].

We compute the evolution of spatially flat, mixed cold and hot dark matter (``MDM") models containing both baryonic matter and two kinds of dark matter. The mean final Zeldovich-Sunyaev $y$ parameter is estimated to be $\bar y=(5.4\pm 2.7)\times 10^{-7}$, with a rms fluctuation of approximately $\bar{\delta y}=(6.0\pm 3.0)\times 10^{-7}$ on arc minute scales. The rate of galaxy formation peaks at an even later epoch ($z\sim 0.3$) than in the standard ($\Omega=1$, $\sigma_8=0.67$) CDM model ($z\sim 0.5$) and, at a redshift of z=4 is nearly a factor of a hundred lower than for the CDM model with the same value of $\sigma_8$. The model yields a cluster mass function which is about a factor of 4 higher than observed but a cluster-cluster correlation length lower by a factor of 2 than what is observed but both are closer to observations than in the COBE normalized CDM model. The one dimensional pairwise velocity dispersion is $605\pm 8$km/s at $1h^{-1}$ separation, lower than that of the CDM model normalized to COBE, but still significantly higher than observations (Davis \& Peebles 1983). A plausible velocity bias $b_v=0.8\pm 0.1$ on this scale will reduce but not remove the discrepancy. Overall, the model is similar both on large and intermediate scales to the standard CDM model normalized to the same value of $\sigma_8$, but the problem with regard to late formation of galaxies is more severe in this model than in that CDM model. Adding hot dark matter significantly improves the ability of COBE normalized CDM scenario to fit existing observations, but the model is in fact not as good as the CDM model with the same $\sigma_8$ and is still probably unsatisfactory with regard to several critical tests.

We challenge a widely accepted assumption of observational cosmology: that successful reconstruction of observed galaxy density fields from measured galaxy velocity fields (or vice versa), using the methods of gravitational instability theory, implies that the observed large-scale structures and large-scale flows were produced by the action of gravity. This assumption is false, in that there exist non-gravitational theories that pass the reconstruction tests and gravitational theories with certain forms of biased galaxy formation that fail them. Gravitational instability theory predicts specific correlations between large-scale velocity and mass density fields, but the same correlations arise in any model where (a) structures in the galaxy distribution grow from homogeneous initial conditions in a way that satisfies the continuity equation, and (b) the present-day velocity field is irrotational and proportional to the time-averaged velocity field. We demonstrate these assertions using analytical arguments and N-body simulations of gravitational and non-gravitational models. We also show examples of gravitational and non-gravitational models that {\it fail} reconstruction tests because galaxy formation is modulated (``biased'') in a way that violates the continuity equation. We discuss the relation between the value of $\Omega$ inferred from velocity-density comparisons and the true cosmological value.

The large-scale bulk flow implied by the measurements of distances to brightest cluster members within 15,000 km/s by Lauer and Postman has the potential to put strong constraints on, or to even rule out, various models of large-scale structure. Using PM simulations 400$h^{-1}$ Mpc on a side of five cosmological models (Standard CDM, Tilted CDM, HDM, $\Omega = 0.3$ CDM, and PBI), we present detailed Monte-Carlo simulations of the peculiar velocity data, which we use to assess the probability that any given model can explain the data. These models are ruled out at confidence levels ranging from 94\% (HDM) to 97\% (PBI). However, these results may be affected by the finite size of the simulation volume. If the observed bulk flow remains unchanged when calculated using a refined distance indicator with smaller scatter, all these models will be ruled out at high confidence levels.