论文部分内容阅读
在本文中我们开发了再电离的半数值模拟算法,在半数值模拟中,宇宙中某个区域是否被电离是由该区域内产生的电离光子数与中性氢原子数的比值决定。以此作为电离条件,半数值模拟可以得到宇宙的电离状态。另外,我们将星系形成的半解析模型引入到再电离的半数值算法中,这个半解析模型是2010年基于MillenniumⅡ数值模拟得到的,MillenniumⅡ是世界上最大的宇宙学模拟之一,它具有再电离计算需要的足够大的尺度和足够高的分辨率。此半解析模型为我们提供了模拟中需要的电离光子产生率,我们假设星际介质中的中性氢原子数正比于MillenniumⅡ中的暗物质密度,再电离半数值算法据此进行模拟。不同于基本的半数值模拟方法,我们将星系形成过程与再电离过程联系起来。我们所应用的星系形成模型对低红移近邻宇宙的观测拟合得很好,同时它也成功的拟合了高红移处的光度函数。但是我们发现半解析模型中的恒星成分不能提供足够的电离光子,因此很难在红移z~6处将宇宙电离。为了解决这个问题,我们对半解析模型稍作修改,假设了含有更多大质量恒星的初始质量函数,这种假设与原始的半解析模型存在一定矛盾,但是它还是能够很好地拟合高红移UV光度函数,因此我们接受了修改后的模型作为最终模型。 在此假设下,再电离发生和结束的时间比原有模型提前,得到了符合观测的再电离历史。我们由此分析了第一代恒星对再电离的贡献。