计算化学公社

标题: 如何从多步计算的chk文件中提取前几步的波函数 [打印本页]

作者
Author: liyuanhe211 时间: 2015-10-7 03:19
标题: 如何从多步计算的chk文件中提取前几步的波函数
假想采用下面的计算路径（简写），实现

%chk = abc.chk
#p opt freq B3LYP/6-31G(d)

--Link1--
%chk = abc.chk
#p freq geom=allcheck int=ultrafine M06-2X/ma-def2-tzvp(假设已定义)

得出的abc.chk文件打开后将是最后一步计算得出的波函数，但考虑到部分波函数分析方法不能接受带弥散函数的波函数，所以希望使用第一步计算得出的chk文件，可以在不重新做freq的情况下得到吗？

作者
Author: sobereva 时间: 2015-10-7 06:56
第二步的时候用%oldchk=abc.chk、%chk=abc2.chk，这样做完第二步后第一步的abc.chk不受影响。

作者
Author: liyuanhe211 时间: 2015-10-7 15:01

sobereva 发表于 2015-10-7 06:56
第二步的时候用%oldchk=abc.chk、%chk=abc2.chk，这样做完第二步后第一步的abc.chk不受影响。

翻译的中文版的G09手册上竟然没有这一内容，简直坑(⊙﹏⊙)b

作者
Author: sobereva 时间: 2015-10-7 16:08

liyuanhe211 发表于 2015-10-7 15:01
翻译的中文版的G09手册上竟然没有这一内容，简直坑(⊙﹏⊙)b

大抵是D.01才有的。得到官网看最新版手册。
中文翻译版对应的是A.02或者B.01版。

作者
Author: 一条君 时间: 2020-7-4 19:53

sobereva 发表于 2015-10-7 06:56
第二步的时候用%oldchk=abc.chk、%chk=abc2.chk，这样做完第二步后第一步的abc.chk不受影响。

请问老师%oldchk方法（1）并没有产生新的chk文件，想看第二步的chk文件怎么办？（2）还导致整个文件夹原来能打开的chk文件出现了can't find gaussian formchk utility，如何解决呢
gjf文件如下，核数与内存与原chk不同，是这个原因吗
原gjf
%nproc=2
%mem=800MB
%chk=E:\Calculate\EC.chk
#p B3LYP/6-311++G** opt freq
接续gjf
%nproc=3
%mem=1000MB
%oldchk=E:\Calculate\EC.chk
#P G4MP2 geom=allcheck

作者
Author: hebrewsnabla 时间: 2020-7-4 20:24

一条君发表于 2020-7-4 19:53
请问老师%oldchk方法（1）并没有产生新的chk文件，想看第二步的chk文件怎么办？（2）还导致整个文件夹原 ...

第二步这样试试
%oldchk=xxx
%chk=yyy
新的chk是yyy

不行的话可能是版本问题。

作者
Author: sobereva 时间: 2020-7-5 01:50

一条君发表于 2020-7-4 19:53
请问老师%oldchk方法（1）并没有产生新的chk文件，想看第二步的chk文件怎么办？（2）还导致整个文件夹原 ...

你第二步没写%chk当然不产生新的chk文件
你的这个计算也显得莫名其妙，本身G4MP2自动就会做opt freq何必还自行做第一步

作者
Author: 一条君 时间: 2020-7-5 15:28
本帖最后由一条君于 2020-7-7 10:44 编辑

sobereva 发表于 2020-7-5 01:50
你第二步没写%chk当然不产生新的chk文件
你的这个计算也显得莫名其妙，本身G4MP2自动就会做opt freq何必 ...

b3lyp是以前练手算的，那请问老师有b3lyp++g**优化的结果下，G4MP2用未优化结构的gjf与b3lyp优化后的结构计算，哪个耗时更短呢（应该是oldchk读取优化后结构吧），谢谢老师

作者
Author: sobereva 时间: 2020-7-6 11:59

一条君发表于 2020-7-5 15:28
第一步是以前练手算的，那请问老师有第一步的结果下，G4MP2用初始还gjf还是oldchk，哪个耗时更短呢（应该 ...

看不懂你的意思
如果只从chk里读取结构，和gjf里直接写结构没有差别

欢迎光临计算化学公社 (http://bbs.keinsci.com/)