2025-11-24T03:43:17.676471

Redshift drift in a universe with structure III: Numerical relativity

Koksbang, Heinesen, Macpherson
Measurements of the cosmic redshift drift - the change in redshift of a source over time - will enable independent detection of cosmological expansion thanks to the immense precision soon reached by new facilities such as the Square Kilometer Array Observatory and the Extremely Large Telescope. We conduct the first ever redshift drift computation in fully relativistic cosmological simulations, with the simulations performed with the Einstein Toolkit. We compute the redshift drift over the full skies of 50 synthetic observers in the simulation. We compare all-sky averages for each observer - and across all observers - to the Einstein-de Sitter (EdS) model which represents the large-scale spatially-averaged spacetime of the simulation. We find that at $z\approx0.2$ the mean redshift drift across the sky for all observers deviates from the EdS prediction at the percent level, reducing to $\sim0.1\%$ by $z\approx 1$. However, fluctuations in the redshift drift across the sky are $\sim 10-30\%$ at $z\approx 0.1$ and a few percent at $z\approx 0.5$. Such fluctuations are large enough to potentially exceed the expected precision of upcoming redshift drift measurements. Additionally, we find that along 0.48% of the light rays the redshift drift becomes temporarily positive at very low redshift of $z\lesssim 0.02$. This occurs despite our simulation data being based on a matter-dominated model universe. By including a cosmological constant, we expect a slower growth of structures than in the leading-order EdS space-time, and this may reduce the anisotropy over the observers' skies, although we generally expect our results to hold as order-of-magnitude estimates. Redshift drift is arguably one of the most important measurements to be made by next-generation telescopes. Our results collectively serve as preparation for interpreting such a measurement in the presence of realistic cosmic structures.
academic

Redshift drift in a universe with structure III: Numerical relativity

基本信息

  • 论文ID: 2404.06242
  • 标题: Redshift drift in a universe with structure III: Numerical relativity
  • 作者: Sofie Marie Koksbang, Asta Heinesen, Hayley J. Macpherson
  • 分类: astro-ph.CO (Cosmology and Nongalactic Astrophysics)
  • 发表时间: 2024年4月(arXiv v2: 2024年8月9日)
  • 论文链接: https://arxiv.org/abs/2404.06242

摘要

本文首次在完全相对论性宇宙学模拟中计算了宇宙红移漂移(redshift drift),使用Einstein Toolkit进行模拟。研究在50个合成观测者的全天空范围内计算红移漂移,并与代表模拟大尺度空间平均时空的Einstein-de Sitter(EdS)模型进行比较。结果显示:在z≈0.2时,所有观测者的平均红移漂移与EdS预测在百分比水平上存在偏差,在z≈1时减少到0.1%。然而,天空中红移漂移的涨落在z≈0.1时为10-30%,在z≈0.5时为几个百分点。这些涨落可能超过即将到来的红移漂移测量的预期精度。此外,在0.48%的光线路径上,红移漂移在z≲0.02的极低红移处暂时变为正值。

研究背景与动机

问题定义

红移漂移是指源天体红移随时间的微小变化,它能够直接探测宇宙膨胀。这是首次能够直接测量宇宙加速或减速的观测量,不同于超新星等间接方法。

研究重要性

  1. 直接探测宇宙膨胀:红移漂移提供了直接测量宇宙膨胀历史的独特机会
  2. 下一代望远镜能力:SKA和ELT等设施即将具备测量红移漂移的精度
  3. 理论准备:需要为解释未来观测做好理论准备,特别是考虑真实宇宙结构的影响

现有方法局限性

  • 以往研究主要基于FLRW模型的解析计算或牛顿N体模拟结合微扰理论
  • 缺乏在完全相对论框架下考虑非线性结构影响的研究
  • 对结构如何影响红移漂移测量的理解不够充分

核心贡献

  1. 首次全相对论计算:在完全相对论性宇宙学模拟中首次计算红移漂移
  2. 系统性统计分析:对50个观测者的38,400条光线进行统计分析
  3. 结构影响量化:量化了宇宙结构对红移漂移测量的影响程度
  4. 观测预测:为SKA等下一代设施的红移漂移观测提供理论预期
  5. 正红移漂移发现:发现在物质主导宇宙中仍可能出现正红移漂移

方法详解

任务定义

计算在包含真实宇宙结构的时空中,光子沿测地线传播时红移的时间变化率,并分析其统计性质和空间分布。

模拟架构

Einstein Toolkit数值相对论模拟

  • 使用3+1 BSSN形式演化Einstein方程
  • 初始条件:z̄=1000的EdS背景上的线性扰动
  • 宇宙学流体:完美尘埃流体(P≪ρ)
  • 边界条件:周期性边界条件
  • 网格分辨率:主要结果使用N=256³,盒子大小L=3072h⁻¹Mpc

光线追踪方法

使用mescaline代码进行后处理:

  • 在模拟域中随机放置50个观测者
  • 每个观测者沿HEALPix像素方向发射768条光线
  • 求解测地线和测地线偏差方程
  • 观测者和源都与流体共动

红移漂移计算

红移漂移公式采用多极分解方法:

dzdτO=EEλOλEdλΠ\frac{dz}{d\tau}\bigg|_O = E_E \int_{\lambda_O}^{\lambda_E} d\lambda \, \Pi

其中积分项简化为:

Π13uμuνRμνeμeνuρuσCρμσν\Pi \approx -\frac{1}{3}u_\mu u_\nu R^{\mu\nu} - e_\mu e_\nu u_\rho u_\sigma C^{\rho\mu\sigma\nu}

包含Ricci张量贡献和Weyl张量贡献两部分。

技术创新点

  1. 完全非线性演化:不依赖微扰理论,包含结构形成的完整非线性效应
  2. 多观测者统计:系统研究不同位置观测者看到的红移漂移差异
  3. 全天空分析:计算每个观测者全天空的红移漂移分布和功率谱
  4. 精确数值实现:在mescaline中实现高精度的红移漂移计算

实验设置

模拟参数

  • 主要模拟:N=256³,L=3072h⁻¹Mpc
  • 收敛性测试:N=128³,L=1536h⁻¹Mpc
  • 宇宙学参数:h=0.45,其他采用Planck默认参数
  • 演化范围:从z̄=1000到z̄≈0

观测设置

  • 观测者数量:50个(主要),10个(高分辨率天空)
  • 视线数量:768条/观测者(Nside=8),12,288条/观测者(Nside=32)
  • 观测时间:δτ₀=30年
  • 红移范围:z≈0到z≈1

评价指标

  • 相对于EdS模型的红移漂移偏差:δz/δzEdS-1
  • 统计量:均值、68.1%和95.4%置信区间
  • 角功率谱分析
  • 天空图分析

实验结果

主要结果

平均红移漂移

  • z>0.03时,平均红移漂移与EdS模型偏差<0.2%
  • z≈1时,95.4%置信区间为±2.5%,68.1%置信区间为±1%
  • 最大涨落在z≈1时超过10%

结构影响的组成分析

  • Ricci贡献倾向于负值
  • Weyl贡献倾向于正值
  • 两者在统计上大致抵消,但沿单条光线不完全抵消

正红移漂移现象

  • 0.48%的光线在z≲0.02时出现正红移漂移
  • 分布在5个观测者的185条光线上
  • 物理机制:低红移时Weyl贡献主导

天空分布分析

各向异性特征

  • z≈0.1时天空涨落达几十个百分点
  • z≈0.5时减少到几个百分点
  • 偶极贡献占主导,比四极大一个数量级

角功率谱

  • z≈0.1时偶极幅度平均为3.7%,四极为1.1%
  • z≈0.5时偶极幅度为0.8%,四极为0.1%
  • 大角度多极占主导

观测者间差异

不同位置观测者的平均红移漂移在z≈0.2时与EdS预测差异达~1%水平,这对高精度测量具有重要意义。

相关工作

理论研究

  • 解析非均匀宇宙学模型(LTB、Szekeres等)
  • 微扰理论计算
  • 各向异性模型研究

数值研究

  • 牛顿N体模拟结合微扰理论
  • 本文是首次全相对论数值研究

观测前景

  • SKA预期精度:1-10%
  • ELT的红移漂移测量能力
  • 本文结果表明结构效应可能达到观测精度水平

结论与讨论

主要结论

  1. 平均效应小:结构对平均红移漂移的影响在宇宙学红移处为亚百分比水平
  2. 涨落显著:个别视线的红移漂移涨落可达10%以上
  3. 观测挑战:涨落水平接近或超过SKA等设施的预期精度
  4. 正漂移现象:即使在物质主导宇宙中也可能观测到正红移漂移

局限性

  1. 无宇宙学常数:模拟中Λ=0,真实宇宙中结构增长会被抑制
  2. 位置漂移忽略:未包含源位置漂移效应(预期<1%贡献)
  3. 分辨率限制:小尺度结构可能影响低红移结果
  4. 统计样本:观测者和视线数量相对有限

未来方向

  1. 包含暗能量:在ΛCDM框架下重复计算
  2. 约束模拟:使用本地环境约束的模拟
  3. 更高统计:增加观测者和视线数量
  4. 观测策略:为实际观测提供选择策略

深度评价

优点

  1. 开创性工作:首次在全相对论框架下计算红移漂移,填补重要空白
  2. 方法严谨:使用成熟的Einstein Toolkit和mescaline工具链
  3. 统计充分:50个观测者、38,400条光线提供良好统计
  4. 结果重要:为即将到来的观测提供关键理论预期
  5. 验证完备:包含LTB模型验证和收敛性测试

不足

  1. 物理假设:忽略宇宙学常数可能高估结构效应
  2. 计算近似:忽略位置漂移和部分高阶项
  3. 分辨率依赖:低红移结果对数值分辨率敏感
  4. 统计限制:某些统计量的误差棒较大

影响力

  1. 理论贡献:为红移漂移理论奠定数值基础
  2. 观测指导:为SKA、ELT等项目提供重要参考
  3. 方法示范:展示了数值相对论在宇宙学观测中的应用
  4. 后续研究:将激发更多相关理论和观测研究

适用场景

  1. 观测准备:SKA Phase 2等高精度红移漂移测量
  2. 理论验证:测试宇宙学模型和暗能量理论
  3. 系统误差:评估结构效应对宇宙学参数测量的影响
  4. 方法发展:为其他观测量的数值相对论计算提供模板

参考文献

论文引用了62篇相关文献,涵盖:

  • 红移漂移的早期理论工作1,2
  • SKA和ELT的观测前景3,4
  • 各种宇宙学模型中的红移漂移计算5-14
  • Einstein Toolkit数值相对论工具42,43
  • 相关的宇宙学模拟工作44-47

这项工作为即将到来的红移漂移观测时代提供了重要的理论准备,首次量化了真实宇宙结构对这一关键观测量的影响,对宇宙学和观测天文学都具有重要意义。