宇宙大爆炸后的三分钟：氦元素如何成为检验引力理论的精密探针？
当爱因斯坦引力不再“够用”：f(R,G,T)引力引入新的作用项
BBN冻结温度之变：模型参数如何扰动早期宇宙的膨胀历史？
四个典型模型逐一验证：哪些参数区间与观测数据相容？
BBN一测就露馅？f(R,G,T)四模型被氦丰度卡住
早期宇宙最狠的地方，不是“想改就改”，而是连氦元素都能把你揪出来。本文用 BBN 和原始氦丰度，直接给扩展引力四个模型做体检，参数一旦飘过头，立刻露馅。
原论文信息如下：
宇宙早期最“诚实”的证据，不是遥远星系，也不是复杂的黑洞图像，而是氦、氘、氦-3、锂-7这些轻元素的原始丰度。原因很简单：它们是在大爆炸后最初几分钟里被“烤”出来的，任何一点点膨胀速率变化、弱相互作用冻结时机变化，都会在最后的元素比例上留下痕迹。换句话说，BBN（Big Bang Nucleosynthesis，大爆炸核合成）就像宇宙给自己留的验血报告，谁想改早期引力，先过它这一关。这篇工作最有意思的地方，不是又写了一个“更一般”的引力函数，而是把理论直接拉到最难糊弄的早期宇宙场景里做约束。BBN对宇宙学模型的脾气很大：只要哈勃膨胀率稍微快一点，弱反应来不及把中子转回质子，中子就会被更多地“保留下来”，最后氦-4 的质量分数就会偏高；反过来，膨胀慢一点，氦又会少一点。氦丰度对扩展历史极其敏感，所以它特别适合拿来卡参数。更关键的是，BBN 的观测数据已经相当精确，原始氦丰度的不确定度已经压到了千分之一的量级，这意味着任何试图修改早期宇宙膨胀历史的引力理论，都必须在这个精度下接受检验。论文选择的 f(R,G,𝒯) 框架，正是当前扩展引力研究中一个非常活跃的方向，因为它同时包含了曲率修正、高阶拓扑项和物质耦合，覆盖面很广。通过 BBN 这把“手术刀”，作者们得以精确解剖这些扩展项在早期宇宙中的具体影响，并给出参数空间的“生死线”。
图1：封面图。本文的核心思路很直接：把扩展引力模型扔进大爆炸核合成的“审讯室”，看它会不会把原始氦丰度弄得太离谱。这个思路的巧妙之处在于，它不需要复杂的数值模拟，而是利用 BBN 中成熟的解析关系，将引力模型的参数直接映射到可观测的氦丰度上，从而得到清晰且严格的约束条件。
在深入技术细节之前，有必要先厘清 BBN 的基本物理图像。大爆炸后约 1 秒到 3 分钟，宇宙的温度从约 100 亿度下降到约 10 亿度，这个过程中质子和中子开始结合成原子核。最先形成的氘核（一个质子加一个中子）非常脆弱，容易被高能光子打碎，直到宇宙冷却到足够低的温度，氘核才能稳定存在，并迅速通过一系列核反应链合成氦-4。氦-4 的最终产量，几乎完全由 BBN 开始时中子和质子的相对数量决定。而这个比例，又由弱相互作用（将质子和中子相互转换）的“冻结”时刻决定。因此，整个 BBN 链条的起点，就是弱相互作用的冻结温度 Tf。任何修改早期宇宙膨胀速率的引力理论，都会通过改变 H(T) 来改变 Tf，进而改变中子-质子比，最终在氦-4 丰度上留下可观测的信号。这就是为什么原始氦丰度能成为检验引力理论的精密探针——它放大了引力模型在早期宇宙中的微小偏差。
先把名字拆开说人话。这里的 f(R,G,𝒯) 引力，就是把标准广义相对论里的爱因斯坦-希尔伯特作用量，扩展成一个同时依赖三个量的函数：Ricci 标量 R、Gauss-Bonnet 不变量 G，以及能量动量张量的迹 𝒯。R 负责普通曲率，G 负责更高阶的几何修正，𝒯 则把“物质本身”也拽进了引力方程里。这个思路的野心很明显：不是只在几何上做文章，而是直接让物质和几何发生更复杂的耦合。为什么需要这样的扩展？因为标准广义相对论虽然在太阳系尺度上通过了大量检验，但在解释宇宙晚期加速膨胀、暗物质和暗能量等观测现象时，却不得不引入神秘的暗能量成分。f(R,G,𝒯) 这类扩展引力理论试图通过修改引力本身的行为，来统一解释这些现象，从而避免引入新的、尚未被直接探测到的物质形式。当然，这种扩展并非没有代价——它引入了更多的自由参数，而这些参数必须受到各种宇宙学观测的严格约束，BBN 就是其中最有力的约束之一。
公式1：f(R,G,𝒯)引力的作用量。这里的 f(R,G,𝒯) 是一个任意函数，Lm 是物质拉格朗日密度，κ2=8πGN。这类写法的关键不是“更复杂”，而是它允许引力响应同时依赖曲率、拓扑型高阶项和物质内容。当 f(R,G,𝒯) = R + 2Λ（Λ为宇宙学常数）时，它就退化为标准的广义相对论加宇宙学常数。因此，这个框架可以看作是广义相对论的一个自然推广，而论文的工作就是通过 BBN 来探测这个推广的“边界”在哪里。
公式2：R、G 与 𝒯 的定义。R 是 Ricci 标量，G 是 Gauss-Bonnet 不变量，𝒯 是能量动量张量的迹。对非专业读者来说，可以把它理解成“几何三件套”：普通曲率、二阶修正、以及物质信息。Gauss-Bonnet 项在四维时空中是一个拓扑不变量，但在修改引力理论中，当它被赋予一个耦合函数时，就会对动力学产生实质影响。𝒯 的引入则更为直接，它意味着引力场的方程不仅取决于时空的弯曲，还取决于物质场的具体性质（如压强和能量密度），这为物质与几何的相互作用提供了新的通道。
论文采用的是平直的 FLRW 宇宙背景，也就是早期宇宙常见的均匀各向同性模型。这个背景下，R 和 G 都能写成哈勃参数 H 及其时间导数的函数。翻译成人话就是：一旦宇宙膨胀速度变了，几何项立刻跟着变，BBN 就能感受到。具体来说，在 FLRW 度规下，Ricci 标量 R 正比于 Ḣ + 2H²，而 Gauss-Bonnet 不变量 G 则正比于 H²(Ḣ + H²)。这意味着，任何修改 H 演化规律的引力项，都会直接反映在 R 和 G 上，进而通过场方程反过来影响 H 的演化。这种自洽的耦合关系，使得 BBN 能够成为一个极其灵敏的探测器——它捕捉的是整个动力学系统的连锁反应。
公式3：FLRW 度规。a(t) 是尺度因子，决定宇宙“长大”的速度。这里用这个背景不是为了显摆数学，而是为了把早期宇宙的膨胀写成可计算的形式。在 BBN 时期，宇宙主要由辐射（光子、中微子、电子-正电子对）主导，其状态方程参数 ω = 1/3，这使得 FLRW 度规下的动力学方程大大简化，为解析求解提供了可能。
公式4：FLRW 背景下的 R 与 G。可以看到，G 不只是“另一个曲率项”，它还显式依赖 H2(Ḣ+H2)，因此对膨胀历史的变化非常敏感。在辐射主导时期，标准广义相对论给出 H ∝ T²，因此 G 的行为也会随之确定。一旦引入 f(R,G,𝒯) 修正，H 的演化偏离标准形式，G 的偏离就会被放大，从而在 BBN 的观测中留下清晰的印记。
公式5：能量动量张量的迹 𝒯=3p−ρ=(3ω−1)ρ。这里 ω 是状态方程参数。把 𝒯 放进引力函数里，本质上就是允许“物质的性质”反过来影响引力方程。对于辐射 (ω=1/3)，𝒯=0；对于非相对论性物质 (ω≈0)，𝒯=-ρ。这意味着，𝒯 项在宇宙的不同演化阶段会扮演不同的角色，从而为模型提供丰富的动力学行为。论文中选取的四个模型，正是通过不同的方式组合 R、G 和 𝒯，来探索这种丰富的可能性。
公式6：连续性方程。它表达的是能量守恒形式下的密度演化。早期宇宙里，ρ 和 p 的变化直接决定膨胀与冷却的节奏。在 f(R,G,𝒯) 理论中，由于物质与几何的耦合，连续性方程本身也可能被修改，但论文中假设标准连续性方程仍然成立，这相当于对模型施加了一个额外的约束，使得分析更加聚焦于引力部分的修正。
如果只看公式，容易觉得这类模型很“理论味”。但它真正的用途很朴素：把标准宇宙学里不够灵活的地方打开一点口子，然后用观测去把口子重新关小。说白了，理论可以大胆，数据负责收拾残局。论文的整个论证过程，就是一次精彩的“收拾残局”示范：它没有试图去证明哪个模型是“正确”的，而是精确地划定了哪些参数区间是“仍然可能”的。这种“证伪”而非“证实”的思路，正是科学方法论在宇宙学前沿研究中的典型体现。
BBN 的主角不是氦，而是冻结温度 Tf。在宇宙温度足够高的时候，弱相互作用足够快，质子和中子可以不停互相转换；当宇宙膨胀冷却到某个临界点，反应速率追不上膨胀速率，弱反应“掉线”，中子-质子比例就被锁死了。这个锁死时刻一旦提前，中子剩得更多，氦就会更多；一旦推后，中子衰变更多，氦就会更少。因此，Tf 是整个 BBN 过程的“节拍器”，它决定了后续所有核合成反应的初始条件。论文的核心工作，就是计算在 f(R,G,𝒯) 引力下，Tf 如何偏离标准值，并将这个偏离与可观测的氦丰度变化联系起来。
公式7：标准辐射主导时期的哈勃率。BBN 发生在辐射主导阶段，因此这里先拿标准 GR 的 H2 作为基准。ρr 是辐射能量密度，MP 是普朗克质量。这个公式表明，在标准宇宙学中，膨胀速率完全由能量密度决定。而在扩展引力中，H2 的表达式中会多出一些由 f(R,G,𝒯) 项贡献的“有效能量密度”项，从而改变膨胀的快慢。
公式8：辐射能量密度 ρr=π2g*T4/30。g* 是有效自由度，早期宇宙里常取约 10 左右。温度越高，辐射越强，宇宙越像一锅滚开的“粒子汤”。在 BBN 时期，g* 主要包括光子、电子、正电子和三种中微子，其具体数值会随着温度变化，因为当温度低于某些粒子的质量时，这些粒子会湮灭，不再贡献自由度。论文中为了简化计算，将 g* 取为一个常数近似值，这在精度要求内是合理的。
公式9：时间与温度的关系。这个式子说明，宇宙越晚期，温度越低，且大体满足 t-1/2 的冷却规律。BBN 的很多约束都来自这里：温度降得快不快，直接影响冻结时机。这个关系是由辐射主导时期的弗里德曼方程和能量守恒方程共同推导出来的，它建立了宇宙年龄和温度之间的直接对应，是连接动力学和热力学的重要桥梁。
这里的关键不是把公式堆满屏，而是把两条时间尺度放在一起比较：一条是宇宙膨胀的时间尺度，一条是弱相互作用的反应时间尺度。前者更快，后者来不及，冻结就发生了。论文后面所有 BBN 约束，实际上都在围绕这个“谁先掉链子”展开。膨胀的时间尺度由哈勃参数 H 的倒数给出，而弱相互作用的反应时间尺度由总反应率 Λ 的倒数给出。冻结条件 H(Tf) ≈ Λ(Tf) 正是这两条时间尺度相等的数学表达。在扩展引力中，H(T) 被修改，因此满足冻结条件的 Tf 也随之改变。
公式10：总反应率由质子转中子和中子转质子两部分组成。这个总率越大，弱反应越能跟上宇宙膨胀；越小，就越容易提前冻结。在早期宇宙的高温高密环境下，这些弱反应主要包括 n + νe ↔ p + e- 和 n + e+ ↔ p + ν̅e 等过程。它们的反应率可以通过量子场论精确计算，是 BBN 理论中已知度最高的部分之一。
公式11：总反应率的近似表达式。Q 是中子-质子质量差，A 是常数。这个式子本身不神秘，真正有用的是它和 H(T) 一起决定了冻结温度。这个近似表达式在 BBN 相关的温度范围内（约 0.7-1 MeV）有很好的精度，它抓住了反应率随温度幂律增长的核心特征，使得冻结条件的求解变得简洁明了。
公式12：Q=mn−mp，约为 1.29×10−3 GeV。中子比质子重一点，所以在热平衡里会更少一些。这个微小的质量差，是决定 BBN 最终产物比例的关键因素之一。在热平衡状态下，中子与质子的数量比由玻尔兹曼因子 exp(-Q/T) 决定，温度越高，两者比例越接近 1:1；温度降低，中子比例指数下降。
公式13：冻结时刻的中子-质子比值。这里的指数项说明，中子比例和温度、寿命都有关，别看只是一个指数，最后能把氦丰度拧出明显变化。这个比值是冻结温度 Tf 的函数，一旦 Tf 确定，中子-质子比就基本确定，后续的氦-4 产量也就随之确定。因此，整个 BBN 对引力模型的约束，最终都可以归结为对 Tf 的约束。
公式14：冻结条件 H(Tf)≈Λ(Tf)≈cqTf5。这就是 BBN 里最重要的门槛条件：膨胀率和反应率一旦相等，冻结就开始了。这个条件将宇宙的膨胀动力学（由引力理论决定）与粒子物理的反应率（由弱相互作用理论决定）联系在了一起，是连接宏观宇宙学和微观粒子物理的桥梁。
公式15：冻结温度 Tf 的解析表达式。它告诉人们，冻结温度受 g*、cq 和普朗克尺度控制。只要引力模型改变了 H(T)，Tf 就会跟着偏移。在标准宇宙学中，代入具体数值可以得到 Tf ≈ 0.8 MeV。而在扩展引力中，Tf 会偏离这个值，偏离的大小和方向取决于具体模型和参数。
论文接着做了一步很关键的近似：把修改引力带来的额外能量密度，视作对标准辐射背景的小修正。这样就能把 Tf 的变化写成一个很直观的相对偏差 |ΔTf/Tf|。这招很实用，因为观测上给出的约束本来就是一个小量，直接拿相对偏差去比，最方便。这个近似成立的前提是，扩展引力的修正项在 BBN 时期远小于标准的辐射能量密度。如果修正项太大，那么整个微扰论的处理方式就会失效，需要直接求解完整的场方程。论文选取的四个模型，其参数范围恰好保证了这种微扰处理的有效性，这也是模型选择的一个隐含条件。
公式16：氦丰度修正与冻结温度偏移的关系。这个式子把“宇宙膨胀变一点点”翻译成“氦丰度变多少”。别小看这个翻译器，它就是 BBN 约束的核心桥梁。这个关系式是通过对 BBN 核反应网络进行线性化分析得到的，它表明氦丰度的相对变化与冻结温度的相对变化成正比，比例系数约为 0.04。这意味着，Tf 每变化 1%，Yp 就会变化约 0.04%。考虑到 Yp 的观测精度已经达到 0.1% 量级，Tf 的可允许变化范围就被限制在了约 2.5% 以内，这是一个非常严格的约束。
公式17：观测给出的氦丰度范围 Yp=0.245±0.003，且 |ΔYp|<10−4。这组数值就是“别乱来”的硬边界。这个观测值是通过对低金属丰度的河外星系 HII 区进行光谱分析得到的，经过了复杂的系统误差修正，是目前宇宙学中最精确的原始元素丰度测量结果之一。它就像一把经过精密校准的卡尺，用来测量早期宇宙的膨胀历史。
公式18：|ΔTf/Tf| 的近似式。只要模型带来的额外密度 ρ 不能太大，冻结温度偏移就能被压在观测上限之内。这个近似式将复杂的场方程积分简化为一个代数关系，使得我们可以快速评估不同模型参数对 Tf 的影响，是论文进行参数扫描和约束的核心工具。
公式19：最终用于约束模型的观测上限。也就是说，理论可以很花，但最后必须老老实实落在 4.7×10−4 这个量级下面。这个数值是结合了 Yp 的观测误差和理论计算的不确定性后得出的一个保守上限，它代表了 BBN 对早期宇宙膨胀速率变化的容忍度。任何导致 |ΔTf/Tf| 超过这个上限的引力模型参数，都会被 BBN 观测所排除。
论文没有只拿一个模型“试水”，而是选了四个代表性的 f(R,G,𝒯) 形式来逐个过筛。这个做法比较务实：如果一个理论框架真的有生命力，不该只在某个特别挑剔的参数点成立，而应该在一批典型构型里都能给出合理区间。四个模型分别把 R、G、𝒯 以不同方式组合，看看谁更容易被 BBN 接受，谁更容易被一票否决。这种“广撒网”的策略，有助于揭示 f(R,G,𝒯) 框架下不同耦合方式对早期宇宙动力学的普遍影响，而不是局限于某一个特例。通过比较四个模型的约束结果，可以更全面地理解 BBN 作为引力理论检验工具的威力。
模型1把 R 设为线性项，G 设为幂律项，𝒯 设为线性耦合。它的结构最像“标准引力 + 一点高阶修正 + 一点物质耦合”，属于比较均衡的试验款。这个模型可以看作是 f(R) 引力和 f(G) 引力的一个简单混合，并额外加入了物质耦合项。它的参数空间包括 α1（R 项的系数，通常取 1 以恢复广义相对论的低能极限）、n（G 项的幂指数）和 γ1（𝒯 项的耦合常数）。论文的目标就是通过 BBN 约束来确定 n 和 γ1 的允许范围。
模型1里额外能量密度 ρ 的表达式已经很说明问题：只要 n 改变，G 的高阶贡献就会被放大或压缩，最后体现在 H 的偏移上。这里的参数 α1、γ1、n 都会影响冻结温度。从表达式可以看出，额外能量密度正比于 H4n，这意味着 n 越大，额外项对 H 的依赖越强，在早期宇宙 H 很大的时候，其影响会被急剧放大。因此，BBN 对 n 的约束会非常敏感，这正是图 2 所展示的结果。
模型1通过把今天的宇宙学参数 H0、ΩDE0 代回去，反推出 γ1 的可行范围。这个步骤很关键，因为它不是纯符号游戏，而是把模型拉回到现实宇宙。通过要求模型在今天的宇宙中也能给出合理的暗能量行为（即 ΩDE0 ≈ 0.7），可以确定 γ1 与 α1 和 n 的关系，从而减少自由参数的数量，使 BBN 的约束更加聚焦。这种“高低能联动”的策略，是宇宙学模型检验中的常用手法。
模型1的 |ΔTf/Tf| 最终写成了 n 的函数。这个式子说明，约束的主战场其实是 n：幂次越高，高阶几何项越容易把偏差推大。从函数形式看，|ΔTf/Tf| 随 n 指数增长，因此存在一个临界值 nc，当 n > nc 时，模型就会被 BBN 排除。论文通过数值计算给出了这个临界值 nc ≈ 0.3721，这是一个非常精确的定量结果。
图2：模型1中冻结温度相对偏差 |ΔTf/Tf| 随参数 n 的变化。红虚线是观测上限，阴影区是允许区域。结论很直白：当 n 增大到约 0.3721 之后，模型开始越线。这个图清晰地展示了 BBN 约束的“一刀切”特性：在 n 小于临界值时，模型与观测完全兼容；一旦超过，就立刻被排除。这种清晰的界限，使得 BBN 成为筛选扩展引力模型的一个高效工具。
模型2换了一个耦合方式，把 R 和 𝒯 直接绑在一起，再让 G 作为幂律修正项。这个结构更像“曲率和物质绑一块，高阶几何单独加戏”。从 BBN 的角度看，这种写法通常更容易把冻结温度推偏，因为物质耦合会直接影响有效引力响应。模型2的形式为 f(R,G,𝒯) = α2R𝒯 + γ2Gn，其中 R𝒯 项意味着时空曲率与物质能量动量张量的迹直接相乘，这是一种比模型1更强的耦合方式。这种耦合在宇宙早期辐射主导时期（𝒯=0）会自动消失，但在物质主导时期（𝒯≠0）会产生显著影响，因此它对宇宙不同演化阶段的区分能力更强。
模型2的核心是 R𝒯 耦合。它比模型1更“激进”，因为曲率和物质不再只是并排出现，而是直接交互。这种交互意味着，物质的分布和运动可以直接改变时空曲率的动力学行为，反之亦然。在 BBN 时期，虽然 𝒯=0 使得 R𝒯 项不直接贡献，但它会通过影响宇宙的整体动力学方程（如弗里德曼方程）来间接改变 H 的演化，从而影响 Tf。
模型2中 ρ 的形式说明，γ2 和 α2 的配合会直接改变 H 的修正强度。只要耦合项稍微放大，BBN 就能感知到。与模型1类似，模型2的额外能量密度也正比于 H4n，因此 n 仍然是决定 BBN 约束强度的关键参数。但不同的是，模型2中还有一个来自 R𝒯 项的贡献，它使得整个约束条件对参数的选择更加敏感。
模型2把 α2 也压到了当前宇宙学参数下进行约束。这样的处理避免了“参数彼此乱飞”的问题，至少保证模型不是纸上谈兵。通过要求模型在今天的宇宙中给出正确的暗能量密度，可以确定 α2 与 γ2 和 n 的关系，从而将自由参数减少到一个（n），使得 BBN 的约束更加清晰和直接。
模型2的 |ΔTf/Tf| 仍然主要受 n 控制，但与模型1相比，耦合结构更复杂，允许区间通常更容易被挤压。从图3可以看出，模型2的 Yp 随 n 的变化曲线比模型1更陡峭，这意味着在相同的 n 值下，模型2对氦丰度的扰动更大，因此其可接受的 n 范围也更窄。
图3：模型2中原始氦丰度 Yp 随参数 n 的变化。曲线越过观测带时，就意味着这个 n 已经不太老实了。BBN 在这里不是“建议”，而是直接判定。与模型1相比，模型2的曲线在更小的 n 值处就穿过了观测带的上限，这表明 R𝒯 耦合确实使得模型对 BBN 约束更加敏感。
模型3更像一个“反过来写”的版本：把 R𝒯 放到分子里，再除以 Gn。这类构造的味道很明显——想让高阶几何项在低 G 区域更敏感。问题也很明显：一旦分母里的 Gn 太强，早期宇宙的修正可能会非常剧烈，所以它往往是最容易被 BBN 盯上的一类。模型3的形式为 f(R,G,𝒯) = α3R𝒯 / Gn + γ3，这种形式在 G 很小（即曲率变化平缓）的区域会产生很大的修正，而在 G 很大（早期宇宙）的区域，修正反而被压制。这种“反直觉”的设计，使得模型3的动力学行为与模型1和2截然不同。
模型3把耦合和高阶项叠得更深。它的自由度虽然多，但也意味着被观测“掐住”的地方更多。由于分母中有 Gn，模型3的额外能量密度表达式会变得更加复杂，它不仅依赖于 H，还可能依赖于 Ḣ，这使得对 Tf 的分析需要更加小心。
模型3中的 ρ 结构更复杂，说明冻结温度偏移会同时受 γ3、α3 和 n 影响。复杂不等于强，复杂有时候只是更容易被数据拆穿。从图4可以看出，模型3的 Yp 随 n 的变化曲线呈现出与模型1和2不同的趋势，其可接受的 n 范围可能更窄，或者对 n 的符号有特殊要求。
模型3通过 α3 的约束把当前宇宙参数重新嵌回模型。这个步骤的价值在于：让模型不是只在抽象符号上成立，而是在今天的宇宙背景下也能自洽。与模型1和2类似，通过要求模型在今天的宇宙中给出正确的暗能量行为，可以确定 α3 与 γ3 和 n 的关系，从而简化 BBN 的分析。
模型3的冻结温度偏差表达式很长，但本质还是一句话：参数一旦让额外修正过大，BBN 就会立刻报警。尽管表达式复杂，但其核心逻辑与模型1和2一致：通过比较 |ΔTf/Tf| 与观测上限，来确定参数 n 的允许范围。
图4：模型3中原始氦丰度 Yp 随参数 n 的变化。相比前两个模型，这一类构造通常更容易被观测上限卡住，说明高阶分母项并不好“蒙混过关”。从图中可以看出，模型3的 Yp 曲线在 n 的某些取值区间内非常陡峭，这意味着模型对 n 的微小变化极其敏感，其可接受的参数空间可能非常有限。
模型4则换成了另一种混合方式：线性 R 项加上 𝒯/Gn，再加一个常数项。它看起来像“把物质耦合和高阶修正打包一起”，实际效果是把约束条件变得更分散。这样的模型有时会比前几个更难调，因为每个项的角色都不一样，谁抢戏太多，BBN 都会看出来。模型4的形式为 f(R,G,𝒯) = α4R + β1𝒯/Gn + γ4，它结合了模型1的线性 R 项和模型3的 𝒯/Gn 结构，是一种混合型模型。这种设计试图在保持广义相对论线性 R 项优点的同时，引入物质与高阶几何的复杂耦合。
模型4把 𝒯 和 Gn 绑在一起，说明物质耦合的强弱会受到高阶几何项调制。在早期宇宙，G 很大，因此 𝒯/Gn 项被压制，模型近似回到广义相对论；而在晚期宇宙，G 变小，该项的影响逐渐显现，可能驱动宇宙的加速膨胀。这种“自动开关”机制是模型4的一个有趣特点。
模型4的 ρ 表达式显示，β1、α4、γ4 与 n 的联合作用会决定偏离程度。这里的约束比“看起来像”更重要，模型是否可行，最终还是看它能不能把偏差压住。由于模型4包含更多的自由参数，其 BBN 约束的分析也更为复杂，但最终仍然可以归结为对 n 等关键参数的约束。
模型4中 β1 的表达式同样把今天的宇宙学参数代入，从而得到可接受范围。这个做法让四个模型都站在同一把尺子上比较。通过要求模型在今天的宇宙中给出正确的暗能量密度，可以确定 β1 与 α4、γ4 和 n 的关系，从而将自由参数减少，使得 BBN 的约束更加聚焦。
模型4的冻结温度偏差依然可以写成参数函数。虽然形式不一样，但逻辑没变：只要模型让宇宙膨胀偏离得太多，BBN 就会把它踢出局。从图5的对比可以看出，模型4的约束结果介于模型1和模型3之间，体现了其混合型模型的特点。
LONGGE AI COMMUNITY
把每天读到的论文，变成长期积累
加入「龙哥读论文」知识星球，持续获取 AI 论文、资讯、开源项目、招聘与研究思路。
加入龙哥读论文微信群：添加微信 kangjinlonghelper，备注“研究方向 + 地点 + 学校/公司 + 昵称”。
微信扫码加入知识星球