液相金纳米粒子非线性光散射方程的研究

  • 格式:pdf
  • 大小:208.43 KB
  • 文档页数:6

第6卷第4期分析测试技术与仪器Volume6Number42000年12月ANALYSISANDTESTINGTECHNOLOGYANDINSTRUMENTSDec.2000研究报告(200~205)液相金纳米粒子非线性光散射方程的研究蒋治良(广西师范大学新技术新材料研究所,广西桂林 541004)摘 要:将激励光与液相金纳米粒子相互作用的运动方程表达为朗之万方程,对金纳米粒子的共振散射、2/3分频共振散射和2倍频共振散射进行了分析,较好地解释了液相金纳米粒子产生的一些非线性散射现象.关键词:液相金纳米粒子;朗之万方程;非线性共振散射中图分类号:O657.32 文献标识码:A 文章编号:1006-3757(2000)04-0200-06 差频、和频、倍频(谐波)等非线性光学原理已用于激光、化学等研究领域[1~5].谐波辐射一直是光学界许多科技人员从事的主要研究内容之一.利用谐波可获得新的激光谱线,可把激光波长推向更短,因此,谐波辐射研究长盛不衰[4].倍频原理在分析化学中也获得了较好的应用效果,建立了一系列高灵敏、高选择性的无机分析新方法[6~10].最近我们采用微波高压液相合成法和光化学法制备了Ag、Au、Ag/Au、Ag/AgCl、Pt等液相纳米粒子,并研究了其共振、倍频和分频散射光谱.从超分子界面能带理论出发,建立了液相纳米粒子的和频、差频和分频原理,较好地解释和预测了液相纳米粒子、超分子、细菌、细胞等的非线性散射现象[11~16].但光(特别是低强度单色光)与液相纳米粒子作用的非线性散射方程的研究尚未见报道.马克思认为,任何一门科学只有它充分应用了数学时才能算作很好地发展了[17].本文企图从光与金纳米粒子相互作用的朗之万(Langevin)方程出发,探讨液相纳米粒子的共振、2/3分频和2倍频散射的数学规律.1 实验部分1.1 仪器与试剂RF-540型荧光分光光度计(日本岛津);U-3400型紫外可见分光光度计(日本日立);H-600型透射片镜(日本电子株式会社);Glanz微波炉(中国顺德,800W,2450MHz).1.2 实验方法量取1.0mL228Lg/mLAu3+溶液于微波反应罐中,再加入0.5mL2.0@10-3mol/L柠檬酸三钠溶液,然后加蒸馏水至10mL,扭紧罐盖,使其混匀.置微波炉中辐照3min,取出冷却后测量其散射光谱和粒径r.2 实验结果当金粒径较小时,溶液呈红色,且柠檬酸三钠的用量较大(柠檬酸在紫外光区有较强吸收).以上两因素均导致散射光大大减弱.因此采用粒径较大、柠檬酸三钠加入量较小的液相金纳米粒子体系进行实验.图1为液相金纳米粒子的吸收光谱,其在537nm和300nm处有二个吸收峰.图2为金纳米粒子的共振散射光谱,其最大共振散射波长为580nm(5.02@1014Hz).当Kex=580nm时,在580nm处产生一个主共振峰,在290nm(2@5.02@1014Hz)处产生一个2倍频共振峰,在870nm(2/3@5.02@1014Hz)处产生一个2/3分频共振峰(图3).3 讨 论图4表示入射光(X)沿X轴入射,X/2分频散射信号以90b垂直检测.当光与金纳米粒子发生非弹性碰撞后,其坐标发生变化.Yomosa把DNA分子视作基转模型[19],并给出激光与DNA作用的扭摆方程:A󰀁25/󰀁t2+Bsin󰀂=0(1) 我们将光强度作用量写为Fo,它是一常矢量;考收稿日期: 2000-08-04; 收到修改稿日期: 2000-10-15.

图1 液相金纳米粒子(r=76nm)的吸收光谱Fig.1 TheabsorptionspectraforAunanoparticleinliquid(r=78nm)a.22.8Lg/mLAu; b.11.4Lg/mLAu; c.5.7Lg/mLAu

图2 液相金纳米粒子(r=76nm)的共振散射光谱Fig.2 ResonancescatteringspectraforAunanoparticleinliquida.22.8Lg/mLAu; b.11.4Lg/mL

Au图3 金粒子的散射光谱Fig.3 ThescatteringspectraforAuparticlea.低灵敏度3档; b.低灵敏度6档虑到金纳米粒子本身的内在随机因素,以及外界环境等的随机因素,引入随机力Bo#(t);并考虑到金纳米粒子受线性阻尼的影响,阻尼记作C(󰀁󰀂/󰀁t).则在方程(1)的基础上,可得出金纳米粒子在入射光、随机力以及阻尼三者作用下的随机动力学方程,即朗之万(Langevin)方程:A󰀂tt+C󰀂t+Bsin󰀂=Bo#(t)+Fo(2)󰀂tt+C󰀂t+bsin󰀂=B#(t)+F(3)式中C=C/A为金纳米粒子所处水溶质环境中的阻尼系数;b=B/A由金纳米粒子体系内在特性决定,反映了金纳米粒子体系的内在特性;B=Bo/A为可调参数;#(t)为随机力;F为入射光强度;󰀂为金纳米粒子相对于平衡位置的角移.式(3)中F=EcosXt,E为入射光场对金纳米粒子中电荷和电偶极矩作用的电场力;t为时间.当无随机力且角移较小时,Bo#(t)=0,sin󰀂U󰀂-(1/6)󰀂3.令X=󰀂,Xc=󰀂t,xcc=󰀂tt,则(2)式化为:Xcc+CXc+bx-(b/6)x3=EcosXt(4)201 第4期 蒋治良:液相金纳米粒子非线性光散射方程的研究

图4 入射光和散射光坐标示意图Fig.4 Diagramofthecoordinateforincidentandscatteringlight3.1 朗之万方程的主共振当散射光振幅较小时,(4)式中的非线性项可以忽略,而化为一普通线性受迫振动方程Xcc+CXc+bX=EcosXt(5)其解为[20]X=Aexp(-C/2)cos(Pt+󰀂)+ Bcos(Xt+H)(6)(6)式中P=(b-C4/4)1/2(7)B=E/[(b-X2)2+C2X2]1/2(8)H=tg-1CX/(b-X2)(9)系统的固有频率Xo=AK/bo系统响应为式(6),第一项为齐次方程通解(自由振动项),第二项为特解.第一项带有决定于起始条件的常数A与󰀂,对于正阻尼C>0,自由振动项随时间衰减,故称为暂态响应.当时间足够长时稳定后所得响应称为稳态响应,仅由特解构成.因此稳态响应与入射光激励项有相同频率,但相位较激励项移动了H,H仅决定于阻尼及X.与X的相对值,与初始条件无关.从(8)式和(9)式可以看出,当XUXo时,即(6)式的稳态响应频率与固有频率相等时,响应振幅会很大,称主共振.这与实验结果(图2和图3)一致,即在金原子团簇的固有频率6.19@1014Hz处,振幅最大.当振幅大时,(4)式中的非线性项x3不能忽略,此时可用多尺度方法求(4)式的解[21].把表示响应解的展开式考虑为多个自变量(或多个尺度)的函数,而不是单个自变量t的函数.即把X看作t和Ent的函数,E为小量.因为考虑主共振时XUXÀ,令X=XÀ+ER,R为常数,为便于求解,将式(4)改写为 x+X20x+2EBx-EAx3= Eecos[(X0+ER)t](10)其中X20=b,2EB=C,EA=b/6,Ee=E.设T0=t,T1=Et,,,则d/dt=(dT0/dt)/T0+(dT1/dt)/T1+,,=D0+ED1+,,, d2/dt2=D20+2xED0D1+,,(11)202 分析测试技术与仪器 第6卷 这里仅选取时间尺度为T0,T1,则式(10)的解用不同时间尺度表示为式(12). x(i;E)=x0(T0,T)+ Ex1(T0+T1)+,,(12)把式(12)代入方程(10),令等式两端的E0和E1的系数相等,则得 D20x0+E20x0=0(13)及 D02x1+X20x1+2D0D1x0+2BD0x0-Ax03= ecos(X0T0+TT1)(14) 式(13)的通解可写为 x0=A(T1)exp(iX0T0)+ Ax(T1)exp(-iX0T0)(15)将式(15)代入式(14),且cos(X0T0+RT1)用复数表示,则得 D20x1+X20x1= -[2iX0(Ac+BA)-3AA2Ax] exp(iX0T0)+AA3exp(3iX0T0)+ (1/2)eexp[i(X0T0+RT0)]+cc(16)式中Ax为A的共轭复数,Ac为A关于T1的导数,cc代表前面各项的共轭复数.使式(16)两边的exp(iXoTo)的系数相等,就得到关于A的方程. 2iXo(Ac+BA)-3AA2Ax- (1/2)eexp(iRT1)=0(17)为便于解方程,将A表示成极式A=(1/2)Cexp(iH)(18)式中C与H都是实函数,设G=RT1-H(19)将式(18)代入式(15),再将结果和式(19)代入式(12),得稳态解的一次近似x=Ccos(XtG)+Q(20)其中Q为解的展开式的余项.从式(20)可以看出响应频率X的共振频率在XUXo处.下面进一步分析,将式(18)代入式(17),分成实部和虚部,并使实部、虚部分别相等,则得 Cc=BC+(e/2X0)sin(RT1-H)CHc=(3/8)(AC3/X0)-(e/2X0)cos(RT1-H)(21)把式(19)代入式(21),可以把式(21)变换成一个自治系统(即不显含T1),则 Cc=BC+(e/2X0)sinG GcC=RC+3RC3/8X0+(e/2X0)cosG(22)对于稳态运动,C与G是常数,T1导数为零,即Cc=Gc=0,则式(22)成为 CB=(e/2X0)sinG CR+3RC3/8X0=-(e/2X0)cosG(23)方程两边平方后相加,得频率响应方程为B2+(R+3RC2/8X0)2=e2/4X20C2(24)响应幅值C与激励幅值e、频率差值R以及非线性项系数等三因素有关,当存在阻尼时,这种关系呈非线性.图5绘出了C与R关系曲线.频率响应曲线的弯曲导致了振幅的多值性,从而导致了跳跃现象.当X远小于X0开始增加,响应的幅值沿AF曲线慢慢增加直到B,在B点只要X有微小增加,就会突然从B点跳跃到C点.随X继续增加,响应幅值沿曲线从C点慢慢减小,再沿DCE曲线慢慢增加.在E点X值的微小减少都会突然发生自发向下跳跃,从E点到F点.X进一步减少,响应幅值继续从F点减小到A点.这表明了大的激励频率,也即光子能量大,不一定引起大的响应.因每一种液相纳米粒子均有其固有的振动频率或泵浦能带.从图5可以看到,激励频率比X0小时,在某一范围,某些条件下,引起响应反而大.这在线性关系203 第4期 蒋治良:液相金纳米粒子非线性光散射方程的研究 中是很难理解的,产生这种现象是系统非线性造成的.

图5 激励光与纳米粒子相互作用系统的主共振频率响应曲线(A=1.0,d=1.0,B=0.5)Fig.5 Thefrequency-responsecurveforprimaryresonanceof theligh-tAunanoparticleinteractionsystem(A=1.0,d=1.0,B=0.5)3.2 朗之万方程的次共振非线性介质的另一特征是分频和倍频共振,当X远离Xo时,式(4)可改写为xcc+X20xc+2EBx-EAx3=EcosXt(25)同理,采用多尺度方法求式(25)的近似解,当3X/2=RE+X0时,其响应的一次近似为x=Ccos(3Xt-G)+2dcosXt+Q(26)式中C和G为常数,d=E/[2(X20-X2)].尽管存在阻尼,在XU(1/3)X0时,其自由振动项(即第一项)并不衰减到零.非线性性质调整了自由振动项共振频率在XU(2/3)X0处.解得其频率响应方程为 B2+(R+3Ad/2X0+3AC2/8X0)2= A2d3/X20C2(27)在XU2X0/3附近,响应幅值与频率差值R的关系如图6.从图6可以看出,XU2X0/3时也能引起共振,此共振称为2/3分频共振.说明激励频率小于系统的固有频率也能引起系统的共振,这是非线性带来的结果.