应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 热点前沿问题讨论 » 【其他】讨论：紧束缚模型（TB）可以用于各种体系吗？计算的结果怎么样？
北京石油化工学院2026年研究生招生接收调剂公告
232/3返回列表上一页123下一页
查看: 3111  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

IPanzer

铁虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
franch(金币+1): 鼓励交流 2011-04-16 20:47:38
这个问题我在以前用matlab算GNR的TB模型里也遇到过，后来正如LZ所述，将所有本征值的虚数部分去掉就行了~但我一直不是太清楚的是去掉本征值的虚数部分是否对结果有影响？
赞一下(2人) 回复此楼
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所述，将所有本征值的虚数部分去掉就行了~但我一直不是太清楚的是去掉本征值的虚数部分是否对结果有影响？

你首先需要确定你的本证值虚数部分是不是为零，是的话就，没问题。在中间你可以把本证值打印出来。确认一下虚数部分是不是零。当然如果你求的本征值虚数不是零，那你的方程就有问题
赞一下(1人) 回复此楼
嘟啦啦嘟啦啦嘟
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》里的方法~
赞一下(1人) 回复此楼
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》里的方法~

我算的虚数全为零，这么小可能是你的矩阵的某个元素没弄对吧
赞一下(1人) 回复此楼
嘟啦啦嘟啦啦嘟
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矩阵，其本征值都是实数
赞一下(2人) 回复此楼
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级波函数时的位置？
赞一下(2人) 回复此楼
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对应的是就是波函数展开成基组时候的系数。
赞一下(2人) 回复此楼
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矩阵，这是没问题的~所以我很奇怪，怀疑是误差造成的~
赞一下(1人) 回复此楼
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的回帖
相关版块跳转 我要订阅楼主 贺仪 的主题更新
232/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +11 111623 2026-04-04 11/550 2026-04-05 11:38 by lqwchd
[考研] 348求调剂 +5 wukira 2026-04-04 5/250 2026-04-05 11:05 by xiayan13521
[考研] 求调剂 +3 电气小神童 2026-04-04 3/150 2026-04-05 10:17 by barlinike
[考研] 359求调剂22408 +3 123456789qw 2026-03-31 3/150 2026-04-05 10:09 by zhq0425
[考研] 调剂求助 +10 想换手机不想解� 2026-04-02 13/650 2026-04-05 09:41 by sam3303
[考研] 求生物学调剂 +14 15172915737 2026-04-01 14/700 2026-04-04 20:13 by babysonlkd
[考研] 321求调剂 +6 认真求上学 2026-04-03 6/300 2026-04-04 19:51 by dongzh2009
[考研] 求调剂 +4 压力??大 2026-04-03 4/200 2026-04-03 21:36 by 啵啵啵0119
[考研] 336求调剂 +8 kiyy 2026-04-01 8/400 2026-04-03 19:41 by lijunpoly
[考研] 求调剂 +4 15064154688 2026-04-03 5/250 2026-04-03 15:07 by zrongyan
[考研] 320求调剂 +5 振—TZ 2026-04-02 5/250 2026-04-03 14:42 by fxue1114
[考研] 319求调剂 +18 太容易1018 2026-04-01 18/900 2026-04-03 11:18 by linyelide
[考研] 交通运输考试264分求工科调剂 +4 jike777 2026-04-02 4/200 2026-04-02 21:53 by zllcz
[考研] 一志愿北京科技材料科学与工程288分，求调剂 +14 是辰啊 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
[考研] 07生物学求调剂 一志愿同济大学359分 +3 LAMC. 2026-03-30 3/150 2026-04-02 10:26 by 18828373951
[考研] 354求调剂 +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[考研] 08生物与医药专硕初试346找调剂 +6 dianeeee 2026-04-01 7/350 2026-04-02 08:23 by guoweigw
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3 kaobao456 2026-03-29 4/200 2026-03-30 23:16 by nBu锂
信息提示
关闭
请填处理意见
关闭