二维对六苯分子团簇中相干激发能量转移过程研究.pdf

二维对六苯分子团簇中相干激发能量转移过程研究*方德音范旭阳卫岸王鹿霞(北京科技大学数理学院,理论物理研究所,北京100083)(2023年 4月 10 日收到;2023年 7月 28 日收到修改稿)激发能量转移是影响对六苯相关器件应用的重要因素之一,研究对六苯团簇在不同外场激发作用下的激子动力学过程及激子相干效应,这对于提高由对六苯及其相关衍生物所组成的分子器件的应用性能具有重要意义.本文基于密度矩阵理论,采用算符的数学平均值近似,应用不同脉冲场激发盘状对六苯分子团簇,对其激子动力学过程及激子相干行为进行了分析.研究发现,当处于强场短脉冲激发时,团簇内产生多激子效应,激子能级出现新的杂化态.在短脉冲激发下,在不同能量区域激发态分布不同,但均表现出明显对称性.由于 H 型团簇最高能级的激子态被优先激发,发现随脉冲场增大,激子相干增强,当脉冲场较大时,在脉冲激发过程中,激子振荡协同性消失,极短时间内大量分子产生激子相干,出现瞬态离域现象.关键词：激发能量转移,多激子效应,激子相干尺寸PACS：73.23.b,78.20.Bh,73.22.LpDOI:10.7498/aps.72.202304761引言近年来,共轭聚集体和分子团簇广泛应用于多种光电器件中,由于电荷及能量转移效率直接或间接影响相关器件的应用,因此其电荷转移和能量转移动力学过程是研究的焦点问题之一16.对六苯(para-sexiphenyl,6P)分子79具有平面骨架结构,可以形成离域的 键,体系能量低,分子稳定性高,其衍生物可以组装成结构各异、功能独特的超分子结构,因此在有机薄膜晶体管1012、场效应晶体管13、有机发光二极管1416、光伏电池17,18等有机光器件19,20中具有广泛的应用.我们知道在弱场激发作用下,分子聚集体或分子团簇内一般形成单激子态,其相干性较好,但能量转移量低;在强场激发作用下,团簇内多个分子被同时激发21,能量转移量高,但由于分子间的相互作用,在团簇内形成非局域激发态,同时需要考虑激发能量转移和静电位移的影响,在外场激发过程中涉及多激子能带、多激子协同衰变和超辐射等过程,其相干行为比较复杂,深入理解其多激子激发及相干转移过程对提高相应分子器件的工作效率具有积极意义.近年来由于分子光刻和金属纳米腔技术的发展22,23,在局域空间内形成强场已经成为可能,强场激发下分子或分子团簇的吸收、荧光和发射谱实验研究2427已经成为当前的热点,但相应的理论研究相对少得多,主要原因是多激子态在分子团簇内形成的激发位形较为复杂,对其准确描述需要耗费大量的计算资源,且多激子态的量子动力学过程比较复杂,结果难与实验比较.本工作基于密度矩阵理论框架,将局域外场激发或表面等离激元局域场激发都化为激光脉冲激发模型,应用算符的数学平均值近似21,28,29,考虑分子上的平均激子占据,*国家自然科学基金(批准号:21961132023,11774026)资助的课题.通信作者.E-mail:2023中国物理学会ChinesePhysicalSocietyhttp:/物理学报ActaPhys.Sin.Vol.72,No.19(2023)197301197301-1不考虑多激子占据的具体位形和表面等离激元引发的量子场效应,大大缩减计算量的同时能够正确反映多激子态的量子动力学过程,得到相干能量转移信息.根据 Kasha 理论,分子聚集体和团簇可分为J 型和 H 型.J 型聚集体中偶极子头尾相连,呈现负激子耦合,聚合物较单体吸收峰红移,出现荧光增强及超辐射现象.而 H 型聚集体中偶极子平行排列,表现正激子耦合,吸收光谱较单体蓝移且出现荧光猝灭现象4,30.近年来,随着具有荧光性质的染料分子出现,荧光型 H 聚集体成为可能,是新一代纳米材料的代表3133.且 H 聚集可产生强相干激发,在非线性光学等领域应用广泛34,因此研究 H 型聚集体的能量及电荷转移行为尤为重要.Efg Eeg对于 H 型 6P 分子团簇的理论研究,主要集中在一维分子链的激发能量转移21及其与无机界面间的电荷转移35,36,但许多有机光器件中 6P 薄膜为直立分子组成的平面团簇,目前对二维结构中多激子行为的研究较少.本工作以具有平面结构的6P 分子团簇作为研究对象,对其内部激发能量转移过程及相干行为进行研究.由于其高阶激发态的能量远大于两倍的第一激发态的能量,即 ,我们将其简化为二能级系统,忽略激子-激子湮灭过程37,38对动力学过程的影响.第 2 节介绍理论模型,第 3 节针对不同激发场下多激子态布居数演变过程和团簇的多激子过程展开讨论,第 4 节对本文进行了总结.2理论模型和动力学方程B+m采用二能级模型对 6P 分子电子态进行描述,引入单分子跃迁算符 :B+m=|memg|,(1)|mg|meB+mBm+BmB+m=|meme|+|mgmg|=1其中 为分子基态,表示分子第一激发态,在单电子模型下满足关系:.这样 6P 分子团簇的哈密顿量可以写为Hmol=m(Em+nVmnB+nBn)B+mBm+m,nJmnB+mBn,(2)Em其中,表示分子团簇中第 m 个分子的激发能,可写为Em=Eme Emg+nJmn(ee,gg)Jmn(gg,gg),(3)Vmn式中,求和部分为分子激发能修正项,表示其余处于基态分子与分子 m 库仑耦合.第 2 项 表示当多个分子被同时激发时,不同分子间的库仑耦合所引起的静电位移对分子激发能的修正,可写为Vmn=12Jmn(gg,gg)Jmn(ee,gg)Jmn(gg,ee)+Jmn(ee,ee),(4)Jmn其中,表示分子共振跃迁所引起的激发能量转移耦合,可写为Jmn=Jmn(eg,ge),(5)这里忽略了其他非共振跃迁过程的影响.外场与分子的耦合项写为Hfield(t)=E(t)mdmB+m+H.c.mRm(t)B+m+H.c.,(6)dmE(t)其中 为第 m 个分子的偶极矩.激光脉冲场 的形式为E(t)=nE(t)ei0t+c.c.,(7)0E(t)其中,n 是激光场的方向矢量,是激发场频率.脉冲场包络 为E(t)=E0exp4ln2 (t tp)2/2p,(8)pE0其中,为脉冲宽度,为脉冲场强度.系统的哈密顿量为H(t)=Hmol+Hfield(t).(9)在密度矩阵理论框架下,应用量子主方程可以描述多电子同时激发的动力学过程.若按照传统方法计算,需要考虑各个激发态在分子团簇内的占据位形.计算不同位形下对应的密度矩阵元的动力学过程,计算量随分子团簇中分子数目的增加呈指数增长,我们采用算符的数学期望值近似21,28,29,只考虑每个分子上电子激发的平均值而不考虑激发态占据的具体位形,研究其动力学过程.系统的量子主方程为t(t)=iH(t),(t)D(t),(10)(t)D(t)其中 为含时密度矩阵算符,等式右边第 1 项为系统内激子相干项.表示该系统耗散,物理学报ActaPhys.Sin.Vol.72,No.19(2023)197301197301-2DO 为耗散超算符,对于任意算符 ,在数学期望值近似下,有O=tr(t)O,(11)DD方程(10)和方程(11)中的耗散超算符 修正为,量子主方程为tO=iH(t),O DO=(Ot)coh+(Ot)diss,(12)其中分子耗散由分子内转换耗散、分子间激发能量转移引起的振动耗散以及分子退相干耗散三部分组成,分别具有以下形式:DICO=12mm(B+mBm,O+2B+mOBm),(13)DEETO=12m,nmn(B+mBmBnB+n,O+2B+mBnOB+nBm),(14)DpdO=m(pd)m(BmB+mOB+mBm+B+mBmOBmB+m).(15)mmn(pd)m这里,为第一激发态的衰变速率,为相邻分子间振动激发能量转移速率,为分子基态激发态间的退相干速率.引入分子平均激发态布居数:Pm(t)=Pme(t)=tr(t)B+mBmB+mBm,(16)Pmg(t)=BmB+m=1 Pme(t),(17)Am(t)Wmn(t)以及辅助量 和 :Am(t)=B+m,(18)Wmn(t)=B+mBn Am(t)An(t).(19)m=Em/vmn=Vmn/jmn=Jmn/rm=Rm/m(pd)mmnm(pd)m引 入 ,由于讨论的动力学过程在皮秒量级,且激子衰变速率 和退相干速率 远小于振动激发能量转移速率 ,所以忽略 项以及 项,得到以下动力学方程:tAm=imAm+ik(1 2Pm)jmkAk+irm(1 2Pm),(20)tPm=2Imkjmk(Wmk+AmAk)+2Im(rmAm)kkmk(Pm Pk),(21)tWmn=i(mn)12kmn(1 Pn)12knm(1 Pm)Wmn ijnm(Pm Pn)+ik=m,njkmWkn(1 2Pm)jnkWmk(1 2Pn)+ijnmAmAm(1 2Pn)ijnm(1 2Pm)AnAn12(kmn+knm)AmAn.(22)m=m+2nmnPni2nkmn(1 Pn)+knmPn.(23)结合方程(3)和方程(4)给出分子激发能的具体形式:m+2nmnPn=m+njmn(ee,gg)Png+jmn(ee,ee)Pne jmn(gg,ee)Pne jmn(gg,gg)Png,(24)nmnPn式中,表示实际激发所引起的库仑位移对分子激发能的校正.团簇处于稳态时激发态的电子占据为P(ss)tot=mP(ss)m,(25)P(ss)m其中 是第 m 个分子达到稳态时的激发态布居.为研究多激子团簇的相干行为,引入相干尺寸39:L(t)=(mn|mn(t)|)2Nmolmn|mn(t)|21,(26)物理学报ActaPhys.Sin.Vol.72,No.19(2023)197301197301-3mnNmolL=1其中,表示密度矩阵内所有对角及非对角元素,为团簇内的分子数.当激子间无相互作用时 ,表明分子团簇内无相干性.3结果与讨论对六苯分子在常温下呈人字形排列9,选取由59 个 6P 分子所组成的圆盘状团簇作为研究模型,mn=nm=如图 1 所示.其偶极矩沿分子长轴方向,表现出H 型聚集体特征4.在脉冲激光场作用下,二维团簇中多个分子被同时激发,形成多激子态,根据相应的量化计算21,圆盘型 6P 分子团簇静电位移及激子能带如图 2 所示,其激子能级主要分布在 45eV 区域,其他参数列在表 1 中,设不同激子间能量转移速率相同 .3.1 光学响应过程apap=E0ppap=5105psVm1p6P 分子团簇的激发和转移过程与外场的能量和脉冲宽度紧密相关.为了深入研究不同脉冲场下,二维 6P 分子团簇的激子能带特征,计算了脉冲场能量 ()一定时,不同脉冲宽度 下,6P 分子团簇的光学响应过程,如图 3 所示.图 3(a)中脉冲场能量较低,发现随脉冲宽度 的减小,相应峰的峰值增大,展宽增强,但在 4.2eV 附近的低能频段,光学响应峰的峰结构表1参数设置Table1.Parametersused.名称/单位数值Nmol分子数 59Emol分子本征态能量 /eV4.02dmol分子偶极矩 /D12.4分子间激发能量转移速率 /(1011s1)20脉冲场频率 E0/随 变化E0脉冲场强度 /(Vm1)105108p脉冲宽度 /fs201000(a)(b)(c)图1二维对六苯团簇结构图(a)俯视图;(b)沿 x 轴侧视图;(c)沿 y 轴侧视图Fig.1.Structureofthetwo-dimensional6Pclusters:(a)Topview;(b)sideviewalongxaxis;(c)sideviewalongyaxis.5.0(b)4.84.6Frenkel激子能带/eV4.44.24.020100-10-20-20-1001020A2468101214相对静电位移/meV(a)图2二维对六苯团簇静电位移示意图(a)分子激发能的相对静电位移;(b)分子激发能无(红色)及有(蓝色)静电位移时的Frenkel 激子能带21Fig.2.Schematicdiagramoftheelectrostaticdisplacementoftwo-dimensional6Pclusters:(a)Relativeelectrostaticshiftofthemolecularexcitationenergies;(b)Frenkelexcitonenergybandwithout(red)andwith(blue)electrostaticshiftofthemolecularexcitationenergies21.物理学报ActaPhys.Sin.Vol.72,No.19(2023)197301197301-4ap=2.5106psVm1p发生变化,由分立的双峰变为有较大展宽的单峰.说明相同脉冲场能量下,脉冲宽度越小,频域内的展宽越大,可以覆盖更多的能级区域,更多的多激子态能级受到激发.图 3(b)给出 强脉冲场下 6P 分子团簇的光学响应过程.与图 3(a)比较可以发现,在脉冲宽度 较大时,强场激发与弱场激发对应的光学响应相差不大,但随着脉冲宽度的减小,强场激发下的光学响应峰出现明显的红移且出现多个共振峰.说明在强场情况下,外场与激子间相互作用产生新的杂化态,表现出明显的多激子效应.为了进一步研究场强对光学响应的影响及其多激子效应,我们给定脉冲宽度,探究不同脉冲场能量下,6P 团簇的光学响应过程.p=1 psp=50 fs图 4(a)对应 的长脉冲场,图 4(b)对应 的短脉冲场.由图 4(a)可发现当脉冲宽度较大时,共振峰结构不随脉冲场能量的增大而改变,峰值对应的能量与图 3(b)中的激子能带结构大致吻合.由图 4(b)可以看到,当脉冲宽度较小时,随着脉冲场能量的增大,高能频段的共振峰数目增多,且峰值红移明显.说明在短脉冲场作用下,4.04.24.44.64.85.000.140.2800.20.400.350.7000.61.200.71.4012Energy/eVtot(ss)(a)1000 fs500 fs200 fs100 fs50 fs20 fstot(ss)(b)4.24.44.64.85.0Energy/eV1000 fs500 fs200 fs100 fs50 fs20 fs00.40.800.71.401.22.4024036051015=2 1011s1P(ss)tot0pap=5 105psVm1ap=2.5 106psVm1图3 时,团簇内激发态布居 随激光 、脉冲宽度 变化(a)弱场 ;(b)强场 P(ss)tot0p=2 1011s1ap=5 105psVm1ap=2.5 106psVm1Fig.3.Excited-statepopulations inclustersvarywithphotonenergy ofthelaserandpulselength at :(a)Inweakfields ;(b)instrongfields .tot(ss)04812Energy/eV4.04.24.44.64.85.05.2(b)p=5T105 psSV/mp=106 psSV/mp=2.5T106 psSV/mp=5T106 psSV/mtot(ss)Energy/eV4.04.24.44.64.85.05.200.51.01.5(a)p=5T105 psSV/mp=106 psSV/mp=2.5T106 psSV/mp=5T106 psSV/m=2 1011s1p=1 psp=50 fsP(ss)tot0ap图4 时,在长脉冲场(a)和短脉冲场(b)下,团簇内激发态布居 随激光 、脉冲能量 变化P(ss)totp=1 psp=50 fs0ap=2 1011s1Fig.4.Excited-statepopulations inclustersinlong (a)andshort (b)pulsedfieldsvarywithphotonenergy ofthelaserandpulseenergy ,.物理学报ActaPhys.Sin.Vol.72,No.19(2023)197301197301-56P 分子团簇的光学响应过程更多的受到脉冲场影响,随着脉冲场强度增大,激子激发程度增加,外场与激子间发生杂化效应,产生了更多的杂化能级.由于 H 型团簇的能级特点,在高能频段电子的布居数增加,多激子效应更为明显.3.2 激子动力学过程0=4.2483 eV4.9243 eVt=100 fs为了进一步研究 6P 分子团簇在高斯脉冲场下的多激子行为,应用宽度为 50fs 的短脉冲激发分子簇,图 5 和图 6 分别给出了 和 时,在 (图(a),200fs(图(b),500fs(图(c)和 1ps(图(d)对六苯分子的激发态布居.可发现在同一激发外场下,6P 团簇的激子局域区域随时间变化不大,对比不同激发态下6P 团簇的激发态布居,可以看到均表现出中心局域的特点,对不同激子态局域特点不尽相同,但由于 6P 团簇为二维对称结构,其激子态分布都具有对称性.为了简化描述并理解团簇内分子激发态动力学过程的特点,在团簇内选取 6 个典型分子,如图 7 所示,研究其动力学特征.0=4.9243 eVH 型团簇的特点是在光激发作用下电子优先占据高能激子态,图 8 给出了 6 个典型分子在 的短脉冲激发下,分子激发态布Pmap=5105ps V/m居 的动力学过程.由图 8(a)可以看出,在脉冲场激发作用下,对六苯团簇的激子布居数迅速增大,约 100fs 左右达到峰值,之后在激发能量转移的影响下激子布居数不断变化,约 6ps 到达稳态.进一步加大脉冲场(图 8(b)(d)时,研究发现,对六苯团簇激发态布居整体变化趋势与 (图 8(a)时的情况相近,但随着脉冲场增大,其达到峰值所需时间减少且出现 Rabi 振荡,Rabi 振幅随脉冲场增强而不断增大.当脉冲场足够大时,在前 100fs 激子振荡的协同性消失,激发局域性被破坏,出现非局域现象.00.0020.0040.0060.0080.010(a)(c)(d)(b)0=4.2483 eV=2 1011s1ap=5 105psVm1p=50 fsPmt=100 fst=200 fst=500 fst=1 ps图5 ,时,弱场短脉冲作用下(,),4 个不同时刻的分子激发态布居 (a);(b);(c);(d)Pmap=5 105psVm1p=50 fs0=4.2483 eV=21011s1t=100 fst=200 fst=500 fst=1 psFig.5.Molecularexcitedstatepopulations atfourdif-ferentmomentsintheweakfieldwithshortpulses(,),:(a);(b);(c);(d).0.010.06(a)(c)(d)(b)0.020.030.040.050=4.9243 eV=2 1011s1ap=5 105psVm1p=50 fsPmt=100 fst=200 fst=500 fst=1 ps图6 ,时,弱场短脉冲作用下(,),4 个不同时刻的分子激发态布居 (a);(b);(c);(d)Pmap=5 105psVm1p=50 fs0=4.9243 eV=21011s1t=100 fst=200 fst=500 fst=1 psFig.6.Molecularexcitedstatepopulations atfourdif-ferentmomentsintheweakfieldwithshortpulses(,),:(a);(b);(c);(d).654231图7二维对六苯团簇选取分子示意图Fig.7.Schematic diagram of selected molecules in two-dimensional6Pclustersclusters.物理学报ActaPhys.Sin.Vol.72,No.19(2023)197301197301-63.3 激子相干为了进一步研究其相干特性,与图 8 对应,图 9给出相应脉冲场条件下二维 6P 团簇激子相干尺寸随时间的变化.可以看到在场强较弱的时候,激子相干尺寸随时间呈现先增后减的趋势,随着脉冲场增强,在脉冲作用时间内(150fs),激子相干尺寸出现一个或多个峰,其峰值位置与图 9(c)和图 9(d)中的波谷对应,说明激子相干与激子布居的非协同振荡紧密相关.L激子相干是非局域激子的一种量子效应,表明激子处于不同分子激发的相干叠加态40.受脉冲激发及分子耗散的影响,对六苯分子的激子相干态不断变化,相应激子局域范围随之改变,当对六苯团簇处于稳态时,分子间的相干能量转移减少,激子态在每个分子上趋于全同.随着脉冲场增大,激子相干效应增强,增大,当脉冲场场强较大时,前 100fs 内大量分子在极短时间内发生激子相干,产生瞬态离域现象.4结论基于二维对六苯团簇的结构特点以及能级分布,将分子系统简化为二能级模型,采用算符的数学平均值近似,从量子主方程出发,研究不同脉冲场激发下圆盘状二维对六苯分子团簇的激子动力学过程,并针对其激子布居及激子相干行为进行讨论.研究发现,应用长脉冲激发二维对六苯团簇时,主要表现为单激子态特征,对激发外场的场强变化不敏感.当外场为强场短脉冲时,多个激子同时被激发,形成多激子态且有新的杂化态出现,在光学响应谱上出现了多个共振峰.强度增大,多激子效应明显,杂化能级增多.在短脉冲激发下,分子团0/fs101102103104123456(b)0.050.100(d)/fs1011021031041234560.060.120.18/fs10110210310400.020.040.06123456(a)(c)0/fs1011021031041234560.060.12=2 1011s1p=50 fs0=4.9243 eVPmap=5105psVm1ap=106psVm1ap=2.5 106psVm1ap=5 106psVm1图8 时,短脉冲(,)场下,分子激发态布居 随时间的演变(a);(b);(c);(d)Pmp=50 fs0=4.9243 eV=21011s1ap=5 105psVm1ap=106psVm1ap=2.5 106psVm1ap=5 106psVm1Fig.8.Molecularexcitedstatepopulations inshortpulsedfields(,)varywithtime,:(a);(b);(c);(d).()/fs1011021031040=1T107 V/m0=2T107 V/m0=5T107 V/m0=10T107 V/m0102030=2 1011s1p=50 fs0=4.9243 eVL图9在 ,短脉冲(,)场下,激子相干尺寸 随时间的演变Lp=50 fs0=4.9243 eV=21011s1Fig.9.Exciton coherence size in short pulsed fields(,)varywithtimeat .物理学报ActaPhys.Sin.Vol.72,No.19(2023)197301197301-7簇在不同时刻的激子态布居均表现为明显的对称性.随着脉冲场强增大,对六苯分子的激子态出现Rabi振荡,当脉冲场达到一定场强时,前 100fs内激子振荡协同性消失,表现非局域特点.从激子相干尺寸角度发现脉冲场增强使相干尺寸出现多个峰,说明激子相干随脉冲场增大而增强,当脉冲场达到一定场强时,脉冲激发过程中大量分子处于激子相干态,出现瞬态离域现象.参考文献 OCarrollDM,PetoukhoffCE,KohlJ,YuBX,CarterCM,GoodmanS2013Polym.Chem.451811ReinekeS,ThomschkeM,LussemB,LeoK2013Rev.Mod.Phys.8512452ZhangCX,LiuRG,MakCH,ZouXL,ShenHH,LeuSY,JiL,HsuHY2018J.Photonics Energy80210013SpanoFC,SilvaC2014Annu.Rev.Phys.Chem.654774BrixnerT,HildnerR,KohlerJ,LambertC,WurthnerF2017Adv.Energy Mater.717002365TamaiY,OhkitaH,BentenH,ItoS2015J.Phys.Chem.Lett.634176ReselR2003Thin Solid Films43317ZojerE,KochN,PuschnigP,MeghdadiF,NikoA,ReselR,Ambrosch-DraxlC,KnupferM,FinkJ,BredasJL,LeisingG2000Phys.Rev.B61165388PalczynskiK,HeimelG,HeydaJ,DzubiellaJ2014Cryst.Growth Des.1437919ChoiA,KwonYN,ChungJW,YunY,ParkJI,LeeYU2020AIP Adv.1002512710ZhangY,WeiXC,ZhangH,ChenX,WangJ2018Appl.Surf.Sci.42745211LiY,WangH,ZhangXH,ZhangQ,WangXS,CaoDF,ShiZS,YanDH,CuiZC2016RSC Adv.6537712KatiyarS,VermaN,JogiJ2022Semicond.Sci.Technol.3702500813YangXX,FengX,XinJH,ZhangPL,WangHB,YanDH2018J.Mater.Chem.C6887914UsluerO,DemicS,KusM,OzelF,SariciftciNS2014J.Lumines.146615SimbrunnerC,Hernandez-SosaG,BaumgartnerE,HesserG,RoitherJ,HeissW,SitterH2009Appl.Phys.Lett.9407350516JiangXX,DaiJG,WangHB,GengYH,YanDH200717Chem.Phys.Lett.446329Winder C,Andreev A,Sitter H,Matt G,Sariciftci N S,MeissnerD2003Synth.Met.13957318QianC,SunJ,KongLA,GouGY,ZhuML,YuanYB,HuangH,GaoYL,YangJL2017Adv.Funct.Mater.27160493319SimbrunnerC2013Semicond.Sci.Technol.2805300120WangL,PlehnT,MayV2020Phys.Rev.B10207540121LeeJ,ZhangQ,ParkS,ChoeA,FanZ,KoH2016ACSAppl.Mater.Interfaces863422LiGC,ZhangQ,MaierSA,LeiD2018Nanophotonics7186523Fernandes J D,Pazin W M,Aroca R F,Macedo W D J,TeixeiraSR,ConstantinoCJL2019Spectroc.Acta A21122124Zhang Y,Meng Q S,Zhang L,Luo Y,Yu Y J,Yang B,ZhangY,EstebanR,AizpuruaJ,LuoY,YangJL,DongZC,HouJG2017Nat.Commun.81522525HanXB,LiF,HeZC,LiuYH,HuHT,WangK,LuPX2022Nanophotonics1160326KramerSN,BrownJ,RiceM,PeteanuLA2021J.Phys.Chem.Lett.12591927WangL,MayV2016Phys.Rev.B9419541328WangL,MayV2018J.Phys.B5106400229KashaM1963Radiat.Res.205530ChakrabortyS,DebnathP,DeyD,BhattacharjeeD,HussainSA2014J.Photochem.Photobiol.A Chem.2935731WuFP,LuYX,MuXLE,ChenZT,LiuSS,ZhouXF,LiuSF,LiZB2020ACS Appl.Mater.Interfaces123238832ZhaoXJ,ZhaoHJ,WangS,FanZW,MaY,YinYM,WangW,XiRM,MengM2021J.Am.Chem.Soc.1432082833Gomez-Sosa G,Beristain M F,Ortega A,Martnez-ViramontesJ,OgawaT,Fernndez-HernndezRC,Tamayo-RiveraL,Reyes-EsquedaJA,IsoshimaT,HaraM2012Opt.Mater.3485634PlehnT,ZiemannD,MayV2018Phys.Chem.Chem.Phys.202687035Plehn T,Ziemann D,May V 2019 J.Phys.Chem.C 1222792536MayV2014J.Chem.Phys.14005410337HaderK,ConsaniC,BrixnerT,EngelV2017Phys.Chem.Chem.Phys.193198938ZhengFL,ChenLP,GaoJB,ZhaoY2021Materials14329139FassioliF,DinshawR,ArpinPC,ScholesGD2014J.R.Soc.Interface112013090140物理学报ActaPhys.Sin.Vol.72,No.19(2023)197301197301-8Coherent excitation energy transfer processes intwo-dimensional para-sexiphenyl molecular clusters*FangDe-YinFanXu-YangWeiAnWangLu-Xia(Institute of Theoretical Physics,School of Mathematics and Physics,University of Science andTechnology Beijing,Beijing 100083,China)(Received10April2023;revisedmanuscriptreceived28July2023)AbstractExcitationenergytransferisoneofthemostimportantfactorsaffectingtheapplicationsofpara-sexiphenyldevices.Thestudyofexcitondynamicsandexcitoncoherenceeffectofpara-sexiphenylclustersunderexternalfield excitation is important in order to improve the performance of molecular devices composed of para-sexiphenylanditsrelatedderivatives.Inthiswork,thetwo-dimensionaldisc-likepara-sexiphenylmolecularclusterisusedastheobjectofstudy.Themolecularsystemissimplifiedintoatwo-levelmodelbasedonitsstructuralfeaturesandenergyleveldistribution.Withintheframeworkofdensitymatrixtheory,theexcitondynamicsandexcitoncoherencebehaviorofdisk-likepara-hexaphenemolecularclustersexcitedbydifferentpulsefieldsareanalyzedthroughusingthemathematicalmeanvalueapproximationoftheoperator.Theresultsshowthatwhenlongpulsesareusedtoexcitepara-sexiphenylclusters,thesingleexcitonstatecharacteristicappearsandisinsensitivetothechangeofexcitedexternalfieldstrength.Whentheclustersaresubjectedtostrongpulsedfieldswithshortpulsewidths,multipleexcitonsareexcitedsimultaneouslyinthecluster,formingmultiple exciton states,with the exciton energy levels shifting toward lower energy and new hybrid statesappearing.Intheopticalresponsespectrum,thereappearmultipleresonancepeaks.Andasthepulsefieldisenhanced,themulti-excitoneffectbecomesapparentandthehybridizationenergylevelincreases.Undershortpulseexcitation,theexcitedstatesaredistributeddifferentlyindifferentenergyregions,butallofthemshowobvioussymmetry.Asthehighest-energyexcitonstatesofH-typeclustersarepreferentiallyexcited,weanalyzetheexcitonstatepopulationandtheexcitoncoherenceevolutionwithtimeinthehigh-energyexcitonstate.Withthepulsefieldincreases,Rabioscillationsappearandtheexcitoncoherenceeffectincreases.Whenthepulsedfieldreachesacertainfieldstrength,theexcitonoscillationcooperativitydisappearsinthefirst100fs,showing the non-local characteristic.The position of the wave trough of the exciton state populationcorrespondstothepeakintheexcitoncoherencesize.Itindicatesthatwhenthepulsefieldisintenseenough,alargenumberofmoleculesareintheexcitoncoherentstateduringthepulsedexcitation,andtransientout-of-domainphenomenaoccur.Keywords:excitatio