收藏本站
《西北大学》 2009年
收藏 | 手机打开
二维码
手机客户端打开本文

激光诱导光化学反应模拟与半经验MRCI及其解析梯度程序化研究

雷依波  
【摘要】: 首先我们概述了分子动力学的理论基础,包括波恩-奥本海默近似(或绝热近似)、半经典近似及两种非绝热动力学近似方法—Ehrenfest动力学和系间窜越(SurfaceHopping近似)。在此基础上我们介绍了一种半经典非绝热动力学模拟方法—半经典的电子-辐射-离子动力学(SERID),此种方法具有以下几个特点: 1.在电子Hamiltonian中通过时间有关的Peierls替代引入了辐射场的矢势A(x,t),即在模型中引入了激光与电子的相互作用。 2.动力学过程采用半经典模型。电子运动用电子波函数描述,而电子波函数的更新需要求解含时Schr(o|¨)dinger方程。核运动轨迹借求解牛顿运动方程获得,求解时采用了一种辛算法(velocity Verlet algorithm),这一算法能够保持能量,动量和几率守恒,并满足Pauli原理。 3.电子结构采用半经验的DFTB(density functional tight-binding method)方法计算。DFTB计算提供了模拟过程中体系的各种能量的变化,电子在轨道之间的跃迁等。 4.是一种直接动力学计算方案,即不需预先构建势能面,能量和力的计算采用了即用即算(on the fly)方案。 在这一章我们还简要剖析了SERID程序(由我们的合作者豆育升教授提供),给出了程序的流程图。并指出了该程序的优点和某些不足之处。 其次我们采用SERID方法模拟了环丁烷光解离生成两个乙烯分子的非绝热动力学过程。模拟结果表明环丁烷光解离是一个两步过程,首先环丁烷一个C-C键断裂形成一个四亚甲基中间体,接着这个中间体解离生成两个乙烯分子。这个结果与以前报道的实验结果一致。为了将模拟结果与以前报道的ab initio结果进行对比,我们采用CASSCF/MRPT2(complete active space self-consistent field/multi-reference secondorder perturbation)方法沿着模拟的反应路径构造了解离过程中基态和低能激发态的势能曲线,还构造了给定C—C键长后各电子态的最小能量路径。这些计算结果表明,环丁烷光解离过程中四亚甲基中间体确实存在,但是却出现在1~3A的势能曲线上,因为此时1~3A的能量比基态的低,而且基态与1~3A的势能曲线发生了两次相交,其中第一次相交对应于四亚甲基中间体的形成,而经过第二次相交后,1~3A的能量比基态的高,从而体系又回到了基态。同样我们采用SERID方法模拟了腺嘌昤非辐射去活化的非绝热动力学过程。应用两个不同波长的激光模拟得到了两条不同的反应路径,分别对应于氨基的平面外扭转和C_2-H键的平面外扭转。为了与以前报道的ab initio结果进行对比,我们采用CASSCF方法沿着两条模拟的反应路径构造了解离过程中基态以及能量较低激发态的势能曲线。结果表明两条反应路径中都存在第一激发态(~1ππ~*(~1L_α))与基态(S_0)的相交,导致电子从激发态回到基态,即活化过程中是从~1ππ~*(~1L_α)到S_0的衰减过程,这与前人的ab initio计算结果不同;他们的计算认为氨基平面外扭转的反应路径是从~1nπ~*到S_0的衰减过程。另外通过对各电子态主组态的分析,我们得到的激发态寿命与实验结果基本一致。 鉴于SERID方法中DFTB的缺点,我们拟采用半经验MRCI(multi-referenceconfiguration interaction)方法来改进电子结构计算。为了将改进方法应用于非绝热动力学中,半经验MRCI的解析梯度也是重要的。因此在第四章中我们概述了半经验MRCI及其解析梯度的理论基础。首先简要介绍了半经验方法的基本原理及发展现状,其中重点介绍了NDDO(neglect of diatomic differential overlap)方法及其参量改进方法。然后重点讨论我们的基于图形酉群的MRCISD算法。我们讨论了酉群不可约表示的基与自旋匹配的组态函数(CSF)的对应关系;给出了Shavitt建议的用于记录CSF的不同行表(DRT)的构造方式;描述了如何搜寻DRT中的Loop,从而计算生成元和生成元乘积的矩阵元以及哈密顿矩阵元。在此基础上简要介绍了最新发展的基于空穴粒子对应的MRCISD算法及其近似算法—双收缩CI(DCCI)的基本原理。在空穴粒子对应的MRCISD算法中,我们重新定义了不同行表(DRT),将空穴空间域外空间的完成Loop和未完成Loop预先计算出来,使得不同行表只包括活性空间的子DRT,这样大大降低了Loop搜索的计算量。作为近似算法,DCCI充分应用了空穴粒子对应的原理,将空穴空间和外空间都收缩为单一步矢,这样大大降低了参与变分的CI系数的数目。在本章的最后,我们介绍了半经验MRCISD解析梯度的理论基础。其中解析梯度的计算可以分为两部分,第一部分计算积分梯度的贡献,称为静态部分;第二部分计算分子轨道系数梯度的贡献,称为响应部分。静态部分需要CI约化密度矩阵与积分梯度相结合,响应部分需要计算Lagrangian矩阵元与求解CPRHF(restricted coupled perturb Hartree-Fock)方程。基于此我们讨论了半经验MRCI及其解析梯度程序的实现过程。半经验MRCI程序将MOPAC 7.0与我们小组发展的XIAN-CI程序包中的基于空穴粒子对应的MRCISD和DCCI程序相结合。从MOPAC 7.0中的MNDO(modified neglect differential overlap)计算得到原子轨道积分及分子轨道系数,完成积分变换得到CI计算所需的分子轨道积分,接着采用MRCISD或DCCI程序计算CI能量。计算过程中我们根据冻结轨道的特点,分开计算冻结轨道和CI轨道对能量的贡献以减少计算量。在半经验MRCI解析梯度程序的实现过程中,首先需要从CI系数和耦合系数确定约化密度矩阵,接着需要将MNDO计算得到的原子轨道积分梯度及分子轨道系数作为输入,计算解析梯度静态部分的贡献;然后计算Lagrangian矩阵元和Z向量,最后计算解析梯度响应部分的贡献。其中CI密度矩阵的计算量最大,直接决定计算瓶颈及计算效率。我们不仅采用直接CI的算法避免了传统CI中需要保存耦合系数或哈密顿矩阵元所带来的计算瓶颈,而且应用空穴粒子对应提高了计算效率。同时针对冻结轨道的特点,应用了许多简化算法以提高计算效率。测试结果显示,我们的半经验MRCI及其解析梯度程序无论从计算规模还是计算效率上都比以前报道的程序优越。
【学位授予单位】:西北大学
【学位级别】:博士
【学位授予年份】:2009
【分类号】:O413.1

手机知网App
【引证文献】
中国硕士学位论文全文数据库 前1条
1 姜富灵;AP与HMX、RDX反应机理的理论研究及UA光致异构化反应的动力学模拟[D];西北大学;2010年
【参考文献】
中国期刊全文数据库 前1条
1 甘正汀 ,苏克和 ,王育邠 ,文振翼;A method to fast determine the coupling coefficients in CI calculation[J];Science in China(Series B);1999年01期
【共引文献】
中国期刊全文数据库 前10条
1 豆育升 ,文振翼 ,刘朝萌 ,王育彬 ,张棣;配体场强场近似的酉群计算及程序化[J];高等学校化学学报;1989年12期
2 朱卡克,甘正汀,苏克和,王育彬,文振翼;CH自由基和NO_2反应研究——Ⅰ.反应的热力学计算[J];化学学报;2000年05期
3 王育彬 ,甘正汀 ,苏克和 ,文振翼;An adjustable contracted CI method[J];Science in China(Series B);1999年06期
4 王育彬,甘正汀,苏克和,文振翼;Configuration-based multi-reference second order perturbation theory[J];Science in China(Series B);2000年06期
5 甘正汀,苏克和,王育邠,文振翼;组态相互作用计算中耦合系数的快速确定[J];中国科学(B辑);1998年05期
6 王育彬,韩慧仙,翟高红,索兵,文振翼;双收缩CI方法及其应用[J];中国科学(B辑 化学);2004年02期
7 宋士吉 ,吴澄;Reverse triple Ⅰ method of fuzzy reasoning[J];Science in China(Series F:Information Sciences);2002年05期
8 ;Approximation and universality of fuzzy Turing machines[J];Science in China(Series F:Information Sciences);2008年10期
9 陈志翔,卢振洋,殷树言,宋永伦;焊缝位姿及焊枪位姿的模型[J];机械工程学报;2003年07期
10 王书宽;;非线性动力学与混沌理论在经济学中的应用前景[J];生产力研究;2006年07期
中国重要会议论文全文数据库 前1条
1 冯承天;;示性类的基本理论[A];数学·物理·力学·高新技术研究进展——2000(8)卷——中国数学力学物理学高新技术交叉研究会第8届学术研讨会论文集[C];2000年
中国博士学位论文全文数据库 前6条
1 苏克和;模型化学精确计算方法研究[D];西北大学;2000年
2 翟高红;H_3O~+(T)的完整势能面计算及空穴-粒子对应在MRCISD中的应用[D];西北大学;2003年
3 程知言;浅表隧道工程多冷源冻结温度、应力、水分场耦合研究[D];中南大学;2003年
4 翟高红;H_3O~+(T)的全域势能面计算及空穴-粒子对应在MRCISD中的应用[D];西北大学;2003年
5 索兵兵;空穴一粒子对应的MRCI程序及其并行化研究[D];西北大学;2005年
6 阮东;群表示论中不可约张量基方法的研究[D];清华大学;1997年
中国硕士学位论文全文数据库 前9条
1 蔡玉平;位移型相变的对称性与序参量[D];河北工业大学;2000年
2 魏俊;G3理论中芯电子相关能及多参考态精确算法研究[D];西北工业大学;2001年
3 宋军锋;关于广义Hamilton系统和Poisson流形上一些问题的研究[D];新疆师范大学;2004年
4 韩慧仙;S_2O分子振动光谱的计算及双收缩CI方法应用[D];西北大学;2004年
5 彭谦;B和Be有关化合物的电子结构及其势能面的理论研究[D];西北大学;2004年
6 于峰峰;可控非完整系统与机电动力系统的对称性和守恒量[D];青岛大学;2006年
7 胡月宏;有关辛李群和微分动力系统的讨论[D];新疆师范大学;2006年
8 侯传燕;Poisson结构和仿射群的拓展及应用[D];新疆师范大学;2005年
9 姜富灵;AP与HMX、RDX反应机理的理论研究及UA光致异构化反应的动力学模拟[D];西北大学;2010年
【同被引文献】
中国期刊全文数据库 前10条
1 王爱英;陈树森;李丽洁;金韶华;孔德斌;;不同量子化方法计算β-HMX分子结构的比较[J];火炸药学报;2007年02期
2 刘子如;施震灏;阴翠梅;赵凤起;;热红联用研究AP与RDX和HMX混合体系的热分解[J];火炸药学报;2007年05期
3 陈煜;刘云飞;谭惠民;;NEPE推进剂的细观力学性能研究[J];火炸药学报;2008年01期
4 戴大地,黎乐民;镱硫属化合物的密度泛函理论研究[J];高等学校化学学报;1997年06期
5 戴大地,黎乐民;几种密度泛函理论公式用于镧系硫属化合物计算的比较[J];高等学校化学学报;1997年07期
6 刘子如,阴翠梅,孔扬辉,赵凤起,罗阳,向海;高氯酸铵的热分解[J];含能材料;2000年02期
7 肖继军,张骥,杨栋,肖鹤鸣;环四甲撑四硝胺(HMX)结构和性质的DFT研究[J];化学物理学报;2002年01期
8 姬广富,肖鹤鸣,董海山;β-HMX晶体结构及其性质的高水平计算研究[J];化学学报;2002年02期
9 肖继军,姬广富,杨栋,肖鹤鸣;环三甲撑三硝胺(RDX)结构和性质的DFT研究[J];结构化学;2002年04期
10 刘子如,阴翠梅,孔扬辉,赵凤起,罗阳,周琥;高氯酸铵与HMX和RDX的相互作用[J];推进技术;2000年06期
中国硕士学位论文全文数据库 前2条
1 陆裕平;RDX的电子结构和热解机理及Ar-N_2,H_2O-H_2O分子间相互作用势的从头算研究[D];四川大学;2005年
2 王新锋;气相RDX热分解及溶剂影响的理论研究[D];中国工程物理研究院;2005年
【二级引证文献】
中国硕士学位论文全文数据库 前1条
1 徐文;RDX和多元硝酸酯的热解及AKⅡ、C2水解和安定机理的理论研究[D];西北大学;2011年
【相似文献】
中国期刊全文数据库 前10条
1 薛鲍;;西北大学专利事务所走出一条适合高校专利事业发展的新路[J];西北大学学报(自然科学版);1992年02期
2 姚远;物理学家岳劼恒[J];西北大学学报(自然科学版);1994年03期
3 朱恪孝,马林安,杨晶;综合大学优势与基础科学研究[J];西北大学学报(自然科学版);1997年06期
4 程馨莹;长江学者:西北大学特聘教授舒德干[J];科技.人才.市场;2000年01期
5 ;分手并不那样痛[J];新知客;2007年12期
6 高立勋;;西北大学与金蝶联合成立中间件技术中心[J];西北大学学报(自然科学版);2008年01期
7 薛鲍;;日本COSMO石油公司持续支持西北大学生命科学学院环保项目[J];西北大学学报(自然科学版);2011年03期
8 姜洪训;西北大学地质系成立四十五周年学术讨论会在西安隆重举行[J];石油与天然气地质;1985年01期
9 赵弘毅;西北大学现代物理研究所介绍[J];西北大学学报(自然科学版);1986年03期
10 ;刊首语[J];西北大学学报(自然科学版);1989年04期
中国重要会议论文全文数据库 前10条
1 ;西北大学毒理学史研究室[A];毒理学史研究文集(第六集)[C];2006年
2 ;西北大学毒理学史研究室简介[A];毒理学史研究文集(第八集)[C];2009年
3 文振翼;王育彬;索兵;甘正汀;;Xian-CI程序及应用[A];中国化学会第九届全国量子化学学术会议暨庆祝徐光宪教授从教六十年论文摘要集[C];2005年
4 ;西北大学生态毒理研究所[A];第一届中国古代毒物学史研讨会论文集[C];2003年
5 胡正海;;从事结构植物学教学和科研工作50年的回顾[A];中国植物学会植物结构与生殖生物学专业委员会、江苏省植物学会2007年学术年会学术报告及研究论文集[C];2007年
6 詹克明;朱幼南;袁建明;蒋晓马;;用密度矩阵方法计算耦合的四核交换体系动态核磁共振(DNMR)线型[A];第六届全国波谱学学术会议论文摘要集[C];1990年
7 张庆营;;非中心力矩阵元的公式(续)[A];第五次核物理会议资料汇编(中册)[C];1982年
8 ;西北大学毒理学史研究室简介[A];全国第二届毒理学史与突发中毒事件处置研讨会论文集[C];2004年
9 幸翀;王楠;刘仁桥;;纳米级散射颗粒溶液的Mueller矩阵测量[A];中国遗传学会第十届全国激光生物学学术会议论文摘要集[C];2009年
10 ;其它[A];中国地质科学院南京地质矿产研究所文集(40)[C];1989年
中国重要报纸全文数据库 前10条
1 姚远;杨钟健出长国立西北大学[N];光明日报;2002年
2 郝京;认真做事 全面发展[N];陕西日报;2008年
3 孟西安 王小娟;西北大学学生党员达14%[N];人民日报;2004年
4 外事;香港新鸿基地产郭氏基金出资780万元资助贫困生[N];各界导报;2007年
5 通讯员 赵亚杰赵芳;西北大学专家考察指导礼县先秦陵园规划工作[N];陇南日报;2007年
6 记者施耀花;我市哲学、社会科学课题进行首次公开招标[N];白银日报;2009年
7 席玮;作家能否做大学教授[N];中国社会科学报;2009年
8 本报记者 王鸣琦;高教大省面临的机遇和挑战[N];陕西日报;2004年
9 记者  杜朋举 实习生 徐丽;启动党的建设共建活动[N];陕西日报;2007年
10 崔延力王鸣琦;西北大学设立辅导员队伍建设长效机制[N];陕西日报;2007年
中国博士学位论文全文数据库 前10条
1 雷依波;激光诱导光化学反应模拟与半经验MRCI及其解析梯度程序化研究[D];西北大学;2009年
2 索兵兵;空穴一粒子对应的MRCI程序及其并行化研究[D];西北大学;2005年
3 刘邦贵;超弦理论[D];西北大学;1988年
4 王萃娟;基于新型吡啶配体的分子构筑及功能配合物的研究[D];西北大学;2008年
5 王翠玲;色胺酮、吴茱萸次碱等天然产物的生物合成及仿生合成研究[D];西北大学;2007年
6 张敏;环糊精对含能化合物的分子识别及其催化研究[D];西北大学;2007年
7 郑继明;单分子电子输运性质的第一性原理研究[D];西北大学;2008年
8 陈学元;晶场对称性方法研究稀土激光晶体的光谱和磁学性能[D];中国科学院研究生院(福建物质结构研究所);1998年
9 王莉衡;基因工程菌E.coli的热动力学研究[D];西北大学;2007年
10 徐哲峰;数论中一些和式的算术性质研究[D];西北大学;2007年
中国硕士学位论文全文数据库 前10条
1 徐刚;西北大学后勤社会化改革研究[D];西北大学;2003年
2 强媛;可逆信息隐藏及无损图像认证技术研究[D];西北大学;2009年
3 张凯;非对易轨形上的投影算子[D];西北大学;2009年
4 张陈俊;gl(1|1)超对称顶角模型的边界关联函数的研究[D];西北大学;2008年
5 张扬;利用密度矩阵元之间的关系构造两体可分准则[D];河北师范大学;2011年
6 刘小川;一类非线性抛物方程解的长时间渐进行为[D];西北大学;2007年
7 乔泽邦;四氮杂杯芳烃和四氮杂环蕃的电子结构研究及新型四氮杂环工番的合成[D];西北大学;2007年
8 张小燕;基于WebGIS的综合风险信息服务平台的研究[D];西北大学;2007年
9 赵艳;P-adic域上的微积分理论和级数理论[D];西北大学;2008年
10 李磊;基于CT断层的颅骨破洞修复技术研究[D];西北大学;2008年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026