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

糖糖糖豆9988 · 发表于 Post on 2022-9-27 15:19:30

sobereva 发表于 2022-9-27 10:08
reverse写IRC里面

老师，再打扰您一下。把reverse和forward分别写在irc里重新计算后，发现两个out文件是一样的，您知道这是怎么回事吗？

wzkchem5 · 发表于 Post on 2022-9-27 15:35:46

忘不掉的回忆发表于 2022-9-27 07:51
有个问题想请教一下老师，
我根据老师的这篇文章（http://bbs.keinsci.com/thread-577-1-1.html）做了两个 ...

（1）很可能是，但是为了完全确定，应该做能量分解分析来确定
（2）如果你只关心能量/自由能，且精度要求不高（比如误差2kcal/mol也可以接受）的话，问题不大；如果对精度要求很高（比如希望误差在1kcal/mol以下），建议做Boltzmann平均；如果要算光谱等，因为能量相差不大的构型的光谱差别可以很大，所以建议对所有构型做Boltzmann平均。
（3）当且仅当计算相互作用能时，计算不同复合物构型所用单体的结构完全一致时可以。因为你考虑了BSSE，所以这样的计算结果会和直接用高精度单点能有一些出入，但我觉得这个出入没有大到让其中一种做法不可信的程度。

wzkchem5 · 发表于 Post on 2022-9-27 15:40:34

dengyh 发表于 2022-9-27 02:40
各位大佬，请教一下这个Orca运行出现“[file orca_mian/run.cpp, line 9893]:Cannot open input file:test. ...

运行orca时，必须先用cd命令，把cmd的当前所在目录切换到test.inp所在目录，才能调用orca。这是因为orca并不默认输入文件和它自己在同一个文件夹，所以如果你没把cmd的当前所在目录切换到test.inp所在目录，那么orca是找不到test.inp的，即使test.inp就在他“眼皮子底下”也看不见。（即便orca说明书可能没有明确说这一点，但是但凡有cmd常识的，不难看出这一点。orca说明书不能负责教用户cmd常识）
此外不建议让test.inp和orca程序处于同一个文件夹，应该在orca程序目录以外的地方新建一个文件夹用来放输入文件。

zcdeng · 发表于 Post on 2022-9-27 19:01:24

wzkchem5 发表于 2022-9-26 00:35
不行。首先HF计算的结果是不能发表的，其次B3LYP只对小体系而言明显比HF慢，对于100原子的体系，B3LYP恐 ...

感谢老师，

忘不掉的回忆 · 发表于 Post on 2022-9-27 19:31:59

wzkchem5 发表于 2022-9-27 15:35
（1）很可能是，但是为了完全确定，应该做能量分解分析来确定
（2）如果你只关心能量/自由能，且精度要 ...

谢谢老师，我是做用量子化学来筛选萃取剂，首先来确定萃取剂与共沸组分的一个稳定构型。
我根据老师的这篇文章（http://sobereva.com/165）做了这五个的构型分布，下图所示，
发现第一个电子能量最低的构型的比例最大，那我将这个构型用来做相互作用能应该可以吧，
谢谢老师

dengyh · 发表于 Post on 2022-9-27 19:32:23

sobereva 发表于 2022-9-27 10:09
你这么运行，test.inp必须放在当前目录下
当前目录是什么在命令行提示符上都写明了

谢谢老师的指点，希望我可以正常应用ORCA程序来解释实际问题

dengyh · 发表于 Post on 2022-9-27 19:33:42

wzkchem5 发表于 2022-9-27 15:40
运行orca时，必须先用cd命令，把cmd的当前所在目录切换到test.inp所在目录，才能调用orca。这是因为orca ...

谢谢老师认真的指导，希望我可以正常应用ORCA程序来解释实际问题，祝你一切顺利

dengyh · 发表于 Post on 2022-9-27 19:34:20

dengyh 发表于 2022-9-27 19:33
谢谢老师认真的指导，希望我可以正常应用ORCA程序来解释实际问题，祝你一切顺利

希望以后不会的多多指教

StormSpirts · 发表于 Post on 2022-9-27 19:51:57

各位老师好。
我在sklearn中定义了LinearRegression模型，想要打印模型参数

model = LinearRegression()
model.fit(X_train_std,y_train)
print(model)

复制代码

输出如下。括号内是空的，print没有打印出模型的具体参数，那么应怎样打印出模型的具体参数呢？

LinearRegression()

复制代码

感谢各位老师。

wzkchem5 · 发表于 Post on 2022-9-27 20:54:56

忘不掉的回忆发表于 2022-9-27 12:31
谢谢老师，我是做用量子化学来筛选萃取剂，首先来确定萃取剂与共沸组分的一个稳定构型。
我根据老师的这 ...

可以，但是如果后续计算发现不同萃取剂的相互作用能差别特别小，必须把相互作用能算得很精确才能解释或预测不同萃取剂之间的差别的话，那么可能需要考虑返回头来做Boltzmann平均

csu-77 · 发表于 Post on 2022-9-27 20:58:30

提交的任务已成功完成，但Maximum Displacement处显示的是No，这个是什么原因导致的呀？出来的数据可信嘛？附件已将计算任务上传

DAP-DAP-OPT-d-p-2.gjf (9.03 KB, 下载次数 Times of downloads: 0)

忘不掉的回忆 · 发表于 Post on 2022-9-27 21:05:24

wzkchem5 发表于 2022-9-27 20:54
可以，但是如果后续计算发现不同萃取剂的相互作用能差别特别小，必须把相互作用能算得很精确才能解释或预 ...

谢谢老师，谢谢老师

wzkchem5 · 发表于 Post on 2022-9-27 21:42:55

底下一行已经写了，因为力已经可以忽略，所以优化完成。
本来收敛的条件就不是四个YES，如果有的YES收敛得极好，可以弥补一或两个NO。因为大多数情况下是四个YES才收敛，就有人讹传说是只有四个YES才是收敛，但凡有一个NO就不是完全收敛，但这个说法其实是错误的。

sobereva · 发表于 Post on 2022-9-27 22:04:42

csu-77 发表于 2022-9-27 20:58
提交的任务已成功完成，但Maximum Displacement处显示的是No，这个是什么原因导致的呀？出来的数据可信嘛？ ...

以后过于初级的问题一律在本帖里问，发主贴时看清楚版头的分类说明，不要乱选“量化理论”分类。而且必须恰当明明标题，不得用诸如“Maximum Displacement No”这种标题，仔细看过http://bbs.keinsci.com/thread-9348-1-1.html后再发帖，否则记违规+删帖+扣分处理

你的问题的回复看16554L

462855 · 发表于 Post on 2022-9-28 09:58:07

求助为什么一个在MS用Dmol3能够正常计算的团簇，转成mol格式后在GaussView中打开会报错？是因为结构太大了吗？

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

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

Maximum Displacement No