碳纳米管组成二维光子晶体的有效介电常数FDTD数值模拟.pdfVIP

维普资讯 计 算 物 理 第 23卷 第 3期 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Vo1．23．No．3 2OO6年 5月 Mav．2006 [文章编号】 1001．246X(2006)03·0375·04 碳纳米管组成二维光子晶体的有效介电 常数 FDTD数值模拟 汤炳书 ， 沈廷根 (1．连云港师范高等专科学校物理系，江苏 连云港 222006；2．江苏大学物理系．江苏 镇江 212003) [摘 要 】 把有序排列的碳纳米管看做二维光子 晶体 ，引入平均场方法 ，先计算单根纳米碳管的有效介 电常 数 ，再用时域有限差分方法 (FDTD)计算碳管阵列的有效介 电常数 ，其结果与 Maxwel1．Gamett(MG)法结果相 比，更接近实验数据 ． [关键词 】 纳米碳管；二维光子晶体 ；时域有限差分法 ；有效介电常数 ；数值模拟 [中图分类号】 Q441．2 [文献标识码】 A O 引言 纳米碳管是 Iijima等人在研究电弧放 电时，首先在碳电极表面发现的…，它具有奇特的光学、电学、力学 性质，依据空间螺旋结构的不同，可以表现为金属，也可表现为半导体，是微观应用领域的理想材料 ．光子 晶体是一种由不同介 电常数的材料间隔周期性排列的人工材料 ，当电磁波入射其上时，在某一频率范 围内 电磁波不能透射过去 ，这很象固体物理学 中晶体 的能带 J．由碳管组成的阵列有效介 电常数实验结果见文 [5]．用 Maxwell—Gamett(MG)近似法可 以计算复合结构的介 电常数 ，Garcia-Vidal等人用该法模拟 了实验结 果 ，发现在长波段 即低频段与实验符合较好 ，因为它忽略了复合结构 中不同物体之间的电场相互作用 ．可 用平面波展开法 ，转移矩阵法 ，时域有限差分法 (Finite—Difference．Time—Domain) 等模拟实验数据 ．时域有 限差分算法(FDTD)是在时域分析电磁波传输 问题的一种有效方法，它直接基于麦克斯韦方程组，具有广泛 的适用性与可操作性 ．本文采用不同近似 ，先计算单根碳管的有效介电常数 e ，然后用 TDFD方法计算碳管 组成二维光子晶体的有效介 电常数 ，计算结果与实验数据符合较好 ．限于篇幅，本文仅计算 TE模式入射的 情况 (即电磁波的电场分量垂直于管轴的方向即横电波)． 1 基本理论与计算公式 1．1 单根碳管的有效介电常数 单根碳管的介电常数在柱坐标系中可以用张量表示为 e(；，，三)：e，， +￡ ( +郑)， (1) ；，，；分别是圆柱坐标的单位矢 ，图1是碳管的截面图．R是外半径 ，r是 内半径 ．e，，(e．)分别是与石墨平面平行的(垂直的)石墨介电常数张量的 2个主要分量 ．在准静态近似下，假定 电场分量沿着管轴方 向，外电势形 式为 (r，声)=V㈩e一，导致单位长度的极化率 一47reo~eR2m甚 专窘 …………图，1……碳……管…截……面………图…一 P= rlR， m =m(上 ／∥)m， =m ／m． Fig．1 Cross-sectionofacarbonnanotube r，R分别是管子的内外半径，e。与 e分别是管内外的介 电常数 ．m是分离变量 中与角度相关的整数 ．引入 [收稿 日期】2OO4—11—09；[修 回日期 】2005—07—11 [基金项 目】江苏省 自然科学基金资助 (BK21X)4059) [作者简介]汤炳书(1963一)，男，江苏丹阳，副教授 ，主要从事计算凝聚态物理研究