近邻和四自旋相互作用对铁电子畴结构影响的模拟.pdfVIP

第23卷第3期 青岛大学学报(自然科学版J V01．23No．3 20l 010 0年9月 JOURNALOF Science Sep．2 QINGDAOUNIVERSITY(NaturalEdition) 文章编号：1006—1037(2010)03—0032—06 doi：10．3969／j．issn．1006—1037．2010．03．009 近邻和四自旋相互作用对铁电体畴结构影响的模拟 崔春玲1，朱化凤2 摘要：在二维Ising模型基础上，采用蒙特卡罗方法研究了二自旋和四自旋相互作用下， 铁电体畴结构的动态演化和极化反转过程。数值模拟结果表明二自旋相互作用下极化反 转过程中显示不出明显的畴纵向生长趋势。而四自旋相互作用下，方型的四自旋相互作 用对于畴的纵向生长起主要作用，T型的四自旋相互作用可以导致反铁电相的出现。 关键词：铁电体；畴结构；蒙特卡罗方法 中图分类号：TQl72．66 文献标志码：A 1 问题背景 当铁电材料由顺温相(大多数情形下为高电相)向低温相(铁电相)转变时，晶体将分成若干小区域，每个 小区域内部极化方向相同，但各个小区域中电偶极子方向不同，这就形成铁电畴n矗]。对铁电体畴结构的研 究是认识铁电体极化反转过程及其相关物理性质的重要手段之一，深刻了解畴结构的形成过程及其动态过 程对于控制和调节铁电材料的特性是非常必要的。铁电体畴结构的极化反转一般可分为四个主要过程：畴 的成核、畴的纵向长大、畴的横向扩张和畴合并[3]。目前对铁电畴反转机理的研究还大多停留在一些初步的 探索设想阶段，缺乏进一步的理论分析和研究。 目前对畴结构的模拟主要是在Ising模型基础上进行的[4’5]，模拟结果与实际的畴结构图案有较大的差 别，而且在极化反转过程中不能够反映出畴的纵向长大过程[6-8]。本工作利用计入四自旋相互作用的Ising 模型格林函数方法进行计算机模拟，直观地显示铁电体的畴结构随外场电压的变化图样，从而对铁电体畴结 构的反转机理进行了一些理论分析和实验研究。这一模型所得结果可以与实验结果定量地符合，而且在足 够大的四自旋相互作用下可以描述一级相变过程。 2 理论模型 子只允许处在两种状态中的任意一种[6]。外场下的Ising模型在仅考虑近邻、次近邻和第三近邻相互作用 时，对于二维自旋网格L×L的系统哈密顿量为 1一 一 H5一专∑J,jS,s；一∑I『』Es7 。 0 i S；，S分别代表第i和『格点上原子或原子团的自旋。其中．，F为交换积分，ep--自旋相互作用系数，它 包括近邻、次近邻和第三近邻相互作用。为方便模拟计算，取二自旋相互作用系数．，#1。岸为自旋偶极矩，E 为外加电场，计算中外加电场E(t)一EoSin(oJt)，它与所应用的频率硼=o．01成正弦变化关系。 ·收稿日期：2010—01—25；修回日期：2010一06—20 课题 ． 作者简介：崔春玲(1973一)，女，讲师。研究方向；铣电体性质理论研究． 万方数据 第3期 崔春玲，等：近邻和四自旋相互作用对铁电体畴结构影响的模拟 33 旋为S，一S。，次近邻四个格点上的自旋为S。--Ss，第三近邻四个格点为S。一s·z。 用Ising模型。其作用系数均为J。