1. 项目概述用机器学习重新定义星系团的“势力范围”在星系天文学里有一个经典问题一个星系团它的“影响范围”到底有多大传统上我们习惯用一些固定的几何边界来回答比如R200即星系团内平均密度达到宇宙临界密度200倍的球体半径。这就像在地图上画一个圈圈内是“城市”圈外是“乡村”。但星系团对成员星系的影响真的会在这个圈上戛然而止吗显然不会。星系团外围是一个复杂的“郊区”星系在进入这个区域时其气体、恒星形成活动、动力学状态就开始发生微妙而持续的变化。这种影响是渐进的、物理的而非几何的。我的这项研究正是试图用机器学习这把“新钥匙”去打开理解这层“郊区”物理的大门。核心思路是与其人为定义一个边界不如让数据自己“说话”。我们训练一个深度神经网络分类器让它学习区分哪些是典型的星系团中心星系哪些是典型的、不受干扰的场星系。我们喂给它的“食物”不是星系的位置而是其内在的物理性质恒星质量占比、气体分数、比恒星形成率、光学颜色、气体和恒星的金属丰度等。如果这个模型学得好那么对于一个处在过渡区的星系模型给出的预测概率或者说“置信度”就会徘徊在0.5附近。通过分析这个预测概率随星系到星系团中心距离的变化我们就能找到一个“过渡区域”——在这里星系的属性既不完全像典型的团成员也不完全像场星系。这个区域就是星系团环境产生可观测影响的真实空间范围。这项工作的价值在于它提供了一种完全数据驱动的方法来量化环境效应。我们不再依赖先验的、可能过于简化的边界定义而是让星系的物理性质本身来揭示环境作用的痕迹。最终我们发现这个机器学习定义的过渡半径与星系团质量存在一个清晰的幂律关系并且其标度行为与暗物质晕的浓度-质量关系惊人地一致。这强烈暗示机器学习捕捉到的正是星系团引力势阱的物理本质如何通过冲压剥离等机制在空间上塑造了星系的性质。2. 核心思路与模型设计从物理问题到机器学习框架2.1 问题重构从边界划分到概率预测传统方法的核心矛盾在于我们用一个静态的、球对称的几何边界如R200去描述一个动态的、非对称的物理过程。星系落入星系团的过程是沿宇宙纤维结构进行的其被“预处理”和“剥离”的效应在远早于穿越R200时就可能开始。因此我们的目标不是找到一个“硬边界”而是刻画一个“软过渡区”。我们将此转化为一个二分类问题给定一个星系的若干物理属性判断它更可能属于“星系团环境”还是“场环境”。这里的“环境”标签在训练初期需要一个人为但合理的定义。我们采用了两种方案基于距离的标签以某个半径如0.5倍R200m内的星系作为“团星系”正样本以足够远处如5倍R200m以外的星系作为“场星系”负样本。中间的星系作为待预测的“灰色地带”。基于动力学子结构的标签利用SUBLINK等算法识别的、真正属于星系团动力学子结构的星系作为正样本。注意标签的定义是模型的“监督信号”其合理性至关重要。基于距离的标签简单直观但可能混入前景/背景星系。基于动力学的标签物理上更干净但依赖于特定的算法和模拟数据。我们的收敛性测试表明只要正负样本在属性空间上有足够区分度模型对标签定义的细节并不敏感这增强了结果的鲁棒性。2.2 特征工程如何让星系“开口说话”模型的好坏一半取决于喂给它的数据。我们精心挑选了15个物理特征覆盖了星系的恒星、气体、动力学三个基本方面并进行了标准化和变换处理以确保模型能高效学习。表1输入特征的缩放与变换核心摘要特征原始单位变换物理类别恒星质量比 (M*/Mtot)—log10(x)恒星气体分数 (Mgas/Mbaryon)—arcsinh(axb)气体比恒星形成率 (sSFR)yr⁻¹arcsinh(axb)气体/恒星颜色 (g-r)mag—恒星气体金属丰度 (Zgas/Z⊙)—arcsinh(axb)气体自旋幅度—log10(x)动力学恒星金属丰度 (Z*/Z⊙)—log10(x)恒星恒星半质量半径 (R*,1/2)kpclog10(x)恒星暗物质质量 (MDM)M⊙log10(x)—暗物质半质量半径 (RDM,1/2)kpclog10(x)—最大速度半径 (Rvmax)kpclog10(x)动力学气体半质量半径 (Rgas,1/2)kpcarcsinh(axb)气体速度大小 (|v|)km s⁻¹—动力学速度弥散 (σv)km s⁻¹—动力学最大速度 (vmax)km s⁻¹—动力学特征选择与变换的考量对数变换 (log10)适用于跨越多个数量级的特征如质量、半径。这能将幂律关系转化为线性关系更符合神经网络的胃口。反双曲正弦变换 (arcsinh)这是本项目的关键技巧之一。对于像气体分数、sSFR这类分布严重偏斜有很多接近零的值也有少数很大的值的特征简单的对数变换会因log(0)而失效。arcsinh(x) ≈ log(2x)当x很大时≈ x当x很小时。我们通过参数a和b将每个特征的分布“对称化”使其更接近高斯分布极大加速了模型训练的收敛速度。物理类别分组我们将特征分为“恒星”、“气体”、“动力学”三组。这不仅有助于后续分析不同物理过程的主导尺度也在特征选择时提供了依据。例如我们的基准模型只使用了6个核心特征恒星质量比、气体分数、sSFR、颜色、气体和恒星金属丰度。这涵盖了星系质量、恒星形成、化学演化等最关键的信息。2.3 模型架构与训练一个轻量但高效的分类器我们没有选择复杂的卷积或图神经网络而是采用了一个结构清晰的多层感知机。在宇宙学模拟数据量巨大但特征维度不高的情况下MLP往往能以更低的计算成本获得优异的表现。网络架构 (6 → 48 → 24 → 12 → 1)这是一个四层全连接网络。输入层6个神经元对应6个核心特征。三个隐藏层神经元数逐层减半这种“漏斗形”设计有助于网络逐层提取和压缩抽象信息。输出层1个神经元通过Sigmoid函数输出一个介于0到1之间的概率值代表该星系属于“场星系”的置信度。防止过拟合的三大“护法”Dropout (p0.10)在每个隐藏层的激活值后以10%的概率随机将其置零。这强迫网络不能过度依赖任何一个神经元必须学习到冗余的、分布式的特征表示是防止过拟合的利器。权重衰减 (L2正则化 λ10⁻⁴)我们使用AdamW优化器它对所有权重施加一个小的L2惩罚项。这相当于在损失函数中增加了一项惩罚过大的权重值促使模型学习更平滑、更简单的函数避免对训练数据中的噪声过于敏感。批量归一化 (BatchNorm)在每个隐藏层的线性变换之后、激活函数之前插入批量归一化层。它会对这一批数据的激活值进行标准化减均值、除以标准差并学习一个缩放和偏移参数。这有重好处稳定内部协方差使得训练过程对学习率不那么敏感提供轻微的正则化效果大大加速训练收敛。训练细节优化器AdamW初始学习率 η0 5×10⁻⁵。学习率调度余弦退火。学习率随训练轮数epoch按余弦函数从初始值衰减到最小值 ηmin 9×10⁻⁶周期 Tmax 200轮。这种调度方式能在训练后期进行更精细的参数微调。早停如果验证集损失在连续5个epoch内没有改善则停止训练并回滚到验证损失最低的模型检查点。这是防止过拟合的最后一道保险。损失函数加权二元交叉熵。由于“团星系”和“场星系”的样本数量通常不均衡我们在损失函数中为样本较少的类别通常是“场星系”赋予更高的权重w N_neg / N_pos确保模型不会简单地偏向多数类。实操心得在宇宙学模拟数据上训练MLP批量归一化是必须的。模拟数据通常来自不同的 halo、不同的宇宙学体积即使经过特征变换其分布也可能存在批次间的差异。BatchNorm 能有效平滑这种差异我实测下来不加BatchNorm的模型收敛速度慢且不稳定验证集损失波动很大。此外学习率不宜设高5e-5是一个安全的起点配合余弦退火能取得稳定且优异的性能。3. 从模型输出到物理洞察分析流程全解析3.1 模型评估与推断训练完成后我们进行两次推断在测试集上评估计算AUROCROC曲线下面积和PR-AUC精确率-召回率曲线下面积等指标定量评估模型的分类性能。一个好的模型其AUROC应接近1这表明模型能很好地区分两类星系。在全数据集上预测这才是重头戏。我们将模拟中所有的星系包括那些没有明确“环境”标签的、位于过渡区的星系输入训练好的最佳模型。对于每个星系i模型会输出一个原始分数z_ilogit。z_i经过Sigmoid函数得到概率p_i。我们同时记录该星系的宿主星系团标识c_i以及它到宿主团中心的归一化距离R_i r_gal / R200m。3.2 构建概率剖面与寻找过渡区接下来我们按宿主星系团将星系分组。对于每个团h我们得到两个向量z_h: 该团所有成员星系模型预测的logit值。ℓ_h: 这些星系归一化距离的对数值log10(R_i)。我们的目标是找到那个“过渡半径”r0即模型认为星系从“团属性”转变为“场属性”的位置。这对应着logit值z从负值团属性跨越零变为正值场属性的径向位置。我们采用两种互补的方法来计算r0逐团中位数法对每个星系团我们将其星系的(log10(R_i), z_i)数据点进行平滑或拟合找到z0对应的log10(R_0)。然后将所有团的R_0按质量分档取中位数得到该质量档的代表性r_med。这种方法保留了星系团之间的个体差异散射。叠加剖面法将同一质量档的所有星系团的星系按log10(R)分到统一的径向 bins 中然后计算每个bin内z值的中位数从而得到一条该质量档的“平均”概率剖面P_stack(r)。这条剖面与z0水平线的交点即定义为r_stack。这种方法能反映整体趋势平滑掉个体噪声。图3展示了叠加剖面法的结果。可以看到对于不同质量的星系团logit剖面都呈现清晰的“S”形在团中心内部r 0.5 R200mlogit为显著的负值模型高度确信这些是团星系在约0.5到1.5倍R200m之间logit快速上升并穿过零点在更远的距离r 2 R200mlogit趋于一个稳定的正值对应场星系。这个S形转变区就是我们寻找的“环境影响过渡区”。3.3 收敛性测试结果是否依赖于标签定义这是一个必须回答的关键问题我们找到的r0会不会只是因为当初画“团”和“场”的标签圈比如0.5和5倍R200m而人为制造出来的为此我们进行了系统的收敛性测试。固定内边界为0.5倍R200m将外边界从1.0倍逐步增加到5.0倍R200m步长0.1用每一组标签训练模型并观察预测的r0如何变化。结果非常振奋人心如图4所示尽管用于定义“场星系”的外边界在不断外推但模型预测的过渡半径r0很快大约在外边界超过3倍R200m后就稳定在一个平台值上其波动范围远小于r0本身的误差棒。这意味着只要用于训练的“团”和“场”样本在属性空间上是足够纯净和分离的模型就能稳健地学习到两者之间的本质区别而这个区别所对应的空间过渡位置是独立于标签的精确距离阈值的。这强有力地证明了r0反映的是星系属性随环境变化的真实物理信号而非我们人为标注的产物。4. 核心发现与物理机制解读4.1 环境影响的空间范围超越R200我们的主要发现如图5左、中面板所示。机器学习预测的过渡半径r0其典型值在0.8 到 1.2倍 R200m之间并且随着宿主星系团质量的增加而轻微增大。这个结果有两层重要含义星系团的影响范围比传统边界更远平均而言星系团环境对星系属性的可观测影响一直延伸到约1.2倍R200m处。这意味着有相当一部分在传统定义上处于“星系团外”的星系其实已经在物理上受到了星系团环境的显著影响。R200m作为一个基于球对称过密度的几何定义并不能作为环境效应的“硬边界”。与 splashback 半径的对比图中我们还叠加了理论预言的 splashback 半径由Diemer 2020和More et al. 2015给出。Splashback半径是暗物质粒子轨道最远点的中位数位置是一个动力学的边界。有趣的是splashback半径随质量减小而我们的r0随质量增大趋势完全相反。这说明星系属性发生转变的位置并非由动力学边界直接决定。它标记的是环境物理过程如气体剥离对星系产生显著影响的区域这个区域与动力学边界相关但并不重合。4.2 本征散射模糊的过渡区图5中的垂直误差棒很大这并非测量误差而是本征散射——它反映了真实的物理多样性。有些星系团的过渡发生得靠内有些则靠外即使在同一个团内不同星系因其轨道、倾角、固有质量等因素开始被“预处理”的位置也不同。图6通过分析不同logit区间的星系径向分布生动地展示了这一点高置信团星系(logit -5紫色曲线)高度集中在r 1倍R200m以内分布尖锐。高置信场星系(logit 5浅绿色曲线)几乎完全分布在r 1倍R200m以外。模糊样本(-0.23 logit 0.23青绿色曲线)分布非常弥散从团内部一直延伸到远郊。这些正是处于“正在转变中”的星系。这个模糊的、弥散的过渡区图像比一个尖锐的边界更符合物理直觉。它告诉我们星系团环境的影响是一个渐进的、概率性的过程。我们的模型成功捕捉到了这种复杂性其输出logit本身就可以作为一个数据驱动的、连续的环境“影强度”指标。4.3 质量依赖性的物理根源浓度-质量关系为什么大质量星系团的过渡半径r0反而更大这背后隐藏着深刻的物理联系。我们发现r0与星系团质量M200m之间存在幂律关系r0 ∝ M200m^α拟合指数α ≈ 0.10见图7。这个指数与宇宙学模拟中建立的暗物质晕浓度-质量关系的指数惊人地一致c(M) ∝ M^{-β}其中β ≈ 0.10。这绝非巧合。低质量暗物质晕形成更早宇宙背景密度更高因此其物质分布更“紧凑”即浓度c更高。高浓度意味着密度剖面在内部更陡峭向外下降得更快。在经典的冲压剥离框架中剥离条件取决于星系所处位置的星系团内介质密度ρ_ICM(r)。对于一个高浓度陡峭密度轮廓的晕ρ_ICM(r)在较小的r/R200处就下降到不足以有效剥离星系气体的阈值。反之对于低浓度平缓轮廓的大质量晕其气体密度在外围下降更慢使得剥离过程可以发生在更大的r/R200处。因此机器学习模型从星系属性中学到的过渡半径r0本质上反比于晕的浓度r0 ∝ 1/c ∝ M^{β}。我们的测量α ≈ β ≈ 0.10完美验证了这一推论。这表明模型并非在学习一个任意的模式而是直接探测到了底层引力势的物理结构通过浓度如何设定了环境过程生效的空间尺度。4.4 不同物理过程的响应尺度我们训练了仅使用特定类别特征的模型如仅用恒星属性、仅用气体属性、仅用动力学属性以探究不同物理过程主导的转变发生在何处。图5右面板和图9展示了不同特征组合预测的r0之比。结果显示动力学属性如速度、自旋定义的过渡区最靠内。这是因为动力学状态如轨道、速度弥散是星系与团引力势直接相互作用的结果在星系深入势阱内部时变化最为剧烈。气体属性和恒星属性定义的过渡区相对更靠外且彼此接近难以完全区分。这与观测相符气体尤其是冷气体作为星系中最“脆弱”的组分会最先受到环境的影响如冲压剥离、潮汐剥离导致恒星形成率下降sSFR变化和颜色变红这些是相对较早的信号。而恒星质量、金属丰度等属性变化可能需要更长时间。这个排序为我们理解星系“被 quenching”停止形成恒星的序列提供了线索可能先是在外围经历气体剥离影响气体属性然后随着深入恒星形成停止影响sSFR、颜色最后动力学状态被完全改造。4.5 与传统边界定义的对比我们同样尝试了直接用传统边界R200m R200c Rsp作为二分类标签即内部为团外部为场来训练模型。结果如图8所示这些模型预测的r0与它们用作标签的边界值并不一致且误差棒巨大。这清晰地表明一个固定的几何或动力学边界并不能完美对应星系属性发生系统性变化的区域。星系在穿越splashback半径之前其属性可能已经开始变化预处理而在R200之内也可能存在一些尚未被完全影响的星系。我们的概率化、数据驱动的方法正是为了捕捉这种物理上的模糊性和渐进性。5. 实操经验、挑战与未来方向5.1 数据准备与特征处理的坑模拟数据的异构性不同宇宙学模拟如IllustrisTNG, EAGLE, SIMBA输出的星系属性其定义、单位、甚至计算方式可能有细微差别。在合并使用或与观测对比前必须进行极其仔细的交叉验证和归一化。我建议为每个模拟数据集单独编写数据读取和预处理管道并保存一份详细的“数据护照”记录每个特征的出处和变换历史。缺失值与异常值模拟数据中也可能存在极端值或无效值例如某些矮星系的气体质量接近零导致气体分数计算溢出。arcsinh变换虽然能处理零值但对于极端负值或由于数值误差产生的异常值仍需警惕。稳健的做法是先进行描述性统计如绘制直方图根据物理意义设定合理的截断范围再进行变换。训练/验证/测试集划分绝不能随机划分因为星系是成团的随机划分会导致同一个星系团的成员同时出现在训练集和测试集造成数据泄露严重高估模型性能。必须按星系团ID进行划分确保用于训练、验证、测试的星系来自完全不同的、独立的星系团集合。5.2 模型训练与调参心得学习率是命门对于宇宙学数据学习率宁小勿大。我从5e-5开始配合余弦退火在大多数任务上都能稳定收敛。如果使用更大的批量大小如512或1024可以适当增大学习率但需密切监控验证集损失。早停的耐心值设置早停的“耐心”参数需要平衡。太短如3可能导致在损失平台期提前停止错过后续下降太长如20则浪费计算资源。我发现在损失曲线进入明显的“平原期”后再观察5-10个epoch是合理的通常设置耐心为5-10。类别权重的计算加权交叉熵中的权重w N_neg / N_pos需要谨慎计算。这里的N_neg和N_pos应该是训练集中的数量而不是全集。并且如果正负样本比例极度失衡如超过1:10可能需要更激进的重采样策略如对少数类过采样或使用Focal Loss等更高级的损失函数。5.3 结果解释的注意事项相关性不等于因果模型找到了属性与环境的关联并预测了一个过渡半径。但这并不意味着在这个半径上“发生”了某个特定事件如冲压剥离启动。它只是表明在这个统计意义上星系团的集体影响开始显著地改变星系属性的整体分布。解释时需要结合具体的物理模型。外推需谨慎我们的模型是在特定模拟如TNG300上训练的。将其直接应用于观测或其他模拟时必须考虑域适应问题。观测数据存在选择效应、测量误差和不同的物理过程如AGN反馈模型不同。一个实用的做法是在目标数据集上对模型进行少量样本的微调或者使用域适应技术。“黑箱”的可解释性虽然我们通过分析不同特征组合的输出来推断物理但深度神经网络的具体决策路径仍是复杂的。可以辅以SHAP或LIME等可解释性AI工具量化每个输入特征对单个星系预测结果的贡献度这能提供更细致的物理洞察。5.4 未来扩展方向引入时间维度当前工作是静态快照。下一步可以引入星系的轨道信息如最近通过 pericenter 的时间或形成历史训练模型预测星系“被影响”的阶段或未来演化轨迹实现更动态的环境诊断。结合空间信息当前模型只用了星系自身的属性。可以尝试图神经网络将星系的局部环境如最近邻星系的距离、所属纤维结构的形态也作为输入捕捉环境效应的非局部性。应用于大规模巡天数据将训练好的模型应用于SDSS、DESI、Euclid等巡天数据。最大的挑战是处理观测误差和选择函数。一个可行的路径是用高保真模拟生成带噪声的、符合观测限制的“模拟观测”数据来训练模型再应用到真实数据。探索更复杂的架构对于更高维的特征集如积分场光谱数据衍生的多维度信息可以尝试注意力机制或Transformer架构让模型自动关注最重要的特征交互。这项研究最让我兴奋的一点是它展示了机器学习不仅仅是一个“黑箱”分类工具。当与深厚的物理洞察结合时它能够揭示出数据中隐藏的、符合物理规律的标度关系和内在联系。我们找到的r0 ∝ M^{0.1}关系就是这种结合产出的“物理可解释”的结晶。它告诉我们在追求预测精度的同时设计合理的实验、提出正的问题、并深入理解模型输出背后的物理才是将AI真正转化为天体物理学发现引擎的关键。
机器学习量化星系团环境影响：从数据驱动边界到物理机制解读
1. 项目概述用机器学习重新定义星系团的“势力范围”在星系天文学里有一个经典问题一个星系团它的“影响范围”到底有多大传统上我们习惯用一些固定的几何边界来回答比如R200即星系团内平均密度达到宇宙临界密度200倍的球体半径。这就像在地图上画一个圈圈内是“城市”圈外是“乡村”。但星系团对成员星系的影响真的会在这个圈上戛然而止吗显然不会。星系团外围是一个复杂的“郊区”星系在进入这个区域时其气体、恒星形成活动、动力学状态就开始发生微妙而持续的变化。这种影响是渐进的、物理的而非几何的。我的这项研究正是试图用机器学习这把“新钥匙”去打开理解这层“郊区”物理的大门。核心思路是与其人为定义一个边界不如让数据自己“说话”。我们训练一个深度神经网络分类器让它学习区分哪些是典型的星系团中心星系哪些是典型的、不受干扰的场星系。我们喂给它的“食物”不是星系的位置而是其内在的物理性质恒星质量占比、气体分数、比恒星形成率、光学颜色、气体和恒星的金属丰度等。如果这个模型学得好那么对于一个处在过渡区的星系模型给出的预测概率或者说“置信度”就会徘徊在0.5附近。通过分析这个预测概率随星系到星系团中心距离的变化我们就能找到一个“过渡区域”——在这里星系的属性既不完全像典型的团成员也不完全像场星系。这个区域就是星系团环境产生可观测影响的真实空间范围。这项工作的价值在于它提供了一种完全数据驱动的方法来量化环境效应。我们不再依赖先验的、可能过于简化的边界定义而是让星系的物理性质本身来揭示环境作用的痕迹。最终我们发现这个机器学习定义的过渡半径与星系团质量存在一个清晰的幂律关系并且其标度行为与暗物质晕的浓度-质量关系惊人地一致。这强烈暗示机器学习捕捉到的正是星系团引力势阱的物理本质如何通过冲压剥离等机制在空间上塑造了星系的性质。2. 核心思路与模型设计从物理问题到机器学习框架2.1 问题重构从边界划分到概率预测传统方法的核心矛盾在于我们用一个静态的、球对称的几何边界如R200去描述一个动态的、非对称的物理过程。星系落入星系团的过程是沿宇宙纤维结构进行的其被“预处理”和“剥离”的效应在远早于穿越R200时就可能开始。因此我们的目标不是找到一个“硬边界”而是刻画一个“软过渡区”。我们将此转化为一个二分类问题给定一个星系的若干物理属性判断它更可能属于“星系团环境”还是“场环境”。这里的“环境”标签在训练初期需要一个人为但合理的定义。我们采用了两种方案基于距离的标签以某个半径如0.5倍R200m内的星系作为“团星系”正样本以足够远处如5倍R200m以外的星系作为“场星系”负样本。中间的星系作为待预测的“灰色地带”。基于动力学子结构的标签利用SUBLINK等算法识别的、真正属于星系团动力学子结构的星系作为正样本。注意标签的定义是模型的“监督信号”其合理性至关重要。基于距离的标签简单直观但可能混入前景/背景星系。基于动力学的标签物理上更干净但依赖于特定的算法和模拟数据。我们的收敛性测试表明只要正负样本在属性空间上有足够区分度模型对标签定义的细节并不敏感这增强了结果的鲁棒性。2.2 特征工程如何让星系“开口说话”模型的好坏一半取决于喂给它的数据。我们精心挑选了15个物理特征覆盖了星系的恒星、气体、动力学三个基本方面并进行了标准化和变换处理以确保模型能高效学习。表1输入特征的缩放与变换核心摘要特征原始单位变换物理类别恒星质量比 (M*/Mtot)—log10(x)恒星气体分数 (Mgas/Mbaryon)—arcsinh(axb)气体比恒星形成率 (sSFR)yr⁻¹arcsinh(axb)气体/恒星颜色 (g-r)mag—恒星气体金属丰度 (Zgas/Z⊙)—arcsinh(axb)气体自旋幅度—log10(x)动力学恒星金属丰度 (Z*/Z⊙)—log10(x)恒星恒星半质量半径 (R*,1/2)kpclog10(x)恒星暗物质质量 (MDM)M⊙log10(x)—暗物质半质量半径 (RDM,1/2)kpclog10(x)—最大速度半径 (Rvmax)kpclog10(x)动力学气体半质量半径 (Rgas,1/2)kpcarcsinh(axb)气体速度大小 (|v|)km s⁻¹—动力学速度弥散 (σv)km s⁻¹—动力学最大速度 (vmax)km s⁻¹—动力学特征选择与变换的考量对数变换 (log10)适用于跨越多个数量级的特征如质量、半径。这能将幂律关系转化为线性关系更符合神经网络的胃口。反双曲正弦变换 (arcsinh)这是本项目的关键技巧之一。对于像气体分数、sSFR这类分布严重偏斜有很多接近零的值也有少数很大的值的特征简单的对数变换会因log(0)而失效。arcsinh(x) ≈ log(2x)当x很大时≈ x当x很小时。我们通过参数a和b将每个特征的分布“对称化”使其更接近高斯分布极大加速了模型训练的收敛速度。物理类别分组我们将特征分为“恒星”、“气体”、“动力学”三组。这不仅有助于后续分析不同物理过程的主导尺度也在特征选择时提供了依据。例如我们的基准模型只使用了6个核心特征恒星质量比、气体分数、sSFR、颜色、气体和恒星金属丰度。这涵盖了星系质量、恒星形成、化学演化等最关键的信息。2.3 模型架构与训练一个轻量但高效的分类器我们没有选择复杂的卷积或图神经网络而是采用了一个结构清晰的多层感知机。在宇宙学模拟数据量巨大但特征维度不高的情况下MLP往往能以更低的计算成本获得优异的表现。网络架构 (6 → 48 → 24 → 12 → 1)这是一个四层全连接网络。输入层6个神经元对应6个核心特征。三个隐藏层神经元数逐层减半这种“漏斗形”设计有助于网络逐层提取和压缩抽象信息。输出层1个神经元通过Sigmoid函数输出一个介于0到1之间的概率值代表该星系属于“场星系”的置信度。防止过拟合的三大“护法”Dropout (p0.10)在每个隐藏层的激活值后以10%的概率随机将其置零。这强迫网络不能过度依赖任何一个神经元必须学习到冗余的、分布式的特征表示是防止过拟合的利器。权重衰减 (L2正则化 λ10⁻⁴)我们使用AdamW优化器它对所有权重施加一个小的L2惩罚项。这相当于在损失函数中增加了一项惩罚过大的权重值促使模型学习更平滑、更简单的函数避免对训练数据中的噪声过于敏感。批量归一化 (BatchNorm)在每个隐藏层的线性变换之后、激活函数之前插入批量归一化层。它会对这一批数据的激活值进行标准化减均值、除以标准差并学习一个缩放和偏移参数。这有重好处稳定内部协方差使得训练过程对学习率不那么敏感提供轻微的正则化效果大大加速训练收敛。训练细节优化器AdamW初始学习率 η0 5×10⁻⁵。学习率调度余弦退火。学习率随训练轮数epoch按余弦函数从初始值衰减到最小值 ηmin 9×10⁻⁶周期 Tmax 200轮。这种调度方式能在训练后期进行更精细的参数微调。早停如果验证集损失在连续5个epoch内没有改善则停止训练并回滚到验证损失最低的模型检查点。这是防止过拟合的最后一道保险。损失函数加权二元交叉熵。由于“团星系”和“场星系”的样本数量通常不均衡我们在损失函数中为样本较少的类别通常是“场星系”赋予更高的权重w N_neg / N_pos确保模型不会简单地偏向多数类。实操心得在宇宙学模拟数据上训练MLP批量归一化是必须的。模拟数据通常来自不同的 halo、不同的宇宙学体积即使经过特征变换其分布也可能存在批次间的差异。BatchNorm 能有效平滑这种差异我实测下来不加BatchNorm的模型收敛速度慢且不稳定验证集损失波动很大。此外学习率不宜设高5e-5是一个安全的起点配合余弦退火能取得稳定且优异的性能。3. 从模型输出到物理洞察分析流程全解析3.1 模型评估与推断训练完成后我们进行两次推断在测试集上评估计算AUROCROC曲线下面积和PR-AUC精确率-召回率曲线下面积等指标定量评估模型的分类性能。一个好的模型其AUROC应接近1这表明模型能很好地区分两类星系。在全数据集上预测这才是重头戏。我们将模拟中所有的星系包括那些没有明确“环境”标签的、位于过渡区的星系输入训练好的最佳模型。对于每个星系i模型会输出一个原始分数z_ilogit。z_i经过Sigmoid函数得到概率p_i。我们同时记录该星系的宿主星系团标识c_i以及它到宿主团中心的归一化距离R_i r_gal / R200m。3.2 构建概率剖面与寻找过渡区接下来我们按宿主星系团将星系分组。对于每个团h我们得到两个向量z_h: 该团所有成员星系模型预测的logit值。ℓ_h: 这些星系归一化距离的对数值log10(R_i)。我们的目标是找到那个“过渡半径”r0即模型认为星系从“团属性”转变为“场属性”的位置。这对应着logit值z从负值团属性跨越零变为正值场属性的径向位置。我们采用两种互补的方法来计算r0逐团中位数法对每个星系团我们将其星系的(log10(R_i), z_i)数据点进行平滑或拟合找到z0对应的log10(R_0)。然后将所有团的R_0按质量分档取中位数得到该质量档的代表性r_med。这种方法保留了星系团之间的个体差异散射。叠加剖面法将同一质量档的所有星系团的星系按log10(R)分到统一的径向 bins 中然后计算每个bin内z值的中位数从而得到一条该质量档的“平均”概率剖面P_stack(r)。这条剖面与z0水平线的交点即定义为r_stack。这种方法能反映整体趋势平滑掉个体噪声。图3展示了叠加剖面法的结果。可以看到对于不同质量的星系团logit剖面都呈现清晰的“S”形在团中心内部r 0.5 R200mlogit为显著的负值模型高度确信这些是团星系在约0.5到1.5倍R200m之间logit快速上升并穿过零点在更远的距离r 2 R200mlogit趋于一个稳定的正值对应场星系。这个S形转变区就是我们寻找的“环境影响过渡区”。3.3 收敛性测试结果是否依赖于标签定义这是一个必须回答的关键问题我们找到的r0会不会只是因为当初画“团”和“场”的标签圈比如0.5和5倍R200m而人为制造出来的为此我们进行了系统的收敛性测试。固定内边界为0.5倍R200m将外边界从1.0倍逐步增加到5.0倍R200m步长0.1用每一组标签训练模型并观察预测的r0如何变化。结果非常振奋人心如图4所示尽管用于定义“场星系”的外边界在不断外推但模型预测的过渡半径r0很快大约在外边界超过3倍R200m后就稳定在一个平台值上其波动范围远小于r0本身的误差棒。这意味着只要用于训练的“团”和“场”样本在属性空间上是足够纯净和分离的模型就能稳健地学习到两者之间的本质区别而这个区别所对应的空间过渡位置是独立于标签的精确距离阈值的。这强有力地证明了r0反映的是星系属性随环境变化的真实物理信号而非我们人为标注的产物。4. 核心发现与物理机制解读4.1 环境影响的空间范围超越R200我们的主要发现如图5左、中面板所示。机器学习预测的过渡半径r0其典型值在0.8 到 1.2倍 R200m之间并且随着宿主星系团质量的增加而轻微增大。这个结果有两层重要含义星系团的影响范围比传统边界更远平均而言星系团环境对星系属性的可观测影响一直延伸到约1.2倍R200m处。这意味着有相当一部分在传统定义上处于“星系团外”的星系其实已经在物理上受到了星系团环境的显著影响。R200m作为一个基于球对称过密度的几何定义并不能作为环境效应的“硬边界”。与 splashback 半径的对比图中我们还叠加了理论预言的 splashback 半径由Diemer 2020和More et al. 2015给出。Splashback半径是暗物质粒子轨道最远点的中位数位置是一个动力学的边界。有趣的是splashback半径随质量减小而我们的r0随质量增大趋势完全相反。这说明星系属性发生转变的位置并非由动力学边界直接决定。它标记的是环境物理过程如气体剥离对星系产生显著影响的区域这个区域与动力学边界相关但并不重合。4.2 本征散射模糊的过渡区图5中的垂直误差棒很大这并非测量误差而是本征散射——它反映了真实的物理多样性。有些星系团的过渡发生得靠内有些则靠外即使在同一个团内不同星系因其轨道、倾角、固有质量等因素开始被“预处理”的位置也不同。图6通过分析不同logit区间的星系径向分布生动地展示了这一点高置信团星系(logit -5紫色曲线)高度集中在r 1倍R200m以内分布尖锐。高置信场星系(logit 5浅绿色曲线)几乎完全分布在r 1倍R200m以外。模糊样本(-0.23 logit 0.23青绿色曲线)分布非常弥散从团内部一直延伸到远郊。这些正是处于“正在转变中”的星系。这个模糊的、弥散的过渡区图像比一个尖锐的边界更符合物理直觉。它告诉我们星系团环境的影响是一个渐进的、概率性的过程。我们的模型成功捕捉到了这种复杂性其输出logit本身就可以作为一个数据驱动的、连续的环境“影强度”指标。4.3 质量依赖性的物理根源浓度-质量关系为什么大质量星系团的过渡半径r0反而更大这背后隐藏着深刻的物理联系。我们发现r0与星系团质量M200m之间存在幂律关系r0 ∝ M200m^α拟合指数α ≈ 0.10见图7。这个指数与宇宙学模拟中建立的暗物质晕浓度-质量关系的指数惊人地一致c(M) ∝ M^{-β}其中β ≈ 0.10。这绝非巧合。低质量暗物质晕形成更早宇宙背景密度更高因此其物质分布更“紧凑”即浓度c更高。高浓度意味着密度剖面在内部更陡峭向外下降得更快。在经典的冲压剥离框架中剥离条件取决于星系所处位置的星系团内介质密度ρ_ICM(r)。对于一个高浓度陡峭密度轮廓的晕ρ_ICM(r)在较小的r/R200处就下降到不足以有效剥离星系气体的阈值。反之对于低浓度平缓轮廓的大质量晕其气体密度在外围下降更慢使得剥离过程可以发生在更大的r/R200处。因此机器学习模型从星系属性中学到的过渡半径r0本质上反比于晕的浓度r0 ∝ 1/c ∝ M^{β}。我们的测量α ≈ β ≈ 0.10完美验证了这一推论。这表明模型并非在学习一个任意的模式而是直接探测到了底层引力势的物理结构通过浓度如何设定了环境过程生效的空间尺度。4.4 不同物理过程的响应尺度我们训练了仅使用特定类别特征的模型如仅用恒星属性、仅用气体属性、仅用动力学属性以探究不同物理过程主导的转变发生在何处。图5右面板和图9展示了不同特征组合预测的r0之比。结果显示动力学属性如速度、自旋定义的过渡区最靠内。这是因为动力学状态如轨道、速度弥散是星系与团引力势直接相互作用的结果在星系深入势阱内部时变化最为剧烈。气体属性和恒星属性定义的过渡区相对更靠外且彼此接近难以完全区分。这与观测相符气体尤其是冷气体作为星系中最“脆弱”的组分会最先受到环境的影响如冲压剥离、潮汐剥离导致恒星形成率下降sSFR变化和颜色变红这些是相对较早的信号。而恒星质量、金属丰度等属性变化可能需要更长时间。这个排序为我们理解星系“被 quenching”停止形成恒星的序列提供了线索可能先是在外围经历气体剥离影响气体属性然后随着深入恒星形成停止影响sSFR、颜色最后动力学状态被完全改造。4.5 与传统边界定义的对比我们同样尝试了直接用传统边界R200m R200c Rsp作为二分类标签即内部为团外部为场来训练模型。结果如图8所示这些模型预测的r0与它们用作标签的边界值并不一致且误差棒巨大。这清晰地表明一个固定的几何或动力学边界并不能完美对应星系属性发生系统性变化的区域。星系在穿越splashback半径之前其属性可能已经开始变化预处理而在R200之内也可能存在一些尚未被完全影响的星系。我们的概率化、数据驱动的方法正是为了捕捉这种物理上的模糊性和渐进性。5. 实操经验、挑战与未来方向5.1 数据准备与特征处理的坑模拟数据的异构性不同宇宙学模拟如IllustrisTNG, EAGLE, SIMBA输出的星系属性其定义、单位、甚至计算方式可能有细微差别。在合并使用或与观测对比前必须进行极其仔细的交叉验证和归一化。我建议为每个模拟数据集单独编写数据读取和预处理管道并保存一份详细的“数据护照”记录每个特征的出处和变换历史。缺失值与异常值模拟数据中也可能存在极端值或无效值例如某些矮星系的气体质量接近零导致气体分数计算溢出。arcsinh变换虽然能处理零值但对于极端负值或由于数值误差产生的异常值仍需警惕。稳健的做法是先进行描述性统计如绘制直方图根据物理意义设定合理的截断范围再进行变换。训练/验证/测试集划分绝不能随机划分因为星系是成团的随机划分会导致同一个星系团的成员同时出现在训练集和测试集造成数据泄露严重高估模型性能。必须按星系团ID进行划分确保用于训练、验证、测试的星系来自完全不同的、独立的星系团集合。5.2 模型训练与调参心得学习率是命门对于宇宙学数据学习率宁小勿大。我从5e-5开始配合余弦退火在大多数任务上都能稳定收敛。如果使用更大的批量大小如512或1024可以适当增大学习率但需密切监控验证集损失。早停的耐心值设置早停的“耐心”参数需要平衡。太短如3可能导致在损失平台期提前停止错过后续下降太长如20则浪费计算资源。我发现在损失曲线进入明显的“平原期”后再观察5-10个epoch是合理的通常设置耐心为5-10。类别权重的计算加权交叉熵中的权重w N_neg / N_pos需要谨慎计算。这里的N_neg和N_pos应该是训练集中的数量而不是全集。并且如果正负样本比例极度失衡如超过1:10可能需要更激进的重采样策略如对少数类过采样或使用Focal Loss等更高级的损失函数。5.3 结果解释的注意事项相关性不等于因果模型找到了属性与环境的关联并预测了一个过渡半径。但这并不意味着在这个半径上“发生”了某个特定事件如冲压剥离启动。它只是表明在这个统计意义上星系团的集体影响开始显著地改变星系属性的整体分布。解释时需要结合具体的物理模型。外推需谨慎我们的模型是在特定模拟如TNG300上训练的。将其直接应用于观测或其他模拟时必须考虑域适应问题。观测数据存在选择效应、测量误差和不同的物理过程如AGN反馈模型不同。一个实用的做法是在目标数据集上对模型进行少量样本的微调或者使用域适应技术。“黑箱”的可解释性虽然我们通过分析不同特征组合的输出来推断物理但深度神经网络的具体决策路径仍是复杂的。可以辅以SHAP或LIME等可解释性AI工具量化每个输入特征对单个星系预测结果的贡献度这能提供更细致的物理洞察。5.4 未来扩展方向引入时间维度当前工作是静态快照。下一步可以引入星系的轨道信息如最近通过 pericenter 的时间或形成历史训练模型预测星系“被影响”的阶段或未来演化轨迹实现更动态的环境诊断。结合空间信息当前模型只用了星系自身的属性。可以尝试图神经网络将星系的局部环境如最近邻星系的距离、所属纤维结构的形态也作为输入捕捉环境效应的非局部性。应用于大规模巡天数据将训练好的模型应用于SDSS、DESI、Euclid等巡天数据。最大的挑战是处理观测误差和选择函数。一个可行的路径是用高保真模拟生成带噪声的、符合观测限制的“模拟观测”数据来训练模型再应用到真实数据。探索更复杂的架构对于更高维的特征集如积分场光谱数据衍生的多维度信息可以尝试注意力机制或Transformer架构让模型自动关注最重要的特征交互。这项研究最让我兴奋的一点是它展示了机器学习不仅仅是一个“黑箱”分类工具。当与深厚的物理洞察结合时它能够揭示出数据中隐藏的、符合物理规律的标度关系和内在联系。我们找到的r0 ∝ M^{0.1}关系就是这种结合产出的“物理可解释”的结晶。它告诉我们在追求预测精度的同时设计合理的实验、提出正的问题、并深入理解模型输出背后的物理才是将AI真正转化为天体物理学发现引擎的关键。
