It has been claimed that the standard model of cosmology (LCDM) cannot easily account for a number of observations on relatively small scales, motivating extensions to the standard model. Here we introduce a new suite of cosmological simulations that systematically explores three plausible extensions: warm dark matter, self-interacting dark matter, and a running of the scalar spectral index of density fluctuations. Current observational constraints are used to specify the additional parameters that come with these extensions. We examine a large range of observable metrics on small scales, including the halo mass function, density and circular velocity profiles, the abundance of satellite subhaloes, and halo concentrations. For any given metric, significant degeneracies can be present between the extensions. In detail, however, the different extensions have quantitatively distinct mass and radial dependencies, suggesting that a multi-probe approach over a range of scales can be used to break the degeneracies. We also demonstrate that the relative effects on the radial density profiles in the different extensions (compared to the standard model) are converged down to significantly smaller radii than are the absolute profiles. We compare the derived cosmological trends with the impact of baryonic physics using the EAGLE and ARTEMIS simulations. Significant degeneracies are also present between baryonic physics and cosmological variations (with both having similar magnitude effects on some observables). Given the inherent uncertainties both in the modelling of galaxy formation physics and extensions to LCDM, a systematic and simultaneous exploration of both is strongly warranted.

Recent analyses of the cosmic microwave background (CMB) and the Lyman-alpha forest indicate a mild preference for a deviation from a power law primordial matter power spectrum (a so-called negative `running'). We use an extension to the BAHAMAS suite of cosmological hydrodynamic simulations to explore the effects that a running scalar spectral index has on large-scale structure (LSS), using Planck CMB constraints to initialize the simulations. We focus on 5 key statistics: i) the non-linear matter power spectrum ii) the halo mass function; iii) the halo two-point auto correlation function; iv) total mass halo density profiles; and v) the halo concentration-mass relation. In terms of the matter power spectrum, we find that a running scalar spectral index affects all k-scales examined in this study, with a negative (positive) running leading to an amplification (suppression) of power. These effects should be easily detectable with upcoming surveys such as LSST and Euclid. In the mass range sampled, a positive running leads to an increase in the mass of galaxy groups and clusters, with the favoured negative running leading to a decrease in mass of lower-mass (M <~ 10^13 M_solar) halos, but an increase for the most massive (M >~ 10^13 M_solar) halos. Changes in the mass are generally confined to 5-10% which, while not insignificant, cannot by itself reconcile the claimed tension between the primary CMB and cluster number counts. We find that running does not significantly affect the shapes of density profiles of matched halos, changing only their amplitude. Finally, we demonstrate that the observed effects on LSS due to a running scalar spectral index are separable from those of baryonic effects to typically a few percent precision.

Dark matter (DM) direct detection experiments aim to place constraints on the DM--nucleon scattering cross-section and the DM particle mass. These constraints depend sensitively on the assumed local DM density and velocity distribution function. While astrophysical observations can inform the former (in a model-dependent way), the latter is not directly accessible with observations. Here we use the high-resolution ARTEMIS cosmological hydrodynamical simulation suite of 42 Milky Way-mass halos to explore the spatial and kinematical distributions of the DM in the solar neighbourhood, and we examine how these quantities are influenced by substructures, baryons, the presence of dark discs, as well as general halo-to-halo scatter (cosmic variance). We also explore the accuracy of the standard Maxwellian approach for modelling the velocity distribution function. We find significant halo-to-halo scatter in the density and velocity functions which, if propagated through the standard halo model for predicting the DM detection limits, implies a significant scatter about the typically quoted limit. We also show that, in general, the Maxwellian approximation works relatively well for simulations that include the important gravitational effects of baryons, but is less accurate for collisionless (DM-only) simulations. Given the significant halo-to-halo scatter in quantities relevant for DM direct detection, we advocate propagating this source of uncertainty through in order to derive conservative DM detection limits.

The standard model of cosmology, the LCDM model, robustly predicts the existence of a multitude of dark matter 'subhaloes' around galaxies like the Milky Way. A wide variety of observations have been proposed to look for the gravitational effects such subhaloes would induce in observable matter. Most of these approaches pertain to the stellar or cool gaseous phases of matter. Here we propose a new approach, which is to search for the perturbations that such dark subhaloes would source in the warm/hot circumgalactic medium (CGM) around normal galaxies. With a combination of analytic theory, carefully-controlled high-resolution idealised simulations, and full cosmological hydrodynamical simulations, we calculate the expected signal and how it depends on important physical parameters (subhalo mass, CGM temperature, and relative velocity). We find that dark subhaloes enhance the local CGM pressure, density, and temperature, in order of decreasing magnitude of the effect. For the pressure, the fluctuations can vary in magnitude from tens of percent (for subhaloes with M_sub=10^10 Msun) to a few percent (for subhaloes with M_sub=10^8 Msun), although this depends strongly on the CGM temperature. The subhaloes also induce fluctuations in the velocity field ranging in magnitude from a few km/s up to 25 km/s. We propose that X-ray, Sunyaev-Zel'dovich effect, radio dispersion measure, and quasar absorption line observations can be used to measure these fluctuations and place constraints on the abundance and distribution of dark subhaloes, thereby placing constraints on the nature of dark matter.

A large body of work based on collisionless cosmological N-body simulations going back over two decades has advanced the idea that collapsed dark matter haloes have simple and approximately universal forms for their mass density and pseudo-phase space density (PPSD) distributions. However, a general consensus on the physical origin of these results has not yet been reached. In the present study, we explore to what extent the apparent universality of these forms holds when we vary the initial conditions (i.e., the primordial power spectrum of density fluctuations) away from the standard CMB-normalised case, but still within the context of LCDM with a fixed expansion history. Using simulations that vary the initial amplitude and shape, we show that the structure of dark matter haloes retains a clear memory of the initial conditions. Specifically, increasing (lowering) the amplitude of fluctuations increases (decreases) the concentration of haloes and, if pushed far enough, the density profiles deviate strongly from the NFW form that is a good approximation for the CMB-normalised case. Although, an Einasto form works well. Rather than being universal, the slope of the PPSD (or pseudo-entropy) profile steepens (flattens) with increasing (decreasing) power spectrum amplitude and can exhibit a strong halo mass dependence. Our results therefore indicate that the previously identified universality of the structure of dark matter haloes is mostly a consequence of adopting a narrow range of (CMB-normalised) initial conditions for the simulations. Our new suite provides a useful test-bench against which physical models for the origin of halo structure can be validated.

The observable properties of galaxies depend on both internal processes and the external environment. In terms of the environmental role, we still do not have a clear picture of the processes driving the transformation of galaxies. The use of proxies for environment (e.g., host halo mass, distance to the N^th nearest neighbour, etc.), as opposed to the real physical conditions (e.g., hot gas density) may bear some responsibility for this. Here we propose a new method that directly links galaxies to their local environments, by using spatial cross-correlations of galaxy catalogues with maps from large-scale structure surveys (e.g., thermal Sunyaev-Zel'dovich [tSZ] effect, diffuse X-ray emission, weak lensing of galaxies or the CMB). We focus here on the quenching of galaxies and its link to local hot gas properties. Maps of galaxy overdensity and quenched fraction excess are constructed from volume-limited SDSS catalogs, which are cross-correlated with tSZ effect and X-ray maps from Planck and ROSAT, respectively. Strong signals out to Mpc scales are detected for most cross-correlations and are compared to predictions from the EAGLE and BAHAMAS cosmological hydrodynamical simulations. The simulations successfully reproduce many, but not all, of the observed power spectra, with an indication that environmental quenching may be too efficient in the simulations. We demonstrate that the cross-correlations are sensitive to both the internal (e.g., AGN and stellar feedback) and external processes (e.g., ram pressure stripping, harassment, strangulation, etc.) responsible for quenching. The methods outlined in this paper can be adapted to other observables and, with upcoming surveys, will provide a stringent test of physical models for environmental transformation.

In this work we consider the impact of spatially-uniform but time-varying dark energy (or `dynamical dark energy', DDE) on large-scale structure in a spatially flat universe, using large cosmological hydrodynamical simulations that form part of the BAHAMAS project. As DDE changes the expansion history of the universe, it impacts the growth of structure. We explore variations in DDE that are constrained to be consistent with the cosmic microwave background. We find that DDE can affect the clustering of matter and haloes at the ~10% level (suppressing it for so-called `freezing' models, while enhancing it for `thawing' models), which should be distinguishable with upcoming large-scale structure surveys. DDE cosmologies can also enhance or suppress the halo mass function (with respect to LCDM) over a wide range of halo masses. The internal properties of haloes are minimally affected by changes in DDE, however. Finally, we show that the impact of baryons and associated feedback processes is largely independent of the change in cosmology and that these processes can be modelled separately to typically better than a few percent accuracy

We introduce the ARTEMIS simulations, a new set of 42 zoomed-in, high-resolution (baryon particle mass of ~ 2x10^4 Msun/h), hydrodynamical simulations of galaxies residing in haloes of Milky Way mass, simulated with the EAGLE galaxy formation code with re-calibrated stellar feedback. In this study, we analyse the structure of stellar haloes, specifically the mass density, surface brightness, metallicity, colour and age radial profiles, finding generally very good agreement with recent observations of local galaxies. The stellar density profiles are well fitted by broken power laws, with inner slopes of ~ -3, outer slopes of ~ -4 and break radii that are typically ~ 20-40 kpc. The break radii generally mark the transition between in situ formation and accretion-driven formation of the halo. The metallicity, colour and age profiles show mild large-scale gradients, particularly when spherically-averaged or viewed along the major axes. Along the minor axes, however, the profiles are nearly flat, in agreement with observations. Overall, the structural properties can be understood by two factors: that in situ stars dominate the inner regions and that they reside in a spatially-flattened distribution that is aligned with the disc. Observations targeting both the major and minor axes of galaxies are thus required to obtain a complete picture of stellar haloes.

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

Female role models reduce the impact on women of stereotype threat, i.e., of being at risk of conforming to a negative stereotype about one's social, gender, or racial group. This can lead women scientists to underperform or to leave their scientific career because of negative stereotypes such as, not being as talented or as interested in science as men. Sadly, history rarely provides role models for women scientists; instead, it often renders these women invisible. In response to this situation, we present a selection of twelve outstanding women who helped to develop nuclear astrophysics.

We study the correlation between stellar mass and stellar metallicity in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) suite of cosmological hydrodynamical simulations. At a given stellar mass, simulated galaxies with lower stellar metallicities show, on average, higher gas fractions, higher specific star formation rates and younger stellar populations. Active galactic nuclei feedback seems to play an important role on the determination of the stellar metallicity at high stellar masses. In general, simulated systems follow a well-defined anticorrelation between stellar metallicity and gas fraction, which does not evolve significantly with redshift. All these trends are consistent with previous findings regarding the metallicity of the star-forming gas in EAGLE.

Recent studies have presented evidence for tension between the constraints on Omega_m and sigma_8 from the cosmic microwave background (CMB) and measurements of large-scale structure (LSS). This tension can potentially be resolved by appealing to extensions of the standard model of cosmology and/or untreated systematic errors in the modelling of LSS, of which baryonic physics has been frequently suggested. We revisit this tension using, for the first time, carefully-calibrated cosmological hydrodynamical simulations, which thus capture the back reaction of the baryons on the total matter distribution. We have extended the BAHAMAS simulations to include a treatment of massive neutrinos, which currently represents the best motivated extension to the standard model. We make synthetic thermal Sunyaev-Zel'dovich effect, weak galaxy lensing, and CMB lensing maps and compare to observed auto- and cross-power spectra from a wide range of recent observational surveys. We conclude that: i) in general there is tension between the primary CMB and LSS when adopting the standard model with minimal neutrino mass; ii) after calibrating feedback processes to match the gas fractions of clusters, the remaining uncertainties in the baryonic physics modelling are insufficient to reconcile this tension; and iii) if one accounts for internal tensions in the Planck CMB dataset (by allowing the lensing amplitude, A_Lens, to vary), invoking a non-minimal neutrino mass, typically of 0.2-0.4 eV, can resolve the tension. This solution is fully consistent with separate constraints from the primary CMB and baryon acoustic oscillations.

We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{\sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.

[Abridged] Typical disc galaxies forming in a LambdaCDM cosmology encounter a violent environment, where they often experience mergers with massive satellites. The fact that disc galaxies are ubiquitous in the local Universe suggests that a quiescent history is not necessary for their formation. Modern cosmological simulations can now obtain relatively realistic populations of disc galaxies, but it still remains to be clarified how discs manage to survive massive mergers. Here we use a suite of high-resolution hydrodynamical simulations set in a LambdaCDM cosmology to elucidate the fate of discs encountering massive mergers. We extract a sample of approximately 100 disc galaxies and follow the changes in their post-merger morphologies, as tracked by their disc-to-total ratios (D/T). We also examine the relations between their present-day morphology, assembly history and gas fractions. We find that approximately half of present-day disc galaxies underwent at least one merger with a satellite of total mass exceeding the host system's stellar mass, a third had mergers with satellites of mass exceeding 3 times the host's stellar mass, and approximately one-sixth had mergers with satellites of mass exceeding 10 times of the host's stellar mass. These mergers lead to a sharp, but often temporary, decrease in the D/T of the hosts, implying that discs are usually disrupted but then quickly re-grow. To do so, high cold gas fractions are required post-merger, as well as a relatively quiescent recent history (over a few Gyrs before z=0). Our results show that discs can form via diverse merger pathways and that quiescent histories are not the dominant mode of disc formation.

We quantify the correlations between gas-phase and stellar metallicities and global properties of galaxies, such as stellar mass, halo mass, age and gas fraction, in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) suite of cosmological hydrodynamical simulations. The slope of the correlation between stellar mass and metallicity of star-forming (SF) gas ($M_*-Z_{\rm SF,gas}$ relation) depends somewhat on resolution, with the higher-resolution run reproducing a steeper slope. This simulation predicts a non-zero metallicity evolution, increasing by $\approx 0.5$ dex at $\sim 10^9 {\rm M}_{\odot}$ since $z = 3$. The simulated relation between stellar mass, metallicity and star formation rate at $z \lesssim 5$ agrees remarkably well with the observed fundamental metallicity relation. At $M_* \lesssim 10^{10.3} {\rm M}_{\odot}$ and fixed stellar mass, higher metallicities are associated with lower specific star formation rates, lower gas fractions and older stellar populations. On the other hand, at higher $M_*$, there is a hint of an inversion of the dependence of metallicity on these parameters. The fundamental parameter that best correlates with the metal content, in the simulations, is the gas fraction. The simulated gas fraction-metallicity relation exhibits small scatter and does not evolve significantly since $z = 3$. In order to better understand the origin of these correlations, we analyse a set of lower resolution simulations in which feedback parameters are varied. We find that the slope of the simulated $M_*-Z_{\rm SF,gas}$ relation is mostly determined by stellar feedback at low stellar masses ($M_* \lesssim 10^{10} {\rm M}_{\odot}$), and at high masses ($M_* \gtrsim 10^{10} {\rm M}_{\odot}$) by the feedback from active galactic nuclei.

The Gaia astrometric mission may offer an unprecedented opportunity to discover new tidal streams in the Galactic halo. To test this, we apply nGC3, a great-circle-cell count method that combines position and proper motion data to identify streams, to ten mock Gaia catalogues of K giants and RR Lyrae stars constructed from cosmological simulations of Milky Way analogues. We analyse two sets of simulations, one using a combination of $N$-body and semi-analytical methods which has extremely high resolution, the other using hydro-dynamical methods, which captures the dynamics of baryons, including the formation of an in situ halo. These ten realisations of plausible Galactic merger histories allow us to assess the potential for the recovery of tidal streams in different Milky Way formation scenarios. We include the Gaia~selection function and observational errors in these mock catalogues. We find that the nGC3 method has a well-defined detection boundary in the space of stream width and projected overdensity, that can be predicted based on direct observables alone. We predict that about 4-13 dwarf galaxy streams can be detected in a typical Milky Way-mass halo with Gaia+nGC3, with an estimated efficiency of $>$80\% inside the detection boundary. The progenitors of these streams are in the mass range of the classical dwarf galaxies and may have been accreted as early as redshift $\sim3$. Finally, we analyse how different possible extensions of the Gaia mission will improve the detection of tidal streams.

We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single {\Lambda}CDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs.

The evolution of the metal content of galaxies and its relations to other global properties [such as total stellar mass (M*), circular velocity, star formation rate (SFR), halo mass, etc.] provides important constraints on models of galaxy formation. Here we examine the evolution of metallicity scaling relations of simulated galaxies in the Galaxies-Intergalactic Medium Interaction Calculation suite of cosmological simulations. We make comparisons to observations of the correlation of gas-phase abundances with M* (the mass-metallicity relation, MZR), as well as with both M* and SFR or gas mass fraction (the so-called 3D fundamental metallicity relations, FMRs). The simulated galaxies follow the observed local MZR and FMRs over an order of magnitude in M*, but overpredict the metallicity of massive galaxies (log M* > 10.5), plausibly due to inefficient feedback in this regime. We discuss the origin of the MZR and FMRs in the context of galactic outflows and gas accretion. We examine the evolution of mass-metallicity relations defined using different elements that probe the three enrichment channels (SNII, SNIa, and AGB stars). Relations based on elements produced mainly by SNII evolve weakly, whereas those based on elements produced preferentially in SNIa/AGB exhibit stronger evolution, due to the longer timescales associated with these channels. Finally, we compare the relations of central and satellite galaxies, finding systematically higher metallicities for satellites, as observed. We show this is due to the removal of the metal poor gas reservoir that normally surrounds galaxies and acts to dilute their gas-phase metallicity (via cooling/accretion onto the disk), but is lost due to ram pressure stripping for satellites.

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

It is shown that Starobinsky-like potentials can be realized in non-geometric flux compactifications of string theory, where the inflaton involves an axion whose shift symmetry can protect UV-corrections to the scalar potential. For that purpose we evaluate the backreacted, uplifted F-term axion-monodromy potential, which interpolates between a quadratic and a Starobinsky-like form. Limitations due to the requirements of having a controlled approximation of the UV theory and of realizing single-field inflation are discussed.

We combine the six high-resolution Aquarius dark matter simulations with a semi-analytic galaxy formation model to investigate the properties of the satellites of Milky Way-like galaxies. We find good correspondence with the observed luminosity function, luminosity-metallicity relation and radial distribution of the Milky Way satellites. The star formation histories of the dwarf galaxies in our model vary widely, in accordance with what is seen observationally. Ram-pressure stripping of hot gas from the satellites leaves a clear imprint of the environment on the characteristics of a dwarf galaxy. We find that the fraction of satellites dominated by old populations of stars matches observations well. However, the internal metallicity distributions of the model satellites appear to be narrower than observed. This may indicate limitations in our treatment of chemical enrichment, which is based on the instantaneous recycling approximation. Our model works best if the dark matter halo of the Milky Way has a mass of ~8 x 10^11 Msun, in agreement with the lower estimates from observations. The galaxy that resembles the Milky Way the most also has the best matching satellite luminosity function, although it does not contain an object as bright as the SMC or LMC. Compared to other semi-analytic models and abundance matching relations we find that central galaxies reside in less massive haloes, but the halo mass-stellar mass relation in our model is consistent both with hydrodynamical simulations and with recent observations.

In the local Universe, galaxies in groups and clusters contain less gas and are less likely to be forming stars than their field counterparts. This effect is not limited to the central group/cluster regions, but is shown by recent observations to persist out to several virial radii. To gain insight into the extent and cause of this large-scale environmental influence, we use a suite of high-resolution cosmological hydrodynamic simulations to analyse galaxies around simulated groups and clusters of a wide range of mass (log M/M_sun = [13.0, 15.2]). In qualitative agreement with the observations, we find a systematic depletion of both hot and cold gas and a decline in the star forming fraction of galaxies as far out as ~ 5 r200 from the host centre. While a substantial fraction of these galaxies are on highly elliptical orbits and are not infalling for the first time (~ 50 per cent at 2 r200, independent of host mass) or are affected by `pre-processing' (less than 10 per cent of galaxies around groups, increasing to ~ 50 per cent around a massive cluster), even a combination of these indirect mechanisms does not fully account for the environmental influence, particularly in the case of the hot gas content. Direct ram pressure interaction with an extended gas `halo' surrounding groups and clusters is shown to be sufficiently strong to strip the hot gas atmospheres of infalling galaxies out to ~ 5 r200. We show that this influence is highly anisotropic, with ram pressure along filaments enhanced by up to a factor of 100 despite significant co-flow of gas and galaxies.

We examine the rotation rates, sizes, and star formation (SF) efficiencies of a representative population of simulated disc galaxies extracted from the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamic simulations. These simulations include efficient, but energetically feasible supernova feedback, but have not been tuned in any way to produce 'realistic' disc galaxies. Yet, they generate a large number of discs, without requiring extremely high resolution. Over the wide galaxy stellar mass range, 9.0 < log10[Mstar (Msun)] < 10.5, the simulations reproduce the observed Tully-Fisher relation, the rotation curves of disc galaxies in bins of stellar mass, the mass-size relation of disc galaxies, the optical rotation to virial circular velocity ratio (Vopt/Vvir), and the SF efficiencies of disc galaxies as inferred from stacked weak lensing and stacked satellite kinematics observations. They also reproduce the specific star formation rates of ~L* galaxies but predict too low levels of star formation for low-mass galaxies, which is plausibly due to the finite resolution of the simulations. At higher stellar masses, log10[Mstar (\Msun)] > 10.6, the simulated galaxies are too concentrated and have too high SF efficiencies. We conjecture that this shortcoming reflects the neglect of feedback from accreting supermassive black holes in these simulations. We conclude that it is possible to generate a representative population of disc galaxies that reproduces many of the observed trends of local disc galaxies using standard numerical hydrodynamic techniques and a plausible implementation of the "subgrid" astrophysical processes thought to be relevant to galaxy formation.

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the global structure and kinematics of stellar spheroids of Milky Way mass disc galaxies. Font et al. have recently demonstrated that these simulations are able to successfully reproduce the satellite luminosity functions and the metallicity and surface brightness profiles of the spheroids of the Milky Way and M31. A key to the success of the simulations is a significant contribution to the spheroid from stars that formed in situ. While the outer halo is dominated by accreted stars, stars formed in the main progenitor of the galaxy dominate at r < ~30 kpc. In the present study we show that this component was primarily formed in a proto-disc at high redshift and was subsequently liberated from the disc by dynamical heating associated with mass accretion. As a consequence of its origin, the in situ component of the spheroid has different kinematics (namely net prograde rotation with respect to the disc) than that of the spheroid component built from the disruption of satellites. In addition, the in situ component has a flattened distribution, that is due in part to its rotation. We make comparisons with measurements of the shape and kinematics of local galaxies, including the Milky Way and M31, and stacked observations of more distant galaxies. We find that the simulated disc galaxies have spheroids of the correct shape (oblate with a median axis ratio of ~0.6 at radii of < ~30 kpc, but note there is significant system-to-system scatter in this quantity) and that the kinematics show evidence for two components (due to in situ vs. accreted), as observed. Our findings therefore add considerable weight to the importance of dissipative processes in the formation of stellar haloes and to the notion of a 'dual stellar halo'.

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the formation of stellar spheroids of Milky Way-mass disc galaxies. The simulations contain accurate treatments of metal-dependent radiative cooling, star formation, supernova feedback, and chemodynamics, and the large volumes that have been simulated yield an unprecedentedly large sample of ~400 simulated L_* disc galaxies. The simulated galaxies are surrounded by low-mass, low-surface brightness stellar haloes that extend out to ~100 kpc and beyond. The diffuse stellar distributions bear a remarkable resemblance to those observed around the Milky Way, M31 and other nearby galaxies, in terms of mass density, surface brightness, and metallicity profiles. We show that in situ star formation typically dominates the stellar spheroids by mass at radii of r < 30 kpc, whereas accretion of stars dominates at larger radii and this change in origin induces a change in slope of the surface brightness and metallicity profiles, which is also present in the observational data. The system-to-system scatter in the in situ mass fractions of the spheroid, however, is large and spans over a factor of 4. Consequently, there is a large degree of scatter in the shape and normalisation of the spheroid density profile within r < 30 kpc (e.g., when fit by a spherical powerlaw profile the indices range from -2.6 to -3.4). We show that the in situ mass fraction of the spheroid is linked to the formation epoch of the system. Dynamically older systems have, on average, larger contributions from in situ star formation, although there is significant system-to-system scatter in this relationship. Thus, in situ star formation likely represents the solution to the longstanding failure of pure accretion-based models to reproduce the observed properties of the inner spheroid.

(Abridged) We study the phase-space distribution of satellite galaxies associated with late-type galaxies in the GIMIC suite of simulations. GIMIC consists of re-simulations of 5 cosmologically representative regions from the Millennium simulation, which have higher resolution and incorporate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions (r ~ 0.1r_200) of the dark matter halo, both in shape and angular momentum, there can be substantial misalignments at larger radii (r ~ r_200). Misalignments of > 45 deg are seen in ~ 30% of our sample. We find that the satellite population aligns with the shape (and angular momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak owing to the mismatch between the disc and dark matter halo. Roughly 20% of the satellite systems with ten bright galaxies within r_200 exhibit a polar spatial alignment with respect to the galaxy --- an orientation reminiscent of the classical satellites of the Milky Way. We find that a small fraction (~10%) of satellite systems show evidence for rotational support which we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). This preference towards co-rotation is stronger in the inner regions of the halo where the most massive satellites accreted at relatively early times are located. We attribute the anisotropic spatial distribution and angular momentum bias of the satellites at z=0 to their directional accretion along the major axes of the dark matter halo. The satellite galaxies have been accreted relatively recently compared to the dark matter mass and have experienced less phase-mixing and relaxation --- the memory of their accretion history can remain intact to z=0.

Extensive photometric and spectroscopic surveys of the Andromeda galaxy (M31) have discovered tidal debris features throughout M31's stellar halo. We present stellar kinematics and metallicities in fields with identified substructure from our on-going SPLASH survey of M31 red giant branch stars with the DEIMOS spectrograph on the Keck II 10-m telescope. Radial velocity criteria are used to isolate members of the kinematically-cold substructures. The substructures are shown to be metal-rich relative to the rest of the dynamically hot stellar population in the fields in which they are found. We calculate the mean metallicity and average surface brightness of the various kinematical components in each field, and show that, on average, higher surface brightness features tend to be more metal-rich than lower surface brightness features. Simulations of stellar halo formation via accretion in a cosmological context are used to illustrate that the observed trend can be explained as a natural consequence of the observed dwarf galaxy mass-metallicity relation. A significant spread in metallicity at a given surface brightness is seen in the data; we show that this is due to time effects, namely the variation in the time since accretion of the tidal streams' progenitor onto the host halo. We show that in this theoretical framework a relationship between the alpha-enhancement and surface brightness of tidal streams is expected, which arises from the varying times of accretion of the progenitor satellites onto the host halo. Thus, measurements of the alpha-enrichment, metallicity, and surface brightness of tidal debris can be used to reconstruct the luminosity and time of accretion onto the host halo of the progenitors of tidal streams.

We use the galaxy stellar mass and halo merger tree information from the semi-analytic model galaxy catalogue of Font et al. (2009) to examine the accretion of galaxies into a large sample of groups and clusters, covering a wide range in halo mass (10E12.9 to 10E15.3 Msun/h), and selected from each of four redshift epochs (z=0, 0.5, 1.0 and 1.5). We find that clusters at all examined redshifts have accreted a significant fraction of their final galaxy populations through galaxy groups. A 10E14.5 Msun/h mass cluster at z=0 has, on average, accreted ~ 40% of its galaxies (Mstellar > 10E9 Msun/h) from halos with masses greater than 10E13 Msun/h. Further, the galaxies which are accreted through groups are more massive, on average, than galaxies accreted through smaller halos or from the field population. We find that at a given epoch, the fraction of galaxies accreted from isolated environments is independent of the final cluster or group mass. In contrast, we find that observing a cluster of the same halo mass at each redshift epoch implies different accretion rates of isolated galaxies, from 5-6 % per Gyr at z=0 to 15% per Gyr at z=1.5. We find that combining the existence of a Butcher Oemler effect at z=0.5 and the observations that galaxies within groups display significant environmental effects with galaxy accretion histories justifies striking conclusions. Namely, that the dominant environmental process must begin to occur in halos of 10E12 -- 10E13 Msun/h, and act over timescales of > 2 Gyrs. This argues in favor of a mechanism like "strangulation", in which the hot halo of a galaxy is stripped upon infalling into a more massive halo . This simple model predicts that by z=1.5 galaxy groups and clusters will display little to no environmental effects.

Current models of galaxy formation predict satellite galaxies in groups and clusters that are redder than observed. We investigate the effect on the colours of satellite galaxies produced by the ram pressure stripping of their hot gaseous atmospheres as the satellites orbit within their parent halo. We incorporate a model of the stripping process based on detailed hydrodynamic simulations within the Durham semi-analytic model of galaxy formation. The simulations show that the environment in groups and clusters is less aggressive than previously assumed. The main uncertainty in the model is the treatment of gas expelled by supernovae. With reasonable assumptions for the stripping of this material, we find that satellite galaxies are able to retain a significant fraction of their hot gas for several Gigayears, thereby replenishing their reservoirs of cold, star forming gas and remaining blue for a relatively long period of time. A bimodal distribution of galaxy colours, similar to that observed in SDSS data, is established and the colours of the satellite galaxies are in good agreement with the data. In addition, our model naturally accounts for the observed dependence of satellite colours on environment, from small groups to high mass clusters.

This paper explores the mapping between the observable properties of a stellar halo in phase- and abundance-space and the parent galaxy's accretion history in terms of the characteristic epoch of accretion and mass and orbits of progenitor objects. The study utilizes a suite of eleven stellar halo models constructed within the context of a standard LCDM cosmology. The results demonstrate that coordinate-space studies are sensitive to the recent (0-8 Gyears ago) merger histories of galaxies (this timescale corresponds to the last few to tens of percent of mass accretion for a Milky-Way-type galaxy). Specifically, the {\it frequency, sky coverage} and {\it fraction of stars} in substructures in the stellar halo as a function of surface brightness are indicators of the importance of recent merging and of the luminosity function of infalling dwarfs. The {\it morphology} of features serves as a guide to the orbital distribution of those dwarfs. Constraints on the earlier merger history (> 8 Gyears ago) can be gleaned from the abundance patterns in halo stars: within our models, dramatic differences in the dominant epoch of accretion or luminosity function of progenitor objects leave clear signatures in the [alpha/Fe] and [Fe/H] distributions of the stellar halo - halos dominated by very early accretion have higher average [alpha/Fe], while those dominated by high luminosity satellites have higher [Fe/H]. This intuition can be applied to reconstruct much about the merger histories of nearby galaxies from current and future data sets.

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

Recent observations have revealed that high surface-brightness, metal-rich debris is present over large regions of the Andromeda (M31) stellar halo. We use a set of numerical models to determine whether extended metal-rich debris is expected to exist in galaxy halos formed in a hierarchical LambdaCDM universe. We identify tidal debris in the simulations according to the current best surface brightness detection limits in M31 and demonstrate that bright features in the outer halo regions usually originate from single satellites, with very little contamination from other sources due to the low probability of tidal streams from two overlapping accretion events. In contrast, high-surface brightness features in the inner halo often originate from multiple progenitors. We also compare the age and metallicity distribution of the debris with the well-mixed stellar halos in which they reside. We find that high-surface brightness tidal debris is produced almost exclusively by relatively high mass progenitors (Mstar ~ 10^7- 10^9 Msun) and thus is expected to be of moderate- to high-metallicity. Similarly, in our models the smooth inner halo is expected to be metal-rich as this region has been built up mainly from massive satellites. Our results imply that the stellar populations of substructure observed around external galaxies with current techniques should not resemble old and metal-poor dwarf spheroidal satellites, nor the underlying component of the stellar halo.

[Abridged] Motivated by upcoming data from astrometric and spectroscopic surveys of the Galaxy, we explore the chemical abundance properties and phase-space distributions in hierarchically-formed stellar halo simulations set in a LambdaCDM Universe. Our sample of Milky-Way type stellar halo simulations result in average metallicities that range from [Fe/H] = -1.3 to -0.9, with the most metal poor halos resulting from accretion histories that lack destructive mergers with massive (metal rich) satellites. Our stellar halo metallicities increase with stellar halo mass. The slope of the [Fe/H]-stellar mass trend mimics that of the satellite galaxies that were destroyed to build the halos, implying that the relation propagates hierarchically. All simulated halos contain a significant fraction of old stellar populations accreted more than 10 Gyr ago and in a few cases, some intermediate age populations exist. In contrast with the Milky Way, many of our simulated stellar halos contain old stellar populations which are metal rich, originating in the early accretion of massive satellites. We suggest that the (metal rich) stellar halo of M31 falls into this category, while the more metal poor halo of the Milky Way is lacking in early massive accretion events. Interestingly, our hierarchically-formed stellar halos often have non-negligible metallicity gradients in both [Fe/H] and [alpha/Fe]. These gradients extend a few tens of kpc, and can be as large as 0.5 dex in [Fe/H] and 0.2 dex in [alpha/Fe]. Finally, we find that chemical abundances can act as a rough substitute for time of accretion of satellite galaxies. We propose a criterion for identifying tidal streams spatially by selecting stars with [alpha/Fe] ratios below solar.

(Abridged) The hierarchical formation scenario for the stellar halo requires the accretion and disruption of dwarf galaxies, yet low-metallicity halo stars are enriched in alpha-elements compared to similar, low-metallicity stars in dwarf spheroidal (dSph) galaxies. We address this primary challenge for the hierarchical formation scenario for the stellar halo by combining chemical evolution modelling with cosmologically-motivated mass accretion histories for the Milky Way dark halo and its satellites. We demonstrate that stellar halo and dwarf galaxy abundance patterns can be explained naturally within the LCDM framework. Our solution relies fundamentally on the LCDM model prediction that the majority of the stars in the stellar halo were formed within a few relatively massive, ~5 x 10^10 Msun, dwarf irregular (dIrr)-size dark matter halos, which were accreted and destroyed ~10 Gyr in the past. These systems necessarily have short-lived, rapid star formation histories, are enriched primarily by Type II supernovae, and host stars with enhanced [a/Fe] abundances. In contrast, dwarf spheroidal galaxies exist within low-mass dark matter hosts of ~10^9 Msun, where supernovae winds are important in setting the intermediate [a/Fe] ratios observed. Our model includes enrichment from Type Ia and Type II supernovae as well as stellar winds, and includes a physically-motivated supernovae feedback prescription calibrated to reproduce the local dwarf galaxy stellar mass - metallicity relation. We use representative examples of the type of dark matter halos we expect to host a destroyed ``stellar halo progenitor'' dwarf, a surviving dIrr, and a surviving dSph galaxy, and show that their derived abundance patterns, stellar masses, and gas masses are consistent with those observed for each type of system.

We present a cosmologically motivated model for the hierachical formation of the stellar halo which includes a semi-analytic treatment of galactic chemical enrichment coupled to numerical simulations that track the orbital evolution and tidal disruption of satellites. A major motivating factor in this investigation is the observed systematic difference between the chemical abundances of stars in satellite galaxies and those in the Milky Way halo. Specifically, for the same [Fe/H] values, stars in neighboring satellite galaxies display significantly lower [alpha/Fe] ratios than stars in the halo. We find that the observed chemical abundance patterns are a natural outcome of the process of hierarchical assembly of the Galaxy. This result follows because the stellar halo in this context is built up from satellite galaxies accreted early on (more than 8-9 Gyr ago) and enriched in alpha-elements produced in Type II supernovae. In contrast, satellites which survive today are typically accreted late (within the last 4-5 Gyr) with nearly solar [alpha/Fe] values as a result of contributions from both Type II and Type Ia supernovae. We use our model to investigate the abundance distribution functions (using both [Fe/H] and [alpha/Fe] ratios) for stars in the halo and in surviving satellites. Our results suggest that the shapes and peaks in the abundance distribution functions provide a direct probe of the accretion histories of galaxies.

We examine the nature of the progenitor of the giant stellar stream in M31 using as constraints new radial velocity measurements of stream red giant stars (presented in the companion paper by Guhathakurta et al. 2005, astro-ph/0406145) along with other M31 data sets available in the literature. We find that the observations are best fit by orbits that are highly eccentric and close to edge-on, with apo- to peri-center ratios of order 25 - 30, and with apocenters at or only slightly beyond the southern edge of the current data. Among these orbits, we are able to find a few that plausibly connect the stream with the northern spur or with the low-surface-brightness feature of similar high metallicity as the stream (originally reported by Ferguson et al. 2002) to the east of M31's center. In the latter case, if the connection is real, then the eastern debris should lie well in front of M31 near the apocenter of the orbit. Both the width of the debris and velocity dispersion measurements imply a rough lower limit on the mass of the progenitor of 10^8 M_sun. We use this limit and our orbits to discuss which of M31's satellites could be plausibly associated with the stream. In addition, we predict that the width of the stream should increase beyond the southern edge of the current data around the apocenter of the orbit and that the line-of-sight velocity dispersion should exhibit significant variations along the stream.

Photoevaporation may provide an explanation for the short lifetimes of disks around young stars. With the exception of neutral oxygen lines, the observed low-velocity forbidden line emission from T Tauri stars can be reproduced by photoevaporating models. The natural formation of a gap in the disk at several AU due to photoevaporation and viscous spreading provides a possible halting mechanism for migrating planets and an explanation for the abundance of observed planets at these radii.

We examine the ability of photoevaporative disk winds to explain the low-velocity components observed in the forbidden line spectra of low-mass T Tauri stars. Using the analytic model of Shu, Johnstone, & Hollenbach (1993) and Hollenbach et al. (1994) as a basis, we examine the characteristics of photoevaporative outflows with hydrodynamic simulations. General results from the simulations agree well with the analytic predictions, although some small differences are present. Most importantly, the flow of material from the disk surface develops at smaller radii than in the analytic approximations and the flow-velocity from the disk surface is only one-third the sound speed. A detailed presentation of observational consequences of the model is given, including predicted line widths, blue-shifts, and integrated luminosities of observable sulfur and nitrogen emission lines. We demonstrate that these predictions are in agreement with current observational data on the low-velocity forbidden line emission of ionized species from T Tauri stars. This is in contrast with magnetic wind models, which systematically under-predict these forbidden line luminosities. However, the present model cannot easily account for the luminosities of neutral oxygen lines in T Tauri stars.

We present an analytical description of the motion in the singular logarithmic potential. This potential plays an important role in the modeling of triaxial systems (like elliptical galaxies) or bars in the centers of galaxy disks. In order to obtain information about the motion near the singularity, we resort to McGehee -type transformations and regularize the vector field. In the axis-symmetric case (b=1), we offer a complete description of the global dynamics. In the non axis-symmetric case (b<1), we prove that all orbits, with the exception of a negligible set, are centrophobic and retrieve numerically partial aspects of the orbital structure.

We examine recent suggestions that substructure in cold dark matter (CDM) halos may be in conflict with the presence of thin, dynamically fragile stellar disks. N-body simulations of an isolated disk/bulge/halo model of the Milky Way that includes several hundred dark matter satellites with masses, densities and orbits derived from high-resolution cosmological CDM simulations indicate that substructure at $z=0$ plays only a minor dynamical role in the heating of the disk over several Gyrs. This is because the orbits of satellites in present-day CDM halos seldom take them near the disk, where their tidal effects are greatest. Unless the effects of substructure are very different at earlier times, our models suggest that substructure might not preclude virialized CDM halos from being acceptable hosts of thin stellar disks like that of the Milky Way.

New molecular (13CO J=3-2) and dust continuum (450 and 850 micron) SCUBA maps of the NGC7129 star forming region are presented, complemented by C18O J=3-2 spectra at several positions within the mapped region. The maps include the Herbig Ae/Be star LkHalpha 234, the far-infrared source NGC 7129 FIRS2 and several other pre-stellar sources embedded within the molecular ridge. The SCUBA maps help us understand the nature of the pre-main sequence stars in this actively star forming region. A deeply embedded submillimeter source, SMM2, not clearly seen in any earlier data set, is shown to be a pre-stellar core or possibly a protostar. The highest continuum peak emission is identified with the deeply embedded source IRS6, a few arcseconds away from LkHalpha 234, and also responsible for both the optical jet and the molecular outflow. The gas and dust masses are found to be consistent, suggesting little or no CO depletion onto grains. The dust emissivity index is lower towards the dense compact sources, beta ~1 - 1.6, and higher, beta ~ 2.0, in the surrounding cloud, implying small size grains in the PDR ridge, whose mantles have been evaporated by the intense UV radiation.