Jiaxin Han | Durham University (original) (raw)

Papers by Jiaxin Han

While various codes exist to systematically and robustly find haloes and subhaloes in cosmologica... more While various codes exist to systematically and robustly find haloes and subhaloes in cosmological simulations , this is the first work to introduce and rigorously test codes that find tidal debris (streams and other unbound substructure) in fully cosmological simulations of structure formation. We use one tracking and three nontracking codes to identify substructure (bound and unbound) in a Milky Way type simulation from the Aquarius suite ) and post-process their output with a common pipeline to determine the properties of these substructures in a uniform way. By using output from a fully cosmological simulation, we also take a step beyond previous studies of tidal debris that have used simple toy models. We find that both tracking and non-tracking codes agree well on the identification of subhaloes and more importantly, the unbound tidal features associated with them. The distributions of basic properties of the total substructure distribution (mass, velocity dispersion, position) are recovered with a scatter of ∼ 20%. Using the tracking code as our reference, we show that the non-tracking codes identify complex tidal debris with purities of ∼ 40%. Analysing the results of the substructure finders, we find that the general distribution of substructures differ significantly from the distribution of bound subhaloes. Most importantly, both bound and unbound substructures together constitute ∼ 18% of the host halo mass, which is a factor of ∼ 2 higher than the fraction in self-bound subhaloes. However, this result is restricted by the remaining challenge to cleanly define when an unbound structure has become part of the host halo. Nevertheless, the more general substructure distribution provides a more complete picture of a halo's accretion history.

We present a maximum-likelihood weak lensing analysis of the mass distribution in optically selec... more We present a maximum-likelihood weak lensing analysis of the mass distribution in optically selected spectroscopic Galaxy Groups (G3Cv5) in the Galaxy And Mass Assembly (GAMA) survey, using background Sloan Digital Sky Survey (SDSS) photometric galaxies. The scaling of halo mass, Mh, with various group observables is investigated. Our main results are: 1) the measured relations of halo mass with group luminosity, virial volume and central galaxy stellar mass, M⋆, agree very well with predictions from mock group catalogues constructed from a GALFORM semi-analytical galaxy formation model implemented in the Millennium ΛCDM N-body simulation; 2) the measured relations of halo mass with velocity dispersion and projected half-abundance radius show weak tension with mock predictions, hinting at problems in the mock galaxy dynamics and their small scale distribution; 3) the median Mh|M⋆ measured from weak lensing depends more sensitively on the lognormal dispersion in M⋆ at fixed Mh than it does on the median M⋆|Mh. Our measurements suggest an intrinsic dispersion of σlog(M⋆)∼0.15; 4) Comparing our mass estimates with those in the catalogue, we find that the G3Cv5 mass can give biased results when used to select subsets of the group sample. Of the various new halo mass estimators that we calibrate using our weak lensing measurements, group luminosity is the best single-proxy estimator of group mass.

The ever increasing size and complexity of data coming from simulations of cosmic structure forma... more The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and V max function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -to a minor extent -the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

We present a study of a comparison of spin distributions of subhaloes found associated with a hos... more We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles. We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as g... more Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the SUSSING MERGER TREES Workshop in which ten different tree-building algorithms were applied to the same set of halo catalogues and their results compared. Although many of these codes were similar in nature, all algorithms produced distinct results. Our main conclusions are that a useful merger-tree code should possess the following features: (i) the use of particle IDs to match haloes between snapshots; (ii) the ability to skip at least one, and preferably more, snapshots in order to recover subhaloes that are temporarily lost during merging; (iii) the ability to cope with (and ideally smooth out) large, temporary flucuations in halo mass. Finally, to enable different groups to communicate effectively, we defined a common terminology that we used when discussing merger trees and we encourage others to adopt the same language. We also specified a minimal output format to record the results.

We present a Revised IRAS-FSC Redshift Catalogue (RIFSCz) of 60,303 galaxies selected at 60 µm fr... more We present a Revised IRAS-FSC Redshift Catalogue (RIFSCz) of 60,303 galaxies selected at 60 µm from the IRAS Faint Source Catalogue (FSC). This revision merges in data from the WISE All-Sky Data Release, the Tenth SDSS Data Release (DR10), the GALEX All-Sky Survey Source Catalog (GASC), the 2MASS Redshift Survey (2MRS) and the Planck Catalogue of Compact Sources (PCCS). The RIFSCz consists of accurate position, ultra-violet (UV), optical, near-, mid-and far-infrared, sub-millimetre (sub-mm) and/or radio identifications, spectroscopic redshift (if available) or photometric redshift (if possible), predicted far-infrared and sub-mm fluxes ranging from 12 to 1380 µm based upon the best-fit infrared template. We also provide stellar masses, star-formation rates and dust masses derived from the optical and infrared template fits, where possible. 56% of the galaxies in the RIFSCz have spectroscopic redshifts and a further 26% have photometric redshifts obtained through the template-fitting method. At S60 > 0.36 Jy, the 90% completeness limit of the FSC, 93% of the sources in the RIFSCz have either spectroscopic or photometric redshifts. An interesting subset of the catalogue is the sources detected by Planck at sub-mm wavelengths. 1200 sources have a detection at better than 5-σ in at least one Planck band and a further 1186 sources have detections at 3-5σ in at least one Planck band.

We report evidence for extended gamma-ray emission from the Virgo, Fornax and Coma clusters based... more 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 bb channel, or 2 − 10 GeV and > 1 TeV annihilating into µ + µ − 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σ 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 bb and µ + µ − channels in all three clusters.

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Si... more We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and v max dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R 200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when v max is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the subsubhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way ... more We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes ). The subhaloes have been identified by ten subhalo finders at up to five different resolution levels ). Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if the analysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. At low resolution levels the shape measurements of the subhalo distributions consisting of the 13 most massive objects identified by the different subhalo finders are strongly scattered. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.

Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in sim... more Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of subhaloes) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced.

We study how to measure the galaxy merger rate from the observed close pair count. Using a highre... more We study how to measure the galaxy merger rate from the observed close pair count. Using a highresolution N-body/SPH cosmological simulation, we find an accurate scaling relation between galaxy pair counts and merger rates down to a stellar mass ratio of about 1:30. The relation explicitly accounts for the dependence on redshift (or time), on pair separation, and on mass of the two galaxies in a pair. With this relation, one can easily obtain the mean merger timescale for a close pair of galaxies. The use of virial masses, instead of stellar masses, is motivated by the fact that the dynamical friction time scale is mainly determined by the dark matter surrounding central and satellite galaxies. This fact can also minimize the error induced by uncertainties in modeling star formation in the simulation. Since the virial mass can be read from the well-established relation between the virial masses and the stellar masses in observation, our scaling relation can be easily applied to observations to obtain the merger rate and merger time scale. For major merger pairs (1:1-1:4) of galaxies above a stellar mass of 4 × 10 10 h −1 M ⊙ at z = 0.1, it takes about 0.31 Gyr to merge for pairs within a projected distance of 20 h −1 kpc with stellar mass ratio of 1:1, while the time taken goes up to 1.6 Gyr for mergers with stellar mass ratio of 1:4. Our results indicate that a single timescale usually used in literature is not accurate to describe mergers with the stellar mass ratio spanning even a narrow range from 1:1 to 1:4.

We present measurements of the velocity dispersion profile (VDP) for galaxy groups in the final d... more We present measurements of the velocity dispersion profile (VDP) for galaxy groups in the final data release of the Sloan Digital Sky Survey (SDSS). For groups of given mass we estimate the redshift-space cross-correlation function (CCF) with respect to a reference galaxy sample, ξ (s) (r p , π), the projected CCF, w p (r p ), and the real-space CCF, ξ cg (r). The VDP is then extracted from the redshift distortion in ξ (s) (r p , π), by comparing ξ (s) (r p , π) with ξ cg (r). We find that the velocity dispersion (VD) within virial radius (R 200 ) shows a roughly flat profile, with a slight increase at radii below ∼ 0.3R 200 for high mass systems. The average VD within the virial radius, σ v , is a strongly increasing function of central galaxy mass. We apply the same methodology to N -body simulations with the concordance Λ cold dark matter cosmology but different values of the density fluctuation parameter σ 8 , and we compare the results to the SDSS results. We show that the σ v − M * relation from the data provides stringent constraints on both σ 8 and σ ms , the dispersion in log M * of central galaxies at fixed halo mass. Our best-fitting model suggests σ 8 = 0.86 ± 0.03 and σ ms = 0.16 ± 0.03. The slightly higher value of σ 8 compared to the WMAP7 result might be due to a smaller matter density parameter assumed in our simulations. Our VD measurements also provide a direct measure of the dark matter halo mass for central galaxies of different luminosities and masses, in good agreement with the results obtained by from stacking the gravitational lensing signals of the SDSS galaxies.

We present a new model to describe the galaxy-dark matter connection across cosmic time, which un... more We present a new model to describe the galaxy-dark matter connection across cosmic time, which unlike the popular subhalo abundance matching technique is self-consistent in that it takes account of the facts that (i) subhalos are accreted at different times, and (ii) the properties of satellite galaxies may evolve after accretion. Using observations of galaxy stellar mass functions out to zsim4z \sim 4zsim4, the conditional stellar mass function at zsim0.1z\sim 0.1zsim0.1 obtained from SDSS galaxy group catalogues, and the two-point correlation function (2PCF) of galaxies at zsim0.1z \sim 0.1zsim0.1 as function of stellar mass, we constrain the relation between galaxies and dark matter halos over the entire cosmic history from zsim4z \sim 4zsim4 to the present. This relation is then used to predict the median assembly histories of different stellar mass components within dark matter halos (central galaxies, satellite galaxies, and halo stars). We also make predictions for the 2PCFs of high-$z$ galaxies as function of stellar mass. Our main findings are the following: (i) Our model reasonably fits all data within the observational uncertainties, indicating that the Lambda\LambdaLambdaCDM concordance cosmology is consistent with a wide variety of data regarding the galaxy population across cosmic time. (ii) ... [abridged]

Astrophysical Journal, 2011

We develop a new method which measures the projected density distribution w_p(r_p)n of photometri... more We develop a new method which measures the projected density distribution w_p(r_p)n of photometric galaxies surrounding a set of spectroscopically-identified galaxies, and simultaneously the projected correlation function w_p(r_p) between the two populations. In this method we are able to divide the photometric galaxies into subsamples in luminosity intervals when redshift information is unavailable, enabling us to measure w_p(r_p)n and w_p(r_p) as a function of not only the luminosity of the spectroscopic galaxy, but also that of the photometric galaxy. Extensive tests show that our method can measure w_p(r_p) in a statistically unbiased way. The accuracy of the measurement depends on the validity of the assumption in the method that the foreground/background galaxies are randomly distributed and thus uncorrelated with those galaxies of interest. Therefore, our method can be applied to the cases where foreground/background galaxies are distributed in large volumes, which is usually valid in real observations. We applied our method to data from SDSS including a sample of 10^5 LRGs at z~0.4 and a sample of about half a million galaxies at z~0.1, both of which are cross-correlated with a deep photometric sample drawn from the SDSS. On large scales, the relative bias factor of galaxies measured from w_p(r_p) at z~0.4 depends on luminosity in a manner similar to what is found at z~0.1, which are usually probed by autocorrelations of spectroscopic samples. On scales smaller than a few Mpc and at both z~0.4 and z~0.1, the photometric galaxies of different luminosities exhibit similar density profiles around spectroscopic galaxies at fixed luminosity and redshift. This provides clear support for the assumption commonly-adopted in HOD models that satellite galaxies of different luminosities are distributed in a similar way, following the dark matter distribution within their host halos.

We develop a new code, the Hierarchical Bound-Tracing (HBT for short) code, to find and trace dar... more We develop a new code, the Hierarchical Bound-Tracing (HBT for short) code, to find and trace dark matter subhaloes in simulations based on the merger hierarchy of dark matter haloes. Application of this code to a recent benchmark test of finding subhaloes demonstrates that HBT stands as one of the best codes to trace the evolutionary history of subhaloes. The success of the code lies in its careful treatment of the complex physical processes associated with the evolution of subhaloes and in its robust unbinding algorithm with an adaptive source subhalo management. We keep a full record of the merger hierarchy of haloes and subhaloes, and allow growth of satellite subhaloes through accretion from its "satellite-of-satellites", hence allowing mergers among satellites. Local accretion of background mass is omitted, while rebinding of stripped mass is allowed. The justification of these treatments is provided by case studies of the lives of individual subhaloes and by the success in finding the complete subhalo catalogue. We compare our result to other popular subhalo finders and show that HBT is able to well resolve subhaloes in high density environment and keep strict physical track of subhaloes' merger history. This code is fully parallelized and freely available upon request to the authors.

Astrophysical Journal, 2011