（原子与分子物理专业论文）双原子分子光解的含时波包动力学研究.pdf

独创声明 本人声明所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研 究成果。据我所知，除了文中特别加以标注和致谢的地方外，论文中不包含其他 人已经发表或撰写过的研究成果，也不包含为获得( 注：如没 、 有其他需要特别声明的，本栏可空) 或其他教育机构的学位或证书使用过的材料。 与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表 示谢意。 学位论文作者签名：气基，、匀易弓 导师签字： 学位论文版权使用授权书 本学位论文作者完全了解堂撞有关保留、使用学位论文的规定，有权保留 并向国家有关部门或机构送交论文的复印件和磁盘，允许论文被查阅和借阅。本 人授权堂撞可以将学位论文的全部或部分内容编入有关数据库进行检索，可以 采用影印、缩印或扫描等复制手段保存、汇编学位论文。( 保密的学位论文在解密 后适用本授权书) 一警锶屹氓聊警 签字日期：2 0 0 8 年4 月8 日签字日期：2 0 0 8 年4 月8 日 山东师范大学硕士学位论文 双原子分子光解的含时波包动力学研究 摘要 分子反应动力学是从原子、分子层次出发研究化学反应微观动态和机理的科 学，它的研究不仅可以阐明分子反应过程中各种瞬态物质的结构、性质和作用， 还能通过对态一态反应动力学，以及对分子相干态之间作用的深入研究来阐明化 学反应的内在规律。光解反应动力学是光化学的一部分，也是分子反应动力学的 一个重要分支。作为分子反应动力学重要研究工具的含时波包动力学理论在分子 物理和场与分子相互作用等研究领域也有着广泛的应用。含时波包法有许多优 点，除了数值计算上的高效率外，该方法还为动力学研究提供了物理意义明确而 直观的图像，它既具有经典的直观，又不乏量子力学的准确。 本论文主要是基于b o m o p p e l l l l e i m e r 近似，应用含时波包动力学方法研究 了双原子分子n a h 光解的吸收截面，讨论了时间传播方法、势能面、初始波 包对光解截面的影响，主要的研究内容和结果如下： 首先在绝热近似下，利用傅立叶格点一哈密顿方法( f g h ) 求解含时薛定谔 方程、劈裂算符方法作为波包传播方案计算了n a h 分子从基态x1 + ( 1 ，”= o ) 跃 迁到激发态b1 n 上的光解截面。为了研究不同的波包传播方案对n a h 分子光解 截面的影响，我们将计算结果与使用c h e b y s h e v 多项式展开方法得到的结果进行 比较，结果表明这两种时间传播方案对总的光吸收截面的计算是等效的。其次， 依据离解阈能的实验值，我们将势能曲线进行了重置，然后将计算结果与 b h a t t a c h 撕e e 等人的结果进行比较，得到了影响光解截面的关键因素。最后选取 了四个不同的振动态( v ”= o 、1 、2 、3 ) 作为初始波包，研究了不同的初始波包对光 解截面的影响。 论文共分为五章。第一章为综述，主要介绍了分子反应动力学的基本理论 及波包动力学方法的发展。第二章对光解动力学进行了简介。第三章介绍了研究 光解动力学所用到的理论和方法。第四章，主要介绍了研究过程和计算结果，并 山东师范大学硕士学位论文 讨论了时间传播方法、势能面、初始波包对光解截面的影响。第五章是本论文的 结束语，简要地对本论文进行了总结。 关键词：含时波包动力学；傅立叶格点一哈密顿方法( f g h ) ；光解截面； n a h 分类号：0 5 6 1 4 4 山东师范大学硕士学位论文 t i m e d e p e n d e n tw a v e p a c k e ta p p r o a c ht 0t h e p h o t o n d i s s o c i a t i o n0 fd i a t o m i cm o l e c u l e s a b s t r a c t m 0 1 e c u l a rr e a c t i o nd y n 锄i c ss t u d i e sm em i c r o c o s 盯l i cd y n 疵c sa n dp 打n c i p l eo f 也ec h e 删c a lr e a c t i o nb a s e do na t o m i ca i l dm o l e c u l a rl e v e l nc 狮e h l c i d a t en o t0 1 1 1 y s 咖c t l l l e s ，c h a r a c t e r sa i l df h l l c t i o n so fa 1 1k i n d so f 仃a l l s i e n ts p e c i e sd 面n gr e a c t i o n ，b u t a l s oi n 删l 】s i cl a wo fc h e r n i c a lr e a c t i o nf b mm er e s e a r c ho ns t a t e t o s t a t er e a c t i o n d y n 锄i c sa n dm ei i l t e r a c t i o no fc o h e r e ms t a t e s p h o t o d i s s o c i a t i o ni s a tm eh e a no f p h o t o c h e m i s 包阱a n di t i sa l s oa ni 瑚p o r t a n te n l b r a i l c l m l e n ti nm o l e c u l a rr e a c t i o n d y n a l l l i c s t h et i m e d e n p e l l d e tw a v e - p a c k e td y n a m i c s 勰廿l et o o lo fm o l e c u l a rr e a c t i o n d y n a i i l i c sr e s e 卸c hh 嬲m 锄ya p p l i c a t i o n si nm o l e c u l a rp h y s i c sa n d 缸l e 丘e l d - m 扰e r n e r a c t i o ns y s t e m s t h et i i i l e - d 印e n d e mw a v e - p a c k e tm e m o di sv e r ya l t r a c t i v e b e s i d e s m en 吼e r i c a le 伍c i e l l c y ，t i i l l e d 印e n d e i l tw a v e - p a c k e tm e t l l o di sc o n c 印t u a l l ys i i n p l e i t p r 0 v i d e sn o to n l yac l a s s i c a l l i k ei n t e 印r e t a t i o n ，b u ta l s om eq u a n t u mp r e c i s i o n t h i st l l e s i ss t l l d i e sm ea b s o r p t i o nc r o s ss e c t i o no fd i a t o m i cm o l e c u l e n a h m o l e c u l e 嘶n gi t sp h o t o d i s s o c i a t i o np r o c e s so nt h eb a s i so f l eb o r i i o p p e n h e i m e r a p p r o x i m a t i o nt h e o 巧t h ei n n u e n c e o ft h ed i f f e r e n tp r o p a g a t i o nm e m o d s ， m e d i f 董醯e n tp o t e n t i a le n e r g ec u r 、，e s ，a sw e l la sm ed i f f e r e n ti n i t i a lw a v e - p a c k e t st om e a b s o 印t i o nc r o s ss e c t i o ni sd i s c u s s e d t h em a i nc o n t c i l t sa n dr e s u l t sa r ep r e s e n t e da s f o l l o w s f i r s t ，m i st h e s i sc a l c u l a t e st h ea b s o 印t i o nc r o s ss e c t i o no fn a hm o l e c u l e 丘o m g r o u n ds t a t ex1 + ( 1 ，”= o ) t 0b1 兀u s i i l gf o u r i e r 鲥dh 锄i l t o i l i a nm e t h o dt 0 s o l v et h et 曲e d 印e n d e n ts c h r 6 d i n g e re q u a t i b na n du s i n gs p l i to p e r a t o rn l e t l l o da sa w a v e p a c k e tp r o p a g a t i o ns c h 锄eo nt h eb a s i so fm ea d i a b a t i ca p p r o x i m a t i o nt h e o 珥 5 山东师范大学硕士学位论文 i i lo r d e rt o s t u d yt h ee a e c t s o fd i f f 色r e n tw a v e p a c k e t p r o p a g a t i o ns c h 锄eo n p h o t o d i s s o c i a t i o nc r o s ss e c t i o no fn a h ，t h er e s u l t so b t a i n e db yp l i t o p e r a t o rm e m o d w a sc o m p a r e dw i mt h o s eo b t a i n e db yt h ec h e b y s h e vp r o p a g a t i o ns c h e m e t h er e s u l t s i n d i c a t et h a tt h et w os c h e m e sa r ee q u i v a l e n ti nc a l c u l a t em et o t a la b s o 印t i o nc r o s s s e c t i o n t h e n ，a c c o r d i n gt ot h ee x p e r i m e n t a lr e s u l to ft h et h r e s h o l de n e r g e ，t h e p o t e n t i a le n e r g ec u r v e sw e r er e a r r a l l g e d ，t h er e s u l t sw eo b t a i n e dw e r ec o m p a i r e dw i t h m er u s u l t s 百v e db yb h a t t a c h 面e e 刃口z ，t h ek e yi n n u e n c e 蠡犯t o ro f p h o t o d i s s o c i a t i o n c r o s ss e c t i o nw a so b t a i n e d l a s t ，f o u rd i 妇陀r e n tv i b r a t i o n a ls t a t e s ( 1 ，”= 0 ，l ，2 ，3 ) w e r e t a k e na si 芏l i t i a lw a v e p a c k e t st oi i l v e s t a g a t em ei i l n u e n c et ot h ep h o t o d i s s o c i a t i o n c r o s ss e c t i o n 啦ss u b r l l i t t e dw o r ki sd i v i d e di n t 0 丘v ec h a p t e r s t h e 丘r s t c h a p t e ri s m e i i l _ 缸- 0 m l c t i o n ，w 1 1 i c hb r i e n yd i s c u s s e st l l eb a s i ct h e o r yo fm 0 1 e c u l a rr e a c t i o nd y n 锄i c s a 1 s op r e s e n t e di st h ed e v e l o p m e n to ft h et i m e - d 印e n d e n tw a v e _ p a c k e tm e t h o d c h a p t e r觚oi n 咖d u c t sm ed e v e l o p m e n ta n ds i g n i f i c a n c eo fp h o t o d i s s o c i a t i o n d y l l 锄i c s hc h 印t e rt l 鹏e ，m em e o r e t i c a la n dc o 砌p u t a t i o n a lm e m o d su s e di n p h o t o d i s s o c i a t i o nd y n 锄i c sa r ep r e s e n t e d t h er e s e a r c hp r o c e s sa r l dt l l er e s u l t sa r e 百v e i lp a n i c u l 砌ya n dt h ei n n u e n c eo fd i f f e r 饥tp r o p a g a t i o nm e t l l o d s 、 m ed i 髓r e n t p o t e n t i a le n e r g ec 眦弋，e s 、m ed i f f e r e n ti n i t i a lw a v e - p a c k e t so nt h ep h o t o d i s s o c i a t i o n c r o s ss e c t i o na r ed i s s c u s s e di nt h ef b l l r t hc h 印t e r t h es 嘲a r y 瓶dm ep r o s p e c ta r e r e p r e s e n t e di nm e1 a s tc h a p t e r k e yw o r d s ：t i i l l e d e p 髓d e n tw a v e p a c k e tm e t h o d ，f o 证e r 鲥dh 锄i l t 砌a nm e t h o d ， p h o t o d i s s o c i a t i o nc r o s ss e c t i o n ，n a h c i a s s i n c a t i o n ：0 5 6 1 。4 6 山东师范大学硕士学位论文 第一章综述帚一早殊尬 1 1 分子反应动力学的基本理论方法 分子反应动力学是从原子、分子层次出发研究化学反应微观动态和机理的科 学。人们已经发展了多种实验和理论方法来研究分子反应动力学过程。在飞秒 ( 1 0 。1 5 s ) 、皮秒( 1 0 1 2 s ) 、阿秒( 1 0 舶s ) 这些超快激光脉冲激光技术出现以前，人们 主要致力于检测产物的各种特性：产物寿命、几何构型、平动、振动、转动各自 由度上的能量分布以及产物的角分布等，以此来间接地了解化学反应的发生过 程。在飞秒( 1 0 。5 s ) 、皮秒( 1 0 。2 s ) 、阿秒( 1 0 。8 s ) 出现以后，人们不仅间接地研究 化学反应过程，而且还可以实时跟踪探测化学反应过程。在理论研究方面，主要 提出四种方法：( 1 ) 经典轨道方法，( 2 ) 半经典方法，( 3 ) 全量子力学方法，( 4 ) 含时量子波包动力学方法。经典轨道方法是研究化学反应常用的方法之一【1 】，该 方法已经成功地应用到大量的化学反应体系中【2 ，3 】。该方法的缺点是：( 1 ) 在进 行全末态分布分析时，需要大量的相空间轨道，计算量将大幅度增加；( 2 ) 该方 法无法准确描述量子效应，如隧道效应、零点能分子的量子运动等。目前使用的 半经典方法主要为高斯波包近似方法【4 ，5 】，该方法包含了量子效应，但要求对不 同的体系做不同的近似，其有效性有待于进一步验证。全量子力学方法主要基于 散射理论中的密耦近似( c l o s e c o u p l i n ga p p r o x i i i l a t i o n ) 【6 】该方法需求解大量的 耦合通道，对于具有丰富振一转耦合能级的大分子体系，精确的数值求解几乎不 可能实现【7 1 。最近几年发展的含时量子波包动力学方法【8 母】主要为求解含时 s h r 6 d i n g e r 方程，该方法对描述化学反应动力学是很有效的，只需知道系统初始 时刻的状态和势能面，系统随时间演化的过程就可以完全确定，采用该方法很容 易研究连续态演化与化学反应问题。 含时量子波包法有许多优点，除了数值计算上的高效率外，该方法还为动力 学提供了物理意义明确而直观的图像，它既具有经典的直观，又不乏量子力学的 准确。另外，含时量子波包法尤其适用于研究体系随时间演化问题。目前，波包 法在化学反应、气一固表面相互作用、光解离、光电离、低温物理、半导体物理 7 山东师范大学硕士学位论文 等领域均有着重要的应用。 1 2 波包动力学的发展及应用 波包动力学是物理化学的重要分支。其理论框架在分子物理及场与物质相互 作用方面都有很多应用。自1 9 2 0 年代至今，波包动力学已经经历了若干个发展 阶段。1 9 2 6 年，为了架起经典力学描述和量子力学描述之间的桥梁，s c h 埔d i n g e r 在量子力学中引入了波包的概念来描述微观粒子运动，提出了含时s c h r 6 d i n g e r 方程， 汤虽l 甲( f ) ) = 日l y ( f ) ) ( 1 1 ) 并指出波包的运动遵循经典的轨迹【1 0 】，标志着波包动力学理论的“开创阶段 。 紧接着，1 9 2 7 年，e 1 1 r e n f e s t 发展了波包的概念，指出在经典极限下，量子力学 的期望值遵循经典力学的规律，这就是著名的i r e i 血s t 定理【1 1 】。此外，w e n t z e l f l 2 】、 融锄e r s 【1 3 1 和b r i l l o u i i l 【1 4 】提出的半经验方法补充了波包的位相等其他一些重要性 质。由于受当时的实验条件的限制，要单独制备一个波包看起来几乎是不可能的， 所以在随后1 9 3 0 年代至1 9 6 0 年代的近3 0 年里，波包的概念一直没有得到实际 的应用，大部分的散射理论都侧重于不含时的量子力学方法。这一时期被称为波 包动力学的“睡眠阶段 。 直到1 9 6 0 年代中期，物理和化学中关于激光和原子、分子相互作用研究的 突破，才使波包及其动力学重新受到重视。计算方法以及计算机技术的飞速发展 为求解含时s c h r 6 d i n g e r 方程提供了可能性。人们发展了各种近似方法来处理含 时s c d i n g e r 方程的空间部分和时间部分，造成了1 9 7 0 年代的“复兴阶段”， 此后波包动力学理论有四个方向的发展：一是，j o n i l e r 等预测了一系列新的时间 相关现象【15 1 ，例如d 锄p e dq u 锄t 啪b e a t s 。二是，w y a t t 等 1 6 1 和后来的h e l l e r 等【1 7 】 发展了在分子结构空间计算波函数l 甲( f ) ) 的方法。尽管存在着计算方法的限制， 但这些早期( g i 甲( f ) ) 计算还是给出了坐标空间中分子波包动力学的模拟方法。三 是，g o n d o n 【1 8 】和h e l l e r 【1 9 】给出了如何将分子波包动力学、时间相关函数和分子 山东师范大学硕士学位论文 c w 光谱联系起来的处理方法。以上三方面都是以早期的h e i s e n b e r g 【2 0 。、 s c h r 6 d i n 9 0 2 2 1 以及l a ) ( 【2 3 】、k u b o 和l n l 锄o t o 【2 5 】的工作为基础发展起来的。最 后一方面的发展是b l o e i n b e r g e i l 和l a r s e n l 2 6 l 开始对纳秒和c w 极限的激光场进行 精确计算。 与此同时，飞秒激光技术的出现及其迅速发展，特别是在物理和化学中的应 用客观上促进了含时波包动力学的发展。7 0 年代用激光锁模( m o d el o c k i n g ) 技 术和染料激光器相结合成功地制备出了皮秒( p s ) 级的脉冲并在化学和生物学领 域中得到了广泛的应用【2 7 1 ，如溶液中的碰撞诱导解离反应，电子转移，分子能 量的内转换以及振动弛豫等。在1 9 8 7 一1 9 8 8 年波包动力学实验出现了“突破阶段”。 众所周知，分子化学键的有效作用范围一般在1 姗以内，而室温下分子的运动 速度大约是1 l c r l l s 。所以分子碰撞造成旧化学键断裂到新化学键生成的过程一般 在1 0 0 2 秒( 即1 p s ) 左右完成。也就是说，过渡态的寿命大约在1 p s 数量级。以 往的实验手段，包括交叉分子束实验也只能观测反应前后的现象。要想更直接地 观察到过渡态，或者说要实现对分子体系从反应物经过渡态到产物连续演化过程 的实时测量和控制，就需要一种具有更高时间分辨的实验技术。s h a l l l 【首先在物 理实验室中制备了短达6 飞秒( f s ) 的脉冲激光【2 8 】，随后，a 1 伽1 0 和s h a p 沛又 制造出了波长连续可调的超短脉冲激光【2 9 1 。1 9 8 7 年，z e w a i l 实现了第一个飞秒 化学实验，他们使用飞秒级的泵浦和探测脉冲激光检测了i c n 和n a 】分子的过 渡态光谱【3 0 1 。这一工作开创了在飞秒尺度上对化学反应进行实时监测，被认为 是奠定了飞秒化学的实验基础。1 9 9 1 年，由s i b b e t t 完成的含钛蓝宝石激光器为 广泛使用超短脉冲激光探测基元反应铺平了道路【3 。由于飞秒脉冲的宽度只有 几个到几十个飞秒，比分子的振动周期还要短，所以它揭开了超快反应过程的神 秘面纱：分子碰撞的详细动态图象已经生动的展现在面前，相应的理论上也有更 进一步的发展。1 9 8 2 1 9 8 3 年，f e i t 与f 1 e c k 【3 2 】和k b s l o 一3 3 1 独立研究出并推广了 f o 嘶e r 网格点方法，从而使波包动力学很快流行起来。此外，“g h t 等【3 4 】发展了 不连续变量表示方法，z h u 等【3 5 】引入了辛算法。这些方法的基本思想都是将时空 划分成若干网格，并用一套函数集合近似代替解析函数。最近十几年，以飞秒时 间尺度上超快化学反应过程为研究对象的飞秒化学已经成为一门十分活跃的新 兴分支学科。 9 山东师范大学硕士学位论文 飞秒激光在最近十几年时间里的快速发展，为人类充分探索微观超快现象 及研究强场物理提供了前所未见的实用手段和发展机遇，但是对于发生在阿秒量 、 级( 1 0 0 8 s ) 的物理现象，如原子内电子的跃迁和弛豫过程，飞秒脉冲就显得无 能为力了。为了研究阿秒时间量级的超快过程中所发生的瞬态现象，就需要探索 产生阿秒脉冲的各种可能的技术和方法。目前人们在产生阿秒脉冲的原理方法上 作了大量的前沿探索工作【3 6 4 2 1 ，其中最有前途的方案是超短脉冲的谐波合成技 术。具体的可行技术有两种：一种是利用超强超短脉冲与惰性气体的非线性效应 产生高次谐波【4 1 】，或者受激拉曼散射从而得到阿秒脉冲的方法：另一种是相位锁 定光学参量或同步飞秒激光产生的可见光亚谐波合成技术【4 3 粕】。最近，前一种方 法，即高次谐波产生阿秒脉冲技术有了突破性的进展。一个由奥地利维也纳技术 大学、加拿大国家研究中心和德国比利斐尔德大学的研究人员组成的国际课题 组，报道了他们产生的单个脉冲6 5 0 a s 的x 射线脉冲【4 5 1 。单个阿秒脉冲的出现在 超快光谱学方面掀开了新的篇章，科学家也将从研究原子分子动力学过程的飞 秒化学时代进入研究原子内电子动力学的阿秒物理时代。随着阿秒脉冲产生技术 及抽运探测技术的进一步发展，能够用以揭示电子动力学过程( 如内壳层动力学 过程) 并控制这些过程的阿秒激光物理无疑将是自然科学的一个重要前沿课题。 波包动力学也被广泛应用于各种化学超快过程和飞秒、皮秒、阿秒脉冲激发的时 间解析光谱的理论研究。 研究波包动力学的一个主要原因是人们期望能实现对化学反应的控制。因为 许多分子内部运动过程和化学反应大都发生在飞秒( 1 0 一5 s ) 、皮秒( 1 0 。1 2 s ) 或者阿 秒( 1 0 。8 s ) 量级，所以超快激光脉冲自问世以来就一直被用来作为实时监测和控制 分子反应动力学过程的工具【4 6 舶】。这种控制是一种“时域”控制，通过调节激光 的强度、时间形状、频率分布以及脉冲之间的延迟时间、相位差等来实现对分子 内部运动的有效控制。分子在超短激光脉冲的作用下，能够表现出许多有趣的现 象：例如超阈值电离 4 7 】、键的硬化或软化4 8 1 、a c s t a r k 效应【4 9 1 、拉比振荡【5 0 5 1 l 等。目前用激光控制化学反应的研究目前在国际上已经掀起了热潮，选键光解离 正是激光控制化学反应的一种手段。通过激发分子到不同的电子激发态或者到同 一电子激发态势能面的不同区域可以达到选择性断键、控制不同反应通道分支比 的目的。利用不同波长的激光可以选择性地断键，可以影响不同通道不同产物的 l o 山东师范大学硕士学位论文 量子产率。 山东师范大学硕士学位论文 2 1 引言 第二章光解动力学简介 束缚态分子吸收一个或多个光子发生解离的过程叫光解。光束的电磁能转变 为分子内部的能量，如果转移的能量超过了最弱键的束缚能，分子将不可避免地 发生解离： ( a b ) + 。枷壳卫b ( a b ) 玛a ( 口) + b ( ) ( 2 1 ) 这里a b 代表母分子，( a b ) 代表激发态络合物，壳国是光子能量，嘞咖雄是吸收 光子的数目，口和表示碎片a 和b 特定的内部量子态。式( 2 1 ) 中第一步表 示母分子吸收光子过程；第二步是激发态络合物解离过程【1 】o 光解反应动力学是光化学的一部分也是分子反应动力学的一个重要分支，它 主要研究分子吸收光子后发生分解反应产生光解碎片的过程以及光解碎片的内 能分布和空间分布。早在1 9 5 0 年n o m s h 和p o r t e r 等人就首次探测到了光解碎片 的内能分布并开创了这一崭新学科。又过了大约1 0 年w i l s o n 【2 】等人第一次测定 了光解产物的平动能分布和角分布。光解反应动力学在过去的几十年里取得了巨 大的发展，它深入到物理化学和化学物理的各个领域，包括环境化学、团簇化学、 气一固表面现象以及大气化学等。到了二十世纪九十年代初，高分辨的h 原子 里德堡时间飞渡光谱仪( h i 曲一r e s 0 1 u t i o nha t o m 玛，d b e r gt i m e o ff 1 i 蛳 s p e c t r o m e t e r ) 发展起来【3 1 ，它可以精确测量光解碎片的平动能分布并能够推测出 光解碎片的内能分布( 包括碎片的振动能分布及转动能分布) ，这使光解反应动力 学又有了新的发展。光解反应动力学实验主要探测两个关系：一是能量分配关系： 即光子能量在光解过程中和光解反应完成后的分配情况；二是光解过程中的矢量 关系：即光解碎片相对于激光电矢量的角分布以及光解碎片的转动角动量等。测 量光解碎片的平动能可以得到重要的动力学信息。如：对a 专b + c 类型的解离， 由动量守恒定理可知，两个反弹碎片的速度是互相关联的。因而，当测定其中一 碎片的速度并将其转变到质心系时，另一部分的速度分布及相应的平动能分布就 山东师范大学硕士学位论文 可完全被确定。如果解离产生了一个原子( 仅有一个电子态) 及一个自由基，则自 由基的内能也可以完全被确定。 在光解理论方面，h e l l 0 4 l 发现在计算总的光吸收截面( 也就是吸收谱) 的 时候，含时方法和非含时方法是等效的。接着k u l a n d e r 和h e l l e r 【5 】利用含时方法 计算得到了微分散射截面。k o s l o 缚6 】等人把快速傅立叶变换( f f t ) 引入到含时量 、 子力学的计算中，使哈密顿算符中的动能部分对波函数作用的计算变得简单而且 节省机时。特别是l i 曲t 和他的合作者【7 ，8 】把离散变量表象( d v r ) 引入到分子反应 动力学研究领域，使含时波包方法得到迅速发展。人们利用含时波包方法成功地 研究了大量分子的光解动力学，特别是小分子n a i 【9 】，h 2 0 【l o 1 1 1 ，0 3 【1 2 ，13 1 ，h 2 s 【1 4 ，1 5 】 以及c h 3 i 【1 6 ，1 7 】。在理论上对光解动力学的反应机理做了深入地研究，模拟出和 实验比较类似的结果，进而给实验工作者提出了很多有益的意见和建议。 光解动力学研究的理论意义在于：作为分子碰撞反应的半过程，它能够控制 反应物的初态和探测产物的终态，对于态一态分辨的分子反应动力学研究具有重 要的指导意义。光解反应是单分子反应的一种，它主要发生在分子的电子激发态。 分子的电子激发态与电子基态在结构、反应活性等方面存在较大的差异，因此研 究分子的光解反应动力学有助于我们了解分子激发态的性质，获得有关复杂分子 势能面的信息，如分子势能面间非绝热交叉的情况等等。 分子光解反应动力学的研究涉及到基元化学反应的核心问题：化学键的断裂 和生成。对单分子吸收一个光子的光解过程的研究可以极大地帮助我们理解分子 断键过程的机理。通过激发分子到不同的电子激发态或者到同一电子激发态势能 面的不同区域可以达到选择性断键、控制不同反应通道分支比的目的。利用不同 波长的激光可以选择性地断键，可以影响不同通道不同产物的量子产率。对激光 控制化学反应的研究目前在国际上已经掀起了热潮，选键光解正是激光控制化学 反应的一种手段。 光解动力学的实际意义在于光解是大气化学中许多重要链反应的发动机，而 太阳提供了足够的紫外辐射来保持这台发动机运转。例如，臭氧层【5 6 】对地球上 的生命来说至关重要，它控制着全球高层大气中臭氧的含量，对于吸收对生物有 害的太阳紫外辐射、保护人类健康和生存环境起着非常重要的作用，是人类必不 可少的保护伞。 1 4 山东师范大学硕士学位论文 此外，光解是许多化学激光产生的起点。光解反应： c h 3 i + 壳缈专c h 3 + i 。 ( 2 3 ) 产生了碘的两个电子态i + ( 2 p l 2 ) 和i ( 2 p 3 2 ) 。利用两个态之间的反转布居，在随后发 生的z p l 2 专2 p 3 2 跃迁过程中可以建立一束激光。光解产物的终态分布往往是反转 的，而这恰恰是激光产生的先决条件。为了系统地发现和控制新的激光光源，毫 无疑问我们应该着重研究光解反应及其相关过程。 2 2 光解离的类型 类。 根据激发态络合物解离过程的不同光解反应可分为直接解离和预解离两大 e 图2 1 分子直接解离示意图 a 直接解离( 如图2 1 ) ：直接解离是指母分子吸收光子后跃迁到排斥型电 子激发态势能面，分子在半个振动周期内即发生解离。其激发态分子寿 命，c 很短暂，大约在1 0 以3 s 秒左右。最显著的特点是母分子的吸收光谱是 连续谱带。 山东师范大学硕：l ：学位论文 1 6 b 预解离( 如图2 2 2 4 ) ：预解离是指处于激发态的母分子不直接发生解离， 而是经过其它一些中间步骤之后发生解离。预解离又分为： 类型i ：振动预解离( 如图2 2 ) ，分子受激到一个电子束缚态的高振动激发 态，经历了一段分子内振动能重新分配的过程后发生分解。解离的 速率般在10 6 l0 1 0 s 一，比电子态直接解离慢得多。 e 图2 2 分子预解离类型i ：振动预解离 类型i i ：电子预解离( 如图2 3 ) ，分子受激到某一上能态，这个能态是一束 缚态，分子不能在这个能态上发生解离。但这个能态与另一排斥型 电子态的势能面交叉，两个电子态之间的耦合使得部分受激的分子 跃迁到排斥型势能面上，则这部分分子可以很快地发生分解。因为 电子排斥型势能面的解离速度极快，所以解离的快慢取决于两个能 态之间耦合的强弱，耦合弱的解离速度慢；耦合强的解离速度快， 但一般小于1 0 s 。 山东师范大学硕士学位论文 e 图2 3 分子预解离类型i i ：电子预解离 光解反应可以根据不同的分类标准分为不同的反应类型，但由于本论文所选 分子势能面的限制，本论文主要涉及到以上三种解离类型，它们是光解离动力学 过程很常见的解离类型。 1 7 山东师范大学硕士学位论文 第三章含时波包动力学基本理论 分子体系的势能面是在b o m o p p e n h e i m e r 近似下获得的。在己知某分子相 应势能面的情况下，使用含时波包方法处理分子动力学问题及其相关的数值方 法，已经比较成熟。处理具体的分子动力学问题时，比如说，计算分子的光解离 问题，计算分子的本征态作为含时问题的初始值，以及波包的传播等等，目前都 有高效精确的数值方法可以参考使用。 含时波包法的核心在于求解含时薛定谔方程，主要可以分为如下几个步骤： 1 计算相关势能面或者利用己知的势能面； 2 利用一定的方法求出体系基态的初始波函数，进而得到体系激发态初始 波包； 3 利用相应的h a m i l t o l l i a n 模型在网格点上计算出传播所用的h a m i l t o n i a l l 矩阵，或使用f f t 技术计算动能算符对波函数的作用； 4 根据模型选择合适的传播波包的方法，传播波包； 5 利用传播后的波包，计算所感兴趣的体系的物理量，完成分子动力学过 程的模拟。 在本章中，就本论文所用到的处理分子动力学问题的数值方法，做一个简单 讨论。 3 1 波恩一奥本海默近似( b o r n o p p e n h e i m e ra p p r o x i m a t i o n ) 任何一个分子体系都是由一定数目的电子和原子核构成的。在不考虑相对论 效应的前提下，体系的运动由薛定谔方程来描述。电子的质量比原子核轻至少3 个数量级以上，原子核运动的速度比电子运动的速度慢很多，因此我们可以采用 波恩一奥本海默近似( b o m o p p e l l h e 曲e ra p p r o x i m a t i o n ) ，即在电子组态中任何 改变都可以认为是瞬时发生的，核的状态发生一些微小的改变，电子都会立刻进 行调整并与之相适应，同时核间的相对运动可以看作是电子运动平均作用的结 1 9 山东师范大学硕士学位论文 设有一个分子体系，采用原子单位制，分子满足定态薛定谔方程 b 去v ；一军扣c 叫y 却 ， 其中y ( r ，) 表示体系的相互作用能， 哪。荟等+ 荟去一丢吾 2 ， 上式中p 和g 用来标记原子核，f 与七用来标记电子。尺朋表示核间距离，表 示电子间距离，名表示核与电子间距离，乙为原子核的质子数。设 少( r ，) = 痧( r ，) z ( r )( 3 3 ) 其中俾，) 和z ( 尺) 分别表示相应电子和核的波函数，r 与，分别对应于核与电 子的位置。将( 3 3 ) 7 式代入( 3 1 ) 式中，并利用 v 2 ( 加) = 2 v 妒缈+ 妒2 缈+ 妒2 矽 ( 3 4 ) 得 一莓去妒；z 一莓丢矽；矽一莓 p 妒p z 一军三押勿+ 矿( r 川觑= e 奶 c 3 s ， 对于一般的分子，有 朋p 1 0 3 1 0 5 ( 3 。6 a ) v ，妒，z 椤；矽 ( 3 6 b ) v ；= v 勿， ( 3 6 c ) 因此( 3 5 ) 式中的第二、第三两项可以忽略。两端同时除以沙，离散变量后得 到 一军参羚幽力= w ( 3 7 ) 一莓# = 励 8 ， 2 0 山东师范大学硕士学位论文 这样，求体系总能量e 的过程可以分为两步：首先，固定核间距r 计算u ；然后 利用( 3 8 ) 式求体系总能量e 。以上就是b o m o p p e l l l l e i m e r 近似的基本思想。 关于核运动与电子运动离散的这一点，已经通过相关的实验得到了证实。但需要 指出，对于较轻的分子体系，这种近似并不十分理想。 3 2 含时薛定谔方程 许多物理和化学现象有着共同的动态学演化模式：体系起源于初始波包，在 外场及自身哈密顿算符的作用下，随时间演化到渐近的终态波包。这些现象包括 化学反应、表面散射、光解离、光电离等等。因此求解含时薛定谔方程具有普遍 的意义。 我们知道，一个分子在外场作用下的哈密顿为： 日= 疋( 尹) + 巧( r ) + 吃( 尹) + 巧州( r ) + w ( 尹，r ) 一d e ( 3 9 ) 上式中夏，尸分别表示核坐标和电子坐标，下标e 和n 表示电子和核，r 是动能 算符，y 表示库仑相互作用。激光场雷= e ( 豆，弦通过偶极算符西与系统耦合。 对于固定的核间距j ，电子基集 l 纯) ) 可以由本征值方程得到： 日d l 纸) = u 。( 天) 1 ) ( 3 1 0 ) 其中 日d = 疋( 尹) + 吃( 尹) + ( r ) + ( 尹，r ) ( 3 1 1 ) 分子的含时波函数可以用这些态来展开： l y ( f ) ) = 眠( 晨，f ) 1 ) ( 3 1 2 ) 在b o m o p p e n h e i m e r 近似下， 瓦( ( 尼f ) i ) ) 巧虬( 天，r ) 1 ) ( 3 1 3 ) 即杀l 纪) 。用( f 左乘s c h r 6 d i n g e r 方程腑昙i 甲( 砩= 疗i 甲( ，) ) ，有 ( i 访昙i 甲( f ) ) = ( i 疗i 甲( f ) ) ( 3 1 4 或者 2 1 山东师范大学硕士学位论文 腕昙莓( i ( 天，r ) i 纯) = ( l 巧+ 也一西雷l 甲( ，) ) ( 3 1 5 ) 再利用正交归一条件( 1 ) = 瓯，。，我们得到了核运动的含时薛定谔方程： 疏导。( 辰，f ) = 兀+ 【，。( j i i ) y 。( j i i ，f ) 一( 缈。l 西吾i 纯) e ( 天，f ) 。( 天，f ) ( 3 1 6 ) u f ” 在无外场的情况下，方程( 3 1 6 ) 描述在势能面u 。( 友) ( 它对应着电子态) 上 的波包成分的时间演化。在b o m o p p e l l l l e i m e r 近似下，这种演化过程完全不依 赖于其它态和它们的可能占有情况。在有外场时，。= 一( 吼l 西引) 耦合着 刀和m 两个势能面。方程( 3 1 6 ) 忽略了相对论效应以及势能面之间的非绝热耦 合即势能面之间的交叉，需要时可以分别在哈密顿的对角元素和非对角元素中加 入。匕述方程可以用一种更简洁明了的形式来表达，即矩阵表示方法： a z 8 t l 2 ： 少。 日1 1 日2 1 日。1 且。、r 。 日2 。0y 2 1 l l ! 日。渺。 ( 3 1 7 ) 其中= 帆( 尼f ) 表示在不同电子态f = 1 ，2 ，z 的原子核波函数。哈密顿算符的 对角元素可以表示为日行= 巧+ u 。( 天) ，非对角元素描述相映两个态之间 的耦合( 包括内部耦合和外部耦合) ，外部耦合指分子与外场之间的耦合而内部 耦合指分子势能面之间的非绝热交叉。 如果只考虑二能级系统及一维问题( 双原子分子( 如h c l ) 的情形) ，可以 取f = 2 以及j = r 。下面主要讨论这种情况。 在半经典理论中，没有进行电磁场的量子化，激光场可以近似表示为： e ( r ，f ) 兰乓( 灭，f ) c o s ( 研) ，这里我们忽略了场的空间依赖性，因为激光场的波长 通常远超过分子的大小。上式中的是激光脉冲的中心频率，e ( r ，f ) 是脉冲的 形状因子，它可以取为高斯函数等形式。 在b o m o p p e n h e 妇e r 近似下，对于被一束外加激光场耦合着的一维两态问 题，有 坦 挖 砣巩也；巩 山东师范大学硕士学位论文 腩旦f 中t ( 尺，1 ：f 丁( 尺) + u 俾) 形 1 f - 职，1 ( 3 1 8 ) l ，z li = ll l j 1 6 ， 研l 2 ( 尺，f ) ll 形 r ( r ) + u 2 ( 尺) l 2 ( 尺，f ) 这里形= “：( r ) e ( r ，f ) 是基态和激发态在激光场下的相互耦合作用，h ：( r ) 是 跃迁偶极矩。