CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 2883  |  回复: 22

IPanzer

铁虫 (初入文坛)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
franch(金币+1): 鼓励交流 2011-04-16 20:47:38
这个问题我在以前用matlab算GNR的TB模型里也遇到过,后来正如LZ所述,将所有本征值的虚数部分去掉就行了~但我一直不是太清楚的是去掉本征值的虚数部分是否对结果有影响?
11楼2011-04-14 10:04:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

★ ★
sunyang1988(金币+2): 谢谢交流 2011-04-16 10:05:43
引用回帖:
Originally posted by IPanzer at 2011-04-14 10:04:43:
这个问题我在以前用matlab算GNR的TB模型里也遇到过,后来正如LZ所述,将所有本征值的虚数部分去掉就行了~但我一直不是太清楚的是去掉本征值的虚数部分是否对结果有影响?

你首先需要确定你的本证值虚数部分是不是为零,是的话就,没问题。在中间你可以把本证值打印出来。确认一下虚数部分是不是零。当然如果你求的本征值虚数不是零,那你的方程就有问题
嘟啦啦嘟啦啦嘟
12楼2011-04-14 14:02:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

IPanzer

铁虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by 贺仪 at 2011-04-14 14:02:25:
你首先需要确定你的本证值虚数部分是不是为零,是的话就,没问题。在中间你可以把本证值打印出来。确认一下虚数部分是不是零。当然如果你求的本征值虚数不是零,那你的方程就有问题

我算出来的本征值的虚数部分都是小数点后3位的量,很接近0,我以前以为是误差造成的~我也考虑的是最近邻近似,直接构建哈密顿矩阵然后求解本征值,实际上就是小木虫上那篇《电子结构计算DIY》里的方法~
13楼2011-04-15 16:52:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)


sunyang1988(金币+1): 谢谢提示 2011-04-16 10:06:02
引用回帖:
Originally posted by IPanzer at 2011-04-15 16:52:17:
我算出来的本征值的虚数部分都是小数点后3位的量,很接近0,我以前以为是误差造成的~我也考虑的是最近邻近似,直接构建哈密顿矩阵然后求解本征值,实际上就是小木虫上那篇《电子结构计算DIY》里的方法~

我算的虚数全为零,这么小可能是你的矩阵的某个元素没弄对吧
嘟啦啦嘟啦啦嘟
14楼2011-04-15 17:10:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2): 谢谢指教 2011-04-16 00:19:05
引用回帖:
Originally posted by IPanzer at 2011-04-15 16:52:17:
我算出来的本征值的虚数部分都是小数点后3位的量,很接近0,我以前以为是误差造成的~我也考虑的是最近邻近似,直接构建哈密顿矩阵然后求解本征值,实际上就是小木虫上那篇《电子结构计算DIY》里的方法~

最好检查下H矩阵。
TB方法得到的H矩阵应该是实对称的
考虑位相因子exp(ikr)后,矩阵是复共轭的,即Hermite矩阵
无论实对称还是Hermite矩阵,其本征值都是实数
15楼2011-04-15 23:30:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
Originally posted by acridine at 2011-04-15 23:30:21:
最好检查下H矩阵。
TB方法得到的H矩阵应该是实对称的
考虑位相因子exp(ikr)后,矩阵是复共轭的,即Hermite矩阵
无论实对称还是Hermite矩阵,其本征值都是实数

前辈你好!我就是在重复别人的结果,单枪匹马,我画是石墨烯的三维能带结构,刚刚才看到你的http://muchong.com/bbs/viewthread.php?tid=1933290,要是在我做之前看到就好了!现在有个问题,希望你能帮忙,你在9楼说:

4.后记
当然除了能级能带之外,使用上述模型还可以分析一些其他东西。如根据特征向量,亦即波函数,可以分析得到一些性质,某个特征波函数主要由哪个原子轨道基函数贡献等等,应该很容易想象到具体方法,这里就不详述了。

我下一步就是求对应一个本证值的的波函数,感觉没头绪。在线性代数里面来说,就是已知本征值怎么求本征向量的问题。为此我们需要选定一个基组,但这个基组应该怎么选呢?或者说产生呢?
嘟啦啦嘟啦啦嘟
16楼2011-04-16 11:03:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夕阳西下

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
franch(金币+1): 鼓励交流 2011-04-16 20:48:47
基组是不是0级波函数时的位置?
17楼2011-04-16 11:43:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
franch(金币+1): 谢谢回帖交流 2011-04-16 20:48:34
引用回帖:
Originally posted by 贺仪 at 2011-04-16 11:03:20:
前辈你好!我就是在重复别人的结果,单枪匹马,我画是石墨烯的三维能带结构,刚刚才看到你的http://muchong.com/bbs/viewthread.php?tid=1933290,要是在我做之前看到就好了!现在有个问题,希望你能 ...

你的这种感觉有点类似“只缘身在此山中”的意思了
你的基组在你选择参数,即矩阵元积分值的时候就已经确定了。而这些矩阵元正是这些基组参与积分后的结果。
详细点说,你的矩阵元Hij对应的就是基组的i第i和j个分量产生的,因此根据各个矩阵元的意义就可以明确基组的具体意义。而你解矩阵得到的特征向量Cn对应的是就是波函数展开成基组时候的系数。
18楼2011-04-16 16:19:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

IPanzer

铁虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by acridine at 2011-04-15 23:30:21:
最好检查下H矩阵。
TB方法得到的H矩阵应该是实对称的
考虑位相因子exp(ikr)后,矩阵是复共轭的,即Hermite矩阵
无论实对称还是Hermite矩阵,其本征值都是实数

非常感谢前辈的回复!
我检查过的,我生成的哈密顿矩阵确实是Hermite矩阵,这是没问题的~所以我很奇怪,怀疑是误差造成的~
19楼2011-04-16 20:10:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
Originally posted by acridine at 2011-04-16 16:19:17:
你的这种感觉有点类似“只缘身在此山中”的意思了
你的基组在你选择参数,即矩阵元积分值的时候就已经确定了。而这些矩阵元正是这些基组参与积分后的结果。
详细点说,你的矩阵元Hij对应的就是基组的i第i和j个 ...

是的 我那个直接取最近邻的为2.7eV,基组就是碳原子的波函数,可是碳原子的波函数具体是什么样子的呢?我求出了Cn后怎么画出来电荷密度呢?谢谢指教!
嘟啦啦嘟啦啦嘟
20楼2011-04-18 18:17:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 贺仪 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见