《2022年银纳米球及其阵列LSPR现象仿真分析》由会员分享，可在线阅读，更多相关《2022年银纳米球及其阵列LSPR现象仿真分析（10页珍藏版）》请在金锄头文库上搜索。
1、银纳米球及其阵列LSPR 现象仿真分析刘国华 , 孙辉, 张维，王程，王育剑，徐凯，岳钊(南开大学信息技术科学学院，天津，300071) 摘要 ：本文采用离散偶极近似DDA: Discrete Dipole Approximation和时域有限差分FDTD: Finite Difference Time Domain 方法对不同半径的银纳米球及其组合阵列的消光现象与耦合作用进行了分析讨论。用DDA方法得到了不同半径银纳米球及其阵列的消光光谱及其消光效率，吸收效率和散射效率等；在此基础上，利用FDTD 方法研究了双纳米球阵列的耦合作用及其电场分布；通过分析比较得到了一些关于纳米颗粒局域外表等离子
2、体共振LSPR：Localized Surface Plasmon Resonance 有意义的变化规律，并给出了它们的折射率灵敏度范围。同时对消光现象和耦合作用的物理机制也进行了分析，对金属纳米颗粒的应用研究具有参考意义。关键词：局域外表等离子体共振；离散偶极近似；时域有限差分；银纳米球Simulation and Analysis of Ag Nanosphere and Array s LSPR Phenomena Liu Guohua ,Sun hui , Zhang Wei,Wang Cheng,Wang yu jian, Xu kai,Yue Zhao(College of Inf
3、ormation Technical Science of Nankai University,Tianjin,300071)Abstract：The Discrete Dipole Approximation(DDA) method and the FDTD (Finite Difference Time Domain) method are used for analyzing different radius of silver nanoparticles and their couples and their arrays. DDA method is used for simulat
4、ing different radius and the array of silver nanoparticles extinction spectra and extinction efficiency absorption efficiency and the scattering efficiency; FDTD method is used for the coupling of two nanoparticles and simulating their around electric field distribution; by analyzing these, we have
5、gotten some important laws of nanoparticlesLocalized Surface Plasmon Resonance(LSPR) phenomenon. It will be used to study and produce the LSPR biosensors.Keyword：Localize Surface Plasmon Resonance，Discrete Dipole Approximation，Finite Difference Time Domain，Ag nanoparticles 精选学习资料 - - - - - - - - - 名
6、师归纳总结 - - - - - - -第 1 页，共 10 页1 引言贵金属纳米颗粒金、银等在紫外- 可见光区域存在很强的吸收光谱带，这是由于入射光子频率与金属传导电子的整体振荡发生共振，纳米颗粒或金属岛会对光子能量产生很强的吸收作用，这就是所谓的局域外表等离子体共振LSPR 现象。吸收峰在光谱上的位置与纳米颗粒的大小、形状、 间距及局域介电环境密切相关。这种现象在十八世纪法拉第研究胶体金时被发现1，目前基于LSPR现象的传感器被广泛研究，用于药物检测，生物检测、细胞标记、定点诊断、分子动力学研究及疾病诊断等方面。这种传感器具有无需标记、实时、无污染、高灵敏度检测、所需的样液少等优点。1908
7、 年，Mie 提出了解释LSPR现象的 Mie 理论2，Mie 理论简单实用， 能很好的解决球形纳米颗粒的消光包括吸收和散射问题， 但是不能解决复杂形状的颗粒以及颗粒之间的互相耦合作用及基底对纳米颗粒消光特性的影响问题。为此，各种数值计算理论和方法相继产生，如有限元法(FEM: Finite Element Method)3，严格耦合波法(RCWA: Rigorous Coupled-Wave Analysis)4，离散偶极近似法(DDA: Discrete Dipole Approximation)5和时域有限差分法(FDTD: Finite Difference Time Domain)6
8、等等。其中 DDA和 FDTD方法在分析LSPR现象中使用较多7，但多数只限于对某一种特例的分析；本文利用两种方法的优势互补，既分析了纳米颗粒的消光现象，也分析了纳米颗粒阵列的耦合作用。2 DDA与 FDTD 方法的比较DDA方法最初是由DeVoe8提出用于分析分子聚合体的光学特性，经过逐步的发展与完善而成为一种研究纳米光学现象的有力工具。到目前为止已经有大量的文献详细讨论了DDA的原理及其应用7,9,10。 Drain 等人在 DDA方法的基础上应用快速傅里叶变换和离散共轭梯度的方法制作了开源软件DDSCAT11。DDSCAT 可以将各种不同形状，不同材质的纳米颗粒离散成偶极子，计算其消光效
9、率因子(Q_ext) ，吸收效率因子(Q_abs) ，散射效率因子(Q_sca) ，其中Q_ext=Q_abs+Q_sca。不仅可以计算在单一频率下的纳米颗粒的消光情况，而且可以计算纳米颗粒在连续频谱下的消光现象。此软件不仅适用于单个纳米颗粒，对周期排列的纳米颗粒阵列也同样适用。缺点是对于不是软件自有形状的颗粒需要重新编程建模，建模的过程不是可视化的，需要对输出结果进行再处理才能描绘颗粒外表电场分布情况。FDTD方12，用变量离散的，含有有限个未知数的差分方程近似的替代连续变量的微分方程。 FDTD是将连续的空间划分为一个个的Yee 元胞，以Yee 元胞为空间电磁场离散单元，将麦克斯韦旋度方程
10、转化为差分方程，结合电脑技术解决电磁学问题。XFDTD6.3.8.4 是 Remcom 公司开发的基于FDTD的电磁学仿真软件，可以应用于分析贵金属纳米颗粒的LSPR现象13。 此软件可以很方便的进行可视化建模，有别于 DDSCAT 中采用的是真实的材料折射率，XFDTD需要用 Debye 模型或者 Lorentz模型来模拟金属在不同频率下的介电常数， 从而带来了误差。而且对于在连续光谱下的消光问题，需要将入射光源设置为高斯脉冲或者修正高斯脉冲来模拟一定波段光源的频谱，而结果通过傅里叶变换将时域的结精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2
11、 页，共 10 页果转化为频域的结果，优点是可以实时的仿真纳米颗粒与光作用的过程，可以很方便的输出纳米颗粒周围的电场分布图。本文将综合运用两种方法的优势，用DDSCAT分析各种纳米颗粒在不同波段的光谱，用XFDTD分析在具体光谱下的纳米颗粒外表电场分布。在分析不同大小的纳米颗粒的消光现象的基础上， 分析消光现象随两个纳米颗粒的间距以及入射光方向而产生的变化，以及纳米颗粒组成的阵列的消光现象。3 纳米颗粒的消光光谱及其电场分布利用 DDSCAT 软件，分析了半径R为 30nm的硅球， 二氧化硅球以及银球在可见光区域内的消光现象。 它们分别代表纳米半导体、绝缘体和金属与入射光之间的相互作用。各折射
12、率参数取自参考文献14 ，仿真结果如图1。300400500600700800-1012345678ExtinctionWavelength /nm Q_ext(Ag) Q_ext(Si) Q_ext(SiO2)图 1 不同材质的纳米球形颗粒的消光光谱Fig.1 The extinction spectra of different materials nanoparticles sphere 从图 1 可以发现二氧化硅的消光趋近于零，硅在波长391nm处存在很弱的的消光峰15，银球在波长371nm处产生很强的消光峰，此时发生了LSPR现象。我们用 XFDTD软件分析了半径为30nm的银和二氧
13、化硅两种球形颗粒的周围电场分布，银的 Modified Debye模型参数取自文献16 ，二氧化硅的折射率为1.47 介电常数，电导率为 10-121014Hz ，沿 x 方向偏振，振幅为1v/m，沿 +z 方向传播。 得到在 - 截面的外表电场分布如图2 所示。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 10 页(a) Silver (b) Silica 图 2 不同材质纳米颗粒的外表电场分布Fig.2 Surface electrical field distribution of different materials nan
14、oparticles 从图 2 中可以看出，在x-y 截面，银球外表的最大电场强度可到达23.99v/m ，且外表电场分布不均匀，说明有共振现象发生；而二氧化硅球的最大电场强度仅为1.546v/m ，且外表电场分布变化不大，说明没有明显的共振现象发生；通过消光光谱和外表电场分布的比照可以看出，只有金属，特别贵金属才有明显的LSPR现象发生。这是由于金属的原子核带正电， 而其周围的价电子带负电，在没有光照的作用下，可以把金属中的价电子看成是均匀正电荷背景中运动的电子气体，这种结构可以看作是一种等离子体，即金属外表等离子体。当特定波长的光入射到颗粒外表时，金属中电子密度分布就会变得不均匀。设想在某
15、一区域电子密度低于平均密度，这样便形成局部的正电荷过剩；这时由于库伦引力的作用，会把邻近的电子吸引到该领域；而由于被吸引的电子具有惯性，又会使该区域聚集过多的负电荷；然后， 由于电子间的排斥作用，使电子再度离开该区域，从而形成价电子相对于正电荷背景的密度起伏振荡，这就是LSPR现象。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 10 页4 半径对球形纳米颗粒消光的影响郭伟杰等人17通过 Mie 理论计算得到， 银纳米颗粒的消光光谱随着颗粒半径R的大小而变化，当半径小于4nm时，消光效率近似为0，当R在 4nm到 56nm之间时，消光效
16、率随波长变化曲线是单峰的形式，而当56nmR500nm,由于纳米颗粒半径的增大，高阶偶极子的作用，消光效率曲线呈现多个波峰的形式。我们用 DDA方法分析了半径R从 5nm到 50nm的银纳米颗粒在光波长300nm到 800nm的消光效率Q_ext 、吸收效率Q_abs 、散射效率Q_sca(Q_ext=Q_abs+Q_sca) 。仿真结果如图 3、图 4 所示。3003504004505005500.00.20.40.60.81.01.21.41.6IntensityWavelength /nmQ_ext Q_abs Q_scaR=5nm300350400450500550012345678I
17、ntensityWavelength /nmQ_ext Q_abs Q_scaR=30nmaR=5nm bR=30nm 图 3 不同大小的银纳米颗粒的消光效率、吸收效率、散射效率Fig.3 The Q_ext,Q_abs,Q_sca of different size silver nanoparticles 02040360390420Wavelength(nm)R(nm) Wavelength(nm) Parabola Fit of Wavelength(nm)01020304050012345678IntensityR/nm Q_ext Q_abs Q_sca01020304050-10
18、0102030405060708090100110Percent(%)R/nmQ_abs/Q_ext Q_sca/Q_exta Peak wavelength bQ_ext、Q_abs 、Q_sca cQ_abs/Q_ext、Q_sca/Q_ext 图 4 消光现象与银纳米颗粒半径之间的关系Fig.4 The relationship between silver extinction phenomenon and radius of nanoparticles 从图 3可以看出当R=5nm时，Q_abs占Q_ext的主导地位， 消光峰在 356nm处， 而当R=40nm时，Q_abs和Q_s
19、ca差不多，消光峰在371nm处，并且此时消光效率大于R=5nm时的消光效率。记录下不同大小的纳米颗粒的消光峰的位置如图4a ，分析在光谱波峰时的Q_ext、Q_abs 、Q_sca。消光峰位置与粒子半径R的拟合曲线为：消光峰波长nm RR2。在同为球形纳米颗粒时，当颗粒半径较小时，消光峰位置波长随颗粒半径变化较小，而当颗粒半径增大时，消光峰位置波长受颗粒半径的影响较大。从图 4b可以看出Q_ext、Q_abs随着粒子半径R的增大先增大后减小；在R=40nm时，Q_ext最大；R=24nm时，Q_abs最大。Q_sca则一直随着颗粒半径而增大。随着R增大，精选学习资料 - - - - - -
20、- - - 名师归纳总结 - - - - - - -第 5 页，共 10 页Q_abs/Q_ext越来越小，Q_sca/Q_ext越来越大18，如图 4c所示。同样，用DDA方法分析了R从 5nm到 50nm之间的银纳米颗粒分别在入射波长为200nm，400nm、600nm、800nm时的Q_abs/Q_ext，Q_sca/Q_ext的变化规律。01020304050-100102030405060708090100110Percent(%)R(nm)Q_abs/Q_ext（200） Q_sca/Q_ext(200nm) Q_abs/Q_ext(400nm) Q_sca/Q_ext(400nm
21、) Q_abs/Q_ext(600nm) Q_sca/Q_ext(600nm) Q_abs/Q_ext(800nm) Q_sca/Q_ext(800nm)图 5 不同大小的银纳米颗粒在波长不同的入射光时Q_abs/Q_ext ，Q_sca/Q_ext Fig.5 The Q_abs / Q_ext, Q_sca / Q_ext of different sizes of silver nanoparticles under different incident wacelength从图中可以看出， 不管是在消光峰处，还是在其它入射波长处，随着颗粒半径R的增大，Q_abs/Q_ext越来越小，Q
22、_sca/Q_ext越来越大；即颗粒越大，散射效率所占的总的消光效率的比例越大。5 纳米颗粒之间的耦合作用为了分析纳米颗粒之间的耦合作用，以半径为30nm的两个银球为例，分析它们相距d分别为 0nm 、5nm 、10nm 、30nm、60nm时的消光现象。分别取入射光平行于和垂直于两颗粒中心轴线两种情况；当两个颗粒相互作用时，其消光现象及外表电场与入射光的方向有关。300350400450500550012345678910ExtinctionWavelength /nm Q_ext(d=0) Q_ext(d=5) Q_ext(d=10) Q_ext(d=30) Q_ext(d=60)3003
23、504004505005506000123456789ExtinctionWavelength/ nm Q_ext(d=0) Q_ext(d=5) Q_ext(d=10) Q_ext(d=30) Q_ext(d=60)a The incident light parallel to the axis (bThe incident light perpendicular to the axis 图 6 消光光谱与纳米颗粒间距及入射光方向的关系Fig. 6 The relationship between the extinction spectra and the nanoparticles d
24、istance,the direction of incident light 可以看出， 当光入射方向平行于颗粒之间连线时，随着颗粒间距的增加，消光峰位置产生红移。而当光线垂直与颗粒之间连线时，随着颗粒间距的增加，消光峰位置产生蓝移，当两个颗粒距离为0nm时，由于高阶偶极子的作用，产生了多个消光峰。当两个颗粒的距离大于颗粒半径时，消光峰的位置与单个银颗粒的消光峰一致。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 10 页1014Hz ，则在两种入射光作用下x-z 截面上的电场分布分别于图7a, b所示。a The incident
25、light parallel to the axis (bThe incident light perpendicular to the axis 图 7 d=10nm时两个颗粒的外表电场分布Fig.7 The surface electric field distribution of two nanoparticles when d = 10nm 由图 7 说明， 当入射光平行于轴线时，最大电场强度为7.879v/m，最大电场强度集中于两个颗粒的两侧，而当入射光垂直于轴线时，最大电场强度为44.16v/m，而且最大电场强度集中于两个颗粒之间。外表电场的大小也与两个纳米颗粒之间的距离有关，当
26、d=0nm 时，外表最大电场强度可以到达129v/m，当两个颗粒距离变大时，最大外表电场强度也减小。这是因为当光线平行于轴线时，其电场方向与入射光方向垂直，即电场方向垂直于两颗粒连线方向， 引起颗粒共振的方向与轴线垂直，此时两个颗粒之间的耦合作用较小。而当入射光垂直于轴线时，入射光电场方向与颗粒轴线平行，引起共振的方向也于轴线平行，在两个颗粒之间由于偶极子的相互作用，产生很强的共振作用，导致其外表电场强度很大，纳米颗粒外表的等离子体是有一定的距离的，所以相隔的越近，作用越剧烈，而当距离相对较大时，无论入射光的方向如何，两个颗粒的耦合作用都较小。用 DDSCAT 分析 R=30nm的两个银纳米颗
27、粒在相隔10nm时在不同介质 真空，n=1；水 ，n=1. 33；丙酮，n=1. 36；二氯甲烷，n=1. 42；及嘧啶，n=1. 51，分析其折射率灵敏度。两个纳米颗粒的消光峰波长折射率灵敏度：当光垂直入射时，消光峰波长nm =+n；当光水平入射时，消光峰波长nm =+ n。即当光垂直入射时，两个纳米颗粒对于外界折射率的变化更敏感。6 银纳米球及其阵列的消光光谱与折射率灵敏度分析与纳米颗粒阵列相比， 单个纳米颗粒在检测应用方面也有很多优点，如检测的极限更高，需要样品的容量更少，可以运用到多通道检测等，不过制作费用也较高19；因此，目前一般采用纳米球光刻NSL：Nanosphere Litho
28、graphy 技术在硅片或玻璃片上制作纳米颗粒的周期排列20。用 DDSCAT 分析方法，对半径为5nm，任意两球间隔d 分别为0nm 及 5nm时的银纳米球的二维周期阵列在可见光范围内的消光光谱进行了分析，其结果如图8 所示。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页，共 10 页3004005006007008000.00.51.01.52.0IntensityWavelength /nmQ_ext Q_abs Q_sca3004005006007008000.00.20.40.60.81.01.21.41.6IntensityWa
29、velength /nm Q_ext Q_abs Q_sca(a)d=0nm (b)d=5nm 图 8 不同间距的银纳米球周期阵列的消光光谱Fig.8 The extinction spectra of different separation distance silver nanoparticles arrays 从上图可以看出，周期排列的银纳米球阵列出现两个消光峰，d=0 的两个消光峰分别在nm 和 nm 处，在 372.92nm 处的消光峰主要是由吸收引起的，nm 处的消光峰比第一个峰强度大很多， 主要是由散射引起的。这是由于当纳米颗粒之间的距离为0时，纳米颗粒的总体形状变得复杂，出现
30、多偶极子的共振，因而出现第二个消光峰。而d=5nm 的第二个消光峰明显比第一个消光峰弱，可以预测当周期阵列型的纳米颗粒间距相对较大时，其消光峰会与单个纳米颗粒的消光峰一致，即为单峰形式。为了分析银纳米球及其阵列的折射率灵敏度，假设待测介质环境分别为真空n=1、水 (n=1. 33) 、丙酮 (n=1. 36) 、二氯甲烷 (n=1. 42) 及嘧啶 (n=1. 51)等，分析比较它们的消光现象由图9 给出。200300400500600700800900100011000.00.51.01.52.0ExtinctionintensityWavelength /nm Q_ext(n=1) Q_e
31、xt(n=1.33) Q_ext(n=1.36) Q_ext(n=1.42) Q_ext(n=1.51)1.01.11.21.31.41.51.6360380400420440460 Peak wavelength 1 Peak wavelength 2Refractive indexPeakwavelength1600650700750800850900950Peak wavelength 2(a)The extinction spectra In different medium (b) The fitting curve between extinction peak position
32、and the external medium refractive index 图 9 银纳米球阵列在不同外部介质下的消光光谱和折射率灵敏度Fig. 9 The xtinction spectra and sensitivity of silver nanoparticles array under different external medium 由图 9(b) 可得两个消光峰的位置波长与外部介质的折射率拟合直线为：消光峰1 波长 nm=214.9+n；消光峰2 波长 nm=116.1+n。在消光峰1 的折射率灵敏度为157.0nm/RIU ， 消光峰 2 的折射率灵敏度为527.7nm
33、/RIU 。 而单个 5nm银球的消光峰波长 nm =+132.5n。第一个消光峰的折射率灵敏度与单个银纳米颗粒相近，而第二个消光峰由于纳米颗粒之间的耦合作用其灵敏度增大。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页，共 10 页7 结论金属纳米颗粒的LSPR 消光现象与纳米颗粒的种类、大小、形状及所处的介质环境有关。本文重点分析了贵金属纳米颗粒LSPR 现象与颗粒大小以及颗粒之间的耦合作用的影响。比照分析了金属球与绝缘体球及半导体球的消光光谱以及外表电场分布，验证了 LSPR 现象产生的原因。通过对不同大小的纳米颗粒的消光、吸收、散射光
34、谱的分析得知，当纳米颗粒半径较小时， 消光峰位置波长随颗粒半径变化不大；而当纳米颗粒半径较大时，消光峰位置波长随颗粒半径变化增大；同时随着颗粒半径的增大，散射效率占总的消光效率的比重增加。同样分析了两个纳米颗粒的耦合作用，发现当入射光垂直与轴线，且纳米颗粒之间间距较小时，耦合的作用较大。 同样纳米颗粒阵列的消光光谱的灵敏度与纳米颗粒的大小以及间距有关，当颗粒间距很小的时候会出现多峰形式的消光谱。这些规律为我们制作基于纳米颗粒LSPR 现象的生物传感器提供了有用的参考依据。参考文献：1 David Thompson. Michael Faraday s Recognition of Ruby G
35、old: the Birth of Modern Nanotechnology ReferencesJ. Gold Bulletin ,2007, 40(4):267-269. 2 Prashant K J, Ivan H E.An nanoparticles target cancerJ.Nanotoday,2002,2:19-29. 3 Jeffrey M, McMahon. Gold nanoparticle dimer plasmonics: finite element method calculations of the electromagnetic enhancement to
36、 surface-enhanced Raman spectroscopyJ. Anal Bioanal Chem,2009, 394:1819 1825. 4 Moharam M G,Gaylord T K. Rigorous coupled-wave analysis of planar-grating diffraction J. Opt. Soc. Am, 1981, 71(7): 811-818. 5 Draine B T. The discrete-dipole approximation andIts application to interstellar graphite gra
37、ins J.The Astrophysical Journal, 1988, 333(22):848872. 6 Yao H M, Li Z, Gong Q H. Coupling-induced excitation of a forbidden surface plasmon mode of a gold nanorod.J. Sci China Ser G-Phys Mech Astron, 2009, 52(8): 1129-1138. 7 Zhou Fei, Li Zhi-Yuan.Quantitative Analysis of Dipole and Quadrupole Exci
38、tation in the Surface Plasmon Resonance of Metal NanoparticlesJ. J. Phys. Chem. C,2008, 112: 20233 20240. 8 DeVoe H. Optical properties of molecular aggregates.I. Classical model of electronic absorption and refractionJ. J. Chem. Phys,1964,41, 393-400. 9 Yang W H,1995, 103:869-875. 10Kelly K L, Laza
39、rides A A,Schatz G C.Computational electromaganetics of metal nanoparticles and their aggregatesJ.Nnao technology, 2001, 12:67-73. 11 Draine B T , Flatau P J. “User Guide for the DiscreteDipole Approximation Code DDSCAT (Version 7.0) ”EB/OL. :/arxiv.org/abs/0809.0337v4, 2000.12Kane Yee. Numerical so
40、lution of initial boundary value problems involving Maxwells equations in isotropic mediaJ.IEEE Transactions on Antennas and Propagation,1966,14: 302 307. 13Zhang Hai-xi,Gu Ying, Gong Qi-huang.A visible-near infrared tunable waveguide based on plasmonic gold nanoshellJ. Chinese Phys. B,2008,17(7):25
41、67-2573 14Edward D Palik. Handbook of Optical Constants of Solids M.Washington DC:Academic press,1985.286-297,350-359. 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页，共 10 页15 Wilcoxon J P, Samara G A.Tailorable, visible light emission from silicon nanocrystalsJ. Applied Physics Letters,1999,74(2
42、1):3164-3166 16 Gai Hongfeng.Modi?ed Debye model parameters of metals applicable for broadband calculationsJ. Applied Optics,2007,46(12):2229:2233. 17 Guo WeiJie,Zhang Yuan. Study of radius effect extinction crosssection of silver-nanopart iclesJ. Journal of Xinjiang University,2008,25(2):187-189. 1
43、8 Eun C C,Chulhong K. Measuring the Optical Absorption Cross Sections of Au Ag Nanocages and Au Nanorods by Photoacoustic ImagingJ. J. Phys. Chem. C, 2009,113(21):9023-9028 19Stuart D A. Biological applications of localised surface plasmonic phenomenaeJ.IEE Proc Nanobiotechnol,2005,152(1):13-32. 20Christy L H, Richard P V. Nanosphere Lithography: A Versatile Nanofabrication Tool for Studies of Size-DependentNanoparticle OpticsJ. J. Phys. Chem. B, 2001, 105: 5599-5611. 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页，共 10 页
