Originally posted by IPanzer at 2011-04-14 10:04:43:
这个问题我在以前用matlab算GNR的TB模型里也遇到过，后来正如LZ所述，将所有本征值的虚数部分去掉就行了~但我一直不是太清楚的是去掉本征值的虚数部分是否对结果有影响？

Originally posted by 贺仪 at 2011-04-14 14:02:25:
你首先需要确定你的本证值虚数部分是不是为零，是的话就，没问题。在中间你可以把本证值打印出来。确认一下虚数部分是不是零。当然如果你求的本征值虚数不是零，那你的方程就有问题

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

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

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

Originally posted by 贺仪 at 2011-04-16 11:03:20:
前辈你好！我就是在重复别人的结果，单枪匹马，我画是石墨烯的三维能带结构，刚刚才看到你的http://muchong.com/bbs/viewthread.php?tid=1933290，要是在我做之前看到就好了！现在有个问题，希望你能 ...

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

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