关于orca计算中读取波函数命令求助

lym563557413 · 发表于 Post on 2022-8-6 18:26:31

我想通过读取cc-pvtz基组下算出的波函数，去计算aug-cc-pvtz基组下的单点能，因为直接算atz基组算不出来，然后查阅资料查到了base这个命令，也尝试着写了一下，但是并没有计算出来，因为我的关键词给的不对。请教大家，读取波函数的关键词怎么用。已经产生的波函数的文件格式为.gbw.xz,不知道这个格式的文件还需要处理吗？下面是我自己给出的命令行，但是不对，我贴出来，有懂的老师同学帮忙解答一下，谢谢。

%pal nprocs 24
end
! dlpno-mp2 sp cc-pvtz cc-pvtz/c tightpno base=129-1-tz.gbw.xz
%maxcore 2500
* xyz 0 1

wzkchem5 · 发表于 Post on 2022-8-6 20:53:57

base不是这么用的。orca关键字不像高斯关键字那样是后面加等号的，仔细看说明书了解orca的关键字怎么写。
读取波函数应该用MORead，自己在手册里搜MORead看怎么用

lym563557413 · 发表于 Post on 2022-8-6 22:11:56

wzkchem5 发表于 2022-8-6 20:53
base不是这么用的。orca关键字不像高斯关键字那样是后面加等号的，仔细看说明书了解orca的关键字怎么写。
...

谢谢大神回复，我确实对这个是小白，请问方便发一份说明书吗？我查了好多都没有查到

ghifi37 · 发表于 Post on 2022-8-6 22:39:00

lym563557413 发表于 2022-8-6 22:11
谢谢大神回复，我确实对这个是小白，请问方便发一份说明书吗？我查了好多都没有查到

https://orcaforum.kofo.mpg.de/ap ... etail&df_id=186

lym563557413 · 发表于 Post on 2022-8-8 11:17:26

ghifi37 发表于 2022-8-6 22:39
https://orcaforum.kofo.mpg.de/app.php/dlext/?view=detail&df_id=186

想继续请教老师，我这个现在知道怎么算了，直接在命令行加了MOread,计算也识别到了关键词，但是却报错说找不到波函数，我的inp文件名gbw文件名是一致的，附3张图如下，第一张是报错的输出文件截图，第二张是输入文件命令行截图，第三张是文件名截图。还希望老师给予帮助解答，多谢。

lym563557413 · 发表于 Post on 2022-8-8 11:18:47

3张图

ghifi37 · 发表于 Post on 2022-8-8 15:26:01

lym563557413 发表于 2022-8-8 11:18
3张图

这个gbw文件是怎么获得的？可以用之前产生gbw的级别再读次gbw看看有没问题，还有问题就是文件坏了

ChemG · 发表于 Post on 2022-8-8 16:14:25

用了！MOREAD，你也要告诉人家你的波函数文件是谁啊。%moinp “file.gbw”。

lym563557413 · 发表于 Post on 2022-8-9 20:03:25

ChemG 发表于 2022-8-8 16:14
用了！MOREAD，你也要告诉人家你的波函数文件是谁啊。%moinp “file.gbw”。

多谢老师，但是我不太会用这个命令，我把我的命令粘贴在下面，但是我这个命令不对，请老师帮忙指正一下，辛苦
%pal nprocs 24
end
! dlpno-mp2 sp cc-pvtz cc-pvtz/c tightpno MORead
%moinp “129-1-tz.gbw”
%maxcore 2500
* xyz 0 1

wzkchem5 · 发表于 Post on 2022-8-9 20:13:04

lym563557413 发表于 2022-8-9 13:03
多谢老师，但是我不太会用这个命令，我把我的命令粘贴在下面，但是我这个命令不对，请老师帮忙指正一下， ...

没错，但是引号必须用英文引号，不能用中文引号

lym563557413 · 发表于 Post on 2022-8-23 15:44:45

wzkchem5 发表于 2022-8-9 20:13
没错，但是引号必须用英文引号，不能用中文引号

老师我之前的问题解决掉了一半，就是我的问题是我想通过读取tz基组的波函数来计算atz基组下的单点能，但是成功读取后，却不能正常完成计算，下面我截了一个图，这里面有warning，最后这个题自己就死掉了，请老师指点一下

wzkchem5 · 发表于 Post on 2022-8-23 15:53:04

lym563557413 发表于 2022-8-23 08:44
老师我之前的问题解决掉了一半，就是我的问题是我想通过读取tz基组的波函数来计算atz基组下的单点能，但 ...

图没有传上来

lym563557413 · 发表于 Post on 2022-8-23 16:44:20

老师，是这样子的

lym563557413 · 发表于 Post on 2022-8-23 16:55:34

wzkchem5 发表于 2022-8-23 15:53
图没有传上来

老师，传到下面了

wzkchem5 · 发表于 Post on 2022-8-23 17:31:06

lym563557413 发表于 2022-8-23 09:44
老师，是这样子的

看起来是基组线性相关问题，因为能量远低于合理值。
试试加大格点，并且把%scf里的sthresh参数调大一些（默认值是1e-8，你先试试1e-7，如果还不行就1e-6）。注意同一篇文章里但凡是相同泛函、相同基组的计算，即便直接算能正常收敛的，都要统一调sthresh和格点来和这个计算统一，不然结果不可比

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

[ORCA] 关于orca计算中读取波函数命令求助