应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】g09中连续计算的命令行如何写？
41/1返回列表
查看: 1407  |  回复: 8
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xtdut

银虫 (小有名气)

[交流] 【求助】g09中连续计算的命令行如何写？

在g09的用户手册和例子中经常看到是一个简单的分子（比如水、乙烷等）优化构型，接着计算单点能什么的。我的问题是：这两步计算是共用一个分子说明部分，还是第一步计算的结果作为第二步计算的分子说明部分？如果我想将第一步计算的结果作为第二步计算的分子说明部分，第二步计算结果作为第三步的分子说明部分，第三步计算结果作为第四步计算的分子说明部分，并且每一步计算中的电荷数和自旋多重度要自己定义，每一步的计算结果都要写在输出文件中（也就是.log中）以方便最后使用gview查看结果，命令行应该如何写？同时如何用gview查看每一步的计算结果也是个问题？另外，听人说g09不稳定，这种连续计算要是中间某一步出错，整个计算就会终止，如何避免？希望大家能够耐心点给出命令行的写法及注意事项，不要一个--link1--命令行就了事，其实这里面有很多说道及注意事项的，用户手册和例子中也没有明确地说明，需要大家给出自己的使用心得与经验，谢谢大家！
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼2009-11-29 09:22:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xtdut(金币+2,VIP+0):谢谢您的详细说明，我为自己的不当言辞向您及各位虫友道歉！ 11-30 10:16
erylingjet(金币+8,VIP+0):感谢交流 12-4 10:57
引用回帖:
Originally posted by xtdut at 2009-11-29 09:22:
在g09的用户手册和例子中经常看到是一个简单的分子（比如水、乙烷等）优化构型，接着计算单点能什么的。我的问题是：这两步计算是共用一个分子说明部分，还是第一步计算的结果作为第二步计算的分子说明部分？如果 ...

您说高斯说明书没有说明您的问题，那小卒我就请您看看，高斯说明书是否真的像您所谓的那样，没有对您的问题进行说明？
引用回帖:
你的第一个问题是：

在g09的用户手册和例子中经常看到是一个简单的分子（比如水、乙烷等）优化构型，接着计算单点能什么的。我的问题是：这两步计算是共用一个分子说明部分，还是第一步计算的结果作为第二步计算的分子说明部分？

好了我们先看看高斯说明书中关于link1的解释。在它的解释中，说明书中举得例子如下
引用回帖:
%Chk=freq
# HF/6-31G(d) Freq

Frequencies at STP

Molecule specification

--Link1--
%Chk=freq
%NoSave
# HF/6-31G(d) Geom=Check Guess=Read Freq=(ReadFC,ReadIsotopes)

Frequencies at 300 K
     
charge and spin
300.0  2.0
Isotope specifications

然后我们就要研究一下这个例子
在第二个任务的命令行中，有两个关键字Geom=Check Guess=Read。那么，我们来看一下这两个关键字代表什么

在高斯说明书中
guess关键字的read参数的解释如下
Requests that the initial guess be read from the checkpoint file (Guess=Read is often specified along with Geom=Checkpoint). This option may be combined with Alter, in which case the orbitals are read from the checkpoint file, projected onto the current basis set, and then the specified alterations are made. Checkpoint is a synonym for Read. The TCheck option says to attempt to read a guess from the checkpoint file, but to generate a new one if necessary.
geom关键字的check参数的解释如下：
Checkpoint（这个参数是check的全称）
Causes the molecule specification (including variables) to be taken from the checkpoint file. Only the charge and multiplicity are read from the input stream. For example, Geom=Checkpoint may be used by a later job step to retrieve the geometry optimized during an earlier job step from the checkpoint file. This action is safe since Gaussian will abort the job if an optimization fails, and consequently subsequent job steps which expect to use the optimized geometry will not be executed. May be combined with the ModRedundant option if you want to retrieve and alter the molecule specification in a checkpoint file using redundant internal coordinate-style modifications.
好了，到此为止，高斯说明书已经对您的第一个问题做出了明确回答：既不是“这两步计算是共用一个分子说明部分”，也不是“第一步计算的结果作为第二步计算的分子说明部分”，而是在第二部计算的时候，从第一步计算得到的chk文件中读取分子构型。
引用回帖:
您的第二个问题是：

如果我想将第一步计算的结果作为第二步计算的分子说明部分，第二步计算结果作为第三步的分子说明部分，第三步计算结果作为第四步计算的分子说明部分，并且每一步计算中的电荷数和自旋多重度要自己定义，每一步的计算结果都要写在输出文件中（也就是.log中）以方便最后使用gview查看结果，命令行应该如何写？

我们再来看看高斯说明书是否明确回答了这个问题：
仍然是参看geom关键字的check参数的解释。从这个解释中我们可以很容易的读出它的含义，只要使用了这个参数，那么每一次计算得到的chk文件都会被下一部计算所读取。那么您仍然按照高斯说明书中对link1举的例子，来写您的多步计算任务不就成了？难道还有别的疑问吗？
你所谓的规定不同的自旋多重度，也在高斯说明书中有明确定义。
geom关键字的checkpoint参数的解释中有很明确的一句话Only the charge and multiplicity are read from the input stream.在allcheck这个参数中也有很明确的一句话Causes the molecule specification (including variables), the charge and multiplicity, and the title section to be taken from the checkpoint file.
这样，答案已经很明确了。如果您使用的是check参数，那么电荷和自旋多重度，高斯请您自己输入。如果您使用allcheck参数，那么电荷和自旋多重度，高斯从上一步计算任务中读取。您看，还有比说明书更明确的回答吗？
引用回帖:
您的第三个问题：
听人说g09不稳定，这种连续计算要是中间某一步出错，整个计算就会终止，如何避免？

我还是从说明书的解释中来直接回答您的问题：
在高斯说明书中关于菜单Process Preferences的End Batch Run on Error这个选项的解释如下：
End Batch Run on Error: Toggles whether to halt batch processing when an error occurs, or to skip to the next job in the batch and keep going. (Default ON). Note: If this feature is active and an error occurs while processing a batch, the batch start entry value is set to the file that caused the error.
看完这句话，我还用多说什么？您的问题并不是高斯是否稳定的问题，而是您采用高斯默认参数的问题。为了解决您的问题，您把End Batch Run on Error这个选项勾掉就行了。还有比高斯说明书更明确的解答吗？
引用回帖:
您最后还发了一通感慨：
希望大家能够耐心点给出命令行的写法及注意事项，不要一个--link1--命令行就了事，其实这里面有很多说道及注意事项的，用户手册和例子中也没有明确地说明，需要大家给出自己的使用心得与经验，谢谢大家！

现在您回答小卒我：您现在觉得高斯说明书没有明确给出您的问题答案吗？事实上，高斯说明书是一个非常全面的高斯指南。事无巨细，面面俱到！我们都是做科学的，没有调查就没有发言权。高斯说明书对某某问题的解答就白纸黑字地写在那里，您却说高斯说明书没有说某某问题，您是不是太武断了些？更让我感觉不解的是：您希望网友们“耐心点给出命令行的写法及注意事项”之前，是不是自己应该耐心点读一读说明书呢？中国古语：“己所不欲勿施于人”，您自己不愿意耐心点读一读说明书，就不要怪罪网友们不愿意“耐心点给出命令行的写法及注意事项”了吧。

[ Last edited by yjcmwgk on 2009-11-30 at 14:12 ]
赞一下(12人) 回复此楼
4楼2009-11-29 12:41:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
引用回帖:
Originally posted by zhangmt at 2009-11-29 16:33:
怕机器空着就写批处理啊。
前一段好多例子啊。好多人有回答。
看http://muchong.com/bbs/viewthread.php?tid=1690945&fpage=1

http://muchong.com/bbs/viewthread.php?tid=1566146 其实这个帖子的第一章第六节就是很全面的高斯批处理说明
赞一下 回复此楼
8楼2009-11-29 16:45:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+2,VIP+0):感谢交流 12-4 10:57
引用回帖:
Originally posted by caohuiming at 2009-11-29 16:20:
要时时监控，发现不合适要kill掉程序，按楼主讲的算这末多步，每步之后还不用分析，我是没做过。

我同意这个说法
不过我们的计算任务往往平行的
比如分别优化一个分子的单重态、叁重态和五重太的构型
那么三个opt放在一起做成批处理文件来先后计算，还是不错的哦
其实我也不建议把一堆任务放在一起进行
除非是平行任务

[ Last edited by yjcmwgk on 2009-11-30 at 16:39 ]
赞一下(12人) 回复此楼
9楼2009-11-30 16:38:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xtdut 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭