过了国庆的星城，早晚的空气中已经带上了几分凉意。
虽然说白天依旧还是那么热，但倒也没有夏天烧烤般的感觉了。
对于星城这座城市来说，一年中仅有的两三周秋季差不多就是这个时候了。
等...
谭蔚泓推了推眼镜，镜片后的目光透出几分期待与笃定：“那就不耽误你时间了，先回去看看材料，我让老周把实验日程表发你邮箱，顺便附上最近三批样品的原始测试数据包——都是未经任何滤波、归一化或人工修正的‘毛数据’，连噪声谱都原样保留着。”
韩川点头应下，指尖无意识摩挲着手机边缘。他忽然想起什么，抬眼问道：“谭院长，你们测迁移率用的是什么构型？顶栅底接触？还是场效应晶体管阵列？”
“顶栅底接触，SiO?作栅介质，厚度300纳米，源漏电极是热蒸发的金，沟道长度50微米。”谭蔚泓语速不疾不徐，却字字清晰，“不过你别被参数吓住——我们更关注的是动态响应。比如同一样品在连续72小时测试中，迁移率曲线出现三次剧烈跳变：第一次在18小时，第二次在41小时，第三次在67小时。每次跳变后，电流-电压曲线的非线性特征都会重构，但又不完全重复。”
韩川微微颔首，没说话，只是把这句话默默记进心里。他注意到谭蔚泓说“不完全重复”时，右手食指在桌面轻轻点了三下，节奏略有差异——像是某种下意识的校准动作。
走出实验室前，袁荷蓓递来一张薄薄的蓝色门禁卡，卡面印着银色分子式C??H??N?O?S?。“这是生物纳米楼B区三层的临时权限，有效期到下周四。你明天来直接刷这张卡就能进洁净间旁边的分析室，设备已经预热好了。”
韩川接过来，指尖触到卡片微凉的塑料表面。他低头看了眼，背面用激光蚀刻着一行极小的字：*No. 2013-BNME-079 —— “青鸾计划”第七批次*
他不动声色地将卡收进衬衫口袋。
回到宿舍已是下午三点。窗外蝉鸣如沸，空气里浮着一层晒得发烫的尘埃光晕。韩川没开空调，只拉开窗帘，任斜阳泼满整张书桌。他打开笔记本电脑，点开邮箱——两封新邮件静静躺在收件箱里：一封来自袁荷蓓，主题为《青鸾计划·基础资料包_v3》；另一封来自谭蔚泓本人，标题只有四个字：《毛数据_0723》。
附件解压后，是三个文件夹：*Raw_IV_curves*、*Noise_Spectra*、*Capacitance_Sweep*。每个文件夹里都塞着上百个文本文件，命名规则统一为“SampleID_Date_Time_Replica.txt”，例如“Q7-LA-20130722_143211_R2.txt”。
韩川点开第一个IV曲线文件，光标停在第17行——那里有一段突兀的电流回滞：在Vg=-20V处，Ids从1.24×10??A骤降至8.31×10??A，持续3.2秒后又弹回原值。而同一时刻的电容数据却平滑如常，毫无扰动。
他皱了皱眉，调出噪声谱文件对比。在10Hz频段，该样本的电压噪声功率谱密度比相邻批次高出1.8个数量级，且呈现典型的1/f2衰减特征——这不该出现在有机半导体中，更像金属接触界面的热涨落异常。
韩川的手指悬在键盘上方，迟迟没有敲下。他忽然想起CERN强子对撞机数据组曾遇到过类似现象：LHCb探测器在特定束流强度下，部分硅像素传感器会周期性出现“暗电流脉冲”，持续时间精确到毫秒级，振幅服从伽马分布。当时德利涅教授带着他们用了整整两周，才确认那是高能质子轰击传感器边缘氧化层，诱发的局部陷阱态瞬态充放电效应。
而眼前这些有机分子薄膜……会不会也存在某种尚未被识别的亚稳态陷阱集群？
他拉开抽屉，取出一支黑色中性笔，在草稿纸上画下第一个模型草图：横轴是时间t，纵轴是迁移率μ。他先描出三条主趋势线——分别对应三批样品的平均迁移率衰减曲线，然后在每条线上刻意标出三处尖锐的“刺状突起”，位置与谭蔚泓提到的72小时跳变点严格对应。
接着，他在图下方写下一行小字：
*假设陷阱态密度ρ(t)随时间演化满足：?ρ/?t = α·μ(t) - β·ρ(t) + γ·δ(t-t?)*
其中δ函数代表外部扰动触发的瞬时陷阱生成，α、β、γ为待拟合参数。
笔尖顿了顿，他又在旁边补了一行：
*但若陷阱本身具有记忆性——即ρ不仅依赖当前μ，还依赖历史μ的积分路径呢？*
这个念头像一道无声的闪电劈开思维。韩川猛地坐直身体，抓起手机拨通袁荷蓓电话。
“老周，你们做迁移率拟合时，用的是标准的SCLC（空间电荷限制电流）模型，还是考虑了载流子捕获动力学的改进版？”
电话那头传来轻微的键盘敲击声，几秒后袁荷蓓答：“标准SCLC。我们试过加入单陷阱项，拟合残差反而更大。”
“有没有保存过原始转移特性曲线的全量数据？不是只导出迁移率数值，而是完整的Ids-Vg扫描点，至少每0.1V一个采样点？”
“有！全都在Raw_IV_curves里，每个文件第一列是Vg，第二列是Ids，第三列是时间戳。”
“好。”韩川语速加快，“再帮我查三件事：第一，所有跳变发生前10分钟，环境温湿度是否记录异常？第二，跳变时刻的源漏偏压是否有微小波动？第三……你们洁净间的氮气纯度监控，最近一周有没有报警记录？”
袁荷蓓明显愣了一下：“氮气纯度？我们用的是99.999%的高纯氮，监控系统没报过警……等等，我查一下后台日志。”
电话陷入短暂沉默。窗外一只知了突然嘶鸣，声音尖锐得近乎破碎。
三十七秒后，袁荷蓓的声音重新响起，带着一丝不易察觉的紧绷：“韩神……七月二十一号凌晨两点十七分，氮气纯度监测仪确实触发过一次0.001%的瞬时跌落，持续2.3秒。系统自动切换了备用气源，所以没生成告警工单。但日志里有记录。”
韩川握着笔的手指缓缓收紧，笔尖在纸上洇开一小团墨渍。
“把那段日志截图发我。另外，把跳变前后五分钟的所有环境参数——温度、湿度、震动频谱、甚至光照强度——全部打包。我要看它们和迁移率跳变的时间偏移量。”
挂断电话，韩川盯着那团墨渍，忽然笑了。不是释然的笑，而是某种近乎锋利的确认——就像数学家终于看清了隐藏在混沌表象下的对称群结构。
他打开MATLAB，新建脚本，第一行敲下：
`clear; clc; close all;`
第二行，他输入：
`addpath('C:\BNME_Data\Raw_IV_curves');`
然后停住。手指悬在回车键上方，迟迟未落。
夕阳正缓缓沉向湘江对岸，把整扇窗染成琥珀色。光带斜切过书桌，在《工程控制论》翻开的那页上投下锐利的阴影——恰巧横亘在“李雅普诺夫稳定性判据”的公式中间，将V?(x) < 0截成两半。
韩川忽然想起今天上午采访结束时，赵琳记者问他的最后一句话：
“很多人说您解题像在迷雾里点灯。但灯亮之后，您看到的究竟是出路，还是更多更深的迷雾？”
他当时笑着答：“迷雾从来不在外面，它就在我们定义问题的方式里。”
此刻，他按下回车键，屏幕泛起幽蓝微光。代码开始运行，光标在命令行里安静闪烁，如同黑暗隧道尽头，一盏刚刚擦亮的矿灯。
三分钟后，第一批数据读取完成。韩川调出Q7-LA-20130722样本的Ids-Vg曲线，将时间戳与氮气纯度跌落时刻对齐。他放大坐标轴，把横轴精度调到毫秒级——
就在纯度跌落起始点+1.7秒的位置，电流曲线上赫然出现一个宽度仅412毫秒的微小凹陷，振幅仅为基线的0.037%，此前所有分析软件都将其判定为仪器噪声而自动滤除。
但韩川知道，这不是噪声。
这是陷阱态集群被氮气中痕量氧分子短暂钝化时，释放出的集体“叹息”。
他关掉MATLAB，打开Word文档，新建一页，标题栏敲下：
**青鸾计划初步建模备忘录 · 韩川**
第一行，他写道：
> 迁移率跳变并非随机失效，而是材料内部陷阱态网络对外界微扰的协同响应。其本质是多尺度耦合动力学：宏观环境扰动（氮气纯度波动）→ 界面化学吸附平衡偏移（O?吸附/脱附）→ 局域电子结构重构（陷阱能级移动）→ 载流子输运路径重配置（迁移率突变）。
第二行，他停顿三秒，删掉“突变”二字，改为：
> ……载流子输运路径**可逆性重配置**。
第三行，他补充：
> 关键验证：若该机制成立，则在跳变发生后12小时内施加相同强度的氮气扰动，应观测到同步的迁移率响应——且响应相位滞后应恒定为1.7±0.2秒。
写完，他合上笔记本，起身走到窗边。
江风裹挟着水汽涌进来，吹散桌上那张画着刺状突起的草稿纸。纸角翻飞，露出背面一行铅笔小字，是昨夜睡前随手记下的：
*“所有看似无序的跳跃，都是系统在更高维度上遵循的确定轨迹。”*
韩川望着远处渐次亮起的路灯，忽然觉得那些光点连起来，很像一组正在自组织的拓扑相变临界点。
手机震了一下。
是袁荷蓓发来的微信，附着一张截图：氮气纯度监测日志。在时间轴2:17:03.1处，那道0.001%的凹陷清晰如刀刻。
下面跟着一行字：
【韩神，刚查完。跳变时刻的洁净间震动频谱——有个频率为3.71Hz的微弱峰，持续了2.8秒。我们昨天检修过通风管道，但没发现松动部件。】
韩川盯着那个数字看了很久。
3.71Hz。
他转身回到书桌前，拉开最底层抽屉，取出一本皮面笔记。翻开泛黄的纸页，某一页上用红笔圈着一个公式：
*ω? = √(k/m) → m ≈ k/(ω?2)*
旁边批注：CERN LHCb冷却系统共振频率反推质量缺陷值。
他拿起计算器，输入：k取典型有机薄膜杨氏模量1.2GPa，ω?=3.71×2π。
结果跳出：m≈0.087克。
不到一枚硬币的重量。
韩川轻轻吐出一口气。
原来迷雾深处，真有路。
只是这条路，要先拆掉自己脑子里的墙。