论文部分内容阅读
本论文发展了有限元映射离散格点表象数值计算方法,模拟了阴离子的光电离过程,解释了气相小分子光剥离电子能谱的实验结果。通过对光电子能谱的研究,我们对分子的电子轨道、振转态和过渡态理论有更深入的了解。主要工作概括如下: (1)提出了可以高效求解薛定谔方程的有限元映射离散格点表象(OFE-DVR)方法。OFE-DVR是一种基于有限元离散格点表象(FE-DVR)的方法,它通过格点映射来优化FE-DVR的格点分布,使之在表征波函数以及求解薛定谔方程时效率更高。为了测试OFE-DVR方法的有效性,我们用它来求解小分子一维及二维薛定谔方程,并与其他被广泛使用的方法(包括FE-DVR、Sinc-DVR、Mapped-Sinc-DVR、MappedFourier等)进行比较。比较结果表明通过映射可以显著提高FE-DVR方法的计算效率。具体来说,在求解一维或者二维分子振动问题中,OFE-DVR的效率远远高于Sinc-DVR,与Mapped-Sinc-DVR的效率几乎相同;在求解涉及到库仑奇点的电子波函数问题中,OFE-DVR的效率远远高于其它几种方法。因此OFE-DVR有望在研究电子动力学问题中发挥重要的作用。 (2)计算Au2H-、Au2D-2的光剥离电子能谱,并解释了相关的实验结果。通过对Au2H-及其同位素取代物的光剥离实验结果的模拟和分析,我们获得了Au2H-、Au2H的电子结构和振转态信息,分析了Au在化学反应中所起的作用。为了计算光剥离电子能谱,我们首先构建了高精度的Au2H-和Au2H从头算势能面。然后利用构建的两个势能面,采用含时量子波包方法计算了Au2H-、Au2D-的光剥离电子能谱,其结果与实验数据吻合得很好。基于理论计算,我们对实验结果中的每一个峰都给出了详细的解释,并确定了Au2H-、Au2D-的绝热激发能和垂直激发能,给出了Au2H和Au2D的三个振动模式。 (3)采用含时量子波包理论计算了FH-2、 FD-2和FHD-的光剥离电子能谱,并与实验结果进行了比较。与计算Au2H-光剥离电子能谱不同的是,计算FH-2光剥离电子能谱主要目的是为了揭示过渡态附近的共振现象对化学反应的影响。为了计算FH-2及其同位素取代物的光剥离电子能谱,我们首先构建了高精度的FH-2从头算势能面。结合原有的高精度的FH2从头算势能面(CSZ势能面),通过用两种不同的方法(FTA方法和FF方法)计算阴离子平衡构型附近的振转态和中性分子本征态之间的弗兰克-康登因子(Frank-Condon Factor)来模拟FH-2及其同位素取代物的光剥离电子能谱。计算结果表明两种不同的方法都能准确地模拟FH-2及其同位素取代物的光剥离电子能谱,并且可以从不同的角度揭示光剥离电子能谱的产生机理,从而让我们对F+H2反应有更全面的理解。另外,我们还把不同势能面(CSZ、SW、FXZ、LWAL)上的计算结果与实验数据进行了全面的比较和分析。结果表明,尽管在这几个势能面上都能得到与实验吻合得很好的理论结果,但是CSZ势能面的计算结果与实验吻合得最好,因此CSZ势能面是目前描述F+H2及其同位素取代物反应最精确的势能面。