企业官网建设流程全解析

1. 模型构建：从标准模型到2HDM2S

粒子物理的标准模型（SM）无疑是二十世纪最成功的理论之一，它精确预言了从电磁相互作用到弱电统一，再到希格斯玻色子的发现等一系列现象。然而，站在今天的视角回望，标准模型更像是一座宏伟但未完工的大厦。它留下了几个根本性的谜题：宇宙中约85%的物质成分——暗物质——其粒子本质是什么？为何中微子拥有如此微小的质量？以及，宇宙早期物质与反物质的不对称性（重子不对称）所需的CP破坏源从何而来？这些“开放性问题”如同大厦图纸上缺失的关键房间，暗示着标准模型之外必然存在新的物理（BSM）。

一个自然而直接的扩展方向，就是希格斯扇区。标准模型中仅有一个希格斯二重态，但这更像是一种“最小化”的审美选择，而非基于第一性原理的强制要求。为什么标量场不能更多？引入额外的希格斯二重态，即双希格斯二重态模型（2HDM），是BSM物理中最受青睐的框架之一。它不仅能自然地提供新的CP破坏源，其丰富的标量粒子谱（如额外的CP偶、CP奇和带电希格斯玻色子）也为探索新物理打开了窗口。

然而，单纯的2HDM本身并不天然包含一个稳定的暗物质候选者。暗物质粒子需要满足几个关键特性：电中性、稳定（长寿命）、非重子性，并且其丰度需与宇宙学观测相符。一个经典的解决方案是引入一个在某种对称性下“惰性”的标量场，该对称性在真空下保持不破缺，从而保证最轻的惰性粒子是稳定的。这就是我们工作的起点：在Type-II 2HDM的基础上，引入两个实标量单态场S和P，并赋予它们一个新的Z‘2对称性（S→-S, P→-P）。这个模型被称为2HDM2S（Two-Higgs-Doublet Model with Two Real Singlets）。

为什么是两个单态，而不是一个？这背后有深刻的动机。首先，单态扩展的2HDM（如N2HDM，即引入一个实单态）已被广泛研究，但近年的研究表明，在满足普朗克卫星测得的暗物质遗迹密度和日益严格的直接探测实验（如LZ实验）约束下，其参数空间正面临巨大压力。增加一个单态，相当于在暗物质“工具箱”里多放了一件工具。两个单态之间的混合可以引入新的相互作用通道，改变暗物质粒子的湮灭截面和与核子的散射截面，从而可能开辟出尚未被实验排除的参数区域。其次，从理论自洽性角度看，两个单态的引入使得标量势的结构更加丰富，可能带来新的真空结构、相变行为，甚至与宇宙早期演化相关联。

我们的核心目标是构建一个完整、自洽的2HDM2S理论框架，并系统性地扫描其参数空间，寻找那些能同时满足所有理论约束（如势能有界性、真空稳定性、微扰幺正性）和实验约束（对撞机物理、暗物质遗迹密度、直接与间接探测）的“幸存”参数点。这不仅仅是一个参数拟合游戏，更是通过系统的“筛选”来理解：什么样的标量势结构、粒子质量谱和耦合强度，才能孕育出一个可行的暗物质候选者，并躲过当今所有实验的追捕。

2. 理论基石：标量势、对称性与粒子谱

构建任何粒子物理模型的第一步，都是写下其拉格朗日量中最重要的部分——标量势。对于2HDM2S模型，我们在Type-II 2HDM的基础上，引入两个实标量单态S和P。为了避免在树级出现味改变中性流（FCNC）这一与实验严重冲突的现象，我们对两个希格斯二重态Φ1和Φ2施加了一个Z2对称性（Φ1→Φ1, Φ2→-Φ2）。这个对称性可以“软破缺”，即通过一个维度为2的项m²₁₂(Φ₁†Φ₂ + h.c.)来轻微破坏，这为我们提供了一个自然的退耦极限。同时，为了保证S和P可以作为暗物质候选者，我们对它们施加了另一个独立的Z‘2对称性（S→-S, P→-P），并且要求这个对称性在真空下保持完好，即S和P的真空期望值（vev）为零。此外，我们假设理论整体是CP守恒的，因此势能中的所有参数都是实数。

基于这些对称性要求，我们可以写出最一般的、可重整化的标量势V。为了清晰起见，我们通常将其分为二次项V₂和四次项V₄。二次项包含了所有场算符的平方项： V₂ = m²₁₁(Φ₁†Φ₁) + m²₂₂(Φ₂†Φ₂) - [m²₁₂(Φ₁†Φ₂) + h.c.] + ½ m²_S S² + ½ m²_P P² - m²_SP S P。

这里，m²₁₁, m²₂₂, m²_S, m²_P是质量参数，m²₁₂是软破缺Z2对称性的参数，而m²_SP是连接两个单态S和P的混合质量项。值得注意的是，m²_SP项的存在是允许的，因为它同时包含了S和P，在Z‘2变换下（S→-S, P→-P）是偶性的。

四次项V₄则更为复杂，包含了所有可能的四次相互作用。我们可以进一步将其分类为几个部分： V₄ = V_N + V_CB + V_HC + V_O。

• V_N（“正规”部分）：这部分可以写成一个关于变量r_i = {|Φ₁|², |Φ₂|², S², P²}的二次型，是分析势能有界性（BFB）的关键。它包含了如(Φ₁†Φ₁)², (Φ₂†Φ₂)², (Φ₁†Φ₁)(Φ₂†Φ₂), S⁴, P⁴, S²P²以及希格斯场与单态场的耦合项如(Φ₁†Φ₁)S²等。
• V_CB（电荷破缺部分）：这部分与一个特定的组合z₁₂ = (Φ₁†Φ₁)(Φ₂†Φ₂) - (Φ₁†Φ₁)(Φ₂†Φ₂)相关，仅包含λ₄耦合常数。它在分析电荷破缺真空时至关重要。
• V_HC（硬性CP破坏部分）：包含(Φ₁†Φ₂)²及其厄米共轭项，与λ₅耦合相关。在CP守恒的假设下，λ₅是实数。
• V_O（“奇数次”部分）：这是模型中最具特色也最棘手的部分，包含了像S³P, SP³, (Φ₁†Φ₁)SP, (Φ₂†Φ₂)SP这样的项（对应耦合λ₁₃, λ₁₄, λ₁₅, λ₁₆）。之所以称其为“奇数次”，是因为在将势能表达为前述双线性形式时，这些项无法简单地纳入一个正定的二次型中，这给势能有界性的分析带来了额外的复杂性。

在自发对称性破缺后，希格斯二重态获得非零的真空期望值：⟨Φ₁⟩ = (0, v₁/√2)^T, ⟨Φ₂⟩ = (0, v₂/√2)^T，其中v = √(v₁² + v₂²) ≈ 246 GeV即为电弱对称性破缺能标。而单态场S和P的vev为零，保持Z‘2对称性。由此，我们可以推导出模型的粒子谱。

模型的粒子内容非常丰富：

1. 带电扇区：一对带电的戈德斯通玻色子G⁺被W玻色子“吃掉”，剩下一个物理的带电希格斯玻色子H⁺。
2. CP奇扇区：一个中性的戈德斯通玻色子G⁰被Z玻色子“吃掉”，剩下一个物理的CP奇希格斯玻色子A。
3. CP偶扇区：这是最有趣的部分。两个希格斯二重态中的CP偶中性分量（y₁, y₂）通过一个2x2的混合矩阵（混合角α₁）对角化，得到两个物理的CP偶希格斯玻色子h₁和h₂。我们将其中质量约为125 GeV的粒子认定为与LHC发现的希格斯粒子一致。同时，两个实单态S和P（它们本身是CP偶的）通过另一个2x2的混合矩阵（混合角α₂）对角化，得到另外两个物理的CP偶标量粒子h₃和h₄。关键在于，由于Z‘2对称性保持，h₃和h₄与标准模型粒子的耦合受到抑制，它们构成了一个“暗物质扇区”。最轻的那个（假设是h₃）如果满足其他条件，就可以成为暗物质候选者。

因此，整个模型包含了四个CP偶希格斯（h₁, h₂, h₃, h₄），一个CP奇希格斯A，以及一个带电希格斯H⁺。模型的自由参数众多，包括这些粒子的质量（m_h1, m_h2, m_h3, m_h4, m_A, m_H±）、混合角（α₁, α₂, tanβ = v₂/v₁���、软破缺参数m²₁₂、单态扇区的质量参数（m²_S, m²_P, m²_SP）以及多个四次耦合常数（λ₆, λ₈, λ₉, λ₁₀, λ₁₂, λ₁₃, λ₁₄, λ₁₆）。我们的任务就是在这样一个高维参数空间中，找到那些物理上可接受的“岛屿”。

3. 理论约束：势能有界性、真空稳定性与微扰幺正性

一个物理上合理的量子场论模型，其标量势必须满足一系列基本的理论自洽性条件。这些条件不是可选的，而是模型得以成立的基石。对于2HDM2S这样复杂的势能，严格推导这些条件是确保后续数值扫描可靠性的关键。

3.1 势能有界性（Boundedness From Below, BFB）

势能有界性要求标量势在任何一个场方向趋于无穷大时，其值必须趋向正无穷。否则，真空将是不稳定的，理论没有意义。对于2HDM2S，挑战主要来自V_O部分，因为它包含了像S³P这样的三次方项，无法直接写成场量平方的二次型。

我们的处理策略是寻找充分条件。我们通过不等式放缩，为V_O部分找到一个“下界”V_O^lower，这个下界可以表达为场量平方（r_i = {|Φ₁|², |Φ₂|², S², P²}）的二次型。例如，对于λ₁₃ S³P项，我们利用不等式 -√(r₃)√(r₄) ≥ -(r₃ + r₄)（当r₃, r₄ ≥ 0时），将其放缩为 -|λ₁₃|(r₃² + r₃ r₄) 的形式。对V_O中所有项进行类似处理，并结合V_N, V_CB, V_HC的类似处理，最终我们将整个势能V的下界表达为一个关于向量(r₁, r₂, r₃, r₄)的二次型：(1/2) Σ_{ij} r_i \bar{A}_{ij} r_j。

这里，矩阵 \bar{A} 的元素由原始耦合常数λ_i以及它们的绝对值|λ_i|组合而成（见原文公式4.15）。于是，问题转化为：确保这个4x4的对称矩阵\bar{A}在所有分量r_i ≥ 0的情况下是正定的。判断一个4x4矩阵的“共正性”（copositivity）比3x3矩阵复杂得多，我们采用了文献中一个包含49个步骤的算法来实现这一检查。这是一个关键点：我们推导出的BFB条件是充分的而非必要的，这意味着满足这些条件的参数点肯定是BFB的，但可能有一些BFB的点被我们排除在外。这在数值扫描中是一种保守但安全的策略。

3.2 全局最小值与真空结构

我们想要的物理真空是电弱对称性破缺、但Z‘2对称性保持的真空（即v₁, v₂ ≠ 0, v_S = v_P = 0）。然而，如此复杂的势能可能拥有多个局域极小值。我们必须确保我们研究的这个真空是全局最小值（即能量最低的真空），否则宇宙将倾向于处在另一个更低能量的状态，我们讨论的真空就是亚稳的，可能在宇宙年龄内发生衰变。

我们系统地分析了所有可能的静止点（即一阶导数为零的点），包括：

• N型真空：即我们想要的真空。
• Ns, Np, Nsp型真空：电弱对称性破缺，但同时S或/和P也获得非零vev，这破坏了Z‘2对称性，暗物质候选者不再稳定。
• 电荷破缺（CB）真空：希格斯场的带电分量获得vev，破坏了电磁U(1)对称性，导致光子获得质量，与实验不符。
• CP破缺真空：希格斯场的CP奇部分获得vev，引入额外的CP破坏源。
• 纯单态真空（S, P, SP）：希格斯场vev为零，电弱对称性未破缺，这与现实严重矛盾。

对于其中一些真空（如纯单态真空、部分混合真空），其势能值V_min可以解析求出（见原文第5节公式5.4-5.18）。我们可以直接比较我们想要的N型真空的势能值（公式5.2）与这些真空的势能值。如果任何一个其他真空的势能值更低，则该参数点被丢弃。

对于另一些无法解析求解的真空（如某些特定的CP破缺或混合真空），我们采用了数值方法（例如，通过变量代换v₁ = v cosβ, v₂ = v sinβ，将方程组化简后数值求解β）。在数值扫描中，对于每一个通过其他约束的参数点，我们都会调用这些数值程序来检查是否存在能量更低的竞争真空。

3.3 微扰幺正性（Perturbative Unitarity）

在量子场论中，散射振幅必须满足幺正性（概率守恒）。在微扰论中，这转化为对散射矩阵的本征值施加约束。对于标量场的高能二体散射过程，幺正性要求其分波振幅的绝对值小于某个常数（通常取1/2）。这直接转化为对模型中四标量耦合常数组合的约束。

我们遵循标准方法，构造所有可能电荷（Q）和超荷（Y）组合下的二标量初态，计算其s波散射振幅矩阵，并要求该矩阵的所有本征值的绝对值小于4π（或8π，取决于约定，我们采用更严格的4π以保证微扰论的有效性）。对于2HDM2S，由于引入了两个新的单态场，散射态的数量比普通的2HDM或N2HDM更多。我们需要计算7个散射矩阵（M^{++}_2, M^{+}_2, M^{+}_1, M^{+}_0, M^{0}_2, M^{0}_1, M^{0}_0），其中M^{+}_1和M^{0}_1是新增的，涉及单态场与希格斯场的散射。

大部分矩阵的本征值可以解析求出（如公式6.3, 6.5, 6.7, 6.9）。最复杂的是中性、零超荷的散射矩阵M^{0}_0，它是一个11x11的矩阵。其本征值中，有6个可以解析得到，另外5个是一个五次多项式的根。在实际计算中，我们更高效地采用了计算矩阵主子式的方法来施加幺正性约束。微扰幺正性约束非常强大，它通常会将许多大耦合常数的参数区域排除在外，确保我们的理论始终处于微扰论有效的范围内，从而使我们的预测是可靠的。

4. 实验与观测约束：对撞机与暗物质

一个模型即使理论自洽，也必须经受住现有实验数据的严格检验。对于2HDM2S，约束来自两大方面：高能对撞机实验和天体物理观测。

4.1 对撞机物理约束

1. 希格斯信号强度：LHC上发现的125 GeV希格斯粒子（在我们的模型中对应h₁或h₂，通常将h₁认定为125 GeV粒子）的产生和衰变信号，必须与ATLAS和CMS的测量结果在2σ置信水平内相符。我们通过计算耦合修饰子（coupling modifiers）来实现这一点。例如，h₁到底夸克、τ轻子、W/Z玻色子的耦合相对于标准模型希格斯耦合的比值，由混合角α₁和β决定：κ_b ≈ sinα₁ / cosβ, κ_τ ≈ sinα₁ / cosβ, κ_V ≈ sin(β - α₁)。实验数据强烈倾向于所谓的“对齐极限”（alignment limit），即sin(β - α₁) ≈ 1，此时h₁的性质非常接近标准模型希格斯。
2. 希格斯不可见衰变：如果暗物质候选者h₃的质量小于m_h₁/2 ≈ 62.5 GeV，那么希格斯粒子h₁可以衰变到一对h₃h₃。这种“不可见衰变”分支比受到LHC的严格限制。目前ATLAS和CMS给出的95%置信水平上限约为BR(h₁ → inv) < 10.7%。
3. 新粒子直接搜索：LHC持续搜索超出标准模型的新粒子。对于2HDM2S，相关的搜索包括：
   • 带电希格斯H⁺：主要衰变道为H⁺ → τ⁺ ν_τ 或 H⁺ → t \bar{b}。根据最新的全局拟合，在Type-II 2HDM中，m_H⁺ ≳ 580 GeV (95% CL)。
   • 重中性希格斯：对A, H (即h₂) 的搜索，通常通过它们衰变到τ⁺τ⁻, b\bar{b}, t\bar{t}或一对规范玻色子等末态来设定质量下限。
   • 双希格斯产生：寻找非标准模型来源的双希格斯玻色子产生过程。 我们使用HiggsTools软件包来系统检查这些约束，它集成了ATLAS��CMS的最新分析结果。
4. 电弱精密测量：高精度的Z玻色子物理和低能实验数据通过所谓的S, T, U参数（或Oblique Parameters）来约束新物理。这些参数描述了新物理粒子对规范玻色子自能图的贡献。我们基于文献中的公式，计算了2HDM2S模型对S, T, U的贡献，并要求其落在最新的全局电弱拟合结果（如PDG或Gfitter小组给出的拟合值）的2σ椭圆内。这主要约束了新标量粒子的质量劈裂，例如，大的m_H⁺ - m_A劈裂会导致对T参数的过大贡献。

4.2 暗物质物理约束

我们将模型实现到micrOMEGAs软件中，用于计算暗物质相关的所有可观测量。

1. 遗迹密度：暗物质候选者（假设为最轻的Z‘2-奇粒子，通常是h₃）在早期宇宙热浴中通过与其他粒子的相互作用保持热平衡。随着宇宙膨胀冷却，当其湮灭速率低于哈勃膨胀速率时，它会“冻结”出来，留下今天的遗迹丰度。我们要求计算出的Ωh²（暗物质密度与临界密度的比值乘以哈勃常数的平方）与普朗克卫星的精确测量值Ωh² = 0.1200 ± 0.0012相符（在3σ范围内）。在初步扫描中，我们有时会采用一个更宽松的范围（如0.09-0.15）以获取更多样本点。
2. 直接探测：暗物质粒子可能与原子核发生弹性散射，被地下探测器（如LZ、XENONnT）捕捉到。散射截面主要由暗物质粒子与核子（质子和中子）的耦合决定。在2HDM2S中，h₃与核子的耦合主要通过t-channel交换h₁和h₂实现。我们计算了暗物质与氙原子核的自旋无关散射截面σ^{SI}_Xe，并与LZ合作组2024年发布的最新排除线进行比较。这是目前最严格的约束之一，许多WIMP模型（包括N2HDM）的参数空间因此被大幅压缩。
3. 间接探测：暗物质粒子在星系中可能相互湮灭，产生高能宇宙射线（如伽马射线、正电子、反质子）。费米卫星（Fermi-LAT）、AMS-02实验和H.E.S.S.等对这类信号进行了搜寻。我们计算了暗物质在今天的银河系晕中的热平均湮灭截面⟨σv⟩，并针对主导的湮灭末态（在我们的模型中主要是b\bar{b}夸克对，或W⁺W⁻/ZZ规范玻色子对），与相应实验的排除上限进行比较。不过，对于满足正确遗迹密度的点，其今天的湮灭截面通常远低于当前间接探测的灵敏度，因此这条约束相对较弱。

一个重要的技术细节：我们的Z‘2对称性作用于两个单态S和P，这意味着h₃和h₄都是Z‘2-奇的。如果m_h₄ > m_h₃，那么h₄可以衰变为h₃加上标准模型粒子（例如通过λ₁₅或λ₁₆耦合项：h₄ → h₃ h₁/h₂，随后h₁/h₂衰变到标准模型粒子）。因此，只有最轻的Z‘2-奇粒子h₃是稳定的暗物质候选者。在micrOMEGAs中，我们使用darkOmegaN例程（而非默认的darkOmega2）来计算包含这种衰变链的遗迹密度。

5. 参数空间扫描策略：从随机漫步到智能演化

2HDM2S模型的参数空间维度高达22维（见公式2.16），且受到复杂非线性理论约束和实验约束的限制。传统的均匀随机扫描在这种高维空间中效率极低，如同大海捞针。我们采用了三种互补的策略来系统性地探索这个空间。

5.1 策略一：无先验的随机扫描

这是最直接的方法。我们固定已知参数：电弱能标v = 246 GeV，以及一个125 GeV的希格斯质量m_h1。其他参数在以下范围内随机均匀抽取：

• 混合角：α₁, α₂ ∈ [-π/2, π/2]；tanβ ∈ [0.3, 10]。
• 标量粒子质量：m_h2, m_h3, m_h4 ∈ [125, 1000] GeV；m_A ∈ [100, 1000] GeV；m_H± ∈ [580, 1000] GeV（基于带电希格斯质量下限）。
• 质量平方参数：m²₁₂, m²_S, m²_P, m²_SP ∈ [±10⁻¹, ±10⁷] GeV²。
• 四次耦合常数：λ₆, λ₈, λ₉, λ₁₀, λ₁₂, λ₁₃, λ₁₄, λ₁₆ ∈ [±10⁻³, ±10¹]。

我们运行了约7000个CPU小时，产生了约12万个随机点。然后，我们逐一施加第3、4节所述的所有理论约束（BFB、全局最小、幺正性、微扰性）和对撞机约束（耦合修饰子、信号强度、HiggsTools）。结果如图2（原文）所示，在α₁-β平面上，只有极少数点（红色和绿色）能通过所有基础约束，而其中能同时满足暗物质遗迹密度要求的点（蓝色）更是凤毛麟角。这直观地展示了传统随机扫描在寻找可行点上的低效性。

5.2 策略二：近对齐极限扫描

LHC数据强烈暗示新物理与电弱扇区的耦合接近“对齐极限”，即sin(β - α₁) ≈ 1。在此极限下，125 GeV的希格斯粒子h₁的性质与标准模型预言几乎一致。因此，我们可以将扫描聚焦在α₁ ≈ -β（对应cos(β - α₁) ≈ 0）的狭窄区域附近，例如让α₁在[-1.1β, -0.9β]范围内随机取值。其他参数的扫描范围与策略一相同。

这种方法效率有所提升（在相似计算时间内获得了约37万个点），因为它直接瞄准了实验上最受青睐的区域。图3（原文）展示了在暗物质质量-遗迹密度平面上，近对齐扫描的结果。可以看到，通过所有基础和对撞机约束的点（绿色）明显增多，但满足遗迹密度要求的点（蓝色）仍然集中在少数几个质量区域，且大部分点（红色）因违反各种约束而被淘汰。更重要的是，当我们把通过遗迹密度筛选的蓝点与LZ直接探测排除线（图4，原文）对比时，几乎所有点都被排除了。这清晰地表明，在严格的直接探测约束下，通过传统随机方法找到“幸存”参数点异常困难。

5.3 策略三：基于进化策略的机器学习优化

面对高维、强约束的优化问题，我们求助于人工智能中的黑盒优化算法。我们采用了协方差矩阵自适应进化策略（CMA-ES）。其核心思想是模拟生物进化：

1. 初始化：在参数空间中随机选择一个点作为“父代”种群的均值，并初始化一个协方差矩阵（描述参数间的关联和搜索步长）。
2. 生成子代：从当前的多变量正态分布中采样，生成一组“子代”候选点。
3. 评估与选择：用一个损失函数L(θ)来评估每个候选点的“好坏”。我们将所有约束（理论+实验）编码为惩罚项C(O)。对于一个约束O，如果其值在允许区间[O_LB, O_UB]内，则C(O)=0；否则，C(O)等于超出边界的距离。损失函数L(θ)就是所有约束惩罚项之和。L(θ)=0的点就是满足所有约束的“好点”。
4. 更新：从子代中选出损失函数最小的几个优秀个体，用它们的分布信息来更新下一代的均值向量和协方差矩阵，使搜索方向朝着损失函数减少的区域移动。
5. 迭代：重复步骤2-4，种群不断“进化”，最终收敛到损失函数最小（理想情况下为0）的区域。

然而，标准的CMA-ES有一个缺点：容易陷入局部最优。一旦找到一个满足约束的区域，它就会集中在该区域采样，而可能错过其他同样好的、但相隔较远的区域。为了解决这个问题，我们引入了新颖性奖励机制。我们记录下所有已找到的“好点”，当算法评估一个新点时，不仅计算其损失函数，还计算它与历史好点在参数空间和可观测量空间（如暗物质质量、散射截面）的“密度距离”。如果一个点位于已有好点密集的区域，它会受到一个额外的惩罚；反之，如果一个点位于稀疏区域，即使它的损失函数略高，也可能因为“新颖”而被保留。这迫使算法去探索更广阔的空间。

我们首先进行了一轮不带新颖性奖励的CMA-ES优化，目标是最小化损失函数（包含遗迹密度约束，但暂未加入直接/间接探测约束）。在约1500个CPU小时内，算法迅速找到了大量满足遗迹密度要求的点（约66万个），如图5（原文）所示。这些点（绿点和蓝点）主要集中在暗物质质量65-300 GeV的区间。然后，我们在损失函数中加入了LZ直接探测和费米卫星等间接探测的约束，并开启了新颖性奖励。这是关键一步。算法开始有目的地避开那些虽然满足遗迹密度但会被直接探测排除的区域，转而探索那些耦合结构特殊、导致核子散射截面被压低的参数空间。最终，我们成功地找到了同时满足所有理论约束、对撞机约束、暗物质遗迹密度和直接/间接探测约束的“幸存”参数点。

6. 结果分析与物理启示

通过结合传统扫描和机器学习优化，我们对2HDM2S模型的可行参数空间有了一个全景式的认识。图5（原文）集中展示了优化后的结果。图中的蓝点代表了模型的“幸存者”，它们通过了所有严苛的筛选。我们可以从中提炼出一些重要的物理模式和启示：

1. 暗物质质量范围：可行的暗物质候选者质量m_h3主要分布在两个区域：一个是低质量区（约62.5 - 80 GeV），另一个是中高质量区（约200 - 600 GeV）。这与湮灭机制密切相关。
2. 低质量区的“希格斯共振”机制：当m_h3 ≈ m_h1/2 ≈ 62.5 GeV时，暗物质通过s-channel交换125 GeV的希格斯玻色子h₁进行湮灭，此时存在共振增强效应。即使暗物质与核子的耦合本身很小，共振处的湮灭截面也足以产生正确的遗迹密度。更重要的是，在共振点附近，暗物质与核子的散射截面（主要通过t-channel的h₁交换）可以被共振效应所抑制，因为此时传递的动量与h₁的质量壳层条件相差很大。这为躲避直接探测提供了天然的机制。我们的扫描确实在m_h3 ~ 65 GeV附近找到了幸存点。
3. 中高质量区的“共湮灭”与“漏斗”机制：在更高的质量区域，简单的希格斯共振不再有效。暗物质湮灭主要通过以下几种渠道：
   • 湮灭到规范玻色子对（WW/ZZ）：当m_h3 > m_W, m_Z时，这个道打开，截面较大。
   • 湮灭到顶夸克对（t\bar{t}）：当m_h3 > m_t时，这个道打开。
   • 共湮灭（Co-annihilation）：如果暗物质粒子h₃的质量与另一个Z‘2-奇粒子（如h₄或A）的质量非常接近（Δm/m ≲ 10%），那么在宇宙早期，h₃和h₄/A可以相互转换并共同参与湮灭过程。这有效地增加了总的湮灭截面，从而在耦合较小的情况下也能获得正确的遗迹密度。共湮灭机制通常有助于降低直接探测截面。
   • A/H漏斗（Funnel）：如果暗物质质量接近某个中性希格斯玻色子（A或h₂）质量的一半，即m_h3 ≈ m_A/2 或 m_h2/2，则可以通过s-channel交换A或h₂发生共振湮灭。这与低质量的h₁共振类似。 机器学习优化找到的中高质量区幸存点，往往对应着这些机制的某种组合，并且其与核子的耦合λ_{h3-Nucleon}被巧妙地压低，以满足LZ的极限。
4. 标量混合角的影响：混合角α₂（控制暗物质扇区h₃和h₄的混合）扮演着关键角色。它直接影响暗物质粒子h₃与标准模型希格斯h₁/h₂的耦合强度。较小的混合意味着h₃更“惰性”，与标准模型的耦合更弱，这有利于降低直接探测截面，但同时也降低了湮灭截面，可能使遗迹密度过高。因此，需要一个精细的平衡。我们的扫描发现，幸存点对应的sinα₂通常不大，但也不为零，表明一定的混合是必要的。
5. 与N2HDM的对比与进步：如前所述，早期的N2HDM（单标量单态扩展）研究显示，在满足遗迹密度和直接探测约束方面已经捉襟见肘。2HDM2S模型通过引入第二个单态，增加了两个关键的自由度：额外的暗物质粒子质量m_h4和混合角α₂。这带来了新的可能性：
   • 新的湮灭和共湮灭通道：h₃ h₃ → h₄ h₄, h₃ h₄ → SM particles等。
   • 更灵活的耦合调节：λ₁₅和λ₁₆耦合项提供了h₃/h₄与希格斯二重态之间新的相互作用顶点，可以更灵活地调节湮灭截面和散射截面之间的关系。
   • 规避直接探测：通过h₃-h₄混合以及新的相互作用，有可能实现破坏性干涉，使得h₃与核子的散射振幅中，不同费曼图（如通过h₁和h₂交换）的贡献相互抵消，从而显著压低散射截面。 正是这些额外的“调节旋钮”，使得2HDM2S在LZ等新一代直接探测实验的严峻挑战下，依然能找到存活的参数空间，而N2HDM的类似区域可能已被完全排除。

7. 总结与展望

通过对2HDM2S模型的全面研究，我们展示了在双希格斯二重态框架下引入两个实标量单态，可以为弱相互作用大质量粒子（WIMP）暗物质提供一个可行的理论平台。这项工作不仅仅是提出一个模型，更完成了一项从理论构建、约束推导到系统性数值扫描的完整“流水线”作业。

理论层面，我们严格推导了该模型势能有界性的充分条件、全局最小值的判定方法以及微扰幺正性约束，这些是任何严肃的BSM模型分析必须跨越的门槛。方法学层面，我们对比了传统随机扫描和基于CMA-ES进化策略的机器学习优化方法。结果清晰地表明，在面对高维、多约束的参数空间时，智能优化算法在寻找可行点方面的效率有数量级的提升，特别是当结合“新颖性奖励”机制来促进探索时。

物理层面，我们的扫描证实，在当今最严格的对撞机约束（特别是125 GeV希格斯性质测量）和暗物质直接探测（LZ实验）的双重压力下，2HDM2S模型仍然拥有存活的参数空间。这些幸存点主要依赖于希格斯共振机制（在m_h3 ≈ 62.5 GeV附近）和共湮灭/漏斗机制（在200-600 GeV质量区间），并通过暗物质扇区内的混合以及新的耦合项来巧妙地压低核子散射截面。

这项研究为未来的探索指明了方向：

1. 对撞机信号：下一步可以深入研究这些幸存参数点在对撞机上的具体信号。例如，寻找 exotic Higgs decays（如h₂ → h₁ h₁, h₂ → h₃ h₃），或带电希格斯H⁺的特殊衰变模式。在LHC的高亮度运行阶段，这些可能是发现新物理的关键通道。
2. 相变与引力波：丰富的标量扇区可能带来强烈的一阶电弱相变，从而产生可探测的随机引力波背景。未来空间引力波探测器（如LISA）的观测数据可以与对撞机搜索形成互补。
3. 中微子质量：本文模型专注于暗物质问题。一个自然的扩展是将此标量扇区与中微子质量生成机制（如跷跷板机制）联系起来，构建一个能同时解决暗物质和中微子质量问题的统一框架。
4. 算法改进：可以尝试将其他更先进的优化算法（如贝叶斯优化、强化学习）应用于此类高能物理模型的参数扫描中，并开发更高效的并行计算和约束预筛选策略。

暗物质的性质是当代物理学最深刻的谜题之一。像2HDM2S这样的具体模型，为我们提供了将宇宙学观测与微观粒子物理联系起来的可检验框架。尽管参数空间日益狭窄，但通过更精密的实验、更深入的理论计算和更智能的数值工具，我们正在一步步缩小暗物质可能藏身的范围。这项工作正是沿着这条道路迈出的坚实一步。