$\Lambda$-Warm Dark Matter (WDM) has been proposed as alternative scenario to $\Lambda$ cold dark matter (CDM), motivated by discrepancies at the scale of dwarf galaxies, with less small-scale power and realized by collisionless particles with energies in the range $1-3$ keV. We present a new approach to constrain the viability of such WDM models using star formation histories of the dwarf spheroidal galaxies (dSphs) in the Local Group. We compare their high time-resolution star formation histories (SFHs) obtained with HST-based color magnitude diagrams with the range of possible collapse redshifts of their dark matter halos expected in CDM and in different WDM scenarios. The collapse redshift is inferred after determining a plausible infall mass of the subhalo. This is based on the current mass of individual dwarf inferred from stellar kinematics combined with results of cosmological simulations providing information on the subhalo evolution. Since WDM subhalos close to the filtering mass scale form significantly later than CDM, we show that they are in the first place difficult to reconcile with a truncation of star formation occurring as early as $z\geq 3$. The Ultra-Faint Dwarfs (UFDs) provide the most stringent constraints. Using 6 UFDs with the best determination of the SFHs, we show that we can exclude a 1 keV warm particle to a 2-$\sigma$ confidence interval consistently with other methods reported in the literature. For some objects the $2$ keV model is also excluded. We discuss the various caveats of the method, most notably the low number of dwarfs with accurately determined star formation histories and the uncertainties in the determination of the infall mass of the subhalos. Our preliminary analysis serves as a pathfinder for future investigations that will combine upcoming accurate SFHs for more local dSphs with direct analysis of WDM cosmological simulations with baryons.

Using a new, high-resolution cosmological hydrodynamic simulation of a Milky Way-type (MW-type) galaxy, we explore how a merger-rich assembly history affects the mass budget of the central supermassive black hole (SMBH). We examine a MW-mass halo at the present epoch whose evolution is characterized by several major mergers to isolate the importance of merger history on black hole accretion. This study is an extension of Bellovary et. al. 2013, which analyzed the accretion of high mass, high redshift galaxies and their central black holes, and found that the gas content of the central black hole reflects what is accreted by the host galaxy halo. In this study, we find that a merger-rich galaxy will have a central SMBH preferentially fed by merger gas. Moreover, we find that nearly 30$\%$ of the accreted mass budget of the SMBH enters the galaxy through the two major mergers in its history, which may account for the increase of merger-gas fueling the SMBH. Through an investigation of the angular momentum of the gas entering the host and its SMBH, we determine that merger gas enters the galaxy with lower angular momentum compared to smooth accretion, partially accounting for the preferential fueling witnessed in the SMBH. In addition, the presence of mergers, particularly major mergers, also helps funnel low angular momentum gas more readily to the center of the galaxy. Our results imply that galaxy mergers play an important role in feeding the SMBH in MW-type galaxies with merger-rich histories. Our results imply that galaxy mergers play an important role in feeding the SMBH in MW-type galaxies with merger-rich histories.

We examine the scalings of galactic outflows with halo mass across a suite of twenty high-resolution cosmological zoom galaxy simulations covering halo masses from 10^9.5 - 10^12 M_sun. These simulations self-consistently generate outflows from the available supernova energy in a manner that successfully reproduces key galaxy observables including the stellar mass-halo mass, Tully-Fisher, and mass-metallicity relations. We quantify the importance of ejective feedback to setting the stellar mass relative to the efficiency of gas accretion and star formation. Ejective feedback is increasingly important as galaxy mass decreases; we find an effective mass loading factor that scales as v_circ^(-2.2), with an amplitude and shape that is invariant with redshift. These scalings are consistent with analytic models for energy-driven wind, based solely on the halo potential. Recycling is common: about half the outflow mass across all galaxy masses is later re-accreted. The recycling timescale is typically about 1 Gyr, virtually independent of halo mass. Recycled material is re-accreted farther out in the disk and with typically about 2-3 times more angular momentum. These results elucidate and quantify how the baryon cycle plausibly regulates star formation and alters the angular momentum distribution of disk material across the halo mass range where most of cosmic star formation occurs.

We present a self-consistent prediction from a large-scale cosmological simulation for the population of `wandering' supermassive black holes (SMBHs) of mass greater than $10^6$ M$_{\odot}$ on long-lived, kpc-scale orbits within Milky Way (MW)-mass galaxies. We extract a sample of MW-mass halos from the Romulus25 cosmological simulation (Tremmel et al. 2017), which is uniquely able to capture the orbital evolution of SMBHs during and following galaxy mergers. We predict that such halos, regardless of recent merger history or morphology, host an average of $5.1 \pm 3.3$ SMBHs, including their central black hole, within 10 kpc from the galactic center and an average of $12.2 \pm 8.4$ SMBHs total within their virial radius, not counting those in satellite halos. Wandering SMBHs exist within their host galaxies for several Gyrs, often accreted by their host halo in the early Universe. We find, with $>4\sigma$ significance, that wandering SMBHs are preferentially found outside of galactic disks.

The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consistent and predict strong radial anisotropy (<beta>~0.7) for Galactocentric radii greater than 10 kpc. Previous observations of beta for the MW's stellar halo claim a detection of an isotropic or tangential "dip" at r~20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of "dips" in beta. We find dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z~1 can result in a present day low (isotropic to tangential) value of beta over a wide range of radii and angular expanse. While all of these mechanisms are plausible drivers for the beta dip observed in the MW, in the simulations, each mechanism has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of beta(r) is required for measuring the mass of the MW halo, we note significant transient dips in beta could cause an overestimate of the halo's mass when using spherical Jeans equation modeling.

We present the first self-consistent prediction for the distribution of formation timescales for close Supermassive Black Hole (SMBH) pairs following galaxy mergers. Using ROMULUS25, the first large-scale cosmological simulation to accurately track the orbital evolution of SMBHs within their host galaxies down to sub-kpc scales, we predict an average formation rate density of close SMBH pairs of 0.013 cMpc^-3 Gyr^-1. We find that it is relatively rare for galaxy mergers to result in the formation of close SMBH pairs with sub-kpc separation and those that do form are often the result of Gyrs of orbital evolution following the galaxy merger. The likelihood and timescale to form a close SMBH pair depends strongly on the mass ratio of the merging galaxies, as well as the presence of dense stellar cores. Low stellar mass ratio mergers with galaxies that lack a dense stellar core are more likely to become tidally disrupted and deposit their SMBH at large radii without any stellar core to aid in their orbital decay, resulting in a population of long-lived 'wandering' SMBHs. Conversely, SMBHs in galaxies that remain embedded within a stellar core form close pairs in much shorter timescales on average. This timescale is a crucial, though often ignored or very simplified, ingredient to models predicting SMBH mergers rates and the connection between SMBH and star formation activity.

Within a Lambda Cold Dark Matter (LCDM) scenario, we use high resolution cosmological simulations spanning over four orders of magnitude in galaxy mass to understand the deficit of dwarf galaxies in observed velocity functions. We measure velocities in as similar a way as possible to observations, including generating mock HI data cubes for our simulated galaxies. We demonstrate that this apples-to-apples comparison yields an "observed" velocity function in agreement with observations, reconciling the large number of low-mass halos expected in a LCDM cosmological model with the low number of observed dwarfs at a given velocity. We then explore the source of the discrepancy between observations and theory, and conclude that the dearth of observed dwarf galaxies is primarily explained by two effects. The first effect is that galactic rotational velocities derived from the HI linewidth severely underestimate the maximum halo velocity. The second effect is that a large fraction of halos at the lowest masses are too faint to be detected by current galaxy surveys. We find that cored dark matter density profiles can contribute to the lower observed velocity of galaxies, but only for galaxies in which the velocity is measured interior to the size of the core (~3 kpc).

We present a novel implementation of supermassive black hole (SMBH) formation, dynamics, and accretion in the massively parallel tree+SPH code, ChaNGa. This approach improves the modeling of SMBHs in fully cosmological simulations, allowing for a more de- tailed analysis of SMBH-galaxy co-evolution throughout cosmic time. Our scheme includes novel, physically motivated models for SMBH formation, dynamics and sinking timescales within galaxies, and SMBH accretion of rotationally supported gas. The sub-grid parameters that regulate star formation (SF) and feedback from SMBHs and SNe are optimized against a comprehensive set of z = 0 galaxy scaling relations using a novel, multi-dimensional parameter search. We have incorporated our new SMBH implementation and parameter optimization into a new set of high resolution, large-scale cosmological simulations called Romulus. We present initial results from our flagship simulation, Romulus25, showing that our SMBH model results in SF efficiency, SMBH masses, and global SF and SMBH accretion histories at high redshift that are consistent with observations. We discuss the importance of SMBH physics in shaping the evolution of massive galaxies and show how SMBH feedback is much more effective at regulating star formation compared to SNe feedback in this regime. Further, we show how each aspect of our SMBH model impacts this evolution compared to more common approaches. Finally, we present a science application of this scheme studying the properties and time evolution of an example dual AGN system, highlighting how our approach allows simulations to better study galaxy interactions and SMBH mergers in the context of galaxy-BH co-evolution.

Using a cosmological $N$-body simulation, we investigate the origin and distribution of stars in the intracluster light (ICL) of a Fornax-like cluster. In a dark matter only simulation we identify a halo which, at $z=0$, has $M_200 \simeq 4.1 \times 10^{13}M_{sun}$ and $r_{200} = 700kpc$, and replace infalling subhalos with models that include spheroid and disc components. As they fall into the cluster, the stars in some of these galaxies are stripped from their hosts, and form the ICL. We consider the separate contributions to the ICL from stars which originate in the haloes and the discs of the galaxies. We find that disc ICL stars are more centrally concentrated than halo ICL stars. The majority of the disc ICL stars are associated with one initially disc-dominated galaxy that falls to the centre of the cluster and is heavily disrupted, producing part of the cD galaxy. At radial distances greater than 200kpc, well beyond the stellar envelope of the cD galaxy, stars formerly from the stellar haloes of galaxies dominate the ICL. Therefore at large distances, the ICL population is dominated by older stars.

We predict the stellar mass-halo mass (SMHM) relationship for dwarf galaxies and their satellites residing in halos down to M$_{halo} =$ 10$^7$ M$_{\odot}$ with 10$^4$ M$_{\odot} <$ M$_{star}$($z=0$) $< 10^8$ M$_{\odot}$, and quantify the predicted scatter in the relation at the low mass end, using cosmological simulations. The galaxies were drawn from a cosmological simulation of dwarf galaxies, run with the N-body + SPH code, ChaNGA, at a high resolution of 60 pc. For M$_{halo} > 10^9$ M$_{\odot}$, the simulated SMHM relationship agrees with literature determinations, including exhibiting a small scatter. However, the scatter in the SMHM relation increases dramatically for lower-mass halos. We find that some of this scatter is due to {\em dark dwarfs}, halos devoid of stars. However, even when only considering well-resolved halos that contain a stellar population, the scatter in stellar mass reaches nearly 1 dex for M$_{halo}$($z=0$) 10$^7$ M$_{\odot}$. Much of this scatter is due to including satellites of the dwarf galaxies that have had their halo masses reduced through tidal stripping. The fraction of dark dwarfs (those that contain no stars) increases substantially with decreasing halo mass. When these dark halos are considered, the true scatter in the SMHM at low masses is even larger. At the faintest end of the SMHM relation probed by our simulations, a galaxy cannot be assigned a unique halo mass based solely on its luminosity. We provide a formula to stochastically populate low-mass halos following our results. Our predicted large scatter at low halo masses increases the slope of the resulting stellar mass function on the ultra-faint dwarf galaxy scales currently being probed by such surveys as the Dark Energy Survey or the Hyper-Suprime Cam Subaru Strategic Program, and in the future by the Large Synoptic Survey Telescope.

We explore for the first time the effect of self-interacting dark matter (SIDM) on the dark matter (DM) and baryonic distribution in massive galaxies formed in hydrodynamical cosmological simulations, including explicit baryonic physics treatment. A novel implementation of Super-Massive Black Hole (SMBH) formation and evolution is used, as in Tremmel et al.(2015, 2016), allowing to explicitly follow SMBH dynamics at the center of galaxies. A high SIDM constant cross-section is chosen, $\sigma$=10 $\rm cm^2/gr$, to amplify differences from CDM models. Milky Way-like galaxies form a shallower DM density profile in SIDM than they do in CDM, with differences already at 20 kpc scales. This demonstrates that even for the most massive spirals the effect of SIDM dominates over the adiabatic contraction due to baryons. Strikingly, the dynamics of SMBHs differs in the SIDM and reference CDM case. SMBHs in massive spirals have sunk to the centre of their host galaxy in both the SIDM and CDM run, while in less massive galaxies about 80$\%$ of the SMBH population is off-centered in the SIDM case, as opposed to the CDM case in which $\sim$90$\%$ of SMBHs have reached their host's centre. SMBHs are found as far as $\sim$9 kpc away from the centre of their host SIDM galaxy. This difference is due to the increased dynamical friction timescale caused by the lower DM density in SIDM galaxies compared to CDM, resulting in 'core stalling'. This pilot work highlights the importance of simulating in a full hydrodynamical context different DM models combined to SMBH physics to study their influence on galaxy formation.

The sources that reionized the universe are still unknown, but likely candidates are faint but numerous galaxies. In this paper we present results from running a high resolution, uniform volume simulation, the Vulcan, to predict the number densities of undetectable, faint galaxies and their escape fractions of ionizing radiation, $f_\mathrm{esc}$, during reionization. Our approach combines a high spatial resolution, a realistic treatment of feedback and hydro processes, a strict threshold for minimum number of resolution elements per galaxy, and a converged measurement of $f_\mathrm{esc}$. We calibrate our physical model using a novel approach to create realistic galaxies at z=0, so the simulation is predictive at high redshifts. With this approach we can (1) robustly predict the evolution of the galaxy UV luminosity function at faint magnitudes down to $M_\mathrm{UV}$~-15, two magnitudes fainter than observations, and (2) estimate $f_\mathrm{esc}$ over a large range of galaxy masses based on the detailed stellar and gas distributions in resolved galaxies. We find steep faint end slopes, implying high number densities of faint galaxies, and the dependence of $f_\mathrm{esc}$ on the UV magnitude of a galaxy, given by the power-law: log $f_\mathrm{esc} = (0.51 \pm 0.04)M_\mathrm{UV} + 7.3 \pm 0.8$, with the faint population having $f_\mathrm{esc}$~35%. Convolving the UV luminosity function with $f_\mathrm{esc}$($M_\mathrm{UV}$), we find an ionizing emissivity that is (1) dominated by the faintest galaxies and (2) reionizes the universe at the appropriate rate, consistent with observational constraints of the ionizing emissivity and the optical depth to the decoupling surface tau_es, without the need for additional sources of ionizing radiation.

We confirm that the object DDO216-A1 is a substantial globular cluster at the center of Local Group galaxy DDO216 (the Pegasus dwarf irregular), using Hubble Space Telescope ACS imaging. By fitting isochrones, we find the cluster metallicity to be [M/H] = -1.6 +/-0.2, for reddening E(B-V) = 0.16 +/-0.02; the best-fit age is 12.3 +/-0.8 Gyr. There are ~30 RR Lyrae variables in the cluster; the magnitude of the fundamental mode pulsators gives a distance modulus of 24.77 +/-0.08 - identical to the host galaxy. The ratio of overtone to fundamental mode variables and their mean periods make DDO216-A1 an Oosterhoff Type I cluster. We find an I-band central surface brightness 20.85 +/-0.17 F814W mag per square arcsecond, a half-light radius of 3.1 arcsec (13.4 pc), and an absolute magnitude M814 = -7.90 +/-0.16 (approximately 10^5 solar masses). King models fit to the cluster give the core radius and concentration index, r_c = 2.1" +/-0.9" and c = 1.24 +/-0.39. The cluster is an "extended" cluster somewhat typical of some dwarf galaxies and the outer halo of the Milky Way. The cluster is projected <30 pc south of the center of DDO216, unusually central compared to most dwarf galaxy globular clusters. Analytical models of dynamical friction and tidal destruction suggest that it probably formed at a larger distance, up to ~1 kpc, and migrated inward. DDO216 has an unexceptional cluster specific frequency, S_N = 10. DDO216 is the lowest-luminosity Local Group galaxy to host a 10^5 solar mass globular cluster, and the only transition-type (dSph/dIrr) in the Local Group with a globular.

We show how the interplay between active galactic nuclei (AGN) and merger history determines whether a galaxy quenches star formation at high redshift. We first simulate, in a full cosmological context, a galaxy of total dynamical mass $10^{12}\,M_{\odot}$ at $z=2$. Then we systematically alter the accretion history of the galaxy by minimally changing the linear overdensity in the initial conditions. This "genetic modification" approach allows the generation of three sets of $\Lambda$CDM initial conditions leading to maximum merger ratios of 1:10, 1:5 and 2:3 respectively. The changes leave the final halo mass, large scale structure and local environment unchanged, providing a controlled numerical experiment. Interaction between the AGN physics and mergers in the three cases lead respectively to a star-forming, temporarily-quenched and permanently-quenched galaxy. However the differences do not primarily lie in the black hole accretion rates, but in the kinetic effects of the merger: the galaxy is resilient against AGN feedback unless its gaseous disk is first disrupted. Typical accretion rates are comparable in the three cases, falling below $0.1\,M_{\odot}$ yr$^{-1}$, equivalent to around $2\%$ of the Eddington rate or $10^{-3}$ times the pre-quenching star formation rate, in agreement with observations. This low level of black hole accretion can be sustained even when there is insufficient dense cold gas for star formation. Conversely, supernova feedback is too distributed to generate outflows in high-mass systems, and cannot maintain quenching over periods longer than the halo gas cooling time.

Cosmological simulations of galaxy formation can produce present-day galaxies with a large range of assembly and star formation histories. A detailed study of the metallicity evolution and star formation history of such simulations can assist in predicting LIGO-detectable compact object binary mergers. Recent simulations of compact binary evolution suggest the compact object merger rate depends sensitively on the progenitor's metallicity. Rare low-metallicity star formation during galaxy assembly can produce more detected compact binaries than typical star formation. Using detailed simulations of galaxy and chemical evolution, we determine how sensitively the compact binary populations of galaxies with similar present-day appearance depend on the details of their assembly. We also demonstrate by concrete example the extent to which dwarf galaxies overabundantly produce compact binary mergers, particularly binary black holes, relative to more massive galaxies. We discuss the implications for transient multimessenger astronomy with compact binary sources.

We study the local dark matter velocity distribution in simulated Milky Way-mass galaxies, generated at high resolution with both dark matter and baryons. We find that the dark matter in the Solar neighborhood is influenced appreciably by the inclusion of baryons, increasing the speed of dark matter particles compared to dark matter-only simulations. The gravitational potential due to the presence of a baryonic disk increases the amount of high velocity dark matter, resulting in velocity distributions which are more similar to the Maxwellian Standard Halo Model than predicted from dark matter-only simulations. Further, the velocity structures present in baryonic simulations possess a greater diversity than expected from dark matter-only simulation. We show the impact on the direct detection experiments LUX, DAMA/Libra, and CoGeNT using our simulated velocity distributions, and explore how resolution and halo mass within the Milky Way's estimated mass range impact the results. A Maxwellian fit to the velocity distribution tends to overpredict the amount of dark matter in the high velocity tail, even with baryons, and thus leads to overly optimistic direct detection bounds on models which are dependent on this region of phase space for an experimental signal. Our work further demonstrates that it is critical to transform simulated velocity distributions to the lab frame of reference, due to the fact that velocity structure in the Solar neighborhood appears when baryons are included. There is more velocity structure present when baryons are included than in dark matter-only simulations. Even when baryons are included, the importance of the velocity structure is not as apparent in the Galactic frame of reference as in the Earth frame.

We present a new zoom-in hydrodynamical simulation, "Erisbh", which follows the cosmological evolution and feedback effects of a supermassive black hole at the center of a Milky Way-type galaxy. ErisBH shares the same initial conditions, resolution, recipes of gas cooling, star formation and feedback, as the close Milky Way-analog "Eris", but it also includes prescriptions for the formation, growth and feedback of supermassive black holes. We find that the galaxy's central black hole grows mainly through mergers with other black holes coming from infalling satellite galaxies. The growth by gas accretion is minimal because very little gas reaches the sub-kiloparsec scales. The final black hole is, at z=0, about 2.6 million solar masses and it sits closely to the position of SgrA* on the MBH-MBulge and MBH-sigma planes, in a location consistent with what observed for pseudobulges. Given the limited growth due to gas accretion, we argue that the mass of the central black hole should be above 10^5 solar masses already at z~8. The effect of AGN feedback on the host galaxy is limited to the very central few hundreds of parsecs. Despite being weak, AGN feedback seems to be responsible for the limited growth of the central bulge with respect to the original Eris, which results in a significantly flatter rotation curve in the inner few kiloparsecs. Moreover, the disk of ErisBH is more prone to instabilities, as its bulge is smaller and its disk larger then Eris. As a result, the disk of ErisBH undergoes a stronger dynamical evolution relative to Eris and around z=0.3 a weak bar grows into a strong bar of a few disk scale lengths in size. The bar triggers a burst of star formation in the inner few hundred parsecs, provides a modest amount of new fuel to the central black hole, and causes the bulge of ErisBH to have, by z=0, a box/peanut morphology.(Abridged)

Supermassive black hole dynamics during galaxy mergers is crucial in determining the rate of black hole mergers and cosmic black hole growth. As simulations achieve higher resolution, it becomes important to assess whether the black hole dynamics is influenced by the treatment of the interstellar medium in different simulation codes. We here compare simulations of black hole growth in galaxy mergers with two codes: the Smoothed Particle Hydrodynamics code Gasoline, and the Adaptive Mesh Refinement code Ramses. We seek to identify predictions of these models that are robust despite differences in hydrodynamic methods and implementations of sub-grid physics. We find that the general behavior is consistent between codes. Black hole accretion is minimal while the galaxies are well-separated (and even as they "fly-by" within 10 kpc at first pericenter). At late stages, when the galaxies pass within a few kpc, tidal torques drive nuclear gas inflow that triggers bursts of black hole accretion accompanied by star formation. We also note quantitative discrepancies that are model-dependent: our Ramses simulations show less star formation and black hole growth, and a smoother gas distribution with larger clumps and filaments, than our Gasoline simulations. We attribute these differences primarily to the sub-grid models for black hole fueling and feedback and gas thermodynamics. The main conclusion is that differences exist quantitatively between codes, and this should be kept in mind when making comparisons with observations. However, reassuringly, both codes capture the same dynamical behaviors in terms of triggering of black hole accretion, star formation, and black hole dynamics.

We use cosmological hydrodynamic simulations to consistently compare the assembly of dwarf galaxies in both $\Lambda$ dominated, Cold (CDM) and Self--Interacting (SIDM) dark matter models. The SIDM model adopts a constant cross section of 2 $cm^{2}/g$, a relatively large value to maximize its effects. These are the first SIDM simulations that are combined with a description of stellar feedback that naturally drives potential fluctuations able to create dark matter cores. Remarkably, SIDM fails to significantly lower the central dark matter density at halo peak velocities V$_{max}$ $<$ 30 Km/s. This is due to the fact that the central regions of very low--mass field halos have relatively low central velocity dispersion and densities, leading to time scales for SIDM collisions greater than a Hubble time. CDM halos with V$_{max}$ $<$ 30 km/s have inefficient star formation, and hence weak supernova feedback. At a fixed 2 cm2/g SIDM cross section, the DM content of very low mass CDM and SIDM halos differs by no more than a factor of two within 100-200pc. At larger halo masses ($\sim$ 10$^{10}$ solar masses), the introduction of baryonic processes creates field dwarf galaxies with dark matter cores and central DM$+$baryon distributions that are effectively indistinguishable between CDM and SIDM. Both models are in broad agreement with observed Local Group field galaxies across the range of masses explored. To significantly differentiate SIDM from CDM at the scale of faint dwarf galaxies, a velocity dependent cross section that rapidly increases to values larger than 2 $cm^{2}/g$ for halos with V$_{max}$ < 25-30 Km/s needs to be introduced.

We introduce a force correction term to better model the dynamical friction (DF) experienced by a supermassive black hole (SMBH) as it orbits within its host galaxy. This new approach accurately follows the orbital decay of a SMBH and drastically improves over commonly used advection methods. The force correction introduced here naturally scales with the force resolution of the simulation and converges as resolution is increased. In controlled experiments we show how the orbital decay of the SMBH closely follows analytical predictions when particle masses are significantly smaller than that of the SMBH. In a cosmological simulation of the assembly of a small galaxy, we show how our method allows for realistic black hole orbits. This approach overcomes the limitations of the advection scheme, where black holes are rapidly and artificially pushed toward the halo center and then forced to merge, regardless of their orbits. We find that SMBHs from merging dwarf galaxies can spend significant time away from the center of the remnant galaxy. Improving the modeling of SMBH orbital decay will help in making robust predictions of the growth, detectability, and merger rates of SMBHs, especially at low galaxy masses or at high redshift.

We use a suite of hydrodynamical simulations of galaxy mergers to compare star formation rate (SFR) and black hole accretion rate (BHAR) for galaxies before the interaction ('stochastic' phase), during the `merger' proper, lasting ~0.2-0.3 Gyr, and in the `remnant' phase. We calculate the bi-variate distribution of SFR and BHAR and define the regions in the SFR-BHAR plane that the three phases occupy. No strong correlation between BHAR and galaxy-wide SFR is found. A possible exception are galaxies with the highest SFR and the highest BHAR. We also bin the data in the same way used in several observational studies, by either measuring the mean SFR for AGN in different luminosity bins, or the mean BHAR for galaxies in bins of SFR. We find that the apparent contradiction or SFR versus BHAR for observed samples of AGN and star forming galaxies is actually caused by binning effects. The two types of samples use different projections of the full bi-variate distribution, and the full information would lead to unambiguous interpretation. We also find that a galaxy can be classified as AGN-dominated up to 1.5 Gyr after the merger-driven starburst took place. Our study is consistent with the suggestion that most low-luminosity AGN hosts do not show morphological disturbances.

Galaxies and the dark matter halos that host them are not spherically symmetric, yet spherical symmetry is a helpful simplifying approximation for idealised calculations and analysis of observational data. The assumption leads to an exact conservation of angular momentum for every particle, making the dynamics unrealistic. But how much does that inaccuracy matter in practice for analyses of stellar distribution functions, collisionless relaxation, or dark matter core-creation? We provide a general answer to this question for a wide class of aspherical systems; specifically, we consider distribution functions that are "maximally stable", i.e. that do not evolve at first order when external potentials (which arise from baryons, large scale tidal fields or infalling substructure) are applied. We show that a spherically-symmetric analysis of such systems gives rise to the false conclusion that the density of particles in phase space is ergodic (a function of energy alone). Using this idea we are able to demonstrate that: (a) observational analyses that falsely assume spherical symmetry are made more accurate by imposing a strong prior preference for near-isotropic velocity dispersions in the centre of spheroids; (b) numerical simulations that use an idealised spherically-symmetric setup can yield misleading results and should be avoided where possible; and (c) triaxial dark matter halos (formed in collisionless cosmological simulations) nearly attain our maximally-stable limit, but their evolution freezes out before reaching it.

We present a new suite of hydrodynamical simulations and use it to study, in detail, black hole and galaxy properties. The high time, spatial and mass resolution, and realistic orbits and mass ratios, down to 1:6 and 1:10, enable us to meaningfully compare star formation rate (SFR) and BH accretion rate (BHAR) timescales, temporal behaviour and relative magnitude. We find that (i) BHAR and galaxy-wide SFR are typically temporally uncorrelated, and have different variability timescales, except during the merger proper, lasting ~0.2-0.3 Gyr. BHAR and nuclear (<100 pc) SFR are better correlated, and their variability are similar. Averaging over time, the merger phase leads typically to an increase by a factor of a few in the BHAR/SFR ratio. (ii) BHAR and nuclear SFR are intrinsically proportional, but the correlation lessens if the long-term SFR is measured. (iii) Galaxies in the remnant phase are the ones most likely to be selected as systems dominated by an active galactic nucleus (AGN), because of the long time spent in this phase. (iv) The timescale over which a given diagnostic probes the SFR has a profound impact on the recovered correlations with BHAR, and on the interpretation of observational data.

We establish a controlled comparison between the properties of galactic stellar halos obtained with hydrodynamical simulations and with `particle tagging'. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints' stars directly onto a suitably defined subset of dark matter particles in a collisionless, dark-matter-only simulation.Our study shows that there are conditions under which particle tagging generates good fits to the hydrodynamical stellar density profiles of a central Milky-Way-like galaxy and its most prominent substructure. Phase-space diffusion processes are crucial to reshaping the distribution of stars in infalling spheroidal systems and hence the final stellar halo. We conclude that the success of any particular tagging scheme hinges on this diffusion being taken into account, at a minimum by making use of `live' tagging schemes, in which particles are regularly tagged throughout the evolution of a galaxy.

We study supermassive black holes (BHs) in merging galaxies, using a suite of hydrodynamical simulations with very high spatial (~10 pc) and temporal (~1 Myr) resolution, where we vary the initial mass ratio, the orbital configuration, and the gas fraction. (i) We address the question of when and why, during a merger, increased BH accretion occurs, quantifying gas inflows and BH accretion rates. (ii) We also quantify the relative effectiveness in inducing AGN activity of merger-related versus secular-related causes, by studying different stages of the encounter: the stochastic (or early) stage, the (proper) merger stage, and the remnant (or late) stage. (iii) We assess which galaxy mergers preferentially enhance BH accretion, finding that the initial mass ratio is the most important factor. (iv) We study the evolution of the BH masses, finding that the BH mass contrast tends to decrease in minor mergers and to increase in major mergers. This effect hints at the existence of a preferential range of mass ratios for BHs in the final pairing stages. (v) In both merging and dynamically quiescent galaxies, the gas accreted by the BH is not necessarily the gas with $low$ angular momentum, but the gas that $loses$ angular momentum.

We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding 10 kpc and galactocentric distances exceeding 20 kpc. Our analysis utilizes Jeans equations to generate two-dimensional acceleration maps throughout the volume; this approach is thoroughly tested on a cosmologically derived N-body+SPH simulation of a MW-like galaxy. We show that the known accelerations (gradients of the gravitational potential) can be successfully recovered in such a realistic system. Leveraging the baryonic gravitational potential derived by Bovy & Rix (2013), we show that the gravitational potential implied by the SDSS observations cannot be explained, assuming Newtonian gravity, by visible matter alone: the gravitational force experienced by stars at galactocentric distances of 20 kpc is as much as three times stronger than what can be attributed to purely visible matter. We also show that the SDSS data provide a strong constraint on the shape of the dark matter halo potential. Within galactocentric distances of 20 kpc, the dark matter halo potential is well described as an oblate halo with axis ratio qDM=0.7+/-0.1; this corresponds to an axis ratio qDM=0.4+/-0.1 for the dark matter density distribution. Because of our precise two-dimensional measurements of the acceleration of the halo stars, we can reject several MOND models as an explanation of the observed behavior.

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $\Lambda$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

The relation of central black hole mass and stellar spheroid velocity dispersion (the M-$\sigma$ relation) is one of the best-known and tightest correlations linking black holes and their host galaxies. There has been much scrutiny concerning the difficulty of obtaining accurate black hole measurements, and rightly so; however, it has been taken for granted that measurements of velocity dispersion are essentially straightforward. We examine five disk galaxies from cosmological SPH simulations and find that line-of-sight effects due to galaxy orientation can affect the measured $\sigma$ by 30%, and consequently black hole mass predictions by up to 1.0 dex. Face-on orientations correspond to systematically lower velocity dispersion measurements, while more edge-on orientations give higher velocity dispersions, due to contamination by disk stars when measuring line of sight quantities. We caution observers that the uncertainty of velocity dispersion measurements is at least 20 km/s, and can be much larger for moderate inclinations. This effect may account for some of the scatter in the locally measured M-$\sigma$ relation, particularly at the low-mass end. We provide a method for correcting observed $\sigma_{\rm los}$ values for inclination effects based on observable quantities.

ChaNGa is an N-body cosmology simulation application implemented using Charm++. In this paper, we present the parallel design of ChaNGa and address many challenges arising due to the high dynamic ranges of clustered datasets. We focus on optimizations based on adaptive techniques for scaling to more than 128K cores. We demonstrate strong scaling on up to 512K cores of Blue Waters evolving 12 and 24 billion particles. We also show strong scaling of highly clustered datasets on up to 128K cores.

We present results from a fully cosmological, very high-resolution, LCDM "zoom-in" simulation of a group of seven field dwarf galaxies with present-day virial masses in the range M_vir=4.4e8-3.6e10 Msun. The simulation includes a blastwave scheme for supernova feedback, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, a scheme for the turbulent diffusion of metals and thermal energy, and a uniform UV background that modifies the ionization and excitation state of the gas. The properties of the simulated dwarfs are strongly modulated by the depth of the gravitational potential well. All three halos with M_vir < 1e9 Msun are devoid of stars, as they never reach the density threshold for star formation of 100 atoms/cc. The other four, M_vir > 1e9 Msun dwarfs have blue colors, low star formation efficiencies, high cold gas to stellar mass ratios, and low stellar metallicities. Their bursty star formation histories are characterized by peak specific star formation rates in excess of 50-100 1/Gyr, far outside the realm of normal, more massive galaxies, and in agreement with observations of extreme emission-line starbursting dwarfs by the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey. Metal-enriched galactic outflows produce sub-solar effective yields and pollute with heavy elements a Mpc-size region of the intergalactic medium, but are not sufficient to completely quench star formation activity and are not ubiquitous in our dwarfs. Within the limited size of the sample, our simulations appear to simultaneously reproduce the observed stellar mass and cold gas content, resolved star formation histories, stellar kinematics, and metallicities of field dwarfs in the Local Volume.

Recent observations have shown that star-forming galaxies like our own Milky Way evolve kinematically into ordered thin disks over the last ~8 billion years since z=1.2, undergoing a process of "disk settling." For the first time, we study the kinematic evolution of a suite of four state of the art "zoom in" hydrodynamic simulations of galaxy formation and evolution in a fully cosmological context and compare with these observations. Until now, robust measurements of the internal kinematics of simulated galaxies were lacking as the simulations suffered from low resolution, overproduction of stars, and overly massive bulges. The current generation of simulations has made great progress in overcoming these difficulties and is ready for a kinematic analysis. We show that simulated galaxies follow the same kinematic trends as real galaxies: they progressively decrease in disordered motions (sigma_g) and increase in ordered rotation (Vrot) with time. The slopes of the relations between both sigma_g and Vrot with redshift are consistent between the simulations and the observations. In addition, the morphologies of the simulated galaxies become less disturbed with time, also consistent with observations, and they both have similarly large scatter. This match between the simulated and observed trends is a significant success for the current generation of simulations, and a first step in determining the physical processes behind disk settling.

We present more results from a fully cosmological LCDM simulation of a group of isolated dwarf galaxies that has been shown to reproduce the observed stellar mass and cold gas content, resolved star formation histories, and metallicities of dwarfs in the Local Volume. Here we investigate the energetics and timetable of the cusp-core transformation. As suggested by previous work, supernova-driven gas outflows remove dark matter (DM) cusps and create kpc-size cores in all systems having a stellar mass above 1e6 Msun. The "DM core mass removal efficiency" -- dark mass ejected per unit stellar mass -- ranges today from a few to a dozen, and increases with decreasing host mass. Because dwarfs form the bulk of their stars prior to redshift 1 and the amount of work required for DM heating and core formation scales approximately as Mvir^{5/3}, the unbinding of the DM cusp starts early and the formation of cored profiles is not as energetically onerous as previously claimed. DM particles in the cusp typically migrate to 2-3 core radii after absorbing a few percent of the energy released by supernovae. The present-day slopes of the inner dark matter mass profiles, Gamma=d log M/dlog R=2.5-3, of the simulated "Bashful" and "Doc" dwarfs are similar to those measured in the luminous Fornax and Sculptor dwarf spheroidals. None of the simulated galaxies has a circular velocity profile exceeding 20 km/s in the inner 1 kpc, implying that supernova feedback is key to solve the "too-big-to-fail" problem for Milky Way subhalos.

We present a detailed comparison between the photometric properties of the bulges of two simulated galaxies and those of a uniform sample of observed galaxies. This analysis shows that the simulated galaxies have bulges with realistic surface brightnesses for their sizes and magnitude. These two field disc galaxies have rotational velocities ~ 100 km/s and were integrated to a redshift of zero in a fully cosmological Lambda cold dark matter context as part of high-resolution smoothed particle hydrodynamic simulations. We performed bulge-disc decompositions of the galaxies using artificial observations, in order to conduct a fair comparison to observations. We also dynamically decomposed the galaxies and compared the star formation histories of the bulges to those of the entire galaxies. These star formation histories showed that the bulges were primarily formed before z = 1 and during periods of rapid star formation. Both galaxies have large amounts of early star formation, which is likely related to the relatively high bulge-to-disc ratios also measured for them. Unlike almost all previous cosmological simulations, the realistically concentrated bulges of these galaxies do not lead to unphysically high rotational velocities, causing them to naturally lie along the observed Tully-Fisher relation.

One of the principal discoveries in modern cosmology is that standard model particles (including baryons, leptons and photons) together comprise only 5% of the mass-energy budget of the Universe. The remaining 95% consists of dark energy and dark matter (DM). Consequently our picture of the universe is known as {\Lambda}CDM, with {\Lambda} denoting dark energy and CDM cold dark matter. {\Lambda}CDM is being challenged by its apparent inability to explain the low density of DM measured at the centre of cosmological systems, ranging from faint dwarf galaxies to massive clusters containing tens of galaxies the size of the Milky Way. But before making conclusions one should carefully include the effect of gas and stars, which were historically seen as merely a passive component during the assembly of galaxies. We now understand that these can in fact significantly alter the DM component, through a coupling based on rapid gravitational potential fluctuations.

We study the dynamical evolution of supermassive black holes (BHs) in merging galaxies on scales of hundreds of kpc to 10 pc, to identify the physical processes that aid or hinder the orbital decay of BHs. We present hydrodynamical simulations of galaxy mergers with a resolution of $\leq$20 pc, chosen to accurately track the motion of the nuclei and provide a realistic environment for the evolution of the BHs. We find that, during the late stages of the merger, tidal shocks inject energy in the nuclei, causing one or both nuclei to be disrupted and leaving their BH `naked', without any bound gas or stars. In many cases, the nucleus that is ultimately disrupted is that of the larger galaxy (`nuclear coup'), as star formation grows a denser nuclear cusp in the smaller galaxy. We supplement our simulations with an analytical estimate of the orbital-decay time required for the BHs to form a binary at unresolved scales, due to dynamical friction. We find that, when a nuclear coup occurs, the time-scale is much shorter than when the secondary's nucleus is disrupted, as the infalling BH is more massive, and it also finds itself in a denser stellar environment.

We compare the central mass distribution of galaxies simulated with three different models of the interstellar medium (ISM) with increasing complexity: primordial (H+He) cooling down to 10^4K, additional cooling via metal lines and to lower temperatures, and molecular hydrogen (H_2) with shielding of atomic and molecular hydrogen, in addition to metal line cooling. In order to analyze the effect of these models, we follow the evolution of four field galaxies with V_peak < 120 km/s to a redshift of zero using high-resolution Smoothed Particle Hydrodynamic simulations in a fully cosmological LCDM context. The spiral galaxies produced in simulations with either primordial cooling or H_2 physics have realistic, rising rotation curves. In contrast, the simulations with metal line cooling and otherwise similar feedback and star formation produced galaxies with the peaked rotation curves typical of most previous LCDM simulations of spiral galaxies. The less-massive bulges and non-peaked rotation curves in the galaxies simulated with primordial cooling or H_2 are linked to changes in the angular momentum distribution of the baryons. These galaxies had smaller amounts of low-angular momentum baryons because of increased gas loss from stellar feedback. When there is only primordial cooling, the star forming gas is hotter and the feedback-heated gas cools more slowly than when metal line cooling is included and so requires less energy to be expelled. When H_2 is included, the accompanying shielding produces large amounts of clumpy, cold gas where H_2 forms. Star formation in clumpy gas results in more concentrated supernova feedback and greater efficiency of mass loss. The higher feedback efficiency causes a decrease of low-angular momentum material. (abridged)

Using a set of zoomed-in cosmological simulations of high-redshift progenitors of massive galaxies, we isolate and trace the history of gas that is accreted by central supermassive black holes. We determine the origins of the accreted gas, in terms of whether it entered the galaxy during a merger event or was smoothly accreted. Furthermore, we designate whether the smoothly accreted gas is accreted via a cold flow or is shocked upon entry into the halo. For moderate-mass (10^6 - 10^7 Msun) black holes at z ~ 4, there is a preference to accrete cold flow gas than gas of shocked or merger origin. However, this result is a consequence of the fact that the entire galaxy has a higher fraction of gas from cold flows. In general, each black hole tends to accrete the same fractions of smooth- and merger-accreted gas as is contained in its host galaxy, suggesting that once gas enters a halo it becomes well-mixed, and its origins are erased. We find that the angular momentum of the gas upon halo entry is a more important factor; black holes preferentially accrete gas that had low angular momentum when it entered the galaxy, regardless of whether it was accreted smoothly or through mergers.

We examine the pressure of the star-forming interstellar medium (ISM) of Milky-Way sized disk galaxies using fully cosmological SPH+N-body, high resolution simulations. These simulations include explicit treatment of metal-line cooling in addition to dust and self-shielding, $\mathrm{H_{2}}$ based star formation. The 4 simulated halos have masses ranging from a few times $10^{10}$ to nearly $10^{12}$ solar masses. Using a kinematic decomposition of these galaxies into present-day bulge and disk components, we find that the typical pressure of the star-forming ISM in the present-day bulge is higher than that in the present-day disk by an order of magnitude. We also find that pressure of the star-forming ISM at high redshift is on average, higher than ISM pressures at low redshift. This explains the why the bulge forms at higher pressures: the disk assembles at lower redshift, when the ISM is lower pressure and the bulge forms at high redshift, when the ISM is at higher pressure. If ISM pressure and IMF variation are tied together as suggested in studies like \cite{Conroy2012}, these results could indicate a time-dependent IMF in Milky-Way like systems, as well as a different IMF in the bulge and the disk.

The cold dark matter (CDM) cosmological model has been remarkably successful in explaining cosmic structure over an enormous span of redshift, but it has faced persistent challenges from observations that probe the innermost regions of dark matter halos and the properties of the Milky Way's dwarf galaxy satellites. We review the current observational and theoretical status of these "small scale controversies." Cosmological simulations that incorporate only gravity and collisionless CDM predict halos with abundant substructure and central densities that are too high to match constraints from galaxy dynamics. The solution could lie in baryonic physics: recent numerical simulations and analytic models suggest that gravitational potential fluctuations tied to efficient supernova feedback can flatten the central cusps of halos in massive galaxies, and a combination of feedback and low star-formation efficiency could explain why most of the dark matter subhalos orbiting the Milky Way do not host visible galaxies. However, it is not clear that this solution can work in the lowest mass galaxies where discrepancies are observed. Alternatively, the small-scale conflicts could be evidence of more complex physics in the dark sector itself. For example, elastic scattering from strong dark matter self-interactions can alter predicted halo mass profiles, leading to good agreement with observations across a wide range of galaxy mass. Gravitational lensing and dynamical perturbations of tidal streams in the stellar halo provide evidence for an abundant population of low mass subhalos in accord with CDM predictions. These observational approaches will get more powerful over the next few years.

Using high resolution cosmological hydrodynamical simulations of Milky Way-massed disk galaxies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright (-15 < M_V < -8) satellite galaxies. These simulations resolve high density regions, comparable to giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H_2 formation and destruction that ties star formation to the presence of shielded, molecular gas. Before infall, supernova feedback from the clumpy, bursty star formation captured by this physically motivated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the massive satellite progenitors (Mvir > 10^9 Msun, Mstar > 10^7 Msun) compared to DM-only simulations. The progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to both heating at reionization and gas loss from initial star forming events, preserving the steep inner density profile predicted by DM-only simulations. After infall, tidal stripping acts to further reduce the central densities of the luminous satellites, particularly those that enter with cored dark matter halos, increasing the discrepancy in the central masses predicted by baryon+DM and DM-only simulations. We show that DM-only simulations, which neglect the baryonic effects described in this work, produce denser satellites with larger central velocities. We provide a simple correction to the central DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations should be used with great caution when interpreting kinematic observations of the Milky Way's dwarf satellites.

We use maximum entropy arguments to derive the phase space distribution of a virialized dark matter halo. Our distribution function gives an improved representation of the end product of violent relaxation. This is achieved by incorporating physically motivated dynamical constraints (specifically on orbital actions) which prevent arbitrary redistribution of energy. We compare the predictions with three high-resolution dark matter simulations of widely varying mass. The numerical distribution function is accurately predicted by our argument, producing an excellent match for the vast majority of particles. The remaining particles constitute the central cusp of the halo (<4% of the dark matter). They can be accounted for within the presented framework once the short dynamical timescales of the centre are taken into account.

We search for evidence of dark matter in the Milky Way by utilizing the stellar number density distribution and kinematics measured by the Sloan Digital Sky Survey (SDSS) to heliocentric distances exceeding ~10 kpc. We employ the cylindrically symmetric form of Jeans equations and focus on the morphology of the resulting acceleration maps, rather than the normalization of the total mass as done in previous, mostly local, studies. Jeans equations are first applied to a mock catalog based on a cosmologically derived N-body + SPH simulation, and the known acceleration (gradient of gravitational potential) is successfully recovered. The same simulation is also used to quantify the impact of dark matter on the total acceleration. We use Galfast, a code designed to quantitatively reproduce SDSS measurements and selection effects, to generate a synthetic stellar catalog. We apply Jeans equations to this catalog and produce two-dimensional maps of stellar acceleration. These maps reveal that in a Newtonian framework, the implied gravitational potential cannot be explained by visible matter alone. The acceleration experienced by stars at galactocentric distances of ~20 kpc is three times larger than what can be explained by purely visible matter. The application of an analytic method for estimating the dark matter halo axis ratio to SDSS data implies an oblate halo with q_DM = 0.47 +/- 0.14 within the same distance range. These techniques can be used to map the dark matter halo to much larger distances from the Galactic center using upcoming deep optical surveys, such as LSST.

We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses ranging from a few times 10^8 to nearly 10^12 solar masses. At z=0 our simulated galaxies have a baryon content and morphology typical of field galaxies. Over a stellar mass range of 2.2 x 10^3 to 4.5 x 10^10 solar masses, we find extremely good agreement between the SHM ratio in simulations and the present-day predictions from the statistical Abundance Matching Technique presented in Moster et al. (2012). This improvement over past simulations is due to a number systematic factors, each decreasing the SHM ratios: 1) gas outflows that reduce the overall SF efficiency but allow for the formation of a cold gas component 2) estimating the stellar masses of simulated galaxies using artificial observations and photometric techniques similar to those used in observations and 3) accounting for a systematic, up to 30 percent overestimate in total halo masses in DM-only simulations, due to the neglect of baryon loss over cosmic times. Our analysis suggests that stellar mass estimates based on photometric magnitudes can underestimate the contribution of old stellar populations to the total stellar mass, leading to stellar mass errors of up to 50 percent for individual galaxies. These results highlight the importance of using proper techniques to compare simulations with observations and reduce the perceived tension between the star formation efficiency in galaxy formation models and in real galaxies.

We examine the properties and evolution of a simulated polar disc galaxy. This galaxy is comprised of two orthogonal discs, one of which contains old stars (old stellar disc), and the other, containing both younger stars and the cold gas (polar disc) of the galaxy. By exploring the shape of the inner region of the dark matter halo, we are able to confirm that the halo shape is a oblate ellipsoid flattened in the direction of the polar disc. We also note that there is a twist in the shape profile, where the innermost 3 kpc of the halo flattens in the direction perpendicular to the old disc, and then aligns with the polar disc out until the virial radius. This result is then compared to the halo shape inferred from the circular velocities of the two discs. We also use the temporal information of the simulation to track the system's evolution, and identify the processes which give rise to this unusual galaxy type. We confirm the proposal that the polar disc galaxy is the result of the last major merger, where the angular moment of the interaction is orthogonal to the angle of the infalling gas. This merger is followed by the resumption of coherent gas infall. We emphasise that the disc is rapidly restored after the major merger and that after this event the galaxy begins to tilt. A significant proportion of the infalling gas comes from filaments. This infalling gas from the filament gives the gas its angular momentum, and, in the case of the polar disc galaxy, the direction of the gas filament does not change before or after the last major merger.

Motivated by the observed connection between molecular hydrogen (H2) and star formation, we present a method for tracking the non-equilibrium abundance and cooling processes of H2 and H2-based star formation in Smoothed Particle Hydrodynamic simulations. The local abundances of H2 are calculated by integrating over the hydrogen chemical network. This calculation includes the gas-phase and dust grain formation of H2, shielding of HI and H2, and photodissociation of H2 by Lyman-Werner radiation from nearby stellar populations. Because this model does not assume equilibrium abundances, it is particularly well suited for simulations that model low-metallicity environments, such as dwarf galaxies and the early Universe. We further introduce an explicit link between star formation and local H2 abundance. This link limits star formation to "star-forming regions," represented by areas with abundant H2. With this implementation, we determine the effect of H2 on star formation in a cosmological simulation of a dwarf galaxy. This simulation is the first cosmological simulation with non-equilibrium H2 abundances to be integrated to a redshift of zero or to include efficient SN feedback. We find that our simulations are consistent with the observed Tully-Fisher, global Kennicutt-Schmidt, and resolved Kennicutt-Schmidt relations. We find that the inclusion of shielding of both the atomic and molecular hydrogen and, to a lesser extent, the additional cooling from H2 at temperatures between 200 and 5000 K increases the amount of cold gas in the galaxies. The changes to the ISM result in an increased amount of cold, dense gas in the disk of the galaxy and the formation of a clumpier interstellar media (ISM). The explicit link between star formation and H2 and the clumpier ISM results in a bluer galaxy with a greater spatial distribution of star formation at a redshift of zero. (abridged)

We examine the evolution of the inner dark matter (DM) and baryonic density profile of a new sample of simulated field galaxies using fully cosmological, Lambda CDM, high resolution SPH + N-Body simulations. These simulations include explicit H2 and metal cooling, star formation (SF) and supernovae (SNe) driven gas outflows. Starting at high redshift, rapid, repeated gas outflows following bursty SF transfer energy to the DM component and significantly flatten the originally `cuspy' central DM mass profile of galaxies with present day stellar masses in the 10^4.5 -- 10^9.8 Msolar range. At z=0, the central slope of the DM density profile of our galaxies (measured between 0.3 and 0.7 kpc from their centre) is well fitted by rhoDM propto r^alpha with alpha \simeq -0.5 + 0.35 log_10(Mstar/10^8Msolar) where Mstar is the stellar mass of the galaxy and 4 < log_10 Mstar < 9.4. These values imply DM profiles flatter than those obtained in DM--only simulations and in close agreement with those inferred in galaxies from the THINGS and LITTLE THINGS survey. Only in very small halos, where by z=0 star formation has converted less than ~ 0.03% of the original baryon abundance into stars, outflows do not flatten the original cuspy DM profile out to radii resolved by our simulations. The mass (DM and baryonic) measured within the inner 500 pc of each simulated galaxy remains nearly constant over four orders of magnitudes in stellar mass for Mstar 10^9 Msolar. This finding is consistent with estimates for faint Local Group dwarfs and field galaxies. These results address one of the outstanding problems faced by the CDM model, namely the strong discrepancy between the original predictions of cuspy DM profiles and the shallower central DM distribution observed in galaxies.

We propose and successfully test against new cosmological simulations a novel analytical description of the physical processes associated with the origin of cored dark matter density profiles. In the simulations, the potential in the central kiloparsec changes on sub-dynamical timescales over the redshift interval 4 > z > 2 as repeated, energetic feedback generates large underdense bubbles of expanding gas from centrally-concentrated bursts of star formation. The model demonstrates how fluctuations in the central potential irreversibly transfer energy into collisionless particles, thus generating a dark matter core. A supply of gas undergoing collapse and rapid expansion is therefore the essential ingredient. The framework, based on a novel impulsive approximation, breaks with the reliance on adiabatic approximations which are inappropriate in the rapidly-changing limit. It shows that both outflows and galactic fountains can give rise to cusp-flattening, even when only a few per cent of the baryons form stars. Dwarf galaxies maintain their core to the present time. The model suggests that constant density dark matter cores will be generated in systems of a wide mass range if central starbursts or AGN phases are sufficiently frequent and energetic.

We investigate the formation of the first massive black holes in high redshift galaxies, with the goal of providing insights to which galaxies do or do not host massive black holes. We adopt a novel approach to forming seed black holes in galaxy halos in cosmological SPH+N-body simulations. The formation of massive black hole seeds is dictated directly by the local gas density, temperature, and metallicity, and motivated by physical models of massive black hole formation. We explore seed black hole populations as a function of halo mass and redshift, and examine how varying the efficiency of massive black hole seed formation affects the relationship between black holes and their hosts. Seed black holes tend to form in halos with mass between 10^7 and 10^9 Msun, and the formation rate is suppressed around z = 5 due to the diffusion of metals throughout the intergalactic medium. We find that the time of massive black hole formation and the occupation fraction of black holes are a function of the host halo mass. By z = 5, halos with mass M_halo > 3 x 10^9 Msun host massive black holes regardless of the efficiency of seed formation, while the occupation fraction for smaller halos increases with black hole formation efficiency. Our simulations explain why massive black holes are found in some bulgeless and dwarf galaxies, but we also predict that their occurrence becomes rarer and rarer in low-mass systems.

We make a direct comparison of the derived dark matter (DM) distributions between hydrodynamical simulations of dwarf galaxies assuming a LCDM cosmology and the observed dwarf galaxies sample from the THINGS survey in terms of (1) the rotation curve shape and (2) the logarithmic inner density slope alpha of mass density profiles. The simulations, which include the effect of baryonic feedback processes, such as gas cooling, star formation, cosmic UV background heating and most importantly physically motivated gas outflows driven by supernovae (SNe), form bulgeless galaxies with DM cores. We show that the stellar and baryonic mass is similar to that inferred from photometric and kinematic methods for galaxies of similar circular velocity. Analyzing the simulations in exactly the same way as the observational sample allows us to address directly the so-called "cusp/core" problem in the LCDM model. We show that the rotation curves of the simulated dwarf galaxies rise less steeply than CDM rotation curves and are consistent with those of the THINGS dwarf galaxies. The mean value of the logarithmic inner density slopes alpha of the simulated galaxies' dark matter density profiles is ~ -0.4 +- 0.1, which shows good agreement with \alpha = -0.29 +- 0.07 of the THINGS dwarf galaxies. The effect of non-circular motions is not significant enough to affect the results. This confirms that the baryonic feedback processes included in the simulations are efficiently able to make the initial cusps with \alpha ~ -1.0 to -1.5 predicted by dark-matter-only simulations shallower, and induce DM halos with a central mass distribution similar to that observed in nearby dwarf galaxies.

We present the analysis of a suite of simulations run with different particle-and grid-based cosmological hydrodynamical codes and compare them with observational data of the Milky Way. This is the first study to make comparisons of properties of galaxies simulated with particle and grid-based codes. Our analysis indicates that there is broad agreement between these different modelling techniques. We study the velocity dispersion - age relation for disc stars at z=0 and find that four of the simulations are more consistent with observations by Holmberg et al. (2008) in which the stellar disc appears to undergo continual/secular heating. Two other simulations are in better agreement with the Quillen & Garnett (2001) observations that suggest a "saturation" in the heating profile for young stars in the disc. None of the simulations have thin discs as old as that of the Milky Way. We also analyse the kinematics of disc stars at the time of their birth for different epochs in the galaxies' evolution and find that in some simulations old stars are born cold within the disc and are subsequently heated, while other simulations possess old stellar populations which are born relatively hot. The models which are in better agreement with observations of the Milky Way's stellar disc undergo significantly lower minor-merger/assembly activity after the last major merger - i.e. once the disc has formed. All of the simulations are significantly "hotter" than the Milky Way disc; on top of the effects of mergers, we find a "floor" in the dispersion that is related to the underlying treatment of the heating and cooling of the interstellar medium, and the low density threshold which such codes use for star formation. This finding has important implications for all studies of disc heating that use hydrodynamical codes.

Using high resolution, fully cosmological smoothed particle hydro-dynamical simulations of dwarf galaxies in a Lambda cold dark matter Universe, we show how baryons attain a final angular momentum distribution which allows pure disc galaxies to form. Blowing out substantial amounts of gas through supernovae and stellar winds, which is well supported observationally, is a key ingredient in forming bulgeless discs. We outline why galactic outflows preferentially remove low angular momentum material, and show that this is a natural result when structure forms in a cold dark matter cosmology. The driving factors are a) the mean angular momentum of accreted material increases with time, b) lower potentials at early times, c) the existence of an extended reservoir of high angular momentum gas which is not within star forming regions, meaning that only gas from the inner region (low angular momentum gas) is expelled and d) the tendency for outflows to follow the path of least resistance which is perpendicular to the disc. We also show that outflows are enhanced during mergers, thus expelling much of the gas which has lost its angular momentum during these events, and preventing the formation of "classical", merger driven bulges in low mass systems. Stars formed prior to such mergers form a diffuse, extended stellar halo component.

We perform SPH+N-body cosmological simulations of massive disk galaxies, including a formalism for black hole seed formation and growth, and find that satellite galaxies containing supermassive black hole seeds are often stripped as they merge with the primary galaxy. These events naturally create a population of "wandering" black holes that are the remnants of stripped satellite cores; galaxies like the Milky Way may host 5 -- 15 of these objects within their halos. The satellites that harbor black hole seeds are comparable to Local Group dwarf galaxies such as the Small and Large Magellanic Clouds; these galaxies are promising candidates to host nearby intermediate mass black holes. Provided that these wandering black holes retain a gaseous accretion disk from their host dwarf galaxy, they give a physical explanation for the origin and observed properties of some recently discovered off-nuclear ultraluminous X-ray sources such as HLX-1.

We have analyzed high resolution N-body simulations of dark matter halos, focusing specifically on the evolution of angular momentum. We find that not only is individual particle angular momentum not conserved, but the angular momentum of radial shells also varies over the age of the Universe by up to factors of a few. We find that torques from external structure are the most likely cause for this distribution shift. Since the model of adiabatic contraction that is often applied to model the effects of galaxy evolution on the dark-matter density profile in a halo assumes angular momentum conservation, this variation implies that there is a fundamental limit on the possible accuracy of the adiabatic contraction model in modeling the response of DM halos to the growth of galaxies.

A sample of very high resolution cosmological disk galaxy simulations is used to investigate the evolution of galaxy disk sizes back to redshift 1 within the Lambda CDM cosmology. Artificial images in the rest frame B band are generated, allowing for a measurement of disk scale lengths using surface brightness profiles as observations would, and avoiding any assumption that light must follow mass as previous models have assumed. We demonstrate that these simulated disks are an excellent match to the observed magnitude - size relation for both local disks, and for disks at z=1 in the magnitude/mass range of overlap. We disentangle the evolution seen in the population as a whole from the evolution of individual disk galaxies. In agreement with observations, our simulated disks undergo roughly 1.5 magnitudes/arcsec^2 of surface brightness dimming since z=1. We find evidence that evolution in the magnitude - size plane varies by mass, such that galaxies with M* > 10^9 M_sun undergo more evolution in size than luminosity, while dwarf galaxies tend to evolve potentially more in luminosity. The disks grow in such a way as to stay on roughly the same stellar mass - size relation with time. Finally, due to an evolving stellar mass - SFR relation, a galaxy at a given stellar mass (or size) at z=1 will reside in a more massive halo and have a higher SFR, and thus a higher luminosity, than a counterpart of the same stellar mass at z=0.

Fully cosmological, high resolution N-Body + SPH simulations are used to investigate the chemical abundance trends of stars in simulated stellar halos as a function of their origin. These simulations employ a physically motivated supernova feedback recipe, as well as metal enrichment, metal cooling and metal diffusion. As presented in an earlier paper, the simulated galaxies in this study are surrounded by stellar halos whose inner regions contain both stars accreted from satellite galaxies and stars formed in situ in the central regions of the main galaxies and later displaced by mergers into their inner halos. The abundance patterns ([Fe/H] and [O/Fe]) of halo stars located within 10 kpc of a solar-like observer are analyzed. We find that for galaxies which have not experienced a recent major merger, in situ stars at the high [Fe/H] end of the metallicity distribution function are more [alpha/Fe]-rich than accreted stars at similar [Fe/H]. This dichotomy in the [O/Fe] of halo stars at a given [Fe/H] results from the different potential wells within which in situ and accreted halo stars form. These results qualitatively match recent observations of local Milky Way halo stars. It may thus be possible for observers to uncover the relative contribution of different physical processes to the formation of stellar halos by observing such trends in the halo populations of the Milky Way, and other local L* galaxies.

We present a detailed analysis of a disc galaxy forming in a high-resolution fully cosmological simulation to investigate the nature of the outer regions of discs and their relevance for the disc formation process. Specifically, we focus on the phenomenon of misaligned disc components and find that the outer disc warp is a consequence of the misalignment between the inner disc and the surrounding hot gaseous halo. As the infalling cold gas sinks toward the centre of the galaxy, it is strongly torqued by the hot gas halo. By the time the fresh gas reaches the central disc-forming region its angular momentum is completely aligned with the spin of the hot gas halo. If the spin of the hot gas halo, in turn, is not aligned with that of the inner disc, a misaligned outer disc forms comprised of newly accreted material. The inner and outer components are misaligned with each other because they respond differently to infalling substructure and accretion. The warped disc feeds the main gas disc due to viscous angular momentum losses, but small amounts of star formation in the warp itself form a low-metallicity thick disc. We show that observations of resolved stellar populations in warped galaxies in the local universe could provide evidence for the presence of these processes and therefore indirectly reveal ongoing gas accretion and the existence of hot gas halos.

Making robust predictions for the phase space distribution of dark matter at the solar neighbourhood is vital for dark matter direct detection experiments. To date, almost all such predictions have been based on simulations that model the dark matter alone. Here, we use three cosmological hydrodynamics simulations of bright, disc dominated galaxies to include the effects of baryonic matter self-consistently for the first time. We find that the addition of baryonic physics drastically alters the dark matter profile in the vicinity of the Solar neighbourhood. A stellar/gas disc, already in place at high redshift, causes merging satellites to be dragged preferentially towards the disc plane where they are torn apart by tides. This results in an accreted dark matter disc that contributes ~0.25 - 1.5 times the non-rotating halo density at the solar position. The dark disc, unlike dark matter streams, is an equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to the Earth significantly boosts WIMP capture in the Earth and Sun, boosts the annual modulation signal, and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP mass to be determined.

Recent focus on the importance of cold, unshocked gas accretion in galaxy formation -- not explicitly included in semi-analytic studies -- motivates the following detailed comparison between two inherently different modelling techniques: direct hydrodynamical simulation and semi-analytic modelling. By analysing the physical assumptions built into the Gasoline simulation, formulae for the emergent behaviour are derived which allow immediate and accurate translation of these assumptions to the Galform semi-analytic model. The simulated halo merger history is then extracted and evolved using these equivalent equations, predicting a strikingly similar galactic system. This exercise demonstrates that it is the initial conditions and physical assumptions which are responsible for the predicted evolution, not the choice of modelling technique. On this level playing field, a previously published Galform model is applied (including additional physics such as chemical enrichment and feedback from active galactic nuclei) which leads to starkly different predictions.

In recent work, we have shown that it is possible to link quantitatively many aspects of damped Lyman alpha (DLA) absorbers in the spectra of quasars to high resolution simulations of galaxy formation. Using runs from the same series of hydrodynamic numerical studies, we consider the expected properties of Lyman alpha absorbers seen in the spectra of high redshift (z>2) gamma ray burst afterglows (GRB-DLAs). If GRBs are associated with the death of massive stars, their afterglows provide insights into otherwise unprobed regions of protogalactic objects, but detailed physical interpretations are currently embryonic. We find that median impact parameters (measured from the potential minimum) are approximately 1 kpc for GRBs compared with 4 kpc for QSO-DLAs. However, an equally important difference is that GRB-DLAs are predominantly associated with halos of mass 10^10<M_vir/M_sol<10^12, an order of magnitude larger than the hosts of QSO-DLAs. Accordingly, there are differences in the stellar properties of hosts. Our simulations accurately predict the form of the GRB-DLA HI column density distribution, producing quantitative agreement for N_HI>10^19 cm^-2, but they somewhat underpredict the incidence of low column densities N_HI<10^19 cm^-2. Line-of-sight neutral gas metallicities predicted by our simulations (10^-2 < Z/Z_sol < 1) are consistent with the modest observational constraints. Because of large internal dispersions in gas metallicities, this agreement is not significantly compromised by imposing a cut-off on the metallicity of stars able to launch GRBs (Z_star<Z_sol/3), confounding claims that the observed metallicity of GRB-DLAs poses a challenge to current GRB models. (Abridged.)

For almost two decades the properties of "dwarf" galaxies have challenged the Cold Dark Matter (CDM) paradigm of galaxy formation. Most observed dwarf galaxies consists of a rotating stellar disc embedded in a massive DM halo with a near constant-density core. Yet, models based on the CDM scenario invariably form galaxies with dense spheroidal stellar "bulges" and steep central DM profiles, as low angular momentum baryons and DM sink to the center of galaxies through accretion and repeated mergers. Processes that decrease the central density of CDM halos have been identified, but have not yet reconciled theory with observations of present day dwarfs. This failure is potentially catastrophic for the CDM model, possibly requiring a different DM particle candidate. This Letter presents new hydrodynamical simulations in a Lambda$CDM framework where analogues of dwarf galaxies, bulgeless and with a shallow central DM profile, are formed. This is achieved by resolving the inhomogeneous interstellar medium, resulting in strong outflows from supernovae explosions which remove low angular momentum gas. This inhibits the formation of bulges and decreases the dark-matter density to less than half within the central kiloparsec. Realistic dwarf galaxies are thus shown to be a natural outcome of galaxy formation in the CDM scenario.

We investigate the formation of the stellar halos of four simulated disk galaxies using high resolution, cosmological SPH + N-Body simulations. These simulations include a self-consistent treatment of all the major physical processes involved in galaxy formation. The simulated galaxies presented here each have a total mass of ~10^12 M_sun, but span a range of merger histories. These simulations allow us to study the competing importance of in-situ star formation (stars formed in the primary galaxy) and accretion of stars from subhalos in the building of stellar halos in a LambdaCDM universe. All four simulated galaxies are surrounded by a stellar halo, whose inner regions (r < 20 kpc) contain both accreted stars, and an in-situ stellar population. The outer regions of the galaxies' halos were assembled through pure accretion and disruption of satellites. Most of the in-situ halo stars formed at high redshift out of smoothly accreted cold gas in the inner 1 kpc of the galaxies' potential wells, possibly as part of their primordial disks. These stars were displaced from their central locations into the halos through a succession of major mergers. We find that the two galaxies with recently quiescent merger histories have a higher fraction of in-situ stars (~20-50%) in their inner halos than the two galaxies with many recent mergers (~5-10% in-situ fraction). Observational studies concentrating on stellar populations in the inner halo of the Milky Way will be the most affected by the presence of in-situ stars with halo kinematics, as we find that their existence in the inner few tens of kpc is a generic feature of galaxy formation.

Using high resolution SPH simulations in a fully cosmological Lambda CDM context we study the formation of a bright disk dominated galaxy that originates from a "wet" major merger at z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly following the last major merger and the present day properties of the final remnant are typical of early type spiral galaxies, with an i band B/D ~0.65, a disk scale length of 7.2 kpc, g-r = 0.5 mag, an HI line width (W_{20}/2) of 238 km/sec and total magnitude i = -22.4. The key ingredients for the formation of a dominant stellar disk component after a major merger are: i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, ii) supernova feedback that is able to partially suppress star formation during mergers and iii) relative fading of the spheroidal component. The gas fraction of the progenitors' disks does not exceed 25% at z<3, emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of disks and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disk galaxies by the high likelihood of interactions and mergers for galaxy sized halos at relatively low z.

The ACS Nearby Galaxy Survey Treasury (ANGST) is a systematic survey to establish a legacy of uniform multi-color photometry of resolved stars for a volume-limited sample of nearby galaxies (D<4 Mpc). The survey volume encompasses 69 galaxies in diverse environments, including close pairs, small & large groups, filaments, and truly isolated regions. The galaxies include a nearly complete range of morphological types spanning a factor of ~10^4 in luminosity and star formation rate. The survey data consists of images taken with ACS on HST, supplemented with archival data and new WFPC2 imaging taken after the failure of ACS. Survey images include wide field tilings covering the full radial extent of each galaxy, and single deep pointings in uncrowded regions of the most massive galaxies in the volume. The new wide field imaging in ANGST reaches median 50% completenesses of m_F475W=28.0 mag, m_F606W=27.3 mag, and m_F814W=27.3 mag, several magnitudes below the tip of the red giant branch (TRGB). The deep fields reach magnitudes sufficient to fully resolve the structure in the red clump. The resulting photometric catalogs are publicly accessible and contain over 34 million photometric measurements of >14 million stars. In this paper we present the details of the sample selection, imaging, data reduction, and the resulting photometric catalogs, along with an analysis of the photometric uncertainties (systematic and random), for both the ACS and WFPC2 imaging. We also present uniformly derived relative distances measured from the apparent magnitude of the TRGB.

We examine the pairing process of supermassive black holes (SMBHs) down to scales of 20-100 pc using a set of N-body/SPH simulations of binary mergers of disk galaxies with mass ratios of 1:4 and 1:10. Our numerical experiments are designed to represent merger events occurring at various cosmic epochs. The initial conditions of the encounters are consistent with the LambdaCDM paradigm of structure formation, and the simulations include the effects of radiative cooling, star formation, and supernovae feedback. We find that the pairing of SMBHs depends sensitively on the amount of baryonic mass preserved in the center of the companion galaxies during the last phases of the merger. In particular, due to the combination of gasdynamics and star formation, we find that a pair of SMBHs can form in 1:10 minor mergers provided that galaxies are relatively gas-rich (gas fractions of 30% of the disk mass) and that the mergers occur at relatively high redshift (z~3), when dynamical friction timescales are shorter. Since 1:10 mergers are most common events during the assembly of galaxies, and mergers are more frequent at high redshift when galaxies are also more gas-rich, our results have positive implications for future gravitational wave experiments such as the Laser Interferometer Space Antenna.

Predicting the local flux of dark matter particles is vital for dark matter direct detection experiments. To date, such predictions have been based on simulations that model the dark matter alone. Here we include the influence of the baryonic matter for the first time. We use two different approaches. Firstly, we use dark matter only simulations to estimate the expected merger history for a Milky Way mass galaxy, and then add a thin stellar disc to measure its effect. Secondly, we use three cosmological hydrodynamic simulations of Milky Way mass galaxies. In both cases, we find that a stellar/gas disc at high redshift (z~1) causes merging satellites to be preferentially dragged towards the disc plane. This results in an accreted dark matter disc that contributes ~0.25 - 1 times the non-rotating halo density at the solar position. An associated thick stellar disc forms with the dark disc and shares a similar velocity distribution. If these accreted stars can be separated from those that formed in situ, future astronomical surveys will be able to infer the properties of the dark disc from these stars. The dark disc, unlike dark matter streams, is an equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to the Earth significantly boosts WIMP capture in the Earth and Sun, increases the likelihood of direct detection at low recoil energy, boosts the annual modulation signal, and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP mass to be determined (see contribution from T. Bruch this volume).

We use high resolution cosmological hydrodynamical simulations to demonstrate that cold flow gas accretion, particularly along filaments, modifies the standard picture of gas accretion and cooling onto galaxy disks. In the standard picture, all gas is initially heated to the virial temperature of the galaxy as it enters the virial radius. Low mass galaxies are instead dominated by accretion of gas that stays well below the virial temperature, and even when a hot halo is able to develop in more massive galaxies there exist dense filaments that penetrate inside of the virial radius and deliver cold gas to the central galaxy. For galaxies up to ~L*, this cold accretion gas is responsible for the star formation in the disk at all times to the present. Even for galaxies at higher masses, cold flows dominate the growth of the disk at early times. Within this modified picture, galaxies are able to accrete a large mass of cold gas, with lower initial gas temperatures leading to shorter cooling times to reach the disk. Although star formation in the disk is mitigated by supernovae feedback, the short cooling times allow for the growth of stellar disks at higher redshifts than predicted by the standard model.

We investigate the population of z=3 damped Lyman alpha absorbers (DLAs) in a recent series of high resolution galaxy formation simulations. The simulations are of interest because they form at z=0 some of the most realistic disk galaxies to date. No free parameters are available in our study: these have been fixed by physical and z=0 observational constraints, and thus our study provides a genuine consistency test. The precise role of DLAs in galaxy formation remains in debate, but they provide a number of strong constraints on the nature of our simulated bound systems at z=3 because of their coupled information on neutral H I densities, kinematics, metallicity and estimates of star formation activity. Our results, without any parameter-tuning, closely match the observed incidence rate and column density distributions of DLAs. Our simulations are the first to reproduce the distribution of metallicities (with a median of Z_{DLA} = Z_{solar}/20) without invoking observationally unsupported mechanisms such as dust biasing. This is especially encouraging given that these simulations have previously been shown to have a realistic 0<z<2 stellar mass-metallicity relation. Additionally, we see a strong positive correlation between sightline metallicity and low-ion velocity width, the normalization and slope of which comes close to matching recent observational results. However, we somewhat underestimate the number of observed high velocity width systems; the severity of this disagreement is comparable to other recent DLA focused studies. By z=0 the majority of the z=3 neutral gas forming the DLAs has been converted into stars, in agreement with rough physical expectations. [Abridged]

Polar Ring Galaxies, such as NGC4650A, are a class of galaxy which have two kinematically distinct components that are inclined by almost 90 degrees to each other. These striking galaxies challenge our understanding of how galaxies form; the origin of their distinct components has remained uncertain, and the subject of much debate. We use high-resolution cosmological simulations of galaxy formation to show that Polar Ring Galaxies are simply an extreme example of the angular moment misalignment that occurs during the hierarchical structure formation characteristic of Cold Dark Matter cosmology. In our model, Polar Ring Galaxies form through the continuous accretion of gas whose angular momentum is misaligned with the central galaxy.

The formation of disk galaxies is one of the most outstanding problems in modern astrophysics and cosmology. We review the progress made by numerical simulations carried out on large parallel supercomputers. Recent progress stems from a combination of increased resolution and improved treatment of the astrophysical processes modeled in the simulations, such as the phenomenological description of the interstellar medium and of the process of star formation. High mass and spatial resolution is a necessary condition in order to obtain large disks comparable with observed spiral galaxies avoiding spurious dissipation of angular momentum. A realistic model of the star formation history. gas-to-stars ratio and the morphology of the stellar and gaseous component is instead controlled by the phenomenological description of the non-gravitational energy budget in the galaxy. We show that simulations of gas collapse within cold dark matter halos including a phenomenological description of supernovae blast-waves allow to obtain stellar disks with nearly exponential surface density profiles as those observed in real disk galaxies, counteracting the tendency of gas collapsing in such halos to form cuspy baryonic profiles. However, the ab-initio formation of a realistic rotationally supported disk galaxy with a pure exponential disk in a fully cosmological simulation is still an open problem. We argue that the suppression of bulge formation is related to the physics of galaxy formation during the merger of the most massive protogalactic lumps at high redshift, where the reionization of the Universe likely plays a key role. A sufficiently high resolution during this early phase of galaxy formation is also crucial to avoid artificial angular momentum loss (Abridged).

We present a new set of multi-million particle SPH simulations of the formation of disk dominated galaxies in a cosmological context. Some of these galaxies are higher resolution versions of the models already described in Governato et al (2007). To correctly compare simulations with observations we create artificial images of our simulations and from them measure photometric Bulge to Disk (B/D) ratios and disk scale lengths. We show how feedback and high force and mass resolution are necessary ingredients to form galaxies that have flatter rotation curves, larger I band disk scale lengths and smaller B/D ratios. A new simulated disk galaxy has an I-band disk scale length of 9.2 kpc and a B/D flux ratio of 0.64 (face on, dust reddened).

We construct mock galaxy catalogues to analyse clustering properties of a Lambda cold dark matter (LCDM) universe within a cosmological dark matter simulation of sufficient resolution to resolve structure down to the scale of dwarfs. We show that there is a strong age-clustering correlation for objects likely to host luminous galaxies, which includes the satellite halo (subhalo) population. Older mock galaxies are significantly more clustered in our catalog, which consists of satellite haloes as well as the central peaks of discrete haloes, selected solely by peak circular velocity. This age dependence is caused mainly by the age-clustering relation for discrete haloes, recently found by Gao et al., acting mostly on field members, combined with the tendency for older mock galaxies to lie within groups and clusters, where galaxy clustering is enhanced. Our results suggest that the clustering age dependence is manifested in real galaxies. At small scales (less than ~5 Mpc/h), the very simple assumption that galaxy colour depends solely on halo age is inconsistent with the strength of the observed clustering colour trends, where red galaxies become increasingly more clustered than blue galaxies toward smaller scales, suggesting that luminosity weighted galaxy ages do not closely trace the assembly epoch of their dark matter hosts. The age dependence is present but is much weaker for satellite haloes lying within groups and clusters than for the global population.

We examine the origin and evolution of the mass-metallicity relationship (MZR, M-Z) for galaxies using high resolution cosmological SPH + N-Body simulations that include a physically motivated description of supernovae feedback and subsequent metal enrichment. We discriminate between two sources that may contribute to the origin of the MZR: 1) metal and baryon loss due to gas outflow, or 2) inefficient star formation at the lowest galaxy masses. Our simulated galaxies reproduce the observed MZR in shape and normalization both at z=0 and z=2. We find that baryon loss occurs due to UV heating before star formation turns on in galaxies with M_baryon < 10^8 M_sun, but that some gas loss due to supernovae induced winds is required to subsequently reproduce the low effective chemical yield observed in low mass galaxies. Despite this, we show that low star formation efficiencies, regulated by supernovae feedback, are primarily responsible for the lower metallicities of low mass galaxies and the overall M-Z trend. We find that the shape of the MZR is relatively constant with redshift, but that its normalization increases with time. Simulations with no energy feedback from supernovae overproduce metals at low galaxy masses by rapidly transforming a large fraction of their gas into stars. Despite the fact that our low mass galaxies have lost a majority of their baryons, they are still the most gas rich objects in our simulations due to their low star formation efficiencies.

We used fully cosmological, high resolution N-body + SPH simulations to follow the formation of disk galaxies with rotational velocities between 135 and 270 km/sec in a Lambda CDM universe. The simulations include gas cooling, star formation, the effects of a uniform UV background and a physically motivated description of feedback from supernovae. The host dark matter halos have a spin and last major merger redshift typical of galaxy sized halos as measured in recent large scale N--Body simulations. The simulated galaxies form rotationally supported disks with realistic exponential scale lengths and fall on both the I-band and baryonic Tully Fisher relations. An extended stellar disk forms inside the Milky Way sized halo immediately after the last major merger. The combination of UV background and SN feedback drastically reduces the number of visible satellites orbiting inside a Milky Way sized halo, bringing it in fair agreement with observations. Our simulations predict that the average age of a primary galaxy's stellar population decreases with mass, because feedback delays star formation in less massive galaxies. Galaxies have stellar masses and current star formation rates as a function of total mass that are in good agreement with observational data. We discuss how both high mass and force resolution and a realistic description of star formation and feedback are important ingredients to match the observed properties of galaxies.

We present an analysis of star formation and feedback recipes appropriate for galactic smoothed particle hydrodynamics simulations. Using an isolated Milky Way-like galaxy, we constrain these recipes based on well-established observational results. Our star formation recipe is based on that of Katz (1992) with the additional inclusion of physically motivated supernova feedback recipes. We propose a new feedback recipe in which type II supernovae are modelled using an analytical treatment of blastwaves. With this feedback mechanism and a tuning of other star formation parameters, the star formation in our isolated Milky Way-like galaxy is constant and follows the slope and normalisation of the observed Schmidt law. In addition, we reproduce the low density cutoff and filamentary structure of star formation observed in disk galaxies. Our final recipe will enable better comparison of cosmological N-body simulations with observations.

Two well studied dwarf galaxies -- NGC 3109 and NGC 6822 -- present some of the strongest observational support for a flat core at the center of galactic dark matter (DM) halos. We use detailed cosmologically motivated numerical models to investigate the systematics and the accuracy of recovering parameters of the galaxies. Some of our models match the observed structure of the two galaxies remarkably well. Our analysis shows that the rotation curves of these two galaxies are instead quite compatible with their DM halos having steep cuspy density profiles. The rotation curves in our models are measured using standard observational techniques. The models reproduce the rotation curves of both galaxies, the disk surface brightness profiles as well as the profile of isophotal ellipticity and position angle. The models are centrally dominated by baryons; however, the dark matter component is globally dominant. The simulated disk mass is marginally consistent with a stellar mass-to-light ratio in agreement with the observed colors. We show that non-circular motions combined with gas pressure support and projection effects results in a large underestimation of the circular velocity in the central $\sim 1$ kpc region, creating the illusion of a constant density core. Although the systematic effects mentioned above are stronger in barred systems, they are also present in axisymetric disks. Our results strongly suggest that there is no contradiction between the observed rotation curves in dwarf galaxies and the cuspy central dark matter density profiles predicted by Cold Dark Matter models.

We used fully cosmological, high resolution N-body + SPH simulations to follow the formation of disk galaxies with a rotational velocity between 140 and 280 Km/sec in a Lambda CDM universe. The simulations include gas cooling, star formation, the effects of a uniform UV background and a physically motivated description of feedback from supernovae. The host dark matter halos have a spin and last major merger redshift typical of galaxy sized halos as measured in recent large scale N-Body simulations. Galaxies formed rotationally supported disks with realistic exponential scale lengths and fall on the I-band and baryonic Tully Fisher relations. The combination of UV background and SN feedback drastically reduced the number of visible satellites orbiting inside a Milky Way sized halo, bringing it fair agreement with observations. Feedback delays SF in small galaxies and more massive ones contain older stellar populations. The current star formation rates as a function of galaxy stellar mass are in good agreement with those measured by the SDSS.

We used fully cosmological, high resolution N-body+SPH simulations to follow the formation of disk galaxies with a rotational velocity between 140 and 280 Km/sec in a Lambda CDM universe. The simulations include gas cooling, star formation (SF), the effects of a uniform UV background and a physically motivated description of feedback from supernovae (SN). Feedback parameters have been chosen to match the star formation rate and interstellar medium (ISM) properties of local galaxies. In cosmological simulations galaxies formed rotationally supported disks with realistic exponential scale lengths and fall on the I-band and baryonic Tully Fisher relations. The combination of UV background and SN feedback drastically reduced the number of visible satellites orbiting inside a Milky Way sized halo, bringing it in fair agreement with observations. Feedback delays SF in small galaxies and more massive ones contain older stellar populations. Here we focus on the SF and feedback implementations. We also briefly discuss how high mass and force resolution and a realistic description of SF and feedback are important ingredients to match the observed properties of galaxies.

We use cosmological LCDM numerical simulations to model the evolution of the substructure population in sixteen dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10^11 to 10^14 Msun at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation: dn/dv=(1/8)(v_cmax/v_cmax,host)^-4. Halo to halo variance in the VDF is a factor of roughly two to four. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for r >~ 0.3r_vir, where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease toward the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias toward more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power law cosmologies of P ~ k^n, and demonstrate that a steeper spectral index, n, results in significantly less substructure.

We use numerical simulations in a Lambda CDM cosmology to model density profiles in a set of 16 dark matter haloes with resolutions of up to 7 million particles within the virial radius. These simulations allow us to follow robustly the formation and evolution of the central cusp over a large mass range of 10^11 to 10^14 M_sun, down to approximately 0.5% of the virial radius, and from redshift 5 to the present. The cusp of the density profile is set at redshifts of 2 or greater and remains remarkably stable to the present time, when considered in non-comoving coordinates. We fit our haloes to a 2 parameter profile where the steepness of the asymptotic cusp is given by gamma, and its radial extent is described by the concentration, c_gamma. In our simulations, we find gamma = 1.4 - 0.08Log(M/M_*) for haloes of 0.01M_* to 1000M_*, with a large scatter of gamma ~ +/-0.3$; and c_gamma = 8*M/M_*^{-0.15}, with a large M/M_* dependent scatter roughly equal to +/- c_gamma. Our redshift zero haloes have inner slope parameters ranging approximately from r^{-1} to r^{-1.5}, with a median of roughly r^{-1.3}. This 2 parameter profile fit works well for all our halo types, whether or not they show evidence of a steep asymptotic cusp. We also model a cluster in power law cosmologies of P ~ k^n (n=0,-1,-2,-2.7). We find larger concentration radii and shallower cusps for steeper n. The minimum resolved radius is well described by the mean interparticle separation. The trend of steeper and more concentrated cusps for smaller $M/M_*$ haloes clearly shows that dwarf sized Lambda CDM haloes have, on average, significantly steeper density profiles within the inner few percent of the virial radius than inferred from recent observations. Code to reproduce this profile can be downloaded from http://www.icc.dur.ac.uk/~reed/profile.html

Using N-body simulations of flat, dark energy dominated cosmologies, we show that galaxies around simulated binary systems resembling the Local Group (LG) have low peculiar velocities, in good agreement with observational data. We have compared results for LG-like systems selected from large, high resolution simulations of three cosmologies: a LCDM model, a LWDM model with a 2 keV warm dark matter candidate and a quintessence model (QCDM) with an equation of state parameter w=-0.6. The Hubble flow is significant colder around LGs selected in a flat, Lambda dominated cosmology than around LGs in open or critical models, showing that a dark energy component manifests itself on the scales of nearby galaxies, cooling galaxy peculiar motions. Flows in the LWDM and QCDM models are marginally colder than in the LCDM one. The results of our simulations have been compared to existing data and a new data set of 28 nearby galaxies with robust distance measures (Cepheids and Surface Brightness Fluctuations). The measured line-of-sight velocity dispersion is sigma = 88 +- 20 km/sec x (R/7 Mpc). The best agreement with observations is found for LGs selected in the $\Lambda$CDM cosmology in environments with -0.1 <delta_rho/rho < 0.6 on scales of 7 Mpc, in agreement with existing observational estimates on the local matter density. These results provide new, independent evidence for the presence of dark energy on scales of few Mpc, corroborating the evidence gathered from observations of distant objects and the early Universe.

We use a multi million particle N-body + SPH simulation to follow the formation of a rich galaxy cluster in a Lambda+CDM cosmology, with the goal of understanding the origin and properties of intracluster stars. The simulation includes gas cooling, star formation, the effects of a uniform UVB and feedback from supernovae. Halos that host galaxies as faint as M_R = -19.0 are resolved by this simulation, which includes 85% of the total galaxy luminosity in a rich cluster. We find that the accumulation of intracluster light (ICL) is an ongoing process, linked to infall and stripping events. The unbound star fraction increases with time and is 20% at z = 0, consistent with observations of galaxy clusters. The surface brightness profile of the cD shows an excess compared to a de Vaucouleur profile near 200 kpc, which is also consistent with observations. Both massive and small galaxies contribute substantially to the formation of the ICL, with stars stripped preferentially from the outer parts of their stellar distributions. Simulated observations of planetary nebulae (PNe) show significant substructure in velocity space. Despite this, individual intracluster PNe might be useful mass tracers if more than 5 fields at a range of radii have measured line-of-sight velocities, where an accurate mass calculation depends more on the number of fields than the number of PNe measured per field. However, the orbits of IC stars are more anisotropic than those of galaxies or dark matter, which leads to a systematic underestimate of cluster mass relative to that calculated with galaxies, if not accounted for in dynamical models. Overall, the properties of ICL formed in a hierarchical scenario are in good agreement with current observations. (Abridged)

We review evidence that the census of Milky Way satellites similar to those known may be incomplete at low latitude due to obscuration and in the outer halo due to a decreasing sensitivity to dwarf satellites with distance. We evaluate the possible impact that incompleteness has on comparisons with substructure models by estimating corrections to the known number of dwarfs using empirical and theoretical models. If we assume that the true distribution of Milky Way satellites is uniform with latitude, then we estimate a 33% incompleteness in the total number of dwarfs due to obscuration at low latitude. Similarly, if we suppose that the radial distribution of Milky Way satellites matches that of M31, or that of the oldest sub-halos or the most massive sub-halos in a simulation, we estimate a total number of Milky Way dwarfs ranging from 1 -- 3 times the known population. Although the true level of incompleteness is quite uncertain, the fact that our extrapolations yield average total numbers of MW dwarfs that are realistically 1.5 -- 2 times the known population, shows that incompleteness needs to be taken seriously when comparing to models of dwarf galaxy formation. Interestingly, the radial distribution of the oldest sub-halos in a Lambda+CDM simulation of a Milky Way-like galaxy possess a close match to the observed distribution of M31's satellites, which suggests that reionization may be an important factor controlling the observability of sub-halos. We also assess the prospects for a new SDSS search for Milky Way satellites to constrain the possible incompleteness in the outer halo.

We simulate the formation of a realistic disk galaxy within the hierarchical scenario of structure formation and study its internal properties to the present epoch. We compare results from a LambdaCDM simulation with a LambdaWDM (2keV) simulation that forms significantly less small scale structure. We show how high mass and force resolution in both the gas and dark matter components play an important role in solving the angular momentum catastrophe claimed from previous simulations of galaxy formation within the hierarchical framework. The stellar material in the disk component has a final specific angular momentum equal to 40% and 90% of that of the dark halo in the LambdaCDM and LambdaWDM models respectively. The LambdaWDM galaxy has a drastically reduced satellite population and a negligible stellar spheroidal component. Encounters with satellites play only a minor role in disturbing the disk. Satellites possess a variety of star formation histories linked to mergers and pericentric passages along their orbit around the primary galaxy. In both cosmologies, the galactic halo retains most of the baryons accreted and builds up a hot gas phase with a substantial X-ray emission. Therefore, while we have been successful in creating a realistic stellar disk in a massive galaxy within the LambdaCDM scenario, energy injection emerges as necessary ingredient to reduce the baryon fraction in galactic halos, independent of the cosmology adopted. (abridged)

We present three-dimensional self-gravitating smoothed-particle hydrodynamics (SPH) simulations of an isothermal gaseous disc interacting with an embedded planet. Discs of varying stability are simulated with planets ranging from 10 Earth-masses to 2 Jupiter-masses. The SPH technique provides the large dynamic range needed to accurately capture the large scale behavior of the disc as well as the small scale interaction of the planet with surrounding material. Most runs used 10^5 gas particles, giving us the spatial resolution required to observe the formation of planets. We find four regions in parameter space: low-mass planets undergo Type I migration; higher-mass planets can form a gap; the gravitational instability mode of planet formation in marginally stable discs can be triggered by embedded planets; discs that are completely unstable can fragment to form many planets. The disc stability is the most important factor in determing which interaction a system will exhibit. For the stable disc cases, our migration and accretion time-scales are shorter and scale differently than previously suggested.

We present three-dimensional SPH calculations of giant planets embedded in gaseous disks. Our findings are collected into a map of parameter space, exhibiting four distinct regions: Type I migration, gap formation, triggered formation of more planets, and wholly unstable disks. For Type I migration of the planet due to secular interactions with the disk material, the migration rate depends linearly on the disk mass, and is independent of the initial planet mass. For more massive disks, the planet can disturb the disk strongly enough to trigger the collapse of gas into additional giant planets. When additional planets form, their interaction as point masses dominates the subsequent behavior of the system. This mechanism allows for the rapid formation of Jupiter-mass and higher planets. Migration due to interaction with the disk can significantly change the orbits of giant planets in gas disks.

We present galaxy circular velocity functions (GCVFs) for 34 low redshift (z < 0.15) clusters identified in the Sloan Digital Sky Survey (SDSS), for fifteen clusters drawn from dark matter simulations of hierarchical structure growth in a LambdaCDM cosmology, and for ~29,000 SDSS field galaxies. We find that the observed and simulated cluster GCVFs take the form of a power law. The cumulative GCVFs of the simulated clusters are very similar across a wide range of cluster masses, provided individual subhalo circular velocities are scaled by the circular velocities of the parent cluster. Once all sources of error are accounted for, the intrinsic scatter in the cumulative, scaled observed cluster GCVF is consistent with the simulations. The slope of the observed cluster GCVF is -2.5, independent of cluster velocity dispersion. The average slope of the simulated GCVFs is somewhat steeper, although formally consistent given the errors. Using our highest resolution cluster, we find that the effects of baryons on galaxy rotation curves is to flatten the simulated cluster GCVF into better agreement with observations. Finally, we find that the field GCVF deviates significantly from a power law, being flatter than the cluster GCVF at circular velocities less than 200 km/s, and steeper at circular velocities greater than 200 km/s.

We use a high resolution $\Lambda$CDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of fifteen, in a 50 $h^{-1}$Mpc volume containing 80 million particles. Our low redshift results allow us to probe low $\sigma$ density fluctuations significantly beyond the range of previous cosmological simulations. The Sheth and Tormen mass function provides an excellent match to all of our data except for redshifts of ten and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 percent for the $10^{10}-10^{11} \msun$ haloes sampled at redshift 15. Our results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence. We provide an empirical fit to our data that corrects for the overprediction of extremely rare objects by the Sheth and Tormen mass function. This overprediction has implications for studies that use the number densities of similarly rare objects as cosmological probes. For example, the number density of high redshift (z $\simeq$ 6) QSOs, which are thought to be hosted by haloes at 5$\sigma$ peaks in the fluctuation field, are likely to be overpredicted by at least a factor of 50%. We test the sensitivity of our results to force accuracy, starting redshift, and halo finding algorithm.

We use a high resolution collisionless simulation of a Virgo--like cluster in a $\Lambda$CDM cosmology to determine the velocity and clustering properties of the diffuse stellar component in the intracluster region at the present epoch. The simulated cluster builds up hierarchically and tidal interactions between member galaxies and the cluster potential produce a diffuse stellar component free-flying in the intracluster medium. Here we adopt an empirical scheme to identify tracers of the stellar component in the simulation and hence study its properties. We find that at $z=0$ the intracluster stellar light is mostly unrelaxed in velocity space and clustered in structures whose typical clustering radii are about 50 kpc at R=400--500 kpc from the cluster center, and predict the radial velocity distribution expected in spectroscopic follow-up surveys. Finally, we compare the spatial clustering in the simulation with the properties of the Virgo intracluster stellar population, as traced by ongoing intracluster planetary nebulae surveys in Virgo. The preliminary results indicate a substantial agreement with the observed clustering properties of the diffuse stellar population in Virgo.

We study the evolution of dark matter satellites orbiting inside more massive haloes using semi-analytical tools coupled with high-resolution N-Body simulations. We select initial satellite sizes, masses, orbital energies, and eccentricities as predicted by hierarchical models of structure formation. Both the satellite and the main halo are described by a Navarro, Frenk & White density profile with various concentrations.

We present a set of gas-dynamical simulations of galaxy groups and clusters aimed at exploring the effect of non-gravitational heating. We use GASOLINE, a parallel Tree+SPH code, to simulate the formation of four cosmic halos with temperature 0.5<T<8 keV. Non-gravitational heating is implemented in two different ways: (1) by imposing a minimum entropy floor at a given redshift, 1<z<5; (2) by gradually heating gas, proportionally to the SN rate expected from semi-analytical modeling of galaxy formation. Our main results are the following. (a) An extra heating energy of about 1 keV per gas particle is required to reproduce the observed Lx-T relation, independent of whether it is provided so as to create an entropy floor of 50-100 keV cm^2, or is modulated in redshift; our SN feedback recipe provides only 1/3 keV/part. (b) The M-T relation is almost unaffected by non-gravitational heating and follows the M T^{3/2} scaling, with a normalization ~40% higher than observed, independent of the heating scheme. The inclusion of cooling in a run of a small group has the effects of increasing T_ew by ~30%, possibly reconciling simulated and observed M-T relations, and of decreasing Lx by ~40%. In spite of the inclusion of SN feedback energy, almost 40% of the gas becomes cold, in excess of current observational estimates. (abridged)

Efforts to understand the deviation of the L--T relation from a simple scaling law valid for clusters and groups have triggered a number of interesting studies on the subject. Techniques and approaches differ widely but most works agree on the important role played by gas cooling and heating sources like AGNs and SNe. Observations set useful constraints on the evolution of the intracluster medium (ICM): a 100KeV/cm^2 entropy floor in the core of groups and about 5--15% of baryons being converted into stars. However, essential details like the nature of the dominant heating mechanism and the quantitative importance of cooling still need to be addressed. I suggest that a new generation of high resolution N-body simulations and a quantitative comparison of results between different approaches is required to improve results and increase our understanding of the problem.

Recent observations suggest that Carina and other nearby dwarf spheroidal galaxies are surrounded by unbound stars tidally stripped by the Milky Way. We run high-resolution N-Body simulations of dwarf galaxies orbiting within the Milky Way halo to determine if such observations can be explained with dark matter potentials as those implied by current structure formation models. We show that tidal forces acting on dwarfs with constant density cores or with cuspy profiles having a low concentration parameter ($c < 5$) lead to flat outer stellar density profiles like that of Carina for a variety of orbital configurations. On the contrary, it is more difficult to remove stars from cuspy dark matter halos with concentrations as high as predicted by CDM models at the mass scale of dwarf galaxies ($c \simgt 10$) and the data can only be reproduced assuming nearly radial orbits. Our simulations show that Carina is losing mass at a fractional rate $< 0.1$ Gyr$^{-1}$ and its mass-to-light ratio could be inflated by at most a factor of 2 due to unbound stars projected along the line of sight. We follow the evolution of the tidal debris within a triaxial clumpy cold dark matter Milky Way halo which causes differential precession and small scale heating of the stellar streams. This renders their use as a dynamical tracer of the Galactic potential practically useless, but does provide a novel test of the nature of the dark matter. Models with warm dark matter (WDM) or fluid dark matter (FDM) produce dwarf halos with lower central densities than CDM and would be consistent with the observed tidal tails even for orbits with eccentricities as low as indicated by current data on nearby dwarf spheroidals. Galactic halos in FDM are smooth and spherical and would be favored by the detection of coherent streams.

We present a new algorithm (PINOCCHIO, PINpointing Orbit-Crossing Collapsed HIerarchical objects) to predict accurately the formation and evolution of individual dark matter haloes in a given realization of an initial linear density field. Compared with the halo population formed in a large (360^3 particles) collisionless simulation of a CDM universe, our method is able to predict to better than 10 per cent statistical quantities such as the mass function, two-point correlation function and progenitor mass function of the haloes. Masses of individual haloes are estimated accurately as well, with errors typically of order 30 per cent in the mass range well resolved by the numerical simulation. These results show that the hierarchical formation of dark matter haloes can be accurately predicted using local approximations to the dynamics when the correlations in the initial density field are properly taken into account. The approach allows one to automatically generate a large ensemble of accurate merging histories of haloes with complete knowledge of their spatial distribution. The construction of the full merger tree for a 256^3 realisation requires a few hours of CPU-time on a personal computer, orders of magnitude faster than the corresponding N-body simulation would take, and without needing any extensive post-processing. The technique can be efficiently used, for instance, for generating the input for galaxy formation modeling.

We investigate how the survival of dark matter satellites inside virialized halos depends on tidal stripping and dynamical friction. We use an analytic approach and then compare the results with N-Body simulations. Both the satellites and the primary halos are similar to cosmological haloes and have NFW density profiles. Satellites can either merge with the primary halo or continue to move on barely perturbed orbits, eventually being disrupted, depending on the relative strength of friction and tidal forces. We propose that their actual fate depends simply on their mass ratio relative to the primary halo.

We present results from high-resolution Tree+SPH simulations of galaxy clusters and groups, aimed at studying the effect of non-gravitational heating on the entropy of the ICM. We simulate three systems, having emission-weighted temperature T=0.6,1 and 3 keV, with spatial resolution better than 1% of the virial radius. We consider the effect of different prescriptions for non-gravitational ICM heating, such as SN energy feedback, as predicted by semi-analytical models of galaxy formation, and two different minimum entropy floors, S_fl=50 and 100 keV cm^2, imposed at z=3. Simulations with only gravitational heating nicely reproduce predictions from self-similar ICM models, while extra heating is shown to break the self-similarity, by a degree which depends on total injected energy and on cluster mass. We use observational results on the excess entropy in central regions of galaxy systems, to constrain the amount of extra-heating required. We find that setting the entropy floor S_fl=50 keV cm^2, which corresponds to an extra heating energy of about 1 keV per particle, is able to reproduce the observed excess of ICM entropy.

We use a cosmological simulation of the Local Group to make quantitative and speculative predictions for direct detection experiments. Cold dark matter (CDM) halos form via a complex series of mergers, accretion events and violent relaxation which precludes the formation of significant caustic features predicted by axially symmetric collapse. The halo density profiles are combined with observational constraints on the galactic mass distribution to constrain the local density of cold dark matter to lie in the range 0.18 <~ rho_CDM(R_solar)/GeV cm^-3 <~ 0.30. In velocity space, coherent streams of dark matter from tidally disrupted halos fill the halo and provide a tracer of the merging hierarchy. The particle velocities within triaxial CDM halos cannot be approximated by a simple Maxwellian distribution and is radially biased at the solar position. The detailed phase space structure within the solar system will depend on the early merger history of the progenitor halos and the importance of major mergers over accretion dominated growth. We follow the formation of a ``Draco'' sized dSph halo of mass 10^8M_solar with several million particles and high force accuracy. Its internal structure and substructure resembles that of galactic or cluster mass halos: the density profile has a singular central cusp and it contains thousands of sub-halos orbiting within its virial radius demonstrating a self-similar nature to collisionless dark matter sub-clustering. The singular cores of substructure halos always survive complete tidal disruption although mass loss is continuous and rapid. Extrapolating wildly to earth mass halos with velocity dispersion of 1 m s^-1 (roughly equal to the free streaming scale for neutralinos) we find that most of the dark matter may remain attached to bound subhalos. (Abridged)

For the last few years the Cold Dark Matter model (ticker: CDM), has been the dominant theory of structure formation. We briefly review the recent advancements and predictions of the model in the field of galaxy clusters. A new set of very high resolution simulations of galaxy clusters show that they have (1) density profiles with central slopes very close to -1.6 and (2) abundance of subhalos similar to the ones observed in real clusters. These results show a remarkably small cluster to cluster variation and a weak dependence from the particular CDM cosmology chosen (LCDM having ~40% less substructure than SCDM). While still a speculative theory with a high prediction/evidence ratio, subject to strong challenges from observational data and competition from other hierarchical theories, we give CDM a rating of ``market outperform'' and of ``long term BUY''.

We present results from high-resolution N-Body/SPH simulations of rotationally supported dwarf irregular galaxies moving on bound orbits in the massive dark matter halo of the Milky Way.The dwarf models span a range in disk surface density and the masses and sizes of their dark halos are consistent with the predictions of cold dark matter cosmogonies. We show that the strong tidal field of the Milky Way determines severe mass loss in their halos and disks and induces bar and bending instabilities that transform low surface brightness dwarfs (LSBs) into dwarf spheroidals (dSphs) and high surface brightness dwarfs (HSBs) into dwarf ellipticals (dEs) in less than 10 Gyr. The final central velocity dispersions of the remnants are in the range 8-30 km/s and their final $v/\sigma$ falls to values $< 0.5$, matching well the kinematics of early-type dwarfs. The transformation requires the orbital time of the dwarf to be $\simlt 3-4$ Gyr, which implies a halo as massive and extended as predicted by hierarchical models of galaxy formation to explain the origin of even the farthest dSph satellites of the Milky Way, Leo I and Leo II. Only dwarfs with central dark matter densities as high as those of Draco and Ursa Minor can survive for 10 Gyr in the proximity of the Milky Way: this is naturally achieved within hierarchical models, where the densest objects should have small orbital times due to their early formation epochs. Part of the gas is stripped and part is funneled to the center due to the bar, generating one strong burst of star formation in HSBs and smaller, multiple bursts in LSBs. Extended low-surface brightness stellar and gaseous streams originate from LSBs and, when projected along the line of sight, can lead to overestimate the mass-to-light ratio of the bound remnant by a factor $\simlt 2$,

N-Body/SPH simulations are used to study the evolution of dwarf irregular galaxies (dIrrs) entering the dark matter halo of the Milky Way or M31 on plunging orbits. We propose a new dynamical mechanism driving the evolution of gas rich, rotationally supported dIrrs, mostly found at the outskirts of the Local Group (LG), into gas free, pressure supported dwarf spheroidals (dSphs) or dwarf ellipticals (dEs), observed to cluster around the two giant spirals. The initial model galaxies are exponential disks embedded in massive dark matter halos and reproduce nearby dIrrs. Repeated tidal shocks at the pericenter of their orbit partially strip their halo and disk and trigger dynamical instabilities that dramatically reshape their stellar component. After only 2-3 orbits low surface brightness (LSB) dIrrs are transformed into dSphs, while high surface brightness (HSB) dIrrs evolve into dEs. This evolutionary mechanism naturally leads to the morphology-density relation observed for LG dwarfs. Dwarfs surrounded by very dense dark matter halos, like the archetypical dIrr GR8, are turned into Draco or Ursa Minor, the faintest and most dark matter dominated among LG dSphs. If disks include a gaseous component, this is both tidally stripped and consumed in periodic bursts of star formation. The resulting star formation histories are in good qualitative agreement with those derived using HST color-magnitude diagrams for local dSphs.

Can N-body simulations reliably determine the structural properties of dark matter halos? Focussing on a Virgo-sized galaxy cluster, we increase the resolution of current ``high resolution simulations'' by almost an order of magnitude to examine the convergence of the important physical quantities. We have 4 million particles within the cluster and force resolution 0.5 kpc/h (0.05% of the virial radius). The central density profile has a logarithmic slope of -1.5, as found in lower resolution studies of the same halo, indicating that the profile has converged to the ``physical'' limit down to scales of a few kpc. Also the abundance of substructure is consistent with that derived from lower resolution runs; on the scales explored, the mass and circular velocity functions are close to power laws of exponents ~ -1.9 and -4. Overmerging appears to be globally unimportant for suhalos with circular velocities > 100 km/s. We can trace most of the cluster progenitors from z=3 to the present; the central object (the dark matter analog of a cD galaxy)is assembled between z=3 and 1 from the merging of a dozen halos with v_circ \sim 300 km/s. The mean circular velocity of the subhalos decreases by ~ 20% over 5 billion years, due to tidal mass loss. The velocity dispersions of halos and dark matter globally agree within 10%, but the halos are spatially anti-biased, and, in the very central region of the cluster, they show positive velocity bias; however, this effect appears to depend on numerical resolution.

We combine semi-analytical methods with a ultra-high resolution simulation of a galaxy cluster (of mass 2.3 10^14h-1Msolar, and 4 10^6 particles within its virial radius) formed in a standard CDM universe to study the spatial distribution and orbital properties of the present-day descendents of Lyman Break Galaxies (LBGs). At the present time only five (out of 12) of halos containing LBGs survive as separate entities inside the cluster virial radius. Their circular velocities are in the range 200 - 550 km/sec. Seven halos merged together to form the central object at the very center of the cluster. Using semi-analytical modeling of galaxy evolution we show that descendents of halos containing LBGs now host giant elliptical galaxies. Galaxy orbits are radial, with a pericenter to apocenter ratio of about 1:5. The orbital eccentricities of LBGs descendents are statistically indistinguishable from those of the average galaxy population inside the cluster, suggesting that the orbits of these galaxies are not significantly affected by dynamical friction decay after the formation of the cluster's main body. In this cluster, possibly due to its early formation time, the descendents of LBGs are contained within the central 60% of the cluster virial radius and have an orbital velocity dispersion lower than the global galaxy population, originating a mild luminosity segregation for the brightest cluster members. Mass estimates based only on LBGs descendents (especially including the central cD) reflect this bias in space and velocity and underestimate the total mass of this well virialized cluster by up to a factor of two compared to estimates using at least 20 cluster members.

We use numerical simulations to examine the substructure within galactic and cluster mass halos that form within a hierarchical universe. Clusters are easily reproduced with a steep mass spectrum of thousands of substructure clumps that closely matches observations. However, the survival of dark matter substructure also occurs on galactic scales, leading to the remarkable result that galaxy halos appear as scaled versions of galaxy clusters. The model predicts that the virialised extent of the Milky Way's halo should contain about 500 satellites with circular velocities larger than Draco and Ursa-Minor i.e. bound masses > 10^8Mo and tidally limited sizes > kpc. The substructure clumps are on orbits that take a large fraction of them through the stellar disk leading to significant resonant and impulsive heating. Their abundance and singular density profiles has important implications for the existence of old thin disks, cold stellar streams, gravitational lensing and indirect/direct detection experiments.

We study the inhomogeneous reionization in a critical density CDM universe due to stellar sources, including Population III objects. The spatial distribution of the sources is obtained from high resolution numerical N-body simulations. We calculate the source properties taking into account a self-consistent treatment of both radiative (ie ionizing and H2 -photodissociating photons) and stellar (ie SN explosions) feedbacks regulated by massive stars. This allows us to describe the topology of the ionized and dissociated regions at various cosmic epochs and derive the evolution of H, He, and H2 filling factors, soft UV background, cosmic star formation rate and the final fate of ionizing objects. The main results are: (i) galaxies reionize the IGM by z~10 (with some uncertainty related to the gas clumping factor), whereas H2 is completely dissociated already by z~25; (ii) reionization is mostly due to the relatively massive objects which collapse via H line cooling, while objects whose formation relies on H2 cooling alone are insufficient to this aim; (iii) the diffuse soft UV background is the major source of radiative feedback effects for z<15; at higher z direct flux from neighboring objects dominates; (iv) the match of the calculated cosmic star formation history with the one observed at lower redshifts suggests that the conversion efficiency of baryons into stars is ~1%; (v) we find that a very large population of dark objects which failed to form stars is present by z~8. We discuss and compare our results with similar previous studies.

The evolution of a small satellite inside a more massive truncated isothermal spherical halo is studied using both the Theory of Linear Response for dynamical friction and N-Body simulations. The analytical approach includes the effects of the gravitational wake, of the tidal deformation and the shift of the barycenter of the primary, so unifying the local versus global interpretation of dynamical friction. Sizes, masses, orbital energies and eccentricities are chosen as expected in hierarchical clustering models. We find that in general the drag force in self-gravitating backgrounds is weaker than in uniform media and that the orbital decay is not accompanied by a significant circularization. We also show that the dynamical friction time scale is weakly dependent on the initial circularity. We provide a fitting formula for the decay time that includes the effect of mass and angular momentum loss. Live satellites with dense cores can survive disruption up to an Hubble time within the primary, notwithstanding the initial choice of orbital parameters. Dwarf spheroidal satellites of the Milky Way, like Sagittarius A and Fornax, have already suffered mass stripping and, with their present masses, the sinking times exceed 10 Gyr even if they are on very eccentric orbits.

We study the evolution of galaxy satellites with high resolution N-body simulations. Satellites are modeled as replicas of typical low and high surface brightness galaxies (LSBs and HSBs). Encounters on high eccentricity orbits (as typical in hierarchical models of galaxy formation) strip LSBs of most of their stars and tend to decrease their surface brightness. On the contrary, bar instability in HSBs leads to substantial loss of angular momentum of the stellar component and to an increase of central surface brightness. In both cases the remnant resembles a spheroidal galaxy with an exponential surface brightness profile. A simple modeling of color evolution and interactions driven star formation gives M/L ratios for the remnants that are roughly consistent with observations. These results suggest an evolutionary scenario for the dwarf galaxies in our Local Group, faint dSphs being the descendents of LSBs and brighter dSphs/dEs being the final state of HSB satellites.

We show that a universe dominated by cold dark matter fails to reproduce the rotation curves of dark matter dominated galaxies, one of the key problems that it was designed to resolve. We perform numerical simulations of the formation of dark matter halos, each containing \gsim 10^6 particles and resolved to 0.003 times the virial radius, allowing an accurate comparison with rotation curve data. A good fit to both galactic and cluster sized halos can be achieved using the density profile rho(r) \propto [(r/r_s)^1.5(1+(r/r_s)^1.5)]^-1, where r_s is a scale radius. This profile has a steeper asymptotic slope, rho(r) \propto r^-1.5, and a sharper turnover than found by lower resolution studies. The central structure of relaxed halos that form within a hierarchical universe has a remarkably small scatter (unrelaxed halos would not host disks). We compare the results with a sample of dark matter dominated, low surface brightness (LSB) galaxies with circular velocities in the range 100-300 km/s. The rotation curves of disks within cold dark matter halos rise too steeply to match these data which require a constant mass density in the central regions. The same conclusion is reached if we compare the scale free shape of observed rotation curves with the simulation data. It is important to confirm these results using stellar rather than HI rotation curves for LSB galaxies. We test the effects of introducing a cut-off in the power spectrum that may occur in a universe dominated by warm dark matter. In this case halos form by a monolithic collapse but the final density profile hardly changes, demonstrating that the merger history does not play a role in determining the halo structure.

Semi-analytical models of galaxy formation based on hierarchical clustering now make a wide range of predictions for observable properties of galaxies at low and high redshift. This article concentrates on 2 aspects: (1) Self-consistent modelling of dust absorption predicts a mean UV extinction A_{UV} ~ 1 mag, depending only weakly on redshift, and similar to observational estimates. (2) The models predict that the Lyman-break galaxies found at z ~ 3 should be strongly clustered with a comoving correlation length r_0 = 4-7 Mpc/h (depending on the cosmology), in good agreement with subsequent observational determinations.

We analyse parallel N-body simulations of three Cold Dark Matter (CDM) universes to study the abundance and clustering of galaxy clusters. The simulations cover a volume comparable to the forthcoming SDSS. We are able to make robust measurements of cluster properties to a redshift larger than unity. We extract halos using two independent, public domain group finders (FOF & HOP) and find consistent results. The correlation function of clusters is in very good agreement with a simple analytic prescription based upon a Lagrangian biasing scheme developed by Mo & White (1996) and the Press-Schechter (PS) formalism for the mass function. The R_0--D_c relation for the open CDM model is in good agreement with the results from the APM Cluster Survey. The SCDM universe shows a robust deviation in the shape and evolution of the mass function when compared with that predicted by the PS formalism. Critical models with a low sigma_8 normalization or small shape parameter Gamma show an excess of massive clusters compared with the PS prediction. When cluster normalized, the SCDM universe at z =1 contains 10 times more clusters with temperatures greater than 7keV, compared with the PS prediction. The agreement between the analytic and N-body mass functions of SCDM can be improved if the value of the delta_c (the extrapolated linear theory threshold for collapse) is revised to be $ delta_c(z) = 1.685[(0.7/sigma_8)(1+z)]^{-0.125}. Our best estimate for the amplitude of fluctuations inferred from the local cluster abundance for SCDM is sigma_{8} = 0.5 \pm 0.04. However, the discrepancy between the temperature function predicted in a critical density universe and that observed at z=0.33 (Henry et al. 1998) remains. (abridged)

The discovery of a population of young galaxies at an epoch when the universe was about one tenth of its current age has shed new light on the question of when and how galaxies formed. Within the context of popular models this is the population of primeval galaxies that built themselves up to the size of present--day galaxies through the process of repeated mergers. But the recent detection of a large concentration of these primeval galaxies (Steidel et al. 1997) appears to be incompatible with hierarchical clustering models, which generally predict that clusters of this size are fully formed later in time. Here we use a combination of two powerful theoretical techniques --semi-analytic modelling and N-body simulations-- to show that such large concentrations should be quite common in a universe dominated by cold dark matter, and that they are the progenitors of the rich galaxy clusters seen today. We predict the clustering properties of primeval galaxies which should, when compared with data that will be collected in the near future, test our understanding of galaxy formation within the framework of a universe dominated by cold dark matter.

We examine the properties of dark matter halos within a rich galaxy cluster using a high resolution simulation that captures the cosmological context of a cold dark matter universe. The mass and force resolution permit the resolution of 150 halos with circular velocities larger than 80 kms within the cluster's virial radius of 2 Mpc. This enables an unprecedented study of the statistical properties of a large sample of dark matter halos evolving in a dense environment. The cumulative fraction of mass attached to these halos varies from 0% at 200 kpc, to 13% at the virial radius. Even at this resolution the overmerging problem persists; halos that pass within 200 kpc of the cluster center are tidally disrupted. Additional substructure is lost at earlier epochs within the massive progenitor halos. The median ratio of apocentric to pericentric radii is 6:1; the orbital distribution is close to isotropic, circular orbits are rare, radial orbits are common. The orbits of halos are unbiased with respect to both position within the cluster and with the orbits of the smooth dark matter background and no velocity bias is detected. The tidal radii of surviving halos are generally well-fit using the simple analytic prediction applied to their orbital pericenters. Halos within clusters have higher concentrations than those in the field. Within the cluster, halo density profiles can be modified by tidal forces and individual encounters with other halos that cause significant mass loss - ``galaxy harassment''. Mergers between halos do not occur inside the clusters virial radius.

We examine the effects of mass resolution and force softening on the density profiles of cold dark matter halos that form within cosmological N-body simulations. As we increase the mass and force resolution, we resolve progenitor halos that collapse at higher redshifts and have very high densities. At our highest resolution we have nearly 3 million particles within the virial radius, several orders of magnitude more than previously used and we can resolve more than one thousand surviving dark matter halos within this single virialised system. The halo profiles become steeper in the central regions and we may not have achieved convergence to a unique slope within the inner 10% of the virialised region. Results from two very high resolution halo simulations yield steep inner density profiles, $\rho(r)\sim r^{-1.4}$. The abundance and properties of arcs formed within this potential will be different from calculations based on lower resolution simulations. The kinematics of disks within such a steep potential may prove problematic for the CDM model when compared with the observed properties of halos on galactic scales.

We take advantage of the biggest cosmological simulation to date for a critical CDM universe in order to test robustness of the cluster mass function on a range of masses much wider than tested before. On the high mass end our results show an excess of hot clusters and a milder evolution compared to the analytical predictions based on the Press & Schechter formula. These features must be properly taken into account in deriving the cluster X-ray temperature function at moderate and high redshifts. On a general basis, the reduced negative evolution in the number of hot clusters could alleviate the discrepancies between the predictions for critical universes and incoming data, which instead seem to favour low Omega_0

We investigate the internal structure of cold dark matter halos using high resolution N-body simulations. As the mass and force resolution are increased, halo density profiles become steeper, asymptoting to a slope of $\sim -4/3$ in the central regions and may not have converged to a unique result. At our highest resolution we have nearly 3 million particles within the virial radius, $R_{200}$, and force softening that is $\sim 0.2% R_{200}$. This resolution has also allowed us to resolve a large part of the overmerging problem - we find over 1000 surviving dark halos orbiting within a single cluster potential. These data have given us unprecedented insights into the dynamics and structure of ``halos within halos'', allowing us for the first time, to compare the distribution of dark matter with observations of galaxies in clusters.

Using cosmological N-body simulations of critical (SCDM) and open (Omega_0=0.3, OCDM) cold dark matter models we have identified dark matter halos which are associated with the progenitors of present day bright early type galaxies. By following their merging history, we show how early type galaxies that formed within massive halos at redshift around 3 are now preferentially residing in clusters and groups. On the other hand, those that formed through major merging events between redshift 1 and the present have not yet been accreted into larger, virialized structures. This result is in agreement with analytical predictions in hierarchical clustering models. CDM models are able to explain both the ancient and uniform population of ellipticals that dominates in clusters together with the more recent and heterogeneous population of field ellipticals. Predictions for the comoving number density of bright early-type galaxies are given, and are shown to be consistent with the observed luminosity function. We predict that the number density of interacting bright binary galaxies, from which the field population of ellipticals may have originated, is proportional to (1+z)**4.2+-0.28 and (1+z)**2.5+-0.42 in SCDM and OCDM respectively. This result is consistent with previous analytical estimates and is discussed together with recent observational constraints.

The dynamics of the Local Group and its environment provide a unique challenge to cosmological models. The velocity field within 5h-1 Mpc of the Local Group (LG) is extremely ``cold''. The deviation from a pure Hubble flow, characterized by the observed radial peculiar velocity dispersion, is measured to be about 60km/s. We compare the local velocity field with similarly defined regions extracted from N-body simulations of Universes dominated by cold dark matter (CDM). This test is able to strongly discriminate between models that have different mean mass densities. We find that neither the Omega=1 (SCDM) nor Omega=0.3 (OCDM) cold dark matter models can produce a single candidate Local Group that is embedded in a region with such small peculiar velocities. For these models, we measure velocity dispersions between 500-700km/s and 150-300km/s respectively, more than twice the observed value. Although both CDM models fail to produce environments similar to those of our Local Group on a scale of a few Mpc, they can give rise to many binary systems that have similar orbital properties as the Milky Way--Andromeda system. The local, gravitationally induced bias of halos in the CDM ``Local Group'' environment, if defined within a sphere of 10 Mpc around each Local Group is about 1.5, independent of Omega. No biasing scheme could reconcile the measured velocity dispersions around Local Groups with the observed one. Identification of binary systems using a halo finder (named Skid (http://www-hpcc.astro.washington.edu/tools/DENMAX for a public version)) based on local density maxima instead of a simple linking algorithm, gives a much more complete sample. We show that a standard ``friend of friends'' algorithm would miss 40% of the LG candidates present in the simulations.

We test the precision with which weak lensing data can provide characteristic cluster mass profiles within Cold Dark Matter scenarios. Using a parallel treecode to simulate volumes as large as 500h-1 Mpc with good resolution, we generate samples of large clusters within a standard CDM model and an open CDM model with Omega_o=0.3. We mock high-quality lensing data by including realistic errors, selecting cluster samples based on velocity dispersion, and fitting profiles within a realistic range in radius. We find that a sample of ten clusters can determine logarithmic profile slopes with 1 sigma errors of about 7%. Increasing the sample size to twenty brings this error down to less than 5%, but this is still insufficient to distinguish the two models. However, measures of cluster profiles obtained with weak lensing do place strong constraints for general CDM--like models of structure formation, and we discuss the optimal strategy for obtaining data samples to use for this purpose.

We study the formation and evolution of small groups of galaxies using new Monte Carlo simulations. These are directly based on the random walk approach to the statistics of condensations collapsing by gravitational instability, and the results are consistent with the Press and Schechter formula. We stress how observational features of groups, such as galaxy membership, depend on the global dynamics of the Universe. This is because the amount of aggregations between infalling galaxies is governed by the ratio of the group age to the time scale for dynamical friction. This ratio strongly depends on cosmology. Thus in low-density models, even if flat, collapse of group sized perturbations followed by fast mergers, would often originate a single massive remnant. This is contrast to observed associations as compact groups of galaxies.

Models of galaxy formation ultimately aim at reproducing the _observed_ properties of galaxies. We report on work in progress to predict luminosities, colours and morphologies of field objects of various masses through N-body + Smoothed Particle Hydrodynamics simulations. We describe our method, illustrate the effects of varying the star formation assumptions, show some preliminary results, draw encouragement from their strengths as well as lessons about what improvements are most needed.

We study the formation and the evolution of galaxy groups in a critical universe, showing the importance of secondary infall for their dynamical evolution. Merging is only slightly accelerated if galaxies have massive halos, because the mass initially associated to the individual galaxies is soon tidally stripped. Stripping is particularly effective for infalling galaxies, which thus avoid merging. We find that, as a rule, merging is effectively terminated when infall becomes dominant. We look for compact groups in our ensemble of simulations, and compare their statistical properties with Hickson's compact groups. We then discuss in terms of the Press & Schechter formalism the statistics of such groups and their evolution in different cosmological scenarios. We show that compact group formation is an ongoing and frequent process in a critical universe. Our model reconciles the contradiction between the observed absence of young merger remnants and the high rate of galaxy interactions expected in compact groups. In open universes, instead, earlier formation of groups and suppression of secondary infall makes it more unlikely that compact groups survive as such until the present. We conclude that the existence of dense and dynamically young groups of galaxies like HCGs points towards a high-density universe.

This paper presents results of a series of numerical N-body experiments that describe the merging of galaxy pairs containing a massive black hole in their core. The aim is to study the orbital evolution of the two black holes through the merger event. It is found that merging does NOT lead to the formation of a CLOSE black hole pair} at the center of the remnant when the two progenitor galaxies have (1) similar masses, (2) different central densities, and (3) relative orbit with non zero angular momentum. Under these conditions, the central black hole of the less dense galaxy settles into an wide orbit, i.e., at a distance from the center comparable to the half mass radius of the remnant, where the time scale for orbital decay by dynamical friction can be several Gyrs. The implications of this result are briefly discussed.