闪光诱导多孔石墨烯COMSOL仿真

1天前浏览10

关键词：COMSOL；多孔石墨烯；EMI屏蔽；闪光灯；光热转换

一、文章简要介绍

无人机雷达、5G手机、可穿戴设备，都怕电磁干扰（EMI）。传统金属屏蔽罩又重又硬，不适合未来移动场景。这篇发表在Nano-Micro Letters 2023（DOI: 10.1007/s40820-023-01157-8）的研究，用闪光灯在几毫秒内把商用聚酰亚胺（PI）膜烧成了多孔石墨烯（FPG）——10×10平方厘米大面积、密度只有0.0354克每立方厘米、方阻低至18欧姆每方，绝对屏蔽效能SSE/t高达1.12×10的5次方分贝平方厘米每克，是目前碳基屏蔽材料里的最优水平。

图1是整体工艺概念：闪光灯宽谱（250到1100纳米）照射PI膜，触发三步协同光效应——先是紫外波段（能量超3.5电子伏）打断C-N键（键能3.04电子伏）制造微裂纹，然后可见-近红外波段的光热效应让膜升温，温度超过1700℃后非晶碳重排成有序石墨烯片层，同时CO2、N2、CO气体逸出撑出多孔结构。

二、仿真步骤

核心问题：闪光灯能不能把PI膜加热到1700℃以上？光靠实验没法直接测膜内温度，我们先用COMSOL Multiphysics 6.0热传导模块做理论验证。控制方程是瞬态热传导方程：密度×热容×温度对时间偏导，减去导热项，等于单位体积热生成率Q。

Q来自吸收的光能，按朗伯-比尔定律衰减：Q(x,y)等于（1减反射率）乘入射光强乘吸收系数，再乘指数衰减项。PI的透射率、反射率、吸收系数用紫外-可见光谱仪实测，材料的热物性参数（密度、热容、导热系数）直接取COMSOL材料库的温度相关数据。

模拟对象是125微米厚的PI膜，闪光能量密度从6到26焦耳每平方厘米扫描。网格划分与边界条件按闪光灯实际照射工况设置，输出膜内沿深度方向的温度分布和表面最高温度。

图2是模拟结果：8、18、22焦耳每平方厘米三种能量密度下的温度分布云图，随能量密度增大，高温区从表面向膜内扩展。

图3是量化曲线：26焦耳每平方厘米时表面最高温度达2300℃，热量有效传递深度约100微米。这个结果说明闪光灯完全能满足石墨烯合成所需的1700℃，理论可行性确认。

三、关键结果与发现

第一个发现：UV是"点火器"，Vis-NIR是"加热器"。实验对比全光谱与滤掉UV的Vis-NIR照射：只有全光谱才能形成多孔石墨烯，因为UV先在PI膜里制造缺陷，提高膜对可见-近红外的吸收，光热效应才得以充分累积——这就是论文说的"协同光效应"。

第二个发现：能量密度决定产物质量。图4的表征显示，能量密度从8升到22焦耳每平方厘米，碳原子占比从75.23%升到95.74%，氧氮含量大幅下降；Raman光谱出现2D峰（ID/IG从1降到0.8、I2D/IG升到0.38），说明形成了4到5层石墨烯；HRTEM看到3.3埃的晶格条纹，XRD 26度处强峰对应002面间距3.4埃，石墨化程度很高。

图4是XPS、FTIR、Raman、HRTEM和SEM的完整表征：C 1s主峰、C-N/C-O/C=O键的消失、石墨烯特征2D峰、3.3埃晶格条纹和约62微米厚的空心柱结构，共同证实了高质量多层石墨烯的生成。

第三个发现：屏蔽性能由吸收主导。图5显示方阻随能量密度从64降到11欧姆每方（22焦耳每平方厘米时为18欧姆每方），膜厚从26增到92微米。图6中SET从12.9升到34分贝，其中吸收项SEA从9.1升到27.8分贝、反射项SER始终低于7分贝——空心柱多孔结构延长了电磁波传播路径，内部多次散射吸收是主机制，这和金属靠反射屏蔽完全不同。

图5是方阻与膜厚随闪光能量密度的变化：导电网络逐步完善，为表面反射和内部吸收提供了电学基础。

图6拆解了屏蔽机制：有效吸收率Aeff从87.5%升到99.8%，反射系数R始终高于吸收系数A，说明先表面反射、再内部散射吸收的复合机制。

第四个发现：柔性耐久+应用落地。图7显示10000次弯折（半径4毫米）后方阻仅增大1.4倍、SE保持0.86倍，仍在20分贝以上。应用到无人机雷达：内部辐射电场降10.8分贝、24吉赫天线增益降12.2分贝；穿戴场景：人体SAR在5.2吉赫从7.7降到1.51瓦每千克（降80.3%）、5.5吉赫从8.5降到1.89（降77.8%），满足IEEE标准。

图7是EMI性能全景：弯折耐久曲线和SSE/t对比图，1.12×10的5次方分贝平方厘米每克的成绩在碳基材料里排第一。

工程启示：COMSOL热传导模拟在工艺开发里干的是"可行性预判"的活——先算清温度场，再定灯功率、脉宽、能量密度这些工艺窗口，比纯试错省太多时间。这套"实验实测光学参数+COMSOL热模拟+工艺参数扫描"的流程，做激光加工、闪光烧结、热退火工艺的同学可以直接复用。

来源：320科技工作室

化学空间对接十亿分子筛药

关键词：虚拟筛选；分子对接；ROCK1；激酶抑制剂；化学空间一、文章简要介绍传统虚拟筛选先枚举几亿个分子再逐个对接，算力开销巨大。这篇文章（Beroza等，2022，NatureCommunications，DOI:10.1038/s41467-022-33981-8）提出一种化学空间对接方法：不枚举完整库，而是把合成砌块直接对接进靶点口袋，用结构信息引导分子生长，实现近十亿分子的可扩展虚拟筛选。方法应用到ROCK1激酶：从858,125,390个虚拟产物中筛出69个分子购买测试，27个抑制活性Ki小于10微摩尔，命中率高达39%，两个先导化合物的对接姿势还得到X射线晶体结构证实。图1展示了80个起始片段在ROCK1结合口袋中的对接姿势：口袋表面为灰色，铰链区药效团为两个绿色球体，浅蓝色连接体作为反应向量用于后续产物库枚举，每个片段代表约一万个分子的库。二、仿真步骤步骤一：构建可搜索化学空间。从EnamineREALSpace双子库（2019年4月访问）出发，包含71,894个合成砌块和102个反应规则，理论最大产物数858,125,390。把每个砌块保留反应后剩余部分并加虚拟连接原子，得到136,835个片段。步骤二：片段对接与筛选。所有片段用FlexX对接软件对接到2ETR结构的ATP结合位点，施加铰链区药效团约束（配体氢键供体对接Glu321主链羰基，或氢键受体对接Met323主链氮），每个片段最多保留10个姿势，共129,125个姿势。用HYDE打分函数评估，选出打分前50,000的姿势导入SeeSAR人工检查，最终挑出500个最优片段及其姿势。筛选标准包括：额外氢键、分子量45到450、cLogP不超过4、连接体几何朝向、扭转能、化学多样性。步骤三：产物库枚举与重对接。按双子反应规则，每个选中片段作为起点全库枚举产物，500个库共生成5,236,824个产物，以片段姿势为模板用FlexX逐个对接，产生23,305,389个姿势，HYDE打分后保留10,391,986个。取每个分子最优姿势约33,000个，用Chemalot做应变能过滤（去掉内应变超过5千卡每摩尔的），剩5,940个分子，再用FRED程序重对接确认，最后k-means聚类成500簇，每簇取FRED打分最好的代表。步骤四：过滤与购买测试。500个簇代表经过两个定性过滤（过度柔性和大疏水基团），77个进入人工目检，最终从Enamine下单69个化合物，用HTRF激酶实验测定Ki值。三、关键结果与发现第一，命中率惊人：69个测试化合物中27个Ki低于10微摩尔，命中率39%，在对接文献中属于最高水平；最强化合物Ki达38纳摩尔，19%的分子达到亚微摩尔级。图2展示了四类化学型中最活跃分子的对接姿势：化合物1（吡唑类，橙色）、化合物16（吡啶酮类，品红）、化合物22（氮杂吲哚类，紫色）、化合物25（吲唑类，黄色）。所有命中都占据铰链区，并与激酶P环发生疏水相互作用，其中两个还与催化赖氨酸相互作用——这两种作用在模板配体中都不存在。第二，X射线证实对接姿势：对最强化合物1和结构新颖的化合物22解析共晶结构（分辨率2.34埃和2.74埃），对接姿势与晶体结构的RMSD分别为0.97埃和2.30埃，铰链区氢键原子偏差仅0.4到1.0埃，证实虚拟筛选预测的结合模式高度可信。图3是化合物1的晶体验证：上图是精修后的2Fo-Fc电子密度（1西格玛等高线），下图是对接姿势（浅绿）与X射线构象（深绿）的叠合，除中间苯环略有倾斜外，整体结合位置高度一致。图4是化合物22的晶体验证：氮杂吲哚的氢键原子仅偏移0.4到1.0埃，远端基团偏差较大主要来自蛋白质柔性（模板结构中Phe87侧链末端相差4.2埃），说明大部分位移源于蛋白构象变化而非对接误差。第三，计算效率大幅提升：传统全库枚举对接的耗时随化学空间规模线性增长，而化学空间对接的耗时只与合成砌块数量相关，比传统方法快多个数量级。图5对比了不同对接方案的计算需求：红黄绿曲线为传统全枚举对接（每分子10秒、1秒、0.1秒），浅蓝色为化学空间对接，文献中大规模对接战役标注为散点，可见本方法在近十亿分子规模下资源消耗远低于传统方案。图6给出了ROCK1结合位点与药效团约束的设置：结合口袋表面与铰链区约束（Glu321羰基或Met323氮），这是整个虚拟筛选流程的对接锚点，保证每个候选分子都形成激酶抑制剂必需的铰链氢键。来源：320科技工作室