Rt-TDDFT揭示冷冻电镜中的保护机制

在模拟之前，首先需要构建不同温度下的平衡结构。为了合理地体现水分子环境对有机分子的影响，本研究构建了一个大小为14*14*14  A³的模型，其中包含1个碳酸乙烯酯（EC）分子和88个水分子。

该结构如图1（a）所示，水分子聚集在EC分子周围，密度约0.94g/cm³，是实验上冷冻样品中低密度无定形冰（LDA）的密度。为了获得合理的无序结构，先将该体系加热到1000 K（在这个过程中没有水分子断键），随后将无序结构分别退火到100 K（冷冻条件，接近液氮温度）和300 K（室温条件）。这样获得的结构的RDFs与Alexander等人的文章结果是一致的（图1（b））。

图 1（a）冷冻条件下模拟样品的结构；红色和白色部分表示水分子，棕色和黄色部分表示EC分子。（b）O原子之间的RDFs的比较。浅蓝色和红色实线表示本研究的结果，而星星和圆点是Alexander等人的结果。（c）电中性EC分子和整个模拟体系的DFT本征能量谱。（d）EC分子的最低能量本征态中的一些的波函数（编号为1、2、3、92、93和94），而那些较高能量本征态（C和O原子的p轨道电子）倾向于与水混合。

建立平衡结构之后，通过rt-TDDFT-MD计算来模拟分子解离过程。由于电子束辐照使一个小分子内产生多次激发的可能性较小，所以该研究只考虑单电子电离情况。

本研究在每个温度下进行了21次NOB模拟，每次模拟的初始条件和时间相同。图2是4次NOB模拟的结果，其中（a）-（d）处于冷冻条件下，（e）-（h）处于室温下。模拟开始时，位于价带的state-1被一个空穴占据。随着空穴跃迁到不同的能级，键长开始变化，4次模拟中均有三种键长显著增加（表明这种键的断裂）：C₁-O₂，C₁-O₃ 和 C₂-C₃（图 2（a）、（c）、（e）和（g））。其中，4次模拟中的断键时间和断键模式不同，证明了模拟的随机性。并且，当空穴从state-1跃迁到state-2 或state-3时，键长没有显著变化。直到跃迁到state-92，键长才开始急剧变化。上述结果表明，断键总是由WFC引起的，且通常发生在两个能量相差比较大的深能级之间。

图 2 左图（a，c，e，g）显示了 EC 分子中键长的变化。在100 K 下（a）和（c）随机数不同，在 300 K 下（e）和（g）随机数不同，导致不同的坍缩。EC 分子模型及其原子编号如（c）所示。右图（b，d，f，h）表示整个系统的绝热态电子占据和能级随时间的变化。线表示能级的变化，红点表示空穴的占据，点的大小与占据的空穴数成正比。

图 3（a）所有模拟中空穴跃迁的路径和分子解离的模式; 红点代表空穴的占据。左图列出了在两个温度下空穴从第一个区域跃迁到另一个区域的概率。右图显示了由于空穴跃迁到第三个区域而产生的分子碎片图像。（b）碎片的类型和发生键断裂的百分比。“Other”是指不太常见的片段。“Break”表示在200-fs模拟期间至少有一个键断裂的概率。

NOB模拟期间，空穴跃迁产生的能量会转移到分子碎片中，因此，碎片具有很大的动能。为了探究碎片的动力学过程，将rt-TDDFT-MD的最终原子位置和速度作为初始条件，进行了2 ps的AIMD模拟。

首先以单个“O”原子碎片的演变来看，在300 K下，“O”原子碎片向远离原始EC位置的方向移动，容易与水分子结合形成水合物，甚至将水分子的H原子吸附到自己身上；在100 K， “O”原子碎片与EC分子和其它碎片保持很近的距离，可以与其它碎片反应形成新的碎片，如图4（a）和（b）所示。

“笼”效应也发生在其它碎片上。例如，在100 K时，碎片"C₂-O₂"和"O₁- C₁-O₃-C₁"将结合形成长链，且具有类似于原始EC的弯曲形状，如图4（c）、（d）和（e）所示。但是，在2 ps AIMD中，长链并没有完全重组成EC。为了估计长链完全重组成EC的可能性，使用NEB方法计算了形成EC所需的反应势垒约1.5 eV。在2 ps AIMD结束时，分子具有800 K的温度。如果该温度保持恒定，则重组为EC的反应可以在大约0.003秒发生。

图 4 碎片的动力学过程。（a）（b）AIMD模拟中单个“O”原子碎片的路径，（a）发生在100 K，（b）发生在300 K；（c）（d）（e）EC分子在冷冻和室温条件下电离2 ps后的动力学。（c）电子束辐射导致键断裂，断键后如（d）；（e）显示了两种温度条件的差异，其中碎片在冷冻过程中由于“笼”效应而重组，在室温下没有“笼”效应而飞散。黄色原子表示EC分子中的O原子，棕色原子表示EC中的C原子，红色原子表示水分子中的O原子，白色原子表示水分子中的氢原子。EB是电子束（electron beam）的缩写。