计算化学公社

标题: 如何用g09如何获得Bq的受力 [打印本页]

作者
Author: alystone 时间: 2018-1-17 15:46
标题: 如何用g09如何获得Bq的受力
请问各位老师，我可以用g03的massage关键词，得到我结构中Bq的受力。但是在G09中删除了massage的此类功能，我想用punch+iop(7/32)来输出，但是也输出不出来。请教有哪位老师知道怎么在g09中得到Bq的受力情况。g03的输入输出文件见附件，以供参考。

谢谢！

作者
Author: sobereva 时间: 2018-1-17 18:14
g09试了，可以跑，结果和你传的g03一样
massage在g09里只是变成了deprecated，但还能用

作者
Author: alystone 时间: 2018-1-17 18:30

sobereva 发表于 2018-1-17 18:14
g09试了，可以跑，结果和你传的g03一样
massage在g09里只是变成了deprecated，但还能用

请问sob老师用的哪个g09版本，输入文件和我的这个输入文件一模一样吗？

作者
Author: sobereva 时间: 2018-1-17 19:04

alystone 发表于 2018-1-17 18:30
请问sob老师用的哪个g09版本，输入文件和我的这个输入文件一模一样吗？

一模一样
G09W D.01

作者
Author: alystone 时间: 2018-1-17 20:30

sobereva 发表于 2018-1-17 19:04
一模一样
G09W D.01

谢谢sob老师。我想是因为刚才那个体系太简单了，我的体系比较大，然后计算的时候出现这样的错误，应该怎么解决。谢谢！具体错误见附件g09计算出的log文件。不加massage关键词。能够算下去，但是结果是错的。加上massage关键词，就这样报错了。

作者
Author: sobereva 时间: 2018-1-17 21:06

alystone 发表于 2018-1-17 20:30
谢谢sob老师。我想是因为刚才那个体系太简单了，我的体系比较大，然后计算的时候出现这样的错误，应该怎 ...

你的赝势没有被读入，搜索 Pseudopotential Parameters，会看到Fe并没有赋予赝势
把背景电荷和赝势字段交换一下看看
massage我自己倒不用。如果之后还有问题，把背景电荷只加一小点，然后逐渐增加，反复尝试找原因。

作者
Author: alystone 时间: 2018-1-17 21:40

sobereva 发表于 2018-1-17 21:06
你的赝势没有被读入，搜索 Pseudopotential Parameters，会看到Fe并没有赋予赝势
把背景电荷和赝势字 ...

谢谢sob老师。的确如你所说，就是位置不对。但是不管怎么换，好像都不对。
情况1：
坐标
空
基组
空
赝势
空
Bq
这种情况下，Bq读不进去。

情况2：
坐标
空
基组
空
Bq
空
赝势
这种情况下，赝势基组读不进去。

请问sob老师，gaussian可以把基组和赝势信息写在外边，用引用文件的方法来避免这种情况吗。或者还有其他什么方法。再次麻烦您了。谢谢！

作者
Author: sobereva 时间: 2018-1-17 23:57

alystone 发表于 2018-1-17 21:40
谢谢sob老师。的确如你所说，就是位置不对。但是不管怎么换，好像都不对。
情况1：
坐标

引用文件并不解决问题，在计算时只不过会把引用的文件进行展开而已。
关于读取顺序，在这里有详细说明
http://sobereva.com/g09/m_input.htm

作者
Author: alystone 时间: 2018-1-18 08:03

sobereva 发表于 2018-1-17 23:57
引用文件并不解决问题，在计算时只不过会把引用的文件进行展开而已。
关于读取顺序，在这里有详细说明
...

谢谢sob老师。

作者
Author: alystone 时间: 2018-1-18 08:10

sobereva 发表于 2018-1-17 23:57
引用文件并不解决问题，在计算时只不过会把引用的文件进行展开而已。
关于读取顺序，在这里有详细说明
...

sob老师您好，如果我用charge代替massage用来添加背景电荷，是否有方法来输出背景电荷的受力。谢谢

作者
Author: sobereva 时间: 2018-1-18 10:09

alystone 发表于 2018-1-18 08:10
sob老师您好，如果我用charge代替massage用来添加背景电荷，是否有方法来输出背景电荷的受力。谢谢

应该没办法

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