24小时热门版块排行榜    

CyRhmU.jpeg
查看: 950  |  回复: 3

kathy2008

木虫 (正式写手)

[交流] 【求助】如何读进去一个12*12矩阵的本征矢量。已有3人参与

本人用huckel 法算一个分子,得到一个多烯烃和三个乙烯分子的本征值和本征矢量,现不知如何读人12*12矩阵的本征矢量(多烯烃)。2*2矩阵的本征矢量(乙烯)。请高手指点。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
nono2009(金币+1):鼓励应助 2010-11-11 07:43:38
引用回帖:
Originally posted by kathy2008 at 2010-11-10 19:59:20:
本人用huckel 法算一个分子,得到一个多烯烃和三个乙烯分子的本征值和本征矢量,现不知如何读人12*12矩阵的本征矢量(多烯烃)。2*2矩阵的本征矢量(乙烯)。请高手指点。

选用基组后,构造出H矩阵,对角化得到你的本征矢量,不就行了嘛?
你要读入干什么?
2楼2010-11-10 21:04:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ykwang

金虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
nono2009(金币+2):鼓励应助 2010-11-12 08:35:22
假定本征矢量是按本征值由低到高存放在文件中,且程序中已声明采用双精度数组A(12,12), B(2,2)来存放要读入的本征矢量,则可用Open语句打开文件(通道号设为8),按下列示例代码把它们读进来:

READ(8,*) ((A(I,J),I=1,12),J=1,12)
READ(8,*) ((B(I,J),I=1,2),J=1,2)
Nothing_Is_Impossible!
3楼2010-11-11 08:00:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kathy2008

木虫 (正式写手)

谢谢

nono2009:建议通使用“引用回复该帖”,不然别人收不到你的回复提示。 2010-11-12 08:35:38
谢谢高手指点,主要是通过这些本征矢量来求轨道混合系数。非常感谢。
4楼2010-11-11 08:40:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kathy2008 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见