(001)应变对正交相Ca2 P0.25 Si0.75能带结构及光学性质的影响.pdf

  • 文档大小:2.33 MB
  • 文档格式:pdf
  • 约 6页
  • 2021-11-23 发布
  • 举报
(001)应变对正交相Ca2 P0.25 Si0.75能带结构及光学性质的影响1 (001)应变对正交相Ca2 P0.25 Si0.75能带结构及光学性质的影响2 (001)应变对正交相Ca2 P0.25 Si0.75能带结构及光学性质的影响3 (001)应变对正交相Ca2 P0.25 Si0.75能带结构及光学性质的影响4 (001)应变对正交相Ca2 P0.25 Si0.75能带结构及光学性质的影响5 (001)应变对正交相Ca2 P0.25 Si0.75能带结构及光学性质的影响6
已阅读完毕,您还可以下载文档进行保存
文档单价:2.00 会员免费
开通会员可免费下载任意文档
  1. 1、本文档共6页,内容下载后可编辑。
  2. 2、本文档内容版权归属内容提供方,所产生的收益全部归内容提供方所有。如果您对本文有版权争议,可选择认领。
  3. 3、本文档由用户上传,本站不保证质量和数量令人满意,可能有诸多瑕疵,付费之前,请仔细先通过免费阅读内容等途径辨别内容交易风险。如存在严重挂羊头卖狗肉之情形,可联系本站下载客服投诉处理。
第23卷第3期2015年6月材料科学与工艺MATERIALSSCIENCE&TECHNOLOGYllVol23llNo3Jun.2015doi:10.11951/j.issn.1005-0299.20150322(001)应变对正交相Ca2P0.25Si0.75能带结构及光学性质的影响岑伟富,杨吟野,范梦慧,杨文帮,姚娟(贵州民族大学理学院,贵阳550025)摘要:为了研究(001)应变对正交相Ca2P0.25Si0.75能带结构及光学性质的影响,采用第一性原理贋势平面波方法对(001)应变下正交相Ca2P0.25Si0.75的能带结构及光学性质进行了模拟计算.计算结果表明:晶格(001)面发生100%~116%张应变时,带隙随着应变增加而减小;在晶格发生88%~100%压应变时,带隙随着张应变的增加而增加;84%~88%压应变时,带隙随着压应变的增加而减小.当施加应变后光学性质发生显著的变化:随着压应变的增加,静态介电常数、折射率逐渐减小,张应变则增大.施加压应变反射向高能方向偏移,施加张应变反射向低能方向偏移,但施加应变对反射区域的影响不显著.施压应变吸收谱、光电导率的变化与介电函数和折射率相反.综上所述,(001)应变改变了Ca2P0.25Si0.75的电子结构和光学常数,是调节Ca2P0.25Si0.75光电传输性能的有效手段.关键词:应变;能带结构;光学性质;第一性原理;Ca2P0.25Si0.75中图分类号:0472.+3文献标志码:A文章编号:1005-0299(2015)03-0112-06Theeffectof(001)strainonenergybandstructureandopticalpropertiesofthesimpleorthorhombicCa2P0.25Si0.75CENWeifu,YANGYinye,FANMenghui,YANGWenbang,YAOJuan(CollegeofScience,GuizhouMinzuUniversity,Guiyang550025)Abstract:Tostudytheeffectof(001)strainonenergybandstructureandopticalpropertiesofthesimpleorthorhombicCa2P0.25Si0.75.TheenergybandstructureandopticalpropertiesofthesimpleorthorhombicCa2P0.25Si0.75⁃⁃bulkunderthe(001)strainhavebeencalculatedbythefirstprinciplepseudopotentialmethodbasedondensityfunctionaltheory(DFT).TheresultsshowthatthesimpleorthorhombicCa2P0.25Si0.75bulkisadirectsemiconductorandthebandgapvaluedecreasedwiththetensilestrainincreasedinrangeof100%~116%andthecompressivestrainincreasedinrangeof84%~88%,butitisincreasedwiththecompressivestrainincreasedinrangeof88%~100%.Theopticalpropertieshaveasignificantlychangebystraining.ThedielectricconstantandtherefractiveindexofthesimpleorthorhombicCa2P0.25Si0.75bulkaredecreasedwhencompressivestrainedincreased,buttheyareincreasedastensilestrainedincreased.ThereflectivityofthesimpleorthorhombicCa2P0.25Si0.75bulkmovedtothehighenergydirectionwiththecompressivestrainincreased,anditmovedtothedirectionoflowenergywiththetensilestrainincreased.TheabsorptionandtheconductivityofthesimpleorthorhombicCa2P0.25Si0.75bulkareincreasedwithcompressivestrainedincreased.ThetransmissionpropertyoftheenergybandstructureandopticalpropertiesofthesimpleorthorhombicCa2P0.25Si0.75bulkarechangedbystrainingonthe(001)surface.ItisausefulmethodformodulatingthephotoelectrictransmissionpropertyofthesimpleorthorhombicCa2P0.25Si0.75bulk.Keywords:straining;energybandstructure;optical⁃properties;firstprinciple;Ca2P0.25Si0.75收稿日期:2014-03-11.基金项目:教育部科学技术研究重点项目(No.210200);贵州省科学技术联合基金项目(No.LKM201130);贵州省优秀科技教育人才省长专项资金项目(No.201174)资助项目.作者简介:岑伟富—(1988),男,硕士研究生.通信作者:杨吟野⁃,Email:Ythin1969@sohu.com.碱土金属硅化物Ca-Si是一种新型的环境友好型半导体材料,Ca和Si反应生成六个中间相,其中Ca2Si是直接带隙半导体.Ca2Si具有两种结构,一种属于正交晶系,为稳定结构[1],禁带宽度为0.36eV,与实验值1.9eV相差较大[2-3];另一种结构属于立方晶系,在外压力作用下为稳定结构,禁带宽度为0.56eV[4].Ca2Si具有优异的光学、电学性质,同时能在硅基上外延生长,在制备和使用过程中对环境无污染长期使用对生命体无害,因此在光电子器件、电子器件及热电器件等领域有广泛的应用前景,是一种潜在的理想半导体材料.Imai等利用第一性原理贋势方法⁃(Thefirstprinciplepseudopotentialmethod)计算得出正交相Ca2Si的能带结构,带隙为0.36eV的直接带隙半导体[2,5].Migas等用全电势线性化缀加平面波方法(FLAPW)和广义梯度近似(GGA)对Ca2Si的基态能带、态密度和介电函数进行计算,得到正交相是一种稳定相.在实验方面,基于电阻率测量实验报道Ca2Si的带隙为1.90eV[6].TakagiN等人应用蒸发加热处理两步法首次制备了单一相的Ca2Si晶体[7-8].杨吟野等利用磁控溅射⁃(MS⁃Magnetronsputtering)技术,在Si(100)衬底上沉积Ca膜,然后在真空退火中获得单一相直接迁移型的Ca-Si膜[9-10].Ca-Si膜在实际生长过程中受到温度、压强等使用条件影响导致其晶格发生应变,从而影响Ca-Si膜的电学性质和光学性质.到目前为此应变对材料性能的影响理论研究的报道较少[11-18],但应变对Ca2P0.25Si0.75电子结构及光学性质的影响研究尚未见报道.Ca2P0.25Si0.75作为新型的半导体材料,有必要深入研究应变对其电子结构及光学性质的影响,揭示应变对其电子结构和光学性质的影响机理.因此,本文采用基于密度泛函理论的第一性原理平面波方法对应变作用下Ca2P0.25Si0.75的电子结构及光学性质进行全面的模拟计算和分析.1计算方法正交晶系Ca2Si空间群为Pnma(No.62),晶格常数为a=0.7667nm,b=0.4779nm,c=0.9002nm每个原胞有12个原子,其中8个Ca原子,4个Si原子[1].由于正交相Ca2Si模型的4个Si原子处于等效位置,P原子的置换不同位置对体系的影响较小,计算采用的正交相相Ca2P0.25Si0.75模型暂不考虑不同位置置换对体系的影响.本文采用基于第一性原理的贋势平面波方法对晶格(001)应变的模拟计算,即将正交相Ca2P0.25Si0.75晶体的晶格常数c在其晶胞优化后所得的晶格平衡值的84%~116%范围内进行线性变化,文中所有的计算由CASTEP(Cambridgeserialtotalenergypackage)软件包完成[19].首先将Ca2P0.25Si0.75的体系采用⁃BFGS(Broyden,Fletcher,GoldfarbandShannon)算法进行几何优化,得到稳定的结构体系.在稳定的结构体系下分别计算不同应变的电子结构和光学性质,体系的价电子波函数用平面波基矢展开并设平面波的截断能为330eV,迭代收敛精度为5×10-6eV,选取广义梯度近似处理交换关联能部分,交换关联势计算采用PBE(PerdewBurkerErnzerhof)提出的广义梯度近似方法(GGA),采用超软贋势(UlterSoftPseudoPotential)计算离子实与电子之间的相互作用,计算总能量在倒易空间中进行,布里渊区积分采用⁃MonkhorstPack方法K点取4×6×3.2结果及讨论2.1能带结构能带结构是分析光电子半导体材料的重要依据,在结构优化的基础上利用广义密度梯度近似(GeneralizedGradientApproximation,GGA)近似处理交换关联泛函,超软贋势处理离子实与价电子之间的相互作用,平面波基组描述体系电子的波函数,计算得到应变条件下Ca2P0.25Si0.75沿布里渊区高对称点方向的能带结构,如图1所示.图1(d)为平衡状态下(未发生应变)的能带结构,由图可知Ca2P0.25Si0.75是带隙值为0.56760eV的直接带隙半导体,价带主要由Si的3p,Ca的3d、4s以及P的3p态电子共同贡献,导带主要由Ca的3d电子贡献.当Ca2P0.25Si0.75材料受到(100)面应力应变时,晶体中原子的相对位置发生变化迫使原子核附近的电子云发生重组,因此晶格常数相应发生变化使得其能带结构发生变化.图1(a)~(c)为施加压应变Ca2P0.25Si0.75的能带结构图,施加压应变改变了晶体中原子的相对位置,从而影响电荷的转移、Ca-Si、Ca-P共价键的键长和键角.随着压应变的增加Ca-Si、Ca-P共价键的键长被压缩,共价性增强,从而导致带隙变宽,这是施加压应变在88%~100%能带变化的主要原因.当应变达到88%时其转变为n型半导体,在84%~88%的压应变范围,带隙逐渐减小,其主要原因是随着应变得增大导致价带分裂.图1(e)~(f)为施加张应变Ca2P0.25Si0.75的能带结构图,随着张应变的增加带隙逐渐变小,当张应变为116%时Ca2P0.25Si0.75转变为金属.计算结果表明了无论施加张应变还是压应变均能改变Ca2P0.25Si0.75的能带结构,验证了应变是调控能带结构的有效手段.·311·第3期岑伟富,等:(001)应变对正交相Ca2P0.25Si0.75能带结构及光学性质的影响3210-1-2-3-4ΓZTYSXUR能带量/eV(a)3210-1-2-3-4ΓZTYSXUR能带量/eV(b)3210-1-2-3-4ΓZTYSXUR能带量/eV(c)3210-1-2-3ΓZTYSXUR能带量/eV(d)3210-1-2-3ΓZTYSXUR能带量/eV(e)3210-1-2-3ΓZTYSXUR能带量/eV(f)图1应变为84%(a)、88%(b)、96%(c)、100%(d)、104%(e)、116%(f)的正交相Ca2P0.25Si0.75的能带结构2.2光学性质2.2.1介电函数介电函数作为沟通带间跃迁微观物理过程与固体电子结构的桥梁,通过介电函数可以得到其他各种光谱信息,Ca2P0.25Si0.75作为新型半导体材料,其光谱是由能级间电子跃迁产生.图2为应变条件下Ca2P0.25Si0.75的介电函数图.3020100-1001234530252015105012345ε1ε284%88%96%100%104%116%84%88%96%100%104%116%E/eVE/eV(a)(b)图2应变为84%、88%、96%、100%、104%、116%的正交相Ca2P0.25Si0.75的介电函数实部(a)和虚部(b)由图2可知,平衡状态下Ca2P0.25Si0.75的静态介电常数为13.81014.当施加应变后静态介电常数发生显著的变化:随着压应变的增加Ca2P0.25Si0.75的静态介电常数逐渐减小;随着张应变的增加Ca2P0.25Si0.75的静态介电常数逐渐增大.应变为84%、88%、96%、104%、116%对应的静态介电常数分别为13.12011、12.97338、13.22993、14.57501、16.18041.介电函数发生变化的主要原因是:在外力作用下晶体中原子发生重排导致晶体场发生变化,导致晶体中原子核周围的电子云发生畸变,使得正负离子中心发生相对位移极化.当施加压应变Ca-Si、Ca-P共价键的键长被压缩,共价性增强,极化削弱,介电常数变小;当施加张应变Ca-Si、Ca-P共价键的键长被拉伸,共价性削弱,极化增强,介电常数增大.介电函数虚部的峰值是由于电子跃迁产生,它反应了Ca2P0.25Si0.75能带结构及其他光谱信息,与未施加应变相比明显的发现:施加张应变,介电函数虚部的峰值及峰数随着·411·材料科学与工艺第23卷张应变增加而增大;施加压应变,介电函数虚部的峰值及峰数随着压应变增加而减小.说明施加张应变电子跃迁增强,施加压应变电子跃迁削弱,与计算得到的能带结构变化一致.2.2.2复折射率由复折射率和介电函数的关系ε1=n2-k2,ε2=2nk得到Ca2P0.25Si0.75的复折射率,图3为Ca2P0.25Si0.75的折射率和消光系数.54321012345η84%88%96%100%104%116%E/eV(a)24681043210κ84%88%96%100%104%116%E/eV(b)图3应变为84%、88%、96%、100%、104%、116%的正交相Ca2P0.25Si0.75的折射率(a)和消光系数(b)由图3可知,平衡状态下Ca2P0.25Si0.75的折射率为3.7162.当施加应变后折射率发生显著的变化:随着压应变的增加Ca2P0.25Si0.75的折射率逐渐减小;随着张应变的增加Ca2P0.25Si0.75的折射率增大.应变为84%、88%、96%、100%、104%、116%对应的折射率分别为3.6221、3.6017、3.6373、3.7162、3.8177、4.0225.当光子能量在约5.769~6.843eV的能量范围,折射趋于零,这与反射率在该能量范围趋于1相对应,表明在此能量范围Ca2P0.25Si0.75呈现很强的反射特性.随着应变从84%变化到116%,折射率的变化规律与介电函数实部的变化规律相似,均是由于原子相对位置变化所致.计算所得折射率n0的值与静态介电常数ε1(0)完全对应.由消光系数图可知,其变化曲线在1eV附近漂移,漂移的能量大小与带隙变化的数值相等.当施加压应变时消光系数k的峰值随着压应变的增加而减小,并且最大峰随着压应变向高能方向偏移;当施加张应变时消光系数k的峰值随着张应变的增加而增加,并且最大峰随着张应变的增加向低能方向偏移.结果表明在此能量范围内施加压应变其光学响应增强,施加张应变其光学响应削弱,同时消光系数在带边表现出强烈的红外吸收特性,并且随着压应变增加向高能方向偏移,随着张应变增加向低能方向偏移.其原因是应变使得Ca2P0.25Si0.75晶体中原子核周围的电子云发生畸变,能带结构和费米面附近的电子态密度分布发生变化所致,对电子跃迁产生影响.2.2.3反射率由复折射率和反射谱的关系R(ω)=(n-1)2+k2(n+1)2+k2得到反射谱.由图4可知,没有施加应变Ca2P0.25Si0.75的反射主要发生在能量为5.260~6.915eV的能量范围;施加84%、88%、96%、104%、116%的应变,反射分别发生在5.966~7.069eV、5.768~7.047eV、5.525~6.981eV、5.282~6.871eV、4.576~6.716eV的能量范围,反射率达到90%以上甚至反射全反射现象.由此可见施加压应变反射向高能方向偏移,施加张应变反射向低能方向偏移,但施加应变对反射区域的影响不显著,计算结果与折射率和吸收系数完全对应.1.00.80.60.40.20246810E/eV反射率84%88%96%100%104%116%图4应变为84%、88%、96%、100%、104%、116%的正交相Ca2P0.25Si0.75的反射谱2.2.4吸收谱由吸收系数和介电函数的关系a(ω)=ω/ncε2(ω)可以得到Ca2P0.25Si0.75的吸收谱如图5所示.无应变时Ca2P0.25Si0.75的吸收系数α(ω)在3.694eV处达到最大峰,最大峰值为160777.57cm-1.吸收系数α(ω)在6.076eV~6.451eV的能量范围吸收为零,与文中反射图(图4)中发生全反射的能量范围相对应.当施加84%、88%、96%、104%、116%的应变,吸收系数α(ω)分别在4.135、4.422、3.871、3.429、2.856eV处达到最大峰,最大峰值分别为171917.76、165056.98、164976.78、·511·第3期岑伟富,等:(001)应变对正交相Ca2P0.25Si0.75能带结构及光学性质的影响156741.98、160663.65cm-1;而分别在6.518~6.694eV、6.385~6.694eV、6.186~6.517eV、5.966~6.407eV、5.393~6.385eV的能量范围系数为零.1.5?1051.0?1055.0?1040246810E/eV最大吸收峰84%88%96%100%104%116%图5应变为84%、88%、96%、100%、104%、116%的正交相Ca2P0.25Si0.75的吸收光谱可见随着压应变的增加最大吸收峰随着增加,吸收系数为零的能量范围在减小,而随着张应变的增加最大吸收峰逐渐减小,吸收系数为零的能量范围逐渐增大.由此可得,施加压应变有利于Ca2P0.25Si0.75对光谱的吸收,而施加张应变则削弱Ca2P0.25Si0.75对光谱的吸收.2.2.5光电导率光电导率描述的是半导体材料电导率随着光照能量的变化而变化的物理参量.图6为Ca2P0.25Si0.75的光电导率实部δ1.从图6可以得出,分别施加84%、88%、96%、100%、104%、116%的应变,δ1在能量分别为2.371、2.017、1.731、1.687、2.392、2.084eV出达到最大值,最大值分别为6.20604、6.48775、6.30715、6.08545、6.48181、7.1172.可见施加压应变可使Ca2P0.25Si0.75的光电导率δ1增加并且向高能方向偏移;施加张应变可使Ca2P0.25Si0.75的光电导率增加但向低能方向偏移,与介电函数虚部和吸收系数进行比较发现,它们的变化趋势基本一致,验证了它们之间的对应关系.电导率87654321012345E/eV84%88%96%100%104%116%图6应变为84%、88%、96%、100%、104%、116%的正交相Ca2P0.25Si0.75的光电导率3结论采用基于密度泛函理论的第一性原理平面波超软贋势方法,模拟研究了应变作用下Ca2P0.25Si0.75(100)面发生晶格应变时的能带结构、介电函数、复折射率、吸收谱、反射谱、光电导率等光学性质.计算得到:Ca2P0.25Si0.75(001)面在晶格发生100%~116%张应变时,Ca2P0.25Si0.75仍然直接带隙半导体,但带隙随着应变增加而减小;Ca2P0.25Si0.75(001)面在晶格发生88%~100%压应变时,Ca2P0.25Si0.75属于直接带隙且禁带宽带随着张应变的增加而增加;84%~88%压应变时,Ca2P0.25Si0.75禁带宽带随着张应变的增加而减小.当施加应变后静态介电常数发生显著的变化:随着压应变的增加Ca2P0.25Si0.75的静态介电常数逐渐减小;随着张应变的增加Ca2P0.25Si0.75的静态介电常数逐渐增大.当施加应变后折射率发生显著的变化:随着压应变的增加Ca2P0.25Si0.75的折射率逐渐减小;随着张应变的增加Ca2P0.25Si0.75的折射率增大.施加压应变反射向高能方向偏移,施加张应变反射向低能方向偏移,但施加应变对反射区域的影响不显著.随着压应变的增加最大吸收峰随着增加,吸收系数为零的能量范围在减小,而随着张应变的增加最大吸收峰逐渐减小,吸收系数为零的能量范围逐渐增大.施加应变可使光电导率增加并且向高能方向偏移.综上所述,应变可以改变Ca2P0.25Si0.75的电子结构和光学常数,是调节Ca2P0.25Si0.75光电传输性能的有效手段.参考文献:[1]LEBÈGUES,ARNAUDB,ALOUANIM.CalculatedquasiparticleandopticalpropertiesoforthorhombicandCa2Si[J].PhysicalReviewB,2005,72(8):1-8.[2]IMAIY,WATANABEA.Energeticsof⁃alkalineearthmetal⁃silicidescalculatedusingafirstprinciplepseudopotentialmethod[J].Intermetallics,2002,10(4):333-341.[3]MATSUIH,KURAMOTOM,ONOT,etal.GrowthofCa2SilayersonMg2Si/Si(111)substrates[J].CrystGrowth,2002,237-239:2121-2124.[4]CHENGJ,YANGYY,LIANGY.SelectivegrowthofCa2SifilminCa-SisystembyR.FMSbyannealing,[C]3M-NANO,2012,9:75-79.[5]IMAIY,WATANABEA,MUKAIDAM.Electronicstructuresof⁃semiconductingalkalineearthmetal⁃silicides[J].JournalofAlloysandCompounds,2003,358(1-2):257-263.·611·材料科学与工艺第23卷[6]MIGASDB,MIGLIOL,SHAPOSHNIKOVVL,etal.Comparativestudyof⁃structural,electronicandopticalpropertiesofCa2Si,Ca2Ge,Ca2SnandCa2Pb[J].Phys.Rev.B,2003,67(20):205203(1-7).[7]TAKAGIN,SATOY,MATSUYAMATetal.GrowthandstructuralpropertiesofMg2SiandCa2⁃Sibulkcrystals[J].AppliedSurfaceScience.2005,244(1-4):330-333.[8]BUSCHC,JUNODP,KATZU,etal.Electrical⁃conductivityoftheintermetalliccompoundsCa2Si,Ca2Sn,Ca2Pb,ZnSb[J].Helv.Phys.Acta.1954,27(33):193-197.[9]YANGYY,XIEQ.Asinglephasesemiconducting⁃⁃Casilicidefilmgrowthbysputteringconditions,annealingtemperatureandannealingtim[J].JournalofMaterialsScience,2009,44(14):3877-3882.[10]YANGYY,XIEQ.SelectivegrowthofCa2SifilmbyannealingfromthedifferentthicknessofCafilmsbysputtering[J].材料科学与工程学报,2009,5(27):675-678.[11]马建立,张鹤鸣,宋健军,等.(001)面任意方向单轴应变硅材料能带结构[J].物理学报,2011,60(2):027101-027106.[12]王冠宇,马建立,张鹤鸣,等.[110]/(0001)单轴应变Si本征载流子浓度模型[J].物理学报,2011,60(7):077105-077111.[13]郭伟国,左红星,孟卫华,等.第三种应变时效与机械波普关联性探讨[J].材料科学与工艺,2012,20(1):128-134.[14]崔冬萌,贾锐,谢泉,等.Ru2Si3在应力作用下的第一性原理研究[J].发光学报,2011,32(9):907-912.[15]WANGLW,YUYW,WEIZJ.MicrostructureoftheHypoeutecticAl-MgAlloySolidifiedat4Gpahighpressure[J].JournalofHarbinInstintuteTechnology(NewSeries),2013,20(5):101-105.[16]胡智向,吴玉喜,顾书林,等.外压调制下ZnO晶体结构与光学性质变化特性的研究[J].量子电子学报,材料科学与工艺,2010,27(5):613-619.[17]宫莹,宗影影,候智敏,等.热变形参数对TC21钛合金微观组织的影响[J].2013,21(5):113-116.[18]闫万珺,张春红,桂放,等.应力作用下β-FeSi2的电子结构及光学性质[J].光学学报,2013,33(7):0716001-0716008.[19]SEGALLMD,LINDANPJD,PROBERTMJ,etal.⁃Firstprinciplessimulation:ideas,illustrationsandtheCASTEPcode[J].JournalofPhysics,2002,14(11):2717-2744.(编辑张积宾)·711·第3期岑伟富,等:(001)应变对正交相Ca2P0.25Si0.75能带结构及光学性质的影响
分享达人
该用户很懒,什么也没介绍
文档单价:2.00 会员免费
开通会员可免费下载任意文档