简单量子化学问题答疑专帖（旧版，已关帖）

xjtu_wcl · 发表于 Post on 2016-6-7 15:41:29

本帖最后由 xjtu_wcl 于 2016-6-7 15:45 编辑

Sob大大
1、请问激发态的极化率（只是需要各向同性的平均极化率）可以用
#p polar B3LYP/aug-cc-pVTZ TD density
这样算吗？
2、我计算了Ar的基态和激发态（上面那个路径是默认root=1，对吗？）发现差很大，基态基本没问题11.15，但激发态这样算出来就有248.695，这个值可信吗？还是说要用更好的DFT方法
3、感觉和SOS方法好混乱，我很纠结激发态有没有极化率了。。。（您推荐个文献吧）谢谢大大

大道之行 · 发表于 Post on 2016-6-7 15:41:30

liyuanhe211 发表于 2016-6-6 23:38
搜一下就第一个，很难吗？

老师我有个问题，我的G09方法中没有mpwpw91只有mpw1pw91.那这个mpw1pw91和mpwpw91是一样的吗？

liyuanhe211 · 发表于 Post on 2016-6-7 16:28:02

大道之行发表于 2016-6-7 15:41
老师我有个问题，我的G09方法中没有mpwpw91只有mpw1pw91.那这个mpw1pw91和mpwpw91是一样的吗？

明明有mpwpw91：

liyuanhe211 · 发表于 Post on 2016-6-7 16:37:13

本帖最后由 liyuanhe211 于 2016-6-7 16:39 编辑

大道之行发表于 2016-6-7 15:31
请教一下前辈们，我用opt+freq计算得到的能量（如图1所示）和用energy计算的能量哪一个是单点能，或者他们 ...

这级别算势垒：“MP2/6-31+G(2df,2p)” 真是只有在量子化学考试题里挑错用才会出现的东西。这么几个原子愿意用RI的话，PWPB95/def2-QZVPP//M06-2X/TZVP一点压力都没有，不知道高到哪里去了

零点能和单点能没关系。
你的截图是频率计算部分。单点能在这上面，做MP2相当于最后的图里MP2=的能量。

活化能用TS的单点能和GS的单点能相减；但热力学能垒（用于算速率的活化自由能）参考下文，不精确的对应于Sum of...Free Energies那一行：

谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327

PUHF和PMP2是spin projected HF/MP2，似乎是用于消除自旋污染的，不过你这里投影不投影数值都一样。

大道之行 · 发表于 Post on 2016-6-7 16:42:44

本帖最后由大道之行于 2016-6-7 16:50 编辑

liyuanhe211 发表于 2016-6-7 16:28
明明有mpwpw91：

老师您看

liyuanhe211 · 发表于 Post on 2016-6-7 17:06:27

大道之行发表于 2016-6-7 16:42
老师您看

不要依赖GaussView，自己输入。是在不愿脱离GV可以你先选个B3LYP，保存输出文件，再打开文件把B3LYP改成你需要的东西

大道之行 · 发表于 Post on 2016-6-7 17:24:22

liyuanhe211 发表于 2016-6-7 17:06
不要依赖GaussView，自己输入。是在不愿脱离GV可以你先选个B3LYP，保存输出文件，再打开文件把B3LYP改成 ...

我之前也看过您们的帖子说尽量自己写输入文件但是每次我都不知道该怎么写，并且我试过参照人家的输入文件写，提交计算之后老是报错，然后我就放弃了。非常感谢老师的解答。

liyuanhe211 · 发表于 Post on 2016-6-7 17:44:12

大道之行发表于 2016-6-7 17:24
我之前也看过您们的帖子说尽量自己写输入文件但是每次我都不知道该怎么写，并且我试过参照人家的输入文件 ...

这个故事的正常版本明明应该是“并且我试过参照人家的输入文件写，提交计算之后老是报错，后来我找到了哪里有错误，然后我就学会了。”

大道之行 · 发表于 Post on 2016-6-7 21:14:43

liyuanhe211 发表于 2016-6-7 17:44
这个故事的正常版本明明应该是“并且我试过参照人家的输入文件写，提交计算之后老是报错，后来我找到了哪 ...

问题是我觉得有错的地方我却不知道怎么改

大道之行 · 发表于 Post on 2016-6-7 21:28:22

liyuanhe211 发表于 2016-6-7 17:43
这个故事的正常版本明明应该是“并且我试过参照人家的输入文件写，提交计算之后老是报错，后来我找到了哪 ...

老师，请允许我再问一个问题，最近我在做一个反应机理方面的计算，我想用亲核指数找出反应物的活性位点，然后以各个反应活性位点作为几条反应的起始反应。请问老师这样做能吗？
另外，我其中计算了一个基元反应的过渡态，找到的过渡态有且只有一个虚频，并且也是连接着产物的，但是对该过渡态的IRC计算结果却不如人意。结果如图1和2.这是不是说明我的过渡态不正确呀？

liyuanhe211 · 发表于 Post on 2016-6-7 23:57:25

大道之行发表于 2016-6-7 21:14
问题是我觉得有错的地方我却不知道怎么改

查。实在查不到问。

liyuanhe211 · 发表于 Post on 2016-6-7 23:59:25

大道之行发表于 2016-6-7 21:28
老师，请允许我再问一个问题，最近我在做一个反应机理方面的计算，我想用亲核指数找出反应物的活性位点， ...

过渡态的震动模式的两边应该分别像反应物和底物。IRC也是。如果你的IRC两边确实接近反应物和底物，则没有问题，只不过是你一边的势垒非常小罢了。如果不是，则要调整初猜重新计算。

sobereva · 发表于 Post on 2016-6-8 04:58:20

大道之行发表于 2016-6-7 15:41
老师我有个问题，我的G09方法中没有mpwpw91只有mpw1pw91.那这个mpw1pw91和mpwpw91是一样的吗？

明显不一样，mpw1pw91是杂化泛函
这些破泛函早不流行了，不用考虑

sobereva · 发表于 Post on 2016-6-8 05:00:49

xjtu_wcl 发表于 2016-6-7 15:41
Sob大大
1、请问激发态的极化率（只是需要各向同性的平均极化率）可以用
#p polar B3LYP/aug-cc-pVTZ TD ...

可能不行。你可以自行用有限场来计算alpha和polar输出的alpha对照来验证下。懒得测的话问高斯客服。

大道之行 · 发表于 Post on 2016-6-8 08:54:04

sobereva 发表于 2016-6-8 04:58
明显不一样，mpw1pw91是杂化泛函
这些破泛函早不流行了，不用考虑

好的，非常感谢老师

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 简单量子化学问题答疑专帖（旧版，已关帖）

评分 Rate

评分 Rate

评分 Rate

浏览过的版块