24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1572  |  回复: 17
当前主题已经存档。
本帖产生 1 个 QC强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

forestwolf9291

木虫 (著名写手)

[交流] 【求助】双原子分子势能曲线已有8人参与

简单的双原子分子,用B3LYP方法做势能曲线
黑色为scan命令扫描键长得到的势能曲线
红色为对应于scan里的键长做单点能计算,用各单点能画出的势能曲线

结果为何不一样呢?

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

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

谢谢张老师瞧得起小卒,小卒就大放厥词啦,哈哈

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+10, QC强帖+1):小卒很厉害啊,呵呵 2010-04-12 14:26
引用回帖:
Originally posted by zhangmt at 2010-04-11 11:00:05:
小卒的回答给我一个提示:
假设那是个H2分子,基态自然是单重态,随着两个原子拉开,共价键断裂,
实际上已经是双自由基了——那么基态按说应该是三重态,但是问题在于,gaussian在计算之初就设置了单重态,也不 ...

看法倒是有,未必合适,哈哈
如果是H2的话,直接用cas做几条势能线的计算(反正双原子分子不存在势能面,只有势能线)分别计算在拉断过程中的S0、S1、T1、S2、T2这几条势能线,算完了就知道了。
---(上面的我敢确定)----确定与不确定的分割线---(下面的有猜的成分)---
小卒猜方法:
cas(2,2)肯定很粗糙,而且只能得到S0、S1、T1。即使仅仅用来计算S0、S1、T1这三个态,还是嫌粗糙。
cas(2,4),仍旧嫌粗糙,但是用它计算出来的S2、T2态已经有了一定的可信度了了
cas(2,6),如果是小卒做这个计算,小卒就用这个。注意此时6个活性轨道绝对不是程序自动选定的那种连续的轨道。此时我们应该把pi 轨道们跳过去。
cas(2,8),哈哈,这就过分了
小卒猜基组:
def2-tzvpp
6-311++g(d,2p2d)(因为体系中只有氢,所以第一个加号和括号中的d,纯属占位之用)
aug-cc-pvqz
以上全部看法均针对氢气。如果是别的分子,另当别论。

其实对于S0、S1、T1态来说,cas(2,6)对活性轨道的选择肯定是有些太大了,小卒估计此时高层轨道的占据数不会超过0.01。但是,为了计算更高的激发态,这些牺牲还是值得的。

[ Last edited by yjcmwgk on 2010-4-11 at 16:54 ]
13楼2010-04-11 16:41:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

sunjibo120

银虫 (小有名气)


forestwolf9291(金币+1):谢谢 2010-04-09 14:46
aylayl08(金币+1):感谢讨论 2010-04-13 19:45
scan命令算的一系列能量值和改变键长算的能量值多多少少会有些不同,但是差距不大,红色曲线的情况我也遇到过,可能是你基态没选择好,也有可能是你原子的选择有问题,有些原子不符合势能曲线,你是做什么原子的?
2楼2010-04-09 12:59:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunjibo120

银虫 (小有名气)

forestwolf9291(金币+1):thanks 2010-04-09 14:49
你用MP2方法再试试
3楼2010-04-09 13:00:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gyli

木虫 (小有名气)

~~~~

★ ★ ★
zhou2009(金币+3): 2010-04-09 13:59
forestwolf9291(金币+2):把每个势能曲线的最低点作为能量零点画图,怎么理解?曲线会变吗? 2010-04-09 14:49
引用回帖:
Originally posted by forestwolf9291 at 2010-04-09 12:40:23:
简单的双原子分子,用B3LYP方法做势能曲线
黑色为scan命令扫描键长得到的势能曲线
红色为对应于scan里的键长做单点能计算,用各单点能画出的势能曲线

结果为何不一样呢?

[img]http://pic.muchong.com/2010 ...

scan并不是简单的挨着个算单点
它是用前一个点的结果作为后一个点的初始猜测
若是取点很密
能减轻部分的计算量
你把每个势能曲线的最低点作为能量零点画图试一下
Gravitation cannot be held responsible for people falling in love.
4楼2010-04-09 13:13:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见