能量供应与生物多样性之间的关系是生态学中的一个长期存在的问题。尽管通常假设多样性随能量可得性单调增加，但在整个生态系统中广泛记录了单峰的物种-能量关系（species–energy relationship），其从第一性原理（first principles）出发的起源仍不清楚。在此，研究人员开发了一个几何框架，将生态可行性（feasibility）明确地重述为能量项。通过将总能量供应视为系统层面的约束并应用于基于能量的网络模型（energy-based network model），研究人员在能量捕获率空间中定义了嵌套的可行性域（feasibility domains），并将可行性概率量化为它们的体积比。研究人员表明，初始化一个可行网络的概率随能量供应单调增加并达到饱和，而维持稳态生物量（steady-state biomass）的概率遵循单峰关系——揭示了一个有限的能量窗口，在此窗口内网络成熟（maturation）最有可能。通过可行性划分（feasibility partitions）将此分析扩展到所有候选子群落，研究人员发现不同群落大小在不同能量水平下最具可行性，且平均多样性本身在中等供应量时达到峰值。总之，这些结果表明，能量约束决定生态网络的多样性不仅仅是通过能量稀缺，而是通过外部能量供应与内部能量交换之间的几何相互作用。
论文解读：《Energetic constraints shape the diversity of feasible ecological networks》
研究背景与意义：
能量与生物多样性之间的关系是生态学的核心问题之一。传统观念常假设生物多样性随外部能量供应（通常等同于初级生产力）单调增加。然而，大量实证观察表明，物种-能量关系（species–energy relationship）往往呈单峰（unimodal）曲线，即在草地、湿地、池塘、藻华和雨林等多种生态系统中，最高生物多样性出现在中等能量水平，而非能量最丰富的区域。这种单峰关系的起源长期以来缺乏从第一性原理（first principles）出发的理论解释。现有理论如共存理论（coexistence theory）侧重于资源竞争和生态位分化，但未将总能量可用性作为明确的系统级约束；消费者-资源理论虽然涉及能量通量和代谢约束，但多局限于资源空间的特定切片；结构方法虽定义了交互作用的几何可行性域（feasibility domain），却难以直接引入系统级能量约束。为此，研究人员开展了一项理论研究，旨在从能量加工与几何可行性的角度，探究系统级能量约束如何塑造生态网络的可行性与多样性，论文发表于《PLOS Computational Biology》。
主要关键技术方法：
研究人员构建了一个基于能量的生态网络模型，将S个种群的生态系统视为开放热力学系统，以单位质量能量捕获率si、维持需求di及能量交换矩阵（交互作用矩阵）M描述生物量动态，并导出广义Lotka-Volterra形式。研究引入总能量供应率Q作为系统级约束，定义了捕获域（capture domain, DC）、初始化域（initialization domain, DI）和成熟域（maturation domain, DM）三个嵌套几何域。可行性概率通过域体积比计算：初始化概率PI= vol(DI)/vol(DC)，成熟概率PM= vol(DM)/vol(DC)。DC体积解析求得；DI（凸多面体）通过V-表示与H-表示分解算法精确计算体积；DM（含二次约束的凸体）通过Cholesky变换标准化为InterPolyQuads类，并利用Hit-and-Run MCMC采样与多相蒙特卡洛方法估计体积（相对误差<3%，S≤12）。为涵盖部分共存，研究定义了可行性划分（feasibility partitions），枚举所有2S个候选群落。参数矩阵通过四步随机矩阵程序生成：采样、效率约束校正、耗散性（Volterra dissipative）对角偏移、整体缩放；最小生物量N0与能量需求d从正态分布抽取并取正。以S=8的100个独立随机网络为主要对象，并通过降采样、参数缩放与连接度变化进行敏感性分析。概率曲线通过功能数据分析（FDA，logit变换）提取总体趋势与变异。
研究结果：
Energetic regimes in a two-population network（两种群网络中的能量 regime）：
以一个随机取参的两种群（S=2）网络为例，随着总能量供应Q增加，网络可行性呈现三个相位。低能区（Q < Qc，临界阈值）下初始化与成熟概率均为零；中能区（Qc< Q < Qopt）两者均上升；高能区（Q > Qopt）初始化概率继续增加至饱和，而成熟概率在Qopt达峰值后下降，呈现单峰响应。几何上，初始化受稳态生物量可行锥与最小能量捕获约束交截控制，Q增大使交截体积比趋近几何极限从而饱和；成熟还受稳态能量捕获约束（sTN* ≤ Q）影响，该约束在高Q时逐步削减可行域，导致最优能量窗口的出现。
Feasibility patterns and drivers for multi-species full communities（多物种完整群落的 feasibility模式与驱动因素）：
在S = 2~8的随机网络群体中，上述模式具有鲁棒性：初始化概率PI随Q单调上升至饱和，存在临界Qc；成熟概率PM随Q呈单峰，存在最优Qopt定义最有利能量窗口。两者均随网络规模S增加而下降，体现“复杂性的能量成本”。即使引入密度依赖的饱和资源吸收动力学，PM在极高Q时趋于非零极限，但中低Q区间的单峰特征依然保持。最小生物量N0升高会提升Qc与Qopt及最大PM；能量交换强度增大则主要提升Qopt与对应PM，而不改变Qc。
Feasibility and average diversity across partial and full communities（部分与完整群落的可行性及平均多样性）：
将分析扩展至S=6网络的所有候选群落（按大小D=1~6分组），各类大小的成熟概率均保持随Q的单峰关系，但最可行能量窗口随群落大小D增大向更高Q偏移：较小群落在较低供应下最可行，较大群落需要更高Q。平均多样性定义为期望群落大小（按候选群落数与对应PM加权），随Q也呈单峰，且峰值能量低于完整群落（D=6）的Qopt。这是因为中等大小群落在组合数量上远多于完整群落，其在中等能量窗口的集体贡献（众多中等可行的中型群落）主导了网络级平均多样性峰值。
讨论部分总结：
将生命斗争视为自由能斗争为研究复杂生命系统提供了统一框架。研究通过引入基于能量的形式体系，将生态动态映射到几何表示，借助高级随机方法估算高维可行性域，揭示了外部能量输入与内部交换共同决定网络可行性的机制。对完整群落的分析显示：初始化存在最小能量门槛并最终饱和；成熟概率呈单峰，揭示最利于持久存在的最优能量供应区，该性质在饱和吸收形式下仍保持。扩展至所有候选群落后，单峰特征依旧，但最可行能量随群落大小系统性上移，故能量增加并不同时提升所有候选群落的可行性，而是沿对数尺度在不同能量范围依次利于不同大小群落。网络级平均多样性也呈单峰，其峰值能量低于完整群落最优能量，因为中等大小群落数量多且在中等能量下集体贡献最大。除了能量供应，网络与群落大小是可行性的关键决定因子，也是能量-多样性关系的中介；多样性限制不仅来自资源稀缺，也来自能量在可行生态网络中可分布的结构约束。研究所给s-空间划分提供了从捕获率到候选群落的静态几何表征，未来可加入动态以描述群落组装/ disassemble 的边界穿越；同时将Q与捕获率s的机理联系（如通过消费者-资源网络）是下一步关键。总体而言，复杂生态群落的持续可能不依赖于最大化能量供应，而依赖于维持在有限能量 regime 内；该框架为理解生命涌现与持续的物理界限提供了基础，并可通向更机理化的生态情景，但应作为基于简化能量表示的基线描述。
研究结论（翻译结论要点）：
能量约束通过外部供应与内部能量交换的几何相互作用（而非仅能量稀缺）决定生态网络多样性；初始化可行性随能量单调饱和，稳态成熟可行性随能量呈单峰，存在有限最优能量窗口；不同大小候选群落的最可行能量水平不同，且网络平均多样性在中等能量供应处达峰，为物种-能量关系中的单峰模式提供了基于第一性原理的几何解释。
