✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,
代码获取、论文复现及科研仿真合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
更多Matlab完整代码及仿真定制内容点击👇
🔥 内容介绍
基于平面波展开法实现二维声子晶体能带计算是固体物理学中一个重要的研究课题。声子晶体是一种具有周期性结构的材料,其能带结构对材料的热学、光学和电子性质有着重要影响。因此,准确地计算声子晶体的能带结构对于理解和预测材料的性质具有重要意义。
平面波展开法是计算固体电子结构和声子结构的常用方法之一。它通过将波函数或位移场展开为平面波的线性组合,将固体的晶格周期性结构考虑在内,从而可以较为准确地描述固体的能带结构。在二维声子晶体的能带计算中,平面波展开法可以很好地适用于描述声子的振动特性和能带结构。
在实际的能带计算中,首先需要建立声子晶体的晶格结构模型。这包括晶格常数、原子位置以及晶格的对称性等信息。然后,通过数值计算的方法,可以求解声子的振动模式和频率,进而得到声子的能带结构。平面波展开法在这一过程中可以提供对声子波函数的有效描述,从而得到准确的能带结构信息。
通过基于平面波展开法的二维声子晶体能带计算,可以得到声子在动量空间中的能量分布情况,从而揭示材料的声子性质。这对于材料的热传导、热电性能等方面具有重要意义。同时,通过对声子能带结构的计算和分析,还可以预测材料的光学性质和电子性质,为材料设计和应用提供重要参考。
总之,基于平面波展开法实现二维声子晶体能带计算是一项重要而复杂的研究课题,它为我们深入理解材料的热学、光学和电子性质提供了重要的理论工具。随着计算方法和计算技术的不断发展,相信这一领域的研究将会取得更加深入和准确的成果,为材料科学和工程技术的发展做出更大的贡献。
📣 部分代码
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%平面波展开法
%二维声子晶体带结构计算
%计算二维正方格子
%散射体立于基体之中
%%%%%%%%%%%%%%%%%%%%%%%%%
clear;clc;tic;epssys=1.0e-6; %设定一个最小量,避免系统截断误差或除零错误
%%%%%%%%%%%%%%%%%%%%%%%%%%
%定义实际的正空间格子基矢
%%%%%%%%%%%%%%%%%%%%%%%%%%
a=0.02;
a1=a*[1 0];
a2=a*[0 1];
%%%%%%%%%%%%%%%%%%%%%%%%%%
%定义晶格的参数
%%%%%%%%%%%%%%%%%%%%%%%%%%
rho1=11600;E1=4.08e10;mju1=1.49e10;lambda1=mju1*(E1-2*mju1)/(3*mju1-E1); %散射体的材料参数
rho2=1300;E2=1.175e5;mju2=4e4;lambda2=mju2*(E2-2*mju2)/(3*mju2-E2); %基体的材料参数
Rc=0.006; %散射体截面半径
Ac=pi*(Rc)^2; %散射体截面面积
Au=a^2; %二维格子原胞面积
⛳️ 运行结果
🔗 参考文献
[1] 赵芳.基于光子晶体相似的二维声子晶体特性研究[D].哈尔滨工程大学,2005.DOI:10.7666/d.y780107.
[2] 甄妮.考虑表/界面效应的二维声子晶体能带结构计算与分析[D].北京交通大学,2013.DOI:10.7666/d.Y2603265.