Search

Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmology Telescope DR5 kinematic Sunyaev-Zel'dovich (kSZ) to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model, `baryonification'. First, using WL only, we compare the $S_8$ constraints using baryonification to a simulation-calibrated halo model, a simulation-based emulator model and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of $S_8$ and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters and shifts $S_8$ to $S_8=0.823^{+0.019}_{-0.020}$, a higher value than attained using the WL-only analysis. We measure the suppression of the non-linear matter power spectrum using WL + kSZ and constrain a mean feedback scenario that is more extreme than the predictions from most hydrodynamical simulations. We constrain the baryon fractions and the gas mass fractions and find them to be generally lower than inferred from X-ray observations and simulation predictions. We conclude that the WL + kSZ measurements provide a new and complementary benchmark for building a coherent picture of the impact of gas around galaxies across observations.

TThe fundamental properties of the postulated dark matter (DM) affect the internal structure of gravitationally-bound structures. In the cold dark matter paradigm, DM particles interact only via gravity. Their distribution is well represented by an Einasto profile with shape parameter $\alpha\approx0.18$, in the smallest dwarf galaxies or the most massive galaxy clusters alike. Conversely, if dark matter particles self-interact via additional forces, we expect the mass density profiles of DM halos to flatten in their central regions, thereby increasing the Einasto shape parameter. We measure the structural properties of the 12 massive X-COP galaxy clusters from observations of their hot gaseous atmosphere using the X-ray observatory XMM-Newton, and of the Sunyaev-Zeldovich effect using the Planck all-sky survey. After removing morphologically disturbed systems, we measure Einasto shape parameters with mean $\langle\alpha\rangle = 0.19 \pm 0.03$ and intrinsic scatter $\sigma_{\alpha}=0.06$, in close agreement with the prediction of the cold dark matter paradigm. We use cosmological hydrodynamical simulations of cluster formation with self-interacting DM (BAHAMAS-SIDM) to determine how the Einasto shape parameter depends on the self-interaction cross section. We use the fitted relation to turn our measurements of $\alpha$ into constraints on the self-interaction cross section, which imply $\sigma/m < 0.19$ cm$^2$/g (95% confidence level) at collision velocity $v_\mathrm{DM-DM}\sim1,000$ km/s. This is lower than the interaction cross-section required for dark matter self-interactions to solve the core-cusp problem in dwarf spheroidal galaxies, unless the cross section is a strong function of velocity., Comment: Re-submitted to A&A after revision

We evaluate the ability of Convolutional Neural Networks (CNNs) to predict galaxy cluster masses in the BAHAMAS hydrodynamical simulations. We train four separate single-channel networks using: stellar mass, soft X-ray flux, bolometric X-ray flux, and the Compton $y$ parameter as observational tracers, respectively. Our training set consists of $\sim$4800 synthetic cluster images generated from the simulation, while an additional $\sim$3200 images form a validation set and a test set, each with 1600 images. In order to mimic real observation, these images also contain uncorrelated structures located within 50 Mpc in front and behind clusters and seen in projection, as well as instrumental systematics including noise and smoothing. In addition to CNNs for all the four observables, we also train a `multi-channel' CNN by combining the four observational tracers. The learning curves of all the five CNNs converge within 1000 epochs. The resulting predictions are especially precise for halo masses in the range $10^{13.25}M_{\odot}

On the scale of galactic haloes, the distribution of matter in the cosmos is affected by energetic, non-gravitational processes; so-called baryonic feedback. A lack of knowledge about the details of how feedback processes redistribute matter is a source of uncertainty for weak-lensing surveys, which accurately probe the clustering of matter in the Universe over a wide range of scales. We develop a cosmology-dependent model for the matter distribution that simultaneously accounts for the clustering of dark matter, gas and stars. We inform our model by comparing it to power spectra measured from the BAHAMAS suite of hydrodynamical simulations. As well as considering matter power spectra, we also consider spectra involving the electron-pressure field, which directly relates to the thermal Sunyaev-Zel'dovich (tSZ) effect. We fit parameters in our model so that it can simultaneously model both matter and pressure data and such that the distribution of gas as inferred from tSZ has influence on the matter spectrum predicted by our model. We present two variants; one that matches the feedback-induced suppression seen in the matter-matter power spectrum at the per-cent level and a second that matches the matter-matter data slightly less well (~2 per cent), but that is able to simultaneously model the matter-electron pressure spectrum at the ~15 per-cent level. We envisage our models being used to simultaneously learn about cosmological parameters and the strength of baryonic feedback using a combination of tSZ and lensing auto- and cross-correlation data., Comment: Main paper: 21 pages, 11 figures; Appendix: 16 pages, 13 figures; matches accepted version - referee report addressed, clarifications and extra text compared to v1, improved model, no significant changes of results

The location of a galaxy cluster's centroid is typically derived from observations of the galactic and/or gas component of the cluster, but these typically deviate from the true centre. This can produce bias when observations are combined to study average cluster properties. Using data from the BAHAMAS cosmological hydrodynamic simulations we study this bias in both two and three dimensions for 2000 clusters over the $10^{13} - 10^{15} ~\mathrm{M_{\odot}}$ mass range. We quantify and model the offset distributions between observationally-motivated centres and the `true' centre of the cluster, which is taken to be the most gravitationally bound particle measured in the simulation. We fit the cumulative distribution function of offsets with an exponential distribution and a Gamma distribution fit well with most of the centroid definitions. The galaxy-based centres can be seen to be divided into a mis-centred group and a well-centred group, with the well-centred group making up about $60\%$ of all the clusters. Gas-based centres are overall less scattered than galaxy-based centres. We also find a cluster-mass dependence of the offset distribution of gas-based centres, with generally larger offsets for smaller mass clusters. We then measure cluster density profiles centred at each choice of the centres and fit them with empirical models. Stacked, mis-centred density profiles fit to the Navarro-Frenk-White dark-matter profile and Komatsu-Seljak gas profile show that recovered shape and size parameters can significantly deviate from the true values. For the galaxy-based centres, this can lead to cluster masses being underestimated by up to $10\%$., Comment: 12 pages, 9 Figures

Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxy's cold gas content and subsequent star formation rate. Conversely, the thermodynamical properties of the ICM in this transition region should also feel the influence of accreting substructure (i.e. galaxies and groups), whose passage can drive shocks. In this paper, we use a suite of cosmological hydrodynamical zoom simulations of a single galaxy cluster, drawn from the nIFTy comparison project, to study how the dynamics of substructure accreted from the cosmic web influences the thermodynamical properties of the ICM in the cluster's outskirts. We demonstrate how features evident in radial profiles of the ICM (e.g. gas density and temperature) can be linked to strong shocks, transient and short-lived in nature, driven by the passage of substructure. The range of astrophysical codes and galaxy formation models in our comparison are broadly consistent in their predictions (e.g. agreeing when and where shocks occur, but differing in how strong shocks will be); this is as we would expect of a process driven by large-scale gravitational dynamics and strong, inefficiently radiating, shocks. This suggests that mapping such shock structures in the ICM in a cluster's outskirts (via e.g. radio synchrotron emission) could provide a complementary measure of its recent merger and accretion history., Comment: 15 pages, 13 figures, version accepted for publication in MNRAS

Weak gravitational lensing depends on the integrated mass along the line of sight. Baryons contribute to the mass distribution of galaxy clusters and the resulting mass estimates from lensing analysis. We use the cosmo-OWLS suite of hydrodynamic simulations to investigate the impact of baryonic processes on the bias and scatter of weak lensing mass estimates of clusters. These estimates are obtained by fitting NFW profiles to mock data using MCMC techniques. In particular, we examine the difference in estimates between dark matter-only runs and those including various prescriptions for baryonic physics. We find no significant difference in the mass bias when baryonic physics is included, though the overall mass estimates are suppressed when feedback from AGN is included. For lowest-mass systems for which a reliable mass can be obtained ($M_{200} \approx 2 \times 10^{14}$ $M_{\odot}$), we find a bias of $\approx -10$ per cent. The magnitude of the bias tends to decrease for higher mass clusters, consistent with no bias for the most massive clusters which have masses comparable to those found in the CLASH and HFF samples. For the lowest mass clusters, the mass bias is particularly sensitive to the fit radii and the limits placed on the concentration prior, rendering reliable mass estimates difficult. The scatter in mass estimates between the dark matter-only and the various baryonic runs is less than between different projections of individual clusters, highlighting the importance of triaxiality., Comment: 11 pages, 6 figures. Accepted by Monthly Notices of the Royal Astronomical Society

Modern cosmological simulations rely heavily on feedback from active galactic nuclei (AGN) in order to stave off overcooling in massive galaxies and galaxy groups and clusters. An important independent test is whether or not the simulations capture the broad demographics of the observed AGN population. Here, we have used the cosmo-OWLS suite of cosmological hydrodynamical simulations to produce realistic synthetic catalogs of X-ray AGN out to $z$=3, with the aim of comparing the catalogs to the observed X-ray AGN population in the XXL survey and other recent surveys. We focused on the unabsorbed X-ray luminosity function (XLF), the Eddington ratio distribution, the black hole mass function, and the projected clustering of X-ray AGN. To compute the unabsorbed XLF of the simulated AGN, we used recent empirically-determined bolometric corrections. We show that the simulated AGN sample accurately reproduces the observed XLF over 3 orders of magnitude in X-ray luminosity in all redshift bins. To compare to the observed Eddington ratio distribution and the clustering of AGN, we produced detailed 'XMM-Newton-detected' catalogs of the simulated AGN. This requires the production of synthetic X-ray images extracted from light cones of the simulations that fold in the relevant instrumental effects of XMM-Newton. We apply a luminosity- and redshift-dependent obscuration function for the AGN and employ the same AGN detection algorithm as used for the real XXL survey. We demonstrate that the detected population of simulated AGN reproduces the observed Eddington ratio distribution and projected clustering from XXL quite well. We conclude that the simulations have a broadly realistic population of AGN and that our synthetic X-ray AGN catalogs should be useful for interpreting additional trends and as a helpful tool for quantifying AGN contamination in galaxy group and cluster X-ray surveys., Comment: 14 pages, accepted for publication in A&A

The evolution of galaxy cluster counts is a powerful probe of several fundamental cosmological parameters. A number of recent studies using this probe have claimed tension with the cosmology preferred by the analysis of the Planck primary CMB data, in the sense that there are fewer clusters observed than predicted based on the primary CMB cosmology. One possible resolution to this problem is systematic errors in the absolute halo mass calibration in cluster studies, which is required to convert the standard theoretical prediction (the halo mass function) into counts as a function of the observable (e.g., X-ray luminosity, Sunyaev-Zel'dovich flux, optical richness). Here we propose an alternative strategy, which is to directly compare predicted and observed cluster counts as a function of the one-dimensional velocity dispersion of the cluster galaxies. We argue that the velocity dispersion of groups/clusters can be theoretically predicted as robustly as mass but, unlike mass, it can also be directly observed, thus circumventing the main systematic bias in traditional cluster counts studies. With the aid of the BAHAMAS suite of cosmological hydrodynamical simulations, we demonstrate the potential of the velocity dispersion counts for discriminating even similar $\Lambda$CDM models. These predictions can be compared with the results from existing redshift surveys such as the highly-complete Galaxy And Mass Assembly (GAMA) survey, and upcoming wide-field spectroscopic surveys such as the Wide Area Vista Extragalactic Survey (WAVES) and the Dark Energy Survey Instrument (DESI)., Comment: 15 pages, 13 figures. Accepted for publication in MNRAS. New section on cosmological forecasts added

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey., Comment: 17 pages, accepted in A&A

Traditionally, galaxy clusters have been expected to retain all the material accreted since their formation epoch. For this reason, their matter content should be representative of the Universe as a whole, and thus their baryon fraction should be close to the Universal baryon fraction. We make use of the sample of the 100 brightest galaxy clusters discovered in the XXL Survey to investigate the fraction of baryons in the form of hot gas and stars in the cluster population. We measure the gas masses of the detected halos and use a mass--temperature relation directly calibrated using weak-lensing measurements for a subset of XXL clusters to estimate the halo mass. We find that the weak-lensing calibrated gas fraction of XXL-100-GC clusters is substantially lower than was found in previous studies using hydrostatic masses. Our best-fit relation between gas fraction and mass reads $f_{\rm gas,500}=0.055_{-0.006}^{+0.007}\left(M_{\rm 500}/10^{14}M_\odot\right)^{0.21_{-0.10}^{+0.11}}$. The baryon budget of galaxy clusters therefore falls short of the Universal baryon fraction by about a factor of two at $r_{\rm 500}$. Our measurements require a hydrostatic bias $1-b=M_X/M_{\rm WL}=0.72_{-0.07}^{+0.08}$ to match the gas fraction obtained using lensing and hydrostatic equilibrium. Comparing our gas fraction measurements with the expectations from numerical simulations, our results favour an extreme feedback scheme in which a significant fraction of the baryons are expelled from the cores of halos. This model is, however, in contrast with the thermodynamical properties of observed halos, which might suggest that weak-lensing masses are overestimated. We note that a mass bias $1-b=0.58$ as required to reconcile Planck CMB and cluster counts should translate into an even lower baryon fraction, which poses a major challenge to our current understanding of galaxy clusters. [Abridged], Comment: 13th paper in the XXL series, A&A in press

We test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at $0.15

We report results for the alignments of galaxies in the EAGLE and cosmo-OWLS simulations as a function of galaxy separation and halo mass. The combination of these hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($10.7

We investigate the evolution of galaxy masses and star formation rates in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These comprise a suite of hydrodynamical simulations in a $\Lambda$CDM cosmogony with subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. The subgrid feedback was calibrated to reproduce the observed present-day galaxy stellar mass function and galaxy sizes. Here we demonstrate that the simulations reproduce the observed growth of the stellar mass density to within 20 per cent. The simulation also tracks the observed evolution of the galaxy stellar mass function out to redshift z = 7, with differences comparable to the plausible uncertainties in the interpretation of the data. Just as with observed galaxies, the specific star formation rates of simulated galaxies are bimodal, with distinct star forming and passive sequences. The specific star formation rates of star forming galaxies are typically 0.2 to 0.4 dex lower than observed, but the evolution of the rates track the observations closely. The unprecedented level of agreement between simulation and data makes EAGLE a powerful resource to understand the physical processes that govern galaxy formation., Comment: Accepted for publication by MNRAS; 18 pages, 7 figures; comments welcome

We investigate the internal structure and density profiles of halos of mass $10^{10}-10^{14}~M_\odot$ in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These follow the formation of galaxies in a $\Lambda$CDM Universe and include a treatment of the baryon physics thought to be relevant. The EAGLE simulations reproduce the observed present-day galaxy stellar mass function, as well as many other properties of the galaxy population as a function of time. We find significant differences between the masses of halos in the EAGLE simulations and in simulations that follow only the dark matter component. Nevertheless, halos are well described by the Navarro-Frenk-White (NFW) density profile at radii larger than ~5% of the virial radius but, closer to the centre, the presence of stars can produce cuspier profiles. Central enhancements in the total mass profile are most important in halos of mass $10^{12}-10^{13}M_\odot$, where the stellar fraction peaks. Over the radial range where they are well resolved, the resulting galaxy rotation curves are in very good agreement with observational data for galaxies with stellar mass $M_*<5\times10^{10}M_\odot$. We present an empirical fitting function that describes the total mass profiles and show that its parameters are strongly correlated with halo mass., Comment: 21 pages, 17 figures, version accepted for publication in MNRAS

We use the "Evolution and assembly of galaxies and their environments" (EAGLE) cosmological simulation to investigate the effect of baryons on the density profiles of rich galaxy clusters. We focus on EAGLE clusters with $M_{200}>10^{14}~M_\odot$ of which we have six examples. The central brightest cluster galaxies (BCGs) in the simulation have steep stellar density profiles, $\rho_*(r) \propto r^{-3}$. Stars dominate the mass density for $r<10~\rm{kpc}$, and, as a result, the total mass density profiles are steeper than the Navarro-Frenk-White (NFW) profile, in remarkable agreement with observations. The dark matter halo itself closely follows the NFW form at all resolved radii ($r\gtrsim3.0~\rm{kpc}$). The EAGLE BCGs have similar surface brightness and line-of-sight velocity dispersion profiles as the BCGs in the sample of Newman et al., which have the most detailed measurements currently available. After subtracting the contribution of the stars to the central density, Newman et al. infer significantly shallower slopes than the NFW value, in contradiction with the EAGLE results. We discuss possible reasons for this discrepancy, and conclude that an inconsistency between the kinematical model adopted by Newman et al. for their BCGs, which assumes isotropic stellar orbits, and the kinematical structure of the EAGLE BCGs, in which the orbital stellar anisotropy varies with radius and tends to be radially biased, could explain at least part of the discrepancy., Comment: 13 pages, 8 figures, accepted for publication in MNRAS

We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation., Comment: Accepted for publication in MNRAS. V2: Minor changes. For images and videos, see http://eagle.strw.leidenuniv.nl/ and http://icc.dur.ac.uk/Eagle/

We employ a large sample of 20171 optically-selected groups and clusters at 0.15 < z < 0.4 in the SDSS to investigate how the stacked stellar mass fraction varies across a wide range of total mass, $M_{500}$. Our study improves upon previous observational studies in a number of important ways, including having a much larger sample size, an explicit inclusion of the intracluster light (ICL) component, and a thorough examination of the accuracy of our total mass estimates via comparisons to simulations and weak lensing observations. We find that the stellar mass fraction depends only weakly on total mass and that the contribution of ICL to the total stellar mass fraction is significant (typically 20-40 per cent). Both of these findings are in excellent accordance with the predictions of cosmological simulations. Under the assumption of a Chabrier (Salpeter) IMF, the derived star formation efficiency ($f_{star}$/$f_{b}$, where $f_b=\Omega_b/\Omega_m$) is relatively low at 8 per cent (14 per cent) and is consistent with the global star formation efficiency of semi-analytic models that reproduce the galaxy stellar mass function. When our measured stellar mass fractions are combined with the observed relation between hot gas mass fraction and total mass from X-ray observations, our results imply that galaxy groups have significantly lower baryon fractions than massive clusters. Ejection of gas due to energetic AGN feedback (most likely at high redshift) provides a plausible mechanism for explaining the trends we observe., Comment: accepted to MNRAS; minor revisions

We present a new catalogue of 55,121 groups and clusters centred on Luminous Red Galaxies from SDSS DR7 in the redshift range 0.15~2.6) that is approximately a factor of two lower than that of the dark matter (as predicted by N-body cosmological simulations) and nearly independent of halo mass. Interestingly, in spite of the difference in shape between the galaxy and dark matter radial distributions, both exhibit a high degree of self-similarity. A self-consistent comparison to several recent semi-analytic models of galaxy formation indicates that: (1) beyond ~0.3 r_500 current models are able to reproduce both the shape and normalisation of the satellite profiles; and (2) within ~0.3 r_500 the predicted profiles are sensitive to the details of the satellite-BCG merger timescale calculation. The former is a direct result of the models being tuned to match the global galaxy luminosity function combined with the assumption that the satellite galaxies do not suffer significant tidal stripping, even though their surrounding dark matter haloes can be removed through this process. Combining our results with measurements of the intracluster light should provide a way to inform theoretical models on the efficacy of the tidal stripping and merging processes., Comment: Resubmitted to MNRAS with minor revisions

We estimate the masses of elliptical galaxies out to five effective radii using planetary nebulae and globular clusters as tracers. A sample of 15 elliptical galaxies with a broad variation in mass is compiled from the literature. A distribution function-maximum likelihood analysis is used to estimate the overall potential slope, normalisation and velocity anisotropy of the tracers. We assume power-law profiles for the potential and tracer density and a constant velocity anisotropy. The derived potential power-law indices lie in between the isothermal and Keplerian regime and vary with mass: there is tentative evidence that the less massive galaxies have steeper potential profiles than the more massive galaxies. We use stellar mass-to-light ratios appropriate for either a Chabrier/KTG (Kroupa, Tout & Gilmore) or Salpeter initial mass function to disentangle the stellar and dark matter components. The fraction of dark matter within five effective radii increases with mass, in agreement with several other studies. We employ simple models to show that a combination of star formation efficiency and baryon extent are able to account for this trend. These models are in good agreement with both our measurements out to five effective radii and recent SLACS measurements within one effective radii when a universal Chabrier/KTG initial mass function is adopted., Comment: 11 pages, 8 figures. Accepted for publication by ApJ

(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., Comment: 19 pages, 17 Figures. Accepted by MNRAS

We investigate the impact of AGN feedback on the gravitational lensing properties of a sample of galaxy clusters with masses in the range 10^14 - 10^15 Msol, using state-of-the-art simulations. Adopting a ray-tracing algorithm, we compute the cross-section of giant arcs from clusters simulated with dark-matter only physics (DM), dark matter plus gas with cooling and star formation (CSF), and dark matter plus gas with cooling, star formation and AGN feedback (CSFBH). Once AGN feedback is included, baryonic physics boosts the strong lensing cross-section by much less than previously estimated using clusters simulated with only cooling and star formation. For a cluster with a virial mass of 7.4 X 10^14 Msol, inclusion of baryonic physics without feedback can boost the cross-section by as much as a factor of 3, in agreement with previous studies, whereas once AGN feedback is included this maximal figure falls to a factor of 2 at most. Typically, clusters simulated with DM and CSFBH physics have similar cross-sections for the production of giant arcs. We also investigate how baryonic physics affects the weak lensing properties of the simulated clusters by fitting NFW profiles to synthetic weak lensing data sets using a Markov Chain Monte Carlo approach, and by performing non-parametric mass reconstructions. Without the inclusion of AGN feedback, measured concentration parameters can be much larger than those obtained with AGN feedback, which are similar to the dark-matter only case., Comment: 12 pages, 7 figures, Accepted to MNRAS

We investigate the physics driving the cosmic star formation (SF) history using the more than fifty large, cosmological, hydrodynamical simulations that together comprise the OverWhelmingly Large Simulations (OWLS) project. We systematically vary the parameters of the model to determine which physical processes are dominant and which aspects of the model are robust. Generically, we find that SF is limited by the build-up of dark matter haloes at high redshift, reaches a broad maximum at intermediate redshift, then decreases as it is quenched by lower cooling rates in hotter and lower density gas, gas exhaustion, and self-regulated feedback from stars and black holes. The higher redshift SF is therefore mostly determined by the cosmological parameters and to a lesser extent by photo-heating from reionization. The location and height of the peak in the SF history, and the steepness of the decline towards the present, depend on the physics and implementation of stellar and black hole feedback. Mass loss from intermediate-mass stars and metal-line cooling both boost the SF rate at late times. Galaxies form stars in a self-regulated fashion at a rate controlled by the balance between, on the one hand, feedback from massive stars and black holes and, on the other hand, gas cooling and accretion. Paradoxically, the SF rate is highly insensitive to the assumed SF law. This can be understood in terms of self-regulation: if the SF efficiency is changed, then galaxies adjust their gas fractions so as to achieve the same rate of production of massive stars. Self-regulated feedback from accreting black holes is required to match the steep decline in the observed SF rate below redshift two, although more extreme feedback from SF, for example in the form of a top-heavy IMF at high gas pressures, can help., Comment: Accepted for publication in MNRAS, 27 pages and 18 figures. Revised version: minor changes

(Abridged) The thermal state of the intracluster medium results from a competition between gas cooling and heating. The heating comes from two distinct sources: gravitational heating from the collapse of the dark matter halo and thermal input from galaxy/black hole formation. However, a long standing problem has been that cosmological simulations based on smoothed particle hydrodynamics (SPH) and Eulerian mesh codes predict different results even when cooling and galaxy/black hole heating are switched off. Clusters formed in SPH simulations show near powerlaw entropy profiles, while those formed in mesh simulations develop a core and do not allow gas to reach such low entropies. Since the cooling rate is closely connected to the minimum entropy of the gas, the differences are of potentially key importance. In this paper, we investigate the origin of this discrepancy. By comparing simulations run using the GADGET-2 SPH code and the FLASH adaptive Eulerian mesh code, we show that the discrepancy arises during the idealised merger of two clusters. The difference is not sensitive to the resolution of our simulations, nor is it is due differences in the gravity solvers, Galilean non-invariance of the mesh code, or an effect of unsuitable artificial viscosity in the SPH code. Instead, we find that the difference is inherent to the treatment of eddies and fluid instabilities. These are suppressed in the SPH simulations, while the cluster mergers generate strong vortices in the mesh simulations that efficiently mix the fluid and erase the low entropy gas. Consequently, particles in the SPH simulations retain a close connection to their initial entropy, while this connection is much weaker in the mesh simulations. We discuss the potentially profound implications of these results., Comment: 19 pages, 10 figures, Accepted for publication in MNRAS

Using a set of high-resolution N-body/SPH cosmological simulations with identical initial conditions but run with different numerical setups, we investigate the influence of baryonic matter on the mass distribution of dark halos when radiative cooling is NOT included. We compare the concentration parameters of about 400 massive halos with virial mass from $10^{13}$ \Msun to $7.1 \times 10^{14}$ \Msun. We find that the concentration parameters for the total mass and dark matter distributions in non radiative simulations are on average larger by ~3% and 10% than those in a pure dark matter simulation. Our results indicate that the total mass density profile is little affected by a hot gas component in the simulations. After carefully excluding the effects of resolutions and spurious two-body heating between dark matter and gas particles, we conclude that the increase of the dark matter concentration parameters is due to interactions between baryons and dark matter. We demonstrate this with the aid of idealized simulations of two-body mergers. The results of individual halos simulated with different mass resolutions show that the gas profiles of densities, temperature and entropy are subjects of mass resolution of SPH particles. In particular, we find that in the inner parts of halos, as the SPH resolution increases the gas density becomes higher but both the entropy and temperature decrease., Comment: 8 pages, 6 figures, 1 table, ApJ in press (v652n1); updated to match with the being published version

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., Comment: 4 pages, 5 figures; invited paper to appear in proceedings of 'Bodenheimer Fest'