Study on the formation and evolution of galaxies by Large-scale numerical simulations on special-purpose computers.

1. 我们在并行GRAPE系统上开发了一种高效的Barnes-Hut树算法并行实现。使用这个计算代码，现在可以使用~ 10^8粒子在小于星系团的空间尺度上进行宇宙结构形成的模拟。以前只报道过一次这样的模拟。通过对超过3000万个粒子的宇宙学模拟，我们以前所未有的准确性和可靠性确定了CDM宇宙学中暗物质晕中心尖端的斜率。斜率在-1.4左右。我们没有观察到斜率趋于-1的任何趋势。这个斜率不依赖于宇宙学模型。我们找到了卫星星系轨道演化的n体模拟与半解析计算不一致的原因。在最近的高精度n体模拟中，卫星星系的演化被证明比使用钱德拉塞卡动力摩擦公式的近似计算结果要慢得多。我们发现这种差异不是由于动力摩擦公式本身，而是由于将母星系的大小作为库仑对数的上限的常见练习。通过采用理论优选值，即卫星到母星系中心的距离，我们解决了这一差异。在卫星内自由度情况下，我们发现其演化速度比解析动力摩擦公式快。

1. We have developed a high-efficiency parallel implementation of the Barnes-Hut tree algorithm on parallel GRAPE system. Using this calculation code, it is now feasible to perform simulations of cosmological structure formation in the spatial scale smaller than clusters of galaxies using 〜 10^8 particles. Only one such simulation has been previously reported.2. By cosmological simulations with more than 30 million particles, we have determined the slope of the central cusp of the dark-matter halo in CDM cosmology with unprecedented accuracy and reliability. The slope is around -1.4. We have not observed any tendency of the slope to approach to ―1. This slope does not depend on the cosmological models.3. We have found the reason of the discrepancy between N-body simulations and semi-analytic calculations of the orbital evolution of satellite galaxies. The evolution of the satellite galaxies in recent high-accuracy N-body simulations turned out to be significantly slower than the result of approximate calculations using Chandrasekhar's dynamical friction formula. We found that the discrepancy is not due to the dynamical friction formula itself, but due to the common exercise of taking the size of the parent galaxy as the upper limit for the Coulomb logarithm. By taking the theoretical preferred value, namely the distance of the satellite from the center of the parent galaxy, we resolved the discrepancy. In the case of the satellite from internal degree of freedom, we found that the evolution is faster than that for analytic dynamical friction formula.