matlab patch函数_MATLAB实现紧束缚近似能带结构画图

这几天原本做电路实验的时候看到FFT的分析,然后查了一下看到了有用matlab来模拟多少个波的叠加可以得到比较好看的三角波、方波之类的,一时心血来潮,就想着要不我也来试试用matlab画出能带图,于是渣渣我又来了,顺便复习一波matlab画三维图的知识,手册走起。


首先,在利用MATLAB进行画图前,要先对公式进行处理,顺便复习就从头再推一次吧。

以下推导以简单立方晶格s态电子为例,后面有MATLAB代码附上(因为我自己学的不是很扎实,所以有很多小白注释,大佬可以跳过)

简单晶格中电子满足薛定谔方程:

equation?tex=%5B-%5Cfrac%7B%5Chbar%7Dm%5Ctriangledown%5E2%2BU%28r%29%5D%5Cpsi%28r%29%3DE%5Cpsi%28r%29%5C%5C

根据紧束缚近似,

equation?tex=%5Cpsi%28r%29 可以写成布洛赫和的形式:

equation?tex=%5Cpsi%28%7Br%7D%29%3D%5Cfrac%7B1%7D%7B%5Csqrt%7BN%7D%7D%5Csum_%7BR_l%7De%5E%7Bik%5Ccdot+R_l%7D+%5Cvarphi_i%28r-R_l%29%5C%5C

其中

equation?tex=%5Cvarphi_i+ 表示孤立原子波动方程的本征态,代入薛定谔方程中:

equation?tex=%5Csum_%7BR_l%27%7D%5B-%5Cfrac%7B%5Chbar%5E2%7D%7B2m%7D%5Ctriangledown%5E2%2BV%28r%29-E%28k%29%5D%5Cfrac%7B1%7D%7B%5Csqrt%7BN%7D%7De%5E%7Bik%5Ccdot+R_l%27%7D%5Cvarphi_i%28r-R_l%27%29%3D0%5C%5C

其中的周期势

equation?tex=V%28r%29%3D%5Csum_%7BR_l%7DU%28r-R_l%29 为孤立原子势的叠加,利用:

equation?tex=%5B-%5Cfrac%7B%5Chbar%5E2%7D%7B2m%7D%5Ctriangledown%5E2%2BU%28r-R_l%29%5D%5Cvarphi_i%28r-R_l%29%3DE_i%5Cvarphi_i%28r-R_l%29%5C%5C

可以得到:

  • 0
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值