折射引力用依赖密度的介电常数取代暗物质
折射引力（RG）是一种现象学修改引力模型。它不引入暗物质粒子，而是通过引入一个依赖重子物质密度ρbar的有效引力介电常数ε(ρbar)来修改泊松方程。当密度低于某个临界值时，该介电常数会增强引力场，从而在宇宙结构中模拟暗物质的行为。
这里涉及三个关键参数。真空介电常数ε0决定了低密度区域中引力的渐近增强幅度。临界密度ρc标记转变发生的位置，指数Q控制转变的陡峭程度。该模型存在协变标量-张量表述，但本文检验的是现象学版本，并指出近期一项分析发现该特定协变模型无法在足够早的时期形成宇宙结构——因此RG仍是一个检验平台，而非已被确立的替代方案。
这关系到结果应如何解读。本文并非在追问RG整体上是否正确，而是在追问星系团运动学能把ε0约束到多紧，以及当重子质量模型改变时这一答案会移动多少。
将74个星系团按速度弥散堆叠成三组以抑制形状系统误差
样本为来自斯隆数字巡天星系团落入区（CIRS）的74个星系团。作者将它们按递增的速度弥散σ堆叠成三组，σ可作为总动力学质量的代理。堆叠是一种有意的权衡：它对偏离球对称和偏离动力学平衡的情形取平均，而这两项假设最可能使单星系团运动学拟合产生偏差。
分析使用MG-MAMPOSSt代码重建引力势，并拟合每个堆叠系统的速度各向异性剖面。只有投影半径R ≥ 0.05 Mpc且R ≤ 1.2 × 平均r200的成员星系进入拟合。内边界剔除动力学由最亮星系团星系主导的区域；外边界剔除可能已不处于平衡、或可能属于周围大尺度结构的星系。作者报告，将上限变动最多10%不会显著改变结果。
另有一组较小的堆叠仅由扰动较小的星系团构成，选取标准为A2 < 0.7。这一选择是本文的内部对照：它检验答案是否取决于系统的弛豫程度。
热气体斜率使ε0移动近一个数量级
星系团内介质主导重子质量预算，其径向剖面用beta模型参数化。作者用β = 2/3和β = 1各运行一次拟合，并对气体和恒星质量剖面参数进行边缘化，使重子不确定度传播到ε0。
对于三个主要堆叠，ε0的中心估计在β = 2/3时为0.04至0.08，在β = 1时为0.19至0.37。这就是核心跨度：同样的星系团、同样的运动学，推断出的真空介电常数相差约五倍，而驱动因素是假设的气体密度斜率。β = 1的情形偏好更大的ε0，在Bin 3中最为明显。
机制在文中表述明确。在远大于气体标度半径的半径处，beta模型密度按ρgas ∝ r^(−3β)下降。对于β = 2/3，密度按r^(−2)下降，累积气体质量随半径近似线性增长；对于β = 1，密度按r^(−3)下降，累积质量仅对数增长。由于RG的引力场完全由重子提供源，较浅的剖面在外围产生更大的引力贡献，拟合则通过移动推断的气体归一化、标度半径和ε0来补偿。
选择弛豫星系团缩小了ε0的气体模型偏移
由A2 < 0.7的星系团构成的两个堆叠在不同气体模型下给出远为一致的图景。它们的中心估计在β = 2/3时为ε0 = 0.04和0.06，在β = 1时为ε0 = 0.04和0.05。改变气体密度斜率所导致的偏移远小于三个主要堆叠。
这是所采用建模框架内的内部一致性结果，作者也如此表述：基于动力学状态的选择可能有助于降低系统不确定度。这并非声称弛豫星系团本质上是更好的引力探针。该比较同时改变了分组方式和重建的重子剖面，因此部分差异可能来自这些变化，而非仅来自弛豫程度。
实际解读是：星系团尺度RG约束对重子建模的敏感度本身依赖于样本。对同一理论、同一巡天的两次分析，可能因堆叠了哪些星系团而在多大程度上依赖气体参数化上有所不同。
牛顿引力拟合几乎不理会RG所敏感的气体斜率
本文首先在牛顿引力下用Navarro-Frenk-White暗物质剖面加独立的气体和星系成分拟合每个堆叠，作为基线。在那里，β = 2/3和β = 1的约束在不确定度范围内吻合良好，χ²差异可忽略。作者得出结论：在牛顿引力中，星系团内介质参数化带来的系统误差相对于统计不确定度是次要的。
与RG运行的对比才是关键。加入暗物质晕吸收了气体剖面原本会控制的许多信息，因此牛顿引力拟合对β基本不敏感。而在RG中，引力场仅由重子分布提供源，同样的β选择会直接传播到ε0。
本文还指出一个关于半径的微妙之处。MG-MAMPOSSt拟合出的r200值与表中列出的平均r200值并不完全一致，因为二者不是等价的估计量：表中数值是对单个基于焦散的测量的平均，而拟合确定的是对完整堆叠星系分布拟合的单一模型的有效动力学尺度。拟合出的r200应被理解为集合星系团的特征动力学半径，而非对平均值的重建。
ρc仅有下限而Q基本无约束
只有ε0被良好定位。转变指数Q完全无界，临界密度仅得到下限log10[ρc/(M⊙/Mpc³)] ≳ 12.3，在所有组中一致。RG参数彼此之间、以及与各向异性和数密度参数之间均无显著相关性，但真空介电常数与重子剖面参数之间存在轻微正相关。
本文解释了原因。运动学探测的大部分径向范围处于低密度区域ρbar ≪ ρc，在那里介电常数趋近ε(ρbar) ≈ ε0。在该区域，动力学约束的是引力场的渐近增强，但几乎不携带关于转变发生在何处或其陡峭程度的信息。排除最内侧50 kpc且未单独建模最亮星系团星系密度剖面，进一步降低了对转变区域的敏感度。
因此，本分析可以检验ε0在不同堆叠间是否一致，但无法检验ρc和Q是否具有普适性，作者也直接说明了这一点。
BIC偏好牛顿引力，但似然值相当
使用BIC = −2 lnLmax + k ln(N)，含NFW暗物质成分的牛顿模型在所有组中均以ΔBIC > 6显著优于仅含重子的RG模型。RG分析有10个自由参数（含重子参数），牛顿分析有9个，N为星系总数。
作者对这意味着什么保持谨慎。两个模型的最大似然值从未有显著差异，因此BIC偏好主要来自对额外参数的惩罚，而非更好的拟合。而且由于ρc和Q仅被弱识别，BIC差异的幅度应谨慎对待：例如将Q固定为3重做分析，两个框架之间没有统计偏好。
这就是模型比较的诚实总结。在当前精度下，数据并未在拟合质量上区分两个框架；它们只在参数经济性上区分，而且仅在弱约束参数保持自由时如此。
ε0与此前星系团和椭圆星系估计重叠，但与盘星系不符
推断的ε0值与已有的RG估计相当。它们与此前分析的两个大质量星系团的界限一致，在所列不确定度内与MACS 1206和Abell S1063均重叠，也与椭圆星系的最佳拟合真空介电常数一致。与盘星系的估计相比，结果显示出约2σ或更大的张力，其显著性取决于组和气体参数化。
本文提醒不要过度解读星系团之间的比较。此前的星系团约束来自两个单独的、大质量的、相对高红移的星系团，使用分辨的气体和恒星质量剖面以及最亮星系团星系的恒星速度弥散剖面。本分析使用较低红移的堆叠系统和由标度关系及各组平均性质重建的有效重子剖面。因此任何比较都混合了RG参数可能的固有变化与星系团质量、径向覆盖和重子建模的差异。
确实成立的一致性特指ε0。此前报告的最强张力涉及ρc和Q，而它们在此处约束不足，无法重新评估。
所有堆叠共有的重子假设无法被平均掉
主要限制在于重子剖面，以及将每个堆叠视为处于动力学平衡的单一球对称系统。气体分布由LX–Mgas标度关系间接推断，而该关系仅用样本中一小部分的分辨剖面校准；恒星成分则通过假设的普适星系与气体质量比与气体挂钩。这些假设为每个堆叠所共有，因此堆叠无法将它们平均掉。
最亮星系团星系和弥散星系团内光均未被显式建模。这一遗漏在RG中比在牛顿拟合中更重要，因为预测的加速度同时取决于包围重子质量Mbar(<r)和进入介电常数的局部密度ρbar(r)。作者还指出，需要模拟来量化堆叠究竟在多大程度上降低了星系团特有的系统误差，以及残余张力是来自重子和动力学建模，还是来自堆叠过程本身。
在当前精度下，本文认为简单的单beta气体参数化已足够，而带有额外弱约束参数的更灵活气体模型尚不必要。该结论特定于本数据集和径向覆盖。
优先事项是控制重子，而非增加星系团数量
本文最具决策相关性的结果不是ε0的数值，而是当气体密度斜率从β = 2/3变为β = 1时，ε0的中心估计从约0.04–0.08移动到约0.19–0.37，而对A2 < 0.7的堆叠这一偏移几乎消失。对于任何跨样本比较RG约束的人，这意味着气体参数化和动力学状态选择不是无关紧要的细节——它们决定了该理论主要参数误差棒的大小。
本文自身的前瞻也遵循同样的逻辑。它指向逐星系团分辨的气体密度和恒星质量剖面、对最亮星系团星系和星系团内光的显式处理，以及对光谱选择函数、成员判定、动力学状态和堆叠过程的独立验证。它指出4MOST eROSITA和4MOST CHANCES等巡天将提供远更多的星系团成员，但明确表示仅靠统计量增加无法控制主导的系统误差。
因此需要关注的阈值不是有多少星系团进入堆叠，而是重子分布是被测量还是被假设，以及堆叠样本是否按动力学状态选取。在这两点都改变之前，星系团尺度的RG约束将持续报告一个不确定度由建模选择而非数据主导的ε0。
本文引用的文献
Testing Refracted Gravity with the kinematics of stacked galaxy clusters
在74个CIRS星系团的三个速度弥散堆叠中，折射引力的真空介电常数ε0随热气体模型在0.04至0.37之间变化；各约束在不确定度内彼此一致，并与此前星系团尺度的估计重叠，而对扰动较小星系团（A2 < 0.7）构建的堆叠，气体模型敏感度缩小。