We present numerical simulations on propagation of Ultra-High Energy Cosmic Rays (UHECRs) above $10^{19}$ eV in a structured extragalactic magnetic field (EGMF) and simulate their arrival distributions at the earth. We use the IRAS PSCz catalogue in order to construct a model of the EGMF and source models of UHECRs, both of which reproduce the local structures observed around the Milky Way. We also consider modifications of UHECR arrival directions by the galactic magnetic field. We follow an inverse process of their propagation from the earth and record the trajectories. This enables us to calculate only trajectories of UHECRs arriving at the earth, which saves the CPU time. From these trajectories and our source models, we construct arrival distributions of UHECRs and calculate the harmonic amplitudes and the two point correlation functions of them. We estimate number density of sources which reproduces the Akeno Ground Air Shower Array (AGASA) observation best. As a result, we find that the most appropriate number density of the sources is $\sim 5 \times 10^{-6}$ Mpc$^{-3}$. This constrains the source candidates of UHECRs. We also demonstrate skymaps of their arrival distribution with the event number expected by future experiments and examine how the EGMF affects their arrival distribution. A main result is diffusion of clustering events which are obtained from calculations in the absence of the EGMF. This tendency allows us to reproduce the observed two point correlation function better.

We find a simple exact solution of 6-dimensional braneworld which captures some essential features of warped flux compactification, including a warped geometry, compactification, a magnetic flux, and one or two 3-brane(s). In this setup we analyze how the Hubble expansion rate on each brane changes when the brane tension changes. It is shown that effective Newton's constant resulting from this analysis agrees with that inferred by simply integrating extra dimensions out. Based on the result, a general formula for effective Newton's constant is conjectured and its application to cosmology with type IIB warped string compactification is discussed.

Recently, the High Resolution Fly's Eye (HiRes) experiment claims that there is no small scale anisotropy in the arrival distribution of ultra-high energy cosmic rays (UHECRs) above $E>10^{19}$ eV contrary to the Akeno Giant Air Shower Array (AGASA) observation. In this paper, we discuss the statistical significance of this discrepancy between the two experiments. We calculate arrival distribution of UHECRs above $10^{19}$ eV predicted by the source models constructed using the Optical Redshift Survey galaxy sample. We apply the new method developed by us for calculating arrival distribution in the presence of the galactic magnetic field. The great advantage of this method is that it enables us to calculate UHECR arrival distribution with lower energy ($\sim 10^{19}$ eV) than previous studies within reasonable time by following only the trajectories of UHECRs actually reaching the earth. It has been realized that the small scale anisotropy observed by the AGASA can be explained with the source number density $\sim 10^{-5 \sim -6}$ Mpc$^{-3}$ assuming weak extragalactic magnetic field ($B \le 1$ nG). We find that the predicted small scale anisotropy for this source number density is also consistent with the current HiRes data. We thus conclude that the statement by the HiRes experiment that they do not find small scale anisotropy in UHECR arrival distribution is not statistically significant at present. We also show future prospect of determining the source number density with increasing amount of observed data.

We present numerical simulations on the propagation of UHE protons with energies of $(10^{19.5}-10^{22})$ eV in extragalactic magnetic fields over 1 Gpc. We use the ORS galaxy sample, which allow us to accurately quantify the contribution of nearby sources to the energy spectrum and the arrival distribution, as a source model. We calculate three observable quantities, cosmic ray spectrum, harmonic amplitude, and two point correlation function from our data of numerical simulations. With these quantities, we compare the results of our numerical calculations with the observation. We show that the three observable quantities including the GZK cutoff of the energy spectrum can be reproduced in the case that the number fraction $\sim 10^{-1.7}$ of the ORS galaxies more luminous than -20.5 mag is selected as UHECR sources. In terms of the source number density, this constraint corresponds to $10^{-6}$ Mpc$^{-3}$. However, since mean number of sources within the GZK sphere is only $\sim 0.5$ in this case, the AGASA 8 events above $10^{20.0}$ eV, which do not constitute the clustered events with each other, can not be reproduced. On the other hand, if the cosmic ray flux measured by the HiRes, which is consistent with the GZK cutoff, is correct and observational features about the arrival distribution of UHECRs are same as the AGASA, our source model can explain both the arrival distribution and the flux at the same time. Thus, we conclude that large fraction of the AGASA 8 events above $10^{20}$ eV might originate in the topdown scenarios, or that the cosmic ray flux measured by the HiRes experiment might be better. We also discuss the origin of UHECRs below $10^{20.0}$ eV through comparisons between the number density of astrophysical source candidates and our result ($\sim 10^{-6}$ Mpc$^{-3}$).

We present a new method for calculating arrival distribution of UHECRs including modifications by the galactic magnetic field. We perform numerical simulations of UHE anti-protons, which are injected isotropically at the earth, in the Galaxy and record the directions of velocities at the earth and outside the Galaxy for all of the trajectories. We then select some of them so that the resultant mapping of the velocity directions outside the Galaxy of the selected trajectories corresponds to a given source location scenario, applying Liouville's theorem. We also consider energy loss processes of UHE protons in the intergalactic space. Applying this method to our source location scenario which is adopted in our recent study and can explain the AGASA observation above 4 \times 10^{19} eV, we calculate the arrival distribution of UHECRs including lower energy (E>10^19 eV) ones. We find that our source model can reproduce the large-scale isotropy and the small-scale anisotropy on UHECR arrival distribution above 10^19 eV observed by the AGASA. We also demonstrate the UHECR arrival distribution above 10^19 eV with the event number expected by future experiments in the next few years. The interesting feature of the resultant arrival distribution is the arrangement of the clustered events in the order of their energies, reflecting the directions of the galactic magnetic field. This is also pointed out by Alvarez-Muniz et al.(2002). This feature will allow us to obtain some kind of information about the composition of UHECRs and the magnetic field with increasing amount of data.

We predict the arrival distribution of UHECRs above $4 \times 10^{19}$ eV with the event number expected by future experiments in the next few years. We perform event simulations with the source model which is adopted in our recent study and can explain the current AGASA observation. At first, we calculate the harmonic amplitude and the two point correlation function for the simulated event sets. We find that significant anisotropy on large angle scale will be observed when $\sim 10^3$ cosmic rays above $4 \times 10^{19}$ eV are detected by future experiments. The statistics of the two point correlation function will also increase. The angle scale at which the events have strong correlation with each other corresponds to deflection angle of UHECR in propagating in the EGMF, which in turn can be determined by the future observations. We further investigate the relation between the number of events clustered at a direction and the distance of their sources. Despite the limited amount of data, we find that the C2 triplet events observed by the AGASA may originate from the source within 100 Mpc. Merger galaxy Arp 299 (NGC 3690 + IC 694) is the best candidate for their source. If data accumulate, the UHECR sources within $\sim 100$ Mpc can be identified from observed event clusterings significantly. This will provide some kinds of information about poorly known parameters which influence the propagation of UHECRs, such as extragalactic and galactic magnetic field, chemical composition of observed cosmic rays. Also, we will reveal their origin with our method to identify the sources of UHECR. Finally, we predict the arrival distribution of UHECRs above $10^{20}$ eV, which is expected to be observed if the current HiRes spectrum is correct, and discuss their statistical features and implications.

If Ultra-High Energy Cosmic Rays (UHECRs) are originated from nearby galaxies, modeling of the distribution of nearby galaxies is important to an accurate estimate the source number density of UHECRs. We investigate uncertainty of the selection function of the Optical Redshift Survey (ORS), which we used to construct a source model of UHECRs. The investigation is based on a comparison of numbe counts of ORS galaxies with those of the spectroscopic sample of the Sloan Digital Sky Survey (SDSS) Early Data Release (EDR). We carefully count galaxies in the same absolute magnitude bin from the two samples. We find a slight systematic overestimate of the ORS counts outside 5000 km s$^{-1}$ by about a factor of 2. We revise the selection function of the ORS assuming that the SDSS counts are correct. Our revision is based on the absorption given in the ORS catalog as well as that computed from Schlegel et al. (1998), which is systematically larger than the former by $A_B \sim 0.1$ mag in the region of low absorption. It is found that introduction of Schlegel et al.'s absorption changes one of the parameters of the ORS selection function by more than 10%. The revision should be taken into account in the future analysis of the source number density of UHECRs based on the ORS. Using the revised selection function, we determine the global structure of the Local Supercluster (LSC) with a source model of UHECRs, that is, a number-density model consisting of a uniform spherical halo and an exponential disk. We find that the revision is insignificant in terms of the structure of the LSC. However, the revised selection function will be useful to other studies such as peculiar velocity and correlation function.

We have performed the detailed numerical simulations on the propagation of the UHE protons in the energy range $E=(10^{19.5} - 10^{22.0}$) eV in the relatively strong extra-galactic magnetic field with strength $B= (10, 100)$ nG within about 40 Mpc. In this case, the deflection angles of UHECRs become so large that the no counterparts problem is simply solved. As for the source distribution, we assumed that it is proportional to the number distribution of galaxies within the GZK sphere. We have found many clusters, which mean the small-scale anisotropy, in our simulations. It has been also shown that the observed energy spectrum is well reproduced in our models without any fine-tuned parameter. We have used the correlation value in order to investigate statistically the similarity between the distribution of arrival directions of UHECRs and that of galaxies. We have found that each correlation value for each parameter set begins to converge when the number of the detected events becomes ${\large O}(10^3)$. Since the expected number counts by the experiment of the next generation such as TA, HiRes, Auger, and EUSO are thought to be the order of $10^3$, we will be able to determine the source distribution and values of the parameters in this study in the very near future. Compared with the AGASA data, the significant anisotropy on the arrival directions of UHECRs are found in the analysis of first and second harmonics. This may originate from the incompleteness of the ORS database. This problem may be also solved if the source distribution is slightly changed. For example, this problem may be solved if we assume that UHECRs come from some of the galaxies such as AGNs and radio galaxies.