化学学报 ›› 2012, Vol. 70 ›› Issue (17): 1869-1876.DOI: 10.6023/A12040139 上一篇    

研究论文

改进的半经典动力学模拟二苯乙烯光致顺反异构化反应

雷依波a,b, 朱超原b, 文振翼a,b, 林聖聖b   

  1. a 合成与天然功能分子化学教育部重点实验室 西北大学化学与材料科学学院 西安 710069;
    b 西北大学现代物理研究所 西安 710069
  • 投稿日期:2012-04-21 发布日期:2012-08-06
  • 通讯作者: 朱超原
  • 基金资助:

    项目受国家自然科学基金(Nos. 21003100, 21033001, 21103136, 21173166)资助.

New Implementation of Semi-classical Dynamic Simulation on the Photoisomerization of cis- and trans-Isomers of Free Stilbene

Lei Yiboa,b, Zhu Chaoyuanb, Wen Zhenyia,b, Lin Sheng-Hsienb   

  1. a Key Laboratory of Synthetic and Natural Functional Molecule Chemistry of Ministry of Education, The College of Chemistry & Materials Science, Shaanxi Key Laboratory of Physico-Inorganic Chemistry, Northwest University, Xi'an 710069;
    b Institute of Modern Physics, Northwest University, Xi'an 710069
  • Received:2012-04-21 Published:2012-08-06
  • Supported by:

    Project supported by the National Natural Science Foundation of China (Nos. 21003100, 21033001, 21103136, 21173166).

发展了一种改进的半经典动力学模拟方法, 并将其程序化用于气相二苯乙烯光致顺反异构化反应的机理研究.新的方法不仅采用e 指数模型改进了原有Zhu-Nakamura 理论中计算电子非绝热跃迁几率的计算方法, 而且将约束哈密顿方法用于限制性分子动力学模拟过程中. 计算结果表明, 采用此方法得到的统计平均的量子产率及反应机理与以前的实验与理论结果吻合较好, 从而可以应用于全量子动力学方法无法进行的大分子体系的动力学研究.

关键词: 改进的半经典动力学模拟, 约束哈密顿系统, Zhu-Nakamura 理论, 二苯乙烯顺反异构化, 二维解析势能面

New implementation of semi-classical trajectory surface hopping dynamic simulation has been developed and applied to the photoisomerization of cis- and trans- isomers on the gas phase. This method not only uses the exponential model to the modification of the originally analytical non-adiabatic transition probability formula, but also involves the constrained Hamiltonian system into the constrained molecular dynamic simulation. Two-dimensional potential energy surfaces of ground S0 and excited S1 states are constructed analytically fitting to ab initio calculations in terms of torsion angle and one dihedral angle around the central ethylenic C=C bond as variables, and the other internal coordinates are all fixed at configuration of one-bond flip conical intersection. The analytical PESs are quite accurate and the mean absolute error is less than 2.4 kcal·mol-1, and much less than 1.0 kcal·mol-1 around conical intersection region. A straight seam line is found on potential energy surfaces that simply separates the cis-area with the trans-area. The constrained Hamiltonian system is employed to run trajectories in the Cartesian coordinate system and surface hopping in terms of the two internal dihedral angles. Typical trajectories are found in which the torsion angle changes monotonically for both cis- to trans- and trans- to cisisomerizations. This is an exact picture of one-bond flip mechanism of photoisomerization around the conical intersection. Quantum yield for trans- to cis- isomerization is simulated as 60.45% in very good agreement with experimental value 55.0%, while quantum yield for cis- to trans- isomerization is simulated as 42.3% in comparison with experimental value 35.0%. As the S1 energy in local minimum of cis-area is higher than that in trans-area, and thus cis- to trans- isomerization is quite possible to access to another Hula-Twist conical intersection. These simulation results demonstrate that the computed cumulative quantum yield and reaction mechanism are consistent with the previously experimental and theoretical results. This means that the present trajectory surface hopping method would be good at the dynamic simulation on the large system with or without constraint Hamiltonian in comparison with the quantum molecular dynamics.

Key words: new implementation of semi-classical dynamic simulation, constrained Hamiltonian system, Zhu-Nakamura theory, isomerization of cis- and trans-stilbene, two dimensional analytical potential energy surfaces