李桃桃 苑航 王兴 张震 郭大龙 朱小龙 闫顺成 赵冬梅 张少锋 许慎跃† 马新文

1)(西安交通大学物理学院,物质非平衡合成与调控教育部重点实验室,西安 710049)

2)(中国科学院近代物理研究所,兰州 730000)

3)(中国科学院大学,北京 100049)

利用反应显微成像谱仪开展了50-keV/u Ne8+ 离子与C3H4 分子碰撞实验,研究了丙二烯(CH2CCH2)和丙炔(CH3CCH)两种同分异构分子形成二价离子并解离产生H+++H 的过程.实验获得了H+ 和的动量,进而通过动量守恒得到第3 个碎片的动量.通过分析3 个碎片的动能及解离的动能释放鉴别出未被探测的碎片为中性H 原子的事件.借助Dalitz 图、Newton 图以及碎片产物的角分布等分析了该通道的动力学机制.结果表明,次序解离是该解离通道的主要机制,在碎裂过程中二价母体离子先解离成H+和,中间体的离子再进一步解离成和H 原子.

1 引言

多原子分子解离碎裂动力学是原子分子物理研究中的热点问题,同时在等离子体物理[1]、行星大气化学[2]、生物组织辐射损伤[3]等领域具有潜在应用价值.尤其是分子的三体碎裂过程中,涉及至少两个化学键的同时或次序断裂,近年来引起了极大关注[4−15].

目前,大多数关于三体解离过程的研究是针对二氧化碳等三原子分子[4−8]或小团簇[9−11]进行.在荷电粒子碰撞导致的分子三体解离方面,Wang等[14]利用较高能量(500 eV)电子碰撞研究了三电荷分子离子的碎裂过程,发现次序解离过程主要发生在基态和较低激发态.Wang 等[15]利用200 eV 电子碰撞研究了二价离子的两体和三体解离通道,通过对比高入射能量下的实验结果,发现不同入射电子能量下离子状态的变化对碎裂机制的影响.高电荷态离子也常常被选作炮弹离子用于分子碎裂动力学的研究.Neumann 等[4]研究了低能(3.2-keV/u)高电荷Ar8+离子碰撞导致离子三体解离过程,通过测量碎片动量区分了直接碎裂、次序解离等不同机制,并发现沉积在系统中的总能量是决定解离机制的关键因素.Yan 等[7]利用56-keV/u 的Ne4+离子研究了三体解离过程,通过分析碎片离子的动量关联和动能分布,鉴别了次序解离和直接解离的解离机制,并根据次序过程中CO2+亚稳态离子解离释放的能量确定了所对应的电子态.

碳氢化合物广泛存在于自然界,且在工业生产中起着重要作用,其解离碎裂机制研究不仅具有重要的基础意义,而且在等离子体物理[1]、星际介质演化[12]等领域具有潜在应用价值.碳氢化合物分子由四个或以上原子组成,其解离过程相比于三原子分子更为复杂,在解离机制鉴别方面增加了新的挑战.丙炔和丙二烯是自然界中存在的具有稳定同分异构体的最小分子,常被选为模型分子研究有机分子解离碎裂机制及同分异构效应.Kusakabe等[16]在电子碰撞导致的C3H4分子解离过程研究中观测到明显的同分异构效应,并发现丙二烯结构异构化至丙炔结构的过程存在.Scully 等[17]在C+离子碰撞导致的解离碎裂研究中同样发现了明显的同分异构效应.Mebel 与Bandrauk[18]和Psciuk等[19]在理论上对C3H4阳离子可能存在的解离路径进行了模拟计算.Yamanouchi 课题组[20−24]利用超快强场激光诱导丙二烯发生库仑爆炸,研究了丙二烯的两体与三体解离过程,并发现了氢迁移过程的存在;Ma 等[25]研究了Ne8+离子碰撞导致的丙二烯三价离子三体碎裂解离过程,发现直接解离过程是三体解离通道的主要机制,次序解离过程的贡献相对较少.

分子碎裂过程中既会产生带正电荷的离子碎片,也可能产生不带电荷的中性碎片,对包含中性碎片的通道的研究是认识分子碎裂机制必不可少的部分.但是,关于丙二烯(CH2CCH2)和丙炔(CH3CCH)解离碎裂中涉及中性碎片过程的研究还未有详细报道.这是由于在实验中探测器往往只能探测到带电粒子,而对中性碎片的直接探测比较困难.在本工作中,利用中国科学院近代物理研究所的反应显微成像谱仪装置,研究了Ne8+诱导丙炔和丙二烯分子三体解离过程.通过碎片离子动能及解离的动能释放分析鉴别出第三个碎片为中性H 原子的事件,进一步结合Newton 图、Dalitz 图及碎片离子角度分布等,对该通道的解离机制进行了详细分析.

2 实验方法

本实验在中国科学院近代物理研究所320 kV高电荷态离子综合实验平台的离子与原子分子碰撞动力学实验终端完成[3,25−28].电子回旋共振离子源(ECR)产生的50-keV/u Ne8+高电荷态离子束先后经过二级磁铁偏转、可调光阑限束及静电偏转后,得到高品质的炮弹离子束,进入反应谱仪(reaction microscope spectroscopy)实验终端,并在反应区与超音速气体喷嘴产生的C3H4分子束交叉碰撞.碰撞产生的离子碎片被垂直于入射炮弹和气体束方向的均匀电场(180 V/cm)加速引出,通过无场漂移区后由反冲离子探测器收集.碰撞之后的散射离子通过静电偏转板与未反应的束流分离后被散射探测器探测,未参与反应的Ne8+离子则被法拉第筒收集.反冲和散射探测器由微通道板和延迟线阳极组成,记录被测离子的二维位置(x,y)和到达时间(t)信息.测量中将炮弹离子到达散射探测器的时间作为触发信号,用于确定碎片离子的绝对飞行时间 (time of flight,TOF).实验选取C3H4的两种同分异构体:丙二烯和丙炔作为研究对象,在相同的实验参数下进行测量.C3H4分子解离末态存在多种碎片离子,在离线数据分析中,通过TOF 谱上离子飞行时间信息区分离子种类.

本研究关注以下两个反应通道:

通过动量计算出离子的动能及解离过程的动能释放(kinetic-energy release,KER,定义为3 个碎片的动能之和).

3 结果与讨论

在实验测量中,H+++H 通道中的中性H 原子不能被直接探测,该通道通过H+和离子的二维飞行时间符合信息进行鉴别,如图1(a)和图1(b)红色椭圆区域所示.这将不可避免地混入三价母体离子解离产生两个质子和离子的事件.由于两种事件所对应的母体离子价态不同,将导致末态碎片离子携带不同的动能,三价母体离子将对应更大的末态碎片动能及KER.为了进一步区分母体离子为二价和三价的事件,首先对末态碎片的动能进行详细分析,其中对于未测量的第3 个碎片(H 原子或质子),首先根据动量守恒得出其动量,进一步计算出其能量.

图1 (a) CH2CCH2 和(b) CH3CCH 二维飞行时间谱[8],其中红色椭圆对应H++的事件;(c) CH2CCH2 和(d) CH3CCH分子的KER-H+能量的二维符合谱;(e) CH2CCH2 和(f) CH3CCH 三体碎裂→H+++H 通道的KER 分布;(g) CH2CCH2 和(h) CH3CCH 对应的H+和未探测粒子动能(H+/H0)的二维符合谱Fig.1.Two-dimensional time-of-flight (TOF) spectra,in which the events of red oval corresponds to the H++ for(a) CH2CCH2 and (b) CH3CCH; two-dimensional coincidence correlation spectra of KER-H+ energy of (c) CH2CCH2 and(d) CH3CCH;the KER distribution of (e) CH2CCH2 and (f) CH3CCH for three-body fragmentation of →H+++H channel;two-dimensional coincidence spectra of kinetic energies between H+ and undetected particle (H+/H0) of (g) CH2CCH2 and(h) CH3CCH.

图1(c)和图1(d)为KER 和质子能量的二维符合谱.从图1(c)和图1(d)可以看出,对于丙二烯和丙炔两种分子,二维符合谱上都存在3 种结构,其中结构I 对应的KER 较小,II 和III 对应的KER 较大.图1(e)和图1(f)分别为CH2CCH2和CH3CCH 分子的KER 分布,即图1(c)和图1(d)在x轴方向的投影.两种分子的KER 都呈双峰结构,其中低能峰位于5.5 eV 左右,来源于二维谱中的结构I,能量较高的峰位于16.5 eV 左右,来源于结构II 和III.根据图1(e)和图1(f)初步判断结构I 来自于二价母体离子解离,结构II 和III 来自于三价母体离子的解离过程.

在图1(g)和图1(h)中,进一步通过H+离子与未探测粒子(H+/H0)动能的二维符合谱对来自不同结构的事件进行分析.对于结构II,H+动能远大于结构I 和III,进一步表明结构II 对应的H+离子来自较高价态的离子解离,即三价母体离子的解离过程.而结构I 和III,H+离子的动能集中在3—6 eV,表明这两结构对应的H+离子来自较低价态的离子解离(二价离子).不同的是,结构III未探测粒子动能(12 eV)远大于结构I (0—2 eV),说明结构III 中未探测粒子来自更高价态(三价母体离子)的解离过程.而结构I 未被探测的粒子的动能集中在0—2 eV,这是由于其不携带电荷,动能不是来自库仑爆炸过程,结合以上分析表明,结构I 对应于→H+++H 的事件,结构II 和III 则对应于→H+++H+的次序解离过程.在后面的讨论中,以图1(c)和图1(d)中的红色实线所包围的区域作为条件,选择出→ H+++H 的事件,深入分析该通道的动力学机制.

反应路径(I)为同步解离过程,路径(II)—(IV)为次序解离过程.在同步碎裂过程(I)中,当C3H42+离子发生库仑爆炸时,中性氢原子由于没有库仑排斥作用,其平均动量远小于另外两个荷电碎片.而H+和由于相互排斥将获得大小相当的动量.对于次序解离路径(II)—(IV),在第一步中所对应的亚稳态中间体离子分别为,,中间体离子存在一定时间之后发生第二步解离.

为了阐明该通道的形成过程,采用广泛用于三体碎裂机制研究的Dalitz 图和Newton 图展示3 个碎片离子的动量关联,进而对解离机理进行详细分析.图2 给出了CH2CCH2和CH3CCH 分子的Dalitz 图.在本研究中,Dalitz 坐标X和Y被定义为

其中a,b,c分别代表3 个碎片离子 H+,C3H2+和H,Pi(i∈{a,b,c}) 表示每个碎片的动量大小.

从Dalitz 图看到明显的宽带状结构,这是由于中间体离子在进一步解离之前发生转动所造成的典型次序解离特征[5],此结构表明直接解离路径(I)不是该通道形成的主要原因.从图2 可看出,大多数事件对应的H 原子的动量比H+和碎片的动量小得多,即.通过动量关系,次序解离路径(II)的贡献可以排除,原因如下:在次序解离路径(II)中,第二步解离时分解为两个质量相同的H+和H 碎片,可以预期中性粒子H 的动量应该与质子H+相当,即pH+≈pH,这与Dalitz 图中呈现的动量关系不符.

图2 (a) CH2CCH2 和 (b) CH3CCH 三体解离 →H+++H 过程的Dalitz 图Fig.2.Dalitz plot for three-body dissociation channel →H+++H of (a) CH2CCH2 and (b) CH3CCH.

为了进一步确认以C3H3+作为中间体离子的次序解离路径(IV)是实验观测到的H+++H 通道的主要贡献,进一步在图3 中给出该通道的Newton 图.如图3 所示,将H+离子的动量定义为1 arb.unit,碎片粒子和 H 的动量矢量归一化到H+离子动量,分别绘制在Newton 图的上、下半平面上.其中图3(a)和图3(g)分别是CH2CCH2和CH3CCH 所有H+++H 事件的Newton 图.在次序解离路径(IV)中,H+离子在第1 步解离中释放,然后离子分解为离子和中性H 原子.在第2 步解离过程中,离子在第1 步获得的动量由和H 碎片分配,分配比等于碎片离子的质量比.从图3(a)和图3(g)可以看到,和H 原子动量分布中心位置与其质量比符合.

在此基础上,以中性 H 碎片的能量范围为条件,给出不同H 原子能量条件下的Newton 图,分段范围为:0—0.5 eV (b),(h);0.5—1.0 eV (c),(i);1.0—2.0 eV (d),(j);2.0—4.0 eV (e),(k);4.0—6.0 eV (f),(l).从分段Newton 图中可以清晰地观察到两个半圆形结构,这是次序解离过程的典型特征.由于在第1 步解离过程中,出射的质子与C 链并不在一条直线上,中间体离子在与质子的库仑作用下会获得一定的角动量而发生旋转,旋转的离子在与质子相互远离的过程中进一步发生第二步解离,产生与中性H 原子,从而形成了Newton 图中上下两个半圆环结构.同时注意到图3(a)与图3(b)和图3(g)与图3(h)并不是典型的半圆型结构,这是由于随着H 原子能量的降低半圆环直径连续减小至0 造成的.

为了进一步确定次序解离通道(IV)的贡献,分析次序解离中两步过程所发生方向的夹角α的分布.对于通道(IV),α定义为第1 步出射的H+的动量(PH+)与在中间态离子质心系中和中性H 的动量之差之间的夹角.对于完全次序解离过程α为正弦分布[5,13].

假定H+++H 通道全部来自于次序解离路径(IV),即H+为第一步出射的粒子,图4(a)和图4(b)分别为CH2CCH2和CH3CCH分子所对应α角分布.从图4(a)和图4(b)可以看出,两个分子的α分布都与正弦分布(蓝色实线)符合.同时,以H 的能量作为条件,在图4(c)和图4(d)中分别给出不同能量段的α分布,分段范围依次为 0—0.5 eV,0.5—1.0 eV,1.0—2.0 eV,2.0—4.0 eV,4.0—6.0 eV.从图4(c)和图4(d)可以看出,对于每一个能量区间,两个分子的α角均服从正弦分布.值得注意的是,在0—0.5 eV 范围内同样满足正弦分布,进一步表明Newton 图3(b)和图3(h)的结构是由于H 原子能量降低时半圆环直径连续减小而形成的连续区域,其同样来自于次序解离路径(IV).α角分布进一步证明次序解离解离通道(IV)是和CH3CCH2+两种同分异构体三体解通道H+++H 的主要贡献.

图3 CH2CCH2 (a)—(f) 和 CH3CCH (g)—(l)三体解离 → H+++H 过程的Newton 图 (a)和(g) 包含所有事件,其他为以第3 个粒子H 的能量为选择条件的Newton 图,相应能量分别为0—0.5 eV (b),(h);0.5—1.0 eV (c),(i);1.0—2.0 eV (d),(j);2.0—4.0 eV (e),(k);4.0—6.0 eV (f),(l)Fig.3.Newton diagrams of CH2CCH2 (a)−(f) and CH3CCH (g)−(l) for three-body fragmentation channel →H+++H.(a) and (g) are Newton diagrams for all the events.The others are for different energy ranges of the neutral H:0−0.5 eV for (b),(h),0.5−1.0 eV for (c),(i);1.0−2.0 eV for (d),(j);2.0−4.0 eV for (e),(k);4.0−6.0 eV for (f),(l).

4 结论

本文研究了丙二烯和丙炔与50-keV/u Ne8+离子碰撞导致的三体解离过程,选取在二维飞行时间谱中观察到的不完全探测通道H++,利用反冲离子的飞行时间和在探测器上的位置信息还原了解离所产生的碎片离子的三维动量,在此基础上根据动量守恒获得第3 个碎片的动量信息.通过KER 及3 个碎片的动能分析,鉴别出未探测碎片为中性H 的事件,从而确定二价离子解离产生H+++H 的三体解离通道.通过3 个碎片离子的动量绘制了该通道的Dalitz 图、Newton 图,详细分析了该通道的发生机制,还给出了第一步和第二步解离方向的夹角(α角)分布,α角呈正弦分布,进一步表明次序解离路径→H++→H+++H 为该通道发生的主要途径.通过对比两种同分异构分子的Dalitz 图、Newton 图和角分布,表明分子结构不是影响此通道的主要因素.

感谢中国科学院近代物理研究所320 千伏高压平台工作人员的技术支持.