We study the role of group infall in the assembly and dynamics of galaxy clusters in $\Lambda$CDM. We select $10$ clusters with virial mass $M_{\rm 200} \sim 10^{14} \, M_{\odot}$ from the cosmological hydrodynamical simulation Illustris and follow their galaxies with stellar mass $M_{\star} \geq 1.5 \times 10^8 \, M_{\odot}$. A median of $\sim 38\%$ of surviving galaxies at $z=0$ are accreted as part of groups and did not infall directly from the field, albeit with significant cluster-to-cluster scatter. The evolution of these galaxy associations is quick, with observational signatures of their common origin eroding rapidly in $1$-$3$ Gyr after infall. Substructure plays a dominant role in fostering the conditions for galaxy mergers to happen, even within the cluster environment. Integrated over time, we identify (per cluster) an average of $17 \pm 6$ mergers that occur in infalling galaxy associations, of which $7 \pm 3$ occur well within the virial radius of their cluster hosts. The number of mergers show large dispersion from cluster to cluster, with our most massive system having $42$ mergers above our mass cut-off. These mergers, which are typically gas rich for dwarfs and a combination of gas rich and gas poor for $M_{\star} \sim 10^{11} \, M_{\odot}$, may contribute significantly within $\Lambda$CDM to the formation of specific morphologies, such as lenticulars (S0) and blue compact dwarfs in groups and clusters.

We study the assembly of globular clusters (GCs) in 9 galaxy clusters using the cosmological simulation Illustris. GCs are tagged to individual galaxies at their infall time. The tidal removal of GCs from their galaxies and the distribution of the GCs within the cluster is later followed self-consistently by the simulation. The method relies on the simple assumption of a single power-law relation between halo mass (M_vir) and mass in GCs (M_GC) as found in observations. We find that the GCs specific frequency $S_N$ as a function of V-band magnitude naturally reproduces the observed "U"-shape, due to the combination of a power law M_GC-M_vir relation and the non-linear M_*-M_vir relation from the simulation. Additional scatter in the $S_N$ values are traced back to galaxies with early infall times due to the evolution in the M_*-M_vir relation with redshift. GCs that have been tidally removed from their galaxies form today the intra-cluster component from which about ~60% were brought in by galaxies that orbit today within the cluster potential. The remaining "orphan" GCs are contributed by satellite galaxies with a wide range of stellar masses that are fully tidally disrupted at z=0. This intra-cluster component is a good dynamical tracer of the dark matter potential. As a consequence of the accreted nature of most intra-cluster GCs, their orbits are fairly radial with a predicted orbital anisotropy \beta >= 0.5. However, local tangential motions may appear as a consequence of localized substructure, providing a possible interpretation to the \beta<0 values suggested in observations of M87.

One possible and natural derivation from the collisionless cold dark matter (CDM) standard cosmological framework is the assumption of the existence of interactions between dark matter (DM) and photons or neutrinos. Such possible interacting dark matter (IDM) model would imply a suppression of small-scale structures due to a large collisional damping effect, even though the weakly interacting massive particle (WIMP) can still be the DM candidate. Because of this, IDM models can help alleviate alleged tensions between standard CDM predictions and observations at small mass scales. In this work, we investigate the properties of DM halo substructure or subhalos formed in a high-resolution cosmological N-body simulation specifically run within these alternative models. We also run its CDM counterpart, which allowed us to compare subhalo properties in both cosmologies. We show that, in the lower mass range covered by our simulation runs, both subhalo concentrations and abundances are systematically lower in IDM compared to the CDM scenario. Yet, as in CDM, we find that median IDM subhalo concentration values increase towards the innermost regions of their hosts for same mass subhalos. Also similarly to CDM, we find IDM subhalos to be more concentrated than field halos of the same mass. Our work has a direct application on studies aimed at the indirect detection of DM where subhalos are expected to boost the DM signal of their host halos significantly. From our results, we conclude that the role of halo substructure in DM searches will be less important in interacting scenarios than in CDM, but is nevertheless far from being negligible.

We study the formation and evolution of ring galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. We use the largest reference model Ref-L100N1504, a cubic cosmological volume of 100 comoving megaparsecs on a side, to identify and characterise these systems through cosmic time. The number density of ring galaxies in EAGLE is in broad agreement with the observations. The vast majority of ring galaxies identified in EAGLE (83 per cent) have an interaction origin, i.e., form when one or more companion galaxies drop-through a disk galaxy. The remainder (17 per cent) have very long-lived ring morphologies (> 2 Gyr) and host strong bars. Ring galaxies are HI rich galaxies, yet display inefficient star formation activity and tend to reside in the green valley particularly at z > 0.5. This inefficiency is mainly due to the low pressure and metallicity of their interstellar medium (ISM) compared with the ISM of similar star-forming galaxies. We find that the interaction(s) is responsible for decreasing the ISM pressure by causing the ISM gas to flow from the inner regions to the outer disk, where the ring feature forms. At a fixed radius, the star formation efficiency of ring galaxies is indistinguishable from their star-forming counterparts, and thus the main reason for their integrated lower efficiency is the different gas surface density profiles. Since galaxy morphologies are not used to tune the parameters in hydrodynamical simulations, the experiment performed here demonstrates the success of the current numerical models in EAGLE.

Using a cosmological dark matter simulation of a galaxy-cluster halo, we follow the temporal evolution of its globular cluster population. To mimic the red and blue globular cluster populations, we select at high redshift $(z\sim 1)$ two sets of particles from individual galactic halos constrained by the fact that, at redshift $z=0$, they have density profiles similar to observed ones. At redshift $z=0$, approximately 60\% of our selected globular clusters were removed from their original halos building up the intra-cluster globular cluster population, while the remaining 40\% are still gravitationally bound to their original galactic halos. Since the blue population is more extended than the red one, the intra-cluster globular cluster population is dominated by blue globular clusters, with a relative fraction that grows from 60\% at redshift $z=0$ up to 83\% for redshift $z\sim 2$. In agreement with observational results for the Virgo galaxy cluster, the blue intra-cluster globular cluster population is more spatially extended than the red one, pointing to a tidally disrupted origin.

We examine the properties of barred disc galaxies in a LCDM cosmological hydrodynamical simulation from the EAGLE project. Our study follows the formation of 269 discs identified at z = 0 in the stellar mass range 10.6 < log Mstr /M < 11. These discs show a wide range of bar strengths, from unbarred discs to weak bars to strongly barred systems (= 20%). Bars in these systems develop after redshift = 1.3, on timescales that depend sen- sitively on the strength of the pattern. Strong bars develop relatively quickly (in a few Gyr, = 10 disc rotation periods) in systems that are disc dominated, gas poor, and have declining rotation curves. Weak bars develop more slowly in systems where the disc is less gravitation- ally important, and are still growing at z = 0. Unbarred galaxies are comparatively gas-rich discs whose rotation speeds do not exceed the maximum circular velocity of the halos they inhabit. Bar lengths compare favourably with observations, ranging from 0.2 to 0.8 times the radius containing 90% of the stars. Bars slow down remarkably quickly as they grow, causing the inner regions of the surrounding dark halo to expand. At z = 0 strong bars have corota- tion radii roughly ten times the bar length. Such slow bars are inconsistent with the few cases where pattern speeds have been measured or inferred observationally, a discrepancy that, if confirmed, might prove a challenge for disc galaxy formation in LCDM.

The Tully-Fisher relation (TFR) links the stellar mass of a disk galaxy, $M_{\rm str}$, to its rotation speed: it is well approximated by a power law, shows little scatter, and evolves weakly with redshift. The relation has been interpreted as reflecting the mass-velocity scaling ($M\propto V^3$) of dark matter halos, but this interpretation has been called into question by abundance-matching (AM) models, which predict the galaxy-halo mass relation to be non-monotonic and rapidy evolving. We study the TFR of luminous spirals and its relation to AM using the EAGLE set of $\Lambda$CDM cosmological simulations. Matching both relations requires disk sizes to satisfy constraints given by the concentration of halos and their response to galaxy assembly. EAGLE galaxies approximately match these constraints and show a tight mass-velocity scaling that compares favourably with the observed TFR. The TFR is degenerate to changes in galaxy formation efficiency and the mass-size relation; simulations that fail to match the galaxy stellar mass function may fit the observed TFR if galaxies follow a different mass-size relation. The small scatter in the simulated TFR results because, at fixed halo mass, galaxy mass and rotation speed correlate strongly, scattering galaxies along the main relation. EAGLE galaxies evolve with lookback time following approximately the prescriptions of AM models and the observed mass-size relation of bright spirals, leading to a weak TFR evolution consistent with observation out to $z=1$. $\Lambda$CDM models that match both the abundance and size of galaxies as a function of stellar mass have no difficulty reproducing the observed TFR and its evolution.

The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in LCDM simulations that match the size and number of galaxies as a function of stellar mass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. Despite the strong feedback, the scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally-supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ~40 km/s seem to deviate from this prediction. We discuss observational biases and modeling uncertainties that may help to explain this disagreement in the context of LCDM models of dwarf galaxy formation.

We use the high-resolution Aquarius cosmological dark matter simulations coupled to the semi-analytic model by Starkenburg et al. (2013) to study the HI content and velocity width properties of field galaxies at the low mass end in the context of $\Lambda$CDM. We compare our predictions to the observed ALFALFA survey HI mass and velocity width functions, and find very good agreement without fine-tuning, when considering central galaxies. Furthermore, the properties of the dark matter halos hosting galaxies, characterised by their peak velocity and circular velocity at 2 radial disk scalelengths overlap perfectly with the inferred values from observations. This suggests that our galaxies are placed in the right dark matter halos, and consequently at face value, we do not find any discrepancy with the predictions from the $\Lambda$CDM model. Our analysis indicates that previous tensions, apparent when using abundance matching models, arise because this technique cannot be straightforwardly applied for objects with masses $M_{vir} < 10^{10} M_{\odot}$.

We use a cosmological simulation of the formation of the Local Group to explore the origin of age and metallicity gradients in dwarf spheroidal galaxies. We find that a number of simulated dwarfs form "outside-in", with an old, metal-poor population that surrounds a younger, more concentrated metal-rich component, reminiscent of dwarf spheroidals like Sculptor or Sextans. We focus on a few examples where stars form in two populations distinct in age in order to elucidate the origin of these gradients. The spatial distributions of the two components reflect their diverse origin; the old stellar component is assembled through mergers, but the young population forms largely in situ. The older component results from a first episode of star formation that begins early but is quickly shut off by the combined effects of stellar feedback and reionization. The younger component forms when a late accretion event adds gas and reignites star formation. The effect of mergers is to disperse the old stellar population, increasing their radius and decreasing their central density relative to the young population. We argue that dwarf-dwarf mergers offer a plausible scenario for the formation of systems with multiple distinct populations and, more generally, for the origin of age and metallicity gradients in dwarf spheroidals.

Using a cosmological N-body numerical simulation of the formation of a galaxy cluster- sized halo, we analyze the temporal evolution of its globular cluster population. We follow the dynamical evolution of 38 galactic dark matter halos orbiting in a galaxy cluster that at redshift z=0 has a virial mass of 1.71 * 10 ^14 Msol h^-1. In order to mimic both "blue" and "red" populations of globular clusters, for each galactic halo we select two different sets of particles at high redshift (z ~ 1), constrained by the condition that, at redshift z=0, their average radial density profiles are similar to the observed profiles. As expected, the general galaxy cluster tidal field removes a significant fraction of the globular cluster populations to feed the intracluster population. On average, halos lost approximately 16% and 29% of their initial red and blue globular cluster populations, respectively. Our results suggest that these fractions strongly depend on the orbital trajectory of the galactic halo, specifically on the number of orbits and on the minimum pericentric distance to the galaxy cluster center that the halo has had. At a given time, these fractions also depend on the current clustercentric distance, just as observations show that the specific frequencyof globular clusters S_N depends on their clustercentric distance.

We study the radial migration of stars driven by recurring multi-arm spiral features in an exponential disk embedded in a dark matter halo. The spiral perturbations redistribute angular momentum within the disk and lead to substantial radial displacements of individual stars, in a manner that largely preserves the circularity of their orbits and that results, after 5 Gyr (~40 full rotations at the disk scalelength), in little radial heating and no appreciable changes to the vertical or radial structure of the disk. Our results clarify a number of issues related to the spatial distribution and kinematics of migrators. In particular, we find that migrators are a heavily biased subset of stars with preferentially low vertical velocity dispersions. This "provenance bias" for migrators is not surprising in hindsight, for stars with small vertical excursions spend more time near the disk plane and thus respond more readily to non-axisymmetric perturbations. We also find that the vertical velocity dispersion of outward migrators always decreases, whereas the opposite holds for inward migrators. To first order, newly arrived migrators simply replace stars that have migrated off to other radii, thus inheriting the vertical bias of the latter. Extreme migrators might therefore be recognized, if present, by the unexpectedly small amplitude of their vertical excursions. Our results show that migration, understood as changes in angular momentum that preserve circularity, can affect strongly the thin disk, but cast doubts on models that envision the Galactic thick disk as a relic of radial migration.

We explore the impact of cosmic reionization on nearby isolated dwarf galaxies using a compilation of SFHs estimated from deep HST data and a cosmological hydrodynamical simulation of the Local Group. The nearby dwarfs show a wide diversity of star formation histories; from ancient systems that have largely completed their star formation $\sim 10$ Gyr ago to young dwarfs that have formed the majority of their stars in the past $\sim 5$ Gyr to two-component systems characterized by the overlap of comparable numbers of old and young stars. Taken as an ensemble, star formation in nearby dwarfs dips to lower-than-average rates at intermediate times ($4<t$/Gyr $<8$), a feature that we trace in the simulation to the effects of cosmic reionization. Reionization heats the gas and drives it out of the shallow potential wells of low mass halos, affecting especially those below a sharp mass threshold that corresponds to a virial temperature of $\sim 2 \times 10^4 $ $\mathrm{K}$ at $z_{\rm reion}$. The loss of baryons leads to a sharp decline in the star forming activity of early-collapsing systems, which, compounded by feedback from early star formation, empties halos of gas and leaves behind systems where a single old stellar component prevails. In halos below the threshold at $z_{\rm reion}$, reionization heating leads to a delay in the onset of star formation that lasts until the halo grows massive enough to allow some of the remaining gas to cool and form stars. Young stellar components therefore dominate in dwarfs whose halos assemble late and thus form few stars before reionization. Two-component systems may be traced to late mergers of individual examples of the two aforementioned cases. The relative dearth of intermediate-age stars in nearby dwarfs might thus be the clearest signature yet identified of the imprint of cosmic reionization on the star formation history of dwarf galaxies.

Counterrotating stars in disk galaxies are a puzzling dynamical feature whose origin has been ascribed to either satellite accretion events or to disk instabilities triggered by deviations from axisymmetry. We use a cosmological simulation of the formation of a disk galaxy to show that counterrotating stellar disk components may arise naturally in hierarchically-clustering scenarios even in the absence of merging. The simulated disk galaxy consists of two coplanar, overlapping stellar components with opposite spins: an inner counterrotating bar-like structure made up mostly of old stars surrounded by an extended, rotationally-supported disk of younger stars. The opposite-spin components originate from material accreted from two distinct filamentary structures which at turn around, when their net spin is acquired, intersect delineating a "V"-like structure. Each filament torques the other in opposite directions; the filament that first drains into the galaxy forms the inner counterrotating bar, while material accreted from the other filament forms the outer disk. Mergers do not play a substantial role and most stars in the galaxy are formed in situ; only 9% of all stars are contributed by accretion events. The formation scenario we describe here implies a significant age difference between the co- and counterrotating components, which may be used to discriminate between competing scenarios for the origin of counterrotating stars in disk galaxies.

We use a cosmological simulation of the formation of the Local Group of Galaxies to identify a mechanism that enables the removal of baryons from low-mass halos without appealing to feedback or reionization. As the Local Group forms, matter bound to it develops a network of filaments and pancakes. This moving web of gas and dark matter drifts and sweeps a large volume, overtaking many halos in the process. The dark matter content of these halos is unaffected but their gas can be efficiently removed by ram-pressure. The loss of gas is especially pronounced in low-mass halos due to their lower binding energy and has a dramatic effect on the star formation history of affected systems. This "cosmic web stripping" may help to explain the scarcity of dwarf galaxies compared with the numerous low-mass halos expected in \Lambda CDM and the large diversity of star formation histories and morphologies characteristic of faint galaxies. Although our results are based on a single high-resolution simulation, it is likely that the hydrodynamical interaction of dwarf galaxies with the cosmic web is a crucial ingredient so far missing from galaxy formation models.

The cold dark matter halo mass function is much steeper than the galaxy stellar mass function on galactic and subgalactic scales. This difference is usually reconciled by assuming that the galaxy formation efficiency drops sharply with decreasing halo mass, so that virtually no dwarf galaxies form in halos less massive than ~ 10^10 M_sun. In turn, this implies that, at any given radius, the dark mass enclosed by a galaxy must exceed a certain minimum. We use rotation curves of dwarf galaxies compiled from the literature to explore whether their enclosed mass is consistent with these constraints. We find that almost one half of the dwarfs in our sample with stellar mass between 10^6-10^7 M_sun are at odds with this restriction: either they live in halos with masses substantially below 10^10 M_sun or there is a mechanism capable of reducing the dark mass enclosed by some of the faintest dwarfs. Neither possibility is easily accommodated within the standard LCDM scenario. Extending galaxy formation to halos well below 10^10 M_sun would lead to large numbers of dwarf galaxies in excess of current estimates; at the same time, the extremely low stellar mass of the systems involved makes it unlikely that baryonic effects can reduce their dark matter content. Resolving this challenge seems to require new insights into dwarf galaxy formation, or perhaps a radical revision of the prevailing paradigm.

We perform controlled numerical experiments to asses the effect of baryon mass loss on the inner structure of large galactic dark matter haloes. This mass expulsion is intended to mimic both the supernovae and AGN feedbacks, as well as the evolution of stellar populations. This study is meant in particular for precursors of massive Early Type Galaxies, wherein strong AGN feedback (often dubbed "QSO mode" in galaxy formation models) has been proposed to remove on a short timescale, of the order of a few dynamical times, a substantial fraction of their baryons. In a previous paper we evaluated the observational consequences (size increase) of this process on the galactic structure (Ragone-Figueroa & Granato 2011). Here we focus on the distribution of dark matter in the galactic region. It is shown that the inner region of the DM halo expands and its density profile flattens by a sizeable amount, with little dependence on the expulsion timescale. We also evaluate the effect of the commonly made approximation of treating the baryonic component as a potential that changes in intensity without any variation in shape. This approximation leads to some underestimates of the halo expansion and its slope flattening. We conclude that cuspy density profiles in ETGs could be difficult to reconcile with an effective AGN (or stellar) feedback during the evolution of these systems.

The Monoceros ring is a collection of stars in nearly-circular orbits at roughly 18 kpc from the Galactic center. It may have originated (i) as the response of the disc to perturbations excited by satellite companions or (ii) from the tidal debris of a disrupted dwarf galaxy. The metallicity of Monoceros stars differs from that of disc stars at comparable Galactocentric distances, an observation that disfavours the first scenario. On the other hand, circular orbits are difficult to accommodate in the tidal-disruption scenario, since it requires a satellite which at the time of disruption was itself in a nearly circular orbit. Such satellite could not have formed at the location of the ring and, given its low mass, dynamical friction is unlikely to have played a major role in its orbital evolution. We search cosmological simulations for low-mass satellites in nearly-circular orbits and find that they result, almost invariably, from orbital changes induced by collisions with more massive satellites: the radius of the circular orbit thus traces the galactocentric distance of the collision. Interestingly, the Sagittarius dwarf, one of the most luminous satellites of the Milky Way, is in a polar orbit that crosses the Galactic plane at roughly the same Galactocentric distance as Monoceros. We use idealized simulations to demonstrate that an encounter with Sagittarius might well have led to the circularization and subsequent tidal demise of the progenitor of the Monoceros ring.

We search for the existence of chemically-distinct stellar components in the solar neighbourhood using published data. Extending earlier work, we show that when the abundances of Fe, alpha elements, and the r-process element Eu are considered, stars separate neatly into two groups that delineate the traditional thin and thick disk of the Milky Way. The group akin to the thin disk is traced by stars with [Fe/H]>-0.7 and alpha/Fe<0.2. The thick disk-like group overlaps the thin disk in [Fe/H] but has higher abundances of alpha elements and Eu. Stars in the range -1.5<[Fe/H]<-0.7 with low [alpha/Fe] ratios, however, seem to belong to a separate, dynamically-cold, non-rotating component that we associate with tidal debris. The kinematically-hot stellar halo dominates the sample for [Fe/H]<-1.5. These results suggest that it may be possible to define the main dynamical components of the solar neighbourhood using only their chemistry, an approach with a number of interesting consequences. The kinematics of thin disk stars is then independent of metallicity: their average rotation speed remains roughly constant in the range -0.7<[Fe/H]<+0.4, a result that argues against radial migration having played a substantial role in the evolution of the thin disk. The velocity dispersion of stars assigned to the thin disk is also independent of [Fe/H], implying that the familiar increase in velocity dispersion with decreasing metallicity is the result of the increased prevalence of the thick disk at lower metallicities, rather than of the sustained operation of a dynamical heating mechanism. The substantial overlap in [Fe/H] and, probably, stellar age, of the various components might affect other reported trends in the properties of stars in the solar neighbourhood.

We develop a new numerical approach to describe the action of ram pressure stripping (RPS) within a semi-analytic model of galaxy formation and evolution which works in combination with non-radiative hydrodynamical simulations of galaxy clusters. The new feature in our method is the use of the gas particles to obtain the kinematical and thermodynamical properties of the intragroup and intracluster medium (ICM). This allows a self-consistent estimation of the RPS experienced by satellite galaxies. We find that the ram pressure (RP) in the central regions of clusters increases approximately one order of magnitude between z = 1 and 0, consistent with the increase in the density of the ICM. The mean RP experienced by galaxies within the virial radius increases with decreasing redshift. In clusters with virial masses M_vir ~10^15 h^-1 M_Sun, over 50 per cent of satellite galaxies have experienced RP ~10^(-11) h^2 dyn cm^-2 or higher for z <= 0.5. In smaller clusters (M_vir ~10^14 h^-1 M_Sun) the mean RP are approximately one order of magnitude lower at all redshifts. RPS has a strong effect on the cold gas content of galaxies for all cluster masses. At z = 0, over 70 per cent of satellite galaxies within the virial radius are completely depleted of cold gas. For the more massive clusters the fraction of depleted galaxies is already established at z ~ 1, whereas for the smaller clusters this fraction increases appreciably between z = 1 and 0. This indicates that the rate at which the cold gas is stripped depends on the virial mass of the host cluster. Compared to our new approach, the use of an analytic profile to describe the ICM results in an overestimation of the RP larger than 50 per cent for z >= 0.5.

Cosmological simulations indicate that cold dark matter (CDM) halos should be triaxial. Verifying observationally this theoretical prediction is, however, less than straightforward because the assembly of galaxies is expected to modify the halo shapes and to render them more axisymmetric. We use a suite of N-body simulations to investigate quantitatively the effect of the growth of a central disk galaxy on the shape of triaxial dark matter halos. As expected, the halo responds to the presence of the disk by becoming more spherical. The net effect depends only weakly on the orientation of the disk relative to the halo principal axes or the timescale of disk assembly, but strongly on the overall gravitational importance of the disk. Our results show that exponential disks whose contribution peaks at less than ~50% of their circular velocity are unable to modify noticeably the shape of the gravitational potential of their surrounding halos. Many dwarf and low surface brightness galaxies are expected to be in this regime, and therefore their detailed kinematics could be used to probe halo triaxiality, one of the basic predictions of the CDM paradigm. We argue that the complex disk kinematics of the dwarf galaxy NGC 2976 might be the reflection of a triaxial halo. Such signatures of halo triaxiality should be common in galaxies where the luminous component is subdominant.

We study the orbital properties of stars in four (published) simulations of thick disks formed by: i) accretion from disrupted satellites, ii) heating of a pre-existing thin disk by a minor merger, iii) radial migration and iv) gas rich mergers. We find that the distribution of orbital eccentricities are predicted to be different for each model: a prominent peak at low eccentricity is expected for the heating, migration and gas-rich merging scenarios, while the eccentricity distribution is broader and shifted towards higher values for the accretion model. These differences can be traced back to whether the bulk of the stars in each case is formed 'in-situ' or is 'accreted', and are robust to the peculiarities of each model. A simple test based on the eccentricity distribution of nearby thick disk stars may thus help elucidate the dominant formation mechanism of the Galactic thick disk.

We use N-body/gasdynamical LambdaCDM cosmological simulations to examine the effect of the assembly of a central galaxy on the shape and mass profile of its dark halo. Two series of simulations are compared; one that follows only the evolution of the dark matter component and a second one where a baryonic component is added. These simulations include radiative cooling but neglect star formation and feedback, leading most baryons to collect at the halo center in a disk which is too small and too massive when compared with typical spiral. This unrealistic model allows us, nevertheless, to gauge the maximum effect that galaxies may have in transforming their dark halos. We find that the shape of the halo becomes more axisymmetric: halos are transformed from triaxial into essentially oblate systems, with well-aligned isopotential contours of roughly constant flattening (c/a ~ 0.85). Halos always contract as a result of galaxy assembly, but the effect is substantially less pronounced than predicted by the "adiabatic contraction" hypothesis. The reduced contraction helps to reconcile LambdaCDM halos with constraints on the dark matter content inside the solar circle and should alleviate the long-standing difficulty of matching simultaneously the scaling properties of galaxy disks and the luminosity function. The halo contraction is also less pronounced than found in earlier simulations, a disagreement that suggests that halo contraction is not solely a function of the initial and final distribution of baryons. Not only how much baryonic mass has been deposited at the center of a halo matters, but also the mode of its deposition. It might prove impossible to predict the halo response without a detailed understanding of a galaxy's assembly history. (Abriged)

Halo stars with unusually high radial velocity ("hypervelocity" stars, or HVS) are thought to be stars unbound to the Milky Way that originate from the gravitational interaction of stellar systems with the supermassive black hole at the Galactic center. We examine the latest HVS compilation and find peculiarities that are unexpected in this black hole-ejection scenario. For example, a large fraction of HVS cluster around the constellation of Leo and share a common travel time of $\sim 100$-200 Myr. Furthermore, their velocities are not really extreme if, as suggested by recent galaxy formation models, the Milky Way is embedded within a $2.5\times 10^{12} h^{-1} M_{\odot}$ dark halo with virial velocity of $\sim 220$ km/s. In this case, the escape velocity at $\sim 50$ kpc would be $\sim 600$ km/s and very few HVS would be truly unbound. We use numerical simulations to show that disrupting dwarf galaxies may contribute halo stars with velocities up to and sometimes exceeding the nominal escape speed of the system. These stars are arranged in a thinly-collimated outgoing ``tidal tail'' stripped from the dwarf during its latest pericentric passage. We speculate that some HVS may therefore be tidal debris from a dwarf recently disrupted near the center of the Galaxy. In this interpretation, the angular clustering of HVS results because from our perspective the tail is seen nearly ``end on'', whereas the common travel time simply reflects the fact that these stars were stripped simultaneously from the dwarf during a single pericentric passage. This proposal is eminently falsifiable, since it makes a number of predictions that are distinct from the black-hole ejection mechanism and that should be testable with improved HVS datasets.

We examine the orbits of satellite galaxies identified in a suite of N-body/gasdynamical simulations of the formation of $L_*$ galaxies in a LCDM universe. Most satellites follow conventional orbits; after turning around, they accrete into their host halo and settle on orbits whose apocentric radii are steadily eroded by dynamical friction. However, a number of outliers are also present, we find that ~1/3 of satellites identified at $z=0$ are on unorthodox orbits, with apocenters that exceed their turnaround radii. This population of satellites on extreme orbits consists typically of the faint member of a satellite pair that has been ejected onto a highly-energetic orbit during its first approach to the primary. Since the concurrent accretion of multiple satellite systems is a defining feature of hierarchical models of galaxy formation, we speculate that this three-body ejection mechanism may be the origin of (i) some of the newly discovered high-speed satellites around M31 (such as Andromeda XIV); (ii) some of the distant fast-receding Local Group members, such as Leo I; and (iii) the oddly isolated dwarf spheroidals Cetus and Tucana in the outskirts of the Local Group. Our results suggest that care must be exercised when using the orbits of the most weakly bound satellites to place constraints on the total mass of the Local Group.

We study the population of satellite galaxies formed in a suite of N-body/gasdynamical simulations of galaxy formation in a LCDM universe. We find little spatial or kinematic bias between the dark matter and the satellite population. The velocity dispersion of the satellites is a good indicator of the virial velocity of the halo: \sigma_{sat}/V_{vir}=0.9 +/- 0.2. Applied to the Milky Way and M31 this gives V_{vir}^{MW}=109 +/- 22$ km/s and V_{vir}^{M31} = 138 +/- 35 km/s, respectively, substantially lower than the rotation speed of their disk components. The detailed kinematics of simulated satellites and dark matter are also in good agreement. By contrast, the stellar halo of the simulated galaxies is kinematically and spatially distinct from the population of surviving satellites. This is because the survival of a satellite depends on mass and on time of accretion; surviving satellites are biased toward low-mass systems that have been recently accreted by the galaxy. Our results support recent proposals for the origin of the systematic differences between stars in the Galactic halo and in Galactic satellites: the elusive ``building blocks'' of the Milky Way stellar halo were on average more massive, and were accreted (and disrupted) earlier than the population of dwarfs that has survived self-bound until the present.

We search for signatures of past accretion events in the Milky Way in the recently published catalogue by Nordstrom et al.(2004), containing accurate spatial and kinematic information as well as metallicities for 13240 nearby stars. To optimize our strategy, we use numerical simulations to characterize the properties of the debris from disrupted satellites. We find that stars with a common progenitor should show distinct correlations between their orbital parameters; in particular, between the apocentre A and pericentre P, as well as their z-angular momentum (L_z). In the APL-space, such stars are expected to cluster around regions of roughly constant eccentricity. The APL space for the Nordstrom catalogue exhibits a wealth of substructure, much of which can be linked to dynamical perturbations induced by spiral arms and the Galactic bar. However, our analysis also reveals a statistically significant excess of stars on orbits of common (moderate) eccentricity, analogous to the pattern expected for merger debris. Besides being dynamically peculiar, the 274 stars in these substructures have very distinct metallicity and age distributions, providing further evidence of their extra-Galactic provenance. It is possible to identify among these stars, three coherent Groups with characteristic ages and metallicities, that, in all likelihood, correspond to the remains of disrupted satellites. [abridged]

(Abridged) We use numerical simulations to investigate the origin and structure of the luminous halos that surround isolated galaxies. These stellar structures extend out to several hundred kpc away from a galaxy, and consist of stars shed by merging subunits during the many accretion events that characterize the hierarchical assembly of galaxies. Such origin suggests that outer luminous halos are ubiquitous and that they should appear as an excess of light over extrapolations of the galaxy's inner profile beyond its traditional luminous radius. The mass profile of the accreted stellar component is well approximated by a model where the logarithmic slope steepens monotonically with radius; from -3 at the luminous edge of the galaxy to -4 or steeper near the virial radius of the system. Such spatial distribution is consistent with that of Galactic and M31 globular clusters, suggesting that many of the globulars were brought in by accretion events, in a manner akin to the classic Searle-Zinn scenario. The outer stellar spheroid is supported by a velocity dispersion tensor with a substantial and radially increasing radial anisotropy. These properties distinguish the stellar halo from the dark matter component, which is more isotropic in velocity space, as well as from some tracers of the outer spheroid such as satellite galaxies. Most stars in the outer halo formed in progenitors that have since merged with the central galaxy; very few stars in the halo are contributed by satellites that survive as self-bound entities at the present. These features are in reasonable agreement with recent observations of the outer halo of the MW, of M31, and of other isolated spirals, and suggest that all of these systems underwent an early period of active merging, as envisioned in hierarchical models of galaxy formation.

We use numerical simulations to investigate the orbital characteristics of tidal debris from satellites whose orbits are dragged into the plane of galactic disks by dynamical friction before disruption. We find that these satellites may deposit a significant fraction of their stars into the disk components of a galaxy, and use our results to motivate the search for accretion relicts in samples of metal-poor disk stars in the vicinity of the Sun. Satellites disrupted on very eccentric orbits coplanar with the disk are expected to shed stars in "trails" of distinct orbital energy and angular momentum during each pericentric passage. To an observer located between the pericenter and apocenter of such orbits, these trails would show as distinct groupings of stars with low vertical velocity and a broad, symmetric, often double-peaked distribution of Galactocentric radial velocities. One group of stars with these characteristics stands out in available compilations of nearby metal-poor stars. These stars have specific angular momenta similar to that of the globular cluster omegaCen, long hypothesized to be the nucleus of a dwarf galaxy disrupted by the Milky Way tidal field. In addition to their kindred kinematics, stars in the omegaCen group share distinct chemical abundance characteristics, and trace a well-defined track in the [alpha/Fe] versus [Fe/H] plane, consistent with simple closed-box enrichment models and a protracted star formation history. The dynamical and chemical coherence of this group suggests that it consists of stars that once belonged to the dwarf that brought omegaCen into the Galaxy. The presence of this and other "tidal relicts" in the solar neighbourhood suggest an extra-Galactic origin for the presence of nearby stars with odd kinematics and chemistry.

We use numerical simulations to investigate the orientation of the angular momentum axis of disk galaxies relative to their surrounding large scale structure. We find that this is closely related to the spatial configuration at turnaround of the material destined to form the galaxy, which is often part of a coherent two-dimensional slab criss-crossed by filaments. The rotation axis is found to align very well with the intermediate principal axis of the inertia momentum tensor at this time. This orientation is approximately preserved during the ensuing collapse, so that the rotation axis of the resulting disk ends up lying on the plane traced by the protogalactic material at turnaround. This suggests a tendency for disks to align themselves so that their rotation axis is perpendicular to the minor axis of the structure defined by surrounding matter. One example of this trend is provided by our own Galaxy, where the Galactic plane is almost at right angles with the supergalactic plane (SGP) drawn by nearby galaxies; indeed, the SGP latitude of the North Galactic Pole is just 6 degrees. We have searched for a similar signature in catalogs of nearby disk galaxies, and find a significant excess of edge-on spirals (for which the orientation of the disk rotation axis may be determined unambiguously) highly inclined relative to the SGP. This result supports the view that disk galaxies acquire their angular momentum as a consequence of early tidal torques acting during the expansion phase of the protogalactic material.

We present a detailed analysis of the dynamical properties of a simulated disk galaxy assembled hierarchically in the Lambda CDM cosmogony. At z=0, two distinct dynamical components are identified on the basis of the orbital parameters of stars in the galaxy: a slowly rotating, centrally concentrated spheroid and a disk-like component largely supported by rotation. These components are also recognized in the surface brightness profile of the galaxy, which can be very well approximated by the superposition of an R^{1/4} spheroid and an exponential disk. However, neither does the dynamically-identified spheroid follow de Vaucouleurs' law nor is the disk purely exponential, a resultwhich calls for caution when estimating the importance of the disk from traditional photometric decomposition techniques. The disk may be further decomposed into a thin, dynamically cold component with stars on nearly circularorbits and a hotter, thicker component with orbital parameters transitional between the thin disk and the spheroid. The spheroid is old, and has essentiallyno stars younger than the time elapsed since the last major accretion event ~8 Gyr ago. The majority of thin disk stars, form after the merging activity is over, although a significant fraction ~15% of thin-disk stars are old enough to predate the last major merger event. This unexpected population of old disk stars consists mainly of the tidal debris of satellites whose orbital plane was coincident with the disk and whose orbits were circularized by dynamical friction prior to full disruption. Our results highlight the role of satellite accretion events in shaping the disk and the spheroidal component and reveal some of the clues to the assembly process of a galaxy preserved in the detailed dynamics of old stellar populations.

We analyze the properties of a disk galaxy simulated with unprecedented numerical resolution in the Lambda CDM cosmogony. The galaxy is assembled through a number of high-redshift mergers followed by a period of quiescent accretion after z~1 which lead to the formation of two distinct dynamical components: a spheroid and a disk. The surface brightness profile is very well approximated by the superposition of an R^{1/4} spheroid and an exponential disk. The surface brightness profile is remarkably similar to that of Sab galaxy UGC615, but the simulated galaxy rotates significantly faster and has a declining rotation curve dominated by the spheroid near the center. The decline in circular velocity is at odds with observation and results from the high concentration of the dark matter and baryonic components, as well as from the relatively high mass-to-light ratio of the stars in the simulation. The simulated galaxy lies ~ 1 mag off the I-band Tully-Fisher relation of late-type spirals, but seems to be in reasonable agreement with Tully-Fisher data on S0 galaxies. The angular momentum of the luminous component is an order of magnitude lower than that of late-type spirals of similar rotation speed. This reflects the dominance of the slowly-rotating, dense spheroidal component. The disk component, on the other hand, has properties rather similar to those of late-type spirals. This suggests that a different form of feedback than adopted in this simulation is required to inhibit the efficient collapse and cooling of gas at high redshift that leads to the formation of the spheroid. Reconciling the properties of disk galaxies with the early collapse and high merging rates characteristic of hierarchical scenarios such as Lambda CDM remains a challenging, yet so far elusive, proposition.

We compute the redshift space power spectrum of two X-ray cluster samples: the X-ray Brightest Abell Cluster Sample (XBACS) and the Brightest Cluster Sample (BCS) using the method developed by Feldman, Kaiser & Peacock. The power spectrums derived for these samples are in agreement with determinations of other optical and X-ray cluster samples. For XBACS we find the largest power spectrum amplitude expected given the high richness of this sample ($R \ge 2$). In the range $0.05 \uk < k < 0.4 \uk$ the power spectrum shows a power law behavior $P(k)\propto k^{n}$ with an index $n\simeq-1.2$. In a similar range $0.04 \uk < k < 0.3 \uk $ BCS power spectrum has a smaller amplitude with index $n\simeq-1.0$. We do not find significant evidence for a peak at $k \simeq 0.05 \uk$ suggesting that claims such of feature detections in some cluster samples could relay on artificial inhomogeneities of the data. We compare our results with power spectrum predictions derived by Moscardini et al. within current cosmological models (LCDM and OCDM). For XBACS we find that both models underestimate the amplitude of the power spectrum but for BCS there is reasonably good agreement at $k\gsim 0.03 \uk$ for both models.

We describe similarity solutions that characterize the collapse of collisional gas onto scalefree perturbations in an Einstein-de Sitter universe. We consider the effects of radiative cooling and derive self-similar solutions under the assumption that the cooling function is a power law of density and temperature, $\Lambda(T,\rho) \propto \rho^{3/2} T$. We use these results to test the ability of Smooth Particle Hydrodynamics (SPH) techniques to follow the collapse and accretion of shocked, rapidly cooling gas in a cosmological context. Our SPH code reproduces the analytical results very well in cases that include or exclude radiative cooling. No substantial deviations from the predicted central mass accretion rates or from the temperature, density, and velocity profiles are observed in well resolved regions inside the shock radius. This test problem lends support to the reliability of SPH techniques to model the complex process of galaxy formation.

We use 3-dimensional SPH/N-BODY simulations to study ram pressure stripping of gas from spiral galaxies orbiting in clusters. We find that the analytic expectation of Gunn & Gott (1972) relating the gravitational restoring force provided by the disk to the ram pressure force, provides a good approximation to the radius that gas will be stripped from a galaxy. However, at small radii it is also important to consider the potential provided by the bulge component. A spiral galaxy passing through the core of a rich cluster such as Coma, will have its gaseous disk truncated to $\sim 4$ kpc, thus losing $\sim 80%$ of its diffuse gas mass. The timescale for this to occur is a fraction of a crossing time $\sim 10^7$ years. Galaxies orbiting within poorer clusters, or inclined to the direction of motion through the intra-cluster medium will lose significantly less gas. We conclude that ram-pressure alone is insufficient to account for the rapid and widespread truncation of star-formation observed in cluster galaxies, or the morphological transformation of Sab's to S0's that is necessary to explain the Butcher-Oemler effect.

We compute the cluster auto-correlation function $\xi_{cc}(r)$ of an X-ray flux limited sample of Abell clusters (XBACs, \cite{ebe}). For the total XBACs sample we find a power-law fit $\xi_{cc}=(r/r_0)^{\gamma}$ with $r_0=21.1$ Mpc h$^{-1}$and $\gamma =-1.9$ consistent with the results of $R \ge 1 $ Abell clusters. We also analyze $\xi_{cc}(r)$ for subsamples defined by different X-ray luminosity thresholds where we find a weak tendency of larger values of $r_0$ with increasing X-ray luminosity although with a low statistical significance. In the different subsamples analyzed we find $21 < r_0 < 35 $ Mpc h$^{-1}$ and $-1.9< \gamma < -1.6$. Our analysis suggests that cluster X-ray luminosities may be used for a reliable confrontation of cluster spatial distribution properties in models and observations.

We study the ability of COBE-normalized CDM models to reproduce observed properties of the distribution of galaxies and clusters using N-body numerical simulations. We analyze the galaxy-galaxy and cluster-galaxy two-point correlation functions, $\xi_{gg}$ and $\xi_{cg}$, in open ($\Omega_{0}=0.4, \Omega_{\Lambda}=0, \sigma_8=0.75$), and flat ($\Omega_{0}=0.3, \Omega_{\Lambda}=0.7, \sigma_8=1.05$) CDM models which both reproduce the observed abundances of rich clusters of galaxies. To compare models with observations we compute projected cross-correlation functions $\omega_{gg}$ and $\omega_{cg}$ to derive the corresponding $\xi_{gg}$ and $\xi_{cg}$. We use target galaxies selected from Las Campanas Redshift Survey, target clusters selected from the APM Cluster Survey and tracer galaxies from the Edinburgh Durham Sky Survey catalog.

We present a hydrodynamical code based on the Smooth Particle Hydrodynamics technique implemented in an AP3M code aimed at solving the hydrodynamical and gravitational equations in a cosmological frame. We analyze the ability of the code to reproduce standard tests and perform numerical simulations to study the formation of galaxies in a typical region of a CDM model. These numerical simulations include gas and dark matter particles and take into account physical processes such as shock waves, radiative cooling, and a simplified model of star formation. Several observed properties of normal galaxies such as $M_{gas}/M_{total}$ ratios, the luminosity function and the Tully-Fisher relation are analyzed within the limits imposed by numerical resolution.