A Consistent Implementation of Cluster Strong Lensing in Cosmological Simulation Light Cones

技术摘要：在宇宙学模拟光锥中实现一致性的星系团强引力透镜成像

问题陈述
虽然星系团强引力透镜效应是精密宇宙学的基石，但稳健的理论预测在历史上一直滞后于高质量的观测结果。这种差距源于两个主要问题：一是标准立方体模拟框与长条形几何结构（光锥所需）之间的几何不匹配；二是容纳大质量星系团所需的巨大体积与解析背景源星系及子结构所需的高分辨率之间的根本矛盾。

目前的预测方法通常依赖于混合方法，即从模拟中提取单个透镜并结合解析或观测到的源族群。这些方法往往无法保持跨红移的相关性，或者忽略了相关及不相关的视线（LoS）结构。因此，目前缺乏能够从模拟输出中一致生成强透镜图像的方法，即无法确保主透镜、源以及所有介入物质都以符合标准冷暗物质范式的方式，取自同一个宇宙学体积。

方法论
作者提出了一种完全基于模拟的程序，直接利用 IllustrisTNG 宇宙学流体动力学模拟的粒子数据生成强透镜图像。该方法由三个主要部分组成：

1. 体积重映射（Volume Remapping）： 为了解决几何约束，作者采用了框重映射技术（使用 BoxRemap 代码），将立方体模拟体积转换为右四棱柱。这种几何形状具有两个短轴（代表天球平面）和一个长轴（代表视线）。这种变换在保持总质量和局部结构的同时，不会在单个光锥内引入重复结构（前提是视场相对于框大小被合理选择）。该重映射一致地应用于模拟的每一个快照（snapshot）。
2. 平面提取与平滑（Plane Extraction and Smoothing）： 从重映射后的框中，作者提取质量和光度信息，以创建一系列对应于模拟快照的透镜平面和源平面。
   • 质量平面： 粒子根据快照红移被划分为径向壳层。通过使用最近邻核函数（32 个邻居，Wendland-C2 核）对暗物质、恒星和气体粒子进行平滑处理，生成质量图。为了处理稀疏区域中的弥散质量并避免샷噪声（shot noise），将该壳层的总质量作为一个均匀薄层添加，并减去宇宙的平均物质密度，从而创建质量对比图。
   • 源平面： 恒星粒子被平滑为特定观测滤波器（例如 JWST f200w）下的通量图。通量基于恒星年龄、初始质量和金属丰度进行估算，经过红移处理至源平面，并转换为视星等。对于粒子数不足的子晕，使用解析的 Sersic 轮廓进行替换，以确保完整性。
3. 多平面射线追踪（Multi-Plane Ray Tracing）： 提取的平面被用于多平面射线追踪算法（使用 GLAMER 代码）。来自红移为 $z_s$ 的源平面光线被向后追踪穿过所有红移为 $z_l < z_s$ 的透镜平面。最终通过堆叠生成的透镜图像来产生理论预测。该流程无需手动注入源或解析轮廓；其宇宙学模型和星系形成模型即作为先验条件。

核心贡献与结果
利用 TNG300-1 模拟，作者生成了以质量较大的星系团（$M_{200m} \ge 5 \times 10^{14} M_\odot$）为中心的光锥，透镜红移为 $z_G \approx 0.38$，延伸至源红移 $z_{max} \approx 4.18$。

• 图像生成： 该流程成功地直接从模拟数据中产生了强透镜特征，包括多重像构型、弧和临界曲线，而无需进行人工源注入。生成的图像展示了典型的特征，如切向弧以及由星系团合并和子结构驱动的复杂临界曲线形态。
• 视线结构的影响： 研究通过对比全光锥射线追踪与仅通过主透镜平面的射线追踪，量化了不相关视线（LoS）结构的影响。
  • 图像位置： 不相关的视线结构可以使透镜像的相对位置发生偏移，偏移量可达数角秒（在极端情况下高达 $\sim 10''$），显著超过了典型的强透镜质量建模残差（$\lesssim 1''$）。
  • 临界曲线形态： 包含视线物质会导致星系团主临界曲线面积产生 $\sim 6\%$ 的散射。在某些情况下，主临界曲线的面积增强了 30%，或减少了 50%。
  • 总临界面积： 在星系团势阱最小值 $100''$ 范围内，当包含视线结构时，在 $z_s = 4$ 时总临界面积变化了 $16^{+20}_{-14}\%$。这种增强既来自于视线天体产生的新的临界曲线，也来自于对现有曲线的修改。

意义与局限性
本文声称提供了首次展示从宇宙学流体动力学模拟光锥中一致生成星系团强透镜图像的范例。这种方法提供了一种更自动化、“免手动”的方法来产生大规模模拟观测样本，可用于研究透镜建模中的系统误差、质量推断偏差以及星系-星系强透镜（GGSL）截面张力。

作者指出了一些局限性：

• 大尺度功率： 该方法不会引入初始模拟框中不存在的大尺度功率模态，这可能会导致对大尺度功率谱采样的偏差。
• 分辨率： 该方法受限于模拟的质量分辨率，难以同时解析高质量星系团和低质量子结构/源。
• 弥散质量： 弥散质量被近似为一个均匀薄层，而不是保留其特定结构，尽管这是为了避免샷噪声所必需的。
• 计算成本： 生成单个光锥需要处理每个快照中模拟框内的所有粒子，这在计算上非常昂贵，并会产生巨大的数据存储需求。
• 外部剪切： 当前实现未包含来自光锥外大尺度结构的外部剪切，不过其预期影响估计较小（$\gamma_{ext} \sim 0.02$）。

作者总结道，纳入完整的光锥信息对于理论预测至关重要，并且在解释观测结果时，必须仔细考虑视线结构，因为不相关的结构可以显著改变图像构型和临界曲线属性。