企业官网建设流程全解析

做一家新闻媒体网站多少钱,网站没有备案可以做seo优化吗,设计笔记页面,邯郸网络作家村铁电材料#xff0c;如BaTiO3#xff0c;PZT和硅掺杂二氧化铪HSO#xff0c;凭借其独特的铁电性、压电性、热释电性以及高介电常数等特性#xff0c;在信息存储领域扮演着关键角色。铁电材料因其固有的非易失性存储能力和高速读写特性#xff0c;成为新一代高性能、低功耗…铁电材料如BaTiO3PZT和硅掺杂二氧化铪HSO凭借其独特的铁电性、压电性、热释电性以及高介电常数等特性在信息存储领域扮演着关键角色。铁电材料因其固有的非易失性存储能力和高速读写特性成为新一代高性能、低功耗存储器如FeRAM的核心基石为物联网设备与边缘计算提供关键的存储支撑。本文将以铁电晶体管为例从其铁电极化原理、铁电晶体管的结构和原理和仿真复现三个方面来进行说明。一、铁电极化成因在说明铁电极化之前先简单说明一下热释电材料的相关内容。热释电材料的特征是其内部存在非零的自发极化强度矢量。这种自发极化是材料晶体结构本身固有的属性源于晶格中正负电荷中心不重合。热释电性产生的必要条件是材料具有唯一确定的、不可被外部电场翻转的自发极化方向即在空间上只有一个可能的极化轴。这类材料的极化通常是“冻结”或难以用电场改变的如GaN材料。铁电材料是热释电材料的一个特殊类别它也具有自发极化但是区别于普通热释电材料的根本特性在于其自发极化矢量具有两个或两个以上的稳定取向方向并且这些不同取向的自发极化可以在施加足够强的外部电场下进行可逆的翻转。图1清晰地展示了铁电材料的标志性特征——电滞回线。横轴为电场强度单位是V/m纵轴为极化强度单位为C/cm2。对于该图所描述的过程归纳为(1) 在初始零电场状态下材料内部虽已存在自发极化的电偶极子但因其方向杂乱无章宏观上对外表现的净极化强度P为零(2) 随着外加电场E的增强这些电偶极子逐渐沿电场方向定向排列(3) 当电场达到足够强度时极化强度达到其最大值即饱和极化强度Ps(4) 若将外加电场强度逐渐减小至零极化强度并不会随之归零而是稳定地保留在一个非零值称之为剩余极化强度±Pr这意味着材料在无外场情况下仍能“记忆”其极化状态(5) 当施加一个足够强的反向电场并超过特定的阈值——矫顽场强度±Ec时材料的极化方向才发生完全翻转。电滞回线正是铁电材料独特极化行为最直观的体现。图1 铁电材料电滞回线[1]2011年德国科学家Böscke及其团队在研究一种特定的薄膜材料——硅掺杂二氧化铪Si:HfO2简称HSO时意外开启了一扇重要的大门。他们观察到经过硅掺杂的HSO薄膜竟然展现出了可被外部电场精确调控的极化特性更关键的是其关键的铁电性能参数——剩余极化强度Pr和矫顽场强度Ec会随着硅掺杂浓度的变化而显著改变。深入的结构分析揭示了这一现象的物理根源在特定的掺杂条件下HSO材料中原子的排列发生了变化。对比没有铁电性的普通HfO2薄膜具有铁电性的HSO薄膜中晶体结构发生了不对称的改变HSO薄膜材料在由单斜相向四方相转变的过程中产生了新的非中心对称的正交相导致了材料内部正负电荷中心的分离从而形成了自发极化并且这种极化状态能够在外电场作用下进行翻转也就是产生了铁电性。至此具有良好硅基CMOS工艺兼容性的氧化铪基铁电材料正式步入了历史舞台。图2 非中心对称铁电性正交相的形成[1]二、铁电晶体管的结构和原理典型的铁电晶体管的结构如图3和4所示。工作原理描述如下如图3当栅极金属施加正向电压使铁电层的电场超过矫顽场 Ec铁电材料的极化状态就会发生改变铁电材料层的极化强度P方向向下撤除栅压后铁电材料的仍然存在向下的剩余极化强度 Pr这会在沟道中吸引电子形成反型层使得晶体管器件在没有栅压的时候也处于开启状态如果要使器件关断则需要继续施加负栅压在沟道中感生空穴把铁电材料极化效应感生出的电子耗尽这就导致了器件的开启电压Vth降低反之如图4当栅极金属继续施加反向电压使铁电层的电场超过矫顽场 -Ec铁电层极化方向翻转为向上撤除栅压后铁电材料的仍然存在向上的剩余极化强度 -Pr这会在沟道吸引空穴使得晶体管器件在没有栅压的时候源漏电流更低继续施加正栅压时要沟道中形成反型层所需的电荷变多使开启电压Vth变得更高。继续增高栅压使铁电材料的电场超过矫顽场 Ec器件回到图3所示的状态如此实现了图3和图4的状态的循环这个过程的源漏电流和栅极电压的曲线就是电滞曲线。图3 对铁电晶体管施加正栅压后再撤出铁电材料处于正极化状态图4 对铁电晶体管施加负栅压后再撤出铁电材料处于负极化状态三、仿真复现使用Silvaco TCAD对铁电晶体管进行基础特性的仿真。使用Victory Device创建一个典型的MOSFET结构代码如图5所示。7~15行代码定义了仿真区域的Line网格19~20行指定了衬底和氧化层的矩形区域其中在氧化层的定义的末尾增加了“lkferro”表明氧化层具有铁电性但根据代码的编写情况此处也可以不指定随后会进行说明22~25行指定了电极的位置虽然也使用了BOX形状的定义方式但是y方向的厚度为0所以电极是线型的27行指定了衬底的p型均匀掺杂浓度29~30行指定了Drain和Source区域的n型高斯掺杂32行保存了结构文件其结果如图6所示。通过这部分的代码定义了基本的MOSFET结构。此时器件还不算是铁电晶体管因为还没有增加铁电的模型和参数。因为是用Line网格生成的结构所以网格也比较粗糙高斯分布的边缘不是很顺滑因此需要对网格进行细化。图5 典型的MOSFET结构的结构创建代码图6 上述仿真代码生成的器件结构如图7通过victorymesh对上述生成的结构进行网格细化。37行载入在上一个Module中生成的结构文件40行声明了网格重构类型为Delaunay三角网格43行提升所有区域的尖角/曲面的采样密度46行设置整个网格区域的最大尺寸为0.2um49~50行指定硅衬底和氧化层区域的网格最大尺寸为0.1um53~57对硅和二氧化硅的界面进行网格细化max.interface.size规定界面处的网格尺寸为1nmmax.interface.distance规定网格逐渐变宽的范围是距离界面处的120nm内 max.interface.distance.size规定渐变网格尺寸最宽为100nmgradinglinear规定网格尺寸会从界面处的1nm线性变宽到100nm60~63行对PN结的边界进行细化定义方法与界面的定义方法类似66行启用了施主掺杂浓度的自适应细化max.impurity.gradient2表示当相邻网格点浓度变化大于等于每微米2个数量级时触发细化使网格尺寸减半min.impurity.size0.1表示网格最小细化到0.1um后就不再进行细化68行保存了细化网格后的结构。运行代码后结果如图8所示可以看到界面和PN结处的网格进行了适当的细化使得PN结的掺杂浓度过渡更加流畅这能使计算准确不容易出现收敛性的问题减少了N掺杂内部和P型衬底的网格减少了重复区域的计算量可以加快仿真求解速度的同时不影响结果准确性。图7 通过Victorymesh进行网格细化图8 经过网格细化后的器件结构随后器件的铁电性质进行了仿真仿真代码如图9所示。74行载入在上一个Module中生成的结构文件76行的末尾添加了”lkferro”对氧化层进行铁电模型的指定使用模型为 朗道-哈拉特尼科夫模型 (Landau-Khalatnikov Model简称LK模型) 需要特别说明的是76行的代码必须要存在于求解电滞曲线的Victorydevice的Module内铁电模型才会生效具体而言就本工程来说如果只有结构生成的Module中添加该行代码而求解电滞曲线的Module不添加该行代码则求解过程不会计入LK模型求解结果是错误的但如果把结构生成的代码和数值求解的代码合并在同一个Module中则只在结构生成的代码中添加lkferro就是可以的78~83行添加了铁电模型的参数和栅极功函数作为变量方便后续构建DOE实验。86~92行是本工程的重点在76增加LK模型后进行数值求解的带入就会包含LK铁电模型的默认参数而86~92行的LK模型参数设置则可以自定模型的参数使参数值更合理或更吻合实验数据。lkferro需要通过interface命令生效并设置各个参数的值和lkferro作用的范围。虽然这样的情况下器件从结构上来说和图3所示的典型铁电晶体管结构并不完全相同但是作用原理是相同的仍然可以通过铁电模型的极化方向的变化来影响半导体沟道中的载流子从而影响电流行为。关于LK模型的详细说明可以参阅Silvaco Victory Device User Guide的”3.9.9 Landau-Khalatnikov Ferroelectric Model”此处简要说明模型中各个参数的主要作用(1) lkremanent即剩余极化强度Pr直接影响铁电材料自由能U的极值点的位置更明显的作用是能显著改变电滞曲线的形状(2) lkcoercive即矫顽场强度Ec控制极化的翻转难度(3) lkoffset即E_offsetLK公式中的Etot E_applied E_polarization E_offset , 能够打破正负极化的对称性(4) lkdelta能影响极化突变的速度调节滞回曲线的过渡陡峭程度(5) lkviscosity即LK模型中的粘滞系数ρ能够控制控制极化趋近稳态的速度。再回到代码作用的解读。94~95行使用probe命令监测了lkferro中可以输出的两个物理量ferro_charge和ferro_state分别是极化电荷和极化状态97行设置了栅极的功函数99行添加了求解的模型其中lknc是启用负电容模式(negative capacitance)强制极化始终处于负电容分支根据user guide在瞬态仿真时Silvaco的求解器可以自动确定界面电荷在电容曲线的哪个分支上因此在瞬态仿真时可以不添加lknc如果使用准静态仿真则需要根据需求考虑是否添加该模型101~107通过两次pulse命令设置了两段幅值大小相同符号相反的脉冲用于实现铁电器件的正负极性的转换109行在method部分设置了求解的限制条件111~125行开始进行瞬态求解使铁电器件从t0栅极为0V开始扫描完整进行两次正负扫描并输出log文件和str文件。图9 电滞曲线的数值求解代码关于101~107行的波形设置和119~122行的瞬态仿真设置这里介绍另一种实现办法。在上面展示的代码中我们让器件执行两次正负极性的变化因为器件的栅极是从0V开始变化的在LK模型中这意味着器件开始处于未极化的状态所以可以观察到初始的0~5V扫描和第二次0~5V的扫描的Drain电流不一致即器件完成了从未极化到有极化的变化并且在第二个扫描周期也不会回到未极化的状态。但如果只关心电滞曲线的核心性质也可以采用第二种设置方法让栅极从开始就是-4V使器件在t0时处于负极化的状态这样只需要扫描一个pulse就可以观察到正负极化的变化虽然看不到从未极化到有极化的变化但是也可以观察到核心的电滞曲线可以节省仿真时间。两种设置的结果对比如图10所示。图10 两种波形条件设置下的仿真结果对比最后通过Victory DOE对所有LK模型中的铁电参数和栅极的功函数进行DOE来观察变化规律结果如图11~16所示。如图11提高剩余极化强度Pr会提高栅极为0时的Drain电流使得电滞曲线的形状显著改变窗口更大高Vth和低Vth的差距变大如图12ferro state为2代表没有极化3代表正极化1代表负极化减小矫顽场绿会让ferro_state更早地改变即铁电的极化方向更容易翻转加大矫顽场蓝会让ferro_state更难改变在如图设置的数值(1e7)下ferro_state无法改变如图13加大偏移电场E_offset蓝会让正负极化整体向正栅压方向变化减小偏移电场E_offset绿会让正负极化整体向负栅压方向变化如图14加大lkdelta蓝会让极化状态的突变程度降低减小lkdelta绿会让极化状态的突变程度增加如图15加大粘滞系数lkviscosity绿会让从极化从无到有的过程变得缓慢同时矫顽场Ec变高减小粘滞系数lkviscosity蓝会让从极化从无到有变得更快同时矫顽场Ec变低如图16提高栅极的功函数会让高Vth和低Vth都变高即让电滞曲线整体向正栅压方向移动。初步了解这些规律后如果有类似器件结构的实验结果就可以较为全面的对曲线进行校正以达到仿真数据和实验数据吻合的目的。图11 LKREMANENT的DOE结果Drain Current vs. Gate Voltage图12 LKCOERCIVE的DOE结果Ferro state vs. Gate Voltage图13 LKOFFSET的DOE结果Ferro Charge vs. Gate Voltage图14 LKDELTA的DOE结果Ferro State vs. Gate Voltage图15 LKVISCOSITY的DOE结果Ferro Charge vs. Gate Voltage图16 Workfunction的DOE结果Drain Current vs. Gate Voltage参考周久人韩根全郝跃. 铁电负电容场效应晶体管[M]. 人民邮电出版社, 2025.本文为半导体器件公众号原创内容未经授权禁止转载侵权必究。欢迎加Q技术交流群638099934名称懒小木 Silvaco TCAD交流群 8719137751群3000人已经满了2群973292092 名称 懒小木Sentaurus TCAD交流群有希望有需要学习器件建模TCAD、SPICE的学生、学者或者朋友可以关注“半导体器件”公众号有系统的课程、专业的答疑累计服务行业内4万余人。