2026-10-01
把 M1/M2 受体图谱装进 68 脑区皮层模型：结构化异质的同步与信息流显著超同质网络，清醒样态中涌现局部睡眠慢波。
主流全脑计算模型有个长期默认：每个皮层区域跑同一套局部动力学，区域之间的差别全部交给结构连接（SC，用弥散磁共振纤维追踪测得的白质连接）去解释。真实皮层明显不是这样，受体密度、基因表达、髓鞘化程度都沿皮层表面呈梯度分布。此前把 T1w:T2w、兴奋抑制比这类区域差异加进模型，确实能改善对实证功能连接的拟合，但分子尺度的异质性对大尺度动力学有多大影响、是否随脑状态变化，一直没有系统答案。
乙酰胆碱（ACh）是切入口。它同时调控皮层兴奋性与神经元适应，清醒与 REM 睡眠时水平高、深睡眠时低；基底前脑投向皮层的胆碱能纤维在空间上并不均匀，小鼠实验里胆碱能投射弱的区域，恰好是 REM 睡眠中冒出局部慢波的区域。巴塞罗那 IDIBAPS 的 Sanchez-Vives 团队联合 Paris-Saclay 的 Destexhe，把这个观察变成了可计算的问题。
框架三层：局部动力学、结构连接、区域异质性。
稳健性两道：用 PET 独立测得的 M1 受体图复跑全部实验；用三个零模型做特异性检验，空间打乱受体图、均值方差匹配的高斯随机图、保留空间自相关的旋转替代图，各 100 次实现。
先复现状态连续谱：平均适应 12 pA 时各区低幅高频活动（主频约 10 Hz）、区间相关低；24 pA 时局部同步簇成形、主频降到约 6 Hz;36 pA 时全网进入约 1 Hz 的 Up/Down 慢波并以行波传播，功能连接与结构连接的耦合从松散变紧，与实证数据里脑状态依赖的耦合变化一致。
异质性的收益集中在信息流。刺激上顶叶 30 ms，测其后 400 ms 窗口内从源区到其余 67 区的传递熵（TE）：
| 网络配置 | 清醒样态 TE(10 pA) | 睡眠样态 TE(50 pA) |
| 转录组异质 | 1.38±0.09 | 0.30±0.07 |
| PET 异质 | 1.28±0.08 | 0.28±0.09 |
| 同质 | 1.23±0.08 | 0.25±0.06 |
同步上，平均适应超过 30 pA 后异质网络的平均 FC 反超同质网络；相位滞后指数（PLI，对零滞后相关不敏感，适合抓行波）在全适应范围都更高，且随 hf 增大单调上升。特异性是混合的：十二组零模型对照六组过 P<0.05，其中四组在睡眠样态，如空间打乱图下睡眠态 z=2.36-2.89、旋转替代图下中间态（hf=1）z=3.32；高斯图在中间态两组都不显著（z 不超过 1.52）。清醒态 TE 与高斯图的差异同样不显著。作者的读法：一般性异质就足以维持宽泛的状态转变，生物结构化的异质在同步态上做选择性增强。
中间态（24 pA）出现全文最抓人的现象：少数区域已进入 Up/Down 慢波，其余区域保持清醒样活动，两种动态在同一个网络里共存。哪个区域先睡由两个弱相关的变量共同决定（r=-0.23,P=0.056）：该区适应水平与 delta 功率（0-4 Hz）正相关，入度连接强度与 delta 功率负相关，即接收兴奋输入越少的区域越容易先睡。零模型下适应与 delta 的相关明显衰减，入度与 delta 的相关向同质基线收拢：受体图的空间组织既放大适应的作用，也削弱结构连接对慢波落点的单独决定权。这与小鼠实验中胆碱能投射弱的区域先出现局部慢波的观察对得上。
对 AI 从业者，这篇是一条旁证。脉冲神经网络与振荡器网络文献里早有单元异质性提升编码与同步的结果（论文引用的 Random heterogeneity outperforms design in network synchronization 是代表），这篇把该原则推到有生物学约束的人脑尺度：大脑没有把异质性平均掉，而是用分子梯度把它组织成全脑通信的增强器。做神经形态硬件或异构单元网络的人，多了一个生物侧的设计论据，不是可直接搬用的配方。
对神经科学侧，这套框架把受体图谱（转录组加 PET）与能产生状态转变的全脑模型装在了一起，并给出局部睡眠的机制解释：适应高、入度低的区域先睡。注意力涣散、睡眠剥夺下的绩效下降、病灶周围慢波都能放进同一框架重算，作者明确指向个体化病变模型。
作者自述三条：胆碱能只建模了适应这一条通路，毒蕈碱受体的大量突触作用未纳入；M1 与 M2 被合并成单一参数，PET 侧只用了 M1；受体图是静态的群体均值，没有神经调制的时间动态。
读下来还要补两点。效应量绝对值小，PLI 差值在 0.001-0.015 量级，论文报的是相对零分布的标准化 z 值；局部慢波的验证停留在模型内部相关，没有对照实证脑电或颅内电极测到的慢波空间分布，这套解释目前是自洽，还不是被外部数据钉死。