We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from the simulations, we find that the inner density slopes of both the total and the dark matter mass distributions can be inferred reasonably well. We compare the density slopes of C-EAGLE clusters with those derived by Newman et al. for 7 massive galaxy clusters in the local Universe. We find that the asymptotic best-fit inner slopes of "generalized" NFW (gNFW) profiles, ${\gamma}_{\rm gNFW}$, of the dark matter haloes of the C-EAGLE clusters are significantly steeper than those inferred by Newman et al. However, the mean mass-weighted dark matter density slopes of the simulated clusters are in good agreement with the Newman et al. estimates. We also find that the estimate of ${\gamma}_{\rm gNFW}$ is very sensitive to the constraints from weak lensing measurements in the outer parts of the cluster and a bias can lead to an underestimate of ${\gamma}_{\rm gNFW}$.

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

We perform a set of non-radiative hydro-dynamical (NHD) simulations of a rich cluster sized dark matter halo from the Phoenix project with 3 different numerical resolutions, to investigate the effect of hydrodynamics alone on the subhalo population in the halo. Compared to dark matter only (DMO) simulations of the same halo, subhaloes are less abundant for relatively massive subhaloes ($M_{sub} > 2.5 \times 10^9h^{-1}M_{\odot}$, or $V_{max} > 70 kms^{-1}$) but more abundant for less massive subhaloes in the NHD simulations. This results in different shapes in the subhalo mass/$V_{max}$ function in two different sets of simulations. At given subhalo mass, the subhaloes less massive than $10^{10} h^{-1}M_{\odot}$ have larger $V_{max}$ in the NHD than DMO simulations, while $V_{max}$ is similar for the subhaloes more massive than the mass value. This is mainly because the progenitors of present day low mass subhaloes have larger concentration parameters in the NHD than DMO simulations. The survival number fraction of the accreted low mass progenitors of the main halo at redshift 2 is about 50 percent higher in the NHD than DMO simulations.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

Kilonovae are optical flashes produced in the aftermath of neutron star-neutron star mergers (NNMs) or neutron star-black hole mergers (NBMs). In this work, we use the Millennium Simulation, combined with a semi-analytic galaxy formation model--GABE (Galaxy Assembly with Binary Evolution) which adopts binary stellar population synthesis models, to explore the cosmic event rate of kilonovae, and the properties of their host galaxies in a cosmological context. We find that model with supernova kick velocity of 0 km/s fits the observation best, in agreement with the exception of some formation channels of binary neutron star. This indicates that NNMs prefer to originate from binary systems with low kick velocities. With V$_{\rm kick}$=0 km/s, the cosmic event rate of NNMs and NBMs at z=0 are 283 Gpc$^{-3}$yr$^{-1}$ and 91 Gpc$^{-3}$yr$^{-1}$, respectively, marginally consistent with the constraint from LVC GWTC-1. For Milky Way-mass galaxies, we predict the NNM rate is $25.7^{+59.6}_{-7.1}$ Myr$^{-1}$, which is also in good agreement with the observed properties of binary neutron stars in the Milky Way. Taking all the NNMs into account in the history of Milky Way-mass galaxies, we find that the averaged r-process elements yield with A>79 in a NNM and NBM event should be 0.01 M$_{\odot}$ to be consistent with observation. We conclude that NGC 4993, the host galaxy of GW170817, is a typical host galaxy for NNMs. However, generally NNMs and NBMs tend to reside in young, blue, star-forming, late-type galaxies, with stellar mass and gaseous metallicity distribution peaking at $10^{10.65}$ M$_{\odot}$ and 8.72-8.85, respectively. By studying kilonovae host galaxies in the cosmological background, it is promising to constrain model details better when we have more events in the forthcoming future. (abridged)

The thermal history of cosmic gas in the Dark Ages remains largely unknown. It is important to quantify the impact of relevant physics on the IGM temperature between $z=10$ and $z \sim 30$, in order to interpret recent and oncoming observations, including results reported by EDGES. We revisit the gas heating due to structure formation shocks in this era, using a set of fixed grid cosmological hydrodynamical simulations performed by three different codes. In all our simulations, the cosmic gas is predicted to be in multiphase state since $z>30$. The gas surrounding high density peaks gradually develops a relation more sharp than $T \propto \rho^{2/3}$, approximately $T \propto \rho^{2}$, from $z=30$ to $z=11$, might due to shock heating. Meanwhile, the gas in void region tends to have a large local mach number, and their thermal state varies significantly from code to code. In the redshift range $11-20$, the mass fraction of gas shock heated above the CMB temperature in our simulations is larger than previous semi-analytical results by a factor of 2 to 8. At $z=15$, the fraction varies from $\sim 19\%$ to $52 \%$ among different codes. Between $z=11$ and $z=20$, the gas temperature $<1/T_{\rm{K}}>_M^{-1}$ is predicted to be $\sim 10-20$ K by two codes, much higher than the adiabatic cooling model and some previous works. However, in our simulations performed by RAMSES, $<1/T_{\rm{K}}>_M^{-1}$ is predicted to be even below the temperature required to explain result of the EDGES. Given the fact that different codes give different predictions, currently, it seems a challenge to make solid prediction on the temperature of gas at $z \sim 17$ in simulations.

In the standard cosmological model, dark matter drives the structure formation and constructs potential wells within which galaxies may form. The baryon fraction in dark halos can reach the universal value (15.7%) in massive clusters and decreases rapidly as the mass of the system decreases. The formation of dwarf galaxies is sensitive both to baryonic processes and the properties of dark matter owing to the shallow potential wells in which they form. In dwarf galaxies in the Local Group, dark matter dominates the mass content even within their optical-light half-radii (r_e ~ 1 kpc). However, recently it has been argued that not all dwarf galaxies are dominated by dark matter. Here we report 19 dwarf galaxies that could consist mainly of baryons up to radii well beyond r_e, at which point they are expected to be dominated by dark matter. Of these, 14 are isolated dwarf galaxies, free from the influence of nearby bright galaxies and high dense environments. This result provides observational evidence that could challenge the formation theory of low-mass galaxies within the framework of standard cosmology. Further observations, in particular deep imaging and spatially-resolved kinematics, are needed to constrain the baryon fraction better in such galaxies.

Dark matter haloes are the basic units of all cosmic structure. They grew by gravitational amplification of weak initial density fluctuations that are still visible on large scales in the cosmic microwave background radiation. Galaxies formed within relatively massive haloes as gas cooled and condensed at their centres, but many hypotheses for the nature of dark matter imply that the halo population should extend to masses many orders of magnitude below those where galaxies can form. Here, we use a novel, multi-zoom technique to create the first consistent simulation of the formation of present-day haloes over the full mass range populated when dark matter is aWeakly Interacting Massive Particle (WIMP) of mass ~100 GeV. The simulation has a dynamic range of 30 orders of magnitude in mass, resolving the internal structure of hundreds of Earth-mass haloes just as well as that of hundreds of rich galaxy clusters. Remarkably, halo density profiles are universal over the entire mass range and are well described by simple two-parameter fitting formulae. Halo mass and concentration are tightly related in a way which depends on cosmology and on the nature of the dark matter. At fixed mass, concentration is independent of local environment for haloes less massive than those of typical galaxies. These results are important for predicting annihilation radiation signals from dark matter, since these should be dominated by contributions from the smallest structures.

Future observations of 21 cm emission from neutral hydrogen survey will become a promising approach to probe the large scale structure of Universe. In this paper, we investigate the impacts of Square Kilometer Array (SKA) 21 cm observation on the estimation of cosmological parameters. We use the simulated data of the baryonic acoustic oscillation (BAO) measurements based on the future SKA experiment with the intensity mapping (IM) technique to do the analysis. For the current observations, we use the latest cosmic microwave background (CMB) observation from {\it Planck} 2018, the optical BAO measurements, and the Type Ia supernovae (SN) observation (Pantheon compilation). We find that the SKA mock data could break the degeneracy between the matter density and the Hubble constant, further improving the cosmological constraints to a great extent. We also find that the constraint on the equation of state parameters of dark energy could be significantly improved by including the SKA mock data into the cosmological global fit.

We investigate the formation of ultra-diffuse galaxies (UDGs) using the Auriga high-resolution cosmological magneto-hydrodynamical simulations of Milky Way-sized galaxies. We identify a sample of $92$ UDGs in the simulations that match a wide range of observables such as sizes, central surface brightness, S\'{e}rsic indices, colors, spatial distribution and abundance. Auriga UDGs have dynamical masses similar to normal dwarfs. In the field, the key to their origin is a strong correlation present in low-mass dark matter haloes between galaxy size and halo spin parameter. Field UDGs form in dark matter haloes with larger spins compared to normal dwarfs in the field, in agreement with previous semi-analytical models. Satellite UDGs, on the other hand, have two different origins: $\sim 55\%$ of them formed as field UDGs before they were accreted; the remaining $\sim 45\%$ were normal field dwarfs that subsequently turned into UDGs as a result of tidal interactions.

The neutral hydrogen (HI) and its 21 cm line are promising probes to the reionization process of the intergalactic medium (IGM). To use this probe effectively, it is imperative to have a good understanding on how the neutral hydrogen traces the underlying matter distribution. Here we study this problem using semi-numerical modeling by combining the HI in the IGM and the HI from halos during the epoch of reionization (EoR), and investigate the evolution and the scale-dependence of the neutral fraction bias as well as the 21 cm line bias. We find that the neutral fraction bias on large scales is negative during reionization, and its absolute value on large scales increases during the early stage of reionization and then decreases during the late stage. During the late stage of reionization, there is a transition scale at which the HI bias transits from negative on large scales to positive on small scales, and this scale increases as the reionization proceeds to the end.

We introduce a new halo/subhalo finder, HIKER (a Halo fInder based on KERnel-shift algorithm), which takes advantage of a machine learning method -- the mean-shift algorithm combined with the Plummer kernel function, to effectively locate density peaks corresponding to halos/subhalos in density field. Based on these density peaks, dark matter halos are identified as spherical overdensity structures, and subhalos are bound substructures with boundaries at their tidal radius. By testing HIKER code with mock halos, we show that HIKER performs excellently in recovering input halo properties. Especially, HIKER has higher accuracy in locating halo/subhalo centres than most halo finders. With cosmological simulations, we further show that HIKER reproduces the abundance of dark matter halos and subhalos quite accurately, and the HIKER halo/subhalo mass functions and $V_{max}$ functions are in good agreement with two widely used halo finders, SUBFIND and AHF.

Low mass galaxies are expected to be dark matter dominated even within their centrals. Recently two observations reported two dwarf galaxies in group environment with very little dark matter in their centrals. We explore the population and origins of dark-matter-deficient galaxies (DMDGs) in two state-of-the-art hydrodynamical simulations, the EAGLE and Illustris projects. For all satellite galaxies with $10^9<M_*<10^{10}$ M$_{\odot}$ in groups with $M_{200}>10^{13}$ M$_{\odot}$, we find that about $2.6\%$ of them in the EAGLE, and $1.5\%$ in the Illustris are DMDGs with dark matter fractions below $50\%$ inside two times half-stellar-mass radii. We demonstrate that DMDGs are highly tidal disrupted galaxies; and because dark matter has higher binding energy than stars, mass loss of the dark matter is much more rapid than stars in DMDGs during tidal interactions. If DMDGs were confirmed in observations, they are expected in current galaxy formation models.

We developed a new semi-analytic galaxy formation model: Galaxy Assembly with Binary Evolution (GABE). For the first time we introduce binary evolution into semi-analytic models of galaxy formation by using Yunnan-II stellar population synthesis model, which includes various binary interactions. When implementing our galaxy formation model onto the merger trees extracted from the Millennium simulation, it can reproduce a large body of observational results. We find that in the local universe the model including binary evolution reduces the luminosity at optical and infrared wavelengths slightly, while increases the luminosity at ultraviolet wavelength significantly, especially in $F_{\rm UV}$ band. The resulting luminosity function does not change very much over SDSS optical bands and infrared band, but the predicted colors are bluer, especially when $F_{\rm UV}$ band is under consideration. The new model allows us to explore the physics of various high energy events related to the remnants of binary stars, e.g. type Ia supernovae, short gamma-ray bursts and gravitational wave events, and their relation with host galaxies in a cosmological context.

Gravitational softening length is one of the key parameters to properly set up a cosmological $N$-body simulation. In this paper, we perform a large suit of high-resolution $N$-body simulations to revise the optimal softening scheme proposed by Power et al. (P03). Our finding is that P03 optimal scheme works well but is over conservative. Using smaller softening lengths than that of P03 can achieve higher spatial resolution and numerically convergent results on both circular velocity and density profiles. However using an over small softening length overpredicts matter density at the inner most region of dark matter haloes. We empirically explore a better optimal softening scheme based on P03 form and find that a small modification works well. This work will be useful for setting up cosmological simulations.

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.

We derive the mass weighted total density slopes within the effective (half-light) radius, $\gamma'$, for more than 2000 nearby galaxies from the SDSS-IV MaNGA survey using Jeans-anisotropic-models applied to IFU observations. Our galaxies span a wide range of the stellar mass ($10^9$ $M_{\rm \odot}< M_* < 10^{12}$ M$_{\odot}$) and the velocity dispersion (30 km/s $< \sigma_v <$ 300 km/s). We find that for galaxies with velocity dispersion $\sigma_v>100$ km/s, the density slope has a mean value $\langle \gamma^{\prime} \rangle = 2.24$ and a dispersion $\sigma_{\gamma}=0.22$, almost independent of velocity dispersion. A clear turn over in the $\gamma'-\sigma_v$ relation is present at $\sigma\sim 100$ km/s, below which the density slope decreases rapidly with $\sigma_v$. Our analysis shows that a large fraction of dwarf galaxies (below $M_* = 10^{10}$ M$_{\odot}$) have total density slopes shallower than 1, which implies that they may reside in cold dark matter halos with shallow density slopes. We compare our results with that of galaxies in hydrodynamical simulations of EAGLE, Illustris and IllustrisTNG projects, and find all simulations predict shallower density slopes for massive galaxies with high $\sigma_v$. Finally, we explore the dependence of $\gamma'$ on the positions of galaxies in halos, namely centrals vs. satellites, and find that for the same velocity dispersion, the amplitude of $\gamma'$ is higher for satellite galaxies by about 0.1.

We compare galaxy morphology predicted by the Illustris hydrodynamical simulation and a semi-analytic model (SAM) grafted in the halo merger trees from the Illustris-Dark matter simulation. Morphology is classified according to the luminous profile and the kinematic bulge-to-total ratio for Illustris galaxy, and the bulge-to-total stellar mass ratio for SAM galaxy. For late-type galaxies in the Illustris catalogue, most of their counterparts in the SAM model have the same type, and the consistency between two models is higher for lower mass galaxies. For early-type (ET) galaxies in Illustris, the consistency is quite low for the counterparts except for most massive galaxies. By comparing in detail the growth histories of some matched galaxy pairs of Milky Way mass in Illustris and the SAM model, we notice two aspects of differences in determining galaxy morphology between the two galaxy formation implementations. First, in the SAM, major merger and frequent minor mergers result in the growth of bulges and turn the galaxy into ET, while bulge formation is not connected to mergers as tightly as in SAM for the Illustris galaxies. In addition, the satellite stellar mass can decrease significantly due to tidal stripping before merging into the central galaxy in Illustris, while it does not decrease in the SAM model. This results in less mergers with large (stellar) mass ratios in the Illustris simulation, and less effect of mergers on shaping galaxy morphology.

We make use of a high-resolution zoom-in hydrodynamical simulation to investigate the impact of filaments on galaxy formation in their residing dark matter haloes. A method based on the density field and the Hoshen-Kopelman algorithm is developed to identify filaments. We show that cold and dense gas preprocessed by dark matter filaments can be further accreted into residing individual low-mass haloes in directions along the filaments. Consequently, comparing with field haloes, gas accretion is very anisotropic for filament haloes. About 30 percent of the accreted gas of a residing filament halo was preprocessed by filaments, leading to two different thermal histories for the gas in filament haloes. Filament haloes have higher baryon and stellar fractions when comparing with their field counterparts. Without including stellar feedback, our results suggest that filaments assist gas cooling and enhance star formation in their residing dark matter haloes at high redshifts (i.e. z=4.0 and 2.5).

The concordance $\Lambda$CDM cosmology predicts tens of satellite galaxies distributed in the inner region ($<40\ \lkpc$) of the Milky Way (MW), yet at most $12$ were discovered at present day, including 3 discovered very recently by Dark Energy Survey (DES) and 5 from other surveys (e.g. PanSTARRS, MagLiteS). We use $5$ ultra-high resolution simulations of MW sized dark matter haloes from the $\aqua$ project, combined with $\galf$ semi-analytical galaxy formation model, to investigate properties of the model satellite galaxy population inside $40\ \lkpc$ of MW sized haloes. On average, in each halo this model predicts about $20$ inner satellite galaxies, among them $5$ are comparable to the classic satellites in the luminosity, these are in stark contrast to the corresponding numbers in observations. We further investigate the survivability of these model inner satellites in the presence of a central stellar disk with a set of ideal simulations. These are done by re-evolving a quarter (30) of the whole $\aqua$ inner satellite galaxies (121) by including a static disk potential in addition to the MW halo. Our finding is that the additional disk completely disrupt 40 percents of these satellites and results in $14$ satellite galaxies within the $40\ \lkpc$ of each $\aqua$ at the end, in reasonably well agreement with observations.

Schwinn et al. (2017) have recently compared the abundance and distribution of massive substructures identified in a gravitational lensing analysis of Abell 2744 by Jauzac et al. (2016) and N-body simulation and found no cluster in {\Lambda}CDM simulation that is similar to Abell 2744. Schwinn et al.(2017) identified the measured projected aperture masses with the actual masses associated with subhaloes in the MXXL N-body simulation. We have used the high resolution Phoenix cluster simulations to show that such an identification is incorrect: the aperture mass is dominated by mass in the body of the cluster that happens to be projected along the line-of-sight to the subhalo. This enhancement varies from factors of a few to factors of more than 100, particularly for subhaloes projected near the centre of the cluster. We calculate aperture masses for subhaloes in our simulation and compare them to the measurements for Abell 2744. We find that the data for Abell 2744 are in excellent agreement with the matched predictions from {\Lambda}CDM. We provide further predictions for aperture mass functions of subhaloes in idealized surveys with varying mass detection thresholds.

We use a volume-limited galaxy sample from the SDSS Data Release 7 to explore the dependence of galactic conformity on the large-scale environment, measured on $\sim$ 4 Mpc scales. We find that the star formation activity of neighbour galaxies depends more strongly on the environment than on the activity of their primary galaxies. In under-dense regions most neighbour galaxies tend to be active, while in over-dense regions neighbour galaxies are mostly passive, regardless of the activity of their primary galaxies. At a given stellar mass, passive primary galaxies reside in higher density regions than active primary galaxies, leading to the apparently strong conformity signal. The dependence of the activity of neighbour galaxies on environment can be explained by the corresponding dependence of the fraction of satellite galaxies. Similar results are found for galaxies in a semi-analytical model, suggesting that no new physics is required to explain the observed large-scale conformity.

We study the radial acceleration relation (RAR) for early-type galaxies (ETGs) in the SDSS MaNGA MPL5 dataset. The complete ETG sample show a slightly offset RAR from the relation reported by McGaugh et al. (2016) at the low-acceleration end; we find that the deviation is due to the fact that the slow rotators show a systematically higher acceleration relation than the McGaugh's RAR, while the fast rotators show a consistent acceleration relation to McGaugh's RAR. There is a 1\sigma significant difference between the acceleration relations of the fast and slow rotators, suggesting that the acceleration relation correlates with the galactic spins, and that the slow rotators may have a different mass distribution compared with fast rotators and late-type galaxies. We suspect that the acceleration relation deviation of slow rotators may be attributed to more galaxy merger events, which would disrupt the original spins and correlated distributions of baryons and dark matter orbits in galaxies.

We introduce a new code for cosmological simulations, PHoToNs, which has features on performing massive cosmological simulations on heterogeneous high performance Computer (HPC) and threads oriented programming. PHoToNs adopts a hybrid scheme to compute gravity force, with the conventional PM to compute the long-range force, the Tree algorithm to compute the short range force, and the direct summation PP to compute the gravity from very close particles. A self-similar space filling Peano-Hilbert curve is used to decompose computing domain. Threads programming is highly used to more flexibly manage the domain communication, PM calculation and synchronization, as well as Dual Tree Traversal on the CPU+MIC platform. The scalability of the PHoToNs performs well and the efficiency of PP kernel achieves 68.6% of peak performance on MIC and 74.4% on CPU platforms. We also test the accuracy of the code against the much used Gadget-2 in the community and found excellent agreement.

The spatial distribution of satellite galaxies encodes rich information of the structure and assembly history of galaxy clusters. In this paper, we select a redMaPPer cluster sample in SDSS Stripe 82 region with $0.1 \le z \le 0.33$, $20<\lambda<100$ and $P_{\rm cen}>0.7$. Using the high-quality weak lensing data from CS82 Survey, we constrain the mass profile of this sample. Then we compare directly the mass density profile with the satellite number density profile. We find that the total mass and number density profiles have the same shape, both well fitted by an NFW profile. The scale radii agree with each other within 1$\sigma$ error ($r_{\rm s,gal}=0.34_{-0.03}^{+0.04}$Mpc vs $r_{\rm s}=0.37_{-0.10}^{+0.15}$Mpc ).

We take advantage of the statistical power of the large-volume dark-matter-only Millennium simulation, combined with a sophisticated semi-analytic galaxy formation model, to explore whether the recently reported $z=3.7$ quiescent galaxy ZF-COSMOS-20115 (ZF; Glazebrook et al. 2017) can be accommodated in current galaxy formation models. In our model, a population of quiescent galaxies (QGs) with stellar masses and star formation rates comparable to those of ZF naturally emerges at redshifts $z<4$. There are two and five ZF analogues at the redshift $3.86$ and $3.58$ in the Millennium simulation volume, respectively. We demonstrate that, while the $z>3.5$ massive QGs are rare (about 2\% of the galaxies with the similar stellar masses), the existing AGN feedback model implemented in the semi-analytic galaxy formation model can successfully explain the formation of the high-redshift QGs as it does on their lower redshift counterparts.

The standard galaxy formation theory assumes that baryons and dark matter are initially well-mixed before becoming segregated due to radiative cooling. We use non-radiative hydrodynamical simulations to explicitly examine this assumption and find that baryons and dark matter can also be segregated because of different physics obeyed by gas and dark matter during the build-up of the halo. As a result, baryons in many haloes do not originate from the same Lagrangian region as the dark matter. When using the fraction of corresponding dark matter and gas particles in the initial conditions (the "paired fraction") as a proxy of the dark matter and gas segregation strength of a halo, on average about $25$ percent of the baryonic and dark matter of the final halo are segregated in the initial conditions. This is at odds with the assumption of the standard galaxy formation model. A consequence of this effect is that the baryons and dark matter of the same halo initially experience different tidal torques and thus their angular momentum vectors are often misaligned. The degree of the misalignment is largely preserved during later halo assembly and can be understood with the tidal torque theory. The result challenges the precision of some semi-analytical approaches which utilize dark matter halo merger trees to infer properties of gas associated to dark matter haloes.

A particular population of galaxies have drawn much interest recently, which are as faint as typical dwarf galaxies but have the sizes as large as $L^*$ galaxies, the so called "ultra-diffuse galaxie" (UDGs). The lack of tidal features of UDGs in dense environments suggests that their host halos are perhaps as massive as that of the Milky Way. On the other hand, galaxy formation efficiency should be much higher in the halos of such masses. Here we use the model galaxy catalog generated by populating two large simulations: the Millennium-II cosmological simulation and Phoenix simulations of 9 big clusters with the semi-analytic galaxy formation model. This model reproduces remarkably well the observed properties of UDGs in the nearby clusters, including the abundance, profile, color, and morphology, etc. We search for UDG candidates using the public data and find 2 UDG candidates in our Local Group and 23 in our Local Volume, in excellent agreement with the model predictions. We demonstrate that UDGs are genuine dwarf galaxies, formed in the halos of $\sim 10^{10}M_{\odot}$. It is the combination of the late formation time and high-spins of the host halos that results in the spatially extended feature of this particular population. The lack of tidal disruption features of UDGs in clusters can also be explained by their late infall-time.

We use a pair of high resolution N-body simulations implementing two dark matter models, namely the standard cold dark matter (CDM) cosmogony and a warm dark matter (WDM) alternative where the dark matter particle is a 1.5keV thermal relic. We combine these simulations with the GALFORM semi-analytical galaxy formation model in order to explore differences between the resulting galaxy populations. We use GALFORM model variants for CDM and WDM that result in the same z=0 galaxy stellar mass function by construction. We find that most of the studied galaxy properties have the same values in these two models, indicating that both dark matter scenarios match current observational data equally well. Even in under-dense regions, where discrepancies in structure formation between CDM and WDM are expected to be most pronounced, the galaxy properties are only slightly different. The only significant difference in the local universe we find is in the galaxy populations of "Local Volumes", regions of radius 1 to 8Mpc around simulated Milky Way analogues. In such regions our WDM model provides a better match to observed local galaxy number counts and is five times more likely than the CDM model to predict sub-regions within them that are as empty as the observed Local Void. Thus, a highly complete census of the Local Volume and future surveys of void regions could provide constraints on the nature of dark matter.

The defining characteristic of the cold dark matter (CDM) hypothesis is the presence of a very large number of low-mass haloes, too small to have made a visible galaxy. Other hypotheses for the nature of the dark matter, such as warm dark matter (WDM), predict a much smaller number of such low-mass haloes. Strong lensing systems offer the possibility of detecting small-mass haloes through the distortions they induce in the lensed image. Here we show that the main contribution to the image distortions comes from haloes along the line of sight rather than subhaloes in the lens as has normally been assumed so far. These interlopers enhance the differences between the predictions of CDM and WDM models. We derive the total perturber mass function, including both subhaloes and interlopers, and show that measurements of approximately 20 strong lens systems with a detection limit of $M_{\rm low}=10^7 h^{-1} M_{\odot}$ would distinguish (at 3 sigma) between CDM and a WDM model consisting of 7 keV sterile neutrinos such as those required to explain the recently detected 3.5 keV X-ray emission line from the centres of galaxies and clusters.

We use the Copernicus Complexio (COCO) high resolution $N$-body simulations to investigate differences in the properties of small-scale structures in the standard cold dark matter (CDM) model and in a model with a cutoff in the initial power spectrum of density fluctuations consistent with both a thermally produced warm dark matter (WDM) particle or a sterile neutrino with mass 7 keV and leptogenesis parameter $L_6=8.7$. The latter corresponds to the "coldest" model with this sterile neutrino mass compatible with the identification of the recently detected 3.5 keV X-ray line as resulting from particle decay. CDM and WDM predict very different number densities of subhaloes with mass $\leq 10^9\,h^{-1}\,M_\odot$ although they predict similar, nearly universal, normalised subhalo radial density distributions. Haloes and subhaloes in both models have cuspy NFW profiles, but WDM subhaloes below the cutoff scale in the power spectrum (corresponding to maximum circular velocities $V_{\mathrm{max}}^{z=0} \leq50~\mathrm{kms}^{-1}$) are less concentrated than their CDM counterparts. We make predictions for observable properties using the GALFORM semi-analytic model of galaxy formation. Both models predict Milky Way satellite luminosity functions consistent with observations, although the WDM model predicts fewer very faint satellites. This model, however, predicts slightly more UV bright galaxies at redshift $z>7$ than CDM, but both are consistent with observations. Gravitational lensing offers the best prospect of distinguishing between the models.

The spatial distribution of the satellite populations of the Milky Way and Andromeda are puzzling in that they are nearly perpendicular to the disks of their central galaxies. To understand the origin of such configurations we study the alignment of the central galaxy, satellite system and dark matter halo in the largest of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation. We find that centrals and their satellite systems tend to be well aligned with their haloes, with a median misalignment angle of $33^{\circ}$ in both cases. While the centrals are better aligned with the inner $10$ kpc halo, the satellite systems are better aligned with the entire halo indicating that satellites preferentially trace the outer halo. The central - satellite alignment is weak (median misalignment angle of $52^{\circ}$) and we find that around $20\%$ of systems have a misalignment angle larger than $78^{\circ}$, which is the value for the Milky Way. The central - satellite alignment is a consequence of the tendency of both components to align with the dark matter halo. As a consequence, when the central is parallel to the satellite system, it also tends to be parallel to the halo. In contrast, if the central is perpendicular to the satellite system, as in the case of the Milky Way and Andromeda, then the central - halo alignment is much weaker. Dispersion-dominated (spheroidal) centrals have a stronger alignment with both their halo and their satellites than rotation-dominated (disk) centrals. We also found that the halo, the central galaxy and the satellite system tend to be aligned with the surrounding large-scale distribution of matter, with the halo being the better aligned of the three.

The cold dark matter (CDM) cosmological model unambigously predicts that a large number of haloes should survive as subhaloes when they are accreted into a larger halo. The CDM model would be ruled out if such substructures were shown not to exist. By contrast, if the dark matter consists of warm particles (WDM), then below a threshold mass that depends on the particle mass far fewer substructures would be present. Finding subhaloes below a certain mass would then rule out warm particle masses below some value. Strong gravitational lensing provides a clean method to measure the subhalo mass function through distortions in the structure of Einstein rings and giant arcs.Using mock lensing observations constructed from high-resolution N-body simulations, we show that measurements of approximately 100 strong lens systems with a detection limit of $10^7 h^{-1} M_{\odot}$ would clearly distinguish CDM from WDM in the case where this consists of 7 keV sterile neutrinos such as those that might be responsible for the 3.5 keV X-ray emission line recently detected in galaxies and clusters.

ARGO-YBJ, located at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l., Tibet, China), is a full coverage air shower array, with an energy threshold of 300 GeV for gamma-ray astronomy. Most of the recorded events are single front showers, satisfying the trigger requirement of at least 20 particles detected in a given time window. However, in 13% of the events, two randomly arriving showers may be recorded in the same time window, and the second one, in generally smaller, does not need to satisfy the trigger condition. These events are called double front shower events. By using these small showers, well under the trigger threshold, the detector primary energy threshold can be lowered to a few tens of GeV. In this paper, the angular resolution that can be achieved with these events is evaluated by a full Monte Carlo simulation. The ARGO-YBJ sensitivity in detecting gamma-ray bursts (GRBs) by using double front shower events is also studied for various cutoff energies, time durations, and zenith angles of GRBs in the field view of ARGO.

We explore the Minkowski functionals of weak lensing convergence map to distinguish between $f(R)$ gravity and the general relativity (GR). The mock weak lensing convergence maps are constructed with a set of high-resolution simulations assuming different gravity models. It is shown that the lensing MFs of $f(R)$ gravity can be considerably different from that of GR because of the environmentally dependent enhancement of structure formation. We also investigate the effect of lensing noise on our results, and find that it is likely to distinguish F5, F6 and GR gravity models with a galaxy survey of $\sim3000$ degree$^2$ and with a background source number density of $n_g=30~{\rm arcmin}^{-2}$, comparable to an upcoming survey dark energy survey (DES). We also find that the $f(R)$ signal can be partially degenerate with the effect of changing cosmology, but combined use of other observations, such as the cosmic microwave background (CMB) data, can help break this degeneracy.

We make use of two suits of ultra high resolution N-body simulations of individual dark matter haloes from the Phoenix and the Aquarius Projects to investigate systematics of assembly history of subhaloes in dark matter haloes differing by a factor of $1000$ in the halo mass. We have found that real progenitors which built up present day subhalo population are relatively more abundant for high mass haloes, in contrast to previous studies claiming a universal form independent of the host halo mass. That is mainly because of repeated counting of the 're-accreted' (progenitors passed through and were later re-accreted to the host more than once) and inclusion of the 'ejected' progenitor population(progenitors were accreted to the host in the past but no longer members at present day) in previous studies. The typical accretion time for all progenitors vary strongly with the host halo mass, which is typical about $z \sim 5$ for the galactic Aquarius and about $z \sim 3$ for the cluster sized Phoenix haloes. Once these progenitors start to orbit their parent haloes, they rapidly lose their original mass but not their identifiers, more than $55$ ($50$) percent of them survive to present day for the Phoenix(Aquarius) haloes. At given redshift, survival fraction of the accreted subhalo is independent of the parent halo mass, whilst the mass-loss of the subhalo is more efficient in high mass haloes. These systematics results in similarity and difference in the subhalo population in dark matter haloes of different masses at present day.

We propose a new cosmological test of gravity, by using the observed mass fraction of X-ray emitting gas in massive galaxy clusters. The cluster gas fraction, believed to be a fair sample of the average baryon fraction in the Universe, is a well-understood observable, which has previously mainly been used to constrain background cosmology. In some modified gravity models, such as $f(R)$ gravity, gas temperature in a massive cluster is determined by the effective mass of that cluster, which can be larger than its true mass. On the other hand, X-ray luminosity is determined by the true gas density, which in both modified gravity and $\Lambda$CDM models depends mainly on $\Omega_{\rm b}/\Omega_{\rm m}$ and hence the true total cluster mass. As a result, the standard practice of combining gas temperatures and X-ray surface brightnesses of clusters to infer their gas fractions can, in modified gravity models, lead to a larger - in $f(R)$ gravity this can be $1/3$ larger - value of $\Omega_{\rm b}/\Omega_{\rm m}$ than that inferred from other observations such as the CMB. A quick calculation shows that the Hu-Sawicki $n=1$ $f(R)$ model with $|\bar{f}_{R0}|=3\sim5\times10^{-5}$ is in tension with the gas fraction data of the 42 clusters analysed by Allen et al. (2008). We also discuss the implications for other modified gravity models.

We investigate the properties of dark matter haloes and subhaloes in an $f(R)$ gravity model with $|f_{R0}|=10^{-6}$, using a very high-resolution N-body simulation. The model is a borderline between being cosmologically interesting and yet still consistent with current data. We find that the halo mass function in this model has a maximum 20% enhancement compared with the $\Lambda$CDM predictions between $z=1$ and $z=0$. Because of the chameleon mechanism which screens the deviation from standard gravity in dense environments, haloes more massive than $10^{13}h^{-1}M_\odot$ in this $f(R)$ model have very similar properties to haloes of similar mass in $\Lambda$CDM, while less massive haloes, such as that of the Milky Way, can have steeper inner density profiles and higher velocity dispersions due to their weaker screening. The halo concentration is remarkably enhanced for low-mass haloes in this model due to a deepening of the total gravitational potential. Contrary to the naive expectation, the halo formation time $z_f$ is later for low-mass haloes in this model, a consequence of these haloes growing faster than their counterparts in $\Lambda$CDM at late times and the definition of $z_f$. Subhaloes, especially those less massive than $10^{11}h^{-1}M_\odot$, are substantially more abundant in this $f(R)$ model for host haloes less massive than $10^{13}h^{-1}M_\odot$. We discuss the implications of these results for the Milky Way satellite abundance problem. Although the overall halo and subhalo properties in this borderline $f(R)$ model are close to their $\Lambda$CDM predictions, our results suggest that studies of the Local Group and astrophysical systems, aided by high-resolution simulations, can be valuable for further tests of it.

We simulate the phase-space distribution of stellar mass in 9 massive Lambda-CDM galaxy clusters by applying the semi-analytic particle tagging method of Cooper et al. to the Phoenix suite of high-resolution N-body simulations (M200 = 7.5 to 33 x 10^14 Msol). The resulting surface brightness (SB) profiles of brightest cluster galaxies (BCGs) match well to observations. On average, stars formed in galaxies accreted by the BCG account for ~90 per cent of its total mass (the remainder is formed in situ). In circular BCG-centred apertures, the superposition of multiple debris clouds (each ~10 per cent of the total BCG mass) from different progenitors can result in an extensive outer diffuse component, qualitatively similar to a 'cD envelope'. These clouds typically originate from tidal stripping at z < 1 and comprise both streams and the extended envelopes of other massive galaxies in the cluster. Stars at very low SB contribute a significant fraction of the total cluster stellar mass budget: in the central 1 Mpc^2 of a z ~ 0.15 cluster imaged at SDSS-like resolution, our fiducial model predicts 80-95 per cent of stellar mass below a SB of mu_V = 26.5 mag arcsec^2 is associated with accreted stars in the envelope of the BCG. The ratio of BCG stellar mass (including this diffuse component) to total cluster stellar mass is ~30 per cent.

We performed a hydrodynamical cosmological simulation of the formation of a Milky Way-like galaxy in a warm dark matter (WDM) cosmology. Smooth and dense filaments, several co-moving mega parsec long, form generically above z 2 in this model. Atomic line cooling allows gas in the centres of these filaments to cool to the base of the cooling function, resulting in a very striking pattern of extended Lyman-limit systems (LLSs). Observations of the correlation function of LLSs might hence provide useful limits on the nature of the dark matter. We argue that the self-shielding of filaments may lead to a thermal instability resulting in star formation. We implement a sub-grid model for this, and find that filaments rather than haloes dominate star formation until z 6. Reionisation decreases the gas density in filaments, and the more usual star formation in haloes dominates below z 6, although star formation in filaments continues until z=2. Fifteen per cent of the stars of the z=0 galaxy formed in filaments. At higher redshift, these stars give galaxies a stringy appearance, which, if observed, might be a strong indication that the dark matter is warm.

Discrepancies between the observed and model-predicted radio flux ratios are seen in a number of quadruply-lensed quasars. The most favored interpretation of these anomalies is that CDM substructures present in lensing galaxies perturb the lens potentials and alter image magnifications and thus flux ratios. So far no consensus has emerged regarding whether or not the predicted CDM substructure abundance fully accounts for the lensing flux anomaly observations. Accurate modeling relies on a realistic lens sample in terms of both the lens environment and internal structures and substructures. In this paper we construct samples of generalised and specific lens potentials, to which we add (rescaled) subhalo populations from the galaxy-scale Aquarius and the cluster-scale Phoenix simulation suites. We further investigate the lensing effects from subhalos of masses several orders of magnitude below the simulation resolution limit. The resulting flux ratio distributions are compared to the currently best available sample of radio lenses. The observed anomalies in B0128+437, B0712+472 and B1555+375 are more likely to be caused by propagation effects or oversimplified lens modeling, signs of which are already seen in the data. Among the quadruple systems that have closely located image triplets/pairs, the anomalous flux ratios of MG0414+0534 can be reproduced by adding CDM subhalos to its macroscopic lens potential, with a probability of 5%-20%; for B0712+472, B1422+231, B1555+375 and B2045+265, these probabilities are only of a few percent. We hence find that CDM substructures are unlikely to be the whole reason for radio flux anomalies. We discuss other possible effects that might also be at work.

Lensing flux-ratio anomalies are most likely caused by gravitational lensing by small-scale dark matter structures. These anomalies offer the prospect of testing a fundamental prediction of the cold dark matter (CDM) cosmological model: the existence of numerous substructures that are too small to host visible galaxies. In two previous studies we found that the number of subhalos in the six high-resolution simulations of CDM galactic halos of the Aquarius project is not sufficient to account for the observed frequency of flux ratio anomalies seen in selected quasars from the CLASS survey. These studies were limited by the small number of halos used, their narrow range of masses (1-2E12 solar masses) and the small range of lens ellipticities considered. We address these shortcomings by investigating the lensing properties of a large sample of halos with a wide range of masses in two sets of high resolution simulations of cosmological volumes and comparing them to a currently best available sample of radio quasars. We find that, as expected, substructures do not change the flux-ratio probability distribution of image pairs and triples with large separations, but they have a significant effect on the distribution at small separations. For such systems, CDM substructures can account for a substantial fraction of the observed flux-ratio anomalies. For large close-pair separation systems, the discrepancies existing between the observed flux ratios and predictions from smooth halo models are attributed to simplifications inherent in these models which do not take account of fine details in the lens mass distributions.

Recent work has suggested that the amplitude of the size mass relation of massive early type galaxies evolves with redshift. Here we use a semi-analytical galaxy formation model to study the size evolution of massive early type galaxies. We find this model is able to reproduce the amplitude of present day amplitude and slope of the relation between size and stellar mass for these galaxies, as well as its evolution. The amplitude of this relation reflects the typical compactness of dark halos at the time when most of the stars are formed. This link between size and star formation epoch is propagated in galaxy mergers. Mergers of high or moderate mass ratio (less than 1:3) become increasingly important with increasing present day stellar mass for galaxies more massive than $10^{11.4}M_{\odot}$. At lower masses, low mass ratio mergers play a more important role. In situ star formation contribute more to the size growth than it does to stellar mass growth. We also find that, for ETGs identified at $z=2$, minor mergers dominate subsequent growth both for stellar mass and in size, consistent with earlier theoretical results.

Well-motivated elementary particle candidates for the dark matter, such as the sterile neutrino, behave as warm dark matter (WDM).For particle masses of order a keV, free streaming produces a cutoff in the linear fluctuation power spectrum at a scale corresponding to dwarf galaxies. We investigate the abundance and structure of WDM haloes and subhaloes on these scales using high resolution cosmological N-body simulations of galactic haloes of mass similar to the Milky Way's. On scales larger than the free-streaming cutoff, the initial conditions have the same power spectrum and phases as one of the cold dark matter (CDM) haloes previously simulated by Springel et al as part of the Virgo consortium Aquarius project. We have simulated four haloes with WDM particle masses in the range 1.4-2.3keV and, for one case, we have carried out further simulations at varying resolution. N-body simulations in which the power spectrum cutoff is resolved are known to undergo artificial fragmentation in filaments producing spurious clumps which, for small masses (<10^7Msun in our case) outnumber genuine haloes. We have developed a robust algorithm to identify these spurious objects and remove them from our halo catalogues. We find that the WDM subhalo mass function is suppressed by well over an order magnitude relative to the CDM case for masses <10^9Msun. Requiring that there should be at least as many subhaloes as there are observed satellites in the Milky Way leads to a conservative lower limit to the (thermal equivalent) WDM particle mass of ~1.5\rmn{keV}. WDM haloes and subhaloes have cuspy density distributions that are well described by NFW or Einasto profiles. Their central densities are lower for lower WDM particle masses and none of the models we have considered suffer from the "too big to fail" problem recently highlighted by Boylan-Kolchin et al.

The emptiness of the Local Void has been put forward as a serious challenge to the current standard paradigm of structure formation in LCDM. We use a high resolution cosmological N-body simulation, the Millennium-II run, combined with a sophisticated semi-analytic galaxy formation model, to explore statistically whether the local void is allowed within our current knowledge of galaxy formation in LCDM. We find that about 15 percent of the Local Group analogue systems (11 of 77) in our simulation are associated with nearby low density regions having size and 'emptiness' similar to those of the observed Local Void. This suggests that, rather than a crisis of the LCDM, the emptiness of the Local Void is indeed a success of the standard LCDM theory. The paucity of faint galaxies in such voids results from a combination of two factors: a lower amplitude of the halo mass function in the voids than in the field, and a lower galaxy formation efficiency in void haloes due to halo assembly bias effects. While the former is the dominated factor, the later also plays a sizable role. The halo assembly bias effect results in a stellar mass fraction 25 percent lower for void galaxies when compared to field galaxies with the same halo mass.

Previous studies indicate that assembly bias effects are stronger for lower mass dark matter haloes. Here we make use of high resolution re-simulations of rich clusters and their surroundings from the Phoenix Project and a large volume cosmological simulation, the Millennium-II run, to quantify assembly bias effects on dwarf-sized dark matter haloes. We find that, in the regions around massive clusters, dwarf-sized haloes ($[10^9,10^{11}]\ms$) form earlier ($\Delta z \sim 2$ in redshift) and possess larger $V_{\rm max}$ ($\sim20%$) than the field galaxies. We find that this environmental dependence is largely caused by tidal interactions between the ejected haloes and their former hosts, while other large scale effects are less important. Finally we assess the effects of assembly bias on dwarf galaxy formation with a sophisticated semi-analytical galaxy formation model. We find that the dwarf galaxies near massive clusters tend to be redder ($\Delta(u-r) = 0.5$) and have three times as much stellar mass compared to the field galaxies with the same halo mass. These features should be seen with observational data.

In this paper we present a catalog which includes 141 bright candidates ($\leq10.27$ mag, V band) showing the infrared (IR) excess at 22 $\mu$m. Of which, 38 stars are known IR excess stars or disk, 23 stars are double or multiple stars and 4 are Be stars. While the remaining more than 70 stars are identified as the 22 $\mu$m excess candidates in our work. The criterion of selecting candidates is $K_s-[22]_{\mu m}$. All these candidates are selected from \emph{WISE} All-sky data cross-correlated with \emph{Hipparcos} Main Catalog and the likelihood-ratio technique is employed. Considering the effect of background, we introduce the \emph{IRAS} 100 $\mu$m level to exclude the high background. We also estimated the coincidence probability of these sources. In addition, we presented the optical to mid-infrared SEDs and optical images of all the candidates, and gave the observed optical spectra of 6 stars with NAOC's 2.16-m telescope. To measure for the dust amount around each star, the fractional luminosity is also provided. We also test whether our method of selecting IR excess stars can be used to search for extra-solar planets, we cross-matched our catalog with known IR-excess stars having planets but none is matched. Finally, we give the fraction of stars showing IR-excess for different spectral type of main-sequence stars.

We explore the cosmological halo-to-halo scatter of the distribution of mass within dark matter halos utilizing a well-resolved statistical sample of clusters from the cosmological Millennium simulation. We find that at any radius, the spherically-averaged dark matter density of a halo (corresponding to the "smooth-component") and its logarithmic slope are well-described by a Gaussian probability distribution. At small radii (within the scale radius), the density distribution is fully determined by the measured Gaussian distribution in halo concentrations. The variance in the radial distribution of mass in dark matter halos is important for the interpretation of direct and indirect dark matter detection efforts. The scatter in mass profiles imparts approximately a 25 percent cosmological uncertainty in the dark matter density at the Solar neighborhood and a factor of ~3 uncertainty in the expected Galactic dark matter annihilation flux. The aggregate effect of halo-to-halo profile scatter leads to a small (few percent) enhancement in dark matter annihilation background if the Gaussian concentration distribution holds for all halo masses versus a 10 percent enhancement under the assumption of a log-normal concentration distribution. The Gaussian nature of the cluster profile scatter implies that the technique of "stacking" halos to improve signal to noise should not suffer from bias.

We study the formation and evolution of Brightest Cluster Galaxies starting from a $z=2$ population of quiescent ellipticals and following them to $z=0$. To this end, we use a suite of nine high-resolution dark matter-only simulations of galaxy clusters in a $\Lambda$CDM universe. We develop a scheme in which simulation particles are weighted to generate realistic and dynamically stable stellar density profiles at $z=2$. Our initial conditions assign a stellar mass to every identified dark halo as expected from abundance matching; assuming there exists a one-to-one relation between the visible properties of galaxies and their host haloes. We set the sizes of the luminous components according to the observed relations for $z\sim2$ massive quiescent galaxies. We study the evolution of the mass-size relation, the fate of satellite galaxies and the mass aggregation of the cluster central. From $z=2$, these galaxies grow on average in size by a factor 5 to 10 of and in mass by 2 to 3. The stellar mass growth rate of the simulated BCGs in our sample is of 1.9 in the range $0.3<z<1.0$ consistent with observations, and of 1.5 in the range $0.0<z<0.3$. Furthermore the satellite galaxies evolve to the present day mass-size relation by $z=0$. Assuming passively evolving stellar populations, we present surface brightness profiles for our cluster centrals which resemble those observed for the cDs in similar mass clusters both at $z=0$ and at $z=1$. This demonstrates that the $\Lambda$CDM cosmology does indeed predict minor and major mergers to occur in galaxy clusters with the frequency and mass ratio distribution required to explain the observed growth in size of passive galaxies since $z=2$. Our experiment shows that Brightest Cluster Galaxies can form through dissipationless mergers of quiescent massive $z=2$ galaxies, without substantial additional star formation.

We have performed a series of numerical experiments to investigate how the primordial thermal velocities of fermionic dark matter particles affect the physical and phase space density profiles of the dark matter haloes into which they collect. The initial particle velocities induce central cores in both profiles, which can be understood in the framework of phase space density theory. We find that the maximum coarse-grained phase space density of the simulated haloes (computed in 6 dimensional phase space using the EnBid code) is very close to the theoretical fine-grained upper bound, while the pseudo phase space density, Q ~ {\rho}/{\sigma}^3, overestimates the maximum phase space density by up to an order of magnitude. The density in the inner regions of the simulated haloes is well described by a 'pseudo-isothermal' profile with a core. We have developed a simple model based on this profile which, given the observed surface brightness profile of a galaxy and its central velocity dispersion, accurately predicts its central phase space density. Applying this model to the dwarf spheroidal satellites of the Milky Way yields values close to 0.5 keV for the mass of a hypothetical thermal warm dark matter particle, assuming the satellite haloes have cores produced by warm dark matter free streaming. Such a small value is in conflict with the lower limit of 1.2 keV set by observations of the Lyman-{\alpha} forest. Thus, if the Milky Way dwarf spheroidal satellites have cores, these are likely due to baryonic processes associated with the forming galaxy, perhaps of the kind proposed by Navarro, Eke and Frenk and seen in recent simulations of galaxy formation in the cold dark matter model.

Cold dark matter models predict the existence of a large number of substructures within dark matter halos. If the cold dark matter consists of weakly interacting massive particles, their annihilation within these substructures could lead to diffuse GeV emission that would dominate over the annihilation signal of the host halo. In this work we search for GeV emission from three nearby galaxy clusters: Coma, Virgo and Fornax. We first remove known extragalactic and galactic diffuse gamma-ray backgrounds and point sources from the Fermi 2-year catalog and find a significant residual diffuse emission in all three clusters. We then investigate whether this emission is due to (i) unresolved point sources; (ii) dark matter annihilation; or (iii) cosmic rays (CR). Using 45 months of Fermi-LAT data we detect several new point sources (not present in the Fermi 2-year point source catalogue) which contaminate the signal previously analyzed by Han et al.(arxiv:1201.1003). Including these and accounting for the effects of undetected point sources, we find no significant detection of extended emission from the three clusters studied. Instead, we determine upper limits on emission due to dark matter annihilation and cosmic rays. For Fornax and Virgo the limits on CR emission are consistent with theoretical models, but for Coma the upper limit is a factor of 2 below the theoretical expectation. Allowing for systematic uncertainties associated with the treatment of CR, the upper limits on the cross section for dark matter annihilation from our clusters are more stringent than those from analyses of dwarf galaxies in the Milky Way. We rule out the thermal cross section for supersymmetric dark matter particles for masses as large as 100 GeV (depending on the annihilation channel).

Analysis of the Pangu N-body simulation validates that the bulk flow of halos follows a Maxwellian distribution which variance is consistent with the prediction of the linear theory of structure formation. We propose that the consistency between the observed bulk velocity and theories should be examined at the effective scale of the radius of a spherical top-hat window function yielding the same smoothed velocity variance in linear theory as the sample window function does. We compared some recently estimated bulk flows from observational samples with the prediction of the \Lambda CDM model we used; some results deviate from expectation at a level of ~ 3\sigma but the discrepancy is not as severe as previously claimed. We show that bulk flow is only weakly correlated with the dipole of the internal mass distribution, the alignment angle between the mass dipole and the bulk flow has a broad distribution peaked at ~ 30-50 deg., and also that the bulk flow shows little dependence on the mass of the halos used in the estimation. In a simulation of box size 1Gpc/h, for a cell of radius 100 Mpc/h the maximal bulk velocity is >500 km/s, dipoles of the environmental mass outside the cell are not tightly aligned with the bulk flow, but are rather located randomly around it with separation angles ~ 20-40 deg. In the fastest cell there is a slightly smaller number of low-mass halos; however halos inside are clustered more strongly at scales > ~ 20 Mpc/h, which might be a significant feature since the correlation between bulk flow and halo clustering actually increases in significance beyond such scales.

The weakly interacting massive particle (WIMP) often serves as a candidate for the cold dark matter, however when produced non-thermally it could behave like warm dark matter. In this paper we study the properties of the $\gamma$-ray emission from annihilation of WIMP dark matter in the halo of our own Milky-Way Galaxy with high resolution $N$-body simulations of a Milky-Way like dark matter halo, assuming different nature of WIMPs. Due to the large free-streaming length in the scenario of warm WIMPs, the substructure contend of the dark matter halo is significantly different from that of the cold WIMP counterpart, resulting in distinct predictions of the $\gamma$-ray signals from the dark matter annihilation. We illustrate these by comparing the predicted $\gamma$-ray signals from the warm WIMP annihilation to that of cold WIMPs. Pronounced differences from the subhalo skymap and statistical properties between two WIMP models are demonstrated. Due to the potentially enhanced cross section of the non-thermal production mechanism in warm WIMP scenario, the Galactic center might be prior for the indirect detection of warm WIMPs to dwarf galaxies, which might be different from the cold dark matter scenario. As a specific example we consider the non-thermally produced neutralino of supersymmetric model and discuss the detectability of warm WIMPs with Fermi $\gamma$-ray telescope.

We report evidence for extended gamma-ray emission from the Virgo, Fornax and Coma clusters based on a maximum-likelihood analysis of the 3-year Fermi-LAT data. For all three clusters, excess emission is observed within three degrees of the center, peaking at the GeV scale. This emission cannot be accounted for by known Fermi sources or by the galactic and extragalactic backgrounds. If interpreted as annihilation emission from supersymmetric dark matter (DM) particles, the data prefer models with a particle mass in the range 20-60 GeV annihilating into the b-bbar channel, or 2-10 GeV and >1 TeV annihilating into mu-mu final states. Our results are consistent with those obtained by Hooper and Linden from a recent analysis of Fermi-LAT data in the region of the Galactic Centre. An extended DM annihilation profile dominated by emission from substructures is preferred over a simple point source model. The significance of DM detection is 4.4 sigma in Virgo and lower in the other two clusters. We also consider the possibility that the excess emission arises from cosmic ray (CR) induced gamma-rays, and infer a CR level within a factor of three of that expected from analytical models. However, the significance of a CR component is lower than the significance of a DM component, and there is no need for such a CR component in the presence of a DM component in the preferred DM mass range. We also set flux and cross-section upper limits for DM annihilation into the b-bbar and mu-mu channels in all three clusters.

[abridged] We introduce the Phoenix Project, a set of $\Lambda$CDM simulations of the dark matter component of nine rich galaxy clusters. Each cluster is simulated at least at two different numerical resolutions. For eight of them, the highest resolution corresponds to $\sim 130$ million particles within the virial radius, while for one this number is over one billion. We study the structure and substructure of these systems and contrast them with six galaxy-sized dark matter haloes from the Aquarius Project, simulated at comparable resolution. This comparison highlights the approximate mass invariance of CDM halo structure and substructure. We find little difference in the spherically-averaged mass, pseudo-phase-space density, and velocity anisotropy profiles of Aquarius and Phoenix haloes. When scaled to the virial properties of the host halo, the abundance and radial distribution of subhaloes are also very similar, despite the fact that Aquarius and Phoenix haloes differ by roughly three decades in virial mass. The most notable difference is that cluster haloes have been assembled more recently and are thus significantly less relaxed than galaxy haloes, which leads to decreased regularity, increased halo-to-halo scatter and sizable deviations from the mean trends. This accentuates the effects of the strong asphericity of individual clusters on surface density profiles, which may vary by up to a factor of three at a given radius, depending on projection. The high apparent concentration reported for some strong-lensing clusters might very well reflect these effects. A more recent assembly also explains why substructure in some Phoenix haloes is slightly more abundant than in Aquarius, especially in the inner regions. Resolved subhaloes nevertheless contribute only $11 \pm 3%$ of the virial mass in Phoenix clusters. .

Recent studies suggest that only three of the twelve brightest satellites of the Milky Way (MW) inhabit dark matter halos with maximum circular velocity, V_max, exceeding 30km/s. This is in apparent contradiction with the LCDM simulations of the Aquarius Project, which suggest that MW-sized halos should have at least 8 subhalos with V_max>30km/s. The absence of luminous satellites in such massive subhalos is thus puzzling and may present a challenge to the LCDM paradigm. We note, however, that the number of massive subhalos depends sensitively on the (poorly-known) virial mass of the Milky Way, and that their scarcity makes estimates of their abundance from a small simulation set like Aquarius uncertain. We use the Millennium Simulation series and the invariance of the scaled subhalo velocity function (i.e., the number of subhalos as a function of \nu, the ratio of subhalo V_max to host halo virial velocity, V_200) to secure improved estimates of the abundance of rare massive subsystems. In the range 0.1<\nu<0.5, N_sub(>\nu) is approximately Poisson-distributed about an average given by <N_sub>=10.2x(\nu/0.15)^(-3.11). This is slightly lower than in Aquarius halos, but consistent with recent results from the Phoenix Project. The probability that a LCDM halo has 3 or fewer subhalos with V_max above some threshold value, V_th, is then straightforward to compute. It decreases steeply both with decreasing V_th and with increasing halo mass. For V_th=30km/s, ~40% of M_halo=10^12 M_sun halos pass the test; fewer than 5% do so for M_halo>= 2x10^12 M_sun; and the probability effectively vanishes for M_halo>= 3x 10^12 M_sun. Rather than a failure of LCDM, the absence of massive subhalos might simply indicate that the Milky Way is less massive than is commonly thought.

The flux-ratio anomalies observed in multiply-lensed quasar images are most plausibly explained as the result of perturbing structures superposed on the underlying smooth matter distribution of the primary lens. The cold dark matter cosmological model predicts that a large number of substructures should survive inside larger halos but, surprisingly, this population alone has been shown to be insufficient to explain the observed distribution of the flux ratios of quasar's multiple images. Other halos (and their own subhalos) projected along the line of sight to the primary lens have been considered as additional source of perturbation. In this work, we use ray tracing through the Millennium II simulation to investigate the importance of projection effects due to halos and subhalos of mass m>1E8 Msun/h and extend our analysis to lower masses, m>1E6 Msun/h, using Monte-Carlo halo distributions. We find that the magnitude of the violation depends strongly on the density profile and concentration of the intervening halos, but clustering plays only a minor role. For a typical lensing geometry (lens at redshift 0.6 and source at redshift 2), background haloes (behind the main lens) are more likely to cause a violation than foreground halos. We conclude that line-of-sight structures can be as important as intrinsic substructures in causing flux-ratio anomalies. The combined effect of perturbing structures within the lens and along the line of sight in the LCDM universe results in a cusp-violation probability of 20-30%. This alleviates the discrepancy between models and current data, but a larger observational sample is required for a stronger test of the theory.

High resolution N-body simulations of galactic cold dark matter haloes indicate that we should expect to find a few satellite galaxies around the Milky Way whose haloes have a maximum circular velocity in excess of 40 kms. Yet, with the exception of the Magellanic Clouds and the Sagittarius dwarf, which likely reside in subhaloes with significantly larger velocities than this, the bright satellites of the Milky Way all appear to reside in subhaloes with maximum circular velocities below 40 kms. As recently highlighted by Boylan-Kolchin et al., this discrepancy implies that the majority of the most massive subhaloes within a cold dark matter galactic halo are much too concentrated to be consistent with the kinematic data for the bright Milky Way satellites. Here we show that no such discrepancy exists if haloes are made of warm, rather than cold dark matter because these haloes are less concentrated on account of their typically later formation epochs. Warm dark matter is one of several possible explanations for the observed kinematics of the satellites.

Recent observations indicate that many if not all galaxies host massive central black holes (BHs). In this paper we explore the influence of supermassive binary black holes (SMBBHs) on their actions as gravitational lenses. When lenses are modelled as singular isothermal ellipsoids, binary black holes change the critical curves and caustics differently as a function of distance. Each black hole can in principle create at least one additional image, which, if observed, provides evidence of black holes. By studying how SMBBHs affect the cumulative distribution of magnification for images created by black holes, we find that the cross section for at least one such additional image to have a magnification larger than $10^{-5}$ is comparable to the cross section for producing multiple-images in singular isothermal lenses. Such additional images may be detectable with high-resolution and large dynamic range maps of multiply-imaged systems from future facilities, such as the SKA. The probability of detecting at least one image (two images) with magnification above $10^{-3}$ is $\sim 0.2 \fBH$ ($\sim 0.05 \fBH$) in a multiply-imaged lens system, where $\fBH$ is the fraction of galaxies housing binary black holes. We also study the effects of SMBBHs on the core images when galaxies have shallower central density profiles (modelled as non-singular isothermal ellipsoids). We find that the cross section of the usually faint core images is further suppressed by SMBBHs. Thus their presence should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.

If the dark matter consists of supersymmetric particles, Gamma-ray observatories such as the Large Area Telescope aboard the Fermi satellite may detect annihilation radiation from the haloes of galaxies and galaxy clusters. Much recent effort has been devoted to searching for this signal around the Milky Way's dwarf satellites. Using a new suite of high-resolution simulations of galaxy cluster haloes (the Phoenix Project), together with the Aquarius simulations of Milky-Way-like galaxy haloes, we show that higher signal-to-noise and equally clean signals are, in fact, predicted to come from nearby rich galaxy clusters. Most of the cluster emission is produced by small subhaloes with masses less than that of the Sun. The large range of mass scales covered by our two sets of simulations allows us to deduce a physically motivated extrapolation to these small (and unresolved) masses. Since tidal effects destroy subhaloes in the dense inner regions of haloes, most cluster emission is then predicted to come from large radii, implying that the nearest and brightest systems should be much more extended than Fermi's angular resolution limit. The most promising targets for detection are clusters such as Coma and Fornax, but detection algorithms must be tuned to the predicted profile of the emission if they are to maximize the chance of finding this weak signal.

We study the population statistics of the surviving subhaloes of LCDM dark matter haloes using a set of very high resolution N-body simulations. These include both simulations of representative regions of the Universe and ultra-high resolution resimulations of individual dark matter haloes. We find that more massive haloes tend to have a larger mass fraction in subhaloes. For example, cluster size haloes typically have 7.5 percent of their mass within R200 in substructures of fractional mass larger than 1e-5, which is 25 percent higher than galactic haloes. There is, however, a large variance in the subhalo mass fraction from halo to halo, whereas the subhalo abundance shows much higher regularity. For dark matter haloes of fixed mass, the subhalo abundance decreases by 30 percent between redshift 2 and 0. The subhalo abundance function correlates with the host halo concentration parameter and formation redshift. However, the intrinsic scatter is not significantly reduced for narrow ranges of concentration parameter or formation redshift, showing that they are not the dominant parameters that determine the subhalo abundance in a halo.

We study the dependence of the galaxy contents within halos on the halo formation time using two galaxy formation models, one being a semianalytic model utilizing the halo assembly history from a high resolution N-body simulation and the other being a smoothed particle hydrodynamics simulation including radiative cooling, star formation, and energy feedback from galactic winds. We confirm the finding by Gao et al. that at fixed mass, the clustering of halos depends on the halo formation time, especially for low-mass halos. This age dependence of halo clustering makes it desirable to study the correlation between the occupation of galaxies within halos and the halo age. We find that, in halos of fixed mass, the number of satellite galaxies has a strong dependence on halo age, with fewer satellites in older halos. The youngest one-third of the halos can have an order of magnitude more satellites than the oldest one-third. For central galaxies, in halos that form earlier, they tend to have more stars and thus appear to be more luminous, and the dependence of their luminosity on halo age is not as strong as that of stellar mass. The results can be understood through the star formation history in halos and the merging of satellites onto central galaxies. The age dependence of the galaxy contents within halos would constitute an important ingredient in a more accurate halo-based model of galaxy clustering.

Lensing flux-ratio anomalies have been frequently observed and taken as evidence for the presence of abundant dark matter substructures in lensing galaxies, as predicted by the cold dark matter (CDM) model of cosmogony. In previous work, we examined the cusp-caustic relations of the multiple images of background quasars lensed by galaxy-scale dark matter haloes, using a suite of high-resolution N-body simulations (the Aquarius simulations). In this work, we extend our previous calculations to incorporate both the baryonic and diffuse dark components in lensing haloes. We include in each lensing simulation: (1) a satellite galaxy population derived from a semi-analytic model applied to the Aquarius haloes, (2) an empirical Milky-Way globular cluster population and (3) satellite streams (diffuse dark component) identified in the simulations. Accounting for these extra components, we confirm our earlier conclusion that the abundance of intrinsic substructures (dark or bright, bound or diffuse) is not sufficient to explain the observed frequency of cusp-caustic violations in the CLASS survey. We conclude that the observed effect could be the result of the small number statistics of CLASS, or intergalactic haloes along the line of sight acting as additional sources of lensing flux anomalies. Another possibility is that this discrepancy signals a failure of the CDM model.

We have implemented a simple model to identify the likely sites of the first stars and galaxies in the high-resolution simulations of the formation of galactic dark matter halos of the Aquarius Project. The first star in a galaxy like the Milky Way formed around redshift z=35; by z=10, the young galaxy contained up to 30000 dark matter haloes capable of forming stars by molecular hydrogen cooling. These minihaloes were strongly clustered and feedback may have severely limited the actual number of Population III stars that formed. By the present day, the remnants of the first stars would be strongly concentrated to the centre of the main halo. If a second generation of long-lived stars formed near the first (the first star relics), we would expect to find half of them within 30kpc/h of the galactic centre and a significant fraction in satellites where they may be singled out by their anomalous metallicity patterns. The first halo in which gas could cool by atomic hydrogen line radiation formed at z=25; by z=10, the number of such first galaxies had increased to ~300. Feedback might have decreased the number of first galaxies at the time they undergo significant star formation, but not the number that survive to the present because near neighbours merge. Half of all the "first galaxies" that form before z=10 merge with the main halo before z~3 and most lose a significant fraction of their mass. However, today there should still be more than 20 remnants orbiting within the central ~30kpc/h of the Milky Way. These satellites have circular velocities of a few kilometers per second or more, comparable to those of known Milky Way dwarfs. They are a promising hunting ground for the remnants of the earliest epoch of star formation.

We use high-resolution Aquarius simulations of Milky Way-sized haloes in the LCDM cosmology to study the effects of dark matter substructures on gravitational lensing. Each halo is resolved with ~ 10^8 particles (at a mass resolution ~ 10^3-4 M_sun/h) within its virial radius. Subhaloes with masses larger than 10^5 M_sun/h are well resolved, an improvement of at least two orders of magnitude over previous lensing studies. We incorporate a baryonic component modelled as a Hernquist profile and account for the response of the dark matter via adiabatic contraction. We focus on the "anomalous" flux ratio problem, in particular on the violation of the cusp-caustic relation due to substructures. We find that subhaloes with masses less than ~ 10^8 M_sun/h play an important role in causing flux anomalies; such low mass subhaloes have been unresolved in previous studies. There is large scatter in the predicted flux ratios between different haloes and between different projections of the same halo. In some cases, the frequency of predicted anomalous flux ratios is comparable to that observed for the radio lenses, although in most cases it is not. The probability for the simulations to reproduce the observed violations of the cusp lenses is about 0.001. We therefore conclude that the amount of substructure in the central regions of the Aquarius haloes is insufficient to explain the observed frequency of violations of the cusp-caustic relation. These conclusions are based purely on our dark matter simulations which ignore the effect of baryons on subhalo survivability.

We perform cosmological hydrodynamic simulations to determine to what extent galaxies lose their gas due to photoheating from an ionizing background. We find that the characteristic mass at which haloes on average have lost half of their baryons is Mc ~ 6.5 x 10^9 Msun/h at z = 0, which corresponds to a circular velocity of 25 km/s. This is significantly lower than the filtering mass obtained by the linear theory, which is often used in semianalytical models of galaxy formation. We demonstrate it is the gas temperature at the virial radius which determines whether a halo can accrete gas. A simple model that follows the merger history of the dark matter progenitors, and where gas accretion is not allowed when this temperature is higher than the virial temperature of the halo, reproduces the results from the simulation remarkably well. This model can be applied to any reionization history, and is easy to incorporate in semianalytical models.

(Abridged) We use two very large cosmological simulations to study how the density profiles of relaxed LCDM dark halos depend on redshift and on halo mass. We confirm that these profiles deviate slightly but systematically from the NFW form and are better approximated by the empirical formula, $d\log \rho/d\log r \propto r^{\alpha}$, first used by Einasto to fit star counts in the Milky Way. The best-fit value of the additional shape parameter, $\alpha$, increases gradually with mass, from alpha ~ 0.16 for present-day galaxy halos to alpha~ 0.3 for the rarest and most massive clusters. Halo concentrations depend only weakly on mass at z=0, and this dependence weakens further at earlier times. At z ~3 the average concentration of relaxed halos does not vary appreciably over the mass range accessible to our simulations (M >=3e11Msun). Furthermore, in our biggest simulation, the average concentration of the most massive}, relaxed halos is constant at c200 ~3.5 to 4 for 0<=z<= 3. These results have important implications for currently planned surveys of distant clusters.

In this paper we use the ``Millennium Simulation'' to re-examine the mass assembly history of dark matter halos and the age dependence of halo clustering. We use eight different definitions of halo formation times to characterize the different aspects of the assembly history of a dark matter halo. We find that these formation times have different dependence on halo mass. While some formation times characterize well the hierarchical nature of halo formation, in the sense that more massive halos have later formation, the trend is reversed for other definitions of the formation time. In particular, the formation times that are likely to be related to the formation of galaxies in dark halos show strong trends of ``down-sizing'', in that lower-mass halos form later. We also investigate how the correlation amplitude of dark matter halos depends on the different formation times. We find that this dependence is quite strong for some definitions of formation time but weak or absent for other definitions. In particular, the correlation amplitude of halos of a given mass is almost independent of their last major merger time. For the definitions that are expected to be more related to the formation of galaxies in dark halos, a significant assembly bias is found only for halos less massive than M_*. We discuss our results in connection to the hierarchical assembly of dark matter halos, the ``archaeological down-sizing'' observed in the galaxy population, and the observed color-dependence of the clustering strength of galaxy groups and clusters.

The first stars in the Universe form when chemically pristine gas heats as it falls into dark matter potential wells, cools radiatively due to the formation of molecular hydrogen, and becomes self-gravitating. We demonstrate with super-computer simulations that their properties depend critically on the currently unknown nature of the dark matter. If the dark matter particles have intrinsic velocities that wipe-out small-scale structure, then the first stars form in filaments with lengths of order the free-streaming scale, which can be about 10^20m (~3kpc, baryonic masses 10^7 solar masses) for realistic "warm dark matter" candidates. Fragmentation of the filaments forms stars with a range of masses which may explain the observed peculiar element abundance pattern of extremely metal-poor stars, while coalescence of fragments and stars during the filament's ultimate collapse may seed the super massive black holes that lurk in the centres of most massive galaxies.

We present a precise estimate of the bulk virial scaling relation of halos formed via hierarchical clustering in an ensemble of simulated cold dark matter cosmologies. The result is insensitive to cosmological parameters, the presence of a trace, dissipationless gas component, and numerical resolution down to a limit of ~1000 particles. The dark matter velocity dispersion scales with total mass as log(sigma_{DM}(M,z)) = log(1082.9 +- 4.0 \kms) + (0.3361 +- 0.0026) log(h(z)M_{200}/10^{15} Msun), with h(z) the dimensionless Hubble parameter. At fixed mass, the velocity dispersion likelihood is nearly log-normal, with scatter sigma_{ln sigma} = 0.0426 +- 0.015, except for a tail to higher dispersions containing 10% of the population that are merger transients. We combine this relation with the halo mass function in LCDM models, and show that a low normalization condition, S_8 = sigma_8(Omega_m/0.3)^{0.35} \= 0.69, favored by recent WMAP and SDSS analysis requires that galaxy and gas specific energies in rich clusters be 50% larger than that of the underlying dark matter. Such large energetic biases are in conflict with the current generation of direct simulations of cluster formation. A higher normalization, S_8 = 0.80, alleviates this tension and implies that the hot gas fraction within r_{500} is (0.71 +- 0.09) h_{70}^{-3/2} Omega_b/\Omega_m, a value consistent with recent Sunyaev-Zel'dovich observations.

We investigate the nature of CO emission from z~6 quasars by combining non-LTE radiative transfer calculations with merger-driven models of z~6 quasar formation that arise naturally in LCDM cosmological simulations. We consider four model quasars formed in 10^12-10^13 M_sun halos from different merging histories. Our main results follow. Owing to massive starbursts and funneling of dense gas into the nuclear regions of merging galaxies, the CO is highly excited and the flux density peaks between J=5-8. The CO morphology of z~6 quasars often exhibits multiple emission peaks which arise from H2 concentrations which have not yet fully coalesced. Quasars at z~6 display a large range of sightline dependent line widths such that the lines are narrowest when the rotating H2 gas associated with the quasar is viewed face-on (when L_B is largest), and broadest when the gas is seen edge-on (when L_B is lowest). Thus for all models selection effects exist such that quasars selected for optical luminosity are preferentially face-on which may result in detected CO line widths narrower than the median. The sightline averaged line width is reflective of the circular velocity (V_c) of the host halo, and ranges from sigma~300-650 km/s. For optically selected QSOs, 10-25% (halo-mass dependant) of sightlines have narrow line widths compatible with the sole CO detection at z>6, J1148+5251. When accounting for both the temporal evolution of CO line widths, as well as the redshift evolution of halo V_c, these models self-consistently account for the CO line widths of both z~2 sub-mm galaxies and QSO's. Finally, the dynamical mass derived from the sightline averaged line widths provides a good estimate of the total mass, and allows for a stellar bulge and SMBH consistent with the local M_BH-M_bulge relation. [abridged]

The discovery of luminous quasars at redshift z ~ 6 indicates the presence of supermassive black holes (SMBHs) of mass ~10^9 Msun when the Universe was less than one billion years old. This finding presents several challenges for theoretical models. Here, we present the first multi-scale simulations that, together with a self-regulated model for the SMBH growth, produce a luminous quasar at z ~ 6.5 in the LCDM paradigm. We follow the hierarchical assembly history of the most massive halo in a ~ 3 Gpc^3 volume, and find that this halo of ~ 8x 10^{12} Msun forming at z ~ 6.5 after several major mergers is able to reproduce a number of observed properties of SDSS J1148+5251, the most distant quasar detected at z =6.42 (Fan et al. 2003). Moreover, the SMBHs grow through gas accretion below the Eddington limit in a self-regulated manner owing to feedback. We find that the progenitors experience significant star formation (up to 10^4 Msun/yr) preceding the major quasar phase such that the stellar mass of the quasar host reaches 10^{12} Msun at z ~ 6.5, consistent with observations of significant metal enrichment in SDSS J1148+5251. Our results provide a viable formation mechanism for z ~ 6 quasars in the standard LCDM cosmology, and demonstrate a common, merger-driven origin for the rarest quasars and the fundamental SMBH-host correlation in a hierarchical Universe.(Abridged)

We use the Millennium Simulation (MS) to study the statistics of LCDM halo concentrations at z = 0. Our results confirm that the average halo concentration declines monotonically with mass; a power-law fits well the concentration-mass relation for over 3 decades in mass, up to the most massive objects to form in a LCDM universe (~ 10^15 h^-1 Msol). This is in clear disagreement with the predictions of the model proposed by Bullock et al. for these rare objects, and agrees better with the original predictions of Navarro, Frenk, & White. The large volume surveyed, together with the unprecedented numerical resolution of the MS, allow us to estimate with confidence the distribution of concentrations and, consequently, the abundance of systems with unusual properties. About one in a hundred cluster haloes (M200 >~ 3x10^14 h^-1 Msol) have concentrations exceeding c200 = 7.5, a result that may be used to interpret the likelihood of unusually strong massive gravitational lenses, such as Abell 1689, in the LCDM cosmogony. A similar fraction (1 in 100) of galaxy-sized haloes (M200 ~ 10^12 h^-1 Msol) have c200 < 4.5, an important constraint on models that attempt to reconcile the rotation curves of low surface-brightness galaxies by appealing to haloes of unexpectedly low concentration. We find that halo concentrations are independent of spin once haloes manifestly out of equilibrium are removed from the sample. Compared to their relaxed brethren, the concentrations of out-of-equilibrium haloes tend to be lower and to have more scatter, while their spins tend to be higher... (continue)

We have performed a large set of high-resolution cosmological simulations using smoothed particle hydrodynamics to study the formation of the first luminous objects in the LCDM cosmology. We follow the collapse of primordial gas clouds in eight early structures and document the scatter in the properties of the first star-forming clouds. Our first objects span formation redshifts from z~10 to z~50 and cover an order of magnitude in halo mas. We find that the physical properties of the central star-forming clouds are very similar in all of the simulated objects despite significant differences in formation redshift and environment. The physical properties of the clouds have little correlation with spin, mass, or assembly history of the host halo. The collapse of protostellar objects at higher redshifts progresses much more rapidly due to the higher densities, which accelerates the formation of molecular hydrogen, enhances initial cooling and shortens the dynamical timescales. The mass of the star-forming clouds cover a broad range, from a few hundred to a few thousand solar masses, and exhibit various morphologies: some of have disk-like structures nearly rotational supported; others form flattened spheroids; still others form bars. All of them develop a single proto-stellar `seed' which does not fragment into multiple objects up to the moment that the central gas becomes optically thick to H2 cooling lines. At this time, the instantaneous mass accretion rate onto the centre varies significantly from object to object, with disk-like structuress have the smallest mass accretion rates. The formation epoch and properties of the star-forming clouds are sensitive to the values of cosmological parameters.

We use a very large simulation of structure growth in a LCDM universe -- the Millennium Simulation -- to study assembly bias, the fact that the large-scale clustering of haloes of given mass varies significantly with their assembly history. We extend earlier work based on the same simulation by superposing results for redshifts from 0 to 3, by defining a less noisy estimator of clustering amplitude, and by considering halo concentration, substructure mass fraction and spin, as well as formation time, as additional parameters. These improvements lead to results with less noise than previous studies and covering a wider range of halo masses and structural properties. We find significant and significantly different assembly bias effects for all the halo properties we consider, although in all cases the dependences on halo mass and on redshift are adequately described as a dependence on equivalent peak height nu(M,z). The nu-dependences for different halo properties differ qualitatively and are not related as might naively be expected given the relations between formation time, concentration, substructure fraction and spin found for the halo population as a whole. These results suggest that it will be difficult to build models for the galaxy populations of dark haloes which can robustly relate the amplitude of large-scale galaxy clustering to that for mass clustering at better than the 10% level.

The clustering of dark halos depends not only on their mass but also on their assembly history, a dependence we term `assembly bias'. Using a galaxy formation model grafted onto the Millennium Simulation of the LCDM cosmogony, we study how assembly bias affects galaxy clustering. We compare the original simulation to `shuffled' versions where the galaxy populations are randomly swapped among halos of similar mass, thus isolating the effects of correlations between assembly history and environment at fixed mass. Such correlations are ignored in the halo occupation distribution models often used populate dark matter simulations with galaxies, but they are significant in our more realistic simulation. Assembly bias enhances 2-point correlations by 10% for galaxies with M_bJ-5logh brighter than -17, but suppresses them by a similar amount for galaxies brighter than -20. When such samples are split by colour, assembly bias is 5% stronger for red galaxies and 5% weaker for blue ones. Halo central galaxies are differently affected by assembly bias than are galaxies of all types. It almost doubles the correlation amplitude for faint red central galaxies. Shuffling galaxies among halos of fixed formation redshift or concentration in addition to fixed mass produces biases which are not much smaller than when mass alone is fixed. Assembly bias must reflect a correlation of environment with aspects of halo assembly which are not encoded in either of these parameters. It induces effects which could compromise precision measurements of cosmological parameters from large galaxy surveys.

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.

Many galaxies display warps, lopsided images, asymmetric rotation curves or other features which suggest that their immediate dynamical environment is neither static nor in equilibrium. In Cold Dark Matter (CDM) theories, such non-equilibrium features are expected in the inner regions of many dark haloes as a result of recent hierarchical growth. We used the excellent statistics provided by the very large Millennium Simulation to study (i) how the distribution of position and velocity asymmetries predicted for halo cores by the concordance LCDM cosmogony depends on halo mass, and (ii) how much of the dark matter in the inner core has been added at relatively recent times. Asymmetries are typically larger in more massive haloes. Thus 20% of cluster halos have density centre separated from barycentre by more than 20% of the virial radius, while only 7% of Milky Way halos have such large asymmetries. About 40% of all cluster halos have a mean core velocity which differs from the barycentre velocity by more than a quarter of the characteristic halo circular velocity, whereas only 10% of Milky Way halos have such large velocity offsets. About 25% of all cluster haloes have acquired more than a quarter of the mass currently in their inner 10kpc through mergers since z=1. The corresponding percentage of Milky Way haloes is 15%. These numbers seem quite compatible with the levels of asymmetry seen in the observable regions of galaxies, but quantitative comparison requires more detailed modelling of the observable components.

We use high-resolution $N$-body simulations to investigate the optical depth of giant arcs with length-to-width ratio larger than 7.5 and 10 in the `standard' $\LCDM$ model with $\sigma_8=0.9$ and $\Omega_{\rm m,0}=0.3$ and a model based on three-year Wilkinson Microwave Anisotropy Probe (WMAP) data. We find that, in dark-matter only simulations, the lensing probability in the three-year WMAP model (with $\sigma_8=0.74$ and $\Omega_{\rm m,0}=0.238)$ decreases by a factor of $\sim 6$ compared with that in the `standard' $\LCDM$ model. The effects of baryonic cooling, star formation and feedbacks are uncertain, but we argue that baryons will only increase the the lensing cross-section by a moderate factor, $\sim 2$. We conclude that the low central value of $\sigma_8$ and $\Omega_{\rm m,0}$ preferred by the WMAP three-year data may be too low to be compatible with observations if conventional assumptions of the background source population are correct.

We simulate the growth of galaxies and their central supermassive black holes by implementing a suite of semi-analytic models on the output of the Millennium Run, a very large simulation of the concordance LCDM cosmogony. Our procedures follow the detailed assembly history of each object and are able to track the evolution of all galaxies more massive than the Small Magellanic Cloud throughout a volume comparable to that of large modern redshift surveys. In this first paper we supplement previous treatments of the growth and activity of central black holes with a new model for `radio' feedback from those AGN that lie at the centre of a quasistatic X-ray emitting atmosphere in a galaxy group or cluster. We show that for energetically and observationally plausible parameters such a model can simultaneously explain: (i) the low observed mass drop-out rate in cooling flows; (ii) the exponential cut-off at the bright end of the galaxy luminosity function; and (iii) the fact that the most massive galaxies tend to be bulge-dominated systems in clusters and to contain systematically older stars than lower mass galaxies. This success occurs because static hot atmospheres form only in the most massive structures, and radio feedback (in contrast, for example, to supernova or starburst feedback) can suppress further cooling and star formation without itself requiring star formation. We discuss possible physical models which might explain the accretion rate scalings required for our phenomenological `radio mode' model to be successful.

Future weak lensing measurements of cosmic shear will reach such high accuracy that second order effects in weak lensing modeling, like the influence of baryons on structure formation, become important. We use a controlled set of state of the art cosmological simulations to quantify this effect by comparing pure N-body dark matter runs with corresponding hydrodynamical simulations, carried out both in non-radiative, and in dissipative form with cooling and star formation. In both hydrodynamical simulations, the clustering of the gas is suppressed while that of dark matter is boosted at scales k>1 h/Mpc. Despite this counterbalance between dark matter and gas, the clustering of the total matter is suppressed by up to 1 percent at 1<k<10 h/Mpc, while for k ~ 20 h/Mpc it is boosted, up to 2 percent in the non-radiative run and 10 percent in the run with star formation. The stellar mass formed in the latter is highly biased relative to the dark matter in the pure N-body simulation. Using our power spectrum measurements to predict the effect of baryons on the weak lensing signal at 100<l<10000, we find that baryons may change the lensing power spectrum by less than 0.5 percent at l<1000, but by 1 to 10 percent at 1000<l<10000. The size of the effect exceeds the predicted accuracy of future lensing power spectrum measurements and will likely be detected. Precise determinations of cosmological parameters with weak lensing, and studies of small-scale fluctuations and clustering, therefore rely on properly including baryonic physics.

In Figure 6 we inadvertently labeled the proxy circular velocity as the virial velocity of the dark matter halo instead of what is actually plotted, the maximum circular velocity of the dark matter halo. The maximum halo circular velocity is a much better estimate of the disk V_c than is V_vir. This confusion influenced the discussion of the Tully-Fisher relation in our paper. In fact, Figure 6 demonstrates that it is possible to simultaneously reproduce both the local Tully-Fisher relation and luminosity function using semi-analytic techniques applied to the standard LCDM cosmology, thus contradicting previous studies of this issue and our own discussion in Section 3.6.

We model gas cooling in high-resolution N-body simulations in order to investigate the formation of the first generation of stars. We follow a region of a LCDM universe especially selected to contain a rich cluster by the present day. The properties of the dark haloes that form in these sub-solar mass resolution simulations are presented in a companion paper by Gao et al. The first gas clouds able to cool by molecular hydrogen line emission collapse at extremely high redshift, z ~ 47, when the mass of the dark halo is 2.4 x 10^5 Msun/h. By z ~ 30, a substantial population of haloes are capable of undergoing molecular hydrogen cooling although their ability to form stars is dependent on the efficiency of feedback processes such as dissociating Lyman-Werner radiation. The mass of the main halo grows extremely rapidly and, by z ~ 36, its virial temperature has reached 10^4K, at which point gas cooling becomes dominated by more effective atomic line processes. By z ~ 30, a small ``group'' of such potential galaxies will have formed unless prevented from doing so by feedback processes. By this redshift, massive (>~ 100 Msun) population III stars are able to ionise gas well beyond their own host halo and neighbouring HII regions can percolate to form an ionized superbubble. Such patches would be too widely separated to contribute significantly to reionisation at this time. The large number density of early cooling haloes in the pre-reionised universe raises the exciting prospect that this ultra-early generation of stars may be observable as gamma-ray bursts or supernovae.

We use a novel technique to simulate the growth of the most massive progenitor of a protocluster region from redshift z~80, when its mass was about 10 Msun until the present day. Our nested sequence of N-body resimulations allows us to study in detail the structure both of the dark matter object itself and of its environment. Our effective resolution is optimal at redshifts of 49, 29, 12, 5 and 0 when the dominant object has mass 1.2e5, 5e7, 2e10, 3e12 and 8e14 Msun/hrespectively, and contains ~1e6 simulation particles within its virial radius. Extended Press-Schechter theory correctly predicts both this rapid growth and the substantial overabundance of massive haloes we find at early times in regions surrounding the dominant object. Although the large-scale structure in these regions differs dramatically from a scaled version of its present-day counterpart, the internal structure of the dominant object is remarkably similar. Molecular hydrogen cooling could start as early as z~49 in this object, while cooling by atomic hydrogen becomes effective at z~39. If the first stars formed in haloes with virial temperature ~2000, their comoving abundance by z=49 would be similar to that of dwarf galaxies today, while their comoving correlation length would be ~2.5Mpc/h.

We use a very large simulation of the concordance LCDM cosmogony to study the clustering of dark matter haloes. For haloes less massive than about 1e13Msun/h, the amplitude of the two-point correlation function on large scales depends strongly on halo formation time. Haloes that assembled at high redshift are substantially more clustered than those that assembled more recently. The effect is a smooth function of halo formation time and its amplitude increases with decreasing halo mass. At 1e11 Msun/h the ``oldest'' 10% of haloes are more than 5 times more strongly correlated than the ``youngest'' 10%. This unexpected result is incompatible with the standard excursion set theory for structure growth, and it contradicts a fundamental assumption of the halo occupation distribution models often used to study galaxy clustering, namely that the galaxy content of a halo of given mass is statistically independent of its larger scale environment.

The cold dark matter model has become the leading theoretical paradigm for the formation of structure in the Universe. Together with the theory of cosmic inflation, this model makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a novel framework for the quantitative physical interpretation of such surveys. This combines the largest simulation of the growth of dark matter structure ever carried out with new techniques for following the formation and evolution of the visible components. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with next generation surveys.

We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular $\Lambda$ Cold Dark Matter cosmogony. We identify voids by requiring them to be regions of space with a mean overdensity of -0.8 or less. Each of the simulations contains thousands of voids. The distribution of void sizes in the different simulations shows good agreement. Voids very clearly correspond to minima in the smoothed initial density field. We find a universal void mass profile of the form $\rho(<r)/\rho(r_{\rm eff}) \propto \exp[(r/r_{\rm eff})^\alpha]$ where $r_{\rm eff}$ is the effective radius of a void and $\alpha\sim 2$. The mass function of haloes in voids is steeper than that of haloes that populate denser regions. The abundances of void haloes seem to evolve somewhat more strongly between redshifts $\sim 1$ and 0 than the global abundances of haloes. (abridged)

We study the spatial distribution of satellite galaxies by assuming that they follow the dark matter distribution. This assumption is supported by semi-analytical studies based on high-resolution numerical simulations. We find that for a Milky-Way type halo, if only a dozen satellite galaxies are observed, then they can lie on a ``great'' disk with an rms height of about 40 kpc. The normal to the plane is roughly isotropic on the sky. These results are consistent with the observed properties of the satellite galaxies in the Milky Way. If, however, the satellite galaxies follow the distribution of substructure selected by present mass, then great disks similar to the one in the Milky Way are rare and difficult to reproduce, in agreement with the conclusion reached by Kroupa et al. (2004).

We use a large cosmological simulation and a variety of high resolution resimulations of individual cluster and galaxy haloes to study the systematics of subhalo populations over ranges of 1000 in halo mass and 1000 in the ratio of subhalo to parent halo mass. The subhalo populations of different haloes are not scaled copies of each other, but vary systematically with halo properties. On average, the amount of substructure increases with halo mass. At fixed mass, it decreases with halo concentration and with halo formation redshift. These trends are comparable in size to the scatter in subhalo abundance between similar haloes. Averaged over all haloes of given mass, the abundance of low mass subhaloes per unit parent halo mass is independendent of parent mass.It is very similar to the abundance per unit mass of low mass haloes in the universe as a whole, once differing boundary definitions for subhaloes and haloes are accounted for. The radial distribution of subhaloes within their parent haloes is substantially less centrally concentrated than that of the dark matter. It varies at most weakly with the mass(or concentration) of the parent halo and not at all with subhalo mass. It does depend on the criteria used to define the subhalo population considered. About 90% of present-day subhaloes were accreted after z=1 and about 70% after z=0.5. Only about 8% the total mass of all haloes accreted at z=1 survives as bound subhaloes at z=0. For haloes accreted at z=2, the survival mass fraction is just 2%. Subhaloes seen near the centre of their parent typically were accreted earlier and retain less of their original mass than those seen near the edge. These strong systematics mean that comparison with galaxies in real clusters is only possible if the formation of the luminous component is modelled appropriately.

The most luminous galaxies in the present Universe are found at the centers of the most massive dark matter haloes, rich galaxy clusters. In the LCDM cosmology, such massive halo cores are present at redshift z=6 with a comoving number density (as a function of mass interior to ~10 kpc) that is comparable to today's value. The identity of the matter in these central regions is, however, predicted to change as major mergers bring together stars and dark matter from initially well separated sub-units. We use N-body simulations to investigate how these mergers push pre-existing matter outwards in the dominant galaxy while preserving the inner density profile of collisionless matter. It appears that the central regions of large galaxies end up dominated by stars formed in a number of dense cores, well before the last major mergers. The density profile of collisionless matter (stars and dark matter combined) in these central regions appears to be stable and to have attractor-like behavior under merging. This suggests that the baryon loading associated with dissipative contraction and star formation may be erased as subsequent mergers drive the mass distribution back to a universal profile. Such suppression of the effects of baryon loading, along with the early assembly of mass concentrations, may help resolve some apparent challenges to the CDM model for structure formation.

We combine 10 high resolution resimulations of cluster--sized dark haloes with semi--analytic galaxy formation modelling in order to compare the number density and velocity dispersion profiles of cluster galaxies with those of dark matter substructures (subhaloes). While the radial distribution of galaxies follows closely that of the dark matter, the distribution of dark matter subhaloes is much less centrally concentrated. The velocity dispersion profiles of galaxies are also very similar to those of the dark matter, while those for subhaloes are biased high, particularly in the inner regions of clusters. We explain how these differences, already clearly visible in earlier work, are a consequence of the formation of galaxies at the centres of dark matter haloes. Galaxies and subhaloes represent different populations and are not directly comparable. Evolution produces a complex and strongly position--dependent relation between galaxies and the subhaloes in which they reside. This relation can be properly modelled only by appropriate physical representation of the galaxy formation process.