We combine measurements of weak gravitational lensing from the CFHTLS-Wide survey, supernovae Ia from CFHT SNLS and CMB anisotropies from WMAP5 to obtain joint constraints on cosmological parameters, in particular, the dark energy equation of state parameter w. We assess the influence of systematics in the data on the results and look for possible correlations with cosmological parameters. We implement an MCMC algorithm to sample the parameter space of a flat CDM model with a dark-energy component of constant w. Systematics in the data are parametrised and included in the analysis. We determine the influence of photometric calibration of SNIa data on cosmological results by calculating the response of the distance modulus to photometric zero-point variations. The weak lensing data set is tested for anomalous field-to-field variations and a systematic shape measurement bias for high-z galaxies. Ignoring photometric uncertainties for SNLS biases cosmological parameters by at most 20% of the statistical errors, using supernovae only; the parameter uncertainties are underestimated by 10%. The weak lensing field-to-field variance pointings is 5%-15% higher than that predicted from N-body simulations. We find no bias of the lensing signal at high redshift, within the framework of a simple model. Assuming a systematic underestimation of the lensing signal at high redshift, the normalisation sigma_8 increases by up to 8%. Combining all three probes we obtain -0.10<1+w<0.06 at 68% confidence (-0.18<1+w<0.12 at 95%), including systematic errors. Systematics in the data increase the error bars by up to 35%; the best-fit values change by less than 0.15sigma. [Abridged]

We explore the radial alignment of subhalos in 2-dimensional projections of cosmological simulations. While most other recent studies focussed on quantifying the signal utilizing the full 3-dimensional spatial information any comparison to observational data has to be done in projection along random lines-of-sight. We have a suite of well resolved host dark matter halos at our disposal ranging from 6 x 10^14 Msun/h down to 6 x 10^13Msun/h. For these host systems we do observe that the major axis of the projected 2D mass distribution of subhalos aligns with its (projected) distance vector to the host's centre. The signal is actually stronger than the observed alignment. However, when considering only the innermost 10-20% of the subhalo's particles for the 2D shape measurement we recover the observed correlation. We further acknowledge that this signal is independent of subhalo mass.

We investigate the structural, kinematical, and chemical properties of globular cluster systems (GCSs) in galaxies of different Hubble types in a self-consistent manner based on high-resolution cosmological N-body simulations combined with semi-analytic models of galaxy and globular cluster (GC) formation. We focus on correlations between the physical properties of GCSs and those of their host galaxies for about 10^5 simulated galaxies located at the centres of dark matter halos (i.e. we do not consider satellite galaxies in sub-halos). Our principal results, which can be tested against observations, are as follows. The majority (about 90%) of GCs currently in halos are formed in low-mass galaxies at redshifts greater than 3 with mean formation redshifts of z = 5.7 (12.7 Gyrs ago) and 4.3 (12.3 Gyrs ago) for metal-poor GCs (MPC) and metal-rich GCs (MRCs), respectively. About 52 % of galaxies with GCs show clear bimodality in their metallicity distribution functions, though less luminous galaxies with M_B fainter than -17 are much less likely to show bimodality owing to little or no MRCs. The number fraction of MRCs does not depend on Hubble type but is generally smaller for less luminous galaxies. The specific frequencies (S_ N) of GCSs are typically higher in ellipticals (S_ N ~ 4.0) than in spirals (S_ N ~ 1.8), and higher again (S_N ~ 5.0) for galaxies located at the centers of clusters of galaxies. The total number of GCs per unit halo mass does not depend strongly on M_B or Hubble type of the host galaxy. The mean metallicities of MPCs and MRCs depend on M_B such that they are higher in more luminous galaxies, though the dependence is significantly weakerfor MPCs.

Recent observations have reported that the Galactic globular clusters (GCs) with unusually extended horizontal-branch (EHB) morphologies show a significantly lower velocity dispersion compared with that of the entire Galactic GC system. We consider that the observed distinctive kinematics of GCs with EHB has valuable information on the formation epochs of GCs and accordingly discuss this observational result based on cosmological N-body simulations with a model of GC formation. We assume that GCs in galaxies were initially formed in low-mass halos at high redshifts and we investigate final kinematics of GCs in their host halos at $z=0$. We find that GCs formed in halos virialized at z>10 show lower velocity dispersions on average than those formed at z>6 for halos with GCs at z=0. We thus suggest that the origin of the observed lower velocity dispersion for the Galactic GCs with EHBs is closely associated with earlier formation epochs (z>10) of halos initially hosting the GCs in the course of the Galaxy formation. Considering that the origin of EHBs can be due to the presence of helium-enhanced second-generation stars in GCs, we discuss the longstanding second parameter problem of GCs in the context of different degrees of chemical pollution in GC-forming gas clouds within low-mass halos virialized at different redshifts.

We investigate formation processes and physical properties of globular cluster systems (GCSs) in galaxies based on high-resolution cosmological simulations with globular clusters. We focus on metal-poor clusters (MPCs) and correlations with their host galaxies by assuming that MPC formation is truncated at a high redshift (z_trun > 6). We find that the correlation between mean metallicities (Z_gc) of MPCs and their host galaxy luminosities (L) flattens from z=z_trun to z=0. We also find that the observed relation (Z_gc ~ L^0.15) in MPCs can be reproduced well in the models with Z_gc ~ L^0.5 at z=z_trun when z_trun ~ 10, if mass-to-light-ratios are assumed to be constant at z=z_trun. However, better agreement with the observed relation is found for models with different mass-to-light-ratios between z=z_trun and z=0. It is also found that the observed color-magnitude relation of luminous MPCs (i.e., ``blue tilts'') may only have a small contribution from the stripped stellar nuclei of dwarf galaxies, which have nuclei masses that correlate with their total mass at z=z_trun. The simulated blue tilts are found to be seen more clearly in more massive galaxies, which reflects the fact that more massive galaxies at z=0 are formed from a larger number of dwarfs with stellar nuclei formed at z>z_trun. The half-number radii (R_e) of GCSs, velocity dispersions of GCSs (sigma), and their host galaxy masses (M_h) are found to be correlated with one another such that R_e ~ M_h^{0.57} and sigma ~ M_h^{0.32}.

We investigate structural properties of old, metal-poor globular clusters (GCs) formed at high redshifts (z>6) and located inside and outside virialized galaxy-scale halos in clusters of galaxies with the total masses of M_CL based on high-resolution cosmological simulations with models of GC formation. We mainly derive the parameter dependences of physical properties of intracluster GCs (ICGCs) based on the results of 14 models. Our principle results are summarized as follows. (1) The projected radial density profiles (Sigma_GC) of ICGCs in clusters with different M_CL can be diverse, though ICGCs have inhomogeneous, asymmetric, and somewhat elongated distributions in most models. If Sigma_GC (R) ~ R^alpha, alpha ranges from -1.5 to -2.5 for GCs in clusters. (2) Although total number of GCs within the central 0.05 Mpc (N_GC,0.05) and 0.2 Mpc (N_GC,0.2) are diverse in different clusters, they can depend weakly on M_CL in such a way that both N_GC,0.05 and N_GC,0.2 are likely to be larger for clusters with larger M_CL. (3) Total number of GCs per cluster masses (specific frequency of GCs for clusters of galaxies) are more likely to be larger in more massive clusters, mainly because a larger number of earlier virialized objects can be located in more massive clusters. (4) Spatial distributions of old GCs in clusters can depend on the truncation epoch of GC formation (z_trun) such that they can be steeper and more compact in the models with higher z_trun. (5) The mean metallicity of ICGCs in a cluster can be smaller than that of GCs within the cluster member galaxy-scale halos by ~ 0.3 in [Fe/H]. Metallicity distribution functions (MDFs) of ICGCs show peak values around [Fe/H] ~ -1.6 and do not have remarkable bimodality.

Using high-resolution numerical simulations, we investigate mass- and luminosity-normalized specific frequencies (T_N and S_N, respectively) of globular cluster systems (GCSs) in order to understand the origin of the observed U-shaped relation between S_N and V-band magnitude (M_V) of their host galaxies. We adopt a biased GC formation scenario in which GC formation is truncated in galaxy halos that are virialized at a later redshift, z_trun. T_N is derived for galaxies with GCs today and converted into S_N for reasonable galaxy mass-to-light-ratios (M/L). We find that T_N depends on halo mass (M_h) in the sense that T_N can be larger in more massive halos with M_h > 10^9 M_sun, if z_trun is as high as 15. We however find that the dependence is too weak to explain the observed S_N-M_V relation and the wide range of S_N in low-mass early-type galaxies with -20.5 < M_V < -16.0 mag for a reasonable constant M/L. The M_V-dependence of S_N for the low-mass galaxies can be well reproduced, if the mass-to-light-ratio M_h/L_V \propto M_h^{\alpha}, where \alpha is as steep as -1. Based on these results, we propose that the origin of the observed U-shaped S_N-M_V relation of GCSs can be understood in terms of the bimodality in the dependence of M_h/L_V on M_h of their host galaxies. We also suggest that the observed large dispersionin S_N in low-mass galaxies is due partly to the large dispersion in T_N.

We explored the clustering properties of Lyman Break Galaxies (LBGs) at z=4 and 5 with an angular two-point correlation function on the basis of the very deep and wide Subaru Deep Field data. We found an apparent dependence of the correlation function slope on UV luminosity for LBGs at both z=4 and 5. More luminous LBGs have a steeper correlation function. To compare these observational results, we constructed numerical mock LBG catalogs based on a semianalytic model of hierarchical clustering combined with high-resolution N-body simulation, carefully mimicking the observational selection effects. The luminosity functions for LBGs predicted by this mock catalog were found to be almost consistent with the observation. Moreover, the overall correlation functions of LBGs were reproduced reasonably well. The observed dependence of the clustering on UV luminosity was not reproduced by the model, unless subsamples of distinct halo mass were considered. That is, LBGs belonging to more massive dark haloes had steeper and larger-amplitude correlation functions. With this model, we found that LBG multiplicity in massive dark halos amplifies the clustering strength at small scales, which steepens the slope of the correlation function. The hierarchical clustering model could therefore be reconciled with the observed luminosity-dependence of the angular correlation function, if there is a tight correlation between UV luminosity and halo mass. Our finding that the slope of the correlation function depends on luminosity could be an indication that massive dark halos hosted multiple bright LBGs (abridged).

We first present the results of numerical simulations on formation processes and physical properties of old globular clusters (GCs) located within clusters of galaxies (``intracluster GCs'') and in between clusters of galaxies (``intercluster GCs''). Our high-resolution cosmological simulations with models of GC formation at high redshifts ($z>6$) show that about 30 % of all GCs in a rich cluster can be ragarded as intracluster GCs that can freely drift being trapped by gravitational potential of the cluster rather than by the cluster member galaxies. The radial surface density profiles of the simulated intracluster GCs are highly likely to be flatter than those of GCs within cluster member galaxies. We also find that about 1% of all GCs formed before $z>6$ are not located within any virialized halos and can be regarded as ``intercluster'' (or ``intergalactic'') GCs. We discuss the dependences of physical properties of intracluster and intercluster GCs on the initial density profiles of GCs within low-mass dark matter halos at high redshifts ($z>6$).

We construct the Numerical Galaxy Catalog ($\nu$GC), based on a semi-analytic model of galaxy formation combined with high-resolution N-body simulations in a $\Lambda$-dominated flat cold dark matter ($\Lambda$CDM) cosmological model. The model includes several essential ingredients for galaxy formation, such as merging histories of dark halos directly taken from N-body simulations, radiative gas cooling, star formation, heating by supernova explosions (supernova feedback), mergers of galaxies, population synthesis, and extinction by internal dust and intervening HI clouds. As the first paper in a series using this model, we focus on basic photometric, structural and kinematical properties of galaxies at present and high redshifts. Two sets of model parameters are examined, strong and weak supernova feedback models, which are in good agreement with observational luminosity functions of local galaxies in a range of observational uncertainty. Both models agree well with many observations such as cold gas mass-to-stellar luminosity ratios of spiral galaxies, HI mass functions, galaxy sizes, faint galaxy number counts and photometric redshift distributions in optical pass-bands, isophotal angular sizes, and cosmic star formation rates. In particular, the strong supernova feedback model is in much better agreement with near-infrared (K'-band) faint galaxy number counts and redshift distribution than the weak feedback model and our previous semi-analytic models based on the extended Press-Schechter formalism. (Abridged)

In this paper, we describe our vectorized and parallelized adaptive mesh refinement (AMR) N-body code with shared time steps, and report its performance on a Fujitsu VPP5000 vector-parallel supercomputer. Our AMR N-body code puts hierarchical meshes recursively where higher resolution is required and the time step of all particles are the same. The parts which are the most difficult to vectorize are loops that access the mesh data and particle data. We vectorized such parts by changing the loop structure, so that the innermost loop steps through the cells instead of the particles in each cell, in other words, by changing the loop order from the depth-first order to the breadth-first order. Mass assignment is also vectorizable using this loop order exchange and splitting the loop into $2^{N_{dim}}$ loops, if the cloud-in-cell scheme is adopted. Here, $N_{dim}$ is the number of dimension. These vectorization schemes which eliminate the unvectorized loops are applicable to parallelization of loops for shared-memory multiprocessors. We also parallelized our code for distributed memory machines. The important part of parallelization is data decomposition. We sorted the hierarchical mesh data by the Morton order, or the recursive N-shaped order, level by level and split and allocated the mesh data to the processors. Particles are allocated to the processor to which the finest refined cells including the particles are also assigned. Our timing analysis using the $\Lambda$-dominated cold dark matter simulations shows that our parallel code speeds up almost ideally up to 32 processors, the largest number of processors in our test.

The mass function of dark halos in a Lambda-dominated cold dark matter (LambdaCDM) universe is investigated. 529 output files from five runs of N-body simulations are analyzed using the friends-of-friends cluster finding algorithm. All the runs use 512^3 particles in the box size of 35 h^{-1}Mpc to 140 h^{-1}Mpc. Mass of particles for 35 h^{-1} Mpc runs is 2.67 times 10^7 h^{-1} M_{solar}. Because of the high mass resolution of our simulations, the multiplicity function in the low-mass range, where the mass is well below the characteristic mass and $\nu = \delta_c / \sigma \lesssim 1.0$, is evaluated in the present work, and is well fitted by the functional form proposed by Sheth & Tormen (ST). However, the maximum value of the multiplicity function from our simulations at $\nu \sim 1$ is smaller, and its low mass tail is shallower when compared with the ST multiplicity function.

It has been suggested that cosmic strings produced at a phase transition in the early universe can be the origin of the extremely high energy cosmic rays (EHCR) observed by AGASA above 10^20 eV. Superheavy cosmic strings with linear mass density of 10^22 g/cm can be indirectly observed through the gravitational lensing effect the distant galaxies. The lensing effect by a long straight object can be characterized by a line of double galaxies or quasars with angular separation of about 5 arcsec. We have searched for aligned double objects from the archived data taken by the Subaru Prime Focus Camera (Suprime-Cam). The Suprime-Cam has a great advantage in observing the wide field of view (30x30 arcmin^2) with high sensitivity (R<26 400s exposure), so it is suitable for this research. In this paper, we describe the result of simulation study for developing the method of searching the objects lensed by cosmic strings, and present the observational result obtained by this method.

We have developed a simulation code with the techniques which enhance both spatial and time resolution of the PM method for which the spatial resolution is restricted by the spacing of structured mesh. The adaptive mesh refinement (AMR) technique subdivides the cells which satisfy the refinement criterion recursively. The hierarchical meshes are maintained by the special data structure and are modified in accordance with the change of particle distribution. In general, as the resolution of the simulation increases, its time step must be shortened and more computational time is required to complete the simulation. Since the AMR enhances the spatial resolution locally, we reduce the time step locally also, instead of shortening it globally. For this purpose we used a technique of hierarchical time steps (HTS) which changes the time step, from particle to particle, depending on the size of the cell in which particles reside. Some test calculations show that our implementation of AMR and HTS is successful. We have performed cosmological simulation runs based on our code and found that many of halo objects have density profiles which are well fitted to the universal profile proposed by Navarro, Frenk, & White (1996) over the entire range of their radius.

We have developed a new parallel tree method which will be called the forest method hereafter. This new method uses the sectional Voronoi tessellation (SVT) for the domain decomposition. The SVT decomposes a whole space into polyhedra and allows their flat borders to move by assigning different weights. The forest method determines these weights based on the load balancing among processors by means of the over-load diffusion (OLD). Moreover, since all the borders are flat, before receiving the data from other processors, each processor can collect enough data to calculate the gravity force with precision. Both the SVT and the OLD are coded in a highly vectorizable manner to accommodate on vector parallel processors. The parallel code based on the forest method with the Message Passing Interface is run on various platforms so that a wide portability is guaranteed. Extensive calculations with 15 processors of Fujitsu VPP300/16R indicate that the code can calculate the gravity force exerted on 10^5 particles in each second for some ideal dark halo. This code is found to enable an N-body simulation with 10^7 or more particles for a wide dynamic range and is therefore a very powerful tool for the study of galaxy formation and large-scale structure in the universe.