给H2+OH化学反应拍一部电影，了解现代化学的神奇之处。
化学给我的第一印象源自电解水：记得初三的时候，黄老师现场给我们演示这个实验。老师把两个通了电的夹子连接上水中的两根电极，在连接上的瞬间，两个电极上顿时都密密麻麻的布满了气泡，气泡不停的往上冒。黄老师用一个气球收集了其中的一种气体，当他把气球口系紧松手后，气球迅速逃离了他的手，飞到天花板上。黄老师说这个气球中是氢气。
他还用一个倒扣的塑料杯子同时收集了一些实验中产生的两种气体，然后倾斜的压在一支笔上。当他把一根点燃的火柴靠近杯口的时候，突然发出了一声巨大的声音。虽然我们事先已经被告知，但坐在第二排的我还是被吓了一跳。当时，我对这个实验有点害怕，但又好奇：为什么在一个电极产生的是氧气，在另一个上产生的是氢气。如果把H2O中的H和O原子重新组合，那么H2和O2应该在同一个地方产生，不是吗？
其实那时候的我并不理解其中具体的过程，黄老师也并未提及他在水里添加了电解质。
最近，跟我的一个侄女讨论了一下甲烷的燃烧。她说这个化学反应很简单，她们初中的时候就已经学明白了，可以用一个简单的化学反应方程式来描述。她还特意强调说反应式配平了。当我追问CH4、O2、CO2和H2O长什么样子的时候，她开始有点怀疑自己对这个反应的认知了。然后，我跟她说，这个反应非常复杂，但她却不相信。
于是我给她看了Nature communications上最近的一篇文章中的图。我跟她介绍说，这些仅是甲烷燃烧初始阶段的一些化学反应。目前，科学家在甲烷燃烧过程中共发现了798种不同的化学反应。听到这些，她彻底否定了自己原先对这个反应的认知，说这个CH4已经不是她原先认识的那个C H4了。
其实，我并不是想要介绍电解水或者甲烷燃烧，而是想通过H2+OH→H2O+H这个看似简单的反应来介绍一下量子化学反应动力学，以此来介绍现代化学中的一些基本概念。虽然这个反应并不常见，但研究这个反应却非常重要。这个反应涉及到的原子种类跟电解水一样，只有氢元素和氧元素。这个反应也参与到了甲烷的燃烧过程。
传统化学和现代化学我们熟悉的化学大都是在试管、烧杯或者烧瓶等容器中进行的，这是传统化学。
而现代化学却截然不同。现代化学可以精细的操控反应物，让单个反应物分子相互碰撞，从而不仅研究化学反应产物，还研究化学反应的中间过程。现代化学的实验通常使用分子束的实验装置，让分子排列着共同朝着一个方向运动。这种实验还需要在超级真空的环境下进行，也就是把所有不相关的分子都赶出门外，让一对分子单独相处然后观察他们之间的相互作用。
这是现代化学的一种研究方法，还有另一种方法，叫作理论计算化学，是用计算机来模拟化学反应。这样的计算一般都比较复杂，往往需要用几十台甚至几百台的计算机共同来完成。
在2011年的时候，现代化学对HD+OH→H2O+D的反应进行了精细的研究。在这个研究中，研究者们让HD分子束和OH分子束在真空环境中进行碰撞，然后查看D产物在空间的分布情况。实验和理论的结果一模一样。通过分析产物的空间分布，研究者们可以推测出HD和OH是如何进行碰撞的，比如在反应过程中，HD和OH之间相对角度是什么样子的，还可以研究HD和OH之间相对速度如何影响反应。
上面介绍的两种现代化学的研究方法，可以让我们非常精确的去了解一个化学反应的详细过程。但从产物的空间分布来了解反应的中间过程是一种推测，而不是直接的观测。要是能够把分子的碰撞过程拍成一个电影，那么我们就可以清楚地知道化学反应的整个过程了。但是，实验上观测这样的微观过程受到很多限制，而理论计算方法可以给我们一个完整的视角。在理论模拟的过程中，我们可以很方便的观察化学反应的整个过程。
我们甚至可以成为量子世界的一名导演，让分子成为演员，给他们穿上不同的衣服，让他们站在不同的高度和位置，然后让他们按照剧本去演绎一个别开生面的景象。这正是接下来我想具体介绍的内容。首先，我需要介绍一下如何成为一个量子世界的导演；然后再介绍我们制作的这部量子世界的电影；最后总结一下我们从中发现了什么重要的现象。
如何成为一名量子世界的导演在我们生活的这个宏观世界里，要成为一名导演，需要经历一个学习的过程，掌握其中的技巧和方法。同样的，在量子世界做导演也需要一些基础技能的学习。我们“拍”这部量子世界的电影，目的是揭示微观世界中分子之间的相互作用，所以我们需要从最基本的原理出发，不使用任何实验参数而只能用基本的物理常量。这样的理论可以很接近的反映宇宙本质，被称为第一性原理。
这样的理论方法十分重要，可以独立揭露物理化学现象的本质。第一性原理中很重要的一条是量子力学，用于描述物质世界中微观粒子的运动规律。就比如在电影中，剧情的发展由人与人的关系、社会法则等来推动，而微观粒子的运动则是由量子力学来描述的。有一个被称作哈密顿量的物理量，它规定了在每个时间各个演员（微观粒子）之间如何“交往”，还有一个法则是“薛定谔方程”，它描述微观世界如何演化发展。
微观粒子的运动特点跟我们日常生活中观察到的物体有很大的区别。日常生活中，我们观察到的物体在某个时间点上只存在于一个位置上。它的运动情况可以准确地用牛顿三大运动定律来描述。然而，微观粒子的运动不具有这个特性。在某个时间点上，微观粒子可以同时存在于空间的多个位置，在不同的位置都有一定的分布概率。了解了这些技能和方法之后，我们对量子世界的特点多了一些了解，可以成为一名正式的量子世界的导演了。
着手拍摄这部量子世界的电影之前，我们需要选剧本、选演员、布置场景。一切就位后，才可以正式拍摄。
一部量子世界的电影研究OH(A)和H2的碰撞反应不仅可以揭示化学反应的本质，也具有重要的现实意义，可以加深对燃烧化学、大气化学和星际化学的认知。比如，很多化学反应中都会产生OH(X)产物。
为了检测OH(X)的产物，经常需要用激光把OH(X)激发到OH(A)这个态，然后检测其产生的激光诱导荧光谱来判断OH(X)的产物。但是，当OH(A)跟某些分子或者原子进行碰撞后，会让OH(A)的电子态激发淬灭，从而降低了荧光谱的强度，造成不准确的检测结果。研究这个反应有助于准确的了解OH(A)淬灭的具体过程。
OH(A)和H2碰撞的结果共有三种可能性：在反应淬灭通道，H2中的H-H键断裂，其中一个H原子和OH结合形成水分子H2O，还留下单独的一个H原子。在非反应淬灭通道，OH(A)和H2碰撞后没有发生化学键的断裂，而是使得OH(A)和H2跳到底下楼层。第三个弹性和非弹性散射通道中，OH和H2在碰撞后依旧停留在原来楼层。
美国宾夕法尼亚大学的Marsha Lester教授课题组对OH(A)和H2碰撞过程中的反应淬灭和非反应淬灭通道的比例进行了详细的研究。他们发现，在他们的实验条件下，反应淬灭通道占了88%，而非反应的淬灭通道仅为12%。随后，不同的理论课题组对这两个通道的比例做了详细的研究。
美国Emory大学的Joel Bowman教授课题组和国内大连化学物理研究所的韩克利教授课题组得出的结论跟实验的一致；而澳大利亚国立大学的Michael Collins教授和国内大连化学物理研究所的张东辉教授在一篇合作的理论文章中，他们得出的结论跟实验恰好相反。这样的不一致让实验研究和理论研究都产生了担忧，不知道具体是什么原因造成的分歧，足足困扰了大家十多年时间。
我们这部电影（研究论文发表在2021年9月份的Nature Chemistry）讲述的是化学反应，参与化学反应的是分子，所以分子是我们的演员。首先，我们来介绍分子演员：它由原子组成，原子又由原子核和电子组成。分子中，原子之间通过化学键相互连接，而化学键是电子在原子核周围形成的一种特定的空间分布，这种空间分布把原子核连接起来，形成化学键。化学反应的过程涉及到旧化学键的断裂和新化学键的形成。
分子中的电子像一团云笼罩在原子核周围。由于电子的质量比原子核小几千或者几万倍，即使原子核处于运动之中，电子总能迅速的调整电子云的分布。相对于电子的迅速调整，原子核似乎是静止不动的。原子核和电子运动的这个特点，让我们能够分别来计算电子和原子核的运动。这就是在理论计算化学中经常采用的“玻恩-奥本海默近似”。
但是在某些特定的原子核位置，电子和原子核运动之间会产生强烈的相互影响，这正是OH(A)+H2碰撞中精彩的部分。具体将在下面详细介绍。
为了拍摄H2+OH化学反应的电影，在电影剧本（哈密顿量H(r,t)）中，我们需要描述四个原子之间的相对位置：H2中两个H之间的距离，OH中O和H之间的距离，还有H2和OH之间的距离，以及H2和OH之间的相对角度。
接下里我们需要来构造拍摄电影的场景，这样的场景和场地布置由电子云的能量决定。在一个固定的H2和OH位置下，都有一系列的电子云分布。当原子核之间的位置发生改变的时候，电子云会随着进行调整，有时候电子云会把原子拉近，而有时候会把原子核推开。我们的场景每一层楼对应的是一种电子云的分布。
但是，这个场景跟我们生活中的楼房有很大的不一样：这里的每一层楼都有六个方向，也就是上面图中描述原子之间相对位置的三个距离和三个角度。每一层楼并不是平的，而是凹凸起伏的。有些地方，两层楼恰好一样高（能量相同）。在图的右侧，OH和H2之间离得比较远，他们之间没有相互作用。当OH(A)和H2逐渐靠近（R逐渐变小），OH(A)和H2的电子云产生重叠，其效果是把OH(A)和H2更进一步的拉近。
但是，在下面的OH(X)和H2逐渐靠近的时候，它们是彼此排斥的。在OH(A)和H2相隔某些距离的时候，第三层楼恰好跟下面两层碰到一起。图中的蓝色和黑色粗线分别标注了第三层与第二层以及第二层和第一层之间的交汇处。在这些交汇处，原子核可以选择受到哪一种电子云的束缚。图中还画了两个具体的OH和H2的相对位置，分别叫做T字形CI和线性CI。其中线性CI的能量低于T字形CI的。
通常，原子核倾向于能量低的电子云，也就是会跑到低楼层。另外，跑到低楼层后，原子核的运动方向可以有两种选择，往左或者往右，它们将使得反应进入不同的产物通道。
到这里，我们总算可以正式拍摄这部量子世界的电影了。第一场拍摄的是第三层楼上的故事，OH(A)和H2两个分子演员分别处于振动和转动的基态。我们首先轻轻的推一下OH(A)和H2，让它们相互靠近。
在相互靠近的时候，它们彼此相互吸引，然后进一步的相互靠近。但是，它们相互靠近需要按照几个特定的角度进行。在视频中，我们可以看到OH(A)和H2在刚开始的时候，OH(A)和H2在空间各个方向的分布具有相同的可能性，这是因为它们都处于转动基态。随着时间流逝，我们知道OH(A)和H2在逐渐靠近。在2100 au的时候，有一部分的OH(A)和H2聚集在了一定的位置。
在整个过程中，T字型的CI位置处自始至终都没有出现聚集。然而，在图中我们介绍过，T字型CI的能量也远远低于初始时候OH(A)和H2的能量，并且在OH(A)和H2靠近的时候，也可以吸引它们相互靠近。但是，在这个过程当中，能量更低的线性CI显然更具有吸引力，在跟T字形CI竞争的过程中明显占了上风，把所有到达这个区域的OH(A)和H2都调整到特定的相对位置。
最后，我们可以完整的播放这部量子世界的电影。
这个电影展示了这个反应的整个过程。我们可以发现很多有意思的物理现象。接下来，我们将对OH(A)和H2的碰撞做进一步的说明。首先，通过介绍反应中的立体动力学，我们发现弹性和非弹性散射通道的贡献，进而发现实验结果的分析过程忽视了这个通道，从而最终解决一个存在了很久的分歧，即反应淬灭和非反应淬灭通道的比例。
最近，我们为OH(A)和H2构造了一个全维的高精度的势能面，考虑了体系的四个电子态对碰撞过程的影响。我们进而在这个势能面上开展了全维的量子碰撞动力学研究，发现了弹性和非弹性散射通道的存在。前几年，英国牛津大学的Mark Brouard教授在实验上也确认了这个通道的存在。
但是，Lester教授在分析他们实验结果的时候，并没有考虑到这个通道的存在，从而意外的把所有弹性和非弹性散射部分归结成了反应淬灭通道，从而造成这个通道占比很大。一旦考虑这个弹性和非弹性散射通道，我们的理论计算结果跟Lester教授的实验结果吻合得非常好。
多年以来，我父亲经常询问我工作的具体内容。我曾尝试用简单的几句话来跟他描述，但似乎都没有成功，因为他过了一段时间后依旧会再次询问。对此我感到十分愧疚。这次恰好需要写一篇文章来介绍我们最近的一个工作，于是想趁次机会多写一点，希望能让我的父亲多了解一点我的工作。文章中一些不太准确的类比，望读者多多包涵。