简单量子化学问题答疑专帖

sai77 · 发表于 Post on 2023-2-7 09:29:45

糖糖糖豆9988 发表于 2023-2-7 09:13
请教一下各位老师，对于如图所示的联环丁烷，C8/C9，以及相应C10/C11是不是等效氢呢？谢谢老师

H8/H9和H10/H11吧？明显不等效啊，H8和H10在环上，远离另一个环丁烷；H9和H11在环下，靠近另外一个环丁烷。
H8的等效氢是H5、H18和H20；H10的等效氢是H22。

糖糖糖豆9988 · 发表于 Post on 2023-2-7 10:10:35

sai77 发表于 2023-2-7 09:29
H8/H9和H10/H11吧？明显不等效啊，H8和H10在环上，远离另一个环丁烷；H9和H11在环下，靠近另外一个环丁烷 ...

谢谢您的回复，所以这个结构里一共有五种H原子？

sai77 · 发表于 Post on 2023-2-7 10:28:11

糖糖糖豆9988 发表于 2023-2-7 10:10
谢谢您的回复，所以这个结构里一共有五种H原子？

是的。
另外这个问题跟量化好像没啥关系，基础有机化学里的核磁共振章节会有介绍。

糖糖糖豆9988 · 发表于 Post on 2023-2-7 11:21:42

请教各位老师，利用G4计算，自由能G(T)直接在图示蓝色部分读取吗？

糖糖糖豆9988 · 发表于 Post on 2023-2-7 19:54:10

请教一下各位老师高斯续算的问题，①过渡态续算：#p opt=(calcfc,ts,noeigen,recalc=4,restart) freq 6-311g(d,p) nosymm int=superfine um062x；②IRC续算：#p irc=(rcfc,LQA,maxpoints=100,restart) um062x/6-311G(d,p) nosymm guess=read geom=allcheck；③扫描：#p opt=(modredundant, restart) 6-311g(d,p) m062x nosymm，或者取来最后一帧结构继续按原关键词接着扫描；问题1：前面三种计算任务，输入文件指定%chk=xxx.chk（中断前的.chk文件），然后加上restart关键词后，会自动读取中断前的.chk文件进行续算？问题2：单点能续算，对于CCSD(T)计算如果没有指定%rwf=xxx.rwf文件，是不是就没办法续算呢？我说的有点啰嗦，谢谢老师

sobereva · 发表于 Post on 2023-2-8 05:19:21

糖糖糖豆9988 发表于 2023-2-7 19:54
请教一下各位老师高斯续算的问题，①过渡态续算：#p opt=(calcfc,ts,noeigen,recalc=4,restart) freq 6-311 ...

会
没法续算

sobereva · 发表于 Post on 2023-2-8 05:19:50

糖糖糖豆9988 发表于 2023-2-7 11:21
请教各位老师，利用G4计算，自由能G(T)直接在图示蓝色部分读取吗？

显然是

糖糖糖豆9988 · 发表于 Post on 2023-2-8 09:16:20

本帖最后由糖糖糖豆9988 于 2023-2-8 10:23 编辑

sobereva 发表于 2023-2-8 05:19
会
没法续算

谢谢社长，所以①restart是包含了guess=read功能吗？②续算重新产生了输出文件，是不是要将两次输出文件合并才是完整的计算结果呀？③利用restart续算，.chk文件是不是会覆盖原来的，即.chk文件只包含新的续算部分？如果是的话，对于过渡态搜寻的续算任务，这个新产生的.chk文件是否包含完整的频率信息呢，还可以在下一步guess=read读取.chk文件并且使用rcfc关键词吗？谢谢社长

糖糖糖豆9988 · 发表于 Post on 2023-2-8 11:08:49

请教各位老师，在利用G4方法进行对称破缺态计算的时候，如下的输入文件正确吗？谢谢老师
%chk=G4_G1_stable.chk
%nprocs=8
#p ub3lyp 6-31G(2df,p) guess=mix stable=opt nosymm

2.7

0 1
坐标

--Link1--
%oldchk=G4_G1_stable.chk
%chk=G4_G1.chk
%nprocs=8
#p G4 guess=read geom=allcheck nosymm

chands · 发表于 Post on 2023-2-8 14:47:38

本帖最后由 chands 于 2023-2-8 14:53 编辑

糖糖糖豆9988 发表于 2023-2-8 11:08
请教各位老师，在利用G4方法进行对称破缺态计算的时候，如下的输入文件正确吗？谢谢老师
%chk=G4_G1_stabl ...

有风险。G4是组合方法，包含好几步，第一步是优化。你这么写，G4的第一步优化可能又把你对称破缺的初猜搞没了。但是原理上可以手动把G4的第一步优化改成对称破缺后计算，见http://sobereva.com/419

或者不用G4，后续用合适的级别算单点能。

糖糖糖豆9988 · 发表于 Post on 2023-2-8 15:42:51

chands 发表于 2023-2-8 14:47
有风险。G4是组合方法，包含好几步，第一步是优化。你这么写，G4的第一步优化可能又把你对称破缺的初猜搞 ...

谢谢老师，在输出文件里显示<S**2>= 1.0124，这能说明最后的计算结构仍然符合双自由基特征吧？

chands · 发表于 Post on 2023-2-8 17:16:05

糖糖糖豆9988 发表于 2023-2-8 15:42
谢谢老师，在输出文件里显示= 1.0124，这能说明最后的计算结构仍然符合双自由基特征吧？

是。

护卫天使 · 发表于 Post on 2023-2-8 18:11:10

本帖最后由护卫天使于 2023-2-8 22:24 编辑

各位老师好，想得到图1中NBO信息，用高斯程序添加关键词得到的输出文件，是不用从图2处获得？为什么二阶扰动理论分析会有好几次，每次分成的unit都不一样，unit会有6个？

sobereva · 发表于 Post on 2023-2-9 02:31:57

护卫天使发表于 2023-2-8 18:11
各位老师好，想得到图1中NBO信息，用高斯程序添加关键词得到的输出文件，是不用从图2处获得？为什么二阶扰 ...

对开壳层波函数会对不同自旋分别做分析，往前翻看提示弄清楚当前的分析对象

波函数分析问题应当发到1L提到的专门的答疑专贴里

sobereva · 发表于 Post on 2023-2-9 02:34:11

糖糖糖豆9988 发表于 2023-2-8 09:16
谢谢社长，所以①restart是包含了guess=read功能吗？②续算重新产生了输出文件，是不是要将两次输出文件 ...

1 是

2 你就看最后输出的就完了

3 搞清楚不同任务chk里都记录什么，看具体任务而定
本来opt=TS freq任务就是最终才做的振动分析，显然如果任务算完了，频率必然是完整的

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

[综合交流] 简单量子化学问题答疑专帖