PPV中极化子动力学方程的龙格-库塔求解.pdfVIP

第28卷 第 5期 计 算 物 理 Vo1．28．No．5 2011年 9月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Sep．，2011 文章编号 ：1001—246X(2011)05-0743-06 PPV 中极化子动力学方程的龙格-库塔求解 赵俊卿， 贾振锋 ， 张天佑， 丁 猛 ， 乔士柱 ， 陈 莹 季燕菊， 张宁玉， 付 刚 (山东建筑大学 理学 院，山东 济南 250101) 摘 要 ：构造求解一阶微分方程组初值问题的八阶龙格．库塔递推公式 ，结合描述有机分子运动 的一维紧束缚模 型，研究PPV原子链 中极化子 的形成及运动．对碳原子数 N=160的PPV原子链 ，由可控步长八阶龙格．库塔公 式求解 2N (2N+1) =102720个方程组成的方程组，用Fortran语言编程计算，得到稳定的极化子结构和运动图 像 ；在场强 E=1x10V ·em-1的电场作用下，极化子沿分子链 的运动速率约为0．2635A ·fs～．计算结果表 明，八阶龙格一库塔方程可 以有效地用于有机分子链 中载流子运动的模拟． 关键词 ：八阶龙格一库塔公式 ；PPV；极化子 ；Fortran 中图分类号 ：O411 文献标识码 ：A 0 引言 很多物理问题的描述最终归结为微分方程组的求解．由于方程组的复杂与庞大，解析求解难以实现，往 往借助于数值计算方法．随着科学的发展，计算方法显得愈加重要．求解微分方程的数值方法多数是建立在 对欧拉 (Euler)方法的改进或应用差分的基础上．龙格．库塔 (Runge．Kutta，简写为 R．K)方法对欧拉方法中的 斜率作 了适 当改进 ，大大减小了积累误差 ，是求解一阶微分方程初值问题的非常实用 的方法 ． 在聚合物材料中，由于电子 一晶格相互作用较强 ，额外的电子或空穴将会诱导晶格畸变 ，从而形成孤子、 极化子 、双极化子等 自局域元激发 ．聚合物 中的载流子通常是这些既带 电荷、又有 晶格畸变 的元激发．这 些元激发可以由求解定态薛定谔方程得到，但这些性质只是研究聚合物性质的第一步，更丰富的内容是研究 其动力学过程．动力学过程以定态元激发为初始条件 ，研究产生的各种元激发的输运和衰变过程．这对理解 聚合物 中的电导、自旋扩散、光吸收、光电导、荧光等现象都是不可缺少的 J． 聚合物材料 PPV(poly—pheny1enevinylene)是世界上第一个报道 、也是20多年来研究最多的有机电致发 光材料之一，具有 良好的化学物理特性 引．物理上，通常将 PPV材料抽象为一维分子链 ，基于一维紧束缚模 型，采用非绝热动力学方法进行研究 ¨．原则上 ，给定哈密顿量 ，任何初态 的演变都容易通过求解含时薛定 谔方程得到．由于电子未必处于晶格位形所对应的电子本征态，因此要把耦合的晶格态和电子态 同时进行求 解 ，允许电子在瞬时本征能级之间发生跃迁． 对于由Ⅳ个原子组成的一维分子链 ，将晶格振动的二阶微分方程转化为一阶微分方程后 ，通过耦合 的 晶格态和电子态确定载流子的动力学状态 ，需要联立求解 由2N(2N+1)个一阶微分方程组成 的方程组 ，计 算量非常大．为了减少累积误差 ，对计算方法提出了很高的要求．一般资料 中对于 R．K方法的介绍只给出低 于四阶的递推公式 ’ 引．为了提高大型计算的精度 ，本文依据 R．K方法的基本原理 ，构造求解一阶微分方 程组初值问题的八阶 R—K递推公式 ，采用可控步长的Fortran程序 ，计算得到 PPV中稳定 的极化子位形和运 动状态 ，取得了较好的效果． 1 八阶 R-K方程 一 般地 ，高阶微分方程总可以化为一阶微分方程组来求解 ，所以只需要研究一阶微分方程初值问题的数 收稿 日期 ：2010—11—19；修 回日期 ：2011—04—09 基金项 目：国家 自然科学基金 及住房与城乡建设部研究计划(2010K4．15)资助项 目 作者简介 ：赵俊卿 (1964一)，女，山东邹平 ，教授 ，博士，主要从事有机半导体光电材料与器件研究，E-mail：zjq@sdjzu．ed